High-Content Analysis of Breast Cancer Using Single-Cell Deep Transfer Learning

High-content analysis has revolutionized cancer drug discovery by identifying substances that alter the phenotype of a cell, which prevents tumor growth and metastasis. The high-resolution biofluorescence images from assays allow precise quantitative measures enabling the distinction of small molecules of a host cell from a tumor. In this work, we are particularly interested in the application of deep neural networks (DNNs), a cutting-edge machine learning method, to the classification of compounds in chemical mechanisms of action (MOAs). Compound classification has been performed using image-based profiling methods sometimes combined with feature reduction methods such as principal component analysis or factor analysis. In this article, we map the input features of each cell to a particular MOA class without using any treatment-level profiles or feature reduction methods. To the best of our knowledge, this is the first application of DNN in this domain, leveraging single-cell information. Furthermore, we use deep transfer learning (DTL) to alleviate the intensive and computational demanding effort of searching the huge parameter's space of a DNN. Results show that using this approach, we obtain a 30% speedup and a 2% accuracy improvement.

Different Electroclinical Manifestations of the Epilepsy Associated with Hamartomas Connecting to the Middle or Posterior Hypothalamus

PURPOSE: The epilepsy associated with hypothalamic hamartomas (HHs) has typical clinical, electrophysiologic, and behavioral manifestations refractory to drug therapy and with unfavorable evolution. It is well known that only sessile lesions produce epilepsy, but no correlation has been established between the different types of sessile hamartomas and the diverse manifestations of the epilepsy. We correlate anatomic details of the hamartoma and the clinical and neurophysiologic manifestations of the associated epilepsy. METHODS: HHs of seven patients with epilepsy (ages 2- 25 years) were classified as to lateralization and connection to the anteroposterior axis of the hypothalamus by using high-resolution brain magnetic resonance imaging. We correlated the anatomic classification with the clinical and neurophysiologic manifestations of the epilepsy as evaluated in long-term (24 h) video-EEG recordings. RESULTS: HHs ranged in size from 0.4 to 2.6 cc, with complete lateralization in six of seven patients. Ictal manifestations showed good correlation with the lobar involvement of ictal/interictal EEGs. These manifestations suggest the existence of two types of cortical involvement, one associated with the temporal lobe, produced by hamartomas connected to the posterior hypothalamus (mamillary bodies), and the other associated with the frontal lobe, seen in lesions connecting to the middle hypothalamus. CONCLUSIONS: A consistent clinical and neurophysiologic pattern of either temporal or frontal lobe cortical secondary involvement was found in the patients of our series. It depends on whether the hamartoma connects to the mamillary bodies (temporal lobe cases) or whether it connects to the medial hypothalamus (frontal lobe cases).

Addressing Sea Surface Salinity retrieval at different spatial resolution. Feasibility study using OCCAM data within SEPS simulator and L2 processor tools

Projecte de recerca elaborat a partir d’una estada a la National Oceanography Centre of Southampton (NOCS), Gran Bretanya, entre maig i juliol del 2006. La possibilitat d’obtenir una estimació precissa de la salinitat marina (SSS) és important per a investigar i predir l’extensió del fenòmen del canvi climàtic. La missió Soil Moisture and Ocean Salinity (SMOS) va ser seleccionada per l’Agència Espacial Europea (ESA) per a obtenir mapes de salinitat de la superfície marina a escala global i amb un temps de revisita petit. Abans del llençament de SMOS es preveu l’anàlisi de la variabilitat horitzontal de la SSS i del potencial de les dades recuperades a partir de mesures de SMOS per a reproduir comportaments oceanogràfics coneguts. L’objectiu de tot plegat és emplenar el buit existent entre les fonts de dades d’entrada/auxiliars fiables i les eines desenvolupades per a simular i processar les dades adquirides segons la configuració de SMOS. El SMOS End-to-end Performance Simulator (SEPS) és un simulador adhoc desenvolupat per la Universitat Politècnica de Catalunya (UPC) per a generar dades segons la configuració de SMOS. Es va utilitzar dades d’entrada a SEPS procedents del projecte Ocean Circulation and Climate Advanced Modeling (OCCAM), utilitzat al NOCS, a diferents resolucions espacials. Modificant SEPS per a poder fer servir com a entrada les dades OCCAM es van obtenir dades de temperatura de brillantor simulades durant un mes amb diferents observacions ascendents que cobrien la zona seleccionada. Les tasques realitzades durant l’estada a NOCS tenien la finalitat de proporcionar una tècnica fiable per a realitzar la calibració externa i per tant cancel•lar el bias, una metodologia per a promitjar temporalment les diferents adquisicions durant les observacions ascendents, i determinar la millor configuració de la funció de cost abans d’explotar i investigar les posibiltats de les dades SEPS/OCCAM per a derivar la SSS recuperada amb patrons d’alta resolució.

