Dense Deformation Field Estimation for Atlas Registration using the Active Contour Framework

In this paper, we propose a new paradigm to carry outthe registration task with a dense deformation fieldderived from the optical flow model and the activecontour method. The proposed framework merges differenttasks such as segmentation, regularization, incorporationof prior knowledge and registration into a singleframework. The active contour model is at the core of ourframework even if it is used in a different way than thestandard approaches. Indeed, active contours are awell-known technique for image segmentation. Thistechnique consists in finding the curve which minimizesan energy functional designed to be minimal when thecurve has reached the object contours. That way, we getaccurate and smooth segmentation results. So far, theactive contour model has been used to segment objectslying in images from boundary-based, region-based orshape-based information. Our registration technique willprofit of all these families of active contours todetermine a dense deformation field defined on the wholeimage. A well-suited application of our model is theatlas registration in medical imaging which consists inautomatically delineating anatomical structures. Wepresent results on 2D synthetic images to show theperformances of our non rigid deformation field based ona natural registration term. We also present registrationresults on real 3D medical data with a large spaceoccupying tumor substantially deforming surroundingstructures, which constitutes a high challenging problem.

A novel method for automated classification of epileptiform activity in the human electroencephalogram-based on independent component analysis.

Diagnosis of several neurological disorders is based on the detection of typical pathological patterns in the electroencephalogram (EEG). This is a time-consuming task requiring significant training and experience. Automatic detection of these EEG patterns would greatly assist in quantitative analysis and interpretation. We present a method, which allows automatic detection of epileptiform events and discrimination of them from eye blinks, and is based on features derived using a novel application of independent component analysis. The algorithm was trained and cross validated using seven EEGs with epileptiform activity. For epileptiform events with compensation for eyeblinks, the sensitivity was 65 +/- 22% at a specificity of 86 +/- 7% (mean +/- SD). With feature extraction by PCA or classification of raw data, specificity reduced to 76 and 74%, respectively, for the same sensitivity. On exactly the same data, the commercially available software Reveal had a maximum sensitivity of 30% and concurrent specificity of 77%. Our algorithm performed well at detecting epileptiform events in this preliminary test and offers a flexible tool that is intended to be generalized to the simultaneous classification of many waveforms in the EEG.

The Effect Of Induced Intraocular Straylight (cataract Simulation) On Sensitivity Measures Of Standard Automated Perimetry, Pulsar Perimetry And Moorfields Motion Displacement Test

Purpose: To investigate the effect of incremental increases in intraocular straylight on threshold measurements made by three modern forms of perimetry: Standard Automated Perimetry (SAP) using Octopus (Dynamic, G-Pattern), Pulsar Perimetry (PP) (TOP, 66 points) and the Moorfields Motion Displacement Test (MDT) (WEBS, 32 points).Methods: Four healthy young observers were recruited (mean age 26yrs [25yrs, 28yrs]), refractive correction [+2 D, -4.25D]). Five white opacity filters (WOF), each scattering light by different amounts were used to create incremental increases in intraocular straylight (IS). Resultant IS values were measured with each WOF and at baseline (no WOF) for each subject using a C-Quant Straylight Meter (Oculus, Wetzlar, Germany). A 25 yr old has an IS value of ~0.85 log(s). An increase of 40% in IS to 1.2log(s) corresponds to the physiological value of a 70yr old. Each WOFs created an increase in IS between 10-150% from baseline, ranging from effects similar to normal aging to those found with considerable cataract. Each subject underwent 6 test sessions over a 2-week period; each session consisted of the 3 perimetric tests using one of the five WOFs and baseline (both instrument and filter were randomised).Results: The reduction in sensitivity from baseline was calculated. A two-way ANOVA on mean change in threshold (where subjects were treated as rows in the block and each increment in fog filters was treated as column) was used to examine the effect of incremental increases in straylight. Both SAP (p<0.001) and Pulsar (p<0.001) were significantly affected by increases in straylight. The MDT (p=0.35) remained comparatively robust to increases in straylight.Conclusions: The Moorfields MDT measurement of threshold is robust to effects of additional straylight as compared to SAP and PP.

