RNA sequencing: advances, challenges and opportunities

Thousands of oscillating LncRNAs in the mouse testis

The long noncoding RNAs (lncRNAs) are involved in numerous fundamental biological processes, including circadian regulation. Although recent studies have revealed insights into the functions of lncRNAs, how the lncRNAs regulate circadian rhythms still requires a deeper investigation. In this study, we generate two datasets of RNA-seq profiles of the mouse (Mus musculus) testis under light-dark (LD) cycle. The first dataset included 18,613 unannotated transcripts measured at 12 time points, each with duplicate samples, under LD conditions; while the second dataset included 21,414 unannotated transcripts measured at six time points, each with three replicates, under desynchronized and control conditions. We identified 5964 testicular lncRNAs in each dataset by BLASTing these transcripts against the known mouse lncRNAs from the NONCODE database. MetaCycle analyses were performed to identify 519, 475, and 494 rhythmically expressed mouse testicular lncRNAs in the 12-time-point dataset, the six-time-point control dataset, and the six-time-point desynchronized dataset, respectively. A comparison of the expression profiles of the lncRNAs under desynchronized and control conditions revealed that 427 rhythmically expressed lncRNAs from the control condition became arrhythmic under the desynchronized condition, suggesting a possible loss of rhythmicity. In contrast, 446 arrhythmic lncRNAs from the control condition became rhythmic under the desynchronized condition, suggesting a possible gain of rhythmicity. Interestingly, 48 lncRNAs were rhythmically expressed under both desynchronized and control conditions. These oscillating lncRNAs were divided into morning lncRNAs, evening lncRNAs, and night lncRNAs based on their time-course expression patterns. We interrogated the promoter regions of these rhythmically expressed mouse testicular lncRNAs to predict their possible regulation by the E-box, D-box, or RORE promoter motifs. GO and KEGG analyses were performed to identify the possible biological functions of these rhythmically expressed mouse testicular lncRNAs. Further, we conducted conservation analyses of the rhythmically expressed mouse testicular lncRNAs with lncRNAs from humans, rats, and zebrafish, and uncovered three mouse testicular lncRNAs conserved across these four species. Finally, we computationally predicted the conserved lncRNA-encoded peptides and their 3D structures from each of the four species. Taken together, our study revealed thousands of rhythmically expressed lncRNAs in the mouse testis, setting the stage for further computational and experimental validations.

RNA-sequencing exploration on SIR2 and SOD genes in Polyalthia longifolia leaf methanolic extracts (PLME) mediated anti-aging effects in Saccharomyces cerevisiae BY611 yeast cells

Polyalthia longifolia is well-known for its abundance of polyphenol content and traditional medicinal uses. Previous research has demonstrated that the methanolic extract of P. longifolia leaves (PLME, 1 mg/mL) possesses anti-aging properties in Saccharomyces cerevisiae BY611 yeast cells. Building on these findings, this study delves deeper into the potential antiaging mechanism of PLME, by analyzing the transcriptional responses of BY611 cells treated with PLME using RNA-sequencing (RNA-seq) technology. The RNA-seq analysis results identified 1691 significantly (padj < 0.05) differentially expressed genes, with 947 upregulated and 744 downregulated genes. Notably, the expression of three important aging-related genes, SIR2, SOD1, and SOD2, showed a significant difference following PLME treatment. The subsequent integration of these targeted genes with GO and KEGG pathway analysis revealed the multifaceted nature of PLME's anti-aging effects in BY611 yeast cells. Enriched GO and KEGG analysis showed that PLME treatment promotes the upregulation of SIR2, SOD1, and SOD2 genes, leading to a boosted cellular antioxidant defense system, reduced oxidative stress, regulated cell metabolism, and maintain genome stability. These collectively increased longevities in PLME-treated BY611 yeast cells and indicate the potential anti-aging action of PLME through the modulation of SIR2 and SOD genes. The present study provided novel insights into the roles of SIR2, SOD1, and SOD2 genes in the anti-aging effects of PLME treatment, offering promising interventions for promoting healthy aging.

Engineering circular RNA for molecular and metabolic reprogramming

The role of messenger RNA (mRNA) in biological systems is extremely versatile. However, it's extremely short half-life poses a fundamental restriction on its application. Moreover, the translation efficiency of mRNA is also limited. On the contrary, circular RNAs, also known as circRNAs, are a common and stable form of RNA found in eukaryotic cells. These molecules are synthesized via back-splicing. Both synthetic circRNAs and certain endogenous circRNAs have the potential to encode proteins, hence suggesting the potential of circRNA as a gene expression machinery. Herein, we aim to summarize all engineering aspects that allow exogenous circular RNA (circRNA) to prolong the time that proteins are expressed from full-length RNA signals. This review presents a systematic engineering approach that have been devised to efficiently assemble circRNAs and evaluate several aspects that have an impact on protein production derived from. We have also reviewed how optimization of the key components of circRNAs, including the topology of vector, 5' and 3' untranslated sections, entrance site of the internal ribosome, and engineered aptamers could be efficiently impacting the translation machinery for molecular and metabolic reprogramming. Collectively, molecular and metabolic reprogramming present a novel way of regulating distinctive cellular features, for instance growth traits to neoplastic cells, and offer new possibilities for therapeutic inventions.

Combined transcriptomic and metabolomic analysis revealed that pH changes affected the expression of carbohydrate and ribosome biogenesis-related genes in Aspergillus niger SICU-33

The process of carbohydrate metabolism and genetic information transfer is an important part of the study on the effects of the external environment on microbial growth and development. As one of the most significant environmental parameters, pH has an important effect on mycelial growth. In this study, the effects of environmental pH on the growth and nutrient composition of Aspergillus niger (A. niger) filaments were determined. The pH values of the medium were 5, 7, and 9, respectively, and the molecular mechanism was further investigated by transcriptomics and metabolomics methods. The results showed that pH 5 and 9 significantly inhibited filament growth and polysaccharide accumulation of A. niger. Further, the mycelium biomass of A. niger and the crude polysaccharide content was higher when the medium's pH was 7. The DEGs related to ribosome biogenesis were the most abundant, and the downregulated expression of genes encoding XRN1, RRM, and RIO1 affected protein translation, modification, and carbohydrate metabolism in fungi. The dynamic changes of pargyline and choline were in response to the oxidative metabolism of A. niger SICU-33. The ribophorin_I enzymes and DL-lactate may be important substances related to pH changes during carbohydrate metabolism of A.niger SICU-33. The results of this study provide useful transcriptomic and metabolomic information for further analyzing the bioinformatic characteristics of A. niger and improving the application in ecological agricultural fermentation.

Machine Learning-Assisted Direct RNA Sequencing with Epigenetic RNA Modification Detection via Quantum Tunneling

RNA sequence information holds immense potential as a drug target for diagnosing various RNA-mediated diseases and viral/bacterial infections. Massively parallel complementary DNA (c-DNA) sequencing helps but results in a loss of valuable information about RNA modifications, which are often associated with cancer evolution. Herein, we report machine learning (ML)-assisted high throughput RNA sequencing with inherent RNA modification detection using a single-molecule, long-read, and label-free quantum tunneling approach. The ML tools achieve classification accuracy as high as 100% in decoding RNA and 98% for decoding both RNA and RNA modifications simultaneously. The relationships between input features and target values have been well examined through Shapley additive explanations. Furthermore, transmission and sensitivity readouts enable the recognition of RNA and its modifications with good selectivity and sensitivity. Our approach represents a starting point for ML-assisted direct RNA sequencing that can potentially decode RNA and its epigenetic modifications at the molecular level.

Screening of reference genes for gene expression study in different tissues from the transcriptome data of the vector leafhopper Psammotettix striatus

Selecting appropriate reference genes is crucial for ensuring the accuracy and reliability of gene expression study using reverse transcription-quantitative PCR (RT-qPCR). To screen the optimal reference genes for analyzing gene expression in different tissues of the vector leafhopper Psammotettix striatus which causes extensive damage to a wide range of crops by vectoring multiple plant pathogenic microorganisms, the transcriptome data from Malpighian tubules (MTs) of P. striatus were mined. Twenty alternative candidate reference genes were initially selected for screening, among which seven genes with diverse Gene Ontology (GO) annotations were choosed as candidate reference genes, i.e., ribosomal protein L7A (RPL7A), ribosomal protein S28 (RPS28), ribosomal protein L22 (RPL22), ribosomal protein LP2 (RPLP2), H3 histone family 3A (H3F3A), elongation factor 1γ (EF-1γ), and elongation factor 1α (EF-1α). Gene expression levels in different tissues of P. striatus adults were examined using RT-qPCR, and their expression stability was analyzed using multiple reference gene screening software. This study revealed EF-1α as the most abundantly expressed gene, while RPL22 exhibited the lowest expression levels. EF-1α showed the most stable expression, whereas RPS28 showed the least stability. Various software tools confirmed EF-1α as the most stable single reference gene, and EF-1α and RPLP2 an optimal combination. This study provides a foundation for future investigation of the transmission of pathogenic microorganisms mediated by the vector leafhoppers, the function of the MTs, the biosynthesis of brochosomes, the coevolutionary processes and nutritional interactions of symbionts and host insects, and the gene expression study of other sap-sucking insects.