Some solutions to obtain very efficient separations in isocratic and gradient modes using small particles size and ultra-high pressure.

The UHPLC strategy which combines sub-2 microm porous particles and ultra-high pressure (>1000 bar) was investigated considering very high resolution criteria in both isocratic and gradient modes, with mobile phase temperatures between 30 and 90 degrees C. In isocratic mode, experimental conditions to reach the maximal efficiency were determined using the kinetic plot representation for DeltaP(max)=1000 bar. It has been first confirmed that the molecular weight of the compounds (MW) was a critical parameter which should be considered in the construction of such curves. With a MW around 1000 g mol(-1), efficiencies as high as 300,000 plates could be theoretically attained using UHPLC at 30 degrees C. By limiting the column length to 450 mm, the maximal plate count was around 100,000. In gradient mode, the longest column does not provide the maximal peak capacity for a given analysis time in UHPLC. This was attributed to the fact that peak capacity is not only related to the plate number but also to column dead time. Therefore, a compromise should be found and a 150 mm column should be preferentially selected for gradient lengths up to 60 min at 30 degrees C, while the columns coupled in series (3x 150 mm) were attractive only for t(grad)>250 min. Compared to 30 degrees C, peak capacities were increased by about 20-30% for a constant gradient length at 90 degrees C and gradient time decreased by 2-fold for an identical peak capacity.

MP2RAGE, a self bias-field corrected sequence for improved segmentation and T-1-mapping at high field

The large spatial inhomogeneity in transmit B, field (B-1(+)) observable in human MR images at hi h static magnetic fields (B-0) severely impairs image quality. To overcome this effect in brain T-1-weighted images the, MPRAGE sequence was modified to generate two different images at different inversion times MP2RAGE By combining the two images in a novel fashion, it was possible to create T-1-weigthed images where the result image was free of proton density contrast, T-2* contrast, reception bias field, and, to first order transmit field inhomogeneity. MP2RAGE sequence parameters were optimized using Bloch equations to maximize contrast-to-noise ratio per unit of time between brain tissues and minimize the effect of B-1(+) variations through space. Images of high anatomical quality and excellent brain tissue differentiation suitable for applications such as segmentation and voxel-based morphometry were obtained at 3 and 7 T. From such T-1-weighted images, acquired within 12 min, high-resolution 3D T-1 maps were routinely calculated at 7 T with sub-millimeter voxel resolution (0.65-0.85 mm isotropic). T-1 maps were validated in phantom experiments. In humans, the T, values obtained at 7 T were 1.15 +/- 0.06 s for white matter (WM) and 1.92 +/- 0.16 s for grey matter (GM), in good agreement with literature values obtained at lower spatial resolution. At 3 T, where whole-brain acquisitions with 1 mm isotropic voxels were acquired in 8 min the T-1 values obtained (0.81 +/- 0.03 S for WM and 1.35 +/- 0.05 for GM) were once again found to be in very good agreement with values in the literature. (C) 2009 Elsevier Inc. All rights reserved.

Targetting of the ventro-intermediate nucleus using ultra-high field (7T) MRI for gamma knife surgery purposes: A pilot in vivo study on healthy subjects.

