A comprehensive assessment of RNA-seq accuracy, reproducibility and information content by the Sequencing Quality Control Consortium

Transipedia.org: k-mer-based exploration of large RNA sequencing datasets and application to cancer data

Indexing techniques relying on k-mers have proven effective in searching for RNA sequences across thousands of RNA-seq libraries, but without enabling direct RNA quantification. We show here that arbitrary RNA sequences can be quantified in seconds through their decomposition into k-mers, with a precision akin to that of conventional RNA quantification methods. Using an index of the Cancer Cell Line Encyclopedia (CCLE) collection consisting of 1019 RNA-seq samples, we show that k-mer indexing offers a powerful means to reveal non-reference sequences, and variant RNAs induced by specific gene alterations, for instance in splicing factors.

Assessing and mitigating batch effects in large-scale omics studies

Batch effects in omics data are notoriously common technical variations unrelated to study objectives, and may result in misleading outcomes if uncorrected, or hinder biomedical discovery if over-corrected. Assessing and mitigating batch effects is crucial for ensuring the reliability and reproducibility of omics data and minimizing the impact of technical variations on biological interpretation. In this review, we highlight the profound negative impact of batch effects and the urgent need to address this challenging problem in large-scale omics studies. We summarize potential sources of batch effects, current progress in evaluating and correcting them, and consortium efforts aiming to tackle them.

NBAtlas: A harmonized single-cell transcriptomic reference atlas of human neuroblastoma tumors

Neuroblastoma, a rare embryonic tumor arising from neural crest development, is responsible for 15% of pediatric cancer-related deaths. Recently, several single-cell transcriptome studies were performed on neuroblastoma patient samples to investigate the cell of origin and tumor heterogeneity. However, these individual studies involved a small number of tumors and cells, limiting the conclusions that could be drawn. To overcome this limitation, we integrated seven single-cell or single-nucleus datasets into a harmonized cell atlas covering 362,991 cells across 61 patients. We use this atlas to decipher the transcriptional landscape of neuroblastoma at single-cell resolution, revealing associations between transcriptomic profiles and clinical outcomes within the tumor compartment. In addition, we characterize the complex immune-cell landscape and uncover considerable heterogeneity among tumor-associated macrophages. Finally, we showcase the utility of our atlas as a resource by expanding it with additional data and using it as a reference for data-driven cell-type annotation.

Enhancing RNA-seq bias mitigation with the Gaussian self-benchmarking framework: towards unbiased sequencing data

BACKGROUND

RNA sequencing is a vital technique for analyzing RNA behavior in cells, but it often suffers from various biases that distort the data. Traditional methods to address these biases are typically empirical and handle them individually, limiting their effectiveness. Our study introduces the Gaussian Self-Benchmarking (GSB) framework, a novel approach that leverages the natural distribution patterns of guanine (G) and cytosine (C) content in RNA to mitigate multiple biases simultaneously. This method is grounded in a theoretical model, organizing k-mers based on their GC content and applying a Gaussian model for alignment to ensure empirical sequencing data closely match their theoretical distribution.

RESULTS

The GSB framework demonstrated superior performance in mitigating sequencing biases compared to existing methods. Testing with synthetic RNA constructs and real human samples showed that the GSB approach not only addresses individual biases more effectively but also manages co-existing biases jointly. The framework's reliance on accurately pre-determined parameters like mean and standard deviation of GC content distribution allows for a more precise representation of RNA samples. This results in improved accuracy and reliability of RNA sequencing data, enhancing our understanding of RNA behavior in health and disease.

CONCLUSIONS

The GSB framework presents a significant advancement in RNA sequencing analysis by providing a well-validated, multi-bias mitigation strategy. It functions independently from previously identified dataset flaws and sets a new standard for unbiased RNA sequencing results. This development enhances the reliability of RNA studies, broadening the potential for scientific breakthroughs in medicine and biology, particularly in genetic disease research and the development of targeted treatments.

Progress in toxicogenomics to protect human health

Toxicogenomics measures molecular features, such as transcripts, proteins, metabolites and epigenomic modifications, to understand and predict the toxicological effects of environmental and pharmaceutical exposures. Transcriptomics has become an integral tool in contemporary toxicology research owing to innovations in gene expression profiling that can provide mechanistic and quantitative information at scale. These data can be used to predict toxicological hazards through the use of transcriptomic biomarkers, network inference analyses, pattern-matching approaches and artificial intelligence. Furthermore, emerging approaches, such as high-throughput dose-response modelling, can leverage toxicogenomic data for human health protection even in the absence of predicting specific hazards. Finally, single-cell transcriptomics and multi-omics provide detailed insights into toxicological mechanisms. Here, we review the progress since the inception of toxicogenomics in applying transcriptomics towards toxicology testing and highlight advances that are transforming risk assessment.

Meta-analysis of RNA interaction profiles of RNA-binding protein using the RBPInper tool

Motivation

Recent RNA-centric experimental methods have significantly expanded our knowledge of proteins with known RNA-binding functions. However, the complete regulatory network and pathways for many of these RNA-binding proteins (RBPs) in different cellular contexts remain unknown. Although critical to understanding the role of RBPs in health and disease, experimentally mapping the RBP-RNA interactomes in every single context is an impossible task due the cost and manpower required. Additionally, identifying relevant RNAs bound by RBPs is challenging due to their diverse binding modes and function.

Results

To address these challenges, we developed RBP interaction mapper RBPInper an integrative framework that discovers global RBP interactome using statistical data fusion. Experiments on splicing factor proline and glutamine rich (SFPQ) datasets revealed cogent global SFPQ interactome. Several biological processes associated with this interactome were previously linked with SFPQ function. Furthermore, we conducted tests using independent dataset to assess the transferability of the SFPQ interactome to another context. The results demonstrated robust utility in generating interactomes that transfers to unseen cellular context. Overall, RBPInper is a fast and user-friendly method that enables a systems-level understanding of RBP functions by integrating multiple molecular datasets. The tool is designed with a focus on simplicity, minimal dependencies, and straightforward input requirements. This intentional design aims to empower everyday biologists, making it easy for them to incorporate the tool into their research.

Availability and implementation

The source code, documentation, and installation instructions as well as results for use case are freely available at https://github.com/AneneLab/RBPInper. A user can easily compile similar datasets for a target RBP.

Proteomic analysis of serum small extracellular vesicles identifies diagnostic biomarkers for neuroblastoma

Background

Neuroblastoma (NB) primarily arises in children who are <10 years of age, and originates from developing sympathetic nervous system, which results in tumors in adrenal glands and/or sympathetic ganglia. The diagnosis of NB involves a combination of laboratory and imaging tests, and biopsies. Small extracellular vesicles (sEVs) have gained attention as potential biomarkers for various types of tumors. Here, we performed proteomic analysis of serum sEVs and identified potential biomarkers for NB.

Methods

Label-free proteomics of serum sEVs were performed in the discovery phase. A bulk RNA-seq dataset of NB tissues was used to analyze the association between genes encoding sEVs proteins and prognosis. Potential biomarkers were validated via multiple reaction monitoring (MRM) or western blot analysis in the validation phase. A public single-cell RNA-seq (scRNA-seq) dataset was integrated to analyze the tissue origin of sEVs harboring biomarkers.

Results

A total of 104 differentially expressed proteins were identified in NB patients with label-free proteomics, and 26 potential biomarkers were validated with MRM analysis. Seven proteins BSG, HSP90AB1, SLC44A1, CHGA, ATP6V0A1, ITGAL and SELL showed the strong ability to distinguish NB patients from healthy controls and non-NB patients as well. Integrated analysis of scRNA-seq and sEVs proteomics revealed that these sEVs-derived biomarkers originated from different cell populations in tumor tissues.

Conclusion

sEVs-based biomarkers may aid the molecular diagnosis of NB, representing an innovative strategy to improve NB detection and management.

