RNA-seq Data Research Articles

A Novel Macrophage Subpopulation Conveys Increased Genetic Risk of Coronary Artery Disease.

Coronary artery disease (CAD), the leading cause of death worldwide, is influenced by both environmental and genetic factors. Although over 250 genetic risk loci have been identified through genome-wide association studies, the specific causal variants and their regulatory mechanisms are still largely unknown, particularly in disease-relevant cell types like macrophages. We utilized single-cell RNA-seq and single-cell multiomics approaches in primary human monocyte-derived macrophages to explore the transcriptional regulatory network involved in a critical pathogenic event of coronary atherosclerosis-the formation of lipid-laden foam cells. The relative genetic contribution to CAD was assessed by partitioning disease heritability across different macrophage subpopulations. Meta-analysis of single-cell RNA-seq data sets from 38 human atherosclerotic samples was conducted to provide high-resolution cross-referencing to macrophage subpopulations in vivo. We identified 18 782 cis-regulatory elements by jointly profiling the gene expression and chromatin accessibility of >5000 macrophages. Integration with CAD genome-wide association study data prioritized 121 CAD-related genetic variants and 56 candidate causal genes. We showed that CAD heritability was not uniformly distributed and was particularly enriched in the gene programs of a novel CD52-hi lipid-handling macrophage subpopulation. These CD52-hi macrophages displayed significantly less lipoprotein accumulation and were also found in human atherosclerotic plaques. We investigated the cis-regulatory effect of a risk variant rs10488763 on FDX1, implicating the recruitment of AP-1 and C/EBP-β in the causal mechanisms at this locus. Our results provide genetic evidence of the divergent roles of macrophage subsets in atherogenesis and highlight lipid-handling macrophages as a key subpopulation through which genetic variants operate to influence disease. These findings provide an unbiased framework for functional fine-mapping of genome-wide association study results using single-cell multiomics and offer new insights into the genotype-environment interactions underlying atherosclerotic disease.

PSMD9 promotes the malignant progression of hepatocellular carcinoma by interacting with c-Cbl to activate EGFR signaling and recycling

BackgroundMounting evidences shows that the ubiquitin‒proteasome pathway plays a pivotal role in tumor progression. The expression of 26S proteasome non-ATPase regulatory subunit 9 (PSMD9) is correlated with recurrence and radiotherapy resistance in several tumor types. However, the role and mechanism of PSMD9 in hepatocellular carcinoma (HCC) progression remain largely unclear.MethodsPSMD9 was identified as a prognosis-related biomarker for HCC based on analysis of clinical characteristics and RNA-seq data from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) and the JP Project of the International Cancer Genome Consortium (ICGC-LIRI-JP). PSMD9 expression was analyzed in cancer tissues and adjacent noncancerous tissues via immunohistochemistry and Western blotting. Multiple in vivo and in vitro experimental techniques (such as CCK-8, colony formation, EdU, and Transwell assays; flow cytometry; Western blotting; quantitative RT-PCR; Coimmunoprecipitation assay and immunofluorescence confocal imaging) were used to assess the functions of PSMD9 in the pathogenesis of HCC.ResultsWe found that the expression of PSMD9 was upregulated and associated with a poor prognosis in HCC patients. PSMD9 promoted HCC cell proliferation, migration, invasion and metastasis. Knockdown of PSMD9 significantly inhibited HCC cell proliferation by inducing G1/S cell cycle arrest and apoptosis. Mechanistically, we demonstrated that PSMD9 promoted HCC cell proliferation and metastasis via direct interaction with the E3 ubiquitin ligase c-Cbl, suppresses EGFR ubiquitination, influenced EGFR endosomal trafficking and degradation and subsequently activated ERK1/2 and Akt signaling. In addition, we showed that PSMD9 knockdown sensitized HCC cells to the tyrosine kinase inhibitor erlotinib in vitro and in vivo.ConclusionsCollectively, our results indicate that PSMD9 drives HCC progression and erlotinib resistance by suppressing c-Cbl mediated EGFR ubiquitination and therefore can be a potential therapeutic target for HCC.

Conserved role of hnRNPL in alternative splicing of epigenetic modifiers enables B cell activation.

