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.

Truncated robust distance for clinical laboratory safety data monitoring and assessment.

Laboratory safety data are routinely collected in clinical studies for safety monitoring and assessment. We have developed a truncated robust multivariate outlier detection method for identifying subjects with clinically relevant abnormal laboratory measurements. The proposed method can be applied to historical clinical data to establish a multivariate decision boundary that can then be used for future clinical trial laboratory safety data monitoring and assessment. Simulations demonstrate that the proposed method has the ability to detect relevant outliers while automatically excluding irrelevant outliers. Two examples from actual clinical studies are used to illustrate the use of this method for identifying clinically relevant outliers.

The Tightness of the Cotton Swabs Meshing Influences the Chances of Getting Conclusive DNA Profiles

Objective: The chance of obtaining a conclusive DNA profile strongly depends on the quantity of biological material that can be recovered from a crime scene sample. Optimizing the collection strategy is therefore of prime interest. A difference in the level of tightness of the cotton meshed around the shaft has been observed between manufacturers and is hypothesized to affect the collection and subsequent release capacity of cotton swabs. Consequently, we compared the performance of cotton swabs from two different suppliers: Applimed SA and DryswabTM. Methods: These swabs were used to recover 50 ml of blood, either pure or diluted (1:1000 and 1:5000), deposited on both smooth and absorbent surfaces. Performance was compared in terms of ease of use, concentration of extracted DNA, and quality of DNA profiles. DNA quantification was obtained by real-time PCR using the QuantifilerTM Human DNA Quantification Kit. Evaluation of DNA profiles was based on profiles obtained using AmpFlSTR® NGM SElectTM PCR Amplification kit. Results: When considering smooth surfaces, recovered DNA was more concentrated when using the DryswabTM than the Applimed SA cotton swab. More precisely, DNA concentrations ranged from 15.7 to 28.8 ng/ml and 6.7 to 21.2 ng/ml, respectively for samples of pure blood. The same trend was observed for the absorbent surface, with 2.0 to 5.0 ng/ml and 0.9 to 1.4 ng/ml, respectively. Conclusion: Our results illustrate that different cotton swabs produce different results in terms of ease of use and quantity of recovered DNA and this should be taken into consideration when choosing which swab to use at both the crime scene and laboratory. More specifically, results from the present study suggest that looser meshing of the cotton fibres

Discovery and refinement of loci associated with lipid levels.

Levels of low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol, triglycerides and total cholesterol are heritable, modifiable risk factors for coronary artery disease. To identify new loci and refine known loci influencing these lipids, we examined 188,577 individuals using genome-wide and custom genotyping arrays. We identify and annotate 157 loci associated with lipid levels at P < 5 × 10(-8), including 62 loci not previously associated with lipid levels in humans. Using dense genotyping in individuals of European, East Asian, South Asian and African ancestry, we narrow association signals in 12 loci. We find that loci associated with blood lipid levels are often associated with cardiovascular and metabolic traits, including coronary artery disease, type 2 diabetes, blood pressure, waist-hip ratio and body mass index. Our results demonstrate the value of using genetic data from individuals of diverse ancestry and provide insights into the biological mechanisms regulating blood lipids to guide future genetic, biological and therapeutic research.

Genetic and phenotypic architecture of metabolic syndrome-associated components in dyslipidemic and normolipidemic subjects: the GEMS Study.

Atherogenic dyslipidemia, manifest by low HDL-cholesterol and high TG levels, is an important component of ATP-III defined metabolic syndrome. Here, we dissected the phenotypic and genetic architecture of these traits by assessing their relationships with other metabolically relevant measures, including plasma adipo-cytokines, highly sensitive C-reactive protein (hsCRP) and LDL particle size, in a large family data set (n=2800) and in an independent set of dyslipidemic cases (n=716) and normolipidemic controls (n=1073). We explored the relationships among these phenotypes using variable clustering and then estimated their genetic heritabilities and cross-trait correlations. In families, four clusters explained 61% of the total variance, with one adiposity-related cluster (including hsCRP), one BP-related cluster, and two lipid-related clusters (HDL-C, TG, adiponectin and LDL particle size; apoB and non-HDL-C). A similar structure was observed in dyslipidemic cases and normolipidemic controls. The genetic correlations in the families largely paralleled the phenotype clustering results, suggesting that common genes having pleiotropic effects contributed to the correlations observed. In summary, our analyses support a model of metabolic syndrome with two major components, body fat and lipids, each with two subcomponents, and quantifies their degree of overlap with each other and with metabolic-syndrome related measures (adipokines, LDL particle size and hsCRP).

