GSAn: an alternative to enrichment analysis for annotating gene sets

Transcriptomic profiling of reward and sensory brain areas in perinatal fentanyl exposed juvenile mice

Use of the synthetic opioid fentanyl increased ~300% in the last decade, including among women of reproductive ages. Adverse neonatal outcomes and long-term behavioral disruptions are associated with perinatal opioid exposure. Our previous work demonstrated that perinatal fentanyl exposed mice displayed enhanced negative affect and somatosensory circuit and behavioral disruptions during adolescence. However, little is known about molecular adaptations across brain regions that underlie these outcomes. We performed RNA sequencing across three reward and two sensory brain areas to study transcriptional programs in perinatal fentanyl exposed juvenile mice. Pregnant dams received 10 μg/ml fentanyl in the drinking water from embryonic day 0 (E0) through gestational periods until weaning at postnatal day 21 (P21). RNA was extracted from nucleus accumbens (NAc), prelimbic cortex (PrL), ventral tegmental area (VTA), somatosensory cortex (S1) and ventrobasal thalamus (VBT) from perinatal fentanyl exposed mice of both sexes at P35. RNA sequencing was performed, followed by analysis of differentially expressed genes (DEGs) and gene co-expression networks. Transcriptome analysis revealed DEGs and gene modules significantly associated with exposure to perinatal fentanyl in a sex-wise manner. The VTA had the most DEGs, while robust gene enrichment occurred in NAc. Genes enriched in mitochondrial respiration were pronounced in NAc and VTA of perinatal fentanyl exposed males, extracellular matrix (ECM) and neuronal migration enrichment were pronounced in NAc and VTA of perinatal fentanyl exposed males, while genes associated with vesicular cycling and synaptic signaling were markedly altered in NAc of perinatal fentanyl exposed female mice. In sensory areas from perinatal fentanyl exposed females, we found alterations in mitochondrial respiration, synaptic and ciliary organization processes. Our findings demonstrate distinct transcriptomes across reward and sensory brain regions, with some showing discordance between sexes. These transcriptome adaptations may underlie structural, functional, and behavioral changes observed in perinatal fentanyl exposed mice.

A systematic review of non-coding RNA genes with differential expression profiles associated with autism spectrum disorders

AIMS

To identify differential expression of shorter non-coding RNA (ncRNA) genes associated with autism spectrum disorders (ASD).

BACKGROUND

ncRNA are functional molecules that derive from non-translated DNA sequence. The HUGO Gene Nomenclature Committee (HGNC) have approved ncRNA gene classes with alignment to the reference human genome. One subset is microRNA (miRNA), which are highly conserved, short RNA molecules that regulate gene expression by direct post-transcriptional repression of messenger RNA. Several miRNA genes are implicated in the development and regulation of the nervous system. Expression of miRNA genes in ASD cohorts have been examined by multiple research groups. Other shorter classes of ncRNA have been examined less. A comprehensive systematic review examining expression of shorter ncRNA gene classes in ASD is timely to inform the direction of research.

METHODS

We extracted data from studies examining ncRNA gene expression in ASD compared with non-ASD controls. We included studies on miRNA, piwi-interacting RNA (piRNA), small NF90 (ILF3) associated RNA (snaR), small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), transfer RNA (tRNA), vault RNA (vtRNA) and Y RNA. The following electronic databases were searched: Cochrane Library, EMBASE, PubMed, Web of Science, PsycINFO, ERIC, AMED and CINAHL for papers published from January 2000 to May 2022. Studies were screened by two independent investigators with a third resolving discrepancies. Data was extracted from eligible papers.

RESULTS

Forty-eight eligible studies were included in our systematic review with the majority examining miRNA gene expression alone. Sixty-four miRNA genes had differential expression in ASD compared to controls as reported in two or more studies, but often in opposing directions. Four miRNA genes had differential expression in the same direction in the same tissue type in at least 3 separate studies. Increased expression was reported in miR-106b-5p, miR-155-5p and miR-146a-5p in blood, post-mortem brain, and across several tissue types, respectively. Decreased expression was reported in miR-328-3p in bloods samples. Seven studies examined differential expression from other classes of ncRNA, including piRNA, snRNA, snoRNA and Y RNA. No individual ncRNA genes were reported in more than one study. Six studies reported differentially expressed snoRNA genes in ASD. A meta-analysis was not possible because of inconsistent methodologies, disparate tissue types examined, and varying forms of data presented.

CONCLUSION

There is limited but promising evidence associating the expression of certain miRNA genes and ASD, although the studies are of variable methodological quality and the results are largely inconsistent. There is emerging evidence associating differential expression of snoRNA genes in ASD. It is not currently possible to say whether the reports of differential expression in ncRNA may relate to ASD aetiology, a response to shared environmental factors linked to ASD such as sleep and nutrition, other molecular functions, human diversity, or chance findings. To improve our understanding of any potential association, we recommend improved and standardised methodologies and reporting of raw data. Further high-quality research is required to shine a light on possible associations, which may yet yield important information.

