Identification of biomarkers predictive of metastasis development in early-stage colorectal cancer using network-based regularization.

Colorectal cancer (CRC) is the third most common cancer and the second most deathly worldwide. It is a very heterogeneous disease that can develop via distinct pathways where metastasis is the primary cause of death. Therefore, it is crucial to understand the molecular mechanisms underlying metastasis. RNA-sequencing is an essential tool used for studying the transcriptional landscape. However, the high-dimensionality of gene expression data makes selecting novel metastatic biomarkers problematic. To distinguish early-stage CRC patients at risk of developing metastasis from those that are not, three types of binary classification approaches were used: (1) classification methods (decision trees, linear and radial kernel support vector machines, logistic regression, and random forest) using differentially expressed genes (DEGs) as input features; (2) regularized logistic regression based on the Elastic Net penalty and the proposed iTwiner-a network-based regularizer accounting for gene correlation information; and (3) classification methods based on the genes pre-selected using regularized logistic regression. Classifiers using the DEGs as features showed similar results, with random forest showing the highest accuracy. Using regularized logistic regression on the full dataset yielded no improvement in the methods' accuracy. Further classification using the pre-selected genes found by different penalty factors, instead of the DEGs, significantly improved the accuracy of the binary classifiers. Moreover, the use of network-based correlation information (iTwiner) for gene selection produced the best classification results and the identification of more stable and robust gene sets. Some are known to be tumor suppressor genes (OPCML-IT2), to be related to resistance to cancer therapies (RAC1P3), or to be involved in several cancer processes such as genome stability (XRCC6P2), tumor growth and metastasis (MIR602) and regulation of gene transcription (NME2P2). We show that the classification of CRC patients based on pre-selected features by regularized logistic regression is a valuable alternative to using DEGs, significantly increasing the models' predictive performance. Moreover, the use of correlation-based penalization for biomarker selection stands as a promising strategy for predicting patients' groups based on RNA-seq data.

Plasma cytokine profiling in sibling pairs discordant for autism spectrum disorder.

OBJECTIVE: Converging lines of evidence point to the existence of immune dysfunction in autism spectrum disorder (ASD), which could directly affect several key neurodevelopmental processes. Previous studies have shown higher cytokine levels in patients with autism compared with matched controls or subjects with other developmental disorders. In the current study, we used plasma-cytokine profiling for 25 discordant sibling pairs to evaluate whether these alterations occur within families with ASD. METHODS: Plasma-cytokine profiling was conducted using an array-based multiplex sandwich ELISA for simultaneous quantitative measurement of 40 unique targets. We also analyzed the correlations between cytokine levels and clinically relevant quantitative traits (Vineland Adaptive Behavior Scale in Autism (VABS) composite score, Social Responsiveness Scale (SRS) total T score, head circumference, and full intelligence quotient (IQ)). In addition, because of the high phenotypic heterogeneity of ASD, we defined four subgroups of subjects (those who were non-verbal, those with gastrointestinal issues, those with regressive autism, and those with a history of allergies), which encompass common and/or recurrent endophenotypes in ASD, and tested the cytokine levels in each group. RESULTS: None of the measured parameters showed significant differences between children with ASD and their related typically developing siblings. However, specific target levels did correlate with quantitative clinical traits, and these were significantly different when the ASD subgroups were analyzed. It is notable that these differences seem to be attributable to a predisposing immunogenetic background, as no other significant differences were noticed between discordant sibling pairs. Interleukin-1ß appears to be the cytokine most involved in quantitative traits and clinical subgroups of ASD. CONCLUSIONS: In the present study, we found a lack of significant differences in plasma-cytokine levels between children with ASD and in their related non-autistic siblings. Thus, our results support the evidence that the immune profiles of children with autism do not differ from their typically developing siblings. However, the significant association of cytokine levels with the quantitative traits and the clinical subgroups analyzed suggests that altered immune responses may affect core feature of ASD.

Clinical-genomic determinants of immune checkpoint blockade response in head and neck squamous cell carcinoma.

