Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences

Nicotinamide phosphoribosyl transferase in mammary gland epithelial cells is required for nicotinamide mononucleotide production in mouse milk

Tissue-specific deficiency of nicotinamide phosphoribosyl transferase (NAMPT), the rate-limiting enzyme of the nicotinamide adenine dinucleotide (NAD+)-salvage pathway, causes a decrease of NAD+ in the tissue, resulting in functional abnormalities. The NAD+-salvage pathway is drastically activated in the mammary gland during lactation, but the significance of this has not been established. To investigate the impact of NAD+ perturbation in the mammary gland, we generated two new lines of mammary gland epithelial-cell-specific Nampt-knockout mice (MGKO). LC-MS/MS analyses confirmed that the levels of NAD+ and its precursor nicotinamide mononucleotide (NMN) were significantly increased in lactating mammary glands. We found that murine milk contained a remarkably high level of NMN. MGKO exhibited a significant decrease in tissue NAD+ and milk NMN levels in the mammary gland during lactation periods. Despite the decline in NAD+ levels, the mammary glands of MGKO appeared to develop normally. Transcriptome analysis revealed that the gene profiles of MGKO were indistinguishable from those of their wild-type counterparts, except for Nampt. Although the NMN levels in milk from MGKO were decreased, the metabolomic profile of milk was otherwise unaltered. The mammary gland also contains adipocytes, but adipocyte-specific deficiency of Nampt did not affect mammary gland NAD+ metabolism or mammary gland development. These results demonstrate that the NAD+ -salvage pathway is activated in mammary epithelial cells during lactation and suggest that this activation is required for production of milk NMN rather than mammary gland development. Our MGKO mice could be a suitable model for exploring the potential roles of NMN in milk.

An association between Dnmt1 and Wnt in the production of oocytes in the whitefly Bemisia tabaci

The function of DNA methylation in insects and the DNA methyltransferase (Dnmt) genes that influence methylation remains uncertain. We used RNA interference to reduce the gene expression of Dnmt1 within the whitefly Bemisia tabaci (Hemiptera:Aleyrodidae; Gennadius), a hemipteran species that relies on Dnmt1 for proper gametogenesis. We then used RNA-seq to test an a priori hypothesis that meiosis-related genetic pathways would be perturbed. We generally did not find an overall effect on meiosis-related pathways. However, we found that genes in the Wnt pathway, genes associated with the entry into meiosis in vertebrates, were differentially expressed. Our results are consistent with Dnmt1 knockdown influencing specific pathways and not causing general transcriptional response. This is a finding that is also seen with other insect species. We also characterised the methylome of B. tabaci and assessed the influence of Dnmt1 knockdown on cytosine methylation. This species has methylome characteristics comparable to other hemipterans regarding overall level, enrichment within gene bodies, and a bimodal distribution of methylated/non-methylated genes. Very little differential methylation was observed, and difference in methylation were not associated with differences in gene expression. The effect on Wnt presents an interesting new candidate pathway for future studies.

lncRNA-gene network analysis reveals the effects of early maternal nutrition on mineral homeostasis and energy metabolism in the fetal liver transcriptome of beef heifers

Maternal nutrition during pregnancy influences fetal development; however, the regulatory markers of fetal programming across different gestational phases remain underexplored in livestock models. Herein, we investigated the regulatory role of long non-coding RNAs (lncRNAs) on fetal liver gene expression, the impacts of maternal vitamin and mineral supplementation, and the rate of maternal body weight gain during the periconceptual period. To this end, crossbred Angus heifers (n=31) were randomly assigned to a 2×2 factorial design to evaluate the main effects of the rate of weight gain (low gain [LG, avg. daily gain of 0.28 kg/day] vs. moderate gain [MG, avg. daily gain of 0.79 kg/day]) and vitamins and minerals supplementation (VTM vs. NoVTM). On day 83±0.27 of gestation, fetuses were collected for morphometric measurements, and fetal liver was collected for transcriptomic and mineral analyses. The maternal diet significantly affected fetal liver development and mineral reserves. Using an RNA-Seq approach, we identified 320 unique differentially expressed genes (DEGs) across all six comparisons (FDR <0.05). Furthermore, lncRNAs were predicted through the FEELnc pipeline, revealing 99 unique differentially expressed lncRNAs (DELs). The over-represented pathways and biological processes (BPs) were associated with energy metabolism, Wnt signaling, CoA carboxylase activity, and fatty acid metabolism. The DEL-regulated BPs were associated with metal ion transport, pyrimidine metabolism, and classical energy metabolism-related glycolytic, gluconeogenic, and TCA cycle pathways. Our findings suggest that lncRNAs regulate mineral homeostasis- and energy metabolism-related gene networks in the fetal liver in response to early maternal nutrition.

Host response to Aplysia Abyssovirus 1 in nervous system and gill

The California sea hare (Aplysia californica) is a model for age associated cognitive decline. Recent researched identified a novel nidovirus, Aplysia Abyssovirus 1, with broad tropism enriched in the Aplysia nervous system. This virus is ubiquitous in wild and maricultured, young and old animals without obvious pathology. Here we re-evaluated gene expression data from several previous studies to investigate differential expression in the nervous system and gill in response to virus and aging as well as the mutational spectrum observed in the viral sequences obtained from these datasets. Viral load and age were highly correlated, indicating persistent infection. Upregulated genes in response to virus were enriched for immune genes and signatures of ER and proteostatic stress, while downregulated genes were enriched for mitochondrial metabolism. Differential expression with respect to age suggested increased iron accumulation and decreased glycolysis, fatty acid metabolism, and proteasome function. Interaction of gene expression trends associated with viral infection and aging suggest that viral infection likely plays a role in aging in the Aplysia nervous system. Mutation analysis of viral RNA identified signatures suggesting ADAR and AID/APOBEC like deaminase act as part of Aplysia anti-viral defense.

1,25-dihydroxyvitamin D3 augments low-dose PMA-based monocyte-to-macrophage differentiation in THP-1 cells

The human monocytic THP-1 cell line is the most routinely employed in vitro model for studying monocyte-to-macrophage differentiation. Despite the wide use of this model, differentiation protocols using phorbol 12-myristate-13-acetate (PMA) or 1,25-dihydroxyvitamin D3 (1,25D3) vary drastically between studies. Given that differences in differentiation protocols have the potential to impact the characteristics of the macrophages produced, we aimed to assess the efficacy of three different THP-1 differentiation protocols by assessing changes in morphology and gene- and cell surface macrophage marker expression. THP-1 cells were differentiated with either 5 nM PMA, 10 nM 1,25D3, or a combination thereof, followed by a rest period. The results indicated that all three protocols significantly increased the expression of the macrophage markers, CD11b (p < 0.001) and CD14 (p < 0.010). Despite this, THP-1 cells exposed to 1,25D3 alone did not adopt the morphological and expression characteristics associated with macrophages. PMA was required to produce these characteristics, which were found to be more pronounced in the presence of 1,25D3. Both PMA- and PMA with 1,25D3-differentiated THP-1 cells were capable of M1 and M2 macrophage polarization, though the gene expression of polarization-associated markers was most pronounced in PMA with 1,25D3-differentiated THP-1 cells. Moreover, the combination of PMA with 1,25D3 appeared to support the process of commitment to a particular polarization state.

Glycolysis inhibition affects proliferation and cytotoxicity of Vγ9Vδ2 T cells expanded for adoptive cell therapy

BACKGROUND AIMS

Vγ9Vδ2 T cells are under investigation as alternative effector cells for adoptive cell therapy (ACT) in cancer. Despite promising in vitro results, anti-tumor efficacies in early clinical studies have been lower than expected, which could be ascribed to the complex interplay of tumor and immune cell metabolism competing for the same nutrients in the tumor microenvironment.

METHODS

To contribute to the scarce knowledge regarding gamma delta T-cell metabolism, we investigated the metabolic phenotype of 25-day-expanded Vγ9Vδ2 T cells and how it is intertwined with functionality.

RESULTS

We found that Vγ9Vδ2 T cells displayed a quiescent metabolism, utilizing both glycolysis and oxidative phosphorylation (OXPHOS) for energy production, as measured in Seahorse assays. Upon T-cell receptor activation, both pathways were upregulated, and inhibition with metabolic inhibitors showed that Vγ9Vδ2 T cells were dependent on glycolysis and the pentose phosphate pathway for proliferation. The dependency on glucose for proliferation was confirmed in glucose-free conditions. Cytotoxicity against malignant melanoma was reduced by glycolysis inhibition but not OXPHOS inhibition.

CONCLUSIONS

These findings lay the groundwork for further studies on manipulation of Vγ9Vδ2 T-cell metabolism for improved ACT outcome.

An optimized workflow for transcriptomic analysis from archival paraformaldehyde-fixed retinal tissues collected by laser capture microdissection

