Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “transcriptome”

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 199 records · Page 11Linked to original sources

FM-GPT: Bayesian fine mapping for phenome-wide transcriptome-wide association studies.

Transcriptome-wide association studies (TWAS) integrate genome wide association studies with expression quantitative trait locus reference panels to identify genes associated with traits of interest. However, linkage disequilibrium and correlated gene expression can induce spurious TWAS signals, motivating fine mapping methods to prioritize putatively causal genes within associated loci. The rapid growth of large-scale phenomic resources (e.g. electronic health records (EHRs)) has shifted genetic studies from single-trait analyses to phenome-wide investigations that jointly evaluate many closely related phenotypes. We introduce FM-GPT (Fine-mapping of causal Genes for Phenome-wide Transcriptome-wide association studies), a novel Bayesian fine mapping method for prioritizing causal genes across multiple correlated phenotypes with potentially mixed outcome types (e.g., binary, count or continuous) in phenome-wide TWAS. FM-GPT performs gene-guided dimension reduction of the phenotypes and reveals pleiotropic or phenotype-specific effects of the identified genes. In simulations, FM-GPT identified true causal genes more accurately than other fine mapping methods while controlling false positives. We applied FM-GPT to two applications using data from UK Biobank: a brain-wide genetic analysis of MRI data derived regional cortical thickness measures and a phenome-wide genetic analysis of clinical phenotypes derived from EHR data. FM-GPT greatly narrowed down the set size of putatively causal genes and identified: 1. genes with pleiotropic effects on regional cortical thickness across the cerebral cortex, including five genes BCAS3, LRRC37A, NOS2P3, ARL17B and UBB on chromosome 17 regulating neuronal morphology and cortical organization; and 2. genes that influence multiple medical conditions across the circulatory, metabolic, digestive, respiratory and genitourinary systems, revealing two major axes of variation among these conditions that point to a potential trade-off in gene regulation between immune and metabolic functions. These results highlight FM-GPT's power to disentangle complex gene-phenotype relationships in large-scale phenome-wide studies, uncovering shared biological mechanisms across diverse human traits and advancing translational and comorbidity research.

Bayesian fine mapping↗

Retinal Transcriptome-Wide Association Study Identifies Novel Alzheimer's Disease Risk Genes.

INTRODUCTION: Alzheimer's disease (AD) is the leading cause of dementia worldwide. The retina shares molecular pathways with the brain, yet no study has systematically linked retinal gene expression to AD risk. METHODS: We performed transcriptome-wide association studies (TWAS) using two independent retinal eQTL panels (Strunz et al., n = 311; EyeGEx, n = 406) and a large meta-analyzed AD genome-wide association study (GWAS) (Bellenguez et al., 111,326 cases, 677,663 controls). Genes were further validated with GWAS in the independent Alzheimer's Disease Sequencing Project (ADSP) using a matched eQTL-panel strategy. RESULTS: We identified 62 AD-associated genes across the two eQTL panels using Bellenguez et al. as the discovery cohort. Of these, 31 were replicated in the ADSP cohort. The findings highlight shared complement-mediated immune dysregulation (CD55, CD46, TREM2) and provide functional transcriptomic evidence to prioritize novel causal drivers of AD pathogenesis, including STYX and the LRRC37 gene family. DISCUSSION: Retinal data capture core AD genetic architecture and reveal novel risk genes, highlighting the retina as a molecularly informative tissue for dementia research.

Alzheimer’s disease↗

MKMC enables reference-free transcriptomic analysis using k-mer representations.

