The RNA-Binding Protein YBX3 Controls Amino Acid Levels by Regulating SLC mRNA Abundance

DNA-binding protein-A promotes kidney ischemia/reperfusion injury and participates in mitochondrial function

DNA-binding protein-A (DbpA; gene: Ybx3) belongs to the cold shock protein family with known functions in cell cycling, transcription, translation, and tight junction communication. In chronic nephritis, DbpA is upregulated. However, its activities in acute injury models, such as kidney ischemia/reperfusion injury (IRI), are unclear. To study this, mice harboring Ybx3+/+, Ybx3+/- or the Ybx3-/- genotype were characterized over 24 months and following experimental kidney IRI. Mitochondrial function, number and integrity were analyzed by mitochondrial stress tests, MitoTracker staining and electron microscopy. Western Blot, immunohistochemistry and flow cytometry were performed to quantify tubular cell damage and immune cell infiltration. DbpA was found to be dispensable for kidney development and tissue homeostasis under healthy conditions. Furthermore, endogenous DbpA protein localizes within mitochondria in primary tubular epithelial cells. Genetic deletion of Ybx3 elevates the mitochondrial membrane potential, lipid uptake and metabolism, oxygen consumption rates and glycolytic activities of tubular epithelial cells. Ybx3-/- mice demonstrated protection from IRI with less immune cell infiltration, endoplasmic reticulum stress and tubular cell damage. A presumed protective mechanism was identified via upregulated antioxidant activities and reduced ferroptosis, when Ybx3 was deleted. Thus, our studies reveal DbpA acts as a mitochondrial protein with profound adverse effects on cell metabolism and highlights a protective effect against IRI when Ybx3 is genetically deleted. Hence, preemptive DbpA targeting in situations with expected IRI, such as kidney transplantation or cardiac surgery, may preserve post-procedure kidney function.

DNA-binding protein-A: a multitool in tubular epithelial cells

Acute kidney injury is still associated with high morbidity and mortality. Reichardt et al. investigated DNA-binding protein-A (Ybx3) in acute kidney injury induced by ischemia-reperfusion injury and found that mice lacking Ybx3 have altered mitochondrial function and increased antioxidant activity, making them more resistant to ischemia-reperfusion injury-acute kidney injury. The study highlights a new role of the multifaceted protein DNA-binding protein-A, which could be potentially therapeutically exploited.

piRNA PROPER Suppresses DUSP1 Translation by Targeting N6-Methyladenosine-Mediated RNA Circularization to Promote Oncogenesis of Prostate Cancer

Genetic and epigenetic alterations occur in many physiological and pathological processes. The existing knowledge regarding the association of PIWI-interacting RNAs (piRNAs) and their genetic variants on risk and progression of prostate cancer (PCa) is limited. In this study, three genome-wide association study datasets are combined, including 85,707 PCa cases and 166,247 controls, to uncover genetic variants in piRNAs. Functional investigations involved manipulating piRNA expression in cellular and mouse models to study its oncogenetic role in PCa. A specific genetic variant, rs17201241 is identified, associated with increased expression of PROPER (piRNA overexpressed in prostate cancer) in tumors and are located within the gene, conferring an increased risk and malignant progression of PCa. Mechanistically, PROPER coupled with YTHDF2 to recognize N6-methyladenosine (m6A) and facilitated RNA-binding protein interactions between EIF2S3 at 5'-untranslated region (UTR) and YTHDF2/YBX3 at 3'-UTR to promote DUSP1 circularization. This m6A-dependent mRNA-looping pattern enhanced DUSP1 degradation and inhibited DUSP1 translation, ultimately reducing DUSP1 expression and promoting PCa metastasis via the p38 mitogen-activated protein kinase (MAPK) signaling pathway. Inhibition of PROPER expression using antagoPROPER effectively suppressed xenograft growth, suggesting its potential as a therapeutic target. Thus, targeting piRNA PROPER-mediated genetic and epigenetic fine control is a promising strategy for the concurrent prevention and treatment of PCa.

Comprehensive analysis of the functions, prognostic and diagnostic values of RNA binding proteins in head and neck squamous cell carcinoma

BACKGROUND

Accumulating evidence has suggested that RNA binding protein (RBP) dysregulation plays an essential role during tumorigenesis. Here, we sought to explore the potential biological functions and clinical significance of RBP and develop diagnostic and prognostic signatures based on RBP in patients with head and neck squamous cell carcinoma (HNSCC).

METHODS

The differently expressed RBPs between HNSCC samples and their normal counterparts were identified using the Limma package. The immunohistochemistry (IHC) images of several RBPs were collected from the Human Protein Atlas database. The diagnostic signature based on RBP was built by LASSO-logistic regression and random forest. The prognostic signature based on RBP was constructed by LASSO and stepwise Cox regression analysis in the training cohort and validated in the validation cohort.

RESULTS

Eighty-four aberrantly expressed RBPs were obtained, comprising 41 up-regulated and 43 down-regulated RBPs. Seven RBP genes (CPEB3, PDCD4, ENDOU, PARP12, DNMT3B, IGF2BP1, EXO1) were identified as diagnostic-related hub genes. They were used to establish a diagnostic RBP signature risk score (DRBPS) model by the coefficients in least absolute shrinkage and selection operator (LASSO)-logistic regression analysis and showed high specificity and sensitivity in the training (area under the receiver operating characteristic curve (AUC) = 0.998), and in all validation cohorts (AUC > 0.95 for all). Similarly, seven RBP genes (MKRN3, ZC3H12D, EIF5A2, AFF3, SIDT1, RBM24, and NR0B1) were identified as prognosis-associated hub genes by LASSO and stepwise multiple Cox regression analyses and were used to construct the prognostic model named as PRBPS. The AUC of the time-dependent receiver operator characteristic curve of the prognostic model was 0.664 at 3 years and 0.635 at 5 years in the training cohort and 0.720, 0.777 in the validation cohort, showing a favorable predictive efficacy for prognosis in HNSCC.

CONCLUSIONS

Our results demonstrate the value of consideration of RBP in the diagnosis and prognosis for HNSCC and provide a novel insight into understanding the potential role of dysregulated RBP in HNSCC.

Identification of tumor-suppressive miRNAs that target amino acid transporter LAT1 and exhibit anti-proliferative effects on cholangiocarcinoma cells

