SEARCH · PubMed Health

Results for “RNAseq”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 19 recordsLinked to original sources

PubMed2024

A PLA2 deletion mutant using CRISPR/Cas9 coupled to RNASeq reveals insect immune genes associated with eicosanoid signaling.

Eicosanoids mediate insect immune responses and synthesized by the catalytic activity of phospholipase A2 (PLA2). A uniquely encoded secretory PLA2 (sPLA2) is associated with immune responses of a lepidopteran insect, Spodoptera exigua. Its deletion mutant was generated using a CRISPR/Cas9 genome editing technology. Both wild and mutant lines were then immune-challenged, and the resulting transcripts were compared with their naïve transcripts by RNASeq using the Illumina-HiSeq platform. In total, 12,878 unigenes were further analyzed by differentially expressed gene tools. Over 69% of the expressed genes in S. exigua larvae are modulated in their expression levels by eicosanoids, recorded from CRISPR/Cas9 mutagenesis against an eicosanoid-synthetic gene, Se-sPLA2. Further, about 36% of the immune-associated genes are controlled by the eicosanoids in S. exigua. Indeed, the deletion mutant suffered significant immunosuppression in both cellular and humoral responses in response to bacterial challenge as well as severely reduced developmental and reproductive potentials.

Animals↗

PubMed2024

RNAseq analysis of oocyte maturation from the germinal vesicle stage to metaphase II in pig and human.

During maturation oocytes at the germinal vesicle (GV) stage progress to metaphase II (MII). However, during in vitro maturation a proportion often fail to progress. To understand these processes, we employed RNA sequencing to examine the transcriptome profile of these three groups of oocytes from the pig. We compared our findings with similar public oocyte data from humans. The transcriptomes in oocytes that failed to progress was similar to those that did. We found in both species, the most upregulated genes in MII oocytes were associated with chromosome segregation and cell cycle processes, while the most down regulated genes were relevant to ribosomal and mitochondrial pathways. Moreover, those genes involved in chromosome segregation during GV to MII transition were conserved in pig and human. We also compared MII and GV oocyte transcriptomes at the isoform transcript level in both species. Several thousands of genes (including DTNBP1, MAPK1, RAB35, GOLGA7, ATP1A1 and ATP2B1) identified as not different in expression at a gene transcript level were found to have differences in isoform transcript levels. Many of these genes were involved in ATPase-dependent or GTPase-dependent intracellular transport in pig and human, respectively. In conclusion, our study suggests the failure to progress to MII in vitro may not be regulated at the level of the genome and that many genes are differentially regulated at the isoform level, particular those involved ATPase- or GTPase-dependent intracellular transport.

Humans↗

PubMed2026

Precision-Based Filtering Facilitates Cross-Referencing of Conventional and Single-Nucleus Transcriptomes to Identify Time- and Temperature-Sensitive Cell Populations.

Transcriptome analysis via RNA sequencing (RNAseq) has become a ubiquitous method of molecular characterization from whole organisms, dissected tissues, and single cells. These experiments continue to provide an extraordinary volume of data describing molecular states and responses to many conditions. However, standard approaches to RNAseq analysis commonly use expression level filters that eliminate potentially useful data in the service of decreasing noise. Here we describe the implementation of a coefficient of variation-based filter for RNAseq gene expression data. This filter prioritizes consistent data across replicates, allowing lowly-expressed genes with low-variation measurements to be retained for downstream analysis. We show, using two independent Arabidopsis RNAseq datasets, that this filter allows for the inclusion of many more transcription factors than even a low-stringency expression level filter. This effect is independent of sequencing depth. We find that these lowly-expressed genes mark specific cell clusters in our single-nucleus (sn)RNAseq dataset and may facilitate future characterization of currently unknown cell types or states. We further characterize communities of co-expressed genes, sampled across the day at two growth temperatures, in relation to snRNAseq cell clusters, finding evidence for a highly photosynthetic cell population, and a cell state marked by high cell division and translation. These methods can be expanded to RNAseq analysis in many systems, facilitating the construction of more detailed models of tissue-specific gene regulatory networks.

