A HIDDEN MARKOV MODEL FOR JOINT ESTIMATION OF GENOTYPE AND COPY NUMBER IN HIGH-THROUGHPUT SNP CHIPS

Amplifications and deletions of chromosomal DNA, as well as copy-neutral loss of heterozygosity have been associated with diseases processes. High-throughput single nucleotide polymorphism (SNP) arrays are useful for making genome-wide estimates of copy number and genotype calls. Because neighboring SNPs in high throughput SNP arrays are likely to have dependent copy number and genotype due to the underlying haplotype structure and linkage disequilibrium, hidden Markov models (HMM) may be useful for improving genotype calls and copy number estimates that do not incorporate information from nearby SNPs. We improve previous approaches that utilize a HMM framework for inference in high throughput SNP arrays by integrating copy number, genotype calls, and the corresponding confidence scores when available. Using simulated data, we demonstrate how confidence scores control smoothing in a probabilistic framework. Software for fitting HMMs to SNP array data is available in the R package ICE.

Animal model for percutaneous creation of traumatic thoracic aortic transection

PURPOSE: This study was conducted to create an animal model for thoracic aortic transection that is suitable for thoracic endograft research. MATERIALS AND METHODS: Percutaneous aortic transection creation was attempted in 12 sheep. A custom collapsible circumferential cutting device was inserted into the proximal descending thoracic aorta via a femoral approach with an 11-F delivery catheter. The device was deployed 2 cm distal to the left subclavian artery origin and rotated 20 times to create aortic transection. Aortic diameters, mean aortic pressures, and heart rates were tested for degrees of difference between measurements before and after the creation of transection. On necropsy, the extent of aortic damage was classified as none, nontransmural, or transmural, and aortic transection was classified as none, partial, or circumferential. RESULTS: On angiography, creation of transmural thoracic aortic transection was successful in 91.7% (11/12) of animals. Aortic transection was circumferential in 54.4% (6/11) of animals and partial in 45.6% (5/11) of animals. Mean aortic diameter was 19.6 +/- 3.4 mm (range 12-24 mm) pre-transection and 25.8 +/- 4.5 mm (range 17.8-33 mm) post-transection (P = .0003). Pre-transection, mean aortic pressure was 79 +/- 13.8 mmHg, and 64.6 +/- 15.8 mmHg 15 min post-transection (P = .041). Pre-transection, mean heart rate was 94.5 +/- 17.2 beats per minute (bpm), and 105.8 +/- 17.2 bpm 15 min post-transection (P = .0057). CONCLUSIONS: Thoracic aortic transection was successfully created percutaneously in most animals. The animals remained in hemodynamically stable condition for as long as 240 minutes after the creation of aortic injury. This percutaneous animal model is straightforward and may be of potential value for future thoracic endograft research.

STUDY OF SPARK IGNITION ENGINE COMBUSTION MODEL FOR THE ANALYSIS OF CYCLIC VARIATION AND COMBUSTION STABILITY AT LEAN OPERATING CONDITIONS

A fundamental combustion model for spark-ignition engine is studied in this report. The model is implemented in SIMULINK to simulate engine outputs (mass fraction burn and in-cylinder pressure) under various engine operation conditions. The combustion model includes a turbulent propagation and eddy burning processes based on literature [1]. The turbulence propagation and eddy burning processes are simulated by zero-dimensional method and the flame is assumed as sphere. To predict pressure, temperature and other in-cylinder variables, a two-zone thermodynamic model is used. The predicted results of this model match well with the engine test data under various engine speeds, loads, spark ignition timings and air fuel mass ratios. The developed model is used to study cyclic variation and combustion stability at lean (or diluted) combustion conditions. Several variation sources are introduced into the combustion model to simulate engine performance observed in experimental data. The relations between combustion stability and the introduced variation amount are analyzed at various lean combustion levels.

Bayesian change-point analysis in linear regression model with scale mixtures of normal distributions

