Using GOstats to test gene lists for GO term association

Diet-Induced Obesity Induces Transcriptomic Changes in Neuroimmunometabolic-Related Genes in the Striatum and Olfactory Bulb

The incidence of obesity has markedly increased globally over the last several decades and is believed to be associated with the easier availability of energy-dense foods, including high-fat foods. The reinforcing hedonic properties of high-fat foods, including olfactory cues, activate reward centers in the brain, motivating eating behavior. Thus, there is a growing interest in the understanding of the genetic changes that occur in the brain that are associated with obesity and eating behavior. This growing interest has paralleled advances in genomic methods that enable transcriptomic-wide analyses. Here, we examined the transcriptomic-level differences in the olfactory bulb and striatum, regions of the brain associated with olfaction and hedonic food-seeking, respectively, in high-fat-diet (HFD)-fed obese mice. To isolate the dietary effects from obesity, we also examined transcriptomic changes in normal-chow-fed and limited-HFD-fed groups, with the latter being pair-fed with an HFD isocaloric to the consumption of the normal-chow-fed mice. Using RNA sequencing, we identified 274 differentially expressed genes (DEGs) in the striatum and 11 in the olfactory bulb of ad libitum HFD-fed mice compared to the chow-fed group, and thirty-eight DEGs in the striatum between the ad libitum HFD and limited-HFD-fed groups. The DEGs in both tissues were associated with inflammation and immune-related pathways, including oxidative stress and immune function, and with mitochondrial dysfunction and reward pathways in the striatum. These results shed light on potential obesity-associated genes in these regions of the brain.

Molecular Mechanisms of Temperature Tolerance Plasticity in an Arthropod

How species thrive in a wide range of environments is a major focus of evolutionary biology. For many species, limited genetic diversity or gene flow among habitats means that phenotypic plasticity must play an important role in their capacity to tolerate environmental heterogeneity and to colonize new habitats. However, we have a limited understanding of the molecular components that govern plasticity in ecologically relevant phenotypes. We examined this hypothesis in a spider species (Stegodyphus dumicola) with extremely low species-wide genetic diversity that nevertheless occupies a broad range of thermal environments. We determined phenotypic responses to temperature stress in individuals from four climatic zones using common garden acclimation experiments to disentangle phenotypic plasticity from genetic adaptations. Simultaneously, we created data sets on multiple molecular modalities: the genome, the transcriptome, the methylome, the metabolome, and the bacterial microbiome to determine associations with phenotypic responses. Analyses of phenotypic and molecular associations reveal that acclimation responses in the transcriptome and metabolome correlate with patterns of phenotypic plasticity in temperature tolerance. Surprisingly, genes whose expression seemed to be involved in plasticity in temperature tolerance were generally highly methylated contradicting the idea that DNA methylation stabilizes gene expression. This suggests that the function of DNA methylation in invertebrates varies not only among species but also among genes. The bacterial microbiome was stable across the acclimation period; combined with our previous demonstrations that the microbiome is temporally stable in wild populations, this is convincing evidence that the microbiome does not facilitate plasticity in temperature tolerance. Our results suggest that population-specific variation in temperature tolerance among acclimation temperatures appears to result from the evolution of plasticity in mainly gene expression.

mTORC1 Signaling in Brain Endothelial Progenitors Contributes to CCM Pathogenesis

BACKGROUND

Cerebral vascular malformations (CCMs) are primarily found within the brain, where they result in increased risk for stroke, seizures, and focal neurological deficits. The unique feature of the brain vasculature is the blood-brain barrier formed by the brain neurovascular unit. Recent studies suggest that loss of CCM genes causes disruptions of blood-brain barrier integrity as the inciting events for CCM development. CCM lesions are proposed to be initially derived from a single clonal expansion of a subset of angiogenic venous capillary endothelial cells (ECs) and respective resident endothelial progenitor cells (EPCs). However, the critical signaling events in the subclass of brain ECs/EPCs for CCM lesion initiation and progression are unclear.

METHODS

Brain EC-specific CCM3-deficient (Pdcd10BECKO) mice were generated by crossing Pdcd10fl/fl mice with Mfsd2a-CreERT2 mice. Single-cell RNA-sequencing analyses were performed by the chromium single-cell platform (10× genomics). Cell clusters were annotated into EC subtypes based on visual inspection and GO analyses. Cerebral vessels were visualized by 2-photon in vivo imaging and tissue immunofluorescence analyses. Regulation of mTOR (mechanistic target of rapamycin) signaling by CCM3 and Cav1 (caveolin-1) was performed by cell biology and biochemical approaches.

RESULTS

Single-cell RNA-sequencing analyses from P10 Pdcd10BECKO mice harboring visible CCM lesions identified upregulated CCM lesion signature and mitotic EC clusters but decreased blood-brain barrier-associated EC clusters. However, a unique EPC cluster with high expression levels of stem cell markers enriched with mTOR signaling was identified from early stages of the P6 Pdcd10BECKO brain. Indeed, mTOR signaling was upregulated in both mouse and human CCM lesions. Genetic deficiency of Raptor (regulatory-associated protein of mTOR), but not of Rictor (rapamycin-insensitive companion of mTOR), prevented CCM lesion formation in the Pdcd10BECKO model. Importantly, the mTORC1 (mTOR complex 1) pharmacological inhibitor rapamycin suppressed EPC proliferation and ameliorated CCM pathogenesis in Pdcd10BECKO mice. Mechanistic studies suggested that Cav1/caveolae increased in CCM3-depleted EPC-mediated intracellular trafficking and complex formation of the mTORC1 signaling proteins.

CONCLUSIONS

CCM3 is critical for maintaining blood-brain barrier integrity and CCM3 loss-induced mTORC1 signaling in brain EPCs initiates and facilitates CCM pathogenesis.

Inorganic Biomaterials Shape the Transcriptome Profile to Induce Endochondral Differentiation

Minerals play a vital role, working synergistically with enzymes and other cofactors to regulate physiological functions including tissue healing and regeneration. The bioactive characteristics of mineral-based nanomaterials can be harnessed to facilitate in situ tissue regeneration by attracting endogenous progenitor and stem cells and subsequently directing tissue-specific differentiation. Here, cellular responses of human mesenchymal stem/stromal cells to traditional bioactive mineral-based nanomaterials, such as hydroxyapatite, whitlockite, silicon-dioxide, and the emerging synthetic 2D nanosilicates are investigated. Transcriptome sequencing is utilized to probe the cellular response and determine the significantly affected signaling pathways due to exposure to these inorganic nanomaterials. Transcriptome profiles of stem cells treated with nanosilicates reveals a stabilized skeletal progenitor state suggestive of endochondral differentiation. This observation is bolstered by enhanced deposition of matrix mineralization in nanosilicate treated stem cells compared to control or other treatments. Specifically, use of 2D nanosilicates directs osteogenic differentiation of stem cells via activation of bone morphogenetic proteins and hypoxia-inducible factor 1-alpha signaling pathway. This study provides  insight into impact of nanomaterials on cellular gene expression profile and predicts downstream effects of nanomaterial induction of endochondral differentiation.

Non-Invasive Quality Control of Organoid Cultures Using Mesofluidic CSTR Bioreactors and High-Content Imaging

Human brain organoids produce anatomically relevant cellular structures and recapitulate key aspects of in vivo brain function, which holds great potential to model neurological diseases and screen therapeutics. However, the long growth time of 3D systems complicates the culturing of brain organoids and results in heterogeneity across samples hampering their applications. We developed an integrated platform to enable robust and long-term culturing of 3D brain organoids. We designed a mesofluidic bioreactor device based on a reaction-diffusion scaling theory, which achieves robust media exchange for sufficient nutrient delivery in long-term culture. We integrated this device with longitudinal tracking and machine learning-based classification tools to enable non-invasive quality control of live organoids. This integrated platform allows for sample pre-selection for downstream molecular analysis. Transcriptome analyses of organoids revealed that our mesofluidic bioreactor promoted organoid development while reducing cell death. Our platform thus offers a generalizable tool to establish reproducible culture standards for 3D cellular systems for a variety of applications beyond brain organoids.

Inorganic ions activate lineage-specific gene regulatory networks

