Leveraging interindividual variability in threat conditioning of inbred mice to model trait anxiety.

Trait anxiety is a major risk factor for stress-induced and anxiety disorders in humans. However, animal models accounting for the interindividual variability in stress vulnerability are largely lacking. Moreover, the pervasive bias of using mostly male animals in preclinical studies poorly reflects the increased prevalence of psychiatric disorders in women. Using the threat imminence continuum theory, we designed and validated an auditory aversive conditioning-based pipeline in both female and male mice. We operationalised trait anxiety by harnessing the naturally occurring variability of defensive freezing responses combined with a model-based clustering strategy. While sustained freezing during prolonged retrieval sessions was identified as an anxiety-endophenotype biomarker in both sexes, females were consistently associated with an increased freezing response. RNA-sequencing of CeA, BLA, ACC, and BNST revealed massive differences in phasic and sustained responders' transcriptomes, correlating with transcriptomic signatures of psychiatric disorders, particularly post-traumatic stress disorder (PTSD). Moreover, we detected significant alterations in the excitation/inhibition balance of principal neurons in the lateral amygdala. These findings provide compelling evidence that trait anxiety in inbred mice can be leveraged to develop translationally relevant preclinical models to investigate mechanisms of stress susceptibility in a sex-specific manner.

Long non-coding RNAs expression and regulation across different brain regions in primates.

Human and non-human primates have strikingly similar genomes, but they strongly differ in many brain-based processes (e.g., behaviour and cognition). While the functions of protein-coding genes have been extensively studied, rather little is known about the role of non-coding RNAs such as long non-coding RNAs (lncRNAs). Here, we predicted lncRNAs and analysed their expression pattern across different brain regions of human and non-human primates (chimpanzee, gorilla, and gibbon). Our analysis identified shared orthologous and non-orthologous lncRNAs, showing striking differences in the genomic features. Differential expression analysis of the shared orthologous lncRNAs from humans and chimpanzees revealed distinct expression patterns in subcortical regions (striatum, hippocampus) and neocortical areas while retaining a homogeneous expression in the cerebellum. Co-expression analysis of lncRNAs and protein-coding genes revealed massive proportions of co-expressed pairs in neocortical regions of humans compared to chimpanzees. Network analysis of co-expressed pairs revealed the distinctive role of the hub-acting orthologous lncRNAs in a region- and species-specific manner. Overall, our study provides novel insight into lncRNA driven gene regulatory landscape, neural regulation, brain evolution, and constitutes a resource for primate's brain lncRNAs.

Characterization of transcriptional profiles associated with stress-induced neuronal activation in Arc-GFP mice.

Chronic stress has become a predominant factor associated with a variety of psychiatric disorders, such as depression and anxiety, in both human and animal models. Although multiple studies have looked at transcriptional changes after social defeat stress, these studies primarily focus on bulk tissues, which might dilute important molecular signatures of social interaction in activated cells. In this study, we employed the Arc-GFP mouse model in conjunction with chronic social defeat (CSD) to selectively isolate activated nuclei (AN) populations in the ventral hippocampus (vHIP) and prefrontal cortex (PFC) of resilient and susceptible animals. Nuclear RNA-seq of susceptible vs. resilient populations revealed distinct transcriptional profiles linked predominantly with neuronal and synaptic regulation mechanisms. In the vHIP, susceptible AN exhibited increased expression of genes related to the cytoskeleton and synaptic organization. At the same time, resilient AN showed upregulation of cell adhesion genes and differential expression of major glutamatergic subunits. In the PFC, susceptible mice exhibited upregulation of synaptotagmins and immediate early genes (IEGs), suggesting a potentially over-amplified neuronal activity state. Our findings provide a novel view of stress-exposed neuronal activation and the molecular response mechanisms in stress-susceptible vs. resilient animals, which may have important implications for understanding mental resilience.

Nanopore adaptive sampling of a metagenomic sample derived from a human monkeypox case.

