Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “transcriptome annotation”

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.

At least 145 records · Page 8Linked to original sources

scGPA: an LLM-assisted workflow for directional virtual gene perturbation analysis from single-cell transcriptomes.

BACKGROUND: Existing virtual perturbation methods can often infer directional changes by comparing predicted post-perturbation expression profiles with control cells. However, workflows that directly return direction-specific downstream candidate genes together with confidence scores, evidence support and interpretable summaries remain limited. We developed scGPA, an LLM-assisted workflow system for directional single-cell virtual gene perturbation analysis. METHODS: scGPA starts from raw single-cell RNA sequencing data and performs quality control, normalization, dimensionality reduction, clustering and cell-group selection. It then constructs cell-group-specific wild-type regulatory networks using repeated subsampling, principal component regression (PCR)/Ridge-based network inference and CP tensor denoising. Based on these networks, scGPA simulates dose-aware virtual knockdown of the target gene and applies signed perturbation propagation to estimate the magnitude and direction of downstream transcriptional responses. LLM assistance is used for marker-based cell-type annotation, evidence-guided candidate prioritization and user-facing biological summarization. RESULTS: We benchmarked scGPA across five public Perturb-seq datasets and compared its performance with GEARS, scGPT and a random baseline. The overall correct prediction rate of scGPA was 23.0%, exceeding those of GEARS (20.7%), scGPT (15.1%) and the random baseline (13.6%). These results indicate that scGPA achieved a higher correct prediction rate than the two comparator models and the random baseline. We subsequently evaluated scGPA using a public osteosarcoma single-cell dataset and performed qRT-PCR validation in 143B osteosarcoma cells. Among genes with significant experimental changes, scGPA achieved a directional concordance of 76.9%. When all tested downstream genes were counted, 37.0% were directionally correct, 51.9% showed no significant change and 11.1% changed in the opposite direction. CONCLUSIONS: scGPA provides a practical workflow system for predicting and prioritizing direction-specific downstream transcriptional responses after target-gene perturbation. By integrating single-cell regulatory network inference, signed virtual perturbation and LLM-assisted interpretation, scGPA supports target-gene function inference and downstream mechanistic investigation from single-cell transcriptomic data.

Single-Cell Gene Expression Analysis↗

Chromosome-level genome assembly of the hemiparasitic Taxillus sutchuenensis (Loranthaceae).

Taxillus sutchuenensis, an ecologically and medicinally important hemiparasitic plant that parasitizes diverse woody hosts, was sequenced to generate a high-quality chromosome-level genome assembly. PacBio HiFi long reads, RNA-seq transcriptome data, and Hi-C data were used to assemble a 406.32 Mb genome anchored onto nine pseudo-chromosomes, with a scaffold N50 of 45.59 Mb. The assembly showed high completeness and accuracy, supported by BUSCO (93.6%) and Merqury QV (70.6) assessments. The LTR Assembly Index (LAI) of 13.98 indicated excellent continuity. A total of 21,795 protein-coding genes were predicted, with 94.46% functionally annotated. Repetitive sequences accounted for 50.05% of the genome, primarily LTR retrotransposons. This genome provides a valuable resource for investigating the evolution, functional genomics, and parasitic mechanisms of hemiparasitic plants.

Genome, Plant↗

Identifying JAK2 and ANXA5 as Key Genes Linking Obstructive Sleep Apnea and Oxidative Stress via Machine Learning and Multilayer Transcriptomic Integration With Functional Validation.