A new, automated, four-colour interphase FISH approach for the simultaneous detection of specific aneuploidies of diagnostic and prognostic significance in high hyperdiploid ALL

In hyperdiploid acute lymphoblastic leukaemia (ALL), the simultaneous occurrence of specific aneuploidies confers a more favourable outcome than hyperdiploidy alone. Interphase (I) FISH complements conventional cytogenetics (CC) through its sensitivity and ability to detect chromosome aberrations in non-dividing cells. To overcome the limits of manual I-FISH, we developed an automated four-colour I-FISH approach and assessed its ability to detect concurrent aneuploidies in ALL. I-FISH was performed using centromeric probes for chromosomes 4, 6, 10 and 17. Parameters established for automatic nucleus selection and signal detection were evaluated (3 controls). Cut-off values were determined (10 controls, 1000 nuclei/case). Combinations of aneuploidies were considered relevant when each aneuploidy was individually significant. Results obtained in 10 ALL patients (1500 nuclei/patient) were compared with those by CC. Various combinations of aneuploidies were identified. All clones detected by CC were observed by I-FISH. I-FISH revealed numerous additional abnormal clones, ranging between 0.1% and 31.6%, based on the large number of nuclei evaluated. Four-colour automated I-FISH permits the identification of concurrent aneuploidies of prognostic significance in hyperdiploid ALL. Large numbers of cells can be analysed rapidly by this method. Owing to its high sensitivity, the method provides a powerful tool for the detection of small abnormal clones at diagnosis and during follow up. Compared to CC, it generates a more detailed cytogenetic picture, the biological and clinical significance of which merits further evaluation. Once optimised for a given set of probes, the system can be easily adapted for other probe combinations.

Simultaneous detection of aneuploidies for chromosomes 4, 6, 10 and 17 by automated four colour I-FISH in high hyperdiploid acute lymphoblastic leukemia : diagnostic assessment, clonal heterogeneity and chromosomal instability in adults

