WebGestalt: an integrated system for exploring gene sets in various biological contexts

Genome data based deep learning identified new genes predicting pharmacological treatment response of attention deficit hyperactivity disorder

Although the efficacy of pharmacy in the treatment of attention deficit/hyperactivity disorder (ADHD) has been well established, the lack of predictors of treatment response poses great challenges for personalized treatment. The current study employed a comprehensive approach, combining genome-wide association analyses (GWAS) and deep learning (DL) methods, to elucidate the genetic underpinnings of pharmacological treatment response in ADHD. Based on genotype data of medication-naïve patients with ADHD who received pharmacological treatments for 12 weeks, the current study performed GWAS using the percentage changes in ADHD-RS score as phenotype. Then, DL models were constructed to predict percentage changes in symptom scores using genetic variants selected based on four different genome-wide P thresholds (E-02, E-03, E-04, E-05) as inputs. The current GWAS results identified two significant loci (rs10880574, P = 2.39E-09; rs2000900, P = 3.31E-09) which implicated two genes, TMEM117 and MYO5B, that were primarily associated with both brain- and gut-related disorders. The convolutional neural network (CNN) model, using variants with genome-wide P values less than E-02 (5516 SNPs), demonstrated the best performance with mean squared error (MSE) equals 0.012 (Accuracy = 0.83; Sensitivity = 0.90; Specificity = 0.75) in the validation dataset, 0.081 in an independent test dataset (Acc = 0.61, Sensitivity = 0.81; Specificity = 0.26). Notably, the variant that contributed most to the CNN model was NKAIN2, an ADHD-related gene, which is also associated with metabolic processes. To conclude, the integration of GWAS and DL methods revealed new genes contribute to ADHD pharmacological treatment responses, and underscored the interplay between neural systems and metabolic processes, potentially providing critical insights into precision treatment. Furthermore, our CNN model exhibited good performance in an independent dataset, encouraged future studies and implied potential clinical applications.

Glymphatic system clearance and Alzheimer's disease risk: a CSF proteome-wide study

BACKGROUND

The emerging evidence of the role of the glymphatic system (GS) in Alzheimer's disease (AD) provides new opportunities for intervention from the earliest stages of the disease. The aim of the study is to evaluate the efficacy of GS in AD to identify new disease biomarkers.

METHODS

We performed a two-stage proteomic study to evaluate the GS health using intravenous gadolinium-based contrast agent (GBCA) with serial T1 3T magnetic resonance imaging (MRI) in individuals with amnestic mild cognitive impairment (aMCI). In Stage 1 (evaluated in the Cohort 1 of aMCI participants (n = 11)), we correlated the levels of 7K cerebrospinal fluid (CSF) proteins (estimated by SOMAscan) with GS health in 78 Freesurfer-segmented brain regions of interest (ROIs).

RESULTS

