Cross recurrence quantification for cover song identification

There is growing evidence that nonlinear time series analysis techniques can be used to successfully characterize, classify, or process signals derived from realworld dynamics even though these are not necessarily deterministic and stationary. In the present study we proceed in this direction by addressing an important problem our modern society is facing, the automatic classification of digital information. In particular, we address the automatic identification of cover songs, i.e. alternative renditions of a previously recorded musical piece. For this purpose we here propose a recurrence quantification analysis measure that allows tracking potentially curved and disrupted traces in cross recurrence plots. We apply this measure to cross recurrence plots constructed from the state space representation of musical descriptor time series extracted from the raw audio signal. We show that our method identifies cover songs with a higher accuracy as compared to previously published techniques. Beyond the particular application proposed here, we discuss how our approach can be useful for the characterization of a variety of signals from different scientific disciplines. We study coupled Rössler dynamics with stochastically modulated mean frequencies as one concrete example to illustrate this point.

Identification of cancer stem cells in Ewing's sarcoma.

Cancer stem cells that display tumor-initiating properties have recently been identified in several distinct types of malignancies, holding promise for more effective therapeutic strategies. However, evidence of such cells in sarcomas, which include some of the most aggressive and therapy-resistant tumors, has not been shown to date. Here, we identify and characterize cancer stem cells in Ewing's sarcoma family tumors (ESFT), a highly aggressive pediatric malignancy believed to be of mesenchymal stem cell (MSC) origin. Using magnetic bead cell separation of primary ESFT, we have isolated a subpopulation of CD133+ tumor cells that display the capacity to initiate and sustain tumor growth through serial transplantation in nonobese diabetic/severe combined immunodeficiency mice, re-establishing at each in vivo passage the parental tumor phenotype and hierarchical cell organization. Consistent with the plasticity of MSCs, in vitro differentiation assays showed that the CD133+ cell population retained the ability to differentiate along adipogenic, osteogenic, and chondrogenic lineages. Quantitative real-time PCR analysis of genes implicated in stem cell maintenance revealed that CD133+ ESFT cells express significantly higher levels of OCT4 and NANOG than their CD133- counterparts. Taken together, our observations provide the first identification of ESFT cancer stem cells and demonstration of their MSC properties, a critical step towards a better biological understanding and rational therapeutic targeting of these tumors.

Visual functions in the healthy and schizophrenic brain : from stimulus control to illusory perceptions using electrical neuroimaging