The multifunctional RNA-binding protein hnRNPL is implicated in antibody class switching but its broader function in B cells is unknown. Here, we show that hnRNPL is essential for B cell activation, germinal center formation,and antibody responses. Upon activation, hnRNPL-deficient B cells show proliferation defects and increased apoptosis. Comparative analysis of RNA-seq data from activated B cells and another eight hnRNPL-depleted cell types reveals common effects on MYC and E2F transcriptional programs required for proliferation. Notably, while individual gene expression changes are cell type specific, several alternative splicing events affecting histone modifiers like KDM6A and SIRT1, are conserved across cell types. Moreover, hnRNPL-deficient B cells show global changes in H3K27me3 and H3K9ac. Epigenetic dysregulation after hnRNPL loss could underlie differential gene expression and upregulation of lncRNAs, and explain common and cell type-specific phenotypes, such as dysfunctional mitochondria and ROS overproduction in mouse B cells. Thus, hnRNPL is essential for the resting-to-activated B cell transition by regulating transcriptional programs and metabolism, at least in part through the alternative splicing of several histone modifiers.

Genome-wide identification and expression analysis of the Eriobotrya japonica TIFY gene family reveals its functional diversity under abiotic stress conditions

BackgroundPlant-specific TIFY proteins are widely found in terrestrial plants and play important roles in plant adversity responses. Although the genome of loquat at the chromosome level has been published, studies on the TIFY family in loquat are lacking. Therefore, the EjTIFY gene family was bioinformatically analyzed by constructing a phylogenetic tree, chromosomal localization, gene structure, and adversity expression profiling in this study.ResultsTwenty-six EjTIFY genes were identified and categorized into four subfamilies (ZML, JAZ, PPD, and TIFY) based on their structural domains. Twenty-four EjTIFY genes were irregularly distributed on 11 of the 17 chromosomes, and the remaining two genes were distributed in fragments. We identified 15 covariate TIFY gene pairs in the loquat genome, 13 of which were involved in large-scale interchromosomal segmental duplication events, and two of which were involved in tandem duplication events. Many abiotic stress cis-elements were widely present in the promoter region. Analysis of the Ka/Ks ratio showed that the paralogous homologs of the EjTIFY family were mainly subjected to purifying selection. Analysis of the RNA-seq data revealed that a total of five differentially expressed genes (DEGs) were expressed in the shoots under gibberellin treatment, whereas only one gene was significantly differentially expressed in the leaves; under both low-temperature and high-temperature stresses, there were significantly differentially expressed genes, and the EjJAZ15 gene was significantly upregulated under both low- and high-temperature stress. RNA-seq and qRT-PCR expression analysis under salt stress conditions revealed that EjJAZ2, EjJAZ4, and EjJAZ9 responded to salt stress in loquat plants, which promoted resistance to salt stress through the JA pathway. The response model of the TIFY genes in the jasmonic acid pathway under salt stress in loquat was systematically summarized.ConclusionsThese results provide a theoretical basis for exploring the characteristics and functions of additional EjTIFY genes in the future. This study also provides a theoretical basis for further research on breeding for salt stress resistance in loquat. RT-qPCR analysis revealed that the expression of one of the three EjTIFY genes increased and the expression of two decreased under salt stress conditions, suggesting that EjTIFY exhibited different expression patterns under salt stress conditions.

RFX6 haploinsufficiency predisposes to diabetes through impaired beta cell function.