Obstructive sleep apnea (OSA) is a common and severe sleep disorder closely associated with oxidative stress (OS). This study aims to identify and validate potential OS-related genes associated with OSA through bioinformatics methods. We successfully identified OS-related differentially expressed genes (OS-DEGs) by combining the limma test, weighted correlation network analysis (WGCNA), and OS-related genes from the GeneCards database. Key genes and potential biological roles were further identified using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), enrichment analysis, protein-protein interaction (PPI) network analysis, Lasso regression analysis, random forest algorithm, and support vector machine recursive feature elimination (SVM-RFE) method. Evaluate and validate the accuracy of key genes through receiver operating characteristic (ROC) curve analysis. The human single-cell RNA sequencing (scRNA-seq) dataset is used for cell classification annotation, analysis of key gene single-cell expression profiles, and virtual gene knockout experiments based on the scTenifoldKnk algorithm. Integrating scRNA-seq sequencing, pseudotime trajectory inference, cell-cell communication analysis, and bulk immune infiltration deconvolution reveals monocyte subtype remodeling in OSA. Finally, the expression levels of key genes in clinical samples were validated using real-time quantitative PCR (RT-qPCR) and Western blotting. A total of 57 common DEGs, indicating significant enrichment in OS, inflammation, and tumor pathways, particularly prominent in the immunometabolism pathway. By integrating DEGs, WGCNA, PPI results, and machine learning methods, key genes Janus kinase 2 (JAK2) and ANXA5 were screened out. JAK2 was significantly upregulated under disease conditions, while ANXA5 was significantly downregulated. ROC curve exhibited high accuracy (area under the curve [AUC] > 0.85). Human scRNA-seq analysis revealed that key genes were predominantly highly expressed in monocytes. Virtual knockout experiments demonstrated that these key genes play a crucial role in regulating immune responses and inflammatory reactions. PPI networks and enrichment analysis verified that downstream genes S100P, ALOX5AP, PROK2, and PADI4 may collaboratively participate in immune response and inflammation regulation. Finally, clinical sample experiment further validated the results of bioinformatics analysis. This study provides new research insights for the diagnosis, mechanism research, and treatment development of OSA in the future by integrating multilayer transcriptomic and machine learning techniques.

Humans↗

Transcriptome analysis of the diseased intervertebral disc tissue in patients with spinal tuberculosis.

OBJECTIVE: To investigate the differential expression genes (DEGs) in spinal tuberculosis using transcriptomics, with the aim of identifying novel therapeutic targets and prognostic indicators for the clinical management of spinal tuberculosis. METHODS: Patients who visited the Department of Orthopedics at the Second Hospital, Lanzhou University from January 2021 to May 2023 were enrolled. Based on the inclusion and exclusion criteria, there were 5 patients in the test group and 5 patients in the control group. Total RNA was extracted and paired-end sequencing was conducted on the sequencing platform. After processing the sequencing data with clean reads and annotating the reference genome, FPKM normalization and differential expression analysis were performed. The DEGs and long non-coding RNAs (LncRNAs) were analyzed for Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment. The cis-regulation of differentially expressed mRNAs (DE mRNAs) by LncRNAs was predicted and analyzed to establish a co-expression network. RESULTS: This study identified 2366 DEGs, with 974 genes significantly upregulated and 1392 genes significantly downregulated. The upregulated genes are associated with cytokine-cytokine receptor interactions, tuberculosis, and TNF-α signaling pathways, primarily enriched in biological processes such as immunity and inflammation. The downregulated genes are related to muscle development, contraction, fungal defense response, and collagen metabolism processes. Analysis of LncRNAs from bone tuberculosis RNA-seq data detected a total of 3652 LncRNAs, with 356 significantly upregulated and 184 significantly downregulated. Further analysis identified 311 significantly different LncRNAs that could cis-regulate 777 target genes, enriched in pathways such as muscle contraction, inflammatory response, and immune response, closely related to bone tuberculosis. There are 51 genes enriched in the immune response pathway regulated by cis-acting LncRNAs. LncRNAs that regulate immune response-related genes, such as upregulated RP11-451G4.2, RP11-701P16.5, AC079767.4, AC017002.1, LINC01094, CTA-384D8.35, and AC092484.1, as well as downregulated RP11-2C24.7, may serve as potential prognostic and therapeutic targets. CONCLUSION: The DE mRNAs and LncRNAs in spinal tuberculosis are both associated with immune regulatory pathways. These pathways promote or inhibit the tuberculosis infection and development at the mechanistic level and play an important role in the process of tuberculosis transferring to bone tissue.

Humans↗

Insights into dill (Anethum graveolens) flavor formation via integrative analysis of chromosomal-scale genome, metabolome and transcriptome.

