A computational framework to improve cross-platform implementation of transcriptomics signatures.

The emergence of next-generation sequencing technologies and computational advances have expanded our understanding of gene expression regulation (i.e., the transcriptome). This has also led to an increased interest in using transcriptomic biomarkers to improve disease diagnosis and stratification, to assess prognosis and predict the response to treatment. Significant progress in identifying transcriptomic signatures for various clinical needs has been made, with large discovery studies accounting for challenges such as patient variability, unwanted batch effects, and data complexities; however, obstacles related to the technical aspects of cross-platform implementation still hinder the successful integration of transcriptomic technologies into standard diagnostic workflows. In this article, we discuss the challenges associated with integrating transcriptomic signatures derived using high-throughput technologies (such as RNA-sequencing) into clinical diagnostic tools using nucleic acid amplification (NAA) techniques. The novelty of the proposed approach lies in our aim to embed constraints related to cross-platform implementation in the process of signature discovery. These constraints could include technical limitations of amplification platform and chemistry, the maximal number of targets imposed by the chosen multiplexing strategy, and the genomic context of identified RNA biomarkers. Finally, we propose to build a computational framework that would integrate these constraints in combination with existing statistical and machine learning models used for signature identification. We envision that this could accelerate the integration of RNA signatures discovered by high-throughput technologies into NAA-based approaches suitable for clinical applications.

Analysis of blood and nasal epithelial transcriptomes to identify mechanisms associated with control of SARS-CoV-2 viral load in the upper respiratory tract.

OBJECTIVES: The amount of SARS-CoV-2 detected in the upper respiratory tract (URT viral load) is a key driver of transmission of infection. Current evidence suggests that mechanisms constraining URT viral load are different from those controlling lower respiratory tract viral load and disease severity. Understanding such mechanisms may help to develop treatments and vaccine strategies to reduce transmission. Combining mathematical modelling of URT viral load dynamics with transcriptome analyses we aimed to identify mechanisms controlling URT viral load. METHODS: COVID-19 patients were recruited in Spain during the first wave of the pandemic. RNA sequencing of peripheral blood and targeted NanoString nCounter transcriptome analysis of nasal epithelium were performed and gene expression analysed in relation to paired URT viral load samples collected within 15 days of symptom onset. Proportions of major immune cells in blood were estimated from transcriptional data using computational differential estimation. Weighted correlation network analysis (adjusted for cell proportions) and fixed transcriptional repertoire analysis were used to identify associations with URT viral load, quantified as standard deviations (z-scores) from an expected trajectory over time. RESULTS: Eighty-two subjects (50% female, median age 54 years (range 3-73)) with COVID-19 were recruited. Paired URT viral load samples were available for 16 blood transcriptome samples, and 17 respiratory epithelial transcriptome samples. Natural Killer (NK) cells were the only blood cell type significantly correlated with URT viral load z-scores (r = -0.62, P = 0.010). Twenty-four blood gene expression modules were significantly correlated with URT viral load z-score, the most significant being a module of genes connected around IFNA14 (Interferon Alpha-14) expression (r = -0.60, P = 1e-10). In fixed repertoire analysis, prostanoid-related gene expression was significantly associated with higher viral load. In nasal epithelium, only GNLY (granulysin) gene expression showed significant negative correlation with viral load. CONCLUSIONS: Correlations between the transcriptional host response and inter-individual variations in SARS-CoV-2 URT viral load, revealed many molecular mechanisms plausibly favouring or constraining viral replication. Existing evidence corroborates many of these mechanisms, including likely roles for NK cells, granulysin, prostanoids and interferon alpha-14. Inhibition of prostanoid production and administration of interferon alpha-14 may be attractive transmission-blocking interventions.

Diagnosis of Multisystem Inflammatory Syndrome in Children by a Whole-Blood Transcriptional Signature.

