Investigating the basis of lineage decisions and developmental trajectories in the dorsal spinal cord through pseudotime analyses.

Dorsal interneurons (dIs) in the spinal cord encode the perception of touch, pain, heat, itchiness and proprioception. Previous studies using genetic strategies in animal models have revealed important insights into dI development, but the molecular details of how dIs arise as distinct populations of neurons remain incomplete. We have developed a resource to investigate dI fate specification by combining a single-cell RNA-Seq atlas of mouse embryonic stem cell-derived dIs with pseudotime analyses. To validate this in silico resource as a useful tool, we used it to first identify genes that are candidates for directing the transition states that lead to distinct dI lineage trajectories, and then validated them using in situ hybridization analyses in the developing mouse spinal cord in vivo. We have also identified an endpoint of the dI5 lineage trajectory and found that dIs become more transcriptionally homogeneous during terminal differentiation. This study introduces a valuable tool for further discovery about the timing of gene expression during dI differentiation and demonstrates its utility in clarifying dI lineage relationships.

Comparative Analysis of Droplet-Based Ultra-High-Throughput Single-Cell RNA-Seq Systems.

Since its establishment in 2009, single-cell RNA sequencing (RNA-seq) has been a major driver behind progress in biomedical research. In developmental biology and stem cell studies, the ability to profile single cells confers particular benefits. Although most studies still focus on individual tissues or organs, the recent development of ultra-high-throughput single-cell RNA-seq has demonstrated potential power in characterizing more complex systems or even the entire body. However, although multiple ultra-high-throughput single-cell RNA-seq systems have attracted attention, no systematic comparison of these systems has been performed. Here, with the same cell line and bioinformatics pipeline, we developed directly comparable datasets for each of three widely used droplet-based ultra-high-throughput single-cell RNA-seq systems, inDrop, Drop-seq, and 10X Genomics Chromium. Although each system is capable of profiling single-cell transcriptomes, their detailed comparison revealed the distinguishing features and suitable applications for each system.

Efficient and Accurate Quantitative Profiling of Alternative Splicing Patterns of Any Complexity on a Laptop.

Alternative splicing (AS) is a widespread process underlying the generation of transcriptomic and proteomic diversity and is frequently misregulated in human disease. Accordingly, an important goal of biomedical research is the development of tools capable of comprehensively, accurately, and efficiently profiling AS. Here, we describe Whippet, an easy-to-use RNA-seq analysis method that rapidly-with hardware requirements compatible with a laptop-models and quantifies AS events of any complexity without loss of accuracy. Using an entropic measure of splicing complexity, Whippet reveals that one-third of human protein coding genes produce transcripts with complex AS events involving co-expression of two or more principal splice isoforms. We observe that high-entropy AS events are more prevalent in tumor relative to matched normal tissues and correlate with increased expression of proto-oncogenic splicing factors. Whippet thus affords the rapid and accurate analysis of AS events of any complexity, and as such will facilitate future biomedical research.

Histological and gene-expression analyses of pyloric sphincter formation during stomach metamorphosis in Xenopus laevis.

During anuran metamorphosis from herbivorous tadpoles to carnivorous frogs, the gastrointestinal (GI) tract undergoes drastic remodeling, such as the formation of the stomach-intestine boundary and the development of the pyloric sphincter at the posterior end of the stomach. However, the morphogenetic process and molecular mechanisms of how the pyloric sphincter is formed during metamorphosis, instead of during embryogenesis as in amniotes, are largely uninvestigated. Using the African clawed frog Xenopus laevis, we histologically examined the development of the pylorus region from embryonic to froglet stages and performed spatiotemporal gene expression analyses. We found that the pyloric sphincter is formed at a flexure within the pyloric region during metamorphic climax, and that the pyloric and duodenal epithelia, which are morphologically indistinguishable before sphincter formation, become clearly demarcated by the sphincter at the end of metamorphosis. Consistent with these morphological changes, expression domains of a stomach marker barx1 and an intestine marker cdx2 overlapped until late metamorphic climax, but became separated after metamorphosis. Despite the absence of the sphincter before metamorphosis, various genes crucial for sphincter formation in amniotes were already expressed in the pylorus region of Xenopus embryos. RNA-sequencing analysis at pre-metamorphic and metamorphic-climax stages suggest unappreciated roles of genes, such as those for retinoic acid signaling and various transcription factors, in suppressing or promoting sphincter formation. These data provide histological and molecular insights into the heterochrony of the pyloric sphincter formation in amniotes and anurans.