RNA sequencing (RNA-seq) coupled with laser capture microdissection (LCM) is a powerful tool for transcriptomic analysis in unfixed fresh-frozen tissues. Fixation of ocular tissues for immunohistochemistry commonly involves the use of paraformaldehyde (PFA) followed by embedding in Optimal Cutting Temperature (OCT) medium for long-term cryopreservation. However, the quality of RNA derived from such archival PFA-fixed/OCT-embedded samples is often compromised, limiting its suitability for transcriptomic studies. In this study, we aimed to develop a methodology to extract high-quality RNA from PFA-fixed canine eyes by utilizing LCM to isolate retinal tissue. We demonstrate the efficacy of an optimized LCM and RNA purification protocol for transcriptomic profiling of PFA-fixed retinal specimens. We compared four pairs of canine retinal tissues, where one eye was subjected to PFA-fixation prior to OCT embedding, while the contralateral eye was embedded fresh frozen (FF) in OCT without fixation. Since the RNA obtained from PFA-fixed retinas were contaminated with genomic DNA, we employed two rounds of DNase I treatment to obtain RNA suitable for RNA-seq. Notably, the quality of sequencing reads and gene sets identified from both PFA-fixed and FF tissues were nearly identical. In summary, our study introduces an optimized workflow for transcriptomic profiling from PFA-fixed archival retina. This refined methodology paves the way for improved transcriptomic analysis of preserved ocular tissue, bridging the gap between optimal sample preservation and high-quality RNA data acquisition.

MCT1-dependent lactate recycling is a metabolic vulnerability in colorectal cancer cells upon acquired resistance to anti-EGFR targeted therapy

Despite the implementation of personalized medicine, patients with metastatic CRC (mCRC) still have a dismal overall survival due to the frequent occurrence of acquired resistance mechanisms thereby leading to clinical relapse. Understanding molecular mechanisms that support acquired resistance to anti-EGFR targeted therapy in mCRC is therefore clinically relevant and key to improving patient outcomes. Here, we observe distinct metabolic changes between cetuximab-resistant CRC cell populations, with in particular an increased glycolytic activity in KRAS-mutant cetuximab-resistant CRC cells (LIM1215 and OXCO2) but not in KRAS-amplified resistant DiFi cells. We show that cetuximab-resistant LIM1215 and OXCO2 cells have the capacity to recycle glycolysis-derived lactate to sustain their growth capacity. This is associated with an upregulation of the lactate importer MCT1 at both transcript and protein levels. Pharmacological inhibition of MCT1, with AR-C155858, reduces the uptake and oxidation of lactate and impairs growth capacity in cetuximab-resistant LIM1215 cells both in vitro and in vivo. This study identifies MCT1-dependent lactate utilization as a clinically actionable, metabolic vulnerability to overcome KRAS-mutant-mediated acquired resistance to anti-EGFR therapy in CRC.

Stem cells tightly regulate dead cell clearance to maintain tissue fitness

Billions of cells are eliminated daily from our bodies1-4. Although macrophages and dendritic cells are dedicated to migrating and engulfing dying cells and debris, many epithelial and mesenchymal tissue cells can digest nearby apoptotic corpses1-4. How these non-motile, non-professional phagocytes sense and eliminate dying cells while maintaining their normal tissue functions is unclear. Here we explore the mechanisms that underlie their multifunctionality by exploiting the cyclical bouts of tissue regeneration and degeneration during hair cycling. We show that hair follicle stem cells transiently unleash phagocytosis at the correct time and place through local molecular triggers that depend on both lipids released by neighbouring apoptotic corpses and retinoids released by healthy counterparts. We trace the heart of this dual ligand requirement to RARγ-RXRα, whose activation enables tight regulation of apoptotic cell clearance genes and provides an effective, tunable mechanism to offset phagocytic duties against the primary stem cell function of preserving tissue integrity during homeostasis. Finally, we provide functional evidence that hair follicle stem cell-mediated phagocytosis is not simply redundant with professional phagocytes but rather has clear benefits to tissue fitness. Our findings have broad implications for other non-motile tissue stem or progenitor cells that encounter cell death in an immune-privileged niche.

Therapeutic potential of human microglia transplantation in a chimeric model of CSF1R-related leukoencephalopathy

Microglia replacement strategies are increasingly being considered for the treatment of primary microgliopathies like adult-onset leukoencephalopathy with axonal spheroids and pigmented glia (ALSP). However, available mouse models fail to recapitulate the diverse neuropathologies and reduced microglia numbers observed in patients. In this study, we generated a xenotolerant mouse model lacking the fms-intronic regulatory element (FIRE) enhancer within Csf1r, which develops nearly all the hallmark pathologies associated with ALSP. Remarkably, transplantation of human induced pluripotent stem cell (iPSC)-derived microglial (iMG) progenitors restores a homeostatic microglial signature and prevents the development of axonal spheroids, white matter abnormalities, reactive astrocytosis, and brain calcifications. Furthermore, transplantation of CRISPR-corrected ALSP-patient-derived iMG reverses pre-existing spheroids, astrogliosis, and calcification pathologies. Together with the accompanying study by Munro and colleagues, our results demonstrate the utility of FIRE mice to model ALSP and provide compelling evidence that iMG transplantation could offer a promising new therapeutic strategy for ALSP and perhaps other microglia-associated neurological disorders.

Transcriptomic landscapes reveal development-related physiological processes in the two-spotted spider mite, Tetranychus urticae

The two-spotted spider mite (Tetranychus urticae Koch, TSSM) is recognized as one of the most problematic spider mite pests. However, the precise gene expression patterns across its key developmental stages remain elusive. Here, we performed a comprehensive transcriptome analysis of TSSM eggs, nymphs and adult females using publicly available RNA sequencing (RNA-seq) data to elucidate the overarching transcriptomic differences between these developmental stages. Principal component analysis and hierarchical clustering analysis unveiled distinct separations among samples across different developmental stages, regardless of their Wolbachia infection status. Differential expression analysis revealed 4,089,2,762, and 1,282 core genes specifically enriched in eggs, nymphs, and adults, respectively. KEGG and GO enrichment analyses showed upregulation of genes in eggs are associated with proteolysis, Wnt signaling pathway, DNA transcription, RNA biosynthetic and metabolic processes, as well as protein folding, sorting, and degradation pathways. Meanwhile, nymphs exhibited increased abundance of genes related to chitin/amino sugar metabolic processes, G protein-coupled receptor signaling pathways, monoatomic ion transport, and neurotransmitter transport pathways. Pathways involving sphingolipid and carbohydrate metabolic processes, proteolysis, lipid transport, and localization were particularly enriched in older females. Altogether, our findings suggest that the egg stage exhibits higher activity in cell differentiation processes, the nymph stage is more involved in chitin development, and the adult stage shows increased metabolic and reproductive activity. This study enhances our understanding of the molecular mechanisms underlying TSSM development and paves the way for further research into the intricate physiological processes of TSSM.

Implications of noncoding regulatory functions in the development of insulinomas

Insulinomas are rare neuroendocrine tumors arising from pancreatic β cells, characterized by aberrant proliferation and altered insulin secretion, leading to glucose homeostasis failure. With the aim of uncovering the role of noncoding regulatory regions and their aberrations in the development of these tumors, we coupled epigenetic and transcriptome profiling with whole-genome sequencing. As a result, we unraveled somatic mutations associated with changes in regulatory functions. Critically, these regions impact insulin secretion, tumor development, and epigenetic modifying genes, including polycomb complex components. Chromatin remodeling is apparent in insulinoma-selective domains shared across patients, containing a specific set of regulatory sequences dominated by the SOX17 binding motif. Moreover, many of these regions are H3K27me3 repressed in β cells, suggesting that tumoral transition involves derepression of polycomb-targeted domains. Our work provides a compendium of aberrant cis-regulatory elements affecting the function and fate of β cells in their progression to insulinomas and a framework to identify coding and noncoding driver mutations.

Cell-type-specific loops linked to RNA polymerase II elongation in human neural differentiation

DNA is folded into higher-order structures that shape and are shaped by genome function. The role of long-range loops in the establishment of new gene expression patterns during cell fate transitions remains poorly understood. Here, we investigate the link between cell-specific loops and RNA polymerase II (RNA Pol II) during neural lineage commitment. We find thousands of loops decommissioned or gained de novo upon differentiation of human induced pluripotent stem cells (hiPSCs) to neural progenitor cells (NPCs) and post-mitotic neurons. During hiPSC-to-NPC and NPC-to-neuron transitions, genes changing from RNA Pol II initiation to elongation are >4-fold more likely to anchor cell-specific loops than repressed genes. Elongated genes exhibit significant mRNA upregulation when connected in cell-specific promoter-enhancer loops but not invariant promoter-enhancer loops or promoter-promoter loops or when unlooped. Genes transitioning from repression to RNA Pol II initiation exhibit a slight mRNA increase independent of loop status. Our data link cell-specific loops and robust RNA Pol II-mediated elongation during neural cell fate transitions.

New potential diagnostic markers for verrucous hyperplasia and verrucous carcinoma based on RNA-sequencing data

