O tratamento de metástases ósseas de fração única vs múltiplas frações com 3D CRT ou SBRT: artigo de revisão sistemática

Introdução – A escolha do tratamento depende de vários fatores, incluindo o estado clínico e prognóstico de cada doente. Estes fatores desempenham um papel importante na escolha da intervenção terapêutica em metástases ósseas. A deteção precoce e o tratamento adequado podem melhorar a qualidade de vida e independência funcional dos doentes. Metodologia – Este artigo pretende realizar uma revisão sistemática da literatura dos últimos 15 anos, identificando os diferentes tipos de fracionamentos (fração única versus múltiplas frações) e técnicas utilizadas em radioterapia no tratamento de metástases ósseas. Resultados – Os recentes avanços na tecnologia e nas técnicas de tratamento de radioterapia ajudam na distribuição de doses altamente conformacionais e com orientação por imagem para uma entrega mais precisa do tratamento. A radioterapia estereotáxica corporal (SBRT, do acrónimo inglês stereotactic body radiotherapy) permite delimitar e aumentar a dose nos tumores a irradiar. No caso das metástases ósseas, os resultados de controlo local do tumor e da dor têm-se revelado promissores. A radioterapia convencional de 8Gyx1, no entanto, continua a ser o tratamento mais indicado nos doentes paliativos. Conclusão – O tratamento de metástases ósseas é complexo e uma abordagem multidisciplinar é sempre necessária. O tratamento deve ser individualizado para se adequar aos sintomas e estado clínico de cada doente.

Radiotherapy Treatment Assessment using DCE-MRI

Abstract

The goal of modern radiotherapy is to precisely deliver a prescribed radiation dose to delineated target volumes that contain a significant amount of tumor cells while sparing the surrounding healthy tissues/organs. Precise delineation of treatment and avoidance volumes is the key for the precision radiation therapy. In recent years, considerable clinical and research efforts have been devoted to integrate MRI into radiotherapy workflow motivated by the superior soft tissue contrast and functional imaging possibility. Dynamic contrast-enhanced MRI (DCE-MRI) is a noninvasive technique that measures properties of tissue microvasculature. Its sensitivity to radiation-induced vascular pharmacokinetic (PK) changes has been preliminary demonstrated. In spite of its great potential, two major challenges have limited DCE-MRI’s clinical application in radiotherapy assessment: the technical limitations of accurate DCE-MRI imaging implementation and the need of novel DCE-MRI data analysis methods for richer functional heterogeneity information.

This study aims at improving current DCE-MRI techniques and developing new DCE-MRI analysis methods for particular radiotherapy assessment. Thus, the study is naturally divided into two parts. The first part focuses on DCE-MRI temporal resolution as one of the key DCE-MRI technical factors, and some improvements regarding DCE-MRI temporal resolution are proposed; the second part explores the potential value of image heterogeneity analysis and multiple PK model combination for therapeutic response assessment, and several novel DCE-MRI data analysis methods are developed.

I. Improvement of DCE-MRI temporal resolution. First, the feasibility of improving DCE-MRI temporal resolution via image undersampling was studied. Specifically, a novel MR image iterative reconstruction algorithm was studied for DCE-MRI reconstruction. This algorithm was built on the recently developed compress sensing (CS) theory. By utilizing a limited k-space acquisition with shorter imaging time, images can be reconstructed in an iterative fashion under the regularization of a newly proposed total generalized variation (TGV) penalty term. In the retrospective study of brain radiosurgery patient DCE-MRI scans under IRB-approval, the clinically obtained image data was selected as reference data, and the simulated accelerated k-space acquisition was generated via undersampling the reference image full k-space with designed sampling grids. Two undersampling strategies were proposed: 1) a radial multi-ray grid with a special angular distribution was adopted to sample each slice of the full k-space; 2) a Cartesian random sampling grid series with spatiotemporal constraints from adjacent frames was adopted to sample the dynamic k-space series at a slice location. Two sets of PK parameters’ maps were generated from the undersampled data and from the fully-sampled data, respectively. Multiple quantitative measurements and statistical studies were performed to evaluate the accuracy of PK maps generated from the undersampled data in reference to the PK maps generated from the fully-sampled data. Results showed that at a simulated acceleration factor of four, PK maps could be faithfully calculated from the DCE images that were reconstructed using undersampled data, and no statistically significant differences were found between the regional PK mean values from undersampled and fully-sampled data sets. DCE-MRI acceleration using the investigated image reconstruction method has been suggested as feasible and promising.

