Skip to main content

ORIGINAL RESEARCH article

Front. Neurosci., 01 February 2024
Sec. Brain Imaging Methods
This article is part of the Research Topic Image Processing Method of Animal MRI and Its Translation to Study Brain Structure and Function Volume II View all 4 articles

From imaging to precision: low cost and accurate determination of stereotactic coordinates for brain surgery Sapajus apella using MRI

  • 1Institute of Biological Sciences, Federal University of Pará, Belém, Brazil
  • 2National Primate Center, Institute Evandro Chagas, Ananindeua, Brazil

The capuchin monkey (Sapajus apella), a New World monkey species, exhibits prominent characteristics that make it an ideal model for neuroscience research. These characteristics include its phylogenetic traits, telencephalization coefficient, anatomical structures and pathways, genetic profile, immune responses, cognitive abilities, and complex behavioral repertoires. Traditionally, methodologies for stereotactic neurosurgery in research models have relied on the use of brain atlases. However, this approach can lead to errors due to the considerable variation in brain size and shape among individual monkeys. To address this issue, we developed a protocol for deriving individual coordinates for each monkey using a straightforward and relatively inexpensive method involving MRI imaging. Our protocol utilizes a specially designed, 3D-printed stereotactic head-holder that is safe to use with an MR magnet, non-invasive placement of fiducial markers, and post-processing with open-source software. This approach enhances MRI data visualization, improves anatomical targeting, and refines the design of neurosurgical experiments. Our technique could also prove beneficial in other areas of neuroscience research that require accurate calculation of stereotaxic coordinates. Furthermore, it could be useful for other nonhuman primate species for which brain atlases are typically unavailable.

1 Introduction

The refinement of neurosurgery techniques for the manipulation of cortical and subcortical targets is crucial for obtaining reliable and reproducible results. This is fundamental to neurology and neuroscience research. Non-human primates (NHPs) serve as invaluable animal models for predicting accurate preclinical results in biomedical research. Despite representing only a small proportion of vertebrate animals, many species are utilized in neuroscience research. Their value as predictive models are further enhanced due to their sophisticated behavior, and the anatomical and functional similarities of their nervous systems to humans (Chudasama and Robbins, 2006; Arnason and Clausen, 2016) (see Table 1).

Table 1
www.frontiersin.org

Table 1. Comparison between different species of non-human primates used in neuroscience research.

Stereotactic neurosurgery methodologies and neurophysiological microelectrode recordings in neuroscience research typically rely on brain atlases. These atlases are based on the assumption that morphological features are consistent across individuals relative to cranial landmarks such as the bregma, interaural line, and infra-orbital ridges. Using a stereotactic device, the atlas is created from a few subjects’ ex-vivo brain histology data to map and target specific brain regions by coordinates (Palazzi et al., 2008; Hardman and Ashwell, 2012). This approach is quite successful in rodents due to low intersubject variability, which across species, is reportedly less than 1 mm (Palazzi et al., 2008; Paxinos and Franklin, 2019). However, in NHPs, the variability of brain volumes is approximately 5-fold larger than in laboratory rodents and can be largely attributed to high intersubject variability (Ose et al., 2022). This significant variability in brain morphology and size presents a unique challenge for the use of brain atlases (Pereira-Pedro et al., 2017; Sansalone et al., 2020). Given the considerable variability exhibited within species, quantitative 3D morphometric analysis showed robust results in geometric allometry brain size and shape variation not associated with size in two New World Monkeys genera. Therefore, it may be inappropriate to use mean values to represent an entire species (Marroig, 2007).

Moreover, brain atlases for NHPs are largely available for the species most commonly used in neuroscience research, such as Rhesus monkeys (Macaca mulatta), Cynomolgus (Macaca fascicularis), and Marmosets (Callitrix jacchus). However, there are scarce references for species that are not among these, such as the robust capuchin monkeys (Sapajus apella), New World species with prominent characteristics for the study of complex neurological functions (Muniz et al., 2021). Species of the genus Sapajus have been successfully used as good models to elucidate cortical and subcortical anatomical connections (Sousa et al., 1991; Rosa et al., 1993; Adams et al., 2000). To the best of our knowledge there is only one atlas available for Sapajus apella, published in 1968 (Manocha et al., 1968). At that time, the genus Cebus comprised only four species: Cebus albifrons, Cebus olivaceus, Cebus capucinus, and Cebus apella. These were later separated into two distinct clades based on evidence from morphological, phylogenetic, and biogeographic studies, giving rise to two genera: the gracile (untufted) capuchins, representing the genus Cebus, and the robust (tufted) capuchins, representing the genus Sapajus (Alfaro et al., 2012).

Robust capuchin monkeys (Sapajus apella) are found in South America and are endemic to the Amazonian region. They represent one of the most widespread primates in the Neotropics. They are not an endangered species and breed well in captivity, easily adapting to human contact (Ryland et al., 1997). Although not traditionally used in neuroscience research, they serve as excellent experimental models due to their ease of handling, adaptability, and reproducibility (Table 2). They exhibit intricated neuronal circuitry, a high encephalization rate, a strong tendency toward exploitation and manipulation, complex social behaviors, and cognitive abilities, such as tool usage (Visalberghi, 1993). These characteristics, along with other cost–benefit related factors, make them a very interesting research model for both basic and applied studies of the nervous system (Muniz et al., 2021).

Table 2
www.frontiersin.org

Table 2. The number of articles using Sapajus apella as an animal model in neuroscience research per decade.

To address the limitations of neurosurgical procedures for NHP species that are not often used in research, we present a low-cost and easily attainable method in this study. This method uses Magnetic Resonance Imaging (MRI) to accurately determine the precise coordinates of brain areas in Sapajus apella. The MRI sessions targeted subcortical areas, but this approach can be applied to a range of procedures that require stereotactic localization of both cortical and subcortical structures in different NHP species. The methodologies and technologies we have developed are highly pertinent in regions where resources are more constrained. Our protocol leverages locally available resources and is tailored to the specific ecological and socio-economic context of Brazil.

2 Materials and methods

2.1 Animal handling and ethics

