Search

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., (© The Author(s) 2023. Published by Oxford University Press.)

Grey matter damage has been established as a key contributor to disability progression in multiple sclerosis. Aside from neuronal loss and axonal transections, which predominate in cortical demyelinated lesions, synaptic alterations have been detected in both demyelinated plaques and normal-appearing grey matter, resulting in functional neuronal damage. The axon initial segment is a key element of neuronal function, responsible for action potential initiation and maintenance of neuronal polarity. Despite several reports of profound axon initial segment alterations in different pathological models, among which experimental auto-immune encephalomyelitis, whether the axon initial segment is affected in multiple sclerosis is still unknown. Using immunohistochemistry, we analysed axon initial segments from control and multiple sclerosis tissue, focusing on layer 5/6 pyramidal neurons in the neocortex and Purkinje cells in the cerebellum and performed analysis on the parameters known to control neuronal excitability, i.e. axon initial segment length and position. We found that the axon initial segment length was increased only in pyramidal neurons of inactive demyelinated lesions, compared with normal appearing grey matter tissue. In contrast, in both cell types, the axon initial segment position was altered, with an increased soma-axon initial segment gap, in both active and inactive demyelinated lesions. In addition, using a computational model, we show that this increased gap between soma and axon initial segment might increase neuronal excitability. Taken together, these results show, for the first time, changes of axon initial segments in multiple sclerosis, in active as well as inactive grey matter lesions in both neocortex and cerebellum, which might alter neuronal function., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Guarantors of Brain.)

Genome-wide association study (GWAS) has been the driving force for identifying association between genetic variants and human phenotypes. Thousands of GWAS summary statistics covering a broad range of human traits and diseases are now publicly available. These GWAS have proven their utility for a range of secondary analyses, including in particular the joint analysis of multiple phenotypes to identify new associated genetic variants. However, although several methods have been proposed, there are very few large-scale applications published so far because of challenges in implementing these methods on real data. Here, we present JASS (Joint Analysis of Summary Statistics), a polyvalent Python package that addresses this need. Our package incorporates recently developed joint tests such as the omnibus approach and various weighted sum of Z -score tests while solving all practical and computational barriers for large-scale multivariate analysis of GWAS summary statistics. This includes data cleaning and harmonization tools, an efficient algorithm for fast derivation of joint statistics, an optimized data management process and a web interface for exploration purposes. Both benchmark analyses and real data applications demonstrated the robustness and strong potential of JASS for the detection of new associated genetic variants. Our package is freely available at https://gitlab.pasteur.fr/statistical-genetics/jass., (© The Author(s) 2019. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

High-throughput measurements of molecular phenotypes provide an unprecedented opportunity to model cellular processes and their impact on disease. These highly structured datasets are usually strongly confounded, creating false positives and reducing power. This has motivated many approaches based on principal components analysis (PCA) to estimate and correct for confounders, which have become indispensable elements of association tests between molecular phenotypes and both genetic and nongenetic factors. Here, we show that these correction approaches induce a bias, and that it persists for large sample sizes and replicates out-of-sample. We prove this theoretically for PCA by deriving an analytic, deterministic, and intuitive bias approximation. We assess other methods with realistic simulations, which show that perturbing any of several basic parameters can cause false positive rate (FPR) inflation. Our experiments show the bias depends on covariate and confounder sparsity, effect sizes, and their correlation. Surprisingly, when the covariate and confounder have [Formula: see text], standard two-step methods all have [Formula: see text]-fold FPR inflation. Our analysis informs best practices for confounder correction in genomic studies, and suggests many false discoveries have been made and replicated in some differential expression analyses., (Copyright © 2019 by the Genetics Society of America.)

The growing number of modalities (e.g. multi-omics, imaging and clinical data) characterizing a given disease provides physicians and statisticians with complementary facets reflecting the disease process but emphasizes the need for novel statistical methods of data analysis able to unify these views. Such data sets are indeed intrinsically structured in blocks, where each block represents a set of variables observed on a group of individuals. Therefore, classical statistical tools cannot be applied without altering their organization, with the risk of information loss. Regularized generalized canonical correlation analysis (RGCCA) and its sparse generalized canonical correlation analysis (SGCCA) counterpart are component-based methods for exploratory analyses of data sets structured in blocks of variables. Rather than operating sequentially on parts of the measurements, the RGCCA/SGCCA-based integrative analysis method aims at summarizing the relevant information between and within the blocks. It processes a priori information defining which blocks are supposed to be linked to one another, thus reflecting hypotheses about the biology underlying the data blocks. It also requires the setting of extra parameters that need to be carefully adjusted.Here, we provide practical guidelines for the use of RGCCA/SGCCA. We also illustrate the flexibility and usefulness of RGCCA/SGCCA on a unique cohort of patients with four genetic subtypes of spinocerebellar ataxia, in which we obtained multiple data sets from brain volumetry and magnetic resonance spectroscopy, and metabolomic and lipidomic analyses. As a first step toward the extraction of multimodal biomarkers, and through the reduction to a few meaningful components and the visualization of relevant variables, we identified possible markers of disease progression.