Verrucous carcinoma (VC) is a rare subtype of squamous cell carcinoma (SCC) characterized by its histological presentation as a low-grade tumor with no potential for metastasis, setting it apart from invasive SCC. However, distinguishing VC from its benign counterpart, verrucous hyperplasia (VH), is challenging due to their clinical and morphological similarities. Despite the importance of accurate diagnosis for determining treatment strategies, diagnosis for of VH and VC relied only on lesion recurrence after resection. To address this challenge, we generated RNA profiling data from tissue samples of VH and VC patients to identify novel diagnostic markers. We analyzed differentially expressed (DE) mRNA and long non-coding RNA (lncRNA) in tissue samples from VH and VC patients. Additionally, ChIP-X Enrichment Analysis 3 (ChEA3) was conducted to identify the top five transcription factors potentially regulating the expression of DE mRNAs in VH and VC. Our analysis of mRNA and lncRNA expression profiles in VH and VC provides insights into the underlying molecular characteristics of these diseases and offers potential new diagnostic markers. The identification of specific DE genes and lncRNAs may enable clinicians to more accurately differentiate between VH and VC, leading to better treatment choices.

Comparative analyses suggest a link between mRNA splicing, stability, and RNA covalent modifications in flowering plants

BACKGROUND

In recent years, covalent modifications on RNA nucleotides have emerged as pivotal moieties influencing the structure, function, and regulatory processes of RNA Polymerase II transcripts such as mRNAs and lncRNAs. However, our understanding of their biological roles and whether these roles are conserved across eukaryotes remains limited.

RESULTS

In this study, we leveraged standard polyadenylation-enriched RNA-sequencing data to identify and characterize RNA modifications that introduce base-pairing errors into cDNA reads. Our investigation incorporated data from three Poaceae (Zea mays, Sorghum bicolor, and Setaria italica), as well as publicly available data from a range of stress and genetic contexts in Sorghum and Arabidopsis thaliana. We uncovered a strong enrichment of RNA covalent modifications (RCMs) deposited on a conserved core set of nuclear mRNAs involved in photosynthesis and translation across these species. However, the cohort of modified transcripts changed based on environmental context and developmental program, a pattern that was also conserved across flowering plants. We determined that RCMs can partly explain accession-level differences in drought tolerance in Sorghum, with stress-associated genes receiving a higher level of RCMs in a drought tolerant accession. To address function, we determined that RCMs are significantly enriched near exon junctions within coding regions, suggesting an association with splicing. Intriguingly, we found that these base-pair disrupting RCMs are associated with stable mRNAs, are highly correlated with protein abundance, and thus likely associated with facilitating translation.

CONCLUSIONS

Our data point to a conserved role for RCMs in mRNA stability and translation across the flowering plant lineage.

Novel Insights into Post-Myocardial Infarction Cardiac Remodeling through Algorithmic Detection of Cell-Type Composition Shifts

Background

Recent advances in single cell sequencing have led to an increased focus on the role of cell-type composition in phenotypic presentation and disease progression. Cell-type composition research in the heart is challenging due to large, frequently multinucleated cardiomyocytes that preclude most single cell approaches from obtaining accurate measurements of cell composition. Our in silico studies reveal that ignoring cell type composition when calculating differentially expressed genes (DEGs) can have significant consequences. For example, a relatively small change in cell abundance of only 10% can result in over 25% of DEGs being false positives.

Methods

We have implemented an algorithmic approach that uses snRNAseq datasets as a reference to accurately calculate cell type compositions from bulk RNAseq datasets through robust data cleaning, gene selection, and multi-sample cross-subject and cross-cell-type deconvolution. We applied our approach to cardiomyocyte-specific α1A adrenergic receptor (CM-α1A-AR) knockout mice. 8-12 week-old mice (either WT or CM-α1A-KO) were subjected to permanent left coronary artery (LCA) ligation or sham surgery (n=4 per group). Transcriptomes from the infarct border zones were collected 3 days later and analyzed using our algorithm to determine cell-type abundances, corrected differential expression calculations using DESeq2, and validated these findings using RNAscope.

Results

Uncorrected DEGs for the CM-α1A-KO X LCA interaction term featured many cell-type specific genes such as Timp4 (fibroblasts) and Aplnr (cardiomyocytes) and overall GO enrichment for terms pertaining to cardiomyocyte differentiation (P=3.1E-4). Using our algorithm, we observe a striking loss of cardiomyocytes and gain in fibroblasts in the α1A-KO + LCA mice that was not recapitulated in WT + LCA animals, although we did observe a similar increase in macrophage abundance in both conditions. This recapitulates prior results that showed a much more severe heart failure phenotype in CM-α1A-KO + LCA mice. Following correction for cell-type, our DEGs now highlight a novel set of genes enriched for GO terms such as cardiac contraction (P=3.7E-5) and actin filament organization (P=6.3E-5).

Conclusions

Our algorithm identifies and corrects for cell-type abundance in bulk RNAseq datasets opening new avenues for research on novel genes and pathways as well as an improved understanding of the role of cardiac cell types in cardiovascular disease.

RNAseq analysis of oocyte maturation from the germinal vesicle stage to metaphase II in pig and human

During maturation oocytes at the germinal vesicle (GV) stage progress to metaphase II (MII). However, during in vitro maturation a proportion often fail to progress. To understand these processes, we employed RNA sequencing to examine the transcriptome profile of these three groups of oocytes from the pig. We compared our findings with similar public oocyte data from humans. The transcriptomes in oocytes that failed to progress was similar to those that did. We found in both species, the most upregulated genes in MII oocytes were associated with chromosome segregation and cell cycle processes, while the most down regulated genes were relevant to ribosomal and mitochondrial pathways. Moreover, those genes involved in chromosome segregation during GV to MII transition were conserved in pig and human. We also compared MII and GV oocyte transcriptomes at the isoform transcript level in both species. Several thousands of genes (including DTNBP1, MAPK1, RAB35, GOLGA7, ATP1A1 and ATP2B1) identified as not different in expression at a gene transcript level were found to have differences in isoform transcript levels. Many of these genes were involved in ATPase-dependent or GTPase-dependent intracellular transport in pig and human, respectively. In conclusion, our study suggests the failure to progress to MII in vitro may not be regulated at the level of the genome and that many genes are differentially regulated at the isoform level, particular those involved ATPase- or GTPase-dependent intracellular transport.

Transcriptional and epigenetic characterization of a new in vitro platform to model the formation of human pharyngeal endoderm

BACKGROUND

The Pharyngeal Endoderm (PE) is an extremely relevant developmental tissue, serving as the progenitor for the esophagus, parathyroids, thyroids, lungs, and thymus. While several studies have highlighted the importance of PE cells, a detailed transcriptional and epigenetic characterization of this important developmental stage is still missing, especially in humans, due to technical and ethical constraints pertaining to its early formation.

RESULTS

Here we fill this knowledge gap by developing an in vitro protocol for the derivation of PE-like cells from human Embryonic Stem Cells (hESCs) and by providing an integrated multi-omics characterization. Our PE-like cells robustly express PE markers and are transcriptionally homogenous and similar to in vivo mouse PE cells. In addition, we define their epigenetic landscape and dynamic changes in response to Retinoic Acid by combining ATAC-Seq and ChIP-Seq of histone modifications. The integration of multiple high-throughput datasets leads to the identification of new putative regulatory regions and to the inference of a Retinoic Acid-centered transcription factor network orchestrating the development of PE-like cells.

CONCLUSIONS

By combining hESCs differentiation with computational genomics, our work reveals the epigenetic dynamics that occur during human PE differentiation, providing a solid resource and foundation for research focused on the development of PE derivatives and the modeling of their developmental defects in genetic syndromes.

Unique features of KGN granulosa-like tumour cells in the regulation of steroidogenic and antioxidant genes

The ovarian KGN granulosa-like tumour cell line is commonly used as a model for human granulosa cells, especially since it produces steroid hormones. To explore this further, we identified genes that were differentially expressed by KGN cells compared to primary human granulosa cells using three public RNA sequence datasets. Of significance, we identified that the expression of the antioxidant gene TXNRD1 (thioredoxin reductase 1) was extremely high in KGN cells. This is ominous since cytochrome P450 enzymes leak electrons and produce reactive oxygen species during the biosynthesis of steroid hormones. Gene Ontology (GO) analysis identified steroid biosynthetic and cholesterol metabolic processes were more active in primary granulosa cells, whilst in KGN cells, DNA processing, chromosome segregation and kinetochore pathways were more prominent. Expression of cytochrome P450 cholesterol side-chain cleavage (CYP11A1) and cytochrome P450 aromatase (CYP19A1), which are important for the biosynthesis of the steroid hormones progesterone and oestrogen, plus their electron transport chain members (FDXR, FDX1, POR) were measured in cultured KGN cells. KGN cells were treated with 1 mM dibutyryl cAMP (dbcAMP) or 10 μM forskolin, with or without siRNA knockdown of TXNRD1. We also examined expression of antioxidant genes, H2O2 production by Amplex Red assay and DNA damage by γH2Ax staining. Significant increases in CYP11A1 and CYP19A1 were observed by either dbcAMP or forskolin treatments. However, no significant changes in H2O2 levels or DNA damage were found. Knockdown of expression of TXNRD1 by siRNA blocked the stimulation of expression of CYP11A1 and CYP19A1 by dbcAMP. Thus, with TXNRD1 playing such a pivotal role in steroidogenesis in the KGN cells and it being so highly overexpressed, we conclude that KGN cells might not be the most appropriate model of primary granulosa cells for studying the interplay between ovarian steroidogenesis, reactive oxygen species and antioxidants.