LsRTDv1, a reference transcript dataset for accurate transcript-specific expression analysis in lettuce.

Accurate quantification of gene and transcript-specific expression, with the underlying knowledge of precise transcript isoforms, is crucial to understanding many biological processes. Analysis of RNA sequencing data has benefited from the development of alignment-free algorithms which enhance the precision and speed of expression analysis. However, such algorithms require a reference transcriptome. Here we generate a reference transcript dataset (LsRTDv1) for lettuce (cv. Saladin), combining long- and short-read sequencing with publicly available transcriptome annotations, and filtering to keep only transcripts with high-confidence splice junctions and transcriptional start and end sites. LsRTDv1 identifies novel genes (mostly long non-coding RNAs) and increases the number of transcript isoforms per gene in the lettuce genome from 1.4 to 2.7. We show that LsRTDv1 significantly increases the mapping rate of RNA-seq data from a lettuce time-series experiment (mock- and Botrytis cinerea-inoculated) and enables detection of genes that are differentially alternatively spliced in response to infection as well as transcript-specific expression changes. LsRTDv1 is a valuable resource for investigation of transcriptional and alternative splicing regulation in lettuce.

Identification of differential m6A RNA methylomes and ALKBH5 as a potential prevention target in the developmental neurotoxicity induced by multiple sevoflurane exposures.

Sevoflurane, as a commonly used inhaled anesthetic for pediatric patients, has been reported that multiple sevoflurane exposures are associated with a greater risk of developing neurocognitive disorder. N6-Methyladenosine (m6A), as the most common mRNA modification in eukaryotes, has emerged as a crucial regulator of brain function in processes involving synaptic plasticity, learning and memory, and neurodevelopment. Nevertheless, the relevance of m6A RNA methylation in the multiple sevoflurane exposure-induced developmental neurotoxicity remains mostly elusive. Herein, we evaluated the genome-wide m6A RNA modification and gene expression in hippocampus of mice that received with multiple sevoflurane exposures using m6A-sequencing (m6A-seq) and RNA-sequencing (RNA-seq). We discovered 19 genes with differences in the m6A methylated modification and differential expression in the hippocampus. Among these genes, we determined that a total of nine differential expressed genes may be closely associated with the occurrence of developmental neurotoxicity induced by multiple sevoflurane exposures. We further found that the alkB homolog 5 (ALKBH5), but not methyltransferase-like 3 (METTL3) and Wilms tumor 1-associated protein (WTAP), were increased in the hippocampus of mice that received with multiple sevoflurane exposures. And the IOX1, as an inhibitor of ALKBH5, significantly improved the learning and memory defects and reduced neuronal damage in the hippocampus of mice induced by multiple sevoflurane exposures. The current study revealed the role of m6A methylated modification and m6A-related regulators in sevoflurane-induced cognitive impairment, which might provide a novel insight into identifying biomarkers and therapeutic strategies for inhaled anesthetic-induced developmental neurotoxicity.

Integrative bioinformatics approach yields a novel gene expression risk model for prognosis and progression prediction in prostate cancer.

