Quantitative gene set analysis generalized for repeated measures, confounder adjustment, and continuous covariates

Tracking early lung cancer metastatic dissemination in TRACERx using ctDNA

Circulating tumour DNA (ctDNA) can be used to detect and profile residual tumour cells persisting after curative intent therapy1. The study of large patient cohorts incorporating longitudinal plasma sampling and extended follow-up is required to determine the role of ctDNA as a phylogenetic biomarker of relapse in early-stage non-small-cell lung cancer (NSCLC). Here we developed ctDNA methods tracking a median of 200 mutations identified in resected NSCLC tissue across 1,069 plasma samples collected from 197 patients enrolled in the TRACERx study2. A lack of preoperative ctDNA detection distinguished biologically indolent lung adenocarcinoma with good clinical outcome. Postoperative plasma analyses were interpreted within the context of standard-of-care radiological surveillance and administration of cytotoxic adjuvant therapy. Landmark analyses of plasma samples collected within 120 days after surgery revealed ctDNA detection in 25% of patients, including 49% of all patients who experienced clinical relapse; 3 to 6 monthly ctDNA surveillance identified impending disease relapse in an additional 20% of landmark-negative patients. We developed a bioinformatic tool (ECLIPSE) for non-invasive tracking of subclonal architecture at low ctDNA levels. ECLIPSE identified patients with polyclonal metastatic dissemination, which was associated with a poor clinical outcome. By measuring subclone cancer cell fractions in preoperative plasma, we found that subclones seeding future metastases were significantly more expanded compared with non-metastatic subclones. Our findings will support (neo)adjuvant trial advances and provide insights into the process of metastatic dissemination using low-ctDNA-level liquid biopsy.

Molecular determinants of response to PD-L1 blockade across tumor types

Immune checkpoint inhibitors targeting the PD-1/PD-L1 axis lead to durable clinical responses in subsets of cancer patients across multiple indications, including non-small cell lung cancer (NSCLC), urothelial carcinoma (UC) and renal cell carcinoma (RCC). Herein, we complement PD-L1 immunohistochemistry (IHC) and tumor mutation burden (TMB) with RNA-seq in 366 patients to identify unifying and indication-specific molecular profiles that can predict response to checkpoint blockade across these tumor types. Multiple machine learning approaches failed to identify a baseline transcriptional signature highly predictive of response across these indications. Signatures described previously for immune checkpoint inhibitors also failed to validate. At the pathway level, significant heterogeneity is observed between indications, in particular within the PD-L1⁺ tumors. mUC and NSCLC are molecularly aligned, with cell cycle and DNA damage repair genes associated with response in PD-L1- tumors. At the gene level, the CDK4/6 inhibitor CDKN2A is identified as a significant transcriptional correlate of response, highlighting the association of non-immune pathways to the outcome of checkpoint blockade. This cross-indication analysis reveals molecular heterogeneity between mUC, NSCLC and RCC tumors, suggesting that indication-specific molecular approaches should be prioritized to formulate treatment strategies.

Longitudinal linear combination test for gene set analysis

Background

Although microarray studies have greatly contributed to recent genetic advances, lack of replication has been a continuing concern in this area. Complex study designs have the potential to address this concern, though they remain undervalued by investigators due to the lack of proper analysis methods. The primary challenge in the analysis of complex microarray study data is handling the correlation structure within data while also dealing with the combination of large number of genetic measurements and small number of subjects that are ubiquitous even in standard microarray studies. Motivated by the lack of available methods for analysis of repeatedly measured phenotypic or transcriptomic data, herein we develop a longitudinal linear combination test (LLCT).

Results