Common lizard primary oviduct cell culture: A model system for the genetic and cellular basis of oviparity and viviparity

Reproduction by egg-laying (oviparity) or live-bearing (viviparity) is a genetically determined trait fundamental to the biology of amniotes. Squamates are an emerging model for the genetics of reproductive mode yet lack cell culture models valuable for exploring molecular mechanisms. Here, we report a novel primary culture model for reproductive biology: cell cultures derived from the oviduct tissues (infundibulum, uterus and vagina) of oviparous and viviparous common lizards (Lacertidae: Zootoca vivipara). We maintained and expanded these cultures for over 100 days, including repeated subculturing and successful revival of cryopreserved cells. Immunocytochemical investigation suggested expression of both epithelial and fibroblast-like proteins, and RNA sequencing of cultured cells as compared to in vivo oviduct tissue showed changes in gene expression in response to the cell culture environment. Despite this, we confirmed the maintenance of distinct gene expression patterns in viviparous and oviparous cells after 60+ days of cell culture, finding 354 differentially expressed genes between viviparous and oviparous cells. Furthermore, we confirmed the expression of 15 viviparity-associated candidate genes in cells maintained for 60+ days in culture. Our study demonstrates the feasibility and utility of oviduct cell culture for molecular analysis of reproductive mode and provides a tool for future genetic experiments.

Angiopoietin-like 4 protects against endothelial dysfunction during bacterial sepsis

Loss of endothelial integrity and vascular leakage are central features of sepsis pathogenesis; however, no effective therapeutic mechanisms for preserving endothelial integrity are available. Here we show that, compared to dermal microvessels, brain microvessels resist infection by Neisseria meningitidis, a bacterial pathogen that causes sepsis and meningitis. By comparing the transcriptional responses to infection in dermal and brain endothelial cells, we identified angiopoietin-like 4 as a key factor produced by the brain endothelium that preserves blood-brain barrier integrity during bacterial sepsis. Conversely, angiopoietin-like 4 is produced at lower levels in the peripheral endothelium. Treatment with recombinant angiopoietin-like 4 reduced vascular leakage, organ failure and death in mouse models of lethal sepsis and N. meningitidis infection. Protection was conferred by a previously uncharacterized domain of angiopoietin-like 4, through binding to the heparan proteoglycan, syndecan-4. These findings reveal a potential strategy to prevent endothelial dysfunction and improve outcomes in patients with sepsis.

Use of equine embryo -derived mesenchymal stromal cells and their extracellular vesicles as a treatment for persistent breeding-induced endometritis in susceptible mares

Persistent breeding induced endometritis (PBIE) is a significant cause of infertility in mares. The development of a safe, universal, readily available therapeutic to manage PBIE and facilitate an optimal uterine environment for embryo development may improve pregnancy rates in susceptible mares. Mesenchymal stromal cells (MSCs) are being used increasingly as a therapeutic mediator for inflammatory conditions such as endometritis, and early gestational tissue provides a unique source of multipotent stem cells for creating MSCs. Extracellular vesicles (EVs) are mediators of cell communication produced by many different cell types. This study utilized embryo-derived mesenchymal stromal cells (EDMSCs) and their EVs as a potential therapeutic modality for PBIE in two groups: a) PBIE-susceptible mares challenged with pooled dead sperm (n=5); and b) client-owned mares diagnosed as susceptible to PBIE (n=37 mares and 40 estrous cycles). Mares pre-treated with intrauterine EDMSCs or their EVs resulted in a significant reduction in the accumulation of intrauterine fluid post-breeding. Nine of 19 (47 %) mares treated with EDMSCs prior to natural breeding and 13 of 20 (65 %) mares treated with EDMSC derived EVs were pregnant after the first cycle and 12 of 18 (67 %) mares treated with EDMSCs, and 15 of 19 (79 %) mares treated with EVs conceived by the end of the breeding season. These preliminary clinical studies are the first reports of the use of EDMSCs or their EVs as a potential intrauterine therapy for the management of PBIE susceptible mares.

SERPINA3 is a marker of cartilage differentiation and is essential for the expression of extracellular matrix genes during early chondrogenesis

Serine proteinase inhibitors (serpins) are a family of structurally similar proteins which regulate many diverse biological processes from blood coagulation to extracellular matrix (ECM) remodelling. Chondrogenesis involves the condensation and differentiation of mesenchymal stem cells (MSCs) into chondrocytes which occurs during early development. Here, and for the first time, we demonstrate that one serpin, SERPINA3 (gene name SERPINA3, protein also known as alpha-1 antichymotrypsin), plays a critical role in chondrogenic differentiation. We observed that SERPINA3 expression was markedly induced at early time points during in vitro chondrogenesis. We examined the expression of SERPINA3 in human cartilage development, identifying significant enrichment of SERPINA3 in developing foetal cartilage compared to total limb, which correlated with well-described markers of cartilage differentiation. When SERPINA3 was silenced using siRNA, cartilage pellets were smaller and contained lower proteoglycan as determined by dimethyl methylene blue assay (DMMB) and safranin-O staining. Consistent with this, RNA sequencing revealed significant downregulation of genes associated with cartilage ECM formation perturbing chondrogenesis. Conversely, SERPINA3 silencing had a negligible effect on the gene expression profile during osteogenesis suggesting the role of SERPINA3 is specific to chondrocyte differentiation. The global effect on cartilage formation led us to investigate the effect of SERPINA3 silencing on the master transcriptional regulator of chondrogenesis, SOX9. Indeed, we observed that SOX9 protein levels were markedly reduced at early time points suggesting a role for SERPINA3 in regulating SOX9 expression and activity. In summary, our data support a non-redundant role for SERPINA3 in enabling chondrogenesis via regulation of SOX9 levels.

Tumor-Intrinsic Activity of Chromobox 2 Remodels the Tumor Microenvironment in High-grade Serous Carcinoma

Chromobox 2 (CBX2), an epigenetic reader and component of polycomb repressor complex 1, is highly expressed in >75% of high-grade serous carcinoma. Increased CBX2 expression is associated with poorer survival, whereas CBX2 knockdown leads to improved chemotherapy sensitivity. In a high-grade serous carcinoma immune-competent murine model, knockdown of CBX2 decreased tumor progression. We sought to explore the impact of modulation of CBX2 on the tumor immune microenvironment (TIME), understanding that the TIME plays a critical role in disease progression and development of therapy resistance. Exploration of existing datasets demonstrated that elevated CBX2 expression significantly correlated with specific immune cell types in the TIME. RNA sequencing and pathway analysis of differentially expressed genes demonstrated immune signature enrichment. Confocal microscopy and co-culture experiments found that modulation of CBX2 leads to increased recruitment and infiltration of macrophages. Flow cytometry of macrophages cultured with CBX2-overexpressing cells showed increased M2-like macrophages and decreased phagocytosis activity. Cbx2 knockdown in the Trp53-null, Brca2-null ID8 syngeneic murine model (ID8 Trp53-/-Brca2-/-) led to decreased tumor progression compared with the control. NanoString immuno-oncology panel analysis suggested that knockdown in Cbx2 shifts immune cell composition, with an increase in macrophages. Multispectral immunohistochemistry (mIHC) further confirmed an increase in macrophage infiltration. Increased CBX2 expression leads to recruitment and polarization of protumor macrophages, and targeting CBX2 may serve to modulate the TIME to enhance the efficacy of immune therapies.

SIGNIFICANCE

CBX2 expression correlates with the TIME. CBX2 modulation shifts the macrophage population, potentially leading to an immunosuppressive microenvironment, highlighting CBX2 as a target to improve efficacy of immunotherapy.

Molecular basis of the unique osmoregulatory strategy in the inshore hagfish, Eptatretus burgeri

Hagfishes are characterized by omo- and iono-conforming nature similar to marine invertebrates. Conventionally, hagfishes had been recognized as the most primitive living vertebrate that retains plesiomorphic features. However, some of the "ancestral" features of hagfishes, such as rudimentary eyes and the lack of vertebrae, have been proven to be deceptive. Similarly, by the principle of maximum parsimony, the unique body fluid regulatory strategy of hagfishes seems to be apomorphic, since the lamprey, another cyclostome, adopts osmo- and iono-regulatory mechanisms as in jawed vertebrates. Although hagfishes are unequivocally important in discussing the origin and evolution of the vertebrate osmoregulatory system, the molecular basis for the body fluid homeostasis in hagfishes has been poorly understood. In the present study, we explored this matter in the inshore hagfish, Eptatretus burgeri, by analyzing the transcriptomes obtained from the gill, kidney, and muscle of the animals acclimated to distinct environmental salinities. Together with the measurement of parameters in the muscular fluid compartment, our data indicate that the hagfish possesses an ability to conduct free amino acid (FAA)-based osmoregulation at a cellular level, which is in coordination with the renal and branchial FAA absorption. We also revealed that the hagfish does possess the orthologs of the known osmoregulatory genes and that the transepithelial movement of inorganic ions in the hagfish gill and kidney is more complex than previously thought. These observations pose a challenge to the conventional view that the physiological features of hagfishes have been inherited from the last common ancestor of the extant vertebrates.

