recount3: summaries and queries for large-scale RNA-seq expression and splicing

Co-expression in tissue-specific gene networks links genes in cancer-susceptibility loci to known somatic driver genes

BACKGROUND

The genetic background of cancer remains complex and challenging to integrate. Many somatic mutations within genes are known to cause and drive cancer, while genome-wide association studies (GWAS) of cancer have revealed many germline risk factors associated with cancer. However, the overlap between known somatic driver genes and positional candidate genes from GWAS loci is surprisingly small. We hypothesised that genes from multiple independent cancer GWAS loci should show tissue-specific co-regulation patterns that converge on cancer-specific driver genes.

RESULTS

We studied recent well-powered GWAS of breast, prostate, colorectal and skin cancer by estimating co-expression between genes and subsequently prioritising genes that show significant co-expression with genes mapping within susceptibility loci from cancer GWAS. We observed that the prioritised genes were strongly enriched for cancer drivers defined by COSMIC, IntOGen and Dietlein et al. The enrichment of known cancer driver genes was most significant when using co-expression networks derived from non-cancer samples of the relevant tissue of origin.

CONCLUSION

We show how genes within risk loci identified by cancer GWAS can be linked to known cancer driver genes through tissue-specific co-expression networks. This provides an important explanation for why seemingly unrelated sets of genes that harbour either germline risk factors or somatic mutations can eventually cause the same type of disease.

Aging-associated Alterations in the Gene Regulatory Network Landscape Associate with Risk, Prognosis and Response to Therapy in Lung Adenocarcinoma

Aging is the primary risk factor for many individual cancer types, including lung adenocarcinoma (LUAD). To understand how aging-related alterations in the regulation of key cellular processes might affect LUAD risk and survival outcomes, we built individual (person)-specific gene regulatory networks integrating gene expression, transcription factor protein-protein interaction, and sequence motif data, using PANDA/LIONESS algorithms, for both non-cancerous lung tissue samples from the Genotype Tissue Expression (GTEx) project and LUAD samples from The Cancer Genome Atlas (TCGA). In GTEx, we found that pathways involved in cell proliferation and immune response are increasingly targeted by regulatory transcription factors with age; these aging-associated alterations are accelerated by tobacco smoking and resemble oncogenic shifts in the regulatory landscape observed in LUAD and suggests that dysregulation of aging pathways might be associated with an increased risk of LUAD. Comparing normal adjacent samples from individuals with LUAD with healthy lung tissue samples from those without LUAD, we found that aging-associated genes show greater aging-biased targeting patterns in younger individuals with LUAD compared to their healthy counterparts of similar age, a pattern suggestive of age acceleration. This implies that an accelerated aging process may be responsible for tumor incidence in younger individuals. Using drug repurposing tool CLUEreg, we found small molecule drugs with potential geroprotective effects that may alter the accelerating aging profiles we found. We also observed that, in contrast to chronological age, a network-informed aging signature was associated with survival and response to chemotherapy in LUAD.

Translatome profiling reveals Itih4 as a novel smooth muscle cell-specific gene in atherosclerosis

AIMS

Vascular smooth muscle cells (SMCs) and their derivatives are key contributors to the development of atherosclerosis. However, studying changes in SMC gene expression in heterogeneous vascular tissues is challenging due to the technical limitations and high cost associated with current approaches. In this paper, we apply translating ribosome affinity purification sequencing to profile SMC-specific gene expression directly from tissue.

METHODS AND RESULTS

To facilitate SMC-specific translatome analysis, we generated SMCTRAP mice, a transgenic mouse line expressing enhanced green fluorescent protein (EGFP)-tagged ribosomal protein L10a (EGFP-L10a) under the control of the SMC-specific αSMA promoter. These mice were further crossed with the atherosclerosis model Ldlr-/-, ApoB100/100 to generate SMCTRAP-AS mice and used to profile atherosclerosis-associated SMCs in thoracic aorta samples of 15-month-old SMCTRAP and SMCTRAP-AS mice. Our analysis of SMCTRAP-AS mice showed that EGFP-L10a expression was localized to SMCs in various tissues, including the aortic wall and plaque. The TRAP fraction demonstrated high enrichment of known SMC-specific genes, confirming the specificity of our approach. We identified several genes, including Cemip, Lum, Mfge8, Spp1, and Serpina3, which are known to be involved in atherosclerosis-induced gene expression. Moreover, we identified several novel genes not previously linked to SMCs in atherosclerosis, such as Anxa4, Cd276, inter-alpha-trypsin inhibitor-4 (Itih4), Myof, Pcdh11x, Rab31, Serpinb6b, Slc35e4, Slc8a3, and Spink5. Among them, we confirmed the SMC-specific expression of Itih4 in atherosclerotic lesions using immunofluorescence staining of mouse aortic roots and spatial transcriptomics of human carotid arteries. Furthermore, our more detailed analysis of Itih4 showed its link to coronary artery disease through the colocalization of genome-wide association studies, splice quantitative trait loci (QTL), and protein QTL signals.

CONCLUSION

We generated a SMC-specific TRAP mouse line to study atherosclerosis and identified Itih4 as a novel SMC-expressed gene in atherosclerotic plaques, warranting further investigation of its putative function in extracellular matrix stability and genetic evidence of causality.

Systematic assessment of long-read RNA-seq methods for transcript identification and quantification

The Long-read RNA-Seq Genome Annotation Assessment Project Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. Using different protocols and sequencing platforms, the consortium generated over 427 million long-read sequences from complementary DNA and direct RNA datasets, encompassing human, mouse and manatee species. Developers utilized these data to address challenges in transcript isoform detection, quantification and de novo transcript detection. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. Incorporating additional orthogonal data and replicate samples is advised when aiming to detect rare and novel transcripts or using reference-free approaches. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis.

The annotation of GBA1 has been concealed by its protein-coding pseudogene GBAP1

Mutations in GBA1 cause Gaucher disease and are the most important genetic risk factor for Parkinson's disease. However, analysis of transcription at this locus is complicated by its highly homologous pseudogene, GBAP1. We show that >50% of short RNA-sequencing reads mapping to GBA1 also map to GBAP1. Thus, we used long-read RNA sequencing in the human brain, which allowed us to accurately quantify expression from both GBA1 and GBAP1. We discovered significant differences in expression compared to short-read data and identify currently unannotated transcripts of both GBA1 and GBAP1. These included protein-coding transcripts from both genes that were translated in human brain, but without the known lysosomal function-yet accounting for almost a third of transcription. Analyzing brain-specific cell types using long-read and single-nucleus RNA sequencing revealed region-specific variations in transcript expression. Overall, these findings suggest nonlysosomal roles for GBA1 and GBAP1 with implications for our understanding of the role of GBA1 in health and disease.

node2vec2rank: Large Scale and Stable Graph Differential Analysis via Multi-Layer Node Embeddings and Ranking

Computational methods in biology can infer large molecular interaction networks from multiple data sources and at different resolutions, creating unprecedented opportunities to explore the mechanisms driving complex biological phenomena. Networks can be built to represent distinct conditions and compared to uncover graph-level differences-such as when comparing patterns of gene-gene interactions that change between biological states. Given the importance of the graph comparison problem, there is a clear and growing need for robust and scalable methods that can identify meaningful differences. We introduce node2vec2rank (n2v2r), a method for graph differential analysis that ranks nodes according to the disparities of their representations in joint latent embedding spaces. Improving upon previous bag-of-features approaches, we take advantage of recent advances in machine learning and statistics to compare graphs in higher-order structures and in a data-driven manner. Formulated as a multi-layer spectral embedding algorithm, n2v2r is computationally efficient, incorporates stability as a key feature, and can provably identify the correct ranking of differences between graphs in an overall procedure that adheres to veridical data science principles. By better adapting to the data, node2vec2rank clearly outperformed the commonly used node degree in finding complex differences in simulated data. In the real-world applications of breast cancer subtype characterization, analysis of cell cycle in single-cell data, and searching for sex differences in lung adenocarcinoma, node2vec2rank found meaningful biological differences enabling the hypothesis generation for therapeutic candidates. Software and analysis pipelines implementing n2v2r and used for the analyses presented here are publicly available.

