S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.

Estimation of inelastic seismic material properties of a surgical sea-bed from multi-component marine seismic data

To date, state-of-the-art seismic material parameter estimates from multi-component sea-bed seismic data are based on the assumption that the sea-bed consists of a fully elastic half-space. In reality, however, the shallow sea-bed generally consists of soft, unconsolidated sediments that are characterized by strong to very strong seismic attenuation. To explore the potential implications, we apply a state-of-the-art elastic decomposition algorithm to synthetic data for a range of canonical sea-bed models consisting of a viscoelastic half-space of varying attenuation. We find that in the presence of strong seismic attenuation, as quantified by Q-values of 10 or less, significant errors arise in the conventional elastic estimation of seismic properties. Tests on synthetic data indicate that these errors can be largely avoided by accounting for the inherent attenuation of the seafloor when estimating the seismic parameters. This can be achieved by replacing the real-valued expressions for the elastic moduli in the governing equations in the parameter estimation by their complex-valued viscoelastic equivalents. The practical application of our parameter procedure yields realistic estimates of the elastic seismic material properties of the shallow sea-bed, while the corresponding Q-estimates seem to be biased towards too low values, particularly for S-waves. Given that the estimation of inelastic material parameters is notoriously difficult, particularly in the immediate vicinity of the sea-bed, this is expected to be of interest and importance for civil and ocean engineering purposes.

Estimation of protein coding density in a corpus of DNA sequence data

A number of experimental methods have been reported for estimating the number of genes in a genome, or the closely related coding density of a genome, defined as the fraction of base pairs in codons. Recently, DNA sequence data representative of the genome as a whole have become available for several organisms, making the problem of estimating coding density amenable to sequence analytic methods. Estimates of coding density for a single genome vary widely, so that methods with characterized error bounds have become increasingly desirable. We present a method to estimate the protein coding density in a corpus of DNA sequence data, in which a ‘coding statistic’ is calculated for a large number of windows of the sequence under study, and the distribution of the statistic is decomposed into two normal distributions, assumed to be the distributions of the coding statistic in the coding and noncoding fractions of the sequence windows. The accuracy of the method is evaluated using known data and application is made to the yeast chromosome III sequence and to C.elegans cosmid sequences. It can also be applied to fragmentary data, for example a collection of short sequences determined in the course of STS mapping.

Temporal diffeomorphic free-form deformation: application to motion and strain estimation from 3D echocardiography

This paper presents a new registration algorithm, called Temporal Di eomorphic Free Form Deformation (TDFFD), and its application to motion and strain quanti cation from a sequence of 3D ultrasound (US) images. The originality of our approach resides in enforcing time consistency by representing the 4D velocity eld as the sum of continuous spatiotemporal B-Spline kernels. The spatiotemporal displacement eld is then recovered through forward Eulerian integration of the non-stationary velocity eld. The strain tensor iscomputed locally using the spatial derivatives of the reconstructed displacement eld. The energy functional considered in this paper weighs two terms: the image similarity and a regularization term. The image similarity metric is the sum of squared di erences between the intensities of each frame and a reference one. Any frame in the sequence can be chosen as reference. The regularization term is based on theincompressibility of myocardial tissue. TDFFD was compared to pairwise 3D FFD and 3D+t FFD, bothon displacement and velocity elds, on a set of synthetic 3D US images with di erent noise levels. TDFFDshowed increased robustness to noise compared to these two state-of-the-art algorithms. TDFFD also proved to be more resistant to a reduced temporal resolution when decimating this synthetic sequence. Finally, this synthetic dataset was used to determine optimal settings of the TDFFD algorithm. Subsequently, TDFFDwas applied to a database of cardiac 3D US images of the left ventricle acquired from 9 healthy volunteers and 13 patients treated by Cardiac Resynchronization Therapy (CRT). On healthy cases, uniform strain patterns were observed over all myocardial segments, as physiologically expected. On all CRT patients, theimprovement in synchrony of regional longitudinal strain correlated with CRT clinical outcome as quanti ed by the reduction of end-systolic left ventricular volume at follow-up (6 and 12 months), showing the potential of the proposed algorithm for the assessment of CRT.

