We also assessed agreement between our study and a recent study by Lisboa and colleagues37, which evaluated DEG in three striatal regions of six OCD and eight control subjects, by comparing DEG across the studies. After identifying a set of genes showing agreement between the two studies using the function ‘fdrThreshold’ (BioNet) to produce a p-value cutoff that corresponds to an FDR of 0.20 for the two striatal regions in our dataset, we determined gene set enrichment for one brain region, the caudate. We then analyzed enriched gene ontology (GO) terms using REVIGO42, which uses semantic similarity43 to cluster similar terms and prioritize terms for semantic interpretation by degree of enrichment. We used a similarity setting of 0.5, favoring a shorter and semantically diverse list of functions, and two default settings (semantic similarity measure SimRel and the database for GO term sizes “whole UniProt”). We used the online version of REVIGO (GO release “go_monthly-termdb.obo-xml.gz” [Jan 2017] and UniProt-to-GO mapping file “goa_uniprot_gcrp.gaf.gz” [15 Mar 2017]).