Amino acid transporter LAT1 is highly upregulated in various cancer types, including cholangiocarcinoma (CHOL), and contributes to the rapid proliferation of cancer cells and disease progression. However, the molecular mechanisms underlying the pathological upregulation of LAT1 remain largely unknown. This study pursued the possibility of miRNA-mediated regulation of the LAT1 expression in CHOL cells. Using online target prediction methods, we extracted five candidate miRNAs commonly predicted to regulate the LAT1 expression. Three of them, miR-194-5p, miR-122-5p, and miR-126-3p, were significantly downregulated in CHOL cancer compared to normal tissues. Correlation analysis revealed weak-to-moderate negative correlations between the expression of these miRNAs and LAT1 mRNA in CHOL cancer tissues. We selected miR-194-5p and miR-122-5p for further analyses and found that both miRNAs functionally target 3'UTR of LAT1 mRNA by a luciferase-based reporter assay. Transfection of the miRNA mimics significantly suppressed the LAT1 expression at mRNA and protein levels and inhibited the proliferation of CHOL cells, with a trend of affecting intracellular amino acids and amino acid-related signaling pathways. This study indicates that the decreased expression of these LAT1-targeting tumor-suppressive miRNAs contributes to the upregulation of LAT1 and the proliferation of CHOL cells, highlighting their potential for developing novel cancer therapeutics and diagnostics.

SARS-CoV-2 RNA stabilizes host mRNAs to elicit immunopathogenesis

SARS-CoV-2 RNA interacts with host factors to suppress interferon responses and simultaneously induces cytokine release to drive the development of severe coronavirus disease 2019 (COVID-19). However, how SARS-CoV-2 hijacks host RNAs to elicit such imbalanced immune responses remains elusive. Here, we analyzed SARS-CoV-2 RNA in situ structures and interactions in infected cells and patient lung samples using RIC-seq. We discovered that SARS-CoV-2 RNA forms 2,095 potential duplexes with the 3' UTRs of 205 host mRNAs to increase their stability by recruiting RNA-binding protein YBX3 in A549 cells. Disrupting the SARS-CoV-2-to-host RNA duplex or knocking down YBX3 decreased host mRNA stability and reduced viral replication. Among SARS-CoV-2-stabilized host targets, NFKBIZ was crucial for promoting cytokine production and reducing interferon responses, probably contributing to cytokine storm induction. Our study uncovers the crucial roles of RNA-RNA interactions in the immunopathogenesis of RNA viruses such as SARS-CoV-2 and provides valuable host targets for drug development.

The YBX3 RNA-binding protein posttranscriptionally controls SLC1A5 mRNA in proliferating and differentiating skeletal muscle cells

In humans, skeletal muscles comprise nearly 40% of total body mass, which is maintained throughout adulthood by a balance of muscle protein synthesis and breakdown. Cellular amino acid (AA) levels are critical for these processes, and mammalian cells contain transporter proteins that import AAs to maintain homeostasis. Until recently, the control of transporter regulation has largely been studied at the transcriptional and posttranslational levels. However, here, we report that the RNA-binding protein YBX3 is critical to sustain intracellular AAs in mouse skeletal muscle cells, which aligns with our recent findings in human cells. We find that YBX3 directly binds the solute carrier (SLC)1A5 AA transporter messenger (m)RNA to posttranscriptionally control SLC1A5 expression during skeletal muscle cell differentiation. YBX3 regulation of SLC1A5 requires the 3' UTR. Additionally, intracellular AAs transported by SLC1A5, either directly or indirectly through coupling to other transporters, are specifically reduced when YBX3 is depleted. Further, we find that reduction of the YBX3 protein reduces proliferation and impairs differentiation in skeletal muscle cells, and that YBX3 and SLC1A5 protein expression increase substantially during skeletal muscle differentiation, independently of their respective mRNA levels. Taken together, our findings suggest that YBX3 regulates AA transport in skeletal muscle cells, and that its expression is critical to maintain skeletal muscle cell proliferation and differentiation.

Establishment and analysis of a novel diagnostic model for systemic juvenile idiopathic arthritis based on machine learning

BACKGROUND

Systemic juvenile idiopathic arthritis (SJIA) is a form of childhood arthritis with clinical features such as fever, lymphadenopathy, arthritis, rash, and serositis. It seriously affects the growth and development of children and has a high rate of disability and mortality. SJIA may result from genetic, infectious, or autoimmune factors since the precise source of the disease is unknown. Our study aims to develop a genetic-based diagnostic model to explore the identification of SJIA at the genetic level.

METHODS

The gene expression dataset of peripheral blood mononuclear cell (PBMC) samples from SJIA was collected from the Gene Expression Omnibus (GEO) database. Then, three GEO datasets (GSE11907-GPL96, GSE8650-GPL96 and GSE13501) were merged and used as a training dataset, which included 125 SJIA samples and 92 health samples. GSE7753 was used as a validation dataset. The limma method was used to screen differentially expressed genes (DEGs). Feature selection was performed using Lasso, random forest (RF)-recursive feature elimination (RFE) and RF classifier.

RESULTS

We finally identified 4 key genes (ALDH1A1, CEACAM1, YBX3 and SLC6A8) that were essential to distinguish SJIA from healthy samples. And we combined the 4 key genes and performed a grid search as well as 10-fold cross-validation with 5 repetitions to finally identify the RF model with optimal mtry. The mean area under the curve (AUC) value for 5-fold cross-validation was greater than 0.95. The model's performance was then assessed once more using the validation dataset, and an AUC value of 0.990 was obtained. All of the above AUC values demonstrated the strong robustness of the SJIA diagnostic model.

CONCLUSIONS

We successfully developed a new SJIA diagnostic model that can be used for a novel aid in the identification of SJIA. In addition, the identification of 4 key genes that may serve as potential biomarkers for SJIA provides new insights to further understand the mechanisms of SJIA.

Single-Cell Transcriptional Analysis of Lamina Propria Lymphocytes in the Jejunum Reveals Innate Lymphoid Cell-like Cells in Pigs

Pigs are the most suitable model to study various therapeutic strategies and drugs for human beings, although knowledge about cell type-specific transcriptomes and heterogeneity is poorly available. Through single-cell RNA sequencing and flow cytometry analysis of the types in the jejunum of pigs, we found that innate lymphoid cells (ILCs) existed in the lamina propria lymphocytes (LPLs) of the jejunum. Then, through flow sorting of live/dead-lineage (Lin)-CD45+ cells and single-cell RNA sequencing, we found that ILCs in the porcine jejunum were mainly ILC3s, with a small number of NK cells, ILC1s, and ILC2s. ILCs coexpressed IL-7Rα, ID2, and other genes and differentially expressed RORC, GATA3, and other genes but did not express the CD3 gene. ILC3s can be divided into four subgroups, and genes such as CXCL8, CXCL2, IL-22, IL-17, and NCR2 are differentially expressed. To further detect and identify ILC3s, we verified the classification of ILCs in the porcine jejunum subgroup and the expression of related hallmark genes at the protein level by flow cytometry. For systematically characterizing ILCs in the porcine intestines, we combined our pig ILC dataset with publicly available human and mice ILC data and identified that the human and pig ILCs shared more common features than did those mouse ILCs in gene signatures and cell states. Our results showed in detail for the first time (to our knowledge) the gene expression of porcine jejunal ILCs, the subtype classification of ILCs, and the markers of various ILCs, which provide a basis for an in-depth exploration of porcine intestinal mucosal immunity.

