Helical tomotherapy (HT) for the treatment of anal canal cancer: preliminary clinical results, and dosimetric comparison between HT and intensity-modulated or 3D conformal radiotherapy

Background: To report a single-center experience in 19 patients (pts) with anal canal cancer treated with helical tomotherapy (HT) and concurrent chemotherapy, and compare the dosimetric results with fixed-field intensitymodulated radiotherapy (IMRT) and 3D conformal radiotherapy (3D RT). Materials and Methods: Between 2007 and 2008, 19 consecutive pts were treated with HT and concurrent CT for anal canal cancer. Median age was 59 years (range, 38−83), and female/male ratio was 14/5. The majority of the pts had T2 or T3 tumours (68.4%), and 52.6% had positive lymph nodes. In all 19 pts, pelvic and inguinal nodes, and tumour irradiation was given using HT upto a median dose of 36 Gy (1.8 Gy/fr) followed by a 1-week gap. A boost dose of 23.4 Gy (1.8 Gy/fr) was delivered to the tumour and involved nodes using 3DRT (n = 12), HT (n = 6), or IMRT (n = 1). Simultaneous integrated boost was used in none of the pts. All but one patient with a T1N0 tumour received concomitant mitomycin/5- fluorouracil (n = 12) or mitomycin/capecitabin (n = 7) CT. Toxicity was scored according to the Common Terminology Criteria for Adverse Events (NCICTCAE v3.0). HT plans and treatments were generated using Tomotherapy, Inc., software and hardware; and 3D or IMRT boost plans with the CMS treatment planning system (TPS), using 6−18 MV photons from a Siemens Primus accelerator. For dosimetric comparison, computed tomography data sets of 10 pts were imported into the TPS, and 3D and 5-field step-andshoot IMRT plans were generated for each case. Plans were optimized with the aim of assessing organs at risk (OAR) and healthy-tissue sparing while enforcing highly conformal target coverage, and evaluated by dose-volume histograms (DVH) of planning target volumes (PTV) and OAR. Results: With a median follow-up of 13 months (range, 3−18), all pts are alive and well; except one patient developing local recurrence at 12 months. No patient developed grade 3 or more acute toxicity. No unplanned treatment interruption was necessary because of toxicity. With 360-degree-of-freedom beam projection, HT showed an advantage over 3D or IMRT plans in terms of dose conformity around the PTV, and dose gradients were steeper outside the PTV, resulting in reduced doses to OARs. Using HT, acute toxicity was acceptable, and seemed to be better than historical standards. Conclusion: We conclude that HT combined with concurrent chemotherapy for anal canal cancer is effective and tolerable. Compared to 3DRT or 5-field IMRT, there is better conformity around the PTV, and OAR sparing.

Respiratory self-navigated postcontrast whole-heart coronary MR angiography: initial experience in patients.

PURPOSE: To assess the diagnostic performance of respiratory self-navigation for whole-heart coronary magnetic resonance (MR) angiography in a patient cohort referred for diagnostic cardiac MR imaging. MATERIALS AND METHODS: Written informed consent was obtained from all participants for this institutional review board-approved study. Self-navigated coronary MR angiography was performed after administration of a contrast agent in 78 patients (mean age, 48.5 years ± 20.7 [standard deviation]; 53 male patients) referred for cardiac MR imaging because of coronary artery disease (n = 40), cardiomyopathy (n = 14), congenital anomaly (n = 17), or "other" (n = 7). Examination duration was recorded, and the image quality for each coronary segment was assessed with consensus reading. Vessel sharpness, length, and diameter were measured. Quantitative values in proximal, middle, and distal segments were compared by using analysis of variance and t tests. A double-blinded comparison with the results of x-ray angiography was performed when such results were available. RESULTS: When patients with different indications for cardiac MR imaging were examined with self-navigated postcontrast coronary MR angiography, whole-heart data sets with 1.15-mm isotropic spatial resolution were acquired in an average of 7.38 minutes ± 1.85. The main and proximal coronary segments could be visualized in 92.3% of cases, while the middle and distal segments could be visualized in 84.0% and 55.8% of cases, respectively. Subjective scores and vessel sharpness were significantly higher in the proximal segments than in the middle and distal segments (P < .05). Anomalies of the coronary arteries could be confirmed or excluded in all cases. Per-vessel sensitivity and specificity for stenosis detection were 64.7% and 85.0%, respectively, in the 31 patients for whom reference standard x-ray coronary angiography results were available. CONCLUSION: The self-navigated coronary MR angiography sequence shows promise for coronary imaging. However, technical improvements are needed to improve image quality, especially in the more distal coronary segments.