Inorganic biomaterials have been shown to direct cellular responses, including cell-cell and cell-matrix interactions. Notably, ions released from these inorganic biomaterials play a vital role in defining cell identity, and promoting tissue-specific functions. However, the effect of inorganic ions on cellular functions have yet to be investigated at the transcriptomic level, representing a critical knowledge gap in the development of next-generation bioactive materials. To address this gap, we investigated the impact of various inorganic ions including silver, copper, titanium, and platinum on human mesenchymal stem cells (hMSCs). Our finding showed that silver and copper induce osteogenic and chondrogenic differentiation respectively, through enrichment of lineage-specific gene expression program. In particular, silver effectively induced Wingless/Integrated (Wnt) and mitogen-activated protein kinase (MAPK) signaling, which are vital for osteogenesis. On the other hand, copper specifically stimulated Transforming growth factor beta (TGFβ) signaling, while suppressing Janus kinase/signal transducers and activators of transcription (JAK-STAT) signaling, thereby promoting chondrogenesis. In contrast, platinum, and tantalum, ions didn't stimulate regenerative responses. Together, our findings highlight the potential of inorganic biomaterials in tissue regeneration strategies, which currently rely largely on growth factors and small molecule therapeutics. STATEMENT OF SIGNIFICANCE: This research emphasizes the critical role of bioactive inorganic ions in controlling lineage-specific gene expression patterns in mesenchymal stem cells, effectively modulating the transcriptome landscape and directing cell fate. The study lays the foundation for a systematic database of biomaterial candidates and their effects on cellular functions, which will ultimately streamline the translation of new biomaterials into clinical applications.

Spodoptera frugiperda Salivary Glucose Oxidase Reduces the Release of Green Leaf Volatiles and Increases Terpene Emission from Maize

The intricate relationships between plants and insects are essential for understanding ecological dynamics. Among these interactions, HIPVs serve as a pivotal defense mechanism. Our findings reveal the highly conserved nature of the GOX gene within the Lepidoptera order, highly expressed in the salivary glands of S. frugiperda, and its role in mediating maize's defense responses. Notably, salivary GOX activity expression significantly decreases subsequent gene knockout. The presence of GOX in the saliva of S. frugiperda significantly modulates the emission of HIPVs during maize consumption. This research delineates that GOX selectively inhibits the emission of certain green leaf volatiles (GLVs) while concurrently enhancing the release of terpene volatiles. This study unveils a novel mechanism whereby S. frugiperda utilizes GOX proteins in OS to modulate volatile emissions from maize, offering fresh perspectives on the adaptive evolution of phytophagous insects and their interactions with their preferred host plants.

SGCP: a spectral self-learning method for clustering genes in co-expression networks

BACKGROUND

A widely used approach for extracting information from gene expression data employs the construction of a gene co-expression network and the subsequent computational detection of gene clusters, called modules. WGCNA and related methods are the de facto standard for module detection. The purpose of this work is to investigate the applicability of more sophisticated algorithms toward the design of an alternative method with enhanced potential for extracting biologically meaningful modules.

RESULTS

We present self-learning gene clustering pipeline (SGCP), a spectral method for detecting modules in gene co-expression networks. SGCP incorporates multiple features that differentiate it from previous work, including a novel step that leverages gene ontology (GO) information in a self-leaning step. Compared with widely used existing frameworks on 12 real gene expression datasets, we show that SGCP yields modules with higher GO enrichment. Moreover, SGCP assigns highest statistical importance to GO terms that are mostly different from those reported by the baselines.

CONCLUSION

Existing frameworks for discovering clusters of genes in gene co-expression networks are based on relatively simple algorithmic components. SGCP relies on newer algorithmic techniques that enable the computation of highly enriched modules with distinctive characteristics, thus contributing a novel alternative tool for gene co-expression analysis.

Facilitated introgression from domestic goat into Alpine ibex at immune loci

Hybridization can result in the transfer of adaptive genetic material from one species to another, known as adaptive introgression. Bottlenecked (and hence genetically depleted) species are expected to be particularly receptive to adaptive introgression, since introgression can introduce new or previously lost adaptive genetic variation. The Alpine ibex (Capra ibex), which recently recovered from near extinction, is known to hybridize with the domestic goat (Capra aegagrus hircus), and signals of introgression previously found at the major histocompatibility complex were suggested to potentially be adaptive. Here, we combine two ancient whole genomes of Alpine ibex with 29 modern Alpine ibex genomes and 31 genomes representing six related Capra species to investigate the genome-wide patterns of introgression and confirm the potential relevance of immune loci. We identified low rates of admixture in modern Alpine ibex through various F statistics and screening for putative introgressed tracts. Further results based on demographic modelling were consistent with introgression to have occurred during the last 300 years, coinciding with the known species bottleneck, and that in each generation, 1-2 out of 100 Alpine ibex had a domestic goat parent. The putatively introgressed haplotypes were enriched at immune-related genes, where the adaptive value of alternative alleles may give individuals with otherwise depleted genetic diversity a selective advantage. While interbreeding with domestic species is a prevalent issue in species conservation, in this specific case, it resulted in putative adaptive introgression. Our findings highlight the complex interplay between hybridization, adaptive evolution, and the potential risks and benefits associated with anthropogenic influences on wild species.

Mechano-inhibition of endocytosis sensitizes cancer cells to Fas-induced Apoptosis

The transmembrane death receptor Fas transduces apoptotic signals upon binding its ligand, FasL. Although Fas is highly expressed in cancer cells, insufficient cell surface Fas expression desensitizes cancer cells to Fas-induced apoptosis. Here, we show that the increase in Fas microaggregate formation on the plasma membrane in response to the inhibition of endocytosis sensitizes cancer cells to Fas-induced apoptosis. We used a clinically accessible Rho-kinase inhibitor, fasudil, that reduces endocytosis dynamics by increasing plasma membrane tension. In combination with exogenous soluble FasL (sFasL), fasudil promoted cancer cell apoptosis, but this collaborative effect was substantially weaker in nonmalignant cells. The combination of sFasL and fasudil prevented glioblastoma cell growth in embryonic stem cell-derived brain organoids and induced tumor regression in a xenograft mouse model. Our results demonstrate that sFasL has strong potential for apoptosis-directed cancer therapy when Fas microaggregate formation is augmented by mechano-inhibition of endocytosis.

Sexual dimorphism in the tardigrade Paramacrobiotus metropolitanus transcriptome

BACKGROUND

In gonochoristic animals, the sex determination pathway induces different morphological and behavioral features that can be observed between sexes, a condition known as sexual dimorphism. While many components of this sex differentiation cascade show high levels of diversity, factors such as the Doublesex-Mab-3-Related Transcription factor (DMRT) are widely conserved across animal taxa. Species of the phylum Tardigrada exhibit remarkable diversity in morphology and behavior between sexes, suggesting a pathway regulating this dimorphism. Despite the wealth of genomic and zoological knowledge accumulated in recent studies, the sexual differences in tardigrades genomes have not been identified. In the present study, we focused on the gonochoristic species Paramacrobiotus metropolitanus and employed omics analyses to unravel the molecular basis of sexual dimorphism.

RESULTS

Transcriptome analysis between sex-identified specimens revealed numerous differentially expressed genes, of which approximately 2,000 male-biased genes were focused on 29 non-male-specific genomic loci. From these regions, we identified two Macrobiotidae family specific DMRT paralogs, which were significantly upregulated in males and lacked sex specific splicing variants. Furthermore, phylogenetic analysis indicated all tardigrade genomes lack the doublesex ortholog, suggesting doublesex emerged after the divergence of Tardigrada. In contrast to sex-specific expression, no evidence of genomic differences between the sexes was found. We also identified several anhydrobiosis genes that exhibit sex-biased expression, suggesting a possible mechanism for protection of sex-specific tissues against extreme stress.

CONCLUSIONS

This study provides a comprehensive analysis for analyzing the genetic differences between sexes in tardigrades. The existence of male-biased, but not male-specific, genomic loci and identification of the family specific male-biased DMRT subfamily provides the foundation for understanding the sex determination cascade. In addition, sex-biased expression of several tardigrade-specific genes which are involved their stress tolerance suggests a potential role in protecting sex-specific tissue and gametes.

Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.

Epigenetic differences between wild and cultivated grapevines highlight the contribution of DNA methylation during crop domestication

