Quality Assessment of Cardiac Magnetic Resonance Images at 1.5/3.0 T: Description and Validation of Standardized Criteria

PURPOSE: Cardiovascular magnetic resonance (CMR) has become a robust and important diagnostic imaging modality in cardiovascular medicine. However,insufficient image quality may compromise its diagnostic accuracy. No standardized criteria are available to assess the quality of CMR studies. We aimed todescribe and validate standardized criteria to evaluate the quality of CMR studies including: a) cine steady-state free precession, b) delayed gadoliniumenhancement, and c) adenosine stress first-pass perfusion. These criteria will serve for the assessment of the image quality in the setting of the Euro-CMR registry.METHOD AND MATERIALS: First, a total of 45 quality criteria were defined (35 qualitative criteria with a score from 0-3, and 10 quantitative criteria). Thequalitative score ranged from 0 to 105. The lower the qualitative score, the better the quality. The quantitative criteria were based on the absolute signal intensity (delayed enhancement) and on the signal increase (perfusion) of the anterior/posterior left ventricular wall after gadolinium injection. These criteria were then applied in 30 patients scanned with a 1.5T system and in 15 patients scanned with a 3.0T system. The examinations were jointly interpreted by 3 CMR experts and 1 study nurse. In these 45 patients the correlation between the results of the quality assessment obtained by the different readers was calculated.RESULTS: On the 1.5T machine, the mean quality score was 3.5. The mean difference between each pair of observers was 0.2 (5.7%) with a mean standarddeviation of 1.4. On the 3.0T machine, the mean quality score was 4.4. The mean difference between each pair of onservers was 0.3 (6.4%) with a meanstandard deviation of 1.6. The quantitative quality assessments between observers were well correlated for the 1.5T machine: R was between 0.78 and 0.99 (pCONCLUSION: The described criteria for the assessment of CMR image quality are robust and have a low inter-observer variability, especially on 1.5T systems.CLINICAL RELEVANCE/APPLICATION: These criteria will allow the standardization of CMR examinations. They will help to improve the overall quality ofexaminations and the comparison between clinical studies.

Validation of the Addiction Severity Index in French-speaking alcoholic patients.

OBJECTIVE: Evaluation of a French translation of the Addiction Severity Index (ASI) in 100 (78 male) alcoholic patients. METHOD: Validity of the instrument was assessed by measuring test-retest and interrater reliability, internal consistency and convergence and discrimination between items and scales. Concurrent validity was assessed by comparing the scores from the ASI with those obtained from three other clinimetric instruments. RESULTS: Test-retest reliability of ASI scores (after a 10-day interval) was good (r = 0.63 to r = 0.95). Interrater reliability was evaluated using six video recordings of patient interviews. Severity ratings assigned by six rates were significantly different (p < .05), but 72% of the ratings assigned by those who viewed the videos were within two points of the interviewer's severity ratings. Cronbach alpha coefficient of internal consistency varied from 0.58 to 0.81 across scales. The average item-to-scale convergent validity (r value) was 0.49 (range 0.0 to 0.84) for composite scores and 0.35 (range 0.00 to 0.68) for severity ratings, whereas discriminant validity was 0.11 on average (range-0.19 to 0.46) for composite scores and 0.12 (range-0.20 to 0.52) for severity ratings. Finally, concurrent validity with the following instruments was assessed: Severity of Alcoholism Dependence Questionnaire (40% shared variance with ASI alcohol scale), Michigan Alcoholism Screening Test (2% shared variance with ASI alcohol scale) and Hamilton Depression Rating Scale (31% shared variance with ASI psychiatric scale). CONCLUSIONS: The Addiction Severity Index covers a large scope of problems encountered among alcoholics and quantifies need for treatment. This French version presents acceptable criteria of reliability and validity.

Development and validation of a postmortem radiological alteration index: the RA-Index.

This study aimed to derive an index quantifying the state of alteration of cadavers by quantifying the presence of gas in the body using postmortem multidetector computed tomography (MDCT) imaging, and to validate the index by defining its sensitivity and specificity. The RA (radiological alteration)-index was derived from postmortem MDCT data from 118 nontraumatically deceased people. To validate the index, 100 additional scanned bodies (50 % traumatically deceased) were retrospectively examined by two independent observers. Presence of gas at 82 sites was assessed by a radiologist, whereas a forensic pathologist only investigated the seven sites used for the RA-index. The RA-index was highly correlated to the overall presence of gas in all 82 sites (R(2) = 0.98 in the derivation set and 0.85 in the validation set). Semiquantitative evaluation of gas presence in each site showed moderate reliability (Cohen's kappa range, 0.41-0.78); nevertheless, the overall RA-index was very reliable (ICC(2,1) = 0.95; 95 % CI 0.92-0.96). Examiner using the RA-index detected heart cavities full of gas with a sensitivity of 100 % (95 % CI 51.7-100) and a specificity of 98.8 % (92.6-99.9). We conclude that determining the presence of gas at seven sites is a valid means to measure the distribution of gas due to cadaveric alteration in the entire body. The RA-index is rapid, easy-to-use, and reliable for nonexperienced users, and it is a valid method to suspect the normal presence of gas from cadaveric alteration. MDCT can be used to screen for gas embolism and to give indications for gas composition analysis (gas chromatography).

Validation of methane measurement using headspace-GC-MS and quantification by a stable isotope-labeled internal standard generated in situ.

A previous study has shown the possibility to identify methane (CH4 ) using headspace-GC-MS and quantify it with a stable isotope as internal standard. The main drawback of the GC-MS methods discussed in literature for CH4 measurement is the absence of a specific internal standard necessary to perform quantification. However, it becomes essential to develop a safer method to limit the manipulation of gaseous CH4 and to precisely control the injected amount of gas for spiking and calibration by comparison with external calibration. To avoid the manipulation of a stable isotope-labeled gas, we have chosen to generate a labeled gas as an internal standard in a vial on the basis of the formation of CH4 by the reaction of Grignard reagent methylmagnesium chloride with deuterated water. This method allows precise measurement of CH4 concentrations in gaseous sample as well as in a solid or a liquid sample after a thermodesorption step in a headspace vial. A full accuracy profile validation of this method is then presented.

Development, design and validation of solid reference samples.

We developed a method of sample preparation using epoxy compound, which was validated in two steps. First, we studied the homogeneity within samples by scanning tubes filled with radioactive epoxy. We found within-sample homogeneity better than 2%. Then, we studied the homogeneity between samples during a 4.5 h dispensing time. The homogeneity between samples was found to be better than 2%. This study demonstrates that we have a validated method, which assures the traceability of epoxy samples.

Identification and validation of copy number variants using SNP genotyping arrays from a large clinical cohort.

BACKGROUND: Genotypes obtained with commercial SNP arrays have been extensively used in many large case-control or population-based cohorts for SNP-based genome-wide association studies for a multitude of traits. Yet, these genotypes capture only a small fraction of the variance of the studied traits. Genomic structural variants (GSV) such as Copy Number Variation (CNV) may account for part of the missing heritability, but their comprehensive detection requires either next-generation arrays or sequencing. Sophisticated algorithms that infer CNVs by combining the intensities from SNP-probes for the two alleles can already be used to extract a partial view of such GSV from existing data sets. RESULTS: Here we present several advances to facilitate the latter approach. First, we introduce a novel CNV detection method based on a Gaussian Mixture Model. Second, we propose a new algorithm, PCA merge, for combining copy-number profiles from many individuals into consensus regions. We applied both our new methods as well as existing ones to data from 5612 individuals from the CoLaus study who were genotyped on Affymetrix 500K arrays. We developed a number of procedures in order to evaluate the performance of the different methods. This includes comparison with previously published CNVs as well as using a replication sample of 239 individuals, genotyped with Illumina 550K arrays. We also established a new evaluation procedure that employs the fact that related individuals are expected to share their CNVs more frequently than randomly selected individuals. The ability to detect both rare and common CNVs provides a valuable resource that will facilitate association studies exploring potential phenotypic associations with CNVs. CONCLUSION: Our new methodologies for CNV detection and their evaluation will help in extracting additional information from the large amount of SNP-genotyping data on various cohorts and use this to explore structural variants and their impact on complex traits.

Photopolymerizable hydrogels for implants: Monte-Carlo modeling and experimental in vitro validation.

Photopolymerization is commonly used in a broad range of bioapplications, such as drug delivery, tissue engineering, and surgical implants, where liquid materials are injected and then hardened by means of illumination to create a solid polymer network. However, photopolymerization using a probe, e.g., needle guiding both the liquid and the curing illumination, has not been thoroughly investigated. We present a Monte Carlo model that takes into account the dynamic absorption and scattering parameters as well as solid-liquid boundaries of the photopolymer to yield the shape and volume of minimally invasively injected, photopolymerized hydrogels. In the first part of the article, our model is validated using a set of well-known poly(ethylene glycol) dimethacrylate hydrogels showing an excellent agreement between simulated and experimental volume-growth-rates. In the second part, in situ experimental results and simulations for photopolymerization in tissue cavities are presented. It was found that a cavity with a volume of 152  mm3 can be photopolymerized from the output of a 0.28-mm2 fiber by adding scattering lipid particles while only a volume of 38  mm3 (25%) was achieved without particles. The proposed model provides a simple and robust method to solve complex photopolymerization problems, where the dimension of the light source is much smaller than the volume of the photopolymerizable hydrogel.

Ruling out coronary heart disease in primary care: external validation of a clinical prediction rule.

BACKGROUND: The Marburg Heart Score (MHS) aims to assist GPs in safely ruling out coronary heart disease (CHD) in patients presenting with chest pain, and to guide management decisions. AIM: To investigate the diagnostic accuracy of the MHS in an independent sample and to evaluate the generalisability to new patients. DESIGN AND SETTING: Cross-sectional diagnostic study with delayed-type reference standard in general practice in Hesse, Germany. METHOD: Fifty-six German GPs recruited 844 males and females aged ≥ 35 years, presenting between July 2009 and February 2010 with chest pain. Baseline data included the items of the MHS. Data on the subsequent course of chest pain, investigations, hospitalisations, and medication were collected over 6 months and were reviewed by an independent expert panel. CHD was the reference condition. Measures of diagnostic accuracy included the area under the receiver operating characteristic curve (AUC), sensitivity, specificity, likelihood ratios, and predictive values. RESULTS: The AUC was 0.84 (95% confidence interval [CI] = 0.80 to 0.88). For a cut-off value of 3, the MHS showed a sensitivity of 89.1% (95% CI = 81.1% to 94.0%), a specificity of 63.5% (95% CI = 60.0% to 66.9%), a positive predictive value of 23.3% (95% CI = 19.2% to 28.0%), and a negative predictive value of 97.9% (95% CI = 96.2% to 98.9%). CONCLUSION: Considering the diagnostic accuracy of the MHS, its generalisability, and ease of application, its use in clinical practice is recommended.

The use of uniaxial accelerometry for the assessment of physical-activity-related energy expenditure: a validation study against whole-body indirect calorimetry.

Assessing the total energy expenditure (TEE) and the levels of physical activity in free-living conditions with non-invasive techniques remains a challenge. The purpose of the present study was to investigate the accuracy of a new uniaxial accelerometer for assessing TEE and physical-activity-related energy expenditure (PAEE) over a 24 h period in a respiratory chamber, and to establish activity levels based on the accelerometry ranges corresponding to the operationally defined metabolic equivalent (MET) categories. In study 1, measurement of the 24 h energy expenditure of seventy-nine Japanese subjects (40 (SD 12) years old) was performed in a large respiratory chamber. During the measurements, the subjects wore a uniaxial accelerometer (Lifecorder; Suzuken Co. Ltd, Nagoya, Japan) on their belt. Two moderate walking exercises of 30 min each were performed on a horizontal treadmill. In study 2, ten male subjects walked at six different speeds and ran at three different speeds on a treadmill for 4 min, with the same accelerometer. O2 consumption was measured during the last minute of each stage and was expressed in MET. The measured TEE was 8447 (SD 1337) kJ/d. The accelerometer significantly underestimated TEE and PAEE (91.9 (SD 5.4) and 92.7 (SD 17.8) % chamber value respectively); however, there was a significant correlation between the two values (r 0.928 and 0.564 respectively; P<0.001). There was a strong correlation between the activity levels and the measured MET while walking (r(2) 0.93; P<0.001). Although TEE and PAEE were systematically underestimated during the 24 h period, the accelerometer assessed energy expenditure well during both the exercise period and the non-structured activities. Individual calibration factors may help to improve the accuracy of TEE estimation, but the average calibration factor for the group is probably sufficient for epidemiological research. This method is also important for assessing the diurnal profile of physical activity.

A comparative study on inter and intralaboratory reproducibility of renin measurement with a conventional enzymatic method and a new chemiluminescent assay of immunoreactive renin.

BACKGROUND: The activity of the renin-angiotensin system is usually evaluated as plasma renin activity (PRA, ngAI/ml per h) but the reproducibility of this enzymatic assay is notoriously scarce. We compared the inter and intralaboratory reproducibilities of PRA with those of a new automated chemiluminescent assay, which allows the direct quantification of immunoreactive renin [chemiluminescent immunoreactive renin (CLIR), microU/ml]. METHODS: Aliquots from six pool plasmas of patients with very low to very high PRA levels were measured in 12 centres with both the enzymatic and the direct assays. The same methods were applied to three control plasma preparations with known renin content. RESULTS: In pool plasmas, mean PRA values ranged from 0.14 +/- 0.08 to 18.9 +/- 4.1 ngAI/ml per h, whereas those of CLIR ranged from 4.2 +/- 1.7 to 436 +/- 47 microU/ml. In control plasmas, mean values of PRA and of CLIR were always within the expected range. Overall, there was a significant correlation between the two methods (r = 0.73, P < 0.01). Similar correlations were found in plasmas subdivided in those with low, intermediate and high PRA. However, the coefficients of variation among laboratories found for PRA were always higher than those of CLIR, ranging from 59.4 to 17.1% for PRA, and from 41.0 to 10.7% for CLIR (P < 0.01). Also, the mean intralaboratory variability was higher for PRA than for CLIR, being respectively, 8.5 and 4.5% (P < 0.01). CONCLUSION: The measurement of renin with the chemiluminescent method is a reliable alternative to PRA, having the advantage of a superior inter and intralaboratory reproducibility.

A Multi-SNP Locus-Association Method Reveals a Substantial Fraction of the Missing Heritability.

There are many known examples of multiple semi-independent associations at individual loci; such associations might arise either because of true allelic heterogeneity or because of imperfect tagging of an unobserved causal variant. This phenomenon is of great importance in monogenic traits but has not yet been systematically investigated and quantified in complex-trait genome-wide association studies (GWASs). Here, we describe a multi-SNP association method that estimates the effect of loci harboring multiple association signals by using GWAS summary statistics. Applying the method to a large anthropometric GWAS meta-analysis (from the Genetic Investigation of Anthropometric Traits consortium study), we show that for height, body mass index (BMI), and waist-to-hip ratio (WHR), 3%, 2%, and 1%, respectively, of additional phenotypic variance can be explained on top of the previously reported 10% (height), 1.5% (BMI), and 1% (WHR). The method also permitted a substantial increase (by up to 50%) in the number of loci that replicate in a discovery-validation design. Specifically, we identified 74 loci at which the multi-SNP, a linear combination of SNPs, explains significantly more variance than does the best individual SNP. A detailed analysis of multi-SNPs shows that most of the additional variability explained is derived from SNPs that are not in linkage disequilibrium with the lead SNP, suggesting a major contribution of allelic heterogeneity to the missing heritability.

Mesurer la stigmatisation perçue chez les personnes souffrant de troubles psychiques : traduction française, validation et adaptation de la Stigma Scale [French translation, validation and adaptation of the Stigma Scale].

AIM: People suffering from mental illness are exposed to stigma. However, only few tools are available to assess stigmatization as perceived from the patient's perspective. The aim of this study is to adapt and validate a French version of the Stigma Scale (King et al., 2007 [8]). This self-report questionnaire has a three-factor structure: discrimination, disclosure and positive aspects of mental illness. Discrimination subscale refers to perceived negative reactions of others. Disclosure subscale refers mainly to managing disclosure to avoid discrimination and finally positive aspects subscale taps into how patients are becoming more accepting, more understanding toward their illness. METHOD: In the first step, internal consistency, convergent validity and test-retest reliability of the French adaptation of the 28-item scale were assessed in a sample of 183 patients. Results of confirmatory factor analyses (CFA) did not confirm the hypothesized structure. In the light of the failed attempts to validate the original version, an alternative 9-item short-form version of the Stigma Scale, maintaining the integrity of the original model, was developed based on results of exploratory factor analyses in the first sample and cross-validated in a new sample of 234 patients. RESULTS: Results of CFA did not confirm that the data fitted well to the three-factor model of the 28-item Stigma Scale (χ(2)/df=2.02, GFI=0.77, AGFI=0.73, RMSEA=0.07, CFI=0.77 and NNFI=0.75). Cronbach's α was excellent for discrimination (0.84) and disclosure (0.83) subscales but poor for potential positive aspects (0.46). External validity was satisfactory. Overall Stigma Scale total score was negatively correlated with the score on Rosenberg's Self-Esteem Scale (r=-0.49), and each subscale was significantly correlated with a visual analogue scale that referred to the specific aspect of stigma (0.43≤|r|≤0.60). Intraclass correlation coefficients between 0.68 and 0.89 indicated good test-retest reliability. The results of the CFA demonstrated that the items chosen for the short version of the Stigma Scale have the expected fit properties (χ(2)/df=1.02, GFI=0.98, AGFI=0.98, RMSEA=0.01, CFI=1.0 and NNFI=1.0). Considering the small number (three) of items in each subscale of the short version of the Stigma Scale, α coefficients for discrimination (0.57), disclosure (0.80) and potential positive aspects subscales (0.62) are considered as good. CONCLUSION: Our results suggest that the 9-item French short version of the Stigma Scale is a useful, reliable and valid self-report questionnaire to assess perceived stigmatization in people suffering from mental illness. The time of completion is really short and questions are well understood and accepted by the patients.

Cross-Cultural Adaptation and Validation of the Spinal Function Sort (SFS) for French- and German-Speaking Patients with Back Complaints

Aim: Functional subjective evaluation through questionnaire is fundamental, but not often realized in patients with back complaints, notably because of lack of validated tools, in accordance with recognized psychometric criteria. The Spinal Function Sort (SFS), developed according to actual standards, was only validated in English. The aim of this study is to translate, adapt and validate the French and German version of the SFS.Method and material: The translation and cross-cultural adaptation were performed following the methodology proposed by the American Association of Orthopedist Surgeon. A total of 344 patients, presenting varied back complaints (especially degenerative and traumatic), took part in this study in a tertiary French- (n=87; mean age 44y; 17 women) and German-speaking (n=257; mean age 41y; 53 women) center. Test-retest reliability was quantified using the intraclass correlation coefficient (ICC) and construct validity was assessed by estimating the Pearson's correlation with the SF-36 physical and mental scales, the Visual Analogue Scale for Pain Intensity (VAS), and subscales of the Hospital Anxiety and Depression Scale (HADS).Results: Respectively for the French and German version, ICC were 0.98 and 0.94. Correlations 0.63 and 0.67 with the SF-36 Physical Functioning subscale; 0.60 and 0.52 with the SF-36 Physical Summary Scale ; -0.33 and -0.51 with the VAS ; -0.08 and 0.25 with the SF-36 Mental Health scale; 0.01 and 0.28 with the SF-36 Mental Summary Scale; -0.26 and -0.42 with the HADS depression; -0.17 and -0.45 with the HADS anxiety.Discussion: For both the French and German version of the SFS, the reliability was excellent. Convergent construct validity with SF-36 physical scales is good, moderated with the VAS. We find out a low correlation with SF-36 mental scales (divergent construct validity). We find out a low correlation with HADS subscales in the French version, and a moderate one in the German version. Selection bias, chronicity of the complaints, as well as cultural differences could explain these results. In conclusion, both the French and German version of the SFS are valid and reliable for evaluation of perceived functional capacity for patients with back complaints.

Simultaneous measurement of regional cerebral blood flow by perfusion CT and stable xenon CT: a validation study.

BACKGROUND AND PURPOSE: Knowledge of cerebral blood flow (CBF) alterations in cases of acute stroke could be valuable in the early management of these cases. Among imaging techniques affording evaluation of cerebral perfusion, perfusion CT studies involve sequential acquisition of cerebral CT sections obtained in an axial mode during the IV administration of iodinated contrast material. They are thus very easy to perform in emergency settings. Perfusion CT values of CBF have proved to be accurate in animals, and perfusion CT affords plausible values in humans. The purpose of this study was to validate perfusion CT studies of CBF by comparison with the results provided by stable xenon CT, which have been reported to be accurate, and to evaluate acquisition and processing modalities of CT data, notably the possible deconvolution methods and the selection of the reference artery. METHODS: Twelve stable xenon CT and perfusion CT cerebral examinations were performed within an interval of a few minutes in patients with various cerebrovascular diseases. CBF maps were obtained from perfusion CT data by deconvolution using singular value decomposition and least mean square methods. The CBF were compared with the stable xenon CT results in multiple regions of interest through linear regression analysis and bilateral t tests for matched variables. RESULTS: Linear regression analysis showed good correlation between perfusion CT and stable xenon CT CBF values (singular value decomposition method: R(2) = 0.79, slope = 0.87; least mean square method: R(2) = 0.67, slope = 0.83). Bilateral t tests for matched variables did not identify a significant difference between the two imaging methods (P >.1). Both deconvolution methods were equivalent (P >.1). The choice of the reference artery is a major concern and has a strong influence on the final perfusion CT CBF map. CONCLUSION: Perfusion CT studies of CBF achieved with adequate acquisition parameters and processing lead to accurate and reliable results.

Développement et validation d'une méthode de calcul indépendant des unités moniteur basée sur les simulations Monte Carlo pour la vérification des traitements IMRT

Résumé : La radiothérapie par modulation d'intensité (IMRT) est une technique de traitement qui utilise des faisceaux dont la fluence de rayonnement est modulée. L'IMRT, largement utilisée dans les pays industrialisés, permet d'atteindre une meilleure homogénéité de la dose à l'intérieur du volume cible et de réduire la dose aux organes à risque. Une méthode usuelle pour réaliser pratiquement la modulation des faisceaux est de sommer de petits faisceaux (segments) qui ont la même incidence. Cette technique est appelée IMRT step-and-shoot. Dans le contexte clinique, il est nécessaire de vérifier les plans de traitement des patients avant la première irradiation. Cette question n'est toujours pas résolue de manière satisfaisante. En effet, un calcul indépendant des unités moniteur (représentatif de la pondération des chaque segment) ne peut pas être réalisé pour les traitements IMRT step-and-shoot, car les poids des segments ne sont pas connus à priori, mais calculés au moment de la planification inverse. Par ailleurs, la vérification des plans de traitement par comparaison avec des mesures prend du temps et ne restitue pas la géométrie exacte du traitement. Dans ce travail, une méthode indépendante de calcul des plans de traitement IMRT step-and-shoot est décrite. Cette méthode est basée sur le code Monte Carlo EGSnrc/BEAMnrc, dont la modélisation de la tête de l'accélérateur linéaire a été validée dans une large gamme de situations. Les segments d'un plan de traitement IMRT sont simulés individuellement dans la géométrie exacte du traitement. Ensuite, les distributions de dose sont converties en dose absorbée dans l'eau par unité moniteur. La dose totale du traitement dans chaque élément de volume du patient (voxel) peut être exprimée comme une équation matricielle linéaire des unités moniteur et de la dose par unité moniteur de chacun des faisceaux. La résolution de cette équation est effectuée par l'inversion d'une matrice à l'aide de l'algorithme dit Non-Negative Least Square fit (NNLS). L'ensemble des voxels contenus dans le volume patient ne pouvant être utilisés dans le calcul pour des raisons de limitations informatiques, plusieurs possibilités de sélection ont été testées. Le meilleur choix consiste à utiliser les voxels contenus dans le Volume Cible de Planification (PTV). La méthode proposée dans ce travail a été testée avec huit cas cliniques représentatifs des traitements habituels de radiothérapie. Les unités moniteur obtenues conduisent à des distributions de dose globale cliniquement équivalentes à celles issues du logiciel de planification des traitements. Ainsi, cette méthode indépendante de calcul des unités moniteur pour l'IMRT step-andshootest validée pour une utilisation clinique. Par analogie, il serait possible d'envisager d'appliquer une méthode similaire pour d'autres modalités de traitement comme par exemple la tomothérapie. Abstract : Intensity Modulated RadioTherapy (IMRT) is a treatment technique that uses modulated beam fluence. IMRT is now widespread in more advanced countries, due to its improvement of dose conformation around target volume, and its ability to lower doses to organs at risk in complex clinical cases. One way to carry out beam modulation is to sum smaller beams (beamlets) with the same incidence. This technique is called step-and-shoot IMRT. In a clinical context, it is necessary to verify treatment plans before the first irradiation. IMRT Plan verification is still an issue for this technique. Independent monitor unit calculation (representative of the weight of each beamlet) can indeed not be performed for IMRT step-and-shoot, because beamlet weights are not known a priori, but calculated by inverse planning. Besides, treatment plan verification by comparison with measured data is time consuming and performed in a simple geometry, usually in a cubic water phantom with all machine angles set to zero. In this work, an independent method for monitor unit calculation for step-and-shoot IMRT is described. This method is based on the Monte Carlo code EGSnrc/BEAMnrc. The Monte Carlo model of the head of the linear accelerator is validated by comparison of simulated and measured dose distributions in a large range of situations. The beamlets of an IMRT treatment plan are calculated individually by Monte Carlo, in the exact geometry of the treatment. Then, the dose distributions of the beamlets are converted in absorbed dose to water per monitor unit. The dose of the whole treatment in each volume element (voxel) can be expressed through a linear matrix equation of the monitor units and dose per monitor unit of every beamlets. This equation is solved by a Non-Negative Least Sqvare fif algorithm (NNLS). However, not every voxels inside the patient volume can be used in order to solve this equation, because of computer limitations. Several ways of voxel selection have been tested and the best choice consists in using voxels inside the Planning Target Volume (PTV). The method presented in this work was tested with eight clinical cases, which were representative of usual radiotherapy treatments. The monitor units obtained lead to clinically equivalent global dose distributions. Thus, this independent monitor unit calculation method for step-and-shoot IMRT is validated and can therefore be used in a clinical routine. It would be possible to consider applying a similar method for other treatment modalities, such as for instance tomotherapy or volumetric modulated arc therapy.