Traditional RNA-seq analysis depends heavily on genome alignment and gene annotation, limiting its utility in non-model organisms and introducing biases that can obscure regulatory complexity. We present MKMC (Multi-sample Kmer Counter), a scalable, reference-free toolkit for RNA-seq analysis that leverages k-mer-based statistics to detect biological variation without requiring alignment. MKMC integrates fast k-mer counting, abundance matrix generation, normalization, dimensionality reduction, and differential analysis into a unified workflow. Across diverse datasets, MKMC recapitulates key biological signals-including sex differences in killifish liver-and matches alignment-based pipelines in differential expression analysis and transcriptomic age prediction. Notably, MKMC detects isoform-specific events missed by traditional methods, one of which we validated using in situ hybridization. These results reveal previously hidden isoform-level regulatory events that contribute to sex- and age-associated transcriptional programs. MKMC offers a robust, extensible alternative to alignment-based approaches, enabling transcriptomic discovery across both model and non-model systems. While we focus here on RNA-seq as a primary application, MKMC is broadly applicable to any k-mer-based analysis of next-generation sequencing data.

MKMC↗

A single-nucleus and spatial transcriptomic atlas of poplar leaves reveals the regulation of leaf polarity and cuticle deposition.

Leaf adaxial-abaxial polarity is fundamental for plant morphogenesis and environmental adaptation through asymmetric cell differentiation. Emerging evidence reveals dorsoventral metabolic gradients act downstream of transcriptional networks to fine-tune cellular specialization. While conserved transcription factors (e.g., HD-ZIP III and KANADI) establish initial polarity, the molecular networks driving position-specific cellular differentiation and their integration with metabolic adaptation remain unclear. Leveraging single-nucleus and spatial transcriptomics, we resolve major cell classes (mesophyll, epidermal, and vascular-associated) and their adaxial-abaxial subtypes, revealing dorsoventral polarity in transcriptional profiles and metabolic pathways. Adaxial cells are enriched in phenylpropanoid/flavonoid biosynthesis, while abaxial cells show preferential activation of stress and hormone signaling. Notably, we identify MYC2 as a key regulator of adaxial cuticle biosynthesis, binding to promoters of lipid biosynthetic and transport genes (e.g., CER10 and LTPG1) and promoting cuticle thickening. Our study uncovers how positional identity shapes transcriptional and metabolic polarity in leaves, with MYC2 emerging as a central regulator coordinating organ-specific adaptations. These findings provide insights into the spatial regulation of plant development and stress resilience, offering potential strategies for engineering stress-tolerant woody crops.

Plant Leaves↗

Single-cell and spatial transcriptomics define a progenitor subpopulation and fibroinflammatory niche at the leading edge of parathyroid carcinoma.

Parathyroid carcinoma (PC) is a rare but clinically aggressive endocrine malignancy with limited treatment options and a poorly defined tumor microenvironment (TME). To elucidate its cellular heterogeneity and spatial architecture, we integrated single-cell and spatial transcriptomic profiling with whole-exome sequencing and multiplex immunohistochemistry on eight parathyroid neoplasm specimens, including PC, parathyroid adenoma, and atypical parathyroid tumor. We identified a distinct progenitor-like endocrine subpopulation (Ca-1) enriched in CDC73-mutant PC, exhibiting stem-like properties, elevated cell cycle activity, and pronounced genomic instability. Spatial mapping revealed that Ca-1 cells preferentially localize at the leading edge, forming a fibroinflammatory niche characterized by the enrichment of inflammatory cancer-associated fibroblasts (iCAFs) and SPP1+ macrophages. Within this niche, the dipeptidyl peptidase 4 (DPP4) is selectively expressed in Ca-1 cells and iCAFs, implicating a potential paracrine axis driving stromal remodeling and immunosuppression. These findings suggest that a spatially organized ecosystem may promote PC progression through TME remodeling and highlight the DPP4-CXCL2 axis as a candidate pathway for future investigation in aggressive parathyroid neoplasms.

Humans↗

Zea mays Drought-Overly Sensitive1/TUBA4 Is Wilty3, and Transcriptome Co-Expression Analysis of Shoot Meristem Mutant Tissues Reveals Wilty2/TUB6:Wi3 Interactions Associated With Stem Vascular Bundle Development.