The molecular and cellular anatomy of a fetal programming defect - the impact of low protein diet on the developing kidney

Low nephron number correlates with the development of hypertension and chronic kidney disease later in life. While intrauterine growth restriction caused by maternal low protein diet (LPD) is thought to be a significant cause of reduced nephron endowment in impoverished communities, its influence on the cellular and molecular processes which drive nephron formation are poorly understood. We conducted a comprehensive characterization of the impact of LPD on kidney development using tomographic and confocal imaging to quantify changes in branching morphogenesis and the cellular and morphological features of nephrogenic niches across development. These analyses were paired with single-cell RNA sequencing to dissect the transcriptional changes that LPD imposes during renal development. Differences in the expression of genes involved in metabolism were identified in most cell types we analyzed, yielding imbalances and shifts in cellular energy production. We further demonstrate that LPD impedes branching morphogenesis and significantly reduces the number of pretubular aggregates - the initial precursors to nephron formation. The most striking observation was that LPD changes the developmental trajectory of nephron progenitor cells, driving the formation of a partially committed cell population which likely reflects a failure of cells to commit to nephron formation and which ultimately reduces endowment. This unique profile of a fetal programming defect demonstrates that low nephron endowment arises from the pleiotropic impact of changes in branching morphogenesis and nephron progenitor cell commitment, the latter of which highlights a critical role for nutrition in regulating the cell fate decisions underpinning nephron endowment.

Significance Statement

While a mother's diet and behavior can negatively impact the number of nephrons in the kidneys of her offspring, the root cellular and molecular drivers of these deficits have not been rigorously explored. In this study we use advanced imaging and gene expression analysis in mouse models to define how a maternal low protein diet, analogous to that of impoverished communities, results in reduced nephron endowment. We find that low protein diet has pleiotropic effects on metabolism and the normal programs of gene expression. These profoundly impact the process of branching morphogenesis necessary to establish niches for nephron generation and change cell behaviors which regulate how and when nephron progenitor cells commit to differentiation.

The 3'‑untranslated region of XB130 regulates its mRNA stability and translational efficiency in non‑small cell lung cancer cells

Silencing XB130 inhibits cell proliferation and epithelial-mesenchymal transition in non-small cell lung cancer (NSCLC), suggesting that downregulating XB130 expression may impede NSCLC progression. However, the molecular mechanism underlying the regulation of XB130 expression remains unclear. In the present study, the role of the 3'-untranslated region (3'-UTR) in the regulation of XB130 expression was investigated. Recombinant psiCHECK-2 vectors with wild-type, truncated, or mutant XB130 3'-UTR were constructed, and the effects of these insertions on reporter gene expression were examined using a dual-luciferase reporter assay and reverse transcription-quantitative PCR. Additionally, candidate proteins that regulated XB130 expression by binding to critical regions of the XB130 3'-UTR were screened for using an RNA pull-down assay, followed by mass spectrometry and western blotting. The results revealed that insertion of the entire XB130 3'-UTR (1,218 bp) enhanced reporter gene expression. Positive regulatory elements were primarily found in nucleotides 113-989 of the 3'-UTR, while negative regulatory elements were found in the 1-112 and 990-1,218 regions of the 3'-UTR. Deletion analyses identified nucleotides 113-230 and 503-660 of the 3'-UTR as two major fragments that likely promote XB130 expression by increasing mRNA stability and translation rate. Additionally, a U-rich element in the 970-1,053 region of the 3'-UTR was identified as a negative regulatory element that inhibited XB130 expression by suppressing translation. Furthermore, seven candidate proteins that potentially regulated XB130 expression by binding to the 113-230, 503-660, and 970-1,053 regions of the 3'-UTR were identified, shedding light on the regulatory mechanism of XB130 expression. Collectively, these results suggested that complex sequence integrations in the mRNA 3'-UTR variably affected XB130 expression in NSCLC cells.

Protein-level mutant p53 reporters identify druggable rare precancerous clones in noncancerous tissues

Detecting and targeting precancerous cells in noncancerous tissues is a major challenge for cancer prevention. Massive stabilization of mutant p53 (mutp53) proteins is a cancer-specific event that could potentially mark precancerous cells, yet in vivo protein-level mutp53 reporters are lacking. Here we developed two transgenic protein-level mutp53 reporters, p53R172H-Akaluc and p53-mCherry, that faithfully mimic the dynamics and function of mutp53 proteins in vivo. Using these reporters, we identified and traced rare precancerous clones in deep noncancerous tissues in various cancer models. In classic mutp53-driven thymic lymphoma models, we found that precancerous clones exhibit broad chromosome number variations, upregulate precancerous stage-specific genes such as Ybx3 and enhance amino acid transport and metabolism. Inhibiting amino acid transporters downstream of Ybx3 at the early but not late stage effectively suppresses tumorigenesis and prolongs survival. Together, these protein-level mutp53 reporters reveal undercharacterized features and vulnerabilities of precancerous cells during early tumorigenesis, paving the way for precision cancer prevention.

MicroRNA profiling of royal jelly extracellular vesicles and their potential role in cell viability and reversing cell apoptosis

MiRNAs are small non-coding RNA molecules that play important regulatory roles in diverse biological processes. Royal jelly, a milky-white substance produced by nurse honeybees (Apis mellifera), is the primary food of queen bees and plays a crucial role in their development. However, little is known about the microRNA (miRNAs) content of royal jelly and their potential functions. In this study, we isolated extracellular vesicles from the royal jelly of 36 samples through sequential centrifugation and targeted nanofiltration and performed high-throughput sequencing to identify and quantify the miRNA content of honeybee royal jelly extracellular vesicles (RJEVs). We found a total of 29 known mature miRNAs and 17 novel miRNAs. Through bioinformatic analysis, we identified several potential target genes of the miRNAs present in royal jelly, including those involved in developmental processes and cell differentiation. To investigate the potential roles of RJEVs in cell viability, RJEVs were supplemented to apoptotic porcine kidney fibroblasts induced by ethanol 6% exposure for 30 min. TUNEL assay showed a significant reduction in the apoptosis percentage after RJEV supplementation when compared with the non-supplemented control group. Moreover, the wound healing assay performed on the apoptotic cells showed a rapid healing capacity of RJEV-supplemented cells compared to the control group. We observed a significant reduction in the expression of the miRNA target genes such as FAM131B, ZEB1, COL5A1, TRIB2, YBX3, MAP2, CTNNA1, and ADAMTS9 suggesting that RJEVs may regulate the target gene expression associated with cellular motility and cell viability. Moreover, RJEVs reduced the expression of apoptotic genes (CASP3, TP53, BAX, and BAK), while significantly increasing the expression of anti-apoptotic genes (BCL2 and BCL-XL). Our findings provide the first comprehensive analysis of the miRNA content of RJEVs and suggest a potential role for these vesicles in the regulation of gene expression and cell survival as well as augmenting cell resurrection or anastasis.