ABSTRACT (FRENCH)Ce travail de thèse basé sur le système visuel chez les sujets sains et chez les patients schizophrènes, s'articule autour de trois articles scientifiques publiés ou en cours de publication. Ces articles traitent des sujets suivants : le premier article présente une nouvelle méthode de traitement des composantes physiques des stimuli (luminance et fréquence spatiale). Le second article montre, à l'aide d'analyses de données EEG, un déficit de la voie magnocellulaire dans le traitement visuel des illusions chez les patients schizophrènes. Ceci est démontré par l'absence de modulation de la composante PI chez les patients schizophrènes contrairement aux sujets sains. Cette absence est induite par des stimuli de type illusion Kanizsa de différentes excentricités. Finalement, le troisième article, également à l'aide de méthodes de neuroimagerie électrique (EEG), montre que le traitement des contours illusoires se trouve dans le complexe latéro-occipital (LOC), à l'aide d'illusion « misaligned gratings ». De plus il révèle que les activités démontrées précédemment dans les aires visuelles primaires sont dues à des inférences « top- down ».Afin de permettre la compréhension de ces trois articles, l'introduction de ce manuscrit présente les concepts essentiels. De plus des méthodes d'analyses de temps-fréquence sont présentées. L'introduction est divisée en quatre parties : la première présente le système visuel depuis les cellules retino-corticales aux deux voix du traitement de l'information en passant par les régions composant le système visuel. La deuxième partie présente la schizophrénie par son diagnostic, ces déficits de bas niveau de traitement des stimuli visuel et ces déficits cognitifs. La troisième partie présente le traitement des contours illusoires et les trois modèles utilisés dans le dernier article. Finalement, les méthodes de traitement des données EEG seront explicitées, y compris les méthodes de temps-fréquences.Les résultats des trois articles sont présentés dans le chapitre éponyme (du même nom). De plus ce chapitre comprendra les résultats obtenus à l'aide des méthodes de temps-fréquenceFinalement, la discussion sera orientée selon trois axes : les méthodes de temps-fréquence ainsi qu'une proposition de traitement de ces données par une méthode statistique indépendante de la référence. La discussion du premier article en montrera la qualité du traitement de ces stimuli. La discussion des deux articles neurophysiologiques, proposera de nouvelles d'expériences afin d'affiner les résultats actuels sur les déficits des schizophrènes. Ceci pourrait permettre d'établir un marqueur biologique fiable de la schizophrénie.ABSTRACT (ENGLISH)This thesis focuses on the visual system in healthy subjects and schizophrenic patients. To address this research, advanced methods of analysis of electroencephalographic (EEG) data were used and developed. This manuscript is comprised of three scientific articles. The first article showed a novel method to control the physical features of visual stimuli (luminance and spatial frequencies). The second article showed, using electrical neuroimaging of EEG, a deficit in spatial processing associated with the dorsal pathway in chronic schizophrenic patients. This deficit was elicited by an absent modulation of the PI component in terms of response strength and topography as well as source estimations. This deficit was orthogonal to the preserved ability to process Kanizsa-type illusory contours. Finally, the third article resolved ongoing debates concerning the neural mechanism mediating illusory contour sensitivity by using electrical neuroimaging to show that the first differentiation of illusory contour presence vs. absence is localized within the lateral occipital complex. This effect was subsequent to modulations due to the orientation of misaligned grating stimuli. Collectively, these results support a model where effects in V1/V2 are mediated by "top-down" modulation from the LOC.To understand these three articles, the Introduction of this thesis presents the major concepts used in these articles. Additionally, a section is devoted to time-frequency analysis methods not presented in the articles themselves. The introduction is divided in four parts. The first part presents three aspects of the visual system: cellular, regional, and its functional interactions. The second part presents an overview of schizophrenia and its sensoiy-cognitive deficits. The third part presents an overview of illusory contour processing and the three models examined in the third article. Finally, advanced analysis methods for EEG are presented, including time- frequency methodology.The Introduction is followed by a synopsis of the main results in the articles as well as those obtained from the time-frequency analyses.Finally, the Discussion chapter is divided along three axes. The first axis discusses the time frequency analysis and proposes a novel statistical approach that is independent of the reference. The second axis contextualizes the first article and discusses the quality of the stimulus control and direction for further improvements. Finally, both neurophysiologic articles are contextualized by proposing future experiments and hypotheses that may serve to improve our understanding of schizophrenia on the one hand and visual functions more generally.

Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Is the observed persistence spurious? A test for fractional integration versus short memory and structural breaks

Although it is commonly accepted that most macroeconomic variables are nonstationary, it is often difficult to identify the source of the non-stationarity. In particular, it is well-known that integrated and short memory models containing trending components that may display sudden changes in their parameters share some statistical properties that make their identification a hard task. The goal of this paper is to extend the classical testing framework for I(1) versus I(0)+ breaks by considering a a more general class of models under the null hypothesis: non-stationary fractionally integrated (FI) processes. A similar identification problem holds in this broader setting which is shown to be a relevant issue from both a statistical and an economic perspective. The proposed test is developed in the time domain and is very simple to compute. The asymptotic properties of the new technique are derived and it is shown by simulation that it is very well-behaved in finite samples. To illustrate the usefulness of the proposed technique, an application using inflation data is also provided.

The impact of generic reference pricing interventions in the statin market

The objective of this study was to evaluate the intended and unintended impact on pharmaceutical use and sales of three public financing reforms applied to the prescription of statins: a Spanish generic reference pricing (RP) system for lovastatin and simvastatin, and two competing policies introduced by the Andalusian Public Health Service (APHS) for all statins, first a maximum consumer price (MCP) and then a so called quality prescribing incentive for general practitioners (MCP plus PI).This study is designed as an observational, retrospective, interrupted time series analysis with comparison series (APHS and the rest of Spain) of 46 monthly drug use and sales ratios from January 2001 to October 2004 for each active ingredient in the group of statins.RP has been effective at reducing the volume of sales growth of the off-patent statins, yet its overall impact on sales of all statins has been relatively modest. The quantity and volume of sales impact heavily depends on regulatory RP details such as when the system is introduced, how often it is updated, and how the reference price is calculated.

Minimum distance estimation of stationary and non-stationary ARFIMA processes