Next generation sequencing technologies to address aberrant mRNA translation in cancer

In this review, we explore the transformative impact of next generation sequencing technologies in the realm of translatomics (the study of how translational machinery acts on a genome-wide scale). Despite the expectation of a direct correlation between mRNA and protein content, the complex regulatory mechanisms that affect this relationship remark the limitations of standard RNA-seq approaches. Then, the review characterizes crucial techniques such as polysome profiling, ribo-seq, trap-seq, proximity-specific ribosome profiling, rnc-seq, tcp-seq, qti-seq and scRibo-seq. All these methods are summarized within the context of cancer research, shedding light on their applications in deciphering aberrant translation in cancer cells. In addition, we encompass databases and bioinformatic tools essential for researchers that want to address translatome analysis in the context of cancer biology.

Microbial metabolite deoxycholic acid-mediated ferroptosis exacerbates high-fat diet-induced colonic inflammation

High-fat diet (HFD) has long been recognized as risk factors for the development and progression of ulcerative colitis (UC), but the exact mechanism remained elusive. Here, HFD increased intestinal deoxycholic acid (DCA) levels, and DCA further exacerbated colonic inflammation. Transcriptome analysis revealed that DCA triggered ferroptosis pathway in colitis mice. Mechanistically, DCA upregulated hypoxia-inducible factor-2α (HIF-2α) and divalent metal transporter-1 (DMT1) expression, causing the ferrous ions accumulation and ferroptosis in intestinal epithelial cells, which was reversed by ferroptosis inhibitor ferrostatin-1. DCA failed to promote colitis and ferroptosis in intestine-specific HIF-2α-null mice. Notably, byak-angelicin inhibited DCA-induced pro-inflammatory and pro-ferroptotic effects through blocking the up-regulation of HIF-2α by DCA. Moreover, fat intake was positively correlated with disease activity in UC patients consuming HFD, with ferroptosis being more pronounced. Collectively, our findings demonstrated that HFD exacerbated colonic inflammation by promoting DCA-mediated ferroptosis, providing new insights into diet-related bile acid dysregulation in UC.

Next-Generation Sequencing-Based Spatial Transcriptomics: A Perspective from Barcoding Chemistry

Gene expression profiling of tissue cells with spatial context is in high demand to reveal cell types, locations, and intercellular or molecular interactions for physiological and pathological studies. With rapid advances in barcoding chemistry and sequencing chemistry, spatially resolved transcriptome (SRT) techniques have emerged to quantify spatial gene expression in tissue samples by correlating transcripts with their spatial locations using diverse strategies. These techniques provide both physical tissue structure and molecular characteristics and are poised to revolutionize many fields, such as developmental biology, neuroscience, oncology, and histopathology. In this context, this Perspective focuses on next-generation sequencing-based SRT methods, particularly highlighting spatial barcoding chemistry. It delves into optically manipulated spatial indexing methods and DNA array-barcoded spatial indexing methods by exploring current advances, challenges, and future development directions in this nascent field.

Transcriptomics confirms IRF1 as a key regulator of pyroptosis in diabetic retinopathy

BACKGROUND

Diabetic retinopathy (DR) is a retinal microvascular complication caused by hyperglycemia, which can lead to visual impairment or blindness. Pyroptosis is a type of inflammation-related programmed cell death, activated by caspase-1, resulting in the maturation of IL-1β and IL-18 and the rupture of the cell membrane. RNA sequencing (RNA-seq) is a high-throughput sequencing technique that reveals the presence and quantity of RNA in the genome at a specific time point, i.e., the transcriptome. RNA-seq can analyze gene expression levels, splicing variants, mutations, fusions, editing and other post-transcriptional modifications, as well as gene expression differences between different samples or conditions. It has been widely used in biological and medical research, clinical diagnosis and new drug development. This study aimed to establish an in vitro model of diabetic retinopathy by culturing human retinal endothelial cells (HREC) with high glucose (30 mmol/L), and to detect their transcriptome expression by RNA-seq, screen for key genes related to pyroptosis, and validate the sequencing results by subsequent experiments.

METHODS

We used RNA-seq to detect the transcriptome expression differences between HREC cells cultured with high glucose and control group, and identified differentially expressed genes by GO/KEGG analysis. We constructed a PPI network and determined the key genes by Cytoscape software and CytoHubba plugin. We validated the expression of related factors by Western Blot, qPCR and ELISA.

RESULTS

We performed GO and KEGG analysis on the RNA-seq data and found differentially expressed genes. We used Cytoscape and CytoHubba plugin to screen out IRF1 as the key gene, and then detected the expression of IRF1 in HREC under high glucose and control group by Western Blot and qPCR. We found that the expression of Caspase-1, GSDMD and IL-1β proteins in HREC under high glucose increased, while the expression of these proteins decreased after the inhibition of IRF1 by siRNA. ELISA showed that the secretion of IL-1β in HREC under high glucose increased, while the inhibition of IRF1 reduced the secretion of IL-1β. These results indicate that IRF1 plays an important role in DR, and provides a new target and strategy for the prevention and treatment of this disease.

ProBac-seq, a bacterial single-cell RNA sequencing methodology using droplet microfluidics and large oligonucleotide probe sets

Methods that measure the transcriptomic state of thousands of individual cells have transformed our understanding of cellular heterogeneity in eukaryotic cells since their introduction in the past decade. While simple and accessible protocols and commercial products are now available for the processing of mammalian cells, these existing technologies are incompatible with use in bacterial samples for several fundamental reasons including the absence of polyadenylation on bacterial messenger RNA, the instability of bacterial transcripts and the incompatibility of bacterial cell morphology with existing methodologies. Recently, we developed ProBac sequencing (ProBac-seq), a method that overcomes these technical difficulties and provides high-quality single-cell gene expression data from thousands of bacterial cells by using messenger RNA-specific probes. Here we provide details for designing large oligonucleotide probe sets for an organism of choice, amplifying probe sets to produce sufficient quantities for repeated experiments, adding unique molecular indexes and poly-A tails to produce finalized probes, in situ probe hybridization and single-cell encapsulation and library preparation. This protocol, from the probe amplification to the library preparation, requires ~7 d to complete. ProBac-seq offers several advantages over other methods by capturing only the desired target sequences and avoiding nondesired transcripts, such as highly abundant ribosomal RNA, thus enriching for signal that better informs on cellular state. The use of multiple probes per gene can detect meaningful single-cell signals from cells expressing transcripts to a lesser degree or those grown in minimal media and other environmentally relevant conditions in which cells are less active. ProBac-seq is also compatible with other organisms that can be profiled by in situ hybridization techniques.

Benchmarking of methods for DNA methylome deconvolution

Defining the number and abundance of different cell types in tissues is important for understanding disease mechanisms as well as for diagnostic and prognostic purposes. Typically, this is achieved by immunohistological analyses, cell sorting, or single-cell RNA-sequencing. Alternatively, cell-specific DNA methylome information can be leveraged to deconvolve cell fractions from a bulk DNA mixture. However, comprehensive benchmarking of deconvolution methods and modalities was not yet performed. Here we evaluate 16 deconvolution algorithms, developed either specifically for DNA methylome data or more generically. We assess the performance of these algorithms, and the effect of normalization methods, while modeling variables that impact deconvolution performance, including cell abundance, cell type similarity, reference panel size, method for methylome profiling (array or sequencing), and technical variation. We observe differences in algorithm performance depending on each these variables, emphasizing the need for tailoring deconvolution analyses. The complexity of the reference, marker selection method, number of marker loci and, for sequencing-based assays, sequencing depth have a marked influence on performance. By developing handles to select the optimal analysis configuration, we provide a valuable source of information for studies aiming to deconvolve array- or sequencing-based methylation data.

Epitranscriptome insights into Riccia fluitans L. (Marchantiophyta) aquatic transition using nanopore direct RNA sequencing

BACKGROUND

Riccia fluitans, an amphibious liverwort, exhibits a fascinating adaptation mechanism to transition between terrestrial and aquatic environments. Utilizing nanopore direct RNA sequencing, we try to capture the complex epitranscriptomic changes undergone in response to land-water transition.

RESULTS

A significant finding is the identification of 45 differentially expressed genes (DEGs), with a split of 33 downregulated in terrestrial forms and 12 upregulated in aquatic forms, indicating a robust transcriptional response to environmental changes. Analysis of N6-methyladenosine (m6A) modifications revealed 173 m6A sites in aquatic and only 27 sites in the terrestrial forms, indicating a significant increase in methylation in the former, which could facilitate rapid adaptation to changing environments. The aquatic form showed a global elongation bias in poly(A) tails, which is associated with increased mRNA stability and efficient translation, enhancing the plant's resilience to water stress. Significant differences in polyadenylation signals were observed between the two forms, with nine transcripts showing notable changes in tail length, suggesting an adaptive mechanism to modulate mRNA stability and translational efficiency in response to environmental conditions. This differential methylation and polyadenylation underline a sophisticated layer of post-transcriptional regulation, enabling Riccia fluitans to fine-tune gene expression in response to its living conditions.