Diversity and Differential Expression of MicroRNAs in the Human Skeletal Muscle with Distinct Fiber Type Composition

The ratio of fast- and slow-twitch fibers in human skeletal muscle is variable and largely determined by genetic factors. In this study, we investigated the contribution of microRNA (miRNA) in skeletal muscle fiber type composition. The study involved biopsy samples of the vastus lateralis muscle from 24 male participants with distinct fiber type ratios. The miRNA study included samples from five endurance athletes and five power athletes with the predominance of slow-twitch (61.6-72.8%) and fast-twitch (69.3-80.7%) fibers, respectively. Total and small RNA were extracted from tissue samples. Total RNA sequencing (N = 24) revealed 352 differentially expressed genes between the groups with the predominance of fast- and slow-twitch muscle fibers. Small RNA sequencing showed upregulation of miR-206, miR-501-3p and miR-185-5p, and downregulation of miR-499a-5p and miR-208-5p in the group of power athletes with fast-twitch fiber predominance. Two miRtronic miRNAs, miR-208b-3p and miR-499a-5p, had strong correlations in expression with their host genes (MYH7 and MYH7B, respectively). Correlations between the expression of miRNAs and their experimentally validated messenger RNA (mRNA) targets were calculated, and 11 miRNA-mRNA interactions with strong negative correlations were identified. Two of them belonged to miR-208b-3p and miR-499a-5p, indicating their regulatory links with the expression of CDKN1A and FOXO4, respectively.

A haplotype-resolved genome assembly of the Nile rat facilitates exploration of the genetic basis of diabetes

Background

The Nile rat (Avicanthis niloticus) is an important animal model because of its robust diurnal rhythm, a cone-rich retina, and a propensity to develop diet-induced diabetes without chemical or genetic modifications. A closer similarity to humans in these aspects, compared to the widely used Mus musculus and Rattus norvegicus models, holds the promise of better translation of research findings to the clinic.

Results

We report a 2.5 Gb, chromosome-level reference genome assembly with fully resolved parental haplotypes, generated with the Vertebrate Genomes Project (VGP). The assembly is highly contiguous, with contig N50 of 11.1 Mb, scaffold N50 of 83 Mb, and 95.2% of the sequence assigned to chromosomes. We used a novel workflow to identify 3613 segmental duplications and quantify duplicated genes. Comparative analyses revealed unique genomic features of the Nile rat, including some that affect genes associated with type 2 diabetes and metabolic dysfunctions. We discuss 14 genes that are heterozygous in the Nile rat or highly diverged from the house mouse.

Conclusions

Our findings reflect the exceptional level of genomic resolution present in this assembly, which will greatly expand the potential of the Nile rat as a model organism.

Transcriptome analysis of immune cells from Behçet's syndrome patients: the importance of IL-17-producing cells and antigen-presenting cells in the pathogenesis of Behçet's syndrome

Background

Behçet’s syndrome (BS) is an immune-mediated disease characterized by recurrent oral ulcers, genital ulcers, uveitis, and skin symptoms. HLA-B51, as well as other genetic polymorphisms, has been reported to be associated with BS; however, the pathogenesis of BS and its relationship to genetic risk factors still remain unclear. To address these points, we performed immunophenotyping and transcriptome analysis of immune cells from BS patients and healthy donors.

Methods

ImmuNexUT is a comprehensive database consisting of RNA sequencing data and eQTL database of immune cell subsets from patients with immune-mediated diseases and healthy donors, and flow cytometry data and transcriptome data from 23 BS patients and 28 healthy donors from the ImmuNexUT study were utilized for this study. Differential gene expression analysis and weighted gene co-expression network analysis (WGCNA) were performed to identify genes associated with BS and clinical features of BS. eQTL database was used to assess the relationship between genetic risk factors of BS with those genes.

Results

The frequency of Th17 cells was increased in BS patients, and transcriptome analysis of Th17 cells suggested the activation of the NFκB pathway in Th17 cells of BS patients. Next, WGCNA was used to group genes into modules with similar expression patterns in each subset. Modules of antigen-presenting cells were associated with BS, and pathway analysis suggested the activation of antigen-presenting cells of BS patients. Further examination of genes in BS-associated modules indicated that the expression of YBX3, a member of a plasmacytoid dendritic cell (pDC) gene module associated with BS, is influenced by a BS risk polymorphism, rs2617170, in pDCs, suggesting that YBX3 may be a key molecule connecting genetic risk factors of BS with disease pathogenesis. Furthermore, pathway analysis of modules associated with HLA-B51 indicated that the association of IL-17-associated pathways in memory CD8⁺ T cells with HLA-B51; therefore, IL-17-producing CD8⁺ T cells, Tc17 cells, may play a critical role in BS.

Conclusions

Various cells including CD4⁺ T cells, CD8⁺ T cells, and antigen-presenting cells are important in the pathogenesis of BS. Tc17 cells and YBX3 may be potential therapeutic targets in BS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13075-022-02867-x.

RNA-Binding Proteins in the Regulation of Adipogenesis and Adipose Function

The obesity epidemic represents a critical public health issue worldwide, as it is a vital risk factor for many diseases, including type 2 diabetes (T2D) and cardiovascular disease. Obesity is a complex disease involving excessive fat accumulation. Proper adipose tissue accumulation and function are highly transcriptional and regulated by many genes. Recent studies have discovered that post-transcriptional regulation, mainly mediated by RNA-binding proteins (RBPs), also plays a crucial role. In the lifetime of RNA, it is bound by various RBPs that determine every step of RNA metabolism, from RNA processing to alternative splicing, nucleus export, rate of translation, and finally decay. In humans, it is predicted that RBPs account for more than 10% of proteins based on the presence of RNA-binding domains. However, only very few RBPs have been studied in adipose tissue. The primary aim of this paper is to provide an overview of RBPs in adipogenesis and adipose function. Specifically, the following best-characterized RBPs will be discussed, including HuR, PSPC1, Sam68, RBM4, Ybx1, Ybx2, IGF2BP2, and KSRP. Characterization of these proteins will increase our understanding of the regulatory mechanisms of RBPs in adipogenesis and provide clues for the etiology and pathology of adipose-tissue-related diseases.

