Incomplete transcripts dominate the Mycobacterium tuberculosis transcriptome

Mycobacterium tuberculosis (Mtb) is a bacterial pathogen that causes tuberculosis (TB), an infectious disease that is responsible for major health and economic costs worldwide1. Mtb encounters diverse environments during its life cycle and responds to these changes largely by reprogramming its transcriptional output2. However, the mechanisms of Mtb transcription and how they are regulated remain poorly understood. Here we use a sequencing method that simultaneously determines both termini of individual RNA molecules in bacterial cells3 to profile the Mtb transcriptome at high resolution. Unexpectedly, we find that most Mtb transcripts are incomplete, with their 5' ends aligned at transcription start sites and 3' ends located 200-500 nucleotides downstream. We show that these short RNAs are mainly associated with paused RNA polymerases (RNAPs) rather than being products of premature termination. We further show that the high propensity of Mtb RNAP to pause early in transcription relies on the binding of the σ-factor. Finally, we show that a translating ribosome promotes transcription elongation, revealing a potential role for transcription-translation coupling in controlling Mtb gene expression. In sum, our findings depict a mycobacterial transcriptome that prominently features incomplete transcripts resulting from RNAP pausing. We propose that the pausing phase constitutes an important transcriptional checkpoint in Mtb that allows the bacterium to adapt to environmental changes and could be exploited for TB therapeutics.

Spatially organized cellular communities form the developing human heart

The heart, which is the first organ to develop, is highly dependent on its form to function1,2. However, how diverse cardiac cell types spatially coordinate to create the complex morphological structures that are crucial for heart function remains unclear. Here we integrated single-cell RNA-sequencing with high-resolution multiplexed error-robust fluorescence in situ hybridization to resolve the identity of the cardiac cell types that develop the human heart. This approach also provided a spatial mapping of individual cells that enables illumination of their organization into cellular communities that form distinct cardiac structures. We discovered that many of these cardiac cell types further specified into subpopulations exclusive to specific communities, which support their specialization according to the cellular ecosystem and anatomical region. In particular, ventricular cardiomyocyte subpopulations displayed an unexpected complex laminar organization across the ventricular wall and formed, with other cell subpopulations, several cellular communities. Interrogating cell-cell interactions within these communities using in vivo conditional genetic mouse models and in vitro human pluripotent stem cell systems revealed multicellular signalling pathways that orchestrate the spatial organization of cardiac cell subpopulations during ventricular wall morphogenesis. These detailed findings into the cellular social interactions and specialization of cardiac cell types constructing and remodelling the human heart offer new insights into structural heart diseases and the engineering of complex multicellular tissues for human heart repair.

Dynamic dysregulation of transcriptomic networks in brainstem autonomic nuclei during hypertension development in the female spontaneously hypertensive rat

Neurogenic hypertension stems from an imbalance in autonomic function that shifts the central cardiovascular control circuits toward a state of dysfunction. Using the female spontaneously hypertensive rat and the normotensive Wistar-Kyoto rat model, we compared the transcriptomic changes in three autonomic nuclei in the brainstem, nucleus of the solitary tract (NTS), caudal ventrolateral medulla, and rostral ventrolateral medulla (RVLM) in a time series at 8, 10, 12, 16, and 24 wk of age, spanning the prehypertensive stage through extended chronic hypertension. RNA-sequencing data were analyzed using an unbiased, dynamic pattern-based approach that uncovered dominant and several subtle differential gene regulatory signatures. Our results showed a persistent dysregulation across all three autonomic nuclei regardless of the stage of hypertension development as well as a cascade of transient dysregulation beginning in the RVLM at the prehypertensive stage that shifts toward the NTS at the hypertension onset. Genes that were persistently dysregulated were heavily enriched for immunological processes such as antigen processing and presentation, the adaptive immune response, and the complement system. Genes with transient dysregulation were also largely region-specific and were annotated for processes that influence neuronal excitability such as synaptic vesicle release, neurotransmitter transport, and an array of neuropeptides and ion channels. Our results demonstrate that neurogenic hypertension is characterized by brainstem region-specific transcriptomic changes that are highly dynamic with significant gene regulatory changes occurring at the hypertension onset as a key time window for dysregulation of homeostatic processes across the autonomic control circuits.NEW & NOTEWORTHY Hypertension is a major disease and is the primary risk factor for cardiovascular complications and stroke. The gene expression changes in the central nervous system circuits driving hypertension are understudied. Here, we show that coordinated and region-specific gene expression changes occur in the brainstem autonomic circuits over time during the development of a high blood pressure phenotype in a rat model of human essential hypertension.

Lipidomic and transcriptomic analysis of the increase in eicosapentaenoic acid under cobalamin deficiency of Schizochytrium sp

Schizochytrium sp. is a heterotrophic microorganism capable of accumulating polyunsaturated fatty acids and has achieved industrial production of docosahexaenoic acid (DHA). It also has the potential for eicosapentaenoic acid (EPA) production. In this study, it was found that the cell growth, lipid synthesis and fatty acid composition of Schizochytrium sp. were significantly affected by the level of cobalamin in the medium, especially with regard to the content of EPA in the fatty acids. The content of EPA in the fatty acids increased 17.91 times, reaching 12.00%, but cell growth and lipid synthesis were significantly inhibited under cobalamin deficiency. The response mechanism for this phenomenon was revealed through combined lipidomic and transcriptomic analysis. Although cell growth was inhibited under cobalamin deficiency, the genes encoding key enzymes in central carbon metabolism were still up-regulated to provide precursors (Acetyl-CoA) and reducing power (NADPH) for the synthesis and accumulation of fatty acids. Moreover, the main lipid subclasses observed during cobalamin deficiency were glycerolipids (including glycerophospholipids), with EPA primarily distributed in them. The genes involved in the biosynthesis of these lipid subclasses were significantly up-regulated, such as the key enzymes in the Kennedy pathway for the synthesis of triglycerides. Thus, this study provided insights into the specific response of Schizochytrium sp. to cobalamin deficiency and identified a subset of new genes that can be engineered for modification.

Cuticular wax biosynthesis in blueberries (Vaccinium corymbosum L.): Transcript and metabolite changes during ripening and storage affect key fruit quality traits

In fruits, cuticular waxes affect fruit quality traits such as surface color at harvest and water loss during postharvest storage. This study investigated the transcriptional regulation of cuticular wax deposition in northern highbush blueberries (Vaccinium corymbosum L.) in relation to fruit water loss and surface color during ripening and postharvest storage, as well as the effects of abscisic acid (ABA)-mediated changes in cuticular wax deposition on these fruit quality traits. Total cuticular wax content (μg∙cm-2) decreased during fruit ripening and increased during postharvest storage. Transcriptome analysis revealed a transcript network for cuticular wax deposition in blueberries. Particularly, five OSC-Likes were identified as putative genes for triterpene alcohol production, with OSC-Like1 and OSC-Like2 encoding mixed amyrin synthases, OSC-Like3 encoding a lupeol synthase, and OSC-Like4 and OSC-Like5 encoding cycloartenol synthases. The expression of three CYP716A-like genes correlated to the accumulation of two triterpene acids oleanolic acid and ursolic acid, the major wax compounds in blueberries. Exogenous ABA application induced the expression of triterpenoid biosynthetic genes and the accumulation of β-amyrin and oleanolic acid, as well as increased the ratio of oleanolic acid to ursolic acid. These changes were associated with reduced fruit water loss. The content of β-diketones was also increased by ABA application, and this increase was associated with increased fruit lightness (measured as L* using CIELAB Color Space by a colorimeter). This study provided key insights on the molecular basis of cuticular wax deposition and its implications on fruit quality traits in blueberries.

Genetic basis and expression of ventral colour in polymorphic common lizards

Colour is an important visual cue that can correlate with sex, behaviour, life history or ecological strategies, and has evolved divergently and convergently across animal lineages. Its genetic basis in non-model organisms is rarely known, but such information is vital for determining the drivers and mechanisms of colour evolution. Leveraging genetic admixture in a rare contact zone between oviparous and viviparous common lizards (Zootoca vivipara), we show that females (N = 558) of the two otherwise morphologically indistinguishable reproductive modes differ in their ventral colouration (from pale to vibrant yellow) and intensity of melanic patterning. We find no association between female colouration and reproductive investment, and no evidence for selection on colour. Using a combination of genetic mapping and transcriptomic evidence, we identified two candidate genes associated with ventral colour differentiation, DGAT2 and PMEL. These are genes known to be involved in carotenoid metabolism and melanin synthesis respectively. Ventral melanic spots were associated with two genomic regions, including a SNP close to protein tyrosine phosphatase (PTP) genes. Using genome re-sequencing data, our results show that fixed coding mutations in the candidate genes cannot account for differences in colouration. Taken together, our findings show that the evolution of ventral colouration and its associations across common lizard lineages is variable. A potential genetic mechanism explaining the flexibility of ventral colouration may be that colouration in common lizards, but also across squamates, is predominantly driven by regulatory genetic variation.