In 2022, a series of human monkeypox cases in multiple countries led to the largest and most widespread outbreak outside the known endemic areas. Setup of proper genomic surveillance is of utmost importance to control such outbreaks. To this end, we performed Nanopore (PromethION P24) and Illumina (NextSeq. 2000) Whole Genome Sequencing (WGS) of a monkeypox sample. Adaptive sampling was applied for in silico depletion of the human host genome, allowing for the enrichment of low abundance viral DNA without a priori knowledge of sample composition. Nanopore sequencing allowed for high viral genome coverage, tracking of sample composition during sequencing, strain determination, and preliminary assessment of mutational pattern. In addition to that, only Nanopore data allowed us to resolve the entire monkeypox virus genome, with respect to two structural variants belonging to the genes OPG015 and OPG208. These SVs in important host range genes seem stable throughout the outbreak and are frequently misassembled and/or misannotated due to the prevalence of short read sequencing or short read first assembly. Ideally, standalone standard Illumina sequencing should not be used for Monkeypox WGS and de novo assembly, since it will obfuscate the structure of the genome, which has an impact on the quality and completeness of the genomes deposited in public databases and thus possibly on the ability to evaluate the complete genetic reason for the host range change of monkeypox in the current pandemic.

A Call for Gene Expression Analysis in Whole Blood of Patients With Rheumatoid Arthritis (RA) as a Biomarker for RA-Associated Interstitial Lung Disease.

OBJECTIVE: Rheumatoid arthritis (RA)-associated interstitial lung disease (ILD) is one of the most common and prognostic organ manifestations of RA. Therefore, to allow effective treatment, it is of crucial importance to diagnose RA-ILD at the earliest possible stage. So far, the gold standard of early detection has been high-resolution computed tomography (HRCT) of the lungs. This procedure involves considerable radiation exposure for the patient and is therefore unsuitable as a routine screening measure for ethical reasons. Here, we propose the analysis of characteristic gene expression patterns as a biomarker to aid in the early detection and initiation of appropriate, possibly antifibrotic, therapy. METHODS: To investigate unique molecular patterns of RA-ILD, whole blood samples were taken from 12 female patients with RA-ILD (n = 7) or RA (n = 5). The RNA was extracted, sequenced by RNA-Seq, and analyzed for characteristic differences in the gene expression patterns between patients with RA-ILD and those with RA without ILD. RESULTS: The differential gene expression analysis revealed 9 significantly upregulated genes in RA-ILD compared to RA without ILD: arginase 1 (ARG1), thymidylate synthetase (TYMS), sortilin 1 (SORT1), marker of proliferation Ki-67 (MKI67), olfactomedin 4 (OLFM4), baculoviral inhibitor of apoptosis repeat containing 5 (BIRC5), membrane spanning 4-domains A4A (MS4A4A), C-type lectin domain family 12 member A (CLEC12A), and the long intergenic nonprotein coding RNA (LINC02967). CONCLUSION: All gene products of these genes (except for LINC02967) are known from the literature to be involved in the pathogenesis of fibrosis. Further, for some, a contribution to the development of pulmonary fibrosis has even been demonstrated in experimental studies. Therefore, the results presented here provide an encouraging perspective for using specific gene expression patterns as biomarkers for the early detection and differential diagnosis of RA-ILD as a routine screening test.

N1-methylation of adenosine (m1A) in ND5 mRNA leads to complex I dysfunction in Alzheimer's disease.

One mechanism of particular interest to regulate mRNA fate post-transcriptionally is mRNA modification. Especially the extent of m1A mRNA methylation is highly discussed due to methodological differences. However, one single m1A site in mitochondrial ND5 mRNA was unanimously reported by different groups. ND5 is a subunit of complex I of the respiratory chain. It is considered essential for the coupling of oxidation and proton transport. Here we demonstrate that this m1A site might be involved in the pathophysiology of Alzheimer's disease (AD). One of the pathological hallmarks of this neurodegenerative disease is mitochondrial dysfunction, mainly induced by Amyloid ß (Aß). Aß mainly disturbs functions of complex I and IV of the respiratory chain. However, the molecular mechanism of complex I dysfunction is still not fully understood. We found enhanced m1A methylation of ND5 mRNA in an AD cell model as well as in AD patients. Formation of this m1A methylation is catalyzed by increased TRMT10C protein levels, leading to translation repression of ND5. As a consequence, here demonstrated for the first time, TRMT10C induced m1A methylation of ND5 mRNA leads to mitochondrial dysfunction. Our findings suggest that this newly identified mechanism might be involved in Aß-induced mitochondrial dysfunction.