Second, for high temporal resolution DCE-MRI, a new PK model fitting method was developed to solve PK parameters for better calculation accuracy and efficiency. This method is based on a derivative-based deformation of the commonly used Tofts PK model, which is presented as an integrative expression. This method also includes an advanced Kolmogorov-Zurbenko (KZ) filter to remove the potential noise effect in data and solve the PK parameter as a linear problem in matrix format. In the computer simulation study, PK parameters representing typical intracranial values were selected as references to simulated DCE-MRI data for different temporal resolution and different data noise level. Results showed that at both high temporal resolutions (<1s) and clinically feasible temporal resolution (~5s), this new method was able to calculate PK parameters more accurate than the current calculation methods at clinically relevant noise levels; at high temporal resolutions, the calculation efficiency of this new method was superior to current methods in an order of 102. In a retrospective of clinical brain DCE-MRI scans, the PK maps derived from the proposed method were comparable with the results from current methods. Based on these results, it can be concluded that this new method can be used for accurate and efficient PK model fitting for high temporal resolution DCE-MRI.

II. Development of DCE-MRI analysis methods for therapeutic response assessment. This part aims at methodology developments in two approaches. The first one is to develop model-free analysis method for DCE-MRI functional heterogeneity evaluation. This approach is inspired by the rationale that radiotherapy-induced functional change could be heterogeneous across the treatment area. The first effort was spent on a translational investigation of classic fractal dimension theory for DCE-MRI therapeutic response assessment. In a small-animal anti-angiogenesis drug therapy experiment, the randomly assigned treatment/control groups received multiple fraction treatments with one pre-treatment and multiple post-treatment high spatiotemporal DCE-MRI scans. In the post-treatment scan two weeks after the start, the investigated Rényi dimensions of the classic PK rate constant map demonstrated significant differences between the treatment and the control groups; when Rényi dimensions were adopted for treatment/control group classification, the achieved accuracy was higher than the accuracy from using conventional PK parameter statistics. Following this pilot work, two novel texture analysis methods were proposed. First, a new technique called Gray Level Local Power Matrix (GLLPM) was developed. It intends to solve the lack of temporal information and poor calculation efficiency of the commonly used Gray Level Co-Occurrence Matrix (GLCOM) techniques. In the same small animal experiment, the dynamic curves of Haralick texture features derived from the GLLPM had an overall better performance than the corresponding curves derived from current GLCOM techniques in treatment/control separation and classification. The second developed method is dynamic Fractal Signature Dissimilarity (FSD) analysis. Inspired by the classic fractal dimension theory, this method measures the dynamics of tumor heterogeneity during the contrast agent uptake in a quantitative fashion on DCE images. In the small animal experiment mentioned before, the selected parameters from dynamic FSD analysis showed significant differences between treatment/control groups as early as after 1 treatment fraction; in contrast, metrics from conventional PK analysis showed significant differences only after 3 treatment fractions. When using dynamic FSD parameters, the treatment/control group classification after 1st treatment fraction was improved than using conventional PK statistics. These results suggest the promising application of this novel method for capturing early therapeutic response.