A modified diffusion-weighted magnetic resonance imaging-based model from the radiologist's perspective: improved performance in determining the surgical resectability of advanced high-grade serous ovarian cancer

BACKGROUND

Complete resection of all visible lesions during primary debulking surgery is associated with the most favorable prognosis in patients with advanced high-grade serous ovarian cancer. An accurate preoperative assessment of resectability is pivotal for tailored management.

OBJECTIVE

This study aimed to assess the potential value of a modified model that integrates the original 8 radiologic criteria of the Memorial Sloan Kettering Cancer Center model with imaging features of the subcapsular or diaphragm and mesenteric lesions depicted on diffusion-weighted magnetic resonance imaging and growth patterns of all lesions for predicting the resectability of advanced high-grade serous ovarian cancer.

STUDY DESIGN

This study included 184 patients with high-grade serous ovarian cancer who underwent preoperative diffusion-weighted magnetic resonance imaging between December 2018 and May 2023 at 2 medical centers. The patient cohort was divided into 3 subsets, namely a study cohort (n=100), an internal validation cohort (n=46), and an external validation cohort (n=38). Preoperative radiologic evaluations were independently conducted by 2 radiologists using both the Memorial Sloan Kettering Cancer Center model and the modified diffusion-weighted magnetic resonance imaging-based model. The morphologic characteristics of the ovarian tumors depicted on magnetic resonance imaging were assessed as either mass-like or infiltrative, and transcriptomic analysis of the primary tumor samples was performed. Univariate and multivariate statistical analyses were performed.

RESULTS

In the study cohort, both the scores derived using the Memorial Sloan Kettering Cancer Center (intraclass correlation coefficients of 0.980 and 0.959, respectively; both P<.001) and modified diffusion-weighted magnetic resonance imaging-based models (intraclass correlation coefficients of 0.962 and 0.940, respectively; both P<.001) demonstrated excellent intra- and interobserver agreement. The Memorial Sloan Kettering Cancer Center model (odds ratio, 1.825; 95% confidence interval, 1.390-2.395; P<.001) and the modified diffusion-weighted magnetic resonance imaging-based model (odds ratio, 1.776; 95% confidence interval, 1.410-2.238; P<.001) independently predicted surgical resectability. The modified diffusion-weighted magnetic resonance imaging-based model demonstrated improved predictive performance with an area under the curve of 0.867 in the study cohort and 0.806 and 0.913 in the internal and external validation cohorts, respectively. Using the modified diffusion-weighted magnetic resonance imaging-based model, patients with scores of 0 to 2, 3 to 4, 5 to 6, 7 to 10, and ≥11 achieved complete tumor debulking rates of 90.3%, 66.7%, 53.3%, 11.8%, and 0%, respectively. Most patients with incomplete tumor debulking had infiltrative tumors, and both the Memorial Sloan Kettering Cancer Center and the modified diffusion-weighted magnetic resonance imaging-based models yielded higher scores. The molecular differences between the 2 morphologic subtypes were identified.

CONCLUSION

When compared with the Memorial Sloan Kettering Cancer Center model, the modified diffusion-weighted magnetic resonance imaging-based model demonstrated enhanced accuracy in the preoperative prediction of resectability for advanced high-grade serous ovarian cancer. Patients with scores of 0 to 6 were eligible for primary debulking surgery.

Transcriptomic changes associated with oral immunotherapy for food allergy

This review summarizes recent advances in characterizing the transcriptional pathways associated with outcomes following Oral Immunotherapy. Recent technological advances including single-cell sequencing are transforming the ways in which the transcriptional landscape is understood. The application of these technologies is still in its infancy in food allergy but here we summarize current understanding of gene expression changes following oral immunotherapy for food allergy and specific signatures underpinning the different clinical outcomes of desensitization and remission (sustained unresponsiveness). T helper 2A cells have been identified as a cell type which correlates with disease activity and is modified by treatment. Molecular features at study entry may differentiate individuals who achieve more positive outcomes during OIT. Recent findings point to T cell anergy and Type 1 interferon pathways as potential mechanisms supporting redirection of the allergen-specific immune response away from allergy towards remission. Despite these developments in our understanding of immune mechanisms following OIT, there are still significant gaps. Additional studies examining immune signatures associated with long term and well-defined clinical outcomes are required to gain a more complete understanding of the pathways leading to remission of allergy, in order to optimize treatments and gain improved outcomes for patients.

Multi-omic analysis of precocious puberty girls: pathway changes and metabolite validation

Objective

Precocious puberty (PP) is a prevalent endocrine disorder affecting the physical and mental wellbeing of children. Identifying the triggering factors of PP has become a central issue. This study seeks to investigate the metabolomic and transcriptomic alterations in PP.

Material and methods

First, 37 school-aged girls diagnosed with PP and 25 age-matched prepubertal control girls were recruited, and the fecal samples were collected for non-targeted metabolomic analysis to screen for differentially expressed metabolites (DEMs). Subsequently, an animal model of PP was constructed by danazol administration to neonatal female rats, and both fecal non-targeted metabolomics and serum next-generation transcriptomic sequencing were performed to screen DEMs and differentially expressed genes (DEGs) in PP. Moreover, the DEM co-existing in clinical and animal models was administrated to PP rats to explore the role of the target metabolite in PP.

Results

A total of 24 DEMs in PP clinical samples and 180 DEMs and 425 DEGs in PP animal samples were identified. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that these DEMs and DEGs were enriched in disease-associated pathways, including fatty acid synthesis, glycerolipid metabolism, pyrimidine metabolism, steroid hormone biosynthesis, progesterone-mediated oocyte maturation, and gonadotropin-releasing hormone (GnRH) signaling pathway, forming a tight DEM-DEG pathway regulatory network. Further DEM validation demonstrated that thymine supplementation delayed the opening of the vagina and development of PP in model rats.

Conclusion

This study reveals that the metabolomic and transcriptomic changes, along with enriched pathways, are implicated in PP based on clinical and animal analyses. The findings may provide new strategies and research avenues for PP treatment.

New insights into the role of immunity and inflammation in diabetic kidney disease in the omics era

Diabetic kidney disease (DKD) is becoming the leading cause of chronic kidney disease, especially in the industrialized world. Despite mounting evidence has demonstrated that immunity and inflammation are highly involved in the pathogenesis and progression of DKD, the underlying mechanisms remain incompletely understood. Substantial molecules, signaling pathways, and cell types participate in DKD inflammation, by integrating into a complex regulatory network. Most of the studies have focused on individual components, without presenting their importance in the global or system-based processes, which largely hinders clinical translation. Besides, conventional technologies failed to monitor the different behaviors of resident renal cells and immune cells, making it difficult to understand their contributions to inflammation in DKD. Recently, the advancement of omics technologies including genomics, epigenomics, transcriptomics, proteomics, and metabolomics has revolutionized biomedical research, which allows an unbiased global analysis of changes in DNA, RNA, proteins, and metabolites in disease settings, even at single-cell and spatial resolutions. They help us to identify critical regulators of inflammation processes and provide an overview of cell heterogeneity in DKD. This review aims to summarize the application of multiple omics in the field of DKD and emphasize the latest evidence on the interplay of inflammation and DKD revealed by these technologies, which will provide new insights into the role of inflammation in the pathogenesis of DKD and lead to the development of novel therapeutic approaches and diagnostic biomarkers.

Ror homolog nhr-23 is essential for both developmental clock and circadian clock in C. elegans

Animals have internal clocks that generate biological rhythms. In mammals, clock genes such as Period form the circadian clock to generate approximately 24-h biological rhythms. In C. elegans, the clock gene homologs constitute the "developmental clock", which has an 8-h period during larval development to determine the timing of molting. Thus, the ancestral circadian clock has been believed to evolve into the oscillator with a shorter period in C. elegans. However, circadian rhythms have also been observed in adult C. elegans, albeit relatively weak. This prompts the question: if the clock gene homologs drive the developmental rhythm with 8-h period, which genes generate the circadian rhythms in C. elegans? In this study, we discovered that nhr-23, a homolog of the mammalian circadian clock gene Ror, is essential for circadian transcriptional rhythms in adult C. elegans. Interestingly, nhr-23 was also known to be essential for the molting clock. The bilaterian ancestral circadian clock genes might have evolved to function over multiple periods depending on developmental contexts rather than a single 8-h period in C. elegans.