NanopoReaTA: a user-friendly tool for nanopore-seq real-time transcriptional analysis.

SUMMARY: Oxford Nanopore Technologies' (ONT) sequencing platform offers an excellent opportunity to perform real-time analysis during sequencing. This feature allows for early insights into experimental data and accelerates a potential decision-making process for further analysis, which can be particularly relevant in the clinical context. Although some tools for the real-time analysis of DNA-sequencing data already exist, there is currently no application available for differential transcriptome data analysis designed for scientists or physicians with limited bioinformatics knowledge. Here, we introduce NanopoReaTA, a user-friendly real-time analysis toolbox for RNA-sequencing data from ONT. Sequencing results from a running or finished experiment are processed through an R Shiny-based graphical user interface with an integrated Nextflow pipeline for whole transcriptome or gene-specific analyses. NanopoReaTA provides visual snapshots of a sequencing run in progress, thus enabling interactive sequencing and rapid decision making that could also be applied to clinical cases. AVAILABILITY AND IMPLEMENTATION: Github https://github.com/AnWiercze/NanopoReaTA; Zenodo https://doi.org/10.5281/zenodo.8099825.

NAPE-PLD deletion in stress-TRAPed neurons results in an anxiogenic phenotype.

Anandamide (AEA) is an endogenous ligand of the cannabinoid CB1 and CB2 receptors, being a component of the endocannabinoid signaling system, which supports the maintenance or regaining of neural homeostasis upon internal and external challenges. AEA is thought to play a protective role against the development of pathological states after prolonged stress exposure, including depression and generalized anxiety disorder. Here, we used the chronic social defeat (CSD) stress as an ethologically valid model of chronic stress in male mice. We characterized a genetically modified mouse line where AEA signaling was reduced by deletion of the gene encoding the AEA synthesizing enzyme N-acyl-phosphatidylethanolamine-hydrolyzing phospholipase D (NAPE-PLD) specifically in neurons activated at the time of CSD stress. One week after the stress, the phenotype was assessed in behavioral tests and by molecular analyses. We found that NAPE-PLD deficiency in neurons activated during the last three days of CSD stress led to an increased anxiety-like behavior. Investigating the molecular mechanisms underlying this phenotype may suggest three main altered pathways to be affected: (i) desensitization of the negative feedback loop of the hypothalamic-pituitary-adrenal axis, (ii) disinhibition of the amygdala by the prefrontal cortex, and (iii) altered neuroplasticity in the hippocampus and prefrontal cortex.

Nuclei on the Rise: When Nuclei-Based Methods Meet Next-Generation Sequencing.

In the last decade, we have witnessed an upsurge in nuclei-based studies, particularly coupled with next-generation sequencing. Such studies aim at understanding the molecular states that exist in heterogeneous cell populations by applying increasingly more affordable sequencing approaches, in addition to optimized methodologies developed to isolate and select nuclei. Although these powerful new methods promise unprecedented insights, it is important to understand and critically consider the associated challenges. Here, we provide a comprehensive overview of the rise of nuclei-based studies and elaborate on their advantages and disadvantages, with a specific focus on their utility for transcriptomic sequencing analyses. Improved designs and appropriate use of the various experimental strategies will result in acquiring biologically accurate and meaningful information.

On cheap entropy-sparsified regression learning.