A new parametric minimum distance time-domain estimator for ARFIMA processes is introduced in this paper. The proposed estimator minimizes the sum of squared correlations of residuals obtained after filtering a series through ARFIMA parameters. The estimator iseasy to compute and is consistent and asymptotically normally distributed for fractionallyintegrated (FI) processes with an integration order d strictly greater than -0.75. Therefore, it can be applied to both stationary and non-stationary processes. Deterministic components are also allowed in the DGP. Furthermore, as a by-product, the estimation procedure provides an immediate check on the adequacy of the specified model. This is so because the criterion function, when evaluated at the estimated values, coincides with the Box-Pierce goodness of fit statistic. Empirical applications and Monte-Carlo simulations supporting the analytical results and showing the good performance of the estimator in finite samples are also provided.

cAMP-dependent chloride secretion mediates tubule enlargement and cyst formation by cultured mammalian collecting duct cells.

Polycystic kidney diseases result from disruption of the genetically defined program that controls the size and geometry of renal tubules. Cysts which frequently arise from the collecting duct (CD) result from cell proliferation and fluid secretion. From mCCD(cl1) cells, a differentiated mouse CD cell line, we isolated a clonal subpopulation (mCCD-N21) that retains morphogenetic capacity. When grown in three-dimensional gels, mCCD-N21 cells formed highly organized tubular structures consisting of a palisade of polarized epithelial cells surrounding a cylindrical lumen. Subsequent addition of cAMP-elevating agents (forskolin or cholera toxin) or of membrane-permeable cAMP analogs (CPT-cAMP) resulted in rapid and progressive dilatation of existing tubules, leading to the formation of cystlike structures. When grown on filters, mCCD-N21 cells exhibited a high transepithelial resistance as well as aldosterone- and/or vasopressin-induced amiloride-sensitive and -insensitive current. The latter was in part inhibited by Na(+)-K(+)-2Cl(-) cotransporter (bumetanide) and chloride channel (NPPB) inhibitors. Real-time PCR analysis confirmed the expression of NKCC1, the ubiquitous Na(+)-K(+)-2Cl(-) cotransporter and cystic fibrosis transmembrane regulator (CFTR) in mCCD-N21 cells. Tubule enlargement and cyst formation were prevented by inhibitors of Na(+)-K(+)-2Cl(-) cotransporters (bumetanide or ethacrynic acid) or CFTR (NPPB or CFTR inhibitor-172). These results further support the notion that cAMP signaling plays a key role in renal cyst formation, at least in part by promoting chloride-driven fluid secretion. This new in vitro model of tubule-to-cyst conversion affords a unique opportunity for investigating the molecular mechanisms that govern the architecture of epithelial tubes, as well as for dissecting the pathophysiological processes underlying cystic kidney diseases.

Comparative morphology of the spermathecae of some species of Chrysomya Robineau-Desvoidy and Cochliomyia Townsend (Diptera, Calliphoridae)

Comparative morphology of the spermathecae of some species of Chrysomya Robineau-Desvoidy and Cochliomyia Townsend (Diptera, Calliphoridae). Little is known about the morphology of the chitinized structures of the spermathecae of the Calliphoridae. In this work, the spermathecae of Chrysomya albiceps Wiedemann, 1819, C. megacephala Fabricius, 1794, Cochliomyia macellaria Fabricius, 1775 and C. hominivorax Coquerel, 1858 are described and illustrated. The occurrence in one species of four spermathecae, an atypical form for blow flies, was recorded for the first time. The analysis of these structures will allow a better understanding of this group as well as provide taxonomic characters for future phylogenetic studies.

Sampling and TDR probe insertion in the determination of the volumetric soil water content