Neurexin 3 is Essential for the Specific Wiring of a Color Pathway in the Mammalian Retina

Precise wiring within sensory systems is critical for the accurate transmission of information. In the visual system, S-cone photoreceptors specialize in detecting short-wavelength light, crucial to color perception and environmental cue detection. S-cones form specific synapses with S-cone bipolar cells (SCBCs), a connection that is remarkably consistent across species. Yet, the molecular mechanisms guiding this specificity remain unexplored. To address this, we used the cone-dominant ground squirrel for deep-sequencing of cone subtype transcriptomes and identified Nrxn3 as an essential molecule for the S-cone to SCBC synapse. Using transgenic mouse models, we further examined the role of Nrxn3 in S-cones and discovered a significant reduction of SCBC connections in the absence of Nrxn3. This finding extends the known functions of neurexins, typically associated with synapse regulation, by highlighting their essential role in a specific synaptic connection for the first time. Moreover, the differentially expressed genes identified here pave the way for further investigations into the unique functions of cone subtypes.

The PAH1-encoded phosphatidate phosphatase of Yarrowia lipolytica differentially affects gene expression and lipid biosynthesis

Yarrowia lipolytica is a model oleaginous yeast with a strong capacity for lipid accumulation, yet its lipid metabolic pathways and regulatory mechanisms remain largely unexplored. The PAH1-encoded phosphatidate (PA) phosphatase governs lipid biosynthesis by its enzymatic activity and regulating the transcription of genes involved in phospholipid biosynthesis. In this work, we examined the effect of the loss of Pah1 (i.e., pah1Δ) on cell metabolism in cells growing in low- and high-glucose media. Multi-omics analyses revealed the global effect of the pah1Δ mutation on lipid and central carbon metabolism. Lipidomics analyses showed that the pah1Δ mutation caused a massive decrease in the masses of triacylglycerol (TAG) and diacylglycerol (DAG), and these effects were independent of glucose concentration in the media. Conversely, phospholipid levels declined in low-glucose media but increased in high-glucose media. The loss of Pah1 affected the expression of genes involved in key pathways of glucose metabolism, such as glycolysis, citric acid cycle, oxidative phosphorylation, and the pentose phosphate pathway, and these effects were more pronounced in high-glucose media. In lipid biosynthesis, the genes catalyzing phosphatidylcholine (PC) synthesis from phosphatidylethanolamine (PE) were upregulated within the CDP-DAG pathway. In contrast, PC synthesis through the Kennedy pathway was downregulated. The ethanolamine branch of the Kennedy pathway that synthesizes PE was also upregulated in pah1Δ. Interestingly, we noted a massive increase in the levels of lysophospholipids, consistent with the upregulation of genes involved in lipid turnover. Overall, this work identified novel regulatory roles of Pah1 in lipid biosynthesis and gene expression.

Hierarchical lncRNA regulatory network in early-onset severe preeclampsia

BACKGROUND

Recent studies have shown that several long non-coding RNAs (lncRNAs) in the placenta are associated with preeclampsia (PE). However, the extent to which lncRNAs may contribute to the pathological progression of PE is unclear.

RESULTS

Here, we report a hierarchical regulatory network involved in early-onset severe PE (EOSPE). We have carried out transcriptome sequencing on the placentae from patients and normal subjects to identify the differentially expressed genes (DEGs), including some lncRNAs (DElncRNAs). We then constructed a high-quality hierarchical regulatory network of lncRNAs, transcription factors (TFs), and target DEGs, containing 1851 lncRNA-TF interactions and 6901 TF-promoter interactions. The lncRNA-to-target regulatory interactions were further validated by the triplex structures between the DElncRNAs and the promoters of the target DEGs. The DElncRNAs in the regulatory network were clustered into 3 clusters, one containing DElncRNAs correlated with the blood pressure, including FLNB-AS1 with targeting 27.89% (869/3116) DEGs in EOSPE. We further demonstrated that FLNB-AS1 could bind the transcription factor JUNB to regulate a series members of the HIF-1 signaling pathway in trophoblast cells.

CONCLUSIONS

Our results suggest that the differential expression of lncRNAs may perturb the lncRNA-TF-DEG hierarchical regulatory network, leading to the dysregulation of many genes involved in EOSPE. Our study provides a new strategy and a valuable resource for studying the mechanism underlying gene dysregulation in EOSPE patients.

High-throughput sensitive screening of small molecule modulators of microexon alternative splicing using dual Nano and Firefly luciferase reporters

Disruption of alternative splicing frequently causes or contributes to human diseases and disorders. Consequently, there is a need for efficient and sensitive reporter assays capable of screening chemical libraries for compounds with efficacy in modulating important splicing events. Here, we describe a screening workflow employing dual Nano and Firefly luciferase alternative splicing reporters that affords efficient, sensitive, and linear detection of small molecule responses. Applying this system to a screen of ~95,000 small molecules identified compounds that stimulate or repress the splicing of neuronal microexons, a class of alternative exons often disrupted in autism and activated in neuroendocrine cancers. One of these compounds rescues the splicing of several analyzed microexons in the cerebral cortex of an autism mouse model haploinsufficient for Srrm4, a major activator of brain microexons. We thus describe a broadly applicable high-throughput screening system for identifying candidate splicing therapeutics, and a resource of small molecule modulators of microexons with potential for further development in correcting aberrant splicing patterns linked to human disorders and disease.

The splicing factor CCAR1 regulates the Fanconi anemia/BRCA pathway

The twenty-three Fanconi anemia (FA) proteins cooperate in the FA/BRCA pathway to repair DNA interstrand cross-links (ICLs). The cell division cycle and apoptosis regulator 1 (CCAR1) protein is also a regulator of ICL repair, though its possible function in the FA/BRCA pathway remains unknown. Here, we demonstrate that CCAR1 plays a unique upstream role in the FA/BRCA pathway and is required for FANCA protein expression in human cells. Interestingly, CCAR1 co-immunoprecipitates with FANCA pre-mRNA and is required for FANCA mRNA processing. Loss of CCAR1 results in retention of a poison exon in the FANCA transcript, thereby leading to reduced FANCA protein expression. A unique domain of CCAR1, the EF hand domain, is required for interaction with the U2AF heterodimer of the spliceosome and for excision of the poison exon. Taken together, CCAR1 is a splicing modulator required for normal splicing of the FANCA mRNA and other mRNAs involved in various cellular pathways.

Association analysis of indel variants and gene expression identifies MDM4 as a novel locus for skeletal muscle hypertrophy and power athlete status

Insertions and deletions (indels) are the second most common type of variation in the human genome. However, limited data on their associations with exercise-related phenotypes have been documented. The aim of the present study was to examine the association between 18,370 indel variants and power athlete status, followed by additional studies in 357,246 individuals. In the discovery phase, the D allele of the MDM4 gene rs35493922 I/D polymorphism was over-represented in power athletes compared with control subjects (P = 7.8 × 10-9) and endurance athletes (P = 0.0012). These findings were replicated in independent cohorts, showing a higher D allele frequency in power athletes compared with control subjects (P = 0.016) and endurance athletes (P = 0.031). Furthermore, the D allele was positively associated (P = 0.0013) with greater fat-free mass in the UK Biobank. MDM4 encodes a protein that inhibits the activity of p53, which induces muscle fibre atrophy. Accordingly, we found that MDM4 expression was significantly higher in the vastus lateralis of power athletes compared with endurance athletes (P = 0.0009) and was positively correlated with the percentage of fast-twitch muscle fibres (P = 0.0062) and the relative area occupied by fast-twitch muscle fibres (P = 0.0086). The association between MDM4 gene expression and an increased proportion of fast-twitch muscle fibres was confirmed in two additional cohorts. Finally, we found that the MDM4 DD genotype was associated with increased MDM4 gene expression in vastus lateralis and greater cross-sectional area of fast-twitch muscle fibres. In conclusion, MDM4 is suggested to be a potential regulator of muscle fibre specification and size, with its indel variant being associated with power athlete status. HIGHLIGHTS: What is the central question of this study? Which indel variants are functional and associated with sport- and exercise-related traits? What is the main finding and its importance? Out of 18,370 tested indels, the MDM4 gene rs35493922 I/D polymorphism was found to be the functional variant (affecting gene expression) and the most significant, with the deletion allele showing associations with power athlete status, fat-free mass and cross-sectional area of fast-twitch muscle fibres. Furthermore, the expression of MDM4 was positively correlated with the percentage of fast-twitch muscle fibres and the relative area occupied by fast-twitch muscle fibres.

Endomucin marks quiescent long-term multi-lineage repopulating hematopoietic stem cells and is essential for their transendothelial migration