Regression learning is one of the long-standing problems in statistics, machine learning, and deep learning (DL). We show that writing this problem as a probabilistic expectation over (unknown) feature probabilities - thus increasing the number of unknown parameters and seemingly making the problem more complex-actually leads to its simplification, and allows incorporating the physical principle of entropy maximization. It helps decompose a very general setting of this learning problem (including discretization, feature selection, and learning multiple piece-wise linear regressions) into an iterative sequence of simple substeps, which are either analytically solvable or cheaply computable through an efficient second-order numerical solver with a sublinear cost scaling. This leads to the computationally cheap and robust non-DL second-order Sparse Probabilistic Approximation for Regression Task Analysis (SPARTAn) algorithm, that can be efficiently applied to problems with millions of feature dimensions on a commodity laptop, when the state-of-the-art learning tools would require supercomputers. SPARTAn is compared to a range of commonly used regression learning tools on synthetic problems and on the prediction of the El Niño Southern Oscillation, the dominant interannual mode of tropical climate variability. The obtained SPARTAn learners provide more predictive, sparse, and physically explainable data descriptions, clearly discerning the important role of ocean temperature variability at the thermocline in the equatorial Pacific. SPARTAn provides an easily interpretable description of the timescales by which these thermocline temperature features evolve and eventually express at the surface, thereby enabling enhanced predictability of the key drivers of the interannual climate.

Acral lamellar ichthyosis with amino acid substitution in the C-terminus of keratin 2.

BACKGROUND: Most cases of hereditary ichthyoses present with generalized scaling and skin dryness. However, in some cases skin involvement is restricted to particular body regions as in acral lamellar ichthyosis. OBJECTIVES: We report on the genetic basis of acral ichthyosis in two families presenting with a similar phenotype. METHODS: Genetic testing was performed by targeted next generation sequencing and whole-exome sequencing. For identity-by-descent analysis, the parents were genotyped and data analysis was performed with the Chromosome Analysis Suite Software. RT-PCR with RNA extracted from skin samples was used to analyse the effect of variants on splicing. RESULTS: Genetic testing identified a few heterozygous variants, but only the variant in KRT2 c.1912 T > C, p.Phe638Leu segregated with the disease and remained the strongest candidate. Pairwise identity-by-descent analysis revealed no indication of family relationship. Phenylalanine 638 is the second last amino acid upstream of the termination codon in the tail of K2, and substitution to leucine is predicted as probably damaging. Assessment of the variant is difficult, in part due to the lack of crystal structures of this region. CONCLUSIONS: Altogether, we show that a type of autosomal dominant acral ichthyosis is most probably caused by an amino acid substitution in the C-terminus of keratin 2.

A Resilience Related Glial-Neurovascular Network Is Transcriptionally Activated after Chronic Social Defeat in Male Mice.

Upon chronic stress, a fraction of individuals shows stress resilience, which can prevent long-term mental dysfunction. The underlying molecular mechanisms are complex and have not yet been fully understood. In this study, we performed a data-driven behavioural stratification together with single-cell transcriptomics of the hippocampus in a mouse model of chronic social defeat stress. Our work revealed that in a sub-group exhibiting molecular responses upon chronic stress, the dorsal hippocampus is particularly involved in neuroimmune responses, angiogenesis, myelination, and neurogenesis, thereby enabling brain restoration and homeostasis after chronic stress. Based on these molecular insights, we applied rapamycin after the stress as a proof-of-concept pharmacological intervention and were able to substantially increase stress resilience. Our findings serve as a data resource and can open new avenues for further understanding of molecular processes underlying stress response and for targeted interventions supporting resilience.

The Big Picture of Neurodegeneration: A Meta Study to Extract the Essential Evidence on Neurodegenerative Diseases in a Network-Based Approach.

The common features of all neurodegenerative diseases, including Alzheimer's disease, Parkinson's disease, Amyotrophic Lateral Sclerosis (ALS), and Huntington's disease, are the accumulation of aggregated and misfolded proteins and the progressive loss of neurons, leading to cognitive decline and locomotive dysfunction. Still, they differ in their ultimate manifestation, the affected brain region, and the kind of proteinopathy. In the last decades, a vast number of processes have been described as associated with neurodegenerative diseases, making it increasingly harder to keep an overview of the big picture forming from all those data. In this meta-study, we analyzed genomic, transcriptomic, proteomic, and epigenomic data of the aforementioned diseases using the data of 234 studies in a network-based approach to study significant general coherences but also specific processes in individual diseases or omics levels. In the analysis part, we focus on only some of the emerging findings, but trust that the meta-study provided here will be a valuable resource for various other researchers focusing on specific processes or genes contributing to the development of neurodegeneration.

Low-Cost Probabilistic 3D Denoising with Applications for Ultra-Low-Radiation Computed Tomography.

