Transcriptome-wide Expression Research Articles

Dysregulation of mRNA expression by hsa-miR-186 overexpression in arsenic-induced skin carcinogenesis.

Dysregulated miRNA expression contributes to development of arsenic-induced cutaneous squamous cell carcinoma (cSCC). hsa-miR-186 (miR-186) is overexpressed in arsenical cSCC tissues as well as in preclinical cell line model of arsenical cSCC. Simultaneous miR-186 overexpression and chronic inorganic trivalent arsenite (iAs; 100nM) exposure transformed human HaCaT cell line preferentially over miR-186 overexpression or iAs exposure alone. Both iAs and miR-186 regulate the expression of wide range of mRNA targets. However, how their interaction impacts the transcriptome-wide mRNA expression landscape ushering in cancer is unknown. We performed longitudinal RNA-seq analysis in passage-matched HaCaT cell clones (±miR-186 overexpression) with simultaneous chronic iAs exposure (0/100nM) at 12 and 29weeks. We determined the impact of each factor and their interaction towards differential gene expression and pathway dysregulation employing two different statistical approaches (t-statistic and 2-factor ANOVA). We show that a core set of pathways are dysregulated deterministically irrespective of the statistical approach chosen, possibly representing necessary changes for transformation. The data suggest that each clonal line could take a unique route to dysregulate this core set of pathways necessary for transformation, highlighting the possible role of stochasticity in cancer development. Evidence is presented to sift the strengths and weaknesses of each statistical methodology in providing biological understanding of events that play crucial roles in carcinogenesis in large datasets with multiple contributing variables.

Transcriptome-wide identification and expression analysis of expansin genes in big and small leaf-morphotypes of Adhatoda vasica

Transcriptome-wide identification and expression analysis of expansin genes in big and small leaf-morphotypes of Adhatoda vasica

Genome- and Transcriptome-Wide Characterization and Expression Analyses of bHLH Transcription Factor Family Reveal Their Relevance to Salt Stress Response in Tomato.

The bHLH (basic helix-loop-helix) transcription factors function as crucial regulators in numerous biological processes including abiotic stress responses and plant development. According to our RNA-seq analysis of tomato seedlings under salt stress, we found that, although the bHLH gene family in tomato has been studied, there are still so many tomato bHLH genes that have not been identified and named, which will hinder the later study of SlbHLHs. In total, 195 SlbHLHs that were unevenly distributed onto 12 chromosomes were identified from the tomato genome and were classified into 27 subfamilies based on their molecular features. The collinearity between SlbHLHs and interrelated orthologs from 10 plants further revealed evolutionary insights into SlbHLHs. Cis-element investigations of SlbHLHs promotors further suggested the potential roles of SlbHLHs in tomato development and stress responses. A total of 30 SlbHLHs were defined as the differentially expressed genes in response to salt stress by RNA-seq. The expression profiles of selected SlbHLHs were varyingly and markedly induced by multiple abiotic stresses and hormone treatments. These results provide valuable information to further understand the significance and intricacy of the bHLH transcription factor family, and lay a foundation for further exploring functions and possible regulatory mechanisms of SlbHLH members in abiotic stress tolerance, which will be significant for the study of tomato stress resistance and agricultural productivity.

Transcriptome-wide RNA 5-methylcytosine profiles of human iPSCs and iPSC-derived cardiomyocytes

Cardiac regenerative therapy has recently progressed by reprogramming somatic cells into induced pluripotent stem cells (iPSCs) and advanced by large-scale differentiation-derived cardiomyocytes (hiPSC-CMs). However, repairing damaged cardiac tissues with hiPSC-CMs remains limited due to immune rejection, cardiac arrhythmias, and concerns over tumor formation after hiPSC-CM transplantation. Despite efforts in profiling epigenomic changes during cardiac differentiation, regulatory mechanisms underlying 5-methylcytosine (m5C) deposition in RNA m5C epitranscriptomic landscape during hiPSC-to-cardiomyocyte differentiation remain unclear. Herein, bisulfite RNA-sequencing analysis was conducted in human pluripotent stem cells (hPSCs) from three independent cellular origins, and their derived cardiomyocytes (hPSC-CM), metabolic-maturation of derived cardiomyocytes (hPSC-CM-lac) and biochemical-enhanced derived cardiomyocytes (hPSC-CM-TDI). Integrated analysis of differentially methylated RNA m5C profiles and transcriptome-wide expression facilitated the identification of m5C sites coupled to the cardiomyocyte differentiation and RNA-dependent regulatory mechanisms of stem cell pluripotency. The RNA m5C profiles in this dataset allow the evaluations of the m5C level and distribution of specific m5C loci and facilitate understanding of the m5C epitranscriptomic landscape in biological functions of hPSC-CM beyond in vivo transplantation barriers.

