A, Acquisition of anatomical image. B, Surface reconstruction was carried out using a software program (FreeSurfer, version 5.3.0; http://surfer.nmr.mgh.harvard.edu/). Local gyrification index (LGI) estimation was performed using the method by Schaer et al.44 C, The atlas by Desikan et al46 was used for parcellating the cortical surface to 68 regions (34 on each hemisphere). D, Association matrices were obtained by calculating the correlations between regional gyrification across individuals within each group separately. E, Binary adjacency matrices were derived from bias estimates using the jackknife procedure for each individual within the 4 study groups. F, Binarization used the selected range of cost densities whereby the resulting graphs were always fully connected and had small-world properties.
The mean (SD) scores of each graph variable were estimated from functional data analysis with the selected density range of 0.05 to 0.25. ARMS-NT indicates at-risk mental state nontransition; ARMS-T, at-risk mental state transition; FEP, first-episode psychosis; and HC, healthy control.
aP < .001 compared with HC.
bP < .001 compared with ARMS-NT.
eMethods. Supplemental Methods
eFigure 1. Parcellated Brain Regions Based on Desikan’s Atlas
eFigure 2. Receiver Operating Characteristic Curve for the Graph-Based Classification Performance
eFigure 3. Permutation Results (Balanced Accuracy and Area Under the Curve)
eTable 1. Pearson’s Correlations Between Graph Variables and IQ (Premorbid) Within Each Group
eTable 2. Raw Local Gyrification Indices (LGI) Across Groups
eTable 3. Classification Performance Measures With Raw LGIs
eTable 4.P Values of Permutation Tests
eResults. Supplemental Results
Customize your JAMA Network experience by selecting one or more topics from the list below.
Das T, Borgwardt S, Hauke DJ, et al. Disorganized Gyrification Network Properties During the Transition to Psychosis. JAMA Psychiatry. Published online April 25, 2018. doi:10.1001/jamapsychiatry.2018.0391
Do patients at high risk of psychosis exhibit disorganized gyrification network properties, and is it possible to predict transition to a first episode of psychosis on the basis of gyrification?
In this cross-sectional magnetic resonance imaging study among 161 individuals in 4 study groups, patients who later develop psychosis exhibit disorganized gyrification network properties even at a clinically high-risk state. Gyrification network measures predict the future outcome of transition with more than 80% accuracy.
Constructing gyrification-based networks from structural magnetic resonance imaging may facilitate individual prediction of future psychosis in those at clinical high risk for psychosis.
There is urgent need to improve the limited prognostic accuracy of clinical instruments to predict psychosis onset in individuals at clinical high risk (CHR) for psychosis. As yet, no reliable biological marker has been established to delineate CHR individuals who will develop psychosis from those who will not.
To investigate abnormalities in a graph-based gyrification connectome in the early stages of psychosis and to test the accuracy of this systems-based approach to predict a transition to psychosis among CHR individuals.
Design, Setting, and Participants
This investigation was a cross-sectional magnetic resonance imaging (MRI) study with follow-up assessment to determine the transition status of CHR individuals. Participants were recruited from a specialized clinic for the early detection of psychosis at the Department of Psychiatry (Universitäre Psychiatrische Kliniken [UPK]), University of Basel, Basel, Switzerland. Participants included individuals in the following 4 study groups: 44 healthy controls (HC group), 63 at-risk mental state (ARMS) individuals without later transition to psychosis (ARMS-NT group), 16 ARMS individuals with later transition to psychosis (ARMS-T group), and 38 antipsychotic-free patients with first-episode psychosis (FEP group). The study dates were November 2008 to November 2014. The dates of analysis were March to November 2017.
Main Outcomes and Measures
Gyrification-based structural covariance networks (connectomes) were constructed to quantify global integration, segregation, and small-worldness. Group differences in network measures were assessed using functional data analysis across a range of network densities. The extremely randomized trees algorithm with repeated 5-fold cross-validation was used to delineate ARMS-T individuals from ARMS-NT individuals. Permutation tests were conducted to assess the significance of classification performance measures.
The 4 study groups comprised 161 participants with mean (SD) ages ranging from 24.0 (4.7) to 25.9 (5.7) years. Small-worldness was reduced in the ARMS-T and FEP groups and was associated with decreased integration and increased segregation in both groups (Hedges g range, 0.666-1.050). Using the connectome properties as features, a good classification performance was obtained (accuracy, 90.49%; balanced accuracy, 81.34%; positive predictive value, 84.47%; negative predictive value, 92.18%; sensitivity, 66.11%; specificity, 96.58%; and area under the curve, 88.30%).
Conclusions and Relevance
These findings suggest that there is poor integration in the coordinated development of cortical folding in patients who develop psychosis. These results further suggest that gyrification-based connectomes might be a promising means to generate systems-based measures from anatomical data to improve individual prediction of a transition to psychosis in CHR individuals.
Predicting psychosis onset in individuals at clinical high risk (CHR) for psychosis1 is essential to administer preventive interventions. However, it is not yet possible to make any personalized prediction of psychosis onset relying only on the initial clinical assessment.2 Therefore, research is striving for reliable brain markers to improve the prediction of psychosis onset in these individuals.3,4
To date, most studies have used structural magnetic resonance imaging (MRI) to investigate the alterations in regional gray matter volume (GMV) in CHR individuals.5 However, most of the available evidence shows that GMV reductions seen in patients occur before the transition to psychosis, during the transition, or in the immediate postonset phase, which may indicate that morphometric methods, such as the measurement of gyrification, are perhaps more sensitive to detect the pathophysiology in the prodromal phase.6 Indeed, altered local gyrification indexes (LGIs) have been reported not only in patients with schizophrenia7 but also in patients with first-episode psychosis (FEP),8 those at genetic risk for schizophrenia,9,10 and CHR individuals,11 as well as in those who later developed an overt psychotic disorder.12-14 However, these localized regional changes fail to quantify the association between concomitant changes in different brain areas.15 The gross morphology of the developing brain undergoes several well-coordinated maturational events to establish neural networks,16 and developmental disturbances in the topological organization of these networks can result in various psychiatric disorders, such as schizophrenia.17,18 Maldevelopment of the structural connectome can be inferred by studying the covariance of morphology using graph theory.15 Such graph-based network studies capture an important aspect of developmental maturation that is crucial for understanding the pathophysiology of psychotic disorders.19,20 Previous studies reported reduced small-worldness of structural brain networks in patients with schizophrenia,21-23 CHR individuals,24 people at increased familial risk for schizophrenia,25-27 and those with subclinical psychotic experiences,28 characterized by increased segregation and decreased interaction of anatomical covariance (see the published reviews18,19,29-31 of network analyses in schizophrenia). Of various regional morphometric properties that can be assessed for structural covariance, gyrification is a compelling marker of early neurodevelopment. The overall pattern of cortical gyrification is established by year 2 of human life and remains consistent for most of the adult life.32 Perinatal complications that affect brain development result in aberrant gyrification.33,34 Notably, experimental introduction of white matter lesions in early life produces localized changes in cortical folding in distant regions that are axonally connected.35 Therefore, aberrations in structural covariance of gyrification can provide an index of the integrity of cortical connectivity in early life. It has previously been shown that the gyrification-based structural connectome is indeed not normal in patients with schizophrenia, with an abnormally segregated pattern of cortical folding in distributed brain regions, especially in those with more severe illness36 and in nonresponders to antipsychotics.37
This study explored the topological organization of gyrification networks in the following 4 study groups: 44 healthy controls (HC group), 63 at-risk mental state (ARMS) individuals without later transition to psychosis (ARMS-NT group), 16 ARMS individuals with later transition to psychosis (ARMS-T group), and 38 antipsychotic-free patients with FEP (FEP group). The study dates were November 2008 to November 2014. The dates of analysis were March to November 2017. Our aims were to investigate whether the network properties of the gyrification-based connectome differed at baseline among the 4 study groups and to test if this information was sufficiently discriminatory to delineate ARMS individuals transitioning vs nontransitioning to psychosis.
We recruited 44 HCs, 79 ARMS individuals, and 38 antipsychotic-free patients with FEP in our specialized clinic for the early detection of psychosis at the Department of Psychiatry (Universitäre Psychiatrische Kliniken [UPK]), University of Basel, Basel, Switzerland. All participants provided written informed consent, and the study was approved by the local ethics committee (Ethikkommission Nordwest-und Zentralschweiz).
The ARMS and FEP statuses were assessed using the Basel Screening Instrument for Psychosis (BSIP),38 the Brief Psychiatric Rating Scale (BPRS),39 the Scale for the Assessment of Negative Symptoms (SANS),40 and the Global Assessment of Functioning (GAF).41 In addition, we recorded current and previous psychotropic medication, as well as nicotine and illegal drug consumption, by using a semistructured interview. The following exclusion criteria were applied to all study groups: history of a psychotic disorder; psychotic symptoms secondary to an organic disorder; substance abuse according to International Statistical Classification of Diseases and Related Health Problems, 10th Revision (ICD-10) research criteria; psychotic symptoms associated with a bipolar disorder, major depression, or a borderline personality disorder; age younger than 18 years; inadequate knowledge of the German language; and IQ less than 70, measured with the Mehrfachwahl-Wortschatz-Test Form B.42
The ARMS status required 1 or more of the following: (1) attenuated psychoticlike symptoms, (2) brief limited intermittent psychotic symptoms (BLIPS), or (3) a first-degree or second-degree relative with a psychotic disorder plus at least 2 additional risk factors for or indicators of beginning psychosis according to the BSIP screening instrument. Inclusion because of attenuated psychotic symptoms required that change in mental state had to be present at least several times a week and for more than 1 week (a score of 2 or 3 on the BPRS hallucination item or a score of 3 or 4 on BPRS items for unusual thought content or suspiciousness). Inclusion because of BLIPS required a score of 4 or higher on the BPRS hallucination item or a score of 5 or higher on BPRS items for unusual thought content, suspiciousness, or conceptual disorganization, with each symptom lasting less than 1 week before resolving spontaneously. After the baseline assessment, ARMS individuals were followed up clinically and received standard psychiatric case management (mean [SD] follow-up, 3.8 [3.2] years). All ARMS individuals were antipsychotic naive, while 29 were taking low-dose antidepressants at the time of imaging (see the eMethods in the Supplement). Sixteen ARMS individuals have transitioned to psychosis (transition rate, 20%; mean [SD] time of transition after MRI, 19.85 [19.60] months).
Patients with FEP were those who met the criteria for a transition to psychosis according to the classification by Yung et al.43 These patients already fulfilled criteria for acute psychotic disorder according to the ICD-10 or the Diagnostic and Statistical Manual of Mental Disorders (Fourth Edition) (DSM-IV) but not yet for schizophrenia. Inclusion required a score of 4 or higher on the BPRS hallucination item or a score of 5 or higher on BPRS items for unusual thought content, suspiciousness, or conceptual disorganization. The symptoms must have occurred at least several times a week and persisted for more than 1 week. All patients with FEP were antipsychotic free at the time of imaging. Six patients with FEP were taking antidepressants (see the eMethods in the Supplement).
Healthy controls were recruited from the same geographic area as the other groups. They had no current psychiatric disorder; no history of a psychiatric illness, head trauma, neurologic illness, serious medical or surgical illness, or substance abuse; and no family history of any psychiatric disorder as assessed by an experienced psychiatrist in a detailed clinical interview.
Magnetic resonance imaging data were obtained in all participants. Details of the data acquisition protocol can be found in the eMethods in the Supplement.
Local gyrification indexes were obtained using the method by Schaer et al,44 with images reconstructed via a software program (FreeSurfer, version 5.3.0; http://surfer.nmr.mgh.harvard.edu/).45 The LGIs were computed for 68 parcellated brain regions according to the atlas by Desikan et al46 (see eFigure 1 in the Supplement). More details can be found in the eMethods in the Supplement.
Multiple analysis of covariance (MANCOVA) was performed. The comparison used LGIs as dependent variables, group as a fixed factor, and age, sex, and intracranial volume as covariates.
We first generated a 68 × 68 correlation matrix based on LGIs for each of the 4 study groups, adjusted for the effect of age, sex, and intracranial volume. We then used a jackknife bias estimation procedure to determine each individual’s contribution to the overall group-level covariance structure, providing an individual-wise 68 × 68 distance matrix for the 4 study groups (see the eMethods in the Supplement for more information).47,48 A graph analysis toolbox49 was used for studying various topological properties. For the weighted connectivity matrices obtained from each individual, a range of network thresholds based on connection density (ie, 0.05-0.25, with interval steps of 0.01) was applied to generate binary undirected adjacency matrices. This choice of range enabled between-group comparisons in topological measures across graphs with comparable number of edges but without inducing disconnection or losing small-worldness. Topological measures were normalized to equivalent values derived from 20 random (“null”) networks with the same degree of distribution. Figure 1 shows a flow diagram of the graph analysis.
In a gyrification-based connectome, the simplest measure of the overall network strength is provided by the average degree of individual nodes. Modular development and regionally selective functional dependency within clusters of high covariance result in a high degree of segregation. In contrast, high integration indicates a coordinated maturational process affecting the entire brain; this effect may result from the presence of certain “central” hub regions whose structure covaries with a large number of other brain regions, leading to widely distributed structural coupling. Such hub regions often show a preferential covariance among each other, making the connectome resilient to pathological processes affecting a single hub region. Assortativity is an index of such resilience. The presence of high segregation in the context of optimum integration gives rise to small-worldness, measured using the small-world index (σ). The σ is the ratio between the normalized clustering coefficient (ϒ) and the normalized characteristic path length (λ). The ϒ is normalized for each node and averaged across the network; hence, it is likely to be influenced disproportionately by low-degree (or weakly connected) nodes. Transitivity, another topological measure, is free of this bias because it is normalized collectively, providing more intuitive information regarding segregated structural covariance. See the eMethods in the Supplement for formal descriptions of these properties.
To perform statistical group comparisons across the range of densities between 0.05 and 0.25, we first constructed curves showing the change in the measures of interest as a function of the density. Functional data analysis was performed with topological measures treated as a function of y = f(x), allowing summation across densities and obviating the need for multiple testing. A 1-way analysis of variance was used, followed by post hoc t tests to compare the density function as obtained using functional data analysis among the 4 study groups. Given that we examined 6 topological measures, Bonferroni-corrected 2-tailed P = .05 divided by 6 (ie, .008) was chosen as the threshold of statistical significance. The unbiased effect size between-group comparisons were estimated using Hedges g.50
Spearman rank correlations were used to test the association between topological measures and clinical variables (BPRS, SANS, and GAF total scores). These correlations were tested for each patient group separately.
The extremely randomized trees algorithm51 (see the eMethods in the Supplement) with all 6 network measures as features was used to delineate ARMS-NT individuals from AMRS-T individuals. We were specifically interested in whether we could improve the positive predictive value (PPV) because clinical instruments have a limited ability to rule in heightened risk of subsequent psychosis in CHR individuals.2 While increasing PPV is important to make clinically viable predictions, reducing the number of false-positives to a minimum is essential, especially because treatments to prevent transition are of uncertain efficacy at present. In line with prior work,52 we thus calculated 2 diagnostic indexes in addition to pretest and posttest probabilities and the positive likelihood ratio. These factors included the predictive summary index (PSI) to quantify the total amount of uncertainty reduced by the gyrification test (PSI = [PPV + NPV] − 1) (where NPV indicates the negative predictive value) and the number needed to predict (NNP) (NNP = 1 / PSI),53 an estimate of the number of patients that need to be examined to correctly predict diagnosis in one person.
The 4 study groups comprised 161 participants with mean (SD) ages ranging from 24.0 (4.7) to 25.9 (5.7) years. The 4 study groups were well matched on age, handedness, cannabis consumption, and antidepressant treatment, but they differed on sex, years of education and premorbid IQ, intracranial volume, cigarettes per day and alcohol consumption, and global functioning. The 3 clinical groups differed on the BPRS and GAF total scores but not on the SANS total score (Table 1). Because premorbid IQ differed between study groups as a feature of the illness, it was not included as a covariate in all group comparisons (see eTable 1 in the Supplement for correlations between IQ and graph variables).
MANCOVA revealed no group effect on regional LGIs (F = 0.9, P = .666). eTable 2 in the Supplement lists the raw LGI values.
We noted statistically significant group differences in the topological properties of λ, σ, and transitivity but not in the overall network strength, assortativity, or ϒ (Table 2 and Figure 2). Post hoc comparisons revealed that small-worldness σ was reduced in the ARMS-T and FEP groups compared with the HC and ARMS-NT groups, with differences between the ARMS-NT and ARMS-T groups showing large effect sizes (Hedges g, >1). Characteristic path length λ was significantly higher in the ARMS-T and FEP groups, indicating reduced integration compared with the HC and ARMS-NT groups. In both cases (for small-worldness σ and characteristic path length λ), there was no difference between the ARMS-T and FEP groups, indicating that these parameters are observable before the transition to psychosis and remain unaltered in the presence of FEP. In addition, we found that transitivity (segregation) was significantly higher in the ARMS-T and FEP groups compared with the HC and ARMS-NT groups. Transitivity was higher in the FEP group compared with the ARMS-T group.
Associations between BPRS, SANS, and GAF total scores and the 6 network measures were tested for each patient group (ARMS-NT, ARMS-T, and FEP) separately. There were no significant correlations between clinical variables and network metrics.
Using network measures as features (see the eResults and eTable 3 in the Supplement for classification performance on the raw LGIs), we obtained a good classification performance (Table 3 and eFigure 2 in the Supplement). With 5000 permutations, both balanced accuracy (classification value, 0.813; permutation mean, 0.499; permutation SD, 0.029) and area under the curve (classification value, 0.883; permutation mean, 0.499; permutation SD, 0.094) reached statistical significance (P < .001) (see eFigure 3 in the Supplement for supporting data on permutation tests). eTable 4 in the Supplement lists details for all permutation tests. Based on this test performance, we estimated a range of diagnostic indexes (Table 3). Notably, the posttest probability of a transition to psychosis increased by more than 60%, up to 83% and 87% compared with the pretest probability, based on standard clinical assessment both in our sample and compared with meta-analytic validity of psychometric CHR interviews.2 We estimated the PSI to be 0.76 and the NNP to be 1.30. These estimates indicate that the use of gyrification connectomes can reduce the uncertainty in predicting a transition to psychosis by 76%, a large gain that is likely to be clinically significant. The NNP of 1.30 also indicates that the overall testing burden is likely to be tolerable to gain the improved precision as observed in this study.
To our knowledge, this is the first gyrification-based connectomic study to predict the transition to psychosis from a CHR state. We had 3 major findings. First, we found that ARMS individuals who later transition to FEP already show abnormalities in the gyrification connectome, indicating subtle neurodevelopmental aberrations in this group. In particular, those who later transition to psychosis demonstrate a highly segregated but poorly integrated pattern in structural covariance of cortical folding compared with nonconverters and healthy individuals. Second, the alterations in the gyrification connectome seen in the ARMS-T group were also present in antipsychotic-free patients with FEP. Third, by using topological measures of the gyrification connectome, we were able to predict the future outcome of a transition to psychosis with 81.34% balanced accuracy.
In search of the neural mechanisms that contribute to emerging psychosis, developmental processes are increasingly being recognized to be crucial for a transition to psychosis.54 Structural covariance based on morphometric measures, especially gyrification, represents synchronized developmental changes.55-57 When developmental lesions are experimentally induced in the white matter, covarying proximal and distal changes occur in gyrification patterns,35 indicating the influence of axonal integrity during early life. Because most cortical folding is complete in fetal life,32 it is likely that the topological abnormalities associated with a transition to psychosis that we report herein had been already present in early childhood. Various perinatal insults and defenses against these insults could result in a pattern of segregated cortical development.33,34 Nevertheless, the lack of regional changes in the degree of cortical folding among the 4 study groups indicates that either these insults were widespread (not localized) but brief or were partially compensated by other factors controlling the morphogenesis.58
As argued by Van Os and Delespaul,59 the number of CHR individuals who need to be treated to prevent one case of full-blown psychotic disorder relates directly to both the PPV of the prognostic test used and the success rate of the preventive treatment. Assuming the latter to be 50%,59 we would need a PPV greater than 50% to obtain a number needed to treat of 4. Given that only 4% to 5% of patients who eventually develop FEP seek help for CHR features even in specialized settings,60,61 it is important that false-negatives are kept to a minimum. To date, approaches like the North American Prodrome Longitudinal Study risk calculator62 and the polyenviromic risk score risk prediction tool63 have not provided a high PPV alongside a high sensitivity. The gyrification connectome offers predictors that appear to achieve both high PPV and moderate sensitivity.
Our results are comparable with the findings of GMV-based transition prediction biomarkers.64-66 While we achieved comparable balanced accuracy and higher PPV, NPV, and specificity, our sensitivity was lower. Low sensitivity may relate to sampling imbalance (63 ARMS-NT vs 16 ARMS-T), although it reflects the true base rates in this population. Furthermore, pathophysiological (ie, GMV vs gyrification67) and methodological differences (support vector machine vs extremely randomized trees algorithm, as well as different cross-validation strategies) may also explain the divergent results. Gyrification connectomes may be well suited as part of sequential testing approaches to improve psychosis prediction.3 Our results need to be validated in larger, enriched samples to ensure their real-world utility.
Several limitations of our study merit comment. It is not yet clear what is the best way of defining individual nodes when constructing graph networks.68 We used neuroanatomically defined boundaries of cortical folding to generate intuitively meaningful values of gyrification that can be applied invariably across different groups being compared and be reliably replicated in future studies. This approach is also likely to provide greater convergence with developmental changes.57 We did not assess perinatal complications and thus could not explore whether developmental problems have driven the topological abnormalities seen in those transitioning to psychosis. Also, we did not consider other developmentally influenced brain measures, such as surface area, that might facilitate improve psychosis prediction.69 Furthermore, we had approximately a 20% transition rate in our sample during almost 4 years of follow-up. Although the transition to psychosis generally plateaus the first 2 years after the initial clinical assessment and the transition rate after this period is likely to be small,70 the duration of the follow-up period should be considered when testing the predictive accuracy of gyrification networks. The lack of a significant association between topological alterations and symptom expression should be considered with caution; disorganized gyrification could also be a marker of poor cognitive or functional outcomes not assessed in this study.
In summary, we provide the first report to date that a transition to psychosis is associated with developmental disruptions in the morphogenesis of cortical folding. This observation makes perturbed neurodevelopment directly relevant to the neurobiology of psychosis onset in a sample with clinically defined ARMS. Gyrification connectomes can potentially be used to provide sample enrichment among CHR individuals to promote prevention of psychosis.
Accepted for Publication: February 4, 2018.
Corresponding Author: André Schmidt, PhD, Department of Psychiatry (Universitäre Psychiatrische Kliniken [UPK]), University of Basel, Wilhelm Klein Strasse 27, 4012 Basel, Switzerland (email@example.com).
Published Online: April 25, 2018. doi:10.1001/jamapsychiatry.2018.0391
Open Access: This article is published under the JN-OA license and is free to read on the day of publication.
Author Contributions: Drs Palaniyappan and Schmidt are joint last authors. Drs Borgwardt and Schmidt had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
Study concept and design: Borgwardt, Lang, Palaniyappan, Schmidt.
Acquisition, analysis, or interpretation of data: Das, Borgwardt, Hauke, Harrisberger, Riecher-Rössler, Palaniyappan, Schmidt.
Drafting of the manuscript: Das, Hauke, Palaniyappan, Schmidt.
Critical revision of the manuscript for important intellectual content: All authors.
Statistical analysis: Das, Hauke, Harrisberger, Palaniyappan, Schmidt.
Obtained funding: Borgwardt, Riecher-Rössler, Schmidt.
Administrative, technical, or material support: Borgwardt, Harrisberger, Palaniyappan, Schmidt.
Study supervision: Borgwardt, Lang, Schmidt.
Conflict of Interest Disclosures: Dr Palaniyappan reported receiving speaker fees from Otsuka Canada in 2017 and reported that the Prevention and Early Intervention Program for Psychoses (PEPP) in London, Ontario, Canada, received investigator-led educational grants from Janssen Canada and Otsuka Canada in 2017, initiated by Dr Palaniyappan. In the last 2 years, Dr Palaniyappan reported holding shares of Shire Inc and GlaxoSmithKline UK in his or his spousal pension funds. No other disclosures were reported.
Funding/Support: Dr Borgwardt was supported by a research grant from the Universitäre Psychiatrische Kliniken [UPK] Basel. Dr Palaniyappan is supported by the Bucke Family Fund, Chrysalis Foundation Fund, Opportunities Fund from the Academic Medical Organization of South Western Ontario, and Canadian Institutes of Health Research Foundation Grant (application 375104).
Role of the Funder/Sponsor: The funding sources had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.