Trans2express - de novo transcriptome assembly pipeline optimized for gene expression analysis

BACKGROUND

As genomes of many eukaryotic species, especially plants, are large and complex, their de novo sequencing and assembly is still a difficult task despite progress in sequencing technologies. An alternative to genome assembly is the assembly of transcriptome, the set of RNA products of the expressed genes. While a bunch of de novo transcriptome assemblers exists, the challenges of transcriptomes (the existence of isoforms, the uneven expression levels across genes) complicates the generation of high-quality assemblies suitable for downstream analyses.

RESULTS

We developed Trans2express - a web-based tool and a pipeline of de novo hybrid transcriptome assembly and postprocessing based on rnaSPAdes with a set of subsequent filtrations. The pipeline was tested on Arabidopsis thaliana cDNA sequencing data obtained using Illumina and Oxford Nanopore Technologies platforms and three non-model plant species. The comparison of structural characteristics of the transcriptome assembly with reference Arabidopsis genome revealed the high quality of assembled transcriptome with 86.1% of Arabidopsis expressed genes assembled as a single contig. We tested the applicability of the transcriptome assembly for gene expression analysis. For both Arabidopsis and non-model species the results showed high congruence of gene expression levels and sets of differentially expressed genes between analyses based on genome and based on the transcriptome assembly.

CONCLUSIONS

We present Trans2express - a protocol for de novo hybrid transcriptome assembly aimed at recovering of a single transcript per gene. We expect this protocol to promote the characterization of transcriptomes and gene expression analysis in non-model plants and web-based tool to be of use to a wide range of plant biologists.

Targeted DNA-seq and RNA-seq of Reference Samples with Short-read and Long-read Sequencing

Next-generation sequencing (NGS) has revolutionized genomic research by enabling high-throughput, cost-effective genome and transcriptome sequencing accelerating personalized medicine for complex diseases, including cancer. Whole genome/transcriptome sequencing (WGS/WTS) provides comprehensive insights, while targeted sequencing is more cost-effective and sensitive. In comparison to short-read sequencing, which still dominates the field due to high speed and cost-effectiveness, long-read sequencing can overcome alignment limitations and better discriminate similar sequences from alternative transcripts or repetitive regions. Hybrid sequencing combines the best strengths of different technologies for a more comprehensive view of genomic/transcriptomic variations. Understanding each technology's strengths and limitations is critical for translating cutting-edge technologies into clinical applications. In this study, we sequenced DNA and RNA libraries of reference samples using various targeted DNA and RNA panels and the whole transcriptome on both short-read and long-read platforms. This study design enables a comprehensive analysis of sequencing technologies, targeting protocols, and library preparation methods. Our expanded profiling landscape establishes a reference point for assessing current sequencing technologies, facilitating informed decision-making in genomic research and precision medicine.

Genomic reproducibility in the bioinformatics era

In biomedical research, validating a scientific discovery hinges on the reproducibility of its experimental results. However, in genomics, the definition and implementation of reproducibility remain imprecise. We argue that genomic reproducibility, defined as the ability of bioinformatics tools to maintain consistent results across technical replicates, is essential for advancing scientific knowledge and medical applications. Initially, we examine different interpretations of reproducibility in genomics to clarify terms. Subsequently, we discuss the impact of bioinformatics tools on genomic reproducibility and explore methods for evaluating these tools regarding their effectiveness in ensuring genomic reproducibility. Finally, we recommend best practices to improve genomic reproducibility.

CleanUpRNAseq: An R/Bioconductor Package for Detecting and Correcting DNA Contamination in RNA-Seq Data

RNA sequencing (RNA-seq) has become a standard method for profiling gene expression, yet genomic DNA (gDNA) contamination carried over to the sequencing library poses a significant challenge to data integrity. Detecting and correcting this contamination is vital for accurate downstream analyses. Particularly, when RNA samples are scarce and invaluable, it becomes essential not only to identify but also to correct gDNA contamination to maximize the data's utility. However, existing tools capable of correcting gDNA contamination are limited and lack thorough evaluation. To fill the gap, we developed CleanUpRNAseq, which offers a comprehensive set of functionalities for identifying and correcting gDNA-contaminated RNA-seq data. Our package offers three correction methods for unstranded RNA-seq data and a dedicated approach for stranded data. Through rigorous validation on published RNA-seq datasets with known levels of gDNA contamination and real-world RNA-seq data, we demonstrate CleanUpRNAseq's efficacy in detecting and correcting detrimental levels of gDNA contamination across diverse library protocols. CleanUpRNAseq thus serves as a valuable tool for post-alignment quality assessment of RNA-seq data and should be integrated into routine workflows for RNA-seq data analysis. Its incorporation into OneStopRNAseq should significantly bolster the accuracy of gene expression quantification and differential expression analysis of RNA-seq data.

LncRNA RP11-818O24.3 promotes hair-follicle recovery via FGF2-PI3K/Akt signal pathway

A previous study indicated that patients with androgenic alopecia (AGA) have significantly reduced levels of LncRNA RP11-818O24.3. This study investigates whether LncRNA RP11-818O24.3 promotes hair-follicle recovery and its possible mechanism. Hair alteration and cutaneous histopathological changes induced by testosterone propionate were observed by H&E and bromodeoxyuridinc (BrdU) stain to evaluate the therapeutic effect of LncRNA RP11-818O24.3 in C57BL/6 J mice. The cellular viability was analyzed in LncRNA RP11-818O24.3-transfected human hair-follicle stem cells (HFSCs) in vitro. The signaling pathways and pro-proliferative factors were investigated by transcriptomic gene sequencing and qRT-PCR. LncRNA RP11-818O24.3 transfection successfully recovered hair growth and hair-follicle cells in AGA mice. In a series of HFSC studies in vitro, LncRNA RP11-818O24.3 transfection greatly promoted cellular proliferation and decreased cellular apoptosis. Transcriptome gene sequencing suggested that the phosphatidylinositol 3-kinase (PI3K)-Akt pathway was upregulated by LncRNA RP11-818O24.3. The qRT-PCR results showed that fibroblast growth factor (FGF)-2 was 14-times upregulated after LncRNA RP11-818O24.3 transfection. Hair-follicle recovery activity of LncRNA RP11-818O24.3 may involve the upregulation of FGF2 and PI3K-Akt to promote follicle stem cell survival. These data not only provide a theoretical basis for AGA development but also reveal a novel therapeutic method for AGA patients.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10616-024-00624-3.

Beyond the base pairs: comparative genome-wide DNA methylation profiling across sequencing technologies

Deoxyribonucleic acid (DNA) methylation plays a key role in gene regulation and is critical for development and human disease. Techniques such as whole-genome bisulfite sequencing (WGBS) and reduced representation bisulfite sequencing (RRBS) allow DNA methylation analysis at the genome scale, with Illumina NovaSeq 6000 and MGI Tech DNBSEQ-T7 being popular due to their efficiency and affordability. However, detailed comparative studies of their performance are not available. In this study, we constructed 60 WGBS and RRBS libraries for two platforms using different types of clinical samples and generated approximately 2.8 terabases of sequencing data. We systematically compared quality control metrics, genomic coverage, CpG methylation levels, intra- and interplatform correlations, and performance in detecting differentially methylated positions. Our results revealed that the DNBSEQ platform exhibited better raw read quality, although base quality recalibration indicated potential overestimation of base quality. The DNBSEQ platform also showed lower sequencing depth and less coverage uniformity in GC-rich regions than did the NovaSeq platform and tended to enrich methylated regions. Overall, both platforms demonstrated robust intra- and interplatform reproducibility for RRBS and WGBS, with NovaSeq performing better for WGBS, highlighting the importance of considering these factors when selecting a platform for bisulfite sequencing.

Transcriptomics and epigenetic data integration learning module on Google Cloud