Transcriptome analysis↗

PubMed2025

RNA sequencing offers new diagnostic opportunities in neurodevelopmental disorders: A systematic review.

PURPOSE: Transcriptomics by way of RNA sequencing (RNAseq) has emerged as a means to increase the diagnostic yield in genetic conditions. In this systematic review, we focus on the contribution of transcriptomics to improve the diagnostic yield in neurodevelopmental disorders. METHODS: We performed a systematic literature search in PubMed until January 2024, including articles describing diagnostic RNAseq on at least 1 individual with a primary neurodevelopmental phenotype. We extracted data on cohort size, phenotype, sample tissue, previously used diagnostic methods, added diagnostic yield of RNAseq, the use of control samples, and technical aspects of the RNA sequencing methodology. RESULTS: A total of 17 articles were eligible for inclusion in the systematic review. We found an average added diagnostic yield of 15.5% through RNA sequencing for individuals with neurodevelopmental disorders. There is heterogeneity in the tissue type, reported quality measures, and the computational pipeline. CONCLUSION: The significantly increased diagnostic yield demonstrates the value of this novel tool in the diagnostic setting of neurodevelopmental disorders. Our results offer an overview of common methodologies for RNAseq and allow us to formulate recommendations for genetic labs and clinicians when implementing RNAseq as a diagnostic tool. Lastly, we provide recommendations for future publications to increase transparency and reproducibility.

Humans↗

PubMed2025

Integrated Multiomics Analyses of the Molecular Landscape of Sarcopenia in Alcohol-Related Liver Disease.

BACKGROUND: Skeletal muscle is a major target for ethanol-induced perturbations, leading to sarcopenia in alcohol-related liver disease (ALD). The complex interactions and pathways involved in adaptive and maladaptive responses to ethanol in skeletal muscle are not well understood. Unlike hypothesis-driven experiments, an integrated multiomics-experimental validation approach provides a comprehensive view of these interactions. METHODS: We performed multiomics analyses with experimental validation to identify novel regulatory mechanisms of sarcopenia in ALD. Studies were done in a comprehensive array of models including ethanol-treated (ET) murine and human-induced pluripotent stem cell-derived myotubes (hiPSCm), skeletal muscle from a mouse model of ALD (mALD) and human patients with alcohol-related cirrhosis and controls. We generated 13 untargeted datasets, including chromatin accessibility (assay for transposase accessible chromatin), RNA sequencing, proteomics, phosphoproteomics, acetylomics and metabolomics, and conducted integrated multiomics analyses using UpSet plots and feature extraction. Key findings were validated using immunoblots, redox measurements (NAD+/NADH ratio), imaging and senescence-associated molecular phenotype (SAMP) assays. Mechanistic studies included mitochondrial-targeted Lactobacillus brevis NADH oxidase (MitoLbNOX) to increase redox ratio and MitoTempo as a mitochondrial free radical scavenger. RESULTS: Multiomics analyses revealed enrichment in mitochondrial oxidative function, protein synthesis and senescence pathways consistent with the known effects of hypoxia-inducible factor 1α (HIF1α) during normoxia. Across preclinical and clinical models, HIF1α targets (n = 32 genes) and signalling genes (n > 100 genes) (n = 3 ATACseq, n = 65 phosphoproteomics, n = 10 acetylomics, n = 6 C2C12 proteomics, n = 106 C2C12 RNAseq, n = 64 hiPSC RNAseq, n = 30 hiPSC proteomics, n = 3 mouse proteomics, n = 25 mouse RNAseq, n = 8 human RNAseq, n = 3 human proteomics) were increased. Stabilization of HIF1α (C2C12, 6hEtOH 0.24 ± 0.09; p = 0.043; mALD 0.32 ± 0.074; p = 0.005; data shown as mean difference ± standard error mean) was accompanied by enrichment in the early transient and late change clusters, -log(p-value) = 1.5-3.8, of the HIF1α signalling pathway. Redox ratio was reduced in ET myotubes (C2C12: 15512 ± 872.1, p < 0.001) and mALD muscle, with decreased expression of electron transport chain components (CI-V, p < 0.05) and Sirt3 (C2C12: 0.067 ± 0.023, p = 0.025; mALD: 0.41 ± 0.12, p = 0.013). Acetylation of mitochondrial proteins was increased in both models (C2C12: 107364 ± 4558, p = 0.03; mALD: 40036 ± 18 987, p = 0.049). Ethanol-induced SAMP was observed across models (P16: C2C12: 0.2845 ± 0.1145, p < 0.05; hiPSCm: 0.2591, p = 0.041). MitoLbNOX treatment reversed redox imbalance, HIF1α stabilization, global acetylation and myostatin expression (p < 0.05). CONCLUSIONS: An integrated multiomics approach, combined with experimental validation, identifies HIF1α stabilization and accelerated post-mitotic senescence as novel mechanisms of sarcopenia in ALD. These findings show the complex molecular interactions leading to mitochondrial dysfunction and progressive sarcopenia in ALD.

