Tximeta: Reference sequence checksums for provenance identification in RNA-seq

Neuron-targeted caveolin-1 overexpression attenuates cognitive loss and pathological transcriptome changes in symptomatic Alzheimer's disease models

Alzheimer's disease (AD) is a devastating neurodegenerative disorder characterized by progressive synaptic loss and cognitive decline. Gene therapy that augments intrinsic neuroprotective pathways offers a promising strategy to mitigate neurodegeneration and prevent further cognitive loss. Caveolin-1 (Cav-1), a membrane lipid raft (MLR) scaffolding protein, regulates multiple pro-growth and pro-survival signaling pathways within plasmalemmal microdomains. Previously, we showed that AAV9-Synapsin-promoted Cav-1 (SynCav1) delivered to presymptomatic AD mice preserved cognitive functions and MLR-associated neurotrophic signaling. However, the therapeutic potential of SynCav1 delivered at the symptomatic stage of AD had not been tested. Therefore, the current study investigated the effect of hippocampal SynCav1 delivery at symptomatic age in two distinct preclinical AD models of amyloid pathology: PSAPP and APPKI mice. Our results demonstrated that SynCav1 delivery to PSAPP and APPKI mice at symptomatic age consistently preserved hippocampal-dependent memory. Transcriptome profiling revealed that PSAPP-SynCav1 mice exhibited a similar transcript profile to age-matched wild-type mice. Gene Ontology enrichment analysis indicated downregulation of neurodegeneration-specific pathways and upregulation of synaptic and cognitive-related pathways in PSAPP-SynCav1 mice. In vitro, SynCav1-transfected mouse primary cortical neurons exhibited increased p-CaMKII and p-CREB expression, suggesting that SynCav1 may protect the CNS by enhancing neuronal and synaptic activity. Furthermore, activity-dependent neuroprotective protein (ADNP) was identified as a potential candidate mediating SynCav1's neuroprotective effects on cognition. Subcellular membrane fractionation revealed that SynCav1 preserved MLR-localized pituitary adenylate cyclase-activating polypeptide type I receptor (PAC1R), a well-known regulator of ADNP expression. Together, these findings highlight SynCav1 as a unique and promising gene therapy candidate in the treatment of AD.

Perceptual and technical barriers in sharing and formatting metadata accompanying omics studies

Metadata, or "data about data," is essential for organizing, understanding, and managing large-scale omics datasets. It enhances data discovery, integration, and interpretation, enabling reproducibility, reusability, and secondary analysis. However, metadata sharing remains hindered by perceptual and technical barriers, including the lack of uniform standards, privacy concerns, study design limitations, insufficient incentives, inadequate infrastructure, and a shortage of trained personnel. These challenges compromise data reliability and obstruct integrative meta-analyses. Addressing these issues requires standardization, education, stronger roles for journals and funding agencies, and improved incentives and infrastructure. Looking ahead, emerging technologies such as artificial intelligence and machine learning may offer promising solutions to automate metadata processes, increasing accuracy and scalability. Fostering a collaborative culture of metadata sharing will maximize the value of omics data, accelerating innovation and scientific discovery.

Investigating endogenous immune-mediated monocyte memory in rheumatoid arthritis

OBJECTIVES

Inflammation triggered by endogenous stimuli that signal cellular stress or tissue injury must be tightly controlled to balance robust protection from intrinsic danger while avoiding catastrophic destruction of healthy tissues. Here, we assess the contribution of innate memory to this balance.

METHODS

Memory evoked by the extracellular matrix protein tenascin-C, a damage-associated, toll-like receptor 4 (TLR4) agonist, was compared to that induced by the pathogenic TLR4 agonist lipopolysaccharide (LPS) by transcriptomic and epigenetic profiling of monocytes from healthy individuals or people wirh rheumatoid arthritis (RA), and tissue macrophages from the RA synovium.

RESULTS

Tenascin-C reprograms monocyte response to subsequent threats, inducing concomitantly suppressed and enhanced responses to rechallenge. Comparative analysis of tenascin-C and LPS revealed common and distinct gene expression signatures, effects controlled transcriptionally and associated with stimulus-specific epigenetic mediators. Altered responses following rechallenge after priming with tenascin-C were not limited to subsequent TLR4 activation but were evident in response to various pathogenic and endogenous stimuli detected by different receptors. In healthy monocytes primed with tenascin-C, rechallenge with stimuli found at high levels in the joints of people with RA resulted in trained responses that were not induced by LPS, including genes associated with chronic inflammation, tissue destruction, altered metabolism, and poor treatment response in RA. The expression of a large subset of these genes was elevated in monocytes from people with RA in the absence of any stimulation and in RA synovial macrophage populations associated with disease flare. Moreover, higher levels of permissive complexes within key epigenetic nodes and increased bivalent modification creating poised loci within endogenously trained genes were observed in RA cells.

CONCLUSIONS

These data highlight how innate reprogramming during 'sterile' inflammatory diseases contributes to chronicity, uncovering pathways unique to endogenous immune triggers that could provide disease-specific points of intervention without engendering global immune suppression.

Cervicovaginal microbiome alters transcriptomic and epigenomic signatures across cervicovaginal epithelial barriers

Background

The cervicovaginal microbiome plays a critical role in women's health, with microbial communities dominated by Lactobacillus species considered optimal. In contrast, the depletion of lactobacilli and the presence of a diverse array of strict and facultative anaerobes, such as Gardnerella vaginalis, have been linked with adverse reproductive outcomes. Despite these associations, the molecular mechanisms by which host-microbial interactions modulate cervical and vaginal epithelial function remains poorly understood.

Results

In this study, we used RNA sequencing to characterize the transcriptional response of cervicovaginal epithelial cells exposed to the culture supernatants of common vaginal bacteria. Our findings revealed that G. vaginalis culture supernatants upregulate genes associated with an activated innate immune response and increased cell death. Conversely, Lactobacillus crispatus culture supernatants induced transcriptional changes indicative of epigenomic modeling in ectocervical epithelial cells. Epigenomic modification by L. crispatus, was confirmed by ATAC-sequencing, which demonstrated reduced chromatin accessibility.

Conclusions

These results provide new insights into host-microbe interactions within the lower reproductive tract and suggests that modulating the vaginal microbiome could offer innovative therapeutic strategies to improve reproductive health.

Dynamic changes in gene expression through aging in Drosophila melanogaster heads

Work in many systems has shown large-scale changes in gene expression during aging. However, many studies employ just 2 arbitrarily chosen timepoints to measure expression and can only observe an increase or a decrease in expression between "young" and "old" animals, failing to capture any dynamic, nonlinear changes that occur throughout the aging process. We used RNA sequencing to measure expression in male head tissue at 15 timepoints through the lifespan of an inbred Drosophila melanogaster strain. We detected >6,000 significant, age-related genes, nearly all of which have been seen in previous Drosophila aging expression studies and that include several known to harbor lifespan-altering mutations. We grouped our gene set into 28 clusters via their temporal expression change, observing a diversity of trajectories; some clusters show a linear change over time, while others show more complex, nonlinear patterns. Notably, reanalysis of our dataset comparing the earliest and latest timepoints-mimicking a 2-timepoint design-revealed fewer differentially expressed genes (around 4,500). Additionally, those genes exhibiting complex expression trajectories in our multitimepoint analysis were most impacted in this reanalysis; their identification, and the inferred change in gene expression with age, was often dependent on the timepoints chosen. Informed by our trajectory-based clusters, we executed a series of gene enrichment analyses, identifying enriched functions/pathways in all clusters, including the commonly seen increase in stress- and immune-related gene expression with age. Finally, we developed a pair of accessible Shiny apps to enable exploration of our differential expression and gene enrichment results.

A comprehensive Schizosaccharomyces pombe atlas of physical transcription factor interactions with proteins and chromatin

Transcription factors (TFs) are key regulators of gene expression, yet many of their targets and modes of action remain unknown. In Schizosaccharomyces pombe, one-third of TFs are solely homology predicted, with few experimentally validated. We created a comprehensive library of 89 endogenously tagged S. pombe TFs, mapping their protein and chromatin interactions using immunoprecipitation-mass spectrometry and chromatin immunoprecipitation sequencing. Our study identified protein interactors for half the TFs, with over a quarter potentially forming stable complexes. We discovered DNA-binding sites for most TFs across 2,027 unique genomic regions, revealing motifs for 38 TFs and uncovering a complex network of extensive TF cross- and autoregulation. Characterization of the largest TF family revealed conserved DNA sequence preferences but diverse binding patterns and identified a repressive heterodimer, Ntu1/Ntu2, linked to perinuclear gene localization. Our TFexplorer webtool makes all data interactively accessible, offering insights into TF interactions and regulatory mechanisms with broad biological relevance.

Extrachromosomal circular DNA as a novel biomarker for the progression of colorectal cancer

BACKGROUND