Regulatory factor X 6 (RFX6) is crucial for pancreatic endocrine development and differentiation. The RFX6 variant p.His293LeufsTer7 is significantly enriched in the Finnish population, with almost 1:250 individuals as a carrier. Importantly, the FinnGen study indicates a high predisposition for heterozygous carriers to develop type 2 and gestational diabetes. However, the precise mechanism of this predisposition remains unknown. To understand the role of this variant in beta cell development and function, we used CRISPR technology to generate allelic series of pluripotent stem cells. We created two isogenic stem cell models: a human embryonic stem cell model; and a patient-derived stem cell model. Both were differentiated into pancreatic islet lineages (stem-cell-derived islets, SC-islets), followed by implantation in immunocompromised NOD-SCID-Gamma mice. Stem cell models of the homozygous variant RFX6-/- predictably failed to generate insulin-secreting pancreatic beta cells, mirroring the phenotype observed in Mitchell-Riley syndrome. Notably, at the pancreatic endocrine stage, there was an upregulation of precursor markers NEUROG3 and SOX9, accompanied by increased apoptosis. Intriguingly, heterozygous RFX6+/- SC-islets exhibited RFX6 haploinsufficiency (54.2% reduction in protein expression), associated with reduced beta cell maturation markers, altered calcium signalling and impaired insulin secretion (62% and 54% reduction in basal and high glucose conditions, respectively). However, RFX6 haploinsufficiency did not have an impact on beta cell number or insulin content. The reduced insulin secretion persisted after in vivo implantation in mice, aligning with the increased risk of variant carriers to develop diabetes. Our allelic series isogenic SC-islet models represent a powerful tool to elucidate specific aetiologies of diabetes in humans, enabling the sensitive detection of aberrations in both beta cell development and function. We highlight the critical role of RFX6 in augmenting and maintaining the pancreatic progenitor pool, with an endocrine roadblock and increased cell death upon its loss. We demonstrate that RFX6 haploinsufficiency does not affect beta cell number or insulin content but does impair function, predisposing heterozygous carriers of loss-of-function variants to diabetes. Ultra-deep bulk RNA-seq data for pancreatic differentiation stages 3, 5 and 7 of H1 RFX6 genotypes are deposited in the Gene Expression Omnibus database with accession code GSE234289. Original western blot images are deposited at Mendeley ( https://data.mendeley.com/datasets/g75drr3mgw/2 ).

GRouNdGAN: GRN-guided simulation of single-cell RNA-seq data using causal generative adversarial networks

We introduce GRouNdGAN, a gene regulatory network (GRN)-guided reference-based causal implicit generative model for simulating single-cell RNA-seq data, in silico perturbation experiments, and benchmarking GRN inference methods. Through the imposition of a user-defined GRN in its architecture, GRouNdGAN simulates steady-state and transient-state single-cell datasets where genes are causally expressed under the control of their regulating transcription factors (TFs). Training on six experimental reference datasets, we show that our model captures non-linear TF-gene dependencies and preserves gene identities, cell trajectories, pseudo-time ordering, and technical and biological noise, with no user manipulation and only implicit parameterization. GRouNdGAN can synthesize cells under new conditions to perform in silico TF knockout experiments. Benchmarking various GRN inference algorithms reveals that GRouNdGAN effectively bridges the existing gap between simulated and biological data benchmarks of GRN inference algorithms, providing gold standard ground truth GRNs and realistic cells corresponding to the biological system of interest.

RNA-Seq transcriptome profiling of immature grain wheat is a technique for understanding comparative modeling of baking quality

Improving the baking quality is a primary challenge in the wheat flour production value chain, as baking quality represents a crucial factor in determining its overall value. In the present study, we conducted a comparative RNA-Seq analysis on the high baking quality mutant “O-64.1.10” genotype and its low baking quality wild type "Omid" cultivar to recognize potential genes associated with bread quality. The cDNA libraries were constructed from immature grains that were 15 days post-anthesis, with an average of 16.24 and 18.97 million paired-end short-read sequences in the mutant and wild-type, respectively. A total number of 733 transcripts with differential expression were identified, 585 genes up-regulated and 188 genes down-regulated in the “O-64.1.10” genotype compared to the “Omid”. In addition, the families of HSF, bZIP, C2C2-Dof, B3-ARF, BES1, C3H, GRF, HB-HD-ZIP, PLATZ, MADS-MIKC, GARP-G2-like, NAC, OFP and TUB were appeared as the key transcription factors with specific expression in the “O-64.1.10” genotype. At the same time, pathways related to baking quality were identified through Kyoto Encyclopedia of Genes and Genomes. Collectively, we found that the endoplasmic network, metabolic pathways, secondary metabolite biosynthesis, hormone signaling pathway, B group vitamins, protein pathways, pathways associated with carbohydrate and fat metabolism, as well as the biosynthesis and metabolism of various amino acids, have a great deal of potential to play a significant role in the baking quality. Ultimately, the RNA-seq results were confirmed using quantitative Reverse Transcription PCR for some hub genes such as alpha-gliadin, low molecular weight glutenin subunit and terpene synthase (gibberellin) and as a resource for future study, 127 EST-SSR primers were generated using RNA-seq data.

Single cells and TRUST4 reveal immunological features of the HFRS transcriptome

The etiology of hemorrhagic fever with renal syndrome (HFRS) is significantly impacted by a variety of immune cells. Nevertheless, the existing techniques for sequencing peripheral blood T cell receptor (TCR) or B cell receptor (BCR) libraries in HFRS are constrained by both limitations and high costs. In this investigation, we utilized the computational tool TRUST4 to generate TCR and BCR libraries utilizing comprehensive RNA-seq data from peripheral blood specimens of HFRS patients. This facilitated the examination of clonality and diversity within immune libraries linked to the condition. Despite previous research on immune cell function, the underlying mechanisms remain intricate, and differential gene expression across immune cell types and cell-to-cell interactions within immune cell clusters have not been thoroughly explored. To address this gap, we performed clustering analysis on 11 cell subsets derived from raw single-cell RNA-seq data, elucidating characteristic changes in cell subset proportions under disease conditions. Additionally, we utilized CellChat, a tool for cell–cell communication analysis, to investigate the impact of MIF family, CD70 family, and GALECTIN family cytokines—known to be involved in cell communication—on immune cell subsets. Furthermore, hdWGCNA analysis identified core genes implicated in HFRS pathogenesis within T cells and B cells. Trajectory analysis revealed that most cell subsets were in a developmental stage, with high expression of transcription factors such as NFKB and JUN in Effector CD8+ T cells, as well as in Naive CD4+ T cells and Naive B cells. Our findings provide a comprehensive understanding of the dynamic changes in immune cells during HFRS pathogenesis, identifying specific V genes and J genes in TCR and BCR that contribute to advancing our knowledge of HFRS. These insights offer potential implications for the diagnosis and treatment of this autoimmune disease.

Integrated single-cell and bulk RNA-Seq analysis enhances prognostic accuracy of PD-1/PD-L1 immunotherapy response in lung adenocarcinoma through necroptotic anoikis gene signatures

In addition to presenting significant diagnostic and treatment challenges, lung adenocarcinoma (LUAD) is the most common form of lung cancer. Using scRNA-Seq and bulk RNA-Seq data, we identify three genes referred to as HMR, FAM83A, and KRT6A these genes are related to necroptotic anoikis-related gene expression. Initial validation, conducted on the GSE50081 dataset, demonstrated the model's ability to categorize LUAD patients into high-risk and low-risk groups with significant survival differences. This model was further applied to predict responses to PD-1/PD-L1 blockade therapies, utilizing the IMvigor210 and GSE78220 cohorts, and showed strong correlation with patient outcomes, highlighting its potential in personalized immunotherapy. Further, LUAD cell lines were analyzed using quantitative PCR (qPCR) and Western blot analysis to confirm their expression levels, further corroborating the model's relevance in LUAD pathophysiology. The mutation landscape of these genes was also explored, revealing their broad implication in various cancer types through a pan-cancer analysis. The study also delved into molecular subclustering, revealing distinct expression profiles and associations with different survival outcomes, emphasizing the model’s utility in precision oncology. Moreover, the diversity of immune cell infiltration, analyzed in relation to the necroptotic anoikis signature, suggested significant implications for immune evasion mechanisms in LUAD. While the findings present a promising stride towards personalized LUAD treatment, especially in immunotherapy, limitations such as the retrospective nature of the datasets and the need for larger sample sizes are acknowledged. Prospective clinical trials and further experimental research are essential to validate these findings and enhance the clinical applicability of our prognostic model.

Comprehensive analysis of competing endogenous RNA networks involved in the regulation of glycolysis in clear cell renal cell carcinoma

Abstract Objectives This study aims to elucidate the role of competing endogenous RNAs (ceRNAs), which are pivotal in the post-transcriptional regulation of cancer cells, in the glycolysis of clear cell renal cell carcinoma (ccRCC). Methods RNA-seq data from ccRCC samples and public datasets were subjected to differential expression analysis to identify the upregulated circular RNAs (circRNAs), long noncoding RNAs (lncRNAs), and mRNAs for constructing a ceRNA network. This network focused on the upregulated mRNAs linked to glycolysis and gluconeogenesis, which were verified at the proteome level using the CPTAC database. The ENCORI database was used to predict the lncRNA/circRNA-miRNA and miRNA-mRNA interactions, which formed a network visualized using Cytoscape. This study further examined the association of selected lncRNAs/circRNAs with biological processes and patient survival and explored lncRNA and mRNA expressions at the single-cell level. Results Three circRNAs and eight lncRNAs were identified to be regulators of four glycolysis pathway genes (HK2, LDHA, PFKP, and ALDOC) via 54 miRNAs. Notably, their expressions were aberrant at the transcriptome and proteome levels. These RNA elements were correlated with malignant and immune pathways, and several lncRNAs showed prognostic significance for ccRCC. Furthermore, single-cell RNA sequencing indicated the abnormal enrichment of these lncRNAs and mRNAs in malignant cells. Conclusions This study identified a ceRNA network that plays a critical role in the glycolysis of ccRCC. These findings highlight the importance of ceRNAs in cancer malignancy and their potential impact on patient prognosis.

The MYB family and their response to abiotic stress in ginger (Zingiber officinale Roscoe)

BackgroundZingiber officinale Roscoe, colloquially known as ginger, is a crop of significant medicinal and culinary value that frequently encounters adversity stemming from inhospitable environmental conditions. The MYB transcription factors have garnered recognition for their pivotal role in orchestrating a multitude of plant biological pathways. Nevertheless, the enumeration and characterization of the MYBs within Z. officinale Roscoe remains unknown. This study embarks on a genome-wide scrutiny of the MYB gene lineage in ginger, with the aim of cataloging all ZoMYB genes implicated in the biosynthesis of gingerols and curcuminoids, and elucidating their potential regulatory mechanisms in counteracting abiotic stress, thereby influencing ginger growth and development.ResultsIn this study, we identified an MYB gene family comprising 231 members in ginger genome. This ensemble comprises 74 singular-repeat MYBs (1R-MYB), 156 double-repeat MYBs (R2R3-MYB), and a solitary triple-repeat MYB (R1R2R3-MYB). Moreover, a comprehensive analysis encompassing the sequence features, conserved protein motifs, phylogenetic relationships, chromosome location, and gene duplication events of the ZoMYBs was conducted. We classified ZoMYBs into 37 groups, congruent with the number of conserved domains and gene structure analysis. Additionally, the expression profiles of ZoMYBs during development and under various stresses, including ABA, cold, drought, heat, and salt, were investigated in ginger utilizing both RNA-seq data and qRT-PCR analysis.ConclusionThis work provides a comprehensive understanding of the MYB family in ginger and lays the foundation for the future investigation of the potential functions of ZoMYB genes in ginger growth, development and abiotic stress tolerance of ginger.

A prediction of the CRNDE role by modulating NF-κB pathway in inflammatory bowel disease (IBD)

Long non-coding RNAs (lncRNAs) regulate multiple pathways and cellular mechanisms. Recent research has emphasized their involvement in the pathogenesis of complex diseases, such as Inflammatory Bowel Disease (IBD) which is characterized by chronic inflammation of the intestines. The two most common types of IBD are ulcerative colitis and Crohn's disease. CRNDE lncRNA was initially detected in colorectal cancer (CRC) and found to be involved in the tumorigenesis pathways. Further studies revealed the role of CRNDE in activating inflammation and promoting the release of inflammatory cytokines. This study utilizes the RNA-seq data analysis and bioinformatics tools to clarify the role of CRNDE in the IBD pathogenesis and confirms its expression in inflamed HT-29 and Caco-2 cell lines and also colonic and blood samples of UC patients and controls ex vivo. Based on our results, CRNDE was significantly upregulated in IBD samples compared to controls in RNA-seq data analysis and Real-time PCR of inflamed HT-29 cell line and colonic biopsies from UC patients. Additionally, predicted that its expression is positively correlated with the pro-inflammatory cytokines production. CRNDE interactions was investigated with several inflammation-related miRNAs and regulatory proteins computationally. Thus, CRNDE upregulation in the colon of IBD patients could be involved in IBD pathogenesis by promoting inflammatory pathways and targeting anti-inflammatory miRNAs.

Ring finger protein 138 inhibits transcription factor C/EBPα protein turnover leading to differentiation arrest in acute myeloid leukemia.

E3 ubiquitin ligase, ring finger protein 138 (RNF138) is involved in several biological processes; however, its role in myeloid differentiation or tumorigenesis remains unclear. RNAseq data from TNMplot showed that RNF138 mRNA levels are highly elevated in acute myeloid leukemia (AML) bone marrow samples as compared with bone marrow of normal volunteers. Here, we show that RNF138 serves as an E3 ligase for the tumor suppressor CCAAT/enhancer binding protein (C/EBPα) and promotes its degradation leading to myeloid differentiation arrest in AML. Wild-type RNF138 physically interacts with C/EBPα and promotes its ubiquitin-dependent proteasome degradation while a mutant RNF-138 deficient in ligase activity though interacts with C/EBPα, fails to down-regulate it. We show that RNF138 depletion enhances endogenous C/EBPα levels in peripheral blood mononuclear cells (PBMCs) isolated from healthy volunteers. Our data further shows that RNF138-mediated degradation of C/EBPα negatively affects its transactivation potential on its target genes. Furthermore, RNF138 overexpression inhibits all-trans-retinoic acid-induced differentiation of HL-60 cells whereas RNF138 RNAi enhances. In line with RNF138 inhibiting C/EBPα protein turnover, we also observed that RNF138 overexpression inhibited β-estradiol (E2)-induced C/EBPα driven granulocytic differentiation in C/EBPα inducible K562-p42C/EBPα-estrogen receptor cells. Furthermore, we also recapitulated these findings in PBMCs isolated from AML patients where depletion of RNF138 increased the expression of myeloid differentiation marker CD11b. These results suggest that RNF138 inhibits myeloid differentiation by targeting C/EBPα for proteasomal degradation and may provide a plausible mechanism for loss of C/EBPα expression often observed in myeloid leukemia. Also, targeting RNF138 may resolve differentiation arrest by restoring C/EBPα expression in AML.

Maternal gastrointestinal nematode infection alters hippocampal neuroimmunity, promotes synaptic plasticity, and improves resistance to direct infection in offspring

The developing brain is vulnerable to maternal bacterial and viral infections which induce strong inflammatory responses in the mother that are mimicked in the offspring brain, resulting in irreversible neurodevelopmental defects, and associated cognitive and behavioural impairments. In contrast, infection during pregnancy and lactation with the immunoregulatory murine intestinal nematode, Heligmosomoides bakeri, upregulates expression of genes associated with long-term potentiation (LTP) of synaptic networks in the brain of neonatal uninfected offspring, and enhances spatial memory in uninfected juvenile offspring. As the hippocampus is involved in spatial navigation and sensitive to immune events during development, here we assessed hippocampal gene expression, LTP, and neuroimmunity in 3-week-old uninfected offspring born to H. bakeri infected mothers. Further, as maternal immunity shapes the developing immune system, we assessed the impact of maternal H. bakeri infection on the ability of offspring to resist direct infection. In response to maternal infection, we found an enhanced propensity to induce LTP at Schaffer collateral synapses, consistent with RNA-seq data indicating accelerated development of glutamatergic synapses in uninfected offspring, relative to those from uninfected mothers. Hippocampal RNA-seq analysis of offspring of infected mothers revealed increased expression of genes associated with neurogenesis, gliogenesis, and myelination. Furthermore, maternal infection improved resistance to direct infection of H. bakeri in offspring, correlated with transfer of parasite-specific IgG1 to their serum. Hippocampal immunohistochemistry and gene expression suggest Th2/Treg biased neuroimmunity in offspring, recapitulating peripheral immunoregulation of H. bakeri infected mothers. These findings indicate maternal H. bakeri infection during pregnancy and lactation alters peripheral and neural immunity in uninfected offspring, in a manner that accelerates neural maturation to promote hippocampal LTP, and upregulates the expression of genes associated with neurogenesis, gliogenesis, and myelination.

PCDHGA10 as a Potential Biomarker of Lung Squamous Cell Carcinoma Based on Bioinformatics and Experimental Verification.

Procalcitonin gamma subfamily A, 10 (PCDHGA10) is a member of the procalcitonin gamma gene cluster, which is associated with neuronal synapse development. However, there are lack of studies on the role and potential prognostic value of PCDHGA10 in lung squamous cell carcinoma (LUSC). We analyzed the RNAseq data of PCDHGA10 to compare their expression differences. Then survival analysis, tumor microenvironment (TME) analysis, and mutation analysis were carried out. Additionally, we performed gene ontology (GO) and kyoto gene encyclopedia (KEGG) enrichment analyses to explore potential signal pathways. PCDHGA10 protein expression was evaluated using immunohistochemistry (IHC) on tissue microarrays (HLugS180Su02). By microarray analysis and database analysis, we found that PCDHGA10 was significantly highly expressed in LUSC. Sufferers with elevated PCDHGA10 levels exhibited a worse prognosis, according to the survival analysis. The PCDHGA10 mutation was also linked to LUSC patient prognosis. Besides, PCDHGA10 was closely related to tumor immune cell infiltration andimmunecheckpoints. In conclusion, PCDHGA10 is expected to become a new molecular marker for LUSC.

Multiplexed bulk and single-cell RNA-seq hybrid enables cost-efficient disease modeling with chimeric organoids

Disease modeling with isogenic Induced Pluripotent Stem Cell (iPSC)-differentiated organoids serves as a powerful technique for studying disease mechanisms. Multiplexed coculture is crucial to mitigate batch effects when studying the genetic effects of disease-causing variants in differentiated iPSCs or organoids, and demultiplexing at the single-cell level can be conveniently achieved by assessing natural genetic barcodes. Here, to enable cost-efficient time-series experimental designs via multiplexed bulk and single-cell RNA-seq of hybrids, we introduce a computational method in our Vireo Suite, Vireo-bulk, to effectively deconvolve pooled bulk RNA-seq data by genotype reference, and thereby quantify donor abundance over the course of differentiation and identify differentially expressed genes among donors. Furthermore, with multiplexed scRNA-seq and bulk RNA-seq, we demonstrate the usefulness and necessity of a pooled design to reveal donor iPSC line heterogeneity during macrophage cell differentiation and to model rare WT1 mutation-driven kidney disease with chimeric organoids. Our work provides an experimental and analytic pipeline for dissecting disease mechanisms with chimeric organoids.

Comparative genomic profiling of transport inhibitor Response1/Auxin signaling F-box (TIR1/AFB) genes in eight Pyrus genomes revealed the intraspecies diversity and stress responsiveness patterns

In the genomics of plants and the phytoecosystem, Pyrus (pear) is among the most nutritious fruits and contains fiber that has great health benefits to humans. It is mostly cultivated in temperate regions and is one of the most cultivated pome fruits globally. Pears are highly subjected to biotic and abiotic stresses that affect their yield. TIR1/AFB proteins act as auxin co-receptors during the signaling of nuclear auxins and play a primary role in development-related regulatory processes and responses to biotic and abiotic stresses. However, this gene family and its members have not been explored in Pyrus genomes, and understanding these genes will help obtain useful insights into stress tolerance and ultimately help maintain a high yield of pears. This study reports a pangenome-wide investigation of TIR1/AFB genes from eight Pyrus genomes: Cuiguan (Pyrus pyrifolia), Shanxi Duli (P. betulifolia), Zhongai 1 [(P. ussuriensis × communis) × spp.], Nijisseiki (P. pyrifolia), Yunhong No.1 (P. pyrifolia), d’Anjou (P. communis), Bartlett v2.0 (P. communis), and Dangshansuli v.1.1 (P. bretschneideri). These genes were randomly distributed on 17 chromosomes in each genome. Based on phylogenetics, the identified TIR1/AFB genes were divided into six groups. Their gene structure and motif pattern showed the intraspecific structural conservation as well as evolutionary patterns of Pyrus TIR1/AFBs. The expansion of this gene family in Pyrus is mainly caused by segmental duplication; however, a few genes showed tandem duplication. Moreover, positive and negative selection pressure equally directed the gene’s duplication process. The GO and PPI analysis showed that Pyrus TIR1/AFB genes are associated with abiotic stress- and development-related signaling pathways. The promoter regions of Pyrus TIR1/AFB genes were enriched in hormone-, light-, development-, and stress-related cis elements. Furthermore, publicly available RNA-seq data analysis showed that DaTIR1/AFBs have varied levels of expression in various tissues and developmental stages, fruit hardening disease conditions, and drought stress conditions. This indicated that DaTIR1/AFB genes might play critical roles in response to biotic and abiotic stresses. The DaTIR1/AFBs have similar protein structures, which show that they are involved in the same function. Hence, this study will broaden our knowledge of the TIR1/AFB gene family in Pyrus, elucidating their contribution to conferring resistance against various environmental stresses, and will also provide valuable insights for future researchers.

Comprehensive assessment of mRNA isoform detection methods for long-read sequencing data

The advancement of Long-Read Sequencing (LRS) techniques has significantly increased the length of sequencing to several kilobases, thereby facilitating the identification of alternative splicing events and isoform expressions. Recently, numerous computational tools for isoform detection using long-read sequencing data have been developed. Nevertheless, there remains a deficiency in comparative studies that systemically evaluate the performance of these tools, which are implemented with different algorithms, under various simulations that encompass potential influencing factors. In this study, we conducted a benchmark analysis of thirteen methods implemented in nine tools capable of identifying isoform structures from long-read RNA-seq data. We evaluated their performances using simulated data, which represented diverse sequencing platforms generated by an in-house simulator, RNA sequins (sequencing spike-ins) data, as well as experimental data. Our findings demonstrate IsoQuant as a highly effective tool for isoform detection with LRS, with Bambu and StringTie2 also exhibiting strong performance. These results offer valuable guidance for future research on alternative splicing analysis and the ongoing improvement of tools for isoform detection using LRS data.

Identification of Molecular Subtypes Associated with PCD in Esophageal Squamous Cell Carcinoma and Analysis of Immune Microenvironment.

Esophageal squamous cell carcinoma (ESCC) is a highly fatal malignancy with increasing incidence, and programmed cell death (PCD) plays an important role in homeostasis. This study aimed to explore the ESCC of heterogeneity based on the PCD signatures for the diagnosis and treatment of patients. The clinical information and RNA-seq data of patients with ESCC and the PCD-related genes set were used to identify PCD signatures.The "limma" package was used to identify the differentially expressed genes (DEGs). "Clusterprofiler" package was used for function enrichment analysis, and the "ConsensusClusterPlus" package was performed for consensus clustering. Finally, the "GSVA" package and the Cibersort algorithm were used for the immune infiltration analysis. We performed differential expression analysis between ESCC and normal samples and identified 1659 DEGs, of which 124 DEGs were PCD genes. Then, the patients were divided into cluster1 and cluster2 based on the expression of 124 PCD genes. There was a significant difference in immune infiltration between the two clusters. The patients in cluster 1 had a higher immune score and more CD56dim natural killer cells, monocytes, activated CD4 T cells, eosinophil, and activated B cells infiltration, while cluster2 had a higher stromal score, more immune regulation, and immune checkpoint genes expression. We identified two clusters based on PCD gene expression and characterized their tumor microenvironment and immune checkpoint difference. Our findings may provide some new insight into the treatment of ESCC.

Transcriptome dynamics of Gossypium purpurascens in response to abiotic stresses by Iso-seq and RNA-seq data

Gossypium purpurascens is a member of the Malvaceae family, holds immense economic significance as a fiber crop worldwide. Abiotic stresses harm cotton crops, reduce yields, and cause economic losses. Generating high-quality reference genomes and large-scale transcriptomic datasets across diverse conditions can offer valuable insights into identifying preferred agronomic traits for crop breeding. The present research used leaf tissues to conduct PacBio Iso-seq and RNA-seq analysis. We carried out an in-depth analysis of DEGs using both correlations with cluster analysis and principal component analysis. Additionally, the study also involved the identification of both lncRNAs and CDS. We have prepared RNA-seq libraries from 75 RNA samples to study the effects of drought, salinity, alkali, and saline-alkali stress, as well as control conditions. A total of 454.06 Gigabytes of transcriptome data were effectively validated through the identification of differentially expressed genes and KEGG and GO analysis. Overwhelmingly, gene expression profiles and full-length transcripts from cotton tissues will aid in understanding the genetic mechanism of abiotic stress tolerance in G. purpurascens.