The domestication process in grapevines has facilitated the fixation of desired traits. Nowadays, vegetative propagation through cuttings enables easier preservation of these genotypes compared to sexual reproduction. Nonetheless, even with vegetative propagation, various phenotypes are often present within the same vineyard due to the accumulation of somatic mutations. These mutations are not the sole factors influencing phenotype. Alongside somatic variations, epigenetic variation has been proposed as a pivotal player in regulating phenotypic variability acquired during domestication. The emergence of these epialleles might have significantly influenced grapevine domestication over time. This study aims to investigate the impact of domestication on methylation patterns in cultivated grapevines. Reduced-representation bisulfite sequencing was conducted on 18 cultivated and wild accessions. Results revealed that cultivated grapevines exhibited higher methylation levels than their wild counterparts. Differential Methylation Analysis between wild and cultivated grapevines identified a total of 9955 differentially methylated cytosines, of which 78% were hypermethylated in cultivated grapevines. Functional analysis shows that core methylated genes (consistently methylated in both wild and cultivated accessions) are associated with stress response and terpenoid/isoprenoid metabolic processes. Meanwhile, genes with differential methylation are linked to protein targeting to the peroxisome, ethylene regulation, histone modifications, and defense response. Collectively, our results highlight the significant roles that epialleles may have played throughout the domestication history of grapevines.

Lack of paternal silencing and ecotype-specific expression in head and body lice hybrids

Paternal genome elimination (PGE) is a non-Mendelian inheritance system, described in numerous arthropod species, in which males develop from fertilized eggs, but their paternally inherited chromosomes are eliminated before or during spermatogenesis. Therefore, PGE males only transmit their maternally inherited set of chromosomes to their offspring. In addition to the elimination of paternal chromosomes, diverse PGE species have also repeatedly evolved the transcriptional silencing of the paternal genome, making males effectively haploid. However, it is unclear if this paternal chromosome silencing is mechanistically linked to the chromosome elimination or has evolved at a later stage, and if so, what drives the haploidization of males under PGE. In order to understand these questions, here we study the human louse, Pediculus humanus, which represents an ideal model system, as it appears to be the only instance of PGE where males eliminate, but not silence their paternal chromosomes, although the latter remains to be shown conclusively. In this study, we analyzed parent-of-origin allele-specific expression patterns in male offspring of crosses between head and body lice ecotypes. We show that hybrid adult males of P. humanus display biparental gene expression, which constitutes the first case of a species with PGE in which genetic activity of paternal chromosomes in the soma is not affected by embryonic silencing or (partial or complete) elimination. We did however also identify a small number of maternally biased genes (potentially imprinted genes), which may be involved in the elimination of paternal chromosomes during spermatogenesis. Finally, we have identified genes that show ecotype-specific expression bias. Given the low genetic diversity between ecotypes, this is suggestive for a role of epigenetic processes in ecotype differences.

Testing the Significance of Ranked Gene Sets in Genome-wide Transcriptome

Background

Popular gene set enrichment analysis approaches assumed that genes in the gene set contributed to the statistics equally. However, the genes in the transcription factors (TFs) derived gene sets, or gene sets constructed by TF targets identified by the ChIP-Seq experiment, have a rank attribute, as each of these genes have been assigned with a p-value which indicates the true or false possibilities of the ownerships of the genes belong to the gene sets.

Objectives

Ignoring the rank information during the enrichment analysis will lead to improper statistical inference. We address this issue by developing of new method to test the significance of ranked gene sets in genome-wide transcriptome profiling data.

Methods

A method was proposed by first creating ranked gene sets and gene lists and then applying weighted Kendall's tau rank correlation statistics to the test. After introducing top-down weights to the genes in the gene set, a new software called "Flaver" was developed.

Results

Theoretical properties of the proposed method were established, and its differences over the GSEA approach were demonstrated when analyzing the transcriptome profiling data across 55 human tissues and 176 human cell-lines. The results indicated that the TFs identified by our method have higher tendency to be differentially expressed across the tissues analyzed than its competitors. It significantly outperforms the well-known gene set enrichment analyzing tools, GOStats (9%) and GSEA (17%), in analyzing well-documented human RNA transcriptome datasets.

Conclusions

The method is outstanding in detecting gene sets of which the gene ranks were correlated with the expression levels of the genes in the transcriptome data.

The interplay between electron transport chain function and iron regulatory factors influences melanin formation in Cryptococcus neoformans

Mitochondrial functions are critical for the ability of the fungal pathogen Cryptococcus neoformans to cause disease. However, mechanistic connections between key functions such as the mitochondrial electron transport chain (ETC) and virulence factor elaboration have yet to be thoroughly characterized. Here, we observed that inhibition of ETC complex III suppressed melanin formation, a major virulence factor. This inhibition was partially overcome by defects in Cir1 or HapX, two transcription factors that regulate iron acquisition and use. In this regard, loss of Cir1 derepresses the expression of laccase genes as a potential mechanism to restore melanin, while HapX may condition melanin formation by controlling oxidative stress. We hypothesize that ETC dysfunction alters redox homeostasis to influence melanin formation. Consistent with this idea, inhibition of growth by hydrogen peroxide was exacerbated in the presence of the melanin substrate L-DOPA. In addition, loss of the mitochondrial chaperone Mrj1, which influences the activity of ETC complex III and reduces ROS accumulation, also partially overcame antimycin A inhibition of melanin. The phenotypic impact of mitochondrial dysfunction was consistent with RNA-Seq analyses of WT cells treated with antimycin A or L-DOPA, or cells lacking Cir1 that revealed influences on transcripts encoding mitochondrial functions (e.g., ETC components and proteins for Fe-S cluster assembly). Overall, these findings reveal mitochondria-nuclear communication via ROS and iron regulators to control virulence factor production in C. neoformans.IMPORTANCEThere is a growing appreciation of the importance of mitochondrial functions and iron homeostasis in the ability of fungal pathogens to sense the vertebrate host environment and cause disease. Many mitochondrial functions such as heme and iron-sulfur cluster biosynthesis, and the electron transport chain (ETC), are dependent on iron. Connections between factors that regulate iron homeostasis and mitochondrial activities are known in model yeasts and are emerging for fungal pathogens. In this study, we identified connections between iron regulatory transcription factors (e.g., Cir1 and HapX) and the activity of complex III of the ETC that influence the formation of melanin, a key virulence factor in the pathogenic fungus Cryptococcus neoformans. This fungus causes meningoencephalitis in immunocompromised people and is a major threat to the HIV/AIDS population. Thus, understanding how mitochondrial functions influence virulence may support new therapeutic approaches to combat diseases caused by C. neoformans and other fungi.

Identification and functional characterization of conserved cis-regulatory elements responsible for early fruit development in cucurbit crops

A number of cis-regulatory elements (CREs) conserved during evolution have been found to be responsible for phenotypic novelty and variation. Cucurbit crops such as cucumber (Cucumis sativus), watermelon (Citrullus lanatus), melon (Cucumis melo), and squash (Cucurbita maxima) develop fruits from an inferior ovary and share some similar biological processes during fruit development. Whether conserved regulatory sequences play critical roles in fruit development of cucurbit crops remains to be explored. In six well-studied cucurbit species, we identified 392,438 conserved noncoding sequences (CNSs), including 82,756 that are specific to cucurbits, by comparative genomics. Genome-wide profiling of accessible chromatin regions (ACRs) and gene expression patterns mapped 20,865 to 43,204 ACRs and their potential target genes for two fruit tissues at two key developmental stages in six cucurbits. Integrated analysis of CNSs and ACRs revealed 4,431 syntenic orthologous CNSs, including 1,687 cucurbit-specific CNSs that overlap with ACRs that are present in all six cucurbit crops and that may regulate the expression of 757 adjacent orthologous genes. CRISPR mutations targeting two CNSs present in the 1,687 cucurbit-specific sequences resulted in substantially altered fruit shape and gene expression patterns of adjacent NAC1 (NAM, ATAF1/2, and CUC2) and EXT-like (EXTENSIN-like) genes, validating the regulatory roles of these CNSs in fruit development. These results not only provide a number of target CREs for cucurbit crop improvement, but also provide insight into the roles of CREs in plant biology and during evolution.

Netrin G1 Ligand is a new stromal immunomodulator that promotes pancreatic cancer

