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.

Alterations of the axon initial segment in multiple sclerosis grey matter.

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.

High prevalence of small intestine bacteria overgrowth and asymptomatic carriage of enteric pathogens in stunted children in Antananarivo, Madagascar.

Environmental Enteric Dysfunction (EED) refers to an incompletely defined syndrome of inflammation, reduced absorptive capacity, and reduced barrier function in the small intestine. It is widespread among children and adults in low- and middle-income countries and is also associated with poor sanitation and certain gut infections possibly resulting in an abnormal gut microbiota, small intestinal bacterial overgrowth (SIBO) and stunting. We investigated bacterial pathogen exposure in stunted and non-stunted children in Antananarivo, Madagascar by collecting fecal samples from 464 children (96 severely stunted, 104 moderately stunted and 264 non-stunted) and the prevalence of SIBO in 109 duodenal aspirates from stunted children (61 from severely stunted and 48 from moderately stunted children). SIBO assessed by both aerobic and anaerobic plating techniques was very high: 85.3% when selecting a threshold of ≥105 CFU/ml of bacteria in the upper intestinal aspirates. Moreover, 58.7% of the children showed more than 106 bacteria/ml in these aspirates. The most prevalent cultivated genera recovered were Streptococcus, Neisseria, Staphylococcus, Rothia, Haemophilus, Pantoea and Branhamella. Feces screening by qPCR showed a high prevalence of bacterial enteropathogens, especially those categorized as being enteroinvasive or causing mucosal disruption, such as Shigella spp., enterotoxigenic Escherichia coli, enteropathogenic E. coli and enteroaggregative E. coli. These pathogens were detected at a similar rate in stunted children and controls, all showing no sign of severe diarrhea the day of inclusion but both living in a highly contaminated environment (slum-dwelling). Interestingly Shigella spp. was the most prevalent enteropathogen found in this study (83.3%) without overrepresentation in stunted children.

The genetic architecture of language functional connectivity.

Language is a unique trait of the human species, of which the genetic architecture remains largely unknown. Through language disorders studies, many candidate genes were identified. However, such complex and multifactorial trait is unlikely to be driven by only few genes and case-control studies, suffering from a lack of power, struggle to uncover significant variants. In parallel, neuroimaging has significantly contributed to the understanding of structural and functional aspects of language in the human brain and the recent availability of large scale cohorts like UK Biobank have made possible to study language via image-derived endophenotypes in the general population. Because of its strong relationship with task-based fMRI (tbfMRI) activations and its easiness of acquisition, resting-state functional MRI (rsfMRI) have been more popularised, making it a good surrogate of functional neuronal processes. Taking advantage of such a synergistic system by aggregating effects across spatially distributed traits, we performed a multivariate genome-wide association study (mvGWAS) between genetic variations and resting-state functional connectivity (FC) of classical brain language areas in the inferior frontal (pars opercularis, triangularis and orbitalis), temporal and inferior parietal lobes (angular and supramarginal gyri), in 32,186 participants from UK Biobank. Twenty genomic loci were found associated with language FCs, out of which three were replicated in an independent replication sample. A locus in 3p11.1, regulating EPHA3 gene expression, is found associated with FCs of the semantic component of the language network, while a locus in 15q14, regulating THBS1 gene expression is found associated with FCs of the perceptual-motor language processing, bringing novel insights into the neurobiology of language.

Susceptibilities of CNS Cells towards Rabies Virus Infection Is Linked to Cellular Innate Immune Responses.

Rabies is caused by neurotropic rabies virus (RABV), contributing to 60,000 human deaths annually. Even though rabies leads to major public health concerns worldwide, we still do not fully understand factors determining RABV tropism and why glial cells are unable to clear RABV from the infected brain. Here, we compare susceptibilities and immune responses of CNS cell types to infection with two RABV strains, Tha and its attenuated variant Th2P-4M, mutated on phospho- (P-protein) and matrix protein (M-protein). We demonstrate that RABV replicates in human stem cell-derived neurons and astrocytes but fails to infect human iPSC-derived microglia. Additionally, we observed major differences in transcription profiles and quantification of intracellular protein levels between antiviral immune responses mediated by neurons, astrocytes (IFNB1, CCL5, CXCL10, IL1B, IL6, and LIF), and microglia (CCL5, CXCL10, ISG15, MX1, and IL6) upon Tha infection. We also show that P- and M-proteins of Tha mediate evasion of NF-κB- and JAK-STAT-controlled antiviral host responses in neuronal cell types in contrast to glial cells, potentially explaining the strong neuron-specific tropism of RABV. Further, Tha-infected astrocytes and microglia protect neurons from Tha infection via a filtrable and transferable agent. Overall, our study provides novel insights into RABV tropism, showing the interest in studying the interplay of CNS cell types during RABV infection.