Endomucin (EMCN) currently represents the only hematopoietic stem cell (HSC) marker expressed by both murine and human HSCs. Here, we report that EMCN+ long-term repopulating HSCs (LT-HSCs; CD150+CD48-LSK) have a higher long-term multi-lineage repopulating capacity compared to EMCN- LT-HSCs. Cell cycle analyses and transcriptional profiling demonstrated that EMCN+ LT-HSCs were more quiescent compared to EMCN- LT-HSCs. Emcn-/- and Emcn+/+ mice displayed comparable steady-state hematopoiesis, as well as frequencies, transcriptional programs, and long-term multi-lineage repopulating capacity of their LT-HSCs. Complementary functional analyses further revealed increased cell cycle entry upon treatment with 5-fluorouracil and reduced granulocyte colony-stimulating factor (GCSF) mobilization of Emcn-/- LT-HSCs, demonstrating that EMCN expression by LT-HSCs associates with quiescence in response to hematopoietic stress and is indispensable for effective LT-HSC mobilization. Transplantation of wild-type bone marrow cells into Emcn-/- or Emcn+/+ recipients demonstrated that EMCN is essential for endothelial cell-dependent maintenance/self-renewal of the LT-HSC pool and sustained blood cell production post-transplant.

A real-world multi-center RNA-seq benchmarking study using the Quartet and MAQC reference materials

Translating RNA-seq into clinical diagnostics requires ensuring the reliability and cross-laboratory consistency of detecting clinically relevant subtle differential expressions, such as those between different disease subtypes or stages. As part of the Quartet project, we present an RNA-seq benchmarking study across 45 laboratories using the Quartet and MAQC reference samples spiked with ERCC controls. Based on multiple types of 'ground truth', we systematically assess the real-world RNA-seq performance and investigate the influencing factors involved in 26 experimental processes and 140 bioinformatics pipelines. Here we show greater inter-laboratory variations in detecting subtle differential expressions among the Quartet samples. Experimental factors including mRNA enrichment and strandedness, and each bioinformatics step, emerge as primary sources of variations in gene expression. We underscore the profound influence of experimental execution, and provide best practice recommendations for experimental designs, strategies for filtering low-expression genes, and the optimal gene annotation and analysis pipelines. In summary, this study lays the foundation for developing and quality control of RNA-seq for clinical diagnostic purposes.

Discovery Phase Agrochemical Predictive Safety Assessment Using High Content In Vitro Data to Estimate an In Vivo Toxicity Point of Departure

Utilization of in vitro (cellular) techniques, like Cell Painting and transcriptomics, could provide powerful tools for agrochemical candidate sorting and selection in the discovery process. However, using these models generates challenges translating in vitro concentrations to the corresponding in vivo exposures. Physiologically based pharmacokinetic (PBPK) modeling provides a framework for quantitative in vitro to in vivo extrapolation (IVIVE). We tested whether in vivo (rat liver) transcriptomic and apical points of departure (PODs) could be accurately predicted from in vitro (rat hepatocyte or human HepaRG) transcriptomic PODs or HepaRG Cell Painting PODs using PBPK modeling. We compared two PBPK models, the ADMET predictor and the httk R package, and found httk to predict the in vivo PODs more accurately. Our findings suggest that a rat liver apical and transcriptomic POD can be estimated utilizing a combination of in vitro transcriptome-based PODs coupled with PBPK modeling for IVIVE. Thus, high content in vitro data can be translated with modest accuracy to in vivo models of ultimate regulatory importance to help select agrochemical analogs in early stage discovery program.

The E3 ubiquitin ligase Nedd4L preserves skeletal muscle stem cell quiescence by inhibiting their activation

Adult stem cells play a critical role in tissue repair and maintenance. In tissues with slow turnover, including skeletal muscle, these cells are maintained in a mitotically quiescent state yet remain poised to re-enter the cell cycle to replenish themselves and regenerate the tissue. Using a panomics approach we show that the PAX7/NEDD4L axis acts against muscle stem cell activation in homeostatic skeletal muscle. Our findings suggest that PAX7 transcriptionally activates the E3 ubiquitin ligase Nedd4L and that the conditional genetic deletion of Nedd4L impairs muscle stem cell quiescence, with an upregulation of cell cycle and myogenic differentiation genes. Loss of Nedd4L in muscle stem cells results in the expression of doublecortin (DCX), which is exclusively expressed during their in vivo activation. Together, these data establish that the ubiquitin proteasome system, mediated by Nedd4L, is a key contributor to the muscle stem cell quiescent state in adult mice.

Adaptive Regularized Tri-Factor Non-Negative Matrix Factorization for Cell Type Deconvolution

Accurate deconvolution of cell types from bulk gene expression is crucial for understanding cellular compositions and uncovering cell-type specific differential expression and physiological states of diseased tissues. Existing deconvolution methods have limitations, such as requiring complete cellular gene expression signatures or neglecting partial biological information. Moreover, these methods often overlook varying cell-type mRNA amounts, leading to biased proportion estimates. Additionally, they do not effectively utilize valuable reference information from external studies, such as means and ranges of population cell-type proportions. To address these challenges, we introduce an Adaptive Regularized Tri-factor non-negative matrix factorization approach for deconvolution (ARTdeConv). We rigorously establish the numerical convergence of our algorithm. Through benchmark simulations, we demonstrate the superior performance of ARTdeConv compared to state-of-the-art semi-reference-based and reference-free methods. In a real-world application, our method accurately estimates cell proportions, as evidenced by the nearly perfect Pearson's correlation between ARTdeConv estimates and flow cytometry measurements in a dataset from a trivalent influenza vaccine study. Moreover, our analysis of ARTdeConv estimates in COVID-19 patients reveals patterns consistent with important immunological phenomena observed in other studies. The proposed method, ARTdeConv, is implemented as an R package and can be accessed on GitHub for researchers and practitioners.

A multiomic atlas of the aging hippocampus reveals molecular changes in response to environmental enrichment

Aging involves the deterioration of organismal function, leading to the emergence of multiple pathologies. Environmental stimuli, including lifestyle, can influence the trajectory of this process and may be used as tools in the pursuit of healthy aging. To evaluate the role of epigenetic mechanisms in this context, we have generated bulk tissue and single cell multi-omic maps of the male mouse dorsal hippocampus in young and old animals exposed to environmental stimulation in the form of enriched environments. We present a molecular atlas of the aging process, highlighting two distinct axes, related to inflammation and to the dysregulation of mRNA metabolism, at the functional RNA and protein level. Additionally, we report the alteration of heterochromatin domains, including the loss of bivalent chromatin and the uncovering of a heterochromatin-switch phenomenon whereby constitutive heterochromatin loss is partially mitigated through gains in facultative heterochromatin. Notably, we observed the multi-omic reversal of a great number of aging-associated alterations in the context of environmental enrichment, which was particularly linked to glial and oligodendrocyte pathways. In conclusion, our work describes the epigenomic landscape of environmental stimulation in the context of aging and reveals how lifestyle intervention can lead to the multi-layered reversal of aging-associated decline.

An epidermal growth factor receptor-targeting immunotoxin based on IgG shows potent antitumor activity against head and neck cancer

The epidermal growth factor receptor (EGFR) is an important target for cancer therapies. Many head and neck cancer (HNC) cells have been reported to overexpress EGFR; therefore, anti-EGFR therapies have been attempted in patients with HNC. However, its clinical efficacy is limited owing to the development of drug resistance. In this study, we developed an EGFR-targeting immunotoxin consisting of a clinically proven anti-EGFR IgG (cetuximab; CTX) and a toxin fragment (LR-LO10) derived from Pseudomonas exotoxin A (PE) using a novel site-specific conjugation technology (peptide-directed photo-crosslinking reaction), as an alternative option. The immunotoxin (CTX-LR-LO10) showed specific binding to EGFR and properties of a typical IgG, such as stability, interactions with receptors of immune cells, and pharmacokinetics, and inhibited protein synthesis via modification of elongation factor-2. Treatment of EGFR-positive HNC cells with the immunotoxin resulted in apoptotic cell death and the inhibition of cell migration and invasion. The efficacy of CTX-LR-LO10 was evaluated in xenograft mouse models, and the immunotoxin exhibited much stronger tumor suppression than CTX or LR-LO10. Transcriptome analyses revealed that the immunotoxins elicited immune responses and altered the expression of genes related to its mechanisms of action. These results support the notion that CTX-LR-LO10 may serve as a new therapeutic agent targeting EGFR-positive cancers.

Fertility decline in Aedes aegypti (Diptera: Culicidae) mosquitoes is associated with reduced maternal transcript deposition and does not depend on female age

Female mosquitoes undergo multiple rounds of reproduction known as gonotrophic cycles (GC). A gonotrophic cycle spans the period from blood meal intake to egg laying. Nutrients from vertebrate host blood are necessary for completing egg development. During oogenesis, a female prepackages mRNA into her oocytes, and these maternal transcripts drive the first 2 h of embryonic development prior to zygotic genome activation. In this study, we profiled transcriptional changes in 1-2 h of Aedes aegypti (Diptera: Culicidae) embryos across 2 GC. We found that homeotic genes which are regulators of embryogenesis are downregulated in embryos from the second gonotrophic cycle. Interestingly, embryos produced by Ae. aegypti females progressively reduced their ability to hatch as the number of GC increased. We show that this fertility decline is due to increased reproductive output and not the mosquitoes' age. Moreover, we found a similar decline in fertility and fecundity across 3 GC in Aedes albopictus. Our results are useful for predicting mosquito population dynamics to inform vector control efforts.