Extrachromosomal circular DNA (eccDNA) has potential in tumor diagnosis, particularly for improving diagnostic accuracy and early cancer detection; however, many challenges remain in its application to clinical practice.

METHODS

We conducted a Circle-Seq analysis on clinical samples at different stages of colorectal cancer progression to examine the dynamic changes of eccDNA during the progression of colorectal cancer. We used breakpoint-specific PCR to verify candidate eccDNAs identified by Circle-Seq. The results were further validated using the AOM/DSS-induced colorectal cancer model.

RESULTS

There was an increase in the abundance of eccDNA with the progression of colorectal cancer. The genes associated with these eccDNA molecules were primarily related to signaling pathways involved in tumor development and metastasis. Our analysis also revealed that eccDNA abundance positively correlates with gene expression, and eccDNA derived from specific genes has potential value for the early diagnosis of tumors.

CONCLUSIONS

This study revealed a connection between eccDNA and colorectal cancer progression and highlights the clinical potential of eccDNA for the early diagnosis of colorectal cancer.

Integrated promoter-capture Hi-C and Hi-C analysis reveals fine-tuned regulation of the 3D chromatin architecture in colorectal cancer

Introduction

Hi-C is a widely used technique for mapping chromosomal interactions within a 3D genomic framework, however, its resolution is often constrained by sequencing depth, making it challenging to detect fine-scale interactions. To overcome this limitation, Promoter-Capture Hi-C (PCHi-C), as it selectively enriches for promoter-associated interactions, was employed. This study integrates PCHi-C and Hi-C datasets from colorectal cancer (CRC) models investigate chromosomal interaction dynamics across various regulatory levels, from cis-regulatory elements to topologically associated domains (TADs). The primary goal is to examine how genomic structural alterations shape the epigenomic landscape in CRC and to assess their potential role in colorectal cancer susceptibility.

Methods

PCHi-C and Hi-C datasets from multiple colorectal cancer (CRC) studies were integrated to enhance the resolution of chromatin interaction mapping. The analysis focused on identifying fine-scale interactions within topologically associated domains (TADs) while incorporating histone modification landscapes (H3K27ac, H3K4me3) and transcriptomic signatures from CRC cell lines and the TCGA database. For experimental validation, ChIP-quantitative PCR was performed at the promoters of target genes using the highly malignant colorectal cell line HT29 and compared it to an embryonic cell line NT2D1.

Results

Our integrated analysis revealed significant genomic structural instability in CRC cells, closely associated with tumor-suppressive transcriptional programs. We identified nine dysregulated genes, including long non-coding RNAs (MALAT1, NEAT1, FTX, and PVT1), small nucleolar RNAs (SNORA26 and SNORA71A), and protein-coding genes (TMPRSS11D, TSPEAR, and DSG4), all of which exhibited a substantial increase in expression in CRC cell lines compared to human embryonic stem cells (hESCs). Additionally, we observed enriched activation-associated histone modifications (H3K27ac and H3K4me3) at the potential enhancer regions of these genes, indicating possible transcriptional activation. ChIP-quantitative PCRs conducted using in the highly malignant CRC cell line HT29, compared to the embryonic cell line NT2D1, further validated these findings, reinforcing the link between altered chromosomal interactions and gene dysregulation in CRC.

Discussion

This study sheds light on the dynamic 3D genome organization in CRC, highlighting critical structural changes associated with disease-associated loci. The identification of nine dysregulated genes points to potential biomarkers for colorectal cancer, with implications for diagnostic and therapeutic strategies. The combination of Hi-C and PCHi-C offers a refined approach for detecting chromosomal interactions at a higher resolution, laying the foundation for future studies on cancer-associated chromatin architecture.

Monocyte-lineage tumor infiltration predicts immunoradiotherapy response in advanced pretreated soft-tissue sarcoma: phase 2 trial results

Immunoradiotherapy holds promise for improving outcomes in patients with advanced solid tumors, including in soft-tissue sarcoma (STS). However, the ideal combination of treatment modalities remains to be determined, and reliable biomarkers to predict which patients will benefit are lacking. Here, we report the results of the STS cohort of the SABR-PDL1 phase II trial that evaluated the anti-PDL1 atezolizumab combined with stereotactic body radiation therapy (SBRT) delivered concurrently with the 2nd cycle to at least one tumor site. Eligible patients received atezolizumab until progression or unmanageable toxicity, with SBRT at 45 Gy in 3 fractions). The primary endpoint was one-year progression-free survival (PFS) rate with success defined as 13 patients achieving 1-year PFS. Sixty-one heavily pretreated patients with STS (median 5 prior lines; 52% men; median age 54 years; 28% leiomyosarcoma) were enrolled across two centers (France, Spain). SBRT was delivered to 55 patients (90%), with the lung being the most commonly irradiated site (50%). After a median follow-up of 45 months, the one-year PFS rate was 8.3% [95% CI: 3.6-18.1]. Median PFS and overall survival were 2.5 and 8.6 months, respectively. Best responses included partial responses (5%) and stable disease (60%). Immune profiling revealed increased immunosuppressive tumor-associated macrophages (e.g., IL4I1, HES1) and monocyte-recruiting chemokines in non-responders. Higher monocyte/lymphocyte ratios (MonoLR) in tumor and blood correlated with progression. PD-L1 status, lymphoid infiltration, and tertiary-lymphoid structures were not predictive. Although the primary endpoint was not met, this study highlights MonoLR imbalance as a potential biomarker to identify STS patients likely to benefit from immunoradiotherapy. EudraCT No. 2015-005464-42; Clinicaltrial.gov number: NCT02992912.

Image-based identification and isolation of micronucleated cells to dissect cellular consequences

Recent advances in isolating cells based on visual phenotypes have transformed our ability to identify the mechanisms and consequences of complex traits. Micronucleus (MN) formation is a frequent outcome of genome instability, triggers extensive changes in genome structure and signaling coincident with MN rupture, and is almost exclusively defined by visual analysis. Automated MN detection in microscopy images has proved challenging, limiting discovery of the mechanisms and consequences of MN. In this study we describe two new MN segmentation modules: a rapid model for classifying micronucleated cells and their rupture status (VCS MN), and a robust model for accurate MN segmentation (MNFinder) from a broad range of cell lines. As proof-of-concept, we define the transcriptome of non-transformed human cells with intact or ruptured MN after chromosome missegregation by combining VCS MN with photoactivation-based cell isolation and RNASeq. Surprisingly, we find that neither MN formation nor rupture triggers a strong unique transcriptional response. Instead, transcriptional changes appear correlated with small increases in aneuploidy in these cell classes. Our MN segmentation modules overcome a significant challenge with reproducible MN quantification, and, joined with visual cell sorting, enable the application of powerful functional genomics assays to a wide-range of questions in MN biology.

Diet-driven transcriptional changes in weaning red abalone (Haliotis rufescens) and its hybrid (H. rufescens [♀] x H. fulgens [♂])

The abalone (Haliotis sp) aquaculture industry is on the rise fueled by its high commercial demand. The use of abalone interspecific hybrids is gaining attention due to their improved growth and tolerance to environmental challenges. However, hybrids may respond differently to dietary inputs compared to their parental species, which could be optimized to achieve maximum growth potential. Accordingly, we previously reported that the postlarvae of a hybrid cross obtained by mating red and green abalone (H. rufescens [♀] x H. fulgens [♂]) presented higher growth and survival rates than red abalone postlarvae when fed with Macrocystis pyrifera (M) and a mixture of Ulva onhoi and Navicula incerta (UN). Here, to gain a deeper insight into the physiological and metabolic responses of both crosses, we report the results of a transcriptome sequencing analysis to unveil the diet-driven physiological changes within each cross. Our findings indicate that the diet induces modulation in the transcriptional machinery related to macromolecules breakdown, especially, the mixed diet UN induced higher expression of lysosomal enzymes, especially in enzymes involved in glycosaminoglycan catabolism, such as arylsulfatases, glucosidases, and alpha-fucosidases, which might indicate a higher ability to break down complex carbohydrates from algae. These transcriptional responses might explain the growth difference between pure and hybrid as well as for abalone fed with M and UN.

Neonatal fungi promote lifelong metabolic health through macrophage-dependent β cell development

Loss of early-life microbial diversity is correlated with diabetes, yet mechanisms by which microbes influence disease remain elusive. We report a critical neonatal window in mice when microbiota disruption results in lifelong metabolic consequences stemming from reduced β cell development. We show evidence for the existence of a similar program in humans and identify specific fungi and bacteria that are sufficient for β cell growth. The microbiota also plays an important role in seeding islet-resident macrophages, and macrophage depletion during development reduces β cells. Candida dubliniensis increases β cells in a macrophage-dependent manner through distinctive cell wall composition and reduces murine diabetes incidence. Provision of C. dubliniensis after β cell ablation or antibiotic treatment improves β cell function. These data identify fungi as critical early-life commensals that promote long-term metabolic health.