Multi-omics (genomics, transcriptomics, epigenomics, proteomics, metabolomics, etc.) research approaches are vital for understanding the hierarchical complexity of human biology and have proven to be extremely valuable in cancer research and precision medicine. Emerging scientific advances in recent years have made high-throughput genome-wide sequencing a central focus in molecular research by allowing for the collective analysis of various kinds of molecular biological data from different types of specimens in a single tissue or even at the level of a single cell. Additionally, with the help of improved computational resources and data mining, researchers are able to integrate data from different multi-omics regimes to identify new prognostic, diagnostic, or predictive biomarkers, uncover novel therapeutic targets, and develop more personalized treatment protocols for patients. For the research community to parse the scientifically and clinically meaningful information out of all the biological data being generated each day more efficiently with less wasted resources, being familiar with and comfortable using advanced analytical tools, such as Google Cloud Platform becomes imperative. This project is an interdisciplinary, cross-organizational effort to provide a guided learning module for integrating transcriptomics and epigenetics data analysis protocols into a comprehensive analysis pipeline for users to implement in their own work, utilizing the cloud computing infrastructure on Google Cloud. The learning module consists of three submodules that guide the user through tutorial examples that illustrate the analysis of RNA-sequence and Reduced-Representation Bisulfite Sequencing data. The examples are in the form of breast cancer case studies, and the data sets were procured from the public repository Gene Expression Omnibus. The first submodule is devoted to transcriptomics analysis with the RNA sequencing data, the second submodule focuses on epigenetics analysis using the DNA methylation data, and the third submodule integrates the two methods for a deeper biological understanding. The modules begin with data collection and preprocessing, with further downstream analysis performed in a Vertex AI Jupyter notebook instance with an R kernel. Analysis results are returned to Google Cloud buckets for storage and visualization, removing the computational strain from local resources. The final product is a start-to-finish tutorial for the researchers with limited experience in multi-omics to integrate transcriptomics and epigenetics data analysis into a comprehensive pipeline to perform their own biological research.This manuscript describes the development of a resource module that is part of a learning platform named ``NIGMS Sandbox for Cloud-based Learning'' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox [16] at the beginning of this Supplement. This module delivers learning materials on the analysis of bulk and single-cell ATAC-seq data in an interactive format that uses appropriate cloud resources for data access and analyses.

HIGHLIGHTS

A real-world multi-center RNA-seq benchmarking study using the Quartet and MAQC reference materials

Translating RNA-seq into clinical diagnostics requires ensuring the reliability and cross-laboratory consistency of detecting clinically relevant subtle differential expressions, such as those between different disease subtypes or stages. As part of the Quartet project, we present an RNA-seq benchmarking study across 45 laboratories using the Quartet and MAQC reference samples spiked with ERCC controls. Based on multiple types of 'ground truth', we systematically assess the real-world RNA-seq performance and investigate the influencing factors involved in 26 experimental processes and 140 bioinformatics pipelines. Here we show greater inter-laboratory variations in detecting subtle differential expressions among the Quartet samples. Experimental factors including mRNA enrichment and strandedness, and each bioinformatics step, emerge as primary sources of variations in gene expression. We underscore the profound influence of experimental execution, and provide best practice recommendations for experimental designs, strategies for filtering low-expression genes, and the optimal gene annotation and analysis pipelines. In summary, this study lays the foundation for developing and quality control of RNA-seq for clinical diagnostic purposes.

Accurate isoform quantification by joint short- and long-read RNA-sequencing

Accurate quantification of transcript isoforms is crucial for understanding gene regulation, functional diversity, and cellular behavior. Existing RNA sequencing methods have significant limitations: short-read (SR) sequencing provides high depth but struggles with isoform deconvolution, whereas long-read (LR) sequencing offers isoform resolution at the cost of lower depth, higher noise, and technical biases. Addressing this gap, we introduce Multi-Platform Aggregation and Quantification of Transcripts (MPAQT), a generative model that combines the complementary strengths of different sequencing platforms to achieve state-of-the-art isoform-resolved transcript quantification, as demonstrated by extensive simulations and experimental benchmarks. By applying MPAQT to an in vitro model of human embryonic stem cell differentiation into cortical neurons, followed by machine learning-based modeling of transcript abundances, we show that untranslated regions (UTRs) are major determinants of isoform proportion and exon usage; this effect is mediated through isoform-specific sequence features embedded in UTRs, which likely interact with RNA-binding proteins that modulate mRNA stability. These findings highlight MPAQT's potential to enhance our understanding of transcriptomic complexity and underline the role of splicing-independent post-transcriptional mechanisms in shaping the isoform and exon usage landscape of the cell.

consICA: an R package for robust reference-free deconvolution of multi-omics data

Motivation

Deciphering molecular signals from omics data helps understanding cellular processes and disease progression. Effective algorithms for extracting these signals are essential, with a strong emphasis on robustness and reproducibility.

Results

R/Bioconductor package consICA implements consensus independent component analysis (ICA)-a data-driven deconvolution method to decompose heterogeneous omics data and extract features suitable for patient stratification and multimodal data integration. The method separates biologically relevant molecular signals from technical effects and provides information about the cellular composition and biological processes. Build-in annotation, survival analysis, and report generation provide useful tools for the interpretation of extracted signals. The implementation of parallel computing in the package ensures efficient analysis using modern multicore systems. The package offers a reproducible and efficient data-driven solution for the analysis of complex molecular profiles, with significant implications for cancer research.

Availability and implementation