CONCLUSIONS

These insights into transcriptome dynamics offer a deeper understanding of plant adaptation strategies at the molecular level, contributing to the broader knowledge of plant biology and evolution. These findings underscore the sophisticated post-transcriptional regulatory strategies Riccia fluitans employs to navigate the challenges of aquatic versus terrestrial living, highlighting the plant's dynamic adaptation to environmental stresses and its utility as a model for studying adaptation mechanisms in amphibious plants.

Detection and identification of single ribonucleotide monophosphates using a dual in-plane nanopore sensor made in a thermoplastic via replication

We report the generation of ∼8 nm dual in-plane pores fabricated in a thermoplastic via nanoimprint lithography (NIL). These pores were connected in series with nanochannels, one of which served as a flight tube to allow the identification of single molecules based on their molecular-dependent apparent mobilities (i.e., dual in-plane nanopore sensor). Two different thermoplastics were investigated including poly(methyl methacrylate), PMMA, and cyclic olefin polymer, COP, as the substrate for the sensor both of which were sealed using a low glass transition cover plate (cyclic olefin co-polymer, COC) that could be thermally fusion bonded to the PMMA or COP substrate at a temperature minimizing nanostructure deformation. Unique to these dual in-plane nanopore sensors was two pores flanking each side of the nanometer flight tube (50 × 50 nm, width × depth) that was 10 μm in length. The utility of this dual in-plane nanopore sensor was evaluated to not only detect, but also identify single ribonucleotide monophosphates (rNMPs) by using the travel time (time-of-flight, ToF), the resistive pulse event amplitude, and the dwell time. In spite of the relatively large size of these in-plane pores (∼8 nm effective diameter), we could detect via resistive pulse sensing (RPS) single rNMP molecules at a mass load of 3.9 fg, which was ascribed to the unique structural features of the nanofluidic network and the use of a thermoplastic with low relative dielectric constants, which resulted in a low RMS noise level in the open pore current. Our data indicated that the identification accuracy of individual rNMPs was high, which was ascribed to an improved chromatographic contribution to the nano-electrophoresis apparent mobility. With the ToF data only, the identification accuracy was 98.3%. However, when incorporating the resistive pulse sensing event amplitude and dwell time in conjunction with the ToF and analyzed via principal component analysis (PCA), the identification accuracy reached 100%. These findings pave the way for the realization of a novel chip-based single-molecule RNA sequencing technology.

The Microglial Transcriptome of Age-Associated Deep Subcortical White Matter Lesions Suggests a Neuroprotective Response to Blood-Brain Barrier Dysfunction

Age-associated deep-subcortical white matter lesions (DSCLs) are an independent risk factor for dementia, displaying high levels of CD68+ microglia. This study aimed to characterize the transcriptomic profile of microglia in DSCLs and surrounding radiologically normal-appearing white matter (NAWM) compared to non-lesional control white matter. CD68+ microglia were isolated from white matter groups (n = 4 cases per group) from the Cognitive Function and Ageing Study neuropathology cohort using immuno-laser capture microdissection. Microarray gene expression profiling, but not RNA-sequencing, was found to be compatible with immuno-LCM-ed post-mortem material in the CFAS cohort and identified significantly differentially expressed genes (DEGs). Functional grouping and pathway analysis were assessed using the Database for Annotation Visualization and Integrated Discovery (DAVID) software, and immunohistochemistry was performed to validate gene expression changes at the protein level. Transcriptomic profiling of microglia in DSCLs compared to non-lesional control white matter identified 181 significant DEGs (93 upregulated and 88 downregulated). Functional clustering analysis in DAVID revealed dysregulation of haptoglobin-haemoglobin binding (Enrichment score 2.5, p = 0.017), confirmed using CD163 immunostaining, suggesting a neuroprotective microglial response to blood-brain barrier dysfunction in DSCLs. In NAWM versus control white matter, microglia exhibited 347 DEGs (209 upregulated, 138 downregulated), with significant dysregulation of protein de-ubiquitination (Enrichment score 5.14, p < 0.001), implying an inability to maintain protein homeostasis in NAWM that may contribute to lesion spread. These findings enhance understanding of microglial transcriptomic changes in ageing white matter pathology, highlighting a neuroprotective adaptation in DSCLs microglia and a potentially lesion-promoting phenotype in NAWM microglia.

LINK-A: unveiling its functional role and clinical significance in human tumors

LINK-A, also recognized as LINC01139, has emerged as a key oncological lncRNA in cancer. LINK-A is upregulated in solid and liquid tumor samples, including breast cancer, ovarian cancer, glioma, non-small-cell lung cancer, and mantle cell lymphoma. Notably, LINK-A is involved in regulating critical cancer-related pathways, such as AKT and HIF1α signaling, and is implicated in a range of oncogenic activities, including cell proliferation, apoptosis, epithelial-mesenchymal transition (EMT), cell invasion and migration, and glycolysis reprogramming. LINK-A's differential expression and its correlation with clinical features enable it to be a promising biomarker for cancer diagnosis, prognosis, and the stratification of tumor progression. Additionally, LINK-A's contribution to the development of resistance to cancer therapies, including AKT inhibitors and immunotherapy, underscores its potential as a therapeutic target. This review provides a comprehensive overview of the available data on LINK-A, focusing on its molecular regulatory pathways and clinical significance. By exploring the multifaceted nature of LINK-A in cancer, the review aims to offer a valuable resource for future research directions, potentially guiding the development of novel therapeutic strategies targeting this lncRNA in cancer treatment.

Bulk RNA-sequencing, single-cell RNA-sequencing analysis, and experimental validation reveal iron metabolism-related genes CISD2 and CYP17A1 are potential diagnostic markers for recurrent pregnancy loss

BACKGROUND

Recurrent pregnancy loss (RPL) is associated with variable causes. Its etiology remains unexplained in about half of the cases, with no effective treatment available. Individuals with RPL have an irregular iron metabolism. In the present study, we identified key genes impacting iron metabolism that could be used for diagnosing and treating RPL.

METHODS

We obtained gene expression profiles from the Gene Expression Omnibus (GEO) database. The Molecular Signatures Database was used to identify 14 gene sets related to iron metabolism, comprising 520 iron metabolism genes. Differential analysis and a weighted gene co-expression network analysis (WGCNA) of gene expression revealed two iron metabolism-related hub genes. Reverse transcriptase-polymerase chain reaction (RT-PCR) and immunohistochemistry were used on clinical samples to confirm our results. The receiver operating characteristic (ROC) analysis and immune infiltration analysis were conducted. In addition, we analyzed the distribution of genes and performed CellChat analysis by single-cell RNA sequencing.

RESULTS

The expression of two hub genes, namely, CDGSH iron sulfur domain 2 (CISD2)and Cytochrome P450 family 17 subfamily A member 1 (CYP17A1), were reduced in RPL, as verified by both qPCR and immunohistochemistry. The Gene Ontology (GO) analysis revealed the genes predominantly engaged in autophagy and iron metabolism. The area under the curve (AUC) demonstrated better diagnostic performance for RPL using CISD2 and CYP17A1. The single-cell transcriptomic analysis of RPL demonstrated that CISD2 is expressed in the majority of cell subpopulations, whereas CYP17A1 is not. The cell cycle analysis revealed highly active natural killer (NK) cells that displayed the highest communications with other cells, including the strongest interaction with macrophages through the migratory inhibitory factor (MIF) pathway.

CONCLUSIONS

Our study suggested that CISD2 and CYP17A1 genes are involved in abnormal iron metabolism, thereby contributing to RPL. These genes could be used as potential diagnostic and therapeutic markers for RPL.

Three-dimensional assessments are necessary to determine the true, spatially-resolved composition of tissues

Methods for spatially resolved cellular profiling using thinly cut sections have enabled in-depth quantitative tissue mapping to study inter-sample and intra-sample differences in normal human anatomy and disease onset and progression. These methods often profile extremely limited regions, which may impact the evaluation of heterogeneity due to tissue sub-sampling. Here, we applied CODA, a deep learning-based tissue mapping platform, to reconstruct the three-dimensional (3D) microanatomy of grossly normal and cancer-containing human pancreas biospecimens obtained from individuals who underwent pancreatic resection. To compare inter- and intra-sample heterogeneity, we assessed bulk and spatially resolved tissue composition in a cohort of two-dimensional (2D) whole slide images (WSIs) and a cohort of thick slabs of pancreas tissue that were digitally reconstructed in 3D from serial sections. To demonstrate the marked under sampling of 2D assessments, we simulated the number of WSIs and tissue microarrays (TMAs) necessary to represent the compositional heterogeneity of 3D data within 10% error to reveal that tens of WSIs and hundreds of TMA cores are sometimes needed. We show that spatial correlation of different pancreatic structures decay significantly within a span of microns, demonstrating that 2D histological sections may not be representative of their neighboring tissues. In sum, we demonstrate that 3D assessments are necessary to accurately assess tissue composition in normal and abnormal specimens and in order to accurately determine neoplastic content. These results emphasize the importance of intra-sample heterogeneity in tissue mapping efforts.