INTRODUCTION: Gamma Knife surgery (GKS) is a non-invasive neurosurgical stereotactic procedure, increasingly used as an alternative to open functional procedures. This includes targeting of the ventro-intermediate nucleus of the thalamus (e.g. Vim) for tremor. We currently perform an indirect targeting, as the Vim is not visible on current 3Tesla MRI acquisitions. Our objective was to enhance anatomic imaging (aiming at refining the precision of anatomic target selection by direct visualisation) in patients treated for tremor with Vim GKS, by using high field 7T MRI. MATERIALS AND METHODSH: Five young healthy subjects were scanned on 3 (T1-w and diffusion tensor imaging) and 7T (high-resolution susceptibility weighted images (SWI)) MRI in Lausanne. All images were further integrated for the first time into the Gamma Plan Software(®) (Elekta Instruments, AB, Sweden) and co-registered (with T1 was a reference). A simulation of targeting of the Vim was done using various methods on the 3T images. Furthermore, a correlation with the position of the found target with the 7T SWI was performed. The atlas of Morel et al. (Zurich, CH) was used to confirm the findings on a detailed analysis inside/outside the Gamma Plan. RESULTS: The use of SWI provided us with a superior resolution and an improved image contrast within the basal ganglia. This allowed visualization and direct delineation of some subgroups of thalamic nuclei in vivo, including the Vim. The position of the target, as assessed on 3T, perfectly matched with the supposed one of the Vim on the SWI. Furthermore, a 3-dimensional model of the Vim-target area was created on the basis of the obtained images. CONCLUSION: This is the first report of the integration of SWI high field MRI into the LGP, aiming at the improvement of targeting validation of the Vim in tremor. The anatomical correlation between the direct visualization on 7T and the current targeting methods on 3T (e.g. quadrilatere of Guyot, histological atlases) seems to show a very good anatomical matching. Further studies are needed to validate this technique, both by improving the accuracy of the targeting of the Vim (potentially also other thalamic nuclei) and to perform clinical assessment.

The feasibility of 350 μm spatial resolution coronary magnetic resonance angiography at 3 T in humans.

PURPOSE: The purposes of this study were to (1) develop a high-resolution 3-T magnetic resonance angiography (MRA) technique with an in-plane resolution approximate to that of multidetector coronary computed tomography (MDCT) and a voxel size of 0.35 × 0.35 × 1.5 mm³ and to (2) investigate the image quality of this technique in healthy participants and preliminarily in patients with known coronary artery disease (CAD). MATERIALS AND METHODS: A 3-T coronary MRA technique optimized for an image acquisition voxel as small as 0.35 × 0.35 × 1.5 mm³ (high-resolution coronary MRA [HRC]) was implemented and the coronary arteries of 22 participants were imaged. These included 11 healthy participants (average age, 28.5 years; 5 men) and 11 participants with CAD (average age, 52.9 years; 5 women) as identified on MDCT. In addition, the 11 healthy participants were imaged using a method with a more common spatial resolution of 0.7 × 1 × 3 mm³ (regular-resolution coronary MRA [RRC]). Qualitative and quantitative comparisons were made between the 2 MRA techniques. RESULTS: Normal vessels and CAD lesions were successfully depicted at 350 × 350 μm² in-plane resolution with adequate signal-to-noise ratio (SNR) and contrast-to-noise ratio. The CAD findings were consistent among MDCT and HRC. The HRC showed a 47% improvement in sharpness despite a reduction in SNR (by 72%) and in contrast-to-noise ratio (by 86%) compared with the regular-resolution coronary MRA. CONCLUSION: This study, as a first step toward substantial improvement in the resolution of coronary MRA, demonstrates the feasibility of obtaining at 3 T a spatial resolution that approximates that of MDCT. The acquisition in-plane pixel dimensions are as small as 350 × 350 μm² with a 1.5-mm slice thickness. Although SNR is lower, the images have improved sharpness, resulting in image quality that allows qualitative identification of disease sites on MRA consistent with MDCT.

MP2RAGE, a self bias-field corrected sequence for improved segmentation and T1-mapping at high field.

The large spatial inhomogeneity in transmit B(1) field (B(1)(+)) observable in human MR images at high static magnetic fields (B(0)) severely impairs image quality. To overcome this effect in brain T(1)-weighted images, the MPRAGE sequence was modified to generate two different images at different inversion times, MP2RAGE. By combining the two images in a novel fashion, it was possible to create T(1)-weighted images where the result image was free of proton density contrast, T(2) contrast, reception bias field, and, to first order, transmit field inhomogeneity. MP2RAGE sequence parameters were optimized using Bloch equations to maximize contrast-to-noise ratio per unit of time between brain tissues and minimize the effect of B(1)(+) variations through space. Images of high anatomical quality and excellent brain tissue differentiation suitable for applications such as segmentation and voxel-based morphometry were obtained at 3 and 7 T. From such T(1)-weighted images, acquired within 12 min, high-resolution 3D T(1) maps were routinely calculated at 7 T with sub-millimeter voxel resolution (0.65-0.85 mm isotropic). T(1) maps were validated in phantom experiments. In humans, the T(1) values obtained at 7 T were 1.15+/-0.06 s for white matter (WM) and 1.92+/-0.16 s for grey matter (GM), in good agreement with literature values obtained at lower spatial resolution. At 3 T, where whole-brain acquisitions with 1 mm isotropic voxels were acquired in 8 min, the T(1) values obtained (0.81+/-0.03 s for WM and 1.35+/-0.05 for GM) were once again found to be in very good agreement with values in the literature.