Understanding pancreatic cancer biology is fundamental for identifying new targets and for developing more effective therapies. In particular, the contribution of the stromal microenvironment to pancreatic cancer tumorigenesis requires further exploration. Here, we report the stromal roles of the synaptic protein Netrin G1 Ligand (NGL-1) in pancreatic cancer, uncovering its pro-tumor functions in cancer-associated fibroblasts and in immune cells. We observed that the stromal expression of NGL-1 inversely correlated with patients' overall survival. Moreover, germline knockout (KO) mice for NGL-1 presented decreased tumor burden, with a microenvironment that is less supportive of tumor growth. Of note, tumors from NGL-1 KO mice produced less immunosuppressive cytokines and displayed an increased percentage of CD8 + T cells than those from control mice, while preserving the physical structure of the tumor microenvironment. These effects were shown to be mediated by NGL-1 in both immune cells and in the local stroma, in a TGF-β-dependent manner. While myeloid cells lacking NGL-1 decreased the production of immunosuppressive cytokines, NGL-1 KO T cells showed increased proliferation rates and overall polyfunctionality compared to control T cells. CAFs lacking NGL-1 were less immunosuppressive than controls, with overall decreased production of pro-tumor cytokines and compromised ability to inhibit CD8 + T cells activation. Mechanistically, these CAFs downregulated components of the TGF-β pathway, AP-1 and NFAT transcription factor families, resulting in a less tumor-supportive phenotype. Finally, targeting NGL-1 genetically or using a functionally antagonistic small peptide phenocopied the effects of chemotherapy, while modulating the immunosuppressive tumor microenvironment (TME), rather than eliminating it. We propose NGL-1 as a new local stroma and immunomodulatory molecule, with pro-tumor roles in pancreatic cancer.

Statement of Significance

Here we uncovered the pro-tumor roles of the synaptic protein NGL-1 in the tumor microenvironment of pancreatic cancer, defining a new target that simultaneously modulates tumor cell, fibroblast, and immune cell functions. This study reports a new pathway where NGL-1 controls TGF-β, AP-1 transcription factor members and NFAT1, modulating the immunosuppressive microenvironment in pancreatic cancer. Our findings highlight NGL-1 as a new stromal immunomodulator in pancreatic cancer.

BugSigDB captures patterns of differential abundance across a broad range of host-associated microbial signatures

The literature of human and other host-associated microbiome studies is expanding rapidly, but systematic comparisons among published results of host-associated microbiome signatures of differential abundance remain difficult. We present BugSigDB, a community-editable database of manually curated microbial signatures from published differential abundance studies accompanied by information on study geography, health outcomes, host body site and experimental, epidemiological and statistical methods using controlled vocabulary. The initial release of the database contains >2,500 manually curated signatures from >600 published studies on three host species, enabling high-throughput analysis of signature similarity, taxon enrichment, co-occurrence and coexclusion and consensus signatures. These data allow assessment of microbiome differential abundance within and across experimental conditions, environments or body sites. Database-wide analysis reveals experimental conditions with the highest level of consistency in signatures reported by independent studies and identifies commonalities among disease-associated signatures, including frequent introgression of oral pathobionts into the gut.

Targeting Colorectal Cancer: Unravelling the Transcriptomic Impact of Cisplatin and High-THC Cannabis Extract

Cisplatin and other platinum-derived chemotherapy drugs have been used for the treatment of cancer for a long time and are often combined with other medications. Unfortunately, tumours often develop resistance to cisplatin, forcing scientists to look for alternatives or synergistic combinations with other drugs. In this work, we attempted to find a potential synergistic effect between cisplatin and cannabinoid delta-9-THC, as well as the high-THC Cannabis sativa extract, for the treatment of HT-29, HCT-116, and LS-174T colorectal cancer cell lines. However, we found that combinations of the high-THC cannabis extract with cisplatin worked antagonistically on the tested colorectal cancer cell lines. To elucidate the mechanisms of drug interactions and the distinct impacts of individual treatments, we conducted a comprehensive transcriptomic analysis of affected pathways within the colorectal cancer cell line HT-29. Our primary objective was to gain a deeper understanding of the underlying molecular mechanisms associated with each treatment modality and their potential interactions. Our findings revealed an antagonistic interaction between cisplatin and high-THC cannabis extract, which could be linked to alterations in gene transcription associated with cell death (BCL2, BAD, caspase 10), DNA repair pathways (Rad52), and cancer pathways related to drug resistance.

Meta-analysis of differential gene expression in lower motor neurons isolated by laser capture microdissection from post-mortem ALS spinal cords

Introduction

ALS is a fatal neurodegenerative disease for which underlying mechanisms are incompletely understood. The motor neuron is a central player in ALS pathogenesis but different transcriptome signatures have been derived from bulk analysis of post-mortem tissue and iPSC-derived motor neurons (iPSC-MNs).

Methods

This study performed a meta-analysis of six gene expression studies (microarray and RNA-seq) in which laser capture microdissection (LCM) was used to isolate lower motor neurons from post-mortem spinal cords of ALS and control (CTL) subjects. Differentially expressed genes (DEGs) with consistent ALS versus CTL expression differences across studies were identified.

Results

The analysis identified 222 ALS-increased DEGs (FDR <0.10, SMD >0.80) and 278 ALS-decreased DEGs (FDR <0.10, SMD < -0.80). ALS-increased DEGs were linked to PI3K-AKT signaling, innate immunity, inflammation, motor neuron differentiation and extracellular matrix. ALS-decreased DEGs were associated with the ubiquitin-proteosome system, microtubules, axon growth, RNA-binding proteins and synaptic membrane. ALS-decreased DEG mRNAs frequently interacted with RNA-binding proteins (e.g., FUS, HuR). The complete set of DEGs (increased and decreased) overlapped significantly with genes near ALS-associated SNP loci (p < 0.01). Transcription factor target motifs with increased proximity to ALS-increased DEGs were identified, most notably DNA elements predicted to interact with forkhead transcription factors (e.g., FOXP1) and motor neuron and pancreas homeobox 1 (MNX1). Some of these DNA elements overlie ALS-associated SNPs within known enhancers and are predicted to have genotype-dependent MNX1 interactions. DEGs were compared to those identified from SOD1-G93A mice and bulk spinal cord segments or iPSC-MNs from ALS patients. There was good correspondence with transcriptome changes from SOD1-G93A mice (r ≤ 0.408) but most DEGs were not differentially expressed in bulk spinal cords or iPSC-MNs and transcriptome-wide effect size correlations were weak (bulk tissue: r ≤ 0.207, iPSC-MN: r ≤ 0.037).

Conclusion

This study defines a robust transcriptome signature from LCM-based motor neuron studies of post-mortem tissue from ALS and CTL subjects. This signature differs from those obtained from analysis of bulk spinal cord segments and iPSC-MNs. Results provide insight into mechanisms underlying gene dysregulation in ALS and highlight connections between these mechanisms, ALS genetics, and motor neuron biology.

Intestinal cell diversity and treatment responses in a parasitic nematode at single cell resolution

BACKGROUND

Parasitic nematodes, significant pathogens for humans, animals, and plants, depend on diverse organ systems for intra-host survival. Understanding the cellular diversity and molecular variations underlying these functions holds promise for developing novel therapeutics, with specific emphasis on the neuromuscular system's functional diversity. The nematode intestine, crucial for anthelmintic therapies, exhibits diverse cellular phenotypes, and unraveling this diversity at the single-cell level is essential for advancing knowledge in anthelmintic research across various organ systems.

RESULTS

Here, using novel single-cell transcriptomics datasets, we delineate cellular diversity within the intestine of adult female Ascaris suum, a parasitic nematode species that infects animals and people. Gene transcripts expressed in individual nuclei of untreated intestinal cells resolved three phenotypic clusters, while lower stringency resolved additional subclusters and more potential diversity. Clusters 1 and 3 phenotypes displayed variable congruence with scRNA phenotypes of C. elegans intestinal cells, whereas the A. suum cluster 2 phenotype was markedly unique. Distinct functional pathway enrichment characterized each A. suum intestinal cell cluster. Cluster 2 was distinctly enriched for Clade III-associated genes, suggesting it evolved within clade III nematodes. Clusters also demonstrated differential transcriptional responsiveness to nematode intestinal toxic treatments, with Cluster 2 displaying the least responses to short-term intra-pseudocoelomic nematode intestinal toxin treatments.

CONCLUSIONS

This investigation presents advances in knowledge related to biological differences among major cell populations of adult A. suum intestinal cells. For the first time, diverse nematode intestinal cell populations were characterized, and associated biological markers of these cells were identified to support tracking of constituent cells under experimental conditions. These advances will promote better understanding of this and other parasitic nematodes of global importance, and will help to guide future anthelmintic treatments.

Cockayne Syndrome Patient iPSC-Derived Brain Organoids and Neurospheres Show Early Transcriptional Dysregulation of Biological Processes Associated with Brain Development and Metabolism