Elevated nuclear TDP-43 induces constitutive exon skipping

BACKGROUND

Cytoplasmic inclusions and loss of nuclear TDP-43 are key pathological features found in several neurodegenerative disorders, suggesting both gain- and loss-of-function mechanisms of disease. To study gain-of-function, TDP-43 overexpression has been used to generate in vitro and in vivo model systems.

METHODS

We analyzed RNA-seq datasets from mouse and human neurons overexpressing TDP-43 to explore species specific splicing patterns. We explored the dynamics between TDP-43 levels and exon repression in vitro. Furthermore we analyzed human brain samples and publicly available RNA datasets to explore the relationship between exon repression and disease.

RESULTS

Our study shows that excessive levels of nuclear TDP-43 protein lead to constitutive exon skipping that is largely species-specific. Furthermore, while aberrant exon skipping is detected in some human brains, it is not correlated with disease, unlike the incorporation of cryptic exons that occurs after loss of TDP-43.

CONCLUSIONS

Our findings emphasize the need for caution in interpreting TDP-43 overexpression data and stress the importance of controlling for exon skipping when generating models of TDP-43 proteinopathy.

Single-cell and spatial transcriptomics analysis of human adrenal aging

OBJECTIVE

The human adrenal cortex comprises three functionally and structurally distinct layers that produce layer-specific steroid hormones. With aging, the human adrenal cortex undergoes functional and structural alteration or "adrenal aging", leading to the unbalanced production of steroid hormones. Given the marked species differences in adrenal biology, the underlying mechanisms of human adrenal aging have not been sufficiently studied. This study was designed to elucidate the mechanisms linking the functional and structural alterations of the human adrenal cortex.

METHODS

We conducted single-cell RNA sequencing and spatial transcriptomics analysis of the aged human adrenal cortex.

RESULTS

The data of this study suggest that the layer-specific alterations of multiple signaling pathways underlie the abnormal layered structure and layer-specific changes in steroidogenic cells. We also highlighted that macrophages mediate age-related adrenocortical cell inflammation and senescence.

CONCLUSIONS

This study is the first detailed analysis of the aged human adrenal cortex at single-cell resolution and helps to elucidate the mechanism of human adrenal aging, thereby leading to a better understanding of the pathophysiology of age-related disorders associated with adrenal aging.

MousiPLIER: A Mouse Pathway-Level Information Extractor Model