Dichotomous metabolic networks govern human ILC2 proliferation and function.

Group 2 innate lymphoid cells (ILC2s) represent innate homologs of type 2 helper T cells (TH2) that participate in immune defense and tissue homeostasis through production of type 2 cytokines. While T lymphocytes metabolically adapt to microenvironmental changes, knowledge of human ILC2 metabolism is limited, and its key regulators are unknown. Here, we show that circulating 'naive' ILC2s have an unexpected metabolic profile with a higher level of oxidative phosphorylation (OXPHOS) than natural killer (NK) cells. Accordingly, ILC2s are severely reduced in individuals with mitochondrial disease (MD) and impaired OXPHOS. Metabolomic and nutrient receptor analysis revealed ILC2 uptake of amino acids to sustain OXPHOS at steady state. Following activation with interleukin-33 (IL-33), ILC2s became highly proliferative, relying on glycolysis and mammalian target of rapamycin (mTOR) to produce IL-13 while continuing to fuel OXPHOS with amino acids to maintain cellular fitness and proliferation. Our results suggest that proliferation and function are metabolically uncoupled in human ILC2s, offering new strategies to target ILC2s in disease settings.

Multitrait GWAS to connect disease variants and biological mechanisms.

Genome-wide association studies (GWASs) have uncovered a wealth of associations between common variants and human phenotypes. Here, we present an integrative analysis of GWAS summary statistics from 36 phenotypes to decipher multitrait genetic architecture and its link with biological mechanisms. Our framework incorporates multitrait association mapping along with an investigation of the breakdown of genetic associations into clusters of variants harboring similar multitrait association profiles. Focusing on two subsets of immunity and metabolism phenotypes, we then demonstrate how genetic variants within clusters can be mapped to biological pathways and disease mechanisms. Finally, for the metabolism set, we investigate the link between gene cluster assignment and the success of drug targets in randomized controlled trials.

Host genetic control of natural killer cell diversity revealed in the Collaborative Cross.

Natural killer (NK) cells are innate effectors armed with cytotoxic and cytokine-secreting capacities whose spontaneous antitumor activity is key to numerous immunotherapeutic strategies. However, current mouse models fail to mirror the extensive immune system variation that exists in the human population which may impact on NK cell-based therapies. We performed a comprehensive profiling of NK cells in the Collaborative Cross (CC), a collection of novel recombinant inbred mouse strains whose genetic diversity matches that of humans, thereby providing a unique and highly diverse small animal model for the study of immune variation. We demonstrate that NK cells from CC strains displayed a breadth of phenotypic and functional variation reminiscent of that reported for humans with regards to cell numbers, key marker expression, and functional capacities. We took advantage of the vast genetic diversity of the CC and identified nine genomic loci through quantitative trait locus mapping driving these phenotypic variations. SNP haplotype patterns and variant effect analyses identified candidate genes associated with lung NK cell numbers, frequencies of CD94+ NK cells, and expression levels of NKp46. Thus, we demonstrate that the CC represents an outstanding resource to study NK cell diversity and its regulation by host genetics.

Immune response profiling of patients with spondyloarthritis reveals signalling networks mediating TNF-blocker function in vivo.