A luminal non-coding RNA-based genomic classifier confirms favourable outcomes in patients with clinically organ-confined bladder cancer treated with radical cystectomy.

To further evaluate a genomic classifier (GC) in a cohort of patients undergoing radical cystectomy (RC), as long non-coding RNA (lncRNA)-based genomic profiling has suggested utility in identifying a distinct tumour subgroup corresponding to a favourable prognosis in patients with bladder cancer. Transcriptome-wide expression profiling using Decipher Bladder was performed on transurethral resection of bladder tumour samples from a cohort of patients with high-grade, clinically organ-confined (cTa-T2N0M0) urothelial carcinoma (UC) who subsequently underwent RC without any neoadjuvant therapy (n = 226). The lncRNA-based luminal favourable status was determined using a previously developed GC. The primary endpoint was overall survival (OS) after RC. Secondary endpoints included cancer-specific mortality and upstaging at RC. In the study, 134 patients were clinical non-muscle-invasive bladder cancer (cTa/Tis/T1) and 92 patients were cT2. We identified 60 patients with luminal favourable subtype, all of which showed robust gene expression patterns associated with less aggressive bladder cancer biology. On multivariate analysis, patients with the luminal favourable subtype (vs without) were significantly associated with lower odds of upstaging to pathological (p)T3+ disease (odds ratio [OR] 0.32, 95% confidence interval [CI] 0.12-0.82; P = 0.02), any upstaging (OR 0.41, 95% CI 0.20-0.83; P = 0.01), and any upstaging and/or pN+ (OR 0.50, 95% CI 0.25-1.00; P = 0.05). Luminal favourable bladder cancer was significantly associated with better OS (hazard ratio 0.33, 95% CI 0.15-0.74; P = 0.007). This study validates the performance of the GC for identifying UCs with a luminal favourable subtype, harbouring less aggressive tumour biology.

Neural network-assisted humanisation of COVID-19 hamster transcriptomic data reveals matching severity states in human disease

Translating findings from animal models to human disease is essential for dissecting disease mechanisms, developing and testing precise therapeutic strategies. The coronavirus disease 2019 (COVID-19) pandemic has highlighted this need, particularly for models showing disease severity-dependent immune responses. Single-cell transcriptomics (scRNAseq) is well poised to reveal similarities and differences between species at the molecular and cellular level with unprecedented resolution. However, computational methods enabling detailed matching are still scarce. Here, we provide a structured scRNAseq-based approach that we applied to scRNAseq from blood leukocytes originating from humans and hamsters affected with moderate or severe COVID-19. Integration of data from patients with COVID-19 with two hamster models that develop moderate (Syrian hamster, Mesocricetus auratus) or severe (Roborovski hamster, Phodopus roborovskii) disease revealed that most cellular states are shared across species. A neural network-based analysis using variational autoencoders quantified the overall transcriptomic similarity across species and severity levels, showing highest similarity between neutrophils of Roborovski hamsters and patients with severe COVID-19, while Syrian hamsters better matched patients with moderate disease, particularly in classical monocytes. We further used transcriptome-wide differential expression analysis to identify which disease stages and cell types display strongest transcriptional changes. Consistently, hamsters' response to COVID-19 was most similar to humans in monocytes and neutrophils. Disease-linked pathways found in all species specifically related to interferon response or inhibition of viral replication. Analysis of candidate genes and signatures supported the results. Our structured neural network-supported workflow could be applied to other diseases, allowing better identification of suitable animal models with similar pathomechanisms across species. This work was supported by German Federal Ministry of Education and Research, (BMBF) grant IDs: 01ZX1304B, 01ZX1604B, 01ZX1906A, 01ZX1906B, 01KI2124, 01IS18026B and German Research Foundation (DFG) grant IDs: 14933180,431232613.