Volumetric soil water content (theta) can be evaluated in the field by direct or indirect methods. Among the direct, the gravimetric method is regarded as highly reliable and thus often preferred. Its main disadvantages are that sampling and laboratory procedures are labor intensive, and that the method is destructive, which makes resampling of a same point impossible. Recently, the time domain reflectometry (TDR) technique has become a widely used indirect, non-destructive method to evaluate theta. In this study, evaluations of the apparent dielectric number of soils (epsilon) and samplings for the gravimetrical determination of the volumetric soil water content (thetaGrav) were carried out at four sites of a Xanthic Ferralsol in Manaus - Brazil. With the obtained epsilon values, theta was estimated using empirical equations (thetaTDR), and compared with thetaGrav derived from disturbed and undisturbed samples. The main objective of this study was the comparison of thetaTDR estimates of horizontally as well as vertically inserted probes with the thetaGrav values determined by disturbed and undisturbed samples. Results showed that thetaTDR estimates of vertically inserted probes and the average of horizontally measured layers were only slightly and insignificantly different. However, significant differences were found between the thetaTDR estimates of different equations and between disturbed and undisturbed samples in the thetaGrav determinations. The use of the theoretical Knight et al. model, which permits an evaluation of the soil volume assessed by TDR probes, is also discussed. It was concluded that the TDR technique, when properly calibrated, permits in situ, nondestructive measurements of q in Xanthic Ferralsols of similar accuracy as the gravimetric method.

Molecular and Functional Characterization of Neuroblastoma-Initiating Cells

Neuroblastoma (NB) is the most common extracranial malignant tumor in young children and arises at any site of the sympathetic nervous system. The disease exhibits a remarkable phenotypic diversity ranging from spontaneous regression to fatal disease. Poor outcome results from a rapidly progressive, metastatic and drug-resistant disease. Recent studies have suggested that solid tumors may arise from a minor population of cancer stem cells (CSCs) with stem cell markers and typical properties such as self-renewal ability, asymmetric division and drug resistance. In this model, CSCs possess the exclusive ability to initiate and maintain the tumor, and to produce distant metastases. Tumor cell subpopulations with stem-like phenotypes have indeed been identified in several cancer including leukemia, breast, brain and colon cancers. CSC hypothesis still needs to be validated in the other cancers including NB.NB originates from neural crest-derived malignant sympatho-adrenal cells. We have identified rare cells that express markers in conformity with neural crest stem cells and their derived lineages within primary NB tissue and cell lines, leading us to postulate the existence of CSCs in NB tumors.In the absence of specific markers to isolate CSCs, we adapted to NB tumor cells the sphere functional assay, based on the ability of stem cells to grow as spheres in non-adherent conditions. By serial passages of spheres from bone marrow NB metastases, a subset of cells was gradually selected and its specific gene expression profile identified by micro-array time-course analysis. The differentially expressed genes in spheres are enriched in genes implicated in development including CD133, ABC-transporters, WNT and NOTCH genes, identified in others solid cancers as CSCs markers, and other new markers, all referred by us as the Neurosphere Expression Profile (NEP). We confirmed the presence of a cell subpopulation expressing a combination of the NEP markers within a few primary NB samples.The tumorigenic potential of NB spheres was assayed by in vivo tumor growth analyses using orthotopic (adrenal glands) implantations of tumor cells into immune-compromised mice. Tumors derived from the sphere cells were significantly more frequent and were detected earlier compared to whole tumor cells. However, NB cells expressing the neurosphere-associated genes and isolated from the bulk tumors did not recapitulate the CSC-like phenotype in the orthotopic model. In addition, the NB sphere cells lost their higher tumorigenic potential when implanted in a subcutaneous heterotopic in vivo model.These results highlighted the complex behavior of CSC functions and led us to consider the stem-like NB cells as a dynamic and heterogeneous cell population influenced by microenvironment signals.Our approach identified for the first time candidate genes that may be associated with NB self-renewal and tumorigenicity and therefore would establish specific functional targets for more effective therapies in aggressive NB.

Breaking the panels: an application to the GDP per capita

This paper proposes a test statistic for the null hypothesis of panel stationarity that allows for the presence of multiple structural breaks. Two different speci¿cations are considered depending on the structural breaks affecting the individual effects and/or the time trend. The model is ¿exible enough to allow the number of breaks and their position to differ across individuals. The test is shown to have an exact limit distribution with a good ¿nite sample performance. Its application to a typical panel data set of real per capita GDP gives support to the trend stationarity of these series

Identification of CT46/HORMAD1, an immunogenic cancer/testis antigen encoding a putative meiosis-related protein.