Sarcopenia↗

PubMed2024

A Transcriptional Signature of Induced Neurons Differentiates Virologically Suppressed People Living With HIV from People Without HIV.

Neurocognitive impairment is a prevalent and important co-morbidity in virologically suppressed people living with HIV (PLWH), yet the underlying mechanisms remain elusive and treatments lacking. Here, we explored for the first time, use of participant-derived directly induced neurons (iNs) to model neuronal biology and injury in PLWH. iNs retain age- and disease-related features of the donors, providing unique opportunities to reveal novel aspects of neurological disorders. We obtained primary dermal fibroblasts from six virologically suppressed PLWH (range: 27 - 64 years, median: 53); 83% Male; 50% White) and seven matched people without HIV (PWOH) (range: 27 - 66, median: 55); 71% Male; 57% White). iNs were generated using transcription factors NGN2 and ASCL1, and validated by immunocytochemistry and single-cell-RNAseq. Transcriptomic analysis using bulk-RNAseq identified 29 significantly differentially expressed genes between iNs from PLWH and PWOH. Of these, 16 genes were downregulated and 13 upregulated in PLWH iNs. Protein-protein interaction network mapping indicates that iNs from PLWH exhibit differences in extracellular matrix organization and synaptic transmission. IFI27 was upregulated in iNs from PLWH, which complements independent post-mortem studies demonstrating elevated IFI27 expression in PLWH-derived brain tissue, indicating that iN generation reconstitutes this pathway. Finally, we observed that expression of the FOXL2NB-FOXL2-LINC01391 genome locus is reduced in iNs from PLWH and negatively correlates with neurocognitive impairment. Thus, we have identified an iN gene signature of HIV through direct reprogramming of skin fibroblasts into neurons revealing novel mechanisms of neurocognitive impairment in PLWH.

HAND↗

PubMed2026

Unraveling a Diagnostic Enigma: A TECPR2 Case Solved Through Multi-Omic Genomics.

TECPR2 is a key regulator of autophagy, encoded by the TECPR2 gene. Pathogenic variants in this gene have been linked to a rare hereditary sensory and autonomic neuropathy with intellectual disability (HSAN9). We report a teenage female with a syndromic intellectual disability disorder associated with neuromuscular abnormalities. Multi-omics analysis including genomics, transcriptomics, and proteomics, together with muscle biopsy from the affected individual, were used in this clinical case. Through trio exome sequencing we identified two heterozygous variants in the TECPR2 gene, NM_014844.4: c.480G>A; p.(Gln160=) and c.2846C>A; p.(Ala949Glu). Both were classified as variants of uncertain significance due to the lack of supporting evidence for pathogenicity. Subsequent long-read sequencing phased the variants and confirmed they were in trans. Additional functional studies using RNAseq and proteomics analyses verified the pathogenicity of the variants. This case study demonstrated the value of a multi-omics assisted analysis, which complemented the traditional phenotype-first approach in reaching a definitive clinical diagnosis.