Wipe sampling procedure coupled to LC-MS/MS analysis for the simultaneous determination of 10 cytotoxic drugs on different surfaces.

A simple wipe sampling procedure was developed for the surface contamination determination of ten cytotoxic drugs: cytarabine, gemcitabine, methotrexate, etoposide phosphate, cyclophosphamide, ifosfamide, irinotecan, doxorubicin, epirubicin and vincristine. Wiping was performed using Whatman filter paper on different surfaces such as stainless steel, polypropylene, polystyrol, glass, latex gloves, computer mouse and coated paperboard. Wiping and desorption procedures were investigated: The same solution containing 20% acetonitrile and 0.1% formic acid in water gave the best results. After ultrasonic desorption and then centrifugation, samples were analysed by a validated liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) in selected reaction monitoring mode. The whole analytical strategy from wipe sampling to LC-MS/MS analysis was evaluated to determine quantitative performance. The lowest limit of quantification of 10 ng per wiping sample (i.e. 0.1 ng cm(-2)) was determined for the ten investigated cytotoxic drugs. Relative standard deviation for intermediate precision was always inferior to 20%. As recovery was dependent on the tested surface for each drug, a correction factor was determined and applied for real samples. The method was then successfully applied at the cytotoxic production unit of the Geneva University Hospitals pharmacy.

Dynamic estimation of three-dimensional cerebrovascular deformation from rotational angiography

Purpose: The objective of this study is to investigate the feasibility of detecting and quantifying 3D cerebrovascular wall motion from a single 3D rotational x-ray angiography (3DRA) acquisition within a clinically acceptable time and computing from the estimated motion field for the further biomechanical modeling of the cerebrovascular wall. Methods: The whole motion cycle of the cerebral vasculature is modeled using a 4D B-spline transformation, which is estimated from a 4D to 2D + t image registration framework. The registration is performed by optimizing a single similarity metric between the entire 2D + t measured projection sequence and the corresponding forward projections of the deformed volume at their exact time instants. The joint use of two acceleration strategies, together with their implementation on graphics processing units, is also proposed so as to reach computation times close to clinical requirements. For further characterizing vessel wall properties, an approximation of the wall thickness changes is obtained through a strain calculation. Results: Evaluation on in silico and in vitro pulsating phantom aneurysms demonstrated an accurate estimation of wall motion curves. In general, the error was below 10% of the maximum pulsation, even in the situation when substantial inhomogeneous intensity pattern was present. Experiments on in vivo data provided realistic aneurysm and vessel wall motion estimates, whereas in regions where motion was neither visible nor anatomically possible, no motion was detected. The use of the acceleration strategies enabled completing the estimation process for one entire cycle in 5-10 min without degrading the overall performance. The strain map extracted from our motion estimation provided a realistic deformation measure of the vessel wall. Conclusions: The authors' technique has demonstrated that it can provide accurate and robust 4D estimates of cerebrovascular wall motion within a clinically acceptable time, although it has to be applied to a larger patient population prior to possible wide application to routine endovascular procedures. In particular, for the first time, this feasibility study has shown that in vivo cerebrovascular motion can be obtained intraprocedurally from a 3DRA acquisition. Results have also shown the potential of performing strain analysis using this imaging modality, thus making possible for the future modeling of biomechanical properties of the vascular wall.

Wall motion estimation in intracranial aneurysms

The quantification of wall motion in cerebral aneurysms is becoming important owing to its potential connection to rupture, and as a way to incorporate the effects of vascular compliance in computational fluid dynamics (CFD) simulations.Most of papers report values obtained with experimental phantoms, simulated images, or animal models, but the information for real patients is limited. In this paper, we have combined non-rigid registration (IR) with signal processing techniques to measure pulsation in real patients from high frame rate digital subtraction angiography (DSA). We have obtained physiological meaningful waveforms with amplitudes in therange 0mm-0.3mm for a population of 18 patients including ruptured and unruptured aneurysms. Statistically significant differences in pulsation were found according to the rupture status, in agreement with differences in biomechanical properties reported in the literature.