Inter-chromosomal insertions at Xq27.1 associated with retinal dystrophy induce dysregulation of LINC00632 and CDR1as/ciRS-7

In two unrelated families with X-linked inherited retinal dystrophy, identification of the causative variants was elusive. Interrogation of the next-generation sequencing (NGS) data revealed a "dark" intergenic region on Xq27.1 with poor coverage. Long-range PCR and DNA walking across this region revealed different inter-chromosomal insertions into the human-specific palindrome on Xq27.1: a 58 kb insertion of 9p24.3 [der(X)dir ins(X;9)(q27.1;p24.3)] in family 1 and a 169 kb insertion of 3p14.2 [der(X)inv ins(X;3)(q27.1;p14.2)] in family 2. To explore the functional consequence of these structural variants in genomic and cellular contexts, induced pluripotent stem cells were derived from affected and control fibroblasts and differentiated to retinal organoids (ROs) and retinal pigment epithelium. Transcriptional dysregulation was evaluated using RNA sequencing (RNA-seq) and RT-qPCR. A downstream long non-coding RNA, LINC00632 (Xq27.1), was upregulated in ROs from both families compared to control samples. In contrast, the circular RNA CDR1as/ciRS-7 (circular RNA sponge for miR-7), spliced from linear LINC00632, was downregulated. To investigate this tissue-specific dysregulation, we interrogated the landscape of the locus using Hi-C and cleavage under targets and tagmentation sequencing (CUT&Tag). This revealed active retinal enhancers within the insertions within a topologically associated domain that also contained the upstream promoter of LINC00632, permitting ectopic contact. Furthermore, CDR1as/ciRS-7 acts as a "sponge" for miR-7, and target genes of miR-7 were also dysregulated in ROs derived from both families. We describe a new genomic mechanism for retinal dystrophy, and our data support a convergent tissue-specific mechanism of altered regulation of LINC00632 and CDR1as/ciRS-7 as a consequence of the insertions within the palindrome on Xq27.1.

A novel role for Friend of GATA1 (FOG-1) in regulating cholesterol transport in murine erythropoiesis

Friend of GATA1 (FOG-1) is an essential transcriptional co-factor of the master erythroid transcription factor GATA1. The knockout of the Zfpm1 gene, coding for FOG-1, results in early embryonic lethality due to anemia in mice, similar to the embryonic lethal phenotype of the Gata1 gene knockout. However, a detailed molecular analysis of the Zfpm1 knockout phenotype in erythropoiesis is presently incomplete. To this end, we used CRISPR/Cas9 to knockout Zfpm1 in mouse erythroleukemic (MEL) cells. Phenotypic characterization of DMSO-induced terminal erythroid differentiation showed that the Zfpm1 knockout MEL cells did not progress past the proerythroblast stage of differentiation. Expression profiling of the Zfpm1 knockout MEL cells by RNAseq showed a lack of up-regulation of erythroid-related gene expression profiles. Bioinformatic analysis highlighted cholesterol transport as a pathway affected in the Zfpm1 knockout cells. Moreover, we show that the cholesterol transporters Abca1 and Ldlr fail to be repressed during erythroid differentiation in Zfpm1 knockout cells, resulting in higher intracellular lipid levels and higher membrane fluidity. We also show that in FOG-1 knockout cells, the nuclear levels of SREBP2, a key transcriptional regulator of cholesterol biosynthesis and transport, are markedly increased. On the basis of these findings we propose that FOG-1 (and, potentially, GATA1) regulate cholesterol homeostasis during erythroid differentiation directly through the down regulation of cholesterol transport genes and indirectly, through the repression of the SREBP2 transcriptional activator of cholesterol homeostasis. Taken together, our work provides a molecular basis for understanding FOG-1 functions in erythropoiesis and reveals a novel role for FOG-1 in cholesterol transport.

Identification of Meibomian gland stem cell populations and mechanisms of aging

Meibomian glands secrete lipid-rich meibum, which prevents tear evaporation. Aging-related Meibomian gland shrinkage may result in part from stem cell exhaustion and is associated with evaporative dry eye disease, a common condition lacking effective treatment. The identities and niche of Meibomian gland stem cells and the signals controlling their activity are poorly defined. Using snRNA-seq, in vivo lineage tracing, ex vivo live imaging, and genetic studies in mice, we identify markers for stem cell populations that maintain distinct regions of the gland and uncover Hedgehog (Hh) signaling as a key regulator of stem cell proliferation. Consistent with this, we show that human Meibomian gland carcinoma exhibits increased Hh signaling. Aged glands display decreased Hh and EGF signaling, deficient innervation, and loss of collagen I in niche fibroblasts, indicating that alterations in both glandular epithelial cells and their surrounding microenvironment contribute to age-related degeneration. These findings suggest new approaches to treat aging-associated Meibomian gland loss.

SPT5 regulates RNA polymerase II stability via Cullin 3-ARMC5 recognition

The stability of RNA polymerase II (Pol II) is tightly regulated during transcriptional elongation for proper control of gene expression. Our recent studies revealed that promoter-proximal Pol II is destabilized via the ubiquitin E3 ligase cullin 3 (CUL3) upon loss of transcription elongation factor SPT5. Here, we investigate how CUL3 recognizes chromatin-bound Pol II as a substrate. Using an unbiased proteomic screening approach, we identify armadillo repeat-containing 5 (ARMC5) as a CUL3 adaptor required for VCP/p97-dependent degradation of SPT5-depleted, chromatin-bound Pol II. Genome-wide analyses indicate that ARMC5 targets promoter-proximal Pol II in a BTB domain-dependent manner. Further biochemical analysis demonstrates that interaction between ARMC5 and Pol II requires the transcriptional cyclin-dependent kinase 9 (CDK9), supporting a phospho-dependent degradation model. We propose that defective, promoter-proximal Pol II that lacks SPT5 is rapidly eliminated from chromatin in a noncanonical early termination pathway that requires CDK9-dependent interaction with the CUL3-ARMC5 ubiquitin ligase complex.

Single-nucleus CUT&RUN elucidates the function of intrinsic and genomics-driven epigenetic heterogeneity in head and neck cancer progression

Interrogating regulatory epigenetic alterations during tumor progression at the resolution of single cells has remained an understudied area of research. Here we developed a highly sensitive single-nucleus CUT&RUN (snCUT&RUN) assay to profile histone modifications in isogenic primary, metastatic, and cisplatin-resistant head and neck squamous cell carcinoma (HNSCC) patient-derived tumor cell lines. We find that the epigenome can be involved in diverse modes to contribute toward HNSCC progression. First, we demonstrate that gene expression changes during HNSCC progression can be comodulated by alterations in both copy number and chromatin activity, driving epigenetic rewiring of cell states. Furthermore, intratumor epigenetic heterogeneity (ITeH) may predispose subclonal populations within the primary tumor to adapt to selective pressures and foster the acquisition of malignant characteristics. In conclusion, snCUT&RUN serves as a valuable addition to the existing toolkit of single-cell epigenomic assays and can be used to dissect the functionality of the epigenome during cancer progression.

Utilization of primary tumor samples for cancer neoantigen discovery

BACKGROUND

The use of tumor-infiltrating T lymphocytes (TIL) that recognize cancer neoantigens has led to lasting remissions in metastatic melanoma and certain cases of metastatic epithelial cancer. For the treatment of the latter, selecting cells for therapy typically involves laborious screening of TIL for recognition of autologous tumor-specific mutations, detected through next-generation sequencing of freshly resected metastatic tumors. Our study explored the feasibility of using archived formalin-fixed, paraffin-embedded (FFPE) primary tumor samples for cancer neoantigen discovery, to potentially expedite this process and reduce the need for resections normally required for tumor sequencing.

METHOD

Whole-exome sequencing was conducted on matched primary and metastatic colorectal cancer samples from 22 patients. The distribution of metastatic tumor mutations that were confirmed as neoantigens through cognate TIL screening was evaluated in the corresponding primary tumors. Mutations unique to primary tumors were screened for recognition by metastasis-derived TIL and circulating T lymphocytes.

RESULTS

We found that 25 (65.8%) of the 38 validated neoantigens identified in metastatic tumors from 18 patients with colorectal cancer were also present in matched primary tumor samples. This included all 12 neoantigens encoded by putative cancer driver genes, which are generally regarded as superior targets for adoptive cell therapy. The detection rate for other neoantigens, representing mutations without an established role in cancer biology, was 50% (13/26). Gene products encoding neoantigens detected in the primary tumors were not more likely to be clonal or broadly distributed among the analyzed metastatic lesions compared with those undetected in the primary tumors. Additionally, we found that mutations detected only in primary tumor samples did not elicit recognition by metastatic tumor-derived TIL but could elicit specific recognition by the autologous circulating memory T cells.

CONCLUSIONS