The second approach of developing novel DCE-MRI methods is to combine PK information from multiple PK models. Currently, the classic Tofts model or its alternative version has been widely adopted for DCE-MRI analysis as a gold-standard approach for therapeutic response assessment. Previously, a shutter-speed (SS) model was proposed to incorporate transcytolemmal water exchange effect into contrast agent concentration quantification. In spite of richer biological assumption, its application in therapeutic response assessment is limited. It might be intriguing to combine the information from the SS model and from the classic Tofts model to explore potential new biological information for treatment assessment. The feasibility of this idea was investigated in the same small animal experiment. The SS model was compared against the Tofts model for therapeutic response assessment using PK parameter regional mean value comparison. Based on the modeled transcytolemmal water exchange rate, a biological subvolume was proposed and was automatically identified using histogram analysis. Within the biological subvolume, the PK rate constant derived from the SS model were proved to be superior to the one from Tofts model in treatment/control separation and classification. Furthermore, novel biomarkers were designed to integrate PK rate constants from these two models. When being evaluated in the biological subvolume, this biomarker was able to reflect significant treatment/control difference in both post-treatment evaluation. These results confirm the potential value of SS model as well as its combination with Tofts model for therapeutic response assessment.

In summary, this study addressed two problems of DCE-MRI application in radiotherapy assessment. In the first part, a method of accelerating DCE-MRI acquisition for better temporal resolution was investigated, and a novel PK model fitting algorithm was proposed for high temporal resolution DCE-MRI. In the second part, two model-free texture analysis methods and a multiple-model analysis method were developed for DCE-MRI therapeutic response assessment. The presented works could benefit the future DCE-MRI routine clinical application in radiotherapy assessment.

Avaliação da eficácia da radioterapia paliativa no controlo da dor em doentes com metástases ósseas

A Radioterapia é um dos tratamentos disponíveis para tratar o cancro seja com intuito curativo ou paliativo. O objectivo da radioterapia paliativa consiste em controlar os sintomas apresentados pelos doentes com doença avançada, como por exemplo a dor provocada pela presença de metástases ósseas, que é o sintoma em estudo neste trabalho. A Radioterapia deve causar o mínimo de efeitos secundários e o mínimo desconforto possível aos doentes, não devendo aumentar a morbilidade, mas antes proporcionar maior qualidade de vida aos doentes. Recentemente, vários estudos investigaram a utilização da radioterapia em fim de vida, proporcionando uma linha orientadora nesta área, a partir das condições clínicas dos doentes submetidos a radioterapia perto do final da vida. A escolha do tema deste trabalho resulta de uma inquietação sentida no dia-a-dia, aquando da realização dos tratamentos de radioterapia a doentes paliativos com dor óssea e com mau estado geral. Exigindo a radioterapia um posicionamento e imobilização precisos e rigorosos pode, por vezes, causar algum desconforto aos doentes. No caso dos doentes com necessidades paliativas, nomeadamente com dor moderada a forte, é importante inferir acerca do real benefício de um tratamento que ao invés de promover alívio, lhes pode provocar mais dor e desconforto, indo contra às boas práticas preconizadas em cuidados paliativos. Considera-se que o tema é relevante para a melhoria dos cuidados prestados ao doente oncológico a realizar tratamentos de radioterapia, na medida em que permite realizar uma reflexão acerca da eficácia deste tratamento em fim de vida e sua aplicabilidade quando estamos perante doentes com prognóstico de vida limitado. O estudo em causa é um estudo descritivo exploratório, tendo sido realizada uma análise sistemática da informação contida nos processos electrónicos dos doentes, disponíveis na rede hospitalar, de acordo com os critérios de inclusão estabelecidos, no sentido de averiguar a eficácia do tratamento de radioterapia no alívio da dor. O estudo permitiu concluir que a radioterapia é um tratamento eficaz no controlo da dor em doentes com metástases ósseas, não vindo alterar as conclusões avançadas por outros estudos e outros autores internacionais. No entanto, é importante a realização de um prognóstico mais preciso na tomada de decisão terapêutica para que doentes que não vão beneficiar do tratamento de radioterapia possam ser referenciados para cuidados paliativos mais precocemente. Quanto à escolha do esquema terapêutico, o esquema de fracção única deveria ser mais comummente utilizado pela sua demonstrada eficácia, como o comprovam os diversos estudos referenciados neste trabalho.