We propose a pipeline for synthetic generation of personalized Computer Tomography (CT) images, with a radiation exposure evaluation and a lifetime attributable risk (LAR) assessment. We perform a patient-specific performance evaluation for a broad range of denoising algorithms (including the most popular deep learning denoising approaches, wavelets-based methods, methods based on Mumford−Shah denoising, etc.), focusing both on accessing the capability to reduce the patient-specific CT-induced LAR and on computational cost scalability. We introduce a parallel Probabilistic Mumford−Shah denoising model (PMS) and show that it markedly-outperforms the compared common denoising methods in denoising quality and cost scaling. In particular, we show that it allows an approximately 22-fold robust patient-specific LAR reduction for infants and a 10-fold LAR reduction for adults. Using a normal laptop, the proposed algorithm for PMS allows cheap and robust (with a multiscale structural similarity index >90%) denoising of very large 2D videos and 3D images (with over 107 voxels) that are subject to ultra-strong noise (Gaussian and non-Gaussian) for signal-to-noise ratios far below 1.0. The code is provided for open access.

Selective targeting of chronic social stress-induced activated neurons identifies neurogenesis-related genes to be associated with resilience in female mice.

Prolonged social stress is a major cause for depression in humans and is associated with a wide range of subsequent pathophysiological changes such as elevated blood pressure. A routinely used model for investigating this kind of stress in mice is the chronic social defeat paradigm where a smaller intruder is exposed to an aggressive inhabitant of a home cage. This model is restricted to males and includes a high proportion of physical stress that might e.g., interfere with immunological aspects of the stress. The prevalence of depression in humans is even higher in women than in men. Therefore, expanding models to female individuals is desirable. We here tested the social instability model as a tool for administering chronic social stress to female C57BL/6J mice and analyzed short-term as well as long-lasting effects. Animals were housed in groups of four and were shuffled two times a week, resulting in a permanent re-structuration of their social hierarchy. While directly after the stress exposure, serum corticosterone was elevated, increased body weight and fat deposits were observed in stressed mice even one year after discontinuation of the stress. At the behavioral level, animals could be stratified into resilient and susceptible animals directly post-stress, but those subgroups were not distinguishable any more in the long-term analysis. To identify molecular contributors to resilience in the here presented social instability induced stress model, Arc-activity dependent trapping of neurons was conducted in Arc-creERT2/sun1sfGFP mice. RNA samples derived from activated nuclei from the ventral hippocampus, a brain region involved in stress-regulation during attacks or explorative behavior of mice, were subjected to a neurogenesis pathway array. While several genes were differentially regulated by stress, in particular, artemin, a neurotrophic factor was upregulated in resilient versus susceptible individuals.

Long-term molecular differences between resilient and susceptible mice after a single traumatic exposure.

BACKGROUND AND PURPOSE: Post-traumatic stress disorder (PTSD) is a heterogeneous disorder induced by trauma, resulting in severe long-term impairments of an individual's mental health. PTSD does not develop in every individual and, thus, some individuals are more resilient. However, the underlying molecular mechanisms are poorly understood. Here, we aimed to elucidate these processes. EXPERIMENTAL APPROACH: We used a single-trauma PTSD model in mice to induce long-term maladaptive behaviours and profiled the mice 4 weeks after trauma into resilient or susceptible individuals. The classification of phenotype was based on individual responses in different behavioural experiments. We analysed microbiome, circulating endocannabinoids, and long-term changes in brain phospholipid and transcript levels. KEY RESULTS: We found many molecular differences between resilient and susceptible individuals across multiple molecular domains, including lipidome, transcriptome and gut microbiome. Some differences were stable even several weeks after the trauma, indicating the long-term impact of traumatic stimuli on the organism's physiology. Furthermore, the integration of these multilayered molecular data revealed that resilient and susceptible individuals have very distinct molecular signatures across various physiological systems. CONCLUSION AND IMPLICATIONS: Trauma induced individual-specific behavioural responses that, in combination with a longitudinal characterisation of mice, could be used to identify distinct sub-phenotypes within the trauma-exposed group. These groups differed significantly not only in their behaviour but also in specific molecular aspects across a variety of tissues and brain regions. This approach may reveal new targets and predictive biomarkers for the pharmacological treatment and prognosis of stress-related disorders. LINKED ARTICLES: This article is part of a themed issue on New discoveries and perspectives in mental and pain disorders. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v179.17/issuetoc.

