- 1Biological Sciences Platform, Sunnybrook Research Institute, Toronto, ON, Canada
- 2Division of Neuroscience, Sirio-Libanês Hospital, São Paulo, Brazil
- 3Division of Neurosurgery, Department of Surgery, University Health Network and University of Toronto, Toronto, ON, Canada
- 4Cerebral Imaging Centre, Douglas Mental Health University Institute, McGill University, Montreal, QC, Canada
- 5Department of Psychiatry, McGill University, Montreal, QC, Canada
- 6Department of Neurology, Division of Functional Neurosurgery of the Institute of Psychiatry, University of São Paulo, Medical School, São Paulo, Brazil
- 7LIM/23, Department of Psychiatry, University of São Paulo, Medical School, São Paulo, Brazil
- 8Department of Biological and Biomedical Engineering, McGill University, Montreal, QC, Canada
Aggressive behaviors comprise verbal and/or physical aggression directed toward oneself, others, or objects and are highly prevalent among psychiatric patients, especially patients diagnosed with autism spectrum disorder and severe intellectual disabilities. Some of these patients are considered refractory to treatment, and functional neurosurgery targeting the amygdala can result in widespread plastic brain changes that might reflect ceasing of some abnormal brain function, offering symptom alleviation. This study investigated cortical thickness changes in refractory aggressive behavior patients that were treated with bilateral amygdala ablation and compared to control patients presenting non-refractory aggressive behavior [three refractory and seven non-refractory patients, all males diagnosed with autism spectrum disorder (ASD) and intellectual disabilities]. The Overt Aggression Scale (OAS) was used to quantify behavior and magnetic resonance imaging was performed to investigate cortical thickness. Before surgery, both groups presented similar total OAS score, however refractory patients presented higher physical aggression against others. After surgery the refractory group showed 88% average reduction of aggressive behavior. Imaging analysis showed that while refractory patients present an overall reduction in cortical thickness compared to non-refractory patients across both timepoints, the local pattern of thickness difference found in areas of the neurocircuitry of aggressive behavior present before surgery is diminished and no longer detected after surgery. These results corroborate the hypotheses on induction of widespread neuronal plasticity following functional neurosurgical procedures resulting in modifications in brain morphology and improvement in behavior. Further studies are necessary to determine the underlying cause of these morphological changes and to better understand and improve treatment options.
Introduction
Aggressive behavior is a multifactorial behavioral syndrome that includes physical aggression and/or verbal aggression and can be directed at oneself (self-injury behavior), to others, and to objects (Corrigan et al., 1993; Yudofsky et al., 2012). It has a high prevalence in patients with Autism Spectrum Disorder and in patients with intellectual disabilities, being a leading cause for patient institutionalization (Brentani et al., 2013; Adler et al., 2015; Gouveia et al., 2019, 2020). Pharmacological therapies targeting the serotonergic and dopaminergic systems (e.g., typical and atypical antipsychotics, as well as selective serotonin reuptake inhibitors), anti-epileptic agents (e.g., valproic acid and carbamazepine), psychostimulants (e.g., amphetamine, methylphenidate), and a large variety of mood stabilizers (e.g., lithium) and tranquilizers (e.g., benzodiazepines) are commonly used to treat aggressive behavior in association with behavioral therapy (Morrissette and Stahl, 2014; Adler et al., 2015; Coleman et al., 2019). Alternative high dose mono and poly-pharmacy therapies can be offered when patients fail to present behavioral improvement with standard therapy (Morrissette and Stahl, 2014). However, these alternative therapies have been associated with symptom deterioration and their long-term side effects are poorly understood. Patients whose behavior does not improve after exhausting all available medication therapies are considered to have refractory aggressive behavior (Adler et al., 2015; Gouveia et al., 2021) and functional neurosurgical procedures may be considered to alleviate symptoms (Sturm et al., 2012; Gouveia et al., 2019, 2021).
The amygdala is a central structure of the neurocircuitry of aggressive behavior that also includes cortical areas (e.g., the orbitofrontal cortex, cingulate cortex) and subcortical structures (e.g., hypothalamus, nucleus accumbens, periaqueductal gray matter) (Siever, 2008; Blair, 2016; Gouveia et al., 2019). An imbalance in serotonin transmission in prefrontal cortical regions along with an increase in limbic activity (especially the amygdala) is believed to be the underlying cause for excessive aggressive behavior (Siever, 2008; Blair, 2016; Gouveia et al., 2019). Previous studies have shown that the bilateral ablation of the amygdala in patients with a wide range of psychiatric disorders associated with severe and/or refractory aggressive behavior can successfully reduce the frequency and/or severity of aggressive behavior and improve patients’ quality of life (Kiloh et al., 1974; Balasubramaniam and Kanaka, 1975; Gouveia et al., 2019, 2021). The marked reduction in aggressive behavior observed in these severe cases following the ablation of the amygdala will affect the neurocircuitry of aggressive behavior, and this circuit-wide restoration of functional equilibrium should be associated with measurable plastic changes (Siever, 2008). To address this hypothesis, in this study we investigated changes in cortical thickness in patients presenting refractory aggressive behavior who were treated with bilateral amygdala radio-frequency ablation for symptom alleviation as compared to patients presenting non-refractory aggressive behavior.
Materials and Methods
Subjects
A total of 10 patients with severe aggressive behavior hallmarked by severe self-injury behavior and extreme aggression toward objects and others were included in this study. All patients were males diagnosed with autism spectrum disorder (ASD) and severe intellectual disability, as qualitatively described by the Diagnostic and Statistical Manual of Mental Disorders, 5th Edition (DSM-5; American Psychiatric Association [APA], 2013). Tests of intelligence quotient (IQ) were not performed as patients were not collaborative, presented high motor agitation and had limited ability to communicate. Of those, three patients were considered to have refractory aggressive behavior (Gouveia et al., 2020, 2021), defined as persistent aggressive behavior despite previous trials of FDA approved antipsychotics targeting the core symptom (Adler et al., 2015), previous trials of single and poly-medication (including antidepressants, mood stabilizers, psycho-stimulants and anticonvulsants) (Morrissette and Stahl, 2014; Gouveia et al., 2019) and insufficient improvement following behavioral therapy. Refractory patients were treated with bilateral amygdala radio-frequency ablation to reduce aggressive behavior (for full description of the surgical procedure see Gouveia et al., 2021). The remaining seven patients were considered non-refractory, as they did not exhaust all possible medication therapies, and surgery was not considered for these patients. Table 1 provides detailed demographics. All procedures were approved by the Research Ethics Board (CAPPesq #742.331; CAAE #31828014.6.3001.5461), and written informed consent was obtained from the patients’ parents.
Questionnaire
The Overt Aggression Scale (OAS) was used to quantify aggressive behavior. The scale is divided in four domains to evaluate all aspects of the behavior and scored from 0 to 16, with 0 representing absence of the behavior and 4 representing the most extreme behavior for each domain (verbal aggression, physical aggression against objects, physical aggression against oneself, and physical aggression against other people) (Silver and Yudofsky, 1991; Yudofsky et al., 2012). No neuropsychological assessments were performed as patients presented with severe intellectual disability and were non-verbal. All questionnaires were reported by the patients’ parents and were performed concomitant with imaging magnetic resonance imaging (MRI) acquisition.
Magnetic Resonance Imaging Acquisition and Processing
Magnetic resonance imaging was performed with patients under deep sedation using a 1.5 Tesla MRI system (Magnetom Espree, Siemens, Germany). T1-weighted structural images were acquired with the following parameters: slice thickness 1.0 mm, no gap, voxel size 1 mm × 0.9 mm × 0.9 mm or 1 mm × 0.5 mm × 0.5 mm, TE/TR 5/300 ms, flip angle 45°, FOV 240 mm. For patients with refractory aggressive behavior, MRI was performed before (on the week preceding surgery for surgical planning) and after bilateral radio-frequency amygdala ablation (3 months for case 2 and 12 months for cases 1 and 3; Supplementary Figure 1B). Non-refractory patients were scanned once. Imaging acquisitions was restricted to the minimum as the cases included here present a higher risk of respiratory depression and death during sedation (Gouveia et al., 2020, 2021). Images were processed as previously described (Bedford et al., 2020; Gouveia et al., 2020). Briefly, images were non-uniformity corrected using the iterativeN4 from minc-bpipe-library1 and subsequently processed using the CIVET pipeline (v2.1.0).2 Images were non-linearly registered to MNI152, and tissue classification was performed using classification priors (gray and white matter and cerebrospinal fluid). Cortical thickness was computed for all patients, at each vertex based on the distance between gray and white matter surfaces using Laplace’s method and blurred with a 40 mm geodesic surface kernel. Quality control of CIVET outputs was performed for all patients. Supplementary Figure 1A shows examples of the CIVET-based segmentations for refractory and non-refractory patients.
Statistical Analysis
For cortical thickness analysis the RMINC (v.1.5.2.0)3 package in R (v.3.4.4)4 was used. The questionnaires’ measures were analyzed using Student’s t-test. Brain wide differences in cortical thickness were evaluated comparing refractory versus non-refractory patients at two timepoints: I. refractory before surgery versus controls and II. refractory after surgery versus controls. To this aim, linear models of thickness by group using individual mean cortical thickness as the covariate to correct for global differences were computed at each vertex. All analyses were corrected for multiple comparisons using False Discovery Rate (FDR) at the pFDRcor < 0.05 threshold.
Results
Only male patients were included in this study. Ages did not differ between patients with refractory (25.66 ± 4.50 years) and non-refractory (18.28 ± 4.60 years) aggressive behavior (t = 1.91, df = 3.39, p = 0.1415). No patient developed post-operative hydrocephalus in the 12-month follow-up and, aside from surgical scars and amygdala lesions, brain morphology and the MRI signal were considered normal and no alterations were detected in relation to the preoperative MRI. Before surgery, no significant differences were found between refractory and non-refractory patients in total OAS score or in the subscales of verbal aggression, physical aggression against objects, and physical aggression against self; however, refractory patients presented significantly higher levels of physical aggression against other people (Table 2). After surgery the refractory group showed a marked reduction in aggressive behavior, showing reduction in total OAS score (88% average) and in all subscales: verbal aggression (66% average), physical aggression against objects (91% average), physical aggression against oneself or himself (91% average), and physical aggression against other people (91% average; see Table 1 for individual improvement rate). Comparing refractory patients after surgery to the non-refractory patients, we observed significant differences in total OAS score and in the subscales of physical aggression against objects, physical aggression against self, and physical aggression against other people (Table 2). No significant difference was found in verbal aggression (Table 2).
Refractory patients presented an overall reduced cortical thickness compared to non-refractory subjects and the bilateral amygdala ablation did not change this morphological characteristic [Before surgery: left (t = −5.42, p < 0.001) and right (t = −5.61, p < 0.001) hemispheres; after surgery: left (t = −8.17, p < 0.001) and right (t = −4.76, p = 0.001) hemispheres). However, the vertex-wise relative cortical thickness analysis identified areas that showed significantly different cortical thickness (both increased/decreased) in refractory patients before surgery when compared to non-refractory patients (pFDRcor < 0.05; Figure 1). Local bilateral increases in cortical thickness were found in the anterior cingulate cortex and posterior cingulate cortex, medial prefrontal cortex, dorsolateral prefrontal cortex, intraparietal sulcus, and orbitofrontal cortex (Figure 1). Local bilateral reduction in cortical thickness was observed in the somatosensory cortex, superior temporal sulcus, premotor cortex, and motor cortex (Figure 1). The local differences diminish following surgery as no significant differences can any longer be found between refractory and non-refractory subjects (Figure 1). While the local thickness differences in the right dorsolateral prefrontal cortex seem to be an exception to that general trend, the fact that no vertex survives multiple comparison corrections suggests that these results may be spurious.
Figure 1. Cortical thickness analysis using CIVET. Cold colors show areas where CT is reduced in patients with refractory aggressive behavior compared with patients with non-refractory aggressive behavior. Warm colors show local cortical thickness increases comparing patients with refractory aggressive behavior to patients with non-refractory aggressive behavior. Please note that there were no significant differences between groups at follow up and colors in the images are uncorrected for multiple comparisons to better illustrate the extent and spatial pattern of the cortical thickness results.
Discussion
This study has investigated changes of cortical thickness and a number of different aggressive behaviors following treatment in patients with refractory aggressive behavior compared to non-refractory patients. Both groups presented similar levels of total aggressive behavior at baseline and after bilateral amygdala radio-frequency ablation the refractory group showed a marked reduction of aggressive behavior now significantly different from the non-refractory group in total OAS and all subscales except for verbal aggression. However, it should be noted that since all patients included in this study present severe intellectual disability and are non-verbal, the range of this subscale is markedly limited with possible scores 0 (absent) or 1 (makes loud noises/shouts angrily). All refractory patients showed score 1 before surgery and after the procedure only one continued with this behavior, with the other two patients scoring zero. Interestingly, this patient that maintained verbal aggressive behavior after surgery (i.e., case 2) presented an average lesion size of 50% of the total volume of the amygdala (reaching parts of the lateral, basal, and central nuclei), while the remaining patients presented almost complete amygdala ablation. The amygdala is a key structure of the limbic system responsible for integrating concurrent inputs from the sensory systems and projecting this information to distinct brain regions to regulate emotional responses. The lateral nucleus of the amygdala in the main area receiving and processing these inputs and along with the basal nucleus, projects highly integrated multimodal information to the central nucleus, the most prominent output area. The high average improvement of the refractory group following amygdala ablation is in line with previous studies showing high rates of success after functional neurosurgery of the amygdala for aggressive behavior (Lee et al., 1998; Langevin et al., 2016; Gouveia et al., 2019, 2021). It is believed that a reduced top-down inhibitory control of the prefrontal cortex over the amygdala would result in hyper-activation of this nucleus and consequent activation of hypothalamic and brainstem regions resulting in lower tolerance of stimuli and greater motor activity. Thus, reducing amygdala activity by neuromodulation treatments is sufficient to restore a homeostatic brain pattern (Siever, 2008; Blair, 2016; Gouveia et al., 2019).
We observed an overall reduced cortical thickness in refractory patients compared to non-refractory patients. It is not possible to determine, however, the underlying cause of this morphological difference, but it is plausible to assume that the chronic use of high dose mono and poly-pharmacy by the refractory group could result in morphometric brain changes, including alterations in cortical thickness. A further cortical thickness analysis of images acquired before surgery revealed significant differences in the pattern of local cortical thickness between refractory and non-refractory patients in areas related to the control of aggressive behavior, such as the orbitofrontal cortex, anterior and posterior cingulate cortices and medial prefrontal cortex (Siever, 2008; Gouveia et al., 2020). Interestingly, no significant local cortical thickness differences can be found any longer after surgery, which corroborates the theory that the functional neurosurgical procedure causes plastic changes and widespread modifications in brain morphology associated with the restoration of brain equilibrium and concomitant improvement in behavior. These characteristic changes are found in all three patients, suggesting that these changes occur within the first months (3 months is the shortest duration since the procedure in patient 2) and remain stable after.
This study provides the rare opportunity to investigate the mechanisms underlying symptom improvement following amygdala ablation in patients that are refractory to all known treatments, with focus on brain-wide cortical remodeling. Although these data allow true causal inferences to be tested, some limitations are to be noted. Refractory patients treated with amygdala ablation are extremely rare, and studying highly aggressive patients that are not collaborative in clinical settings is a very challenging task; nevertheless, the small sample size prevents a greater generalization of the findings. All patients included here required image acquisition to be performed under sedation; however, these patients present a high risk for respiratory depression and death during sedation as a result of the pharmacological regime used to control aggressive behavior. Thus, non-refractory patients (included in the control group) were scanned only once, and no longitudinal analysis could be performed, limiting our knowledge about possible cortical thickness changes in patients treated only with medication. In this sense, this study does not attempt to compare cortical thickness changes in patients subjected to distinct treatments for aggressive behavior but aims to demonstrate that the pattern of cortical thickness observed in non-refractory patients after amygdala lesion (when patients return to show beneficial responses to pharmacological therapy) is similar to the one found in patients that do respond to medication. It is important to highlight that only ASD patients were included in this study as it has been shown that this patient population presents different brain morphology than healthy controls with greater cortical thickness (Bedford et al., 2020).
This study showed evidence that bilateral amygdala radio-frequency ablation induces profound neuronal plasticity along the neurocircuitry of aggressive behavior resulting in brain morphological changes that could be integral to the behavioral improvement observed. Further studies are necessary to determine the underlying cause of these morphological changes and to better understand and improve treatment options.
Data Availability Statement
All data generated and analyzed during this study are available from the corresponding author upon reasonable request.
Ethics Statement
The studies involving human participants were reviewed and approved by the Comissão de Ética para Análise de Projetos de Pesquisa do Hospital das Clinicas da Faculdade de Medicina da Universidade de São Paulo (CAPPesq HC FMUSP). Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin (CAPPesq #742.331; CAAE #31828014.6.3001.5461).
Author Contributions
FG, JG, and RCRM designed the work. FG, RCRM, EF, RMCBM, and HB acquired the data. FG, JG, and GD analyzed the data. FG, JG, GD, and MC interpreted the data. FG and JG drafted the manuscript. GD, EF, RCRM, HB, MC, and RMCBM critically revised the manuscript. All authors approved the final version of the manuscript.
Funding
This work was supported by grants FAPESP #11/08575-7 (RCRM), #13/20602-5 (FG), and #17/10466-8 (FG) from the Brazilian Government. This research was enabled in part by support provided by SciNet (www.scinethpc.ca) and Compute Canada (www.computecanada.ca).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We wish to thank the research assistants and staff of Hospital Sirio-Libanes, Brazil, and of Hospital das Clinicas da Faculdade de Medicina da Universidade de São Paulo, Brazil.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnhum.2021.653631/full#supplementary-material
Supplementary Figure 1 | (A) Cortical thickness segmentation in patients treated with amygdala lesion (I,II) and in controls (III,IV). (B) Magnetic resonance imaging of patients treated with amygdala ablation (Cases 1–3, images I–III, respectively) before and after surgery in axial (top row) and coronal (bottom row) planes.
Footnotes
- ^ https://github.com/CobraLab/minc-bpipe-library
- ^ http://www.bic.mni.mcgill.ca/ServicesSoftware/CIVET
- ^ https://github.com/Mouse-Imaging-Centre/RMINC
- ^ https://www.r-project.org/
References
Adler, B. A., Wink, L. K., Early, M., Shaffer, R., Minshawi, N., McDougle, C. J., et al. (2015). Drug-refractory aggression, self-injurious behavior, and severe tantrums in autism spectrum disorders: a chart review study. Autism 19, 102–106. doi: 10.1177/1362361314524641
American Psychiatric Association [APA] (2013). Diagnostic and Statistical Manual of Mental Disorders (DSM-5®). Virginia: American Psychiatric Association.
Balasubramaniam, V., and Kanaka, T. S. (1975). Amygdalotomy and hypothalamotomy - a comparative study. Confin. Neurol. 37, 195–201. doi: 10.1159/000102740
Bedford, S. A., Park, M. T. M., Devenyi, G. A., Tullo, S., Germann, J., Patel, R., et al. (2020). Large-scale analyses of the relationship between sex, age and intelligence quotient heterogeneity and cortical morphometry in autism spectrum disorder. Mol. Psychiatry 25, 614–628. doi: 10.1038/s41380-019-0420-6
Blair, R. J. R. (2016). The neurobiology of impulsive aggression. J. Child Adolesc. Psychopharmacol. 26, 4–9.
Brentani, H., de Paula, C. S., Bordini, D., Rolim, D., Sato, F., Portolese, J., et al. (2013). Autism spectrum disorders: an overview on diagnosis and treatment. Rev. Bras. Psiquiatr. 35, S62–S72.
Coleman, D. M., Adams, J. B., Anderson, A. L., and Frye, R. E. (2019). Rating of the effectiveness of 26 psychiatric and seizure medications for autism spectrum disorder: results of a national survey. J. Child Adolesc. Psychopharmacol. 29, 107–123. doi: 10.1089/cap.2018.0121
Corrigan, P. W., Yudofsky, S. C., and Silver, J. M. (1993). Pharmacological and behavioral treatments for aggressive psychiatric inpatients. Hosp. Community Psychiatr. 44, 125–133. doi: 10.1176/ps.44.2.125
Gouveia, F. V., Germann, J., de Morais, R., Fonoff, E. T., Hamani, C., Alho, E. J., et al. (2021). Longitudinal changes after amygdala surgery for intractable aggressive behavior: clinical, imaging genetics, and deformation-based morphometry study-A case series. Neurosurgery 88, E158–E169.
Gouveia, F. V., Germann, J., Devenyi, G. A., Morais, R. M. C. B., Santos, A. P. M., Fonoff, E. T., et al. (2020). Refractoriness of aggressive behaviour to pharmacological treatment: cortical thickness analysis in autism spectrum disorder. BJPsych Open 6:e85.
Gouveia, F. V., Hamani, C., Fonoff, E. T., Brentani, H., Alho, E. J. L., de Morais, R. M. C. B., et al. (2019). Amygdala and hypothalamus: historical overview with focus on aggression. Neurosurgery 85, 11–30. doi: 10.1093/neuros/nyy635
Kiloh, L. G., Gye, R. S., Rushworth, R. G., Bell, D. S., and White, R. T. (1974). Stereotactic amygdaloidotomy for aggressive behaviour. J. Neurol. Neurosurg. Psychiatr. 37, 437–444. doi: 10.1136/jnnp.37.4.437
Langevin, J.-P., Chen, J. W. Y., Koek, R. J., Sultzer, D. L., Mandelkern, M. A., Schwartz, H. N., et al. (2016). Deep brain stimulation of the basolateral amygdala: targeting technique and electrodiagnostic findings. Brain Sci. 6:28. doi: 10.3390/brainsci6030028
Lee, G. P., Bechara, A., Adolphs, R., Arena, J., Meador, K. J., Loring, D. W., et al. (1998). Clinical and physiological effects of stereotaxic bilateral amygdalotomy for intractable aggression. J. Neuropsychiatry Clin. Neurosci. 10, 413–420. doi: 10.1176/jnp.10.4.413
Morrissette, D. A., and Stahl, S. M. (2014). Treating the violent patient with psychosis or impulsivity utilizing antipsychotic polypharmacy and high-dose monotherapy. CNS Spectr. 19, 439–448. doi: 10.1017/s1092852914000388
Siever, L. J. (2008). Neurobiology of aggression and violence. Am. J. Psychiatry 165, 429–442. doi: 10.1176/appi.ajp.2008.07111774
Silver, J. M., and Yudofsky, S. C. (1991). The overt aggression scale: overview and guiding principles. J. Neuropsychiatry Clin. Neurosci. 3, S22–S29.
Sturm, V., Fricke, O., Bührle, C. P., Lenartz, D., Maarouf, M., Treuer, H., et al. (2012). DBS in the basolateral amygdala improves symptoms of autism and related self-injurious behavior: a case report and hypothesis on the pathogenesis of the disorder. Front. Hum. Neurosci. 6:341. doi: 10.3389/fnhum.2012.00341
Keywords: aggressive behavior, amygdala, cortical thickness (CT), magnetic resonance imaging, neurosurgery, self-injury behavior
Citation: Gouveia FV, Germann J, Devenyi GA, Fonoff ET, Morais RMCB, Brentani H, Chakravarty MM and Martinez RCR (2021) Bilateral Amygdala Radio-Frequency Ablation for Refractory Aggressive Behavior Alters Local Cortical Thickness to a Pattern Found in Non-refractory Patients. Front. Hum. Neurosci. 15:653631. doi: 10.3389/fnhum.2021.653631
Received: 14 January 2021; Accepted: 03 May 2021;
Published: 09 June 2021.
Edited by:
Darin D. Dougherty, Massachusetts General Hospital and Harvard Medical School, United StatesReviewed by:
Hao Yan, Peking University, ChinaBrendon M. Nacewicz, School of Medicine and Public Health, University of Wisconsin-Madison, United States
Copyright © 2021 Gouveia, Germann, Devenyi, Fonoff, Morais, Brentani, Chakravarty and Martinez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Flavia Venetucci Gouveia, fvenetucci@gmail.com; Raquel C. R. Martinez, quelmartinez@yahoo.com.br
†These authors have contributed equally to this work and share first authorship