Void fraction measurement and signal analysis from multiple-electrode impedance sensors

Void fraction sensors are important instruments not only for monitoring two-phase flow, but for furnishing an important parameter for obtaining flow map pattern and two-phase flow heat transfer coefficient as well. This work presents the experimental results obtained with the analysis of two axially spaced multiple-electrode impedance sensors tested in an upward air-water two-phase flow in a vertical tube for void fraction measurements. An electronic circuit was developed for signal generation and post-treatment of each sensor signal. By phase shifting the electrodes supplying the signal, it was possible to establish a rotating electric field sweeping across the test section. The fundamental principle of using a multiple-electrode configuration is based on reducing signal sensitivity to the non-uniform cross-section void fraction distribution problem. Static calibration curves were obtained for both sensors, and dynamic signal analyses for bubbly, slug, and turbulent churn flows were carried out. Flow parameters such as Taylor bubble velocity and length were obtained by using cross-correlation techniques. As an application of the void fraction tested, vertical flow pattern identification could be established by using the probability density function technique for void fractions ranging from 0% to nearly 70%.

A multiple time scale survival model with a cure fraction

Many recent survival studies propose modeling data with a cure fraction, i.e., data in which part of the population is not susceptible to the event of interest. This event may occur more than once for the same individual (recurrent event). We then have a scenario of recurrent event data in the presence of a cure fraction, which may appear in various areas such as oncology, finance, industries, among others. This paper proposes a multiple time scale survival model to analyze recurrent events using a cure fraction. The objective is analyzing the efficiency of certain interventions so that the studied event will not happen again in terms of covariates and censoring. All estimates were obtained using a sampling-based approach, which allows information to be input beforehand with lower computational effort. Simulations were done based on a clinical scenario in order to observe some frequentist properties of the estimation procedure in the presence of small and moderate sample sizes. An application of a well-known set of real mammary tumor data is provided.

Lossless filter for multiple repeats with bounded edit distance

Background: Identifying local similarity between two or more sequences, or identifying repeats occurring at least twice in a sequence, is an essential part in the analysis of biological sequences and of their phylogenetic relationship. Finding such fragments while allowing for a certain number of insertions, deletions, and substitutions, is however known to be a computationally expensive task, and consequently exact methods can usually not be applied in practice. Results: The filter TUIUIU that we introduce in this paper provides a possible solution to this problem. It can be used as a preprocessing step to any multiple alignment or repeats inference method, eliminating a possibly large fraction of the input that is guaranteed not to contain any approximate repeat. It consists in the verification of several strong necessary conditions that can be checked in a fast way. We implemented three versions of the filter. The first is simply a straightforward extension to the case of multiple sequences of an application of conditions already existing in the literature. The second uses a stronger condition which, as our results show, enable to filter sensibly more with negligible (if any) additional time. The third version uses an additional condition and pushes the sensibility of the filter even further with a non negligible additional time in many circumstances; our experiments show that it is particularly useful with large error rates. The latter version was applied as a preprocessing of a multiple alignment tool, obtaining an overall time (filter plus alignment) on average 63 and at best 530 times smaller than before (direct alignment), with in most cases a better quality alignment. Conclusion: To the best of our knowledge, TUIUIU is the first filter designed for multiple repeats and for dealing with error rates greater than 10% of the repeats length.

An isolated perfused human placental lobule model for multiple indicator dilution studies

