Search

Abstract Background Single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) have led to groundbreaking advancements in life sciences. To develop bioinformatics tools for scRNA-seq and SRT data and perform unbiased benchmarks, data simulation has been widely adopted by providing explicit ground truth and generating customized datasets. However, the performance of simulation methods under multiple scenarios has not been comprehensively assessed, making it challenging to choose suitable methods without practical guidelines. Results We systematically evaluated 49 simulation methods developed for scRNA-seq and/or SRT data in terms of accuracy, functionality, scalability, and usability using 152 reference datasets derived from 24 platforms. SRTsim, scDesign3, ZINB-WaVE, and scDesign2 have the best accuracy performance across various platforms. Unexpectedly, some methods tailored to scRNA-seq data have potential compatibility for simulating SRT data. Lun, SPARSim, and scDesign3-tree outperform other methods under corresponding simulation scenarios. Phenopath, Lun, Simple, and MFA yield high scalability scores but they cannot generate realistic simulated data. Users should consider the trade-offs between method accuracy and scalability (or functionality) when making decisions. Additionally, execution errors are mainly caused by failed parameter estimations and appearance of missing or infinite values in calculations. We provide practical guidelines for method selection, a standard pipeline Simpipe ( https://github.com/duohongrui/simpipe ; https://doi.org/10.5281/zenodo.11178409 ), and an online tool Simsite ( https://www.ciblab.net/software/simshiny/ ) for data simulation. Conclusions No method performs best on all criteria, thus a good-yet-not-the-best method is recommended if it solves problems effectively and reasonably. Our comprehensive work provides crucial insights for developers on modeling gene expression data and fosters the simulation process for users.

A phylogenetic tree can reflect the evolutionary relationships between species or gene families, and they play a critical role in modern biological research. In this review, we summarize common methods for constructing phylogenetic trees, including distance methods, maximum parsimony, maximum likelihood, Bayesian inference, and tree-integration methods (supermatrix and supertree). Here we discuss the advantages, shortcomings, and applications of each method and offer relevant codes to construct phylogenetic trees from molecular data using packages and algorithms in R. This review aims to provide comprehensive guidance and reference for researchers seeking to construct phylogenetic trees while also promoting further development and innovation in this field. By offering a clear and concise overview of the different methods available, we hope to enable researchers to select the most appropriate approach for their specific research questions and datasets.

Erythritol has shown excellent insecticidal performance against a wide range of insect species, but the molecular mechanism by which it causes insect mortality and sterility is not fully understood. The mortality and sterility of Drosophila melanogaster were assessed after feeding with 1M erythritol for 72 h and 96 h, and gene expression profiles were further compared through RNA sequencing. Enrichment analysis of GO and KEGG revealed that expressions of the adipokinetic hormone gene (Akh), amylase gene (Amyrel), α-glucosidase gene (Mal-B1/2, Mal-A1-4, Mal-A7/8), and triglyceride lipase gene (Bmm) were significantly up-regulated, while insulin-like peptide genes (Dilp2, Dilp3 and Dilp5) were dramatically down-regulated. Seventeen genes associated with eggshell assembly, including Dec-1 (down 315-fold), Vm26Ab (down 2014-fold) and Vm34Ca (down 6034-fold), were significantly down-regulated or even showed no expression. However, there were no significant differences in the expression of three diuretic hormone genes (DH44, DH31, CAPA) and eight aquaporin genes (Drip, Big brain, AQP, Eglp1, Eglp2, Eglp3, Eglp4 and Prip) involved in osmolality regulation (all p value > 0.05). We concluded that erythritol, a competitive inhibitor of α-glucosidase, severely reduced substrates and enzyme binding, inhibiting effective carbohydrate hydrolysis in the midgut and eventually causing death due to energy deprivation. It was clear that Drosophila melanogaster did not die from the osmolality of the hemolymph. Our findings elucidate the molecular mechanism underlying the mortality and sterility in Drosophila melanogaster induced by erythritol feeding. It also provides an important theoretical basis for the application of erythritol as an environmentally friendly pesticide.

Insect-associated microorganisms play important roles in the health and development of insects. This study aimed to investigate the similarities and differences in bacterial community structure and composition between the larval gut of Osmia excavata, nest soil, and brood provision from the nest tube. We sequenced larvae gut and their environments’ microorganisms of O. excavata from four locations based on full-length 16S rRNA gene amplicons. The results showed 156, 280, and 366 bacterial OTUs from gut, brood provision, and nest soil, respectively, and three groups shared 131 bacterial OTUs. In the gut, the top two dominant bacteria were Sodalis praecaptivus (68.99%), Lactobacillus micheneri (17.95%). In the brood provision, the top two dominant bacteria were S. praecaptivus (26.66%), Acinetobacter nectaris (13.05%), and in the nest soil, the two most abundant bacteria were Gaiella occulta (4.33%), Vicinamibacter silvestris (3.88%). There were significant differences in diversity between the brood provision groups and the nest soil groups, respectively. Three of the four locations did not differ for gut microbial diversity. Bacteria similar to other solitary bees also existed in the gut of the larvae. Results indicated when the habitat environments were similar, the bacterial community diversity of the gut of O. excavata was similar, despite significant differences among brood provisions and soils, respectively.

Asthma is a common chronic airway disease worldwide. Due to its clinical and genetic heterogeneity, the cellular and molecular processes in asthma are highly complex and relatively unknown. To discover novel biomarkers and the molecular mechanisms underlying asthma, several studies have been conducted by focusing on gene expression patterns in epithelium through microarray analysis. However, few robust specific biomarkers were identified and some inconsistent results were observed. Therefore, it is imperative to conduct a robust analysis to solve these problems. Herein, an integrated gene expression analysis of ten independent, publicly available microarray data of bronchial epithelial cells from 348 asthmatic patients and 208 healthy controls was performed. As a result, 78 up- and 75 down-regulated genes were identified in bronchial epithelium of asthmatics. Comprehensive functional enrichment and pathway analysis revealed that response to chemical stimulus, extracellular region, pathways in cancer, and arachidonic acid metabolism were the four most significantly enriched terms. In the protein-protein interaction network, three main communities associated with cytoskeleton, response to lipid, and regulation of response to stimulus were established, and the most highly ranked 6 hub genes (up-regulated CD44, KRT6A, CEACAM5, SERPINB2, and down-regulated LTF and MUC5B) were identified and should be considered as new biomarkers. Pathway cross-talk analysis highlights that signaling pathways mediated by IL-4/13 and transcription factor HIF-1α and FOXA1 play crucial roles in the pathogenesis of asthma. Interestingly, three chemicals, polyphenol catechin, antibiotic lomefloxacin, and natural alkaloid boldine, were predicted and may be potential drugs for asthma treatment. Taken together, our findings shed new light on the common molecular pathogenesis mechanisms of asthma and provide theoretical support for further clinical therapeutic studies.