Transcriptomic Responses to Koi Herpesvirus in Isolated Blood Leukocytes from Infected Common Carp

Koi herpesvirus (KHV, CyHV-3) causes severe economic losses in carp farms. Its eradication is challenging due to the establishment of latency in blood leukocytes and other tissues. To understand the molecular mechanisms leading to KHV infection in leukocytes, common carp were bath-exposed to KHV at 17 °C. After confirming the presence of viral transcripts in blood leukocytes at ten days post infection, RNA-Seq was performed on peripheral blood leukocytes on the Illumina NovaSeq. KHV infection triggered a robust immune response mediated by pattern recognition receptors, mainly toll-like receptors (tlr2, tlr5, tlr7, and tlr13), urokinase plasminogen activator surface receptor-like, galectin proteins, and lipid mediators such as leukotriene B4 receptor 1. Enriched pathways showed increased mitochondria oxidative phosphorylation and the activation of signalling pathways such as mitogen-activated protein kinases (MAPKs) and vascular endothelial growth factor (VEGF). KHV-infected leukocytes showed low production of reactive oxygen species (ROS) and glutathione metabolism, high iron export and phagocytosis activity, and low autophagy. Macrophage polarization was deduced from the up-regulation of genes such as arginase non-hepatic 1-like, macrophage mannose receptor-1, crem, il-10, and il-13 receptors, while markers for cytotoxic T cells were observed to be down-regulated. Further work is required to characterise these leukocyte subsets and the molecular events leading to KHV latency in blood leukocytes.

Refining the transcriptional landscapes for distinct clades of virulent phages infecting Pseudomonas aeruginosa

The introduction of high-throughput sequencing has resulted in a surge of available bacteriophage genomes, unveiling their tremendous genomic diversity. However, our current understanding of the complex transcriptional mechanisms that dictate their gene expression during infection is limited to a handful of model phages. Here, we applied ONT-cappable-seq to reveal the transcriptional architecture of six different clades of virulent phages infecting Pseudomonas aeruginosa. This long-read microbial transcriptomics approach is tailored to globally map transcription start and termination sites, transcription units, and putative RNA-based regulators on dense phage genomes. Specifically, the full-length transcriptomes of LUZ19, LUZ24, 14-1, YuA, PAK_P3, and giant phage phiKZ during early, middle, and late infection were collectively charted. Beyond pinpointing traditional promoter and terminator elements and transcription units, these transcriptional profiles provide insights in transcriptional attenuation and splicing events and allow straightforward validation of Group I intron activity. In addition, ONT-cappable-seq data can guide genome-wide discovery of novel regulatory element candidates, including noncoding RNAs and riboswitches. This work substantially expands the number of annotated phage-encoded transcriptional elements identified to date, shedding light on the intricate and diverse gene expression regulation mechanisms in Pseudomonas phages, which can ultimately be sourced as tools for biotechnological applications in phage and bacterial engineering.

A multiomics dataset for the study of RNA modifications in human macrophage differentiation and polarisation

RNA modifications have emerged as central regulators of gene expression programs. Amongst RNA modifications are N6-methyladenosine (m6A) and RNA 5-hydroxymethylcytosine (5hmC). While m6A is established as a versatile regulator of RNA metabolism, the functions of RNA 5hmC are unclear. Despite some evidence linking RNA modifications to immunity, their implications in gene expression control in macrophage development and functions remain unclear. Here we present a multi-omics dataset capturing different layers of the gene expression programs driving macrophage differentiation and polarisation. We obtained mRNA-Seq, m6A-IP-Seq, 5hmC-IP-Seq, Polyribo-Seq and LC-MS/MS data from monocytes and resting-, pro- and anti-inflammatory-like macrophages. We present technical validation showing high quality and correlation between samples for all datasets, and evidence of biological consistency of modelled macrophages at the transcriptomic, epitranscriptomic, translational and proteomic levels. This multi-omics dataset provides a resource for the study of RNA m6A and 5hmC in the context of macrophage biology and spans the gene expression process from transcripts to proteins.

Hybridization and gene expression: Beyond differentially expressed genes

Gene expression has a key role in reproductive isolation, and studies of hybrid gene expression have identified mechanisms causing hybrid sterility. Here, we review the evidence for altered gene expression following hybridization and outline the mechanisms shown to contribute to altered gene expression in hybrids. Transgressive gene expression, transcending that of both parental species, is pervasive in early generation sterile hybrids, but also frequently observed in viable, fertile hybrids. We highlight studies showing that hybridization can result in transgressive gene expression, also in established hybrid lineages or species. Such extreme patterns of gene expression in stabilized hybrid taxa suggest that altered hybrid gene expression may result in hybridization-derived evolutionary novelty. We also conclude that while patterns of misexpression in hybrids are well documented, the understanding of the mechanisms causing misexpression is lagging. We argue that jointly assessing differences in cell composition and cell-specific changes in gene expression in hybrids, in addition to assessing changes in chromatin and methylation, will significantly advance our understanding of the basis of altered gene expression. Moreover, uncovering to what extent evolution of gene expression results in altered expression for individual genes, or entire networks of genes, will advance our understanding of how selection moulds gene expression. Finally, we argue that jointly studying the dual roles of altered hybrid gene expression, serving both as a mechanism for reproductive isolation and as a substrate for hybrid ecological adaptation, will lead to significant advances in our understanding of the evolution of gene expression.

Systems immunology of transcriptional responses to viral infection identifies conserved antiviral pathways across macaques and humans

Viral pandemics and epidemics pose a significant global threat. While macaque models of viral disease are routinely used, it remains unclear how conserved antiviral responses are between macaques and humans. Therefore, we conducted a cross-species analysis of transcriptomic data from over 6,088 blood samples from macaques and humans infected with one of 31 viruses. Our findings demonstrate that irrespective of primate or viral species, there are conserved antiviral responses that are consistent across infection phase (acute, chronic, or latent) and viral genome type (DNA or RNA viruses). Leveraging longitudinal data from experimental challenges, we identify virus-specific response kinetics such as host responses to Coronaviridae and Orthomyxoviridae infections peaking 1-3 days earlier than responses to Filoviridae and Arenaviridae viral infections. Our results underscore macaque studies as a powerful tool for understanding viral pathogenesis and immune responses that translate to humans, with implications for viral therapeutic development and pandemic preparedness.

Statistical method scDEED for detecting dubious 2D single-cell embeddings and optimizing t-SNE and UMAP hyperparameters

Two-dimensional (2D) embedding methods are crucial for single-cell data visualization. Popular methods such as t-distributed stochastic neighbor embedding (t-SNE) and uniform manifold approximation and projection (UMAP) are commonly used for visualizing cell clusters; however, it is well known that t-SNE and UMAP's 2D embeddings might not reliably inform the similarities among cell clusters. Motivated by this challenge, we present a statistical method, scDEED, for detecting dubious cell embeddings output by a 2D-embedding method. By calculating a reliability score for every cell embedding based on the similarity between the cell's 2D-embedding neighbors and pre-embedding neighbors, scDEED identifies the cell embeddings with low reliability scores as dubious and those with high reliability scores as trustworthy. Moreover, by minimizing the number of dubious cell embeddings, scDEED provides intuitive guidance for optimizing the hyperparameters of an embedding method. We show the effectiveness of scDEED on multiple datasets for detecting dubious cell embeddings and optimizing the hyperparameters of t-SNE and UMAP.

The Future of Targeted Therapy for Leiomyosarcoma

Leiomyosarcoma (LMS) is an aggressive subtype of soft tissue sarcoma that arises from smooth muscle cells, most commonly in the uterus and retroperitoneum. LMS is a heterogeneous disease with diverse clinical and molecular characteristics that have yet to be fully understood. Molecular profiling has uncovered possible targets amenable to treatment, though this has yet to translate into approved targeted therapies in LMS. This review will explore historic and recent findings from molecular profiling, highlight promising avenues of current investigation, and suggest possible future strategies to move toward the goal of molecularly matched treatment of LMS. We focus on targeting the DNA damage response, the macrophage-rich micro-environment, the PI3K/mTOR pathway, epigenetic regulators, and telomere biology.