Guidelines for the use and interpretation of assays for monitoring autophagy.

In 2008 we published the first set of guidelines for standardizing research in autophagy. Since then, research on this topic has continued to accelerate, and many new scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Accordingly, it is important to update these guidelines for monitoring autophagy in different organisms. Various reviews have described the range of assays that have been used for this purpose. Nevertheless, there continues to be confusion regarding acceptable methods to measure autophagy, especially in multicellular eukaryotes. A key point that needs to be emphasized is that there is a difference between measurements that monitor the numbers or volume of autophagic elements (e.g., autophagosomes or autolysosomes) at any stage of the autophagic process vs. those that measure flux through the autophagy pathway (i.e., the complete process); thus, a block in macroautophagy that results in autophagosome accumulation needs to be differentiated from stimuli that result in increased autophagic activity, defined as increased autophagy induction coupled with increased delivery to, and degradation within, lysosomes (in most higher eukaryotes and some protists such as Dictyostelium) or the vacuole (in plants and fungi). In other words, it is especially important that investigators new to the field understand that the appearance of more autophagosomes does not necessarily equate with more autophagy. In fact, in many cases, autophagosomes accumulate because of a block in trafficking to lysosomes without a concomitant change in autophagosome biogenesis, whereas an increase in autolysosomes may reflect a reduction in degradative activity. Here, we present a set of guidelines for the selection and interpretation of methods for use by investigators who aim to examine macroautophagy and related processes, as well as for reviewers who need to provide realistic and reasonable critiques of papers that are focused on these processes. These guidelines are not meant to be a formulaic set of rules, because the appropriate assays depend in part on the question being asked and the system being used. In addition, we emphasize that no individual assay is guaranteed to be the most appropriate one in every situation, and we strongly recommend the use of multiple assays to monitor autophagy. In these guidelines, we consider these various methods of assessing autophagy and what information can, or cannot, be obtained from them. Finally, by discussing the merits and limits of particular autophagy assays, we hope to encourage technical innovation in the field.

Covalent assembly of a soluble T cell receptor-peptide-major histocompatibility class I complex.

We used stepwise photochemical cross-linking for specifically assembling soluble and covalent complexes made of a T-cell antigen receptor (TCR) and a class I molecule of the major histocompatibility complex (MHC) bound to an antigenic peptide. For that purpose, we have produced in myeloma cells a single-chain Fv construct of a TCR specific for a photoreactive H-2Kd-peptide complex. Photochemical cross-linking of this TCR single-chain Fv with a soluble form of the photoreactive H-2Kd-peptide ligand resulted in the formation of a ternary covalent complex. We have characterized the soluble ternary complex and showed that it reacted with antibodies specific for epitopes located either on the native TCR or on the Kd molecules. By preventing the fast dissociation kinetics observed with most T cell receptors, this approach provides a means of preparing soluble TCR-peptide-MHC complexes on large-scale levels.

Common variants associated with plasma triglycerides and risk for coronary artery disease.

Triglycerides are transported in plasma by specific triglyceride-rich lipoproteins; in epidemiological studies, increased triglyceride levels correlate with higher risk for coronary artery disease (CAD). However, it is unclear whether this association reflects causal processes. We used 185 common variants recently mapped for plasma lipids (P < 5 × 10(-8) for each) to examine the role of triglycerides in risk for CAD. First, we highlight loci associated with both low-density lipoprotein cholesterol (LDL-C) and triglyceride levels, and we show that the direction and magnitude of the associations with both traits are factors in determining CAD risk. Second, we consider loci with only a strong association with triglycerides and show that these loci are also associated with CAD. Finally, in a model accounting for effects on LDL-C and/or high-density lipoprotein cholesterol (HDL-C) levels, the strength of a polymorphism's effect on triglyceride levels is correlated with the magnitude of its effect on CAD risk. These results suggest that triglyceride-rich lipoproteins causally influence risk for CAD.

The role of transcription factories-mediated interchromosomal contacts in the organization of nuclear architecture.