The participants were seven adult male capuchin monkeys (Sapajus apella) from a colony at the National Primates Center primate facilities in Ananindeua, Pará, Brazil. The animals were 21.14 ± 5.17 years old and weighed 4.3 ± 0.99 kg (Table 3). They were housed in pairs in standard cages, each measuring 2.5 × 2.0 × 2.5 m. Each cage had a communication window that allowed free movement in a total space of 2.5 × 4.0 × 2.5 m. The cages were regularly enriched with toys and objects to facilitate climbing and were cleaned regularly. The lighting followed a regular day/night cycle provided by natural light (latitude 01° 21′ 56” S; longitude 48° 22′ 20″). The animals’ diet consisted of specific chow for laboratory animals, specifically for non-human primates, and was offered daily. Fruits and natural juice were also part of the diet. Water was available ad libitum. All procedures were performed according to the guidelines of Directive 2010/63/EU of the European Union and under the approval of the Ethics Committee for the Use of Animals from the Evandro Chagas Institute (CEUA/IEC), protocol numbers 45/2016 and 37/2018.

Table 3
www.frontiersin.org

Table 3. Characteristics of the experimental subjects.

2.2 Preparation for imaging procedures

The imaging procedures were performed with the animals anesthetized by using 9 mg/kg of pethidine and 1.2 mg/kg of midazolam, followed by propofol at 2 mg/kg maintained by infusion at the rate of 0.4 mg/kg/min (Galante et al., 2019). The animals were transported to the imaging facilities in portable stainless steel containment cages (0.80 m × 0.90 m × 0.80 m), with constant anesthetic monitoring by a veterinarian and two handling technicians. For proper positioning of the animals to enter the MRI machine, they were placed in a non-metallic (polylactic acid plastic/PLA) stereotactic apparatus (a 3D printed copy of the original stereotactic apparatus—Narishige®, see Figures 1, 2 for some details). Correct placement of ear bars into the acoustic meatus, and placement of the mouth adaptor and orbital bars were performed to guarantee alignment of the orbitomeatal plane. Small fish oil capsules were glued to the head at the locations of the infraorbital foramen, and the tip of each ear bar was coated with cotton embedded in a fish oil solution, to increase signal in these areas and facilitate further localization. The animal in the apparatus was then positioned on an MRI bed-sled that fit into the bore of the MRI machine. An 8-channel receive coil was placed over this stereotaxic system (see Figure 3A). To enhance flexibility and ensure animal safety, the orbital bars and the edges of the ear bar were coated with industrial silicone (Figure 3B). After the image acquisition, the anesthetic supply was suspended, and the animal was monitored until it fully recovered.

Figure 1
www.frontiersin.org

Figure 1. (A) Diagram of the 3D-printed, non-metallic stereotactic apparatus. The ear bars (1), orbital bars (2), and mouth adaptor (3) are similar to those used in the surgical apparatus. (B) Schematic drawing of the animal positioned in the apparatus for image acquisition. The ear bars are positioned in the right and left external acoustic meatus, the eyepieces are on the right and left inferior orbital margins, and the mouth adaptor is securely fixed to the palatus.

Figure 2
www.frontiersin.org

Figure 2. Non-metallic apparatus for MRI acquisition. The frame, depicted in dark gray, was constructed from utility engineering materials, including PVC pipes and elbows. The components in lighter gray were produced using a Sethi3D S3 printer, which employs Fused Deposition Modeling (FDM) technology. Polylactic acid (PLA) was chosen for printing due to its high resolution, structural rigidity, and minimal shrinkage. The 3D models of each component can be accessed in the. STL file (refer to Supplementary material S2). Assembly was made easy by straightforward hardware fittings. The alignment of the stereotaxic device was executed with precision using a square and level.

Figure 3
www.frontiersin.org

Figure 3. (A) This image showcases the stereotactic apparatus equipped with a standard receive coil for MRI examinations. For this study, a phantom was used to enhance the control and precision of the image acquisition. (B) Close-up of the stereotaxic apparatus with a phantom near the ears and orbital bars coated with industrial silicone.

2.3 Assessment of MRI compatibility

To ensure that the stereotactic device does not cause image loss, an MRI was also acquired using a phantom framed and frameless. Phantom volumes were calculated using rendering and parameters calculations of SLICER 3D© software (see Supplementary material S3). To define the SNR (Signal-to-Noise Ratio) values of each image, the average signal intensity was calculated and then divided by the standard deviation of noise reference. The average signal intensity was measured using ROI with a 10-point square in the center of each image. The reference noise, ROI with a 10-point square in the top left corner of each image.

2.4 Image acquisition and data processing

MR anatomical images were acquired in a SIGNA Creator 1,5 T using a 3D T1 Cube (voxel/pixel ratio = 1). Whole-brain images were acquired in a 3D volume with a 256x256x256mm matrix; a slice thickness of 0.6 mm, NEX = 1, TR = 500 ms, and a minimum echo time. Images were uploaded in DICOM (Digital Imaging and Communications in Medicine) format. Post-processing was performed in SLICER 3D© software. This software provides a tridimensional reconstruction of the animal’s brain through MRI sections uploaded in DICOM files. These reconstructions were used for brain volume rendering and calculation, followed by the placement of anatomical (infraorbital foramina, acoustic meatus) and fiducial markers (middle point of the interaural line) that were used as references for aligning the stereotactic planes. For a step-by-step description of the procedures on the software, see Supplementary material S3. The anteroposterior zero plane (AP-0) was represented by a straight line passing through the width of the brain at the center of the external auditory meatus. The origin (zero point) of the stereotactic coordinates corresponded to the midpoint between the two ear-bars, with the horizontal plane parallel to the orbitomeatal (OM) plane as determined by the two ear-bars and two eye-bars. The axial, coronal, and sagittal planes were then aligned perpendicularly at the zero point (Figure 4).

Figure 4
www.frontiersin.org

Figure 4. Image data processing for determination of stereotactic coordinates. (A) Image rendering and placement of the reference planes. (B–D) Placement of anatomical and fiducial markers and alignment of the coronal, axial, and sagittal planned perpendicularly and aligned with the reference points. (B) Coronal image passing through the intermeatus line (green line, y). (C) Axial image passing through the transverse axis (red line, z). (D) Sagittal image passing through the longitudinal fissure (yellow line, x) (A, Accustic meatus; IF, infraorbital foramen; P0, the origin).

2.5 Determination of the stereotaxic coordinates