Impaired bone morphogenetic protein (BMP) signaling pathways disrupt decidualization in endometriosis

Endometriosis is linked to increased infertility and pregnancy complications due to defective endometrial decidualization. We hypothesized that identification of altered signaling pathways during decidualization could identify the underlying cause of infertility and pregnancy complications. Our study reveals that transforming growth factor β (TGFβ) pathways are impaired in the endometrium of individuals with endometriosis, leading to defective decidualization. Through detailed transcriptomic analyses, we discovered abnormalities in TGFβ signaling pathways and key regulators, such as SMAD4, in the endometrium of affected individuals. We also observed compromised activity of bone morphogenetic proteins (BMP), a subset of the TGFβ family, that control endometrial receptivity. Using 3-dimensional models of endometrial stromal and epithelial assembloids, we showed that exogenous BMP2 improved decidual marker expression in individuals with endometriosis. Our findings reveal dysfunction of BMP/SMAD signaling in the endometrium of individuals with endometriosis, explaining decidualization defects and subsequent pregnancy complications in these individuals.

Tetraspanin CD82 Correlates with and May Regulate S100A7 Expression in Oral Cancer

Many metastatic cancers with poor prognoses correlate to downregulated CD82, but exceptions exist. Understanding the context of this correlation is essential to CD82 as a prognostic biomarker and therapeutic target. Oral squamous cell carcinoma (OSCC) constitutes over 90% of oral cancer. We aimed to uncover the function and mechanism of CD82 in OSCC. We investigated CD82 in human OSCC cell lines, tissues, and healthy controls using the CRISPR-Cas9 gene knockout, transcriptomics, proteomics, etc. CD82 expression is elevated in CAL 27 cells. Knockout CD82 altered over 300 genes and proteins and inhibited cell migration. Furthermore, CD82 expression correlates with S100 proteins in CAL 27, CD82KO, SCC-25, and S-G cells and some OSCC tissues. The 37-50 kDa CD82 protein in CAL 27 cells is upregulated, glycosylated, and truncated. CD82 correlates with S100 proteins and may regulate their expression and cell migration. The truncated CD82 explains the invasive metastasis and poor outcome of the CAL 27 donor. OSCC with upregulated truncated CD82 and S100A7 may represent a distinct subtype with a poor prognosis. Differing alternatives from wild-type CD82 may elucidate the contradictory functions and pave the way for CD82 as a prognostic biomarker and therapeutic target.

Corrigendum: Advances and potential of omics studies for understanding the development of food allergy

[This corrects the article DOI: 10.3389/falgy.2023.1149008.].

A 3D Bioprinted Human Neurovascular Unit Model of Glioblastoma Tumor Growth

A 3D bioprinted neurovascular unit (NVU) model is developed to study glioblastoma (GBM) tumor growth in a brain-like microenvironment. The NVU model includes human primary astrocytes, pericytes and brain microvascular endothelial cells, and patient-derived glioblastoma cells (JHH-520) are used for this study. Fluorescence reporters are used with confocal high content imaging to quantitate real-time microvascular network formation and tumor growth. Extensive validation of the NVU-GBM model includes immunostaining for brain relevant cellular markers and extracellular matrix components; single cell RNA sequencing (scRNAseq) to establish physiologically relevant transcriptomics changes; and secretion of NVU and GBM-relevant cytokines. The scRNAseq reveals changes in gene expression and cytokines secretion associated with wound healing/angiogenesis, including the appearance of an endothelial mesenchymal transition cell population. The NVU-GBM model is used to test 18 chemotherapeutics and anti-cancer drugs to assess the pharmacological relevance of the model and robustness for high throughput screening.

Multi-Omics Analysis Reveals That Anthocyanin Degradation and Phytohormone Changes Regulate Red Color Fading in Rapeseed (Brassica napus L.) Petals

Flower color is an important trait for the ornamental value of colored rapeseed (Brassica napus L.), as the plant is becoming more popular. However, the color fading of red petals of rapeseed is a problem for its utilization. Unfortunately, the mechanism for the process of color fading in rapeseed is unknown. In the current study, a red flower line, Zhehuhong, was used as plant material to analyze the alterations in its morphological and physiological characteristics, including pigment and phytohormone content, 2 d before flowering (T1), at flowering (T2), and 2 d after flowering (T3). Further, metabolomics and transcriptomics analyses were also performed to reveal the molecular regulation of petal fading. The results show that epidermal cells changed from spherical and tightly arranged to totally collapsed from T1 to T3, according to both paraffin section and scanning electron microscope observation. The pH value and all pigment content except flavonoids decreased significantly during petal fading. The anthocyanin content was reduced by 60.3% at T3 compared to T1. The content of three phytohormones, 1-aminocyclopropanecarboxylic acid, melatonin, and salicylic acid, increased significantly by 2.2, 1.1, and 30.3 times, respectively, from T1 to T3. However, auxin, abscisic acid, and jasmonic acid content decreased from T1 to T3. The result of metabolomics analysis shows that the content of six detected anthocyanin components (cyanidin, peonidin, pelargonidin, delphinidin, petunidin, and malvidin) and their derivatives mainly exhibited a decreasing trend, which was in accordance with the trend of decreasing anthocyanin. Transcriptomics analysis showed downregulation of genes involved in flavonol, flavonoid, and anthocyanin biosynthesis. Furthermore, genes regulating anthocyanin biosynthesis were preferentially expressed at early stages, indicating that the degradation of anthocyanin is the main issue during color fading. The corresponding gene-encoding phytohormone biosynthesis and signaling, JASMONATE-ZIM-DOMAIN PROTEIN, was deactivated to repress anthocyanin biosynthesis, resulting in fading petal color. The results clearly suggest that anthocyanin degradation and phytohormone regulation play essential roles in petal color fading in rapeseed, which is a useful insight for the breeding of colored rapeseed.

Establishment and characterization of turtle liver organoids provides a potential model to decode their unique adaptations

Painted turtles are remarkable for their freeze tolerance and supercooling ability along with their associated resilience to hypoxia/anoxia and oxidative stress, rendering them an ideal biomedical model for hypoxia-induced injuries (including strokes), tissue cooling during surgeries, and organ cryopreservation. Yet, such research is hindered by their seasonal reproduction and slow maturation. Here we developed and characterized adult stem cell-derived turtle liver organoids (3D self-assembled in vitro structures) from painted, snapping, and spiny softshell turtles spanning ~175My of evolution, with a subset cryopreserved. This development is, to the best of our knowledge, a first for this vertebrate Order, and complements the only other non-avian reptile organoids from snake venom glands. Preliminary characterization, including morphological, transcriptomic, and proteomic analyses, revealed organoids enriched in cholangiocytes. Deriving organoids from distant turtles and life stages demonstrates that our techniques are broadly applicable to chelonians, permitting the development of functional genomic tools currently lacking in herpetological research. Such platform could potentially support studies including genome-to-phenome mapping, gene function, genome architecture, and adaptive responses to climate change, with implications for ecological, evolutionary, and biomedical research.

Mechanism of 2'-fucosyllactose degradation by human-associated Akkermansia

Among the first microorganisms to colonize the human gut of breastfed infants are bacteria capable of fermenting human milk oligosaccharides (HMOs). One of the most abundant HMOs, 2'-fucosyllactose (2'-FL), may specifically drive bacterial colonization of the intestine. Recently, differential growth has been observed across multiple species of Akkermansia on various HMOs including 2'-FL. In culture, we found growth of two species, A. muciniphila MucT and A. biwaensis CSUN-19,on HMOs corresponded to a decrease in the levels of 2'-FL and an increase in lactose, indicating that the first step in 2'-FL catabolism is the cleavage of fucose. Using phylogenetic analysis and transcriptional profiling, we found that the number and expression of fucosidase genes from two glycoside hydrolase (GH) families, GH29 and GH95, vary between these two species. During the mid-log phase of growth, the expression of several GH29 genes was increased by 2'-FL in both species, whereas the GH95 genes were induced only in A. muciniphila. We further show that one putative fucosidase and a β-galactosidase from A. biwaensis are involved in the breakdown of 2'-FL. Our findings indicate that the plasticity of GHs of human-associated Akkermansia sp. enables access to additional growth substrates present in HMOs, including 2'-FL. Our work highlights the potential for Akkermansia to influence the development of the gut microbiota early in life and expands the known metabolic capabilities of this important human symbiont.IMPORTANCEAkkermansia are mucin-degrading specialists widely distributed in the human population. Akkermansia biwaensis has recently been observed to have enhanced growth relative to other human-associated Akkermansia on multiple human milk oligosaccharides (HMOs). However, the mechanisms for enhanced growth are not understood. Here, we characterized the phylogenetic diversity and function of select genes involved in the growth of A. biwaensis on 2'-fucosyllactose (2'-FL), a dominant HMO. Specifically, we demonstrate that two genes in a genomic locus, a putative β-galactosidase and α-fucosidase, are likely responsible for the enhanced growth on 2'-FL. The functional characterization of A. biwaensis growth on 2'-FL delineates the significance of a single genomic locus that may facilitate enhanced colonization and functional activity of select Akkermansia early in life.