One major challenge in multiple sclerosis is to understand the cellular and molecular mechanisms leading to disease severity progression. The recently demonstrated correlation between disease severity and remyelination emphasizes the importance of identifying factors leading to a favourable outcome. Why remyelination fails or succeeds in multiple sclerosis patients remains largely unknown, mainly because remyelination has never been studied within a humanized pathological context that would recapitulate major events in plaque formation such as infiltration of inflammatory cells. Therefore, we developed a new paradigm by grafting healthy donor or multiple sclerosis patient lymphocytes in the demyelinated lesion of nude mice spinal cord. We show that lymphocytes play a major role in remyelination whose efficacy is significantly decreased in mice grafted with multiple sclerosis lymphocytes compared to those grafted with healthy donors lymphocytes. Mechanistically, we demonstrated in vitro that lymphocyte-derived mediators influenced differentiation of oligodendrocyte precursor cells through a crosstalk with microglial cells. Among mice grafted with lymphocytes from different patients, we observed diverse remyelination patterns reproducing for the first time the heterogeneity observed in multiple sclerosis patients. Comparing lymphocyte secretory profile from patients exhibiting high and low remyelination ability, we identified novel molecules involved in oligodendrocyte precursor cell differentiation and validated CCL19 as a target to improve remyelination. Specifically, exogenous CCL19 abolished oligodendrocyte precursor cell differentiation observed in patients with high remyelination pattern. Multiple sclerosis lymphocytes exhibit intrinsic capacities to coordinate myelin repair and further investigation on patients with high remyelination capacities will provide new pro-regenerative strategies., (© The Author (2017). Published by Oxford University Press on behalf of the Guarantors of Brain.)

Regularized generalized canonical correlation analysis (RGCCA) is a generalization of regularized canonical correlation analysis to 3 or more sets of variables. RGCCA is a component-based approach which aims to study the relationships between several sets of variables. The quality and interpretability of the RGCCA components are likely to be affected by the usefulness and relevance of the variables in each block. Therefore, it is an important issue to identify within each block which subsets of significant variables are active in the relationships between blocks. In this paper, RGCCA is extended to address the issue of variable selection. Specifically, sparse generalized canonical correlation analysis (SGCCA) is proposed to combine RGCCA with an [Formula: see text]-penalty in a unified framework. Within this framework, blocks are not necessarily fully connected, which makes SGCCA a flexible method for analyzing a wide variety of practical problems. Finally, the versatility and usefulness of SGCCA are illustrated on a simulated dataset and on a 3-block dataset which combine gene expression, comparative genomic hybridization, and a qualitative phenotype measured on a set of 53 children with glioma. SGCCA is available on CRAN as part of the RGCCA package., (© The Author 2014. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Motivation: In statistical bioinformatics research, different optimization mechanisms potentially lead to 'over-optimism' in published papers. So far, however, a systematic critical study concerning the various sources underlying this over-optimism is lacking., Results: We present an empirical study on over-optimism using high-dimensional classification as example. Specifically, we consider a 'promising' new classification algorithm, namely linear discriminant analysis incorporating prior knowledge on gene functional groups through an appropriate shrinkage of the within-group covariance matrix. While this approach yields poor results in terms of error rate, we quantitatively demonstrate that it can artificially seem superior to existing approaches if we 'fish for significance'. The investigated sources of over-optimism include the optimization of datasets, of settings, of competing methods and, most importantly, of the method's characteristics. We conclude that, if the improvement of a quantitative criterion such as the error rate is the main contribution of a paper, the superiority of new algorithms should always be demonstrated on independent validation data., Availability: The R codes and relevant data can be downloaded from http://www.ibe.med.uni-muenchen.de/organisation/mitarbeiter/020_professuren/boulesteix/overoptimism/, such that the study is completely reproducible.