Using numerical simulations, we investigate the underlying physical effects responsible for the overall organization of chromosomal territories in interphase nuclei. In particular, we address the following three questions: (i) why are chromosomal territories with relatively high transcriptional activity on average, closer to the centre of cell's nucleus than those with the lower activity? (ii) Why are actively transcribed genes usually located at the periphery of their chromosomal territories? (iii) Why are pair-wise contacts between active and inactive genes less frequent than those involving only active or only inactive genes? We show that transcription factories-mediated contacts between active genes belonging to different chromosomal territories are instrumental for all these features of nuclear organization to emerge spontaneously due to entropic effects arising when chromatin fibres are highly crowded.

Cellular and molecular mechanisms underlying the effects of sleep loss on hippocampal synaptic plasticity

SUMMARY : The function of sleep for the organism is one of the most persistent and perplexing questions in biology. Current findings lead to the conclusion that sleep is primarily for the brain. In particular, a role for sleep in cognitive aspects of brain function is supported by behavioral evidence both in humans and animals. However, in spite of remarkable advancement in the understanding of the mechanisms underlying sleep generation and regulation, it has been proven difficult to determine the neurobiological mechanisms underlying the beneficial effect of sleep, and the detrimental impact of sleep loss, on learning and memory processes. In my thesis, I present results that lead to several critical steps forward in the link between sleep and cognitive function. My major result is the molecular identification and physiological analysis of a protein, the NR2A subunit of NMDA receptor (NMDAR), that confers sensitivity to sleep loss to the hippocampus, a brain structure classically involved in mnemonic processes. Specifically, I used a novel behavioral approach to achieve sleep deprivation in adult C57BL6/J mice, yet minimizing the impact of secondary factors associated with the procedure,.such as stress. By using in vitro electrophysiological analysis, I show, for the first time, that sleep loss dramatically affects bidirectional plasticity at CA3 to CA1 synapses in the hippocampus, a well established cellular model of learning and memory. 4-6 hours of sleep loss elevate the modification threshold for bidirectional synaptic plasticity (MT), thereby promoting long-term depression of CA3 to CA 1 synaptic strength after stimulation in the theta frequency range (5 Hz), and rendering long-term potentiation induction.more difficult. Remarkably, 3 hours of recovery sleep, after the deprivation, reset the MT at control values, thus re-establishing the normal proneness of synapses to undergo long-term plastic changes. At the molecular level, these functional changes are paralleled by a change in the NMDAR subunit composition. In particular, the expression of the NR2A subunit protein of NMDAR at CA3 to CA1 synapses is selectively and rapidly increased by sleep deprivation, whereas recovery sleep reset NR2A synaptic content to control levels. By using an array of genetic, pharmacological and computational approaches, I demonstrate here an obligatory role for NR2A-containing NMDARs in conveying the effect of sleep loss on CA3 to CAl MT. Moreover, I show that a genetic deletion of the NR2A subunit fully preserves hippocampal plasticity from the impact of sleep loss, whereas it does not alter sleepwake behavior and homeostatic response to sleep deprivation. As to the mechanism underlying the effects of the NR2A subunit on hippocampal synaptic plasticity, I show that the increased NR2A expression after sleep loss distinctly affects the contribution of synaptic and more slowly recruited NMDAR pools activated during plasticity-induction protocols. This study represents a major step forward in understanding the mechanistic basis underlying sleep's role for the brain. By showing that sleep and sleep loss affect neuronal plasticity by regulating the expression and function of a synaptic neurotransmitter receptor, I propose that an important aspect of sleep function could consist in maintaining and regulating protein redistribution and ion channel trafficking at central synapses. These findings provide a novel starting point for investigations into the connections between sleep and learning, and they may open novel ways for pharmacological control over hippocampal .function during periods of sleep restriction. RÉSUMÉ DU PROJET La fonction du sommeil pour l'organisme est une des questions les plus persistantes et difficiles dans la biologie. Les découvertes actuelles mènent à la conclusion que le sommeil est essentiel pour le cerveau. En particulier, le rôle du sommeil dans les aspects cognitifs est soutenu par des études comportementales tant chez les humains que chez les animaux. Cependant, malgré l'avancement remarquable dans la compréhension des mécanismes sous-tendant la génération et la régulation du sommeil, les mécanismes neurobiologiques qui pourraient expliquer l'effet favorable du sommeil sur l'apprentissage et la mémoire ne sont pas encore clairs. Dans ma thèse, je présente des résultats qui aident à clarifier le lien entre le sommeil et la fonction cognitive. Mon résultat le plus significatif est l'identification moléculaire et l'analyse physiologique d'une protéine, la sous-unité NR2A du récepteur NMDA, qui rend l'hippocampe sensible à la perte de sommeil. Dans cette étude, nous avons utilisé une nouvelle approche expérimentale qui nous a permis d'induire une privation de sommeil chez les souris C57BL6/J adultes, en minimisant l'impact de facteurs confondants comme, par exemple, le stress. En utilisant les techniques de l'électrophysiologie in vitro, j'ai démontré, pour la première fois, que la perte de sommeil est responsable d'affecter radicalement la plasticité bidirectionnelle au niveau des synapses CA3-CA1 de l'hippocampe. Cela correspond à un mécanisme cellulaire de l'apprentissage et de la mémoire bien établi. En particulier, 4-6 heures de privation de sommeil élèvent le seuil de modification pour la plasticité synaptique bidirectionnelle (SM). Comme conséquence, la dépression à long terme de la transmission synaptique est induite par la stimulation des fibres afférentes dans la bande de fréquences thêta (5 Hz), alors que la potentialisation à long terme devient plus difficile. D'autre part, 3 heures de sommeil de récupération sont suffisant pour rétablir le SM aux valeurs contrôles. Au niveau moléculaire, les changements de la plasticité synaptiques sont associés à une altération de la composition du récepteur NMDA. En particulier, l'expression synaptique de la protéine NR2A du récepteur NMDA est rapidement augmentée de manière sélective par la privation de sommeil, alors que le sommeil de récupération rétablit l'expression de la protéine au niveau contrôle. En utilisant des approches génétiques, pharmacologiques et computationnelles, j'ai démontré que les récepteurs NMDA qui expriment la sous-unité NR2A sont responsables de l'effet de la privation de sommeil sur le SM. De plus, nous avons prouvé qu'une délétion génétique de la sous-unité NR2A préserve complètement la plasticité synaptique hippocampale de l'impact de la perte de sommeil, alors que cette manipulation ne change pas les mécanismes de régulation homéostatique du sommeil. En ce qui concerne les mécanismes, j'ai .découvert que l'augmentation de l'expression de la sous-unité NR2A au niveau synaptique modifie les propriétés de la réponse du récepteur NMDA aux protocoles de stimulations utilisés pour induire la plasticité. Cette étude représente un pas en avant important dans la compréhension de la base mécaniste sous-tendant le rôle du sommeil pour le cerveau. En montrant que le sommeil et la perte de sommeil affectent la plasticité neuronale en régulant l'expression et la fonction d'un récepteur de la neurotransmission, je propose qu'un aspect important de la fonction du sommeil puisse être finalisé au règlement de la redistribution des protéines et du tracking des récepteurs aux synapses centraux. Ces découvertes fournissent un point de départ pour mieux comprendre les liens entre le sommeil et l'apprentissage, et d'ailleurs, ils peuvent ouvrir des voies pour des traitements pharmacologiques dans le .but de préserver la fonction hippocampale pendant les périodes de restriction de sommeil.