Short-Term Blockade of Pro-Inflammatory Alarmin S100A9 Favorably Modulates Left Ventricle Proteome and Related Signaling Pathways Involved in Post-Myocardial Infarction Recovery

Prognosis after myocardial infarction (MI) varies greatly depending on the extent of damaged area and the management of biological processes during recovery. Reportedly, the inhibition of the pro-inflammatory S100A9 reduces myocardial damage after MI. We hypothesize that a S100A9 blockade induces changes of major signaling pathways implicated in post-MI healing. Mass spectrometry-based proteomics and gene analyses of infarcted mice left ventricle were performed. The S100A9 blocker (ABR-23890) was given for 3 days after coronary ligation. At 3 and 7 days post-MI, ventricle samples were analyzed versus control and Sham-operated mice. Blockade of S100A9 modulated the expressed proteins involved in five biological processes: leukocyte cell–cell adhesion, regulation of the muscle cell apoptotic process, regulation of the intrinsic apoptotic signaling pathway, sarcomere organization and cardiac muscle hypertrophy. The blocker induced regulation of 36 proteins interacting with or targeted by the cellular tumor antigen p53, prevented myocardial compensatory hypertrophy, and reduced cardiac markers of post-ischemic stress. The blockade effect was prominent at day 7 post-MI when the quantitative features of the ventricle proteome were closer to controls. Blockade of S100A9 restores key biological processes altered post-MI. These processes could be valuable new pharmacological targets for the treatment of ischemic heart. Mass spectrometry data are available via ProteomeXchange with identifier PXD033683.

Identification of crucial circRNAs in skeletal muscle during chicken embryonic development

Background

Chicken provides humans with a large amount of animal protein every year, in which skeletal muscle plays a leading role. The embryonic skeletal muscle development determines the number of muscle fibers and will affect the muscle production of chickens. CircRNAs are involved in a variety of important biological processes, including muscle development. However, studies on circRNAs in the chicken embryo muscle development are still lacking.

Results

In the study, we collected chicken leg muscles at 14 and 20-day embryo ages both in the fast- and slow-growing groups for RNA-seq. We identified 245 and 440 differentially expressed (DE) circRNAs in the comparison group F14vsF20 and S14vsS20 respectively. GO enrichment analysis for the host genes of DE circRNAs showed that biological process (BP) terms in the top 20 related to growth in F14vsF20 were found such as positive regulation of transcription involved in G1/S phase of mitotic cell cycle, multicellular organismal macromolecule metabolic process, and multicellular organismal metabolic process. In group S14vsS20, we also found some BP terms associated with growth in the top 20 including actomyosin structure organization, actin cytoskeleton organization and myofibril assembly. A total of 7 significantly enriched pathways were obtained, containing Adherens junction and Tight junction. Further analysis of those pathways found three crucial host genes MYH9, YBX3, IGF1R in both fast- and slow-growing groups, three important host genes CTNNA3, AFDN and CREBBP only in the fast-growing group, and six host genes FGFR2, ACTN2, COL1A2, CDC42, DOCK1 and MYL3 only in the slow-growing group. In addition, circRNA-miRNA network also revealed some key regulation pairs such as novel_circ_0007646-miR-1625-5p, novel_circ_0007646-miR-1680-5p, novel_circ_0008913-miR-148b-5p, novel_circ_0008906-miR-148b-5p and novel_circ_0001640-miR-1759-3p.

Conclusions

Comprehensive analysis of circRNAs and their targets would contribute to a better understanding of the molecular mechanisms in poultry skeletal muscle and it also plays an important guiding role in the next research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08588-4.

A novel 27-bp indel in the intron region of the YBX3 gene is associated with growth traits in chickens

1. The DNA/RNA binding protein YBX3 is associated with gene transcription, DNA repair, and the progression of various diseases and is highly conserved from bacteria to humans.2. The following experiment found a 27-bp insertion/deletion polymorphism in the intron region of the YBX3 gene through resequencing. In cross-designed, F₂ resource groups, the indel was significantly associated with broiler weight and body size at 0, 2, 4, 6, 8, 10 and 12 weeks of age and several other traits (semi evisceration weight (SEW), evisceration weight (EW), semi evisceration rate (SER), evisceration rate (ER), head weight (HW), claw weight (CLW), wing weight (DWW), gizzard weight (GW), pancreas weight (PW), chest muscle weight (CMW), leg weight (LW), leg muscle weight (LMW), shedding weight (SW), carcass weight (CW) and pectoral area (PA)) (P<0.05).3. The insertion-insertion (II) genotype was significantly associated with the greatest growth traits and carcass traits, whereas the values associated with the insertion-deletion (ID) genotype were the lowest in the F₂ reciprocal cross chickens.4. The mutation sites were genotyped in 3611 individuals from 13 different chicken breeds and cross-designed F₂ resource groups. The II genotype is the most important in commercial broilers, and the I allele frequency observed in these breeds was relatively high. However, there is still considerable potential in breeding dual-purpose chickens and commercial laying hens.5. The mRNA expression of the YBX3 gene in tissues from different breeds and developmental stages demonstrated that the 27-bp indel may affect the entire development process of poultry by affecting muscle development. These findings are beneficial for elucidating the function of the YBX3 gene and facilitating enhanced reproduction in the chicken industry.

Loss of CLDN5 in podocytes deregulates WIF1 to activate WNT signaling and contributes to kidney disease

Although mature podocytes lack tight junctions, tight junction integral membrane protein claudin-5 (CLDN5) is predominantly expressed on plasma membranes of podocytes under normal conditions. Using podocyte-specific Cldn5 knockout mice, we identify CLDN5 as a crucial regulator of podocyte function and reveal that Cldn5 deletion exacerbates podocyte injury and proteinuria in a diabetic nephropathy mouse model. Mechanistically, CLDN5 deletion reduces ZO1 expression and induces nuclear translocation of ZONAB, followed by transcriptional downregulation of WNT inhibitory factor-1 (WIF1) expression, which leads to activation of WNT signaling pathway. Podocyte-derived WIF1 also plays paracrine roles in tubular epithelial cells, as evidenced by the finding that animals with podocyte-specific deletion of Cldn5 or Wif1 have worse kidney fibrosis after unilateral ureteral obstruction than littermate controls. Systemic delivery of WIF1 suppresses the progression of diabetic nephropathy and ureteral obstruction-induced renal fibrosis. These findings establish a function for podocyte CLDN5 in restricting WNT signaling in kidney.