After the image rendering and placement of the reference planes, 14 potential surgical targets were placed in easily evident and well-defined brain structures: corpus callosum (genus – G, and splenium – S), anterior horn of lateral ventricles (left – A1 and right – A2), posterior horn of the lateral ventricles (left – P1 and right – P2), substantia nigra (four points distributed laterally through the nucleus length – SN1-4) and caudate nucleus (four points distributed anteroposteriorly through the nucleus length – C1-4) (Figures 5A,B). The coordinates of each target were calculated automatically by the 3D Slicer software. All coordinates calculated by the program were visually confirmed by two experimenters who measured the target points manually using the software’s rule (Figure 5C).

Figure 5
www.frontiersin.org

Figure 5. Determination of the stereotaxic coordinates. (A) Placement of surgical targets through basal ganglia structures (SN, substantia nigra; CD, caudate nucleus). (B) Placement of surgical targets through lateral ventricles and corpus callosum (HLV, horns of lateral ventricles; CC, corpus callosum; S, splenium; G, genus). (C) Visual confirmation of coordinates determination. In this case, SN1 coordinate was manually measured as 5,45 mm in medio-lateral (ML) which confirms the surgical coordinates previously calculated.

2.6 Statistical analysis

We performed statistical analysis using GraphPad Prism1 8.0 Software. First, we tested the normality (Kolmogorov–Smirnov) and homogeneity (Levine test) of the data. A parametric statistical analysis using ANOVA (One-way, α = 0.05, Tukey post-test) was performed for each axis (AP, DV, ML) to detect differences between coordinates. Pearson correlation coefficients were calculated to verify possible correlations between stereotaxic coordinates, brain volumes, and cephalic perimeters. The coefficient of variation was calculated for axis coordinates and brain volumes. The confidence interval was set to 95% (p < 0.05).

For MRI compatibility comparisons, we tested the normality test (Shapiro–Wilk). A non-parametric statistical analysis using the Mann–Whitney test (U, α = 0.05) was performed for SNRs values to detect any device susceptibility on MRI acquisition.

3 Results

3.1 Assessment of MRI compatibility

The same MRI acquisition protocol was conducted with a phantom in the stereotactic system and solely phantom to assess MRI compatibility comparisons (Table 4 and Supplementary material S1). There is no significant difference in SNR values (p = 1). This confirmed that there was no interference in MRI acquisition, which validates the use of the machined non-metallic stereotactic apparatus. Note there is no difference induced by the device (Figure 6).

Table 4
www.frontiersin.org

Table 4. Volume and average SNRs for the scans using the phantom with and without the stereotactic frame.

Figure 6
www.frontiersin.org

Figure 6. Temporal SNR for phantom frameless, and phantom framed. The same protocol was used for both conditions.

3.2 Normality test

The normality test identified one of the animals (AM-AOR/Velhinho) as an outlier. Given that the MRI images taken for this animal showed a significant movement artifact, which made it difficult to place the markers and surgical targets, we excluded this subject from the sample. The results of the normality tests (Kolmogorov–Smirnov) are presented in Table 5.

Table 5
www.frontiersin.org

Table 5. Normality test results for all coordinates in each of the three axes (mm).

3.3 Individual calculation of stereotaxic coordinates

The stereotactic coordinate calculations were performed individually according to the procedures described above. The individual coordinates for each of the surgical targets, and for the three anatomical axes are presented in Table 6. ANOVA tests comparing the coordinates for all animals in each target did not show significant differences between the measurements obtained for the subjects. The averages were significantly different enough to reject the null hypothesis of coordinates equality between monkeys at each axis (AP [p = 0.98]; DV [p = 0.19]; ML [p = 0.76]).

Table 6
www.frontiersin.org

Table 6. Individual coordinate values were calculated for each surgical target in each subject (mm).

3.4 Surgical confirmation of the coordinates

The surgical planning yielded robust results in terms of experimental accuracy. The surgical procedure, aimed at the basal ganglia, was informed by coordinates determined through MRI Slicer 3D processing and demonstrated efficiency. The conventional apparatus Narishige® was used for surgical manipulation. Our surgical protocol involved unilateral injection of 6-hydroxydopamine into the substantia nigra and caudate nucleus to induce a reduction in dopaminergic neurons (see Figure 7).

Figure 7
www.frontiersin.org

Figure 7. This photomicrograph depicts the mesencephalic region (substantia nigra) of a Sapajus apella subject, highlighting Tyrosine Hydroxylase (TH) immunohistochemistry. (A) At a lower magnification, the left side of the brain exhibits a lesion induced by our surgical protocol, which involved the injection of a neurotoxin. In contrast, the right side remained unlesioned, serving as a control. (B) A higher magnification of the area shown in Figure A indicates a notable reduction in TH-positive neurons on the lesioned left side. (C) A higher magnification of the unlesioned right side, as shown in Figure A, displays numerous TH-positive neurons. This data has not yet been published.

4 Correlation analysis of brain coordinates and volumes

We observed no significant correlation between brain volumes and stereotactic coordinates on either axis. This suggests that the variability in coordinates does not mirror changes in overall brain volume. In other words, variations in brain volume might be indicative of nonlinear differences in skull size and shape. To demonstrate this, we arranged the data on brain volumes and coordinates in descending order and emphasized the coordinates corresponding to each animal’s brain volume. For clearer visualization, we assigned different colors to each subject’s brain, creating a coordinate matrix for each animal (see Figure 8).

Figure 8
www.frontiersin.org

Figure 8. This color matrix displays individual coordinates in descending order. It is important to note that for each subject, the descending order of brain volume does not necessarily correspond to the arrangement of the coordinates in the same manner. This discrepancy reveals nonlinear variances in brain shapes, as illustrated by the 3D renderings of entire brain volumes from each MRI scan (right).

4.1 Coefficient of variation

The CV for brain volumes was 9.4%. The COV for all coordinates in different axes can be seen in Figure 9. This suggests a non-linear variation, with a higher COV value observed for AP coordinates.

Figure 9
www.frontiersin.org

Figure 9. (A) Brain volume rendering in descending order (1–6). Brain shape overlapping comparison (7). (B) COV distribution of the surgical targets in anteroposterior (AP), dorsoventral (DV), and mediolateral (ML) axes.

5 Discussion