Cockayne syndrome (CS) is a rare hereditary autosomal recessive disorder primarily caused by mutations in Cockayne syndrome protein A (CSA) or B (CSB). While many of the functions of CSB have been at least partially elucidated, little is known about the actual developmental dysregulation in this devasting disorder. Of particular interest is the regulation of cerebral development as the most debilitating symptoms are of neurological nature. We generated neurospheres and cerebral organoids utilizing Cockayne syndrome B protein (CSB)-deficient induced pluripotent stem cells derived from two patients with distinct severity levels of CS and healthy controls. The transcriptome of both developmental timepoints was explored using RNA-Seq and bioinformatic analysis to identify dysregulated biological processes common to both patients with CS in comparison to the control. CSB-deficient neurospheres displayed upregulation of the VEGFA-VEGFR2 signalling pathway, vesicle-mediated transport and head development. CSB-deficient cerebral organoids exhibited downregulation of brain development, neuron projection development and synaptic signalling. We further identified the upregulation of steroid biosynthesis as common to both timepoints, in particular the upregulation of the cholesterol biosynthesis branch. Our results provide insights into the neurodevelopmental dysregulation in patients with CS and strengthen the theory that CS is not only a neurodegenerative but also a neurodevelopmental disorder.

HIHISIV: a database of gene expression in HIV and SIV host immune response

In the battle of the host against lentiviral pathogenesis, the immune response is crucial. However, several questions remain unanswered about the interaction with different viruses and their influence on disease progression. The simian immunodeficiency virus (SIV) infecting nonhuman primates (NHP) is widely used as a model for the study of the human immunodeficiency virus (HIV) both because they are evolutionarily linked and because they share physiological and anatomical similarities that are largely explored to understand the disease progression. The HIHISIV database was developed to support researchers to integrate and evaluate the large number of transcriptional data associated with the presence/absence of the pathogen (SIV or HIV) and the host response (NHP and human). The datasets are composed of microarray and RNA-Seq gene expression data that were selected, curated, analyzed, enriched, and stored in a relational database. Six query templates comprise the main data analysis functions and the resulting information can be downloaded. The HIHISIV database, available at  https://hihisiv.github.io , provides accurate resources for browsing and visualizing results and for more robust analyses of pre-existing data in transcriptome repositories.

A Critical Role for IFN-β Signaling for IFN-κ Induction in Keratinocytes

Background/Purpose

Cutaneous lupus erythematosus (CLE) affects up to 70% of patients with systemic lupus erythematosus (SLE), and type I interferons (IFNs) are important promoters of SLE and CLE. Our previous work identified IFN-kappa (IFN-κ), a keratinocyte-produced type I IFN, as upregulated in non-lesional and lesional lupus skin and as a critical regulator for enhanced UVB-mediated cell death in SLE keratinocytes. Importantly, the molecular mechanisms governing regulation of IFN-κ expression have been relatively unexplored. Thus, this study sought to identify critical regulators of IFN-κ and identified a novel role for IFN-beta (IFN-β).

Methods

Human N/TERT keratinocytes were treated with the RNA mimic poly (I:C) or 50 mJ/cm2 ultraviolet B (UVB), followed by mRNA expression quantification by RT-qPCR in the presence or absence neutralizing antibody to the type I IFN receptor (IFNAR). IFNB and STAT1 knockout (KO) keratinocytes were generated using CRISPR/Cas9.

Results

Time courses of poly(I:C) and UVB treatment revealed a differential expression of IFNB, which was upregulated between 3-6 hours and IFNK, which was upregulated 24 hours after stimulation. Intriguingly, only IFNK expression was substantially abrogated by neutralizing antibodies to IFNAR, suggesting that IFNK upregulation required type I IFN signaling for induction. Indeed, deletion of IFNB abrogated IFNK expression. Further exploration confirmed a role for type I IFN-triggered STAT1 activation.

Conclusion

Collectively, our work describes a novel mechanistic paradigm in keratinocytes in which initial IFN-κ induction in response to poly(I:C) and UVB is IFNβ1-dependent, thus describing IFNK as both an IFN gene and an interferon-stimulated gene.

Single-cell RNA-seq revealed heterogeneous responses and functional differentiation of hemocytes against white spot syndrome virus infection in Litopenaeus vannamei

Shrimp hemocytes are the vital immune cells participating in innate immune response to defend against viruses. However, the lack of specific molecular markers for shrimp hemocyte hindered the insightful understanding of their functional clusters and differential roles in combating microbial infections. In this study, we used single-cell RNA sequencing to map the transcriptomic landscape of hemocytes from the white spot syndrome virus (WSSV)-infected Litopenaeus vannamei and conjointly analyzed with our previous published single-cell RNA sequencing technology data from the healthy hemocytes. A total of 16 transcriptionally distinct cell clusters were identified, which occupied different proportions in healthy and WSSV-infected hemocytes and exerted differential roles in antiviral immune response. Following mapping of the sequencing data to the WSSV genome, we found that all types of hemocytes could be invaded by WSSV virions, especially the cluster 8, which showed the highest transcriptional levels of WSSV genes and exhibited a cell type-specific antiviral response to the viral infection. Further evaluation of the cell clusters revealed the delicate dynamic balance between hemocyte immune response and viral infestation. Unsupervised pseudo-time analysis of hemocytes showed that the hemocytes in immune-resting state could be significantly activated upon WSSV infection and then functionally differentiated to different hemocyte subsets. Collectively, our results revealed the differential responses of shrimp hemocytes and the process of immune-functional differentiation post-WSSV infection, providing essential resource for the systematic insight into the synergistic immune response mechanism against viral infection among hemocyte subtypes.

IMPORTANCE

Current knowledge of shrimp hemocyte classification mainly comes from morphology, which hinder in-depth characterization of cell lineage development, functional differentiation, and different immune response of hemocyte types during pathogenic infections. Here, single-cell RNA sequencing was used for mapping hemocytes during white spot syndrome virus (WSSV) infection in Litopenaeus vannamei, identifying 16 cell clusters and evaluating their potential antiviral functional characteristics. We have described the dynamic balance between viral infestation and hemocyte immunity. And the functional differentiation of hemocytes under WSSV stimulation was further characterized. Our results provided a comprehensive transcriptional landscape and revealed the heterogeneous immune response in shrimp hemocytes during WSSV infection.

Unraveling athletic performance: Transcriptomics and external load monitoring in handball competition

OBJECTIVE

This study aims to comprehend the impact of handball practice on sub-elite athletes by investigating transcriptomic changes that occur during a match. The primary focus encompasses a dual objective: firstly, to identify and characterize these transcriptomic alterations, and secondly, to establish correlations between internal factors (gene expression), and external loads measured through Electronic Performance and Tracking Systems (EPTS variables). Ultimately, this comprehensive analysis seeks to evaluate both acute and chronic responses to exercise within the context of handball training.

METHODS

The study included sixteen elite male athletes from the FC Barcelona handball second team. Blood samples were extracted at three different time points: before the match at baseline levels (T1), immediately upon completion (T2), and 24 hours after completion (T3). Differential gene expression, Gene Ontology Term and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were conducted in two comparisons: Comparison 1 (T1 vs T2) and Comparison 2 (T1 vs T3). Further, the correlation between gene expression levels and training variables (external load) was conducted.

RESULTS

In T1 vs T2, 3717 of the 14632 genes detected were differentially expressed (adjusted p-value < 0.05), and enrichment of terms related to the immune system, mitochondria, and metabolic processes was found. Further, significant linear correlations were obtained between High-Speed running (HSR) and high-intensity variables such as acceleration ACC and deceleration DEC values with amino acids, and inflammatory and oxidative environment-related pathways, both in chronic and acute response.

CONCLUSIONS

This research highlights the effects of external workload on elite athletes during a handball match and throughout the season. The study identifies deregulation in the immune system, mitochondrial functions, and various metabolic pathways during the match. Additionally, it establishes correlations between the external load and pathways associated with amino acids, inflammation, oxidative environment, and regulation. These findings offer insights into the immediate and chronic responses of athletes to physical effort.

Decoding the leaf apical meristem of Guarea glabra Vahl (Meliaceae): insight into the evolution of indeterminate pinnate leaves