Humans↗

PubMed2017

Nonpromoter methylation of the CDKN2A gene with active transcription is associated with improved locoregional control in laryngeal squamous cell carcinoma.

We previously reported a novel association between CDKN2A nonpromoter methylation and transcription (ARF/INK4a) in human papillomavirus associated oropharyngeal tumors. In this study we assessed whether nonpromoter CDKN2A methylation in laryngeal squamous cell carcinomas (LXSCC) conferred a similar association with transcription that predicted patient outcome. We compared DNA methylation and ARF/INK4a RNA expression levels for the CDKN2A locus using the Illumina HumanMethylation27 beadchip and RT-PCR in 43 LXSCC tumor samples collected from a prospective study of head and neck cancer patients treated at Montefiore Medical Center (MMC). Validation was performed using RNAseq data on 111 LXSCC tumor samples from the Cancer Genome Atlas (TCGA). The clinical relevance of combined nonpromoter CDKN2A methylation and transcription was assessed by multivariate Cox regression for locoregional recurrence on a subset of 69 LXSCC patients with complete clinicopathologic data from the MMC and TCGA cohorts. We found evidence of CDKN2A nonpromoter hypermethylation in a third of LXSCC from our MMC cohort, which was significantly associated with increased ARF and INK4a RNA expression (Wilcoxon rank-sum, P = 0.007 and 0.003, respectively). A similar association was confirmed in TCGA samples (Wilcoxon rank-sum test P < 0.0001 for ARF and INK4a). Patients with CDKN2A hypermethylation or high ARF/INK4a expression were significantly less likely to develop a locoregional recurrence compared to those with neither of the features, independent of other clinicopatholgic risk factors (adjusted hazard ratio=0.21, 95% confidence interval:0.05-0.81). These results support the conclusion that CDKN2A nonpromoter methylation is associated with increased ARF and INK4a RNA expression, and improved locoregional control in LXSCC.

Carcinoma, Squamous Cell↗

PubMed2025

Transposable elements are dysregulated in brains of individuals with major depressive disorder.

Transposable elements (TEs) are repetitive DNA sequences capable of being transcribed and re-integrated, or transposed, into distinct loci throughout the genome. While thought to be largely transcriptionally silenced in brain, TE transcription is increasingly recognized as dynamic and involved in human health and disease states, including in disorders of the brain. In this study, we annotated TE transcripts in publicly available RNA sequencing (RNAseq) of postmortem human brain tissue to investigate the expression profile of TE transcripts in individuals with major depressive disorder (MDD) compared to healthy controls. Our findings reveal a robust impact to TE transcript expression in the brains of subjects with MDD relative to controls. This work points to the aberrant transcription of cortical TEs as a potentially overlooked molecular signature of MDD.

Humans↗

PubMed2026

Genomic and transcriptomic quality control for an autologous iPSC-derived cell therapy for Parkinson's disease.

Toward development of an autologous, induced pluripotent stem cell (iPSC)-based cell therapy for Parkinson's disease (PD), we demonstrate successful, reproducible genomic and transcriptomic qualification of patient-derived dopaminergic neuron precursor cells (DANPCs) across multiple donors. Our analysis includes whole-genome sequencing data from fibroblasts, iPSCs, and DANPCs and the development of NeuriTest, an RNAseq-based bioinformatic analysis of DANPCs designed to predict cell quality based on empirical animal data. Autologous cell therapies are immune matched to the patient, potentially augmenting durability of benefit compared to allogeneic cells while negating the need for immunosuppression and accompanying side effects. Patient-specific iPSCs are an autologous cell source that can be differentiated to dopaminergic neurons, the cell type lost in PD. We report here our preclinical manufacturing strategy and results demonstrating efficacy in a PD rodent model and safety in a 9-month GLP toxicology study.

Parkinson’s disease↗

PubMed2025

Integrated genomic analysis of NF1-associated peripheral nerve sheath tumors: an updated biorepository dataset.