Proteomic analysis of sialoliths from calcified, lipid and mixed groups as a source of potential biomarkers of deposit formation in the salivary glands

Salivary stones, also known as sialoliths, are formed in a pathological situation in the salivary glands. So far, neither the mechanism of their formation nor the factors predisposing to their formation are known despite several hypotheses. While they do not directly threaten human life, they significantly deteriorate the patient's quality of life. Although this is not a typical research material, attempts are made to apply various analytical tools to characterise sialoliths and search for the biomarkers in their proteomes. In this work, we used mass spectrometry and SWATH-MS qualitative and quantitative analysis to investigate the composition and select proteins that may contribute to solid deposits in the salivary glands. Twenty sialoliths, previously characterized spectroscopically and divided into the following groups: calcified (CAL), lipid (LIP) and mixed (MIX), were used for the study. Proteins unique for each of the groups were found, including: for the CAL group among them, e.g. proteins from the S100 group (S100 A8/A12 and P), mucin 7 (MUC7), keratins (KRT1/2/4/5/13), elastase (ELANE) or stomatin (STOM); proteins for the LIP group-transthyretin (TTR), lactotransferrin (LTF), matrix Gla protein (MPG), submandibular gland androgen-regulated protein 3 (SMR3A); mixed stones had the fewest unique proteins. Bacterial proteins present in sialoliths have also been identified. The analysis of the results indicates the possible role of bacterial infections, disturbances in calcium metabolism and neutrophil extracellular traps (NETs) in the formation of sialoliths.

Distinct Cellular Origins and Differentiation Process Account for Distinct Oncogenic and Clinical Behaviors of Leiomyosarcomas

In leiomyosarcoma (LMS), a very aggressive disease, a relatively transcriptionally uniform subgroup of well-differentiated tumors has been described and is associated with poor survival. The question raised how differentiation and tumor progression, two apparently antagonist processes, coexist and allow tumor malignancy. We first identified the most transcriptionally homogeneous LMS subgroup in three independent cohorts, which we named 'hLMS'. The integration of multi-omics data and functional analysis suggests that hLMS originate from vascular smooth muscle cells and show that hLMS transcriptional program reflects both modulations of smooth muscle contraction activity controlled by MYOCD/SRF regulatory network and activation of the cell cycle activity controlled by E2F/RB1 pathway. We propose that the phenotypic plasticity of vascular smooth muscle cells coupled with MYOCD/SRF pathway amplification, essential for hLMS survival, concomitant with PTEN absence and RB1 alteration, could explain how hLMS balance this uncommon interplay between differentiation and aggressiveness.

X-chromosome inactivation patterns depend on age and tissue but not conception method in humans

Female somatic X-chromosome inactivation (XCI) balances the X-linked transcriptional dosages between the sexes, randomly silencing the maternal or paternal X chromosome in each cell of 46,XX females. Skewed XCI toward one parental X has been observed in association with ageing and in some female carriers of X-linked diseases. To address the problem of non-random XCI, we quantified the XCI skew in different biological samples of naturally conceived females of different age groups and girls conceived after in vitro fertilization (IVF). Generally, XCI skew differed between saliva, blood, and buccal swabs, while saliva and blood had the most similar XCI patterns in individual females. XCI skew increased with age in saliva, but not in other tissues. We showed no significant differences in the XCI patterns in tissues of naturally conceived and IVF females. The gene expression profile of the placenta and umbilical cord blood was determined depending on the XCI pattern. The increased XCI skewing in the placental tissue was associated with the differential expression of several genes out of 40 considered herein. Notably, skewed XCI patterns (> 80:20) were identified with significantly increased expression levels of four genes: CD44, KDM6A, PHLDA2, and ZRSR2. The differences in gene expression patterns between samples with random and non-random XCI may shed new light on factors contributing to the XCI pattern outcome and indicate new paths in future research on the phenomenon of XCI skewing.

Three Microbial Musketeers of the Seas: Shewanella baltica, Aliivibrio fischeri and Vibrio harveyi, and Their Adaptation to Different Salinity Probed by a Proteomic Approach

Osmotic changes are common challenges for marine microorganisms. Bacteria have developed numerous ways of dealing with this stress, including reprogramming of global cellular processes. However, specific molecular adaptation mechanisms to osmotic stress have mainly been investigated in terrestrial model bacteria. In this work, we aimed to elucidate the basis of adjustment to prolonged salinity challenges at the proteome level in marine bacteria. The objects of our studies were three representatives of bacteria inhabiting various marine environments, Shewanella baltica, Vibrio harveyi and Aliivibrio fischeri. The proteomic studies were performed with bacteria cultivated in increased and decreased salinity, followed by proteolytic digestion of samples which were then subjected to liquid chromatography with tandem mass spectrometry analysis. We show that bacteria adjust at all levels of their biological processes, from DNA topology through gene expression regulation and proteasome assembly, to transport and cellular metabolism. The finding that many similar adaptation strategies were observed for both low- and high-salinity conditions is particularly striking. The results show that adaptation to salinity challenge involves the accumulation of DNA-binding proteins and increased polyamine uptake. We hypothesize that their function is to coat and protect the nucleoid to counteract adverse changes in DNA topology due to ionic shifts.