Plant vasculature is essential for the transport of water, nutrients, and signaling molecules across organs, while also providing critical mechanical support for growth and development. Disruptions in vascular bundle formation can therefore lead to severe physiological and developmental defects. In maize, ethyl methanesulfonate (EMS)-induced dominant nonallelic Wilty mutants exhibit a pronounced wilting phenotype even under well-watered conditions, indicating underlying defects in vascular function. In this study, we characterized the Wi3 mutant, identified as ZmDrought-Overly-Sensitive1/DOS1, and compared it with the previously described Wi2 mutant to uncover shared mechanisms underlying their phenotypes. We provide evidence, by bulk segregant resequencing linkage disequilibrium of SNPs adjacent to the causal Wilty SNPs in respective ß- and α-tubulin genes, for the personal communication from Gerry Neuffer that Wi2/ß-tub6 provenance is from ACR-related stock, whereas Wi3/α-tub4 allele is from Mo17, not B73 as claimed by the authors who cloned Dos1. Histochemical staining and Fourier-transform infrared (FTIR) spectroscopy of vascular bundles in Wi3 indicated apparent alterations in cellulose and lignin content consistent with those observed in Wi2. Transcriptome analysis of shoot meristems further indicated that similar sets of genes and pathways are differentially expressed in both mutants, suggesting convergence on common biological pathways. Using bulk-segregant whole-genome resequencing, we identified alpha-tubulin4 (TUA4) as the causal gene in Wi3 (ZmDOS1), harboring a C-to-T substitution within the N-terminal GTPase-binding domain. This mutation results in a glutamic acid196-to-lysine substitution. Given that α- and β-tubulin subunits heterodimerize, and in many plants and animal mutant alleles are dominant-negative gains-of-function, we infer Wi2, Wi3, and likely Wi4, based on very similar FTIR biophysical difference spectra, may act as effectors of vascular bundle cell wall deposition, potentially involving vesicle trafficking as recently shown for asymmetric cell divisions in maize stomatal development. Together, these findings highlight the functional interdependence of tubulin subunits and provide a plausible mechanistic framework for the striking biophysical, transcriptomic, and phenotypic similarities observed between Wi2, Wi3/ZmDOS1, and Wi4 mutants.

bulk segregant analysis↗

Transcriptomic landscape of NK cell-related genes in hepatocellular carcinoma: associations with prognosis and therapeutic response.

Hepatocellular carcinoma (HCC) is a highly aggressive and heterogeneous malignancy, in which natural killer (NK) cells play a crucial role in tumor progression and immune surveillance. This study aimed to characterize the transcriptomic landscape of NK cell-associated genes (NAGs) and explore their associations with clinical outcomes and therapeutic responses in HCC. Using transcriptomic data from The Cancer Genome Atlas (TCGA), we identified key NAGs through comprehensive statistical analyses. Patients were stratified into distinct risk groups based on NAG expression profiles. Low-risk patients demonstrated better survival, higher immune infiltration, and greater predicted sensitivity to immunotherapy, whereas high-risk patients were associated with reduced chemotherapy responsiveness. These findings contribute to a deeper understanding of the immunogenomic features of HCC and provide a basis for developing personalized therapeutic approaches centered on NK cell-related mechanisms.

Hepatocellular carcinoma↗

Integrated transcriptomic analysis of mRNA and miRNA in Brown adipose tissue of the greater horseshoe bats during hibernation.

