Mantle cell lymphoma (MCL) is a rare, incurable, and aggressive B-cell non-Hodgkin lymphoma (NHL). Early MCL diagnosis and treatment is critical and puzzling due to inter/intra-tumoral heterogeneity and limited understanding of the underlying molecular mechanisms. We developed and applied a multifaceted analysis of selected publicly available transcriptomic data of well-defined MCL stages, integrating network-based methods for pathway enrichment analysis, co-expression module alignment, drug repurposing, and prediction of effective drug combinations. We demonstrate the \"butterfly effect\" emerging from a small set of initially differentially expressed genes, rapidly expanding into numerous deregulated cellular processes, signaling pathways, and core machineries as MCL becomes aggressive. We explore pathogenicity-related signaling circuits by detecting common co-expression modules in MCL stages, pointing out, among others, the role of VEGFA and SPARC proteins in MCL progression and recommend further study of precise drug combinations. Our findings highlight the benefit that can be leveraged by such an approach for better understanding pathobiology and identifying high-priority novel diagnostic and prognostic biomarkers, drug targets, and efficacious combination therapies against MCL that should be further validated for their clinical impact.

Despite its prevalence, preeclampsia (PE) remains unclear as to its etiology. Here, we aimed to investigate the mechanisms regulating differences in the gene expression of zinc-finger protein 516 (ZNF516) in the placenta. The expression of the placental ZNF516 gene and its association with critical clinical markers were verified, and a rigorous correlation analysis was conducted. With a dual-luciferase reporter gene assay, microRNA targeting the ZNF516 gene was predicted and confirmed. Finally, the molecular processes associated with ZNF516 were explored via microarray and bioinformatic analyses. In hypoxic conditions, miR-371-5p expression was reduced, resulting in ZNF516 expression being induced. Moreover, ZNF516 was shown to hinder trophoblast cell migration and invasion while enhancing trophoblast cell death in various in vitro cellular assays, such as cell counting kit-8, colony formation, wound healing, and Transwell assays. Our findings reveal a new regulatory network facilitated by ZNF516. ZNF516 overexpression inhibits trophoblast growth, movement, and penetration, potentially causing problems with placenta formation with the help of miR-371-5p suppression.

Abnormal miRNA expression has been associated with breast cancer. Knowing miRNA and its target genes gives a better understanding of the biological mechanism behind the development of breast cancer. Here, we evaluated the potential prognostic and predictive values of miRNAs in breast cancer development by analyzing Malay women with breast cancer expression profiles. Seven differentially expressed miRNAs (DEMs) were subjected to miRNA‒target interaction network analysis (MTIN). A comprehensive MTIN was developed by integrating the information on miRNA and target gene interactions from five independent databases, including DIANA-TarBase, miRTarBase, miRNet, miRDB, and DIANA-microT. To understand the role of miRNAs in the progress of breast cancer, functional enrichment analysis of the miRNA target genes was conducted, followed by survival analysis to assess the prognostic values of the miRNAs and their target genes. In total, 1416 interactions were discovered among seven DEMs and 1274 target genes with a confidence score (CS) > 0.8. The overall survival analysis of the three most DEMs revealed a significant association of miR-27b-3p with poor prognosis in the TCGA breast cancer patient cohort. Further functional analysis of 606 miR-27b-3p target genes revealed their involvement in cancer-related processes and pathways, including the progesterone receptor signaling pathway, PI3K-Akt pathway, and EGFR transactivation. Notably, six high-confidence target genes (BTG2, DNAJC13, GRB2, GSK3B, KRAS, and UBR5) were discovered to be associated with worse overall survival in breast cancer patients, underscoring their essential roles in breast cancer development. Thus, we suggest that miR-27b-3p has significant potential as a biomarker for detecting breast cancer and can provide valuable understanding regarding the molecular mechanisms of the disease.