BACKGROUNDRecurrent and/or metastatic (R/M) head and neck squamous cell carcinoma (HNSCC) is generally an incurable disease, with patients experiencing median survival of under 10 months and significant morbidity. While immune checkpoint blockade (ICB) drugs are effective in approximately 20% of patients, the remaining experience limited clinical benefit and are exposed to potential adverse effects and financial costs. Clinically approved biomarkers, such as tumor mutational burden (TMB), have a modest predictive value in HNSCC.METHODSWe analyzed clinical and genomic features, generated using whole-exome sequencing, in 133 ICB-treated patients with R/M HNSCC, of whom 69 had virus-associated and 64 had non-virus-associated tumors.RESULTSHierarchical clustering of genomic data revealed 6 molecular subtypes characterized by a wide range of objective response rates and survival after ICB therapy. The prognostic importance of these 6 subtypes was validated in an external cohort. A random forest-based predictive model, using several clinical and genomic features, predicted progression-free survival (PFS), overall survival (OS), and response with greater accuracy than did a model based on TMB alone. Recursive partitioning analysis identified 3 features (systemic inflammatory response index, TMB, and smoking signature) that classified patients into risk groups with accurate discrimination of PFS and OS.CONCLUSIONThese findings shed light on the immunogenomic characteristics of HNSCC tumors that drive differential responses to ICB and identify a clinical-genomic classifier that outperformed the current clinically approved biomarker of TMB. This validated predictive tool may help with clinical risk stratification in patients with R/M HNSCC for whom ICB is being considered.FUNDINGFundación Alfonso Martín Escudero, NIH R01 DE027738, US Department of Defense CA210784, The Geoffrey Beene Cancer Research Center, The MSKCC Population Science Research Program, the Jayme Flowers Fund, the Sebastian Nativo Fund, and the NIH/NCI Cancer Center Support Grant P30 CA008748.

Genetic and microenvironmental intra-tumor heterogeneity impacts colorectal cancer evolution and metastatic development.

Colorectal cancer (CRC) is a highly diverse disease, where different genomic instability pathways shape genetic clonal diversity and tumor microenvironment. Although intra-tumor heterogeneity has been characterized in primary tumors, its origin and consequences in CRC outcome is not fully understood. Therefore, we assessed intra- and inter-tumor heterogeneity of a prospective cohort of 136 CRC samples. We demonstrate that CRC diversity is forged by asynchronous forms of molecular alterations, where mutational and chromosomal instability collectively boost CRC genetic and microenvironment intra-tumor heterogeneity. We were able to depict predictor signatures of cancer-related genes that can foresee heterogeneity levels across the different tumor consensus molecular subtypes (CMS) and primary tumor location. Finally, we show that high genetic and microenvironment heterogeneity are associated with lower metastatic potential, whereas late-emerging copy number variations favor metastasis development and polyclonal seeding. This study provides an exhaustive portrait of the interplay between genetic and microenvironment intra-tumor heterogeneity across CMS subtypes, depicting molecular events with predictive value of CRC progression and metastasis development.

Spatiotemporal evolution of the clear cell renal cell carcinoma microenvironment links intra-tumoral heterogeneity to immune escape.

BACKGROUND: Intratumoral heterogeneity (ITH) is a hallmark of clear cell renal cell carcinoma (ccRCC) that reflects the trajectory of evolution and influences clinical prognosis. Here, we seek to elucidate how ITH and tumor evolution during immune checkpoint inhibitor (ICI) treatment can lead to therapy resistance. METHODS: Here, we completed a single-arm pilot study to examine the safety and feasibility of neoadjuvant nivolumab in patients with localized RCC. Primary endpoints were safety and feasibility of neoadjuvant nivolumab. Then, we spatiotemporally profiled the genomic and immunophenotypic characteristics of 29 ccRCC patients, including pre- and post-therapy samples from 17 ICI-treated patients. Deep multi-regional whole-exome and transcriptome sequencing were performed on 29 patients at different time points before and after ICI therapy. T cell repertoire was also monitored from tissue and peripheral blood collected from a subset of patients to study T cell clonal expansion during ICI therapy. RESULTS: Angiogenesis, lymphocytic infiltration, and myeloid infiltration varied significantly across regions of the same patient, potentially confounding their utility as biomarkers of ICI response. Elevated ITH associated with a constellation of both genomic features (HLA LOH, CDKN2A/B loss) and microenvironmental features, including elevated myeloid expression, reduced peripheral T cell receptor (TCR) diversity, and putative neoantigen depletion. Hypothesizing that ITH may itself play a role in shaping ICI response, we derived a transcriptomic signature associated with neoantigen depletion that strongly associated with response to ICI and targeted therapy treatment in several independent clinical trial cohorts. CONCLUSIONS: These results argue that genetic and immune heterogeneity jointly co-evolve and influence response to ICI in ccRCC. Our findings have implications for future biomarker development for ICI response across ccRCC and other solid tumors and highlight important features of tumor evolution under ICI treatment. TRIAL REGISTRATION: The study was registered on ClinicalTrial.gov (NCT02595918) on November 4, 2015.

