Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2021 May;17(5):e10267.
doi: 10.15252/msb.202110267.

Behavioral fingerprints predict insecticide and anthelmintic mode of action

Adam McDermott-Rouse  1   2 , Eleni Minga  1   2 , Ida Barlow  1   2 , Luigi Feriani  1   2 , Philippa H Harlow  3 , Anthony J Flemming  3 , André E X Brown  1   2
Affiliations

Behavioral fingerprints predict insecticide and anthelmintic mode of action

Adam McDermott-Rouse et al. Mol Syst Biol. 2021 May.

Abstract

Novel invertebrate-killing compounds are required in agriculture and medicine to overcome resistance to existing treatments. Because insecticides and anthelmintics are discovered in phenotypic screens, a crucial step in the discovery process is determining the mode of action of hits. Visible whole-organism symptoms are combined with molecular and physiological data to determine mode of action. However, manual symptomology is laborious and requires symptoms that are strong enough to see by eye. Here, we use high-throughput imaging and quantitative phenotyping to measure Caenorhabditis elegans behavioral responses to compounds and train a classifier that predicts mode of action with an accuracy of 88% for a set of ten common modes of action. We also classify compounds within each mode of action to discover substructure that is not captured in broad mode-of-action labels. High-throughput imaging and automated phenotyping could therefore accelerate mode-of-action discovery in invertebrate-targeting compound development and help to refine mode-of-action categories.

Keywords: C. elegans; anthelmintics; computational ethology; pesticide; phenotypic screen.

PubMed Disclaimer

Conflict of interest statement

Research grant support was provided by the Biotechnology and Biological Sciences Research Council of the UK in partnership with Syngenta UK. AJF and PHH are employees of Syngenta UK. AEXB has consulted with Syngenta UK.

Figures

Figure 1
Figure 1. Insecticides affect phenotypes in multiple behavioral dimensions

  1. We image an entire 96‐well plate with a megapixel camera array with enough resolution and high enough frame rate to track, segment, and estimate the posture of C. elegans over time.

  2. All compound doses in the speed/tail curvature space, with points and lines showing the mean and standard deviation of biological dose replicates. On average, 12 biological replicates were collected per compound dose, together with 601 DMSO replicates, across at least 3 different tracking days for each condition. Several compounds, including the serotonin receptor antagonist mianserin (blue), glutamate‐gated chloride channel activator emamectin benzoate (purple), and vesicular acetylcholine transporter inhibitor SY1713 (red), have a strong effect on the worms' behavioral phenotype. They can be distinguished from the DMSO control (black) and from each other based on speed and tail curvature alone. Not all compounds are well separated in these two dimensions (gray points). Inset images are samples that show postural differences.

  3. Sample worm skeletons over time show the effect of the compounds highlighted in (B) on motion.

  4. Number of features significantly different from the DMSO control at a false discovery rate of 1% for each compound, grouped by mode of action. The pre‐stimulus, blue light stimulus, and post‐stimulus data are shown separately (a total of 3,020 features are tested for each assay period). The percentage of significantly different features is highest for the blue light stimulus recording.

Figure 2
Figure 2. Clustering and dose response of behavioral fingerprints

  1. Hierarchical clustering of behavioral fingerprints highlights structure in the responses to different compounds. Each row of the heat map represents the mean dose fingerprint of a specific compound described by 256 pre‐selected features from each blue light condition. Clear clusters can be observed for some compound classes, e.g., AChE inhibitors, vAchT inhibitors, GluCl agonists, and mAchR agonists. Low doses and low potency compounds from different classes cluster together around the DMSO averages at the center‐top part of the heat map.

  2. Cluster purity as a function of the hierarchical cluster distance shows that the degree of mode‐of‐action clustering (red) is greater than expected by chance for random clusters (gray).

  3. Compounds in the same class can have different dose–response curves. (upper) The three mitochondrial inhibitors cyazofamid, rotenone, and SY1048 all decrease angular velocity, but the concentration at which their effect is measurable is not conserved across compounds. (lower) Different spiroindolines affect body curvature differently. SY1786's dose–response curve is non‐monotonic. The central band and box limits show the median and quartiles of the distribution of the biological replicates for each compound dose (on average 12 wells per dose and 601 DMSO wells), while the whiskers extend to 1.5 IQRs beyond the lower and upper quartile. The P‐values reported in the legend represent the significance of the drug dose effect and were estimated using linear mixed models with tracking day as random effect and drug dose as fixed effect. The positions of these compounds in the heat map in (B) are marked using the color bar on the right side of the heat map.