lncRNA-miRNA-mRNA network in kidney transcriptome of Labeo rohita under hypersaline environment

The present study describes the kidney transcriptome of Labeo rohita, a freshwater fish, exposed to gradually increased salinity concentrations (2, 4, 6 and 8ppt). A total of 10.25 Gbps data was generated, and a suite of bioinformatics tools, including FEELnc, CPC2 and BLASTn were employed for identification of long non-coding RNAs (lncRNAs) and micro RNAs (miRNAs). Our analysis revealed a total of 170, 118, 99, and 269 differentially expressed lncRNA and 120, 118, 99, and 124 differentially expressed miRNAs in 2, 4, 6 and 8 ppt treatment groups respectively. Two competing endogenous RNA (ceRNA) networks were constructed i.e. A* ceRNA network with up-regulated lncRNAs and mRNAs, down-regulated miRNAs; and B* ceRNA network vice versa. 2ppt group had 131 and 83 lncRNA-miRNA-mRNA pairs in A* and B* networks, respectively. 4ppt group featured 163 pairs in A* network and 191 in B* network, while the 6ppt had 103 and 105 pairs. 8ppt group included 192 and 174 pairs. These networks illuminate the intricate RNA interactions in freshwater fish to varying salinity conditions.

A First Omics Data Integration Approach in Hass Avocados to Evaluate Rootstock-Scion Interactions: From Aerial and Root Plant Growth to Fruit Development

Grafting, the careful selection of rootstocks and scions, has played a crucial role maintaining Chilean avocado fruit quality standards in a scenario in which climate change and drought-related issues have considerably decreased avocado fruit production in the last fifteen years. The historical use of seedling rootstocks in Chile has experienced a recent shift towards clonal rootstocks, driven by the potential to produce more consistent and predictable crops. This research aims to compare Hass avocado plants grafted on Mexicola seedling and Dusa® clonal rootstocks in a soilless and protected system using (i) a differential expression analysis of root and leaf samples and (ii) a fruit transcriptomic and metabolomic integration analysis to improve our understanding of rootstock-scion interaction and its impact on avocado tree performance and fruit quality. The results demonstrated that no significant transcriptomic and metabolomic differences were identified at fruit level in the ready-to-eat (RTE) stage for Hass avocado fruit from both rootstocks. However, Hass avocados grafted on the clonal rootstock showed greater aerial growth and slightly increased fruit size than the seedling rootstock due to the enrichment of cell wall-remodeling genes as revealed in leaves and fruit at harvest stage.

The Biology and Genomics of Human Hair Follicles: A Focus on Androgenetic Alopecia

Androgenetic alopecia is a highly prevalent condition mainly affecting men. This complex trait is related to aging and genetics; however, multiple other factors, for example, lifestyle, are also involved. Despite its prevalence, the underlying biology of androgenetic alopecia remains elusive, and thus advances in its treatment have been hindered. Herein, we review the functional anatomy of hair follicles and the cell signaling events that play a role in follicle cycling. We also discuss the pathology of androgenetic alopecia and the known molecular mechanisms underlying this condition. Additionally, we describe studies comparing the transcriptional differences in hair follicles between balding and non-balding scalp regions. Given the genetic contribution, we also discuss the most significant risk variants found to be associated with androgenetic alopecia. A more comprehensive understanding of this pathology may be generated through using multi-omics approaches.

Sex differences in physiological response to increased neuronal excitability in a knockin mouse model of pediatric epilepsy

BACKGROUND

Epilepsy is a common neurological disease; however, few if any of the currently marketed antiseizure medications prevent or cure epilepsy. Discovery of pathological processes in the early stages of epileptogenesis has been challenging given the common use of preclinical models that induce seizures in physiologically normal animals. Moreover, despite known sex dimorphism in neurological diseases, females are rarely included in preclinical epilepsy models.

METHODS

We characterized sex differences in mice carrying a pathogenic knockin variant (p.N1768D) in the Scn8a gene that causes spontaneous tonic-clonic seizures (TCs) at ∼3 months of age and found that heterozygous females are more resilient than males in mortality and morbidity. To investigate the cellular mechanisms that underlie female resilience, we utilized blood-brain barrier (BBB) and hippocampal transcriptomic analyses in heterozygous mice before seizure onset (pre-TC) and in mice that experienced ∼20 TCs (post-TC).

RESULTS

In the pre-TC latent phase, both sexes exhibited leaky BBB; however, patterns of gene expression were sexually dimorphic. Females exhibited enhanced oxidative phosphorylation and protein biogenesis, while males activated gliosis and CREB signaling. After seizure onset (chronic phase), females exhibited a metabolic switch to lipid metabolism, while males exhibited increased gliosis and BBB dysfunction and a strong activation of neuroinflammatory pathways.

CONCLUSION

The results underscore the central role of oxidative stress and BBB permeability in the early stages of epileptogenesis, as well as sex dimorphism in response to increasing neuronal hyperexcitability. Our results also highlight the need to include both sexes in preclinical studies to effectively translate results of drug efficacy studies.

Benchmarking RNA-Seq Aligners at Base-Level and Junction Base-Level Resolution Using the Arabidopsis thaliana Genome

The utmost goal of selecting an RNA-Seq alignment software is to perform accurate alignments with a robust algorithm, which is capable of detecting the various intricacies underlying read-mapping procedures and beyond. Most alignment software tools are typically pre-tuned with human or prokaryotic data, and therefore may not be suitable for applications to other organisms, such as plants. The rapidly growing plant RNA-Seq databases call for the assessment of the alignment tools on curated plant data, which will aid the calibration of these tools for applications to plant transcriptomic data. We therefore focused here on benchmarking RNA-Seq read alignment tools, using simulated data derived from the model organism Arabidopsis thaliana. We assessed the performance of five popular RNA-Seq alignment tools that are currently available, based on their usage (citation count). By introducing annotated single nucleotide polymorphisms (SNPs) from The Arabidopsis Information Resource (TAIR), we recorded alignment accuracy at both base-level and junction base-level resolutions for each alignment tool. In addition to assessing the performance of the alignment tools at their default settings, accuracies were also recorded by varying the values of numerous parameters, including the confidence threshold and the level of SNP introduction. The performances of the aligners were found consistent under various testing conditions at the base-level accuracy; however, the junction base-level assessment produced varying results depending upon the applied algorithm. At the read base-level assessment, the overall performance of the aligner STAR was superior to other aligners, with the overall accuracy reaching over 90% under different test conditions. On the other hand, at the junction base-level assessment, SubRead emerged as the most promising aligner, with an overall accuracy over 80% under most test conditions.

Delayed maturation and migration of excitatory neurons in the juvenile mouse paralaminar amygdala

The human amygdala paralaminar nucleus (PL) contains many immature excitatory neurons that undergo prolonged maturation from birth to adulthood. We describe a previously unidentified homologous PL region in mice that contains immature excitatory neurons and has previously been considered part of the amygdala intercalated cell clusters or ventral endopiriform cortex. Mouse PL neurons are born embryonically, not from postnatal neurogenesis, despite a subset retaining immature molecular and morphological features in adults. During juvenile-adolescent ages (P21-P35), the majority of PL neurons undergo molecular, structural, and physiological maturation, and a subset of excitatory PL neurons migrate into the adjacent endopiriform cortex. Alongside these changes, PL neurons develop responses to aversive and appetitive olfactory stimuli. The presence of this homologous region in both humans and mice points to the significance of this conserved mechanism of neuronal maturation and migration during adolescence, a key time period for amygdala circuit maturation and related behavioral changes.

Single-nucleus RNA sequencing demonstrates an autosomal dominant Alzheimer's disease profile and possible mechanisms of disease protection