High prevalence of peripheral atherosclerosis in a rapidly developing country.

Cardiovascular disease is rapidly increasing in developing countries experiencing epidemiological transition. We investigated the prevalence of peripheral atherosclerosis in a rapidly developing country and compared our findings with data previously reported in Western populations. A cardiovascular risk factor survey was conducted in 1067 individuals aged 25-64 randomly selected from the general population of Seychelles. High-resolution ultrasonography of the right and left carotid and femoral arteries was performed in a random subgroup of 503 subjects (245 men and 258 women). In each of the four arteries, arterial wall thickness (in plaque-free segments) and atherosclerotic plaques (i.e. focal wall thickening at least 1.0 mm thick) were measured separately. The prevalence of peripheral atherosclerosis was high in this population. For instance, at least one plaque > or =1.0 mm was found in, respectively, 34.9 and 27.5% of men and women aged 25-34 and at least one plaque > or =2.5 mm was found in, respectively, 58.2 and 36.9% of men and women aged 55-64. With reference to data found in the literature, the prevalence of carotid atherosclerosis appeared to be significantly higher in Seychelles than in Western populations. This study provides further evidence for the importance of cardiovascular disease in developing countries. Determinants should be identified and relevant prevention and control programs implemented.

Psychogenic seizures and frontal disconnection: EEG synchronisation study.

Objective Psychogenic non-epileptic seizures (PNES) are paroxysmal events that, in contrast to epileptic seizures, are related to psychological causes without the presence of epileptiform EEG changes. Recent models suggest a multifactorial basis for PNES. A potentially paramount, but currently poorly understood factor is the interplay between psychiatric features and a specific vulnerability of the brain leading to a clinical picture that resembles epilepsy. Hypothesising that functional cerebral network abnormalities may predispose to the clinical phenotype, the authors undertook a characterisation of the functional connectivity in PNES patients. Methods The authors analysed the whole-head surface topography of multivariate phase synchronisation (MPS) in interictal high-density EEG of 13 PNES patients as compared with 13 age- and sex-matched controls. MPS mapping reduces the wealth of dynamic data obtained from high-density EEG to easily readable synchronisation maps, which provide an unbiased overview of any changes in functional connectivity associated with distributed cortical abnormalities. The authors computed MPS maps for both Laplacian and common-average-reference EEGs. Results In a between-group comparison, only patchy, non-uniform changes in MPS survived conservative statistical testing. However, against the background of these unimpressive group results, the authors found widespread inverse correlations between individual PNES frequency and MPS within the prefrontal and parietal cortices. Interpretation PNES appears to be associated with decreased prefrontal and parietal synchronisation, possibly reflecting dysfunction of networks within these regions.

Ultra-High Field (7 Tesla) MRI for Gamma Knife Surgery Targeting of the Ventro-Intermediate Nucleus: A Pilot in Vivo Study