Improving drug response prediction via integrating gene relationships with deep learning

Predicting the drug response of cancer cell lines is crucial for advancing personalized cancer treatment, yet remains challenging due to tumor heterogeneity and individual diversity. In this study, we present a deep learning-based framework named Deep neural network Integrating Prior Knowledge (DIPK) (DIPK), which adopts self-supervised techniques to integrate multiple valuable information, including gene interaction relationships, gene expression profiles and molecular topologies, to enhance prediction accuracy and robustness. We demonstrated the superior performance of DIPK compared to existing methods on both known and novel cells and drugs, underscoring the importance of gene interaction relationships in drug response prediction. In addition, DIPK extends its applicability to single-cell RNA sequencing data, showcasing its capability for single-cell-level response prediction and cell identification. Further, we assess the applicability of DIPK on clinical data. DIPK accurately predicted a higher response to paclitaxel in the pathological complete response (pCR) group compared to the residual disease group, affirming the better response of the pCR group to the chemotherapy compound. We believe that the integration of DIPK into clinical decision-making processes has the potential to enhance individualized treatment strategies for cancer patients.

Unlocking melanoma Suppression: Insights from Plasma-Induced potent miRNAs through PI3K-AKT-ZEB1 axis

INTRODUCTION

Melanoma is a rare but highly malignant form of skin cancer. Although recent targeted and immune-based therapies have improved survival rates by 10-15%, effective melanoma treatment remains challenging. Therefore, novel, combinatorial therapy options such as non-thermal atmospheric pressure plasma (NTP) are being investigated to inhibit and prevent chemoresistance. Although several studies have reported the apoptotic and inhibitory effects of reactive oxygen species produced by NTP in the context of melanoma, the intricate molecular network that determines the role of microRNAs (miRNAs) in regulating NTP-mediated cell death remains unexplored.

OBJECTIVES

This study aimed to explore the molecular mechanisms and miRNA networks regulated by NTP-induced oxidative stress in melanoma cells.

METHODS

Melanoma cells were exposed to NTP and then subjected to high-throughput miRNA sequencing to identify NTP-regulated miRNAs. Various biological processes and underlying molecular mechanisms were assessed using Alamar Blue, propidium iodide (PI) uptake, cell migration, and clonogenic assays followed by qRT-PCR and flow cytometry.

RESULTS

NTP exposure for 3 min was sufficient to modulate the expression of several miRNAs, inhibiting cell growth. Persistent NTP exposure for 5 min increased differential miRNA regulation, PI uptake, and the expression of genes involved in cell cycle arrest and death. qPCR confirmed that miR-200b-3p and miR-215-5p upregulation contributed to decreased cell viability and migration. Mechanistically, inhibiting miR-200b-3p and miR-215-5p in SK-2 cells enhancedZEB1, PI3K, and AKT expression, increasing cell proliferation and viability.

CONCLUSION

This study demonstrated that NTP exposure for 5 min results in the differential regulation of miRNAs related to the PI3K-AKT-ZEB1 axis and cell cycle dysregulation to facilitate melanoma suppression.

GERWR: Identifying the Key Pathogenicity- Associated sRNAs of Magnaporthe Oryzae Infection in Rice Based on Graph Embedding and Random Walk With Restart

Rice blast, caused by Magnaporthe oryzae(M.oryzae), is a destructive rice disease that reduces rice yield by 10% to 30% annually. It also affects other cereal crops such as barley, wheat, rye, millet, sorghum, and maize. Small RNAs (sRNAs) play an essential regulatory role in fungus-plant interaction during the fungal invasion, but studies on pathogenic sRNAs during the fungal invasion of plants based on multi-omics data integration are rare. This paper proposes a novel approach called Graph Embedding combined with Random Walk with Restart (GERWR) to identify pathogenic sRNAs based on multi-omics data integration during M.oryzae invasion. By constructing a multi-omics network (MRMO), we identified 29 pathogenic sRNAs of rice blast fungus. Further analysis revealed that these sRNAs regulate rice genes in a many-to-many relationship, playing a significant regulatory role in the pathogenesis of rice blast disease. This paper explores the pathogenic factors of rice blast disease from the perspective of multi-omics data analysis, revealing the inherent connection between pathogenic factors of different omics. It has essential scientific significance for studying the pathogenic mechanism of rice blast fungus, the rice blast fungus-rice model system, and the pathogen-host interaction in related fields.

Non-severe thermal burn injuries induce long-lasting downregulation of gene expression in cortical excitatory neurons and microglia

Burn injuries are devastating traumas, often leading to life-long consequences that extend beyond the observable burn scar. In the context of the nervous system, burn injury patients commonly develop chronic neurological disorders and have been suggested to have impaired motor cortex function, but the long-lasting impact on neurons and glia in the brain is unknown. Using a mouse model of non-severe burn injury, excitatory and inhibitory neurons in the primary motor cortex were labelled with fluorescent proteins using adeno-associated viruses (AAVs). A total of 5 weeks following the burn injury, virus labelled excitatory and inhibitory neurons were isolated using fluorescence-activated cell sorting (FACS). In addition, microglia and astrocytes from the remaining cortical tissue caudal to the motor cortex were immunolabelled and isolated with FACS. Whole transcriptome RNA-sequencing was used to identify any long-lasting changes to gene expression in the different cell types. RNA-seq analysis showed changes to the expression of a small number of genes with known functions in excitatory neurons and microglia, but not in inhibitory neurons or astrocytes. Specifically, genes related to GABA-A receptors in excitatory neurons and several cellular functions in microglia were found to be downregulated in burn injured mice. These findings suggest that non-severe burn injuries lead to long lasting transcriptomic changes in the brain, but only in specific cell types. Our findings provide a broad overview of the long-lasting impact of burn injuries on the central nervous system which may help identify potential therapeutic targets to prevent neurological dysfunction in burn patients.

Endophytic fungi as a potential source of anti-cancer drug

Endophytes are considered one of the major sources of bioactive compounds used in different aspects of health care including cancer treatment. When colonized, they either synthesize these bioactive compounds as a part of their secondary metabolite production or augment the host plant machinery in synthesising such bioactive compounds. Hence, the study of endophytes has drawn the attention of the scientific community in the last few decades. Among the endophytes, endophytic fungi constitute a major portion of endophytic microbiota. This review deals with a plethora of anti-cancer compounds derived from endophytic fungi, highlighting alkaloids, lignans, terpenes, polyketides, polyphenols, quinones, xanthenes, tetralones, peptides, and spirobisnaphthalenes. Further, this review emphasizes modern methodologies, particularly omics-based techniques, asymmetric dihydroxylation, and biotic elicitors, showcasing the dynamic and evolving landscape of research in this field and describing the potential of endophytic fungi as a source of anticancer drugs in the future.

Comparison of RNA-Seq and microarray in the prediction of protein expression and survival prediction

Gene expression profiling using RNA-sequencing (RNA-seq) and microarray technologies is widely used in cancer research to identify biomarkers for clinical endpoint prediction. We compared the performance of these two methods in predicting protein expression and clinical endpoints using The Cancer Genome Atlas (TCGA) datasets of lung cancer, colorectal cancer, renal cancer, breast cancer, endometrial cancer, and ovarian cancer. We calculated the correlation coefficients between gene expression measured by RNA-seq or microarray and protein expression measured by reverse phase protein array (RPPA). In addition, after selecting the top 103 survival-related genes, we compared the random forest survival prediction model performance across test platforms and cancer types. Both RNA-seq and microarray data were retrieved from TCGA dataset. Most genes showed similar correlation coefficients between RNA-seq and microarray, but 16 genes exhibited significant differences between the two methods. The BAX gene was recurrently found in colorectal cancer, renal cancer, and ovarian cancer, and the PIK3CA gene belonged to renal cancer and breast cancer. Furthermore, the survival prediction model using microarray was better than the RNA-seq model in colorectal cancer, renal cancer, and lung cancer, but the RNA-seq model was better in ovarian and endometrial cancer. Our results showed good correlation between mRNA levels and protein measured by RPPA. While RNA-seq and microarray performance were similar, some genes showed differences, and further clinical significance should be evaluated. Additionally, our survival prediction model results were controversial.

GAM-MDR: probing miRNA-drug resistance using a graph autoencoder based on random path masking