Highly penetrant autosomal dominant Alzheimer's disease (ADAD) comprises a distinct disease entity as compared to the far more prevalent form of AD in which common variants collectively contribute to risk. The downstream pathways that distinguish these AD forms in specific cell types have not been deeply explored. We compared single-nucleus transcriptomes among a set of 27 cases divided among PSEN1-E280A ADAD carriers, sporadic AD, and controls. Autophagy genes and chaperones clearly defined the PSEN1-E280A cases compared to sporadic AD. Spatial transcriptomics validated the activation of chaperone-mediated autophagy genes in PSEN1-E280A. The PSEN1-E280A case in which much of the brain was spared neurofibrillary pathology and harbored a homozygous APOE3-Christchurch variant revealed possible explanations for protection from AD pathology including overexpression of LRP1 in astrocytes, increased expression of FKBP1B, and decreased PSEN1 expression in neurons. The unique cellular responses in ADAD and sporadic AD require consideration when designing clinical trials.

Higher Maternal Body Mass Index Is Associated With Lower Placental Expression of EPYC: A Genome-Wide Transcriptomic Study

CONTEXT

Elevated body mass index (BMI) in pregnancy is associated with adverse maternal and fetal outcomes. The placental transcriptome may elucidate molecular mechanisms underlying these associations.

OBJECTIVE

We examined the association of first-trimester maternal BMI with the placental transcriptome in the Gen3G prospective cohort.

METHODS

We enrolled participants at 5 to 16 weeks of gestation and measured height and weight. We collected placenta samples at delivery. We performed whole-genome RNA sequencing using Illumina HiSeq 4000 and aligned RNA sequences based on the GTEx v8 pipeline. We conducted differential gene expression analysis of over 15 000 genes from 450 placental samples and reported the change in normalized gene expression per 1-unit increase in log2 BMI (kg/m2) as a continuous variable using Limma Voom. We adjusted models for maternal age, fetal sex, gestational age at delivery, gravidity, and surrogate variables accounting for technical variability. We compared participants with BMI of 18.5 to 24.9 mg/kg2 (N = 257) vs those with obesity (BMI ≥30 kg/m2, N = 82) in secondary analyses.

RESULTS

Participants' mean ± SD age was 28.2 ± 4.4 years and BMI was 25.4 ± 5.5 kg/m2 in early pregnancy. Higher maternal BMI was associated with lower placental expression of EPYC (slope = -1.94, false discovery rate [FDR]-adjusted P = 7.3 × 10-6 for continuous BMI; log2 fold change = -1.35, FDR-adjusted P = 3.4 × 10-3 for BMI ≥30 vs BMI 18.5-24.9 kg/m2) and with higher placental expression of IGFBP6, CHRDL1, and CXCL13 after adjustment for covariates and accounting for multiple testing (FDR < 0.05).

CONCLUSION

Our genome-wide transcriptomic study revealed novel genes potentially implicated in placental biologic response to higher maternal BMI in early pregnancy.

Wolfberry genome database: integrated genomic datasets for studying molecular biology

Wolfberry, also known as goji berry or Lycium barbarum, is a highly valued fruit with significant health benefits and nutritional value. For more efficient and comprehensive usage of published L. barbarum genomic data, we established the Wolfberry database. The utility of the Wolfberry Genome Database (WGDB) is highlighted through the Genome browser, which enables the user to explore the L. barbarum genome, browse specific chromosomes, and access gene sequences. Gene annotation features provide comprehensive information about gene functions, locations, expression profiles, pathway involvement, protein domains, and regulatory transcription factors. The transcriptome feature allows the user to explore gene expression patterns using transcripts per kilobase million (TPM) and fragments per kilobase per million mapped reads (FPKM) metrics. The Metabolism pathway page provides insights into metabolic pathways and the involvement of the selected genes. In addition to the database content, we also introduce six analysis tools developed for the WGDB. These tools offer functionalities for gene function prediction, nucleotide and amino acid BLAST analysis, protein domain analysis, GO annotation, and gene expression pattern analysis. The WGDB is freely accessible at https://cosbi7.ee.ncku.edu.tw/Wolfberry/. Overall, WGDB serves as a valuable resource for researchers interested in the genomics and transcriptomics of L. barbarum. Its user-friendly web interface and comprehensive data facilitate the exploration of gene functions, regulatory mechanisms, and metabolic pathways, ultimately contributing to a deeper understanding of wolfberry and its potential applications in agronomy and nutrition.

Global Transcriptome Analysis of the Peach (Prunus persica) in the Interaction System of Fruit-Chitosan-Monilinia fructicola

The peach (Prunus persica L.) is one of the most important stone-fruit crops worldwide. Nevertheless, successful peach fruit production is seriously reduced by losses due to Monilinia fructicola the causal agent of brown rot. Chitosan has a broad spectrum of antimicrobial properties and may also act as an elicitor that activate defense responses in plants. As little is known about the elicitation potential of chitosan in peach fruits and its impact at their transcriptional-level profiles, the aim of this study was to uncover using RNA-seq the induced responses regulated by the action of chitosan in fruit-chitosan-M. fructicola interaction. Samples were obtained from fruits treated with chitosan or inoculated with M. fructicola, as well from fruits pre-treated with chitosan and thereafter inoculated with the fungus. Chitosan was found to delay the postharvest decay of fruits, and expression profiles showed that its defense-priming effects were mainly evident after the pathogen challenge, driven particularly by modulations of differentially expressed genes (DEGs) related to cell-wall modifications, pathogen perception, and signal transduction, preventing the spread of fungus. In contrast, as the compatible interaction of fruits with M. fructicola was challenged, a shift towards defense responses was triggered with a delay, which was insufficient to limit fungal expansion, whereas DEGs involved in particular processes have facilitated early pathogen colonization. Physiological indicators of peach fruits were also measured. Additionally, expression profiles of particular M. fructicola genes highlight the direct antimicrobial activity of chitosan against the fungus. Overall, the results clarify the possible mechanisms of chitosan-mediated tolerance to M. fructicola and set new foundations for the potential employment of chitosan in the control of brown rot in peaches.

Corrigendum: Identification of ACKR4 as an immune checkpoint in pulmonary arterial hypertension

[This corrects the article DOI: 10.3389/fimmu.2023.1153573.].

Intrahepatic Transcriptomics Differentiate Advanced Fibrosis and Clinical Outcomes in Adults With Fontan Circulation

BACKGROUND

The molecular mechanisms underlying Fontan-associated liver disease (FALD) remain largely unknown.

OBJECTIVES

This study aimed to assess intrahepatic transcriptomic differences among patients with FALD according to the degree of liver fibrosis and clinical outcomes.

METHODS

This retrospective cohort study included adults with the Fontan circulation. Baseline clinical, laboratory, imaging, and hemodynamic data as well as a composite clinical outcome (CCO) were extracted from medical records. Patients were classified into early or advanced fibrosis. RNA was isolated from formalin-fixed paraffin-embedded liver biopsy samples; RNA libraries were constructed with the use of an rRNA depletion method and sequenced on an Illumina Novaseq 6000. Differential gene expression and gene ontology analyses were performed with the use of DESeq2 and Metascape.

RESULTS

A total of 106 patients (48% male, median age 31 years [IQR: 11.3 years]) were included. Those with advanced fibrosis had higher B-type natriuretic peptide levels and Fontan, mean pulmonary artery, and capillary wedge pressures. The CCO was present in 23 patients (22%) and was not predicted by advanced liver fibrosis, right ventricular morphology, presence of aortopulmonary collaterals, or Fontan pressures on multivariable analysis. Samples with advanced fibrosis had 228 upregulated genes compared with early fibrosis. Samples with the CCO had 894 upregulated genes compared with those without the CCO. A total of 136 upregulated genes were identified in both comparisons and were enriched in cellular response to cytokine stimulus or oxidative stress, VEGFA-VEGFR2 signaling pathway, TGF-β signaling pathway, and vasculature development.

CONCLUSIONS

Patients with FALD and advanced fibrosis or the CCO exhibited upregulated genes related to inflammation, congestion, and angiogenesis.

Characterization of non-invasive oropharyngeal samples and nucleic acid isolation for molecular diagnostics