Aim: Gamma Knife surgery (GKS) is a non-invasive neurosurgical stereotactic procedure, increasingly used as an alternative to open functional procedures. This includes the targeting of the ventro-intermediate (Vim) nucleus of the thalamus for tremor. We currently perform an indirect targeting, using the "quadrilatere of Guyot," as the Vim nucleus is not visible on current 3 Tesla (T) MRI acquisitions. The primary objective of the current study was to enhance anatomic imaging for Vim GKS using high-field (7 T) MRI, with the aim of refining the visualization and precision of anatomical targeting. Method: Five young healthy subjects (mean age 23 years) were scanned both on 3 and 7 T MRI in Lausanne University Hospital (CHUV) and Center for Biomedical Imaging (CIBM). Classical T1-weighted MPRAGE, T2 CISS sequences (replacing former ventriculography) and diffusion tensor imaging were acquired at 3T. We obtained high-resolution susceptibility weighted images (SWI) at 7T for the visualization of thalamic subparts. SWI was further integrated for the first time into Leksell Gamma Plan® (LGP) software and co-registered with the 3T images. A simulation of targeting of the Vim was done using the "quadrilatere of Guyot" methodology on the 3T images. Furthermore, a correlation with the position of the found target on SWI was performed. The atlas of Morel et al. was used to confirm the findings on a detailed computer analysis outside LGP. Also, 3T and 7T MRI of one patient undergoing GKS Vim thalamotomy, were obtained before and 2 years after the procedure, and studied similarly. Results: The use of SWI provided a superior resolution and improved image contrast within the central gray matter. This allowed visualization and direct delineation of groups of thalamic nuclei in vivo, including the Vim. The position of the target, as assessed with the "quadrilatere of Guyot" method on 3 T, perfectly matched with the supposed one of the Vim on the SWI. Furthermore, a 3-dimensional model of the Vim target area was created on the basis of 3T and 7T images. Conclusion: This is the first report of the integration of SWI high-field MRI into the LGP in healthy subjects and in one patient treated GKS Vim thalamotomy. This approach aims at the improvement of targeting validation and further direct targeting of the Vim in tremor. The anatomical correlation between the direct visualization on 7T and the current targeting methods on 3T seems to show a very good anatomical matching.

A post-processing pipeline to reconstruct the developing fetal brain using low-resolution MRI