In our study, we described a technique in which MRI imaging was used to obtain stereotactic coordinates for different surgical targets in Sapajus apella, a gyrencephalic New World species that is a promising model for many approaches in neuroscience research. The method described here represents an important experimental refinement when compared to conventional methods (see MRI processing guide in Supplementary material S1). In primate experiments, atlas-based stereotaxy is commonly used. Errors in targeting are generally a result of either morphologic differences between the individual macaque brains, or human errors in identifying the brain landmarks (Subramanian et al., 2005). These errors in targeting and the consequent misplaced targeting are rarely reported but should be included in future publications.

Moreover, the availability of brain atlases for different species of non-human primates is scarce, and the available atlases are either entirely based on a single animal or lack adequate 3D data (Paxinos et al., 2000). As mentioned, the only brain atlas available for capuchin monkeys was published in 1968 (Manocha et al., 1968) for Cebus apella. At that time, the genus Cebus included Tufted and Untufted Capuchin monkeys, which posteriorly. These were later separated into two different genera after phylogenetic studies have confirmed a huge diversity, justifying the creation of a new taxonomic division, the genus Sapajus (Lynch Alfaro et al., 2012). Additionally, the aforementioned publication does not cover the integrity of the brain regions and displays only images of the basal telencephalon in the coronal plane without coverage of more superior or cortical structures.

Therefore, the need for the development of our protocol primarily stems from the lack of available stereotactic atlases for Sapajus apella. The literature on the use of Sapajus apella as an animal model for experimental neurosurgeries is quite scarce. Despite being a promising model for the study of both normal and abnormal nervous systems (Hecht et al., 2021; Muniz et al., 2021; Miss et al., 2022), capuchin monkeys are not among the most frequently used species in non-human primate biomedical research. This is likely due to practical reasons that make other NHP species more suitable for laboratory use (Bert et al., 2012; Lankau et al., 2014).

Unlike rodents, whose brain structures exhibit relatively constant size and location values, the brains of NHPs can be highly variable (Bentley et al., 2018). Primates are altricial animals, experiencing a significant postnatal brain growth spurt. Therefore, variations in environmental factors such as maternal care, food, and environmental stimuli can induce a wide range of phenotypes, reflecting on the development of the size and shape of the skull (Ikonomidou et al., 2001; De Vareilles et al., 2022) (see Table 6; Figure 8). Thus, the use of MRI for targeting significantly improves surgical accuracy as it addresses the issue of intraspecies variability in brain structure.

Based on the classical Cartesian plane, frame-based systems for neurosurgery have the advantage of proven clinical accuracy and a high degree of mechanical stability, making them the gold standard. However, there are many frameless methods available that appear to offer flexibility (Fukaya et al., 2010; Patil, 2010; Sutherland et al., 2013). Frameless systems demand greater complexity and material expenditure. In addition to requiring precise hand coordination, natural tremors may influence the error margin and surgical results (Bradac et al., 2017).

Convection-enhanced delivery (CED) is a frameless system of drug delivery to the brain through intraparenchymal microcatheters with broader achievement using fewer infusions (Bradac et al., 2017; Yazdan-Shahmorad et al., 2018). In addition to the need for compatible and expensive materials for implantation in the parenchyma, it constitutes an invasive surgical method with the potential for encephalitic states (Khateeb et al., 2019). Surgical methods must be efficient and minimize animal risk. MRI-guided surgeries are preferable to intraoperative plain film roentgenography due to associated radiation exposure issues (Damilakis et al., 2010). Beyond radiation safety, MRI-guided surgeries, as a preoperative refinement, enable virtual fiducial placement on software platforms for virtual surgical planning, which is crucial for minimizing the error rate (Farrell et al., 2014).

Our method, which uses MRI, provides superior soft tissue contrast and allows us to view the structures of the monkey brain in vivo. The procedure we followed for stereotactic calculations was both simple and inexpensive. Although we are aware of the existence of complex scripts and software reported in previous studies, they typically employ high-cost computational frameworks and involve several pre-processing steps (Subramanian et al., 2005; Frey et al., 2011), yet they offer similar accuracy to our method. Some protocols involve the use of fiducial markers attached to the skull or in subcortical regions, which are more invasive and can trigger unintended inflammatory responses (Asahi et al., 2003). Our method was designed to avoid invasive fixation on the stereotaxic frame when positioning the animal for MRI scanning, similar to the protocol used by Saunders (Saunders et al., 1990).

Other studies have previously applied similar protocols using Old World NHPs, primarily Rhesus (Macaca mulatta). For New World monkeys, the most commonly used species is the common marmoset (Callithrix jacchus), a New World species with a lissencephalic brain (Hikishima et al., 2011). An interesting study by Seidlitz (Seidlitz et al., 2018) provided a population-averaged MRI standard template with anatomical references for NHP brains, taking into account individual asymmetry. This included surface reconstructions and transformations to previously published digital brain atlases for Rhesus. We did not find a similar study for Sapajus apella.

6 Error margin

Although we applied a highly developed imaging method, the protocols described here are not without potential errors. We calculated a small error margin in our study by comparing the software-generated coordinates with visually controlled coordinates, which involved confirming the coordinates using the software’s manual tool. Based on these comparisons, the error margin for each coordinate in the three axes was calculated by taking the average of the differences between the software-calculated and visually-calculated coordinates, along with their standard deviation. A mean error of 0.02 ± 0.01 in the anterior–posterior, 0.09 ± 0.11 in the dorsoventral, and 0.01 ± 0.03 mm in the mediolateral was observed (see Supplementary material S1). This can occur due to factors such as MRI resolution, potential uncontrolled movement artifacts during image acquisition, and the examiner’s performance during visual control. Nonetheless, compared to the size and volume of the targeted anatomical structures, the error margin obtained here is negligible and may not be associated with a surgical error. On the other hand, the intraspecies variation observed here ranged between 4.04 ± 1.03 in the anteroposterior, 3.88 ± 0.58 in the dorsoventral, and 1.42 ± 0.62 mm in the mediolateral (see Table 4), which represents a significant distance that could lead to a targeting error in neurosurgical approaches.

6.1 Correlation between brain volume and stereotaxic coordinates