Hibernation enables animals survive harsh environments by conserving energy through reduced metabolism and body temperature. Brown adipose tissue (BAT) plays a critical role in non-shivering thermogenesis, crucial for warming up during arousal phase. The greater horseshoe bats (Rhinolophus nippon) are typical hibernators and non-shivering thermogenesis in BAT tissue may persist throughout the arousal process in bats. This study examines gene expression and regulatory changes in BAT of these bats across active, hibernation, and arousal phases using transcriptome and miRNA sequencing. A total of 2721 differentially expressed mRNAs and 268 differentially expressed miRNAs were identified. The results reveal that the BAT transcriptome undergoes state-dependent remodeling throughout the hibernation process. The most pronounced divergence occurs between the active phase and torpor, involving cell cycle arrest, immunosuppression, thermogenic signal desensitization, and upregulation of lipid metabolism and autophagy pathways, reflecting the coordinated adaptation of energy conservation and thermogenic reserve. In contrast, transcriptional alterations between torpor and arousal are extremely limited, indicating that torpid BAT is already pre-primed for thermogenesis and requires only modest transcriptional adjustments to activate heat production. Notably, although body temperature recovers to active-phase levels during arousal, the molecular signature of BAT remains highly similar to that of the torpid state. Furthermore, the core thermogenic gene UCP1 showed no significant expression differences across the three groups. In conclusion, this study systematically delineates the miRNA-mRNA regulatory landscape of bat BAT across the hibernation process, and deepens our understanding of the thermoregulatory mechanisms underlying mammalian hibernation.

BAT↗

Linking MRI radiomics to transcriptomics-based radiosensitivity in lower-grade glioma: A radiogenomic framework.

BACKGROUND: RSI is a transcriptomics-based biomarker associated with radiotherapy outcomes, but its clinical application is constrained by the requirement for tumor tissue and RNA sequencing. This study investigates whether MRI-derived radiomic features can reflect RSI-defined intrinsic radiosensitivity in lower-grade glioma.This addresses a critical gap arising from the limited availability of matched imaging and genomic data in routine clinical practice. METHODS: MRI-derived radiomic features were extracted from FLAIR images of lower-grade glioma patients obtained from TCIA and matched with transcriptomic data from TCGA. A total of 107 patients with both MRI and RNA sequencing data were included in the radiogenomic analysis. Radiomic features were ranked using a Borda-based ensemble feature selection strategy. Five supervised machine-learning classifiers were trained to predict RSI-based radiosensitivity classification, and model interpretability was assessed using SHAP within radiogenomic framework. RESULTS: Classification performance increased with feature number and stabilized at compact subset of 13 radiomic features. Logistic regression showed stable performance with an AUC of 0.82 (95 % CI: 0.71-0.93). SHAP analysis indicated that heterogeneity-related texture features were dominant contributors to model predictions, with many associated with the RR phenotype, while others were linked to the RS phenotype. CONCLUSION: An MRI-based radiomic signature enables non-invasive prediction of RSI-defined radiosensitivity in lower-grade glioma. Rather than offering an immediately deployable clinical tool, this study establishes a proof-of-concept radiogenomic framework demonstrating that intrinsic radiosensitivity, traditionally assessed through invasive molecular assays, can be approximated using quantitative imaging features. These findings highlight the potential of imaging-based radiosensitivity assessment and provide a foundation for future radiogenomic investigations.

Lower-grade glioma↗

Multi-omics and spatial transcriptomics reveal that S100A10 drives CD8+ T-cell exhaustion and immune evasion in hepatocellular carcinoma through cPLA2-5-LOX-mediated arachidonic acid metabolism and ferroptosis.