Multiuser detection in a dynamic environment — Part II: Joint user identification and parameter estimation

The problem of jointly estimating the number, the identities, and the data of active users in a time-varying multiuser environment was examined in a companion paper (IEEE Trans. Information Theory, vol. 53, no. 9, September 2007), at whose core was the use of the theory of finite random sets on countable spaces. Here we extend that theory to encompass the more general problem of estimating unknown continuous parameters of the active-user signals. This problem is solved here by applying the theory of random finite sets constructed on hybrid spaces. We doso deriving Bayesian recursions that describe the evolution withtime of a posteriori densities of the unknown parameters and data.Unlike in the above cited paper, wherein one could evaluate theexact multiuser set posterior density, here the continuous-parameter Bayesian recursions do not admit closed-form expressions. To circumvent this difficulty, we develop numerical approximationsfor the receivers that are based on Sequential Monte Carlo (SMC)methods (“particle filtering”). Simulation results, referring to acode-divisin multiple-access (CDMA) system, are presented toillustrate the theory.

Sequential estimation of multipath MIMO-OFDM channels

Wireless “MIMO” systems, employing multiple transmit and receive antennas, promise a significant increase of channel capacity, while orthogonal frequency-division multiplexing (OFDM) is attracting a good deal of attention due to its robustness to multipath fading. Thus, the combination of both techniques is an attractive proposition for radio transmission. The goal of this paper is the description and analysis of a new and novel pilot-aided estimator of multipath block-fading channels. Typical models leading to estimation algorithms assume the number of multipath components and delays to be constant (and often known), while their amplitudes are allowed to vary with time. Our estimator is focused instead on the more realistic assumption that the number of channel taps is also unknown and varies with time following a known probabilistic model. The estimation problem arising from these assumptions is solved using Random-Set Theory (RST), whereby one regards the multipath-channel response as a single set-valued random entity.Within this framework, Bayesian recursive equations determine the evolution with time of the channel estimator. Due to the lack of a closed form for the solution of Bayesian equations, a (Rao–Blackwellized) particle filter (RBPF) implementation ofthe channel estimator is advocated. Since the resulting estimator exhibits a complexity which grows exponentially with the number of multipath components, a simplified version is also introduced. Simulation results describing the performance of our channel estimator demonstrate its effectiveness.

Simultaneous detection of aneuploidies for chromosomes 4, 6, 10 and 17 by automated four colour I-FISH in high hyperdiploid acute lymphoblastic leukemia : diagnostic assessment, clonal heterogeneity and chromosomal instability in adults