It is known that brain growth is accentuated during early life development. After a highly dynamic early postnatal period, brain growth stabilizes in adulthood, when interindividual differences become well established. These differences result from varying and nonlinear growth patterns in different brain structures throughout the lifespan (Coupé et al., 2017; Danielsen et al., 2020), leading to interindividual differences in brain volume, shape, development, and proportion between brain structures.

In our study, we found no correlation between total brain volume and stereotaxic coordinates, which reinforces the fact that differences in brain development are nonlinear within species. It has been shown that the same subcortical structure may exhibit different stereotaxic locations throughout the lifespan, for example, the midsagittal area of the corpus callosum. Watson et al. (2022) observed sustained growth patterns of the corpus callosum in Sapajus apella monkeys in the splenium and genu regions. Regional differences are apparent in these growth patterns, with the genu increasing considerably during the first 6 years of life, and the splenium showing a greater increase between seven and 18 years (Vannucci et al., 2017). Considering our sample of adult monkeys, we observed higher variations for the coordinates P1, P2, and S, the most posterior structures analyzed in this study, which may reflect the patterns reported by Watson et al. (see Figure 9).

Based on the ratio of AC–PC (anterior commissure to posterior commissure) analysis in selected groups of humans and NHPs, there are higher variances in the AP axis. The mean ± SD for AC–PC distances were 28.3 ± 1.6 mm, 12.3 ± 0.8 mm, and 13.8 ± 0.7 mm for the human, Cynomolgus, and Rhesus groups, respectively (Fiandaca et al., 2011). This correlates with our results about the higher CV in AP coordinates and also supports the difficulty of localizing deep brain targets such as striatal brain structures (Yin et al., 2009).

Morphological brain variations can be related to different environmental contexts, which involve capuchins’ extractive foraging behaviors, including motor processing and spatial ability for prey capturing, tool usage, etc. (Janson and Boinski, 1992). Furthermore, Sapajus apella monkeys exhibit sex differences in behavior, with differences in cortical plasticity in brain organization (Hecht et al., 2021). It has been suggested that the morphology of the corpus callosum in capuchins is influenced by sex and handedness (Phillips et al., 2007), which could be linked to differences in cortical anatomy (Phillips and Sherwood, 2008).

Previous studies have calculated coefficients of variation (COV) of brain volume for different species, having reported a higher COV of brain volume in New World monkey species than that observed in rodents: 2.3% for mice (Ma et al., 2008), 3.2% for rats (Hasegawa et al., 2010), compared to 6.6% in marmosets (Hayashi et al., 2021). The significant primate’s COV correlates with the acceptable limits of the classical stereotactic approach. Our data suggest a volume brain COV of 9,4% for Sapajus apella adult population.

6.2 Advantages and limitations

This simple, cost-effective protocol with minimal pre-processing steps for primates’ MRI-based stereotactic calculations has the potential to assist experimenters working with different animal species in refining their neurosurgical methods and improving targeting (Figure 7). The most significant advantages of our protocol lie in the ease of determining any potential brain target, both cortical and subcortical, and in considering intraspecies anatomical differences. This method, which involves using a non-metallic replica of the surgical holder during imaging and referencing anatomical landmarks, can overcome potential errors induced by anatomical asymmetries of the subjects’ heads. For instance, substantial differences in the diameters of the external auditory meatus between the right and left sides can result in a discrepancy between two coordinate systems, because the midpoint of the ear bars does not correspond to the actual interaural zero. The validation by predicting the error margin met the required safety standards, leading us to believe that our protocol will have promising applications.

We believe that an open-source paradigm is a major trend in the development of medical imaging applications. 3D Slicer has the advantage of having an intuitive work interface, an internal DICOM image convertor, and the possibility of filing results through a scene system that can make changes on previous reconstructions, besides combining different scenes to compare studies.

Despite its simplicity and consistency, there are some limitations to this alternative protocol. First, there are technical limitations such as the quality of the MRI and post-processing software tools. For instance, the manual ruler for visual control depends on the operator, and the positioning of fiducial markers relies on anatomical knowledge (Gonzalo Domínguez et al., 2016). The resolution of the MRI must be considered and can be improved by adjusting the scanning sequence and extending the scanning time. Another significant limitation lies in determining the coordinates of areas where cytoarchitectural criteria are crucial for defining boundaries. For example, many neocortical areas cannot be confidently identified by MRI imaging, which only shows the differences between gray and white matter.

7 Conclusion

This technique addresses the issue of significant intraspecies variability in brain size and shape, a major source of bias in atlas-based studies for stereotactic brain surgery in NHPs. Furthermore, for species that are infrequently used in neuroscience research, the lack of available brain atlases can be overcome. The use of open-source and user-friendly software for image processing makes this method cost-effective and readily accessible for use in various experimental settings. The adoption of anatomical landmarks, along with the pre-placement of the animal in a stereotactic frame for imaging, may help to avoid stereotactic errors induced by the anatomical asymmetry of the animal’s head. Additionally, preoperative MRI scanning can be beneficial in identifying preexisting brain abnormalities that meet exclusion criteria for behavioral or neuroanatomic studies. In summary, our method may be highly beneficial to areas of neuroscience research that depend on accurate stereotaxic apparatus, with minimal procedural complexity.

We hope our protocol serves as a model for similar environments, proving that significant scientific contributions can emerge from regions traditionally underrepresented in global research landscapes. This aligns with the growing recognition of the importance of contextually relevant research in developing countries, which can address local challenges and contribute to the global scientific community.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The animal study was approved by Ethics Committee for the Use of Animals from the Evandro Chagas Institute (CEUA/IEC). The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

LP: Conceptualization, Data curation, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing. LL: Data curation, Formal analysis, Software, Writing – review & editing. JM: Conceptualization, Project administration, Writing – review & editing. CB: Data curation, Methodology, Software, Writing – review & editing. BG: Data curation, Writing – review & editing. LK: Conceptualization, Funding acquisition, Supervision, Writing – review & editing, Writing – original draft.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by the Horizon 2020 Research and Innovation Program of the European Union (n° 767092). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. publications costs by PROPESP/UFPA (PAPQ).

Acknowledgments

We would like to express our gratitude to the staff of the MEDICAL clinic in Castanhal-PA. We extend our thanks to the technicians of the National Center of Primates in Ananindeua-PA for their skillful assistance with MRI imaging acquisition, and collaborators of Eduardo Oswaldo-Cruz Neurophysiology Laboratory of the Federal University of Pará (UFPA). We also thank the M.sc Arthur Gonsales da Silva for the SNRs measurements script.

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.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2024.1324669/full#supplementary-material