Molecular diagnostics is an increasingly important clinical tool, especially in routine sampling. We evaluated two non-invasive methods (oral swabs and mouthwashes) for sampling nucleic acids from the oral/pharyngeal area. We created a workflow from sample collection (n = 59) to RT-qPCR based analysis. The samples were further characterized in terms of their cellular composition as well as the purity, degradation and microbial content of the derived DNA/RNA. We determined the optimal housekeeping genes applicable for these types of samples. The cellular composition indicated that mouthwashes contained more immune cells and bacteria. Even though the protocol was not specifically optimized to extract bacterial RNA it was possible to derive microbial RNA, from both sampling methods. Optimizing the protocol allowed us to generate stable quantities of DNA/RNA. DNA/RNA purity parameters were not significantly different between the two sampling methods. Even though integrity analysis demonstrated a high level of degradation of RNA, corresponding parameters confirmed their sequencing potential. RT-qPCR analysis determined TATA-Box Binding Protein as the most favorable housekeeping gene. In summary, we have developed a robust method suitable for multiple downstream diagnostic techniques. This protocol can be used as a foundation for further research endeavors focusing on developing molecular diagnostics for the oropharyngeal cavity.

Longitudinal Transcriptomic, Proteomic, and Metabolomic Response of Citrus sinensis to Diaphorina citri Inoculation of Candidatus Liberibacter asiaticus

Huanglongbing (HLB) is a fatal citrus disease that is currently threatening citrus varieties worldwide. One putative causative agent, Candidatus Liberibacter asiaticus (CLas), is vectored by Diaphorina citri, known as the Asian citrus psyllid (ACP). Understanding the details of CLas infection in HLB disease has been hindered by its Candidatus nature and the inability to confidently detect it in diseased trees during the asymptomatic stage. To identify early changes in citrus metabolism in response to inoculation of CLas using its natural psyllid vector, leaves from Madam Vinous sweet orange (Citrus sinensis (L.) Osbeck) trees were exposed to CLas-positive ACP or CLas-negative ACP and longitudinally analyzed using transcriptomics (RNA sequencing), proteomics (liquid chromatography-tandem mass spectrometry; data available in Dryad: 10.25338/B83H1Z), and metabolomics (proton nuclear magnetic resonance). At 4 weeks postexposure (wpe) to psyllids, the initial HLB plant response was primarily to the ACP and, to a lesser extent, the presence or absence of CLas. Additionally, analysis of 4, 8, 12, and 16 wpe identified 17 genes and one protein as consistently differentially expressed between leaves exposed to CLas-positive ACP versus CLas-negative ACP. This study informs identification of early detection molecular targets and contributes to a broader understanding of vector-transmitted plant pathogen interactions.

RNA-seq analysis of parasitism by Intoshia linei (Orthonectida) reveals protein effectors of defence, communication, feeding and growth

Orthonectida is a group of multicellular endoparasites of a wide range of marine invertebrates. Their parasitic stage is a multinuclear shapeless plasmodium infiltrating host tissues. The development of the following worm-like sexual generation takes place within the cytoplasm of the plasmodium. The existence of the plasmodial stage and the development of a sexual stage within the plasmodium are unique features to Bilateria. However, the molecular mechanisms that maintain this peculiar organism, and hence enable parasitism in orthonectids, are unknown. Here, we present the first-ever RNA-seq analysis of the plasmodium, aimed at the identification and characterization of the plasmodium-specific protein-coding genes and corresponding hypothetical proteins that distinguish the parasitic plasmodium stage from the sexual stage of the orthonectid Intoshia linei Giard, 1877, parasite of nemertean Lineus ruber Müller, 1774. We discovered 119 plasmodium-specific proteins, 82 of which have inferred functions based on known domains. Thirty-five of the detected proteins are orphans, at least part of which may reflect the unique evolutionary adaptations of orthonectids to parasitism. Some of the identified proteins are known effector molecules of other endoparasites suggesting convergence. Our data indicate that the plasmodium-specific proteins might be involved in the plasmodium defense against the host, host-parasite communication, feeding and nutrient uptake, growth within the host, and support of the sexual stage development. These molecular processes in orthonectids have not been described before, and the particular protein effectors remained unknown until now.

Multi-omics resources for the Australian southern stuttering frog (Mixophyes australis) reveal assorted antimicrobial peptides

The number of genome-level resources for non-model species continues to rapidly expand. However, frog species remain underrepresented, with up to 90% of frog genera having no genomic or transcriptomic data. Here, we assemble the first genomic and transcriptomic resources for the recently described southern stuttering frog (Mixophyes australis). The southern stuttering frog is ground-dwelling, inhabiting naturally vegetated riverbanks in south-eastern Australia. Using PacBio HiFi long-read sequencing and Hi-C scaffolding, we generated a high-quality genome assembly, with a scaffold N50 of 369.3 Mb and 95.1% of the genome contained in twelve scaffolds. Using this assembly, we identified the mitochondrial genome, and assembled six tissue-specific transcriptomes. We also bioinformatically characterised novel sequences of two families of antimicrobial peptides (AMPs) in the southern stuttering frog, the cathelicidins and β-defensins. While traditional peptidomic approaches to peptide discovery have typically identified one or two AMPs in a frog species from skin secretions, our bioinformatic approach discovered 12 cathelicidins and two β-defensins that were expressed in a range of tissues. We investigated the novelty of the peptides and found diverse predicted activities. Our bioinformatic approach highlights the benefits of multi-omics resources in peptide discovery and contributes valuable genomic resources in an under-represented taxon.

A transcriptomics-based drug repositioning approach to identify drugs with similar activities for the treatment of muscle pathologies in spinal muscular atrophy (SMA) models

Spinal muscular atrophy (SMA) is a genetic neuromuscular disorder caused by the reduction of survival of motor neuron (SMN) protein levels. Although three SMN-augmentation therapies are clinically approved that significantly slow down disease progression, they are unfortunately not cures. Thus, complementary SMN-independent therapies that can target key SMA pathologies and that can support the clinically approved SMN-dependent drugs are the forefront of therapeutic development. We have previously demonstrated that prednisolone, a synthetic glucocorticoid (GC) improved muscle health and survival in severe Smn-/-;SMN2 and intermediate Smn2B/- SMA mice. However, long-term administration of prednisolone can promote myopathy. We thus wanted to identify genes and pathways targeted by prednisolone in skeletal muscle to discover clinically approved drugs that are predicted to emulate prednisolone's activities. Using an RNA-sequencing, bioinformatics, and drug repositioning pipeline on skeletal muscle from symptomatic prednisolone-treated and untreated Smn-/-; SMN2 SMA and Smn+/-; SMN2 healthy mice, we identified molecular targets linked to prednisolone's ameliorative effects and a list of 580 drug candidates with similar predicted activities. Two of these candidates, metformin and oxandrolone, were further investigated in SMA cellular and animal models, which highlighted that these compounds do not have the same ameliorative effects on SMA phenotypes as prednisolone; however, a number of other important drug targets remain. Overall, our work further supports the usefulness of prednisolone's potential as a second-generation therapy for SMA, identifies a list of potential SMA drug treatments and highlights improvements for future transcriptomic-based drug repositioning studies in SMA.

Identification of unique genomic signatures in patients with fibromyalgia and chronic pain

Fibromyalgia (FM) is a chronic pain syndrome characterized by widespread pain. The pathophysiology of fibromyalgia is not clearly understood and there are no specific biomarkers available for accurate diagnosis. Here we define genomic signatures using high throughput RNA sequencing on 96 fibromyalgia and 93 control cases. Our findings revealed three major fibromyalgia-associated expression signatures. The first group included 43 patients with a signature enriched for gene expression associated with extracellular matrix and downregulation of RhoGDI signaling pathway. The second group included 30 patients and showed a profound reduction in the expression of inflammatory mediators with an increased expression of genes involved in the CLEAR signaling pathway. These results suggest defective tissue homeostasis associated with the extra-cellular matrix and cellular program that regulates lysosomal biogenesis and participates in macromolecule clearance in fibromyalgia. The third group of 17 FM patients showed overexpression of pathways that control acute inflammation and dysfunction of the global transcriptional process. The result of this study indicates that FM is a heterogeneous and complex disease. Further elucidation of these pathways will lead to the development of accurate diagnostic markers, and effective therapeutic options for fibromyalgia.

Whole blood RNA signatures in tuberculosis patients receiving H56:IC31 vaccine as adjunctive therapy

Introduction

Therapeutic vaccination in tuberculosis (TB) represents a Host Directed Therapy strategy which enhances immune responses in order to improve clinical outcomes and shorten TB treatment. Previously, we have shown that the subunit H56:IC31 vaccine induced both humoral and cellular immune responses when administered to TB patients adjunctive to standard TB treatment (TBCOX2 study, NCT02503839). Here we present the longitudinal whole blood gene expression patterns in H56:IC31 vaccinated TB patients compared to controls receiving standard TB treatment only.