INTRODUCTION: Dill (Anethum graveolens) is a significant medicinal herb belonging to the Apiaceae family. Owing to its high levels of volatile organic compounds (VOCs), dill is commonly utilized for essential oil extraction and medicine purpose. However, the biosynthesis of the crucial VOC in dill remains obscure. OBJECTIVES: Identify the key VOCs related to the flavor formation in dill and dissect the regulatory mechanism of their synthesis. METHODS: The dill chromosomal-level genome was constructed by PacBio HiFi, Hi-C, and BGISEQ second generation sequencing and assembly. The VOCs in dill leaves were identified through GC-MS. The potential mechanism involved in regulating the VOC accumulation in dill flavor formation was analyzed by multi-omics analysis. RESULTS: A 1.17 Gb chromosome-scale genome of dill with a contig N50 of 10.78 Mb was constructed. A total of 46,538 genes were annotated across 11 assembled chromosomes. Comparative genomics analysis suggested that transposable element insertions, especially LTR-Gypsy, have contributed to the evolution and expansion of the dill genome. The flavor formation of dill was mainly attributed to terpenoids, especially α-phellandrene, β-ocimene, and o-cymene. The contribution of expansion and replication of terpenoid synthesis pathway genes, especially terpene synthase (TPS), to the abundant terpenoid production of dill was identified. Differential gene expression patterns observed at various developmental stages and tissues provided key candidate genes for the regulation of terpenoid synthesis, as well as transcription factors. The different accumulation of esters and aromatics also affected the flavor formation of dill. The key genes implicated in the synthesis of anethole, namely AIS and AMT were further identified. CONCLUSION: This study constructed the chromosome level genome and identified the main VOCs and related key genes in flavor formation of dill, shedding lights on our understanding of terpenoid biosynthesis but also offered guidance for future genetic research on molecular breeding in Anethum graveolens.

Transcriptome↗

ceRNA network of lncRNAs and mRNAs in OSF-to-OSCC progression: Diagnostic biomarkers and functional pathways.

BACKGROUND: Oral submucous fibrosis (OSF) is a chronic potentially malignant disorder that can progress to oral squamous cell carcinoma (OSCC). Although dysregulated non-coding RNAs have been implicated in oral carcinogenesis, the competing endogenous RNA (ceRNA)-mediated regulatory mechanisms underlying OSF-to-OSCC progression remain poorly understood. This study aimed to identify candidate regulatory molecules and construct a putative lncRNA-miRNA-mRNA network associated with malignant transformation. METHODS: Publicly available microarray datasets (GSE117973 and GSE125866) were analyzed to identify differentially expressed genes between OSF and OSCC. Differentially expressed transcripts were classified into mRNAs and lncRNAs based on public transcript annotations. Highly correlated lncRNA-mRNA pairs were identified using Pearson correlation analysis and integrated with multiMiR-supported miRNA-mRNA interactions obtained from public databases to construct a putative ceRNA regulatory network. Functional characterization focused on apoptosis, epithelial-mesenchymal transition (EMT), and immune checkpoint-related pathways. Receiver operating characteristic (ROC) analysis was performed to evaluate diagnostic performance, and selected biomarkers were externally validated using The Cancer Genome Atlas (TCGA) OSCC cohort. RESULTS: Integrated transcriptomic analysis identified several dysregulated mRNAs and lncRNAs associated with OSF-to-OSCC progression. Network analysis highlighted TBC1D3B, RREB1, TEAD3, SREBF1, TMEM41B, FOXK2, and KIAA1958 as prominent hub genes within the putative regulatory network. Functional analyses demonstrated significant associations with apoptosis-, EMT-, and immune checkpoint-related genes, suggesting potential involvement in multiple biological processes contributing to malignant transformation. Several hub genes exhibited strong diagnostic performance, with ROC analysis yielding AUC values ranging from 0.891 to 1.000, indicating excellent discrimination between OSF and OSCC samples. External validation using TCGA further supported the relevance of the identified biomarkers in OSCC. CONCLUSIONS: This study provides a comprehensive transcriptomic framework describing putative lncRNA-miRNA-mRNA regulatory interactions associated with OSF progression to OSCC. The identified hub genes and regulatory networks represent candidate biomarkers for early detection and provide a foundation for future mechanistic and experimental validation. As the proposed ceRNA interactions are computationally inferred, further biological validation is required before clinical application.