LLCT is a two-step method to analyze multiple longitudinal phenotypes when there is high dimensionality in response and/or explanatory variables. Alternating between calculating within-subjects and between-subjects variations in two steps, LLCT examines if the maximum possible correlation between a linear combination of the time trends and a linear combination of the predictors given by the gene expressions is statistically significant. A generalization of this method can handle family-based study designs when the subjects are not independent. This method is also applicable to time-course microarray, with the ability to identify gene sets that exhibit significantly different expression patterns over time. Based on the results from a simulation study, LLCT outperformed its alternative: pathway analysis via regression. LLCT was shown to be very powerful in the analysis of large gene sets even when the sample size is small.

Conclusions

This self-contained pathway analysis method is applicable to a wide range of longitudinal genomics, proteomics, metabolomics (OMICS) data, allows adjusting for potentially time-dependent covariates and works well with unbalanced and incomplete data. An important potential application of this method could be time-course linkage of OMICS, an attractive possibility for future genetic researchers.

Availability: R package of LLCT is available at: https://github.com/its-likeli-jeff/LLCT

Longitudinal profiling of human blood transcriptome in healthy and lupus pregnancy

Healthy and uncomplicated lupus pregnancies exhibit early and sustained transcriptional modulation of lupus-related pathways. This might contribute to fetal tolerance while predisposing pregnant women to certain infections. Failure to modulate these pathways is associated with lupus pregnancy complications.

Systemic lupus erythematosus carries an increased risk of pregnancy complications, including preeclampsia and fetal adverse outcomes. To identify the underlying molecular mechanisms, we longitudinally profiled the blood transcriptome of 92 lupus patients and 43 healthy women during pregnancy and postpartum and performed multicolor flow cytometry in a subset of them. We also profiled 25 healthy women undergoing assisted reproductive technology to monitor transcriptional changes around embryo implantation. Sustained down-regulation of multiple immune signatures, including interferon and plasma cells, was observed during healthy pregnancy. These changes appeared early after embryo implantation and were mirrored in uncomplicated lupus pregnancies. Patients with preeclampsia displayed early up-regulation of neutrophil signatures that correlated with expansion of immature neutrophils. Lupus pregnancies with fetal complications carried the highest interferon and plasma cell signatures as well as activated CD4⁺ T cell counts. Thus, blood immunomonitoring reveals that both healthy and uncomplicated lupus pregnancies exhibit early and sustained transcriptional modulation of lupus-related signatures, and a lack thereof associates with adverse outcomes.

Gene set meta-analysis with Quantitative Set Analysis for Gene Expression (QuSAGE)

Small sample sizes combined with high person-to-person variability can make it difficult to detect significant gene expression changes from transcriptional profiling studies. Subtle, but coordinated, gene expression changes may be detected using gene set analysis approaches. Meta-analysis is another approach to increase the power to detect biologically relevant changes by integrating information from multiple studies. Here, we present a framework that combines both approaches and allows for meta-analysis of gene sets. QuSAGE meta-analysis extends our previously published QuSAGE framework, which offers several advantages for gene set analysis, including fully accounting for gene-gene correlations and quantifying gene set activity as a full probability density function. Application of QuSAGE meta-analysis to influenza vaccination response shows it can detect significant activity that is not apparent in individual studies.

Phantom: investigating heterogeneous gene sets in time-course data

Motivation

Gene set analysis is a powerful tool to study the coordinative change of time-course data. However, most existing methods only model the overall change of a gene set, yet completely overlooked heterogeneous time-dependent changes within sub-sets of genes.

Results

We have developed a novel statistical method, Phantom, to investigate gene set heterogeneity. Phantom employs the principle of multi-objective optimization to assess the heterogeneity inside a gene set, which also accounts for the temporal dependency in time-course data. Phantom improves the performance of gene set based methods to detect biological changes across time.

Availability and implementation

Phantom webpage can be accessed at: http://www.baylorhealth.edu/Phantom . R package of Phantom is available at https://cran.r-project.org/web/packages/phantom/index.html .

Contact

jinghua.gu@bswhealth.org.

Supplementary information

Supplementary data are available at Bioinformatics online.