OBJECTIVES: Antitumour necrosis factor (TNF) therapy has revolutionised treatment of several chronic inflammatory diseases, including spondyloarthritis (SpA). However, TNF inhibitors (TNFi) are not effective in all patients and the biological basis for treatment failure remains unknown. We have analysed induced immune responses to define the mechanism of action of TNF blockers in SpA and to identify immunological correlates of responsiveness to TNFi. METHODS: Immune responses to microbial and pathway-specific stimuli were analysed in peripheral blood samples from 80 patients with axial SpA before and after TNFi treatment, using highly standardised whole-blood stimulation assays. Cytokines and chemokines were measured in a Clinical Laboratory Improvement Amendments (CLIA)-certified laboratory, and gene expression was monitored using nCounter assays. RESULTS: Anti-TNF therapy induced profound changes in patients' innate immune responses. TNFi action was selective, and had only minor effects on Th1/Th17 immunity. Modular transcriptional repertoire analysis identified prostaglandin E2 synthesis and signalling, leucocyte recirculation, macrophage polarisation, dectin and interleukin (IL)-1 signalling, as well as the nuclear factor kappa B (NF-kB) transcription factor family as key pathways targeted by TNF blockers in vivo. Analysis of induced immune responses before treatment initiation revealed that expression of molecules associated with leucocyte adhesion and invasion, chemotaxis and IL-1 signalling are correlated with therapeutic responses to anti-TNF. CONCLUSIONS: We show that TNFi target multiple immune cell pathways that cooperate to resolve inflammation. We propose that immune response profiling provides new insight into the biology of TNF-blocker action in patients and can identify signalling pathways associated with therapeutic responses to biological therapies.

JASS: command line and web interface for the joint analysis of GWAS results.

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.

Exploration Deficits Under Ecological Conditions as a Marker of Apathy in Frontotemporal Dementia.

Apathy is one of the six clinical criteria for the behavioral variant of frontotemporal dementia (bvFTD), and it is almost universal in this disease. Although its consequences in everyday life are debilitating, its underlying mechanisms are poorly known, its assessment is biased by subjectivity and its care management is very limited. In this context, we have developed "ECOCAPTURE," a method aimed at providing quantifiable and objective signature(s) of apathy in order to assess it and identify its precise underlying mechanisms. ECOCAPTURE consists of the observation and recording of the patient's behavior when the participant is being submitted to a multiple-phase scenario reproducing a brief real-life situation. It is performed in a functional exploration platform transformed into a fully furnished waiting room equipped with a video and sensor-based data acquisition system. This multimodal method allowed video-based behavior analyses according to predefined behavioral categories (exploration behavior, sustained activities or inactivity) and actigraphy analyses from a 3D accelerometer. The data obtained were also correlated with behavioral/cognitive tests and scales assessing global cognitive efficiency, apathy, cognitive disinhibition, frontal syndrome, depression and anxiety. Here, bvFTD patients (n = 14) were compared to healthy participants (n = 14) during the very first minutes of the scenario, when the participants discovered the room and were encouraged to explore it. We showed that, in the context of facing a new environment, healthy participants first explored it and then engaged in sustained activities. By contrast, bvFTD patients were mostly inactive and eventually explored this new place, but in a more irregular and less efficient mode than normal subjects. This exploration deficit was correlated with apathy, disinhibition and cognitive and behavioral dysexecutive syndromes. These findings led us to discuss the presumed underlying mechanisms responsible for the exploration deficit (an inability to self-initiate actions, to integrate reward valuation and to inhibit involuntary behavior). Altogether, these results pave the way for simple and objective assessment of behavioral changes that represents a critical step for the evaluation of disease progression and efficacy of treatment in bvFTD.

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.

Adjusting for Principal Components of Molecular Phenotypes Induces Replicating False Positives.

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.

Continuation of Nesterov's Smoothing for Regression With Structured Sparsity in High-Dimensional Neuroimaging.