GENCODE: the reference human genome annotation for The ENCODE Project.

The GENCODE Consortium aims to identify all gene features in the human genome using a combination of computational analysis, manual annotation, and experimental validation. Since the first public release of this annotation data set, few new protein-coding loci have been added, yet the number of alternative splicing transcripts annotated has steadily increased. The GENCODE 7 release contains 20,687 protein-coding and 9640 long noncoding RNA loci and has 33,977 coding transcripts not represented in UCSC genes and RefSeq. It also has the most comprehensive annotation of long noncoding RNA (lncRNA) loci publicly available with the predominant transcript form consisting of two exons. We have examined the completeness of the transcript annotation and found that 35% of transcriptional start sites are supported by CAGE clusters and 62% of protein-coding genes have annotated polyA sites. Over one-third of GENCODE protein-coding genes are supported by peptide hits derived from mass spectrometry spectra submitted to Peptide Atlas. New models derived from the Illumina Body Map 2.0 RNA-seq data identify 3689 new loci not currently in GENCODE, of which 3127 consist of two exon models indicating that they are possibly unannotated long noncoding loci. GENCODE 7 is publicly available from gencodegenes.org and via the Ensembl and UCSC Genome Browsers.

Méthodes sismiques actives adaptées à l'étude de milieux hétérogènes : application aux glissements de terrain