Our findings indicate that primary FFPE tumor-derived screening libraries could be used to discover most neoantigens present in metastatic tumors requiring treatment. Furthermore, this approach can reveal additional neoantigens not present in resected metastatic tumors, prompting further research to understand their clinical relevance as potential therapeutic targets.

Whole blood transcriptome profile identifies motor neurone disease RNA biomarker signatures

Blood-based biomarkers for motor neuron disease are needed for better diagnosis, progression prediction, and clinical trial monitoring. We used whole blood-derived total RNA and performed whole transcriptome analysis to compare the gene expression profiles in (motor neurone disease) MND patients to the control subjects. We compared 42 MND patients to 42 aged and sex-matched healthy controls and described the whole transcriptome profile characteristic for MND. In addition to the formal differential analysis, we performed functional annotation of the genomics data and identified the molecular pathways that are differentially regulated in MND patients. We identified 12,972 genes differentially expressed in the blood of MND patients compared to age and sex-matched controls. Functional genomic annotation identified activation of the pathways related to neurodegeneration, RNA transcription, RNA splicing and extracellular matrix reorganisation. Blood-based whole transcriptomic analysis can reliably differentiate MND patients from controls and can provide useful information for the clinical management of the disease and clinical trials.

Response eQTLs, chromatin accessibility, and 3D chromatin structure in chondrocytes provide mechanistic insight into osteoarthritis risk

Osteoarthritis (OA) poses a significant healthcare burden with limited treatment options. While genome-wide association studies (GWASs) have identified over 100 OA-associated loci, translating these findings into therapeutic targets remains challenging. To address this gap, we mapped gene expression, chromatin accessibility, and 3D chromatin structure in primary human articular chondrocytes in both resting and OA-mimicking conditions. We identified thousands of differentially expressed genes, including those associated with differences in sex and age. RNA sequencing in chondrocytes from 101 donors across two conditions uncovered 3,782 unique eGenes, including 420 that exhibited strong and significant condition-specific effects. Colocalization with OA GWAS signals revealed 13 putative OA risk genes, 6 of which have not been previously identified. Chromatin accessibility and 3D chromatin structure provided insights into the mechanisms and conditional specificity of these variants. Our findings shed light on OA pathogenesis and highlight potential targets for therapeutic development.

Clec12a controls colitis by tempering inflammation and restricting expansion of specific commensals

Microbiota composition regulates colitis severity, yet the innate immune mechanisms that control commensal communities and prevent disease remain unclear. We show that the innate immune receptor, Clec12a, impacts colitis severity by regulating microbiota composition. Transplantation of microbiota from a Clec12a-/- animal is sufficient to worsen colitis in wild-type mice. Clec12a-/- mice have expanded Faecalibaculum rodentium, and treatment with F. rodentium similarly exacerbates disease. However, Clec12a-/- animals are resistant to colitis development when rederived into an 11-member community, underscoring the role of specific species. Colitis in Clec12a-/- mice is dependent on monocytes, and cytokine and sequencing analysis in Clec12a-/- macrophages and serum shows enhanced inflammation with a reduction in phagocytic genes. F. rodentium specifically binds to Clec12a, and Clec12a-/--deficient macrophages are impaired in their ability to phagocytose F. rodentium. Thus, Clec12a contributes to an innate-immune-surveillance mechanism that controls the expansion of potentially harmful commensals while limiting inflammation.

Galaxy as a gateway to bioinformatics: Multi-Interface Galaxy Hands-on Training Suite (MIGHTS) for scRNA-seq

BACKGROUND

Bioinformatics is fundamental to biomedical sciences, but its mastery presents a steep learning curve for bench biologists and clinicians. Learning to code while analyzing data is difficult. The curve may be flattened by separating these two aspects and providing intermediate steps for budding bioinformaticians. Single-cell analysis is in great demand from biologists and biomedical scientists, as evidenced by the proliferation of training events, materials, and collaborative global efforts like the Human Cell Atlas. However, iterative analyses lacking reinstantiation, coupled with unstandardized pipelines, have made effective single-cell training a moving target.

FINDINGS

To address these challenges, we present a Multi-Interface Galaxy Hands-on Training Suite (MIGHTS) for single-cell RNA sequencing (scRNA-seq) analysis, which offers parallel analytical methods using a graphical interface (buttons) or code. With clear, interoperable materials, MIGHTS facilitates smooth transitions between environments. Bridging the biologist-programmer gap, MIGHTS emphasizes interdisciplinary communication for effective learning at all levels. Real-world data analysis in MIGHTS promotes critical thinking and best practices, while FAIR data principles ensure validation of results. MIGHTS is freely available, hosted on the Galaxy Training Network, and leverages Galaxy interfaces for analyses in both settings. Given the ongoing popularity of Python-based (Scanpy) and R-based (Seurat & Monocle) scRNA-seq analyses, MIGHTS enables analyses using both.

CONCLUSIONS

MIGHTS consists of 11 tutorials, including recordings, slide decks, and interactive visualizations, and a demonstrated track record of sustainability via regular updates and community collaborations. Parallel pathways in MIGHTS enable concurrent training of scientists at any programming level, addressing the heterogeneous needs of novice bioinformaticians.

Safety of baricitinib in vaccinated patients with severe and critical COVID-19 sub study of the randomised Bari-SolidAct trial

BACKGROUND

The Bari-SolidAct randomized controlled trial compared baricitinib with placebo in patients with severe COVID-19. A post hoc analysis revealed a higher incidence of serious adverse events (SAEs) among SARS-CoV-2-vaccinated participants who had received baricitinib. This sub-study aimed to investigate whether vaccination influences the safety profile of baricitinib in patients with severe COVID-19.

METHODS

Biobanked samples from 146 participants (55 vaccinated vs. 91 unvaccinated) were analysed longitudinally for inflammation markers, humoral responses, tissue viral loads, and plasma viral antigens on days 1, 3, and 8. High-dimensional analyses, including RNA sequencing and flow cytometry, were performed on available samples. Mediation analyses were used to assess relationships between SAEs, baseline-adjusted biomarkers, and treatment-vaccination status.

FINDINGS

Vaccinated participants were older, more frequently hospitalized, had more comorbidities, and exhibited higher nasopharyngeal viral loads. Baricitinib treatment did not affect antibody responses or viral clearance, but reduced markers of T-cell and monocyte activation compared to placebo (sCD25, sCD14, sCD163, sTIM-3). Age, baseline levels of plasma viral antigen, and several inflammatory markers, as well as IL-2, IL-6, Neopterin, CXCL16, sCD14, and suPAR on day 8 were associated with the occurrence of SAEs. However, mediation analyses of markers linked to SAEs, baricitinib treatment, or vaccination status did not reveal statistically significant interactions between vaccination status and SAEs.

INTERPRETATION

This sub-study did not identify any virus- or host-related biomarkers significantly associated with the interaction between SARS-CoV-2 vaccination status and the safety of baricitinib. However, caution should be exercised due to the moderate sample size.

FUNDING

EU Horizon 2020 (grant number 101015736).

Metabolomic Response to Non-Steroidal Anti-Inflammatory Drugs

Non-steroidal anti-inflammatory drugs (NSAIDs) are popular choices for the mitigation of pain and inflammation; however, they are accompanied by side effects in the gastrointestinal and cardiovascular systems. We compared the effects of naproxen, a traditional NSAID, and celecoxib, a cyclooxygenase - 2 (Cox-2) inhibitor, in humans. Our findings showed a decrease in tryptophan and kynurenine levels in plasma of volunteers treated with naproxen. We further validated this result in mice. Additionally, we find that the depression of tryptophan was independent of both Cox-1 and Cox-2 inhibition, but rather was due to the displacement of bound tryptophan by naproxen. Supplementation of tryptophan in naproxen-treated mice rescued fecal blood loss and inflammatory gene expression driven by IL-1β in the heart.

WNT signaling contributes to the extrahepatic bile duct proliferative response to obstruction in mice

Biliary obstruction and cholangiocyte hyperproliferation are important features of cholangiopathies affecting the large extrahepatic bile duct (EHBD). The mechanisms underlying obstruction-induced cholangiocyte proliferation in the EHBD remain poorly understood. Developmental pathways, including WNT signaling, are implicated in regulating injury responses in many tissues, including the liver. To investigate the contribution of WNT signaling to obstruction-induced cholangiocyte proliferation in the EHBD, we used complementary in vivo and in vitro models with pharmacologic interventions and transcriptomic analyses. To model obstruction, we used bile duct ligation (BDL) in mice. Human and mouse biliary organoids and mouse biliary explants were used to investigate the effects of WNT activation and inhibition in vitro. We observed an upregulation of WNT ligand expression associated with increased biliary proliferation following obstruction. Cholangiocytes were identified as both WNT ligand-expressing and WNT-responsive cells. Inhibition of WNT signaling decreased cholangiocyte proliferation in vivo and in vitro, while activation increased proliferation. WNT effects on cholangiocyte proliferation were β-catenin dependent, and we showed a direct effect of WNT7B on cholangiocyte growth. Our studies suggested that cholangiocyte-derived WNT ligands can activate WNT signaling to induce proliferation after obstructive injury. These findings implicate the WNT pathway in injury-induced cholangiocyte proliferation within the EHBD.