Identification of genetic variants using bar-coded multiplexed sequencing.

We developed a generalized framework for multiplexed resequencing of targeted human genome regions on the Illumina Genome Analyzer using degenerate indexed DNA bar codes ligated to fragmented DNA before sequencing. Using this method, we simultaneously sequenced the DNA of multiple HapMap individuals at several Encyclopedia of DNA Elements (ENCODE) regions. We then evaluated the use of Bayes factors for discovering and genotyping polymorphisms. For polymorphisms that were either previously identified within the Single Nucleotide Polymorphism database (dbSNP) or visually evident upon re-inspection of archived ENCODE traces, we observed a false positive rate of 11.3% using strict thresholds for predicting variants and 69.6% for lax thresholds. Conversely, false negative rates were 10.8-90.8%, with false negatives at stricter cut-offs occurring at lower coverage (<10 aligned reads). These results suggest that >90% of genetic variants are discoverable using multiplexed sequencing provided sufficient coverage at the polymorphic base.

HERVs establish a distinct molecular subtype in stage II/III colorectal cancer with poor outcome.

Colorectal cancer (CRC) is one of the most lethal malignancies. The extreme heterogeneity in survival rate is driving the need for new prognostic biomarkers. Human endogenous retroviruses (hERVs) have been suggested to influence tumor progression, oncogenesis and elicit an immune response. We examined multiple next-generation sequencing (NGS)-derived biomarkers in 114 CRC patients with paired whole-exome and whole-transcriptome sequencing (WES and WTS, respectively). First, we demonstrate that the median expression of hERVs can serve as a potential biomarker for prognosis, relapse, and resistance to chemotherapy in stage II and III CRC. We show that hERV expression and CD8+ tumor-infiltrating T-lymphocytes (TILs) synergistically stratify overall and relapse-free survival (OS and RFS): the median OS of the CD8-/hERV+ subgroup was 29.8 months compared with 37.5 months for other subgroups (HR = 4.4, log-rank P < 0.001). Combing NGS-based biomarkers (hERV/CD8 status) with clinicopathological factors provided a better prediction of patient survival compared to clinicopathological factors alone. Moreover, we explored the association between genomic and transcriptomic features of tumors with high hERV expression and establish this subtype as distinct from previously described consensus molecular subtypes of CRC. Overall, our results underscore a previously unknown role for hERVs in leading to a more aggressive subtype of CRC.

Single-cell sequencing links multiregional immune landscapes and tissue-resident T cells in ccRCC to tumor topology and therapy efficacy.

Clear cell renal cell carcinomas (ccRCCs) are highly immune infiltrated, but the effect of immune heterogeneity on clinical outcome in ccRCC has not been fully characterized. Here we perform paired single-cell RNA (scRNA) and T cell receptor (TCR) sequencing of 167,283 cells from multiple tumor regions, lymph node, normal kidney, and peripheral blood of two immune checkpoint blockade (ICB)-naïve and four ICB-treated patients to map the ccRCC immune landscape. We detect extensive heterogeneity within and between patients, with enrichment of CD8A+ tissue-resident T cells in a patient responsive to ICB and tumor-associated macrophages (TAMs) in a resistant patient. A TCR trajectory framework suggests distinct T cell differentiation pathways between patients responding and resistant to ICB. Finally, scRNA-derived signatures of tissue-resident T cells and TAMs are associated with response to ICB and targeted therapies across multiple independent cohorts. Our study establishes a multimodal interrogation of the cellular programs underlying therapeutic efficacy in ccRCC.