In seed plants, growth of shoots and roots is indeterminate, while leaves are typically determinate organs that cease to grow after a certain developmental stage. This is due to the characteristics of the leaf meristem, where cell proliferation activity is retained only for a limited period. However, several plants exhibit indeterminacy in their leaves, exemplified by the pinnate compound leaves of Guarea and Chisocheton genera in the Meliaceae family. In these plants, the leaf meristem at the tip of the leaf retains meristematic activity and produces leaflets over years, resulting in a single leaf that resembles a twig. The molecular mechanism underlying the indeterminate leaf meristem of these plants has not been examined. In this research, we used Guarea glabra as a model to investigate the development of indeterminate pinnate leaves. Transcriptome analyses revealed that the gene expression profile in leaf apex tissue differed from that in the shoot apex. However, a class 1 KNOTTED-LIKE HOMEOBOX (KNOX1) gene which is lost in Brassicaceae was highly expressed in both tissues. We established an in situ hybridisation system for this species using Technovit 9100 to analyse the spatial expression patterns of genes. We revealed that the leaf meristematic region of G. glabra expresses KNOX1, LEAFY and ANGUSTIFORIA3 simultaneously, suggesting the involvement of these genes in the indeterminacy of the leaf meristem.

A single-nuclei paired multiomic analysis of the human midbrain reveals age- and Parkinson's disease-associated glial changes

Age is the primary risk factor for Parkinson's disease (PD), but how aging changes the expression and regulatory landscape of the brain remains unclear. Here we present a single-nuclei multiomic study profiling shared gene expression and chromatin accessibility of young, aged and PD postmortem midbrain samples. Combined multiomic analysis along a pseudopathogenesis trajectory reveals that all glial cell types are affected by age, but microglia and oligodendrocytes are further altered in PD. We present evidence for a disease-associated oligodendrocyte subtype and identify genes lost over the aging and disease process, including CARNS1, that may predispose healthy cells to develop a disease-associated phenotype. Surprisingly, we found that chromatin accessibility changed little over aging or PD within the same cell types. Peak-gene association patterns, however, are substantially altered during aging and PD, identifying cell-type-specific chromosomal loci that contain PD-associated single-nucleotide polymorphisms. Our study suggests a previously undescribed role for oligodendrocytes in aging and PD.

The mutation of NSUN5 R295C promotes preeclampsia by impairing decidualization through downregulating IL-11Rα

Preeclampsia (PE) is a pregnancy-specific hypertensive disorder that severely impairs maternal and fetal health. However, its pathogenesis remains elusive. NOP2/Sun5 (NSUN5) is an RNA methyltransferase. This study discovered a significant correlation between rs77133388 of NSUN5 and PE in a cohort of 868 severe PE patients and 982 healthy controls. To further explore this association, the researchers generated single-base mutant mice (NSUN5 R295C) at rs77133388. The pregnant NSUN5 R295C mice exhibited PE symptoms. Additionally, compared to the controls, the decidual area of the placenta was significantly reduced in NSUN5 R295C mice, and their decidualization was impaired with a significantly decrease in polyploid cell numbers after artificially induced decidualization. The study also found a decrease in phosphorylated JAK2, STAT3, and IL-11Rα, Cyclin D3 expression in NSUN5 R295C mice. Overall, these findings suggest that NSUN5 mutation potentially alters decidualization through the IL-11Rα/JAK2/STAT3/Cyclin D3 pathway, ultimately impairing placental development and contributing to PE occurrence.

Drug-target identification in COVID-19 disease mechanisms using computational systems biology approaches

Introduction

The COVID-19 Disease Map project is a large-scale community effort uniting 277 scientists from 130 Institutions around the globe. We use high-quality, mechanistic content describing SARS-CoV-2-host interactions and develop interoperable bioinformatic pipelines for novel target identification and drug repurposing.

Methods

Extensive community work allowed an impressive step forward in building interfaces between Systems Biology tools and platforms. Our framework can link biomolecules from omics data analysis and computational modelling to dysregulated pathways in a cell-, tissue- or patient-specific manner. Drug repurposing using text mining and AI-assisted analysis identified potential drugs, chemicals and microRNAs that could target the identified key factors.

Results

Results revealed drugs already tested for anti-COVID-19 efficacy, providing a mechanistic context for their mode of action, and drugs already in clinical trials for treating other diseases, never tested against COVID-19.

Discussion

The key advance is that the proposed framework is versatile and expandable, offering a significant upgrade in the arsenal for virus-host interactions and other complex pathologies.

vissE: a versatile tool to identify and visualise higher-order molecular phenotypes from functional enrichment analysis

Functional analysis of high throughput experiments using pathway analysis is now ubiquitous. Though powerful, these methods often produce thousands of redundant results owing to knowledgebase redundancies upstream. This scale of results hinders extensive exploration by biologists and can lead to investigator biases due to previous knowledge and expectations. To address this issue, we present vissE, a flexible network-based analysis and visualisation tool that organises information into semantic categories and provides various visualisation modules to characterise them with respect to the underlying data, thus providing a comprehensive view of the biological system. We demonstrate vissE's versatility by applying it to three different technologies: bulk, single-cell and spatial transcriptomics. Applying vissE to a factor analysis of a breast cancer spatial transcriptomic data, we identified stromal phenotypes that support tumour dissemination. Its adaptability allows vissE to enhance all existing gene-set enrichment and pathway analysis workflows, empowering biologists during molecular discovery.

MpTGA, together with MpNPR, regulates sexual reproduction and independently affects oil body formation in Marchantia polymorpha

In angiosperms, basic leucine-zipper (bZIP) TGACG-motif-binding (TGA) transcription factors (TFs) regulate developmental and stress-related processes, the latter often involving NON EXPRESSOR OF PATHOGENESIS-RELATED GENES (NPR) coregulator interactions. To gain insight into their functions in an early diverging land-plant lineage, the single MpTGA and sole MpNPR genes were investigated in the liverwort Marchantia polymorpha. We generated Marchantia MpTGA and MpNPR knockout and overexpression mutants and conducted morphological, transcriptomic and expression studies. Furthermore, we investigated MpTGA interactions with wild-type and mutagenized MpNPR and expanded our analyses including TGA TFs from two streptophyte algae. Mptga mutants fail to induce the switch from vegetative to reproductive development and lack gametangiophore formation. MpTGA and MpNPR proteins interact and Mpnpr mutant analysis reveals a novel coregulatory NPR role in sexual reproduction. Additionally, MpTGA acts independently of MpNPR as a repressor of oil body (OB) formation and can thereby affect herbivory. The single MpTGA TF exerts a dual role in sexual reproduction and OB formation in Marchantia. Common activities of MpTGA/MpNPR in sexual development suggest that coregulatory interactions were established after emergence of land-plant-specific NPR genes and contributed to the diversification of TGA TF functions during land-plant evolution.

Single-cell multi-omics profiling of human preimplantation embryos identifies cytoskeletal defects during embryonic arrest

Human in vitro fertilized embryos exhibit low developmental capabilities, and the mechanisms that underlie embryonic arrest remain unclear. Here using a single-cell multi-omics sequencing approach, we simultaneously analysed alterations in the transcriptome, chromatin accessibility and the DNA methylome in human embryonic arrest due to unexplained reasons. Arrested embryos displayed transcriptome disorders, including a distorted microtubule cytoskeleton, increased genomic instability and impaired glycolysis, which were coordinated with multiple epigenetic reprogramming defects. We identified Aurora A kinase (AURKA) repression as a cause of embryonic arrest. Mechanistically, arrested embryos induced through AURKA inhibition resembled the reprogramming abnormalities of natural embryonic arrest in terms of the transcriptome, the DNA methylome, chromatin accessibility and H3K4me3 modifications. Mitosis-independent sequential activation of the zygotic genome in arrested embryos showed that YY1 contributed to human major zygotic genome activation. Collectively, our study decodes the reprogramming abnormalities and mechanisms of human embryonic arrest and the key regulators of zygotic genome activation.

Cellular heterogeneity in red and melanized focal muscle changes in farmed Atlantic salmon (Salmo salar) visualized by spatial transcriptomics

Spatial transcriptomics is a technique that provides insight into gene expression profiles in tissue sections while retaining structural information. We have employed this method to study the pathological conditions related to red and melanized focal changes in farmed Atlantic salmon (Salmo salar). Our findings support a model where similar molecular mechanisms are involved in both red and melanized filet discolorations and genes associated with several relevant pathways show distinct expression patterns in both sample types. Interestingly, there appears to be significant cellular heterogeneity in the foci investigated when looking at gene expression patterns. Some of the genes that show differential spatial expression are involved in cellular processes such as hypoxia and immune responses, providing new insight into the nature of muscle melanization in Atlantic salmon.

The Current Situation and Development Prospect of Whole-Genome Screening