Claudin-5 is a tight junction integral membrane protein, but it is also expressed in mature podocytes which lack tight junctions. Here the authors report that podocyte claudin-5 regulates WNT signaling activity by modulating WIF1 expression, and its downregulation contributes to kidney disease progression in mice.

Role of Y-Box Binding Proteins in Ontogenesis

Y-box binding proteins (YB proteins) are multifunctional DNA/RNA-binding proteins capable of regulating gene expression at multiple levels. At present, the most studied function of these proteins is the regulation of protein synthesis. Special attention in this review has been paid to the role of YB proteins in the control of mRNA translation and stability at the earliest stages of organism formation, from fertilization to gastrulation. Furthermore, the functions of YB proteins in the formation of germ cells, in which they accumulate in large amounts, are summarized. The review then discusses the contribution of YB proteins to the regulation of gene expression during the differentiation of various types of somatic cells. Finally, future directions in the study of YB proteins and their role in ontogenesis are considered.

Single-cell analysis identifies dynamic gene expression networks that govern B cell development and transformation

Integration of external signals and B-lymphoid transcription factor activities organise B cell lineage commitment through alternating cycles of proliferation and differentiation, producing a diverse repertoire of mature B cells. We use single-cell transcriptomics/proteomics to identify differentially expressed gene networks across B cell development and correlate these networks with subtypes of B cell leukemia. Here we show unique transcriptional signatures that refine the pre-B cell expansion stages into pre-BCR-dependent and pre-BCR-independent proliferative phases. These changes correlate with reciprocal changes in expression of the transcription factor EBF1 and the RNA binding protein YBX3, that are defining features of the pre-BCR-dependent stage. Using pseudotime analysis, we further characterize the expression kinetics of different biological modalities across B cell development, including transcription factors, cytokines, chemokines, and their associated receptors. Our findings demonstrate the underlying heterogeneity of developing B cells and characterise developmental nodes linked to B cell transformation.

Transcriptional Alterations in Dorsolateral Prefrontal Cortex and Nucleus Accumbens Implicate Neuroinflammation and Synaptic Remodeling in Opioid Use Disorder

BACKGROUND

Prevalence rates of opioid use disorder (OUD) have increased dramatically, accompanied by a surge of overdose deaths. While opioid dependence has been extensively studied in preclinical models, an understanding of the biological alterations that occur in the brains of people who chronically use opioids and who are diagnosed with OUD remains limited. To address this limitation, RNA sequencing was conducted on the dorsolateral prefrontal cortex and nucleus accumbens, regions heavily implicated in OUD, from postmortem brains in subjects with OUD.

METHODS

We performed RNA sequencing on the dorsolateral prefrontal cortex and nucleus accumbens from unaffected comparison subjects (n = 20) and subjects diagnosed with OUD (n = 20). Our transcriptomic analyses identified differentially expressed transcripts and investigated the transcriptional coherence between brain regions using rank-rank hypergeometric orderlap. Weighted gene coexpression analyses identified OUD-specific modules and gene networks. Integrative analyses between differentially expressed transcripts and genome-wide association study datasets using linkage disequilibrium scores assessed the genetic liability of psychiatric-related phenotypes in OUD.

RESULTS

Rank-rank hypergeometric overlap analyses revealed extensive overlap in transcripts between the dorsolateral prefrontal cortex and nucleus accumbens in OUD, related to synaptic remodeling and neuroinflammation. Identified transcripts were enriched for factors that control proinflammatory cytokine, chondroitin sulfate, and extracellular matrix signaling. Cell-type deconvolution implicated a role for microglia as a potential driver for opioid-induced neuroplasticity. Linkage disequilibrium score analysis suggested genetic liabilities for risky behavior, attention-deficit/hyperactivity disorder, and depression in subjects with OUD.

CONCLUSIONS

Overall, our findings suggest connections between the brain's immune system and opioid dependence in the human brain.

Characterization of Tumor-Associated Macrophages and the Immune Microenvironment in Limited-Stage Neuroendocrine-High and -Low Small Cell Lung Cancer

Simple Summary

To date, the therapeutic strategy and guidelines in small cell lung cancer (SCLC) are based on cancer cell-related attributes with no biomarker used in the clinical practice. In the present study, using RNAseq and IHC, we aim to characterize in the frontline the latest biomarkers of tumor-associated macrophages (TAMs), myeloid-derived suppressor cells (MDSC) and related critical elements, regulating the anti-tumor immune response. Accordingly, we extensively evaluated the TME associations in primary tumors and matched lymph node metastases in different tumor compartments (stroma and tumor nests) and neuroendocrine (NE) subtypes in limited-stage SCLC. We show the RNA gene enrichment of the most critical molecular pathways based on the Gene Ontology (GO) iteration system using thorough bioinformatics analysis to identify new molecular targets in distinct NE subtypes.

Abstract

This study aims to characterize tumor-infiltrating macrophages (TAMs), myeloid-derived suppressor cells (MDSC), and the related molecular milieu regulating anti-tumor immunity in limited-stage neuroendocrine (NE)-high and NE-low small cell lung cancer. Primary tumors and matched lymph node (LN) metastases of 32 resected, early-stage SCLC patients were analyzed by immunohistochemistry (IHC) with antibodies against pan-macrophage marker CD68, M2-macrophage marker CD163, and MDSC marker CD33. Area-adjusted cell counting on TMAs showed that TAMs are the most abundant cell type in the TME, and their number in tumor nests exceeds the number of CD3 + T-cells (64% vs. 38% in NE-low and 71% vs. 18% in NE-high). Furthermore, the ratio of CD163-expressing M2-polarized TAMs in tumor nests was significantly higher in NE-low vs. NE-high tumors (70% vs. 31%). TAM density shows a strong positive correlation with CD45 and CD3 in tumor nests, but not in the stroma. fGSEA analysis on a targeted RNAseq oncological panel of 2560 genes showed that NE-high tumors exhibited increased enrichment in pathways related to cell proliferation, whereas in NE-low tumors, immune response pathways were significantly upregulated. Interestingly, we identified a subset of NE-high tumors representing an immune-oasis phenotype, but with a different gene expression profile compared to NE-low tumors. In contrast, we found that a limited subgroup of NE-low tumors is immune-deserted and express distinct cellular pathways from NE-high tumors. Furthermore, we identified potential molecular targets based on our expression data in NE-low and immune-oasis tumor subsets, including CD70, ANXA1, ITGB6, TP63, IFI27, YBX3 and CXCR2.