In this thesis, we consider Bayesian inference on the detection of variance change-point models with scale mixtures of normal (for short SMN) distributions. This class of distributions is symmetric and thick-tailed and includes as special cases: Gaussian, Student-t, contaminated normal, and slash distributions. The proposed models provide greater flexibility to analyze a lot of practical data, which often show heavy-tail and may not satisfy the normal assumption. As to the Bayesian analysis, we specify some prior distributions for the unknown parameters in the variance change-point models with the SMN distributions. Due to the complexity of the joint posterior distribution, we propose an efficient Gibbs-type with Metropolis- Hastings sampling algorithm for posterior Bayesian inference. Thereafter, following the idea of [1], we consider the problems of the single and multiple change-point detections. The performance of the proposed procedures is illustrated and analyzed by simulation studies. A real application to the closing price data of U.S. stock market has been analyzed for illustrative purposes.

Development of a Mathematical Model for the Calculation of the Tool Appropriation Delay depending on the Tool Inventory

Compliance with punctual delivery under the high pressure of costs can be implemented through the optimization of the in-house tool supply. Within the Transfer Project 13 of the Collaborative Research Centre 489 using the example of the forging industry, a mathematical model was developed which determines the minimum inventory of forging tools required for production, considering the tool appropriation delay.

Development of a novel material flow simulation model for the integration of spatial and process relevant information

Simulation techniques are almost indispensable in the analysis of complex systems. Materials- and related information flow processes in logistics often possess such complexity. Further problem arise as the processes change over time and pose a Big Data problem as well. To cope with these issues adaptive simulations are more and more frequently used. This paper presents a few relevant advanced simulation models and intro-duces a novel model structure, which unifies modelling of geometrical relations and time processes. This way the process structure and their geometric relations can be handled in a well understandable and transparent way. Capabilities and applicability of the model is also presented via a demonstrational example.

Flow-R, a model for susceptibility mapping of debris flows and other gravitational hazards at a regional scale