Increased expression of urokinase-type plasminogen activator mRNA determines adverse prognosis in ErbB2-positive primary breast cancer.

PURPOSE: To evaluate and validate mRNA expression markers capable of identifying patients with ErbB2-positive breast cancer associated with distant metastasis and reduced survival. PATIENTS AND METHODS: Expression of 60 genes involved in breast cancer biology was assessed by quantitative real-time PCR (qrt-PCR) in 317 primary breast cancer patients and correlated with clinical outcome data. Results were validated subsequently using two previously published and publicly available microarray data sets with different patient populations comprising 295 and 286 breast cancer samples, respectively. RESULTS: Of the 60 genes measured by qrt-PCR, urokinase-type plasminogen activator (uPA or PLAU) mRNA expression was the most significant marker associated with distant metastasis-free survival (MFS) by univariate Cox analysis in patients with ErbB2-positive tumors and an independent factor in multivariate analysis. Subsequent validation in two microarray data sets confirmed the prognostic value of uPA in ErbB2-positive tumors by both univariate and multivariate analysis. uPA mRNA expression was not significantly associated with MFS in ErbB2-negative tumors. Kaplan-Meier analysis showed in all three study populations that patients with ErbB2-positive/uPA-positive tumors exhibited significantly reduced MFS (hazard ratios [HR], 4.3; 95% CI, 1.6 to 11.8; HR, 2.7; 95% CI, 1.2 to 6.2; and, HR, 2.8; 95% CI, 1.1 to 7.1; all P < .02) as compared with the group with ErbB2-positive/uPA-negative tumors who exhibited similar outcome to those with ErbB2-negative tumors, irrespective of uPA status. CONCLUSION: After evaluation of 898 breast cancer patients, uPA mRNA expression emerged as a powerful prognostic indicator in ErbB2-positive tumors. These results were consistent among three independent study populations assayed by different techniques, including qrt-PCR and two microarray platforms.

News from compositions, the R package

The R-package “compositions”is a tool for advanced compositional analysis. Its basicfunctionality has seen some conceptual improvement, containing now some facilitiesto work with and represent ilr bases built from balances, and an elaborated subsys-tem for dealing with several kinds of irregular data: (rounded or structural) zeroes,incomplete observations and outliers. The general approach to these irregularities isbased on subcompositions: for an irregular datum, one can distinguish a “regular” sub-composition (where all parts are actually observed and the datum behaves typically)and a “problematic” subcomposition (with those unobserved, zero or rounded parts, orelse where the datum shows an erratic or atypical behaviour). Systematic classificationschemes are proposed for both outliers and missing values (including zeros) focusing onthe nature of irregularities in the datum subcomposition(s).To compute statistics with values missing at random and structural zeros, a projectionapproach is implemented: a given datum contributes to the estimation of the desiredparameters only on the subcompositon where it was observed. For data sets withvalues below the detection limit, two different approaches are provided: the well-knownimputation technique, and also the projection approach.To compute statistics in the presence of outliers, robust statistics are adapted to thecharacteristics of compositional data, based on the minimum covariance determinantapproach. The outlier classification is based on four different models of outlier occur-rence and Monte-Carlo-based tests for their characterization. Furthermore the packageprovides special plots helping to understand the nature of outliers in the dataset.Keywords: coda-dendrogram, lost values, MAR, missing data, MCD estimator,robustness, rounded zeros

Predicting spatial patterns of plant biodiversity : from species to communities