Prostate cancer (PCa), a prevalent malignancy among elderly males, exhibits a notable rate of advancement, even when subjected to conventional androgen deprivation therapy or chemotherapy. An effective progression prediction model would prove invaluable in identifying patients with a higher progression risk. Using bioinformatics strategies, we integrated diverse data sets of PCa to construct a novel risk model predicated on gene expression and progression-free survival (PFS). The accuracy of the model was assessed through validation using an independent data set. Eight genes were discerned as independent prognostic factors and included in the prediction model. Patients assigned to the high-risk cohort demonstrated a diminished PFS, and the areas under the curve of our model in the validation set for 1-year, 3-year, and 5-year PFS were 0.9325, 0.9041 and 0.9070, respectively. Additionally, through the application of single-cell RNA sequencing to two castration-related prostate cancer (CRPC) samples and two hormone-related prostate cancer (HSPC) samples, we discovered that luminal cells within CRPC exhibited an elevated risk score. Subsequent molecular biology experiments corroborated our findings, illustrating heightened SYK expression levels within tumour tissues and its contribution to cancer cell migration. We found that the knockdown of SYK could inhibit migration in PCa cells. Our progression-related risk model demonstrated the potential prognostic value of SYK and indicated its potential as a target for future diagnosis and treatment strategies in PCa management.

Transcriptome disclosure of hormones inducing stigma exsertion in Nicotiana tabacum by corolla shortening.

BACKGROUND: Stigma exsertion is an essential agricultural trait that can promote cross-pollination to improve hybrid seed production efficiency. However, the molecular mechanism controlling stigma exsertion remains unknown. RESULTS: In this study, the Nicotiana tabacum cv. K326 and its two homonuclear-heteroplasmic lines, MSK326 (male-sterile) and MSK326SE (male-sterile and stigma exserted), were used to investigate the mechanism of tobacco stigma exsertion. A comparison of the flowers between the three lines showed that the stigma exsertion of MSK326SE was mainly due to corolla shortening. Therefore, the corollas of the three lines were sampled and presented for RNA-seq analysis, which found 338 candidate genes that may cause corolla shortening. These genes were equally expressed in K326 and MSK326, but differentially expressed in MSK326SE. Among these 338 genes, 15 were involved in hormone synthesis or signal transduction pathways. Consistently, the content of auxin, dihydrozeatin, gibberellin, and jasmonic acid was significantly decreased in the MSK326SE corolla, whereas abscisic acid levels were significantly increased. Additionally, seven genes involved in cell division, cell cycle, or cell expansion were identified. Protein-protein interaction network analysis identified 45 nodes and 79 protein interactions, and the largest module contained 20 nodes and 52 protein interactions, mainly involved in the hormone signal transduction and pathogen defensive pathways. Furthermore, a putative hub gene coding a serine/threonine-protein kinase was identified for the network. CONCLUSIONS: Our results suggest that hormones may play a key role in regulating tobacco stigma exsertion induced by corolla shortening.

Genome-wide identification and in-silico expression analysis of CCO gene family in sunflower (Helianthus annnus) against abiotic stress.

Carotenoid cleavage oxygenases (CCOs) enzymes play an important role in plant growth and development by producing a wide array of apocarotenoids and their derivatives. These compounds are vital for colouring flowers and fruits and synthesizing plant hormones such as abscisic acid and strigolactones. Despite their importance, the gene family responsible for CCO enzymes in sunflowers has not been identified. In this study, we identify the CCO genes of the sunflower plant to fill this knowledge gap. Phylogenetic and synteny analysis indicated that the Helianthus annnus CCO (HaCCO) genes were conserved in different plant species and they could be divided into three subgroups based on their conserved domains. Analysis using MEME tool and multiple sequence alignment identified conserved motifs in the HaCCO gene sequence. Cis-regulatory elements (CREs) analysis of the HaCCO genes indicated the presence of various responsive elements related to plant hormones, development, and responses to both biotic and abiotic stresses. This implies that these genes may respond to plant hormones, developmental cues, and drought stress, offering potential applications in the development of more resistant crops. Genes belonging to the 9-cis-epoxy carotenoid dioxygenases (NCED) subgroups predominantly exhibited chloroplast localization, whereas the genes found in other groups are primarily localized in the cytoplasm. These 21 identified HaCCOs were regulated by 60 miRNAs, indicating the crucial role of microRNAs in gene regulation in sunflowers. Gene expression analysis under drought stress revealed significant up-regulation of HaNCED16 and HaNCED19, genes that are pivotal in ABA hormone biosynthesis. During organ-specific gene expression analysis, HaCCD12 and HaCCD20 genes exhibit higher activity in leaves, indicating a potential role in leaf pigmentation. This study provides a foundation for future research on the regulation and functions of the CCO gene family in sunflower and beyond. There is potential for developing molecular markers that could be employed in breeding programs to create new sunflower lines resistant to biotic and abiotic stresses.