Development of a fixed module repertoire for the analysis and interpretation of blood transcriptome data

As the capacity for generating large-scale molecular profiling data continues to grow, the ability to extract meaningful biological knowledge from it remains a limitation. Here, we describe the development of a new fixed repertoire of transcriptional modules, BloodGen3, that is designed to serve as a stable reusable framework for the analysis and interpretation of blood transcriptome data. The construction of this repertoire is based on co-clustering patterns observed across sixteen immunological and physiological states encompassing 985 blood transcriptome profiles. Interpretation is supported by customized resources, including module-level analysis workflows, fingerprint grid plot visualizations, interactive web applications and an extensive annotation framework comprising functional profiling reports and reference transcriptional profiles. Taken together, this well-characterized and well-supported transcriptional module repertoire can be employed for the interpretation and benchmarking of blood transcriptome profiles within and across patient cohorts. Blood transcriptome fingerprints for the 16 reference cohorts can be accessed interactively via:  https://drinchai.shinyapps.io/BloodGen3Module/ .

Virus-Host Interaction Gets Curiouser and Curiouser. PART II: Functional Transcriptomics of the E. coli DksA-Deficient Cell upon Phage P1vir Infection

The virus–host interaction requires a complex interplay between the phage strategy of reprogramming the host machinery to produce and release progeny virions, and the host defense against infection. Using RNA sequencing, we investigated the phage–host interaction to resolve the phenomenon of improved lytic development of P1vir phage in a DksA-deficient E. coli host. Expression of the ant1 and kilA P1vir genes in the wild-type host was the highest among all and most probably leads to phage virulence. Interestingly, in a DksA-deficient host, P1vir genes encoding lysozyme and holin are downregulated, while antiholins are upregulated. Gene expression of RepA, a protein necessary for replication initiating at the phage oriR region, is increased in the dksA mutant; this is also true for phage genes responsible for viral morphogenesis and architecture. Still, it seems that P1vir is taking control of the bacterial protein, sugar, and lipid metabolism in both, the wild type and dksA⁻ hosts. Generally, bacterial hosts are reacting by activating their SOS response or upregulating the heat shock proteins. However, only DksA-deficient cells upregulate their sulfur metabolism and downregulate proteolysis upon P1vir infection. We conclude that P1vir development is enhanced in the dksA mutant due to several improvements, including replication and virion assembly, as well as a less efficient lysis.

Activation of a neural stem cell transcriptional program in parenchymal astrocytes

Adult neural stem cells, located in discrete brain regions, generate new neurons throughout life. These stem cells are specialized astrocytes, but astrocytes in other brain regions do not generate neurons under physiological conditions. After stroke, however, striatal astrocytes undergo neurogenesis in mice, triggered by decreased Notch signaling. We used single-cell RNA sequencing to characterize neurogenesis by Notch-depleted striatal astrocytes in vivo. Striatal astrocytes were located upstream of neural stem cells in the neuronal lineage. As astrocytes initiated neurogenesis, they became transcriptionally very similar to subventricular zone stem cells, progressing through a near-identical neurogenic program. Surprisingly, in the non-neurogenic cortex, Notch-depleted astrocytes also initiated neurogenesis. Yet, these cortical astrocytes, and many striatal ones, stalled before entering transit-amplifying divisions. Infusion of epidermal growth factor enabled stalled striatal astrocytes to resume neurogenesis. We conclude that parenchymal astrocytes are latent neural stem cells and that targeted interventions can guide them through their neuronal differentiation.

A modular framework for the development of targeted Covid-19 blood transcript profiling panels

BACKGROUND

Covid-19 morbidity and mortality are associated with a dysregulated immune response. Tools are needed to enhance existing immune profiling capabilities in affected patients. Here we aimed to develop an approach to support the design of targeted blood transcriptome panels for profiling the immune response to SARS-CoV-2 infection.

METHODS

We designed a pool of candidates based on a pre-existing and well-characterized repertoire of blood transcriptional modules. Available Covid-19 blood transcriptome data was also used to guide this process. Further selection steps relied on expert curation. Additionally, we developed several custom web applications to support the evaluation of candidates.

RESULTS

As a proof of principle, we designed three targeted blood transcript panels, each with a different translational connotation: immunological relevance, therapeutic development relevance and SARS biology relevance.

CONCLUSION

Altogether the work presented here may contribute to the future expansion of immune profiling capabilities via targeted profiling of blood transcript abundance in Covid-19 patients.