Introduction. Development of the fetal brain surfacewith concomitant gyrification is one of the majormaturational processes of the human brain. Firstdelineated by postmortem studies or by ultrasound, MRIhas recently become a powerful tool for studying in vivothe structural correlates of brain maturation. However,the quantitative measurement of fetal brain developmentis a major challenge because of the movement of the fetusinside the amniotic cavity, the poor spatial resolution,the partial volume effect and the changing appearance ofthe developing brain. Today extensive efforts are made todeal with the âeurooepost-acquisitionâeuro reconstruction ofhigh-resolution 3D fetal volumes based on severalacquisitions with lower resolution (Rousseau, F., 2006;Jiang, S., 2007). We here propose a framework devoted tothe segmentation of the basal ganglia, the gray-whitetissue segmentation, and in turn the 3D corticalreconstruction of the fetal brain. Method. Prenatal MRimaging was performed with a 1-T system (GE MedicalSystems, Milwaukee) using single shot fast spin echo(ssFSE) sequences in fetuses aged from 29 to 32gestational weeks (slice thickness 5.4mm, in planespatial resolution 1.09mm). For each fetus, 6 axialvolumes shifted by 1 mm were acquired (about 1 min pervolume). First, each volume is manually segmented toextract fetal brain from surrounding fetal and maternaltissues. Inhomogeneity intensity correction and linearintensity normalization are then performed. A highspatial resolution image of isotropic voxel size of 1.09mm is created for each fetus as previously published byothers (Rousseau, F., 2006). B-splines are used for thescattered data interpolation (Lee, 1997). Then, basalganglia segmentation is performed on this superreconstructed volume using active contour framework witha Level Set implementation (Bach Cuadra, M., 2010). Oncebasal ganglia are removed from the image, brain tissuesegmentation is performed (Bach Cuadra, M., 2009). Theresulting white matter image is then binarized andfurther given as an input in the Freesurfer software(http://surfer.nmr.mgh.harvard.edu/) to provide accuratethree-dimensional reconstructions of the fetal brain.Results. High-resolution images of the cerebral fetalbrain, as obtained from the low-resolution acquired MRI,are presented for 4 subjects of age ranging from 29 to 32GA. An example is depicted in Figure 1. Accuracy in theautomated basal ganglia segmentation is compared withmanual segmentation using measurement of Dice similarity(DSI), with values above 0.7 considering to be a verygood agreement. In our sample we observed DSI valuesbetween 0.785 and 0.856. We further show the results ofgray-white matter segmentation overlaid on thehigh-resolution gray-scale images. The results arevisually checked for accuracy using the same principlesas commonly accepted in adult neuroimaging. Preliminary3D cortical reconstructions of the fetal brain are shownin Figure 2. Conclusion. We hereby present a completepipeline for the automated extraction of accuratethree-dimensional cortical surface of the fetal brain.These results are preliminary but promising, with theultimate goal to provide âeurooemovieâeuro of the normal gyraldevelopment. In turn, a precise knowledge of the normalfetal brain development will allow the quantification ofsubtle and early but clinically relevant deviations.Moreover, a precise understanding of the gyraldevelopment process may help to build hypotheses tounderstand the pathogenesis of several neurodevelopmentalconditions in which gyrification have been shown to bealtered (e.g. schizophrenia, autismâeuro¦). References.Rousseau, F. (2006), 'Registration-Based Approach forReconstruction of High-Resolution In Utero Fetal MR Brainimages', IEEE Transactions on Medical Imaging, vol. 13,no. 9, pp. 1072-1081. Jiang, S. (2007), 'MRI of MovingSubjects Using Multislice Snapshot Images With VolumeReconstruction (SVR): Application to Fetal, Neonatal, andAdult Brain Studies', IEEE Transactions on MedicalImaging, vol. 26, no. 7, pp. 967-980. Lee, S. (1997),'Scattered data interpolation with multilevel B-splines',IEEE Transactions on Visualization and Computer Graphics,vol. 3, no. 3, pp. 228-244. Bach Cuadra, M. (2010),'Central and Cortical Gray Mater Segmentation of MagneticResonance Images of the Fetal Brain', ISMRM Conference.Bach Cuadra, M. (2009), 'Brain tissue segmentation offetal MR images', MICCAI.

In vivo assessment of myelination by phase imaging at high magnetic field.

The present study evaluated the potential of using the phase of T2* weighted MR images to characterize myelination during brain development and pathology in rodents at 9.4 T. Phase contrast correlated with myelin content assessed by histology and suggests that most contrast between white and cortical gray matter is modulated by myelin. Ex vivo experiments showed that gray-white matter phase contrast remains unchanged after iron extraction. In dysmyelinated shiverer mice, phase imaging correlated strongly with myelin staining, showing reduced contrast between white and gray matter when compared to healthy controls. We conclude that high-resolution phase images, acquired at high field, allow assessment of myelination and dysmyelination.

Crystallization and resolution of the lipoxygenase of Pseudomonas aeruginosa 42A2 and phylogenetic study of the subfamilies of the lipoxygenases

Lipoxygenases are non-heme iron enzymes essential in eukaryotes, where they catalyze the formation of the fatty acid hydroperoxides that are required by a large diversity of biological and pathological processes. In prokaryotes, most of them totally lacking in polyunsaturated fatty acids, the possible biological roles oflipoxygenases have remained obscure. In this study, it is reported the crystallization of a lipoxygenase of Pseudomonas aeruginosa (Pa_LOX), the first from a prokaryote. High resolution data has been acquired which is expected to yield structural clues to the questions adressed. Besides, a preliminar phylogenetic analysis using 14 sequences has confirmed the existence of this subfamily of bacterial lipoxygenases, on one side, and a greater diversity than in the corresponding eukaryotic ones, on the other. Finally, an evolutionary study of bacteriallipoxygenases on the same set of lipoxygenases, show a selection pressure of a basically purifying or neutral character except for a single aminoacid, which would have been selected after a positive selection event.

Crystallization and resolution of the lipoxygenase of Pseudomonas aeruginosa 42A2 and phylogenetic study of the subfamilies of the lipoxygenases

Lipoxygenases are non-heme iron enzymes essential in eukaryotes, where they catalyze the formation of the fatty acid hydroperoxides that are required by a large diversity of biological and pathological processes. In prokaryotes, most of them totally lacking in polyunsaturated fatty acids, the possible biological roles oflipoxygenases have remained obscure. In this study, it is reported the crystallization of a lipoxygenase of Pseudomonas aeruginosa (Pa_LOX), the first from a prokaryote. High resolution data has been acquired which is expected to yield structural clues to the questions adressed. Besides, a preliminar phylogenetic analysis using 14 sequences has confirmed the existence of this subfamily of bacterial lipoxygenases, on one side, and a greater diversity than in the corresponding eukaryotic ones, on the other. Finally, an evolutionary study of bacteriallipoxygenases on the same set of lipoxygenases, show a selection pressure of a basically purifying or neutral character except for a single aminoacid, which would have been selected after a positive selection event.