References

Adams, M. M., Hof, P. R., Gattass, R., Webster, M. J., and Ungerleider, L. G. (2000). Visual cortical projections and chemoarchitecture of macaque monkey pulvinar. J. Comp. Neurol. 419, 377–393. doi: 10.1002/(SICI)1096-9861(20000410)419:3<377::AID-CNE9>3.0.CO;2-E

Crossref Full Text | Google Scholar

Alfaro, J. W. L., Silva Jr, J. D. S. E., and Rylands, A. B. (2012). How different are robust and gracile capuchin monkeys? An argument for the use of S apajus and C ebus. Am. J. Primatol. 74, 273–286. doi: 10.1002/ajp.22007

PubMed Abstract | Crossref Full Text | Google Scholar

Arnason, G., and Clausen, J. (2016). On balance: weighing harms and benefits in fundamental neurological research using nonhuman primates. Med. Health Care Philos. 19, 229–237. doi: 10.1007/s11019-015-9663-4

PubMed Abstract | Crossref Full Text | Google Scholar

Asahi, T., Tamura, R., Eifuku, S., Hayashi, N., Endo, S., Nishijo, H., et al. (2003). A method for accurate determination of stereotaxic coordinates in single-unit recording studies in monkeys by high-resolution three-dimensional magnetic resonance imaging. Neurosci. Res. 47, 255–260. doi: 10.1016/S0168-0102(03)00202-5

PubMed Abstract | Crossref Full Text | Google Scholar

Bentley, J. N., Khalsa, S. S., Kobylarek, M., Schroeder, K. E., Chen, K., Bergin, I. L., et al. (2018). A simple, inexpensive method for subcortical stereotactic targeting in nonhuman primates. J. Neurosci. Methods 305, 89–97. doi: 10.1016/j.jneumeth.2018.05.007

PubMed Abstract | Crossref Full Text | Google Scholar

Bert, A., Abbott, D. H., Nakamura, K., and Fuchs, E. (2012). The marmoset monkey: a multi-purpose preclinical and translational model of human biology and disease. Drug Discov. Today 17, 1160–1165. doi: 10.1016/j.drudis.2012.06.009

Crossref Full Text | Google Scholar

Bradac, O., Steklacova, A., Nebrenska, K., Vrana, J., De Lacy, P., and Benes, V. (2017). Accuracy of Vario guide frameless stereotactic system against frame-based stereotaxy: prospective, randomized, single-center study. World Neurosurg. 104, 831–840. doi: 10.1016/j.wneu.2017.04.104

PubMed Abstract | Crossref Full Text | Google Scholar

Chudasama, Y., and Robbins, T. (2006). Functions of frontostriatal systems in cognition: comparative neuropsychopharmacological studies in rats, monkeys and humans. Biol. Psychol. 73, 19–38. doi: 10.1016/j.biopsycho.2006.01.005

PubMed Abstract | Crossref Full Text | Google Scholar

Coupé, P., Catheline, G., Lanuza, E., Manjón, J. V., and Initiative, A. S. D. N. (2017). Towards a unified analysis of brain maturation and aging across the entire lifespan: a MRI analysis. Hum. Brain Mapp. 38, 5501–5518. doi: 10.1002/hbm.23743

PubMed Abstract | Crossref Full Text | Google Scholar

Damilakis, J., Adams, J. E., Guglielmi, G., and Link, T. M. (2010). Radiation exposure in X-ray-based imaging techniques used in osteoporosis. Eur. Radiol. 20, 2707–2714. doi: 10.1007/s00330-010-1845-0

PubMed Abstract | Crossref Full Text | Google Scholar

Danielsen, V. M., Vidal-Piñeiro, D., Mowinckel, A. M., Sederevicius, D., Fjell, A. M., Walhovd, K. B., et al. (2020). Lifespan trajectories of relative corpus callosum thickness: regional differences and cognitive relevance. Cortex 130, 127–141. doi: 10.1016/j.cortex.2020.05.020

PubMed Abstract | Crossref Full Text | Google Scholar

De Vareilles, H., Rivière, D., Sun, Z.-Y., Fischer, C., Leroy, F., Neumane, S., et al. (2022). Shape variability of the central sulcus in the developing brain: a longitudinal descriptive and predictive study in preterm infants. Neuro Image 251:118837. doi: 10.1016/j.neuroimage.2021.118837

Crossref Full Text | Google Scholar

Farrell, B. B., Franco, P. B., and Tucker, M. R. (2014). Virtual surgical planning in orthognathic surgery. Oral Maxillofac Surg Clin 26, 459–473. doi: 10.1016/j.coms.2014.08.011

Crossref Full Text | Google Scholar

Fiandaca, M. S., Salegio, E. A., Yin, D., Richardson, R. M., Valles, F. E., Larson, P. S., et al. (2011). Human/nonhuman primate AC–PC ratio—considerations for translational brain measurements. J. Neurosci. Methods 196, 124–130. doi: 10.1016/j.jneumeth.2010.12.023

PubMed Abstract | Crossref Full Text | Google Scholar

Frey, S., Pandya, D. N., Chakravarty, M. M., Bailey, L., Petrides, M., and Collins, D. L. (2011). An MRI based average macaque monkey stereotaxic atlas and space (MNI monkey space). NeuroImage 55, 1435–1442. doi: 10.1016/j.neuroimage.2011.01.040

Crossref Full Text | Google Scholar

Fukaya, C., Sumi, K., Otaka, T., Obuchi, T., Kano, T., Kobayashi, K., et al. (2010). Nexframe frameless stereotaxy with multitract microrecording: accuracy evaluated by frame-based stereotactic X-ray. Stereotact. Funct. Neurosurg. 88, 163–168. doi: 10.1159/000313868

PubMed Abstract | Crossref Full Text | Google Scholar

Galante, R., Carvalho, E. R., Muniz, J. A., Castro, P. H., Gris, V. N., Amora, D. S., et al. (2019). Comparison between total intravenous anesthesia with propofol and intermittent bolus of tiletamine-zolazepam in capuchin monkey (Sapajus apella). Pesquisa Veterinária Brasileira 39, 271–277. doi: 10.1590/1678-5150-pvb-5847