Résumé Les glissements de terrain représentent un des principaux risques naturels dans les régions montagneuses. En Suisse, chaque année les glissements de terrains causent des dégâts qui affectent les infrastructures et ont des coûts financiers importants. Une bonne compréhension des mécanismes des glissements peut permettre d'atténuer leur impact. Celle-ci passe notamment par la connaissance de la structure interne du glissement, la détermination de son volume et de son ou ses plans de glissement. Dans un glissement de terrain, la désorganisation et la présence de fractures dans le matériel déplacé engendre un changement des paramètres physiques et en particulier une diminution des vitesses de propagation des ondes sismiques ainsi que de la densité du matériel. Les méthodes sismiques sont de ce fait bien adaptées à l'étude des glissements de terrain. Parmi les méthodes sismiques, l'analyse de la dispersion des ondes de surface est une méthode simple à mettre en oeuvre. Elle présente l'avantage d'estimer les variations des vitesses de cisaillement avec la profondeur sans avoir spécifiquement recours à l'utilisation d'une source d'onde S et de géophones horizontaux. Sa mise en oeuvre en trois étapes implique la mesure de la dispersion des ondes de surface sur des réseaux étendus, la détermination des courbes de dispersion pour finir par l'inversion de ces courbes. Les modèles de vitesse obtenus à partir de cette procédure ne sont valides que lorsque les milieux explorés ne présentent pas de variations latérales. En pratique cette hypothèse est rarement vérifiée, notamment pour un glissement de terrain dans lequel les couches remaniées sont susceptibles de présenter de fortes hétérogénéités latérales. Pour évaluer la possibilité de déterminer des courbes de dispersion à partir de réseaux de faible extension des mesures testes ont été effectuées sur un site (Arnex, VD) équipé d'un forage. Un profil sismique de 190 m de long a été implanté dans une vallée creusée dans du calcaire et remplie par des dépôts glacio-lacustres d'une trentaine de mètres d'épaisseur. Les données acquises le long de ce profil ont confirmé que la présence de variations latérales sous le réseau de géophones affecte l'allure des courbes de dispersion jusqu'à parfois empêcher leur détermination. Pour utiliser l'analyse de la dispersion des ondes de surface sur des sites présentant des variations latérales, notre approche consiste à déterminer les courbes de dispersions pour une série de réseaux de faible extension, à inverser chacune des courbes et à interpoler les différents modèles de vitesse obtenus. Le choix de la position ainsi que de l'extension des différents réseaux de géophones est important. Il tient compte de la localisation des hétérogénéités détectées à partir de l'analyse de sismique réfraction, mais également d'anomalies d'amplitudes observées sur des cartes qui représentent dans le domaine position de tir - position du récepteur, l'amplitude mesurée pour différentes fréquences. La procédure proposée par Lin et Lin (2007) s'est avérée être une méthode efficace permettant de déterminer des courbes de dispersion à partir de réseaux de faible extension. Elle consiste à construire à partir d'un réseau de géophones et de plusieurs positions de tir un enregistrement temps-déports qui tient compte d'une large gamme de distances source-récepteur. Au moment d'assembler les différentes données une correction de phase est appliquée pour tenir compte des hétérogénéités situées entre les différents points de tir. Pour évaluer cette correction nous suggérons de calculer pour deux tir successif la densité spectrale croisée des traces de même offset: Sur le site d'Arnex, 22 courbes de dispersions ont été déterminées pour de réseaux de géophones de 10 m d'extension. Nous avons également profité du forage pour acquérir un profil de sismique verticale en ondes S. Le modèle de vitesse S déduit de l'interprétation du profil de sismique verticale est utilisé comme information à priori lors l'inversion des différentes courbes de dispersion. Finalement, le modèle en deux dimension qui a été établi grâce à l'analyse de la dispersion des ondes de surface met en évidence une structure tabulaire à trois couches dont les limites coïncident bien avec les limites lithologiques observées dans le forage. Dans celui-ci des argiles limoneuses associées à une vitesse de propagation des ondes S de l'ordre de 175 m/s surmontent vers 9 m de profondeur des dépôts de moraine argilo-sableuse caractérisés par des vitesses de propagation des ondes S de l'ordre de 300 m/s jusqu'à 14 m de profondeur et supérieur ou égal à 400 m/s entre 14 et 20 m de profondeur. Le glissement de la Grande Combe (Ballaigues, VD) se produit à l'intérieur du remplissage quaternaire d'une combe creusée dans des calcaires Portlandien. Comme dans le cas du site d'Arnex les dépôts quaternaires correspondent à des dépôts glacio-lacustres. Dans la partie supérieure la surface de glissement a été localisée à une vingtaine de mètres de profondeur au niveau de l'interface qui sépare des dépôts de moraine jurassienne et des dépôts glacio-lacustres. Au pied du glissement 14 courbes de dispersions ont été déterminées sur des réseaux de 10 m d'extension le long d'un profil de 144 m. Les courbes obtenues sont discontinues et définies pour un domaine de fréquence de 7 à 35 Hz. Grâce à l'utilisation de distances source-récepteur entre 8 et 72 m, 2 à 4 modes de propagation ont été identifiés pour chacune des courbes. Lors de l'inversion des courbes de dispersion la prise en compte des différents modes de propagation a permis d'étendre la profondeur d'investigation jusqu'à une vingtaine de mètres de profondeur. Le modèle en deux dimensions permet de distinguer 4 couches (Vs1 < 175 m/s, 175 m/s < Vs2 < 225 m/s, 225 m/s < Vs3 < 400 m/s et Vs4 >.400 m/s) qui présentent des variations d'épaisseur. Des profils de sismiques réflexion en ondes S acquis avec une source construite dans le cadre de ce travail, complètent et corroborent le modèle établi à partir de l'analyse de la dispersion des ondes de surface. Un réflecteur localisé entre 5 et 10 m de profondeur et associé à une vitesse de sommation de 180 m/s souligne notamment la géométrie de l'interface qui sépare la deuxième de la troisième couche du modèle établi à partir de l'analyse de la dispersion des ondes de surface. Abstract Landslides are one of the main natural hazards in mountainous regions. In Switzerland, landslides cause damages every year that impact infrastructures and have important financial costs. In depth understanding of sliding mechanisms may help limiting their impact. In particular, this can be achieved through a better knowledge of the internal structure of the landslide, the determination of its volume and its sliding surface or surfaces In a landslide, the disorganization and the presence of fractures in the displaced material generate a change of the physical parameters and in particular a decrease of the seismic velocities and of the material density. Therefoe, seismic methods are well adapted to the study of landslides. Among seismic methods, surface-wave dispersion analysis is a easy to implement. Through it, shearwave velocity variations with depth can be estimated without having to resort to an S-wave source and to horizontal geophones. Its 3-step implementation implies measurement of surface-wave dispersion with long arrays, determination of the dispersion curves and finally inversion of these curves. Velocity models obtained through this approach are only valid when the investigated medium does not include lateral variations. In practice, this assumption is seldom correct, in particular for landslides in which reshaped layers likely include strong lateral heterogeneities. To assess the possibility of determining dispersion curves from short array lengths we carried out tests measurements on a site (Arnex, VD) that includes a borehole. A 190 m long seismic profile was acquired in a valley carved into limestone and filled with 30 m of glacio-lacustrine sediments. The data acquired along this profile confirmed that the presence of lateral variations under the geophone array influences the dispersion-curve shape so much that it sometimes preventes the dispersion curves determination. Our approach to use the analysis of surface-wave dispersion on sites that include lateral variations consists in obtaining dispersion curves for a series of short length arrays; inverting each so obtained curve and interpolating the different obtained velocity model. The choice of the location as well as the geophone array length is important. It takes into account the location of the heterogeneities that are revealed by the seismic refraction interpretation of the data but also, the location of signal amplitude anomalies observed on maps that represent, for a given frequency, the measured amplitude in the shot position - receiver position domain. The procedure proposed by Lin and Lin (2007) turned out to be an efficient one to determine dispersion curves using short extension arrays. It consists in building a time-offset from an array of geophones with a wide offset range by gathering seismograms acquired with different source-to-receiver offsets. When assembling the different data, a phase correction is applied in order to reduce static phase error induced by lateral variation. To evaluate this correction, we suggest to calculate, for two successive shots, the cross power spectral density of common offset traces. On the Arnex site, 22 curves were determined with 10m in length geophone-arrays. We also took advantage of the borehole to acquire a S-wave vertical seismic profile. The S-wave velocity depth model derived from the vertical seismic profile interpretation is used as prior information in the inversion of the dispersion-curves. Finally a 2D velocity model was established from the analysis of the different dispersion curves. It reveals a 3-layer structure in good agreement with the observed lithologies in the borehole. In it a clay layer with a shear-wave of 175 m/s shear-wave velocity overlies a clayey-sandy till layer at 9 m depth that is characterized down to 14 m by a 300 m/s S-wave velocity; these deposits have a S-wave velocity of 400 m/s between depths of 14 to 20 m. The La Grand Combe landslide (Ballaigues, VD) occurs inside the Quaternary filling of a valley carved into Portlandien limestone. As at the Arnex site, the Quaternary deposits correspond to glaciolacustrine sediments. In the upper part of the landslide, the sliding surface is located at a depth of about 20 m that coincides with the discontinuity between Jurassian till and glacio-lacustrine deposits. At the toe of the landslide, we defined 14 dispersion curves along a 144 m long profile using 10 m long geophone arrays. The obtained curves are discontinuous and defined within a frequency range of 7 to 35 Hz. The use of a wide range of offsets (from 8 to 72 m) enabled us to determine 2 to 4 mode of propagation for each dispersion curve. Taking these higher modes into consideration for dispersion curve inversion allowed us to reach an investigation depth of about 20 m. A four layer 2D model was derived (Vs1< 175 m/s, 175 m/s <Vs2< 225 m/s, 225 m/s < Vs3 < 400 m/s, Vs4> 400 m/s) with variable layer thicknesses. S-wave seismic reflection profiles acquired with a source built as part of this work complete and the velocity model revealed by surface-wave analysis. In particular, reflector at a depth of 5 to 10 m associated with a 180 m/s stacking velocity image the geometry of the discontinuity between the second and third layer of the model derived from the surface-wave dispersion analysis.