Signatures of inflammation and impending multiple organ dysfunction in the hyperacute phase of trauma: A prospective cohort study

BACKGROUND

Severe trauma induces a widespread response of the immune system. This "genomic storm" can lead to poor outcomes, including Multiple Organ Dysfunction Syndrome (MODS). MODS carries a high mortality and morbidity rate and adversely affects long-term health outcomes. Contemporary management of MODS is entirely supportive, and no specific therapeutics have been shown to be effective in reducing incidence or severity. The pathogenesis of MODS remains unclear, and several models are proposed, such as excessive inflammation, a second-hit insult, or an imbalance between pro- and anti-inflammatory pathways. We postulated that the hyperacute window after trauma may hold the key to understanding how the genomic storm is initiated and may lead to a new understanding of the pathogenesis of MODS.

METHODS AND FINDINGS

We performed whole blood transcriptome and flow cytometry analyses on a total of 70 critically injured patients (Injury Severity Score [ISS] ≥ 25) at The Royal London Hospital in the hyperacute time period within 2 hours of injury. We compared transcriptome findings in 36 critically injured patients with those of 6 patients with minor injuries (ISS ≤ 4). We then performed flow cytometry analyses in 34 critically injured patients and compared findings with those of 9 healthy volunteers. Immediately after injury, only 1,239 gene transcripts (4%) were differentially expressed in critically injured patients. By 24 hours after injury, 6,294 transcripts (21%) were differentially expressed compared to the hyperacute window. Only 202 (16%) genes differentially expressed in the hyperacute window were still expressed in the same direction at 24 hours postinjury. Pathway analysis showed principally up-regulation of pattern recognition and innate inflammatory pathways, with down-regulation of adaptive responses. Immune deconvolution, flow cytometry, and modular analysis suggested a central role for neutrophils and Natural Killer (NK) cells, with underexpression of T- and B cell responses. In the transcriptome cohort, 20 critically injured patients later developed MODS. Compared with the 16 patients who did not develop MODS (NoMODS), maximal differential expression was seen within the hyperacute window. In MODS versus NoMODS, 363 genes were differentially expressed on admission, compared to only 33 at 24 hours postinjury. MODS transcripts differentially expressed in the hyperacute window showed enrichment among diseases and biological functions associated with cell survival and organismal death rather than inflammatory pathways. There was differential up-regulation of NK cell signalling pathways and markers in patients who would later develop MODS, with down-regulation of neutrophil deconvolution markers. This study is limited by its sample size, precluding more detailed analyses of drivers of the hyperacute response and different MODS phenotypes, and requires validation in other critically injured cohorts.

CONCLUSIONS

In this study, we showed how the hyperacute postinjury time window contained a focused, specific signature of the response to critical injury that led to widespread genomic activation. A transcriptomic signature for later development of MODS was present in this hyperacute window; it showed a strong signal for cell death and survival pathways and implicated NK cells and neutrophil populations in this differential response.

Personalized Immunomonitoring Uncovers Molecular Networks that Stratify Lupus Patients

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by loss of tolerance to nucleic acids and highly diverse clinical manifestations. To assess its molecular heterogeneity, we longitudinally profiled the blood transcriptome of 158 pediatric patients. Using mixed models accounting for repeated measurements, demographics, treatment, disease activity (DA), and nephritis class, we confirmed a prevalent IFN signature and identified a plasmablast signature as the most robust biomarker of DA. We detected gradual enrichment of neutrophil transcripts during progression to active nephritis and distinct signatures in response to treatment in different nephritis subclasses. Importantly, personalized immunomonitoring uncovered individual correlates of disease activity that enabled patient stratification into seven groups, supported by patient genotypes. Our study uncovers the molecular heterogeneity of SLE and provides an explanation for the failure of clinical trials. This approach may improve trial design and implementation of tailored therapies in genetically and clinically complex autoimmune diseases. PAPERCLIP.