MicroRNAs (miRNAs) are found ubiquitously in biological cells and play a pivotal role in regulating the expression of numerous target genes. Therapies centered around miRNAs are emerging as a promising strategy for disease treatment, aiming to intervene in disease progression by modulating abnormal miRNA expressions. The accurate prediction of miRNA-drug resistance (MDR) is crucial for the success of miRNA therapies. Computational models based on deep learning have demonstrated exceptional performance in predicting potential MDRs. However, their effectiveness can be compromised by errors in the data acquisition process, leading to inaccurate node representations. To address this challenge, we introduce the GAM-MDR model, which combines the graph autoencoder (GAE) with random path masking techniques to precisely predict potential MDRs. The reliability and effectiveness of the GAM-MDR model are mainly reflected in two aspects. Firstly, it efficiently extracts the representations of miRNA and drug nodes in the miRNA-drug network. Secondly, our designed random path masking strategy efficiently reconstructs critical paths in the network, thereby reducing the adverse impact of noisy data. To our knowledge, this is the first time that a random path masking strategy has been integrated into a GAE to infer MDRs. Our method was subjected to multiple validations on public datasets and yielded promising results. We are optimistic that our model could offer valuable insights for miRNA therapeutic strategies and deepen the understanding of the regulatory mechanisms of miRNAs. Our data and code are publicly available at GitHub:https://github.com/ZZCrazy00/GAM-MDR.

Classifying breast cancer using multi-view graph neural network based on multi-omics data

Introduction: As the evaluation indices, cancer grading and subtyping have diverse clinical, pathological, and molecular characteristics with prognostic and therapeutic implications. Although researchers have begun to study cancer differentiation and subtype prediction, most of relevant methods are based on traditional machine learning and rely on single omics data. It is necessary to explore a deep learning algorithm that integrates multi-omics data to achieve classification prediction of cancer differentiation and subtypes. Methods: This paper proposes a multi-omics data fusion algorithm based on a multi-view graph neural network (MVGNN) for predicting cancer differentiation and subtype classification. The model framework consists of a graph convolutional network (GCN) module for learning features from different omics data and an attention module for integrating multi-omics data. Three different types of omics data are used. For each type of omics data, feature selection is performed using methods such as the chi-square test and minimum redundancy maximum relevance (mRMR). Weighted patient similarity networks are constructed based on the selected omics features, and GCN is trained using omics features and corresponding similarity networks. Finally, an attention module integrates different types of omics features and performs the final cancer classification prediction. Results: To validate the cancer classification predictive performance of the MVGNN model, we conducted experimental comparisons with traditional machine learning models and currently popular methods based on integrating multi-omics data using 5-fold cross-validation. Additionally, we performed comparative experiments on cancer differentiation and its subtypes based on single omics data, two omics data, and three omics data. Discussion: This paper proposed the MVGNN model and it performed well in cancer classification prediction based on multiple omics data.

Peripheral inflammation-induced changes in songbird brain gene expression: 3' mRNA transcriptomic approach

Species-specific neural inflammation can be induced by profound immune signalling from periphery to brain. Recent advances in transcriptomics offer cost-effective approaches to study this regulation. In a population of captive zebra finch (Taeniopygia guttata), we compare the differential gene expression patterns in lipopolysaccharide (LPS)-triggered peripheral inflammation revealed by RNA-seq and QuantSeq. The RNA-seq approach identified more differentially expressed genes but failed to detect any inflammatory markers. In contrast, QuantSeq results identified specific expression changes in the genes regulating inflammation. Next, we adopted QuantSeq to relate peripheral and brain transcriptomes. We identified subtle changes in the brain gene expression during the peripheral inflammation (e.g. up-regulation in AVD-like and ACOD1 expression) and detected co-structure between the peripheral and brain inflammation. Our results suggest benefits of the 3'end transcriptomics for association studies between peripheral and neural inflammation in genetically heterogeneous models and identify potential targets for the future brain research in birds.

A Circular RNA Derived from the Pumilio 1 Gene Could Regulate PTEN in Human Cumulus Cells

CircRNAs are a class of non-coding RNAs able to regulate gene expression at multiple levels. Their involvement in physiological processes, as well as their altered regulation in different human diseases, both tumoral and non-tumoral, is well documented. However, little is known about their involvement in female reproduction. This study aims to identify circRNAs potentially involved in reproductive women's health. Candidate circRNAs expressed in ovary and sponging miRNAs, already known to be expressed in the ovary, were selected by a computational approach. Using real time PCR, we verified their expression and identified circPUM1 as the most interesting candidate circRNA for further analyses. We assessed the expression of circPUM1 and its linear counterpart in all the follicle compartments and, using a computational and experimental approach, identified circPUM1 direct and indirect targets, miRNAs and mRNAs, respectively, in cumulus cells. We found that both circPUM1 and its mRNA host gene are co-expressed in all the follicle compartments and proposed circPUM1 as a potential regulator of PTEN, finding a strong positive correlation between circPUM1 and PTEN mRNA. These results suggest a possible regulation of PTEN by circPUM1 in cumulus cells and point out the important role of circRNA inside the pathways related to follicle growth and oocyte maturation.

Uncovering transcriptomic biomarkers for enhanced diagnosis of methamphetamine use disorder: a comprehensive review

Introduction

Methamphetamine use disorder (MUD) is a chronic relapsing disorder characterized by compulsive Methamphetamine (MA) use despite its detrimental effects on physical, psychological, and social well-being. The development of MUD is a complex process that involves the interplay of genetic, epigenetic, and environmental factors. The treatment of MUD remains a significant challenge, with no FDA-approved pharmacotherapies currently available. Current diagnostic criteria for MUD rely primarily on self-reporting and behavioral assessments, which have inherent limitations owing to their subjective nature. This lack of objective biomarkers and unidimensional approaches may not fully capture the unique features and consequences of MA addiction.

Methods

We performed a literature search for this review using the Boolean search in the PubMed database.

Results

This review explores existing technologies for identifying transcriptomic biomarkers for MUD diagnosis. We examined non-invasive tissues and scrutinized transcriptomic biomarkers relevant to MUD. Additionally, we investigated transcriptomic biomarkers identified for diagnosing, predicting, and monitoring MUD in non-invasive tissues.

Discussion

Developing and validating non-invasive MUD biomarkers could address these limitations, foster more precise and reliable diagnostic approaches, and ultimately enhance the quality of care for individuals with MA addiction.

Transcriptomic and Metabolomic Analysis of Liver Cirrhosis

BACKGROUND

Liver cirrhosis is one of the leading causes of decreased life expectancy worldwide. However, the molecular mechanisms underlying liver cirrhosis remain unclear. In this study, we performed a comprehensive analysis using transcriptome and metabolome sequencing to explore the genes, pathways, and interactions associated with liver cirrhosis.

METHODS

We performed transcriptome and metabolome sequencing of blood samples from patients with cirrhosis and healthy controls (1:1 matched for sex and age). We validated the differentially expressed microRNA (miRNA) and mRNAs using real-time quantitative polymerase chain reaction.

RESULTS

For transcriptome analysis, we screened for differentially expressed miRNAs and mRNAs, analyzed mRNAs to identify possible core genes and pathways, and performed coanalysis of miRNA and mRNA sequencing results. In terms of the metabolome, we screened five pathways that were substantially enriched in the differential metabolites. Next, we identified the metabolites with the most pronounced differences among these five metabolic pathways. We performed receiver operating characteristic (ROC) curve analysis of these five metabolites to determine their diagnostic efficacy for cirrhosis. Finally, we explored possible links between the transcriptome and metabolome.

CONCLUSION

Based on sequencing and bioinformatics, we identified miRNAs and genes that were differentially expressed in the blood of patients with liver cirrhosis. By exploring pathways and disease-specific networks, we identified unique biological mechanisms. In terms of metabolomes, we identified novel biomarkers and explored their diagnostic efficacy. We identified possible common pathways in the transcriptome and metabolome that could serve as candidates for further studies.

Analyzing RNA posttranscriptional modifications to decipher the epitranscriptomic code

The discovery of RNA silencing has revealed that non-protein-coding sequences (ncRNAs) can cover essential roles in regulatory networks and their malfunction may result in severe consequences on human health. These findings have prompted a general reassessment of the significance of RNA as a key player in cellular processes. This reassessment, however, will not be complete without a greater understanding of the distribution and function of the over 170 variants of the canonical ribonucleotides, which contribute to the breathtaking structural diversity of natural RNA. This review surveys the analytical approaches employed for the identification, characterization, and detection of RNA posttranscriptional modifications (rPTMs). The merits of analyzing individual units after exhaustive hydrolysis of the initial biopolymer are outlined together with those of identifying their position in the sequence of parent strands. Approaches based on next generation sequencing and mass spectrometry technologies are covered in depth to provide a comprehensive view of their respective merits. Deciphering the epitranscriptomic code will require not only mapping the location of rPTMs in the various classes of RNAs, but also assessing the variations of expression levels under different experimental conditions. The fact that no individual platform is currently capable of meeting all such demands implies that it will be essential to capitalize on complementary approaches to obtain the desired information. For this reason, the review strived to cover the broadest possible range of techniques to provide readers with the fundamental elements necessary to make informed choices and design the most effective possible strategy to accomplish the task at hand.

From cancer big data to treatment: Artificial intelligence in cancer research