Metabolomic Response to Non-Steroidal Anti-Inflammatory Drugs

Non-steroidal anti-inflammatory drugs (NSAIDs) are popular choices for the mitigation of pain and inflammation; however, they are accompanied by side effects in the gastrointestinal and cardiovascular systems. We compared the effects of naproxen, a traditional NSAID, and celecoxib, a cyclooxygenase -2 (Cox-2) inhibitor, in humans. Our findings showed a decrease in tryptophan and kynurenine levels in plasma of volunteers treated with naproxen. We further validated this result in mice. Additionally, we find that the depression of tryptophan was independent of both Cox-1 and Cox-2 inhibition, but rather was due to the displacement of bound tryptophan by naproxen. Supplementation of tryptophan in naproxen-treated mice rescued fecal blood loss and inflammatory gene expression driven by IL-1β in the heart.

Connecting verified databases with clinical practice and the patient' s experience through omnichannel communication

INTRODUCTION

Patient-reported outcomes (PRO) collect data directly from patients. These data are utilized in clinical practice, helping decision-making. Studies emphasize the importance of omnichannel communication (WhatsApp, e-mail, SMS) with healthcare professionals and patients. Omnichannel communication enables the integration of different communication channels to improve the end-client experience. In addition to the means of communication, the daily practice of professionals requires different activities that can be performed in distinct systems. The existence of various separate systems for other activities in medical practice may result in complexities and bottlenecks in their use by healthcare professionals and patients.

OBJECTIVE

To present the Digital Health Ecosystem (DHE) that unifies scientific research with medical practice in omnichannel communication and mechanisms to verify the authenticity and integrity of the data collected and stored.

METHODOLOGY

The system requirements and needs were met utilizing the Iconix development methodology. Microsoft Dot Net was used to develop software. Usability, usefulness and user satisfaction with the system were measured using the Post-Study System Usability Questionnaire (PSSUQ).

RESULTS

Omnichannel communication was utilized to contact patients and healthcare professionals autonomously. A single system enabled the carrying out of patientreported outcome data collection, telemedicine, image storage, and notes from patient consultations. The data was collected through structured questionnaires via link and chatbot. The functionalities created in the HDE allowed the integrity and authenticity verification of the data collected and stored.

CONCLUSION

Personalized omnichannel communication via links and chatbots using WhatsApp, E-mail, and SMS accelerates autonomous interaction with patients and healthcare professionals. In addition, the structured and non-structured data were stored in the EHD and able to be verified for integrity and authenticity.

A SIRT7-dependent acetylation switch regulates early B cell differentiation and lineage commitment through Pax5

B lymphopoiesis is orchestrated by lineage-specific transcription factors. In B cell progenitors, lineage commitment is mediated by Pax5, which is commonly mutated in B cell acute lymphoblastic leukemia. Despite its essential role in immunity, the mechanisms regulating Pax5 function remain largely unknown. Here, we found that the NAD+-dependent enzyme SIRT7 coordinates B cell development through deacetylation of Pax5 at K198, which promotes Pax5 protein stability and transcriptional activity. Neither Pax5K198 deacetylated nor acetylated mimics rescued B cell differentiation in Pax5-/- pro-B cells, suggesting that B cell development requires Pax5 dynamic deacetylation. The Pax5K198 deacetylation mimic restored lineage commitment in Pax5-/- pro-B cells and B cell differentiation in Sirt7-/- pro-B cells, suggesting the uncoupling of differentiation from lineage commitment. The SIRT7-Pax5 interplay was conserved in B cell acute lymphoblastic leukemia, where SIRT7 expression correlated with good prognosis. Our findings reveal a crucial mechanism for B lymphopoiesis and highlight the relevance of sirtuins in immune function.

Modulation of UPF1 catalytic activity upon interaction of SARS-CoV-2 Nucleocapsid protein with factors involved in nonsense mediated-mRNA decay

The RNA genome of the SARS-CoV-2 virus encodes for four structural proteins, 16 non-structural proteins and nine putative accessory factors. A high throughput analysis of interactions between human and SARS-CoV-2 proteins identified multiple interactions of the structural Nucleocapsid (N) protein with RNA processing factors. The N-protein, which is responsible for packaging of the viral genomic RNA was found to interact with two RNA helicases, UPF1 and MOV10 that are involved in nonsense-mediated mRNA decay (NMD). Using a combination of biochemical and biophysical methods, we investigated the interaction of the SARS-CoV-2 N-protein with NMD factors at a molecular level. Our studies led us to identify the core NMD factor, UPF2, as an interactor of N. The viral N-protein engages UPF2 in multipartite interactions and can negate the stimulatory effect of UPF2 on UPF1 catalytic activity. N also inhibits UPF1 ATPase and unwinding activities by competing in binding to the RNA substrate. We further investigate the functional implications of inhibition of UPF1 catalytic activity by N in mammalian cells. The interplay of SARS-CoV-2 N with human UPF1 and UPF2 does not affect decay of host cell NMD targets but might play a role in stabilizing the viral RNA genome.

GGNBP2 regulates MDA5 sensing triggered by self double-stranded RNA following loss of ADAR1 editing

Adenosine-to-inosine (A-to-I) editing of double-stranded RNA (dsRNA) by ADAR1 is an essential modifier of the immunogenicity of cellular dsRNA. The role of MDA5 in sensing unedited cellular dsRNA and the downstream activation of type I interferon (IFN) signaling are well established. However, we have an incomplete understanding of pathways that modify the response to unedited dsRNA. We performed a genome-wide CRISPR screen and showed that GGNBP2, CNOT10, and CNOT11 interact and regulate sensing of unedited cellular dsRNA. We found that GGNBP2 acts between dsRNA transcription and its cytoplasmic sensing by MDA5. GGNBP2 loss prevented induction of type I IFN and autoinflammation after the loss of ADAR1 editing activity by modifying the subcellular distribution of endogenous A-to-I editing substrates and reducing cytoplasmic dsRNA load. These findings reveal previously undescribed pathways to modify diseases associated with ADAR mutations and may be determinants of response or resistance to small-molecule ADAR1 inhibitors.

Solid tumor immunotherapy using NKG2D-based adaptor CAR T cells

NKG2D ligands (NKG2DLs) are broadly expressed in cancer. To target these, we describe an adaptor chimeric antigen receptor (CAR) termed NKG2D/Dap10-12. Herein, T cells are engineered to co-express NKG2D with a fusion protein that comprises Dap10 joined to a Dap12 endodomain. NKG2D/Dap10-12 T cells elicit compelling efficacy, eradicating or controlling NKG2DL-expressing tumors in several established xenograft models. Importantly, durable responses, long-term survival, and rejection of tumor re-challenge are reproducibly achieved. Efficacy is markedly superior to a clinical stage CAR analog, comprising an NKG2D-CD3ζ fusion. Structure-function analysis using an extended CAR panel demonstrates that potency is dependent on membrane proximity of signaling units, high NKG2D cell surface expression, adaptor structure, provision of exogenous Dap10, and inclusion of one rather than three immune tyrosine activation motifs per signaling unit. Potent therapeutic impact of NKG2D/Dap10-12 T cells is also underpinned by enhanced oxidative phosphorylation, reduced senescence, and transcriptomic re-programming for increased ribosomal biogenesis.

Epigenome Reprogramming Through H3K27 and H3K4 Trimethylation as a Resistance Mechanism to DNA Methylation Inhibition in BRAFV600E-Mutated Colorectal Cancer

PURPOSE

BRAFV600E-mutated colorectal cancer exhibits a strong correlation with DNA hypermethylation, suggesting that this subgroup of tumors presents unique epigenomic phenotypes. Nonetheless, 5-azacitidine, which inhibits DNA methyltransferase activity, is not efficacious in BRAFV600E colorectal cancer in vivo.

EXPERIMENTAL DESIGN

We randomized and treated mice implanted with patient-derived tumor xenografts harboring BRAFV600E mutation with control, 5-azacitidine, vemurafenib (BRAF inhibitor), or the combination. Comprehensive epigenomic profiling was conducted on control and 5-azacitidine-treated tumor samples, including DNA methylation, histone modifications, chromatin accessibility, and gene expression. Combinations of epigenetic agents were explored in preclinical BRAFV600E colorectal cancer models.

RESULTS