SummarySimultaneous detection of aneuploidies for chromosomes 4, 6,10 and 17 by automated four color l-FISH in high hyperdiploid acute lymphoblastic leukemia: diagnostic assessment, clonal heterogeneity and chromosomal instability in adultsAnna Talamo BlandinService de Génétique Médicale, Unité de Cytogénétique du Cancer, CHUVAcute lymphoblastic leukemia (ALL) is a malignant hemopathy characterized by the accumulation of the immature lymphoid cells in the bone marrow and, most often, in the peripheral blood. ALL is a heterogeneous disease with distinct biological and prognostic entities. At diagnosis, cytogenetic and molecular findings constitute important and independent prognostic factors. High hyperdiploidy with 51-67 chromosomes (HeH), one of the largest cytogenetic subsets of ALL, in childhood particularly, is generally associated with a relatively favorable outcome. Chromosome gain is nonrandom, extracopies of some chromosome occurring more frequently than those of others. Concurrent presence of trisomy for chromosomes 4, 10 and 17 confers an especially good prognosis. The first aim of our work was to develop an automated four color interphase fluorescence in situ hybridization (l-FISH) methodology and to assess its ability to detect concurrent aneuploidies 4, 6, 10 and 17 in 10 ALL patients. Various combinations of aneuploidies were identified. All clones detected by conventional cytogenetics were also observed by l-FISH. However, in all patients, l-FISH revealed numerous additional abnormal clones, leading to a high level of clonal heterogeneity. Our second aim has been to investigate the nature and origin of this clonal heterogeneity and to test for the presence of chromosome instability (CIN) in HeH ALL at initial presentation. Ten HeH ALL and 10 non-HeH ALL patients were analysed by four colour l-FISH and numerical CIN values were determined for all four chromosomes together and for each chromosome and patient group, an original approach in ALL. CIN values in HeH ALL proved to be much higher than#iose in non-HeH ALL, suggesting that numerical CIN may be at the origin of the high level of clonal heterogeneity revealed by l-FISH. Our third aim has been to study the evolution of these cytogenetic features during the course of the disease in 10 HeH ALL patients. Clonal heterogeneity was also observed again during disease progression, particularly at relapse. Clones detected at initial presentation generally reappeared in relapse, in most cases with newly generated ones. A significant correlation between the number of abnormal clones and CIN suggested that the higher the instability, the larger the number of abnormal clones. Whereas clonal heterogeneity and its evolution most probably result from underlying chromosome instability, operating processes remain conjectural.RésuméLa leucémie lymphoblastique aiguë (LLA) est une hémopathie maligne qui résulte de l'accumulationde cellules lymphoïdes immatures dans la moelle osseuse, et, le plus souvent, dans le sangpériphérique également. La LLA est une affection hétérogène au sein de laquelle se distinguentplusieurs entités biologiques et pronostiques. Les données cytogénétiques et moléculaires font partieintégrante du diagnostic et jouent un rôle essentiel dans l'évaluation du pronostic. L'hyperdiploïdieélevée à 51-­67 chromosomes (HeH), relativement fréquente, en particulier chez l'enfant, s'associe àun pronostic favorable. Le gain de chromosomes ne relève pas du hasard, certains chromosomesétant plus fréquemment impliqués que d'autres. La présence simultanée des trisomies 4, 6, et 17s'associe à un pronostic particulièrement bon. Le premier but du travail a été de développer uneméthode d'analyse automatique par hybridation in situ fluorescente interphasique (I-­FISH) à 4couleurs et de tester sa capacité à identifier la présence simultanée d'aneuploïdies 4, 6, 10 et 17 dans10 cas de LLA. Différentes combinaisons d'aneuploïdies ont été identifiées. Tous les clones détectéspar cytogénétique conventionnelle l'ont été par I-­FISH. Or, chez tous les patients, l'I-­FISH a révélé denombreux clones anormaux additionnels générant un degré élevé d'hétérogénéité clonale. Notredeuxième but a été d'investiguer la nature et l'origine de cette hétérogénéité et de tester la présenced'instabilité chromosomique (CIN) chez les patients avec une LLA HeH en presentation initiale. DixLLA HeH et 10 LLA non-­HeH ont été analysées par I-­FISH et les valeurs de CIN numérique ont étédéterminées pour les 4 chromosomes ensemble et pour chaque chromosome et groupe de patients,approche originale dans la LLA. Ces valeurs étant beaucoup plus élevées dans la LLA HeH que dansla LLA non-­HeH, elles favorisent l'hypothèse selon laquelle la CIN serait à l'origine de l'hétérogénéitéclonale révélée par I-­FISH. Le troisième but de notre travail a été d'étudier l'évolution de cescaractéristiques cytogénétiques au cours de la maladie dans 10 cas de LLA HeH. L'hétérogénéitéclonale a été retrouvée lors de la progression de la maladie, en particulier en rechute, où les clonesanormaux détectés en présentation initiale réapparaissent, généralement accompagnés de clonesnouveaux. La corrélation existant entre nombre de clones anormaux et valeurs de CIN suggère queplus l'instabilité est élevée, plus le nombre de clones anormaux est grand. Bien que l'hétérogénéitéclonale et son évolution résultent très probablement de l'instabilité chromosomique, les processus àl'oeuvre ne sont pas entièrement élucidés.

Sequential estimation of time-varying multipath channel for MIMO-OFDM systems