High-throughput gene expression profiling measures individual gene expression across conditions. However, genes are regulated in complex networks, not as individual entities, limiting the interpretability of gene expression data. Machine learning models that incorporate prior biological knowledge are a powerful tool to extract meaningful biology from gene expression data. Pathway-level information extractor (PLIER) is an unsupervised machine learning method that defines biological pathways by leveraging the vast amount of published transcriptomic data. PLIER converts gene expression data into known pathway gene sets, termed latent variables (LVs), to substantially reduce data dimensionality and improve interpretability. In the current study, we trained the first mouse PLIER model on 190,111 mouse brain RNA-sequencing samples, the greatest amount of training data ever used by PLIER. We then validated the mousiPLIER approach in a study of microglia and astrocyte gene expression across mouse brain aging. mousiPLIER identified biological pathways that are significantly associated with aging, including one latent variable (LV41) corresponding to striatal signal. To gain further insight into the genes contained in LV41, we performed k-means clustering on the training data to identify studies that respond strongly to LV41. We found that the variable was relevant to striatum and aging across the scientific literature. Finally, we built a Web server (http://mousiplier.greenelab.com/) for users to easily explore the learned latent variables. Taken together, this study defines mousiPLIER as a method to uncover meaningful biological processes in mouse brain transcriptomic studies.

Flexible modeling of regulatory networks improves transcription factor activity estimation

Transcriptional regulation plays a crucial role in determining cell fate and disease, yet inferring the key regulators from gene expression data remains a significant challenge. Existing methods for estimating transcription factor (TF) activity often rely on static TF-gene interaction databases and cannot adapt to changes in regulatory mechanisms across different cell types and disease conditions. Here, we present a new algorithm - Transcriptional Inference using Gene Expression and Regulatory data (TIGER) - that overcomes these limitations by flexibly modeling activation and inhibition events, up-weighting essential edges, shrinking irrelevant edges towards zero through a sparse Bayesian prior, and simultaneously estimating both TF activity levels and changes in the underlying regulatory network. When applied to yeast and cancer TF knock-out datasets, TIGER outperforms comparable methods in terms of prediction accuracy. Moreover, our application of TIGER to tissue- and cell-type-specific RNA-seq data demonstrates its ability to uncover differences in regulatory mechanisms. Collectively, our findings highlight the utility of modeling context-specific regulation when inferring transcription factor activities.

Illuminating the function of the orphan transporter, SLC22A10, in humans and other primates

SLC22A10 is an orphan transporter with unknown substrates and function. The goal of this study is to elucidate its substrate specificity and functional characteristics. In contrast to orthologs from great apes, human SLC22A10, tagged with green fluorescent protein, is not expressed on the plasma membrane. Cells expressing great ape SLC22A10 orthologs exhibit significant accumulation of estradiol-17β-glucuronide, unlike those expressing human SLC22A10. Sequence alignments reveal a proline at position 220 in humans, which is a leucine in great apes. Replacing proline with leucine in SLC22A10-P220L restores plasma membrane localization and uptake function. Neanderthal and Denisovan genomes show proline at position 220, akin to modern humans, indicating functional loss during hominin evolution. Human SLC22A10 is a unitary pseudogene due to a fixed missense mutation, P220, while in great apes, its orthologs transport sex steroid conjugates. Characterizing SLC22A10 across species sheds light on its biological role, influencing organism development and steroid homeostasis.

Signature reversion of three disease-associated gene signatures prioritizes cancer drug repurposing candidates

Drug repurposing is promising because approving a drug for a new indication requires fewer resources than approving a new drug. Signature reversion detects drug perturbations most inversely related to the disease-associated gene signature to identify drugs that may reverse that signature. We assessed the performance and biological relevance of three approaches for constructing disease-associated gene signatures (i.e., limma, DESeq2, and MultiPLIER) and prioritized the resulting drug repurposing candidates for four low-survival human cancers. Our results were enriched for candidates that had been used in clinical trials or performed well in the PRISM drug screen. Additionally, we found that pamidronate and nimodipine, drugs predicted to be efficacious against the brain tumor glioblastoma (GBM), inhibited the growth of a GBM cell line and cells isolated from a patient-derived xenograft (PDX). Our results demonstrate that by applying multiple disease-associated gene signature methods, we prioritized several drug repurposing candidates for low-survival cancers.

Steroids-producing nodules: a two-layered adrenocortical nodular structure as a precursor lesion of cortisol-producing adenoma

BACKGROUND

The human adrenal cortex consists of three functionally and structurally distinct layers; zona glomerulosa, zona fasciculata (zF), and zona reticularis (zR), and produces adrenal steroid hormones in a layer-specific manner; aldosterone, cortisol, and adrenal androgens, respectively. Cortisol-producing adenomas (CPAs) occur mostly as a result of somatic mutations associated with the protein kinase A pathway. However, how CPAs develop after adrenocortical cells acquire genetic mutations, remains poorly understood.

METHODS

We conducted integrated approaches combining the detailed histopathologic studies with genetic, RNA-sequencing, and spatially resolved transcriptome (SRT) analyses for the adrenal cortices adjacent to human adrenocortical tumours.

FINDINGS

Histopathological analysis revealed an adrenocortical nodular structure that exhibits the two-layered zF- and zR-like structure. The nodular structures harbour GNAS somatic mutations, known as a driver mutation of CPAs, and confer cell proliferative and autonomous steroidogenic capacities, which we termed steroids-producing nodules (SPNs). RNA-sequencing coupled with SRT analysis suggests that the expansion of the zF-like structure contributes to the formation of CPAs, whereas the zR-like structure is characterised by a macrophage-mediated immune response.

INTERPRETATION

We postulate that CPAs arise from a precursor lesion, SPNs, where two distinct cell populations might contribute differently to adrenocortical tumorigenesis. Our data also provide clues to the molecular mechanisms underlying the layered structures of human adrenocortical tissues.

FUNDING

KAKENHI, The Uehara Memorial Foundation, Daiwa Securities Health Foundation, Kaibara Morikazu Medical Science Promotion Foundation, Secom Science and Technology Foundation, ONO Medical Research Foundation, and Japan Foundation for Applied Enzymology.

Dopamine signaling enriched striatal gene set predicts striatal dopamine synthesis and physiological activity in vivo

The polygenic architecture of schizophrenia implicates several molecular pathways involved in synaptic function. However, it is unclear how polygenic risk funnels through these pathways to translate into syndromic illness. Using tensor decomposition, we analyze gene co-expression in the caudate nucleus, hippocampus, and dorsolateral prefrontal cortex of post-mortem brain samples from 358 individuals. We identify a set of genes predominantly expressed in the caudate nucleus and associated with both clinical state and genetic risk for schizophrenia that shows dopaminergic selectivity. A higher polygenic risk score for schizophrenia parsed by this set of genes predicts greater dopamine synthesis in the striatum and greater striatal activation during reward anticipation. These results translate dopamine-linked genetic risk variation into in vivo neurochemical and hemodynamic phenotypes in the striatum that have long been implicated in the pathophysiology of schizophrenia.

Heterogenous nuclear ribonucleoprotein D-like controls endothelial cell functions

HnRNPs are ubiquitously expressed RNA-binding proteins, tightly controlling posttranscriptional gene regulation. Consequently, hnRNP networks are essential for cellular homeostasis and their dysregulation is associated with cancer and other diseases. However, the physiological function of hnRNPs in non-cancerous cell systems are poorly understood. We analyzed the importance of HNRNPDL in endothelial cell functions. Knockdown of HNRNPDL led to impaired proliferation, migration and sprouting of spheroids. Transcriptome analysis identified cyclin D1 (CCND1) and tropomyosin 4 (TPM4) as targets of HNRNPDL, reflecting the phenotypic changes after knockdown. Our findings underline the importance of HNRNPDL for the homeostasis of physiological processes in endothelial cells.

Tumor-agnostic transcriptome-based classifier identifies spatial infiltration patterns of CD8+T cells in the tumor microenvironment and predicts clinical outcome in early-phase and late-phase clinical trials

BACKGROUND

The immune status of a patient's tumor microenvironment (TME) may guide therapeutic interventions with cancer immunotherapy and help identify potential resistance mechanisms. Currently, patients' immune status is mostly classified based on CD8+tumor-infiltrating lymphocytes. An unmet need exists for comparable and reliable precision immunophenotyping tools that would facilitate clinical treatment-relevant decision-making and the understanding of how to overcome resistance mechanisms.

METHODS

We systematically analyzed the CD8 immunophenotype of 2023 patients from 14 phase I-III clinical trials using immunohistochemistry (IHC) and additionally profiled gene expression by RNA-sequencing (RNA-seq). CD8 immunophenotypes were classified by pathologists into CD8-desert, CD8-excluded or CD8-inflamed tumors using CD8 IHC staining in epithelial and stromal areas of the tumor. Using regularized logistic regression, we developed an RNA-seq-based classifier as a surrogate to the IHC-based spatial classification of CD8+tumor-infiltrating lymphocytes in the TME.

RESULTS

The CD8 immunophenotype and associated gene expression patterns varied across indications as well as across primary and metastatic lesions. Melanoma and kidney cancers were among the strongest inflamed indications, while CD8-desert phenotypes were most abundant in liver metastases across all tumor types. A good correspondence between the transcriptome and the IHC-based evaluation enabled us to develop a 92-gene classifier that accurately predicted the IHC-based CD8 immunophenotype in primary and metastatic samples (area under the curve inflamed=0.846; excluded=0.712; desert=0.855). The newly developed classifier was prognostic in The Cancer Genome Atlas (TCGA) data and predictive in lung cancer: patients with predicted CD8-inflamed tumors showed prolonged overall survival (OS) versus patients with CD8-desert tumors (HR 0.88; 95% CI 0.80 to 0.97) across TCGA, and longer OS on immune checkpoint inhibitor administration (phase III OAK study) in non-small-cell lung cancer (HR 0.75; 95% CI 0.58 to 0.97).

CONCLUSIONS

We provide a new precision immunophenotyping tool based on gene expression that reflects the spatial infiltration patterns of CD8+ lymphocytes in tumors. The classifier enables multiplex analyses and is easy to apply for retrospective, reverse translation approaches as well as for prospective patient enrichment to optimize the response to cancer immunotherapy.

Rummagene: massive mining of gene sets from supporting materials of biomedical research publications

Many biomedical research publications contain gene sets in their supporting tables, and these sets are currently not available for search and reuse. By crawling PubMed Central, the Rummagene server provides access to hundreds of thousands of such mammalian gene sets. So far, we scanned 5,448,589 articles to find 121,237 articles that contain 642,389 gene sets. These sets are served for enrichment analysis, free text, and table title search. Investigating statistical patterns within the Rummagene database, we demonstrate that Rummagene can be used for transcription factor and kinase enrichment analyses, and for gene function predictions. By combining gene set similarity with abstract similarity, Rummagene can find surprising relationships between biological processes, concepts, and named entities. Overall, Rummagene brings to surface the ability to search a massive collection of published biomedical datasets that are currently buried and inaccessible. The Rummagene web application is available at https://rummagene.com .

RummaGEO: Automatic Mining of Human and Mouse Gene Sets from GEO

The Gene Expression Omnibus (GEO) is a major open biomedical research repository for transcriptomics and other omics datasets. It currently contains millions of gene expression samples from tens of thousands of studies collected by many biomedical research laboratories from around the world. While users of the GEO repository can search the metadata describing studies for locating relevant datasets, there are currently no methods or resources that facilitate global search of GEO at the data level. To address this shortcoming, we developed RummaGEO, a webserver application that enables gene expression signature search of a large collection of human and mouse RNA-seq studies deposited into GEO. To develop the search engine, we performed offline automatic identification of sample conditions from the uniformly aligned GEO studies available from ARCHS4. We then computed differential expression signatures to extract gene sets from these studies. In total, RummaGEO currently contains 135,264 human and 158,062 mouse gene sets extracted from 23,395 GEO studies. Next, we analyzed the contents of the RummaGEO database to identify statistical patterns and perform various global analyses. The contents of the RummaGEO database are provided as a web-server search engine with signature search, PubMed search, and metadata search functionalities. Overall, RummaGEO provides an unprecedented resource for the biomedical research community enabling hypothesis generation for many future studies. The RummaGEO search engine is available from: https://rummageo.com/.

Benchmark of cellular deconvolution methods using a multi-assay reference dataset from postmortem human prefrontal cortex

Background

Cellular deconvolution of bulk RNA-sequencing (RNA-seq) data using single cell or nuclei RNA-seq (sc/snRNA-seq) reference data is an important strategy for estimating cell type composition in heterogeneous tissues, such as human brain. Computational methods for deconvolution have been developed and benchmarked against simulated data, pseudobulked sc/snRNA-seq data, or immunohistochemistry reference data. A major limitation in developing improved deconvolution algorithms has been the lack of integrated datasets with orthogonal measurements of gene expression and estimates of cell type proportions on the same tissue sample. Deconvolution algorithm performance has not yet been evaluated across different RNA extraction methods (cytosolic, nuclear, or whole cell RNA), different library preparation types (mRNA enrichment vs. ribosomal RNA depletion), or with matched single cell reference datasets.

Results

A rich multi-assay dataset was generated in postmortem human dorsolateral prefrontal cortex (DLPFC) from 22 tissue blocks. Assays included spatially-resolved transcriptomics, snRNA-seq, bulk RNA-seq (across six library/extraction RNA-seq combinations), and RNAScope/Immunofluorescence (RNAScope/IF) for six broad cell types. The Mean Ratio method, implemented in the DeconvoBuddies R package, was developed for selecting cell type marker genes. Six computational deconvolution algorithms were evaluated in DLPFC and predicted cell type proportions were compared to orthogonal RNAScope/IF measurements.

Conclusions

Bisque and hspe were the most accurate methods, were robust to differences in RNA library types and extractions. This multi-assay dataset showed that cell size differences, marker genes differentially quantified across RNA libraries, and cell composition variability in reference snRNA-seq impact the accuracy of current deconvolution methods.

Large-scale RNA-seq mining reveals ciclopirox triggers TDP-43 cryptic exons

Nuclear clearance and cytoplasmic aggregation of TDP-43 in neurons, initially identified in ALS-FTD, are hallmark pathological features observed across a spectrum of neurodegenerative diseases. We previously found that TDP-43 loss-of-function leads to the transcriptome-wide inclusion of deleterious cryptic exons in brains and biofluids post-mortem as well as during the presymptomatic stage of ALS-FTD, but upstream mechanisms that lead to TDP-43 dysregulation remain unclear. Here, we developed a web-based resource (SnapMine) to determine the levels of TDP-43 cryptic exon inclusion across hundreds of thousands of publicly available RNA sequencing datasets. We established cryptic exon inclusion across a variety of human cells and tissues to provide ground truth references for future studies on TDP-43 dysregulation. We then explored studies that were entirely unrelated to TDP-43 or neurodegeneration and found that ciclopirox olamine (CPX), an FDA-approved antifungal, can trigger the inclusion of TDP-43-associated cryptic exons in a variety of mouse and human primary cells. CPX induction of cryptic exon occurs via heavy metal toxicity and oxidative stress, suggesting that similar vulnerabilities could play a role in neurodegeneration. Our work demonstrates how diverse datasets can be linked through common biological features and underscores that public archives of sequencing data represent a vastly underutilized resource with tremendous potential for uncovering novel insights into complex biological mechanisms and diseases.

A transfer learning framework to elucidate the clinical relevance of altered proximal tubule cell states in kidney disease

The application of single-cell technologies in clinical nephrology remains elusive. We generated an atlas of transcriptionally defined cell types and cell states of human kidney disease by integrating single-cell signatures reported in the literature with newly generated signatures obtained from 5 patients with acute kidney injury. We used this information to develop kidney-specific cell-level information ExtractoR (K-CLIER), a transfer learning approach specifically tailored to evaluate the role of cell types/states on bulk RNAseq data. We validated the K-CLIER as a reliable computational framework to obtain a dimensionality reduction and to link clinical data with single-cell signatures. By applying K-CLIER on cohorts of patients with different kidney diseases, we identified the most relevant cell types associated with fibrosis and disease progression. This analysis highlighted the central role of altered proximal tubule cells in chronic kidney disease. Our study introduces a new strategy to exploit the power of single-cell technologies toward clinical applications.

Genotype prediction of 336,463 samples from public expression data

Tens of thousands of RNA-sequencing experiments comprising hundreds of thousands of individual samples have now been performed. These data represent a broad range of experimental conditions, sequencing technologies, and hypotheses under study. The Recount project has aggregated and uniformly processed hundreds of thousands of publicly available RNA-seq samples. Most of these samples only include RNA expression measurements; genotype data for these same samples would enable a wide range of analyses including variant prioritization, eQTL analysis, and studies of allele specific expression. Here, we developed a statistical model based on the existing reference and alternative read counts from the RNA-seq experiments available through Recount3 to predict genotypes at autosomal biallelic loci in coding regions. We demonstrate the accuracy of our model using large-scale studies that measured both gene expression and genotype genome-wide. We show that our predictive model is highly accurate with 99.5% overall accuracy, 99.6% major allele accuracy, and 90.4% minor allele accuracy. Our model is robust to tissue and study effects, provided the coverage is high enough. We applied this model to genotype all the samples in Recount 3 and provide the largest ready-to-use expression repository containing genotype information. We illustrate that the predicted genotype from RNA-seq data is sufficient to unravel the underlying population structure of samples in Recount3 using Principal Component Analysis.

Polycomb-mediated silencing of miR-8 is required for maintenance of intestinal stemness in Drosophila melanogaster

Balancing maintenance of self-renewal and differentiation is a key property of adult stem cells. The epigenetic mechanisms controlling this balance remain largely unknown. Herein, we report that the Polycomb Repressive Complex 2 (PRC2) is required for maintenance of the intestinal stem cell (ISC) pool in the adult female Drosophila melanogaster. We show that loss of PRC2 activity in ISCs by RNAi-mediated knockdown or genetic ablation of the enzymatic subunit Enhancer of zeste, E(z), results in loss of stemness and precocious differentiation of enteroblasts to enterocytes. Mechanistically, we have identified the microRNA miR-8 as a critical target of E(z)/PRC2-mediated tri-methylation of histone H3 at Lys27 (H3K27me3) and uncovered a dynamic relationship between E(z), miR-8 and Notch signaling in controlling stemness versus differentiation of ISCs. Collectively, these findings uncover a hitherto unrecognized epigenetic layer in the regulation of stem cell specification that safeguards intestinal homeostasis.

Noninvasive Detection of Neuroendocrine Prostate Cancer through Targeted Cell-free DNA Methylation

Castration-resistant prostate cancer (CRPC) is a heterogeneous disease associated with phenotypic subtypes that drive therapy response and outcome differences. Histologic transformation to castration-resistant neuroendocrine prostate cancer (CRPC-NE) is associated with distinct epigenetic alterations, including changes in DNA methylation. The current diagnosis of CRPC-NE is challenging and relies on metastatic biopsy. We developed a targeted DNA methylation assay to detect CRPC-NE using plasma cell-free DNA (cfDNA). The assay quantifies tumor content and provides a phenotype evidence score that captures diverse CRPC phenotypes, leveraging regions to inform transcriptional state. We tested the design in independent clinical cohorts (n = 222 plasma samples) and qualified it achieving an AUC > 0.93 for detecting pathology-confirmed CRPC-NE (n = 136). Methylation-defined cfDNA tumor content was associated with clinical outcomes in two prospective phase II clinical trials geared towards aggressive variant CRPC and CRPC-NE. These data support the application of targeted DNA methylation for CRPC-NE detection and patient stratification.

SIGNIFICANCE

Neuroendocrine prostate cancer is an aggressive subtype of treatment-resistant prostate cancer. Early detection is important, but the diagnosis currently relies on metastatic biopsy. We describe the development and validation of a plasma cell-free DNA targeted methylation panel that can quantify tumor fraction and identify patients with neuroendocrine prostate cancer noninvasively. This article is featured in Selected Articles from This Issue, p. 384.

Germline determinants of aberrant signaling pathways in cancer

Cancer is a complex disease influenced by a heterogeneous landscape of both germline genetic variants and somatic aberrations. While there is growing evidence suggesting an interplay between germline and somatic variants, and a substantial number of somatic aberrations in specific pathways are now recognized as hallmarks in many well-known forms of cancer, the interaction landscape between germline variants and the aberration of those pathways in cancer remains largely unexplored. Utilizing over 8500 human samples across 33 cancer types characterized by TCGA and considering binary traits defined using a large collection of somatic aberration profiles across ten well-known oncogenic signaling pathways, we conducted a series of GWAS and identified genome-wide and suggestive associations involving 276 SNPs. Among these, 94 SNPs revealed cis-eQTL links with cancer-related genes or with genes functionally correlated with the corresponding traits' oncogenic pathways. GWAS summary statistics for all tested traits were then used to construct a set of polygenic scores employing a customized computational strategy. Polygenic scores for 24 traits demonstrated significant performance and were validated using data from PCAWG and CCLE datasets. These scores showed prognostic value for clinical variables and exhibited significant effectiveness in classifying patients into specific cancer subtypes or stratifying patients with cancer-specific aggressive phenotypes. Overall, we demonstrate that germline genetics can describe patients' genetic liability to develop specific cancer molecular and clinical profiles.

Autonomous transposons tune their sequences to ensure somatic suppression

Transposable elements (TEs) are a major constituent of human genes, occupying approximately half of the intronic space. During pre-messenger RNA synthesis, intronic TEs are transcribed along with their host genes but rarely contribute to the final mRNA product because they are spliced out together with the intron and rapidly degraded. Paradoxically, TEs are an abundant source of RNA-processing signals through which they can create new introns1, and also functional2 or non-functional chimeric transcripts3. The rarity of these events implies the existence of a resilient splicing code that is able to suppress TE exonization without compromising host pre-mRNA processing. Here we show that SAFB proteins protect genome integrity by preventing retrotransposition of L1 elements while maintaining splicing integrity, via prevention of the exonization of previously integrated TEs. This unique dual role is possible because of L1's conserved adenosine-rich coding sequences that are bound by SAFB proteins. The suppressive activity of SAFB extends to tissue-specific, giant protein-coding cassette exons, nested genes and Tigger DNA transposons. Moreover, SAFB also suppresses LTR/ERV elements in species in which they are still active, such as mice and flies. A significant subset of splicing events suppressed by SAFB in somatic cells are activated in the testis, coinciding with low SAFB expression in postmeiotic spermatids. Reminiscent of the division of labour between innate and adaptive immune systems that fight external pathogens, our results uncover SAFB proteins as an RNA-based, pattern-guided, non-adaptive defence system against TEs in the soma, complementing the RNA-based, adaptive Piwi-interacting RNA pathway of the germline.

AutoTransOP: translating omics signatures without orthologue requirements using deep learning

The development of therapeutics and vaccines for human diseases requires a systematic understanding of human biology. Although animal and in vitro culture models can elucidate some disease mechanisms, they typically fail to adequately recapitulate human biology as evidenced by the predominant likelihood of clinical trial failure. To address this problem, we developed AutoTransOP, a neural network autoencoder framework, to map omics profiles from designated species or cellular contexts into a global latent space, from which germane information for different contexts can be identified without the typically imposed requirement of matched orthologues. This approach was found in general to perform at least as well as current alternative methods in identifying animal/culture-specific molecular features predictive of other contexts-most importantly without requiring homology matching. For an especially challenging test case, we successfully applied our framework to a set of inter-species vaccine serology studies, where 1-to-1 mapping between human and non-human primate features does not exist.

Aggregation of recount3 RNA-seq data improves inference of consensus and tissue-specific gene co-expression networks

Background

Gene co-expression networks (GCNs) describe relationships among expressed genes key to maintaining cellular identity and homeostasis. However, the small sample size of typical RNA-seq experiments which is several orders of magnitude fewer than the number of genes is too low to infer GCNs reliably. recount3, a publicly available dataset comprised of 316,443 uniformly processed human RNA-seq samples, provides an opportunity to improve power for accurate network reconstruction and obtain biological insight from the resulting networks.

Results

We compared alternate aggregation strategies to identify an optimal workflow for GCN inference by data aggregation and inferred three consensus networks: a universal network, a non-cancer network, and a cancer network in addition to 27 tissue context-specific networks. Central network genes from our consensus networks were enriched for evolutionarily constrained genes and ubiquitous biological pathways, whereas central context-specific network genes included tissue-specific transcription factors and factorization based on the hubs led to clustering of related tissue contexts. We discovered that annotations corresponding to context-specific networks inferred from aggregated data were enriched for trait heritability beyond known functional genomic annotations and were significantly more enriched when we aggregated over a larger number of samples.

Conclusion

This study outlines best practices for network GCN inference and evaluation by data aggregation. We recommend estimating and regressing confounders in each data set before aggregation and prioritizing large sample size studies for GCN reconstruction. Increased statistical power in inferring context-specific networks enabled the derivation of variant annotations that were enriched for concordant trait heritability independent of functional genomic annotations that are context-agnostic. While we observed strictly increasing held-out log-likelihood with data aggregation, we noted diminishing marginal improvements. Future directions aimed at alternate methods for estimating confounders and integrating orthogonal information from modalities such as Hi-C and ChIP-seq can further improve GCN inference.

CMAUP database update 2024: extended functional and association information of useful plants for biomedical research

Knowledge of the collective activities of individual plants together with the derived clinical effects and targeted disease associations is useful for plant-based biomedical research. To provide the information in complement to the established databases, we introduced a major update of CMAUP database, previously featured in NAR. This update includes (i) human transcriptomic changes overlapping with 1152 targets of 5765 individual plants, covering 74 diseases from 20 027 patient samples; (ii) clinical information for 185 individual plants in 691 clinical trials; (iii) drug development information for 4694 drug-producing plants with metabolites developed into approved or clinical trial drugs; (iv) plant and human disease associations (428 737 associations by target, 220 935 reversion of transcriptomic changes, 764 and 154121 associations by clinical trials of individual plants and plant ingredients); (v) the location of individual plants in the phylogenetic tree for navigating taxonomic neighbors, (vi) DNA barcodes of 3949 plants, (vii) predicted human oral bioavailability of plant ingredients by the established SwissADME and HobPre algorithm, (viii) 21-107% increase of CMAUP data over the previous version to cover 60 222 chemical ingredients, 7865 plants, 758 targets, 1399 diseases, 238 KEGG human pathways, 3013 gene ontologies and 1203 disease ontologies. CMAUP update version is freely accessible at https://bidd.group/CMAUP/index.html.

MAJIQlopedia: an encyclopedia of RNA splicing variations in human tissues and cancer

Quantification of RNA splicing variations based on RNA-Sequencing can reveal tissue- and disease-specific splicing patterns. To study such splicing variations, we introduce MAJIQlopedia, an encyclopedia of splicing variations that encompasses 86 human tissues and 41 cancer datasets. MAJIQlopedia reports annotated and unannotated splicing events for a total of 486 175 alternative splice junctions in normal tissues and 338 317 alternative splice junctions in cancer. This database, available at https://majiq.biociphers.org/majiqlopedia/, includes a user-friendly interface that provides graphical representations of junction usage quantification for each junction across all tissue or cancer types. To demonstrate case usage of MAJIQlopedia, we review splicing variations in genes WT1, MAPT and BIN1, which all have known tissue or cancer-specific splicing variations. We also use MAJIQlopedia to highlight novel splicing variations in FDX1 and MEGF9 in normal tissues, and we uncover a novel exon inclusion event in RPS6KA6 that only occurs in two cancer types. Users can download the database, request the addition of data to the webtool, or install a MAJIQlopedia server to integrate proprietary data. MAJIQlopedia can serve as a reference database for researchers seeking to understand what splicing variations exist in genes of interest, and those looking to understand tissue- or cancer-specific splice isoform usage.

NCBI GEO: archive for gene expression and epigenomics data sets: 23-year update

The Gene Expression Omnibus (GEO) is an international public repository that archives gene expression and epigenomics data sets generated by next-generation sequencing and microarray technologies. Data are typically submitted to GEO by researchers in compliance with widespread journal and funder mandates to make generated data publicly accessible. The resource handles raw data files, processed data files and descriptive metadata for over 200 000 studies and 6.5 million samples, all of which are indexed, searchable and downloadable. Additionally, GEO offers web-based tools that facilitate analysis and visualization of differential gene expression. This article presents the current status and recent advancements in GEO, including the generation of consistently computed gene expression count matrices for thousands of RNA-seq studies, and new interactive graphical plots in GEO2R that help users identify differentially expressed genes and assess data set quality. The GEO repository is built and maintained by the National Center for Biotechnology Information (NCBI), a division of the National Library of Medicine (NLM), and is publicly accessible at https://www.ncbi.nlm.nih.gov/geo/.

The landscape of SETBP1 gene expression and transcription factor activity across human tissues

The SET binding protein 1 (SETBP1) gene encodes a transcription factor (TF) involved in various cellular processes. Variants in SETBP1 can result in three different diseases determined by the introduction (germline vs. somatic) and location of the variant. Germline variants cause the ultra-rare pediatric Schinzel Giedion Syndrome (SGS) and SETBP1 haploinsufficiency disorder (SETBP1-HD), characterized by severe multisystemic abnormalities with neurodegeneration or a less severe brain phenotype accompanied by hypotonia and strabismus, respectively. Somatic variants in SETBP1 are associated with hematological malignancies and cancer development in other tissues in adults. To better understand the tissue-specific mechanisms involving SETBP1, we analyzed publicly available RNA-sequencing (RNA-seq) data from the Genotype-Tissue Expression (GTEx) project. We found SETBP1 and its known target genes were widely expressed across 31 adult human tissues. K-means clustering identified three distinct expression patterns of SETBP1 targets across tissues. Functional enrichment analysis (FEA) of each cluster revealed gene sets related to transcriptional regulation, DNA binding, and mitochondrial function. TF activity analysis of SETBP1 and its target TFs revealed tissue-specific TF activity, underscoring the role of tissue context-driven regulation and suggesting its impact in SETBP1-associated disease. In addition to uncovering tissue-specific molecular signatures of SETBP1 expression and TF activity, we provide a Shiny web application to facilitate exploring TF activity across human tissues for 758 TFs. This study provides insight into the landscape of SETBP1 expression and TF activity across 31 non-diseased human tissues and reveals tissue-specific expression and activity of SETBP1 and its targets. In conjunction with the web application we constructed, our framework enables researchers to generate hypotheses related to the role tissue backgrounds play with respect to gene expression and TF activity in different disease contexts.

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single-cell RNA-sequencing datasets

Deconvolution of cell mixtures in "bulk" transcriptomic samples from homogenate human tissue is important for understanding disease pathologies. However, several experimental and computational challenges impede transcriptomics-based deconvolution approaches using single-cell/nucleus RNA-seq reference atlases. Cells from the brain and blood have substantially different sizes, total mRNA, and transcriptional activities, and existing approaches may quantify total mRNA instead of cell type proportions. Further, standards are lacking for the use of cell reference atlases and integrative analyses of single-cell and spatial transcriptomics data. We discuss how to approach these key challenges with orthogonal "gold standard" datasets for evaluating deconvolution methods.

Pan-cancer analysis of antibody-drug conjugate targets and putative predictors of treatment response

BACKGROUND

Antibody-drug conjugates (ADCs) are a rapidly expanding class of compounds in oncology. Our goal was to assess the expression of ADC targets and potential downstream determining factors of activity across pan-cancer and normal tissues.

MATERIALS AND METHODS

ADCs in clinical trials (n = 121) were identified through ClinicalTrials.gov, corresponding to 54 targets. Genes potentially implicated in treatment response were identified in the literature. Gene expression from The Cancer Genome Atlas (9000+ cancers of 31 cancer types), the Genotype-Tissue Expression database (n = 19,000 samples from 31 normal tissue types), and the TNMplot.com (n = 12,494 unmatched primary and metastatic samples) were used in this analysis. To compare relative expression across and within tumour types we used pooled normal tissues as reference.

RESULTS

For most ADC targets, mRNA levels correlated with protein expression. Pan-cancer target expression distributions identified appealing cancer types for each ADC development. Co-expression of multiple targets was common and suggested opportunities for ADC combinations. Expression levels of genes potentially implicated in ADC response downstream of the target might provide additional information (e.g. TOP1 was highly expressed in many tumour types, including breast and lung cancers). Metastatic compared to primary tissues overexpressed some ADCs targets. Single sample "targetgram" plots were generated to visualise the expression of potentially competing ADC targets and resistance/sensitivity markers highlighting high inter-patient heterogeneity. Off-cancer target expression only partially explains adverse events, while expression of determinants of payload activity explained more of the observed toxicities.

CONCLUSION

Our findings draw attention to new therapeutic opportunities for ADCs that can be tested in the clinic and our web platform (https://tnmplot.com) can assist in prioritising upcoming ADC targets for clinical development.

Comparative RNA-sequencing analysis of the prostate in a mouse model of benign prostatic hyperplasia with bladder outlet obstruction

In ageing men, benign prostatic hyperplasia (BPH) is a chronic disease that leads to progressive lower urinary tract symptoms (LUTS) caused by obstruction of the bladder outlet (BOO). Patients with LUTS (such as increased frequency and urgency of urination) and complications of BOO (such as hydronephrosis and bladder stones) are at risk of serious health problems. BPH causes a rapidly rising burden of LUTS far exceeding that of other urological conditions. Treatment outcomes are unsatisfactory for BPH largely due to the lacking of fully understanding of the pathogenesis. Hormonal imbalances related to androgen and oestrogen can cause BPH, but the exact mechanism is still unknown, even the animal model is not fully understood. Additionally, there are no large-scale data to explain this mechanism. A BPH mouse model was established using mixed slow-release pellets of testosterone (T) and estradiol (E2), and we measured gene expression in mouse prostate tissue using RNA-seq, verified the results using qRT‒PCR, and used bioinformatics methods to analyse the differentially expressed genes (DEGs).

SpliceWiz: interactive analysis and visualization of alternative splicing in R

Alternative splicing (AS) is a crucial mechanism for regulating gene expression and isoform diversity in eukaryotes. However, the analysis and visualization of AS events from RNA sequencing data remains challenging. Most tools require a certain level of computer literacy and the available means of visualizing AS events, such as coverage and sashimi plots, have limitations and can be misleading. To address these issues, we present SpliceWiz, an R package with an interactive Shiny interface that allows easy and efficient AS analysis and visualization at scale. A novel normalization algorithm is implemented to aggregate splicing levels within sample groups, thereby allowing group differences in splicing levels to be accurately visualized. The tool also offers downstream gene ontology enrichment analysis, highlighting ASEs belonging to functional pathways of interest. SpliceWiz is optimized for speed and efficiency and introduces a new file format for coverage data storage that is more efficient than BigWig. Alignment files are processed orders of magnitude faster than other R-based AS analysis tools and on par with command-line tools. Overall, SpliceWiz streamlines AS analysis, enabling reliable identification of functionally relevant AS events for further characterization. SpliceWiz is a Bioconductor package and is also available on GitHub (https://github.com/alexchwong/SpliceWiz).

N-of-one differential gene expression without control samples using a deep generative model

Differential analysis of bulk RNA-seq data often suffers from lack of good controls. Here, we present a generative model that replaces controls, trained solely on healthy tissues. The unsupervised model learns a low-dimensional representation and can identify the closest normal representation for a given disease sample. This enables control-free, single-sample differential expression analysis. In breast cancer, we demonstrate how our approach selects marker genes and outperforms a state-of-the-art method. Furthermore, significant genes identified by the model are enriched in driver genes across cancers. Our results show that the in silico closest normal provides a more favorable comparison than control samples.

Statistical learning quantifies transposable element-mediated cis-regulation

BACKGROUND

Transposable elements (TEs) have colonized the genomes of most metazoans, and many TE-embedded sequences function as cis-regulatory elements (CREs) for genes involved in a wide range of biological processes from early embryogenesis to innate immune responses. Because of their repetitive nature, TEs have the potential to form CRE platforms enabling the coordinated and genome-wide regulation of protein-coding genes by only a handful of trans-acting transcription factors (TFs).

RESULTS

Here, we directly test this hypothesis through mathematical modeling and demonstrate that differences in expression at protein-coding genes alone are sufficient to estimate the magnitude and significance of TE-contributed cis-regulatory activities, even in contexts where TE-derived transcription fails to do so. We leverage hundreds of overexpression experiments and estimate that, overall, gene expression is influenced by TE-embedded CREs situated within approximately 500 kb of promoters. Focusing on the cis-regulatory potential of TEs within the gene regulatory network of human embryonic stem cells, we find that pluripotency-specific and evolutionarily young TE subfamilies can be reactivated by TFs involved in post-implantation embryogenesis. Finally, we show that TE subfamilies can be split into truly regulatorily active versus inactive fractions based on additional information such as matched epigenomic data, observing that TF binding may better predict TE cis-regulatory activity than differences in histone marks.

CONCLUSION

Our results suggest that TE-embedded CREs contribute to gene regulation during and beyond gastrulation. On a methodological level, we provide a statistical tool that infers TE-dependent cis-regulation from RNA-seq data alone, thus facilitating the study of TEs in the next-generation sequencing era.

Comparison of gene expression in living and postmortem human brain

Molecular mechanisms of neuropsychiatric disorders are challenging to study in human brain. For decades, the preferred model has been to study postmortem human brain samples despite the limitations they entail. A recent study generated RNA sequencing data from biopsies of prefrontal cortex from living patients with Parkinson's Disease and compared gene expression to postmortem tissue samples, from which they found vast differences between the two. This led the authors to question the utility of postmortem human brain studies. Through re-analysis of the same data, we unexpectedly found that the living brain tissue samples were of much lower quality than the postmortem samples across multiple standard metrics. We also performed simulations that illustrate the effects of ignoring RNA degradation in differential gene expression analyses, showing the effects can be substantial and of similar magnitude to what the authors find. For these reasons, we believe the authors' conclusions are unjustified. To the contrary, while opportunities to study gene expression in the living brain are welcome, evidence that this eclipses the value of postmortem analyses is not apparent.

Expression and splicing mediate distinct biological signals

BACKGROUND

Through alternative splicing, most human genes produce multiple isoforms in a cell-, tissue-, and disease-specific manner. Numerous studies show that alternative splicing is essential for development, diseases, and their treatments. Despite these important examples, the extent and biological relevance of splicing are currently unknown.

RESULTS

To solve this problem, we developed pairedGSEA and used it to profile transcriptional changes in 100 representative RNA-seq datasets. Our systematic analysis demonstrates that changes in splicing, on average, contribute to 48.1% of the biological signal in expression analyses. Gene-set enrichment analysis furthermore indicates that expression and splicing both convey shared and distinct biological signals.

CONCLUSIONS

These findings establish alternative splicing as a major regulator of the human condition and suggest that most contemporary RNA-seq studies likely miss out on critical biological insights. We anticipate our results will contribute to the transition from a gene-centric to an isoform-centric research paradigm.

The landscape of SETBP1 gene expression and transcription factor activity across human tissues

Background

The SET binding protein 1 (SETBP1) gene encodes a transcription factor (TF) involved in various cellular processes. Distinct SETBP1 variants have been linked to three different diseases. Germline variants cause the ultra-rare pediatric Schinzel Giedion Syndrome (SGS) and SETBP1 haploinsufficiency disorder (SETBP1-HD), characterized by severe multisystemic abnormalities with neurodegeneration or a less severe brain phenotype accompanied by hypotonia and strabismus, respectively. Somatic variants in SETBP1 are associated with hematological malignancies and cancer development in other tissues in adults.

Results

To better understand the tissue-specific mechanisms involving SETBP1, we analyzed publicly available RNA-sequencing data from the Genotype-Tissue Expression (GTEx) project. We found SETBP1, and its known target genes were widely expressed across 31 adult human tissues. K-means clustering identified three distinct expression patterns of SETBP1 targets across tissues. Functional enrichment analysis (FEA) of each cluster revealed gene sets related to transcription regulation, DNA binding, and mitochondrial function. TF activity analysis of SETBP1 and its target TFs revealed tissue-specific TF activity, underscoring the role of tissue context-driven regulation and suggesting its impact in SETBP1-associated disease. In addition to uncovering tissue-specific molecular signatures of SETBP1 expression and TF activity, we provide a Shiny web application to facilitate exploring TF activity across human tissues for 758 TFs.

Conclusions

This study provides insight into the landscape of SETBP1 expression and TF activity across 31 non-diseased human tissues and reveals tissue-specific expression and activity of SETBP1 and its targets. In conjunction with the web application we constructed, our framework enables researchers to generate hypotheses related to the role tissue backgrounds play with respect to gene expression and TF activity in different disease contexts.

Integration of 168,000 samples reveals global patterns of the human gut microbiome

Understanding the factors that shape variation in the human microbiome is a major goal of research in biology. While other genomics fields have used large, pre-compiled compendia to extract systematic insights requiring otherwise impractical sample sizes, there has been no comparable resource for the 16S rRNA sequencing data commonly used to quantify microbiome composition. To help close this gap, we have assembled a set of 168,484 publicly available human gut microbiome samples, processed with a single pipeline and combined into the largest unified microbiome dataset to date. We use this resource, which is freely available at microbiomap.org, to shed light on global variation in the human gut microbiome. We find that Firmicutes, particularly Bacilli and Clostridia, are almost universally present in the human gut. At the same time, the relative abundance of the 65 most common microbial genera differ between at least two world regions. We also show that gut microbiomes in undersampled world regions, such as Central and Southern Asia, differ significantly from the more thoroughly characterized microbiomes of Europe and Northern America. Moreover, humans in these overlooked regions likely harbor hundreds of taxa that have not yet been discovered due to this undersampling, highlighting the need for diversity in microbiome studies. We anticipate that this new compendium can serve the community and enable advanced applied and methodological research.

Functional Profiling of Soft Tissue Sarcoma Using Mechanistic Models

Soft tissue sarcoma is an umbrella term for a group of rare cancers that are difficult to treat. In addition to surgery, neoadjuvant chemotherapy has shown the potential to downstage tumors and prevent micrometastases. However, finding effective therapeutic targets remains a research challenge. Here, a previously developed computational approach called mechanistic models of signaling pathways has been employed to unravel the impact of observed changes at the gene expression level on the ultimate functional behavior of cells. In the context of such a mechanistic model, RNA-Seq counts sourced from the Recount3 resource, from The Cancer Genome Atlas (TCGA) Sarcoma project, and non-diseased sarcomagenic tissues from the Genotype-Tissue Expression (GTEx) project were utilized to investigate signal transduction activity through signaling pathways. This approach provides a precise view of the relationship between sarcoma patient survival and the signaling landscape in tumors and their environment. Despite the distinct regulatory alterations observed in each sarcoma subtype, this study identified 13 signaling circuits, or elementary sub-pathways triggering specific cell functions, present across all subtypes, belonging to eight signaling pathways, which served as predictors for patient survival. Additionally, nine signaling circuits from five signaling pathways that highlighted the modifications tumor samples underwent in comparison to normal tissues were found. These results describe the protective role of the immune system, suggesting an anti-tumorigenic effect in the tumor microenvironment, in the process of tumor cell detachment and migration, or the dysregulation of ion homeostasis. Also, the analysis of signaling circuit intermediary proteins suggests multiple strategies for therapy.

Gene regulatory Networks Reveal Sex Difference in Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) has been observed to have significant sex differences in incidence, prognosis, and response to therapy. However, the molecular mechanisms responsible for these disparities have not been investigated extensively. Sample-specific gene regulatory network methods were used to analyze RNA sequencing data from non-cancerous human lung samples from The Genotype Tissue Expression Project (GTEx) and lung adenocarcinoma primary tumor samples from The Cancer Genome Atlas (TCGA); results were validated on independent data. We observe that genes associated with key biological pathways including cell proliferation, immune response and drug metabolism are differentially regulated between males and females in both healthy lung tissue, as well as in tumor, and that these regulatory differences are further perturbed by tobacco smoking. We also uncovered significant sex bias in transcription factor targeting patterns of clinically actionable oncogenes and tumor suppressor genes, including AKT2 and KRAS. Using differentially regulated genes between healthy and tumor samples in conjunction with a drug repurposing tool, we identified several small-molecule drugs that might have sex-biased efficacy as cancer therapeutics and further validated this observation using an independent cell line database. These findings underscore the importance of including sex as a biological variable and considering gene regulatory processes in developing strategies for disease prevention and management.

Dopamine and schizophrenia from bench to bedside: Discovery of a striatal co-expression risk gene set that predicts in vivo measures of striatal function

Schizophrenia (SCZ) is characterized by a polygenic risk architecture implicating diverse molecular pathways important for synaptic function. However, how polygenic risk funnels through these pathways to translate into syndromic illness is unanswered. To evaluate biologically meaningful pathways of risk, we used tensor decomposition to characterize gene co-expression in post-mortem brain (of neurotypicals: N=154; patients with SCZ: N=84; and GTEX samples N=120) from caudate nucleus (CN), hippocampus (HP), and dorsolateral prefrontal cortex (DLPFC). We identified a CN-predominant gene set showing dopaminergic selectivity that was enriched for genes associated with clinical state and for genes associated with SCZ risk. Parsing polygenic risk score for SCZ based on this specific gene set (parsed-PRS), we found that greater pathway-specific SCZ risk predicted greater in vivo striatal dopamine synthesis capacity measured by [ 18 F]-FDOPA PET in three independent cohorts of neurotypicals and patients (total N=235) and greater fMRI striatal activation during reward anticipation in two additional independent neurotypical cohorts (total N=141). These results reveal a 'bench to bedside' translation of dopamine-linked genetic risk variation in driving in vivo striatal neurochemical and hemodynamic phenotypes that have long been implicated in the pathophysiology of SCZ.

Fluid flow shear stress and tissue remodeling-an orthodontic perspective: evidence synthesis and differential gene expression network analysis

Introduction: This study aimed to identify and analyze in vitro studies investigating the biological effect of fluid-flow shear stress (FSS) on cells found in the periodontal ligament and bone tissue. Method: We followed the PRISMA guideline for systematic reviews. A PubMed search strategy was developed, studies were selected according to predefined eligibility criteria, and the risk of bias was assessed. Relevant data related to cell source, applied FSS, and locus-specific expression were extracted. Based on this evidence synthesis and, as an original part of this work, analysis of differential gene expression using over-representation and network-analysis was performed. Five relevant publicly available gene expression datasets were analyzed using gene set enrichment analysis (GSEA). Result: A total of 6,974 articles were identified. Titles and abstracts were screened, and 218 articles were selected for full-text assessment. Finally, 120 articles were included in this study. Sample size determination and statistical analysis related to methodological quality and the ethical statement item in reporting quality were most frequently identified as high risk of bias. The analyzed studies mostly used custom-made fluid-flow apparatuses (61.7%). FSS was most frequently applied for 0.5 h, 1 h, or 2 h, whereas FSS magnitudes ranged from 6 to 20 dyn/cm2 depending on cell type and flow profile. Fluid-flow frequencies of 1 Hz in human cells and 1 and 5 Hz in mouse cells were mostly applied. FSS upregulated genes/metabolites responsible for tissue formation (AKT1, alkaline phosphatase, BGLAP, BMP2, Ca2+, COL1A1, CTNNB1, GJA1, MAPK1/MAPK3, PDPN, RUNX2, SPP1, TNFRSF11B, VEGFA, WNT3A) and inflammation (nitric oxide, PGE-2, PGI-2, PTGS1, PTGS2). Protein-protein interaction networks were constructed and analyzed using over-representation analysis and GSEA to identify shared signaling pathways. Conclusion: To our knowledge, this is the first review giving a comprehensive overview and discussion of methodological technical details regarding fluid flow application in 2D cell culture in vitro experimental conditions. Therefore, it is not only providing valuable information about cellular molecular events and their quantitative and qualitative analysis, but also confirming the reproducibility of previously published results.

Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates

SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17β-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.

Integrative omics approaches to advance rare disease diagnostics

Over the past decade high-throughput DNA sequencing approaches, namely whole exome and whole genome sequencing became a standard procedure in Mendelian disease diagnostics. Implementation of these technologies greatly facilitated diagnostics and shifted the analysis paradigm from variant identification to prioritisation and evaluation. The diagnostic rates vary widely depending on the cohort size, heterogeneity and disease and range from around 30% to 50% leaving the majority of patients undiagnosed. Advances in omics technologies and computational analysis provide an opportunity to increase these unfavourable rates by providing evidence for disease-causing variant validation and prioritisation. This review aims to provide an overview of the current application of several omics technologies including RNA-sequencing, proteomics, metabolomics and DNA-methylation profiling for diagnostics of rare genetic diseases in general and inborn errors of metabolism in particular.

Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates

SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17β-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.