Proteomic insights into extracellular vesicles in ALS for therapeutic potential of Ropinirole and biomarker discovery

BACKGROUND

Extracellular vesicles (EVs) hold the potential for elucidating the pathogenesis of amyotrophic lateral sclerosis (ALS) and serve as biomarkers. Notably, the comparative and longitudinal alterations in the protein profiles of EVs in serum (sEVs) and cerebrospinal fluid (CSF; cEVs) of sporadic ALS (SALS) patients remain uncharted. Ropinirole hydrochloride (ROPI; dopamine D2 receptor [D2R] agonist), a new anti-ALS drug candidate identified through induced pluripotent stem cell (iPSC)-based drug discovery, has been suggested to inhibit ALS disease progression in the Ropinirole Hydrochloride Remedy for Amyotrophic Lateral Sclerosis (ROPALS) trial, but its mechanism of action is not well understood. Therefore, we tried to reveal longitudinal changes with disease progression and the effects of ROPI on protein profiles of EVs.

METHODS

We collected serum and CSF at fixed intervals from ten controls and from 20 SALS patients participating in the ROPALS trial. Comprehensive proteomic analysis of EVs, extracted from these samples, was conducted using liquid chromatography/mass spectrometer (LC/MS). Furthermore, we generated iPSC-derived astrocytes (iPasts) and performed RNA sequencing on astrocytes with or without ROPI treatment.

RESULTS

The findings revealed notable disparities yet high congruity in sEVs and cEVs protein profiles concerning disease status, time and ROPI administration. In SALS, both sEVs and cEVs presented elevated levels of inflammation-related proteins but reduced levels associated with unfolded protein response (UPR). These results mirrored the longitudinal changes after disease onset and correlated with the revised ALS Functional Rating Scale (ALSFRS-R) at sampling time, suggesting a link to the onset and progression of SALS. ROPI appeared to counteract these changes, attenuating inflammation-related protein levels and boosting those tied to UPR in SALS, proposing an anti-ALS impact on EV protein profiles. Reverse translational research using iPasts indicated that these changes may partly reflect the DRD2-dependent neuroinflammatory inhibitory effects of ROPI. We have also identified biomarkers that predict diagnosis and disease progression by machine learning-driven biomarker search.

CONCLUSIONS

Despite the limited sample size, this study pioneers in reporting time-series proteomic alterations in serum and CSF EVs from SALS patients, offering comprehensive insights into SALS pathogenesis, ROPI-induced changes, and potential prognostic and diagnostic biomarkers.

Molecular characterization of the salivary adenoid cystic carcinoma immune landscape by anatomic subsites

Adenoid cystic carcinoma (AdCC) is a slow-growing salivary gland malignancy that relapses frequently. AdCCs of the submandibular gland exhibit unique differences in prognosis and treatment response to adjuvant radiotherapy compared to other sites, yet the role of tumor anatomic subsite on gene expression and tumor immune microenvironment (TIME) composition remains unclear. We used 87 samples, including 48 samples (27 AdCC and 21 normal salivary gland tissue samples) from 4 publicly available AdCC RNA sequencing datasets, a validation set of 33 minor gland AdCCs, and 39 samples from an in-house cohort (30 AdCC and 9 normal salivary gland samples). RNA sequencing data were used for single sample gene set enrichment analysis and TIME deconvolution. Quantitative PCR and multiplex immunofluorescence were performed on the in-house cohort. Wilcoxon rank-sum, nonparametric equality-of-medians tests and linear regression models were used to evaluate tumor subsite differences. AdCCs of different anatomic subsites including parotid, submandibular, sublingual, and minor salivary glands differed with respect to expression of several key tumorigenic pathways. Among the three major salivary glands, the reactive oxygen species (ROS)/nuclear factor erythroid 2-related factor 2 (NRF2) pathway signature was significantly underexpressed in AdCC of submandibular compared to parotid and sublingual glands while this association was not observed among normal glands. Additionally, the NRF2 pathway, whose expression was associated with favorable overall survival, was overexpressed in AdCCs of parotid gland compared to minor and submandibular glands. The TIME deconvolution identified differences in CD4+ T cell populations between AdCC of major and minor glands and natural killer (NK) cells among AdCC of minor, submandibular, and parotid glands while plasma cells were enriched in normal submandibular glands compared to other normal gland controls. Our data reveal key molecular differences in AdCC of different anatomic subsites. The ROS and NRF2 pathways are underexpressed in submandibular and minor AdCCs compared to parotid gland AdCCs, and NRF2 pathway expression is associated with favorable overall survival. The CD4+ T, NK, and plasma cell populations also vary by tumor subsites, suggesting that the observed submandibular AdCC tumor-intrinsic pathway differences may be responsible for influencing the TIME composition and survival differences.

The Pgb1 locus controls glycogen aggregation in astrocytes of the aged hippocampus without impacting cognitive function

In aged humans and mice, aggregates of hypobranched glycogen molecules called polyglucosan bodies (PGBs) accumulate in hippocampal astrocytes. PGBs are known to drive cognitive decline in neurological diseases but remain largely unstudied in the context of typical brain aging. Here, we show that PGBs arise in autophagy-dysregulated astrocytes of the aged C57BL/6J mouse hippocampus. To map the genetic cause of age-related PGB accumulation, we quantified PGB burden in 32 fully sequenced BXD-recombinant inbred mouse strains, which display a 400-fold variation in hippocampal PGB burden at 16-18 months of age. A major modifier locus was mapped to chromosome 1 at 72-75 Mb, which we defined as the Pgb1 locus. To evaluate candidate genes and downstream mechanisms by which Pgb1 controls the aggregation of glycogen, extensive hippocampal transcriptomic and proteomic datasets were produced for aged mice of the BXD family. We utilized these datasets to identify Smarcal1 and Usp37 as potential regulators of PGB accumulation. To assess the effect of PGB burden on age-related cognitive decline, we performed phenome-wide association scans, transcriptomic analyses as well as conditioned fear memory and Y-maze testing. Importantly, we did not find any evidence suggesting a negative impact of PGBs on cognition. Taken together, our study demonstrates that the Pgb1 locus controls glycogen aggregation in astrocytes of the aged hippocampus without affecting age-related cognitive decline.

RNA splicing variants of the novel long non-coding RNA, CyKILR, possess divergent biological functions in non-small cell lung cancer

The CDKN2A gene, responsible for encoding the tumor suppressors p16(INK4A) and p14(ARF), is frequently inactivated in non-small cell lung cancer (NSCLC). In this study, an uncharacterized long non-coding RNA (lncRNA) (ENSG00000267053) on chromosome 19p13.12 was found to be overexpressed in NSCLC cells with an active CDKN2A gene. This lncRNA, named Cy clin-Dependent K inase I nhibitor 2A-regulated l nc R NA (CyKILR), also correlated with the STK11 gene-coded tumor suppressor Liver kinase B1 (LKB1). CyKILR displayed two splice variants, CyKILRa (with exon 3) and CyKILRb (without exon 3), which are synergistically regulated by CDKN2A and STK11 as knockdown of both tumor suppressor genes led to a significant loss of exon 3 inclusion in mature CyKILR RNA. CyKILRa localized to the nucleus, and its downregulation using antisense RNA oligonucleotides enhanced cellular proliferation, migration, clonogenic survival, and tumor incidence. In contrast, CyKILRb localized to the cytoplasm, and downregulation of CyKILRb using siRNA reduced cell proliferation, migration, clonogenic survival, and tumor incidence. Transcriptomics analyses revealed enhancement of apoptotic pathways with concomitant suppression of key cell cycle pathways by CyKILRa demonstrating its tumor-suppressive role, while CyKILRb inhibited tumor suppressor microRNAs, indicating an oncogenic nature. These findings elucidate the intricate roles of lncRNAs in cell signaling and tumorigenesis.

Venom gland organogenesis in the common house spider

Venom is a remarkable innovation found across the animal kingdom, yet the evolutionary origins of venom systems in various groups, including spiders, remain enigmatic. Here, we investigated the organogenesis of the venom apparatus in the common house spider, Parasteatoda tepidariorum. The venom apparatus consists of a pair of secretory glands, each connected to an opening at the fang tip by a duct that runs through the chelicerae. We performed bulk RNA-seq to identify venom gland-specific markers and assayed their expression using RNA in situ hybridisation experiments on whole-mount time-series. These revealed that the gland primordium emerges during embryonic stage 13 at the chelicera tip, progresses proximally by the end of embryonic development and extends into the prosoma post-eclosion. The initiation of expression of an important toxin component in late postembryos marks the activation of venom-secreting cells. Our selected markers also exhibited distinct expression patterns in adult venom glands: sage and the toxin marker were expressed in the secretory epithelium, forkhead and sum-1 in the surrounding muscle layer, while Distal-less was predominantly expressed at the gland extremities. Our study provides the first comprehensive analysis of venom gland morphogenesis in spiders, offering key insights into their evolution and development.

IL-13 decreases susceptibility to airway epithelial SARS-CoV-2 infection but increases disease severity in vivo