Understanding the distribution and composition of species assemblages and being able to predict them in space and time are highly important tasks io investigate the fate of biodiversity in the current global changes context. Species distribution models are tools that have proven useful to predict the potential distribution of species by relating their occurrences to environmental variables. Species assemblages can then be predicted by combining the prediction of individual species models. In the first part of my thesis, I tested the importance of new environmental predictors to improve species distribution prediction. I showed that edaphic variables, above all soil pH and nitrogen content could be important in species distribution models. In a second chapter, I tested the influence of different resolution of predictors on the predictive ability of species distribution models. I showed that fine resolution predictors could ameliorate the models for some species by giving a better estimation of the micro-topographic condition that species tolerate, but that fine resolution predictors for climatic factors still need to be ameliorated. The second goal of my thesis was to test the ability of empirical models to predict species assemblages' characteristics such as species richness or functional attributes. I showed that species richness could be modelled efficiently and that the resulting prediction gave a more realistic estimate of the number of species than when obtaining it by stacking outputs of single species distribution models. Regarding the prediction of functional characteristics (plant height, leaf surface, seed mass) of plant assemblages, mean and extreme values of functional traits were better predictable than indices reflecting the diversity of traits in the community. This approach proved interesting to understand which environmental conditions influence particular aspects of the vegetation functioning. It could also be useful to predict climate change impacts on the vegetation. In the last part of my thesis, I studied the capacity of stacked species distribution models to predict the plant assemblages. I showed that this method tended to over-predict the number of species and that the composition of the community was not predicted exactly either. Finally, I combined the results of macro- ecological models obtained in the preceding chapters with stacked species distribution models and showed that this approach reduced significantly the number of species predicted and that the prediction of the composition is also ameliorated in some cases. These results showed that this method is promising. It needs now to be tested on further data sets. - Comprendre la manière dont les plantes se répartissent dans l'environnement et s'organisent en communauté est une question primordiale dans le contexte actuel de changements globaux. Cette connaissance peut nous aider à sauvegarder la diversité des espèces et les écosystèmes. Des méthodes statistiques nous permettent de prédire la distribution des espèces de plantes dans l'espace géographique et dans le temps. Ces modèles de distribution d'espèces, relient les occurrences d'une espèce avec des variables environnementales pour décrire sa distribution potentielle. Cette méthode a fait ses preuves pour ce qui est de la prédiction d'espèces individuelles. Plus récemment plusieurs tentatives de cumul de modèles d'espèces individuelles ont été réalisées afin de prédire la composition des communautés végétales. Le premier objectif de mon travail est d'améliorer les modèles de distribution en testant l'importance de nouvelles variables prédictives. Parmi différentes variables édaphiques, le pH et la teneur en azote du sol se sont avérés des facteurs non négligeables pour prédire la distribution des plantes. Je démontre aussi dans un second chapitre que les prédicteurs environnementaux à fine résolution permettent de refléter les conditions micro-topographiques subies par les plantes mais qu'ils doivent encore être améliorés avant de pouvoir être employés de manière efficace dans les modèles. Le deuxième objectif de ce travail consistait à étudier le développement de modèles prédictifs pour des attributs des communautés végétales tels que, par exemple, la richesse en espèces rencontrée à chaque point. Je démontre qu'il est possible de prédire par ce biais des valeurs de richesse spécifiques plus réalistes qu'en sommant les prédictions obtenues précédemment pour des espèces individuelles. J'ai également prédit dans l'espace et dans le temps des caractéristiques de la végétation telles que sa hauteur moyenne, minimale et maximale. Cette approche peut être utile pour comprendre quels facteurs environnementaux promeuvent différents types de végétation ainsi que pour évaluer les changements à attendre au niveau de la végétation dans le futur sous différents régimes de changements climatiques. Dans une troisième partie de ma thèse, j'ai exploré la possibilité de prédire les assemblages de plantes premièrement en cumulant les prédictions obtenues à partir de modèles individuels pour chaque espèce. Cette méthode a le défaut de prédire trop d'espèces par rapport à ce qui est observé en réalité. J'ai finalement employé le modèle de richesse en espèce développé précédemment pour contraindre les résultats du modèle d'assemblage de plantes. Cela a permis l'amélioration des modèles en réduisant la sur-prédiction et en améliorant la prédiction de la composition en espèces. Cette méthode semble prometteuse mais de nouveaux tests sont nécessaires pour bien évaluer ses capacités.

New tools for epidemiology: a space odyssey

Geographical Information Systems (GIS) facilitate access to epidemiological data through visualization and may be consulted for the development of mathematical models and analysis by spatial statistics. Variables such as land-cover, land-use, elevations, surface temperatures, rainfall etc. emanating from earth-observing satellites, complement GIS as this information allows the analysis of disease distribution based on environmental characteristics. The strength of this approach issues from the specific environmental requirements of those causative infectious agents, which depend on intermediate hosts for their transmission. The distribution of these diseases is restricted, both by the environmental requirements of their intermediate hosts/vectors and by the ambient temperature inside these hosts, which effectively govern the speed of maturation of the parasite. This paper discusses the current capabilities with regard to satellite data collection in terms of resolution (spatial, temporal and spectral) of the sensor instruments on board drawing attention to the utility of computer-based models of the Earth for epidemiological research. Virtual globes, available from Google and other commercial firms, are superior to conventional maps as they do not only show geographical and man-made features, but also allow instant import of data-sets of specific interest, e.g. environmental parameters, demographic information etc., from the Internet.

Leprosy frequency in the world, 1999-2010