The package is implemented in R and available under MIT license at Bioconductor (https://bioconductor.org/packages/consICA) or at GitHub (https://github.com/biomod-lih/consICA).

Quartet RNA reference materials improve the quality of transcriptomic data through ratio-based profiling

Certified RNA reference materials are indispensable for assessing the reliability of RNA sequencing to detect intrinsically small biological differences in clinical settings, such as molecular subtyping of diseases. As part of the Quartet Project for quality control and data integration of multi-omics profiling, we established four RNA reference materials derived from immortalized B-lymphoblastoid cell lines from four members of a monozygotic twin family. Additionally, we constructed ratio-based transcriptome-wide reference datasets between two samples, providing cross-platform and cross-laboratory 'ground truth'. Investigation of the intrinsically subtle biological differences among the Quartet samples enables sensitive assessment of cross-batch integration of transcriptomic measurements at the ratio level. The Quartet RNA reference materials, combined with the ratio-based reference datasets, can serve as unique resources for assessing and improving the quality of transcriptomic data in clinical and biological settings.

Reliable biological and multi-omics research through biometrology

Metrology is the science of measurement and its applications, whereas biometrology is the science of biological measurement and its applications. Biometrology aims to achieve accuracy and consistency of biological measurements by focusing on the development of metrological traceability, biological reference measurement procedures, and reference materials. Irreproducibility of biological and multi-omics research results from different laboratories, platforms, and analysis methods is hampering the translation of research into clinical uses and can often be attributed to the lack of biologists' attention to the general principles of metrology. In this paper, the progresses of biometrology including metrology on nucleic acid, protein, and cell measurements and its impacts on the improvement of reliability and comparability in biological research are reviewed. Challenges in obtaining more reliable biological and multi-omics measurements due to the lack of primary reference measurement procedures and new standards for biological reference materials faced by biometrology are discussed. In the future, in addition to establishing reliable reference measurement procedures, developing reference materials from single or multiple parameters to multi-omics scale should be emphasized. Thinking in way of biometrology is warranted for facilitating the translation of high-throughput omics research into clinical practices.

Multi-omics data integration using ratio-based quantitative profiling with Quartet reference materials

Characterization and integration of the genome, epigenome, transcriptome, proteome and metabolome of different datasets is difficult owing to a lack of ground truth. Here we develop and characterize suites of publicly available multi-omics reference materials of matched DNA, RNA, protein and metabolites derived from immortalized cell lines from a family quartet of parents and monozygotic twin daughters. These references provide built-in truth defined by relationships among the family members and the information flow from DNA to RNA to protein. We demonstrate how using a ratio-based profiling approach that scales the absolute feature values of a study sample relative to those of a concurrently measured common reference sample produces reproducible and comparable data suitable for integration across batches, labs, platforms and omics types. Our study identifies reference-free 'absolute' feature quantification as the root cause of irreproducibility in multi-omics measurement and data integration and establishes the advantages of ratio-based multi-omics profiling with common reference materials.

Asporin and CD109, expressed in the injured neonatal spinal cord, attenuate axonal re-growth in vitro

Axonal regeneration is restricted in adults and causes irreversible motor dysfunction following spinal cord injury (SCI). In contrast, neonates have prominent regenerative potential and can restore their neural function. Although the distinct cellular responses in neonates have been studied, how they contribute to neural recovery remains unclear. To assess whether the secreted molecules in neonatal SCI can enhance neural regeneration, we re-analyzed the previously performed single-nucleus RNA-seq (snRNA-seq) and focused on Asporin and Cd109, the highly expressed genes in the injured neonatal spinal cord. In the present study, we showed that both these molecules were expressed in the injured spinal cords of adults and neonates. We treated the cortical neurons with recombinant Asporin or CD109 to observe their direct effects on neurons in vitro. We demonstrated that these molecules enhance neurite outgrowth in neurons. However, these molecules did not enhance re-growth of severed axons. Our results suggest that Asporin and CD109 influence neurites at the lesion site, rather than promoting axon regeneration, to restore neural function in neonates after SCI.

Distinct serum exosomal miRNA profiles detected in acute and asymptomatic dengue infections: A community-based study in Baiyun District, Guangzhou

Background

In recent years, research on exosomal miRNAs has provided new insights into exploring the mechanism of viral infection and disease prevention. This study aimed to investigate the serum exosomal miRNA expression profile of dengue-infected individuals through a community survey of dengue virus (DENV) infection.

Methods

A seroprevalence study of 1253 healthy persons was first conducted to ascertain the DENV infection status in Baiyun District, Guangzhou. A total of 18 serum samples, including 6 healthy controls (HC), 6 asymptomatic DENV infections (AsymptDI), and 6 confirmed dengue fever patients (AcuteDI), were collected for exosome isolation and then sRNA sequencing. Through bioinformatics analysis, we discovered distinct serum exosomal miRNA profiles among the different groups and identified differentially expressed miRNAs (DEMs). These findings were further validated by qRT-PCR.

Results

The community survey of DENV infection indicated that the DENV IgG antibody positivity rate among the population was 11.97 % in the study area, with asymptomatic infected individuals accounting for 93.06 % of the anti-DENV IgG positives. The age and Guangzhou household registration were associated with DENV IgG antibody positivity by logistic regression analysis. Distinct miRNA profiles were observed between healthy individuals and DENV infections. A total of 1854 miRNAs were identified in 18 serum exosome samples from the initial analysis of the sequencing data. Comparative analysis revealed 23 DEMs comprising 5 upregulated and 18 downregulated miRNAs in the DENV-infected group (mergedDI). In comparison to AcuteDI, 18 upregulated miRNAs were identified in AsymptDI. Moreover, functional enrichment of the predicted target genes of DEMs indicated that these miRNAs were involved in biological processes and pathways related to cell adhesion, focal adhesion, endocytosis, and ECM-receptor interaction. Eight DEMs were validated by qRT-PCR.

Conclusion

The Baiyun District of Guangzhou exhibits a notable proportion of asymptomatic DENV infections as suggested in other research, highlighting the need for enhanced monitoring and screening of asymptomatic persons and the elderly. Differential miRNA expression among healthy, symptomatic and asymptomatic DENV-infected individuals suggests their potential as biomarkers for distinguishing DENV infection and offers new avenues of investigating the mechanisms underlying DENV asymptomatic infections.

Comprehensive assessment of mRNA isoform detection methods for long-read sequencing data

The advancement of Long-Read Sequencing (LRS) techniques has significantly increased the length of sequencing to several kilobases, thereby facilitating the identification of alternative splicing events and isoform expressions. Recently, numerous computational tools for isoform detection using long-read sequencing data have been developed. Nevertheless, there remains a deficiency in comparative studies that systemically evaluate the performance of these tools, which are implemented with different algorithms, under various simulations that encompass potential influencing factors. In this study, we conducted a benchmark analysis of thirteen methods implemented in nine tools capable of identifying isoform structures from long-read RNA-seq data. We evaluated their performances using simulated data, which represented diverse sequencing platforms generated by an in-house simulator, RNA sequins (sequencing spike-ins) data, as well as experimental data. Our findings demonstrate IsoQuant as a highly effective tool for isoform detection with LRS, with Bambu and StringTie2 also exhibiting strong performance. These results offer valuable guidance for future research on alternative splicing analysis and the ongoing improvement of tools for isoform detection using LRS data.

A human neural crest model reveals the developmental impact of neuroblastoma-associated chromosomal aberrations

Early childhood tumours arise from transformed embryonic cells, which often carry large copy number alterations (CNA). However, it remains unclear how CNAs contribute to embryonic tumourigenesis due to a lack of suitable models. Here we employ female human embryonic stem cell (hESC) differentiation and single-cell transcriptome and epigenome analysis to assess the effects of chromosome 17q/1q gains, which are prevalent in the embryonal tumour neuroblastoma (NB). We show that CNAs impair the specification of trunk neural crest (NC) cells and their sympathoadrenal derivatives, the putative cells-of-origin of NB. This effect is exacerbated upon overexpression of MYCN, whose amplification co-occurs with CNAs in NB. Moreover, CNAs potentiate the pro-tumourigenic effects of MYCN and mutant NC cells resemble NB cells in tumours. These changes correlate with a stepwise aberration of developmental transcription factor networks. Together, our results sketch a mechanistic framework for the CNA-driven initiation of embryonal tumours.

Normalization of RNA-Seq data using adaptive trimmed mean with multi-reference

The normalization of RNA sequencing data is a primary step for downstream analysis. The most popular method used for the normalization is the trimmed mean of M values (TMM) and DESeq. The TMM tries to trim away extreme log fold changes of the data to normalize the raw read counts based on the remaining non-deferentially expressed genes. However, the major problem with the TMM is that the values of trimming factor M are heuristic. This paper tries to estimate the adaptive value of M in TMM based on Jaeckel's Estimator, and each sample acts as a reference to find the scale factor of each sample. The presented approach is validated on SEQC, MAQC2, MAQC3, PICKRELL and two simulated datasets with two-group and three-group conditions by varying the percentage of differential expression and the number of replicates. The performance of the present approach is compared with various state-of-the-art methods, and it is better in terms of area under the receiver operating characteristic curve and differential expression.

Histone lysine demethylase 4 family proteins maintain the transcriptional program and adrenergic cellular state of MYCN-amplified neuroblastoma

Neuroblastoma with MYCN amplification (MNA) is a high-risk disease that has a poor survival rate. Neuroblastoma displays cellular heterogeneity, including more differentiated (adrenergic) and more primitive (mesenchymal) cellular states. Here, we demonstrate that MYCN oncoprotein promotes a cellular state switch in mesenchymal cells to an adrenergic state, accompanied by induction of histone lysine demethylase 4 family members (KDM4A-C) that act in concert to control the expression of MYCN and adrenergic core regulatory circulatory (CRC) transcription factors. Pharmacologic inhibition of KDM4 blocks expression of MYCN and the adrenergic CRC transcriptome with genome-wide induction of transcriptionally repressive H3K9me3, resulting in potent anticancer activity against neuroblastomas with MNA by inducing neuroblastic differentiation and apoptosis. Furthermore, a short-term KDM4 inhibition in combination with conventional, cytotoxic chemotherapy results in complete tumor responses of xenografts with MNA. Thus, KDM4 blockade may serve as a transformative strategy to target the adrenergic CRC dependencies in MNA neuroblastomas.

Comparison of RNA-Seq and microarray in the prediction of protein expression and survival prediction

Gene expression profiling using RNA-sequencing (RNA-seq) and microarray technologies is widely used in cancer research to identify biomarkers for clinical endpoint prediction. We compared the performance of these two methods in predicting protein expression and clinical endpoints using The Cancer Genome Atlas (TCGA) datasets of lung cancer, colorectal cancer, renal cancer, breast cancer, endometrial cancer, and ovarian cancer. We calculated the correlation coefficients between gene expression measured by RNA-seq or microarray and protein expression measured by reverse phase protein array (RPPA). In addition, after selecting the top 103 survival-related genes, we compared the random forest survival prediction model performance across test platforms and cancer types. Both RNA-seq and microarray data were retrieved from TCGA dataset. Most genes showed similar correlation coefficients between RNA-seq and microarray, but 16 genes exhibited significant differences between the two methods. The BAX gene was recurrently found in colorectal cancer, renal cancer, and ovarian cancer, and the PIK3CA gene belonged to renal cancer and breast cancer. Furthermore, the survival prediction model using microarray was better than the RNA-seq model in colorectal cancer, renal cancer, and lung cancer, but the RNA-seq model was better in ovarian and endometrial cancer. Our results showed good correlation between mRNA levels and protein measured by RPPA. While RNA-seq and microarray performance were similar, some genes showed differences, and further clinical significance should be evaluated. Additionally, our survival prediction model results were controversial.

Integrative analysis of neuroblastoma by single-cell RNA sequencing identifies the NECTIN2-TIGIT axis as a target for immunotherapy

Pediatric patients with high-risk neuroblastoma have poor survival rates and urgently need more effective treatment options with less side effects. Since novel and improved immunotherapies may fill this need, we dissect the immunoregulatory interactions in neuroblastoma by single-cell RNA-sequencing of 24 tumors (10 pre- and 14 post-chemotherapy, including 5 pairs) to identify strategies for optimizing immunotherapy efficacy. Neuroblastomas are infiltrated by natural killer (NK), T and B cells, and immunosuppressive myeloid populations. NK cells show reduced cytotoxicity and T cells have a dysfunctional profile. Interaction analysis reveals a vast immunoregulatory network and identifies NECTIN2-TIGIT as a crucial immune checkpoint. Combined blockade of TIGIT and PD-L1 significantly reduces neuroblastoma growth, with complete responses (CR) in vivo. Moreover, addition of TIGIT+PD-L1 blockade to standard relapse treatment in a chemotherapy-resistant Th-ALKF1174L/MYCN 129/SvJ syngeneic model induces CR. In conclusion, our integrative analysis provides promising targets and a rationale for immunotherapeutic combination strategies.

A type 1 immunity-restricted promoter of the IL-33 receptor gene directs antiviral T-cell responses

The pleiotropic alarmin interleukin-33 (IL-33) drives type 1, type 2 and regulatory T-cell responses via its receptor ST2. Subset-specific differences in ST2 expression intensity and dynamics suggest that transcriptional regulation is key in orchestrating the context-dependent activity of IL-33-ST2 signaling in T-cell immunity. Here, we identify a previously unrecognized alternative promoter in mice and humans that is located far upstream of the curated ST2-coding gene and drives ST2 expression in type 1 immunity. Mice lacking this promoter exhibit a selective loss of ST2 expression in type 1- but not type 2-biased T cells, resulting in impaired expansion of cytotoxic T cells (CTLs) and T-helper 1 cells upon viral infection. T-cell-intrinsic IL-33 signaling via type 1 promoter-driven ST2 is critical to generate a clonally diverse population of antiviral short-lived effector CTLs. Thus, lineage-specific alternative promoter usage directs alarmin responsiveness in T-cell subsets and offers opportunities for immune cell-specific targeting of the IL-33-ST2 axis in infections and inflammatory diseases.

Early-stage idiopathic Parkinson's disease is associated with reduced circular RNA expression

Neurodegeneration in Parkinson's disease (PD) precedes diagnosis by years. Early neurodegeneration may be reflected in RNA levels and measurable as a biomarker. Here, we present the largest quantification of whole blood linear and circular RNAs (circRNA) in early-stage idiopathic PD, using RNA sequencing data from two cohorts (PPMI = 259 PD, 161 Controls; ICICLE-PD = 48 PD, 48 Controls). We identified a replicable increase in TMEM252 and LMNB1 gene expression in PD. We identified novel differences in the expression of circRNAs from ESYT2, BMS1P1 and CCDC9, and replicated trends of previously reported circRNAs. Overall, using circRNA as a diagnostic biomarker in PD did not show any clear improvement over linear RNA, minimising its potential clinical utility. More interestingly, we observed a general reduction in circRNA expression in both PD cohorts, accompanied by an increase in RNASEL expression. This imbalance implicates the activation of an innate antiviral immune response and suggests a previously unknown aspect of circRNA regulation in PD.

A uniform data processing pipeline enables harmonized nanoparticle protein corona analysis across proteomics core facilities

Protein corona, a layer of biomolecules primarily comprising proteins, forms dynamically on nanoparticles in biological fluids and is crucial for predicting nanomedicine safety and efficacy. The protein composition of the corona layer is typically analyzed using liquid chromatography-mass spectrometry (LC-MS/MS). Our recent study, involving identical samples analyzed by 17 proteomics facilities, highlighted significant data variability, with only 1.8% of proteins consistently identified across these centers. Here, we implement an aggregated database search unifying parameters such as variable modifications, enzyme specificity, number of allowed missed cleavages and a stringent 1% false discovery rate at the protein and peptide levels. Such uniform search dramatically harmonizes the proteomics data, increasing the reproducibility and the percentage of consistency-identified unique proteins across distinct cores. Specifically, out of the 717 quantified proteins, 253 (35.3%) are shared among the top 5 facilities (and 16.2% among top 11 facilities). Furthermore, we note that reduction and alkylation are important steps in protein corona sample processing and as expected, omitting these steps reduces the number of total quantified peptides by around 20%. These findings underscore the need for standardized procedures in protein corona analysis, which is vital for advancing clinical applications of nanoscale biotechnologies.

Transcriptomic Profiles of Pectoralis major Muscles Affected by Spaghetti Meat and Woody Breast in Broiler Chickens

Spaghetti meat (SM) and woody breast (WB) are breast muscle myopathies of broiler chickens, characterized by separation of myofibers and by fibrosis, respectively. This study sought to investigate the transcriptomic profiles of breast muscles affected by SM and WB. Targeted sampling was conducted on a flock to obtain 10 WB, 10 SM, and 10 Normal Pectoralis major muscle samples from 37-day-old male chickens. Total RNA was extracted, cDNA was used for pair-end sequencing, and differentially expressed genes (DEGs) were determined by a false discovery rate of <0.1 and a >1.5-fold change. Principal component and heatmap cluster analyses showed that the SM and WB samples clustered together. No DEGs were observed between SM and WB fillets, while a total of 4018 and 2323 DEGs were found when comparing SM and WB, respectively, against Normal samples. In both the SM and WB samples, Gene Ontology terms associated with extracellular environment and immune response were enriched. The KEGG analysis showed enrichment of cytokine-cytokine receptor interaction and extracellular matrix-receptor interaction pathways in both myopathies. Although SM and WB are macroscopically different, the similar transcriptomic profiles suggest that these conditions may share a common pathogenesis. This is the first study to compare the transcriptomes of SM and WB, and it showed that, while both myopathies had profiles different from the normal breast muscle, SM and WB were similar, with comparable enriched metabolic pathways and processes despite presenting markedly different macroscopic features.

MYCN and SNRPD3 cooperate to maintain a balance of alternative splicing events that drives neuroblastoma progression

Many of the pro-tumorigenic functions of the oncogene MYCN are attributed to its regulation of global gene expression programs. Alternative splicing is another important regulator of gene expression and has been implicated in neuroblastoma development, however, the molecular mechanisms remain unknown. We found that MYCN up-regulated the expression of the core spliceosomal protein, SNRPD3, in models of neuroblastoma initiation and progression. High mRNA expression of SNRPD3 in human neuroblastoma tissues was a strong, independent prognostic factor for poor patient outcome. Repression of SNRPD3 expression correlated with loss of colony formation in vitro and reduced tumorigenicity in vivo. The effect of SNRPD3 on cell viability was in part dependent on MYCN as an oncogenic co-factor. RNA-sequencing revealed a global increase in the number of genes being differentially spliced when MYCN was overexpressed. Surprisingly, depletion of SNRPD3 in the presence of overexpressed MYCN further increased differential splicing, particularly of cell cycle regulators, such as BIRC5 and CDK10. MYCN directly bound SNRPD3, and the protein arginine methyltransferase, PRMT5, consequently increasing SNRPD3 methylation. Indeed, the PRMT5 inhibitor, JNJ-64619178, reduced cell viability and SNRPD3 methylation in neuroblastoma cells with high SNRPD3 and MYCN expression. Our findings demonstrate a functional relationship between MYCN and SNRPD3, which maintains the fidelity of MYCN-driven alternative splicing in the narrow range required for neuroblastoma cell growth. SNRPD3 methylation and its protein-protein interface with MYCN represent novel therapeutic targets. Hypothetical model for SNRPD3 as a co-factor for MYCN oncogenesis. SNRPD3 and MYCN participate in a regulatory loop to balance splicing fidelity in neuroblastoma cells. First MYCN transactivates SNRPD3 to lead to high-level expression. Second, SNRPD3 and MYCN form a protein complex involving PRMT5. Third, this leads to balanced alterative splicing (AS) activitiy that is favorable to neuroblastoma. Together this forms as a therapeutic vulnerability where SNRPD3 perturbation or PRMT5 inhibitors are selectively toxic to neuroblastoma by conditionally disturbing splicing activity.

PCR-based analytics of gene therapies using adeno-associated virus vectors: Considerations for cGMP method development

The field of gene therapy has evolved and improved so that today the treatment of thousands of genetic diseases is now possible. An integral aspect of the drug development process is generating analytical methods to be used throughout clinical and commercial manufacturing. Enumeration and identification assays using genetic testing are critical to ensure the safety, efficacy, and stability of many active pharmaceutical ingredients. While nucleic acid-based methods are already reliable and rapid, there are unique biological, technological, and regulatory aspects in gene therapies that must be considered. This review surveys aspects of method development and validation using nucleic acid-based testing of gene therapies by focusing on adeno-associated virus (AAV) vectors and their co-transfection factors. Key differences between quantitative PCR and droplet digital technologies are discussed to show how improvements can be made while still adhering to regulatory guidance. Example validation parameters for AAV genome titers are described to demonstrate the scope of analytical development. Finally, several areas for improving analytical testing are presented to inspire future innovation, including next-generation sequencing and artificial intelligence. Reviewing the broad characteristics of gene therapy assessment serves as an introduction for new researchers, while clarifying processes for professionals already involved in pharmaceutical manufacturing.

Ambiguous genes due to aligners and their impact on RNA-seq data analysis

The main scope of the study is ambiguous genes, i.e. genes whose expression is difficult to estimate from the data produced by next-generation sequencing technologies. We focused on the RNA sequencing (RNA-Seq) type of experiment performed on the Illumina platform. It is crucial to identify such genes and understand the cause of their difficulty, as these genes may be involved in some diseases. By giving misleading results, they could contribute to a misunderstanding of the cause of certain diseases, which could lead to inappropriate treatment. We thought that the ambiguous genes would be difficult to map because of their complex structure. So we looked at RNA-seq analysis using different mappers to find genes that would have different measurements from the aligners. We were able to identify such genes using a generalized linear model with two factors: mappers and groups introduced by the experiment. A large proportion of ambiguous genes are pseudogenes. High sequence similarity of pseudogenes to functional genes may indicate problems in alignment procedures. In addition, predictive analysis verified the performance of difficult genes in classification. The effectiveness of classifying samples into specific groups was compared, including the expression of difficult and not difficult genes as covariates. In almost all cases considered, ambiguous genes have less predictive power.

Targeting TRIM24 promotes neuroblastoma differentiation and decreases tumorigenicity via LSD1/CoREST complex

PURPOSE

High-risk neuroblastoma (NB) still has an unfavorable prognosis and inducing NB differentiation is a potential strategy in clinical treatment, yet underlying mechanisms are still elusive. Here we identify TRIM24 as an important regulator of NB differentiation.

METHODS

Multiple datasets and clinical specimens were analyzed to define the role of TRIM24 in NB. The effects of TRIM24 on differentiation and growth of NB were determined by cell morphology, spheres formation, soft agar assay, and subcutaneous xenograft in nude mice. RNA-Seq and qRT-PCR were used to identify genes and pathways involved. Mass spectrometry and co-immunoprecipitation were used to explore the interaction of proteins.

RESULTS

Trim24 is highly expressed in spontaneous NB in TH-MYCN transgenic mice and clinical NB specimens. It is associated with poor NB differentiation and unfavorable prognostic. Knockout of TRIM24 in neuroblastoma cells promotes cell differentiation, reduces cell stemness, and inhibits colony formation in soft agar and subcutaneous xenograft tumor growth in nude mice. Mechanistically, TRIM24 knockout alters genes and pathways related to neural differentiation and development by suppressing LSD1/CoREST complex formation. Besides, TRIM24 knockout activates the retinoic acid pathway. Targeting TRIM24 in combination with retinoic acid (RA) synergistically promotes NB cell differentiation and inhibits cell viability.

CONCLUSION

Our findings demonstrate that TRIM24 is critical for NB differentiation and suggest that TRIM24 is a promising therapeutic target in combination with RA in NB differentiation therapy.

A fluorescent reporter model for the visualization and characterization of TDC

TDC are hematopoietic cells that combine dendritic cell (DC) and conventional T-cell markers and functional properties. They were identified in secondary lymphoid organs (SLOs) of naïve mice as cells expressing CD11c, major histocompatibility molecules (MHC)-II, and the T-cell receptor (TCR). Despite thorough characterization, a physiological role for TDC remains to be determined. Unfortunately, using CD11c as a marker for TDC has the caveat of its upregulation on different cells, including T cells, upon activation. Here, we took advantage of Zbtb46-GFP reporter mice to explore the frequency and localization of TDC in different tissues at steady state and upon viral infection. RNA sequencing analysis confirmed that TDC sorted from Zbtb46-GFP mice have a gene signature that is distinct from conventional T cells and DC. In addition, this reporter model allowed for identification of TDC in situ not only in SLOs but also in the liver and lung of naïve mice. Interestingly, we found that TDC numbers in the SLOs increased upon viral infection, suggesting that TDC might play a role during viral infections. In conclusion, we propose a visualization strategy that might shed light on the physiological role of TDC in several pathological contexts, including infection and cancer.

Transcriptional Profiling of SARS-CoV-2-Infected Calu-3 Cells Reveals Immune-Related Signaling Pathways

The COVID-19 disease, caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), emerged in late 2019 and rapidly spread worldwide, becoming a pandemic that infected millions of people and caused significant deaths. COVID-19 continues to be a major threat, and there is a need to deepen our understanding of the virus and its mechanisms of infection. To study the cellular responses to SARS-CoV-2 infection, we performed an RNA sequencing of infected vs. uninfected Calu-3 cells. Total RNA was extracted from infected (0.5 MOI) and control Calu-3 cells and converted to cDNA. Sequencing was performed, and the obtained reads were quality-analyzed and pre-processed. Differential expression was assessed with the EdgeR package, and functional enrichment was performed in EnrichR for Gene Ontology, KEGG pathways, and WikiPathways. A total of 1040 differentially expressed genes were found in infected vs. uninfected Calu-3 cells, of which 695 were up-regulated and 345 were down-regulated. Functional enrichment analyses revealed the predominant up-regulation of genes related to innate immune response, response to virus, inflammation, cell proliferation, and apoptosis. These transcriptional changes following SARS-CoV-2 infection may reflect a cellular response to the infection and help to elucidate COVID-19 pathogenesis, in addition to revealing potential biomarkers and drug targets.