SummarySimultaneous detection of aneuploidies for chromosomes 4, 6,10 and 17 by automated four color l-FISH in high hyperdiploid acute lymphoblastic leukemia: diagnostic assessment, clonal heterogeneity and chromosomal instability in adultsAnna Talamo BlandinService de Génétique Médicale, Unité de Cytogénétique du Cancer, CHUVAcute lymphoblastic leukemia (ALL) is a malignant hemopathy characterized by the accumulation of the immature lymphoid cells in the bone marrow and, most often, in the peripheral blood. ALL is a heterogeneous disease with distinct biological and prognostic entities. At diagnosis, cytogenetic and molecular findings constitute important and independent prognostic factors. High hyperdiploidy with 51-67 chromosomes (HeH), one of the largest cytogenetic subsets of ALL, in childhood particularly, is generally associated with a relatively favorable outcome. Chromosome gain is nonrandom, extracopies of some chromosome occurring more frequently than those of others. Concurrent presence of trisomy for chromosomes 4, 10 and 17 confers an especially good prognosis. The first aim of our work was to develop an automated four color interphase fluorescence in situ hybridization (l-FISH) methodology and to assess its ability to detect concurrent aneuploidies 4, 6, 10 and 17 in 10 ALL patients. Various combinations of aneuploidies were identified. All clones detected by conventional cytogenetics were also observed by l-FISH. However, in all patients, l-FISH revealed numerous additional abnormal clones, leading to a high level of clonal heterogeneity. Our second aim has been to investigate the nature and origin of this clonal heterogeneity and to test for the presence of chromosome instability (CIN) in HeH ALL at initial presentation. Ten HeH ALL and 10 non-HeH ALL patients were analysed by four colour l-FISH and numerical CIN values were determined for all four chromosomes together and for each chromosome and patient group, an original approach in ALL. CIN values in HeH ALL proved to be much higher than#iose in non-HeH ALL, suggesting that numerical CIN may be at the origin of the high level of clonal heterogeneity revealed by l-FISH. Our third aim has been to study the evolution of these cytogenetic features during the course of the disease in 10 HeH ALL patients. Clonal heterogeneity was also observed again during disease progression, particularly at relapse. Clones detected at initial presentation generally reappeared in relapse, in most cases with newly generated ones. A significant correlation between the number of abnormal clones and CIN suggested that the higher the instability, the larger the number of abnormal clones. Whereas clonal heterogeneity and its evolution most probably result from underlying chromosome instability, operating processes remain conjectural.RésuméLa leucémie lymphoblastique aiguë (LLA) est une hémopathie maligne qui résulte de l'accumulationde cellules lymphoïdes immatures dans la moelle osseuse, et, le plus souvent, dans le sangpériphérique également. La LLA est une affection hétérogène au sein de laquelle se distinguentplusieurs entités biologiques et pronostiques. Les données cytogénétiques et moléculaires font partieintégrante du diagnostic et jouent un rôle essentiel dans l'évaluation du pronostic. L'hyperdiploïdieélevée à 51-­67 chromosomes (HeH), relativement fréquente, en particulier chez l'enfant, s'associe àun pronostic favorable. Le gain de chromosomes ne relève pas du hasard, certains chromosomesétant plus fréquemment impliqués que d'autres. La présence simultanée des trisomies 4, 6, et 17s'associe à un pronostic particulièrement bon. Le premier but du travail a été de développer uneméthode d'analyse automatique par hybridation in situ fluorescente interphasique (I-­FISH) à 4couleurs et de tester sa capacité à identifier la présence simultanée d'aneuploïdies 4, 6, 10 et 17 dans10 cas de LLA. Différentes combinaisons d'aneuploïdies ont été identifiées. Tous les clones détectéspar cytogénétique conventionnelle l'ont été par I-­FISH. Or, chez tous les patients, l'I-­FISH a révélé denombreux clones anormaux additionnels générant un degré élevé d'hétérogénéité clonale. Notredeuxième but a été d'investiguer la nature et l'origine de cette hétérogénéité et de tester la présenced'instabilité chromosomique (CIN) chez les patients avec une LLA HeH en presentation initiale. DixLLA HeH et 10 LLA non-­HeH ont été analysées par I-­FISH et les valeurs de CIN numérique ont étédéterminées pour les 4 chromosomes ensemble et pour chaque chromosome et groupe de patients,approche originale dans la LLA. Ces valeurs étant beaucoup plus élevées dans la LLA HeH que dansla LLA non-­HeH, elles favorisent l'hypothèse selon laquelle la CIN serait à l'origine de l'hétérogénéitéclonale révélée par I-­FISH. Le troisième but de notre travail a été d'étudier l'évolution de cescaractéristiques cytogénétiques au cours de la maladie dans 10 cas de LLA HeH. L'hétérogénéitéclonale a été retrouvée lors de la progression de la maladie, en particulier en rechute, où les clonesanormaux détectés en présentation initiale réapparaissent, généralement accompagnés de clonesnouveaux. La corrélation existant entre nombre de clones anormaux et valeurs de CIN suggère queplus l'instabilité est élevée, plus le nombre de clones anormaux est grand. Bien que l'hétérogénéitéclonale et son évolution résultent très probablement de l'instabilité chromosomique, les processus àl'oeuvre ne sont pas entièrement élucidés.

Examining the Agreement of Structural Loss Measured With Heidelberg Retinal Tomograhy and Functional Loss With Standard Automated Perimetry (SAP; Octopus), Pulsar Perimetry (PP), and Moorfields Motion Displacement Test (MDT) Perimetry

