Covariance-based analyses of biological pathways

Distinct spatial immune microlandscapes are independently associated with outcomes in triple-negative breast cancer

The utility of spatial immunobiomarker quantitation in prognostication and therapeutic prediction is actively being investigated in triple-negative breast cancer (TNBC). Here, with high-plex quantitative digital spatial profiling, we map and quantitate intraepithelial and adjacent stromal tumor immune protein microenvironments in systemic treatment-naïve (female only) TNBC to assess the spatial context in immunobiomarker-based prediction of outcome. Immune protein profiles of CD45-rich and CD68-rich stromal microenvironments differ significantly. While they typically mirror adjacent, intraepithelial microenvironments, this is not uniformly true. In two TNBC cohorts, intraepithelial CD40 or HLA-DR enrichment associates with better outcomes, independently of stromal immune protein profiles or stromal TILs and other established prognostic variables. In contrast, intraepithelial or stromal microenvironment enrichment with IDO1 associates with improved survival irrespective of its spatial location. Antigen-presenting and T-cell activation states are inferred from eigenprotein scores. Such scores within the intraepithelial compartment interact with PD-L1 and IDO1 in ways that suggest prognostic and/or therapeutic potential. This characterization of the intrinsic spatial immunobiology of treatment-naïve TNBC highlights the importance of spatial microenvironments for biomarker quantitation to resolve intrinsic prognostic and predictive immune features and ultimately inform therapeutic strategies for clinically actionable immune biomarkers.

Homogeneity tests of covariance matrices with high-dimensional longitudinal data

This paper deals with the detection and identification of changepoints among covariances of high-dimensional longitudinal data, where the number of features is greater than both the sample size and the number of repeated measurements. The proposed methods are applicable under general temporal-spatial dependence. A new test statistic is introduced for changepoint detection, and its asymptotic distribution is established. If a changepoint is detected, an estimate of the location is provided. The rate of convergence of the estimator is shown to depend on the data dimension, sample size, and signal-to-noise ratio. Binary segmentation is used to estimate the locations of possibly multiple changepoints, and the corresponding estimator is shown to be consistent under mild conditions. Simulation studies provide the empirical size and power of the proposed test and the accuracy of the changepoint estimator. An application to a time-course microarray dataset identifies gene sets with significant gene interaction changes over time.

Stepwise approach to SNP-set analysis illustrated with the Metabochip and colorectal cancer in Japanese Americans of the Multiethnic Cohort

BACKGROUND

Common variants have explained less than the amount of heritability expected for complex diseases, which has led to interest in less-common variants and more powerful approaches to the analysis of whole-genome scans. Because of low frequency (low statistical power), less-common variants are best analyzed using SNP-set methods such as gene-set or pathway-based analyses. However, there is as yet no clear consensus regarding how to focus in on potential risk variants following set-based analyses. We used a stepwise, telescoping approach to analyze common- and rare-variant data from the Illumina Metabochip array to assess genomic association with colorectal cancer (CRC) in the Japanese sub-population of the Multiethnic Cohort (676 cases, 7180 controls). We started with pathway analysis of SNPs that are in genes and pathways having known mechanistic roles in colorectal cancer, then focused on genes within the pathways that evidenced association with CRC, and finally assessed individual SNPs within the genes that evidenced association. Pathway SNPs downloaded from the dbSNP database were cross-matched with Metabochip SNPs and analyzed using the logistic kernel machine regression approach (logistic SNP-set kernel-machine association test, or sequence kernel association test; SKAT) and related methods.

RESULTS

The TGF-β and WNT pathways were associated with all CRC, and the WNT pathway was associated with colon cancer. Individual genes demonstrating the strongest associations were TGFBR2 in the TGF-β pathway and SMAD7 (which is involved in both the TGF-β and WNT pathways). As partial validation of our approach, a known CRC risk variant in SMAD7 (in both the TGF-β and WNT pathways: rs11874392) was associated with CRC risk in our data. We also detected two novel candidate CRC risk variants (rs13075948 and rs17025857) in TGFBR2, a gene known to be associated with CRC risk.

CONCLUSIONS

A stepwise, telescoping approach identified some potentially novel risk variants associated with colorectal cancer, so it may be a useful method for following up on results of set-based SNP analyses. Further work is required to assess the statistical characteristics of the approach, and additional applications should aid in better clarifying its utility.