RNA, Long Noncoding↗

[DNA arrays: technological aspects and applications].

The Human Genome Project has allowed considerable progress in the construction of physical and genetic maps and the identification of genes involved in human sicknesses. The accelerated accumulation of biological information and knowledge is due in large part to the sequencing projects of other organisms, which in fact paved the way for the Human Genome Project. In parallel, recently developed techniques which take advantage of genomic sequences allow large scale molecular analyses resulting in the functional annotation of many of the proteins represented by these genes. This is the goal of functional genomics. These progresses are at the origin of the present revolution in biomedical research. DNA microarrays are playing a dominant role compared to the other developing technologies since they are relatively easy to make and use and are applicable to numerous scientific inquiries. They allow the simultaneous analysis of several thousands of genes in biological samples from sick or healthy tissues, at the genome or transcriptome level. The data obtained is expected to result in major advances in the health sciences. In addition to an improved understanding of the complex molecular interaction networks of healthy cells and tissues, a more precise genetic characterization of the molecular mechanisms involved in pathology should result in the identification of new therapeutic targets and the development of new medicines. The genetic profiles thus obtained should also permit the definition of new pathologic subclasses not recognizable by traditional clinical factors, as well as new markers for susceptibility to certain illnesses, and new prognostic markers or methods of predicting responses to treatment. In this article, we present the different approaches and potential applications of DNA microarray technology, in particular as applied to cancer research.

Chromosome Mapping↗

Identification of a putative RocS homolog through phenotypic profiling of uncharacterized essential genes in Streptococcus mutans.

Genome-wide viability catalogs produced by transposon sequencing (Tn-seq) and CRISPR interference (CRISPRi) have successfully mapped the essential genome of Streptococcus mutans . In this study, we combined predictive bioinformatics, conditional CRISPRi transcriptional silencing, transmission electron microscopy, transcriptomics, and genetic suppressor screens to investigate nine poorly characterized essential genes in S. mutans . From this screen, phenotypic and genetic analyses identified SMU_393 as a functional homolog of the pneumococcal chromosome segregation factor, RocS. Depletion of SMU_393 resulted in abnormal cell widening, hypersensitivity to DNA damage, and a significant subpopulation of anucleate cells. These phenotypes were bypassed by a spontaneous surface-exposed missense mutation ( dnaA Q197E ) within the AAA+ ATPase domain of the replication initiator. Together, this study refines annotations within the S. mutans essential genome and provides genetic insights into streptococcal chromosome segregation and cell cycle control.

Journal Article↗

Analysis of bovine mammary gland EST and functional annotation of the Bos taurus gene index.

Functional genomic studies of the mammary gland require an appropriate collection of cDNA sequences to assess gene expression patterns from the different developmental and operational states of underlying cell types. To better capture the range of gene expression, a normalized cDNA library was constructed from pooled bovine mammary tissues, and 23,202 expressed sequence tags (EST) were produced and deposited into GenBank. Assembly of these EST with sequences in the Bos taurus Gene Index (BtGI) helped to form 5751 of the current 23,883 tentative consensus (TC) sequences. The majority (87%) of these 5751 assemblies contained only one to three mammary-derived EST. In contrast, 18% of the mammary EST assembled with TC sequences corresponding to 12 genes. These results suggest library normalization was only partially effective, because the reduction in EST for genes abundantly transcribed during lactation could be attributed to pooling. For better assessment of novel content in the mammary library and to add to existing annotation of all bovine sequence elements, gene ontology assignments, and comparative sequence analyses against human genome sequence, human and rodent gene indices, and an index of orthologous alignments of genes across eukaryotes (TOGA) were performed, and results were added to existing BtGI annotation. Over 35,000 of the bovine elements significantly matched human genome sequence, and the positions of some alignments (3%) were unique relative to those using human expressed sequences. Because 3445 TC sequences had no significant match with any data set, mammary-derived cDNA clones representing 23 of these elements were analyzed further for expression and novelty. Only one clone met criteria suggesting the corresponding gene was a divergent ortholog or expressed sequence unique to cattle. These results demonstrate that bovine sequence expression data serve as a resource for characterizing mammalian transcriptomes and identifying those genes potentially unique to ruminants.