High-throughput genetic screening is useful for discovering critical genes or gene sequences that trigger specific cell functions and/or phenotypes. Loss-of-function genetic screening is mainly achieved through RNA interference (RNAi), CRISPR knock-out (CRISPRko), and CRISPR interference (CRISPRi) technologies. Gain-of-function genetic screening mainly depends on the overexpression of a cDNA library and CRISPR activation (CRISPRa). Base editing can perform both gain- and loss-of-function genetic screening. This review discusses genetic screening techniques based on Cas9 nuclease, including Cas9-mediated genome knock-out and dCas9-based gene activation and interference. We compare these methods with previous genetic screening techniques based on RNAi and cDNA library overexpression and propose future prospects and applications for CRISPR screening.

Curation of causal interactions mediated by genes associated with autism accelerates the understanding of gene-phenotype relationships underlying neurodevelopmental disorders

Autism spectrum disorder (ASD) comprises a large group of neurodevelopmental conditions featuring, over a wide range of severity and combinations, a core set of manifestations (restricted sociality, stereotyped behavior and language impairment) alongside various comorbidities. Common and rare variants in several hundreds of genes and regulatory regions have been implicated in the molecular pathogenesis of ASD along a range of causation evidence strength. Despite significant progress in elucidating the impact of few paradigmatic individual loci, such sheer complexity in the genetic architecture underlying ASD as a whole has hampered the identification of convergent actionable hubs hypothesized to relay between the vastness of risk alleles and the core phenotypes. In turn this has limited the development of strategies that can revert or ameliorate this condition, calling for a systems-level approach to probe the cross-talk of cooperating genes in terms of causal interaction networks in order to make convergences experimentally tractable and reveal their clinical actionability. As a first step in this direction, we have captured from the scientific literature information on the causal links between the genes whose variants have been associated with ASD and the whole human proteome. This information has been annotated in a computer readable format in the SIGNOR database and is made freely available in the resource website. To link this information to cell functions and phenotypes, we have developed graph algorithms that estimate the functional distance of any protein in the SIGNOR causal interactome to phenotypes and pathways. The main novelty of our approach resides in the possibility to explore the mechanistic links connecting the suggested gene-phenotype relations.

Epigenetic changes in sperm are associated with paternal and child quantitative autistic traits in an autism-enriched cohort

There is a need to consider paternal contributions to autism spectrum disorder (ASD) more strongly. Autism etiology is complex, and heritability is not explained by genetics alone. Understanding paternal gametic epigenetic contributions to autism could help fill this knowledge gap. In the present study, we explored whether paternal autistic traits, and the sperm epigenome, were associated with autistic traits in children at 36 months enrolled in the Early Autism Risk Longitudinal Investigation (EARLI) cohort. EARLI is a pregnancy cohort that recruited and enrolled pregnant women in the first half of pregnancy who already had a child with ASD. After maternal enrollment, EARLI fathers were approached and asked to provide a semen specimen. Participants were included in the present study if they had genotyping, sperm methylation data, and Social Responsiveness Scale (SRS) score data available. Using the CHARM array, we performed genome-scale methylation analyses on DNA from semen samples contributed by EARLI fathers. The SRS-a 65-item questionnaire measuring social communication deficits on a quantitative scale-was used to evaluate autistic traits in EARLI fathers (n = 45) and children (n = 31). We identified 94 significant child SRS-associated differentially methylated regions (DMRs), and 14 significant paternal SRS-associated DMRs (fwer p < 0.05). Many child SRS-associated DMRs were annotated to genes implicated in ASD and neurodevelopment. Six DMRs overlapped across the two outcomes (fwer p < 0.1), and, 16 DMRs overlapped with previous child autistic trait findings at 12 months of age (fwer p < 0.05). Child SRS-associated DMRs contained CpG sites independently found to be differentially methylated in postmortem brains of individuals with and without autism. These findings suggest paternal germline methylation is associated with autistic traits in 3-year-old offspring. These prospective results for autism-associated traits, in a cohort with a family history of ASD, highlight the potential importance of sperm epigenetic mechanisms in autism.

Human disease-specific cell signatures in non-lesional tissue in Multiple Sclerosis detected by single-cell and spatial transcriptomics

Recent investigations of cell type changes in Multiple Sclerosis (MS) using single-cell profiling methods have focused on active lesional and peri-lesional brain tissue, and have implicated a number of peripheral and central nervous system cell types. However, an important question is the extent to which so-called "normal-appearing" non-lesional tissue in individuals with MS accumulate changes over the lifespan. Here, we compared post-mortem non-lesional brain tissue from donors with a pathological or clinical diagnosis of MS from the Religious Orders Study or Rush Memory and Aging Project (ROSMAP) cohorts to age- and sex-matched brains from persons without MS (controls). We profiled three brain regions using single-nucleus RNA-seq: dorsolateral prefrontal cortex (DLPFC), normal appearing white matter (NAWM) and the pulvinar in thalamus (PULV), from 15 control individuals, 8 individuals with MS, and 5 individuals with other detrimental pathologies accompanied by demyelination, resulting in a total of 78 samples. We identified region- and cell type-specific differences in non-lesional samples from individuals diagnosed with MS and/or exhibiting secondary demyelination with other neurological conditions, as compared to control donors. These differences included lower proportions of oligodendrocytes with expression of myelination related genes MOBP, MBP, PLP1, as well as higher proportions of CRYAB+ oligodendrocytes in all three brain regions. Among microglial signatures, we identified subgroups that were higher in both demyelination (TMEM163+/ERC2+), as well as those that were specifically higher in MS donors (HIF1A+/SPP1+) and specifically in donors with secondary demyelination (SOCS6+/MYO1E+), in both white and grey matter. To validate our findings, we generated Visium spatial transcriptomics data on matched tissue from 13 donors, and recapitulated our observations of gene expression differences in oligodendrocytes and microglia. Finally, we show that some of the differences observed between control and MS donors in NAWM mirror those previously reported between control WM and active lesions in MS donors. Overall, our investigation sheds additional light on cell type- and disease-specific differences present even in non-lesional white and grey matter tissue, highlighting widespread cellular signatures that may be associated with downstream pathological changes.

Mapping protein states and interactions across the tree of life with co-fractionation mass spectrometry

We present CFdb, a harmonized resource of interaction proteomics data from 411 co-fractionation mass spectrometry (CF-MS) datasets spanning 21,703 fractions. Meta-analysis of this resource charts protein abundance, phosphorylation, and interactions throughout the tree of life, including a reference map of the human interactome. We show how large-scale CF-MS data can enhance analyses of individual CF-MS datasets, and exemplify this strategy by mapping the honey bee interactome.

AAGGG repeat expansions trigger RFC1-independent synaptic dysregulation in human CANVAS Neurons

Cerebellar ataxia with neuropathy and vestibular areflexia syndrome (CANVAS) is a late onset, recessively inherited neurodegenerative disorder caused by biallelic, non-reference pentameric AAGGG(CCCTT) repeat expansions within the second intron of replication factor complex subunit 1 (RFC1). To investigate how these repeats cause disease, we generated CANVAS patient induced pluripotent stem cell (iPSC) derived neurons (iNeurons) and utilized calcium imaging and transcriptomic analysis to define repeat-elicited gain-of-function and loss-of-function contributions to neuronal toxicity. AAGGG repeat expansions do not alter neuronal RFC1 splicing, expression, or DNA repair pathway functions. In reporter assays, AAGGG repeats are translated into pentapeptide repeat proteins that selectively accumulate in CANVAS patient brains. However, neither these proteins nor repeat RNA foci were detected in iNeurons, and overexpression of these repeats in isolation did not induce neuronal toxicity. CANVAS iNeurons exhibit defects in neuronal development and diminished synaptic connectivity that is rescued by CRISPR deletion of a single expanded allele. These phenotypic deficits were not replicated by knockdown of RFC1 in control neurons and were not rescued by ectopic expression of RFC1. These findings support a repeat-dependent but RFC1-independent cause of neuronal dysfunction in CANVAS, with important implications for therapeutic development in this currently untreatable condition.

Single-cell immune profiling reveals markers of emergency myelopoiesis that distinguish severe from mild respiratory syncytial virus disease in infants