MicroRNAs (miRNAs) are a subset of small non-coding single-stranded RNA molecules involved in the regulation of post-transcriptional gene expression of a variety of transcript targets. Therefore altered miRNA expression may result in the dysregulation of key genes and biological pathways that has been reported with the onset and progression of neurodegenerative diseases, such as Amyotrophic lateral sclerosis (ALS). ALS is marked by a progressive degeneration of motor neurons (MNs) present in the spinal cord, brain stem and motor cortex. Although the pathomechanism underlying molecular interactions of ALS remains poorly understood, alterations in RNA metabolism, including dysregulation of miRNA expression in familial as well as sporadic forms are still scarcely studied. In this study, we performed combined transcriptomic data and miRNA profiling in MN samples of the same samples of iPSC-derived MNs from SOD1- and TARDBP (TDP-43 protein)-mutant-ALS patients and healthy controls. We report a global upregulation of mature miRNAs, and suggest that differentially expressed (DE) miRNAs have a significant impact on mRNA-level in SOD1-, but not in TARDBP-linked ALS. Furthermore, in SOD1-ALS we identified dysregulated miRNAs such as miR-124-3p, miR-19b-3p and miR-218 and their potential targets previously implicated in important functional process and pathogenic pathways underlying ALS. These miRNAs may play key roles in the neuronal development and cell survival related functions in SOD1-ALS. Altogether, we provide evidence of miRNA regulated genes expression mainly in SOD1 rather than TDP43-ALS.

Batch effects are notoriously common technical variations in multiomics data and may result in misleading outcomes if uncorrected or over-corrected. A plethora of batch-effect correction algorithms are proposed to facilitate data integration. However, their respective advantages and limitations are not adequately assessed in terms of omics types, the performance metrics, and the application scenarios.
As part of the Quartet Project for quality control and data integration of multiomics profiling, we comprehensively assess the performance of seven batch effect correction algorithms based on different performance metrics of clinical relevance, i.e., the accuracy of identifying differentially expressed features, the robustness of predictive models, and the ability of accurately clustering cross-batch samples into their own donors. The ratio-based method, i.e., by scaling absolute feature values of study samples relative to those of concurrently profiled reference material(s), is found to be much more effective and broadly applicable than others, especially when batch effects are completely confounded with biological factors of study interests. We further provide practical guidelines for implementing the ratio based approach in increasingly large-scale multiomics studies.
Multiomics measurements are prone to batch effects, which can be effectively corrected using ratio-based scaling of the multiomics data. Our study lays the foundation for eliminating batch effects at a ratio scale.

Circular RNAs (circRNAs) are important regulators of diverse biological processes of plants. However, the evolution and potential functions of circRNAs during winter dormancy and spring bud flushing of tea plant is largely unknown. Using RNA-seq data, a total of 1184 circRNAs were identified in the winter dormant and spring bud flushing leaf samples of tea plants in two different cultivars exhibiting different duration of winter dormancy. A total of 156 circRNAs are found to be differentially expressed and the weighted gene co-expression network (WGCNA) analysis revealed that 22 and 20 differentially expressed-circRNAs (DE-circRNAs) positively correlated with the flushing and dormant leaf traits, respectively, in both the tea cultivars used. Some transcription factors (TFs) viz. MYB, WRKY, ERF, bHLH and several genes related to secondary metabolite biosynthetic pathways are found to co-express with circRNAs. DE-circRNAs also predicted to interact with miRNAs and can regulate phytohormone biosynthesis and various signalling pathways in tea plant. This study uncovers the potential roles of circRNAs to determine winter dormancy and spring bud flushing conditions in tea plants.

The response of Spodoptera frugiperda genes toward insecticides is crucial for guiding insecticide use. The regulation of the S. frugiperda genes via long noncoding RNAs (lncRNAs) under insecticide treatment should be investigated. In this study, 452 differentially expressed lncRNAs were identified by analyzing RNA-sequencing data of S. frugiperda under 23 pesticide treatments. We found 59 and 43 differentially expressed lncRNAs that could regulate detoxification-related cytochrome P450 and UDP-glucuronosyltransferase genes, respectively. Furthermore, the target genes of differentially expressed lncRNAs were enriched in Pfam, including chitin bind 4 and gene ontology terms such as structural constituent of the cuticle, revealing their potential mechanism of action on the growth inhibition of S. frugiperda larvae. Insecticide-specific expression of lncRNAs highlights the properties and commonalities of different insecticide-induced lncRNA regulatory mechanisms. To conclude, the results of this study provide new insights and perspectives on the use of 23 insecticides via lncRNA regulation of mRNAs.