Dex-Benchmark: datasets and code to evaluate algorithms for transcriptomics data analysis

Many tools and algorithms are available for analyzing transcriptomics data. These include algorithms for performing sequence alignment, data normalization and imputation, clustering, identifying differentially expressed genes, and performing gene set enrichment analysis. To make the best choice about which tools to use, objective benchmarks can be developed to compare the quality of different algorithms to extract biological knowledge maximally and accurately from these data. The Dexamethasone Benchmark (Dex-Benchmark) resource aims to fill this need by providing the community with datasets and code templates for benchmarking different gene expression analysis tools and algorithms. The resource provides access to a collection of curated RNA-seq, L1000, and ChIP-seq data from dexamethasone treatment as well as genetic perturbations of its known targets. In addition, the website provides Jupyter Notebooks that use these pre-processed curated datasets to demonstrate how to benchmark the different steps in gene expression analysis. By comparing two independent data sources and data types with some expected concordance, we can assess which tools and algorithms best recover such associations. To demonstrate the usefulness of the resource for discovering novel drug targets, we applied it to optimize data processing strategies for the chemical perturbations and CRISPR single gene knockouts from the L1000 transcriptomics data from the Library of Integrated Network Cellular Signatures (LINCS) program, with a focus on understudied proteins from the Illuminating the Druggable Genome (IDG) program. Overall, the Dex-Benchmark resource can be utilized to assess the quality of transcriptomics and other related bioinformatics data analysis workflows. The resource is available from: https://maayanlab.github.io/dex-benchmark.