Purpose: To examine the relationship of functional measurements with structural measures. Methods: 146 eyes of 83 test subjects underwent Heidelberg Retinal Tomography (HRTIII) (disc area<2.43, mphsd<40), and perimetry testing with Octopus (SAP; Dynamic), Pulsar (PP; TOP) and Moorfields MDT (ESTA). Glaucoma was defined as progressive structural or functional loss (20 eyes). Perimetry test points were grouped into 6 sectors based on the estimated optic nerve head angle into which the associated nerve fiber bundle enters (Garway-Heath map). Perimetry summary measures (PSM) (MD SAP/ MD PP/ PTD MDT) were calculated from the average total deviation of each measured threshold from the normal for each sector. We calculated the 95% significance level of the sectorial PSM from the respective normative data. We calculated the percentage agreement with group1 (G1), healthy on HRT and within normal perimetric limits, and group 2 (G2), abnormal on HRT and outside normal perimetric limits. We also examined the relationship of PSM and rim area (RA) in those sectors classified as abnormal by MRA (Moorfields Regression Analysis) of HRT. Results: The mean age was 65 (range= [37, 89]). The global sensitivity versus specificity of each instrument in detecting glaucomatous eyes was: MDT 80% vs. 88%, SAP 80% vs. 80%, PP 70% vs. 89% and HRT 80% vs. 79%. Highest percentage agreement of HRT (respectively G1, G2, sector) with PSM were MDT (89%, 57%, nasal superior), SAP (83%, 74%, temporal superior), PP (74%, 63%, nasal superior). Globally percentage agreement (respectively G1, G2) was MDT (92%, 28%), SAP (87%, 40%) and PP (77%, 49%). Linear regression showed there was no significant trend globally associating RA and PSM. However, sectorally the supero-nasal sector had a statistically significant (p<0.001) trend with each instrument, the associated r2 coefficients are (MDT 0.38 SAP 0.56 and PP 0.39). Conclusions: There were no significant differences in global sensitivity or specificity between instruments. Structure-function relationships varied significantly between instruments and were consistently strongest supero-nasally. Further studies are required to investigate these relationships in detail.

Automated auditory mismatch negativity paradigm improves coma prognostic accuracy after cardiac arrest and therapeutic hypothermia.

PURPOSE: EEG and somatosensory evoked potential are highly predictive of poor outcome after cardiac arrest; their accuracy for good recovery is however low. We evaluated whether addition of an automated mismatch negativity-based auditory discrimination paradigm (ADP) to EEG and somatosensory evoked potential improves prediction of awakening. METHODS: EEG and ADP were prospectively recorded in 30 adults during therapeutic hypothermia and in normothermia. We studied the progression of auditory discrimination on single-trial multivariate analyses from therapeutic hypothermia to normothermia, and its correlation to outcome at 3 months, assessed with cerebral performance categories. RESULTS: At 3 months, 18 of 30 patients (60%) survived; 5 had severe neurologic impairment (cerebral performance categories = 3) and 13 had good recovery (cerebral performance categories = 1-2). All 10 subjects showing improvements of auditory discrimination from therapeutic hypothermia to normothermia regained consciousness: ADP was 100% predictive for awakening. The addition of ADP significantly improved mortality prediction (area under the curve, 0.77 for standard model including clinical examination, EEG, somatosensory evoked potential, versus 0.86 after adding ADP, P = 0.02). CONCLUSIONS: This automated ADP significantly improves early coma prognostic accuracy after cardiac arrest and therapeutic hypothermia. The progression of auditory discrimination is strongly predictive of favorable recovery and appears complementary to existing prognosticators of poor outcome. Before routine implementation, validation on larger cohorts is warranted.

Automated quantitative pupillometry for the prognostication of coma after cardiac arrest.