Targeting SLC1A5 and SLC3A2/SLC7A5 as a Potential Strategy to Strengthen Anti-Tumor Immunity in the Tumor Microenvironment

Cancer cells are metabolically vigorous and are superior in the uptake of nutrients and in the release of the tumor microenvironment (TME)-specific metabolites. They create an acidic, hypoxic, and nutrient-depleted TME that makes it difficult for the cytotoxic immune cells to adapt to the metabolically hostile environment. Since a robust metabolism in immune cells is required for optimal anti-tumor effector functions, the challenges caused by the TME result in severe defects in the invasion and destruction of the established tumors. There have been many recent developments in NK and T cell-mediated immunotherapy, such as engineering them to express chimeric antigen receptors (CARs) to enhance tumor-recognition and infiltration. However, to defeat the tumor and overcome the limitations of the TME, it is essential to fortify these novel therapies by improving the metabolism of the immune cells. One potential strategy to enhance the metabolic fitness of immune cells is to upregulate the expression of nutrient transporters, specifically glucose and amino acid transporters. In particular, the amino acid transporters SLC1A5 and SLC7A5 as well as the ancillary subunit SLC3A2, which are required for efficient uptake of glutamine and leucine respectively, could strengthen the metabolic capabilities and effector functions of tumor-directed CAR-NK and T cells. In addition to enabling the influx and efflux of essential amino acids through the plasma membrane and within subcellular compartments such as the lysosome and the mitochondria, accumulating evidence has demonstrated that the amino acid transporters participate in sensing amino acid levels and thereby activate mTORC1, a master metabolic regulator that promotes cell metabolism, and induce the expression of c-Myc, a transcription factor essential for cell growth and proliferation. In this review, we discuss the regulatory pathways of these amino acid transporters and how we can take advantage of these processes to strengthen immunotherapy against cancer.

Towards the mechanism(s) of YB-3 synthesis regulation by YB-1

Y-box binding proteins are members of the family of proteins containing the evolutionarily conserved cold shock domain. Their cellular functions are quite diverse, including transcription and translation regulation, participation in pre-mRNA splicing, mRNA stabilization and packaging into mRNPs, involvement in DNA repair, and some others. To date, we know little about the plausible functional interchangeability of Y-box binding proteins. Our previous finding was that in YB-1-null HEK293T cells the synthesis of YB-3 is enhanced, thus enabling YB-3 to interact with a larger set of mRNAs and compensate for the YB-1 absence. We suggested the existence of a mechanism of YB-3 synthesis regulation by its paralog, YB-1. Here we demonstrate that YB-1 participates in the translational control and stabilization of YB-3 mRNA through untranslated regions of YB-3 mRNA.

Identification of novel mRNA isoforms associated with meat tenderness using RNA sequencing data in beef cattle

The Warner-Bratzler shear force (WBSF) and myofibrillar fragmentation index (MFI) are complementary methodologies used to measure beef tenderness. Longissimus thoracis samples from the 20 most extreme bulls (out of 80 bulls set) for WBSF (tender (n = 10) and tough (n = 10)) and MFI (high (n = 10) and low (n = 10)) traits were collected to perform transcriptomic analysis using RNA-Sequencing. All analysis were performed through CLC Genomics Workbench. A total of 39 and 27 transcripts for WBSF and MFI phenotypes were DE, respectively. The possible DE novel mRNA isoforms, for WBSF and MFI traits, are myosin encoders (e.g. MYL1 and MYL6). In addition, we identified potential mRNA isoforms related to genes affecting the speed fibers degradation during the meat aging process. The DE novel transcripts are transcripted by genes with biological functions related to oxidative process, energy production and striated muscle contraction. The results suggest that the identified mRNA isoforms could be used as potential candidate to select animals in order to improve meat tenderness.

Circulating Coding and Long Non-Coding RNAs as Potential Biomarkers of Idiopathic Pulmonary Fibrosis

BACKGROUND

Idiopathic Pulmonary Fibrosis (IPF) is a chronic degenerative disease with a median survival of 2-5 years after diagnosis. Therefore, IPF patient identification represents an important and challenging clinical issue. Current research is still searching for novel reliable non-invasive biomarkers. Therefore, we explored the potential use of long non-coding RNAs (lncRNAs) and mRNAs as biomarkers for IPF.

METHODS

We first performed a whole transcriptome analysis using microarray (n = 14: 7 Control, 7 IPF), followed by the validation of selected transcripts through qPCRs in an independent cohort of 95 subjects (n = 95: 45 Control, 50 IPF). Diagnostic performance and transcript correlation with functional/clinical data were also analyzed.

RESULTS

1059 differentially expressed transcripts were identified. We confirmed the downregulation of FOXF1 adjacent non-coding developmental regulatory RNA (FENDRR) lncRNA, hsa_circ_0001924 circularRNA, utrophin (UTRN) and Y-box binding protein 3 (YBX3) mRNAs. The two analyzed non-coding RNAs correlated with Forced Vital Capacity (FVC)% and Diffusing Capacity of the Lung for carbon monoxide (DLCO)% functional data, while coding RNAs correlated with smock exposure. All analyzed transcripts showed excellent performance in IPF identification with Area Under the Curve values above 0.87; the most outstanding one was YBX3: AUROC 0.944, CI 95% = 0.895-0.992, sensitivity = 90%, specificity = 88.9%, p-value = 1.02 × 10^-13.

CONCLUSIONS

This study has identified specific transcript signatures in IPF suggesting that validated transcripts and microarray data could be useful for the potential future identification of RNA molecules as non-invasive biomarkers for IPF.

Development and validation of a prognostic model for kidney renal clear cell carcinoma based on RNA binding protein expression

Dysregulated expression of RNA-binding proteins (RBPs) is strongly associated with the development and progression of multiple tumors. However, little is known about the role of RBPs in kidney renal clear cell carcinoma (KIRC). In this study, we examined RBP expression profiles using The Cancer Genome Atlas database and identified 133 RBPs that were differentially expressed in KIRC and non-tumor tissues. We then systematically analyzed the potential biological functions of these RBPs and established PPIs. Based on Lasso regression and Cox survival analyses, we constructed a risk model that could independently and accurately predict prognosis based on seven RBPs (NOL12, PABPC1L, RNASE2, RPL22L1, RBM47, OASL, and YBX3). Survival times were shorter in patients with high risk scores for cohorts stratified by different characteristics. Gene set enrichment analysis was also performed to further understand functional differences between high- and low-risk groups. Finally, we developed a clinical nomogram with a concordance index of 0.792 for estimating 3- and 5-year survival probabilities. Our results demonstrate that this risk model could potentially improve individualized diagnostic and therapeutic strategies.