BACKGROUND: To identify a diagnostic blood transcriptomic signature that distinguishes multisystem inflammatory syndrome in children (MIS-C) from Kawasaki disease (KD), bacterial infections, and viral infections. METHODS: Children presenting with MIS-C to participating hospitals in the United Kingdom and the European Union between April 2020 and April 2021 were prospectively recruited. Whole-blood RNA Sequencing was performed, contrasting the transcriptomes of children with MIS-C (n = 38) to those from children with KD (n = 136), definite bacterial (DB; n = 188) and viral infections (DV; n = 138). Genes significantly differentially expressed (SDE) between MIS-C and comparator groups were identified. Feature selection was used to identify genes that optimally distinguish MIS-C from other diseases, which were subsequently translated into RT-qPCR assays and evaluated in an independent validation set comprising MIS-C (n = 37), KD (n = 19), DB (n = 56), DV (n = 43), and COVID-19 (n = 39). RESULTS: In the discovery set, 5696 genes were SDE between MIS-C and combined comparator disease groups. Five genes were identified as potential MIS-C diagnostic biomarkers (HSPBAP1, VPS37C, TGFB1, MX2, and TRBV11-2), achieving an AUC of 96.8% (95% CI: 94.6%-98.9%) in the discovery set, and were translated into RT-qPCR assays. The RT-qPCR 5-gene signature achieved an AUC of 93.2% (95% CI: 88.3%-97.7%) in the independent validation set when distinguishing MIS-C from KD, DB, and DV. CONCLUSIONS: MIS-C can be distinguished from KD, DB, and DV groups using a 5-gene blood RNA expression signature. The small number of genes in the signature and good performance in both discovery and validation sets should enable the development of a diagnostic test for MIS-C.

Immunoglobulin, glucocorticoid, or combination therapy for multisystem inflammatory syndrome in children: a propensity-weighted cohort study.

Background: Multisystem inflammatory syndrome in children (MIS-C), a hyperinflammatory condition associated with SARS-CoV-2 infection, has emerged as a serious illness in children worldwide. Immunoglobulin or glucocorticoids, or both, are currently recommended treatments. Methods: The Best Available Treatment Study evaluated immunomodulatory treatments for MIS-C in an international observational cohort. Analysis of the first 614 patients was previously reported. In this propensity-weighted cohort study, clinical and outcome data from children with suspected or proven MIS-C were collected onto a web-based Research Electronic Data Capture database. After excluding neonates and incomplete or duplicate records, inverse probability weighting was used to compare primary treatments with intravenous immunoglobulin, intravenous immunoglobulin plus glucocorticoids, or glucocorticoids alone, using intravenous immunoglobulin as the reference treatment. Primary outcomes were a composite of inotropic or ventilator support from the second day after treatment initiation, or death, and time to improvement on an ordinal clinical severity scale. Secondary outcomes included treatment escalation, clinical deterioration, fever, and coronary artery aneurysm occurrence and resolution. This study is registered with the ISRCTN registry, ISRCTN69546370. Findings: We enrolled 2101 children (aged 0 months to 19 years) with clinically diagnosed MIS-C from 39 countries between June 14, 2020, and April 25, 2022, and, following exclusions, 2009 patients were included for analysis (median age 8·0 years [IQR 4·2-11·4], 1191 [59·3%] male and 818 [40·7%] female, and 825 [41·1%] White). 680 (33·8%) patients received primary treatment with intravenous immunoglobulin, 698 (34·7%) with intravenous immunoglobulin plus glucocorticoids, 487 (24·2%) with glucocorticoids alone; 59 (2·9%) patients received other combinations, including biologicals, and 85 (4·2%) patients received no immunomodulators. There were no significant differences between treatments for primary outcomes for the 1586 patients with complete baseline and outcome data that were considered for primary analysis. Adjusted odds ratios for ventilation, inotropic support, or death were 1·09 (95% CI 0·75-1·58; corrected p value=1·00) for intravenous immunoglobulin plus glucocorticoids and 0·93 (0·58-1·47; corrected p value=1·00) for glucocorticoids alone, versus intravenous immunoglobulin alone. Adjusted average hazard ratios for time to improvement were 1·04 (95% CI 0·91-1·20; corrected p value=1·00) for intravenous immunoglobulin plus glucocorticoids, and 0·84 (0·70-1·00; corrected p value=0·22) for glucocorticoids alone, versus intravenous immunoglobulin alone. Treatment escalation was less frequent for intravenous immunoglobulin plus glucocorticoids (OR 0·15 [95% CI 0·11-0·20]; p<0·0001) and glucocorticoids alone (0·68 [0·50-0·93]; p=0·014) versus intravenous immunoglobulin alone. Persistent fever (from day 2 onward) was less common with intravenous immunoglobulin plus glucocorticoids compared with either intravenous immunoglobulin alone (OR 0·50 [95% CI 0·38-0·67]; p<0·0001) or glucocorticoids alone (0·63 [0·45-0·88]; p=0·0058). Coronary artery aneurysm occurrence and resolution did not differ significantly between treatment groups. Interpretation: Recovery rates, including occurrence and resolution of coronary artery aneurysms, were similar for primary treatment with intravenous immunoglobulin when compared to glucocorticoids or intravenous immunoglobulin plus glucocorticoids. Initial treatment with glucocorticoids appears to be a safe alternative to immunoglobulin or combined therapy, and might be advantageous in view of the cost and limited availability of intravenous immunoglobulin in many countries. Funding: Imperial College London, the European Union's Horizon 2020, Wellcome Trust, the Medical Research Foundation, UK National Institute for Health and Care Research, and National Institutes of Health.

