Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “scRNA-seq”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

50 recordsLinked to original sources
PubMed2026

ScRNA-seq analysis reveals the effects of nitrite stress on the endocrine system of the eyestalk in Litopenaeus vannamei.

Nitrite is a harmful substance generated in Litopenaeus vannamei farming systems, largely originating from the inadequate breakdown of surplus feed and shrimp feces. Its accumulation in the water can affect the growth and physiological functions of shrimp, damage the immune system, and even cause mass mortality, thus becoming a key environmental factor restricting the green development of the industry. Under nitrite stress, the eyestalk, as an important neuroendocrine regulatory center in crustaceans, participates in the stress adaptation of the organism and exerts a protective effect by regulating energy metabolism and immune function. However, the molecular regulatory mechanism of the eyestalk in response to nitrite stress remains unclear. In this study, single-cell RNA sequencing (scRNA-seq) technology was used to analyze the heterogeneity of eyestalk cells in L. vannamei under nitrite stress. A total of 18, 394 high-quality cells were obtained, and six major cell subpopulations, including Neurosecretory cell, Motor neuron, Sensory neuron, Interneuron, Neurogliocyte, and Support cell, were identified. Differential expression analysis identified 839 differentially expressed genes, and different cell types showed distinct specific responses to nitrite stress. Functional enrichment analysis indicated that pathways such as glycolysis, oxidative phosphorylation, ribosome function, and endoplasmic reticulum protein processing were significantly activated, while signal transduction and DNA repair-related pathways were inhibited. Further analysis revealed that nitrite stress could induce mitochondrial function changes and trigger oxidative stress, thereby affecting the neuroendocrine system function of the eyestalk. This study provided insights into transcriptomic responses of the eyestalk to nitrite stress at the single-cell level, laying a theoretical foundation for the management of aquaculture environments.

Animals↗
PubMed2026

Integrated single-cell and bulk transcriptomic analysis identifies a novel senescent fibroblast subtype associated with poor prognosis in acral melanoma.

BACKGROUND: Acral melanoma (AM) exhibits significant intratumoral heterogeneity, but its tumor microenvironment (TME) and immune regulation remain unclear. This study aims to dissect TME heterogeneity and establish a prognostic model based on key cell subpopulations. METHODS: We collected AM single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data from the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA). Unsupervised clustering, CellChat, and Scissor analysis were performed to characterize cellular heterogeneity, cell-cell communication, and prognosis-related cell subpopulations. Kaplan-Meier analysis was used to assess the prognostic value of key genes, which were further validated by multiplex immunohistochemistry (mIHC). RESULTS: In AM, Mel_C2, C7, and C9 with high SEMA6A and KIT expression were strongly linked to poor prognosis. We further identified a senescent fibroblast subpopulation (sCAF_CDKN2A) characterized by high fibroblast senescence signature (FSS) scores. Integrating Scissor analysis of fibroblast subtypes with bulk prognostic data, we identified COL3A1, VCAN, and KIT as prognosis-associated genes upregulated in poor-outcome-related fibroblast subsets. Cell-cell communication analysis revealed that sCAF_CDKN2A engages in an immunosuppressive network, interacting with regulatory T cells (Tregs) via MIF signaling and receiving signals from exhausted CD8+ T cells through PPIA-BSG interactions. Using transcription factor expression patterns from these fibroblast subtypes, we constructed a prognostic model that effectively stratified patients into distinct risk groups with significant differences in overall survival (OS). mIHC confirmed significantly higher protein levels of SEMA6A and COL3A1 in tumor tissues compared to matched normal tissues. CONCLUSIONS: We established a novel prognostic model for AM and identified sCAF_CDKN2A as an immunosuppressive senescent fibroblast subpopulation driving poor prognosis.

Acral melanoma↗
PubMed2026

Deciphering CD8+ T cell exhaustion in human cancers through single-cell and spatial transcriptomics.