The Quartet Data Portal: integration of community-wide resources for multiomics quality control

The Quartet Data Portal facilitates community access to well-characterized reference materials, reference datasets, and related resources established based on a family of four individuals with identical twins from the Quartet Project. Users can request DNA, RNA, protein, and metabolite reference materials, as well as datasets generated across omics, platforms, labs, protocols, and batches. Reproducible analysis tools allow for objective performance assessment of user-submitted data, while interactive visualization tools support rapid exploration of reference datasets. A closed-loop "distribution-collection-evaluation-integration" workflow enables updates and integration of community-contributed multiomics data. Ultimately, this portal helps promote the advancement of reference datasets and multiomics quality control.

Cerebrospinal fluid immune cells appear similar across neuropathic and non-neuropathic pain conditions

Background

Microglia have been implicated in the pathophysiology of neuropathic pain. Here, we sought to investigate whether cerebrospinal fluid (CSF) might be used as a proxy-measure of microglial activation in human participants.

Methods

We preformed fluorescence-activated cell sorting (FACS) of CSF immune cell populations derived from individuals who experienced pain with neuropathic features. We sorted CD4+, CD8+ T cells and monocytes and analyzed their transcriptome using RNA sequencing. We also performed Cellular Indexing of Transcriptomes and Epitopes (CITE) sequencing to characterize the expression of all CSF immune cells in a patient with postherpetic neuralgia and in a patient with neuropathic pain after failed back surgery.