Pain Sensitivity Genes as Therapeutic Targets in Knee Osteoarthritis: A Comprehensive Analysis.

Pain sensitivity is a significant factor in knee osteoarthritis (KOA), influencing patient outcomes and complicating treatment. Genetic differences, particularly in pain-sensing genes (PSRGs), are known to contribute to the variability in pain experiences among KOA patients. This study aims to systematically analyze PSRGs in KOA to better understand their role and potential as therapeutic targets. We utilized bulk RNA-seq data from the GSE114007 and GSE169077 datasets to identify differentially expressed genes, with 20 genes found to be significantly altered. Key PSRGs, including PENK, NGF, HOXD1, and TRPA1, were identified using LASSO, SVM, and random forest algorithms. Further, KEGG and GO enrichment analyses revealed pathways such as "Neuroactive ligand-receptor interaction" and "ECM-receptor interaction," which were validated through external datasets. Single-cell RNA-seq analysis from GSE152805, GSE133449, and GSE104782 datasets demonstrated the heterogeneity and dynamic expression of PSRGs across different cell subpopulations in synovium, meniscus, and cartilage samples. UMAP and pseudotime analyses were used to visualize spatial distribution and developmental trajectories of these genes. The findings emphasize the critical roles of PSRGs in KOA, highlighting their potential as therapeutic targets and suggesting that integrating genetic information into clinical practice could significantly improve pain management and treatment strategies for KOA.

SISTER OF FCA physically associates with SKB1 to regulate flowering time in Arabidopsis thaliana.

BACKGROUND: Proper flowering time is important for the growth and development of plants, and both too early and too late flowering impose strong negative influences on plant adaptation and seed yield. Thus, it is vitally important to study the mechanism underlying flowering time control in plants. In a previous study by the authors, genome-wide association analysis was used to screen the candidate gene SISTER OF FCA (SSF) that regulates FLOWERING LOCUS C (FLC), a central gene encoding a flowering suppressor in Arabidopsis thaliana. RESULTS: SSF physically interacts with Protein arginine methyltransferase 5 (PRMT5, SKB1). Subcellular co-localization analysis showed that SSF and SKB1 interact in the nucleus. Genetically, SSF and SKB1 exist in the same regulatory pathway that controls FLC expression. Furthermore, RNA-sequencing analysis showed that both SSF and SKB1 regulate certain common pathways. CONCLUSIONS: This study shows that PRMT5 interacts with SSF, thus controlling FLC expression and facilitating flowering time control.

Exercise training upregulates CD55 to suppress complement-mediated synaptic phagocytosis in Parkinson's disease.

The primary pathological change in Parkinson's disease (PD) is the progressive degeneration of dopaminergic neurons in the substantia nigra. Additionally, excessive microglial activation and synaptic loss are also typical features observed in PD samples. Exercise trainings have been proven to improve PD symptoms, delay the disease progression as well as affect excessive microglial synaptic phagocytosis. In this study, we established a mouse model of PD by injecting mouse-derived α-synuclein preformed fibrils (M-α-syn PFFs) into the substantia nigra, and demonstrated that treadmill exercise inhibits microglial activation and synaptic phagocytosis in striatum. Using RNA-Seq and proteomics, we also found that PD involves excessive activation of the complement pathway which is closely related to over-activation of microglia and abnormal synaptic function. More importantly, exercise training can inhibit complement levels and complement-mediated microglial phagocytosis of synapses. It is probably triggered by CD55, as we observed that CD55 in the striatum significantly increased after exercise training and up-regulation of that molecule rescued motor deficits of PD mice, accompanied with reduced microglial synaptic phagocytosis in the striatum. This research elucidated the interplay among microglia, complement, and synapses, and analyzed the effects of exercise training on these factors. Our work also suggested CD55 as a complement-relevant candidate molecule for developing therapeutic strategies of PD.