We describe a method for multiple indicator dilution studies in the isolated perfused human placental lobule developed to investigate the relationships between changes in pressure and flow and solute clearance. A peripheral lobule of a human placenta is perfused with a tissue culture-based medium and the perfusate oxygen tension, arterial and venous pressures, pH and perfusion temperature continuously monitored by a computerized system. Flow rates are readily changed. Bolus injections of vascular, extracellular and water space markers, and study compounds can be made into either maternal or fetal circulations, and precisely timed outflow fractions can be collected with computer-controlled fraction collectors, allowing simultaneous determination of concentration-time profiles of each marker. (C) 1997 Elsevier Science Inc.

Use of multiple channel pH monitoring for evaluation of ultra-distal reflux in patients after fundoplication for treatment of Barrett`s esophagus

Dysplasia and esophageal adenocarcinoma may arise in patients with Barrett`s esophagus after fundoplication esophageal pH monitoring showing no acid in esophagus. This suggests the need to develop methodology to evaluate the occurrence of ultra-distal reflux (1 cm above the LES). The objective of the study was to compare acid exposition in three different levels: 5 cm above the upper border of the LES, 1 cm above the LES and in the intrasphincteric region. Eleven patients with Barrett`s esophagus after Nissen fundoplication with no clinical, endoscopic and radiologic evidence of reflux were selected. Four-channel pH monitoring took place: channel A, 5 cm above the upper border of the LES; channel B, 1 cm above the LES; channel C, intrasphincteric; channel D, intragastric. The results of channels A, B and C were compared. There was significant increase in number of reflux episodes and a higher fraction of time with pH <4.0 in channel B compared to channel A. There was significant decrease in fraction of time with pH <4.0 in channel B compared to channel C. Two cases of esophageal adenocarcinoma were diagnosed in the studied patients. The region 1 cm above the upper border of the LES is more exposed to acid than the region 5 cm above the upper border of the LES, although this exposure occurred in reduced levels. The region 1 cm above the upper border of the LES is less exposed to acid than the intrasphincteric region.

Distribution kinetics of solutes in the isolated in-situ perfused rat head using the multiple indicator dilution technique and a physiological two-barrier model

The purpose of this study was to determine the pharmacokinetics of [C-14]diclofenac, [C-14]salicylate and [H-3]clonidine using a single pass rat head perfusion preparation. The head was perfused with 3-[N-morpholino] propane-sulfonic acid-buffered Ringer's solution. Tc-99m-red blood cells and a drug were injected in a bolus into the internal carotid artery and collected from the posterior facial vein over 28 min. A two-barrier stochastic organ model was used to estimate the statistical moments of the solutes. Plasma, interstitial and cellular distribution volumes for the solutes ranged from 1.0 mL (diclofenac) to 1.6 mL (salicylate), 2.0 mL (diclofenac) to 4.2 mL (water) and 3.9 mL (salicylate) to 20.9 mL (diclofenac), respectively. A comparison of these volumes to water indicated some exclusion of the drugs from the interstitial space and salicylate from the cellular space. Permeability-surface area (PS) products calculated from plasma to interstitial fluid permeation clearances (CLPI) (range 0.02-0.40 mL s(-1)) and fractions of solute unbound in the perfusate were in the order: diclofenac>salicylate >clonidine>sucrose (from 41.8 to 0.10 mL s(-1)). The slow efflux of diclofenac, compared with clonidine and salicylate, may be related to its low average unbound fraction in the cells. This work accounts for the tail of disposition curves in describing pharmacokinetics in the head.

Why do females mate with multiple males ? The sexually selected sperm hypothesis

One debated issues in evolutionary biology is, why in many species females mate with multiple males. Several hypotheses have been put forward, yet the benefits of multiple mating (here defined as mating with several males) remain unclear in many cases. The sperm sexual selection (SSS) hypothesis has been developed to account for the widespread occurrence of multiple mating in females. It argues that multiple mating by females may rapidly spread, when initially a small fraction of the females mate multiply, and if there is a heritable difference among males in one or several of the four characteristics: (1) the quantity of sperm they produce; (2) the success of their sperm in reaching and fertilizing an egg; (3) their ability to displace the sperm that females stored during previous mating; and (4) their ability to prevent any other male from subsequently introducing sperm (e.g., differential efficiency of mating plugs).

On-pump beating heart coronary surgery for high risk patients requiring emergency multiple coronary artery bypass grafting.

BACKGROUND: Cardiopulmonary bypass (CPB) with aortic cross-clamping and cardioplegic arrest remains the method of choice for patients requiring standard myocardial revascularization. Therefore, very high-risk patients presenting with acute coronary syndrome, unstable angina, onset of cardiac decompensation and requiring emergency multiple myocardial revascularization, can have a poor outcome. The on-pump beating heart technique can reduce the mortality and the morbidity in such a selected group of patients and this report describes our clinical experience. METHODS: Out of 290 patients operated for CABG from January 2005 to January 2006, 25 (8.6%) selected high-risk patients suffering from life threatening coronary syndrome (mean age 69 +/- 7 years) and requiring emergency multiple myocardial revascularization, underwent on-pump beating heart surgery. The mean pre-operative left ventricle ejection fraction (LVEF) was 27 +/- 8%. The majority of them (88%) suffered of tri-vessel coronary disease and 6 (24%) had a left main stump disease. Nine patients (35%) were on severe cardiac failure and seven among them (28%) received a pre-operative intra-aortic balloon pump. The pre-operative EuroScore rate was equal or above 8 in 18 patients (73%). RESULTS: All patients underwent on-pump-beating heart coronary revascularization. The mean number of graft/patient was 2.9 +/- 0.6 and the internal mammary artery was used in 23 patients (92%). The mean CPB time was 84 +/- 19 minutes. Two patients died during the recovery stay in the intensive care unit, and there were no postoperative myocardial infarctions between the survivors. Eight patients suffered of transitorily renal failure and 1 patient developed a sternal wound infection. The mean hospital stay was 12 +/- 7 days. The follow-up was complete for all 23 patients survived at surgery and the mean follow-up time was 14 +/- 5 months. One patient died during the follow-up for cardiac arrest and 2 patients required an implantable cardiac defibrillator. One year after surgery they all had a standard trans-thoracic echocardiogram showing a mean LVEF rate of 36 +/- 11.8%. CONCLUSION: Standard on-pump arrested heart coronary surgery has higher mortality and morbidity in emergencies. The on-pump beating heart myocardial revascularization seems to be a valid alternative for the restricted and selected cohort of patients suffering from life threatening coronary syndrome and requiring multiple emergency CABG.

A bi-dimensional genome scan for prolificacy traits in pigs shows the existence of multiple epistatic QTL

Background: Prolificacy is the most important trait influencing the reproductive efficiency of pig production systems. The low heritability and sex-limited expression of prolificacy have hindered to some extent the improvement of this trait through artificial selection. Moreover, the relative contributions of additive, dominant and epistatic QTL to the genetic variance of pig prolificacy remain to be defined. In this work, we have undertaken this issue by performing one-dimensional and bi-dimensional genome scans for number of piglets born alive (NBA) and total number of piglets born (TNB) in a three generation Iberian by Meishan F2 intercross. Results: The one-dimensional genome scan for NBA and TNB revealed the existence of two genome-wide highly significant QTL located on SSC13 (P < 0.001) and SSC17 (P < 0.01) with effects on both traits. This relative paucity of significant results contrasted very strongly with the wide array of highly significant epistatic QTL that emerged in the bi-dimensional genome-wide scan analysis. As much as 18 epistatic QTL were found for NBA (four at P < 0.01 and five at P < 0.05) and TNB (three at P < 0.01 and six at P < 0.05), respectively. These epistatic QTL were distributed in multiple genomic regions, which covered 13 of the 18 pig autosomes, and they had small individual effects that ranged between 3 to 4% of the phenotypic variance. Different patterns of interactions (a × a, a × d, d × a and d × d) were found amongst the epistatic QTL pairs identified in the current work.Conclusions: The complex inheritance of prolificacy traits in pigs has been evidenced by identifying multiple additive (SSC13 and SSC17), dominant and epistatic QTL in an Iberian × Meishan F2 intercross. Our results demonstrate that a significant fraction of the phenotypic variance of swine prolificacy traits can be attributed to first-order gene-by-gene interactions emphasizing that the phenotypic effects of alleles might be strongly modulated by the genetic background where they segregate.

A pipeline approach with spatial information for segmenting multiple sclerosis lesions on brain magnetic resonance imaging

Background: Conventional magnetic resonance imaging (MRI) techniques are highly sensitive to detect multiple sclerosis (MS) plaques, enabling a quantitative assessment of inflammatory activity and lesion load. In quantitative analyses of focal lesions, manual or semi-automated segmentations have been widely used to compute the total number of lesions and the total lesion volume. These techniques, however, are both challenging and time-consuming, being also prone to intra-observer and inter-observer variability.Aim: To develop an automated approach to segment brain tissues and MS lesions from brain MRI images. The goal is to reduce the user interaction and to provide an objective tool that eliminates the inter- and intra-observer variability.Methods: Based on the recent methods developed by Souplet et al. and de Boer et al., we propose a novel pipeline which includes the following steps: bias correction, skull stripping, atlas registration, tissue classification, and lesion segmentation. After the initial pre-processing steps, a MRI scan is automatically segmented into 4 classes: white matter (WM), grey matter (GM), cerebrospinal fluid (CSF) and partial volume. An expectation maximisation method which fits a multivariate Gaussian mixture model to T1-w, T2-w and PD-w images is used for this purpose. Based on the obtained tissue masks and using the estimated GM mean and variance, we apply an intensity threshold to the FLAIR image, which provides the lesion segmentation. With the aim of improving this initial result, spatial information coming from the neighbouring tissue labels is used to refine the final lesion segmentation.Results:The experimental evaluation was performed using real data sets of 1.5T and the corresponding ground truth annotations provided by expert radiologists. The following values were obtained: 64% of true positive (TP) fraction, 80% of false positive (FP) fraction, and an average surface distance of 7.89 mm. The results of our approach were quantitatively compared to our implementations of the works of Souplet et al. and de Boer et al., obtaining higher TP and lower FP values.Conclusion: Promising MS lesion segmentation results have been obtained in terms of TP. However, the high number of FP which is still a well-known problem of all the automated MS lesion segmentation approaches has to be improved in order to use them for the standard clinical practice. Our future work will focus on tackling this issue.

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.

An allele-specific, stochastic gene expression process controls the expression of multiple Ly49 family genes and generates a diverse, MHC-specific NK cell receptor repertoire.

Mouse NK cells express MHC class I-specific inhibitory Ly49 receptors. Since these receptors display distinct ligand specificities and are clonally distributed, their expression generates a diverse NK cell receptor repertoire specific for MHC class I molecules. We have previously found that the Dd (or Dk)-specific Ly49A receptor is usually expressed from a single allele. However, a small fraction of short-term NK cell clones expressed both Ly49A alleles, suggesting that the two Ly49A alleles are independently and randomly expressed. Here we show that the genes for two additional Ly49 receptors (Ly49C and Ly49G2) are also expressed in a (predominantly) mono-allelic fashion. Since single NK cells can co-express multiple Ly49 receptors, we also investigated whether mono-allelic expression from within the tightly linked Ly49 gene cluster is coordinate or independent. Our clonal analysis suggests that the expression of alleles of distinct Ly49 genes is not coordinate. Thus Ly49 alleles are apparently independently and randomly chosen for stable expression, a process that directly restricts the number of Ly49 receptors expressed per single NK cell. We propose that the Ly49 receptor repertoire specific for MHC class I is generated by an allele-specific, stochastic gene expression process that acts on the entire Ly49 gene cluster.