Crossref Full Text | Google Scholar

Gonzalo Domínguez, M., Hernández, C., Ruisoto, P., Juanes, J. A., Prats, A., and Hernández, T. (2016). Morphological and volumetric assessment of cerebral ventricular system with 3D slicer software. J. Med. Syst. 40, 1–8. doi: 10.1007/s10916-016-0510-9

Crossref Full Text | Google Scholar

Hardman, C. D., and Ashwell, K. W. (2012). Stereotaxic and chemoarchitectural atlas of the brain of the common marmoset (Callithrix jacchus). Boca Raton, FL: CRC Press.

Google Scholar

Hasegawa, M., Kida, I., and Wada, H. (2010). A volumetric analysis of the brain and hippocampus of rats rendered perinatal hypothyroid. Neurosci. Lett. 479, 240–244. doi: 10.1016/j.neulet.2010.05.070

PubMed Abstract | Crossref Full Text | Google Scholar

Hayashi, T., Hou, Y., Glasser, M. F., Autio, J. A., Knoblauch, K., Inoue-Murayama, M., et al. (2021). The nonhuman primate neuroimaging and neuroanatomy project. NeuroImage 229:117726. doi: 10.1016/j.neuroimage.2021.117726

PubMed Abstract | Crossref Full Text | Google Scholar

Hecht, E. E., Reilly, O. T., Benítez, M. E., Phillips, K. A., and Brosnan, S. F. (2021). Sex differences in the brains of capuchin monkeys (Sapajus [Cebus] apella). J. Comp. Neurol. 529, 327–339. doi: 10.1002/cne.24950

PubMed Abstract | Crossref Full Text | Google Scholar

Hikishima, K., Quallo, M., Komaki, Y., Yamada, M., Kawai, K., Momoshima, S., et al. (2011). Population-averaged standard template brain atlas for the common marmoset (Callithrix jacchus). NeuroImage 54, 2741–2749. doi: 10.1016/j.neuroimage.2010.10.061

PubMed Abstract | Crossref Full Text | Google Scholar

Ikonomidou, C., Bittigau, P., Koch, C., Genz, K., Hoerster, F., Felderhoff-Mueser, U., et al. (2001). Neurotransmitters and apoptosis in the developing brain. Biochem. Pharmacol. 62, 401–405. doi: 10.1016/S0006-2952(01)00696-7

Crossref Full Text | Google Scholar

Janson, C. H., and Boinski, S. (1992). Morphological and behavioral adaptations for foraging in generalist primates: the case of the cebines. Am. J. Phys. Anthropol. 88, 483–498. doi: 10.1002/ajpa.1330880405

PubMed Abstract | Crossref Full Text | Google Scholar

Khateeb, K., Griggs, D. J., Sabes, P. N., and Yazdan-Shahmorad, A. (2019). Convection enhanced delivery of optogenetic adeno-associated viral vector to the cortex of rhesus macaque under guidance of online MRI images. JoVE (J. Vis. Exp.) 147:e59232. doi: 10.3791/59232

Crossref Full Text | Google Scholar

Lankau, E. W., Turner, P. V., Mullan, R. J., and Galland, G. G. (2014). Use of nonhuman primates in research in North America. J. Am. Assoc. Lab. Anim. Sci. 53, 278–282. doi: 10.30802/AALAS-JAALAS-13-000107

PubMed Abstract | Crossref Full Text | Google Scholar

Lynch Alfaro, J. W., Boubli, J. P., Olson, L. E., Di Fiore, A., Wilson, B., Gutiérrez-Espeleta, G. A., et al. (2012). Explosive Pleistocene range expansion leads to widespread Amazonian sympatry between robust and gracile capuchin monkeys. J. Biogeogr. 39, 272–288. doi: 10.1111/j.1365-2699.2011.02609.x

Crossref Full Text | Google Scholar

Ma, Y., Smith, D., Hof, P. R., Foerster, B., Hamilton, S., Blackband, S. J., et al. (2008). In vivo 3D digital atlas database of the adult C57BL/6J mouse brain by magnetic resonance microscopy. Front. Neuroanat. 2:175. doi: 10.3389/neuro.05.001.2008

Crossref Full Text | Google Scholar

Manocha, S. L., Shantha, T. R., and Bourne, G. H. (1968). A stereotaxic atlas of the brain of the Cebus monkey (Cebus apella). Oxford University, Great Britain Clarendon Press.

Google Scholar

Marroig, G. (2007). When size makes a difference: allometry, life-history and morphological evolution of capuchins (Cebus) and squirrels (Saimiri) monkeys (Cebinae, Platyrrhini). BMC Evol. Biol. 7, 20–26. doi: 10.1186/1471-2148-7-20

Crossref Full Text | Google Scholar

Miss, F. M., Meunier, H., and Burkart, J. M. (2022). Primate origins of corepresentation and cooperative flexibility: a comparative study with common marmosets (Callithrix jacchus), brown capuchins (Sapajus apella), and Tonkean macaques (Macaca tonkeana). J. Comp. Psychol. 136, 199–212. doi: 10.1037/com0000315

PubMed Abstract | Crossref Full Text | Google Scholar

Muniz, J. A. C. P., Leal, L. C. P., Bahia, C. P., and Krejcová, L. V. (2021). Sapajus apella as a model for the development of novel therapeutic approaches for Parkinson's disease. Revista Pan-Amazônica de Saúde 12:11. doi: 10.5123/S2176-6223202100514

Crossref Full Text | Google Scholar

Ose, T., Autio, J. A., Ohno, M., Frey, S., Uematsu, A., Kawasaki, A., et al. (2022). Anatomical variability, multi-modal coordinate systems, and precision targeting in the marmoset brain. Neuro Image 250:118965. doi: 10.1016/j.neuroimage.2022.118965

Crossref Full Text | Google Scholar

Palazzi, X., Bordier, N., Palazzi, X., and Bordier, N. (2008). The marmoset brain in stereotaxic coordinates New York, NY: Springer.

Google Scholar

Patil, A. A. (2010). A modified stereotactic frame as an instrument holder for frameless stereotaxis. Surg. Neurol. Int. 1:62. doi: 10.4103/2152-7806.70957

