Abstract
Disorders of the medial temporal lobe (MTL) adversely affect visual working memory (vWM) performance, including feature binding. It is unclear whether these impairments generalize across visual dimensions or are specifically spatial. To address this issue, we compared performance in two tasks of 13 epilepsy patients, who had undergone a temporal lobectomy, and 15 healthy controls. In the vWM task, participants recalled the color of one of two polygons, previously displayed side by side. At recall, a location or shape probe identified the target. In the perceptual task, participants estimated the centroid of three visible disks. Patients recalled the target color less accurately than healthy controls because they frequently swapped the nontarget with the target color. Moreover, healthy controls and right temporal lobectomy patients made more swap errors following shape than space probes. Left temporal lobectomy patients, showed the opposite pattern of errors instead. Patients and controls performed similarly in the perceptual task. We conclude that left MTL damage impairs spatial binding in vWM, and that this impairment does not reflect a perceptual or attentional deficit.
Significance Statement
This study examined color recall in temporal lobectomy patients and healthy controls, to determine whether patients show differential impairments binding color and shape versus color and location of memorized objects. Left temporal lobectomy patients were less accurate recalling color, especially when the target object was identified by the location, rather than the shape it had in the initial display. We found no group difference in a task, which required estimating the centroid of three circles, indicating that the memory impairment was not accounted by perceptual or attentional difficulties. Our findings indicate that lateralized medial temporal circuits are crucial for binding visual features to the location where they had appeared, thus ensuring the primacy of space in organizing declarative memories.
Introduction
The role of the medial temporal lobe (MTL) in episodic memory is well established (Squire, 2009). Despite initial reports of preserved immediate memory span in temporal lobectomy patients (Drachman and Arbit, 1966), later studies found that MTL lesions also engender substantial working memory (WM) deficits (Olton et al., 1979; Holdstock et al., 1995; Hannula et al., 2006). Which WM processes are specifically supported by the MTL is, nevertheless, a matter of ongoing investigations.
An early, seminal model suggested that visual WM (vWM) contains few discrete “slots,” each used to store one and only one object with high fidelity (Luck and Vogel, 1997; Zhang and Luck, 2008). Despite its simplicity, the slot model makes nontrivial predictions. First, the complexity of memorized objects should not affect recall accuracy. Second, recalling feature conjunctions should carry no additional cost over remembering features, since features are stored ipso facto as parts of an object into vWM. Wheeler and Treisman (2002) found, instead, that simple objects were recalled more accurately than complex ones, and that recall accuracy was equalized for displays containing the same number of color features rather than the same number of objects. They concluded that memory limitations reflect feature rather than object-based storage mechanisms. Moreover, observers were worse at detecting changes of feature conjunctions than features, indicating that conjunctions are stored or recalled less efficiently than features. Later studies confirmed that changes in feature conjunctions are poorly detected (Allen et al., 2014), leading to the suggestion that dimensionally specific registers store features, while an “episodic buffer” is dedicated to holding bound object representations in vWM (Baddeley et al., 2011). The need for binding processes follows logically from the alternative model of vWM which proposes that visual features are stored in dimensionally specific, limited resolution stores (Alvarez and Cavanagh, 2004; Smyrnis et al., 2005; Bays et al., 2009). Clearly, if different feature dimensions are stored separately, then a binding process is required to ensure that features belonging to the same object, but different feature dimensions, are identified as such (Wheeler and Treisman, 2002; Smyrnis et al., 2005).
While the idea that conjunctive binding is required to preserve object identity is not unanimously shared (Luck and Vogel, 2013), there are several proposals regarding its nature. Treisman and Zhang (2006) concluded that binding is automatic, established initially by the features’ shared location, but then becomes location independent. Schneegans and Bays (2017) proposed instead that location information is always required, because features from different visual dimensions are stored in separate retinotopic maps.
Investigators examining the neurologic underpinnings of declarative memory largely embrace the idea that space plays a crucial role in indexing declarative memories. Animal and patient studies (Chalfonte et al., 1996; Brown and Aggleton, 2001; Eacott and Gaffan, 2005; Piekema et al., 2006; Ranganath, 2010; Libby et al., 2014) documented a functional parcellation of the MTL with separate structures representing, respectively, the environmental layout, the objects within it, as well as binding the latter to the former. According to this view, space is crucial for recalling events, but not for binding object features. Olson et al. (2006) for example, reported that patients with postanoxic or postencephalitic MTL pathology had impaired object-location binding in a WM task. This impairment was unaccounted by either diminished recognition or spatial memory. In animals, MTL lesions are also followed by dissociated impairments in object recognition and recall of object-location conjunctions, suggesting that feature and spatial binding depend on distinct MTL processes (Meunier et al., 1993; Murray and Mishkin, 1998; Malkova and Mishkin, 2003). Studies in temporal lobe epilepsy (TLE) patients reported deficits in spatial recall, spatial binding and visual recognition (Abrahams et al., 1997; Bohbot et al., 1998; Stepankova et al., 2004); however, it remains unclear whether these impairments should be attributed to diminished spatial precision (Kolarik et al., 2016) or a binding deficit (Zokaei et al., 2019) and whether binding impairments are dimensionally general (Hannula et al., 2006; Pertzov et al., 2013) or specific.
To examine these issues, we tested TLE patients who had undergone temporal lobectomies and healthy controls with two tasks used in a previous investigation of a stroke patient (Dundon et al., 2018). The first requires recalling the color of one of two polygons identified by either a location or a shape probe, thus directly pitting spatial versus nonspatial binding. The second probes participants’ ability to estimate the average position of three visible disks. In this task healthy participants show a pseudo-neglect pattern of leftward static errors (Baud-Bovy and Soechting, 2001), which suggests that centroid estimation is sensitive to attentional biases (Dundon et al., 2018). The centroid task was therefore used to highlight the presence of unilateral neglect, which can follow lesions of the nondominant parahippocampal cortex (Mort et al., 2003) as well as the integrity of spatial perception and attention (Drew et al., 2010).
Materials and Methods
The aim of the present study was to compare nonspatial and spatial binding performance in TLE patients with medically refractory epilepsy who had undergone temporal lobectomy and healthy controls. Recruitment and testing took place at the Neuropsychology section of the Department of Neurosciences of King Faisal Specialist Hospital and Research Centre in Riyadh, Saudi Arabia. The experimental protocol was approved by the local Institutional Review Board. Participants gave written informed consent before engaging in any experimental procedure.
Study participants
Over a month period, an opportunity sample of patients attending the Neurology and Neurosurgical Clinics were invited to participate in the study. All had been diagnosed with TLE on the basis of clinical presentation and instrumental diagnostic procedures, inclusive of ambulatory EEG and neuroimaging, and after failing medical therapy, had undergone surgical treatment. All patients had normal or corrected to normal visual acuity. Those with an estimated full-scale IQ of <75, as assessed with an Arabic version of the Wechsler Abbreviated Scale of Intelligence, Second Edition (WASI-II; Al-Joudi et al., 2019), were excluded, as well as those with a history of traumatic brain injury or major psychiatric disorders. Patients who suffered a seizure in the preceding 48 h had the testing session postponed. Thirteen patients took part in this study. The study’s neurosurgeon identified the anatomic structures involved by the resection on the basis of the surgical record and postsurgical MRI scans.
Fifteen healthy participants were concurrently recruited from the local community as controls. Potential participants were excluded if they had a history of a major neurologic or psychiatric disorder or an uncorrected visual impairment or an estimated IQ <75.
Testing protocol
Testing took place in a quiet, dimly lit room. Participants sat comfortably at a distance of ∼70 cm from a MacBook Pro set at a resolution of 1680 × 1050 pixels. Custom-coded MATLAB (MathWorks) scripts used a set of freely available routines (Brainard, 1997) to control the timing of the displays. Two tasks were conducted using the computerized set-up.
Cued color recall
In each trial, an equilateral triangle and a square, whose side lengths are 1.29° and 0.92°, respectively, appear side-to-side in the lower half of the screen, as shown in Figure 1A. The shapes are centered at an eccentricity of 2.27° along the main screen diagonals and remain visible for 2.0 s. The two shapes are of different colors, either red, blue, or green. The sample display is followed by a 0.2-s-long pattern mask and a 1.3-s blank screen. The recall screen contains three colored rectangles, 0.53° wide and 1.59° high, whose lower edges are aligned 1.39° above the screen center and spaced horizontally 4.81° apart. A bright cross (location probe) or the outline of one of the polygons (shape probe) identifies the target. The location probe, which also includes a dark cross, appear at the locations occupied by the two shapes. The shape cue appears 3.0° below the screen center. Participants report the target color by placing the cursor over the corresponding colored rectangle and clicking the mouse button. The mouse click prompts the beginning of a new trial after a 1.0-s delay during which the screen is blank. Participants practiced the task over 16 trials and then completed two blocks of 96 trials each, including both shape and location cued recalls. Trial order was randomized, minimizing participants’ ability to predict whether a shape or location probe would follow the sample display. To ensure that patients had not forgotten the task instructions, at the end of each block they were asked to describe what they had done.
Centroid estimation
This task assesses the accuracy and precision of estimates of the average location of three visible white disks and is illustrated in Figure 1B. Each disk’s diameter is 0.27°. Participants place a crosshair-shaped cursor at the estimated centroid location and click the mouse. Following a 1.0-s interval, a novel set of disks appears and the procedure is repeated. Disks can occupy any three of seven canonical locations, including the screen center and the vertices of a virtual concentric hexagon, with a side length of 3.67°. All permutations of three canonical locations, less any resulting in a collinear configuration, are used as test arrays. Each possible permutation appears twice, for a total of 64 trials. Pseudorandom, zero mean, independent circular Gaussian distributions, with a SD of 0.6°, are sampled to jitter each disk’s position. Before testing, instructions were read to the participants. The centroid was defined as the point where the triangle, whose vertices coincided with the disks’ locations, would balance in the horizontal plane (Baud-Bovy and Soechting, 2001). Participants completed 25 practice trials, without feed-back, followed by two blocks of 64 trials each.
Neuropsychological tests
Three neuropsychological instruments were used to assess participants: (1) Hopkins Verbal Learning Test–Revised (HVLT-R), (2) Brief Visuospatial Memory Test–Revised (BVMT-R), and (3) Color Trails Test (CTT). The Arabic version of these tests were recently validated (Al-Joudi et al., 2019).
Data analysis
In the recall task, participants could either report (1) the color of the target, correct response; (2) the color of the nontarget item, that is make a swap error; or (3) the color absent from the sample, that is make a generic error. We approached the group level hypothesis testing as a metanalysis of prevalence data, treating each participant as a separate study. All inferential analysis presented in this study is Bayesian.
Group differences in recall error rates
Group level effects were analyzed with mixed Bayesian ANOVAs using the JASP software (JASP Team, 2021; https://jasp-stats.org/). The between group variable coded whether the participant was (1) a healthy control, (2) a patient following left, or (3) right temporal lobectomy, respectively. The within group variables were probe dimension, i.e., whether a shape or space probe was used to cue recall, and block order, i.e., first or second block. The error rates were normalized with a Friedman–Tukey double arcsine transformation (Barendregt et al., 2013). The transformation stabilizes error rates variances and allows the use of parametric methods to compute group statistics:
Analysis of centroid task
The analysis of the centroid task was conducted by fitting the following regression model to each participants’ reports:
Results
Participants’ demographic, clinical, and neuropsychometric characteristics
Table 1 reports the demographic and clinical characteristics of the left and right temporal lobectomy patients and healthy controls. The groups were matched on age, gender, and educational level. Both patient groups showed a lower full-scale IQ that healthy controls; however, all of the variables were more likely to reflect a null effect than a group difference. Table 2 details the sex, education level, and neuropsychometric performance of the 13 patients. Both raw scores as well the values normalized on the basis of a reference sample of healthy controls, whose first language is Arabic (Al-Joudi et al., 2019), are shown. We examined whether patients showed a material specific pattern of lateralized deficits (Saling, 2009) by comparing performance of the left and right temporal lobectomy patients on the HVMT and the BVMT with Bayesian independent samples t tests (JASP Team, 2021; https://jasp-stats.org/). There was moderate strength evidence for left temporal lobectomy patients having worse delayed verbal recall on the HVLT than controls (BF10 = 7.99), while there was anecdotal evidence for no group difference in the delayed visuo-spatial recall (BF10 = 0.46). In Table 3, the MTL structures affected by the surgical excisions are listed, patient by patient, while Extended Data Table 3-1 shows representative postsurgical axial, sagittal, and coronal multimodal MRI slices for each patient.
Extended Data Table 3-1
Postsurgical MRI scans for 13 patients. The images were obtained with T1 weighted, T2 weighted, fluid-attenuated inversion recovery, and gradient recalled echo sequences. ERC, entorhinal cortex; PHC, parahippocampal cortex; PRC, perirhinal cortex; Hipp, hippocampus; AMG, amygdala; TP, temporal pole; antSTG, anterior superior temporal gyrus; antMTG, anterior middle temporal gyrus. Download Table 3-1, DOCX file.
Cued color recall performance
In the cued recall task, participants completed two blocks of 96 trials each. Mixed effect Bayesian ANOVAs were used to examine the influence of three factors on recall: (1) group, namely, whether participants were controls, patients following left and right temporal lobectomy, respectively; (2) block; and (3) probe dimension. Generic and swap errors were analyzed separately.
Extended Data Table 4-1 reports the result of the ANOVA model comparison for swap errors. The model with the highest posterior probability included group, probe dimension, and block, as well as the interaction of group by probe dimension. Table 4 summarizes the evidence for each predictor. There was moderate evidence for an effect of probe dimension and block. There was very strong evidence for an effect of group and strong evidence for an interaction of group by probe dimension.
Extended Data Table 4-1
Best model comparisons for swap errors. The table presents each of the model comparisons from the Bayesian ANOVA. The within factors are block (B) and probe dimension (D). The between factor is group G. P(M) is the a priori model probability, P(M|d) is the posterior model probability. BFM is the Bayes factor of the model, BF10 is the Bayes factor of the model relative to the best one. The best model contained all three factors and the interaction of group by probe dimension. Download Table 4-1, DOC file.
Extended Data Table 5-1 summarizes the results of the model comparisons for generic errors. The best model was the one which included all three factors, but none of the interactions. Table 5 summarizes the analysis of the effects. There was anecdotal evidence for the block and probe dimension affecting the proportion of generic errors. On the other hand, there was anecdotal evidence for a null effect of group.
Extended Data Table 5-1
Best model comparisons for generic errors. The table presents the model comparisons obtained from a Bayesian ANOVA. The within factors are block (B) and probe dimension (D). The between factor is group G. P(M) is the a priori model probability, P(M|d) is the posterior model probability. BFM is the Bayes factor of the model, BF10 is the Bayes factor of the model relative to the best one. The best model included the three main factors, namely, block (B), probe dimension (D), and group (G). Download Table 5-1, DOC file.
Figure 2 shows the mean proportion of swap and generic errors following space and shape probe, respectively. Controls made fewest swap and generic errors. Controls and patients with right temporal lobectomies made more swap errors following shape than space probes. Patients with left temporal lobectomies made most swap errors and three of them made more swap errors following space than shape probes. Participants made more generic errors following space than shape probe. Participants also made more swap and generic errors in the first than second block (data not shown).
Centroid estimation performance
We examined how group and screen coordinates affected three indices of performance in the centroid estimation task: (1) static offsets, (2) accuracy, and (3) precision (see Materials and Methods). For static offsets the model with the highest posterior probability included the effect of screen coordinates only. In fact, there was strong evidence (BFinc = 32.13) that the horizontal and vertical offsets differed. While there was no appreciable horizontal bias, m = 0.0°, 95%CI = [−0.04°, 0.05°], participants showed an upward bias, m = 0.15°, 95%CI = [0.07°, 0.23°]. There was moderate evidence in favor of the null and against both an effect of group (BFinc = 0.22), and its interaction with coordinate (BFinc = 0.28). For accuracy, the model with highest posterior probability included group. However, there was anecdotal evidence for a null effect of group (BFinc = 0.76), with moderate evidence in favor of a null effect of coordinate (BFinc = 0.26) and its interaction with group (BFinc = 0.28). For precision, namely, the variance of the variable errors, the null model had the highest posterior probability. There was anecdotal evidence for a null effect of screen coordinate (BFinc = 0.4), anecdotal evidence for a null effect of group (BFinc = 0.53), and moderate evidence for a null interaction of group and dimension (BFinc = 0.15).
Discussion
In this study, we compared performance of healthy controls and temporal lobectomy patients in two tasks, one probing conjunctive binding in vWM, the other perceptual, spatial averaging of disks patterns. The vWM task required the recall of a target’s color, where the target was identified by either a location or a shape probe. The task was thus designed to determine dimensional specificity of WM conjunctive binding. Controls were more accurate than patients overall. Moreover, they made fewer swap errors following space than shape probes, while left temporal lobectomy patients made more swap errors following space than shape probes. There was no evidence of group differences in static biases, accuracy and precision of perceptual centroid estimates. The implication of these findings for the organization of binding and spatial processes in vWM is discussed in the next paragraph, following a brief overview of prior evidence.
MTL lesions specifically disrupt spatial binding in vWM
Previous studies addressed whether MTL pathology is associated with impairments in vWM binding. The ability to recall shape-color, shape-location or item-item conjunctions has been found to be diminished in patients with anoxic/ischemic or infectious pathology involving the MTL as well as neurodegenerative disorders, suggesting an impairment in conjunctive and relational binding (Hannula et al., 2006; Olson et al., 2006; Parra et al., 2009). van Geldorp et al. (2014) compared patients, who had undergone anterior temporal lobectomies for medically refractory epilepsy, and healthy controls’ performance in four match-to-sample tasks. The tasks were difficult and required participants to remember three separate frames presented sequentially. Each frame contained the picture of a face and a building, which differed in location and color. A cue, presented before the sample, indicated whether participants should only remember the identity of the items or also their location (spatial binding condition), color (color binding condition), or the item they had been presented with (relational binding condition). Overall, recall was less accurate in patients compared with controls, and particularly so in the relational binding condition. Recall performance in the spatial and color binding conditions were equally affected, suggesting that spatial and nonspatial WM binding were not differentially compromised in these patients. Zokaei et al. (2019) found that patients who had undergone a temporal lobectomy made more swap errors when recalling the location of fractal patterns, compared with controls. Since neither fractal recognition nor memory for locations were found to be appreciably impaired in these patients, it was inferred that they suffered a primary binding deficit. Using a similar paradigm, Pertzov et al. (2013) documented both a spatial as well as a nonspatial binding deficit in individuals recovering from autoimmune encephalitis, suggesting that binding impairments because of MTL dysfunction are not dimensionally specific. Braun et al. (2011) concluded instead that TLE patients, who had undergone a right temporal lobectomy, were only impaired when the vWM task required spatial binding but performed similarly to healthy controls when it required binding of nonspatial features.
Our own findings contribute new, crucial evidence for understanding the role of MTL in vWM binding by confirming the association between MTL pathology and spatial binding impairments, unaccounted for by impairments of either spatial vision or feature memory. We found an increase in the proportion of swap errors in the TLE group. Crucially, while healthy participants made significantly more swap errors following the shape probe, left temporal lobectomy patients made more swap errors following the space probe. In healthy participants the results are therefore in keeping with the hypothesis that binding of nonspatial features is mediated by the features’ shared location (Treisman and Zhang, 2006; Schneegans and Bays, 2017). In fact, the likelihood of swap errors should be greater following shape compared with space probes, because in the latter case both target shape and color need to be bound to the target location before they can be bound to each other (Schneegans and Bays, 2017). On the other hand, the observation that some patients made significantly more swap errors following space probes than controls may indicate that patients gained the ability to bind nonspatial features directly, without the mediation of a shared location, allowing them to achieve higher accuracies following shape than space probes. Whether this inference is warranted remains to be established. Regardless, the group level pattern of dimensionally specific binding impairments observed in left TLE patients replicates a previous observation in a stroke patient with bilateral MTL damage, found to be impaired only in vWM tasks requiring spatial binding, but not those requiring nonspatial binding (Dundon et al., 2018). These observations thus allow us to draw the following conclusion: MTL pathology can be associated with WM binding impairments that are spatially specific and reverse the spatial advantage characteristic of healthy controls. If processes underlying spatial binding in vWM are independent from processes devoted to binding of nonspatial features, then the role of space in organizing vWM may extend beyond providing a common index for the conjunctive binding of visual features.
Hemispheric lateralization and binding
Bohbot et al. (1998) found that patients who had undergone thermocoagulation of structures within the right, but not the left, MTL were more impaired in a number of spatial WM tasks than those who had not undergone surgery, suggesting that right MTL structures may play an overarching role in spatial memory. Braun et al. (2011) compared patients with right temporal lobectomy and healthy controls’ performance on a number of single feature and feature conjunction recall tasks and concluded that these patients are specifically impaired in spatial binding. However, the tasks employed memory samples of different complexity to probe recall of features and conjunctions, respectively, thus introducing a load confound in the comparison. Our own results are in keeping with the idea that left rather than right MTL structures are specifically involved in spatial binding, since patients who had undergone left temporal lobectomies showed greater spatial binding impairments than patients with right temporal lobectomies. While our results need to be interpreted cautiously, given the small sample size, they agree with the conclusion drawn by Kessels et al. (2004), who found that patients who had undergone left, but not right temporal lobectomies were impaired in spatial binding, confirming lateralization effects previously observed by the same group in a sample of patients with cerebrovascular pathology (Kessels et al., 2002). Spiers et al. (2001) found that left temporal lesion specifically affect object location binding while right temporal lesions affect spatial memory more generally. However, the existing literature remains inconclusive to the relation between laterality and WM spatial binding. A distinct view of lateralization of binding impairments is that the latter reflect attention deficits which follow cortical strokes, especially involving the nondominant right hemisphere. For example. Cohen-Dallal et al. (2021) reported that patients with unilateral attentional neglect show delay-dependent decrements in spatial binding performance. In light of the fact that parahippocampal damage is associated with unilateral neglect (Mort et al., 2003), Cohen-Dallal et al. (2021)’s finding raises the possibility that lateralized attentional deficits may also contribute to the binding impairments observed in our study. However, performance in a centroid estimation task did not show group differences in lateralized static biases, suggesting that lesion in our sample was not associated with unilateral neglect. Moreover, patients showed neither diminished accuracy nor lower precision in the centroid task compared with controls, indicating that spatial perception and attention were not compromised (Drew et al., 2010).
A further concern is that uncontrolled verbal strategies may have confounded the interaction of lesion laterality and probe dimension. We found in pilot studies that healthy participants maintain a spatial advantage in WM binding under condition of articulatory suppression and therefore concluded that the spatial advantage in WM binding does not depend on verbal strategies (d’Avossa G, Dundon NM, Katshu MZUH, unpublished data). However, we cannot rule out the possibility that left temporal lobectomy patients used a verbal strategy and thus managed to selectively improve binding of shape and color.
It is important to note, with regard to the issue of localization and the nature of cognitive impairments encountered in TLE patients, that group level results belie substantial interindividual differences (see Fig. 2). Previous electrophysiological recordings from the left MTL indicated that trial by trial changes in the amplitude of low-frequency oscillations, obtained during encoding episodes, predict subsequent recall accuracy of object/place conjunctions in TLE patients (Miller et al., 2018), in keeping with our own conclusion that spatial binding is dependent on left lateralized processes. Interestingly, this was not the case in all suggesting that spatial binding processes are not universally left lateralized in these patients. Unfortunately, the study did not report whether the lateralization of spatial binding processes was affected by the laterality of the seizure focus, precluding firmer conclusion regarding the relation between the two. Other studies have, however, demonstrated anomalous lateralization in high proportion of TLE patients, as memory processes often shift to the contralesional hemisphere both preoperatively (Bellgowan et al., 1998; Golby et al., 2002; Janszky et al., 2004) and postoperatively (Sidhu et al., 2016). These findings may provide one possible interpretation of the interindividual differences highlighted above, namely, that neural plasticity in some TLE patients modifies the lateralization of memory processes usually encountered in healthy controls. An alternative explanation, initially born out of observations in nonhuman primates with focal lesions (Browning et al., 2010; Croxson et al., 2012), is that effects of MTL functional and structural damage in postsurgical TLE patients may be attenuated by nonlateralized recruitment of neocortical areas (Sidhu et al., 2013). To determine whether either or both of these hypotheses can account for interindividual differences in spatial binding impairments and to what extent other factors, like age of seizure onset, severity and frequency of seizures, neuropsychiatric complications, and antiepileptic medications, known to affect neural plasticity and degree of cognitive impairment (Bell et al., 2011), also contribute to hemispheric lateralization in TLE patients will require new experimental evidence.
Despite the potential confound listed above, the present study indicates that following left temporal lobectomy, vWM binding is diminished in a dimensionally specific manner in the absence of appreciable perceptual, attentional, or visual-spatial memory deficits.
Footnotes
The authors declare no competing financial interests.
M.F.A. was supported by a doctoral scholarship from the Organ Transplant Centre at King Faisal Specialist Hospital and Research Centre in Riyadh, Saudi Arabia.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.