Integrative analysis of transcriptome-wide association study and mRNA expression profile identified risk genes for bipolar disorder

ObjectiveBipolar disorder (BD) is a debilitating neuropsychiatric disorder, which is associated with genetic variation through “vast but mixed” Genome-Wide Association Studies (GWAS). Transcriptome-Wide Association Study (TWAS) is more effective in explaining genetic factors that influence complex diseases and can help identifying risk genes more reliably. So, this study aims to identify potential BD risk genes in pedigrees with TWAS. MethodsWe conducted a TWAS analysis with expression quantitative trait loci (eQTL) analysis on extended BD pedigrees, and the BD genome-wide association study (GWAS) summary data acquired from the Psychiatric Genomics Consortium (PGC). Furthermore, the BD-associated genes identified by TWAS were validated by mRNA expression profiles from the Gene Expression Omnibus (GEO) Datasets (GSE23848 and GSE46416). Functional enrichment and annotation analysis were implemented by RStudio (version 4.2.0). ResultsTWAS identified 362 genes with P value < 0.05, and 18 genes remain significant after Bonferroni correction, such as SEMA3G (PTWAS=1.07 × 10-11), ALOX5AP (PTWAS=3.12 × 10-8), and PLEC (PTWAS=1.27 × 10-7). Further 6 overlapped genes were detected in integrative analysis, such as UQCRB (PTWAS=0.0020, PmRNA=0.0000), TMPRSS9 (PTWAS=0.0405, PmRNA=0.0032), and SNX10 (PTWAS=0.0104, PmRNA=0.0015). Using genes identified by TWAS, Gene Ontology (GO) enrichment analysis identified 40 significant GO terms, such as mitochondrial ATP synthesis coupled electron transport, mitochondrial respiratory, aerobic electron transport chain, oxidative phosphorylation, mitochondrial membrane proteins, and ubiquinone activity. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway enrichment analysis identified significant 15 pathways for BD, such as Oxidative phosphorylation, endocannabinoids signaling, neurodegeneration, and reactive oxide species. ConclusionsWe found a set of BD-associated genes and pathways, validating the important role of neurodevelopmental abnormalities, inflammatory responses, and mitochondrial dysfunction in the pathology of BD, offering novel information for comprehending the genetic basis of BD.

Systematic selection of suitable reference genes for quantitative real-time PCR normalization studies of gene expression in Lutjanuserythropterus

Quantitative real-time PCR (qRT-PCR) has been widely employed for the study of gene expression in fish, and accurate normalization is crucial. In this study, we aimed to identify the most stably expressed genes in various tissues, different developmental stages, and within astaxanthin treatment groups in Lutjanuserythropterus. Twelve candidate genes (EEF1A, CYB5R3, DLD, IDH3A, MRPL17, MRPL43, NDUFS7, PABPC1, PAGR1, PFDN2, PSMC3, and RAB10) were examined via qRT-PCR. We employed geNorm and NormFinder to assess their stability. The results revealed that RAB10 and PFDN2 exhibited relatively stable expression patterns across different tissue and astaxanthin treatment groups, while NDUFS7 and MRPL17 proved to be the most reliable reference gene combinations across various developmental stages. The stability of these selected genes was further validated by assessing the expression of two target genes, CRADD and CAPNS1, across developmental stages, reinforcing the reliability of NDUFS7 as it closely aligned with transcriptome-wide expression patterns at these stages. The present results will help researchers to obtain more accurate results in future qRT-PCR analysis in L.erythropterus.

Transcriptome-wide expression landscape and starch synthesis pathway co-expression network in sorghum.