Comparative transcriptome analysis reveals the redirection of metabolic flux from cell growth to astaxanthin biosynthesis in Yarrowia lipolytica.

Engineering Yarrowia lipolytica to produce astaxanthin provides a promising route. Here, Y. lipolytica M2 producing a titer of 181 mg/L astaxanthin was isolated by iterative atmospheric and room-temperature plasma mutagenesis and diphenylamine-mediated screening. Interestingly, a negative correlation was observed between cell biomass and astaxanthin production. To reveal the underlying mechanism, RNA-seq analysis of transcriptional changes was performed in high producer M2 and reference strain M1, and a total of 1379 differentially expressed genes were obtained. Data analysis revealed that carbon flux was elevated through lipid metabolism, acetyl-CoA and mevalonate supply, but restrained through central carbon metabolism in strain M2. Moreover, upregulation of other pathways such as ATP-binding cassette transporter and thiamine pyrophosphate possibly provided more cofactors for carotenoid hydroxylase and relieved cell membrane stress caused by astaxanthin insertion. These results suggest that balancing cell growth and astaxanthin production may be important to promote efficient biosynthesis of astaxanthin in Y. lipolytica.

Integration of QTL and transcriptome approaches for the identification of genes involved in tomato response to nitrogen deficiency.

Optimizing plant nitrogen (N) usage and inhibiting N leaching loss in the soil-crop system is crucial to maintaining crop yield and reducing environmental pollution. This study aimed at identifying quantitative trait loci (QTLs) and differentially expressed genes (DEGs) between two N treatments in order to list candidate genes related to nitrogen-related contrasting traits in tomato varieties. We characterized a genetic diversity core-collection (CC) and a multi-parental advanced generation intercross (MAGIC) tomato population grown in a greenhouse under two nitrogen levels and assessed several N-related traits and mapped QTLs. Transcriptome response under the two N conditions was also investigated through RNA sequencing of fruit and leaves in four parents of the MAGIC population. Significant differences in response to N input reduction were observed at the phenotypic level for biomass and N-related traits. Twenty-seven QTLs were detected for three target traits (leaf N content, leaf nitrogen balance index, and petiole NO3- content), 10 and six in the low and high N condition, respectively, while 19 QTLs were identified for plasticity traits. At the transcriptome level, 4752 and 2405 DEGs were detected between the two N conditions in leaves and fruits, respectively, among which 3628 (50.6%) in leaves and 1717 (71.4%) in fruit were genotype specific. When considering all the genotypes, 1677 DEGs were shared between organs or tissues. Finally, we integrated DEG and QTL analyses to identify the most promising candidate genes. The results highlighted a complex genetic architecture of N homeostasis in tomato and novel putative genes useful for breeding tomato varieties requiring less N input.

Transcriptome-wide profiling for melanocytes derived from newborn and adult human epidermis with enhanced proliferation.

The senescence-associated protein p16INK4A acts as a limiter element in cell-cycle progression. The loss of p16INK4A function is causally related to cellular immortalization. The increase in p16INK4A levels with advancing age was demonstrated in melanocytes. However, the characteristic difference between young and senescent melanocytes affecting immortalization of melanocytes remains unclear. In this study, we generated 10 different cell lines in total from newborn (NB) and adult (AD) primary normal human epidermal melanocytes (NHEM) using four different methods, transduction of CDK4R24C and cyclin D1 (K4D), K4D with TERT (K4DT), SV40 T-antigen (SV40T), and HPV16 E6 and E7 (E6/E7) and performed whole transcriptome sequencing analysis (RNA-Seq) to elucidate the differences of genome-wide expression profiles among cell lines. The analysis data revealed distinct differences in expression pattern between cell lines from NB and AD although no distinct biological differences were detected in analyses such as comparison of cell morphology, evaluation of cell proliferation, and cell cycle profiles. This study may provide useful in vitro models to benefit the understanding of skin-related diseases.