BACKGROUND: Sedation and therapeutic hypothermia (TH) delay neurological responses and might reduce the accuracy of clinical examination to predict outcome after cardiac arrest (CA). We examined the accuracy of quantitative pupillary light reactivity (PLR), using an automated infrared pupillometry, to predict outcome of post-CA coma in comparison to standard PLR, EEG, and somato-sensory evoked potentials (SSEP). METHODS: We prospectively studied over a 1-year period (June 2012-June 2013) 50 consecutive comatose CA patients treated with TH (33 °C, 24 h). Quantitative PLR (expressed as the % of pupillary response to a calibrated light stimulus) and standard PLR were measured at day 1 (TH and sedation; on average 16 h after CA) and day 2 (normothermia, off sedation: on average 46 h after CA). Neurological outcome was assessed at 90 days with Cerebral Performance Categories (CPC), dichotomized as good (CPC 1-2) versus poor (CPC 3-5). Predictive performance was analyzed using area under the ROC curves (AUC). RESULTS: Patients with good outcome [n = 23 (46 %)] had higher quantitative PLR than those with poor outcome [n = 27; 16 (range 9-23) vs. 10 (1-30) % at day 1, and 20 (13-39) vs. 11 (1-55) % at day 2, both p < 0.001]. Best cut-off for outcome prediction of quantitative PLR was <13 %. The AUC to predict poor outcome was higher for quantitative than for standard PLR at both time points (day 1, 0.79 vs. 0.56, p = 0.005; day 2, 0.81 vs. 0.64, p = 0.006). Prognostic accuracy of quantitative PLR was comparable to that of EEG and SSEP (0.81 vs. 0.80 and 0.73, respectively, both p > 0.20). CONCLUSIONS: Quantitative PLR is more accurate than standard PLR in predicting outcome of post-anoxic coma, irrespective of temperature and sedation, and has comparable prognostic accuracy than EEG and SSEP.

Atlas-based segmentation of pathological MR brain images using a model of lesion growth.

We propose a method for brain atlas deformation in the presence of large space-occupying tumors, based on an a priori model of lesion growth that assumes radial expansion of the lesion from its starting point. Our approach involves three steps. First, an affine registration brings the atlas and the patient into global correspondence. Then, the seeding of a synthetic tumor into the brain atlas provides a template for the lesion. The last step is the deformation of the seeded atlas, combining a method derived from optical flow principles and a model of lesion growth. Results show that a good registration is performed and that the method can be applied to automatic segmentation of structures and substructures in brains with gross deformation, with important medical applications in neurosurgery, radiosurgery, and radiotherapy.

Model-based Segmentation and Fusion of 3D Computed Tomography and 3D Ultrasound of the Eye for Radiotherapy Planning

Computed Tomography (CT) represents the standard imaging modality for tumor volume delineation for radiotherapy treatment planning of retinoblastoma despite some inherent limitations. CT scan is very useful in providing information on physical density for dose calculation and morphological volumetric information but presents a low sensitivity in assessing the tumor viability. On the other hand, 3D ultrasound (US) allows a highly accurate definition of the tumor volume thanks to its high spatial resolution but it is not currently integrated in the treatment planning but used only for diagnosis and follow-up. Our ultimate goal is an automatic segmentation of gross tumor volume (GTV) in the 3D US, the segmentation of the organs at risk (OAR) in the CT and the registration of both modalities. In this paper, we present some preliminary results in this direction. We present 3D active contour-based segmentation of the eye ball and the lens in CT images; the presented approach incorporates the prior knowledge of the anatomy by using a 3D geometrical eye model. The automated segmentation results are validated by comparing with manual segmentations. Then, we present two approaches for the fusion of 3D CT and US images: (i) landmark-based transformation, and (ii) object-based transformation that makes use of eye ball contour information on CT and US images.

Automated analysis of background EEG and reactivity during therapeutic hypothermia in comatose patients after cardiac arrest.

Visual analysis of electroencephalography (EEG) background and reactivity during therapeutic hypothermia provides important outcome information, but is time-consuming and not always consistent between reviewers. Automated EEG analysis may help quantify the brain damage. Forty-six comatose patients in therapeutic hypothermia, after cardiac arrest, were included in the study. EEG background was quantified with burst-suppression ratio (BSR) and approximate entropy, both used to monitor anesthesia. Reactivity was detected through change in the power spectrum of signal before and after stimulation. Automatic results obtained almost perfect agreement (discontinuity) to substantial agreement (background reactivity) with a visual score from EEG-certified neurologists. Burst-suppression ratio was more suited to distinguish continuous EEG background from burst-suppression than approximate entropy in this specific population. Automatic EEG background and reactivity measures were significantly related to good and poor outcome. We conclude that quantitative EEG measurements can provide promising information regarding current state of the patient and clinical outcome, but further work is needed before routine application in a clinical setting.