A total of seven different proteins were significantly associated with GS health (p-value < 6.4 × 10-4). The stronger correlations were identified for NSUN6, GRAAK, OLFML3, ACTN2, RUXF, SHPS1 and TIM-4. A pathway enrichment analysis revealed that the proteins associated with GS health were mainly implicated in neurodegenerative processes, immunity and inflammation. In Stage 2, we validated these proteomic results in a new cohort of aMCI participants (with and without evidence of AD pathology in CSF (aMCI(-) and aMCI/AD( +); n = 22 and 7, respectively) and healthy controls (n = 10). Proteomic prediction models were generated in each ROI. These were compared with demographic-only models for identifying participants with aMCI(-) and aMCI/AD( +) vs controls. This analysis was repeated to determine if the models could identify those with aMCI/AD( +) from both aMCI(-) and controls. The proteomic models were found to outperform the demographic-only models.

CONCLUSIONS

Our study identifies proteins linked with GS health and involved the immune system in aMCI participants.

The immunomodulatory role of the MAFB gene in hepatocellular carcinoma and its impact on biological activities

OBJECTIVE

The transcription factor MAFB is part of the MAF family and is known to promote hepatocellular carcinoma (HCC) by upregulating cyclin D1. However, its role in HCC immunity and prognosis remains unclear. This study explores the biological function, prognostic significance, and immune impact of MAFB in HCC.

METHODS

Immunohistochemistry was used to analyze MAFB expression in HCC and adjacent non-tumor tissues. RT-qPCR and Western blotting measured MAFB levels in HCC cell lines. Specific siRNA was used to knockdown MAFB in HCCLM3 and MHCC97H cells, followed by assays to evaluate cell proliferation, migration, and colony formation. Data from the TCGA database and online tools TIMER and TISDB were used to assess the relationship between MAFB and immune responses. A prognostic model based on MAFB-related immune genes was established, and drug sensitivity analysis was performed.

RESULTS

MAFB was significantly overexpressed in HCC tissues. Knockdown of MAFB in HCC cell lines reduced their proliferation and migration abilities. The risk model based on MAFB-related immune genes effectively predicted patient prognosis, supported by ROC curves. Gene set enrichment analysis indicated that MAFB is involved in immune-related pathways. Several drugs were identified as potentially sensitive to MAFB expression levels.

CONCLUSION

MAFB plays a significant role in HCC development and immune regulation. The prognostic model combining MAFB-related immune genes provides valuable insights for predicting patient outcomes and identifying potential therapeutic targets.

Animal and vegetal materials of mouse oocytes segregate at first zygotic cleavage: a simple mechanism that makes the two-cell blastomeres differ reciprocally from the start

Recent advances in embryology have shown that the sister blastomeres of two-cell mouse and human embryos differ reciprocally in potency. An open question is whether the blastomeres became different as opposed to originating as different. Here we wanted to test two relevant but conflicting models: one proposing that each blastomere contains both animal and vegetal materials in balanced proportions because the plane of first cleavage runs close to the animal-vegetal axis of the fertilized oocyte (meridional cleavage); and the other model proposing that each blastomere contains variable proportions of animal and vegetal materials because the plane of the first cleavage can vary - up to an equatorial orientation - depending on the topology of fertilization. Therefore, we imposed the fertilization site in three distinct regions of mouse oocytes (animal pole, vegetal pole, equator) via ICSI. After the first zygotic cleavage, the sister blastomeres were dissociated and subjected to single-cell transcriptome analysis, keeping track of the original pair associations. Non-supervised hierarchical clustering revealed that the frequency of correct pair matches varied with the fertilization site (vegetal pole > animal pole > equator), thereby, challenging the first model of balanced partitioning. However, the inter-blastomere differences had similar signatures of gene ontology across the three groups, thereby, also challenging the competing model of variable partitioning. These conflicting observations could be reconciled if animal and vegetal materials were partitioned at the first cleavage: an event considered improbable and possibly deleterious in mammals. We tested this occurrence by keeping the fertilized oocytes immobilized from the time of ICSI until the first cleavage. Image analysis revealed that cleavage took place preferentially along the short (i.e. equatorial) diameter of the oocyte, thereby partitioning the animal and vegetal materials into the two-cell blastomeres. Our results point to a simple mechanism by which the two sister blastomeres start out as different, rather than becoming different.

Gold(I) N-heterocyclic carbene complexes show strong proapoptotic, antioxidant and anti-inflammatory effects in A2780 and endothelial cells

A series of eight gold(I) N-heterocyclic carbene (NHC) complexes [Au(IMes)(Ln)] based on 1,3-bis(2,4,6-trimethylphenyl)imidazole-2-ylidene (IMes) and 7-azaindole derivatives (HLn), where n = 1-8 for HL1 = 5-fluoro-7-azaindole, HL2 = 5-bromo-7-azaindole, HL3 = 3-chloro-7-azaindole, HL4 = 3-iodo-7-azaindole, HL5 = 5-bromo-3-chloro-7-azaindole, HL6 = 5-bromo-3-iodo-7-azaindole, HL7 = 4-chloro-2-methyl-7-azaindole and HL8 = 7-azaindole, was prepared, characterised and studied for their in vitro anti-cancer and anti-inflammatory effects. The complexes showed significant cytotoxicity on human ovarian cancer cell lines (A2780, IC50 ≈ 8-19 μM and A2780R, IC50 ≈ 8-19 μM) and lowered toxicity in normal HaCat and MRC-5 cells. Cellular effects of the selected complexes 1 and 7 were evaluated in A2780 cells using flow cytometry. Moreover, the time-dependent cellular uptake in A2780 cells, a shotgun proteomic analysis, an ESI-MS study of hydrolysis and interactions with l-cysteine and reduced glutathione (GSH) were performed. Complexes 1 and 7 revealed remarkable anti-inflammatory effects via inhibition of NF-κB activity in human endothelial cells.

Genetic insights into MIS-C Post-COVID-19 in Kuwaiti children: investigating monogenic factors

Background

Multisystem inflammatory syndrome in children (MIS-C) is a severe complication arising from SARS-CoV-2 infection, with indications that rare inborn errors of immunity may play a role in its pathogenesis. Recent studies suggest that genetic predispositions, particularly monogenic forms, could significantly influence the immune responses to SARS-CoV-2 in MIS-C.

Methods

We analysed 24 children under 12 years old, all of whom met the criteria provided by the World Health Organization, 2020 for MIS-C diagnosis, from the Paediatric COVID-19 Registry in Kuwait (PCR-Q8). Demographic and clinical data were collected from medical records, and exome sequencing was performed on the children and their parents to identify rare exonic variants. These variants were prioritized using two approaches: a candidate genes approach employing trio segregation analysis, and a candidate variants approach using a gene panel informed by previous studies on MIS-C-related genetic variants and datasets of differentially expressed genes in MIS-C patients.

Results

The candidate genes approach identified 53 unique genes in 20 of the 24 probands, including DDX60 and TMEM154, which were also differentially expressed between MIS-C and control groups. The candidate variants approach identified 33 rare, predicted deleterious heterozygous variants across 19 unique genes in 19 of the 24 probands, including both previously described and novel candidate variants for MIS-C. Pathway analysis of the identified genes from both approaches revealed significant involvement in immune response, viral defence, and inflammatory pathways.

Conclusion

This study underscores the monogenic susceptibility to MIS-C, enhancing the evidence base through comprehensive genetic analysis. The findings highlight the critical role of genetic predispositions in MIS-C and suggest that further functional genomics work is necessary to explore the mechanistic contributions of these genes, facilitating the development of targeted diagnostic strategies.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Response strategies to acute and chronic environmental stress in the arctic breeding Lapland longspur (Calcarius lapponicus)

The potentially devastating effects of climate change have raised awareness of the need to understand how the biology of wild animals is influenced by extreme-weather events. We investigate how a wild arctic-breeding bird, the Lapland longspur (Calcarius lapponicus), responds to different environmental perturbations and its coping strategies. We explore the transcriptomic response to environmental adversity during the transition from arrival at the breeding grounds to incubation on the Arctic tundra. The effects of an extremely cold spring on arrival and a severe storm during incubation are examined through RNA-seq analysis of pertinent tissues sampled across the breeding cycle. The stress response, circadian rhythms, reproduction, and metabolism are all affected. A key gene of the Hypothalamic-Pituitary-Adrenal axis, FKBP5, was significantly up-regulated in hypothalamus. The genome assembly and gene expression profiles provide comprehensive resources for future studies. Our findings on different coping strategies to chronic and acute stressors will contribute to understanding the interplay between changing environments and genomic regulation.

Exploration of key mechanisms underlying the therapeutic effects of AMD3100 on attenuating lipopolysaccharide-induced acute lung injury in mice

Context

AMD3100, a CXCR4 antagonist, has beneficial effects immaculate in the treatment of acute lung injury (ALI).

Objective

ALI is a severe inflammatory condition associated with poor prognosis and limited treatment options. AMD3100, has therapeutic effects that reduce ALI. Our study explored the regulatory mechanisms of AMD3100 in alleviating the injury of lipopolysaccharide (LPS)-induced ALI in mice.

Materials and Methods

Male ICR mice were randomly divided into control, LPS-treated, AMD3100-treated, and LPS + AMD3100-treatment groups. The histological changes of lung tissues from different groups were evaluated using hematoxylin and eosin staining. Lung injury was measured by ELISA and lung wet/dry ratio. Moreover, lung tissues from the four groups were subjected to transcriptome sequencing followed by differential expression, functional enrichment, protein-protein interaction (PPI) networks, and transcription factor analyses. The validation of mRNAs and protein levels were conducted with qRT-PCR and ELISA.

Results

Hematoxylin and eosin staining combined with the concentration of IL-1 and IL1-β and lung wet/dry ratios revealed that AMD3100 reduced the level of LPS-induced lung injury. Analysis of the transcriptome sequencing data identified 294 differentially expressed genes in the LPS-induced ALI mouse model. Based on the PPI network and module analysis, hub targets of AMD3100, such as Cxcl10 and Cxcl9, were identified in module 1, and hub targets, such as Cxcl12 and Cxcl1, were identified in module 2. Cxcl10 and Cxcl9 are involved in the Toll-like receptor signaling pathway, and Cxcl12 and Cxcl1 arae enriched in the nuclear factor-kappa B signaling pathway. Cxcl19, Cxcl10, and Cxcl1 are targeted by transcription factors like NF-κB. The validation of mRNAs and protein levels conducted by PCR and ELISA supported our transcriptome data.

Conclusions

Our findings indicate that AMD3100 may exhibit a therapeutic effect on LPS-induced ALI in mice by modulating multiple chemokines to inhibit the Toll-like receptor/nuclear factor-kappa B signaling pathway.

Heat shock protein 27 downregulation attenuates isoprenaline-induced myocardial fibrosis and diastolic dysfunction by modulating the endothelial-mesenchymal transition

Heart failure (HF), an end-stage clinical syndrome secondary to cardiac impairment, significantly affects patients' quality of life and long-term prognosis. Myocardial fibrosis leads to systolic and diastolic dysfunction, and promotes the progression of HF. Several studies involving the modulation of myocardial fibrosis have been conducted in an effort to improve cardiac function. Heat shock protein 27 (HSP27) is a small chaperone protein that is overexpressed in cellular stress states. HSP27 modulates epithelial-mesenchymal transition, playing a crucial role in the pathology of several fibrotic diseases. However, its association with myocardial fibrosis regulation is unknown. This study aimed to investigate the mechanisms by which HSP27 regulates myocardial fibrosis. We created cardiac-specific HSP25 (the murine ortholog of human HSP27) knockout mice and found that HSP25 knockdown inhibited endothelial-mesenchymal transition (EndMT), attenuated myocardial fibrosis, and ameliorated diastolic dysfunction in isoproterenol-induced HF mice via echocardiography, histology, and western bloting. In vitro, HSP27 knockdown attenuated transforming growth factor beta-induced EndMT, whereas HSP27 overexpression promoted EndMT. Furthermore, the SMAD3/SNAIL1 pathway was found to be crucial for HSP27-mediated EndMT regulation. As an essential molecule in EndMT regulation and myocardial fibrosis modulation, HSP27 may hold promise as a therapeutic target for patients with HF.

Using clusterProfiler to characterize multiomics data

With the advent of multiomics, software capable of multidimensional enrichment analysis has become increasingly crucial for uncovering gene set variations in biological processes and disease pathways. This is essential for elucidating disease mechanisms and identifying potential therapeutic targets. clusterProfiler stands out for its comprehensive utilization of databases and advanced visualization features. Importantly, clusterProfiler supports various biological knowledge, including Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, through performing over-representation and gene set enrichment analyses. A key feature is that clusterProfiler allows users to choose from various graphical outputs to visualize results, enhancing interpretability. This protocol describes innovative ways in which clusterProfiler has been used for integrating metabolomics and metagenomics analyses, identifying and characterizing transcription factors under stress conditions, and annotating cells in single-cell studies. In all cases, the computational steps can be completed within ~2 min. clusterProfiler is released through the Bioconductor project and can be accessed via https://bioconductor.org/packages/clusterProfiler/ .

Molecular insights into Solanum sisymbriifolium's resistance against Globodera pallida via RNA-seq

BACKGROUND

The presence of potato cyst nematodes (PCN) causes a significant risk to potato crops globally, leading to reduced yields and economic losses. While the plant Solanum sisymbriifolium is known for its resistance to PCN and can be used as a trap crop, the molecular mechanisms behind this resistance remain poorly understood. In this study, genes differentially expressed were identified in control and infected plants during the early stages of the S. sisymbriifolium - G. pallida interaction.

RESULTS

Gene expression profiles were characterized for two S. sisymbriifolium cultivars, Melody and Sis6001, uninfected and infected by G. pallida. The comparative transcriptome analysis revealed a total of 4,087 and 2,043 differentially expressed genes (DEGs) in response to nematode infection in the cultivars Melody and Sis6001, respectively. Gene ontology (GO) enrichment analysis provided insights into the response of the plant to nematode infection, indicating an activation of the plant metabolism, oxidative stress leading to defence mechanism activation, and modification of the plant cell wall. Genes associated with the jasmonic and salicylic acid pathways were also found to be differentially expressed, suggesting their involvement in the plant's defence response. In addition, the analysis of NBS-LRR domain-containing transcripts that play an important role in hypersensitive response and programmed cell death led to the identification of ten transcripts that had no annotations from the databases, with emphasis on TRINITY_DN52667_C1_G1, found to be upregulated in both cultivars.

CONCLUSIONS

These findings represent an important step towards understanding the molecular basis underlying plant resistance to nematodes and facilitating the development of more effective control strategies against PCN.

Astrocyte regional specialization is shaped by postnatal development

Astrocytes are an abundant class of glial cells with critical roles in neural circuit assembly and function. Though many studies have uncovered significant molecular distinctions between astrocytes from different brain regions, how this regionalization unfolds over development is not fully understood. We used single-nucleus RNA sequencing to characterize the molecular diversity of brain cells across six developmental stages and four brain regions in the mouse and marmoset brain. Our analysis of over 170,000 single astrocyte nuclei revealed striking regional heterogeneity among astrocytes, particularly between telencephalic and diencephalic regions, at all developmental time points surveyed in both species. At the stages sampled, most of the region patterning was private to astrocytes and not shared with neurons or other glial types. Though astrocytes were already regionally patterned in late embryonic stages, this region-specific astrocyte gene expression signature changed dramatically over postnatal development, and its composition suggests that regional astrocytes further specialize postnatally to support their local neuronal circuits. Comparing across species, we found divergence in the expression of astrocytic region- and age-differentially expressed genes and the timing of astrocyte maturation relative to birth between mouse and marmoset, as well as hundreds of species differentially expressed genes. Finally, we used expansion microscopy to show that astrocyte morphology is largely conserved across gray matter regions of prefrontal cortex, striatum, and thalamus in the mouse, despite substantial molecular divergence.

Olduvai domain expression downregulates mitochondrial pathways: implications for human brain evolution and neoteny

Encoded by the NBPF gene family, Olduvai (formerly DUF1220) protein domains have undergone the largest human lineage-specific copy number expansion of any coding region in the genome. Olduvai copy number shows a linear relationship with several brain size-related measures and cortical neuron number among primates and with normal and disease-associated (micro- and macrocephaly) variation in brain size in human populations. While Olduvai domains have been shown to promote proliferation of neural stem cells, the mechanism underlying such effects has remained unclear. Here, we investigate the function of Olduvai by transcriptome and proteome analyses of cells overexpressing NBPF1, a gene encoding 7 Olduvai domains. Our results from both RNAseq and mass spectrometry approaches suggest a potential downregulation of mitochondria. In our proteomics study, a Gene Ontology (GO) enrichment analysis for the downregulated proteins revealed a striking overrepresentation of the biological process related to the mitochondrial electron transport chain (p value: 1.81e-11) and identified deregulation of the NADH dehydrogenase activity (p value: 2.43e-11) as the primary molecular function. We verify the reduction of apparent mitochondria via live-cell imaging experiments. Given these and previous Olduvai findings, we suggest that the Olduvai-mediated, dosage-dependent reduction in available energy via mitochondrial downregulation may have resulted in a developmental slowdown such that the neurogenic window among primates, and most extremely in humans, was expanded over a greater time interval, allowing for production of greater numbers of neurons and a larger brain. We further suggest that such a slowdown may extend to other developmental processes that also exhibit neotenic features.

mulea: An R package for enrichment analysis using multiple ontologies and empirical false discovery rate

Traditional gene set enrichment analyses are typically limited to a few ontologies and do not account for the interdependence of gene sets or terms, resulting in overcorrected p-values. To address these challenges, we introduce mulea, an R package offering comprehensive overrepresentation and functional enrichment analysis. mulea employs a progressive empirical false discovery rate (eFDR) method, specifically designed for interconnected biological data, to accurately identify significant terms within diverse ontologies. mulea expands beyond traditional tools by incorporating a wide range of ontologies, encompassing Gene Ontology, pathways, regulatory elements, genomic locations, and protein domains. This flexibility enables researchers to tailor enrichment analysis to their specific questions, such as identifying enriched transcriptional regulators in gene expression data or overrepresented protein domains in protein sets. To facilitate seamless analysis, mulea provides gene sets (in standardised GMT format) for 27 model organisms, covering 22 ontology types from 16 databases and various identifiers resulting in almost 900 files. Additionally, the muleaData ExperimentData Bioconductor package simplifies access to these pre-defined ontologies. Finally, mulea's architecture allows for easy integration of user-defined ontologies, or GMT files from external sources (e.g., MSigDB or Enrichr), expanding its applicability across diverse research areas. mulea is distributed as a CRAN R package downloadable from https://cran.r-project.org/web/packages/mulea/ and https://github.com/ELTEbioinformatics/mulea . It offers researchers a powerful and flexible toolkit for functional enrichment analysis, addressing limitations of traditional tools with its progressive eFDR and by supporting a variety of ontologies. Overall, mulea fosters the exploration of diverse biological questions across various model organisms.

Stromal Expression Profiling Reveals Immune-Driven Adaption to Malignancy in Canine Melanoma Subtypes

Canine mucosal melanoma (CMM) is the most common oral malignancy in dogs and is significantly more aggressive than its cutaneous counterpart (CCM), yet the reasons for this disparity remain unclear. Cancer-associated stroma (CAS) plays a crucial role in tumour progression, but a detailed understanding of CAS in canine melanoma is missing. To assess stromal reprogramming, we analysed CAS from 21 CMM, 14 CCM and normal stroma from 10 skin and 9 oral mucosa samples by laser-capture microdissection followed by RNA sequencing. Results were assessed in relation to subtypes, prognostic factors including mitotic count (MC), ulceration, necrosis, pigmentation and immune cell infiltration (CD3, CD20 and CD68), scored using immunohistochemistry and RNA in situ hybridisation. Stromal reprogramming was evident in both subtypes but significantly more pronounced in CMM. Immune-excluded tumours exhibited higher MC than desert/cold ones. MC strongly correlated with genes associated with B-cells, T-helper cells and CTLA4 in CCM, suggesting CAS reprogramming to depend on tumour malignancy. Finally, we identify an immune-suppressive stromal signature in a subset of CMM characterised by the downregulation of key immune checkpoints and pathways. Together, these findings provide a solid foundation for understanding the role of CAS in canine melanoma, specific to cutaneous and mucosal subtypes.

Systemic Sclerosis Dermal Fibroblast Exosomes Trigger Type 1 Interferon Responses in Keratinocytes via a TBK/JAK/STAT Signaling Axis

OBJECTIVE

Activation of type I interferon (IFN) response has been shown to correlate with disease activity in systemic sclerosis (SSc). It is currently unknown whether the tissue-specific type I IFN activation is a consequence of the response observed in blood or rather its source. Exosomes from SSc fibroblasts were recently shown to activate macrophages in vitro. Here, we aimed to determine the source of type I IFN signature in SSc skin biopsies and the potential role of exosomes from SSc dermal fibroblasts in the process.

METHODS

Skin biopsies were obtained from the forearms of healthy patients and of those with SSc and processed for dermal fibroblasts and keratinocytes. Exosomes were isolated from healthy and SSc dermal fibroblast supernatants by ultracentrifugation and added to human skin keratinocytes. Keratinocyte transcriptome was analyzed by RNA sequencing (RNA-seq) analysis. TANK-binding kinase (TBK) and JAK were inhibited using a small molecule inhibitor (GSK8612) and tofacitinib, respectively.

RESULTS

SSc skin biopsies showed the highest levels of type I IFN response in the epidermal layer. RNA-seq analysis of keratinocytes transcriptome following exposure to dermal fibroblast exosomes showed strong up-regulation of IFN signature genes induced by SSc exosomes compared to healthy control. Inhibition of TBK or JAK activity suppressed the up-regulation of the IFN signature induced by SSc exosomes.

CONCLUSION

IFN activation of SSc keratinocytes is dependent on their crosstalk with dermal fibroblasts and inducible by extracellular exosomes. Our data indicate that SSc fibroblast exosomes contribute to the type I IFN activation in SSc skin through activation of pattern recognition receptors upstream of TBK.

Transcriptomic and proteomic spatial profiling of pediatric and adult diffuse midline glioma H3 K27-Altered

Diffuse midline glioma, H3 K27-altered (DMG) are highly aggressive malignancies of the central nervous system (CNS) that primarily affect the pediatric population. Large scale spatial transcriptomic studies have implicated that tumor microenvironmental landscape plays an important role in determining the phenotypic differences in tumor presentation and clinical course, however, data connecting overall transcriptomic changes to the protein level is lacking. The NanoString GeoMx™ Digital Spatial Profiler platform was used to determine the spatial transcriptomic and proteomic landscape in a cohort of both pediatric and adult H3 K27-altered DMG biopsy samples. Three fluorescently labeled antibodies targeting immune cells (CD45), epithelial cells (PanCK), tumor cells (H3 K27M) and a nucleic acid stain (SYTO-13) were used to establish regions of interest (ROI) for genomic and proteomic analysis. We found genetic alterations within the tumor which can be delineated across patient age and spatial location. We show that the H3 K27M mutation itself has a profound impact on tumor cells transcriptomics and interestingly we found limited fidelity between overall transcriptome and proteome. Our data also validate a previously described genomic signature at the proteomic level and reveal a special shift in the signature based on the local TME composition.

Hyperactivating EZH2 to augment H3K27me3 levels in regulatory T cells enhances immune suppression by driving early effector differentiation

The immunosuppressive function of regulatory T (Treg) cells is essential for maintaining immune homeostasis. Enhancer of zeste homolog 2 (EZH2), a histone H3 lysine 27 (H3K27) methyltransferase, plays a key role in maintaining Treg cell function upon CD28 co-stimulation, and Ezh2 deletion in Treg cells causes autoimmunity. Here, we assess whether increasing H3K27me3 levels, by using an Ezh2Y641F gain-of-function mutation, will improve Treg cell function. We find that Treg cells expressing Ezh2Y641F display an effector Treg phenotype, are poised for improved homing to organ tissues, and can accelerate remission from autoimmunity. The H3K27me3 landscape and transcriptome of naive Ezh2Y641F Treg cells exhibit a redistribution of H3K27me3 modifications that recapitulates the gene expression profile of activated Ezh2WT Treg cells after CD28 co-stimulation. Altogether, increased H3K27me3 levels promote the differentiation of effector Treg cells that can better suppress autoimmunity.

Linked regulation of genome integrity and senescence-associated inflammation by p53

Genomic instability and inflammation are distinct hallmarks of aging, but the connection between them is poorly understood. Understanding their interrelationship will help unravel new mechanisms and therapeutic targets of aging and age-associated diseases. Here we report a novel mechanism directly linking genomic instability and inflammation in senescent cells through a mitochondria-regulated molecular circuit driven by p53 and cytoplasmic chromatin fragments (CCF). We show, through activation or inactivation of p53 by genetic and pharmacologic approaches, that p53 suppresses CCF accumulation and the downstream inflammatory senescence-associated secretory phenotype (SASP), without affecting cell cycle arrest. p53 activation suppressed CCF formation by promoting DNA repair, and this is reflected in maintenance of genomic integrity, particularly in subtelomeric regions, as shown by single cell genome resequencing. Activation of p53 in aged mice by pharmacological inhibition of MDM2 reversed signatures of aging, including age- and senescence-associated transcriptomic signatures of inflammation and age-associated accumulation of monocytes and macrophages in liver. Remarkably, mitochondria in senescent cells suppressed p53 activity by promoting CCF formation and thereby restricting ATM-dependent nuclear DNA damage signaling. These data provide evidence for a mitochondria-regulated p53 signaling circuit in senescent cells that controls DNA repair, genome integrity, and senescence- and age-associated inflammation. This pathway is immunomodulatory in mice and a potential target for healthy aging interventions by small molecules already shown to activate p53.

Behavioral, neurotransmitter and transcriptomic analyses in male and female Fmr1 KO mice

Introduction

Fragile X syndrome is an inherited X-linked disorder associated with intellectual disabilities that begin in childhood and last a lifetime. The symptoms overlap with autism spectrum disorder, and the syndrome predominantly affects males. Consequently, FXS research tends to favor analysis of social behaviors in males, leaving a gap in our understanding of other behavioral traits, especially in females.

Methods

We used a mouse model of FXS to analyze developmental, behavioral, neurochemical, and transcriptomic profiles in males and females.

Results

Our behavioral assays demonstrated locomotor hyperactivity, motor impulsivity, increased "approach" behavior in an approach-avoidance assay, and deficits in nest building behavior. Analysis of brain neurotransmitter content revealed deficits in striatal GABA, glutamate, and serotonin content. RNA sequencing of the ventral striatum unveiled expression changes associated with neurotransmission as well as motivation and substance use pathways. Sex differences were identified in nest building behavior, striatal neurotransmitter content, and ventral striatal gene expression.

Discussion

In summary, our study identified sex differences in specific behavioral, neurotransmitter, and gene expression phenotypes and gene set enrichment analysis identified significant enrichment of pathways associated with motivation and drug reward.

Integration of estimated regional gene expression with neuroimaging and clinical phenotypes at biobank scale

An understanding of human brain individuality requires the integration of data on brain organization across people and brain regions, molecular and systems scales, as well as healthy and clinical states. Here, we help advance this understanding by leveraging methods from computational genomics to integrate large-scale genomic, transcriptomic, neuroimaging, and electronic-health record data sets. We estimated genetically regulated gene expression (gr-expression) of 18,647 genes, across 10 cortical and subcortical regions of 45,549 people from the UK Biobank. First, we showed that patterns of estimated gr-expression reflect known genetic-ancestry relationships, regional identities, as well as inter-regional correlation structure of directly assayed gene expression. Second, we performed transcriptome-wide association studies (TWAS) to discover 1,065 associations between individual variation in gr-expression and gray-matter volumes across people and brain regions. We benchmarked these associations against results from genome-wide association studies (GWAS) of the same sample and found hundreds of novel associations relative to these GWAS. Third, we integrated our results with clinical associations of gr-expression from the Vanderbilt Biobank. This integration allowed us to link genes, via gr-expression, to neuroimaging and clinical phenotypes. Fourth, we identified associations of polygenic gr-expression with structural and functional MRI phenotypes in the Human Connectome Project (HCP), a small neuroimaging-genomic data set with high-quality functional imaging data. Finally, we showed that estimates of gr-expression and magnitudes of TWAS were generally replicable and that the p-values of TWAS were replicable in large samples. Collectively, our results provide a powerful new resource for integrating gr-expression with population genetics of brain organization and disease.

Mitochondrial complex I promotes kidney cancer metastasis

Most kidney cancers are metabolically dysfunctional1-4, but how this dysfunction affects cancer progression in humans is unknown. We infused 13C-labelled nutrients in over 80 patients with kidney cancer during surgical tumour resection. Labelling from [U-13C]glucose varies across subtypes, indicating that the kidney environment alone cannot account for all tumour metabolic reprogramming. Compared with the adjacent kidney, clear cell renal cell carcinomas (ccRCCs) display suppressed labelling of tricarboxylic acid (TCA) cycle intermediates in vivo and in ex vivo organotypic cultures, indicating that suppressed labelling is tissue intrinsic. [1,2-13C]acetate and [U-13C]glutamine infusions in patients, coupled with measurements of respiration in isolated human kidney and tumour mitochondria, reveal lower electron transport chain activity in ccRCCs that contributes to decreased oxidative and enhanced reductive TCA cycle labelling. However, ccRCC metastases unexpectedly have enhanced TCA cycle labelling compared with that of primary ccRCCs, indicating a divergent metabolic program during metastasis in patients. In mice, stimulating respiration or NADH recycling in kidney cancer cells is sufficient to promote metastasis, whereas inhibiting electron transport chain complex I decreases metastasis. These findings in humans and mice indicate that metabolic properties and liabilities evolve during kidney cancer progression, and that mitochondrial function is limiting for metastasis but not growth at the original site.

Robust assessment of sample preparation protocols for proteomics of cells and tissues

In proteomic studies, the reliability and reproducibility of results hinge on well-executed protein extraction and digestion protocols. Here, we systematically compared three established digestion methods for macrophages, namely filter-assisted sample preparation (FASP), in-solution, and in-gel digestion protocols. We also compared lyophilization and manual lysis for liver tissue protein extraction, each of them tested using either sodium deoxycholate (SDC)- or RIPA-based lysis buffer. For the macrophage cell line, FASP using passivated filter units outperformed the other tested methods regarding the number of identified peptides and proteins. However, a careful standardization has shown that all three methods can yield robust results across a wide range of starting material (even starting with 1 μg of proteins). Importantly, inter and intra-day coefficients of variance (CVs) were determined for all sample preparation protocols. Thus, the median inter-day CVs for in-solution, in-gel and FASP protocols were respectively 10, 8 and 9%, very similar to the median CVs obtained for the intra-day analysis (9, 8 and 8%, respectively). Moreover, FASP digestion presented 80% of proteins with a CV lower than 25%, followed closely by in-gel digestion (78%) and in-solution sample preparation (72%) protocols. For tissue proteomics, both manual lysis and lyophilization presented similar proteome coverage and reproducibility, but the efficiency of protein extraction depended on the lysis buffer used, with RIPA buffer showing better results. In conclusion, although each sample preparation method has its own particularity, they are all suited for successful proteomic experiments if a careful standardization of the sample preparation workflow is carried out.

Genome-Wide Analysis of Genetic Diversity and Selection Signatures in Zaobei Beef Cattle

This investigation provides a comprehensive analysis of genomic diversity and selection signatures in Zaobei beef cattle, an indigenous breed known for its adaptation to hot and humid climates and superior meat quality. Whole-genome resequencing was conducted on 23 Zaobei cattle, compared with 46 Simmental cattle to highlight genetic distinctions. Population structure analysis confirmed the genetic uniqueness of Zaobei cattle. Using methods such as DASDC v1.01, XPEHH, and θπ ratio, we identified 230, 232, and 221 genes through DASDC, including hard sweeps, soft sweeps, and linkage sweeps, respectively. Coincidentally, 109 genes were identified when using XPEHH and θπ ratio methods. Together, these analyses revealed eight positive selection genes (ARHGAP15, ZNF618, USH2A, PDZRN4, SPATA6, ROR2, KCNIP3, and VWA3B), which are linked to critical traits such as heat stress adaptation, fertility, and meat quality. Moreover, functional enrichment analyses showed pathways related to autophagy, immune response, energy metabolism, and muscle development. The comprehensive genomic insights gained from this study provide valuable knowledge for breeding programs aimed at enhancing the beneficial traits in Zaobei cattle.

Identification and validation of Rab GTPases RAB13 as biomarkers for peritoneal metastasis and immune cell infiltration in colorectal cancer patients

Background

As one of the most common cancer, colorectal cancer (CRC) is with high morbidity and mortality. Peritoneal metastasis (PM) is a fatal state of CRC, and few patients may benefit from traditional therapies. There is a complex interaction between PM and immune cell infiltration. Therefore, we aimed to determine biomarkers associated with colorectal cancer peritoneal metastasis (CRCPM) and their relationship with immune cell infiltration.

Methods

By informatic analysis, differently expressed genes (DEGs) were selected and hub genes were screened out. RAB13, one of the hub genes, was identificated from public databases and validated in CRC tissues. The ESTIMATE, CEBERSORT and TIMER algorithms were applied to analyze the correlation between RAB13 and immune infiltration in CRC. RAB13's expression in different cells were analyzed at the single-cell level in scRNA-Seq. The Gene Set Enrichment Analysis (GSEA) was performed for RAB13 enrichment and further confirmed. Using oncoPredict algorithm, RAB13's impact on drug sensitivity was evaluated.

Results

High RAB13 expression was identified in public databases and led to a poor prognosis. RAB13 was found to be positively correlated with the macrophages and other immune cells infiltration and from scRNA-Seq, RAB13 was found to be located in CRC cells and macrophages. GSEA revealed that high RAB13 expression enriched in a various of biological signaling, and oncoPredict algorithm showed that RAB13 expression was correlated with paclitaxel sensitivity.

Conclusion

Our study indicated clinical role of RAB13 in CRC-PM, suggesting its potential as a therapeutic target in the future.

A Comprehensive Review of Bioinformatics Tools for Genomic Biomarker Discovery Driving Precision Oncology

The rapid advancement of high-throughput technologies, particularly next-generation sequencing (NGS), has revolutionized cancer research by enabling the investigation of genetic variations such as SNPs, copy number variations, gene expression, and protein levels. These technologies have elevated the significance of precision oncology, creating a demand for biomarker identification and validation. This review explores the complex interplay of oncology, cancer biology, and bioinformatics tools, highlighting the challenges in statistical learning, experimental validation, data processing, and quality control that underpin this transformative field. This review outlines the methodologies and applications of bioinformatics tools in cancer genomics research, encompassing tools for data structuring, pathway analysis, network analysis, tools for analyzing biomarker signatures, somatic variant interpretation, genomic data analysis, and visualization tools. Open-source tools and repositories like The Cancer Genome Atlas (TCGA), Genomic Data Commons (GDC), cBioPortal, UCSC Genome Browser, Array Express, and Gene Expression Omnibus (GEO) have emerged to streamline cancer omics data analysis. Bioinformatics has significantly impacted cancer research, uncovering novel biomarkers, driver mutations, oncogenic pathways, and therapeutic targets. Integrating multi-omics data, network analysis, and advanced ML will be pivotal in future biomarker discovery and patient prognosis prediction.

CREB-regulated transcription during glycogen synthesis in astrocytes

Glycogen storage, conversion and utilization in astrocytes play an important role in brain energy metabolism. The conversion of glycogen to lactate through glycolysis occurs through the coordinated activities of various enzymes and inhibition of this process can impair different brain processes including formation of long-lasting memories. To replenish depleted glycogen stores, astrocytes undergo glycogen synthesis, a cellular process that has been shown to require transcription and translation during specific stimulation paradigms. However, the detail nuclear signaling mechanisms and transcriptional regulation during glycogen synthesis in astrocytes remains to be explored. In this report, we study the molecular mechanisms of vasoactive intestinal peptide (VIP)-induced glycogen synthesis in astrocytes. VIP is a potent neuropeptide that triggers glycogenolysis followed by glycogen synthesis in astrocytes. We show evidence that VIP-induced glycogen synthesis requires CREB-mediated transcription that is calcium dependent and requires conventional Protein Kinase C but not Protein Kinase A. In parallel to CREB activation, we demonstrate that VIP also triggers nuclear accumulation of the CREB coactivator CRTC2 in astrocytic nuclei. Transcriptome profiles of VIP-induced astrocytes identified robust CREB transcription, including a subset of genes linked to glucose and glycogen metabolism. Finally, we demonstrate that VIP-induced glycogen synthesis shares similar as well as distinct molecular signatures with glucose-induced glycogen synthesis, including the requirement of CREB-mediated transcription. Overall, our data demonstrates the importance of CREB-mediated transcription in astrocytes during stimulus-driven glycogenesis.

Influence of antipsychotic drugs on microRNA expression in schizophrenia patients - A systematic review

Schizophrenia (SCZ) is a severe psychiatric disorder with unclear pathophysiology. Moreover, there is no specific biological marker to help clinicians to define a diagnosis, and medication is decided according to the psychiatrist's experience. In this scenario, microRNAs (miRNAs), which are small noncoding RNA molecules that regulate several genes, emerge as potential peripheral biomarkers to help not only the evaluation of the disease state but also the treatment response. Here, we systematically reviewed indexed literature and evaluated follow-up studies investigating the changes in miRNA expression due to antipsychotic treatment. We also assessed target genes and performed pathway enrichment analysis of miRNAs listed in this systematic review. A total of 11 studies were selected according to research criteria, and we observed that 28 miRNAs play a relevant role in schizophrenia pathogenesis or response to antipsychotic treatment, seven of those of extreme interest as possible biomarkers either for condition or treatment. Predicted targets of the miRNAs reviewed here were previously associated with schizophrenia in genome-wide studies, and pathway analysis showed enrichment for genes related to neural processes. With this review, we expect to highlight the importance of miRNAs in schizophrenia pathogenesis and its treatment and point out interesting miRNAs to future studies.

Identifying Universal Fish Biomarker Genes in Response to PCB126 Exposure by Comparative Transcriptomic Analyses

Water pollution remains a major environmental concern, with increased toxic by-products being released into water bodies. Many of these chemical contaminants persist in the environment and bio-accumulate in aquatic organisms. At present, toxicological tests are mostly based on laboratory tests, and effective methods for monitoring wild aquatic environments remain lacking. In the present study, we used a well-characterized toxic chemical, 3,3',4,4',5-polychlorinated biphenyl (PCB126), as an example to try to identify common biomarker genes to be used for predictive toxicity of this toxic substance. First, we used two laboratory fish models, the zebrafish (Danio rerio) and medaka (Oryzias latipes), to expose PCB126 to obtain liver transcriptomic data by RNA-seq. Comparative transcriptomic analyses indicated generally conserved and concerted changes from the two species, thus validating the transcriptomic data for biomarker gene selection. Based on the common up- and downregulated genes in the two species, we selected nine biomarker genes to further test in other fish species. The first validation experiment was carried out using the third fish species, Mozambique tilapia (Oreochromis mossambicus), and essentially, all these biomarker genes were validated for consistent responses with the two laboratory fish models. Finally, to develop universal PCR primers suitable for potentially all teleost fish species, we designed degenerate primers and tested them in the three fish species as well as in another fish species without a genomic sequence available: guppy (Poecilia reticulata). We found all the biomarker genes showed consistent response to PCB126 exposure in at least 50% of the species. Thus, our study provides a promising strategy to identify common biomarker genes to be used for teleost fish analyses. By using degenerate PCR primers and analyzing multiple biomarker genes, it is possible to develop diagnostic PCR arrays to predict water contamination from any wild fish species sampled in different water bodies.

Co-analysis of methylation platforms for signatures of biological aging in the domestic dog reveals previously unexplored confounding factors

Chronological age reveals the number of years an individual has lived since birth. By contrast, biological age varies between individuals of the same chronological age at a rate reflective of physiological decline. Differing rates of physiological decline are related to longevity and result from genetics, environment, behavior, and disease. The creation of methylation biological age predictors is a long-standing challenge in aging research due to the lack of individual pre-mortem longevity data. The consistent differences in longevity between domestic dog breeds enable the construction of biological age estimators which can, in turn, be contrasted with methylation measurements to elucidate mechanisms of biological aging. We draw on three flagship methylation studies using distinct measurement platforms and tissues to assess the feasibility of creating biological age methylation clocks in the dog. We expand epigenetic clock building strategies to accommodate phylogenetic relationships between individuals, thus controlling for the use of breed standard metrics. We observe that biological age methylation clocks are affected by population stratification and require heavy parameterization to achieve effective predictions. Finally, we observe that methylation-related markers reflecting biological age signals are rare and do not colocalize between datasets.

WebGestalt 2024: faster gene set analysis and new support for metabolomics and multi-omics

Enrichment analysis, crucial for interpreting genomic, transcriptomic, and proteomic data, is expanding into metabolomics. Furthermore, there is a rising demand for integrated enrichment analysis that combines data from different studies and omics platforms, as seen in meta-analysis and multi-omics research. To address these growing needs, we have updated WebGestalt to include enrichment analysis capabilities for both metabolites and multiple input lists of analytes. We have also significantly increased analysis speed, revamped the user interface, and introduced new pathway visualizations to accommodate these updates. Notably, the adoption of a Rust backend reduced gene set enrichment analysis time by 95% from 270.64 to 12.41 s and network topology-based analysis by 89% from 159.59 to 17.31 s in our evaluation. This performance improvement is also accessible in both the R package and a newly introduced Python package. Additionally, we have updated the data in the WebGestalt database to reflect the current status of each source and have expanded our collection of pathways, networks, and gene signatures. The 2024 WebGestalt update represents a significant leap forward, offering new support for metabolomics, streamlined multi-omics analysis capabilities, and remarkable performance enhancements. Discover these updates and more at https://www.webgestalt.org.

Age, sex and Alzheimer's disease: a longitudinal study of 3xTg-AD mice reveals sex-specific disease trajectories and inflammatory responses mirrored in postmortem brains from Alzheimer's patients

BACKGROUND

Aging and sex are major risk factors for developing late-onset Alzheimer's disease. Compared to men, women experience worse neuropathological burden and cognitive decline despite living longer with the disease. Similarly, male 3xTg-AD mice, developed to model Alzheimer's disease, no longer consistently exhibit standard Alzheimer's neuropathology yet experience higher rates of mortality - providing a unique opportunity to further elucidate this dichotomy. We hypothesized that sex differences in the biological aging process yield distinct pathological and molecular Alzheimer's disease signatures in males and females, which could be harnessed for therapeutic and biomarker development.

METHODS

We aged male and female, 3xTg-AD and B6129 control mice across their respective lifespans (n = 3-8 mice per sex, strain, and age group) and longitudinally assessed neuropathological hallmarks of Alzheimer's disease, markers of hepatic inflammation, splenic mass and morphology, as well as plasma cytokine levels. We conducted RNA sequencing analysis on bulk brain tissue and examined differentially expressed genes (DEGs) between 3xTg-AD and B6129 samples and across ages in each sex. We also examined DEGs between clinical Alzheimer's and control parahippocampal gyrus brain tissue samples from the Mount Sinai Brain Bank study in each sex.

RESULTS

3xTg-AD females significantly outlived 3xTg-AD males and exhibited progressive Alzheimer's neuropathology, while 3xTg-AD males demonstrated progressive hepatic inflammation, splenomegaly, circulating inflammatory proteins, and minimal Alzheimer's neuropathological hallmarks. Instead, 3xTg-AD males experienced an accelerated upregulation of immune-related gene expression in the brain relative to females. Our clinical investigations revealed that individuals with Alzheimer's disease develop similar sex-specific alterations in neuronal and immune function. In diseased males of both species, we observed greater upregulation of complement-related gene expression, and lipopolysaccharide was predicted as the top upstream regulator of DEGs.

CONCLUSIONS

Our data demonstrate that chronic inflammation and complement activation are associated with increased mortality, indicating that age-related changes in immune response contribute to sex differences in Alzheimer's disease trajectories. We provide evidence that aging and transgene-driven disease progression trigger a widespread inflammatory response in 3xTg-AD males, which mimics the impact of lipopolysaccharide stimulation despite the absence of infection.

Human serum proteomics reveals a molecular signature after one night of sleep deprivation

Study Objectives

Sleep deprivation is highly prevalent and caused by conditions such as night shift work or illnesses like obstructive sleep apnea. Compromised sleep affects cardiovascular-, immune-, and neuronal systems. Recently, we published human serum proteome changes after a simulated night shift. This pilot proteomic study aimed to further explore changes in human blood serum after 6 hours of sleep deprivation at night.

Methods

Human blood serum samples from eight self-declared healthy females were analyzed using Orbitrap Eclipse mass spectrometry (MS-MS) and high-pressure liquid chromatography. We used a within-participant design, in which the samples were taken after 6 hours of sleep at night and after 6 hours of sleep deprivation the following night. Systems biological databases and bioinformatic software were used to analyze the data and comparative analysis were done with other published sleep-related proteomic datasets.

Results

Out of 494 proteins, 66 were found to be differentially expressed proteins (DEPs) after 6 hours of sleep deprivation. Functional enrichment analysis revealed the associations of these DEPs with several biological functions related to the altered regulation of cellular processes such as platelet degranulation and blood coagulation, as well as associations with different curated gene sets.

Conclusions

This study presents serum proteomic changes after 6 hours of sleep deprivation, supports previous findings showing that short sleep deprivation affects several biological processes, and reveals a molecular signature of proteins related to pathological conditions such as altered coagulation and platelet function, impaired lipid and immune function, and cell proliferation. Data are available via ProteomeXchange with identifier PXD045729. This paper is part of the Genetic and other molecular underpinnings of sleep, sleep disorders, and circadian rhythms including translational approaches Collection.

A chromosome-level genome assembly of the Asian house martin implies potential genes associated with the feathered-foot trait

The presence of feathers is a vital characteristic among birds, yet most modern birds had no feather on their feet. The discoveries of feathers on the hind limbs of basal birds and dinosaurs have sparked an interest in the evolutionary origin and genetic mechanism of feathered feet. However, the majority of studies investigating the genes associated with this trait focused on domestic populations. Understanding the genetic mechanism underpinned feathered-foot development in wild birds is still in its infancy. Here, we assembled a chromosome-level genome of the Asian house martin (Delichon dasypus) using the long-read High Fidelity sequencing approach to initiate the search for genes associated with its feathered feet. We employed the whole-genome alignment of D. dasypus with other swallow species to identify high-SNP regions and chromosomal inversions in the D. dasypus genome. After filtering out variations unrelated to D. dasypus evolution, we found six genes related to feather development near the high-SNP regions. We also detected three feather development genes in chromosomal inversions between the Asian house martin and the barn swallow genomes. We discussed their association with the wingless/integrated (WNT), bone morphogenetic protein, and fibroblast growth factor pathways and their potential roles in feathered-foot development. Future studies are encouraged to utilize the D. dasypus genome to explore the evolutionary process of the feathered-foot trait in avian species. This endeavor will shed light on the evolutionary path of feathers in birds.

Tracking of Nascent Matrix Deposition during Muscle Stem Cell Activation across Lifespan Using Engineered Hydrogels

Adult stem cells occupy a niche that contributes to their function, but how stem cells rebuild their microenvironment after injury remains an open-ended question. Herein, biomaterial-based systems and metabolic labeling are utilized to evaluate how skeletal muscle stem cells deposit extracellular matrix. Muscle stem cells and committed myoblasts are observed to generate less nascent matrix than muscle resident fibro-adipogenic progenitors. When cultured on substrates that matched the stiffness of physiological uninjured and injured muscles, muscle stem cells increased nascent matrix deposition with activation kinetics. Reducing the ability to deposit nascent matrix by an inhibitor of vesicle trafficking (Exo-1) attenuated muscle stem cell function and mimicked impairments observed from muscle stem cells isolated from old muscles. Old muscle stem cells are observed to deposit less nascent matrix than young muscle stem cells, which is rescued with therapeutic supplementation of insulin-like growth factors. These results highlight the role of nascent matrix production with muscle stem cell activation.

Detecting the effect of genetic diversity on brain composition in an Alzheimer's disease mouse model

Alzheimer's disease (AD) is broadly characterized by neurodegeneration, pathology accumulation, and cognitive decline. There is considerable variation in the progression of clinical symptoms and pathology in humans, highlighting the importance of genetic diversity in the study of AD. To address this, we analyze cell composition and amyloid-beta deposition of 6- and 14-month-old AD-BXD mouse brains. We utilize the analytical QUINT workflow- a suite of software designed to support atlas-based quantification, which we expand to deliver a highly effective method for registering and quantifying cell and pathology changes in diverse disease models. In applying the expanded QUINT workflow, we quantify near-global age-related increases in microglia, astrocytes, and amyloid-beta, and we identify strain-specific regional variation in neuron load. To understand how individual differences in cell composition affect the interpretation of bulk gene expression in AD, we combine hippocampal immunohistochemistry analyses with bulk RNA-sequencing data. This approach allows us to categorize genes whose expression changes in response to AD in a cell and/or pathology load-dependent manner. Ultimately, our study demonstrates the use of the QUINT workflow to standardize the quantification of immunohistochemistry data in diverse mice, - providing valuable insights into regional variation in cellular load and amyloid deposition in the AD-BXD model.

DExplore: An Online Tool for Detecting Differentially Expressed Genes from mRNA Microarray Experiments

Microarray experiments, a mainstay in gene expression analysis for nearly two decades, pose challenges due to their complexity. To address this, we introduce DExplore, a user-friendly web application enabling researchers to detect differentially expressed genes using data from NCBI's GEO. Developed with R, Shiny, and Bioconductor, DExplore integrates WebGestalt for functional enrichment analysis. It also provides visualization plots for enhanced result interpretation. With a Docker image for local execution, DExplore accommodates unpublished data. To illustrate its utility, we showcase two case studies on cancer cells treated with chemotherapeutic drugs. DExplore streamlines microarray data analysis, empowering molecular biologists to focus on genes of biological significance.

Homeobox regulator Wilms Tumour 1 is displaced by androgen receptor at cis-regulatory elements in the endometrium of PCOS patients

Decidualisation, the process whereby endometrial stromal cells undergo morphological and functional transformation in preparation for trophoblast invasion, is often disrupted in women with polycystic ovary syndrome (PCOS) resulting in complications with pregnancy and/or infertility. The transcription factor Wilms tumour suppressor 1 (WT1) is a key regulator of the decidualization process, which is reduced in patients with PCOS, a complex condition characterized by increased expression of androgen receptor in endometrial cells and high presence of circulating androgens. Using genome-wide chromatin immunoprecipitation approaches on primary human endometrial stromal cells, we identify key genes regulated by WT1 during decidualization, including homeobox transcription factors which are important for regulating cell differentiation. Furthermore, we found that AR in PCOS patients binds to the same DNA regions as WT1 in samples from healthy endometrium, suggesting dysregulation of genes important to decidualisation pathways in PCOS endometrium due to competitive binding between WT1 and AR. Integrating RNA-seq and H3K4me3 and H3K27ac ChIP-seq metadata with our WT1/AR data, we identified a number of key genes involved in immune response and angiogenesis pathways that are dysregulated in PCOS patients. This is likely due to epigenetic alterations at distal enhancer regions allowing AR to recruit cofactors such as MAGEA11, and demonstrates the consequences of AR disruption of WT1 in PCOS endometrium.

SRRM2 splicing factor modulates cell fate in early development

Embryo development is an orchestrated process that relies on tight regulation of gene expression to guide cell differentiation and fate decisions. The Srrm2 splicing factor has recently been implicated in developmental disorders and diseases, but its role in early mammalian development remains unexplored. Here, we show that Srrm2 dosage is critical for maintaining embryonic stem cell pluripotency and cell identity. Srrm2 heterozygosity promotes loss of stemness, characterised by the coexistence of cells expressing naive and formative pluripotency markers, together with extensive changes in gene expression, including genes regulated by serum-response transcription factor (SRF) and differentiation-related genes. Depletion of Srrm2 by RNA interference in embryonic stem cells shows that the earliest effects of Srrm2 heterozygosity are specific alternative splicing events on a small number of genes, followed by expression changes in metabolism and differentiation-related genes. Our findings unveil molecular and cellular roles of Srrm2 in stemness and lineage commitment, shedding light on the roles of splicing regulators in early embryogenesis, developmental diseases and tumorigenesis.

SIEVE: One-stop differential expression, variability, and skewness analyses using RNA-Seq data

Motivation

RNA-Seq data analysis is commonly biased towards detecting differentially expressed genes and insufficiently conveys the complexity of gene expression changes between biological conditions. This bias arises because discrete models of RNA-Seq count data cannot fully characterize the mean, variance, and skewness of gene expression distribution using independent model parameters. A unified framework that simultaneously tests for differential expression, variability, and skewness is needed to realize the full potential of RNA-Seq data analysis in a systems biology context.

Results

We present SIEVE, a statistical methodology that provides the desired unified framework. SIEVE embraces a compositional data analysis framework that transforms discrete RNA-Seq counts to a continuous form with a distribution that is well-fitted by a skew-normal distribution. Simulation results show that SIEVE controls the false discovery rate and probability of Type II error better than existing methods for differential expression analysis. Analysis of the Mayo RNA-Seq dataset for Alzheimer's disease using SIEVE reveals that a gene set with significant expression difference in mean, standard deviation and skewness between the control and the Alzheimer's disease group strongly predicts a subject's disease state. Furthermore, functional enrichment analysis shows that relying solely on differentially expressed genes detects only a segment of a much broader spectrum of biological aspects associated with Alzheimer's disease. The latter aspects can only be revealed using genes that show differential variability and skewness. Thus, SIEVE enables fresh perspectives for understanding the intricate changes in gene expression that occur in complex diseases.

Availability

The SIEVE R package and source codes are available at https://github.com/Divo-Lee/SIEVE .

MAGP1 maintains tumorigenicity and angiogenesis of laryngeal cancer by activating Wnt/β-catenin/MMP7 pathway

Microfibril-associated glycoprotein-1 (MAGP1), a crucial extracellular matrix protein, contributes to the initiation and progression of different cancers. However, the role of MAGP1 in laryngeal cancer is not clear. The purpose of this study was to investigate the clinical significance and biological function of MAGP1 in laryngeal cancer. MAGP1 was upregulated in public databases and laryngeal cancer tissues, and high MAGP1 expression led to a poor prognosis and was identified as an independent prognostic marker. Knocking-down MAGP1 inhibited laryngeal cancer cell growth and metastasis. According to gene set enrichment analysis, high MAGP1 expression revealed enrichment in Wnt/β-catenin signaling and knocking-down MAGP1 in laryngeal cancer cells also caused degradation, de-activation, re-location and loss of stability of β-catenin. Additionally, we observed MAGP1 in laryngeal cancer cells inhibits angiogenesis in an MMP7-dependent way. In conclusion, our study suggests a clinical role of MAGP1 in laryngeal cancer, signifying its potential as a therapeutic target in the future.

Increased EZH2 function in regulatory T cells promotes their capacity to suppress autoimmunity by driving effector differentiation prior to activation

The immunosuppressive function of regulatory T (Treg) cells is essential for maintaining immune homeostasis. Enhancer of zeste homolog 2 (EZH2), a histone H3 lysine 27 (H3K27) methyltransferase, plays a key role in maintaining Treg cell function upon CD28 co-stimulation, and Ezh2 deletion in Treg cells causes autoimmunity. Here we assessed whether increased EZH2 activity in Treg cells would improve Treg cell function. Using an Ezh2 gain-of-function mutation, Ezh2 Y641F , we found that Treg cells expressing Ezh2 Y641F displayed an increased effector Treg phenotype and were poised for improved homing to organ tissues. Expression of Ezh2 Y641F in Treg cells led to more rapid remission from autoimmunity. H3K27me3 profiling and transcriptomic analysis revealed a redistribution of H3K27me3, which prompted a gene expression profile in naïve Ezh2 Y641F Treg cells that recapitulated aspects of CD28-activated Ezh2 WT Treg cells. Altogether, increased EZH2 activity promotes the differentiation of effector Treg cells that can better suppress autoimmunity.

Highlights

EZH2 function promotes effector differentiation of Treg cells.EZH2 function promotes Treg cell migration to organ tissues.EZH2 function in Treg cells improves remission from autoimmunity.EZH2 function poises naïve Treg cells to adopt a CD28-activated phenotype.

Transcriptomic and Proteomic Spatial Profiling of Pediatric and Adult Diffuse Midline Glioma H3 K27-Altered, Reveals Region Specific Differences and Limited Overlap between mRNA and Protein

Diffuse midline glioma, H3 K27-altered (DMG-Alt) are highly aggressive malignancies of the central nervous system (CNS) that primarily affect the pediatric population. Large scale spatial transcriptomic studies have implicated that tumor microenvironmental landscape plays an important role in determining the phenotypic differences in tumor presentation and clinical course, however, data connecting overall transcriptomic changes to the protein level is lacking. The NanoString GeoMx™ Digital Spatial Profiler platform was used to determine the spatial transcriptomic and proteomic landscape in a cohort of both pediatric and adult H3 K27-altered DMG biopsy samples. Three fluorescently labeled antibodies targeting immune cells (CD45), epithelial cells (PanCK), tumor cells (H3 K27M) and a nucleic acid stain (SYTO-13) were used to establish regions of interest (ROI) for genomic and proteomic analysis. We found genetic alterations within the tumor which can be delineated across patient age and spatial location. We show that the H3 K27M mutation itself has a profound impact on tumor cells transcriptomics and interestingly we found limited fidelity between overall transcriptome and proteome. Our data also validate the previously described OPC like precursor signature at the proteomic level and reveal a special shift in the signature based on the local TME composition.

Genetic contribution of reproductive traits to risk of uterine leiomyomata: a large-scale, genome-wide, cross-trait analysis

BACKGROUND

Although phenotypic associations between female reproductive characteristics and uterine leiomyomata have long been observed in epidemiologic investigations, the shared genetic architecture underlying these complex phenotypes remains unclear.

OBJECTIVE

We aimed to investigate the shared genetic basis, pleiotropic effects, and potential causal relationships underlying reproductive traits (age at menarche, age at natural menopause, and age at first birth) and uterine leiomyomata.

STUDY DESIGN

With the use of large-scale, genome-wide association studies conducted among women of European ancestry for age at menarche (n=329,345), age at natural menopause (n=201,323), age at first birth (n=418,758), and uterine leiomyomata (ncases/ncontrols=35,474/267,505), we performed a comprehensive, genome-wide, cross-trait analysis to examine systematically the common genetic influences between reproductive traits and uterine leiomyomata.

RESULTS

Significant global genetic correlations were identified between uterine leiomyomata and age at menarche (rg, -0.17; P=3.65×10-10), age at natural menopause (rg, 0.23; P=3.26×10-07), and age at first birth (rg, -0.16; P=1.96×10-06). Thirteen genomic regions were further revealed as contributing significant local correlations (P<.05/2353) to age at natural menopause and uterine leiomyomata. A cross-trait meta-analysis identified 23 shared loci, 3 of which were novel. A transcriptome-wide association study found 15 shared genes that target tissues of the digestive, exo- or endocrine, nervous, and cardiovascular systems. Mendelian randomization suggested causal relationships between a genetically predicted older age at menarche (odds ratio, 0.88; 95% confidence interval, 0.85-0.92; P=1.50×10-10) or older age at first birth (odds ratio, 0.95; 95% confidence interval, 0.90-0.99; P=.02) and a reduced risk for uterine leiomyomata and between a genetically predicted older age at natural menopause and an increased risk for uterine leiomyomata (odds ratio, 1.08; 95% confidence interval, 1.06-1.09; P=2.30×10-27). No causal association in the reverse direction was found.

CONCLUSION

Our work highlights that there are substantial shared genetic influences and putative causal links that underlie reproductive traits and uterine leiomyomata. The findings suggest that early identification of female reproductive risk factors may facilitate the initiation of strategies to modify potential uterine leiomyomata risk.

Experimentally Induced Peripheral Venous Congestion Exacerbates Inflammation, Oxidative Stress, and Neurohormonal and Endothelial Cell Activation in Patients With Systolic Heart Failure

BACKGROUND

Venous congestion (VC) is a hallmark of symptomatic heart failure (HF) requiring hospitalization; however, its role in the pathogenesis of HF progression remains unclear. We investigated whether peripheral VC exacerbates inflammation, oxidative stress and neurohormonal and endothelial cell (EC) activation in patients with HF with reduced ejection fraction (HFrEF).

METHODS AND RESULTS

Two matched groups of patients with HFrEF and with no peripheral VC vs without recent HF hospitalization were studied. We modeled peripheral VC by inflating a cuff around the dominant arm, targeting ∼ 30 mmHg increase in venous pressure (venous stress test [VST]). Blood and ECs were sampled before and after 90 minutes of VST. We studied 44 patients (age 53 ± 12 years, 32% female). Circulating endothelin-1, tumor necrosis factor-α, interleukin-6, isoprostane, angiotensin II (ang-2), angiopoietin-2, vascular cell adhesion molecule-1, and CD146 significantly increased after the VST. Enhanced endothelin-1 and angiopoietin-2 responses to the VST were present in patients with vs without recent hospitalization and were prospectively associated with incident HF-related events; 6698 messenger ribonucleic acid (mRNA probe sets were differentially expressed in ECs after VST.

CONCLUSIONS

Experimental VC exacerbates inflammation, oxidative stress, neurohormonal and EC activation and promotes unfavorable transcriptome remodeling in ECs of patients with HFrEF. A distinct biological sensitivity to VC appears to be associated with high risk for HF progression.

Bulk RNA sequencing for analysis of post COVID-19 condition in adolescents and young adults

BACKGROUND

Post COVID-19 condition (PCC) is a complication of SARS-COV-2 infection and can lead to long-term disability.

METHODS

The present study was designed to analyse the gene expression patterns of PCC through bulk RNA sequencing of whole blood and to explore the potential molecular mechanisms of PCC. Whole blood was collected from 80 participants enrolled in a prospective cohort study following SARS-CoV-2 infected and non-infected individuals for 6 months after recruitment and was used for bulk RNA sequencing. Identification of differentially expressed genes (DEG), pathway enrichment and immune cell deconvolution was performed to explore potential biological pathways involved in PCC.

RESULTS

We have found 13 differentially expressed genes associated with PCC. Enriched pathways were related to interferon-signalling and anti-viral immune processes.

CONCLUSION

The PCC transcriptome is characterized by a modest overexpression of interferon-stimulated genes, pointing to a subtle ongoing inflammatory response.

Loss of GTF2I promotes neuronal apoptosis and synaptic reduction in human cellular models of neurodevelopment

Individuals with Williams syndrome (WS), a neurodevelopmental disorder caused by hemizygous loss of 26-28 genes at 7q11.23, characteristically portray a hypersocial phenotype. Copy-number variations and mutations in one of these genes, GTF2I, are associated with altered sociality and are proposed to underlie hypersociality in WS. However, the contribution of GTF2I to human neurodevelopment remains poorly understood. Here, human cellular models of neurodevelopment, including neural progenitors, neurons, and three-dimensional cortical organoids, are differentiated from CRISPR-Cas9-edited GTF2I-knockout (GTF2I-KO) pluripotent stem cells to investigate the role of GTF2I in human neurodevelopment. GTF2I-KO progenitors exhibit increased proliferation and cell-cycle alterations. Cortical organoids and neurons demonstrate increased cell death and synaptic dysregulation, including synaptic structural dysfunction and decreased electrophysiological activity on a multielectrode array. Our findings suggest that changes in synaptic circuit integrity may be a prominent mediator of the link between alterations in GTF2I and variation in the phenotypic expression of human sociality.

Identification of expression profiles and prognostic value of RFCs in colorectal cancer

Colorectal cancer (CRC) ranks among the most prevalent cancers globally, with its incidence closely tied to DNA damage. The Replication Factor C (RFC) complexes comprises five protein subunits: RFC1, RFC2, RFC3, RFC4, and RFC5. These RFC complexes play crucial roles in DNA replication, repair pathways, activities post DNA damage, and ATP-dependent processes during DNA synthesis. However, the impact of RFC complexes proteins on CRC prognosis remains unclear. To explore this, we employed a computational analysis approach, utilizing platforms such as the DepMap portal, GEPIA, DAVID Bioinformatics for KEGG pathway analysis, Human Protein Atlas (HPA), STRING, and TIMER. Our results indicate that the mRNA levels of RFC1 and RFC5 were the least expressed among CRC cell lines compared to other RFC complex subunits. Notably, low RFC1 and RFC5 expression was correlated with poor prognosis in terms of CRC patients' overall survival (OS). Immunohistochemical results from the Human Protein Atlas demonstrated medium staining for RFC1, RFC2, and RFC5 in CRC tissues. Furthermore, the low expression of RFC1 and RFC5 showed a significant correlation with high expression levels of miR-26a-5p and miR-636, impacting cell proliferation through mismatch repair, DNA replication, and the nucleotide excision repair pathway. Although the precise functions of RFC1 in cancer are still unknown, our findings suggest that the small-molecule single target, CHEMBL430483, and multiple target molecules could be potential treatments for CRC. In conclusion, the elevated expression of miR-26a-5p and miR-636 targeting RFC1 and RFC5 expression holds promise as a potential biomarker for early-stage CRC detection. These insights provide novel directions and strategies for CRC therapies.

Extensive dynamic changes in the human transcriptome and its circadian organization during prolonged bed rest

Physiological and molecular processes including the transcriptome change across the 24-h day, driven by molecular circadian clocks and behavioral and systemic factors. It is not known how the temporal organization of the human transcriptome responds to a long-lasting challenge. This may, however, provide insights into adaptation, disease, and recovery. We investigated the human 24-h time series transcriptome in 20 individuals during a 90-day constant bed rest protocol. We show that the protocol affected 91% of the transcriptome with 76% of the transcriptome still affected after 10 days of recovery. Dimensionality-reduction approaches revealed that many affected transcripts were associated with mRNA translation and immune function. The number, amplitude, and phase of rhythmic transcripts, including clock genes, varied significantly across the challenge. These findings of long-lasting changes in the temporal organization of the transcriptome have implications for understanding the mechanisms underlying health consequences of conditions such as microgravity and bed rest.