The development of susceptibility maps for debris flows is of primary importance due to population pressure in hazardous zones. However, hazard assessment by process-based modelling at a regional scale is difficult due to the complex nature of the phenomenon, the variability of local controlling factors, and the uncertainty in modelling parameters. A regional assessment must consider a simplified approach that is not highly parameter dependant and that can provide zonation with minimum data requirements. A distributed empirical model has thus been developed for regional susceptibility assessments using essentially a digital elevation model (DEM). The model is called Flow-R for Flow path assessment of gravitational hazards at a Regional scale (available free of charge under http://www.flow-r.org) and has been successfully applied to different case studies in various countries with variable data quality. It provides a substantial basis for a preliminary susceptibility assessment at a regional scale. The model was also found relevant to assess other natural hazards such as rockfall, snow avalanches and floods. The model allows for automatic source area delineation, given user criteria, and for the assessment of the propagation extent based on various spreading algorithms and simple frictional laws. We developed a new spreading algorithm, an improved version of Holmgren's direction algorithm, that is less sensitive to small variations of the DEM and that is avoiding over-channelization, and so produces more realistic extents. The choices of the datasets and the algorithms are open to the user, which makes it compliant for various applications and dataset availability. Amongst the possible datasets, the DEM is the only one that is really needed for both the source area delineation and the propagation assessment; its quality is of major importance for the results accuracy. We consider a 10 m DEM resolution as a good compromise between processing time and quality of results. However, valuable results have still been obtained on the basis of lower quality DEMs with 25 m resolution.

Papain-induced in vitro disc degeneration model for the study of injectable nucleus pulposus therapy

BACKGROUND CONTEXT Proteolytic enzyme digestion of the intervertebral disc (IVD) offers a method to simulate a condition of disc degeneration for the study of cell-scaffold constructs in the degenerated disc. PURPOSE To characterize an in vitro disc degeneration model (DDM) of different severities of glycosaminoglycans (GAG) and water loss by using papain, and to determine the initial response of the human mesenchymal stem cells (MSCs) introduced into this DDM. STUDY DESIGN Disc degeneration model of a bovine disc explant with an end plate was induced by the injection of papain at various concentrations. Labeled MSCs were later introduced in this model. METHODS Phosphate-buffered saline (PBS control) or papain in various concentrations (3, 15, 30, 60, and 150 U/mL) were injected into the bovine caudal IVD explants. Ten days after the injection, GAG content of the discs was evaluated by dimethylmethylene blue assay and cell viability was determined by live/dead staining together with confocal microscopy. Overall matrix composition was evaluated by histology, and water content was visualized by magnetic resonance imaging. Compressive and torsional stiffness of the DDM were also recorded. In the second part, MSCs were labeled with a fluorescence cell membrane tracker and injected into the nucleus of the DDM or a PBS control. Mesenchymal stem cell viability and distribution were evaluated by confocal microscopy. RESULTS A large drop of GAG and water content of the bovine disc were obtained by injecting >30 U/mL papain. Magnetic resonance imaging showed Grade II, III, and IV disc degeneration by injecting 30, 60, and 150 U/mL papain. A cavity in the center of the disc could facilitate later injection of the nucleus pulposus tissue engineering construct while retaining an intact annulus fibrosus. The remaining disc cell viability was not affected. Mesenchymal stem cells injected into the protease-treated DDM disc showed significantly higher cell viability than when injected into the PBS-injected control disc. CONCLUSIONS By varying the concentration of papain for injection, an increasing amount of GAG and water loss could be induced to simulate the different severities of disc degeneration. MSC suspension introduced into the disc has a very low short-term survival. However, it should be clear that this bovine IVD DDM does not reflect a clinical situation but offers exciting possibilities to test novel tissue engineering protocols.

A papain-induced disc degeneration model for the assessment of thermo-reversible hydrogel-cells therapeutic approach

Nucleus pulposus (NP) regeneration by the application of injectable cell-embedded hydrogels is an appealing approach for tissue engineering. We investigated a thermo-reversible hydrogel (TR-HG), based on a modified polysaccharide with a thermo-reversible polyamide [poly(N-isopropylacrylamide), pNIPAM], which is made to behave as a liquid at room temperature and hardens at > 32 °C. In order to test the hydrogel, a papain-induced bovine caudal disc degeneration model (PDDM), creating a cavity in the NP, was employed. Human mesenchymal stem cells (hMSCs) or autologous bovine NP cells (bNPCs) were seeded in TR-HG; hMSCs were additionally preconditioned with rhGDF-5 for 7 days. Then, TR-HG was reversed to a fluid and the cell suspension injected into the PDDM and kept under static loading for 7 days. Experimental design was: (D1) fresh disc control + PBS injection; (D2) PDDM + PBS injection; (D3) PDDM + TR-HG (material control); (D4) PDDM + TR-HG + bNPCs; (D5) PDDM + TR-HG + hMSCs. Magnetic resonance imaging performed before and after loading, on days 9 and 16, allowed imaging of the hydrogel-filled PDDM and assessment of disc height and volume changes. In gel-injected discs the NP region showed a major drop in volume and disc height during culture under static load. The RT–PCR results of injected hMSCs showed significant upregulation of ACAN, COL2A1, VCAN and SOX9 during culture in the disc cavity, whereas the gene expression profile of NP cells remained unchanged. The cell viability of injected cells (NPCs or hMSCs) was maintained at over 86% in 3D culture and dropped to ~72% after organ culture. Our results underline the need for load-bearing hydrogels that are also cyto-compatible.

Subcutaneous infection model facilitates treatment assessment of secondary Alveolar echinococcosis in mice.

Alveolar echinococcosis (AE) in humans is a parasitic disease characterized by severe damage to the liver and occasionally other organs. AE is caused by infection with the metacestode (larval) stage of the fox tapeworm Echinococcus multilocularis, usually infecting small rodents as natural intermediate hosts. Conventionally, human AE is chemotherapeutically treated with mebendazole or albendazole. There is, however still the need for improved chemotherapeutical options. Primary in vivo studies on drugs of interest are commonly performed in small laboratory animals such as mice and Mongolian jirds, and in most cases, a secondary infection model is used, whereby E. multilocularis metacestodes are directly injected into the peritoneal cavity or into the liver. Disadvantages of this methodological approach include risk of injury to organs during the inoculation and, most notably, a limitation in the macroscopic (visible) assessment of treatment efficacy. Thus, in order to monitor the efficacy of chemotherapeutical treatment, animals have to be euthanized and the parasite tissue dissected. In the present study, mice were infected with E. multilocularis metacestodes through the subcutaneous route and were then subjected to chemotherapy employing albendazole. Serological responses to infection were comparatively assessed in mice infected by the conventional intraperitoneal route. We demonstrate that the subcutaneous infection model for secondary AE facilitates the assessment of the progress of infection and drug treatment in the live animal.