Periphilin self-association underpins epigenetic silencing by the HUSH complex

Transcription of integrated DNA from viruses or transposable elements is tightly regulated to prevent pathogenesis. The Human Silencing Hub (HUSH), composed of Periphilin, TASOR and MPP8, silences transcriptionally active viral and endogenous transgenes. HUSH recruits effectors that alter the epigenetic landscape and chromatin structure, but how HUSH recognizes target loci and represses their expression remains unclear. We identify the physicochemical properties of Periphilin necessary for HUSH assembly and silencing. A disordered N-terminal domain (NTD) and structured C-terminal domain are essential for silencing. A crystal structure of the Periphilin-TASOR minimal core complex shows Periphilin forms an α-helical homodimer, bound by a single TASOR molecule. The NTD forms insoluble aggregates through an arginine/tyrosine-rich sequence reminiscent of low-complexity regions from self-associating RNA-binding proteins. Residues required for TASOR binding and aggregation were required for HUSH-dependent silencing and genome-wide deposition of repressive mark H3K9me3. The NTD was functionally complemented by low-complexity regions from certain RNA-binding proteins and proteins that form condensates or fibrils. Our work suggests the associative properties of Periphilin promote HUSH aggregation at target loci.

Comprehensive Analysis of LincRNAs in Classical and Basal-Like Subtypes of Pancreatic Cancer

Pancreatic ductal adenocarcinomas (PDAC) belong to the deadliest malignancies in the western world. Mutations in TP53 and KRAS genes along with some other frequent polymorphisms occur almost universally and are major drivers of tumour initiation. However, these mutations cannot explain the heterogeneity in therapeutic responses and differences in overall survival observed in PDAC patients. Thus, recent classifications of PDAC tumour samples have leveraged transcriptome-wide gene expression data to account for epigenetic, transcriptional and post-transcriptional mechanisms that may contribute to this deadly disease. Intriguingly, long intervening RNAs (lincRNAs) are a special class of long non-coding RNAs (lncRNAs) that can control gene expression programs on multiple levels thereby contributing to cancer progression. However, their subtype-specific expression and function as well as molecular interactions in PDAC are not fully understood yet. In this study, we systematically investigated the expression of lincRNAs in pancreatic cancer and its molecular subtypes using publicly available data from large-scale studies. We identified 27 deregulated lincRNAs that showed a significant different expression pattern in PDAC subtypes suggesting context-dependent roles. We further analyzed these lincRNAs regarding their common expression patterns. Moreover, we inferred clues on their functions based on correlation analyses and predicted interactions with RNA-binding proteins, microRNAs, and mRNAs. In summary, we identified several PDAC-associated lincRNAs of prognostic relevance and potential context-dependent functions and molecular interactions. Hence, our study provides a valuable resource for future investigations to decipher the role of lincRNAs in pancreatic cancer.

Y-Box Binding Proteins in mRNP Assembly, Translation, and Stability Control

Y-box binding proteins (YB proteins) are DNA/RNA-binding proteins belonging to a large family of proteins with the cold shock domain. Functionally, these proteins are known to be the most diverse, although the literature hardly offers any molecular mechanisms governing their activities in the cell, tissue, or the whole organism. This review describes the involvement of YB proteins in RNA-dependent processes, such as mRNA packaging into mRNPs, mRNA translation, and mRNA stabilization. In addition, recent data on the structural peculiarities of YB proteins underlying their interactions with nucleic acids are discussed.

A Targeted Multi-omic Analysis Approach Measures Protein Expression and Low-Abundance Transcripts on the Single-Cell Level

High-throughput single-cell RNA sequencing (scRNA-seq) has become a frequently used tool to assess immune cell heterogeneity. Recently, the combined measurement of RNA and protein expression was developed, commonly known as cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq). Acquisition of protein expression data along with transcriptome data resolves some of the limitations inherent to only assessing transcripts but also nearly doubles the sequencing read depth required per single cell. Furthermore, there is still a paucity of analysis tools to visualize combined transcript-protein datasets. Here, we describe a targeted transcriptomics approach that combines an analysis of over 400 genes with simultaneous measurement of over 40 proteins on 2 × 104 cells in a single experiment. This targeted approach requires only about one-tenth of the read depth compared to a whole-transcriptome approach while retaining high sensitivity for low abundance transcripts. To analyze these multi-omic datasets, we adapted one-dimensional soli expression by nonlinear stochastic embedding (One-SENSE) for intuitive visualization of protein-transcript relationships on a single-cell level.

Integrated Analysis of the Functions and Prognostic Values of RNA Binding Proteins in Lung Squamous Cell Carcinoma

Lung cancer is the leading cause of cancer-related deaths worldwide. Dysregulation of RNA binding proteins (RBPs) has been found in a variety of cancers and is related to oncogenesis and progression. However, the functions of RBPs in lung squamous cell carcinoma (LUSC) remain unclear. In this study, we obtained gene expression data and corresponding clinical information for LUSC from The Cancer Genome Atlas (TCGA) database, identified aberrantly expressed RBPs between tumors and normal tissue, and conducted a series of bioinformatics analyses to explore the expression and prognostic value of these RBPs. A total of 300 aberrantly expressed RBPs were obtained, comprising 59 downregulated and 241 upregulated RBPs. Functional enrichment analysis indicated that the differentially expressed RBPs were mainly associated with mRNA metabolic processes, RNA processing, RNA modification, regulation of translation, the TGF-beta signaling pathway, and the Toll-like receptor signaling pathway. Nine RBP genes (A1CF, EIF2B5, LSM1, LSM7, MBNL2, RSRC1, TRMU, TTF2, and ZCCHC5) were identified as prognosis-associated hub genes by univariate, least absolute shrinkage and selection operator (LASSO), Kaplan–Meier survival, and multivariate Cox regression analyses, and were used to construct the prognostic model. Further analysis demonstrated that high risk scores for patients were significantly related to poor overall survival according to the model. The area under the time-dependent receiver operator characteristic curve of the prognostic model was 0.712 at 3 years and 0.696 at 5 years. We also developed a nomogram based on nine RBP genes, with internal validation in the TCGA cohort, which showed a favorable predictive efficacy for prognosis in LUSC. Our results provide novel insights into the pathogenesis of LUSC. The nine-RBP gene signature showed predictive value for LUSC prognosis, with potential applications in clinical decision-making and individualized treatment.