A profound reduction of DNA methylation levels upon 5-azacitidine treatment was confirmed, however, transcriptional repression was not relieved. This study unbiasedly explored the adaptive engagement of other epigenomic modifications upon 5-azacitidine treatment. A loss of histone acetylation and a gain of histone methylations, including H3K27 and H3K4 trimethylation, were observed around these hypomethylated regions, suggesting the involvement of polycomb repressive complex (PRC) activity around the genome with loss of DNA methylation, therefore maintaining the repression of key tumor-suppressor genes. Combined inhibition of PRC activity through EZH2 inhibition with 5-azacitidine treatment additively improved efficacies in BRAFV600E colorectal cancer cells.

CONCLUSIONS

In conclusion, DNA hypomethylation by 5-azacitidine exhibits a close association with H3K27me3 and PRC activity in BRAFV600E colorectal cancer, and simultaneous blockade of DNA methyltransferase and EZH2 holds promise as a potential therapeutic strategy for patients with BRAFV600E-mutated colorectal cancer.

Response splicing QTLs in primary human chondrocytes identifies putative osteoarthritis risk genes

Osteoarthritis affects millions worldwide, yet effective treatments remain elusive due to poorly understood molecular mechanisms. While genome-wide association studies (GWAS) have identified over 100 OA-associated loci, identifying the genes impacted at each locus remains challenging. Several studies have mapped expression quantitative trait loci (eQTL) in chondrocytes and colocalized them with OA GWAS variants to identify putative OA risk genes; however, the degree to which genetic variants influence OA risk via alternative splicing has not been explored. We investigated the role of alternative splicing in OA pathogenesis using RNA-seq data from 101 human chondrocyte samples treated with PBS (control) or fibronectin fragment (FN-f), an OA trigger. We identified 590 differentially spliced genes between conditions, with FN-f inducing splicing events similar to those in primary OA tissue. We used CRISPR/Cas9 to mimic an SNRNP70 splicing event observed in OA and FN-f-treated chondrocytes and found that it induced an OA-like expression pattern. Integration with genotyping data revealed 7,188 splicing quantitative trait loci (sQTL) affecting 3,056 genes. While many sQTLs were shared, we identified 738 and 343 condition-specific sQTLs for control and FN-f, respectively. We identified 15 RNA binding proteins whose binding sites were enriched at sQTL splice junctions and found that expression of those RNA binding proteins correlated with exon inclusion. Colocalization with OA GWAS identified 6 putative risk genes, including a novel candidate, PBRM1. Our study highlights the significant impact of alternative splicing in OA and provides potential therapeutic targets for future research.

Continuous infusion of resolvin D2 in combination with Angiotensin-II show contrary effects on blood pressure and intracardiac artery remodeling

Specialized pro-resolving mediators (SPMs) are key effectors of resolution of inflammation. This is highly relevant for cardiac and vessel remodeling, where the net inflammatory response contributes to determine disease outcome. Herein, we used a mice model of angiotensin (Ang)-II-induced hypertension to study the effect of the SPM Resolvin D2 (RvD2), on hypertension and cardiac remodeling. By using subcutaneous osmotic minipumps, mice were treated with PBS or Ang-II in combination with or without RvD2 for two weeks. Mice receiving RvD2 gained less blood pressure increase compared to Ang-II alone. Surprisingly, however, examination of intracardiac arteries revealed that RvD2 treatment in combination with Ang-II exacerbated Ang-II-induced fibrosis. Measures of vascular smooth muscle cell dedifferentiation correlated with the level of vascular remodeling, indicating that this dedifferentiation, including increased proliferation and migration, is a contributing factor. RNA sequencing of left ventricle cardiac tissue supported these findings as pathways related to cell proliferation and cell differentiation were upregulated in mice treated with Ang-II in combination with RvD2. Additionally, the RNA sequencing also showed upregulation of pathways related to SPM metabolism. In line with this, Mass spectrometry analysis of lipid mediators showed reduced cardiac levels of the arachidonic acid derived metabolite leukotriene E4 in RvD2 treated mice. Our study suggests that continuous infusion through osmotic minipumps should not be the recommended route of RvD2 administration in future studies.

Response eQTLs, chromatin accessibility, and 3D chromatin structure in chondrocytes provide mechanistic insight into osteoarthritis risk

Osteoarthritis (OA) poses a significant healthcare burden with limited treatment options. While genome-wide association studies (GWAS) have identified over 100 OA-associated loci, translating these findings into therapeutic targets remains challenging. Integrating expression quantitative trait loci (eQTL), 3D chromatin structure, and other genomic approaches with OA GWAS data offers a promising approach to elucidate disease mechanisms; however, comprehensive eQTL maps in OA-relevant tissues and conditions remain scarce. We mapped gene expression, chromatin accessibility, and 3D chromatin structure in primary human articular chondrocytes in both resting and OA-mimicking conditions. We identified thousands of differentially expressed genes, including those associated with differences in sex and age. RNA-seq in chondrocytes from 101 donors across two conditions uncovered 3782 unique eGenes, including 420 that exhibited strong and significant condition-specific effects. Colocalization with OA GWAS signals revealed 13 putative OA risk genes, 10 of which have not been previously identified. Chromatin accessibility and 3D chromatin structure provided insights into the mechanisms and conditional specificity of these variants. Our findings shed light on OA pathogenesis and highlight potential targets for therapeutic development.

Isoform-level analyses of 6 cancers uncover extensive genetic risk mechanisms undetected at the gene-level

Integrating genome-wide association study (GWAS) and transcriptomic datasets can help identify potential mediators for germline genetic risk of cancer. However, traditional methods have been largely unsuccessful because of an overreliance on total gene expression. These approaches overlook alternative splicing, which can produce multiple isoforms from the same gene, each with potentially different effects on cancer risk. Here, we integrate genetic and multi-tissue isoform-level gene expression data from the Genotype Tissue-Expression Project (GTEx, N = 108-574) with publicly available European-ancestry GWAS summary statistics (all N > 20,000 cases) to identify both isoform- and gene-level risk associations with six cancers (breast, endometrial, colorectal, lung, ovarian, prostate) and six related cancer subtype classifications (N = 12 total). Compared to traditional methods leveraging total gene expression, directly modeling isoform expression through transcriptome-wide association studies (isoTWAS) substantially increases discovery of transcriptomic mechanisms underlying genetic associations. Using the same RNA-seq datasets, isoTWAS identified 164% more significant unique gene associations compared to TWAS (6,163 and 2,336, respectively), with isoTWAS-prioritized genes enriched 4-fold for evolutionarily-constrained genes (P = 6.1 × 10-13). isoTWAS tags transcriptomic associations at 52% more independent GWAS loci compared to TWAS across the six cancers. Additionally, isoform expression mediates an estimated 63% greater proportion of cancer risk SNP heritability compared to gene expression when evaluating cis-genetic influence on isoform expression. We highlight several notable isoTWAS associations that demonstrate GWAS colocalization at the isoform level but not at the gene level, including, CLPTM1L (lung cancer), LAMC1 (colorectal), and BABAM1 (breast). These results underscore the critical importance of modeling isoform-level expression to maximize discovery of genetic risk mechanisms for cancers.

Integration of long-read sequencing, DNA methylation and gene expression reveals heterogeneity in Y chromosome segment lengths in phenotypic males with 46,XX testicular disorder/difference of sex development

BACKGROUND

46,XX testicular disorder/difference of sex development (46,XX DSD) is a rare congenital condition, characterized by a combination of the typical female sex chromosome constitution, 46,XX, and a variable male phenotype. In the majority of individuals with 46,XX DSD, a Y chromosome segment containing the sex-determining region gene (SRY) has been translocated to the paternal X chromosome. However, the precise genomic content of the translocated segment and the genome-wide effects remain elusive.

METHODS

We performed long-read DNA sequencing, RNA sequencing and DNA methylation analyses on blood samples from 46,XX DSD (n = 11), male controls (46,XY; variable cohort sizes) and female controls (46,XX; variable cohort sizes), in addition to RNA sequencing and DNA methylation analysis on blood samples from males with Klinefelter syndrome (47,XXY, n = 22). We also performed clinical measurements on all 46,XX DSD and a subset of 46,XY (n = 10).

RESULTS

We identified variation in the translocated Y chromosome segments, enabling subcategorization into 46,XX DSD (1) lacking Y chromosome material (n = 1), (2) with short Yp arms (breakpoint at 2.7-2.8 Mb, n = 2), (3) with medium Yp arms (breakpoint at 7.3 Mb, n = 1), and (4) with long Yp arms (n = 7), including deletions of AMELY, TBLY1 and in some cases PRKY. We also identified variable expression of the X-Y homologues PRKY and PRKX. The Y-chromosomal transcriptome and methylome reflected the Y chromosome segment lengths, while changes to autosomal and X-chromosomal regions indicated global effects. Furthermore, transcriptional changes tentatively correlated with phenotypic traits of 46,XX DSD, including reduced height, lean mass and testicular size.

CONCLUSION

This study refines our understanding of the genetic composition in 46,XX DSD, describing the translocated Y chromosome segment in more detail than previously and linking variability herein to genome-wide changes in the transcriptome and methylome.