Predictive models can be used on high-dimensional brain images to decode cognitive states or diagnosis/prognosis of a clinical condition/evolution. Spatial regularization through structured sparsity offers new perspectives in this context and reduces the risk of overfitting the model while providing interpretable neuroimaging signatures by forcing the solution to adhere to domain-specific constraints. Total variation (TV) is a promising candidate for structured penalization: it enforces spatial smoothness of the solution while segmenting predictive regions from the background. We consider the problem of minimizing the sum of a smooth convex loss, a non-smooth convex penalty (whose proximal operator is known) and a wide range of possible complex, non-smooth convex structured penalties such as TV or overlapping group Lasso. Existing solvers are either limited in the functions they can minimize or in their practical capacity to scale to high-dimensional imaging data. Nesterov's smoothing technique can be used to minimize a large number of non-smooth convex structured penalties. However, reasonable precision requires a small smoothing parameter, which slows down the convergence speed to unacceptable levels. To benefit from the versatility of Nesterov's smoothing technique, we propose a first order continuation algorithm, CONESTA, which automatically generates a sequence of decreasing smoothing parameters. The generated sequence maintains the optimal convergence speed toward any globally desired precision. Our main contributions are: gap to probe the current distance to the global optimum in order to adapt the smoothing parameter and the To propose an expression of the duality convergence speed. This expression is applicable to many penalties and can be used with other solvers than CONESTA. We also propose an expression for the particular smoothing parameter that minimizes the number of iterations required to reach a given precision. Furthermore, we provide a convergence proof and its rate, which is an improvement over classical proximal gradient smoothing methods. We demonstrate on both simulated and high-dimensional structural neuroimaging data that CONESTA significantly outperforms many state-of-the-art solvers in regard to convergence speed and precision.

A strategy for multimodal data integration: application to biomarkers identification in spinocerebellar ataxia.

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.

Covariate selection for association screening in multiphenotype genetic studies.

Testing for associations in big data faces the problem of multiple comparisons, wherein true signals are difficult to detect on the background of all associations queried. This difficulty is particularly salient in human genetic association studies, in which phenotypic variation is often driven by numerous variants of small effect. The current strategy to improve power to identify these weak associations consists of applying standard marginal statistical approaches and increasing study sample sizes. Although successful, this approach does not leverage the environmental and genetic factors shared among the multiple phenotypes collected in contemporary cohorts. Here we developed covariates for multiphenotype studies (CMS), an approach that improves power when correlated phenotypes are measured on the same samples. Our analyses of real and simulated data provide direct evidence that correlated phenotypes can be used to achieve increases in power to levels often surpassing the power gained by a twofold increase in sample size.

Modeling a linkage between blood transcriptional expression and activity in brain regions to infer the phenotype of schizophrenia patients.

Hundreds of genetic loci participate to schizophrenia liability. It is also known that impaired cerebral connectivity is directly related to the cognitive and affective disturbances in schizophrenia. How genetic susceptibility and brain neural networks interact to specify a pathological phenotype in schizophrenia remains elusive. Imaging genetics, highlighting brain variations, has proven effective to establish links between vulnerability loci and associated clinical traits. As previous imaging genetics works in schizophrenia have essentially focused on structural DNA variants, these findings could be blurred by epigenetic mechanisms taking place during gene expression. We explored the meaningful links between genetic data from peripheral blood tissues on one hand, and regional brain reactivity to emotion task assayed by blood oxygen level-dependent functional magnetic resonance imaging on the other hand, in schizophrenia patients and matched healthy volunteers. We applied Sparse Generalized Canonical Correlation Analysis to identify joint signals between two blocks of variables: (i) the transcriptional expression of 33 candidate genes, and (ii) the blood oxygen level-dependent activity in 16 region of interest. Results suggested that peripheral transcriptional expression is related to brain imaging variations through a sequential pathway, ending with the schizophrenia phenotype. Generalization of such an approach to larger data sets should thus help in outlining the pathways involved in psychiatric illnesses such as schizophrenia. IMAGING: SEARCHING FOR LINKS TO AID DIAGNOSIS: Researchers explore links between the expression of genes associated with schizophrenia in blood cells and variations in brain activity during emotion processing. El Chérif Ibrahim and Eric Fakra at Aix-Marseille Université, France, and colleagues have developed a method to relate the expression levels of 33 schizophrenia susceptibility genes in blood cells and functional magnetic resonance imaging (fMRI) data obtained as individuals carry out a task that triggers emotional responses. Although they found no significant differences in the expression of genes between the 26 patients with schizophrenia and 26 healthy controls they examined, variations in activity in the superior temporal gyrus were strongly linked to schizophrenia-associated gene expression and presence of disease. Similar analyses of larger data sets will shed further light on the relationship between peripheral molecular changes and disease-related behaviors and ultimately, aid the diagnosis of neuropsychiatric disease.