Exhausted CD8+ T cells (Tex) within the tumor microenvironment (TME) represents a critical barrier limiting anti-tumor immune responses. Tex cells are characterized by upregulated inhibitory immune checkpoint receptors, reduced cytotoxicity, and functional heterogeneity. Their genomic features and regulatory networks remain poorly defined, and only a minority of patients respond to immune checkpoint blockade (ICB) therapy. Single-cell RNA sequencing (scRNA-seq), through high-resolution transcriptomic profiling, has revealed diverse Tex subpopulations, identified subpopulation-specific marker genes and regulatory pathways. Spatial transcriptomics has further mapped the spatial distribution of Tex and their interaction networks with immune cells, tumor cells, and stromal cells, elucidating the impact of spatial heterogeneity on Tex functionality. Current studies indicate that the exhausted state of Tex is dynamic and modifiable, with functional differences among subpopulations closely associated with tumor progression and therapeutic response. However, the genomic characteristics, epigenetic regulation, and spatial interaction mechanisms of Tex require further exploration. This review summarizes recent advances in high-resolution omics technologies for precisely dissecting Tex heterogeneity, functional features, and interactions with other cells. It emphasizes the central value of optimizing Tex-targeted tumor immunotherapy strategies, providing theoretical foundations and directional guidance for developing more effective anti-tumor immunotherapies.

Humans↗
PubMed2026

Dynamic evolution of chaperone-mediated autophagy is associated with tumor microenvironment remodeling and prognostic stratification in lung adenocarcinoma: insights from single-cell transcriptomics, ensemble machine learning, and experimental validation.

BACKGROUND: Lung adenocarcinoma (LUAD) shows prognostic heterogeneity, and tumor-node-metastasis (TNM) staging is limited for individualized management. Chaperone-mediated autophagy (CMA) maintains proteostasis, but its role during adenocarcinoma in situ (AIS)-minimally invasive adenocarcinoma (MIA)-invasive adenocarcinoma (IAC) progression remains unclear. METHODS: Single-cell RNA sequencing (scRNA-seq) data from GSE189357 and bulk transcriptomes from The Cancer Genome Atlas (TCGA)-LUAD and Gene Expression Omnibus (GEO) cohorts were integrated. CMA activity, cell-cell communication, weighted gene co-expression network analysis (WGCNA), tumor-normal differential expression, machine-learning survival modeling, tumor microenvironment (TME) features, drug sensitivity, and EPC1 function were analyzed. RESULTS: CMA-high tumor epithelial cells increased from AIS (58.1%) to MIA (65.7%) but declined in IAC (44.4%; p < 0.001). CMA-low cells preferentially received fibroblast-derived extracellular matrix cues. A CMA-negatively correlated module identified 69 core genes. Random survival forest (RSF) performed best among 117 machine-learning combinations (mean concordance index > 0.873). High-risk patients had worse survival across cohorts, and the risk score was independently associated with overall survival (hazard ratio = 16.013, 95% confidence interval: 9.579-26.768, p < 0.001). High-risk tumors showed proliferative activation and M0 macrophage enrichment, whereas low-risk tumors showed stronger immune-related signaling. EPC1 overexpression suppressed malignant phenotypes in A549 cells. CONCLUSION: CMA dynamics are associated with stromal and immune remodeling during LUAD progression. A CMA-based model provides robust prognostic stratification and may offer a basis for future TME-guided studies.

Chaperone-mediated autophagy↗
PubMed2026

ChIP-seq profiling identifies diapause-regulated H3K27me3 targets in the fat body of Culex pipiens.