Immune evasion in hepatocellular carcinoma (HCC) represents a major biological barrier limiting the efficacy of immunotherapy, yet its molecular basis remains incompletely understood. Increasing evidence indicates that tumor metabolic reprogramming and ferroptosis-related signaling play critical roles in shaping an immunosuppressive tumor microenvironment (TME); however, the specific regulatory factors involved remain unclear. This study aims to systematically elucidate the functional role of S100 calcium-binding protein A10 (S100A10) in immune evasion in HCC, with a particular focus on the molecular mechanisms by which S100A10 regulates CD8+ T-cell exhaustion through arachidonic acid (AA) metabolism and ferroptosis, as well as its potential therapeutic implications. To this end, data from The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) cohort are integrated to analyze the expression patterns of S100A10, its prognostic value, and its association with the immune microenvironment. S100A10 overexpression and knockout models are established in HCCLM3 and MHCC97L cell lines, and S100A10-mediated metabolic pathway reprogramming is characterized using transcriptomic profiling, untargeted metabolomics, and ferroptosis-related functional assays. In parallel, single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics are employed to delineate the cell-type specificity and spatial distribution of S100A10. Furthermore, human CD8+ T-cell co-culture systems and orthotopic mouse HCC models are used to evaluate the impact of S100A10 on immune function and responsiveness to anti-programmed cell death protein 1 (anti-PD-1) therapy. The results demonstrate that S100A10 is significantly upregulated in HCC and is closely associated with poor prognosis and an immunosuppressive state. Mechanistically, S100A10 activates cytosolic phospholipase A2-arachidonate 5-lipoxygenase (cPLA2-5-LOX)-mediated AA oxidative metabolism, leading to the accumulation of lipid peroxidation products and ferroptosis-associated signals, thereby driving CD8+ T-cell exhaustion and promoting immune evasion. Significantly, inhibition of S100A10 reshapes the tumor immune microenvironment (TIME) and enhances the therapeutic efficacy of anti-PD-1 treatment. Collectively, these findings identify S100A10 as a critical regulator of metabolic-immune coupling in HCC and provide a theoretical basis for combinatorial strategies targeting metabolism and immunotherapy.

Arachidonic acid metabolism↗

Single-nucleus transcriptomics reveals cell type-specific remodeling and epilepsy-associated microglia.

Temporal lobe epilepsy (TLE) is the most common acquired epilepsy, causing refractory seizures and cognitive deficits. We performed single-nucleus RNA sequencing on hippocampal tissue from mice 3 and 6 weeks following pilocarpine-induced status epilepticus, a robust model of TLE. Epilepsy samples showed reductions in Cck and Lamp5-Lhx6 interneuron subclusters, alongside increases in Cajal-Retzius cells, dentate granule (DG) cell precursors, and a mature DG cell subcluster. Among glia, an astrocyte subcluster and a markedly expanded microglia sublcuster were increased. We term this microglia population epilepsy-associated microglia (EAM). The transcriptomic profile of EAM overlaps with microglia described in models of Alzheimer's disease and traumatic brain injury, including enrichment of Myo1e and Igf1. EAM display amoeboid morphology, can be found in clumps around pyramidal and granule cell body layers, and exhibit enlarged vesicles and mitochondria. Cell-cell interaction analysis predicts DG cells as their primary interaction partners. This dataset defines transcriptomic programs underlying key cellular alterations in TLE, enabling mechanistic dissection of epileptogenesis.

TLE↗

Genomic and transcriptomic features of relapsed small cell lung cancer.

BACKGROUND: Relapsed small cell lung cancer is characterized by treatment resistance and poor outcomes. Genomic and transcriptomic alterations in relapsed SCLC have not been characterized well. We comprehensively profiled relapsed SCLC samples along with patient-matched treatment-naive samples, when available, using whole-exome (WES), whole-genome (WGS), and RNA-sequencing (RNA-seq) to describe the molecular landscape of relapsed SCLC. Our goal is to identify potential novel pathways for additional functional validation and eventually novel therapeutic options. METHODS: We analyzed 54 relapsed and 27 treatment-naive SCLC samples using WES (with 26 patient-matched paired samples). A subset of the samples was also analyzed by WGS (n=28) and RNA-seq (n=31). Differences in mutational signatures, gene expression, structural variants, splicing, and neoantigen profiles at diagnosis and relapse were investigated. RESULTS: Relapsed SCLC samples demonstrated mutation signatures characteristic of platinum and APOBEC mutagenesis. Furthermore, these samples were characterized by MYC, MYCL and MYCN amplifications. Both treatment-naive and relapsed SCLC samples showed high prevalence of mutation-associated neoantigens (median= 86 in treatment-naive and 90 in relapsed SCLC; p=0.8) and TP53 was the most frequently altered gene to result in a neoantigen (48% of analyzed samples). Potential mechanisms of immune evasion, including amplification of CD24, overexpression of IDO1, increased M2 macrophage presence, and upregulation of HLA-E were also observed in relapse samples. Differences in alternative splicing patterns were observed between treatment-naive and relapsed small cell samples. Retained intron events were significantly enriched in treatment-naive samples and affected genes involved in DNA repair, metabolism, and WNT and MYC pathways. CONCLUSIONS: This study highlights the genomic and transcriptomic features of relapsed SCLC. These samples were characterized by genomic instability, WNT and MYC dysregulation, and splicing aberrations. Additional studies targeting the splicing machinery, WNT signaling, and immune evasion pathways could identify novel therapeutic vulnerabilities in SCLC.