In recent years, developing the idea of "cancer big data" has emerged as a result of the significant expansion of various fields such as clinical research, genomics, proteomics and public health records. Advances in omics technologies are making a significant contribution to cancer big data in biomedicine and disease diagnosis. The increasingly availability of extensive cancer big data has set the stage for the development of multimodal artificial intelligence (AI) frameworks. These frameworks aim to analyze high-dimensional multi-omics data, extracting meaningful information that is challenging to obtain manually. Although interpretability and data quality remain critical challenges, these methods hold great promise for advancing our understanding of cancer biology and improving patient care and clinical outcomes. Here, we provide an overview of cancer big data and explore the applications of both traditional machine learning and deep learning approaches in cancer genomic and proteomic studies. We briefly discuss the challenges and potential of AI techniques in the integrated analysis of omics data, as well as the future direction of personalized treatment options in cancer.

Revolutionizing Malaria Vector Control: The Importance of Accurate Species Identification through Enhanced Molecular Capacity

Many factors, such as the resistance to pesticides and a lack of knowledge of the morphology and molecular structure of malaria vectors, have made it more challenging to eradicate malaria in numerous malaria-endemic areas of the globe. The primary goal of this review is to discuss malaria vector control methods and the significance of identifying species in vector control initiatives. This was accomplished by reviewing methods of molecular identification of malaria vectors and genetic marker classification in relation to their use for species identification. Due to its specificity and consistency, molecular identification is preferred over morphological identification of malaria vectors. Enhanced molecular capacity for species identification will improve mosquito characterization, leading to accurate control strategies/treatment targeting specific mosquito species, and thus will contribute to malaria eradication. It is crucial for disease epidemiology and surveillance to accurately identify the Plasmodium spp. that are causing malaria in patients. The capacity for disease surveillance will be significantly increased by the development of more accurate, precise, automated, and high-throughput diagnostic techniques. In conclusion, although morphological identification is quick and achievable at a reduced cost, molecular identification is preferred for specificity and sensitivity. To achieve the targeted malaria elimination goal, proper identification of vectors using accurate techniques for effective control measures should be prioritized.

Integrating genomics and genome editing for orphan crop improvement: a bridge between orphan crops and modern agriculture system

Domestication of orphan crops could be explored by editing their genomes. Genome editing has a lot of promise for enhancing agricultural output, and there is a lot of interest in furthering breeding in orphan crops, which are sometimes plagued with unwanted traits that resemble wild cousins. Consequently, applying model crop knowledge to orphan crops allows for the rapid generation of targeted allelic diversity and innovative breeding germplasm. We explain how plant breeders could employ genome editing as a novel platform to accelerate the domestication of semi-domesticated or wild plants, resulting in a more diversified base for future food and fodder supplies. This review emphasizes both the practicality of the strategy and the need to invest in research that advances our understanding of plant genomes, genes, and cellular systems. Planting more of these abandoned orphan crops could help alleviate food scarcities in the challenge of future climate crises.

Ambiguous genes due to aligners and their impact on RNA-seq data analysis

The main scope of the study is ambiguous genes, i.e. genes whose expression is difficult to estimate from the data produced by next-generation sequencing technologies. We focused on the RNA sequencing (RNA-Seq) type of experiment performed on the Illumina platform. It is crucial to identify such genes and understand the cause of their difficulty, as these genes may be involved in some diseases. By giving misleading results, they could contribute to a misunderstanding of the cause of certain diseases, which could lead to inappropriate treatment. We thought that the ambiguous genes would be difficult to map because of their complex structure. So we looked at RNA-seq analysis using different mappers to find genes that would have different measurements from the aligners. We were able to identify such genes using a generalized linear model with two factors: mappers and groups introduced by the experiment. A large proportion of ambiguous genes are pseudogenes. High sequence similarity of pseudogenes to functional genes may indicate problems in alignment procedures. In addition, predictive analysis verified the performance of difficult genes in classification. The effectiveness of classifying samples into specific groups was compared, including the expression of difficult and not difficult genes as covariates. In almost all cases considered, ambiguous genes have less predictive power.

Hydrated lime promoted the polysaccharide content and affected the transcriptomes of Lentinula edodes during brown film formation

Brown film formation, a unique developmental stage in the life cycle of Lentinula edodes, is essential for the subsequent development of fruiting bodies in L. edodes cultivation. The pH of mushroom growth substrates are usually adjusted with hydrated lime, yet the effects of hydrated lime on cultivating L. edodes and the molecular mechanisms associated with the effects have not been studied systemically. We cultivated L. edodes on substrates supplemented with 0% (CK), 1% (T1), 3% (T2), and 5% (T3) hydrated lime (Ca (OH)2), and applied transcriptomics and qRT-PCR to study gene expression on the brown film formation stage. Hydrated lime increased polysaccharide contents in L. edodes, especially in T2, where the 5.3% polysaccharide content was approximately 1.5 times higher than in the CK. The addition of hydrated lime in the substrate promoted laccase, lignin peroxidase and manganese peroxidase activities, implying that hydrated lime improved the ability of L. edodes to decompose lignin and provide nutrition for its growth and development. Among the annotated 9,913 genes, compared to the control, 47 genes were up-regulated and 52 genes down-regulated in T1; 73 genes were up-regulated and 44 were down-regulated in T2; and 125 genes were up-regulated and 65 genes were down-regulated in T3. Differentially expressed genes (DEGs) were enriched in the amino acid metabolism, lipid metabolism and carbohydrate metabolism related pathways. The carbohydrate-active enzyme genes up-regulated in the hydrated lime treatments were mostly glycosyl hydrolase genes. The results will facilitate future optimization of L. edodes cultivation techniques and possibly shortening the production cycle.

Virulence factors of the genus Fusarium with targets in plants

Fusarium spp. comprise various species of filamentous fungi that cause severe diseases in plant crops of both agricultural and forestry interest. These plant pathogens produce a wide range of molecules with diverse chemical structures and biological activities. Genetic functional analyses of some of these compounds have shown their role as virulence factors (VF). However, their mode of action and contributions to the infection process for many of these molecules are still unknown. This review aims to analyze the state of the art in Fusarium VF, emphasizing their biological targets on the plant hosts. It also addresses the current experimental approaches to improve our understanding of their role in virulence and suggests relevant research questions that remain to be answered with a greater focus on species of agroeconomic importance. In this review, a total of 37 confirmed VF are described, including 22 proteinaceous and 15 non-proteinaceous molecules, mainly from Fusarium oxysporum and Fusarium graminearum and, to a lesser extent, in Fusarium verticillioides and Fusarium solani.

Transcriptome profiles of Streptomyces clavuligerus strains producing different titers of clavulanic acid

Streptomyces clavuligerus NRRL 3585 is a native producer of clavulanic acid (CA), a clinically used β-lactamase inhibitor, and is widely used as an industrial strain for the production of antibiotics. Selective random mutagenesis has successfully generated the improved CA-producing S. clavuligerus mutant strains as well as the strain with the loss of CA biosynthesis. To understand the molecular mechanisms associated with the improved CA-production potential, genome-scale RNA-sequencing-based transcriptional data were obtained for the wild-type S. clavuligerus strain and its three mutant strains. Total RNA samples for each strain were collected across four different growth stages, and all 32 sequencing data points exhibited an average Phred score of 36. The high-quality genome-scale transcriptional profile of S. clavuligerus strains with varied CA biosynthetic potential provides valuable insights and new opportunities for discovering efficient metabolic engineering strategies for the development of improved industrial strains.

Radical-Induced Dissociation for Oligonucleotide Sequencing by TiO2/ZnAl-Layered Double Oxide-Assisted Laser Desorption/Ionization Mass Spectrometry

De novo sequencing of oligonucleotides remains challenging, especially for oligonucleotides with post-transcriptional or synthetic modifications. Mass spectrometry (MS) sequencing can reliably detect and locate all of the modification sites in oligonucleotides via m/z variance. However, current MS-based sequencing methods exhibit complex spectra and low ion abundance and usually require coupled instrumentation. Herein, we demonstrate a method of oligonucleotide sequencing using TiO2/ZnAl-layered double oxide (LDO)-assisted laser desorption/ionization (LDI)-MS based on radical-induced dissociation (RID). ·CH2OH radicals can be produced on the surface of a TiO2/ZnAl-LDO matrix via ultraviolet light, inducing an attack on the active site of the oligonucleotide phosphate skeleton to create typical "a-, a-B-, c·-, d-, w-, and y"-type fragments. Compared with the spectra obtained via collision-based methods, such as collision-induced dissociation and higher-energy collisional dissociation, the LDI-MS spectra based on RID exhibit single-charged signals, fewer types of fragments, and a lower proportion of unknown noise peaks. We demonstrate full sequence coverage for a 6-mer 2'-O-methyl-modified oligonucleotide and a 21-mer small interfering RNA and show that RID can sequence oligonucleotides with modifications. Importantly, the mechanism responsible for the RID of the oligonucleotide phosphate skeleton was investigated through offline experiments, demonstrating consistent results with density functional theory calculations.

O. Improved meta-analysis pipeline ameliorates distinctive gene regulators of diabetic vasculopathy in human endothelial cell (hECs) RNA-Seq data