Avian pathogenic E. coli (APEC) can cause localized and systemic diseases in poultry, threatening human health via meat or egg contamination and resulting in considerable economic losses to the poultry industry globally. Increasing evidence shows circRNAs were widely involved in various biological processes. However, the role of circRNAs in the host response against APEC infection, especially correlated with the regulation of RIP2, remains unclear. Herein, the RNAseq technology was used to identify the circRNA expression profiles in the overexpression of RIP2 macrophages with or without APEC infection. A total of 256 and 287 differentially expressed (DE) circRNAs were identified in the overexpression of RIP2 group (oeRIP2) vs. the wild-type group (WT) and oeRIP2 + APEC vs. APEC, respectively, whose parental genes were involved in MAPK signalling pathway, Wnt signalling pathway, focal adhesion, tight junction, and VEGF signalling pathways. Specifically, the key circRNAs, such as 5:814443-825127, 10:18922360-18928461, 2:8746306-8750639, and 2:124177751-124184063 might play a critical role in APEC infection and the regulation of RIP2. As a whole, these findings will facilitate understanding the molecular mechanism underlying circRNAs, especially related to the regulation of the RIP2 gene. Meanwhile, the study may offer new ideas to improve host immune and inflammatory response against APEC infection.

A miRNA-mRNA combination analysis was performed on the longissimus dorsi muscle of adult Queshan Black and Large White pigs by RNA-seq technology to reveal the molecular mechanism affecting pork quality traits. The sequencing results showed that 39 miRNAs were differentially expressed between Queshan Black and Large White pigs, which targeted 5234 mRNAs, and 15 differentially expressed miRNAs targeted 86 differentially expressed mRNAs. The qRT-PCR results showed that miRNAs showed similar expression patterns to RNA-seq. The GO analysis indicated that differentially expressed miRNAs with differential target mRNAs were primarily involved in biological processes such as phospholipase activity, MAP-kinase scaffold activity, lipase activity, and regulation of the extent of cell growth. The KEGG analysis also revealed that such mRNAs were significantly enriched in the ECM-receptor interaction, sphingolipid metabolism, apoptosis, PI3K-Akt signaling pathway, and AMPK signaling pathway. In addition, software predictions showed that 17 (13 of which were upregulated and four were downregulated) of 39 differentially expressed miRNAs targeted 118 negatively correlated expression mRNAs. The upregulated miRNAs contained 103 negatively correlated target mRNAs, whereas the downregulated miRNAs contained 15 negatively correlated target mRNAs. The GO analysis showed that such mRNAs were primarily involved in MAP-kinase scaffold activity, myoblast development, and peptidyl-lysine methylation, and the KEGG analysis showed significant enrichment in ECM-receptor interaction and focal adhesion. The functional enrichment analysis of miRNA target genes revealed that miR-328 was screened out as a key miRNA, and preliminary functional validation was performed. Moreover, the overexpressed miR-328 could affect the expression of proliferation-related genes, such as CDK2, CDK4, CCNB1, CCND1, CCNE1, and PCNA. These results indicated that miR-328 may regulate fat deposition and affect meat quality by influencing related pathways. This study revealed that the miRNA-mRNA regulatory axis affects fat deposition and skeletal muscle development, which provides a theoretical basis for further study on the molecular mechanism of meat quality.

There is growing evidences indicating that long intergenic ncRNAs (lincRNAs) play key roles in plant development and stress responses. To research tomato lincRNA functions during the interaction between tomato and Ralstonia solanacearum, RNA-seq data of tomato plants inoculated with R. solanacearum was analyzed. In this study, 315 possible lincRNAs were identified from RNA-seq data. Then 23 differentially expressed lincRNAs between tomato plants inoculated with R. solanacearum and control were identified and a total of 171 possible target genes for these differentially expressed lincRNAs were predicted. Through GO and KEGG analysis, we found that lincRNA might be involved in jasmonic acid and ethylene signaling pathways to respond to tomato bacterial wilt infection. Furthermore, lincRNA may also be involved in regulating the expression of AGO protein. Subsequently, analysis of expression patterns between differentially expressed lincRNAs and adjacent mRNAs by qRT-PCR revealed that part of lincRNAs and their possible target genes exhibited positive correlation. Taken together, these results suggest that lincRNAs play potential roles in tomato against R. solanacearum infection and will provide fundamental information about the lincRNA-based plant defense mechanisms.