In this paper, we introduce a pilot-aided multipath channel estimator for Multiple-Input Multiple-Output (MIMO) Orthogonal Frequency Division Multiplexing (OFDM) systems. Typical estimation algorithms assume the number of multipath components and delays to be known and constant, while theiramplitudes may vary in time. In this work, we focus on the more realistic assumption that also the number of channel taps is unknown and time-varying. The estimation problem arising from this assumption is solved using Random Set Theory (RST), which is a probability theory of finite sets. Due to the lack of a closed form of the optimal filter, a Rao-Blackwellized Particle Filter (RBPF) implementation of the channel estimator is derived. Simulation results demonstrate the estimator effectiveness.

One case, 3 rare simultaneous findings: intramyocardial bronchogenic cyst, P.H558R variant of SCN5A gene, and granular cell tumor of the esophagus.

We describe the sudden death of a 42-year-old white man. The decedent was a healthy young man with a short clinical history of chest pain, fatigue, dizziness, and pyrosis. Two weeks before his death, he underwent medical evaluation for the aforementioned symptoms. Electrocardiogram, chest x-ray, and serum troponin were all within normal limits. Gastroesophageal reflux disease was suspected, and the decedent was treated with omeprazole. Medicolegal autopsy disclosed an incidental intramyocardial bronchogenic cyst and p.H558R variant of the SCN5A gene. The cyst was located between the epicardium and myocardium of the posterior face of the left superior ventricular wall, adjacent to the base of the heart. An incidental granular cell tumor of the esophagus was also identified, which was likely unrelated to death.

Simple estimation of incident HIV infection rates in notification cohorts based on window periods of algorithms for evaluation of line-immunoassay result patterns.

BACKGROUND: Tests for recent infections (TRIs) are important for HIV surveillance. We have shown that a patient's antibody pattern in a confirmatory line immunoassay (Inno-Lia) also yields information on time since infection. We have published algorithms which, with a certain sensitivity and specificity, distinguish between incident (< = 12 months) and older infection. In order to use these algorithms like other TRIs, i.e., based on their windows, we now determined their window periods. METHODS: We classified Inno-Lia results of 527 treatment-naïve patients with HIV-1 infection < = 12 months according to incidence by 25 algorithms. The time after which all infections were ruled older, i.e. the algorithm's window, was determined by linear regression of the proportion ruled incident in dependence of time since infection. Window-based incident infection rates (IIR) were determined utilizing the relationship 'Prevalence = Incidence x Duration' in four annual cohorts of HIV-1 notifications. Results were compared to performance-based IIR also derived from Inno-Lia results, but utilizing the relationship 'incident = true incident + false incident' and also to the IIR derived from the BED incidence assay. RESULTS: Window periods varied between 45.8 and 130.1 days and correlated well with the algorithms' diagnostic sensitivity (R(2) = 0.962; P<0.0001). Among the 25 algorithms, the mean window-based IIR among the 748 notifications of 2005/06 was 0.457 compared to 0.453 obtained for performance-based IIR with a model not correcting for selection bias. Evaluation of BED results using a window of 153 days yielded an IIR of 0.669. Window-based IIR and performance-based IIR increased by 22.4% and respectively 30.6% in 2008, while 2009 and 2010 showed a return to baseline for both methods. CONCLUSIONS: IIR estimations by window- and performance-based evaluations of Inno-Lia algorithm results were similar and can be used together to assess IIR changes between annual HIV notification cohorts.

Technical note: Preliminary estimation of rockfall runout zones

Rockfall propagation areas can be determined using a simple geometric rule known as shadow angle or energy line method based on a simple Coulomb frictional model implemented in the CONEFALL computer program. Runout zones are estimated from a digital terrain model (DTM) and a grid file containing the cells representing rockfall potential source areas. The cells of the DTM that are lowest in altitude and located within a cone centered on a rockfall source cell belong to the potential propagation area associated with that grid cell. In addition, the CONEFALL method allows estimation of mean and maximum velocities and energies of blocks in the rockfall propagation areas. Previous studies indicate that the slope angle cone ranges from 27° to 37° depending on the assumptions made, i.e. slope morphology, probability of reaching a point, maximum run-out, field observations. Different solutions based on previous work and an example of an actual rockfall event are presented here.