beeRapp: an R shiny app for automated high-throughput explorative analysis of multivariate behavioral data.

Summary: Animal behavioral studies typically generate high-dimensional datasets consisting of multiple correlated outcome measures across distinct or related behavioral domains. Here, we introduce the BEhavioral Explorative analysis R shiny APP (beeRapp) that facilitates explorative and inferential analysis of behavioral data in a high-throughput fashion. By employing an intuitive and user-friendly graphical user interface, beeRapp empowers behavioral scientists without programming and data science expertise to perform clustering, dimensionality reduction, correlational and inferential statistics and produce up to thousands of high-quality output plots visualizing results in a standardized and automated way. Availability and implementation: The code and data underlying this article are available at https://github.com/anmabu/beeRapp.

The gut microbiota instructs the hepatic endothelial cell transcriptome.

The gut microbiota affects remote organ functions but its impact on organotypic endothelial cell (EC) transcriptomes remains unexplored. The liver endothelium encounters microbiota-derived signals and metabolites via the portal circulation. To pinpoint how gut commensals affect the hepatic sinusoidal endothelium, a magnetic cell sorting protocol, combined with fluorescence-activated cell sorting, was used to isolate hepatic sinusoidal ECs from germ-free (GF) and conventionally raised (CONV-R) mice for transcriptome analysis by RNA sequencing. This resulted in a comprehensive map of microbiota-regulated hepatic EC-specific transcriptome profiles. Gene Ontology analysis revealed that several functional processes in the hepatic endothelium were affected. The absence of microbiota influenced the expression of genes involved in cholesterol flux and angiogenesis. Specifically, genes functioning in hepatic endothelial sphingosine metabolism and the sphingosine-1-phosphate pathway showed drastically increased expression in the GF state. Our analyses reveal a prominent role for the microbiota in shaping the transcriptional landscape of the hepatic endothelium.

IntelliPy: a GUI for analyzing IntelliCage data.

SUMMARY: The IntelliCage systems offer the possibility to conduct long-term behavioral experiments on mice in social groups without human intervention. Although this setup provides new findings, only about 150 studies with the IntelliCage system have been published in the last two decades, which is also caused by the challenging problems of processing and handling the large and heterogeneous amounts of captured data. This application note introduces the Python-GUI IntelliPy, especially designed for users not very experienced in using programming languages. IntelliPy allows users to quickly analyze the IntelliCage output in a user-friendly way, thus making the systems more accessible to a broader audience. AVAILABILITY AND IMPLEMENTATION: https://github.com/NiRuff/IntelliPy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

N6 -Methyladenosine Modification in Chronic Stress Response Due to Social Hierarchy Positioning of Mice.

Appropriately responding to stressful events is essential for maintaining health and well-being of any organism. Concerning social stress, the response is not always as straightforward as reacting to physical stressors, e.g., extreme heat, and thus has to be balanced subtly. Particularly, regulatory mechanisms contributing to gaining resilience in the face of mild social stress are not fully deciphered yet. We employed an intrinsic social hierarchy stress paradigm in mice of both sexes to identify critical factors for potential coping strategies. While global transcriptomic changes could not be observed in male mice, several genes previously reported to be involved in synaptic plasticity, learning, and anxiety-like behavior were differentially regulated in female mice. Moreover, changes in N6-methyladenosine (m6A)-modification of mRNA occurred associated with corticosterone level in both sexes with, e.g., increased global amount in submissive female mice. In accordance with this, METTL14 and WTAP, subunits of the methyltransferase complex, showed elevated levels in submissive female mice. N6-adenosyl-methylation is the most prominent type of mRNA methylation and plays a crucial role in processes such as metabolism, but also response to physical stress. Our findings underpin its essential role by also providing a link to social stress evoked by hierarchy building within same-sex groups. As recently, search for small molecule modifiers for the respective class of RNA modifying enzymes has started, this might even lead to new therapeutic approaches against stress disorders.