Animals↗

SpliceHarmonization: an integrated method for identifying RNA splicing events in therapeutics for splicing modulation.

MOTIVATION: Splicing, a critical co-transcriptional process in eukaryotes, enhances transcriptome diversity by generating isoforms specific to cell types, tissues, or developmental stages. Recent advancements in splicing modulators have opened new avenues for targeting previously undruggable genes by inducing significant perturbations in splicing events. These developments underscore the need for comprehensive methods to accurately identify and compare splicing events. While several tools have been developed to detect local splice variants, inconsistencies across methods remain a significant challenge. To address this, we present SpliceHarmonization, an integrated approach that combines the strengths of rMATS, LeafCutter, and MAJIQ, enabling robust and reliable splicing analysis with event type annotations. RESULTS: In a comprehensive evaluation using diverse simulated datasets, SpliceHarmonization streamlined and standardized the outputs from three detection methods into a unified format, thereby improving splicing detection with event type annotation and outperforming individual methods. By integrating the outputs from rMATS, LeafCutter, and MAJIQ, our approach not only enhanced identification of a wide range of splicing events but also effectively mitigated method-specific discrepancies. This integration led to an accuracy exceeding 0.8 and a recall of up to 0.5, with an observed increase in AUC of up to 10%. Furthermore, SpliceHarmonization demonstrated high sensitivity in detecting low-abundance and complex splicing events, providing annotations including genomic coordinates and event type. AVAILABILITY AND IMPLEMENTATION: SpliceHarmonization is available at https://github.com/interactivereport/SpliceHarmonization.

RNA Splicing↗

The UTRs of Leishmania donovani vary in length and are enriched in potential regulatory structures.

Leishmania spp. regulate gene expression largely post-transcriptionally, yet untranslated regions (UTRs) remain poorly delineated. We generated high-quality genome and transcriptome datasets for Leishmania donovani strain 1S2D (Ld1S) by combining PacBio HiFi de novo assembly with Oxford Nanopore direct RNA sequencing of promastigotes and axenic amastigotes. The genome assembly consists of 65 scaffolds totaling ~33.3 Mb. Structural comparisons to LdBPK282A1 revealed numerous rearrangements, including some reshuffling genes among polycistronic transcription units and validated by polycistronic reads from RNA sequencing. Promastigote and amastigote RNA sequencing produced 469,010 and 46,729 monocistronic reads containing a spliced-leader and a polyA tail sequences, defining 8,479 transcripts and supporting 7,415 of the 7,969 annotated protein coding genes, as well as 604 putative long non-coding RNAs. We annotated UTRs for 4,921 genes and observed that putative RNA G-quadruplexes were markedly enriched in UTRs. We also noted that 31.9% and 11.5% were expressed into multiple isoforms in promastigotes and amastigotes, respectively. Collectively, these data provide a comprehensive annotation of L. donovani genes and their UTRs and reveal widespread and stage-specific UTR length polymorphisms, and, overall, points to an important role of 3' UTR in post-transcriptional regulation in L. donovani.

Journal Article↗

High-Content CRISPR Screening: Methods and Applications.

Clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 screening has become a central technology in functional genomics, enabling genome-scale interrogation via pooled perturbations. Early CRISPR screens employed survival or simple phenotypic readouts to identify essential genes and drug resistance mechanisms. However, as biological questions have shifted toward understanding regulatory networks, cellular heterogeneity, and context-dependent gene functions, there has been increasing demand for screening strategies capable of capturing complex cellular phenotypes beyond cell fitness. Recent advances in single-cell sequencing, high-content imaging, and spatial transcriptomics have expanded the resolution of CRISPR screening by enabling multidimensional phenotypic characterization following genetic perturbation. By integrating pooled perturbations with diverse readouts, these approaches systematically map targeted gene edits to transcriptional states, cellular phenotypes, and microenvironmental contexts. Meanwhile, innovations in library design, delivery, and computational pipelines have further improved the robustness and interpretability of high-content screening platforms. This review synthesizes the methodological evolution of CRISPR screening, emphasizing advances in perturbation strategies, delivery systems, and multimodal readouts. Representative applications spanning oncology, immunotherapy, developmental biology, neurobiology, and infectious diseases are delineated to demonstrate refined gene network annotations. Additionally, existing technical bottlenecks, such as scalability, cost constraints, and in vivo limitations, are critically assessed. Finally, future directions are proposed to facilitate the development of precise medicine.

CRISPR screening↗

Structural variant discovery and diagnostic impact in rare diseases from short-read and long-read sequencing.

Rare diseases collectively affect 1 in 10 individuals, yet current genetic testing fails to identify a causal variant for most cases. At present, cytogenetic methods and/or sequencing approaches such as exome (ES) or short-read genome sequencing (srGS) represent the state-of-the-art for comprehensive clinical discovery of sequence and structural variants (SVs), including copy number variants, balanced SVs, complex SVs, and tandem repeats (TRs). Recently, long-read genome sequencing (lrGS), coupled with multiomics data, has presented great promise to resolve variation in genomic regions recalcitrant to characterization by srGS such as highly repetitive simple repeat sequences and segmental duplications. However, there are few guidelines to enable clinical interpretation of genetic variation in these highly repetitive genomic regions, and the enthusiasm of the field in adopting lrGS has made it difficult to assess the true added diagnostic yield of this technology due to widely variable and inconsistently applied analytic pipelines and variable degrees of pre-screening by ES or srGS. Here, we investigated the contribution of SVs to rare diseases using srGS as a front-line strategy when paired with highly sensitive SV discovery and evaluate the added diagnostic yield of incorporating lrGS for a subset of cases. Our srGS analysis encompassed 1,462 families (3,450 individuals) recruited through the Broad Institute Center for Mendelian Genetics and the Genomics Research to Elucidate the Genetics of Rare Diseases (GREGoR) programs. Diagnostic SVs were identified in 5.4% of cases (79/1,462), of which 80% were uniquely detectable by srGS compared to standard cytogenetic techniques. For 96 families (including 10 families with a heterozygous variant observed in a known recessive gene of clinical relevance), we performed lrGS with methylation profiling, as well as long-read transcriptomic analyses in a subset of 20 trios. Analyses with lrGS yielded over 25,000 SVs per genome, 63% of which were not captured by srGS, along with an additional ~200 rare SNV/indels per genome not previously captured and 12 differentially methylated regions per genome. Among these, we identified only one diagnostic variant not interpreted by srGS, an apparently mosaic de novo SNV in CASK that was absent in the srGS callset due to allelic imbalance. No new diagnoses were supported by long-read transcriptomics or episignatures. In this well characterized rare disease cohort, the added diagnostic yield was thus 1.04% (1/96 families). Following a systematic literature review of prior lrGS studies, we find that most reported diagnoses were detectable by srGS and that our added diagnostic yield is consistent with those prior studies. These studies emphasize the significant impact of comprehensive SV discovery in rare disease cases and further demonstrate the power for increased discovery of novel genomic variation and episignatures from lrGS. Nonetheless, they also serve to temper expectations of dramatic diagnostic advances in rare disease patients until there is more extensive annotation of the functional and clinical impact of all coding and noncoding variation uniquely accessible to lrGS with extensive reference databases spanning highly repetitive genomic sequencing that could be enabled by this transformative technology.

Journal Article↗

An Annotated Biobank of Triple-Negative Breast Cancer Patient-Derived Xenografts Features Treatment-Naïve and Longitudinal Samples during Neoadjuvant Chemotherapy.

