The identification of predictors of treatment response is crucial for improving treatment outcome for children with anxiety disorders. Machine learning methods provide opportunities to identify combinations of factors that contribute to risk prediction models.

A machine learning approach was applied to predict anxiety disorder remission in a large sample of 2114 anxious youth (5–18 years). Potential predictors included demographic, clinical, parental, and treatment variables with data obtained pre-treatment, post-treatment, and at least one follow-up.

All machine learning models performed similarly for remission outcomes, with AUC between 0.67 and 0.69. There was significant alignment between the factors that contributed to the models predicting two target outcomes: remission of all anxiety disorders and the primary anxiety disorder. Children who were older, had multiple anxiety disorders, comorbid depression, comorbid externalising disorders, received group treatment and therapy delivered by a more experienced therapist, and who had a parent with higher anxiety and depression symptoms, were more likely than other children to still meet criteria for anxiety disorders at the completion of therapy. In both models, the absence of a social anxiety disorder and being treated by a therapist with less experience contributed to the model predicting a higher likelihood of remission.

These findings underscore the utility of prediction models that may indicate which children are more likely to remit or are more at risk of non-remission following CBT for childhood anxiety.

Anxiety disorders are common, and cognitive–behavioural therapy (CBT) is a first-line treatment. Candidate gene studies have suggested a genetic basis to treatment response, but findings have been inconsistent.

To perform the first genome-wide association study (GWAS) of psychological treatment response in children with anxiety disorders (n = 980).

Presence and severity of anxiety was assessed using semi-structured interview at baseline, on completion of treatment (post-treatment), and 3 to 12 months after treatment completion (follow-up). DNA was genotyped using the Illumina Human Core Exome-12v1.0 array. Linear mixed models were used to test associations between genetic variants and response (change in symptom severity) immediately post-treatment and at 6-month follow-up.

No variants passed a genome-wide significance threshold (P=5×10–8) in either analysis. Four variants met criteria for suggestive significance (P<5×10–6) in association with response post-treatment, and three variants in the 6-month follow-up analysis.

This is the first genome-wide therapygenetic study. It suggests no common variants of very high effect underlie response to CBT. Future investigations should maximise power to detect single-variant and polygenic effects by using larger, more homogeneous cohorts.