Treatments available to prevent progression of virus-induced lung diseases, including coronavirus disease 2019 (COVID-19) are of limited benefit once respiratory failure occurs. The efficacy of approved and emerging cytokine signaling-modulating antibodies is variable and is affected by disease course and patient-specific inflammation patterns. Therefore, understanding the role of inflammation on the viral infectious cycle is critical for effective use of cytokine-modulating agents. We investigated the role of the type 2 cytokine IL-13 on SARS-CoV-2 binding/entry, replication, and host response in primary HAE cells in vitro and in a model of mouse-adapted SARS-CoV-2 infection in vivo. IL-13 protected airway epithelial cells from SARS-CoV-2 infection in vitro by decreasing the abundance of ACE2-expressing ciliated cells rather than by neutralization in the airway surface liquid or by interferon-mediated antiviral effects. In contrast, IL-13 worsened disease severity in mice; the effects were mediated by eicosanoid signaling and were abolished in mice deficient in the phospholipase A2 enzyme PLA2G2D. We conclude that IL-13-induced inflammation differentially affects multiple steps of COVID-19 pathogenesis. IL-13-induced inflammation may be protective against initial SARS-CoV-2 airway epithelial infection; however, it enhances disease progression in vivo. Blockade of IL-13 and/or eicosanoid signaling may be protective against progression to severe respiratory virus-induced lung disease.

Multi-omic Analyses Shed Light on The Genetic Control of High-altitude Adaptation in Sheep

Sheep were domesticated in the Fertile Crescent and then spread globally, where they have been encountering various environmental conditions. The Tibetan sheep has adapted to high altitudes on the Qinghai-Tibet Plateau over the past 3000 years. To explore genomic variants associated with high-altitude adaptation in Tibetan sheep, we analyzed Illumina short-reads of 994 whole genomes representing ∼ 60 sheep breeds/populations at varied altitudes, PacBio High fidelity (HiFi) reads of 13 breeds, and 96 transcriptomes from 12 sheep organs. Association testing between the inhabited altitudes and 34,298,967 variants was conducted to investigate the genetic architecture of altitude adaptation. Highly accurate HiFi reads were used to complement the current ovine reference assembly at the most significantly associated β-globin locus and to validate the presence of two haplotypes A and B among 13 sheep breeds. The haplotype A carried two homologous gene clusters: (1) HBE1, HBE2, HBB-like, and HBBC, and (2) HBE1-like, HBE2-like, HBB-like, and HBB; while the haplotype B lacked the first cluster. The high-altitude sheep showed highly frequent or nearly fixed haplotype A, while the low-altitude sheep dominated by haplotype B. We further demonstrated that sheep with haplotype A had an increased hemoglobin-O2 affinity compared with those carrying haplotype B. Another highly associated genomic region contained the EGLN1 gene which showed varied expression between high-altitude and low-altitude sheep. Our results provide evidence that the rapid adaptive evolution of advantageous alleles play an important role in facilitating the environmental adaptation of Tibetan sheep.

Comparative Genomics and the Salivary Transcriptome of the Redbanded Stink Bug Shed Light on Its High Damage Potential to Soybean

The redbanded stink bug, Piezodorus guildinii (Westwood) (Hemiptera: Pentatomidae), is a significant soybean pest in the Americas, which inflicts more physical damage on soybean than other native stink bugs. Studies suggest that its heightened impact is attributed to the aggressive digestive properties of its saliva. Despite its agricultural importance, the factors driving its greater ability to degrade plant tissues have remained unexplored in a genomic evolutionary context. In this study, we hypothesized that lineage-specific gene family expansions have increased the copy number of digestive genes expressed in the salivary glands. To investigate this, we annotated a previously published genome assembly of the redbanded stink bug, performed a comparative genomic analysis on 11 hemipteran species, and reconstructed patterns of gene duplication, gain, and loss in the redbanded stink bug. We also performed RNA-seq on the redbanded stink bug's salivary tissues, along with the rest of the body without salivary glands. We identified hundreds of differentially expressed salivary genes, including a subset lost in other stink bug lineages, but retained and expressed in the redbanded stink bug's salivary glands. These genes were significantly enriched with protein families involved in proteolysis, potentially explaining the redbanded stink bug's heightened damage to soybeans. Contrary to our hypothesis, we found no support for an enrichment of duplicated digestive genes that are also differentially expressed in the salivary glands of the redbanded stink bug. Nonetheless, these results provide insight into the evolution of this important crop pest, establishing a link between its genomic history and its agriculturally important physiology.

SLAM/SAP signaling regulates discrete γδ T cell developmental checkpoints and shapes the innate-like γδ TCR repertoire

During thymic development, most γδ T cells acquire innate-like characteristics that are critical for their function in tumor surveillance, infectious disease, and tissue repair. The mechanisms, however, that regulate γδ T cell developmental programming remain unclear. Recently, we demonstrated that the SLAM-SAP signaling pathway regulates the development and function of multiple innate-like γδ T cell subsets. Here, we used a single-cell proteogenomics approach to identify SAP-dependent developmental checkpoints and to define the SAP-dependent γδ TCR repertoire. SAP deficiency resulted in both a significant loss of an immature Gzma + Blk + Etv5 + Tox2 + γδT17 precursor population, and a significant increase in Cd4 + Cd8+ Rorc + Ptcra + Rag1 + thymic γδ T cells. SAP-dependent diversion of embryonic day 17 thymic γδ T cell clonotypes into the αβ T cell developmental pathway was associated with a decreased frequency of mature clonotypes in neonatal thymus, and an altered γδ TCR repertoire in the periphery. Finally, we identify TRGV4/TRAV13-4(DV7)-expressing T cells as a novel, SAP-dependent Vγ4 γδT1 subset. Together, the data suggest that SAP-dependent γδ/αβ T cell lineage commitment regulates γδ T cell developmental programming and shapes the γδ TCR repertoire.

Genetic and functional traits limit the success of colonisation by arbuscular mycorrhizal fungi in a tomato wild relative

To understand whether domestication had an impact on susceptibility and responsiveness to arbuscular mycorrhizal fungi (AMF) in tomato (Solanum lycopersicum), we investigated two tomato cultivars ("M82" and "Moneymaker") and a panel of wild relatives including S. neorickii, S. habrochaites and S. pennellii encompassing the whole Lycopersicon clade. Most genotypes revealed good AM colonisation levels when inoculated with the AMF Funneliformis mosseae. By contrast, both S. pennellii accessions analysed showed a very low colonisation, but with normal arbuscule morphology, and a negative response in terms of root and shoot biomass. This behaviour was independent of fungal identity and environmental conditions. Genomic and transcriptomic analyses revealed in S. pennellii the lack of genes identified within QTLs for AM colonisation, a limited transcriptional reprogramming upon mycorrhization and a differential regulation of strigolactones and AM-related genes compared to tomato. Donor plants experiments indicated that the AMF could represent a cost for S. pennellii: F. mosseae could extensively colonise the root only when it was part of a mycorrhizal network, but a higher mycorrhization led to a higher inhibition of plant growth. These results suggest that genetics and functional traits of S. pennellii are responsible for the limited extent of AMF colonisation.

Erg251 has complex and pleiotropic effects on sterol composition, azole susceptibility, filamentation, and stress response phenotypes

Ergosterol is essential for fungal cell membrane integrity and growth, and numerous antifungal drugs target ergosterol. Inactivation or modification of ergosterol biosynthetic genes can lead to changes in antifungal drug susceptibility, filamentation and stress response. Here, we found that the ergosterol biosynthesis gene ERG251 is a hotspot for point mutations during adaptation to antifungal drug stress within two distinct genetic backgrounds of Candida albicans. Heterozygous point mutations led to single allele dysfunction of ERG251 and resulted in azole tolerance in both genetic backgrounds. This is the first known example of point mutations causing azole tolerance in C. albicans. Importantly, single allele dysfunction of ERG251 in combination with recurrent chromosome aneuploidies resulted in bona fide azole resistance. Homozygous deletions of ERG251 caused increased fitness in low concentrations of fluconazole and decreased fitness in rich medium, especially at low initial cell density. Homozygous deletions of ERG251 resulted in accumulation of ergosterol intermediates consistent with the fitness defect in rich medium. Dysfunction of ERG251, together with FLC exposure, resulted in decreased accumulation of the toxic sterol (14-ɑ-methylergosta-8,24(28)-dien-3β,6α-diol) and increased accumulation of non-toxic alternative sterols. The altered sterol composition of the ERG251 mutants had pleiotropic effects on transcription, filamentation, and stress responses including cell membrane, osmotic and oxidative stress. Interestingly, while dysfunction of ERG251 resulted in azole tolerance, it also led to transcriptional upregulation of ZRT2, a membrane-bound Zinc transporter, in the presence of FLC, and overexpression of ZRT2 is sufficient to increase azole tolerance in wild-type C. albicans. Finally, in a murine model of systemic infection, homozygous deletion of ERG251 resulted in decreased virulence while the heterozygous deletion mutants maintain their pathogenicity. Overall, this study demonstrates that single allele dysfunction of ERG251 is a recurrent and effective mechanism of acquired azole tolerance. We propose that altered sterol composition resulting from ERG251 dysfunction mediates azole tolerance as well as pleiotropic effects on stress response, filamentation and virulence.