Culex pipiens, a principal vector of significant arboviruses, survives winter through diapause, a hormonally controlled inactive phase that enhances endurance under severe cold circumstances. Recent data suggests that epigenetic processes, namely histone post-translational modifications (hPTMs), play a crucial role in regulating seasonal dormancy. Prior studies from our laboratory indicated a decrease in the methylation of Histone 3 (H3K27me3) in diapausing fat body tissue, associated with elevated expression of the histone demethylase UTX. Nonetheless, the precise genomic areas impacted by these chromatin alterations remained unidentified. We used chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) to delineate the genome-wide distribution of H3K27me3 across fat body chromatin in diapausing (D) and non-diapausing (ND) female Cx. pipiens. Notably, the higher signal at transcription start sites (TSSs) reflects localized redistribution rather than a global decrease, as diapausing fat bodies retain less H3K27me3 overall but concentrate it at promoters. To investigate the functional significance of these chromatin alterations, we confirmed a number of target loci via ChIP-qPCR and assessed gene expression with qRT-PCR. We identified many critical genes that were markedly increased in diapausing mosquitoes, exhibiting an inverse relation to H3K27me3 enrichment. Our data demonstrates different H3K27me3 chromatin landscapes between diapausing and non-diapausing Cx. pipiens, corroborating a hypothesis of selective, locus-specific repression in the non-diapause state and its targeted removal during diapause to permit activation of dormancy-associated genes. These results suggest that chromatin remodeling is a core driver of the diapause switch.

Animals↗
PubMed2026

A contextual activity score (CAS) for inferring ADAR-associated transcriptional activity across RNA-seq, single-cell, and spatial transcriptomics.

BACKGROUND AND OBJECTIVE: Adenosine-to-inosine RNA editing, catalyzed by Adenosine Deaminases Acting on RNA (ADARs), is a widespread modification involved in neural function, immune regulation, and cancer. The Alu Editing Index (AEI) is the standard metric to estimate ADAR activity but requires raw sequencing reads and is poorly suited for single-cell and spatial transcriptomic data. This study aimed to develop an alternative framework for inferring ADAR-associated transcriptional activity from gene expression data across diverse transcriptomic technologies. METHODS: We developed the Contextual Activity Score (CAS), a framework based on transcriptional signatures from ADAR perturbation experiments. Context-specific signatures were generated for human neurons, mouse neurons, and cancer models to infer ADAR1 and ADAR2 activity. CAS was computed from normalized gene expression matrices using regulon-based enrichment analysis. Performance was evaluated by comparing with the Alu Editing Index across bulk RNA sequencing datasets, simulated sequencing depths, and library preparation protocols. RESULTS: CAS showed strong concordance with the Alu Editing Index across multiple datasets, while remaining robust to reduced sequencing depth and different library protocols. Unlike the Alu Editing Index, CAS can be applied to single-cell and spatial transcriptomic data and enables the independent assessment of ADAR2 activity. In cancer and neuronal contexts, CAS captured biologically meaningful variations in ADAR-associated transcriptional activity at sample, cell-type, and spatial levels. CONCLUSION: CAS provides a scalable approach applicable across multiple RNA-seq protocols for estimating ADAR-associated transcriptional activity using gene expression data. This method, implemented in an open-source R package for broad adoption, expands the ability to study ADAR-associated transcriptional activity across transcriptomic modalities where direct editing quantification is challenging, such as single-cell and spatial transcriptomics.

Adenosine Deaminase↗
PubMed2026

Functional role and regulatory network of miR-22-3p in chicken hepatic lipid metabolism.

