Across scientific disciplines, there is a rapidly growing recognition of the need for more statistically robust, transparent approaches to data visualization. Complementary to this, many scientists have called for plotting tools that accurately and transparently convey key aspects of statistical effects and raw data with minimal distortion. Previously common approaches, such as plotting conditional mean or median barplots together with error-bars have been criticized for distorting effect size, hiding underlying patterns in the raw data, and obscuring the assumptions upon which the most commonly used statistical tests are based. Here we describe a data visualization approach which overcomes these issues, providing maximal statistical information while preserving the desired ‘inference at a glance’ nature of barplots and other similar visualization devices. These “raincloud plots” can visualize raw data, probability density, and key summary statistics such as median, mean, and relevant confidence intervals in an appealing and flexible format with minimal redundancy. In this tutorial paper, we provide basic demonstrations of the strength of raincloud plots and similar approaches, outline potential modifications for their optimal use, and provide open-source code for their streamlined implementation in R, Python and Matlab (https://github.com/RainCloudPlots/RainCloudPlots). Readers can investigate the R and Python tutorials interactively in the browser using Binder by Project Jupyter.

Significance: Adolescence is a period of human brain growth and high incidence of mental health disorders. Here, we show consistently in two MRI cohorts that human brain changes in adolescence were concentrated on the more densely connected hubs of the connectome (i.e., association cortical regions that mediated efficient connectivity throughout the human brain structural network). Hubs were less myelinated at 14 y but had faster rates of myelination and cortical shrinkage in the 14- to 24-y period. This topologically focused process of cortical consolidation was associated with expression of genes enriched for normal synaptic and myelin-related processes and risk of schizophrenia. Consolidation of anatomical network hubs could be important for normal and clinically disordered adolescent brain development.

Across scientific disciplines, there is a rapidly growing recognition of the need for more statistically robust, transparent approaches to data visualization. Complementary to this, many scientists have called for plotting tools that accurately and transparently convey key aspects of statistical effects and raw data with minimal distortion. Previously common approaches, such as plotting conditional mean or median barplots together with error-bars have been criticized for distorting effect size, hiding underlying patterns in the raw data, and obscuring the assumptions upon which the most commonly used statistical tests are based. Here we describe a data visualization approach which overcomes these issues, providing maximal statistical information while preserving the desired ‘inference at a glance’ nature of barplots and other similar visualization devices. These “raincloud plots” can visualize raw data, probability density, and key summary statistics such as median, mean, and relevant confidence intervals in an appealing and flexible format with minimal redundancy. In this tutorial paper, we provide basic demonstrations of the strength of raincloud plots and similar approaches, outline potential modifications for their optimal use, and provide open-source code for their streamlined implementation in R, Python and Matlab (https://github.com/RainCloudPlots/RainCloudPlots). Readers can investigate the R and Python tutorials interactively in the browser using Binder by Project Jupyter.

Schizophrenia has been conceived as a disorder of brain connectivity, but it is unclear how this network phenotype is related to the underlying genetics. We used morphometric similarity analysis of MRI data as a marker of interareal cortical connectivity in three prior case–control studies of psychosis: in total, n = 185 cases and n = 227 controls. Psychosis was associated with globally reduced morphometric similarity in all three studies. There was also a replicable pattern of case–control differences in regional morphometric similarity, which was significantly reduced in patients in frontal and temporal cortical areas but increased in parietal cortex. Using prior brain-wide gene expression data, we found that the cortical map of case–control differences in morphometric similarity was spatially correlated with cortical expression of a weighted combination of genes enriched for neurobiologically relevant ontology terms and pathways. In addition, genes that were normally overexpressed in cortical areas with reduced morphometric similarity were significantly up-regulated in three prior post mortem studies of schizophrenia. We propose that this combined analysis of neuroimaging and transcriptional data provides insight into how previously implicated genes and proteins as well as a number of unreported genes in their topological vicinity on the protein interaction network may drive structural brain network changes mediating the genetic risk of schizophrenia.

Motivated by prior data on local cortical shrinkage and intracortical myelination, we predicted age-related changes in topological organization of cortical structural networks during adolescence. We estimated structural correlation from magnetic resonance imaging measures of cortical thickness at 308 regions in a sample of N = 297 healthy participants, aged 14–24 years. We used a novel sliding-window analysis to measure age-related changes in network attributes globally, locally and in the context of several community partitions of the network. We found that the strength of structural correlation generally decreased as a function of age. Association cortical regions demonstrated a sharp decrease in nodal degree (hubness) from 14 years, reaching a minimum at approximately 19 years, and then levelling off or even slightly increasing until 24 years. Greater and more prolonged age-related changes in degree of cortical regions within the brain network were associated with faster rates of adolescent cortical myelination and shrinkage. The brain regions that demonstrated the greatest age-related changes were concentrated within prefrontal modules. We conclude that human adolescence is associated with biologically plausible changes in structural imaging markers of brain network organization, consistent with the concept of tuning or consolidating anatomical connectivity between frontal cortex and the rest of the connectome.