Transcripts with ESTs derived exclusively or predominantly from testis, and not from other normal tissues, are likely to be products of genes with testis-restricted expression, and are thus potential cancer/testis (CT) antigen genes. A list of 371 genes with such characteristics was compiled by analyzing publicly available EST databases. RT-PCR analysis of normal and tumor tissues was performed to validate an initial selection of 20 of these genes. Several new CT and CT-like genes were identified. One of these, CT46/HORMAD1, is expressed strongly in testis and weakly in placenta; the highest level of expression in other tissues is <1% of testicular expression. The CT46/HORMAD1 gene was expressed in 31% (34/109) of the carcinomas examined, with 11% (12/109) showing expression levels >10% of the testicular level of expression. CT46/HORMAD1 is a single-copy gene on chromosome 1q21.3, encoding a putative protein of 394 aa. Conserved protein domain analysis identified a HORMA domain involved in chromatin binding. The CT46/HORMAD1 protein was found to be homologous to the prototype HORMA domain-containing protein, Hop1, a yeast meiosis-specific protein, as well as to asy1, a meiotic synaptic mutant protein in Arabidopsis thaliana.

Increased connexin43 expression in human saphenous veins in culture is associated with intimal hyperplasia.

OBJECTIVE: Intimal hyperplasia is a vascular remodelling process that occurs after a vascular injury. The mechanisms involved in intimal hyperplasia are proliferation, dedifferentiation, and migration of medial smooth muscle cells towards the subintimal space. We postulated that gap junctions, which coordinate physiologic processes such as cell growth and differentiation, might participate in the development of intimal hyperplasia. Connexin43 (Cx43) expression levels may be altered in intimal hyperplasia, and we therefore evaluated the regulated expression of Cx43 in human saphenous veins in culture in the presence or not of fluvastatin, an inhibitor of 3-hydroxy-3-methylglutaryl-coenzyme A reductase activity. METHODS: Segments of harvested human saphenous veins, obtained at the time of bypass graft, were opened longitudinally with the luminal surface uppermost and maintained in culture for 14 days. Vein fragments were then processed for histologic examination, neointimal thickness measurements, immunocytochemistry, RNA, and proteins analysis. RESULTS: Of the four connexins (Cx37, 40, 43, and 45), we focused on Cx43 and Cx40, which we found by real-time polymerase chain reaction to be expressed in the saphenous vein because they are the predominant connexins expressed by smooth muscle cells and endothelial cells. After 14 days of culture, histomorphometric analysis showed a significant increase in the intimal thickness as observed during the process of intimal hyperplasia. A time-course analysis revealed a progressive upregulation of Cx43 to reach a maximal increase of sixfold to eightfold at both transcript and protein levels after 14 days in culture. In contrast, the expression of Cx40, abundantly expressed in the endothelial cells, was not altered. Immunofluorescence showed a large increase in Cx43 within smooth muscle cell membranes of the media layer. The development of intimal hyperplasia in vitro was decreased in presence of fluvastatin and was associated with reduced Cx43 expression. CONCLUSIONS: These data show that Cx43 is increased in vitro during the process of intimal hyperplasia and that fluvastatin could prevent this induction, supporting a critical role for Cx43-mediated gap-junctional communication in the human vein during the development of intimal hyperplasia. CLINICAL RELEVANCE: Stenosis due to intimal hyperplasia is the most common cause of failure of venous bypass grafts. To better understand the development of intimal hyperplasia, we used an ex vivo organ culture model to study saphenous veins harvested from patients undergoing a lower limb bypass surgery. In this model, the morphologic and functional integrity of the vessel wall is maintained and significant intimal hyperplasia development occurs after 14 days in culture. We have postulated that gap junctions, which coordinate physiologic processes such as cell growth and differentiation, may participate in the development of intimal hyperplasia. Indeed, intimal hyperplasia consists of proliferation and migration of smooth muscle cells into the subendothelial space. Intercellular communication is responsible for the direct transfer of ions and small molecules from one cell to the other through gap-junction channels found at cell-cell appositions. No study to date has evaluated whether gap junctional communication is involved in the process of intimal hyperplasia in humans. This assertion was investigated by using the aforementioned organ culture model of intimal hyperplasia in human saphenous veins, and our data support a critical role for Cx43-mediated gap junctional communication in human vein during the development of intimal hyperplasia.

Breaking the panels: an application to the GDP per capita

This paper proposes a test statistic for the null hypothesis of panel stationarity that allows for the presence of multiple structural breaks. Two different speci¿cations are considered depending on the structural breaks affecting the individual effects and/or the time trend. The model is ¿exible enough to allow the number of breaks and their position to differ across individuals. The test is shown to have an exact limit distribution with a good ¿nite sample performance. Its application to a typical panel data set of real per capita GDP gives support to the trend stationarity of these series