Genome-Wide Identification and Characterization of MYB Gene Family and Analysis of Its Sex-Biased Expression Pattern in Spinacia oleracea L.

The members of the myeloblastosis (MYB) family of transcription factors (TFs) participate in a variety of biological regulatory processes in plants, such as circadian rhythm, metabolism, and flower development. However, the characterization of MYB genes across the genomes of spinach Spinacia oleracea L. has not been reported. Here, we identified 140 MYB genes in spinach and described their characteristics using bioinformatics approaches. Among the MYB genes, 54 were 1R-MYB, 80 were 2R-MYB, 5 were 3R-MYB, and 1 was 4R-MYB. Almost all MYB genes were located in the 0-30 Mb region of autosomes; however, the 20 MYB genes were enriched at both ends of the sex chromosome (chromosome 4). Based on phylogeny, conserved motifs, and the structure of genes, 2R-MYB exhibited higher conservation relative to 1R-MYB genes. Tandem duplication and collinearity of spinach MYB genes drive their evolution, enabling the functional diversification of spinach genes. Subcellular localization prediction indicated that spinach MYB genes were mainly located in the nucleus. Cis-acting element analysis confirmed that MYB genes were involved in various processes of spinach growth and development, such as circadian rhythm, cell differentiation, and reproduction through hormone synthesis. Furthermore, through the transcriptome data analysis of male and female flower organs at five different periods, ten candidate genes showed biased expression in spinach males, suggesting that these genes might be related to the development of spinach anthers. Collectively, this study provides useful information for further investigating the function of MYB TFs and novel insights into the regulation of sex determination in spinach.

Transcriptomic Analyses Reveal Insights into the Shared Regulatory Network of Phenolic Compounds and Steviol Glycosides in Stevia rebaudiana.

Stevia rebaudiana (Bertoni) is a highly valuable crop for the steviol glycoside content in its leaves, which are no-calorie sweeteners hundreds of times more potent than sucrose. The presence of health-promoting phenolic compounds, particularly flavonoids, in the leaf of S. rebaudiana adds further nutritional value to this crop. Although all these secondary metabolites are highly desirable in S. rebaudiana leaves, the genes regulating the biosynthesis of phenolic compounds and the shared gene network between the regulation of biosynthesis of steviol glycosides and phenolic compounds still need to be investigated in this species. To identify putative candidate genes involved in the synergistic regulation of steviol glycosides and phenolic compounds, four genotypes with different contents of these compounds were selected for a pairwise comparison RNA-seq analysis, yielding 1136 differentially expressed genes. Genes that highly correlate with both steviol glycosides and phenolic compound accumulation in the four genotypes of S. rebaudiana were identified using the weighted gene co-expression network analysis. The presence of UDP-glycosyltransferases 76G1, 76H1, 85C1, and 91A1, and several genes associated with the phenylpropanoid pathway, including peroxidase, caffeoyl-CoA O-methyltransferase, and malonyl-coenzyme A:anthocyanin 3-O-glucoside-6″-O-malonyltransferase, along with 21 transcription factors like SCL3, WRK11, and MYB111, implied an extensive and synergistic regulatory network involved in enhancing the production of such compounds in S. rebaudiana leaves. In conclusion, this work identified a variety of putative candidate genes involved in the biosynthesis and regulation of particular steviol glycosides and phenolic compounds that will be useful in gene editing strategies for increasing and steering the production of such compounds in S. rebaudiana as well as in other species.

Genome-Wide Identification and Characterization of U-Box Gene Family Members and Analysis of Their Expression Patterns in Phaseolus vulgaris L. under Cold Stress.