Journal Article↗

jsPCA: fast, scalable, and interpretable identification of spatial domains and variable genes across multi-slice and multi-sample spatial transcriptomics data.

MOTIVATION: Spatial transcriptomics technologies record genome-wide measurements of gene expression with high spatial resolution. These technologies generate large and high-dimensional datasets requiring efficient automated methods for their analysis. We introduce joint spatial PCA (jsPCA), a novel, fast, scalable and interpretable method for the automatic identification of spatial domains and variable genes in multi-slice and multi-sample spatial transcriptomics data. RESULTS: jsPCA relies on a simple mathematical formulation of a spatial covariance defined as the product of the gene expression covariance with the spatial autocorrelation. The principal components of this spatial covariance yield a biologically meaningful low-dimensional representation. From this representation, spatial domains are derived by simple clustering and spatially variable genes are identified directly from the principal component coefficients. A joint representation of multiple slices and samples without spatial alignment is obtained by computing common principal components via joint diagonalization. By leveraging data sparsity and non-convex manifold optimization, jsPCA leads to computing time in the order of seconds to minutes, substantially outperforming state-of-the-art approaches. We benchmarked jsPCA against 10 state-of-the-art methods on two reference databases. Our approach demonstrated excellent performance, comparable or better than state-of-the-art methods, while being much faster, interpretable, and scalable to very large datasets.

Journal Article↗

Admission whole-blood transcriptomic characterization of a neutrophil-predominant systemic immune response in patients with acute traumatic brain injury.