Enormous gene expression data generated through next-generation sequencing (NGS) technologies are accessible to the scientific community via public repositories. The data harboured in these repositories are foundational for data integrative studies enabling large-scale data analysis whose potential is yet to be fully realized. Prudent integration of individual gene expression data i.e. RNA-Seq datasets is remarkably challenging as it encompasses an assortment and series of data analysis steps that requires to be accomplished before arriving at meaningful insights on biological interrogations. These insights are at all times latent within the data and are not usually revealed from the modest individual data analysis owing to the limited number of biological samples in individual studies. Nevertheless, a sensibly designed meta-analysis of select individual studies would not only maximize the sample size of the analysis but also significantly improves the statistical power of analysis thereby revealing the latent insights. In the present study, a custom-built meta-analysis pipeline is presented for the integration of multiple datasets from different origins. As a case study, we have tested with the integration of two relevant datasets pertaining to diabetic vasculopathy retrieved from the open source domain. We report the meta-analysis ameliorated distinctive and latent gene regulators of diabetic vasculopathy and uncovered a total of 975 i.e. 930 up-regulated and 45 down-regulated gene signatures. Further investigation revealed a subset of 14 DEGs including CTLA4, CALR, G0S2, CALCR, OMA1, and DNAJC3 as latent i.e. novel as these signatures have not been reported earlier. Moreover, downstream investigations including enrichment analysis, and protein-protein interaction (PPI) network analysis of DEGs revealed durable disease association signifying their potential as novel transcriptomic biomarkers of diabetic vasculopathy. While the meta-analysis of individual whole transcriptomic datasets for diabetic vasculopathy is exclusive to our comprehension, however, the novel meta-analysis pipeline could very well be extended to study the mechanistic links of DEGs in other disease conditions.

Nitrate improves aluminium resistance through SLAH-mediated citrate exudation from roots

Aluminium (Al) toxicity is one of the major constraint for crop production in acidic soil, and the inappropriate utilization of nitrogen fertilizer can accelerate soil acidification. Despite previous studies investigating the regulation of nitrogen forms in Al toxicity of plants, the underlying mechanism, particularly at the molecular level, remains unclear. This study aims to uncover the potentially regulatory mechanism of nitrate (NO3 - ) in the Al resistance of maize and Arabidopsis. NO3 - conservatively improves Al resistance in maize and Arabidopsis, with nitrate-elevated citrate synthesis and exudation potentially playing critical roles in excluding Al from the root symplast. ZmSLAH2 in maize and AtSLAH1 in Arabidopsis are essential for the regulation of citrate exudation and NO3 - -promoted Al resistance, with ZmMYB81 directly targeting the ZmSLAH2 promoter to activate its activity. Additionally, NO3 - transport is necessary for NO3 - -promoted Al resistance, with ZmNRT1.1A and AtNRT1.1 potentially playing vital roles. The suppression of NO3 - transport in roots by ammonium (NH4 + ) may inhibit NO3 - -promoted Al resistance. This study provides novel insights into the understanding of the crucial role of NO3 - -mediated signalling in the Al resistance of plants and offers guidance for nitrogen fertilization on acid soils.

Antifungal mechanism of volatile compounds emitted by Actinomycetota Paenarthrobacter ureafaciens from a disease-suppressive soil on Saccharomyces cerevisiae

Increasing evidence suggests that in disease-suppressive soils, microbial volatile compounds (mVCs) released from bacteria may inhibit the growth of plant-pathogenic fungi. However, the antifungal activities and molecular responses of fungi to different mVCs remain largely undescribed. In this study, we first evaluated the responses of pathogenic fungi to treatment with mVCs from Paenarthrobacter ureafaciens. Then, we utilized the well-characterized fungal model organism Saccharomyces cerevisiae to study the potential mechanistic effects of the mVCs. Our data showed that exposure to P. ureafaciens mVCs leads to reduced growth of several pathogenic fungi, and in yeast cells, mVC exposure prompts the accumulation of reactive oxygen species. Further experiments with S. cerevisiae deletion mutants indicated that Slt2/Mpk1 and Hog1 MAPKs play major roles in the yeast response to P. ureafaciens mVCs. Transcriptomic analysis revealed that exposure to mVCs was associated with 1,030 differentially expressed genes (DEGs) in yeast. According to gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses, many of these DEGs are involved in mitochondrial dysfunction, cell integrity, mitophagy, cellular metabolism, and iron uptake. Genes encoding antimicrobial proteins were also significantly altered in the yeast after exposure to mVCs. These findings suggest that oxidative damage and mitochondrial dysfunction are major contributors to the fungal toxicity of mVCs. Furthermore, our data showed that cell wall, antioxidant, and antimicrobial defenses are induced in yeast exposed to mVCs. Thus, our findings expand upon previous research by delineating the transcriptional responses of the fungal model. IMPORTANCE Since the use of bacteria-emitted volatile compounds in phytopathogen control is of considerable interest, it is important to understand the molecular mechanisms by which fungi may adapt to microbial volatile compounds (mVCs). Paenarthrobacter ureafaciens is an isolated bacterium from disease-suppressive soil that belongs to the Actinomycetota phylum. P. ureafaciens mVCs showed a potent antifungal effect on phytopathogens, which may contribute to disease suppression in soil. However, our knowledge about the antifungal mechanism of mVCs is limited. This study has proven that mVCs are toxic to fungi due to oxidative stress and mitochondrial dysfunction. To deal with mVC toxicity, antioxidants and physical defenses are required. Furthermore, iron uptake and CAP proteins are required for antimicrobial defense, which is necessary for fungi to deal with the thread from mVCs. This study provides essential foundational knowledge regarding the molecular responses of fungi to inhibitory mVCs.

New insights into transcriptome variation during cattle adipocyte adipogenesis by direct RNA sequencing

We performed direct RNA sequencing (DRS) together with PCR-amplified cDNA long and short read sequencing for cattle adipocyte at different stages. We proved that the DRS was with advantages to avoid artificial transcripts and questionable exitrons. Totally, we obtained 68,124 transcripts with information of alternative splicing, poly (A) length and mRNA modification. The number of transcripts for adipogenesis was expanded by alternative splicing, which lead regulation mechanisms far more complex than ever known. We detected 891 differentially expressed genes (DEGs). However, 62.78% transcripts of DEGs were not significantly differentially expressed, and 248 transcripts showed opposite changing directions with their genes. The poly (A) tail became globally shorter in differentiated adipocyte than in primary adipocyte, and had a weak negative correlation with gene/transcript expression. Moreover, the study of different mRNA modifications implied their potential roles in gene expression and alternative splicing. Overall, our study promoted better understanding of adipogenesis mechanisms in cattle adipocytes.

A urine extracellular vesicle lncRNA classifier for high-grade prostate cancer and increased risk of progression: A multi-center study

To construct a urine extracellular vesicle long non-coding RNA (lncRNA) classifier that can detect high-grade prostate cancer (PCa) of grade group 2 or greater and estimate the risk of progression during active surveillance, we identify high-grade PCa-specific lncRNAs by combined analyses of cohorts from TAHSY, TCGA, and the GEO database. We develop and validate a 3-lncRNA diagnostic model (Clnc, being made of AC015987.1, CTD-2589M5.4, RP11-363E6.3) that can detect high-grade PCa. Clnc shows higher accuracy than prostate cancer antigen 3 (PCA3), multiparametric magnetic resonance imaging (mpMRI), and two risk calculators (Prostate Cancer Prevention Trial [PCPT]-RC 2.0 and European Randomized Study of Screening for Prostate Cancer [ERSPC]-RC) in the training cohort (n = 350), two independent cohorts (n = 232; n = 251), and TCGA cohort (n = 499). In the prospective active surveillance cohort (n = 182), Clnc at diagnosis remains a powerful independent predictor for overall active surveillance progression. Thus, Clnc is a potential biomarker for high-grade PCa and can also serve as a biomarker for improved selection of candidates for active surveillance.

Harnessing large language models (LLMs) for candidate gene prioritization and selection

BACKGROUND

Feature selection is a critical step for translating advances afforded by systems-scale molecular profiling into actionable clinical insights. While data-driven methods are commonly utilized for selecting candidate genes, knowledge-driven methods must contend with the challenge of efficiently sifting through extensive volumes of biomedical information. This work aimed to assess the utility of large language models (LLMs) for knowledge-driven gene prioritization and selection.

METHODS

In this proof of concept, we focused on 11 blood transcriptional modules associated with an Erythroid cells signature. We evaluated four leading LLMs across multiple tasks. Next, we established a workflow leveraging LLMs. The steps consisted of: (1) Selecting one of the 11 modules; (2) Identifying functional convergences among constituent genes using the LLMs; (3) Scoring candidate genes across six criteria capturing the gene's biological and clinical relevance; (4) Prioritizing candidate genes and summarizing justifications; (5) Fact-checking justifications and identifying supporting references; (6) Selecting a top candidate gene based on validated scoring justifications; and (7) Factoring in transcriptome profiling data to finalize the selection of the top candidate gene.

RESULTS

