Integrating multiomics and prior knowledge: a study of the Graphnet penalty impact.

MOTIVATION: In the field of oncology, statistical models are used for the discovery of candidate factors that influence the development of the pathology or its outcome. These statistical models can be designed in a multiblock framework to study the relationship between different multiomic data, and variable selection is often achieved by imposing constraints on the model parameters. A priori graph constraints have been used in the literature as a way to improve feature selection in the model, yielding more interpretability. However, it is still unclear how these graphs interact with the models and how they impact the feature selection. Additionally, with the availability of different graphs encoding different information, one can wonder how the choice of the graph meaningfully impacts the results obtained. RESULTS: We proposed to study the graph penalty impact on a multiblock model. Specifically, we used the SGCCA as the multiblock framework. We studied the effect of the penalty on the model using the TCGA-LGG dataset. Our findings are 3-fold. We showed that the graph penalty increases the number of selected genes from this dataset, while selecting genes already identified in other works as pertinent biomarkers in the pathology. We demonstrated that using different graphs leads to different though consistent results, but that graph density is the main factor influencing the obtained results. Finally, we showed that the graph penalty increases the performance of the survival prediction from the model-derived components and the interpretability of the results. AVAILABILITY AND IMPLEMENTATION: Source code is freely available at https://github.com/neurospin/netSGCCA.

Multiway generalized canonical correlation analysis.

Regularized generalized canonical correlation analysis (RGCCA) is a general multiblock data analysis framework that encompasses several important multivariate analysis methods such as principal component analysis, partial least squares regression, and several versions of generalized canonical correlation analysis. In this article, we extend RGCCA to the case where at least one block has a tensor structure. This method is called multiway generalized canonical correlation analysis (MGCCA). Convergence properties of the MGCCA algorithm are studied, and computation of higher-level components are discussed. The usefulness of MGCCA is shown on simulation and on the analysis of a cognitive study in human infants using electroencephalography (EEG).

Cortical signatures in behaviorally clustered autistic traits subgroups: a population-based study.

Extensive heterogeneity in autism spectrum disorder (ASD) has hindered the characterization of consistent biomarkers, which has led to widespread negative results. Isolating homogenized subtypes could provide insight into underlying biological mechanisms and an overall better understanding of ASD. A total of 1093 participants from the population-based "Healthy Brain Network" cohort (Child Mind Institute in the New York City area, USA) were selected based on score availability in behaviors relevant to ASD, aged 6-18 and IQ >= 70. All participants underwent an unsupervised clustering analysis on behavioral dimensions to reveal subgroups with ASD traits, identified by the presence of social deficits. Analysis revealed three socially impaired ASD traits subgroups: (1) high in emotionally dysfunctional traits, (2) high in ADHD-like traits, and (3) high in anxiety and depressive symptoms. 527 subjects had good quality structural MRI T1 data. Site effects on cortical features were adjusted using the ComBat method. Neuroimaging analyses compared cortical thickness, gyrification, and surface area, and were controlled for age, gender, and IQ, and corrected for multiple comparisons. Structural neuroimaging analyses contrasting one combined heterogeneous ASD traits group against controls did not yield any significant differences. Unique cortical signatures, however, were observed within each of the three individual ASD traits subgroups versus controls. These observations provide evidence of ASD traits subtypes, and confirm the necessity of applying dimensional approaches to extract meaningful differences, thus reducing heterogeneity and paving the way to better understanding ASD traits.

Clinical, molecular, and radiomic profile of gliomas with FGFR3-TACC3 fusions.

BACKGROUND: Actionable fibroblast growth factor receptor 3 (FGFR3)-transforming acidic coiled-coil protein 3 fusions (F3T3) are found in approximately 3% of gliomas, but their characteristics and prognostic significance are still poorly defined. Our goal was to characterize the clinical, radiological, and molecular profile of F3T3 positive diffuse gliomas. METHODS: We screened F3T3 fusion by real-time (RT)-PCR and FGFR3 immunohistochemistry in a large series of gliomas, characterized for main genetic alterations, histology, and clinical evolution. We performed a radiological and radiomic case control study, using an exploratory and a validation cohort. RESULTS: We screened 1162 diffuse gliomas (951 unselected cases and 211 preselected for FGFR3 protein immunopositivity), identifying 80 F3T3 positive gliomas. F3T3 was mutually exclusive with IDH mutation (P < 0.001) and EGFR amplification (P = 0.01), defining a distinct molecular cluster associated with CDK4 (P = 0.04) and MDM2 amplification (P = 0.03). F3T3 fusion was associated with longer survival for the whole series and for glioblastomas (median overall survival was 31.1 vs 19.9 mo, P = 0.02) and was an independent predictor of better outcome on multivariate analysis.F3T3 positive gliomas had specific MRI features, affecting preferentially insula and temporal lobe, and with poorly defined tumor margins. F3T3 fusion was correctly predicted by radiomics analysis on both the exploratory (area under the curve [AUC] = 0.87) and the validation MRI (AUC = 0.75) cohort. Using Cox proportional hazards models, radiomics predicted survival with a high C-index (0.75, SD 0.04), while the model combining clinical, genetic, and radiomic data showed the highest C-index (0.81, SD 0.04). CONCLUSION: F3T3 positive gliomas have distinct molecular and radiological features, and better outcome.

A constrained singular value decomposition method that integrates sparsity and orthogonality.

We propose a new sparsification method for the singular value decomposition-called the constrained singular value decomposition (CSVD)-that can incorporate multiple constraints such as sparsification and orthogonality for the left and right singular vectors. The CSVD can combine different constraints because it implements each constraint as a projection onto a convex set, and because it integrates these constraints as projections onto the intersection of multiple convex sets. We show that, with appropriate sparsification constants, the algorithm is guaranteed to converge to a stable point. We also propose and analyze the convergence of an efficient algorithm for the specific case of the projection onto the balls defined by the norms L1 and L2. We illustrate the CSVD and compare it to the standard singular value decomposition and to a non-orthogonal related sparsification method with: 1) a simulated example, 2) a small set of face images (corresponding to a configuration with a number of variables much larger than the number of observations), and 3) a psychometric application with a large number of observations and a small number of variables. The companion R-package, csvd, that implements the algorithms described in this paper, along with reproducible examples, are available for download from https://github.com/vguillemot/csvd.