Next-generation molecular diagnostics: Leveraging digital technologies to enhance multiplexing in real-time PCR.

Real-time polymerase chain reaction (qPCR) enables accurate detection and quantification of nucleic acids and has become a fundamental tool in biological sciences, bioengineering and medicine. By combining multiple primer sets in one reaction, it is possible to detect several DNA or RNA targets simultaneously, a process called multiplex PCR (mPCR) which is key to attaining optimal throughput, cost-effectiveness and efficiency in molecular diagnostics, particularly in infectious diseases. Multiple solutions have been devised to increase multiplexing in qPCR, including single-well techniques, using target-specific fluorescent oligonucleotide probes, and spatial multiplexing, where segregation of the sample enables parallel amplification of multiple targets. However, these solutions are mostly limited to three or four targets, or highly sophisticated and expensive instrumentation. There is a need for innovations that will push forward the multiplexing field in qPCR, enabling for a next generation of diagnostic tools which could accommodate high throughput in an affordable manner. To this end, the use of machine learning (ML) algorithms (data-driven solutions) has recently emerged to leverage information contained in amplification and melting curves (AC and MC, respectively) - two of the most standard bio-signals emitted during qPCR - for accurate classification of multiple nucleic acid targets in a single reaction. Therefore, this review aims to demonstrate and illustrate that data-driven solutions can be successfully coupled with state-of-the-art and common qPCR platforms using a variety of amplification chemistries to enhance multiplexing in qPCR. Further, because both ACs and MCs can be predicted from sequence data using thermodynamic databases, it has also become possible to use computer simulation to rationalize and optimize the design of mPCR assays where target detection is supported by data-driven technologies. Thus, this review also discusses recent work converging towards the development of an end-to-end framework where knowledge-based and data-driven software solutions are integrated to streamline assay design, and increase the accuracy of target detection and quantification in the multiplex setting. We envision that concerted efforts by academic and industry scientists will help advance these technologies, to a point where they become mature and robust enough to bring about major improvements in the detection of nucleic acids across many fields.

Characterisation of the blood RNA host response underpinning severity in COVID-19 patients.

Infection with SARS-CoV-2 has highly variable clinical manifestations, ranging from asymptomatic infection through to life-threatening disease. Host whole blood transcriptomics can offer unique insights into the biological processes underpinning infection and disease, as well as severity. We performed whole blood RNA Sequencing of individuals with varying degrees of COVID-19 severity. We used differential expression analysis and pathway enrichment analysis to explore how the blood transcriptome differs between individuals with mild, moderate, and severe COVID-19, performing pairwise comparisons between groups. Increasing COVID-19 severity was characterised by an abundance of inflammatory immune response genes and pathways, including many related to neutrophils and macrophages, in addition to an upregulation of immunoglobulin genes. In this study, for the first time, we show how immunomodulatory treatments commonly administered to COVID-19 patients greatly alter the transcriptome. Our insights into COVID-19 severity reveal the role of immune dysregulation in the progression to severe disease and highlight the need for further research exploring the interplay between SARS-CoV-2 and the inflammatory immune response.