Butterfly eggs prime anti-herbivore defense in an annual but not perennial Arabidopsis species

MAIN CONCLUSION

Unlike Arabidopsis thaliana, defenses of Arabidopsis lyrata against Pieris brassicae larval feeding are not primable by P. brassicae eggs. Thus, egg primability of plant anti-herbivore defenses is not phylogenetically conserved in the genus Arabidopsis. While plant anti-herbivore defenses of the annual species Arabidopsis thaliana were shown to be primable by Pieris brassicae eggs, the primability of the phylogenetically closely related perennial Arabidopsis lyrata has not yet been investigated. Previous studies revealed that closely related wild Brassicaceae plant species, the annual Brassica nigra and the perennial B. oleracea, exhibit an egg-primable defense trait, even though they have different life spans. Here, we tested whether P. brassicae eggs prime anti-herbivore defenses of the perennial A. lyrata. We exposed A. lyrata to P. brassicae eggs and larval feeding and assessed their primability by (i) determining the biomass of P. brassicae larvae after feeding on plants with and without prior P. brassicae egg deposition and (ii) investigating the plant transcriptomic response after egg deposition and/or larval feeding. For comparison, these studies were also conducted with A. thaliana. Consistent with previous findings, A. thaliana's response to prior P. brassicae egg deposition negatively affected conspecific larvae feeding upon A. thaliana. However, this was not observed in A. lyrata. Arabidopsis thaliana responded to P. brassicae eggs with strong transcriptional reprogramming, whereas A. lyrata responses to eggs were negligible. In response to larval feeding, A. lyrata exhibited a greater transcriptome change compared to A. thaliana. Among the strongly feeding-induced A. lyrata genes were those that are egg-primed in feeding-induced A. thaliana, i.e., CAX3, PR1, PR5, and PDF1.4. These results suggest that A. lyrata has evolved a robust feeding response that is independent from prior egg exposure.

HSP60 chaperone deficiency disrupts the mitochondrial matrix proteome and dysregulates cholesterol synthesis

OBJECTIVE

Mitochondrial proteostasis is critical for cellular function. The molecular chaperone HSP60 is essential for cell function and dysregulation of HSP60 expression has been implicated in cancer and diabetes. The few reported patients carrying HSP60 gene variants show neurodevelopmental delay and brain hypomyelination. Hsp60 interacts with more than 260 mitochondrial proteins but the mitochondrial proteins and functions affected by HSP60 deficiency are poorly characterized.

METHODS

We studied two model systems for HSP60 deficiency: (1) engineered HEK cells carrying an inducible dominant negative HSP60 mutant protein, (2) zebrafish HSP60 knockout larvae. Both systems were analyzed by RNASeq, proteomics, and targeted metabolomics, and several functional assays relevant for the respective model. In addition, skin fibroblasts from patients with disease-associated HSP60 variants were analyzed by proteomics.

RESULTS

We show that HSP60 deficiency leads to a differentially downregulated mitochondrial matrix proteome, transcriptional activation of stress responses, and dysregulated cholesterol biosynthesis. This leads to lipid accumulation in zebrafish knockout larvae.

CONCLUSIONS

Our data provide a compendium of the effects of HSP60 deficiency on the mitochondrial matrix proteome. We show that HSP60 is a master regulator and modulator of mitochondrial functions and metabolic pathways. HSP60 dysfunction also affects cellular metabolism and disrupts the integrated stress response. The effect on cholesterol synthesis explains the effect of HSP60 dysfunction on myelination observed in patients carrying genetic variants of HSP60.

A global transcriptional atlas of the effect of acute sleep deprivation in the mouse frontal cortex

Sleep deprivation (SD) has negative effects on brain and body function. Sleep problems are prevalent in a variety of disorders, including neurodevelopmental and psychiatric conditions. Thus, understanding the molecular consequences of SD is of fundamental importance in biology. In this study, we present the first simultaneous bulk and single-nuclear RNA sequencing characterization of the effects of SD in the male mouse frontal cortex. We show that SD predominantly affects glutamatergic neurons, specifically in layers 4 and 5, and produces isoform switching of over 1500 genes, particularly those involved in splicing and RNA binding. At both the global and cell-type specific level, SD has a large repressive effect on transcription, downregulating thousands of genes and transcripts. As a resource we provide extensive characterizations of cell-types, genes, transcripts, and pathways affected by SD. We also provide publicly available tutorials aimed at allowing readers adapt analyses performed in this study to their own datasets.

ICOS limits memory-like properties and function of exhausted PD-1 + CD8 T cells

During persistent antigen stimulation, PD-1 + CD8 T cells are maintained by progenitor exhausted PD-1 + TCF-1 + CD8 T cells (Tpex). Tpex respond to PD-1 blockade, and regulation of Tpex differentiation into more functional Tex is of major interest for cancer immunotherapies. Tpex express high levels of Inducible Costimulator (ICOS), but the role of ICOS for PD-1 + CD8 T cell responses has not been addressed. In chronic infection, ICOS-deficiency increased both number and quality of virus-specific CD8 T cells, with accumulation of effector-like Tex due to enhanced survival. Mechanistically, loss of ICOS signaling potentiated FoxO1 activity and memory-like features of Tpex. In mice with established chronic infection, ICOS-Ligand blockade resulted in expansion of effector-like Tex and reduction in viral load. In a mouse model of hepatocellular carcinoma, ICOS inhibition improved cytokine production by tumor-specific PD-1 + CD8 T cells and delayed tumor growth. Overall, we show that ICOS limits CD8 T cell responses during chronic antigen exposure.

Differential biological effect of low doses of ionizing radiation depending on the radiosensitivity in a cell line model

PURPOSE

Exposure to low doses (LD) of ionizing radiation (IR), such as the ones employed in computed tomography (CT) examination, can be associated with cancer risk. However, cancer development could depend on individual radiosensitivity. In the present study, we evaluated the differences in the response to a CT-scan radiation dose of 20 mGy in two lymphoblastoid cell lines with different radiosensitivity.

MATERIALS AND METHODS

Several parameters were studied: gene expression, DNA damage, and its repair, as well as cell viability, proliferation, and death. Results were compared with those after a medium dose of 500 mGy.

RESULTS

After 20 mGy of IR, the radiosensitive (RS) cell line showed an increase in DNA damage, and higher cell proliferation and apoptosis, whereas the radioresistant (RR) cell line was insensitive to this LD. Interestingly, the RR cell line showed a higher expression of an antioxidant gene, which could be used by the cells as a protective mechanism. After a dose of 500 mGy, both cell lines were affected by IR but with significant differences. The RS cells presented an increase in DNA damage and apoptosis, but a decrease in cell proliferation and cell viability, as well as less antioxidant response.

CONCLUSIONS

A differential biological effect was observed between two cell lines with different radiosensitivity, and these differences are especially interesting after a CT scan dose. If this is confirmed by further studies, one could think that individuals with radiosensitivity-related genetic variants may be more vulnerable to long-term effects of IR, potentially increasing cancer risk after LD exposure.

MTGR1 is required to maintain small intestinal stem cell populations

Undifferentiated intestinal stem cells (ISCs) are crucial for maintaining homeostasis and resolving injury. Lgr5+ cells in the crypt base constantly divide, pushing daughter cells upward along the crypt axis where they differentiate into specialized cell types. Coordinated execution of complex transcriptional programs is necessary to allow for the maintenance of undifferentiated stem cells while permitting differentiation of the wide array of intestinal cells necessary for homeostasis. Previously, members of the myeloid translocation gene (MTG) family have been identified as transcriptional co-repressors that regulate stem cell maintenance and differentiation programs in multiple organ systems, including the intestine. One MTG family member, myeloid translocation gene related 1 (MTGR1), has been recognized as a crucial regulator of secretory cell differentiation and response to injury. However, whether MTGR1 contributes to the function of ISCs has not yet been examined. Here, using Mtgr1-/- mice, we have assessed the effects of MTGR1 loss specifically in ISC biology. Interestingly, loss of MTGR1 increased the total number of cells expressing Lgr5, the canonical marker of cycling ISCs, suggesting higher overall stem cell numbers. However, expanded transcriptomic and functional analyses revealed deficiencies in Mtgr1-null ISCs, including deregulated ISC-associated transcriptional programs. Ex vivo, intestinal organoids established from Mtgr1-null mice were unable to survive and expand due to aberrant differentiation and loss of stem and proliferative cells. Together, these results indicate that the role of MTGR1 in intestinal differentiation is likely stem cell intrinsic and identify a novel role for MTGR1 in maintaining ISC function.

Benchmarking transcriptome deconvolution methods for estimating tissue- and cell-type-specific extracellular vesicle abundances