Adaptive human immunity drives remyelination in a mouse model of demyelination.

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.

Transcriptome Analysis of Peripheral Blood in Chronic Inflammatory Demyelinating Polyradiculoneuropathy Patients Identifies TNFR1 and TLR Pathways in the IVIg Response.

We have studied the response to intravenous immunoglobulins (IVIg) by a transcriptomic approach in 11 chronic inflammatory demyelinating polyradiculoneuropathy (CIDP) patients (CIDP duration = 6 [0.83-6.5] years). RNA was extracted from cells in whole blood collected before and 3 weeks after IVIg treatment, and hybridized on Illumina chips. After RNA quality controls, gene expression was analyzed using statistical tests fitted for microarrays (R software, limma package), and a pathway analysis was performed using DAVID software. We identified 52 genes with expression that varied significantly after IVIg (fold change [FC] > 1.2, P < 0.001, false discovery rate [FDR] <0.05). Among these 52 genes, 7 were related to immunity, 3 were related to the tumor necrosis factor (TNF)-α receptor 1 (TNFR1) pathway (inhibitor of caspase-activated DNase (ICAD): FC = 1.8, P = 1.7E-7, FDR = 0.004; p21 protein-activated kinase 2 [PAK2]: FC = 1.66, P = 2.6E-5, FDR = 0.03; TNF-α-induced protein 8-like protein 1 [TNFAIP8L1]: P = 1.00E-05, FDR = 0.026), and 2 were related to Toll-like receptors (TLRs), especially TLRs 7 and 9, and were implicated in autoimmunity. These genes were UNC93B1 (FC = 1.6, P = 2E-5, FDR = 0.03), which transports TLRs 7 and 9 to the endolysosomes, and RNF216 (FC = 1.5, P = 1E-05, FDR = 0.03), which promotes TLR 9 degradation. Pathway analysis showed that the TNFR1 pathway was significantly lessened by IVIg (enrichment score = 24, Fischer exact test = 0.003). TNF-α gene expression was higher in responder patients than in nonresponders; however, it decreased after IVIg in responders (P = 0.04), but remained stable in nonresponders. Our data suggest the actions of IVIg on the TNFR1 pathway and an original mechanism involving innate immunity through TLRs in CIDP pathophysiology and the response to IVIg. We conclude that responder patients have stronger inflammatory activity that is lessened by IVIg.

Clinical-genetic model predicts incident impulse control disorders in Parkinson's disease.

OBJECTIVES: Impulse control disorders (ICD) are commonly associated with dopamine replacement therapy (DRT) in patients with Parkinson's disease (PD). Our aims were to estimate ICD heritability and to predict ICD by a candidate genetic multivariable panel in patients with PD. METHODS: Data from de novo patients with PD, drug-naïve and free of ICD behaviour at baseline, were obtained from the Parkinson's Progression Markers Initiative cohort. Incident ICD behaviour was defined as positive score on the Questionnaire for Impulsive-Compulsive Disorders in PD. ICD heritability was estimated by restricted maximum likelihood analysis on whole exome sequencing data. 13 candidate variants were selected from the DRD2, DRD3, DAT1, COMT, DDC, GRIN2B, ADRA2C, SERT, TPH2, HTR2A, OPRK1 and OPRM1 genes. ICD prediction was evaluated by the area under the curve (AUC) of receiver operating characteristic (ROC) curves. RESULTS: Among 276 patients with PD included in the analysis, 86% started DRT, 40% were on dopamine agonists (DA), 19% reported incident ICD behaviour during follow-up. We found heritability of this symptom to be 57%. Adding genotypes from the 13 candidate variants significantly increased ICD predictability (AUC=76%, 95% CI (70% to 83%)) compared to prediction based on clinical variables only (AUC=65%, 95% CI (58% to 73%), p=0.002). The clinical-genetic prediction model reached highest accuracy in patients initiating DA therapy (AUC=87%, 95% CI (80% to 93%)). OPRK1, HTR2A and DDC genotypes were the strongest genetic predictive factors. CONCLUSIONS: Our results show that adding a candidate genetic panel increases ICD predictability, suggesting potential for developing clinical-genetic models to identify patients with PD at increased risk of ICD development and guide DRT management.