Figure EV1
Figure EV1. The smoothing and balancing procedure using bootstrapped averages reduces the effect of outliers and balances the classes increasing classification accuracy
One of the sparsely populated classes (GluCl agonists with 5 compounds) and the most well‐populated class (AChE inhibitors with 10 compounds) are shown before and after the smoothing and balancing procedure in the 3 top PCA components.
Figure 3
Figure 3. Classifiers trained on behavioral fingerprints can predict the mode of action of unseen test compounds

  1. Toy data illustrating the potential benefit of normalization in correcting for potency differences within mode‐of‐action classes. Following normalization, each behavioral fingerprint exists on a hypersphere in the phenotype space regardless of effect size in the original space. Nonlinear dose–response curves will not collapse perfectly following normalization, which is a linear transformation.

  2. The confusion matrix obtained through cross‐validation for the best performing feature set (1,024 features) and logistic regression classifier following feature selection and hyperparameter tuning on the training data.

  3. The confusion matrix for the classifier trained in (B) applied to previously unseen test compounds without any further tuning.

  4. The novelty score assigned to novel test compounds with a mode of action not seen during training compared to the novelty score of compounds from the test set in (C). Novel compounds tend to have higher novelty scores than compounds from previously seen modes of action. The non‐novel compounds with high novelty scores include the two incorrectly classified test compounds (in red box).

Figure EV2
Figure EV2. Illustration of the effect of normalization with real data
The normalization of samples to unit L2 norm brings compounds with different potencies closer together and helps separate the classes in the phenotypic space.

  1. PCA of the data from 2 different classes with 3 compounds each before normalization (data simply standardized).

  2. PCA of the normalized data.

Figure 4
Figure 4. Mode of action can be resolved within compound classes

  1. The confusion matrix showing cross‐validation performance of a classifier trained to distinguish serotonin receptor antagonists from each other. Ritanserin is distinguishable from all other compounds, and the two structurally similar antidepressants (mianserin and methiothepin) are somewhat mutually confused by the classifier but are distinguishable from the non‐antidepressants.

  2. The confusion matrix for the mitochondrial inhibitors also shows some substructure: complex II and III inhibitors (cyazofamid, antimycin, cyenopyrafen) are phenotypically similar, and distinct from the complex I inhibitors.

Figure EV3
Figure EV3
Confusion matrices showing cross‐validation performance of a classifier trained to distinguish compounds of the same class from each other.
Figure EV4
Figure EV4. Randomization of columns using liquid handling robot
To minimize any position‐induced bias, we use an Opentrons liquid handling robot to randomly shuffle the position of the compounds in the imaging plates. We programmed the robot to keep a record of the randomized shuffle and use this log to create the correct well‐compound mapping.
Figure EV5
Figure EV5
Flow chart showing the data pre‐processing steps.
Figure EV6
Figure EV6
Flow chart showing the pipeline for the clustering of average drug doses.
Figure EV7
Figure EV7
Flow chart showing the steps followed for the tuning and training of the classifier using the training dataset and the prediction of the mode of action in the test set.
Figure EV8
Figure EV8
Flow chart showing the steps of the novelty detection method.
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

See all "Cited by" articles

References

    1. Anderson DJ, Perona P (2014) Toward a science of computational ethology. Neuron 84: 18–31 - PubMed
    1. Benjamini Y, Yekutieli D (2001) The control of the false discovery rate in multiple testing under dependency. Ann Stat 29: 1165–1188
    1. Berman GJ (2018) Measuring behavior across scales. BMC Biol 16: 23 - PMC - PubMed
    1. Berman GJ, Choi DM, Bialek W, Shaevitz JW (2014) Mapping the stereotyped behaviour of freely moving fruit flies. J R Soc Interface 11: 20140672 - PMC - PubMed
    1. Brenner S (1974) The genetics of Caenorhabditis elegans . Genetics 77: 71–94 - PMC - PubMed

Publication types