PubMed Abstract | Crossref Full Text | Google Scholar

Paxinos, G., and Franklin, K. B. (2019). Paxinos and Franklin's the mouse brain in stereotaxic coordinates. Cambridge: Cambridge Academic press.

Google Scholar

Paxinos, G., Huang, X.-F., and Toga, A. W. (2000). The rhesus monkey brain in stereotaxic coordinates. San Diego, USA: Academic Press.

Google Scholar

Pereira-Pedro, A. S., Rilling, J. K., Chen, X., Preuss, T. M., and Bruner, E. (2017). Midsagittal brain variation among non-human primates: insights into evolutionary expansion of the human precuneus. Brain Behav. Evol. 90, 255–263. doi: 10.1159/000481085

PubMed Abstract | Crossref Full Text | Google Scholar

Phillips, K. A., and Sherwood, C. C. (2008). Cortical development in brown capuchin monkeys: a structural MRI study. NeuroImage 43, 657–664. doi: 10.1016/j.neuroimage.2008.08.031

PubMed Abstract | Crossref Full Text | Google Scholar

Phillips, K. A., Sherwood, C. C., and Lilak, A. L. (2007). Corpus callosum morphology in capuchin monkeys is influenced by sex and handedness. PLoS One 2:e792. doi: 10.1371/journal.pone.0000792

PubMed Abstract | Crossref Full Text | Google Scholar

Rosa, M. G., Soares, J. G., Fiorani, M., and Gattass, R. (1993). Cortical afferents of visual area MT in the Cebus monkey: possible homologies between new and Old World monkeys. Vis. Neurosci. 10, 827–855. doi: 10.1017/S0952523800006064

PubMed Abstract | Crossref Full Text | Google Scholar

Ryland, A. B., Mittermeier, R. A., and Rodriguez-Luna, E. (1997). Conservation of neotropical primates: threatened species and an analysis of primate diversity by country and region. Folia Primatol. 68, 134–160. doi: 10.1159/000157243

Crossref Full Text | Google Scholar

Sansalone, G., Allen, K., Ledogar, J., Ledogar, S., Mitchell, D., Profico, A., et al. (2020). Variation in the strength of allometry drives rates of evolution in primate brain shape. Proc. R. Soc. B 287:20200807. doi: 10.1098/rspb.2020.0807

PubMed Abstract | Crossref Full Text | Google Scholar

Saunders, R., Aigner, T., and Frank, J. (1990). Magnetic resonance imaging of the rhesus monkey brain: use for stereotactic neurosurgery. Exp. Brain Res. 81, 443–446. doi: 10.1007/BF00228139

PubMed Abstract | Crossref Full Text | Google Scholar

Seidlitz, J., Sponheim, C., Glen, D., Frank, Q. Y., Saleem, K. S., Leopold, D. A., et al. (2018). A population MRI brain template and analysis tools for the macaque. NeuroImage 170, 121–131. doi: 10.1016/j.neuroimage.2017.04.063

PubMed Abstract | Crossref Full Text | Google Scholar

Sousa, A. P., Piñon, M. C. G., Gattass, R., and Rosa, M. G. (1991). Topographic organization of cortical input to striate cortex in the Cebus monkey: a fluorescent tracer study. J. Comp. Neurol. 308, 665–682. doi: 10.1002/cne.903080411

PubMed Abstract | Crossref Full Text | Google Scholar

Subramanian, T., Deogaonkar, M., Brummer, M., and Bakay, R. (2005). MRI guidance improves accuracy of stereotaxic targeting for cell transplantation in parkinsonian monkeys. Exp. Neurol. 193, 172–180. doi: 10.1016/j.expneurol.2004.11.032

PubMed Abstract | Crossref Full Text | Google Scholar

Sutherland, G. R., Wolfsberger, S., Lama, S., and Zarei-Nia, K. (2013). The evolution of neuro arm. Neurosurgery 72, A27–A32. doi: 10.1227/NEU.0b013e318270da19

PubMed Abstract | Crossref Full Text | Google Scholar

Vannucci, R. C., Barron, T. F., and Vannucci, S. J. (2017). Development of the corpus callosum: an MRI study. Dev. Neurosci. 39, 97–106. doi: 10.1159/000453031

Crossref Full Text | Google Scholar

Visalberghi, E. (1993). “Capuchin monkeys: a window into tool use in apes and humans” in Tools, language and cognition in human evolution, 138–150.

Google Scholar

Watson, C. M., Sherwood, C. C., and Phillips, K. A. (2022). Myelin characteristics of the corpus callosum in capuchin monkeys (Sapajus [Cebus] apella) across the lifespan. Sci. Rep. 12:8786. doi: 10.1038/s41598-022-12893-z

PubMed Abstract | Crossref Full Text | Google Scholar

Yazdan-Shahmorad, A., Silversmith, D. B., Kharazia, V., and Sabes, P. N. (2018). Targeted cortical reorganization using optogenetics in non-human primates. elife 7:e31034. doi: 10.7554/eLife.31034

PubMed Abstract | Crossref Full Text | Google Scholar

Yin, D., Valles, F. E., Fiandaca, M. S., Forsayeth, J., Larson, P., Starr, P., et al. (2009). Striatal volume differences between non-human and human primates. J. Neurosci. Methods 176, 200–205. doi: 10.1016/j.jneumeth.2008.08.027

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: non-human primates, magnetic resonance imaging, brain, neurosurgery, stereotaxic coordinates

Citation: Pedrosa LRR, Leal LCP, Muniz JAPC, Bastos CdO, Gomes BD and Krejcová LV (2024) From imaging to precision: low cost and accurate determination of stereotactic coordinates for brain surgery Sapajus apella using MRI. Front. Neurosci. 18:1324669. doi: 10.3389/fnins.2024.1324669

Received: 19 October 2023; Accepted: 17 January 2024;
Published: 01 February 2024.

Edited by:

Shengxiang Liang, Fujian University of Traditional Chinese Medicine, China

Reviewed by:

Ashik Banstola, University of Otago, New Zealand
David Schaeffer, University of Pittsburgh, United States

Copyright © 2024 Pedrosa, Leal, Muniz, Bastos, Gomes and Krejcová. 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: Lane V. Krejcová, lane@ufpa.br

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.