In 1991, the World Health Organization (WHO) committed to reducing the prevalence of leprosy to below 1 in 10,000 inhabitants by 2000. Significant improvements in leprosy control have occurred, but leprosy remains a public health problem in many countries due to its high incidence and rate of transmission. This paper reviews data published by the WHO in the years 2000, 2005 and 2010. These data sets included 148 countries or territories that reported to the WHO at least once. Only four countries reported higher prevalence rates in 2010 than in 2000 and eight reported higher case detection rate (CDR) in 2009 than in 1999. Prevalence rate reductions were greater for the first five-year period examined, while CDR reductions were greater in the second five-year period. Thirty-six countries and territories reported at least one prevalence value higher than 1 per 10,000 inhabitants and 32 reported at least one CDR value higher than 9 per 100,000 inhabitants. A total of 39 countries fit at least one of these criteria and all were located in tropical regions.

Synchrotron-based X-ray tomographic microscopy for rock physics investigations

Synchrotron radiation X-ray tomographic microscopy is a nondestructive method providing ultra-high-resolution 3D digital images of rock microstructures. We describe this method and, to demonstrate its wide applicability, we present 3D images of very different rock types: Berea sandstone, Fontainebleau sandstone, dolomite, calcitic dolomite, and three-phase magmatic glasses. For some samples, full and partial saturation scenarios are considered using oil, water, and air. The rock images precisely reveal the 3D rock microstructure, the pore space morphology, and the interfaces between fluids saturating the same pore. We provide the raw image data sets as online supplementary material, along with laboratory data describing the rock properties. By making these data sets available to other research groups, we aim to stimulate work based on digital rock images of high quality and high resolution. We also discuss and suggest possible applications and research directions that can be pursued on the basis of our data.

Measuring and Explaining the Increase of Travel Distance : A Multilevel Analysis Using Repeated Cross Sectional Travel Surveys

We want to shed some light on the development of person mobility by analysing the repeated cross-sectional data of the four National Travel Surveys (NTS) that were conducted in Germany since the mid seventies. The above mentioned driving forces operate on different levels of the system that generates the spatial behaviour we observe: Travel demand is derived from the needs and desires of individuals to participate in spatially separated activities. Individuals organise their lives in an interactive process within the context they live in, using given infrastructure. Essential determinants of their demand are the individual's socio-demographic characteristics, but also the opportunities and constraints defined by the household and the environment are relevant for the behaviour which ultimately can be realised. In order to fully capture the context which determines individual behaviour, the (nested) hierarchy of persons within households within spatial settings has to be considered. The data we will use for our analysis contains information on these three levels. With the analysis of this micro-data we attempt to improve our understanding of the afore subsumed macro developments. In addition we will investigate the prediction power of a few classic sociodemographic variables for the daily travel distance of individuals in the four NTS data sets, with a focus on the evolution of this predictive power. The additional task to correctly measure distances travelled by means of the NTS is threatened by the fact that although these surveys measure the same variables, different sampling designs and data collection procedures were used. So the aim of the analysis is also to detect variables whose control corrects for the known measurement error, as a prerequisite to apply appropriate models in order to better understand the development of individual travel behaviour in a multilevel context. This task is complicated by the fact that variables that inform on survey procedures and outcomes are only provided with the data set for 2002 (see Infas and DIW Berlin, 2003).

Nonparametric estimation of Value-at-Risk

A method to estimate an extreme quantile that requires no distributional assumptions is presented. The approach is based on transformed kernel estimation of the cumulative distribution function (cdf). The proposed method consists of a double transformation kernel estimation. We derive optimal bandwidth selection methods that have a direct expression for the smoothing parameter. The bandwidth can accommodate to the given quantile level. The procedure is useful for large data sets and improves quantile estimation compared to other methods in heavy tailed distributions. Implementation is straightforward and R programs are available.

alr approach for replacing values below the detection limit

All of the imputation techniques usually applied for replacing values below thedetection limit in compositional data sets have adverse effects on the variability. In thiswork we propose a modification of the EM algorithm that is applied using the additivelog-ratio transformation. This new strategy is applied to a compositional data set and theresults are compared with the usual imputation techniques

The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models.

Gene expression data from microarrays are being applied to predict preclinical and clinical endpoints, but the reliability of these predictions has not been established. In the MAQC-II project, 36 independent teams analyzed six microarray data sets to generate predictive models for classifying a sample with respect to one of 13 endpoints indicative of lung or liver toxicity in rodents, or of breast cancer, multiple myeloma or neuroblastoma in humans. In total, >30,000 models were built using many combinations of analytical methods. The teams generated predictive models without knowing the biological meaning of some of the endpoints and, to mimic clinical reality, tested the models on data that had not been used for training. We found that model performance depended largely on the endpoint and team proficiency and that different approaches generated models of similar performance. The conclusions and recommendations from MAQC-II should be useful for regulatory agencies, study committees and independent investigators that evaluate methods for global gene expression analysis.