Of the four LLMs evaluated, OpenAI's GPT-4 and Anthropic's Claude demonstrated the best performance and were chosen for the implementation of the candidate gene prioritization and selection workflow. This workflow was run in parallel for each of the 11 erythroid cell modules by participants in a data mining workshop. Module M9.2 served as an illustrative use case. The 30 candidate genes forming this module were assessed, and the top five scoring genes were identified as BCL2L1, ALAS2, SLC4A1, CA1, and FECH. Researchers carefully fact-checked the summarized scoring justifications, after which the LLMs were prompted to select a top candidate based on this information. GPT-4 initially chose BCL2L1, while Claude selected ALAS2. When transcriptional profiling data from three reference datasets were provided for additional context, GPT-4 revised its initial choice to ALAS2, whereas Claude reaffirmed its original selection for this module.

CONCLUSIONS

Taken together, our findings highlight the ability of LLMs to prioritize candidate genes with minimal human intervention. This suggests the potential of this technology to boost productivity, especially for tasks that require leveraging extensive biomedical knowledge.

Enhanced DNA Entropy-Driven Circuit by Locked Nucleic Acids and Simulation-Guided Localization

Signal amplification methods based on DNA molecular interactions are promising tools for detecting various biomarkers in low abundance. The entropy-driven circuit (EDC), as an enzyme-free signal amplification method, has been used in detecting and imaging a variety of biomarkers. The localization strategy can effectively increase the local concentration of the DNA reaction modules to improve the signal amplification effect. However, the localization strategy may also amplify the leak reaction of the EDC, and effective signal amplification can be limited by the unclear structure-function relationship. Herein, we utilized locked nucleic acid (LNA) modification to enhance the stability of the localized entropy-driven circuit (LEDC), which suppressed a 94.6% leak signal. The coarse-grained model molecular simulation was used to guide the structure design of the LEDC, and the influence of critical factors such as the localized distance and spacer length was analyzed at the molecular level to obtain the best reaction performance. The sensitivities of miR-21 and miR-141 detected by a simulation-guided optimal LEDC probe were 17.45 and 65 pM, 1345 and 521 times higher than free-EDC, respectively. The LEDC was further employed for the fluorescence imaging of miRNA in cancer cells, showing excellent specificity and sensitivity. This work utilizes LNA and molecular simulations to comprehensively improve the performance of a localized DNA signal amplification circuit, providing an advanced DNA probe design strategy for biosensing and imaging as well as valuable information for the designers of DNA-based probes.

Targeting Desulfovibrio vulgaris flagellin-induced NAIP/NLRC4 inflammasome activation in macrophages attenuates ulcerative colitis

INTRODUCTION

The perturbations of gut microbiota could interact with excessively activated immune responses and play key roles in the etiopathogenesis of ulcerative colitis (UC). Desulfovibrio, the most predominant sulfate reducing bacteria (SRB) resided in the human gut, was observed to overgrow in patients with UC. The interactions between specific gut microbiota and drugs and their impacts on UC treatment have not been demonstrated well.

OBJECTIVES

This study aimed to elucidate whether Desulfovibrio vulgaris (D. vulgaris, DSV) and its flagellin could activate nucleotide-binding oligomerization domain-like receptors (NLR) family of apoptosis inhibitory proteins (NAIP) / NLR family caspase activation and recruitment domain-containing protein 4 (NLRC4) inflammasome and promote colitis, and further evaluate the efficacy of eugeniin targeting the interaction interface of D. vulgaris flagellin (DVF) and NAIP to attenuate UC.

METHODS

The abundance of DSV and the occurrence of macrophage pyroptosis in human UC tissues were investigated. Colitis in mice was established by dextran sulfate sodium (DSS) and gavaged with DSV or its purified flagellin. NAIP/NLRC4 inflammasome activation and macrophage pyroptosis were evaluated in vivo and in vitro. The effects of eugeniin on blocking the interaction of DVF and NAIP/NLRC4 and relieving colitis were also assessed.

RESULTS

The abundance of DSV increased in the feces of patients with UC and was found to be associated with disease activity. DSV and its flagellin facilitated DSS-induced colitis in mice. Mechanistically, RNA sequencing showed that gene expression associated with inflammasome complex and pyroptosis was upregulated after DVF treatment in macrophages. DVF was further demonstrated to induce significant macrophage pyroptosis in vitro, depending on NAIP/NLRC4 inflammasome activation. Furthermore, eugeniin was screened as an inhibitor of the interface between DVF and NAIP and successfully alleviated the proinflammatory effect of DVF in colitis.

CONCLUSION

Targeting DVF-induced NAIP/NLRC4 inflammasome activation and macrophage pyroptosis ameliorates UC. This finding is of great significance for exploring the gut microbiota-host interactions in UC development and providing new insights for precise treatment.

An Integrative ATAC-Seq and RNA-Seq Analysis of the Endometrial Tissues of Meishan and Duroc Pigs

Meishan pigs are a well-known indigenous pig breed in China characterized by a high fertility. Notably, the number of endometrial grands is significantly higher in Meishan pigs than Duroc pigs. The characteristics of the endometrial tissue are related to litter size. Therefore, we used the assay for transposase-accessible chromatin with sequencing (ATAC-seq) and RNA-sequencing (RNA-seq) to analyze the mechanisms underlying the differences in fecundity between the breeds. We detected the key transcription factors, including Double homeobox (Dux), Ladybird-like homeobox gene 2 (LBX2), and LIM homeobox 8 (Lhx8), with potentially pivotal roles in the regulation of the genes related to endometrial development. We identified the differentially expressed genes between the breeds, including SOX17, ANXA4, DLX3, DMRT1, FLNB, IRF6, CBFA2T2, TFCP2L1, EFNA5, SLIT2, and CYFIP2, with roles in epithelial cell differentiation, fertility, and ovulation. Interestingly, ANXA4, CBFA2T2, and TFCP2L1, which were upregulated in the Meishan pigs in the RNA-seq analysis, were identified again by the integration of the ATAC-seq and RNA-seq data. Moreover, we identified genes in the cancer or immune pathways, FoxO signaling, Wnt signaling, and phospholipase D signaling pathways. These ATAC-seq and RNA-seq analyses revealed the accessible chromatin and potential mechanisms underlying the differences in the endometrial tissues between the two types of pigs.

Early Pharmacodynamic Changes Measured Using RNA Sequencing of Peripheral Blood from Patients in a Phase I Study with Mitazalimab, a Potent CD40 Agonistic Monoclonal Antibody

CD40-targeting therapies can enhance the dendritic cell priming of tumor-specific T cells and repolarize intratumoral macrophages to alleviate the tumoral immunosuppressive environment and remodel the extracellular matrix. Mitazalimab is a potent agonistic CD40 monoclonal IgG1 antibody currently under clinical development. This study used RNA sequencing of blood samples from a subset of patients from a Phase I trial with mitazalimab (NCT02829099) to assess peripheral pharmacodynamic activity. We found that mitazalimab induced transient peripheral transcriptomic alterations (at 600 µg/kg and 900 µg/kg dose administered intravenously), which mainly were attributed to immune activation. In particular, the transcriptomic alterations showed a reduction in effector cells (e.g., CD8+ T cells and natural killer cells) and B cells peripherally with the remaining cells (e.g., dendritic cells, monocytes, B cells, and natural killer cells) showing transcription profiles consistent with activation. Lastly, distinct patient subgroups based on the pattern of transcriptomic alterations could be identified. In summary, the data presented herein reinforce the anticipated mode of action of mitazalimab and support its ongoing clinical development.

Reproducible Bioinformatics Analysis Workflows for Detecting IGH Gene Fusions in B-Cell Acute Lymphoblastic Leukaemia Patients

B-cell acute lymphoblastic leukaemia (B-ALL) is characterised by diverse genomic alterations, the most frequent being gene fusions detected via transcriptomic analysis (mRNA-seq). Due to its hypervariable nature, gene fusions involving the Immunoglobulin Heavy Chain (IGH) locus can be difficult to detect with standard gene fusion calling algorithms and significant computational resources and analysis times are required. We aimed to optimize a gene fusion calling workflow to achieve best-case sensitivity for IGH gene fusion detection. Using Nextflow, we developed a simplified workflow containing the algorithms FusionCatcher, Arriba, and STAR-Fusion. We analysed samples from 35 patients harbouring IGH fusions (IGH::CRLF2 n = 17, IGH::DUX4 n = 15, IGH::EPOR n = 3) and assessed the detection rates for each caller, before optimizing the parameters to enhance sensitivity for IGH fusions. Initial results showed that FusionCatcher and Arriba outperformed STAR-Fusion (85-89% vs. 29% of IGH fusions reported). We found that extensive filtering in STAR-Fusion hindered IGH reporting. By adjusting specific filtering steps (e.g., read support, fusion fragments per million total reads), we achieved a 94% reporting rate for IGH fusions with STAR-Fusion. This analysis highlights the importance of filtering optimization for IGH gene fusion events, offering alternative workflows for difficult-to-detect high-risk B-ALL subtypes.