Results

Immune cell numbers and phenotypes were not obviously different between individuals regardless of the etiology of their pain. This was true when examining our own dataset, as well as when comparing it to previously published single-cell RNA sequencing data of human CSF. In all instances, CSF monocytes showed expression of myeloid cell markers commonly associated with microglia ( P2RY12, TMEM119 and OLFML3), which will make it difficult to ascertain the origin of CSF proteins: do they derive directly from circulating CSF monocytes or could some originate in spinal cord microglia in the parenchyma?

Conclusions

We conclude that it will not be straightforward to use CSF as a biomarker for microglial function in humans.

miR-124 and VAMP3 Act Antagonistically in Human Neuroblastoma

Neuroblastoma (NB) is the most common extracranial solid tumor that affects developing nerve cells in the fetus, infants, and children. miR-124 is a microRNA (miRNA) enriched in neuronal tissues, and VAMP3 (vesicle-associated membrane protein 3) has been reported to be an miR-124 target, although the relationship between NB and miR-124 or VAMP3 is unknown. Our current work identified that miR-124 levels are high in NB cases and that elevated miR-124 correlates with worse NB outcomes. Conversely, depressed VAMP3 correlates with worse NB outcomes. To investigate the mechanisms by which miR-124 and VAMP3 regulate NB, we altered miR-124 or VAMP3 expression in human NB cells and observed that increased miR-124 and reduced VAMP3 stimulated cell proliferation and suppressed apoptosis, while increased VAMP3 had the opposite effects. Genome-wide mRNA expression analyses identified gene and pathway changes which might explain the NB cell phenotypes. Together, our studies suggest that miR-124 and VAMP3 could be potential new markers of NB and targets of NB treatments.

Merging short and stranded long reads improves transcript assembly

Long-read RNA sequencing has arisen as a counterpart to short-read sequencing, with the potential to capture full-length isoforms, albeit at the cost of lower depth. Yet this potential is not fully realized due to inherent limitations of current long-read assembly methods and underdeveloped approaches to integrate short-read data. Here, we critically compare the existing methods and develop a new integrative approach to characterize a particularly challenging pool of low-abundance long noncoding RNA (lncRNA) transcripts from short- and long-read sequencing in two distinct cell lines. Our analysis reveals severe limitations in each of the sequencing platforms. For short-read assemblies, coverage declines at transcript termini resulting in ambiguous ends, and uneven low coverage results in segmentation of a single transcript into multiple transcripts. Conversely, long-read sequencing libraries lack depth and strand-of-origin information in cDNA-based methods, culminating in erroneous assembly and quantitation of transcripts. We also discover a cDNA synthesis artifact in long-read datasets that markedly impacts the identity and quantitation of assembled transcripts. Towards remediating these problems, we develop a computational pipeline to "strand" long-read cDNA libraries that rectifies inaccurate mapping and assembly of long-read transcripts. Leveraging the strengths of each platform and our computational stranding, we also present and benchmark a hybrid assembly approach that drastically increases the sensitivity and accuracy of full-length transcript assembly on the correct strand and improves detection of biological features of the transcriptome. When applied to a challenging set of under-annotated and cell-type variable lncRNA, our method resolves the segmentation problem of short-read sequencing and the depth problem of long-read sequencing, resulting in the assembly of coherent transcripts with precise 5' and 3' ends. Our workflow can be applied to existing datasets for superior demarcation of transcript ends and refined isoform structure, which can enable better differential gene expression analyses and molecular manipulations of transcripts.

Shifts in isoform usage underlie transcriptional differences in regulatory T cells in type 1 diabetes

Genome-wide association studies have identified numerous loci with allelic associations to Type 1 Diabetes (T1D) risk. Most disease-associated variants are enriched in regulatory sequences active in lymphoid cell types, suggesting that lymphocyte gene expression is altered in T1D. Here we assay gene expression between T1D cases and healthy controls in two autoimmunity-relevant lymphocyte cell types, memory CD4+/CD25+ regulatory T cells (Treg) and memory CD4+/CD25- T cells, using a splicing event-based approach to characterize tissue-specific transcriptomes. Limited differences in isoform usage between T1D cases and controls are observed in memory CD4+/CD25- T-cells. In Tregs, 402 genes demonstrate differences in isoform usage between cases and controls, particularly RNA recognition and splicing factor genes. Many of these genes are regulated by the variable inclusion of exons that can trigger nonsense mediated decay. Our results suggest that dysregulation of gene expression, through shifts in alternative splicing in Tregs, contributes to T1D pathophysiology.