UNLABELLED: Triple-negative breast cancer (TNBC) that fails to respond to neoadjuvant chemotherapy (NACT) can be lethal. Developing effective strategies to eradicate chemoresistant disease requires experimental models that recapitulate the heterogeneity characteristic of TNBC. To that end, we established a biobank of 92 orthotopic patient-derived xenograft (PDX) models of TNBC from the tumors of 75 patients enrolled in A Robust TNBC Evaluation fraMework to Improve Survival clinical trial (ARTEMIS, NCT02276443), including 12 longitudinal sets generated from serial patient biopsies collected throughout NACT treatment and from metastatic disease. Models were established from both chemosensitive and chemoresistant tumors, and nearly 30% of the PDX models were capable of metastasizing to the lungs. Comprehensive molecular profiling demonstrated conservation of genomes and transcriptomes between patient and corresponding PDX tumors, with representation of all major transcriptional subtypes. Transcriptional changes observed in the longitudinal PDX models highlighted dysregulation in pathways associated with DNA integrity, extracellular matrix interactions, the ubiquitin-proteasome system, epigenetics, and inflammatory signaling. These alterations revealed a complex network of adaptations associated with chemoresistance. Overall, this PDX biobank provides a valuable tool for tackling the most pressing issues facing the clinical management of TNBC. SIGNIFICANCE: The development of a patient-derived xenograft biobank that comprehensively captures the genomic and transcriptional diversity of triple-negative breast cancer promises to be a robust resource to investigate and overcome chemoresistance and metastasis.

Animals↗

Moderate expression and activity of flocculins underlie the characteristic flocculation phenotype of Saccharomyces pastorianus.

Flocculation is a key technological trait in lager brewing, governing fermentation performance, yeast recovery, and beer quality. In the allo-aneuploid hybrid yeast Saccharomyces pastorianus, the genetic basis of flocculation remains poorly resolved due to its complex dual sub-genome architecture. Here, we systematically re-annotated and functionally characterized the complete FLO gene repertoire of the Group II strain CBS 1483. Thirteen FLO genes were identified, including allelic variants and a previously uncharacterized adhesin, Flo12, containing a Hyphal_reg_CWP domain instead of the canonical PA14 lectin-binding domain. Structural modeling revealed strong conservation of Ca²+-binding residues in PA14 domains, alongside repeat-region diversification likely contributing to functional variability. Using optogenetic expression in a FLO-null background, we demonstrated that SpcI-FLO9-1 and SpcI-FLO9-2_1 are the strongest drivers of flocculation, exhibiting NewFlo-like sugar sensitivity. Transcriptomic analysis during 17°P wort fermentation showed dynamic induction of these genes coinciding with flocculation onset. Surprisingly, deletion of both loci in CBS 1483 did not abolish but only delayed sedimentation in wort, accompanied by improved maltose utilization and attenuation. These findings reveal functional redundancy and compensatory mechanisms within the FLO network of lager yeast, highlighting the genetic complexity underlying flocculation, and providing a molecular framework to inform yeast selection, strain development, and optimization of the lager fermentation processes.IMPORTANCEFlocculation, the process by which yeast cells aggregate and settle, is essential for producing clear, high-quality lager beer, and for efficient yeast recovery during brewing. However, the genetic basis of this trait in lager yeast has remained poorly understood because these strains possess unusually complex hybrid genomes. In this study, we systematically identified and characterized the complete set of flocculation genes in the industrial lager yeast Saccharomyces pastorianus CBS 1483. We demonstrated that lager yeast flocculation is not controlled by a single dominant gene, but instead emerges from the combined action of several moderately active adhesion proteins that are expressed at low levels during fermentation. Surprisingly, deleting the two strongest candidate genes only delayed, rather than eliminated, sedimentation, revealing a robust compensatory network that preserves brewing performance. These findings refine the current understanding of yeast flocculation and provide a molecular framework for developing brewing strains with improved fermentation efficiency, product consistency, and flavor quality.

Saccharomyces pastorianus↗

Enhanced chromatin compaction is associated with de novo expression of a nuclear microprotein, global loss of H3 acetylation and local transcriptional changes in retinal rod photoreceptors.