The gene expression landscape across different tissues and developmental stages reflects their biological functions and evolutionary patterns. Integrative and comprehensive analyses of all transcriptomic data in an organism are instrumental to obtaining a comprehensive picture of gene expression landscape. Such studies are still very limited in sorghum, which limits the discovery of the genetic basis underlying complex agricultural traits in sorghum. We characterized the genome-wide expression landscape for sorghum using 873 RNA-sequencing (RNA-seq) datasets representing 19 tissues. Our integrative analysis of these RNA-seq data provides the most comprehensive transcriptomic atlas for sorghum, which will be valuable for the sorghum research community for functional characterizations of sorghum genes. Based on the transcriptome atlas, we identified 595 housekeeping genes (HKGs) and 2080 tissue-specific expression genes (TEGs) for the 19 tissues. We identified different gene features between HKGs and TEGs, and we found that HKGs have experienced stronger selective constraints than TEGs. Furthermore, we built a transcriptome-wide co-expression network (TW-CEN) comprising 35 modules with each module enriched in specific Gene Ontology terms. High-connectivity genes in TW-CEN tend to express at high levels while undergoing intensive selective pressure. We also built global and seed-preferential co-expression networks of starch synthesis pathways, which indicated that photosynthesis and microtubule-based movement play important roles in starch synthesis. The global transcriptome atlas of sorghum generated by this study provides an important functional genomics resource for trait discovery and insight into starch synthesis regulation in sorghum.

Synergistic response to climate stressors in coral is associated with genotypic variation in baseline expression.

As environments are rapidly reshaped due to climate change, phenotypic plasticity plays an important role in the ability of organisms to persist and is considered an especially important acclimatization mechanism for long-lived sessile organisms such as reef-building corals. Often, this ability of a single genotype to display multiple phenotypes depending on the environment is modulated by changes in gene expression, which can vary in response to environmental changes via two mechanisms: baseline expression and expression plasticity. We used transcriptome-wide expression profiling of eleven genotypes of common-gardened Acropora cervicornis to explore genotypic variation in the expression response to thermal and acidification stress, both individually and in combination. We show that the combination of these two stressors elicits a synergistic gene expression response, and that both baseline expression and expression plasticity in response to stress show genotypic variation. Additionally, we demonstrate that frontloading of a large module of coexpressed genes is associated with greater retention of algal symbionts under combined stress. These results illustrate that variation in the gene expression response of individuals to climate change stressors can persist even when individuals have shared environmental histories, affecting their performance under future climate change scenarios.

Abstract 7568: Molecular effects of neo-adjuvant hormone therapy for locally advanced prostate cancer: An ancillary study of the ARNEO phase 2 trial