Registration of congenital anomalies in Switzerland by EUROCAT.

Since 1988 the epidemiological surveillance of congenital anomalies (malformations, chromosomal aberrations, metabolic diseases, hereditary diseases, neurosensorial defects, etc.) is carried out by the Swiss registry of EUROCAT (European Registry of Congenital Anomalies and Twins). Several Swiss cantons collaborate through their own local registry, transmitting data to the central registry in Lausanne. We present the main objectives and methods of registration and give the global prevalence rates for the main malformations for 1996 and the period 1993-1996.

Performance of an automated multiplex immunofluorescence assay for detection of Chlamydia trachomatis immunoglobulin G.

Chlamydia serology is indicated to investigate etiology of miscarriage, infertility, pelvic inflammatory disease, and ectopic pregnancy. Here, we assessed the reliability of a new automated-multiplex immunofluorescence assay (InoDiag test) to detect specific anti-C. trachomatis immunoglobulin G. Considering immunofluorescence assay (IF) as gold standard, InoDiag tests exhibited similar sensitivities (65.5%) but better specificities (95.1%-98%) than enzyme-linked immunosorbent assays (ELISAs). InoDiag tests demonstrated similar or lower cross-reactivity rates when compared to ELISA or IF.

Population-based registration of phenotypic anatomic and familial data for melanoma: results from a Swiss multicentric study

CONTEXT AND OBJECTIVES: A multicentric study was set up to assess the feasibility for Swiss cancer registries of actively retrieving 3 additional variables of epidemiological and a etiological relevance for melanoma, and of potential use for the evaluation of prevention campaigns. MATERIAL AND METHODS: The skin type, family history of melanoma and precise anatomical site were retrieved for melanoma cases registered in 5 Swiss cantons (Neuchâtel, St-Gall and Appenzell, Vaud and Wallis) over 3 to 6 consecutive years (1995-2002). Data were obtained via a short questionnaire administered by the physicians - mostly dermatologists - who originally excised the lesions. As the detailed body site was routinely collected in Ticino, data from this Cancer Registry were included in the body site analysis. Relative melanoma density (RMD) was computed by the ratio of observed to expected numbers of melanomas allowing for body site surface areas, and further adjusted for site-specific melanocyte density. RESULTS: Of the 1,645 questionnaires sent, 1,420 (86.3%) were returned. The detailed cutaneous site and skin type were reliably obtained for 84.7% and 78.7% of questionnaires, and family history was known in 76% of instances. Prevalence of sun-sensitive subjects and patients with melanoma affected first-degree relatives, two target groups for early detection and surveillance campaigns were 54.1% and 3.4%, respectively. After translation into the 4th digit of the International Classification of Diseases for Oncology, the anatomical site codes from printed (original information) and pictorial support (body chart from the questionnaire) concurred for 94.6% of lesions. Discrepancies occurred mostly for lesions on the upper, outer part of the shoulder for which the clinician's textual description was "shoulder blade". This differential misclassification suggests under-estimation by about 10% of melanomas of the upper limbs and an over-estimation of 5% for truncal melanomas. Sites of highest melanoma risk were the face, the shoulder and the upper arm for sexes, the back for men and the leg for women. Three major features of this series were: (1) an unexpectedly high RMD for the face in women (6.2 vs 4.2 in men), (2) the absence of a male predominance for melanomas on the ears, and (3) for the upper limbs, a steady gradient of increasing melanoma density with increasing proximity to the trunk, regardless of sex. DISCUSSION AND CONCLUSION: The feasibility of retrieving the skin type, the precise anatomical location and family history of melanoma in a reliable manner was demonstrated thanks to the collaboration of Swiss dermatologists. Use of a schematic body drawing improves the quality of the anatomical site data and facilitate the reporting task of doctors. Age and sex patterns of RMD paralleled general indicators of sun exposure and behaviour, except for the hand (RMD=0.2). These Swiss results support some site or sun exposure specificity in the aetiology of melanoma.