Whereas most infants infected with respiratory syncytial virus (RSV) show no or only mild symptoms, an estimated 3 million children under five are hospitalized annually due to RSV disease. This study aimed to investigate biological mechanisms and associated biomarkers underlying RSV disease heterogeneity in young infants, enabling the potential to objectively categorize RSV-infected infants according to their medical needs. Immunophenotypic and functional profiling demonstrated the emergence of immature and progenitor-like neutrophils, proliferative monocytes (HLA-DRLow , Ki67+), impaired antigen-presenting function, downregulation of T cell response and low abundance of HLA-DRLow B cells in severe RSV disease. HLA-DRLow monocytes were found as a hallmark of RSV-infected infants requiring hospitalization. Complementary transcriptomics identified genes associated with disease severity and pointed to the emergency myelopoiesis response. These results shed new light on mechanisms underlying the pathogenesis and development of severe RSV disease and identified potential new candidate biomarkers for patient stratification.

Higher expression of programmed cell death 4 (PDCD4) in acute myeloid leukemia is associated with better prognosis after chemotherapy

Acute myeloid leukemia (AML) is a common heterogeneous malignancy. Novel molecular markers aid diagnosis, patient sub-categorization, and optimal clinical decisions. Here, we explored the prognostic implications associated with the expression of the programmed cell death (PDCD) family of molecules in AML patients. Based on the findings from the TCGA and OHSU cohorts, we observed that the mRNA abundance of PDCD4 is significantly higher compared to other molecules within the PDCD family. Furthermore, high expression of PDCD4 was associated with predicted long-term patient survival in diagnosed AML patients. In the chemotherapy group, patients with high PDCD4 expression showed a tendency toward longer overall survival (OS) (P = 0.0266) and event-free survival (EFS) (P = 0.0008). High PDCD4 levels served as a favorable independent predictor for both OS and EFS in AML patients. However, subgroup analyses in the hematopoietic stem cell transplantation (HSCT) group revealed no significant difference in OS or EFS between individuals with high and low PDCD4 expression. Furthermore, in the low PDCD4 expression group, AML patients who underwent HSCT experienced improved survival outcomes (P = 0.0015), helping to mitigate the unfavorable prognosis associated with PDCD4 downregulation. Conversely, in the high PDCD4 expression group, HSCT offered a notable short-term survival advantage, while patients with high PDCD4 expression responded favorably to long-term survival through chemotherapy. Biological function enrichment showed that the expression of PDCD4 was correlated with complement and coagulation cascades, cell receptor signaling pathways, and cholesterol metabolism. The findings from this study will aid in better categorizing heterogeneous AML patients and guiding more appropriate clinical decision-making.

Intergenerational effects of early life-stage temperature modulation on gene expression and DNA methylation in Atlantic cod (Gadus morhua)

After suffering several collapses, the cod farming industry is now in the process of trying to re-establish itself. We have used material from Norway's National Cod Breeding Program to study how different early life-stage temperature regimes affect DNA methylation and gene expression. Long-term effects were detected by sampling fish several weeks after the end of differential treatments, and offspring from the different exposure groups was also sampled. Many overlapping genes were found between the different exposure groups and generations, coupled with genes associated with differential CpG methylation levels. Genes involved in muscle fibre development, general metabolic processes and formation of deformities were significantly affected, and genes relevant for intergenerational transfer of epigenetic marks were also detected. We believe the use of environmental cues can be a useful strategy for improving the production of Atlantic cod.

Sis2 regulates yeast replicative lifespan in a dose-dependent manner

Application of microfluidic platforms facilitated high-precision measurements of yeast replicative lifespan (RLS); however, comparative quantification of lifespan across strain libraries has been missing. Here we microfluidically measure the RLS of 307 yeast strains, each deleted for a single gene. Despite previous reports of extended lifespan in these strains, we found that 56% of them did not actually live longer than the wild-type; while the remaining 44% showed extended lifespans, the degree of extension was often different from what was previously reported. Deletion of SIS2 gene led to the largest RLS increase observed. Sis2 regulated yeast lifespan in a dose-dependent manner, implying a role for the coenzyme A biosynthesis pathway in lifespan regulation. Introduction of the human PPCDC gene in the sis2Δ background neutralized the lifespan extension. RNA-seq experiments revealed transcriptional increases in cell-cycle machinery components in sis2Δ background. High-precision lifespan measurement will be essential to elucidate the gene network governing lifespan.

An Artificial Intelligence Prediction Model of Insulin Sensitivity, Insulin Resistance, and Diabetes Using Genes Obtained through Differential Expression

Insulin is a powerful pleiotropic hormone that affects processes such as cell growth, energy expenditure, and carbohydrate, lipid, and protein metabolism. The molecular mechanisms by which insulin regulates muscle metabolism and the underlying defects that cause insulin resistance have not been fully elucidated. This study aimed to perform a microarray data analysis to find differentially expressed genes. The analysis has been based on the data of a study deposited in Gene Expression Omnibus (GEO) with the identifier "GSE22309". The selected data contain samples from three types of patients after taking insulin treatment: patients with diabetes (DB), patients with insulin sensitivity (IS), and patients with insulin resistance (IR). Through an analysis of omics data, 20 genes were found to be differentially expressed (DEG) between the three possible comparisons obtained (DB vs. IS, DB vs. IR, and IS vs. IR); these data sets have been used to develop predictive models through machine learning (ML) techniques to classify patients with respect to the three categories mentioned previously. All the ML techniques present an accuracy superior to 80%, reaching almost 90% when unifying IR and DB categories.

Sctensor detects many-to-many cell-cell interactions from single cell RNA-sequencing data

BACKGROUND

Complex biological systems are described as a multitude of cell-cell interactions (CCIs). Recent single-cell RNA-sequencing studies focus on CCIs based on ligand-receptor (L-R) gene co-expression but the analytical methods are not appropriate to detect many-to-many CCIs.

RESULTS

In this work, we propose scTensor, a novel method for extracting representative triadic relationships (or hypergraphs), which include ligand-expression, receptor-expression, and related L-R pairs.

CONCLUSIONS

Through extensive studies with simulated and empirical datasets, we have shown that scTensor can detect some hypergraphs that cannot be detected using conventional CCI detection methods, especially when they include many-to-many relationships. scTensor is implemented as a freely available R/Bioconductor package.

Proteomics and Metabolomics Reveal that an Abundant α-Glucosidase Drives Sorghum Fermentability for Beer Brewing

Sorghum (Sorghum bicolor), a grass native to Africa, is a popular alternative to barley for brewing beer. The importance of sorghum to beer brewing is increasing because it is a naturally gluten-free cereal, and climate change is expected to cause a reduction in the production of barley over the coming decades. However, there are challenges associated with the use of sorghum instead of barley in beer brewing. Here, we used proteomics and metabolomics to gain insights into the sorghum brewing process to advise processes for efficient beer production from sorghum. We found that during malting, sorghum synthesizes the amylases and proteases necessary for brewing. Proteomics revealed that mashing with sorghum malt required higher temperatures than barley malt for efficient protein solubilization. Both α- and β-amylase were considerably less abundant in sorghum wort than in barley wort, correlating with lower maltose concentrations in sorghum wort. However, metabolomics revealed higher glucose concentrations in sorghum wort than in barley wort, consistent with the presence of an abundant α-glucosidase detected by proteomics in sorghum malt. Our results indicate that sorghum can be a viable grain for industrial fermented beverage production, but that its use requires careful process optimization for efficient production of fermentable wort and high-quality beer.

COX17 acetylation via MOF-KANSL complex promotes mitochondrial integrity and function

Reversible acetylation of mitochondrial proteins is a regulatory mechanism central to adaptive metabolic responses. Yet, how such functionally relevant protein acetylation is achieved remains unexplored. Here we reveal an unprecedented role of the MYST family lysine acetyltransferase MOF in energy metabolism via mitochondrial protein acetylation. Loss of MOF-KANSL complex members leads to mitochondrial defects including fragmentation, reduced cristae density and impaired mitochondrial electron transport chain complex IV integrity in primary mouse embryonic fibroblasts. We demonstrate COX17, a complex IV assembly factor, as a bona fide acetylation target of MOF. Loss of COX17 or expression of its non-acetylatable mutant phenocopies the mitochondrial defects observed upon MOF depletion. The acetylation-mimetic COX17 rescues these defects and maintains complex IV activity even in the absence of MOF, suggesting an activatory role of mitochondrial electron transport chain protein acetylation. Fibroblasts from patients with MOF syndrome who have intellectual disability also revealed respiratory defects that could be restored by alternative oxidase, acetylation-mimetic COX17 or mitochondrially targeted MOF. Overall, our findings highlight the critical role of MOF-KANSL complex in mitochondrial physiology and provide new insights into MOF syndrome.