A fast and globally optimal solution for RNA-seq quantification

Alignment-based RNA-seq quantification methods typically involve a time-consuming alignment process prior to estimating transcript abundances. In contrast, alignment-free RNA-seq quantification methods bypass this step, resulting in significant speed improvements. Existing alignment-free methods rely on the Expectation-Maximization (EM) algorithm for estimating transcript abundances. However, EM algorithms only guarantee locally optimal solutions, leaving room for further accuracy improvement by finding a globally optimal solution. In this study, we present TQSLE, the first alignment-free RNA-seq quantification method that provides a globally optimal solution for transcript abundances estimation. TQSLE adopts a two-step approach: first, it constructs a k-mer frequency matrix A for the reference transcriptome and a k-mer frequency vector b for the RNA-seq reads; then, it directly estimates transcript abundances by solving the linear equation ATAx = ATb. We evaluated the performance of TQSLE using simulated and real RNA-seq data sets and observed that, despite comparable speed to other alignment-free methods, TQSLE outperforms them in terms of accuracy. TQSLE is freely available at https://github.com/yhg926/TQSLE.

Target Genes of c-MYC and MYCN with Prognostic Power in Neuroblastoma Exhibit Different Expressions during Sympathoadrenal Development

Deregulation of the MYC family of transcription factors c-MYC (encoded by MYC), MYCN, and MYCL is prevalent in most human cancers, with an impact on tumor initiation and progression, as well as response to therapy. In neuroblastoma (NB), amplification of the MYCN oncogene and over-expression of MYC characterize approximately 40% and 10% of all high-risk NB cases, respectively. However, the mechanism and stage of neural crest development in which MYCN and c-MYC contribute to the onset and/or progression of NB are not yet fully understood. Here, we hypothesized that subtle differences in the expression of MYCN and/or c-MYC targets could more accurately stratify NB patients in different risk groups rather than using the expression of either MYC gene alone. We employed an integrative approach using the transcriptome of 498 NB patients from the SEQC cohort and previously defined c-MYC and MYCN target genes to model a multigene transcriptional risk score. Our findings demonstrate that defined sets of c-MYC and MYCN targets with significant prognostic value, effectively stratify NB patients into different groups with varying overall survival probabilities. In particular, patients exhibiting a high-risk signature score present unfavorable clinical parameters, including increased clinical risk, higher INSS stage, MYCN amplification, and disease progression. Notably, target genes with prognostic value differ between c-MYC and MYCN, exhibiting distinct expression patterns in the developing sympathoadrenal system. Genes associated with poor outcomes are mainly found in sympathoblasts rather than in chromaffin cells during the sympathoadrenal development.

Analysis of High-Risk Neuroblastoma Transcriptome Reveals Gene Co-Expression Signatures and Functional Features

Neuroblastoma represents a neoplastic expansion of neural crest cells in the developing sympathetic nervous system and is childhood's most common extracranial solid tumor. The heterogeneity of gene expression in different types of cancer is well-documented, and genetic features of neuroblastoma have been described by classification, development stage, malignancy, and progression of tumors. Here, we aim to analyze RNA sequencing datasets, publicly available in the GDC data portal, of neuroblastoma tumor samples from various patients and compare them with normal adrenal gland tissue from the GTEx data portal to elucidate the gene expression profile and regulation networks they share. Our results from the differential expression, weighted correlation network, and functional enrichment analyses that we performed with the count data from neuroblastoma and standard normal gland samples indicate that the analysis of transcriptome data from 58 patients diagnosed with high-risk neuroblastoma shares the expression pattern of 104 genes. More importantly, our analyses identify the co-expression relationship and the role of these genes in multiple biological processes and signaling pathways strongly associated with this disease phenotype. Our approach proposes a group of genes and their biological functions to be further investigated as essential molecules and possible therapeutic targets of neuroblastoma regardless of the etiology of individual tumors.

Quartet protein reference materials and datasets for multi-platform assessment of label-free proteomics

BACKGROUND

Quantitative proteomics is an indispensable tool in life science research. However, there is a lack of reference materials for evaluating the reproducibility of label-free liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based measurements among different instruments and laboratories.

RESULTS

Here, we develop the Quartet standard as a proteome reference material with built-in truths, and distribute the same aliquots to 15 laboratories with nine conventional LC-MS/MS platforms across six cities in China. Relative abundance of over 12,000 proteins on 816 mass spectrometry files are obtained and compared for reproducibility among the instruments and laboratories to ultimately generate proteomics benchmark datasets. There is a wide dynamic range of proteomes spanning about 7 orders of magnitude, and the injection order has marked effects on quantitative instead of qualitative characteristics.

CONCLUSION

Overall, the Quartet offers valuable standard materials and data resources for improving the quality control of proteomic analyses as well as the reproducibility and reliability of research findings.

Correcting batch effects in large-scale multiomics studies using a reference-material-based ratio method

BACKGROUND

Batch effects are notoriously common technical variations in multiomics data and may result in misleading outcomes if uncorrected or over-corrected. A plethora of batch-effect correction algorithms are proposed to facilitate data integration. However, their respective advantages and limitations are not adequately assessed in terms of omics types, the performance metrics, and the application scenarios.

RESULTS

As part of the Quartet Project for quality control and data integration of multiomics profiling, we comprehensively assess the performance of seven batch effect correction algorithms based on different performance metrics of clinical relevance, i.e., the accuracy of identifying differentially expressed features, the robustness of predictive models, and the ability of accurately clustering cross-batch samples into their own donors. The ratio-based method, i.e., by scaling absolute feature values of study samples relative to those of concurrently profiled reference material(s), is found to be much more effective and broadly applicable than others, especially when batch effects are completely confounded with biological factors of study interests. We further provide practical guidelines for implementing the ratio based approach in increasingly large-scale multiomics studies.

CONCLUSIONS

Multiomics measurements are prone to batch effects, which can be effectively corrected using ratio-based scaling of the multiomics data. Our study lays the foundation for eliminating batch effects at a ratio scale.

Allele-specific quantification of human leukocyte antigen transcript isoforms by nanopore sequencing

Introduction

While tens of thousands of HLA alleles have been identified by DNA sequencing, the contribution of alternative splicing to HLA diversity is not well characterized. In this study, we sought to determine if long-read sequencing could be used to accurately quantify allele-specific HLA transcripts in primary human lymphocytes.

Methods

cDNA libraries were prepared from peripheral blood lymphocytes from 12 donors and sequenced by nanopore long-read sequencing. HLA reads were aligned to donor-specific reference sequences based on the known type of each donor. Allele-specific exon utilization was calculated as the proportion of reads aligning to each allele containing known exons, and transcript isoforms were quantified based on patterns of exon utilization within individual reads.

Results

Splice variants were rare among class I HLA genes (median exon retention rate 99%-100%), except for several HLA-C alleles with exon 5 spliced out of up to 15% of reads. Splice variants were also rare among class II HLA genes (median exon retention rate 98%-100%), except for HLA-DQB1. Consistent with previous work, exon 5 of HLA-DQB1 was spliced out in alleles with a mutated splice acceptor site at rs28688207. Surprisingly, a 28% loss of exon 5 was also observed in HLA-DQB1 alleles with an intact splice acceptor site at rs28688207.

Discussion

We describe a simple bioinformatic workflow to quantify allele-specific expression of HLA transcript isoforms. Further studies are warranted to characterize the repertoire of HLA transcripts expressed in different cell types and tissues across diverse populations.