Immunoregulation in human malaria: the challenge of understanding asymptomatic infection

Asymptomatic Plasmodium infection carriers represent a major threat to malaria control worldwide as they are silent natural reservoirs and do not seek medical care. There are no standard criteria for asymptomaticPlasmodium infection; therefore, its diagnosis relies on the presence of the parasite during a specific period of symptomless infection. The antiparasitic immune response can result in reducedPlasmodium sp. load with control of disease manifestations, which leads to asymptomatic infection. Both the innate and adaptive immune responses seem to play major roles in asymptomatic Plasmodiuminfection; T regulatory cell activity (through the production of interleukin-10 and transforming growth factor-β) and B-cells (with a broad antibody response) both play prominent roles. Furthermore, molecules involved in the haem detoxification pathway (such as haptoglobin and haeme oxygenase-1) and iron metabolism (ferritin and activated c-Jun N-terminal kinase) have emerged in recent years as potential biomarkers and thus are helping to unravel the immune response underlying asymptomatic Plasmodium infection. The acquisition of large data sets and the use of robust statistical tools, including network analysis, associated with well-designed malaria studies will likely help elucidate the immune mechanisms responsible for asymptomatic infection.

Diffusion spectrum magnetic resonance imaging (DSI) tractography of crossing fibers.

MRI tractography is the mapping of neural fiber pathways based on diffusion MRI of tissue diffusion anisotropy. Tractography based on diffusion tensor imaging (DTI) cannot directly image multiple fiber orientations within a single voxel. To address this limitation, diffusion spectrum MRI (DSI) and related methods were developed to image complex distributions of intravoxel fiber orientation. Here we demonstrate that tractography based on DSI has the capacity to image crossing fibers in neural tissue. DSI was performed in formalin-fixed brains of adult macaque and in the brains of healthy human subjects. Fiber tract solutions were constructed by a streamline procedure, following directions of maximum diffusion at every point, and analyzed in an interactive visualization environment (TrackVis). We report that DSI tractography accurately shows the known anatomic fiber crossings in optic chiasm, centrum semiovale, and brainstem; fiber intersections in gray matter, including cerebellar folia and the caudate nucleus; and radial fiber architecture in cerebral cortex. In contrast, none of these examples of fiber crossing and complex structure was identified by DTI analysis of the same data sets. These findings indicate that DSI tractography is able to image crossing fibers in neural tissue, an essential step toward non-invasive imaging of connectional neuroanatomy.

Arqueología de la producción agrícola. Aplicación de técnicas de análisis estadístico y geocomputacional a un caso de estudio del IV-II milenio ANE en la Cuenca del Río Grande (Perú)

Muchas investigaciones arqueobotánicas, desde un enfoque cualitativo-descriptivo, limitan el propio campo de estudio a los análisis de presencia/ausencia y/o de frecuencia de taxones a partir de su recuento en los conjuntos vegetales. De esa manera, los datos proporcionados resultan ser inconcluyentes y no fiables para la reconstrucción del paleoambiente, la determinación de la dieta alimenticia y de la práctica económica realizada (recolección VS agricultura), y totalmente insuficientes para determinar los cambios históricos ocurridos en los procesos productivos. Por lo que concierne el Perú, desde los primeros estudios con referencia a restos vegetales recuperados en yacimientos arqueológicos, principalmente de la costa, se documenta el importante papel que han desarrollado las especies vegetales en la vida de las comunidades pre-hispánicas. No obstante la excepcional abundancia y óptima preservación de este tipo de material (botánico) en muchos de los yacimientos arqueológicos de esta región, gracias a las extremas condiciones climáticas y ambientales sobre todo de sus áridas zonas costeras, los estudios arqueobotánicos desarrollados hasta el momento son muy escasos y las limitaciones análiticas que presentan en su mayoría reflejan la poca importancia dada a las investigaciones arqueobotánicas. En el presente trabajo desarrollamos y aplicamos una metodología analítica de tipo cuantitativa para el estudio de los macrorestos vegetales procedentes de un yacimiento de la Costa sur del Perú. Con ello pretendemos obtener datos representativos y objetivos de los conjuntos analizados, cuyo procesado lleve a una exhaustiva y correcta interpretación de la información.