Quantification of myocardial blood flow with 82Rb positron emission tomography: clinical validation with 15O-water.

PURPOSE: Quantification of myocardial blood flow (MBF) with generator-produced (82)Rb is an attractive alternative for centres without an on-site cyclotron. Our aim was to validate (82)Rb-measured MBF in relation to that measured using (15)O-water, as a tracer 100% of which can be extracted from the circulation even at high flow rates, in healthy control subject and patients with mild coronary artery disease (CAD). METHODS: MBF was measured at rest and during adenosine-induced hyperaemia with (82)Rb and (15)O-water PET in 33 participants (22 control subjects, aged 30 ± 13 years; 11 CAD patients without transmural infarction, aged 60 ± 13 years). A one-tissue compartment (82)Rb model with ventricular spillover correction was used. The (82)Rb flow-dependent extraction rate was derived from (15)O-water measurements in a subset of 11 control subjects. Myocardial flow reserve (MFR) was defined as the hyperaemic/rest MBF. Pearson's correlation r, Bland-Altman 95% limits of agreement (LoA), and Lin's concordance correlation ρ (c) (measuring both precision and accuracy) were used. RESULTS: Over the entire MBF range (0.66-4.7 ml/min/g), concordance was excellent for MBF (r = 0.90, [(82)Rb-(15)O-water] mean difference ± SD = 0.04 ± 0.66 ml/min/g, LoA = -1.26 to 1.33 ml/min/g, ρ(c) = 0.88) and MFR (range 1.79-5.81, r = 0.83, mean difference = 0.14 ± 0.58, LoA = -0.99 to 1.28, ρ(c) = 0.82). Hyperaemic MBF was reduced in CAD patients compared with the subset of 11 control subjects (2.53 ± 0.74 vs. 3.62 ± 0.68 ml/min/g, p = 0.002, for (15)O-water; 2.53 ± 1.01 vs. 3.82 ± 1.21 ml/min/g, p = 0.013, for (82)Rb) and this was paralleled by a lower MFR (2.65 ± 0.62 vs. 3.79 ± 0.98, p = 0.004, for (15)O-water; 2.85 ± 0.91 vs. 3.88 ± 0.91, p = 0.012, for (82)Rb). Myocardial perfusion was homogeneous in 1,114 of 1,122 segments (99.3%) and there were no differences in MBF among the coronary artery territories (p > 0.31). CONCLUSION: Quantification of MBF with (82)Rb with a newly derived correction for the nonlinear extraction function was validated against MBF measured using (15)O-water in control subjects and patients with mild CAD, where it was found to be accurate at high flow rates. (82)Rb-derived MBF estimates seem robust for clinical research, advancing a step further towards its implementation in clinical routine.