Abstract
Neurotransmission imaging studies require normative data for the statistical assessment of neurophysiologic dysfunctions. 2′-Methoxyphenyl-(N-2′-pyridinyl)-p-18F-fluoro-benzamidoethylpiperazine (18F-MPPF) is a specific serotonin 5-HT1A antagonist PET tracer recently characterized, modeled, and used for clinical research to explore abnormalities in the serotoninergic system. Our study reports, to our knowledge, the first large normative imaging database of 18F-MPPF binding potential (BP) over aging, for both males and females. Methods: Fifty-three healthy volunteers (27 females, 26 males; age, 20–70 y) were selected to undergo structural MRI and single-injection 18F-MPPF multiframe dynamic PET. 18F-MPPF BP values were computed using a nonlinear modeling method with tissue reference. The statistical assessment of the effect of age and sex was performed both at the anatomic structure level, using regions of interest drawn manually on individual MR images, and at the voxel level, using normalized BP parametric images in different statistical parametric mapping designs. Results: A negative linear correlation between age and 18F-MPPF binding (3.6% decrease by decade) was found in females but not in males and involved most of the limbic and paralimbic regions; on the other hand, males in their 30s showed decreased binding in most cerebral regions. Conclusion: A comparison of males and females revealed higher BP values independent of age in females in the right hemisphere and a different evolution of BP over aging. These results confirm the necessity of a database for further statistical analysis in individuals or groups with pathology.
In the mammalian nervous system, serotonin (5-HT) mediates a large variety of physiologic responses (development, pain, sleep, mood, eating, memory, attention), behaviors (stress, aggression, panic, sexual behavior), or neuropsychiatric problems (depression, sleep disturbance, eating disorders, anxiety, suicide, schizophrenia, obsessive compulsive disorder) through one of the widest range of receptors known for any neurotransmitter (1). Most of these receptors have been characterized; yet, to date, the 5-HT1A subtype remains the best known in terms of structure, pharmacology, and action.
The recent development of selective PET ligands for 5-HT receptors permits the in vivo exploration of the serotoninergic system in humans (2). In particular, two specific radioligands for 5-HT1A receptors have been used with PET. The first to be developed, 11C-WAY100635, is a high-affinity selective antagonist of 5-HT1A receptors. Binding site extension was delineated (3) and quantified (4). These studies confirmed that 5-HT1A receptors are mainly located in the limbic areas and in the raphe nuclei but are also present in the paralimbic and neocortical regions.
The fluoro analog 2′-methoxyphenyl-(N-2′-pyridinyl)-p-18F-fluoro-benzamidoethylpiperazine (18F-MPPF) was synthesized shortly after (5) and, when compared with 11C-WAY100635, offers an affinity more similar to that of 5-HT for 5-HT1A receptors (inhibition constant = 3.3 nmol/L) (6) and a longer half-life (107 min). In vivo and ex vivo distribution studies in rat (7), cat (5,8), and humans (9) confirmed the selective binding of the molecule in cerebral regions rich in 5-HT1A receptors, such as the dorsal raphe nuclei, the limbic and paralimbic regions, and, with less intensity, the neocortex. Quantification studies using a 3-compartmental simplified reference tissue model (SRTM) confirmed that binding potential (BP) values were correlated with the maximum number of available binding sites (Bmax) and, therefore, could be considered as a reliable index of local receptor concentration (10).
Although these studies emphasized the potential use of both ligands for quantitative studies of the serotoninergic 5-HT1A system in vivo, small groups of subjects were examined, which lead to a rather high variability of BP values. Certainly, significant variations in the quantification of 5-HT1A receptor binding have been described previously using 11C-WAY100635 (11–13). These variations were observed especially in small structures, such as the raphe nuclei or the amygdala, which are of particular interest for the study of neurologic diseases.
The aim of our study was to obtain data for 18F-MPPF PET in a large group of healthy subjects (males and females) to test for the possible influence of age and sex on the variability of 18F-MPPF PET BP. Moreover, our goal was also to provide a large control database appropriate for statistical parametric mapping (SPM) to 18F-MPPF PET studies in 5-HT–related neuropsychiatric disorders (such as depression, eating disorders, Alzheimer, migraine or epilepsy) in which patients of wide age range or of specific sex are recruited.
MATERIALS AND METHODS
Subjects
Fifty-three healthy subjects (26 males [mean age ± SD, 41.7 ± 13.5 y; range, 20–68 y] and 27 females [mean age ± SD, 41.7 ± 14.1 y; range, 19–70 y]) were selected to participate in the study. Male and female control subjects were recruited to obtain a homogeneous age distribution between 19 and 70 y, with approximately 6 male and 6 female subjects for each decade of age. Subjects gave their written consent to participate in the study, which was approved by the local ethical committee (Centre Léon Bérard, Lyon, France) in accordance with the Helsinki Declaration and French rules protecting persons. Subjects were interviewed by a medical doctor to verify that they did not meet any of the following exclusion criteria: (i) sign or history of neurologic, psychiatric, cardiovascular, pleuro-pulmonary or hematologic disease; (ii) ongoing neuroleptic, antiparkinsonian α-methyl-dopa, β-blocker, MAOI-A or -B (MAOI = monoamine oxidase inhibitor), tricyclic antidepressant, or thymoregulator treatment; (iii) pregnancy; (iv) hormone replacement therapy; (v) consumption of recreational drugs (canabis, exctasy); (vi) contraindication to MRI; (vii) brain MRI lesion. Moreover, each subject was asked for his or her consumption of tobacco or alcohol. Heavy consumers of alcohol or tobacco were excluded. Five subjects smoking <10 cigarettes/d and 9 subjects taking 1–2 wine glass units/d (1 wine glass unit = 125 mL of wine according to the Medical Council of Alcohol) were included in the study. Second, all subjects were evaluated for depression using the General Health Questionnaire (GHQ-28) by Goldberg and Hillier (14). Only subjects with a score below the threshold for depression (<7) were included in the study.
MRI
A 3-dimensional multiplanar reconstruction anatomic MRI scan was performed on each subject, yielding a volume of 130–170 slices with cubic voxels of 1 mm3 with 256 × 256 voxels in each transverse plane.
PET Scan of 18F-MPPF
The methodology of tracer production, scan acquisition, and data processing has been described previously (15,16). Briefly, 18F-MPPF was acquired in 35 frames with a CTI Exact HR+ scanner for 60 min after the injection of 153–250 MBq of 18F-MPPF (mean ± SD, 192.7 ± 23.8 MBq). Sinograms were normalized, attenuated, and scatter corrected and then subsequently reconstructed with filtered backprojection (Hanning filter: cutoff, 0.5 cycle/voxel). This yielded to a dynamic study of 35 volumes of 128 × 128 × 63 with a voxel size of 2.04 × 2.04 × 2.42 mm3. The binding parameters of the tracer were estimated according to the 3-compartment SRTM previously applied to 11C-WAY100635 (17). This model is based on the analytic solution of the compartment model that is used to estimate 3 parameters without the use of an arterial sampling input function: R1 (ratio of plasma to brain transport constant in the region of interest [ROI] and in the reference region), k2 (tracer’s efflux in the vascular system), and BP (ratio of available receptor density to receptor affinity, BP = Bmax/Kd, where Kd is the dissociation constant). The cerebellum was taken as the reference region. The SRTM has been evaluated with a correlation study between parameters obtained with a complete nonlinear 3-compartement model (10) and the SRTM BP calculated in 20 ROIs.
The images were processed with the SRTM 2 different ways: one based on mean ROI activity curves (ROI analysis) and the other based on voxel level and SPM (parametric image analysis).
ROI Analysis.
Individual 18F-MPPF static images (summed from 0 to 60 min) were resliced in a transverse orientation parallel to the bihippocampal plane to optimally outline the mesial and lateral temporal regions. The same rotation parameters were applied to the dynamic acquisition. Independently, the MR image was coregistered with the 18F-MPPF static image by an automated method based on Mutual Information criteria (SPM; Welcome Department of Cognitive Neurology, London, U.K.). Anatomic ROIs were manually drawn on the coregistered MR image by one of the authors. Four hundred ROIs were drawn and regrouped into anatomic volumes of interest (VOIs) to describe left and right amygdala, hippocampi, parahippocampal gyri, temporal poles, insula, anterior and posterior cingulate gyri, orbitofrontal cortices, inferior, middle, and superior temporal gyri, subcentral operculi, inferior and superior parietal lobules, frontal and occipital cortices, and thalami. As raphe nuclei cannot be identified on MRI, this region was outlined directly on the static 18F-MPPF image by thresholding the activity at 80% of the local maximum in the brainstem. This region was a posteriori visualized on the MR image to check for its proper location in the periacqueductal gray matter. All ROIs were then applied to dynamic PET to extract time–activity curves and to obtain mean regional values of R1, k2, and BP for each ROI.
Parametric Image Analysis.
For each subject, the MR image was coregistered with the static PET image in its acquisition orientation. On the coregistered MR image, a large ROI was outlined in the cerebellum and was used as a unique reference region for the simplified model applied to each voxel. From individual voxel time–activity curves, the parametric images of R1, k2, and BP were determined. Individual parametric images were then transformed into a standard space using the nonlinear transformation matrix derived from the spatial normalization of the individual’s MR image to the T1 MRI default template (MNI template of the ICBM Project), and smoothed using an 8 × 8 × 8 mm3 full width at half maximum isotropic gaussian kernel that takes into account interindividual anatomy variability and improves the sensitivity of the statistical analysis.
Statistical Inference.
For the ROI analysis, 2 separate statistical analyses for men and women were performed on mean regional BP values considered as independent measures.
First, we performed a linear regression of BP with age for each anatomic ROI with left and right BP values averaged. The r2 coefficient was taken as the index of correlation with a standard Fisher z-score and P for significant assessment. To evaluate the impact of partial-volume effect (PVE) on BP values, we calculated a linear regression of BP versus ROI volumes.
We then performed an ANOVA for each anatomic ROI using laterality (left and right hemisphere) as an independent factor and the decade of age (20s, 30s, 40s, 50s, and 60s) as a group factor.
SPM99 was used for the voxel-based analysis on the normalized smoothed parametric BP images of the 53 control subjects. Three linear models were used to assess the variability due to sex or age between groups. In analysis I, we applied a 2-sample t test between males and females with individual BP values taken as covariates for interindividual adjustment. In analysis II, the design was an analysis of covariance (ANCOVA) with age considered as a covariate, sex as a group factor, and individual global BP values taken as covariates for interindividual adjustment. To perform analysis III, we used an ANOVA with sex and decade of age as group factors and individual global BP as a covariate for adjustment. Statistical parametric maps of the t statistic (SPM{t}) were calculated with a threshold of P = 0.001 uncorrected at the voxel level. Significant clusters were selected at a corrected cluster level of P < 0.05 determined from a joint probability of peak height and cluster size (18).
RESULTS
Relationship Between BP and ROI Volumes
No significant correlation between BP and ROI volumes was found for any regions except in the right insula (P < 0.04) of females and in the raphe nuclei (P < 0.035) of males. However, in these regions, BP variations were inversely correlated with ROI volumes.
Normative Data for ROI Analysis
Table 1 reports BP values for 18 anatomic regions in men and women. As there was no significant difference between left and right BP, the left and right values were mixed to determine the mean BP and the SD in each ROI for each decade of age. Minimal BP values were found in the thalamus and maxima in the hippocampus, with a highest mean of 1.27 ± 0.07 for women in their 30s and 1.28 ± 0.12 for men in their 20s. Minima were found in women in their 60s and in the males in their 30s. BP mean values ranged between 0.5 and 1.2 in limbic areas, between 0.5 and 0.6 in the lateral temporal cortex, between 0.3 and 0.6 in other neocortical regions, and between 0.3 ± 0.1 and 0.6 ± 0.1 in the raphe nuclei. Scatter plots of BP values versus age are illustrated in Figure 1 for all 18 ROIs.
Sex Difference in 18F-MPPF Binding
For the SPM analysis I, the unpaired t test comparing BP images in males and females did not reveal any difference in 18F-MPPF binding. However, when considering age as a covariate (analysis II), BP appeared significantly greater for women in several areas located in the limbic, paralimbic, and frontal regions (Table 2; Fig. 2A). The size of significant clusters was bigger in the right than in the left hemisphere (Fig. 2B). For men, BP was significantly greater in a single cluster, which appeared to be related to a side effect due to the normalization procedure (Fig. 2C).
Linear Variation of BP with Age
Females.
The SPM analysis II, exploring the linear regression between age and 18F-MPPF binding, showed a significant BP decrease with age involving most cortical areas (Fig. 2D). The most significantly correlated areas were found in the superior and middle temporal gyri bilaterally, the left middle frontal gyrus, the insula, the temporooccipital region, and the cingulate gyrus (Table 3; Fig. 2E), sparing only the left amygdala and the orbitofrontal cortex. The inverse contrast did not reveal any significant BP increase.
The ROI analysis, illustrated on Figure 1, showed a significant BP decrease with age in all ROIs (hippocampus, parahippocampal gyrus, insula, anterior and posterior cingulate gyrus, orbitofrontal and frontal gyrus, temporal, parietal, and occipital lobes, and raphe nucleus), except in amygdala and temporal pole. Global BP values decreased significantly with age in the female population with a slope of 3.6% per decade and a mean reduction in 18F-MPPF BP of 15% from 20 to 60 y.
Males.
SPM analysis II revealed a significant BP decrease with age in small clusters of voxels located in the superior frontal gyrus and in the precuneus (Table 3). The inverse contrast showed a significant BP increase in the left middle frontal gyrus and in the rectus gyrus bilaterally (Fig. 2F).
The ROI analysis showed a significant positive linear correlation between age and BP in amygdala and temporal pole. Conversely, a significant negative linear correlation was found in the superior temporal gyri, in the opercular region, in the inferior and superior parietal cortices, and in the occipital lobe. Global BP values were not correlated with age (Fig. 1).
Evolution in 18F-MPPF BP with Decades
From the ANOVA of ROI results, where decade of age and lateralization were considered as independent factors, we did not find any significant influence of the lateralization on regional BP values. Conversely, BP was primarily influenced by the age factor in a manner different for male and female groups.
Females.
The statistical analysis of ROI data showed no influence of the decade of age on the variance in the temporal pole and in the inferior temporal, superior parietal, and occipital gyri. In the remaining regions, where the decade factor had a significant influence, the post hoc analysis revealed a significant decrease between subjects in their 20s and their 30s in amygdala, between subjects in their 40s and their 50s in amygdala and operculum, and between subjects in their 50s and their 60s in all ROIs (mean decrease of 17%) except the posterior cingulate, middle temporal, and opercular cortices.
The SPM analysis III (Table 4), revealed a significant BP decrease independent of global variations between subjects in their 40s and their 50s in the right superior temporal gyrus (Fig. 3A) and a significant BP increase between subjects in their 50s and their 60s in the left fusiform gyrus (Fig. 3B). Other comparisons did not show any significant clusters.
Males.
The statistical analysis of ROI data showed a significant influence of decade of age on BP in all regions except in the raphe nucleus, and in the inferior temporal gyrus. In regions where the decade of age had a significant influence on variance, the post hoc analysis revealed (i) a significant BP decrease between the 20–29 and the 30–39 age group involving all ROIs except the amygdala and the temporal pole, and (ii) a significant BP increase between the 30–39 y and the 40–49 age group in all ROIs except in the cingulate, superior parietal and occipital gyri. No difference was found either between subjects in their 40s and 50s or between those in their 50s and 60s.
The SPM analysis III (Table 4) showed (i) a decrease in BP between the 20–29 y and 30–39 y age group in the inferior parietal gyrus and precuneus bilaterally, in the left cingulate gyrus (Fig. 3C), in the left supplementary motor area, in the medial and lateral superior frontal gyrus (Fig. 3D), and (ii) an increase in BP between the 30–39 and the 40–49 age group in the right middle frontal gyrus (Fig. 3E), and left superior frontal gyrus (Fig. 3F).
DISCUSSION
Use of a Large Database
To our knowledge, this study is the first report exploring 5-HT1A receptor binding with 18F-MPPF in a large population including both male and female subjects. Previous PET studies have used another antagonist of 5-HT1A, 11C-WAY100635, in either smaller populations (n = <25) (19–21) or in a large group of only male volunteers (n = 61) (13). Our database is large enough to perform statistical inference of 18F-MPPF binding modifications to 5-HT1A receptors in the context of pathologies or specific cerebral functions involving the serotoninergic system, or to study the effect of drugs on the serotoninergic system. Accordingly, as patients began to be compared with control populations, we could assess significant bilateral abnormalities in 18F-MPPF binding to 5-HT1A receptors using our database in epileptic patients with temporal lobe epilepsy (15,16).
Moreover, the large sampling of female and male subjects over aging permits the study of pathologies such as depression, Alzheimer’ disease, or anorexia, which may differentially affect men versus women or younger versus older patients.
This study continues our previous work (10) in which we quantified a whole set of 18F-MPPF transport, binding, and metabolization parameters. This quantification implied arterial blood sampling and a multiple injection protocol, which was not appropriate for large-scale studies and justified the use of a simplified protocol. In the present work, the SRTM, with the cerebellum as reference, provides 3 parameters related to binding (BP), washout (k2), and delivery (R1) with sufficient stability to be computed as parametric images suitable for statistical voxel-by-voxel analysis. Furthermore, from our spatially normalized static PET images we were able to create a 18F-MPPF template in the MNI stereotactic space, which allows future spatial normalization procedures in the absence of individual anatomic MR images.
Methodologic Issues
A limitation in our database is the small number of subjects included in the 60–69 y-old age group. Indeed, the different size of our population samples impairs the statistical power for the detection of specific differences between each group and might have an impact on our results. However, SDs for males and females in their 60s were of similar magnitude as for other age groups (Table 1). Nevertheless, to obtain a direct-paired comparison of parametric data between patients of the 60s decade (like patients with Alzheimer’s disease) and controls, additional healthy volunteers of this age group should be included.
Another methodologic issue relates to the drawing of ROIs that were manually traced directly on the individual MR images. The uncertainty of the delineation in the cortical and subcortical regions is not suspected. The problem might be located in the raphe nucleus. Given that this structure is not really a nucleus in the cytoarchitectural sense of the term, a delineation should be arbitrarily chosen either on a fixed volume embedding the raphe in the brainstem with an a priori knowledge of its position or with an arbitrary delineation based on the 18F-MPPF activity threshold from a static equilibrium 18F-MPPF image. We chose this second approach, fixing a threshold at 80% of the local maximum activity located in the brainstem. However, this delineation would not be suitable in patients presenting a severely decreased 18F-MPPF binding in the raphe nucleus. In such cases, a parametric image without a priori delineation of the raphe would be more appropriate.
The question of volume delineation of structure contributes to the complication of the PVE and spillover due to the limited resolution of the PET scanner. The correction method based on geometric matrix transfer (22) can be applied. In this study we considered uncorrected results, yet we verified that the variation of volume over subjects could not explain the variation of BP. As in Tauscher et al. (21), we found that ROI volumes were not correlated with BP changes except in 2 regions—that is, the insula in women and the raphe nucleus in men. However, in these 2 regions, linear regression revealed a negative correlation between volume and BP. Therefore, the hypothesis that BP reductions were attributed to a PVE induced by ROI volume reduction was ruled out. This does not mean that no PVE affected our measures but that PVE cannot explain BP reduction in our study.
Sex Difference in 18F-MPPF Binding
Sex differences are observed in several serotonin-regulated behaviors (sexual, aggression) or mental disorders (eating disorders, anxiety disorders, depression). As an example, the rate of mood disorders was shown to be almost double in females versus that in males (23).
To our knowledge, only 2 PET studies explored male/females differences in 5-HT1A binding using 11C-WAY100635 (20,21). Parsey et al. (20) found higher values of distribution volume in specific compartments (V3) for women compared with men. When considering the SRTM BP (i.e., V3 normalized to the nonspecific distribution volume), this sex effect disappeared as in Tauscher et al. (21). In our study, where BP was calculated with the SRTM, the 2-sample t test comparing parametric images of BP in male versus female groups consistently did not show any significant difference. However, once the effect of age had been removed by an ANCOVA, significantly higher BP values in females than in males were shown in limbic, paralimbic, and frontal regions, areas known to be particularly rich in 5-HT1A receptors. It is noteworthy that Tauscher et al. found a “sex × age” interaction in some regions, but that did not survive a correction for multiple comparisons. In our study, statistical differences between women and men are more likely to be attributed to a higher sensitivity resulting from a larger number of subjects (about twice as in Tauscher et al. or Parsey et al. (20)) or to the different properties of 18F-MPPF versus 11C-WAY100635 relative to 5-HT1A receptors.
The higher BP values found in women can be interpreted as a higher density of 5-HT1A receptors in women compared with men; on the other hand, given that 18F-MPPF is sensitive to endogenous variations of 5-HT concentration (24), higher BP values in women could also be interpreted as a higher availability of 5-HT1A receptor sites due to a lower concentration of 5-HT. BP variations could also be interpreted as a Kd variation of 18F-MPPF binding to 5-HT1A. However, 5-HT synthesis has not yet been assessed in vivo in male and female subjects, and, in addition, postmortem studies in humans using 8-hydroxy-2-(di-n-propylamino)tetralin failed to find any sex effect on Bmax (25,26).
Variation of BP with Age
In our study, we showed that aging influenced the BP values in our healthy subjects. This is an expected result as normal aging leads to changes in most neurotransmission systems as a consequence of presynaptic reduction of the neurotransmitter (which could be due to neuronal loss, decrease in synthesis, or increase of metabolism degradation), decrease in the density of postsynaptic receptors, or deficiency in signal transduction mechanisms. Although some studies in humans fail to find any effect of aging on 5-HT1A receptor binding in healthy subjects (27), other postmortem studies in humans (26) show a substantial reduction of binding with age (up to 50%) attributed to a decrease in receptor density with no change in affinity (26). These studies, however, do not provide any information on sex.
Our results revealed a sex-specific effect of age on BP. We found a linear decrease of BP with age in almost all regions in women. In men, changes in BP were region specific with no correlation between BP and age in most limbic and paralimbic regions, with either an increase or a decrease of BP with age in the remaining regions. Among PET studies using 11C-WAY100635 to explore the linear relation between age and BP, the results are contradictory. Some studies found a BP decrease with age but no effect of sex (21) or a BP decrease in men but not in women (19) or no correlation between BP in men (13,28) or in men and women (20). These results, along with ours, most likely reflect differences in population sizes, age range (oldest subjects, 55–80 y), camera field of view, or affinity of the radioligand for 5-HT1A receptors.
However, using comparable methodologies, postmortem studies in humans also give inconsistent results. Although Matsubara et al. (29) found a negative correlation between the 5-HT1A receptor density (Bmax) and age in male but not in female controls, Palego et al. (25) described both significant increases and decreases of Bmax with age in women, whereas men showed a decrease of affinity with age in the occipital regions.
Analysis of BP changes by decade of age adds to the understanding of the results regarding the linear changes with age. This analysis reveals specific patterns of BP changes over the decade of age. In women, our data suggest that there is a widespread decrease of BP with age, which worsens after 50 y. In men, the analysis by decade reveals that BP decreases between 20 and 30 y but that this decrease in compensated by an increase between 30 and 40 y for almost all regions. This would explain the absence of a linear correlation between age and BP in men. Interestingly, the SPM analysis shows that the decrease between 20 and 30 y is even more significant in the parietal, frontal, and precuneus regions, which are precisely the areas were a linear negative correlation between BP and age was identified. These specific patterns of BP changes with age in men are difficult to interpret, primarily as linear changes would be expected as a consequence of neuronal loss during aging. Unfortunately, the only study reporting such specific patterns of BP changes over the decade of age (30) explored 5-HT2A receptors and pooled the results from male and female subjects. Interestingly, however, on the comparison of these data with those of Palego et al. (25) and Dillon et al. (26), the authors suggested that aging effects on the serotoninergic system seem to depend on age groups and that the decrease in 5-HT1A receptor density is more pronounced in younger men and older women. Whether regulation processes occur in men during midlife is difficult to state from our data but might explain differences in the incidence of psychiatric or anxiety disorders in older male and female subjects.
CONCLUSION
Our data revealed that 18F-MPPF binding to 5-HT1A receptors is different in men and women independently of age and that changes of BP with age are sex specific, which corroborates the different epidemiology of psychiatric diseases in elderly. Our results justify the use of a large normative database to eliminate age and sex effects from the statistical analysis of parameters of interest in the study of patient populations with PET and 18F-MPPF.
Acknowledgments
The authors thank the entire medical and radiochemistry staff of the CERMEP for their technical and human assistance. The authors thank Professor Jean-François Pujol for the initiative and his support to this research.
Footnotes
Received Feb. 14, 2005; revision accepted Aug. 1, 2005.
For correspondence or reprints contact: Nicolas Costes, MSc, CERMEP, 59 boulevard Pinel, 69677 Bron Cedex, France.
E-mail: costes{at}cermep.fr