Extracellular vesicles (EVs) contain cell-derived lipids, proteins and RNAs; however, determining the tissue- and cell-type-specific EV abundances in body fluids remains a significant hurdle for our understanding of EV biology. While tissue- and cell-type-specific EV abundances can be estimated by matching the EV's transcriptome to a tissue's/cell type's expression signature using deconvolutional methods, a comparative assessment of deconvolution methods' performance on EV transcriptome data is currently lacking. We benchmarked 11 deconvolution methods using data from four cell lines and their EVs, in silico mixtures, 118 human plasma and 88 urine EVs. We identified deconvolution methods that estimated cell type-specific abundances of pure and in silico mixed cell line-derived EV samples with high accuracy. Using data from two urine EV cohorts with different EV isolation procedures, four deconvolution methods produced highly similar results. The three methods were also concordant in their tissue- and cell-type-specific plasma EV abundance estimates. We identified driving factors for deconvolution accuracy and highlighted the importance of implementing biological knowledge in creating the tissue/cell type signature. Overall, our analyses demonstrate that the deconvolution algorithms DWLS and CIBERSORTx produce highly similar and accurate estimates of tissue- and cell-type-specific EV abundances in biological fluids.

TnP and AHR-CYP1A1 Signaling Crosstalk in an Injury-Induced Zebrafish Inflammation Model

Aryl Hydrocarbon Receptor (AHR) signaling is crucial for regulating the biotransformation of xenobiotics and physiological processes like inflammation and immunity. Meanwhile, Thalassophryne nattereri Peptide (TnP), a promising anti-inflammatory candidate from toadfish venom, demonstrates therapeutic effects through immunomodulation. However, its influence on AHR signaling remains unexplored. This study aimed to elucidate TnP's molecular mechanisms on the AHR-cytochrome P450, family 1 (CYP1) pathway upon injury-induced inflammation in wild-type (WT) and Ahr2-knockdown (KD) zebrafish larvae through transcriptomic analysis and Cyp1a reporters. TnP, while unable to directly activate AHR, potentiated AHR activation by the high-affinity ligand 6-Formylindolo [3,2-b]carbazole (FICZ), implying a role as a CYP1A inhibitor, confirmed by in vitro studies. This interplay suggests TnP's ability to modulate the AHR-CYP1 complex, prompting investigations into its influence on biotransformation pathways and injury-induced inflammation. Here, the inflammation model alone resulted in a significant response on the transcriptome, with most differentially expressed genes (DEGs) being upregulated across the groups. Ahr2-KD resulted in an overall greater number of DEGs, as did treatment with the higher dose of TnP in both WT and KD embryos. Genes related to oxidative stress and inflammatory response were the most apparent under inflamed conditions for both WT and KD groups, e.g., Tnfrsf1a, Irf1b, and Mmp9. TnP, specifically, induces the expression of Hspa5, Hsp90aa1.2, Cxcr3.3, and Mpeg1.2. Overall, this study suggests an interplay between TnP and the AHR-CYP1 pathway, stressing the inflammatory modulation through AHR-dependent mechanisms. Altogether, these results may offer new avenues in novel therapeutic strategies, such as based on natural bioactive molecules, harnessing AHR modulation for targeted and sustained drug effects in inflammatory conditions.

Transcriptional dynamics of sleep deprivation and subsequent recovery sleep in the male mouse cortex

Sleep is an essential, tightly regulated biological function. Sleep is also a homeostatic process, with the need to sleep increasing as a function of being awake. Acute sleep deprivation (SD) increases sleep need, and subsequent recovery sleep (RS) discharges it. SD is known to alter brain gene expression in rodents, but it remains unclear which changes are linked to sleep homeostasis, SD-related impairments, or non-sleep-specific effects. To investigate this question, we analyzed RNA-seq data from adult wild-type male mice subjected to 3 and 5-6 hours of SD and 2 and 6 hours of RS after SD. We hypothesized molecular changes associated with sleep homeostasis mirror sleep pressure dynamics as defined by brain electrical activity, peaking at 5-6 hours of SD, and are no longer differentially expressed after 2 hours of RS. We report 5-6 hours of SD produces the largest effect on gene expression, affecting approximately half of the cortical transcriptome, with most differentially expressed genes (DEGs) downregulated. The majority of DEGs normalize after 2 hours of RS and are involved in redox metabolism, chromatin regulation, and DNA damage/repair. Additionally, RS affects gene expression related to mitochondrial metabolism and Wnt-signaling, potentially contributing to its restorative effects. DEGs associated with cholesterol metabolism and stress response do not normalize within 6 hours and may be non-sleep-specific. Finally, DEGs involved in insulin signaling, MAPK signaling, and RNA-binding may mediate the impairing effects of SD. Overall, our results offer insight into the molecular mechanisms underlying sleep homeostasis and the broader effects of SD.

Circular RNA Profile in Atherosclerotic Disease: Regulation during ST-Elevated Myocardial Infarction

Circular (circ) RNAs are non-coding RNAs with important functions in the nervous system, cardiovascular system, and cancer. Their role in atherosclerosis and myocardial infarction (MI) remains poorly described. We aim to investigate the potential circRNAs in immune cells during atherogenesis and examine the most regulated during MI and the modulation by interleukin (IL)-6 receptor inhibition by tocilizumab. Wild-type (WT) and ApoE-/- mice were fed an atherogenic diet for 10 weeks, and the circRNA profile was analyzed by circRNA microarray. Whole blood from patients with ST-elevated MI (STEMI) and randomized to tocilizumab (n = 21) or placebo (n = 19) was collected at admission, 3-7 days, and at 6 months, in addition to samples from healthy controls (n = 13). Primers for human circRNA were designed, and circRNA levels were measured using RT-qPCR. mRNA regulation of predicted circRNA targets was investigated by RNA sequencing. The expression of 867 circRNAs differed between atherogenic and WT mice. In STEMI patients, circUBAC2 was significantly lower than in healthy controls. CircANKRD42 and circUBAC2 levels were inversely correlated with troponin T, and for circUBAC2, an inverse correlation was also seen with final infarct size at 6 months. The predicted mRNA targets for circUBAC2 and circANKRD42 were investigated and altered levels of transcripts involved in the regulation of inflammatory/immune cells, apoptosis, and mitochondrial function were found. Finally, tocilizumab induced an up-regulation of circANKRD42 and circUBAC2 3-7 days after percutaneous coronary intervention. CircRNA levels were dysregulated in STEMI, potentially influencing the immune system, apoptosis, and mitochondrial function.

Identifying a gene signature of metastatic potential by linking pre-metastatic state to ultimate metastatic fate

Identifying the key molecular pathways that enable metastasis by analyzing the eventual metastatic tumor is challenging because the state of the founder subclone likely changes following metastatic colonization. To address this challenge, we labeled primary mouse pancreatic ductal adenocarcinoma (PDAC) subclones with DNA barcodes to characterize their pre-metastatic state using ATAC-seq and RNA-seq and determine their relative in vivo metastatic potential prospectively. We identified a gene signature separating metastasis-high and metastasis-low subclones orthogonal to the normal-to-PDAC and classical-to-basal axes. The metastasis-high subclones feature activation of IL-1 pathway genes and high NF-κB and Zeb/Snail family activity and the metastasis-low subclones feature activation of neuroendocrine, motility, and Wnt pathway genes and high CDX2 and HOXA13 activity. In a functional screen, we validated novel mediators of PDAC metastasis in the IL-1 pathway, including the NF-κB targets Fos and Il23a, and beyond the IL-1 pathway including Myo1b and Tmem40. We scored human PDAC tumors for our signature of metastatic potential from mouse and found that metastases have higher scores than primary tumors. Moreover, primary tumors with higher scores are associated with worse prognosis. We also found that our metastatic potential signature is enriched in other human carcinomas, suggesting that it is conserved across epithelial malignancies. This work establishes a strategy for linking cancer cell state to future behavior, reveals novel functional regulators of PDAC metastasis, and establishes a method for scoring human carcinomas based on metastatic potential.

Identification of Meibomian gland stem cell populations and mechanisms of aging

Meibomian glands secrete lipid-rich meibum, which prevents tear evaporation. Aging-related Meibomian gland shrinkage may result in part from stem cell exhaustion and is associated with evaporative dry eye disease, a common condition lacking effective treatment. The identities and niche of Meibomian gland stem cells and the signals controlling their activity are poorly defined. Using snRNA-seq, in vivo lineage tracing, ex vivo live imaging, and genetic studies in mice, we identified markers for stem cell populations that maintain distinct regions of the gland and uncovered Hh signaling as a key regulator of stem cell proliferation. Consistent with this, human Meibomian gland carcinoma exhibited increased Hh signaling. Aged glands displayed decreased Hh and EGF signaling, deficient innervation, and loss of collagen I in niche fibroblasts, indicating that alterations in both glandular epithelial cells and their surrounding microenvironment contribute to age-related degeneration. These findings suggest new approaches to treat aging-associated Meibomian gland loss.

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.