The common bean (Phaseolus vulgaris L.) is an economically important food crop grown worldwide; however, its production is affected by various environmental stresses, including cold, heat, and drought stress. The plant U-box (PUB) protein family participates in various biological processes and stress responses, but the gene function and expression patterns of its members in the common bean remain unclear. Here, we systematically identified 63 U-box genes, including 8 tandem genes and 55 non-tandem genes, in the common bean. These PvPUB genes were unevenly distributed across 11 chromosomes, with chromosome 2 holding the most members of the PUB family, containing 10 PUB genes. The analysis of the phylogenetic tree classified the 63 PUB genes into three groups. Moreover, transcriptome analysis based on cold-tolerant and cold-sensitive varieties identified 4 differentially expressed PvPUB genes, suggesting their roles in cold tolerance. Taken together, this study serves as a valuable resource for exploring the functional aspects of the common bean U-box gene family and offers crucial theoretical support for the development of new cold-tolerant common bean varieties.

Exploring the impact of TGF-ß family gene mutations and expression on skin wound healing and tissue repair.

Transforming Growth Factor-Beta (TGF-ß) signalling pathway is of paramount importance in the processes of wound healing, epidermal integrity maintenance and development of skin cancer. The objective of this research endeavour was to clarify the impact of gene mutations and variations in expression within TGF-ß family on mechanisms of tissue repair, as well as to identify potential targets for therapeutic purposes in non-melanoma skin cancer (NMSC). The methods utilized in this study involved obtaining RNA-seq data from 224 NMSC patients and paired normal skin tissues from the PRJNA320473 and PRJEB27606 databases. The purpose of the differential gene expression analysis was to identify genes whose expression had changed significantly. In order to evaluate the effects and interrelationships of identified gene variants, structural analysis with AlphaFold and PDB data and network analysis with the STRING database were both utilized. Critical gene expression was externally validated through the utilization of the GEPIA database. Tumour tissues exhibited a notable upregulation of genes associated with the TGF-ß pathway, specifically MMP1, MMP3, MMP9, EGF, COL3A1 and COL1A2, in comparison with normal tissues. As indicated by the central node status of these genes in the network analysis, they play a crucial role in the progression of NMSCs. The results of the structural analysis suggested that mutations might cause functional disruptions. External validation of the upregulation confirmed the expression trends and emphasized the biomarker potential of the upregulated genes. In conclusion, this research offered thorough examination of molecular modifications that occur in TGF-ß family genes, which are linked to cutaneous wound healing and NMSC. The modified expression of the identified hub genes may represent innovative targets for therapeutic intervention.

Scellpam: an R package/C++ library to perform parallel partitioning around medoids on scRNAseq data sets.

BACKGROUND: Partitioning around medoids (PAM) is one of the most widely used and successful clustering method in many fields. One of its key advantages is that it only requires a distance or a dissimilarity between the individuals, and the fact that cluster centers are actual points in the data set means they can be taken as reliable representatives of their classes. However, its wider application is hampered by the large amount of memory needed to store the distance matrix (quadratic on the number of individuals) and also by the high computational cost of computing such distance matrix and, less importantly, by the cost of the clustering algorithm itself. RESULTS: Therefore, new software has been provided that addresses these issues. This software, provided under GPL license and usable as either an R package or a C++ library, calculates in parallel the distance matrix for different distances/dissimilarities ([Formula: see text], [Formula: see text], Pearson, cosine and weighted Euclidean) and also implements a parallel fast version of PAM (FASTPAM1) using any data type to reduce memory usage. Moreover, the parallel implementation uses all the cores available in modern computers which greatly reduces the execution time. Besides its general application, the software is especially useful for processing data of single cell experiments. It has been tested in problems including clustering of single cell experiments with up to 289,000 cells with the expression of about 29,000 genes per cell. CONCLUSIONS: Comparisons with other current packages in terms of execution time have been made. The method greatly outperforms the available R packages for distance matrix calculation and also improves the packages that implement the PAM itself. The software is available as an R package at https://CRAN.R-project.org/package=scellpam and as C++ libraries at https://github.com/JdMDE/jmatlib and https://github.com/JdMDE/ppamlib The package is useful for single cell RNA-seq studies but it is also applicable in other contexts where clustering of large data sets is required.