Methods

The H56:IC31 group (N=11) and Control group (N=7) underwent first-line TB treatment for 182 days. The H56:IC31 group received 5 micrograms of the H56:IC31 vaccine (Statens Serum Institut; SSI, Valneva Austria GmbH) intramuscularly at day 84 and day 140. Total RNA was extracted from whole blood samples collected in PAXgene tubes on days 0, 84, 98, 140, 154, 182 and 238. The expression level of 183 immune-related genes was measured by high-throughput microfluidic qPCR (Biomark HD system, Standard BioTools).

Results

The targeted gene expression profiling unveiled the upregulation of modules such as interferon (IFN) signalling genes, pattern recognition receptors and small nucleotide guanosine triphosphate (GTP)-ases in the vaccinated group compared to controls two weeks after administration of the first H56:IC31 vaccine. Additionally, the longitudinal analysis of the Adolescent Cohort Study-Correlation of Risk (ACS-COR) signature showed a progressive downregulation in both study arms towards the end of TB treatment, in congruence with reported treatment responses and clinical improvements. Still, two months after the end of TB treatment, vaccinated patients, and especially those developing both cellular and humoral vaccine responses, showed a lower expression of the ACS-COR genes compared to controls.

Discussion

Our data report gene expression patterns following H56:IC31 vaccination which might be interpreted as a lower risk of relapse in therapeutically vaccinated patients. Further studies are needed to conclude if these gene expression patterns could be used as prognostic biosignatures for therapeutic TB vaccine responses.

Exploring the impact of the PNPLA3 I148M variant on primary human hepatic stellate cells using 3D extracellular matrix models

BACKGROUND & AIMS

The PNPLA3 rs738409 C>G (encoding for I148M) variant is a risk locus for the fibrogenic progression of chronic liver diseases, a process driven by hepatic stellate cells (HSCs). We investigated how the PNPLA3 I148M variant affects HSC biology using transcriptomic data and validated findings in 3D-culture models.

METHODS

RNA sequencing was performed on 2D-cultured primary human HSCs and liver biopsies of individuals with obesity, genotyped for the PNPLA3 I148M variant. Data were validated in wild-type (WT) or PNPLA3 I148M variant-carrying HSCs cultured on 3D extracellular matrix (ECM) scaffolds from human healthy and cirrhotic livers, with/without TGFB1 or cytosporone B (Csn-B) treatment.

RESULTS

Transcriptomic analyses of liver biopsies and HSCs highlighted shared PNPLA3 I148M-driven dysregulated pathways related to mitochondrial function, antioxidant response, ECM remodelling and TGFB1 signalling. Analogous pathways were dysregulated in WT/PNPLA3-I148M HSCs cultured in 3D liver scaffolds. Mitochondrial dysfunction in PNPLA3-I148M cells was linked to respiratory chain complex IV insufficiency. Antioxidant capacity was lower in PNPLA3-I148M HSCs, while reactive oxygen species secretion was increased in PNPLA3-I148M HSCs and higher in bioengineered cirrhotic vs. healthy scaffolds. TGFB1 signalling followed the same trend. In PNPLA3-I148M cells, expression and activation of the endogenous TGFB1 inhibitor NR4A1 were decreased: treatment with the Csn-B agonist increased total NR4A1 in HSCs cultured in healthy but not in cirrhotic 3D scaffolds. NR4A1 regulation by TGFB1/Csn-B was linked to Akt signalling in PNPLA3-WT HSCs and to Erk signalling in PNPLA3-I148M HSCs.

CONCLUSION

HSCs carrying the PNPLA3 I148M variant have impaired mitochondrial function, antioxidant responses, and increased TGFB1 signalling, which dampens antifibrotic NR4A1 activity. These features are exacerbated by cirrhotic ECM, highlighting the dual impact of the PNPLA3 I148M variant and the fibrotic microenvironment in progressive chronic liver diseases.

IMPACT AND IMPLICATIONS

Hepatic stellate cells (HSCs) play a key role in the fibrogenic process associated with chronic liver disease. The PNPLA3 genetic mutation has been linked with increased risk of fibrogenesis, but its role in HSCs requires further investigation. Here, by using comparative transcriptomics and a novel 3D in vitro model, we demonstrate the impact of the PNPLA3 genetic mutation on primary human HSCs' behaviour, and we show that it affects the cell's mitochondrial function and antioxidant response, as well as the antifibrotic gene NR4A1. Our publicly available transcriptomic data, 3D platform and our findings on NR4A1 could facilitate the discovery of targets to develop more effective treatments for chronic liver diseases.

Differential responses of selectively bred mussels (Perna canaliculus) to heat stress-survival, immunology, gene expression and microbiome diversity

New Zealand's green-lipped mussel (Perna canaliculus) is an ecologically and economically important species. Marine heatwaves are increasing in frequency around NZ's coastline, and these events are correlated with increased stress and mortality of some aquaculture species. This study aimed to identify general biomarkers of heat stress in P. canaliculus and to assess whether responses differed between genetically distinct selectively bred mussels. We exposed three families of selectively bred mussels (families A, B and C) to three seawater temperature regimes in the laboratory: 1) a "control" treatment (ambient 12°C), 2) a 26°C heat challenge with a subsequent recovery period, and 3) a sustained 26°C heat challenge with no recovery. We investigated whether the survival, immune response (hemocyte concentration and viability, oxidative stress and total antioxidant capacity), hemocyte gene expression and gill microbiome differed between the families during the temperature challenges. In the sustained heat-stress treatment, family A had the highest survival rate (42% compared with 25% and 5% for families C and B, respectively). Gene expression levels significantly shifted during thermal stress and differed between families, with family A more dissimilar than families B and C. Family C had substantially more genes impacted by temperature treatment and timepoint than the other families, while family B had very little genes/pathways that responded to thermal stress. Genes related to heat shock proteins and immune responses (e.g., AIF1, CTSC, TOLL8, CASP9, FNTA, AHCY, CRYAB, PPIF) were upregulated in all families during heat stress. Microbiome species-richness differed between families before and during heat-stress, with family A having a distinctly different microbiome flora than the other families. Microbial diversity changed similarly in all families exposed to prolonged heat-stress, with species of Vibrio and Campylobacter increasing in these mussels. Our study highlights the use of non-lethal sampling of hemocytes as a diagnostic tool to explore the immune response and gene expression of selectively bred mussels, to predict their response to ocean warming. This approach can identify potential thermotolerant candidates for further selective breeding, which may increase the resilience of the mussel aquaculture industry in a warming ocean.

High resolution spatial profiling of kidney injury and repair using RNA hybridization-based in situ sequencing

Emerging spatially resolved transcriptomics technologies allow for the measurement of gene expression in situ at cellular resolution. We apply direct RNA hybridization-based in situ sequencing (dRNA HybISS, Cartana part of 10xGenomics) to compare male and female healthy mouse kidneys and the male kidney injury and repair timecourse. A pre-selected panel of 200 genes is used to identify cell state dynamics patterns during injury and repair. We develop a new computational pipeline, CellScopes, for the rapid analysis, multi-omic integration and visualization of spatially resolved transcriptomic datasets. The resulting dataset allows us to resolve 13 kidney cell types within distinct kidney niches, dynamic alterations in cell state over the course of injury and repair and cell-cell interactions between leukocytes and kidney parenchyma. At late timepoints after injury, C3+ leukocytes are enriched near pro-inflammatory, failed-repair proximal tubule cells. Integration of snRNA-seq dataset from the same injury and repair samples also allows us to impute the spatial localization of genes not directly measured by dRNA HybISS.

Exploring JC polyomavirus sequences and human gene expression in brain tissue of patients with progressive multifocal leukoencephalopathy

Progressive multifocal leukoencephalopathy (PML) is a rare neurological condition associated with reactivation of dormant JC polyomavirus (JCPyV). In this study, we characterized gene expression and JCPyV rearrangements in PML brain tissue. Infection of white matter astrocytes and oligodendrocytes as well as occasional brain cortex neurons was shown. PML brain harbored exclusively rearranged JCPyV variants. Viral transcripts covered the whole genome on both strands. Strong differential expression of human genes associated with neuroinflammation, blood-brain-barrier permeability and neurodegenerative diseases was shown. Pathway analysis revealed wide immune activation in PML brain. The study provides novel insights into the pathogenesis of PML.