Abstract Intensified hormone therapy in combination with radiotherapy improves outcome for very high-risk locally advanced prostate cancer. The impact of hormone therapy with surgery, however, is uncertain. In order to understand this better, we ran an ancillary translational study on ARNEO (NCT03080116), a Phase 2 clinical trial, which randomized 89 locally-advanced prostate cancer patients to 12 weeks LHRHa +/- apalutamide (APA) followed by prostatectomy. We performed clinical-grade transcriptome-wide expression array analysis (Veracyte, Inc. San Diego, CA) of FFPE diagnostic core biopsies (N=89) and tumors removed at prostatectomy (N=32 LHRHa, N=23 LHRHa + APA). We also performed RNASeq on fresh frozen biopsies of prostatectomy tumors (N=42 LHRHa, N=45 LHRHa + apalutamide) and LHRHa-naive control cases (N=61) matched for clinical and pathological risk factors. Pre-surgery serum PSA decreased by a median of 98% (IQR:90.2-99.3) after hormone treatment. This decrease was better associated with tumor AR signaling (r2=0.11) versus pre-treatment serum PSA (r2=0.0005), or between pre-treatment serum PSA with diagnostic biopsies (r2=0.008). We confirmed APA more potently suppressed AR signaling than LHRHa alone, with no detectable difference in direction of effect for differentially regulated genes/pathways. LHRHa + APA had less residual tumor (median Residual Cancer Burden, RCB=0.48cm3) than in LHRHa (median RCB=1.7cm3) but in cases with sufficient tumor, we observed residual active AR signaling pathways at prostatectomy in both treatment arms. To orthogonally confirm this, we selected 28 ERG positive cases (by IHC pre-treatment) and found retained ERG expression in 12/13 on LHRHa and 10/15 on LHRHa + APA. Using pathway analysis in the RNAseq data, we identified LHRHa +/- APA induced down-regulation of AR, metabolic pathways, as well as DNA repair, and up-regulation of hallmark NOTCH, Hedgehog, TGF-beta, WNT and immune signaling (adjusted p-value cutoff=0.001). Comparing differentially expressed genes on RNASeq between treatment-naïve versus LHRHa +/- APA treated tumors, we observed STING1 as significantly higher in LHRHa +/- APA vs control (Wilcoxon rank sum test p-value &amp;lt;0.001). We found that expression of STING1 was positively correlated with hallmark interferon alpha response (r=0.7, rs=0.6) that was correlated with apoptosis (r=0.5, rs=0.5). Whilst we show that adding APA to LHRHa increases suppression of AR activity and ERG expression, we find that 12 weeks is insufficient to abrogate AR signaling in most tumors, despite ~98% decrease in serum PSA. We identify novel processes modulated in humans by potent suppression of androgen activity, including upregulation of STING1-type 1 interferon response as a potential key modulator of treatment effect. Citation Format: Mazlina Ismail, Wout Devlies, Marina Parry, Gaetan Devos, Stefanie Friedrich, Charles Parker, Larissa Mendes, Elai Davicioni, Mark Linch, Frank Claessens, Steven Joniau, Gerhardt Attard. Molecular effects of neo-adjuvant hormone therapy for locally advanced prostate cancer: An ancillary study of the ARNEO phase 2 trial [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7568.

Chronic arsenic exposure induces malignant transformation of human HaCaT cells through both deterministic and stochastic changes in transcriptome expression

Biological processes are inherently stochastic, i.e., are partially driven by hard to predict random probabilistic processes. Carcinogenesis is driven both by stochastic and deterministic (predictable non-random) changes. However, very few studies systematically examine the contribution of stochastic events leading to cancer development. In differential gene expression studies, the established data analysis paradigms incentivize expression changes that are uniformly different across the experimental versus control groups, introducing preferential inclusion of deterministic changes at the expense of stochastic processes that might also play a crucial role in the process of carcinogenesis. In this study, we applied simple computational techniques to quantify: (i) The impact of chronic arsenic (iAs) exposure as well as passaging time on stochastic gene expression and (ii) Which genes were expressed deterministically and which were expressed stochastically at each of the three stages of cancer development. Using biological coefficient of variation as an empirical measure of stochasticity we demonstrate that chronic iAs exposure consistently suppressed passaging related stochastic gene expression at multiple time points tested, selecting for a homogenous cell population that undergo transformation. Employing multiple balanced removal of outlier data, we show that chronic iAs exposure induced deterministic and stochastic changes in the expression of unique set of genes, that populate largely unique biological pathways. Together, our data unequivocally demonstrate that both deterministic and stochastic changes in transcriptome-wide expression are critical in driving biological processes, pathways and networks towards clonal selection, carcinogenesis, and tumor heterogeneity.

Transcriptome-Wide Identification and Expression Analysis of bHLH Family Genes in Iris domestica under Drought and Cu Stress.

The role of bHLH transcription factors in plant response to abiotic stress and regulation of flavonoid metabolism is well documented. However, to date, the bHLH transcription factor family in Iris domestica remains unreported, impeding further research on flavonoid metabolism in this plant. To address this knowledge gap, we employed bioinformatics to identify 39 IdbHLH genes and characterised their phylogenetic relationships and gene expression patterns under both drought and copper stress conditions. Our evolutionary tree analysis classified the 39 IdbHLHs into 17 subfamilies. Expression pattern analysis revealed that different IdbHLH transcription factors had distinct expression trends in various organs, suggesting that they might be involved in diverse biological processes. We found that IdbHLH36 was highly expressed in all organs (Transcripts Per Million (TPM) > 10), while only 12 IdbHLH genes in the rhizome and four in the root were significantly upregulated under drought stress. Of these, four genes (IdbHLH05, -37, -38, -39) were co-upregulated in both the rhizome and root, indicating their potential role in drought resistance. With regards to copper stress, we found that only 12 genes were upregulated. Further co-expression analysis revealed that most bHLH genes were significantly correlated with key enzyme genes involved in isoflavone biosynthesis. Thereinto, IdbHLH06 showed a significant positive correlation with IdC4H1 and Id4CL1 (p < 0.05). Furthermore, a transient expression assay confirmed that the IdbHLH06 protein was localised in the nucleus. Our findings provide new insights into the molecular basis and regulatory mechanisms of bHLH transcription factors in isoflavone biosynthesis in I. domestica.

The role of androgen response pathway in association with tumor biology and response to neoadjuvant immune-checkpoint inhibitors (ICI) in muscle-invasive urothelial bladder carcinoma (MIBC).

662 Background: Emerging findings from ICI trials are suggestive of differential outcomes in females compared with males. Sex-specific differences in tumor biology, microbiome, and outcomes in MIBC are well-established. Conversely, limited data are available regarding the role of the androgen response pathway in explaining subtype differences between males and females and, most importantly, response to ICI. Methods: We retrospectively evaluated the biomarker findings of the transurethral resection of the bladder tumors (TURBT) and baseline clinical characteristics of patients (pts) with cT2-4N0M0 MIBC who received neoadjuvant pembrolizumab and radical cystectomy (RC) in the PURE-01 study. Transcriptome-wide expression profiling was performed on 102 TURBT samples. We focused on androgen receptor (AR) gene expression and the hallmark androgen response pathway expression quantifying 117 different genes related to the androgen response pathway as previously reported (PMID: 26771021). Pt and tumor characteristics were compared between subgroups using χ2 tests and two-sided Wilcoxon rank-sum tests. Kaplan-Meier curves and log rank tests evaluated gene expression and signature scores in association with the event-free survival (EFS), calculated from start of pembrolizumab treatment. Signature scores were analyzed as continuous values and were split by median and quartile values. Results: Pts received pembrolizumab and RC from 2017 to 2022. There were 87 (85.3%) males and 15 (14.7%) females. AR gene expression was not significantly different between males and females (p = 0.5), whereas androgen response signature was significantly higher in male (p = 0.03). The AR gene expression and androgen response signature scores were more highly expressed in luminal tumors vs rest of molecular subtypes (p = 0.005 and p &lt; 0.001). Despite AR gene expression was not associated with ypT0N0 response (p = 0.425), hallmark androgen response signature scores were significantly lower in ypT0N0 responders (p = 0.033). Higher hallmark androgen response signature scores were statistically significantly associated with shorter EFS (p = 0.028), whereas AR gene expression was not (p = 0.152). Conclusions: For the first time to our knowledge, we reported the androgen response pathway expression to be a potential biomarker of ICI benefit in MIBC. This signature was also associated with a luminal subtype, linked to inferior responsiveness to ICI. If confirmed at a larger scale these analyses will strengthen the need for sex-specific approaches and potentially new combination therapies in MIBC.

Molecularly defined and spatially resolved cell atlas of the whole mouse brain

In mammalian brains, millions to billions of cells form complex interaction networks to enable a wide range of functions. The enormous diversity and intricate organization of cells have impeded our understanding of the molecular and cellular basis of brain function. Recent advances in spatially resolved single-cell transcriptomics have enabled systematic mapping of the spatial organization of molecularly defined cell types in complex tissues1–3, including several brain regions (for example, refs. 1–11). However, a comprehensive cell atlas of the whole brain is still missing. Here we imaged a panel of more than 1,100 genes in approximately 10 million cells across the entire adult mouse brains using multiplexed error-robust fluorescence in situ hybridization12 and performed spatially resolved, single-cell expression profiling at the whole-transcriptome scale by integrating multiplexed error-robust fluorescence in situ hybridization and single-cell RNA sequencing data. Using this approach, we generated a comprehensive cell atlas of more than 5,000 transcriptionally distinct cell clusters, belonging to more than 300 major cell types, in the whole mouse brain with high molecular and spatial resolution. Registration of this atlas to the mouse brain common coordinate framework allowed systematic quantifications of the cell-type composition and organization in individual brain regions. We further identified spatial modules characterized by distinct cell-type compositions and spatial gradients featuring gradual changes of cells. Finally, this high-resolution spatial map of cells, each with a transcriptome-wide expression profile, allowed us to infer cell-type-specific interactions between hundreds of cell-type pairs and predict molecular (ligand–receptor) basis and functional implications of these cell–cell interactions. These results provide rich insights into the molecular and cellular architecture of the brain and a foundation for functional investigations of neural circuits and their dysfunction in health and disease.

Genome- and Transcriptome-wide Association Studies to Discover Candidate Genes for Diverse Root Phenotypes in Cultivated Rice

Root system architecture plays a crucial role in nutrient and water absorption during rice production. Genetic improvement of the rice root system requires elucidating its genetic control. Genome-wide association studies (GWASs) have identified genomic regions responsible for rice root phenotypes. However, candidate gene prioritization around the peak region often suffers from low statistical power and resolution. Transcriptomics enables other statistical mappings, such as transcriptome-wide association study (TWAS) and expression GWAS (eGWAS), which improve candidate gene identification by leveraging the natural variation of the expression profiles. To explore the genes responsible for root phenotypes, we conducted GWAS, TWAS, and eGWAS for 12 root phenotypes in 57 rice accessions using 427,751 single nucleotide polymorphisms (SNPs) and the expression profiles of 16,901 genes expressed in the roots. The GWAS identified three significant peaks, of which the most significant peak responsible for seven root phenotypes (crown root length, crown root surface area, number of crown root tips, lateral root length, lateral root surface area, lateral root volume, and number of lateral root tips) was detected at 6,199,732 bp on chromosome 8. In the most significant GWAS peak region, OsENT1 was prioritized as the most plausible candidate gene because its expression profile was strongly negatively correlated with the seven root phenotypes. In addition to OsENT1, OsEXPA31, OsSPL14, OsDEP1, and OsDEC1 were identified as candidate genes responsible for root phenotypes using TWAS. Furthermore, a cis-eGWAS peak SNP was detected for OsDjA6, which showed the eighth strongest association with lateral root volume in the TWAS. The cis-eGWAS peak SNP for OsDjA6 was in strong linkage disequilibrium (LD) with a GWAS peak SNP on the same chromosome for lateral root volume and in perfect LD with another SNP variant in a putative cis-element at the 518 bp upstream of the gene. These candidate genes provide new insights into the molecular breeding of root system architecture.

Integrating spatial and single-cell transcriptomics data using deep generative models with SpatialScope

The rapid emergence of spatial transcriptomics (ST) technologies is revolutionizing our understanding of tissue spatial architecture and biology. Although current ST methods, whether based on next-generation sequencing (seq-based approaches) or fluorescence in situ hybridization (image-based approaches), offer valuable insights, they face limitations either in cellular resolution or transcriptome-wide profiling. To address these limitations, we present SpatialScope, a unified approach integrating scRNA-seq reference data and ST data using deep generative models. With innovation in model and algorithm designs, SpatialScope not only enhances seq-based ST data to achieve single-cell resolution, but also accurately infers transcriptome-wide expression levels for image-based ST data. We demonstrate SpatialScope’s utility through simulation studies and real data analysis from both seq-based and image-based ST approaches. SpatialScope provides spatial characterization of tissue structures at transcriptome-wide single-cell resolution, facilitating downstream analysis, including detecting cellular communication through ligand-receptor interactions, localizing cellular subtypes, and identifying spatially differentially expressed genes.

Spatial transcriptomics and the anatomical pathologist: Molecular meets morphology.

In recent years anatomical pathology has been revolutionised by the incorporation of molecular findings into routine diagnostic practice, and in some diseases the presence of specific molecular alterations are nowessential for diagnosis. Spatial transcriptomics describes a group of technologies that provide up to transcriptome-wide expression profiling while preserving the spatial origin of the data, with many of these technologies able to provide these data using a single tissue section. Spatial transcriptomics allows expression profiling of highly specific areas within a tissue section potentially to subcellular resolution, and allows correlation of expression data with morphology, tissue type and location relative to other structures. While largely still research laboratory-based, several spatial transcriptomics methods have now achieved compatibility with formalin-fixed paraffin-embedded tissue (FFPE), allowing their use in diagnostic tissue samples, and with further development potentially leading to their incorporation in routine anatomical pathology practice. This mini review provides an overview of spatial transcriptomics methods, with an emphasis on platforms compatible with FFPE tissue, approaches to assess the data and potential applications in anatomical pathology practice.

Vesical Imaging-Reporting and Data System use predicting the outcome of neoadjuvant pembrolizumab in muscle-invasive bladder cancer.

To evaluate the predictive capability of the pre- and post-pembrolizumab Vesical Imaging-Reporting and Data System (VI-RADS) to identify ypT0N0 or ypT≤1N0 response in muscle-invasive bladder cancer (MIBC) within the PURE-01 trial (ClinicalTrials.gov identifier: NCT02736266). Patients were staged with bladder multiparametric magnetic resonance imaging (mpMRI) before and after treatment (three cycles of pembrolizumab) prior to radical cystectomy (RC). Logistic regression models were used to analyse the pre- and post- pembrolizumab VI-RADS against ypT≤1N0 and ypT0N0 response. The VI-RADS scores were dichotomised between 0 and 3 (0 = no evidence of disease) and 4-5. Event-free survival (EFS) and overall survival (OS) analyses were performed. Comprehensive genomic profiling and transcriptome-wide expression profiling data were matched with the VI-RADS scores. In total, 110 patients underwent centrally reviewed scans (N = 220 mpMRI), treated between February 2017 and July 2020. Both pre- and post-pembrolizumab VI-RADS 0-3 scores were the only significant covariates that predicted the ypT≤1N0 endpoint in multivariable analyses, and the strongest effect was seen with post-pembrolizumab VI-RADS 0-3 predicting the ypT≤1N0 response (P < 0.001). The area under the curve for this model was 0.90. Post-pembrolizumab VI-RADS 0-3 also predicted a longer EFS (P < 0.001) and OS (P = 0.044). The scores of several gene signatures from baseline tumours differed between the pre-pembrolizumab VI-RADS 0-3 and 4-5 categories. Post-pembrolizumab VI-RADS scores are strongly associated with pathological downstaging and survival. VI-RADS scores were also characterised by distinct biomarker features. These results indicate that the VI-RADS is emerging as an important tool for designing next-generation trials for MIBC.

Consequences of light spectra for pigment composition and gene expression in the cryptophyte Rhodomonas salina.

Algae with a more diverse suite of pigments can, in principle, exploit a broader swath of the light spectrum through chromatic acclimation, the ability to maximize light capture via plasticity of pigment composition. We grew Rhodomonas salina in wide-spectrum, red, green, and blue environments and measured how pigment composition differed. We also measured expression of key light-capture and photosynthesis-related genes and performed a transcriptome-wide expression analysis. We observed the highest concentration of phycoerythrin in green light, consistent with chromatic acclimation. Other pigments showed trends inconsistent with chromatic acclimation, possibly due to feedback loops among pigments or high-energy light acclimation. Expression of some photosynthesis-related genes was sensitive to spectrum, although expression of most was not. The phycoerythrin α-subunit was expressed two-orders of magnitude greater than the β-subunit even though the peptides are needed in an equimolar ratio. Expression of genes related to chlorophyll-binding and phycoerythrin concentration were correlated, indicating a potential synthesis relationship. Pigment concentrations and expression of related genes were generally uncorrelated, implying post-transcriptional regulation of pigments. Overall, most differentially expressed genes were not related to photosynthesis; thus, examining associations between light spectrum and other organismal functions, including sexual reproduction and glycolysis, may be important.