Neurofibromatosis type 1 (NF1) is an inherited neurocutaneous condition that predisposes to the development of peripheral nerve sheath tumors (PNST) including cutaneous neurofibromas (CNF), plexiform neurofibromas (PNF), atypical neurofibromatous neoplasms of uncertain biologic potential (ANNUBP), and malignant peripheral nerve sheath tumors (MPNST). The Johns Hopkins NF1 biospecimen repository promotes the successful advancement of therapeutic developments for NF1-associated PNST through acquisition and genomic analysis of human tumor specimens. RNA sequencing (RNAseq) and whole exome sequencing (WES) data were generated from 73 and 114 primary human tumor samples, respectively. These pre-processed data, standardized for immediate computational analysis, are accessible through the NF Data Portal, allowing immediate interrogation. This dataset combines new and previously released samples, offering a comprehensive view of the entire cohort sequenced. As a dedicated effort to systematically bank tumor samples from people with NF1, in collaboration with molecular geneticists and computational biologists, the Johns Hopkins NF1 biospecimen repository offers access to tissue samples and genomic data to promote the advancement of NF1-related tumor biologic insights and therapies.

Humans↗

PubMed2024

Decorin suppresses tumor lymphangiogenesis: A mechanism to curtail cancer progression.

The complex interplay between malignant cells and the cellular and molecular components of the tumor stroma is a key aspect of cancer growth and development. These tumor-host interactions are often affected by soluble bioactive molecules such as proteoglycans. Decorin, an archetypical small leucine-rich proteoglycan primarily expressed by stromal cells, affects cancer growth in its soluble form by interacting with several receptor tyrosine kinases (RTK). Overall, decorin leads to a context-dependent and protracted cessation of oncogenic RTK activity by attenuating their ability to drive a prosurvival program and to sustain a proangiogenic network. Through an unbiased transcriptomic analysis using deep RNAseq, we identified that decorin down-regulated a cluster of tumor-associated genes involved in lymphatic vessel (LV) development when systemically delivered to mice harboring breast carcinoma allografts. We found that Lyve1 and Podoplanin, two established markers of LVs, were markedly suppressed at both the mRNA and protein levels, and this suppression correlated with a significant reduction in tumor LVs. We further identified that soluble decorin, but not its homologous proteoglycan biglycan, inhibited LV sprouting in an ex vivo 3D model of lymphangiogenesis. Mechanistically, we found that decorin interacted with vascular endothelial growth factor receptor 3 (VEGFR3), the main lymphatic RTK, and its activity was required for the decorin-mediated block of lymphangiogenesis. Finally, we identified that Lyve1 was in part degraded via decorin-evoked autophagy in a nutrient- and energy-independent manner. These findings implicate decorin as a biological factor with antilymphangiogenic activity and provide a potential therapeutic agent for curtailing breast cancer growth and metastasis.

Decorin↗

PubMed2025

scACCorDiON: a clustering approach for explainable patient level cell-cell communication graph analysis.

MOTIVATION: Combining single-cell sequencing with ligand-receptor (LR) analysis paves the way for the characterization of cell communication events in complex tissues. In particular, directed weighted graphs naturally represent cell-cell communication events. However, current computational methods cannot yet analyze sample-specific cell-cell communication events, as measured in single-cell data produced in large patient cohorts. Cohort-based cell-cell communication analysis presents many challenges, such as the nonlinear nature of cell-cell communication and the high variability given by the patient-specific single-cell RNAseq datasets. RESULTS: Here, we present scACCorDiON (single-cell Analysis of Cell-Cell Communication in Disease clusters using Optimal transport in Directed Networks), an optimal transport algorithm exploring node distances on the Markov Chain as the ground metric between directed weighted graphs. Benchmarking indicates that scACCorDiON performs a better clustering of samples according to their disease status than competing methods that use undirected graphs. We provide a case study of pancreas adenocarcinoma, where scACCorDion detects a sub-cluster of disease samples associated with changes in the tumor microenvironment. Our study case corroborates that clusters provide a robust and explainable representation of cell-cell communication events and that the expression of detected LR pairs is predictive of pancreatic cancer survival. AVAILABILITY AND IMPLEMENTATION: The code of scACCorDiON is available at https://scaccordion.readthedocs.io/en/latest/. and https://doi.org/10.5281/zenodo.15267648. The survival analysis package can be found at https://github.com/CostaLab/scACCorDiON.su.