We have limited understanding of how aging alters gene expression and remodels cellular architecture in post-mitotic neurons. The inverted nuclear organization of mouse rod photoreceptors provides a unique model to gain mechanistic insights into age-associated decline in neuronal function. We have generated and integrated multi-omic datasets including 3D-genome topology, histone modifications, chromatin accessibility, DNA methylation and transcriptome of rod photoreceptors from young- and aged-mice. We show that aging drives global chromatin compaction, with regional alterations enriched at active chromatin. Epigenomic and transcriptional changes broadly correlate with chromatin dynamics as validated by high resolution microscopy. We uncover a megabase-sized genomic region with multi-level alterations, including de novo transcription of Gm7239, which encodes a functional microprotein carrying histone acetyltransferase-inhibitor domain. Overexpression of Gm7239 is associated with global loss of histone H3 acetylation, highlighting a potential new axis of genomic regulation in aging. Finally, we identify multiple significant local transcriptional alterations in non-annotated regions and genes associated with age-related macular degeneration. Our studies link age-related chromatin landscape changes with gene expression that may influence rod function and vulnerability to diseases.

Journal Article↗

A single-cell transcriptomic atlas of the pigtail macaque placenta in late gestation.

The placenta is a complex organ with multiple immune and non-immune cell types that promote fetal tolerance and facilitate the transfer of nutrients and oxygen. The nonhuman primate (NHP) is a key experimental model for studying human pregnancy complications, in part due to similarities in placental structure, which makes it essential to understand how single-cell populations compare across the human and NHP maternal-fetal interface. We constructed a single-cell RNA-Seq (scRNA-Seq) atlas of the placenta from the pigtail macaque ( Macaca nemestrina ) in the third trimester, comprising three different tissues at the maternal-fetal interface: the chorionic villi (placental disc), chorioamniotic membranes, and the maternal decidua. Each tissue was separately dissociated into single cells and processed through the 10X Genomics and Seurat pipeline, followed by aggregation, unsupervised clustering, and cluster annotation. Next, we determined the maternal-fetal origins of cell populations and analyzed single-cell RNA trajectory, Gene Ontology enrichment, and cell-cell communication. Single-cell populations in the pigtail macaque were strikingly similar in their identity and frequency to those found in the human placenta, including cells from trophoblast, stromal cell, immune, and macrophage lineages. An advantage of our approach was the deep sequencing of three tissues at the maternal-fetal interface, which yielded a rich diversity of common and rare single-cell populations. The third-trimester pigtail macaque single-cell atlas enables the identification of cellular subclusters analogous to those in humans and provides a powerful resource for understanding experimental perturbations on the NHP placenta.

Journal Article↗

Receptor-defined targeting of a genomically unique melanoma-enriched noncanonical antigen.

Effective T cell-based immunotherapies require functional receptors that can be engineered and redeployed to recognize tumor-restricted antigens. Noncanonical peptides arising from transcription outside annotated protein-coding regions expand the antigenic landscape of cancer; however, systematic strategies to biologically prioritize and functionally validate such targets remain underdeveloped. Here, we integrated de novo transcript analysis, exon-resolved quantification, RNA in situ hybridization, and immunopeptidomics to identify melanoma-associated noncanonical transcripts and advance candidates through receptor-level validation. Among three recurrent melanoma-associated transcripts, EVA003 emerged as a lead target based on its distinct repeat-enriched genomic architecture, consistent tumor-enriched exon-level expression across independent datasets, and a genomically unique immunogenic core sequence. We demonstrate endogenous presentation of EVA003-derived peptides on HLA-A*03:01 and detect specific reactivity in patient-derived tumor-infiltrating lymphocytes. Single-cell transcriptomic profiling identified a dominant peptide-reactive clonotype, enabling isolation of a naturally occurring T cell receptor. Transfer of this receptor into healthy donor T cells conferred antigen-dependent activation and cytotoxicity against both peptide-pulsed targets and melanoma cells expressing EVA003 endogenously. Together, these findings establish a biologically informed strategy for prioritizing noncanonical tumor antigens and demonstrate that genomically unique, tumor-enriched noncanonical peptides can be presented to molecularly defined receptors capable of mediating cancer cell killing. These findings support the integration of prioritized noncanonical antigens into engineered T cell therapeutic strategies.

Humans↗