BACKGROUND: Acute traumatic brain injury (TBI) is accompanied by systemic immune responses, but their whole-blood transcriptomic features at hospital arrival remain incompletely characterized. We aimed to characterize these features in patients with acute TBI compared with healthy controls. METHODS: In this single-center prospective observational study, we performed whole-blood RNA sequencing on hospital-arrival samples from 42 patients with acute TBI and 21 healthy controls. Analyses included differential expression (limma-voom; FDR < 0.05, |log2FC| > 0.7), functional enrichment, Ingenuity Pathway Analysis, CIBERSORTx LM22 deconvolution, and per-sample neutrophil degranulation signature scoring. RESULTS: Differential expression analysis identified 996 upregulated and 863 downregulated genes, with marked upregulation of inflammation-, innate immunity-, and neutrophil-related genes including DUSP1, HMGB2, MMP9, and S100A8. Canonical pathways with positive IPA z-scores included Neutrophil degranulation, Neutrophil Extracellular Trap Signaling Pathway, and Toll-like Receptor Signaling; upstream regulators included TNF, IL1B, IFNG, and STAT3. Deconvolution identified 7 of 22 differing subsets (q < 0.05), with relatively higher myeloid and lower lymphoid fractions in TBI. The Neutrophil degranulation signature score correlated with Injury Severity Score within TBI (Spearman &#x3c1; = +0.55; q < 0.001). CONCLUSIONS: Admission whole-blood transcriptomics characterized a neutrophil-predominant systemic transcriptional response in patients with acute TBI. This response was also evident among patients without major extracranial injury and was associated with total ISS. However, because the study lacked an appropriately matched non-TBI trauma comparator, the findings should be interpreted as a descriptive characterization of a systemic injury response accompanying TBI and do not establish a TBI-specific molecular signature or mechanism.

gene expression↗

SpaceBar enables clone tracing in spatial transcriptomic data.

We report a cellular barcoding strategy, SpaceBar, that enables simultaneous clone tracing and spatial transcriptomics profiling. Our approach uses a library of 96 synthetic barcode sequences that can be robustly detected by imaging based spatial transcriptomics (seqFISH), delivered such that each cell is labeled with a combination of barcodes. We used these barcodes to label melanoma cells in a tumor xenograft model and profiled both clone identity and spatial gene expression in situ. We developed a gene scoring metric that quantifies how strongly gene expression is driven by intrinsic cellular cues or extrinsic environmental signals. Our framework distinguishes between clonal dynamics and environmentally-driven transcriptional regulation in complex tissue contexts.

Journal Article↗

Single-cell-scale spatial transcriptome of the developing and adult mouse ovary.

Mammalian ovary development is essential for female fertility, involving the complex spatial patterning of diverse cell types to establish the finite reserve of ovarian follicles. While single-cell transcriptome analyses have provided important insights into the mechanisms driving specification and developmental trajectories of ovarian cells, they disrupt this crucial spatial context. To overcome this limitation, we used 10X Genomics Visium HD spatial transcriptomics to analyze the developing mouse ovary while maintaining its native cellular architecture. We captured all ovarian cell types at eight key fetal and postnatal timepoints, generating a near single cell resolution library of spatial gene expression across ovarian development. This comprehensive dataset allows analysis of dynamic transcriptional signatures associated with unique spatial patterning throughout development, including the establishment of cortex and medulla and assembly of ovarian follicles in each region. This dataset represents a fundamental resource for the investigation of regulatory mechanisms driving spatial patterning of the ovary and opens new avenues to explore the spatial determinants of female fertility and reproductive longevity.

Journal Article↗

Comprehensive Transcriptome Annotation of Thousands of HIV-1 Genomes.

Alternative splicing in HIV-1 has been a central focus of decades of research, uncovering key mechanisms of viral gene regulation, immune evasion, and therapeutic response - yet, no reference resource has existed to support transcriptome-wide analysis, limiting adoption of modern computational methods. We present HIV Atlas (https://ccb.jhu.edu/HIV_Atlas), the first reference-quality annotation of HIV-1 and SIV transcriptional diversity. We manually curated transcriptomes for HIV-1HXB2 and SIVmac239 and developed Vira, an automated annotation-transfer method specifically designed to address unique challenges of viral genome biology, to generate high-quality annotations for 2,077 complete HIV-1 genomes. Using the resources presented in our work, we evaluated conservation of splice sites, revealing near-perfect preservation of major donors and acceptors. Furthermore, using several public datasets, we demonstrate how HIV Atlas enhances methodology, improves the quality and novelty of results, and opens novel avenues for research, supporting more accurate and comprehensive analyses of bulk, single-cell, and spatial RNA-seq in HIV-1 studies.

Journal Article↗

The mouse brain transcriptome by SAGE: differences in gene expression between P30 brains of the partial trisomy 16 mouse model of Down syndrome (Ts65Dn) and normals.

Trisomy 21, or Down syndrome (DS), is the most common genetic cause of mental retardation. Changes in the neuropathology, neurochemistry, neurophysiology, and neuropharmacology of DS patients' brains indicate that there is probably abnormal development and maintenance of central nervous system structure and function. The segmental trisomy mouse (Ts65Dn) is a model of DS that shows analogous neurobehavioral defects. We have studied the global gene expression profiles of normal and Ts65Dn male and normal female mice brains (P30) using the serial analysis of gene expression (SAGE) technique. From the combined sample we collected a total of 152,791 RNA tags and observed 45,856 unique tags in the mouse brain transcriptome. There are 14 ribosomal protein genes (nine under expressed) among the 330 statistically significant differences between normal male and Ts65Dn male brains, which possibly implies abnormal ribosomal biogenesis in the development and maintenance of DS phenotypes. This study contributes to the establishment of a mouse brain transcriptome and provides the first overall analysis of the differences in gene expression in aneuploid versus normal mammalian brain cells.

Animals↗