Humans↗

PubMed2026

De Novo Genome Sequence Assembly of the Algal Endosymbiont Micractinium conductrix Derived From Its Host Paramecium bursaria 186b.

Endosymbiosis is a major driver of evolutionary innovation and underpins the function of diverse ecosystems. The origins and evolution of endosymbiosis are challenging to study experimentally due to the short-lived culturability of many microbial strains derived from endosymbiotic interactions. The facultative endosymbiosis between the ciliate, Paramecium bursaria, and the green alga, Micractinium conductrix (Chlorellaceae, Trebouxiophyceae), is ecologically widespread and has emerged as a powerful lab-tractable model system. This endosymbiosis is founded upon a reciprocal nutrient exchange, but each of the species can be cultured independently enabling quantification of symbiotic fitness effects, new partnerships to be generated in the lab, and co-associations to be subject to experimental evolution. To date, evolve-and-resequence approaches have been limited due to a lack of high-quality genome assemblies enabling gene variants to be identified. Here, we report a near telomere-to-telomere genome assembly for M. conductrix 186b, using a range of sequencing technologies. Comparative analysis shows that this is one of the most complete Chlorellaceae algal genome assemblies available to date. To aid accurate gene calling and annotation, we conducted both RNAseq and Iso-Seq transcriptome sequencing experiments. Collectively, these 'omics datasets will facilitate: (i) comparative genomics studies of endosymbiont evolution, (ii) evolve-and-resequence experiments, (iii) genome-scale metabolic modeling studies, and (iv) identification of targets for genetic modification experiments and biotechnological applications.

Symbiosis↗

PubMed2026

Chromosome-Level Reference Genome of the Desert Night Lizard Xantusia vigilis.

We present a reference-quality genome assembly for the desert night lizard (Xantusia vigilis). The night lizards (Xantusiidae) are a family of small-bodied lizards found in North America (Xantusia), Central America (Lepidophyma), and Cuba (Cricosaura). The night lizard family has an independent evolutionary history of at least 80 million years from its sister taxa within Scincoidea. The Xantusiids have several unique ecological, behavioral and evolutionary characteristics. For instance, the family contains the only squamate species that form diploid, unisexual, parthenogenic lineages. In addition, most night lizards are viviparous and form stable kin groups that are maintained over multiple years, an unusual life history strategy among lizards. Combining PacBio long-read sequencing, Hi-C, and RNAseq data we developed a reference-quality genome for the desert night lizard, X. vigilis. We assembled a complete mitochondrion and ~ 2.2 Gb nuclear genome, with 20 scaffolds that correlate in size to the X. vigilis karyotype. In addition, we found that X. vigilis chromosome 1 aligns with gene content of both of macrochromosome 1 and microchromosome 9 from a genome assembly of a species in the sister family Cordylidae (Hemicordylus capensis).

Xantusia↗

PubMed2026

Reference genome of the Californian trapdoor spider Aptostichus stephencolberti Bond 2008 (Araneae: Mygalomorphae: Euctenizidae).