Candidate Agtr2 influenced genes and pathways identified by expression profiling in the developing brain of Agtr2(-/y) mice.

Intellectual disability (ID) is a common developmental disability observed in 1 to 3% of the human population. A possible role for the Angiotensin II type 2 receptor (AGTR2) in brain function, affecting learning, memory, and behavior, has been suggested in humans and rodents. Mice lacking the Agtr2 gene (Agtr2(-/y)) showed significant impairment in their spatial memory and exhibited abnormal dendritic spine morphology. To identify Agtr2 influenced genes and pathways, we performed whole genome microarray analysis on RNA isolated from brains of Agtr2(-/y) and control male mice at embryonic day 15 (E15) and postnatal day one (P1). The gene expression profiles of the Agtr2(-/y) brain samples were significantly different when compared to profiles of the age-matched control brains. We identified 62 differently expressed genes (p< or =0.005) at E15 and in P1 brains of the Agtr2(-/y) mice. We verified the differential expression of several of these genes in brain samples using quantitative RT-PCR. Differentially expressed genes encode molecules involved in multiple cellular processes including microtubule functions associated with dendritic spine morphology. This study provides insight into Agtr2 influenced candidate genes and suggests that expression dysregulation of these genes may modulate Agtr2 actions in the brain that influences learning and memory.

AT2 Receptor-Interacting Proteins ATIPs in the Brain.

A complete renin-angiotensin system (RAS) is locally expressed in the brain and fulfills important functions. Angiotensin II, the major biologically active peptide of the RAS, acts via binding to two main receptor subtypes designated AT1 and AT2. The present paper focuses on AT2 receptors, which have been reported to have neuroprotective effects on stroke, degenerative diseases, and cognitive functions. Our group has identified a family of AT2 receptor interacting proteins (ATIPs) comprising three major members (ATIP1, ATIP3, and ATIP4) with different intracellular localization. Of interest, all ATIP members are expressed in brain tissues and carry a conserved domain able to interact with the AT2 receptor intracellular tail, suggesting a role in AT2-mediated brain functions. We summarize here current knowledge on the ATIP family of proteins, and we present new experimental evidence showing interaction defects between ATIP1 and two mutant forms of the AT2 receptor identified in cases of mental retardation. These studies point to a functional role of the AT2/ATIP1 axis in cognition.

Identification of a common variant affecting human episodic memory performance using a pooled genome-wide association approach: a case study of disease gene identification.

Genome-wide association studies (GWAS) are an important tool for discovering novel genes associated with disease or traits. Careful design of case-control groups greatly facilitates the efficacy of these studies. Here we describe a pooled GWAS study undertaken to find novel genes associated with human episodic memory performance. A genomic locus for the WW and C2 domain-containing 1 protein, KIBRA (also known as WWC1), was found to be associated with memory performance in three cognitively normal cohorts from Switzerland and the USA. This result was further supported by correlation of KIBRA genotype and differences in hippocampal activation as measured by functional magnetic resonance imaging (fMRI). These findings provide an excellent example of the application of GWAS using a pooled genomic DNA approach to successfully identify a locus with strong effects on human memory.

Hippocampal gene expression changes during age-related cognitive decline.

As humans age, cognitive performance decreases differentially across individuals. This age-related decline in otherwise healthy individuals is likely due to the interaction of multiple factors including genetics and environment. We hypothesized that altered spatial memory performance in genetically similar mice could be in part due to differential gene expression patterns in the hippocampus. To investigate this we utilized Morris water maze (MWM) testing in a group of young (3 months) and aged (24 months) C57BL/J male mice. Two sub-groups were identified in the aged animals; one in which MWM performance was not significantly different when compared to the young animals (aged-unimpaired; "AU") and one in which performance was significantly different by 1.5 standard deviations from the mean (aged-impaired; "AI"). One week after testing was completed the entire hippocampus was collected from six each of AU, AI and young mice and their gene expression profiles were compared using Affymetrix microarrays. Benjamini and Hochberg FDR correction at p<0.05 identified 18 genes differentially expressed between the AI and AU mice. The correlation between behavioral deficits and gene expression patterning allows a better understanding of how altered gene expression in the hippocampus contributes to accelerated age-related cognitive decline and delineates between gene expression changes associated with normal aging vs. memory performance.