Although microRNA-22-3p (miR-22-3p) is abundantly expressed in the avian liver, its epigenetic role in lipid homeostasis remains largely uncharacterized. To elucidate its in vivo function, 14-day-old female Qingyuan Partridge chickens were intravenously injected with lentiviral vectors to establish miR-22-3p overexpression and knockdown models. Phenotypic analysis demonstrated that miR-22-3p knockdown significantly elevated hepatic triglyceride (TG) levels (p&#xa0;<&#xa0;0.05) and drove marked steatosis, whereas its overexpression reduced TG content. Transcriptome sequencing (RNA-Seq) revealed profound metabolic remodeling, identifying 23 core lipid-associated genes (e.g., ELOVL6, FADS2, ACSBG2, and PTGIS) heavily enriched in steroid biosynthesis, fatty acid metabolism, and elongation pathways. In conclusion, miR-22-3p functions as a bidirectional epigenetic rheostat that negatively regulates hepatic lipid deposition by orchestrating a multilayered polygenic network, providing novel molecular targets for mitigating avian metabolic disorders and optimizing production traits in indigenous poultry breeds.

Animals↗
PubMed2026

Exercise-associated epigenetic remodeling and TCR repertoire dynamics in Lynch syndrome carriers.

Lynch syndrome (LS) carriers are at elevated cancer risk. Emerging evidence suggests that exercise may serve as a non-pharmacologic preventive strategy, yet the epigenetic and immunological mechanisms underlying its protective effects in this population remain unclear. Here, we perform integrative multi-omics profiling of DNA methylation, gene expression, and the T cell receptor (TCR) repertoire in LS carriers undergoing a 52-week aerobic cycling intervention. We identify compartment-specific DNA methylation changes, including innate immune activation in cfDNA and oncogenic pathway repression in tissue. Integrative transcriptomic analysis highlights ISL1 as a key exercise-repressed, epigenetically regulated gene, and identifies FLCN as a colorectal cancer (CRC)-associated methylation target. TCR analysis reveals an exercise-associated increase in systemic repertoire diversity and tissue-specific clonal convergence, thus suggesting antigen-driven recruitment. Collectively, these findings uncover epigenetic and immune remodeling as potential mechanisms of exercise-mediated protection in LS.

Lynch syndrome↗
PubMed2026

Transcriptomic analysis provides molecular insights into the innate immune defense of Mactra veneriformis against Vibrio alginolyticus infection.

Mactra veneriformis is an economically important bivalve mollusc in China, but its aquaculture is frequently threatened by Vibrio infections, particularly Vibrio alginolyticus. To investigate the molecular immune response of M. veneriformis to V. alginolyticus, we performed RNA-seq analysis of hepatopancreatic tissues collected at 48&#xa0;h post-infection, the peak mortality time point, with PBS-injected individuals used as controls. Infection with V. alginolyticus caused severe histopathological damage in the hepatopancreas and resulted in a cumulative mortality of 53.3% over 14 d, compared with 3.3% in the control group. Transcriptomic analysis identified 2623 differentially expressed genes (DEGs), including 1585 significantly up-regulated genes and 1038 down-regulated genes. KEGG enrichment analysis demonstrated that DEGs were significantly enriched in immune related and metabolism pathways, including the JAK-STAT signaling pathway, RIG-I-like receptor (RLR) signaling pathway, and cytochrome P450 (CYP450) signaling pathway. Collectively, these findings revealed candidate immune related genes (tlr3, tlr5, myd88, nfkb1, il-17d, and ifi44l), a putative TLR-MyD88-NF-&#x3ba;B signaling axis, and KEGG signaling pathways, including JAK-STAT, RLR and CYP450, that may be involved in the innate immune response of M. veneriformis to V. alginolyticus infection. These results provide a transcriptomic basis for understanding host-pathogen interactions in this species and highlight candidate genes and pathways for future functional validation and potential application in disease-resistance breeding.

Animals↗
PubMed2026

Transcriptomic insights into thermal stress reveal physiological trade-off between thermal stress adaptation and reproductive investment in Spodoptera litura.

Spodoptera litura, a highly polyphagous lepidopteran pest, poses a major threat to agricultural productivity due to its remarkable adaptability to diverse environmental conditions. Although heat stress is known to trigger transcriptional reprogramming in insects, the molecular mechanisms underlying thermal stress responses in S. litura remain poorly understood. In the present study, fourth-instar larvae were exposed to acute heat stress (44&#xa0;&#xb0;C) and compared with control conditions (27&#xa0;&#xb1;&#xa0;1&#xa0;&#xb0;C) to investigate heat-induced transcriptional alterations affecting physiology and reproduction. High-quality RNA-Seq data achieved more than 80% mapping efficiency, with a total of 15,782 transcripts were identified. Transcriptome analysis of S. litura larvae showed 323 differentially expressed genes (DEGs), of which 262 genes were significantly upregulated and 61 were downregulated in heat-stressed larvae compared to the control group. The DEGs were associated with stress response, reproduction, signalling, proteostasis, detoxification, oxidative stress, metabolism, development, and chromatin regulation. Heat shock proteins genes, including HSP70, HSP90, and HSP27, together with co-chaperones such as TRET-1, STIP1, and Starvin, were strongly upregulated, indicating enhanced cellular protection against protein damage and oxidative stress under heat stress. Conversely, key reproductive and cell cycle-related genes, including BARR, CAPD2, FEO, CDK2 and MORULA, were significantly downregulated, suggesting reproductive impairment and developmental arrest. RT-qPCR validation corroborated the RNA-Seq findings, demonstrating a heat-induced physiological trade-off that prioritizes survival over reproduction. Consistent with these molecular responses, heat-stressed insects exhibited marked reproductive impairment, including significant reductions in gonadosomatic index, eupyrene sperm bundle count, mating frequency, mating success, female calling behaviour, copulation duration, fecundity, and egg fertility. Collectively, these findings provide comprehensive insights into the molecular basis of thermal adaptation in S. litura and demonstrate that acute heat stress compromises reproductive fitness while activating conserved stress-response pathways that promote short-term survival.

Animals↗
PubMed2026

Single-cell RNA sequencing provides further insights into the immunostimulatory action of freeze-dried Lactiplantibacillus plantarum on Penaeus vannamei shrimp.

Immunostimulation through dietary interventions opened new avenues in developing disease control and prevention tools for shrimp aquaculture. We have previously shown that feeding with freeze-dried Lactiplantibacillus plantarum (LAB) increased disease resistance of Penaeus vannamei against both Vibrio parahaemolyticus and white spot syndrome virus (WSSV) based on bulk RNA sequencing of shrimp gills. This tissue participates in ion transport and serves as a first line of defense against environmental stressors and pathogenic infections. However, characterization of their cell composition and functions remains limited. Here, we implemented a single-cell RNA sequencing approach to further gather insights into how feeding with freeze-dried LAB modulates host immunity which may not be evident with bulk RNA sequencing approach. A total of five clusters with unique transcriptional signatures were identified, corresponding to pillar cells, septal cells, and sessile hemocytes. Pseudo-bulk analyses at global- and cluster-levels showed differential expression of genes related to host immunity and metabolism. We further revealed how overall transcriptomic changes are not exclusively caused by gene expression changes but may also be driven by cell population dynamics. This study highlighted how single-cell RNA sequencing approach may shed light on the mechanisms of action of immunostimulants which may be masked in bulk transcriptome analyses.

Animals↗
PubMed2026

Integrative analysis of transcriptome and DNA methylome dynamics during caudal fin regeneration in silver pomfret (Pampus argenteus).

Caudal fin regeneration in teleost fish is a complex, multi-stage process involving coordinated molecular and cellular changes. While the role of epigenetic regulation particularly DNA methylation has been studied in model freshwater species such as zebrafish, its contribution to regeneration in marine teleosts remains largely unexplored. In this study, we integrated transcriptomic and DNA methylomic data to characterize the temporal dynamics of gene expression and methylation during caudal fin regeneration in the silver pomfret (Pampus argenteus). Using RNA-sequencing and reduced representation bisulfite sequencing (RRBS) at three biologically critical time points 1, 3, and 7&#xa0;days post-amputation (dpa), we characterized the spatiotemporal molecular landscape of caudal fin regeneration. These time points capture the key transitional phases of wound healing and inflammation (1 dpa), blastema formation and progenitor proliferation (3 dpa), and regenerative outgrowth with tissue remodeling (7 dpa), enabling robust detection of the major molecular programs underlying epimorphic regeneration. Concurrently, CG-methylome analysis identified thousands of dynamically changing differentially methylated regions (DMRs). A strong global inverse correlation was observed between promoter methylation and gene expression. Integrative analysis pinpointed key regeneration genes (fgf20a, msxb, sox9b) whose expression was associated with dynamic methylation changes in their promoters or gene bodies. We conclude that DNA methylation is a dynamic and key regulatory layer that acts in concert with transcriptional reprogramming to coordinate tissue regeneration, providing new insights into the epigenetic mechanisms underlying complex regenerative processes in teleosts.

Animals↗
PubMed2026

Upscaling Genotyping by Amplicon Sequencing With GBAS-GUI.

Genotyping by amplicon sequencing (GBAS) is a relatively low-cost approach for generating genotypic data compared with established genomic methods, making it highly scalable and particularly suitable for large-scale genetic monitoring projects. However, most existing analytical pipelines are either marker-specific, insufficiently scalable, or lacking efficient data management systems for the long-term integration of genotypic information, limiting the full potential of GBAS. Here, we address this gap by introducing GBAS-GUI (https://github.com/sonnenbe-dot/GBAS-GUI), a pipeline capable of generating GBAS-based genotypic data for a wide variety of loci at scale. GBAS-GUI integrates a graphical user interface with multiple checkpoints to improve accessibility and robustness. It implements multiprocessing architecture and a relational database that links genotypic data with associated sample metadata to enhance scalability and data management. The pipeline further enables marker screening through automated calculation of polymorphism information content (PIC) and implements a strategy to recover homologous genotypic information from paralogous loci with non-overlapping amplicon length ranges. Using multiple empirical datasets, we demonstrate substantial improvements in processing speed, database management and handling artefacts related to co-amplification of unspecific regions and duplicates of the same genomic region. We further show that incorporating the full sequence information captured by an amplicon increases marker information content beyond what is achievable with length-based genotyping alone and expands the analytical versatility of GBAS. Overall, GBAS-GUI provides a robust, scalable and versatile framework that unlocks the potential of GBAS for large-scale population genetic and phylogeographic studies.

Genotyping Techniques↗
PubMed2026

Transcriptomic insights into temperature regulation of proliferation and differentiation in skeletal muscle cells of Nibea albiflora.

Myogenesis involves sequential stages of muscle satellite cell activation, myoblast proliferation, differentiation, and fusion into multinucleated myotubes. Teleost muscle exhibits indeterminate growth and is highly sensitive to environmental temperature, yet the underlying mechanisms by which temperature regulate proliferation and differentiation remain poorly understood. In this study, we established a primary skeletal muscle cell culture from the yellow drum (Nibea albiflora), an economically important marine fish, and integrated morphological observations with comparative transcriptomics analysis to characterize cellular and molecular responses at 28&#xa0;&#xb0;C and 20&#xa0;&#xb0;C during both proliferation and differentiation stages. Phenotypic analysis revealed that 28&#xa0;&#xb0;C significantly enhanced both myoblast proliferation and myogenic differentiation ability compared with 20&#xa0;&#xb0;C. Transcriptomic profiling revealed that at 28&#xa0;&#xb0;C, differentiation upregulated extracellular matrix(ECM) organization, calcium signaling, and sarcomere assembly, while proliferation enhanced focal adhesion, growth factor signaling, and lipid metabolism. At 20&#xa0;&#xb0;C, differentiation was characterized by glutathione metabolism, and ferroptosis, while proliferation involved cytokine-cytokine receptor interaction and negative regulation of signal transduction. Core myogenic regulatory factors (MRFs), particularly myogenin, were strongly upregulated at 28&#xa0;&#xb0;C during the differentiation stage, serving as an internal control. Based on these findings, we propose a testable model that elevated temperature coordinates Ca2+-dependent MRF activation with ECM-integrin signaling to drive sarcomere assembly and muscle growth. Key differentially expressed genes (DEGs) regulating myogenesis in N. albiflora were also identified. This study provides a mechanistic framework for temperature adaptation in teleost skeletal muscle and identifies candidate genes for functional validation and marker-assisted selection, as well as a rationale for temperature management strategies to improve aquaculture yield of N. albiflora.

Animals↗
PubMed2026

In vivo genome-wide CRISPR screens identify FOXR1 as a suppressor of CD8+ T cell antitumor immunity.

T cell dysfunction critically limits the efficacy of T cell-based immunotherapies in solid tumors, yet the intrinsic regulators of T cell dysfunction remain incompletely understood. Through an in vivo genome-wide CRISPR screen in tumor-infiltrating CD8+ T cells, we identified Forkhead Box R1 (FOXR1) as a potent transcriptional suppressor of CD8+ T cell effector functions. Genetic ablation of FOXR1 significantly enhanced cytokine production and cytotoxic capacity in both murine and human CD8+ T cells, whereas its overexpression impaired T cell activation and effector molecule expression. Mechanistically, multiomics integration of RNA-seq, CUT&Tag-seq, and ATAC-seq revealed that FOXR1 binds directly to promoter regions of key effector genes, including IL2, GZMB, and PRF1, and represses their expression. Importantly, FOXR1 deletion in human anti-CD19 CAR T cells improved their efficacy against solid tumors, demonstrating that FOXR1 is a checkpoint of T cell effector function and targeting FOXR1 is a promising strategy to enhance CAR T cell efficacy against solid tumors.

Animals↗
PubMed2026

Comparison of paralog identification methods and their impact on species tree topologies in target capture phylogenomics within the Sindora clade (Detarioideae: Leguminosae).

Target capture is a common method of generating high throughput DNA sequencing data for phylogenetic reconstruction of species relationships, for which single copy genes are usually most informative. However, a pervasive problem with target capture is that putatively single copy genes may in fact be paralogs resulting from gene duplication, which are problematic for phylogenetic inference because their evolutionary history may differ from the divergence history of species. Here, we use as a case study a target enrichment dataset of 88 species of Detarioideae (Leguminosae) with a focus on the Sindora clade to examine approaches for handling paralogs, including the built-in paralog handling functions in HybPiper and CAPTUS, plus subsequent steps using Putative Paralog Detection and the tree-based Yang & Smith orthology inference approach. We compare the paralogs flagged using these methods and verify their performance with BLAST mapping against a reference genome sequence of Sindora glabra, and then subsequently compare the species tree topologies produced across these methods. Our comparisons of paralogs flagged across the Sindora clade show that the Putative Paralog Detection pipeline was the most accurate in identifying paralogs in terms of its similarity to the BLAST mapping, followed by the built-in paralog identification function of CAPTUS. However, the results we recovered for the Detarioideae subfamily suggest that the largest differences in species tree topology resulted from the use of paralog-filtered alignments (such as with the Putative Paralog Detection pipeline and the Yang & Smith orthology inference approaches) rather than just by removing the sequences of identified paralogous genes. This was the true for HybPiper-assembled datasets but was not seen in CAPTUS-assembled datasets. In all comparisons, the topological differences caused by different paralog handling methods tended to be confined to clades where processes such as hybridisation and introgression are prevalent. Our study provides a roadmap to establish the best approach to identify, eliminate or separate paralogs in the absence of a chromosomally contiguous reference genome for a study group, and highlights the importance of careful data inspection and processing in addition to understanding the extent of paralogy and paralog characteristics (e.g. sequence divergence between copies) for their study group.

Phylogeny↗
PubMed2026

Translational reprogramming of TGF-&#x3b2; signaling via TRMT61A-mediated tRNA m1A drives prostatic fibrosis and hyperplasia.

Dysregulation of the epitranscriptomic landscape is closely linked to pathological proliferation, but its specific role in benign prostatic hyperplasia (BPH) remains unclear. Here, we identify the tRNA methyltransferase TRMT61A as a critical driver of BPH progression. We found that TRMT61A and global N1-methyladenosine (m1A) levels are aberrantly upregulated in human BPH tissues. Functionally, TRMT61A knockdown potently suppresses prostate cell proliferation and reduces stromal fibrosis, inducing G1 cell cycle arrest and reversing pathological remodeling both in vitro and in vivo. By integrating ribosome profiling (Ribo-seq) and tRNA-seq, we observed that TRMT61A drives translational reprogramming. TRMT61A preserves the stability of specific tRNA isoacceptors (e.g., tRNA-Leu-CAA), which is required for the efficient decoding of mRNAs containing m1A-dependent codons. Consequently, TRMT61A selectively promotes the translational elongation of the key receptor TGF&#x3b2;R1. This amplifies downstream TGF-&#x3b2;/SMAD signaling and drives epithelial-mesenchymal transition (EMT) without affecting mRNA transcription. In summary, our study reveals how TRMT61A drives BPH progression through TGF&#x3b2;R1 translation, highlighting the therapeutic potential of targeting epitranscriptomic pathways to reverse prostatic hyperplasia and fibrosis.

Male↗