We present a reference genome assembly for the trapdoor spider Aptostichus stephencolberti. This species, described in 2008, is endemic to the highly fragmented coastal dune habitats of Northern California from Monterey to the San Francisco Bay Area. Trapdoor spiders are ideal taxa for landscape scale genomic studies owing to their extreme site fidelity and limited dispersal capabilities; these same characteristics make them prone to extinction. Genomic studies of species like A. stephencolberti can reveal novel areas of endemism and high conservation value that may not be evident in species with wider ranges and greater dispersal capabilities. As part of the California Conservation Genomics Project, we constructed the A. stephencolberti reference genome from high quality long-read sequences, scaffolded with proximity ligation Omni-C data. The primary assembly comprises 551 scaffolds spanning 3.63 Gbp, a scaffold N50 of 62.2 Mbp and BUSCO completeness of 95.6%. We estimate 52 chromosomes yet find no (TTAGG)n telomer repeats. Expanding the telomeric repeat search finds an ancestral loss of the repeat from all spiders. Automated annotation using the NCBI refseq pipeline and RNAseq data from whole adults finds 14,067 genes with a BUSCO annotation completeness of 95.56%. Repeat annotation identified 77% of the genome to be interspersed repeats. This resource, the first for family Euctenizidae will facilitate future study and resulting conservation actions of A. stephencolberti and other Aptostichus sp. populations associated with the rapidly changing California coastal dune ecosystem.

Aptostichus stephencolberti↗

PubMed2026

Mutation accumulation in a hybrid parthenogenetic vertebrate.

Asexual lineages are thought to experience elevated extinction rates compared with sexual species, yet direct evidence for the underlying genetic causes remains scarce. Muller's ratchet predicts that the absence of recombination in asexual organisms facilitates the accumulation of deleterious mutations, thereby reducing long-term fitness. Here, we test this hypothesis in the hybrid-origin, parthenogenetic whiptail lizard Aspidoscelis tesselatus by integrating short-read RNAseq and long-read IsoSeq data from both the asexual lineage and its parental sexual species. We reconstructed phased transcripts for A. tesselatus to quantify mutation accumulation relative to the parental sexual species. Comparative analyses revealed elevated ω ratios in both parental genomic complements (subgenomes) of the parthenogenetic lineage, consistent with accelerated accumulation of nonsynonymous mutations. Structural variant analyses identified multiple indels in expressed transcripts predicted to disrupt protein domains. Functional annotation indicated that genes affected by both single-nucleotide variants and indels were enriched for roles in chromatin organization, apoptosis regulation, and transcriptional control. While both parental subgenomes showed similar evolutionary patterns, the maternal complement exhibited more structural and missense mutations than the paternal complement. Together, these results provide evidence that mutations accumulate in asexual A. tesselatus in genes involved in core cellular functions, supporting theoretical predictions that Muller's ratchet contributes to mutation accumulation in asexual lineages.

Animals↗

PubMed2023

GBMdeconvoluteR accurately infers proportions of neoplastic and immune cell populations from bulk glioblastoma transcriptomics data.

BACKGROUND: Characterizing and quantifying cell types within glioblastoma (GBM) tumors at scale will facilitate a better understanding of the association between the cellular landscape and tumor phenotypes or clinical correlates. We aimed to develop a tool that deconvolutes immune and neoplastic cells within the GBM tumor microenvironment from bulk RNA sequencing data. METHODS: We developed an IDH wild-type (IDHwt) GBM-specific single immune cell reference consisting of B cells, T-cells, NK-cells, microglia, tumor associated macrophages, monocytes, mast and DC cells. We used this alongside an existing neoplastic single cell-type reference for astrocyte-like, oligodendrocyte- and neuronal progenitor-like and mesenchymal GBM cancer cells to create both marker and gene signature matrix-based deconvolution tools. We applied single-cell resolution imaging mass cytometry (IMC) to ten IDHwt GBM samples, five paired primary and recurrent tumors, to determine which deconvolution approach performed best. RESULTS: Marker-based deconvolution using GBM-tissue specific markers was most accurate for both immune cells and cancer cells, so we packaged this approach as GBMdeconvoluteR. We applied GBMdeconvoluteR to bulk GBM RNAseq data from The Cancer Genome Atlas and recapitulated recent findings from multi-omics single cell studies with regards associations between mesenchymal GBM cancer cells and both lymphoid and myeloid cells. Furthermore, we expanded upon this to show that these associations are stronger in patients with worse prognosis. CONCLUSIONS: GBMdeconvoluteR accurately quantifies immune and neoplastic cell proportions in IDHwt GBM bulk RNA sequencing data and is accessible here: https://gbmdeconvoluter.leeds.ac.uk.

Humans↗