@article{BorongCannistraciConrad2017, author = {Borong, Shao and Cannistraci, Carlo and Conrad, Tim}, title = {Epithelial Mesenchymal Transition Network-based Feature Engineering in Lung Adenocarcinoma Prognosis Prediction Using Multiple Omic Data}, volume = {3}, journal = {Genomics and Computational Biology}, number = {3}, doi = {http://dx.doi.org/10.18547/gcb.2017.vol3.iss3.e57}, pages = {1 -- 13}, year = {2017}, abstract = {Epithelial mesenchymal transition (EMT) process has been shown as highly relevant to cancer prognosis. However, although different biological network-based biomarker identification methods have been proposed to predict cancer prognosis, EMT network has not been directly used for this purpose. In this study, we constructed an EMT regulatory network consisting of 87 molecules and tried to select features that are useful for prognosis prediction in Lung Adenocarcinoma (LUAD). To incorporate multiple molecular profiles, we obtained four types of molecular data including mRNA-Seq, copy number alteration (CNA), DNA methylation, and miRNA-Seq data from The Cancer Genome Atlas. The data were mapped to the EMT network in three alternative ways: mRNA-Seq and miRNA-Seq, DNA methylation, and CNA and miRNA-Seq. Each mapping was employed to extract five different sets of features using discretization and network-based biomarker identification methods. Each feature set was then used to predict prognosis with SVM and logistic regression classifiers. We measured the prediction accuracy with AUC and AUPR values using 10 times 10-fold cross validation. For a more comprehensive evaluation, we also measured the prediction accuracies of clinical features, EMT plus clinical features, randomly picked 87 molecules from each data mapping, and using all molecules from each data type. Counter-intuitively, EMT features do not always outperform randomly selected features and the prediction accuracies of the five feature sets are mostly not significantly different. Clinical features are shown to give the highest prediction accuracies. In addition, the prediction accuracies of both EMT features and random features are comparable as using all features (more than 17,000) from each data type.}, language = {en} } @article{ConradGenzelCvetkovicetal.2017, author = {Conrad, Tim and Genzel, Martin and Cvetkovic, Nada and Wulkow, Niklas and Vybiral, Jan and Kutyniok, Gitta and Sch{\"u}tte, Christof}, title = {Sparse Proteomics Analysis - a compressed sensing-based approach for feature selection and classification of high-dimensional proteomics mass spectrometry data}, volume = {18}, journal = {BMC Bioinformatics}, number = {160}, doi = {10.1186/s12859-017-1565-4}, pages = {1 -- 20}, year = {2017}, abstract = {Motivation: High-throughput proteomics techniques, such as mass spectrometry (MS)-based approaches, produce very high-dimensional data-sets. In a clinical setting one is often interested how MS spectra dier between patients of different classes, for example spectra from healthy patients vs. spectra from patients having a particular disease. Machine learning algorithms are needed to (a) identify these discriminating features and (b) classify unknown spectra based on this feature set. Since the acquired data is usually noisy, the algorithms should be robust to noise and outliers, and the identied feature set should be as small as possible. Results: We present a new algorithm, Sparse Proteomics Analysis (SPA), based on the theory of Compressed Sensing that allows to identify a minimal discriminating set of features from mass spectrometry data-sets. We show how our method performs on artificial and real-world data-sets.}, language = {en} } @article{SeeberConradHoppeetal.2017, author = {Seeber, L. and Conrad, Tim and Hoppe, Christian and Obermeier, Patrick and Chen, X. and Karsch, K. and Muehlhans, S. and Tief, Franziska and Boettcher, Sindy and Diedrich, S. and Schweiger, Brunhilde and Rath, Barbara}, title = {Educating parents about the vaccination status of their children: A user-centered mobile application}, volume = {5}, journal = {Preventive Medicine Reports}, doi = {10.1016/j.pmedr.2017.01.002}, pages = {241 -- 250}, year = {2017}, abstract = {Parents are often uncertain about the vaccination status of their children. In times of vaccine hesitancy, vaccination programs could benefit from active patient participation. The Vaccination App (VAccApp) was developed by the Vienna Vaccine Safety Initiative, enabling parents to learn about the vaccination status of their children, including 25 different routine, special indication and travel vaccines listed in the WHO Immunization Certificate of Vaccination (WHO-ICV). Between 2012 and 2014, the VAccApp was validated in a hospital-based quality management program in Berlin, Germany, in collaboration with the Robert Koch Institute. Parents of 178 children were asked to transfer the immunization data of their children from the WHO-ICV into the VAccApp. The respective WHO-ICV was photocopied for independent, professional data entry (gold standard). Demonstrating the status quo in vaccine information reporting, a Recall Group of 278 parents underwent structured interviews for verbal immunization histories, without the respective WHO-ICV. Only 9\% of the Recall Group were able to provide a complete vaccination status; on average 39\% of the questions were answered correctly. Using the WHO-ICV with the help of the VAccApp resulted in 62\% of parents providing a complete vaccination status; on average 95\% of the questions were answered correctly. After using the VAccApp, parents were more likely to remember key aspects of the vaccination history. User-friendly mobile applications empower parents to take a closer look at the vaccination record, thereby taking an active role in providing accurate vaccination histories. Parents may become motivated to ask informed questions and to keep vaccinations up-to-date.}, language = {en} } @misc{AmbellanLameckervonTycowiczetal.2019, author = {Ambellan, Felix and Lamecker, Hans and von Tycowicz, Christoph and Zachow, Stefan}, title = {Statistical Shape Models - Understanding and Mastering Variation in Anatomy}, issn = {1438-0064}, doi = {10.1007/978-3-030-19385-0_5}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-72699}, year = {2019}, abstract = {In our chapter we are describing how to reconstruct three-dimensional anatomy from medical image data and how to build Statistical 3D Shape Models out of many such reconstructions yielding a new kind of anatomy that not only allows quantitative analysis of anatomical variation but also a visual exploration and educational visualization. Future digital anatomy atlases will not only show a static (average) anatomy but also its normal or pathological variation in three or even four dimensions, hence, illustrating growth and/or disease progression. Statistical Shape Models (SSMs) are geometric models that describe a collection of semantically similar objects in a very compact way. SSMs represent an average shape of many three-dimensional objects as well as their variation in shape. The creation of SSMs requires a correspondence mapping, which can be achieved e.g. by parameterization with a respective sampling. If a corresponding parameterization over all shapes can be established, variation between individual shape characteristics can be mathematically investigated. We will explain what Statistical Shape Models are and how they are constructed. Extensions of Statistical Shape Models will be motivated for articulated coupled structures. In addition to shape also the appearance of objects will be integrated into the concept. Appearance is a visual feature independent of shape that depends on observers or imaging techniques. Typical appearances are for instance the color and intensity of a visual surface of an object under particular lighting conditions, or measurements of material properties with computed tomography (CT) or magnetic resonance imaging (MRI). A combination of (articulated) statistical shape models with statistical models of appearance lead to articulated Statistical Shape and Appearance Models (a-SSAMs).After giving various examples of SSMs for human organs, skeletal structures, faces, and bodies, we will shortly describe clinical applications where such models have been successfully employed. Statistical Shape Models are the foundation for the analysis of anatomical cohort data, where characteristic shapes are correlated to demographic or epidemiologic data. SSMs consisting of several thousands of objects offer, in combination with statistical methods ormachine learning techniques, the possibility to identify characteristic clusters, thus being the foundation for advanced diagnostic disease scoring.}, language = {en} } @misc{AmbellanTackEhlkeetal.2019, author = {Ambellan, Felix and Tack, Alexander and Ehlke, Moritz and Zachow, Stefan}, title = {Automated Segmentation of Knee Bone and Cartilage combining Statistical Shape Knowledge and Convolutional Neural Networks: Data from the Osteoarthritis Initiative}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-72704}, year = {2019}, abstract = {We present a method for the automated segmentation of knee bones and cartilage from magnetic resonance imaging (MRI) that combines a priori knowledge of anatomical shape with Convolutional Neural Networks (CNNs).The proposed approach incorporates 3D Statistical Shape Models (SSMs) as well as 2D and 3D CNNs to achieve a robust and accurate segmentation of even highly pathological knee structures.The shape models and neural networks employed are trained using data from the Osteoarthritis Initiative (OAI) and the MICCAI grand challenge "Segmentation of Knee Images 2010" (SKI10), respectively. We evaluate our method on 40 validation and 50 submission datasets from the SKI10 challenge.For the first time, an accuracy equivalent to the inter-observer variability of human readers is achieved in this challenge.Moreover, the quality of the proposed method is thoroughly assessed using various measures for data from the OAI, i.e. 507 manual segmentations of bone and cartilage, and 88 additional manual segmentations of cartilage. Our method yields sub-voxel accuracy for both OAI datasets. We make the 507 manual segmentations as well as our experimental setup publicly available to further aid research in the field of medical image segmentation.In conclusion, combining localized classification via CNNs with statistical anatomical knowledge via SSMs results in a state-of-the-art segmentation method for knee bones and cartilage from MRI data.}, language = {en} } @article{ObermeierHeimBiereetal.2018, author = {Obermeier, Patrick and Heim, A. and Biere, Barbara and Hage, E. and Alchikh, Maren and Conrad, Tim and Schweiger, Brunhilde and Rath, Barbara}, title = {Clinical characteristics and disease severity associated with adenovirus infections in infants and children - discovery of a novel adenovirus, HAdV-D80}, journal = {Clinical Infectious Diseases}, year = {2018}, language = {en} } @article{KorkSpiesConradetal.2018, author = {Kork, F. and Spies, Claudia and Conrad, Tim and Weiss, B. and Roenneberg, T. and Wernecke, K.-D. and Balzer, Felix}, title = {Associations of postoperative mortality with the time of day, week and year}, journal = {Anaesthesia}, year = {2018}, language = {en} } @article{GelssKlusEisertetal.2019, author = {Gelß, Patrick and Klus, Stefan and Eisert, Jens and Sch{\"u}tte, Christof}, title = {Multidimensional Approximation of Nonlinear Dynamical Systems}, volume = {14}, journal = {Journal of Computational and Nonlinear Dynamics}, number = {6}, doi = {10.1115/1.4043148}, year = {2019}, abstract = {A key task in the field of modeling and analyzing nonlinear dynamical systems is the recovery of unknown governing equations from measurement data only. There is a wide range of application areas for this important instance of system identification, ranging from industrial engineering and acoustic signal processing to stock market models. In order to find appropriate representations of underlying dynamical systems, various data-driven methods have been proposed by different communities. However, if the given data sets are high-dimensional, then these methods typically suffer from the curse of dimensionality. To significantly reduce the computational costs and storage consumption, we propose the method multidimensional approximation of nonlinear dynamical systems (MANDy) which combines data-driven methods with tensor network decompositions. The efficiency of the introduced approach will be illustrated with the aid of several high-dimensional nonlinear dynamical systems.}, language = {en} } @incollection{AmbellanLameckervonTycowiczetal.2019, author = {Ambellan, Felix and Lamecker, Hans and von Tycowicz, Christoph and Zachow, Stefan}, title = {Statistical Shape Models - Understanding and Mastering Variation in Anatomy}, volume = {3}, booktitle = {Biomedical Visualisation}, number = {1156}, editor = {Rea, Paul M.}, edition = {1}, publisher = {Springer Nature Switzerland AG}, isbn = {978-3-030-19384-3}, doi = {10.1007/978-3-030-19385-0_5}, pages = {67 -- 84}, year = {2019}, abstract = {In our chapter we are describing how to reconstruct three-dimensional anatomy from medical image data and how to build Statistical 3D Shape Models out of many such reconstructions yielding a new kind of anatomy that not only allows quantitative analysis of anatomical variation but also a visual exploration and educational visualization. Future digital anatomy atlases will not only show a static (average) anatomy but also its normal or pathological variation in three or even four dimensions, hence, illustrating growth and/or disease progression. Statistical Shape Models (SSMs) are geometric models that describe a collection of semantically similar objects in a very compact way. SSMs represent an average shape of many three-dimensional objects as well as their variation in shape. The creation of SSMs requires a correspondence mapping, which can be achieved e.g. by parameterization with a respective sampling. If a corresponding parameterization over all shapes can be established, variation between individual shape characteristics can be mathematically investigated. We will explain what Statistical Shape Models are and how they are constructed. Extensions of Statistical Shape Models will be motivated for articulated coupled structures. In addition to shape also the appearance of objects will be integrated into the concept. Appearance is a visual feature independent of shape that depends on observers or imaging techniques. Typical appearances are for instance the color and intensity of a visual surface of an object under particular lighting conditions, or measurements of material properties with computed tomography (CT) or magnetic resonance imaging (MRI). A combination of (articulated) statistical shape models with statistical models of appearance lead to articulated Statistical Shape and Appearance Models (a-SSAMs).After giving various examples of SSMs for human organs, skeletal structures, faces, and bodies, we will shortly describe clinical applications where such models have been successfully employed. Statistical Shape Models are the foundation for the analysis of anatomical cohort data, where characteristic shapes are correlated to demographic or epidemiologic data. SSMs consisting of several thousands of objects offer, in combination with statistical methods ormachine learning techniques, the possibility to identify characteristic clusters, thus being the foundation for advanced diagnostic disease scoring.}, language = {en} } @article{HildebrandtBrueningSchmidtetal.2019, author = {Hildebrandt, Thomas and Bruening, Jan Joris and Schmidt, Nora Laura and Lamecker, Hans and Heppt, Werner and Zachow, Stefan and Goubergrits, Leonid}, title = {The Healthy Nasal Cavity - Characteristics of Morphology and Related Airflow Based on a Statistical Shape Model Viewed from a Surgeon's Perspective}, volume = {35}, journal = {Facial Plastic Surgery}, number = {1}, doi = {10.1055/s-0039-1677721}, pages = {9 -- 13}, year = {2019}, abstract = {Functional surgery on the nasal framework requires referential criteria to objectively assess nasal breathing for indication and follow-up. Thismotivated us to generate amean geometry of the nasal cavity based on a statistical shape model. In this study, the authors could demonstrate that the introduced nasal cavity's mean geometry features characteristics of the inner shape and airflow, which are commonly observed in symptom-free subjects. Therefore, the mean geometry might serve as a reference-like model when one considers qualitative aspects. However, to facilitate quantitative considerations and statistical inference, further research is necessary. Additionally, the authorswere able to obtain details about the importance of the isthmus nasi and the inferior turbinate for the intranasal airstream.}, language = {en} } @article{HildebrandtBrueningLameckeretal.2019, author = {Hildebrandt, Thomas and Bruening, Jan Joris and Lamecker, Hans and Zachow, Stefan and Heppt, Werner and Schmidt, Nora and Goubergrits, Leonid}, title = {Digital Analysis of Nasal Airflow Facilitating Decision Support in Rhinosurgery}, volume = {35}, journal = {Facial Plastic Surgery}, number = {1}, doi = {10.1055/s-0039-1677720}, pages = {1 -- 8}, year = {2019}, abstract = {Successful functional surgery on the nasal framework requires reliable and comprehensive diagnosis. In this regard, the authors introduce a new methodology: Digital Analysis of Nasal Airflow (diANA). It is based on computational fluid dynamics, a statistical shape model of the healthy nasal cavity and rhinologic expertise. diANA necessitates an anonymized tomographic dataset of the paranasal sinuses including the complete nasal cavity and, when available, clinical information. The principle of diANA is to compare the morphology and the respective airflow of an individual nose with those of a reference. This enablesmorphometric aberrations and consecutive flow field anomalies to localize and quantify within a patient's nasal cavity. Finally, an elaborated expert opinion with instructive visualizations is provided. Using diANA might support surgeons in decision-making, avoiding unnecessary surgery, gaining more precision, and target-orientation for indicated operations.}, language = {en} } @article{HettichSchierjottSchillingetal.2018, author = {Hettich, G. and Schierjott, R. A. and Schilling, C. and Maas, A. and Ramm, Heiko and Bindernagel, Matthias and Lamecker, Hans and Grupp, T. M.}, title = {Validation of a Statistical Shape Model for Acetabular Bone Defect Analysis}, journal = {ISTA 2018 London Abstract Book}, year = {2018}, abstract = {Acetabular bone defects are still challenging to quantify. Numerous classification schemes have been proposed to categorize the diverse kinds of defects. However, these classification schemes are mainly descriptive and hence it remains difficult to apply them in pre-clinical testing, implant development and pre-operative planning. By reconstructing the native situation of a defect pelvis using a Statistical Shape Model (SSM), a more quantitative analysis of the bone defects could be performed. The aim of this study is to develop such a SSM and to validate its accuracy using relevant clinical scenarios and parameters.}, language = {en} } @article{MirelesConrad2018, author = {Mireles, Victor and Conrad, Tim}, title = {Reusable building blocks in biological systems}, volume = {15}, journal = {Journal of the Royal Society Interface}, number = {149}, doi = {10.1098/rsif.2018.0595}, pages = {1 -- 10}, year = {2018}, abstract = {One of the most widely recognized features of biological systems is their modularity. The modules that constitute biological systems are said to be redeployed and combined across several conditions, thus acting as building blocks. In this work, we analyse to what extent are these building blocks reusable as compared with those found in randomized versions of a system. We develop a notion of decompositions of systems into phenotypic building blocks, which allows them to overlap while maximizing the number of times a building block is reused across several conditions. Different biological systems present building blocks whose reusability ranges from single use (e.g. condition specific) to constitutive, although their average reusability is not always higher than random equivalents of the system. These decompositions reveal a distinct distribution of building block sizes in real biological systems. This distribution stems, in part, from the peculiar usage pattern of the elements of biological systems, and constitutes a new angle to study the evolution of modularity.}, language = {en} } @article{MaConradAlchikhetal.2018, author = {Ma, Xiaolin and Conrad, Tim and Alchikh, Maren and Reiche, J. and Schweiger, Brunhilde and Rath, Barbara}, title = {Can we distinguish respiratory viral infections based on clinical features? A prospective pediatric cohort compared to systematic literature review}, volume = {28}, journal = {Medical Virology}, number = {5}, doi = {10.1002/rmv.1997}, pages = {1 -- 12}, year = {2018}, abstract = {Studies have shown that the predictive value of "clinical diagnoses" of influenza and other respiratory viral infections is low, especially in children. In routine care, pediatricians often resort to clinical diagnoses, even in the absence of robust evidence-based criteria. We used a dual approach to identify clinical characteristics that may help to differentiate infections with common pathogens including influenza, respiratory syncytial virus, adenovirus, metapneumovirus, rhinovirus, bocavirus-1, coronaviruses, or parainfluenza virus: (a) systematic review and meta-analysis of 47 clinical studies published in Medline (June 1996 to March 2017, PROSPERO registration number: CRD42017059557) comprising 49 858 individuals and (b) data-driven analysis of an inception cohort of 6073 children with ILI (aged 0-18 years, 56\% male, December 2009 to March 2015) examined at the point of care in addition to blinded PCR testing. We determined pooled odds ratios for the literature analysis and compared these to odds ratios based on the clinical cohort dataset. This combined analysis suggested significant associations between influenza and fever or headache, as well as between respiratory syncytial virus infection and cough, dyspnea, and wheezing. Similarly, literature and cohort data agreed on significant associations between HMPV infection and cough, as well as adenovirus infection and fever. Importantly, none of the abovementioned features were unique to any particular pathogen but were also observed in association with other respiratory viruses. In summary, our "real-world" dataset confirmed published literature trends, but no individual feature allows any particular type of viral infection to be ruled in or ruled out. For the time being, laboratory confirmation remains essential. More research is needed to develop scientifically validated decision models to inform best practice guidelines and targeted diagnostic algorithms.}, language = {en} } @article{AlchikhConradHoppeetal.2018, author = {Alchikh, Maren and Conrad, Tim and Hoppe, Christian and Ma, Xiaolin and Broberg, Eeva K. and Penttinen, P. and Reiche, J. and Biere, Barbara and Schweiger, Brunhilde and Rath, Barbara}, title = {Are we missing respiratory viral infections in infants and children? Comparison of a hospital-based quality management system with standard of care. Clinical Microbiology and Infection}, journal = {Clinical Microbiology and Infection}, number = {06/18}, doi = {10.1016/j.cmi.2018.05.023}, pages = {1 -- 1}, year = {2018}, language = {en} } @inproceedings{JayrannejadConrad2017, author = {Jayrannejad, Fahrnaz and Conrad, Tim}, title = {Better Interpretable Models for Proteomics Data Analysis Using rule-based Mining}, booktitle = {Springer Lecture Notes in Artificial Intelligence}, pages = {studi}, year = {2017}, abstract = {Recent advances in -omics technology has yielded in large data-sets in many areas of biology, such as mass spectrometry based proteomics. However, analyzing this data is still a challenging task mainly due to the very high dimensionality and high noise content of the data. One of the main objectives of the analysis is the identification of relevant patterns (or features) which can be used for classification of new samples to healthy or diseased. So, a method is required to find easily interpretable models from this data. To gain the above mentioned goal, we have adapted the disjunctive association rule mining algorithm, TitanicOR, to identify emerging patterns from our mass spectrometry proteomics data-sets. Comparison to five state-of-the-art methods shows that our method is advantageous them in terms of identifying the inter-dependency between the features and the TP-rate and precision of the features selected. We further demonstrate the applicability of our algorithm to one previously published clinical data-set.}, language = {en} } @inproceedings{TackZachow2019, author = {Tack, Alexander and Zachow, Stefan}, title = {Accurate Automated Volumetry of Cartilage of the Knee using Convolutional Neural Networks: Data from the Osteoarthritis Initiative}, booktitle = {IEEE 16th International Symposium on Biomedical Imaging (ISBI 2019)}, doi = {10.1109/ISBI.2019.8759201}, pages = {40 -- 43}, year = {2019}, abstract = {Volumetry of cartilage of the knee is needed for knee osteoarthritis (KOA) assessment. It is typically performed manually in a tedious and subjective process. We developed a method for an automated, segmentation-based quantification of cartilage volume by employing 3D Convolutional Neural Networks (CNNs). CNNs were trained in a supervised manner using magnetic resonance imaging data and cartilage volumetry readings performed by clinical experts for 1378 subjects provided by the Osteoarthritis Initiative. It was shown that 3D CNNs are able to achieve volume measures comparable to the magnitude of variation between expert readings and the real in vivo situation. In the future, accurate automated cartilage volumetry might support both, diagnosis of KOA as well as longitudinal analysis of KOA progression.}, language = {en} } @misc{AmbellanTackEhlkeetal.2019, author = {Ambellan, Felix and Tack, Alexander and Ehlke, Moritz and Zachow, Stefan}, title = {Automated Segmentation of Knee Bone and Cartilage combining Statistical Shape Knowledge and Convolutional Neural Networks: Data from the Osteoarthritis Initiative (Supplementary Material)}, volume = {52}, journal = {Medical Image Analysis}, number = {2}, doi = {10.12752/4.ATEZ.1.0}, pages = {109 -- 118}, year = {2019}, abstract = {We present a method for the automated segmentation of knee bones and cartilage from magnetic resonance imaging that combines a priori knowledge of anatomical shape with Convolutional Neural Networks (CNNs). The proposed approach incorporates 3D Statistical Shape Models (SSMs) as well as 2D and 3D CNNs to achieve a robust and accurate segmentation of even highly pathological knee structures. The shape models and neural networks employed are trained using data of the Osteoarthritis Initiative (OAI) and the MICCAI grand challenge "Segmentation of Knee Images 2010" (SKI10), respectively. We evaluate our method on 40 validation and 50 submission datasets of the SKI10 challenge. For the first time, an accuracy equivalent to the inter-observer variability of human readers has been achieved in this challenge. Moreover, the quality of the proposed method is thoroughly assessed using various measures for data from the OAI, i.e. 507 manual segmentations of bone and cartilage, and 88 additional manual segmentations of cartilage. Our method yields sub-voxel accuracy for both OAI datasets. We made the 507 manual segmentations as well as our experimental setup publicly available to further aid research in the field of medical image segmentation. In conclusion, combining statistical anatomical knowledge via SSMs with the localized classification via CNNs results in a state-of-the-art segmentation method for knee bones and cartilage from MRI data.}, language = {en} } @article{KoltaiSchuette2018, author = {Koltai, P{\´e}ter and Sch{\"u}tte, Christof}, title = {A multiscale perturbation expansion approach for Markov state modeling of non-stationary molecular dynamics}, volume = {16}, journal = {SIAM J. Multiscale Model. Simul.}, number = {4}, publisher = {SIAM}, doi = {10.1137/17M1146403}, pages = {1455 -- 1485}, year = {2018}, abstract = {We investigate metastable dynamical systems subject to non-stationary forcing as they appear in molecular dynamics for systems driven by external fields. We show, that if the strength of the forcing is inversely proportional to the length of the slow metastable time scales of the unforced system, then the effective behavior of the forced system on slow time scales can be described by a low-dimensional reduced master equation. Our construction is explicit and uses the multiscale perturbation expansion method called two-timing, or method of multiple scales. The reduced master equation—a Markov state model—can be assembled by constructing two equilibrium Markov state models; one for the unforced system, and one for a slightly perturbed one.}, language = {en} } @article{AmbellanTackEhlkeetal.2019, author = {Ambellan, Felix and Tack, Alexander and Ehlke, Moritz and Zachow, Stefan}, title = {Automated Segmentation of Knee Bone and Cartilage combining Statistical Shape Knowledge and Convolutional Neural Networks: Data from the Osteoarthritis Initiative}, volume = {52}, journal = {Medical Image Analysis}, number = {2}, doi = {10.1016/j.media.2018.11.009}, pages = {109 -- 118}, year = {2019}, abstract = {We present a method for the automated segmentation of knee bones and cartilage from magnetic resonance imaging that combines a priori knowledge of anatomical shape with Convolutional Neural Networks (CNNs). The proposed approach incorporates 3D Statistical Shape Models (SSMs) as well as 2D and 3D CNNs to achieve a robust and accurate segmentation of even highly pathological knee structures. The shape models and neural networks employed are trained using data of the Osteoarthritis Initiative (OAI) and the MICCAI grand challenge "Segmentation of Knee Images 2010" (SKI10), respectively. We evaluate our method on 40 validation and 50 submission datasets of the SKI10 challenge. For the first time, an accuracy equivalent to the inter-observer variability of human readers has been achieved in this challenge. Moreover, the quality of the proposed method is thoroughly assessed using various measures for data from the OAI, i.e. 507 manual segmentations of bone and cartilage, and 88 additional manual segmentations of cartilage. Our method yields sub-voxel accuracy for both OAI datasets. We made the 507 manual segmentations as well as our experimental setup publicly available to further aid research in the field of medical image segmentation. In conclusion, combining statistical anatomical knowledge via SSMs with the localized classification via CNNs results in a state-of-the-art segmentation method for knee bones and cartilage from MRI data.}, language = {en} } @article{ShaoBjaanaesHellandetal.2019, author = {Shao, Borong and Bjaanaes, Maria and Helland, Aslaug and Sch{\"u}tte, Christof and Conrad, Tim}, title = {EMT network-based feature selection improves prognosis prediction in lung adenocarcinoma}, volume = {14}, journal = {PLOS ONE}, number = {1}, doi = {10.1371/journal.pone.0204186}, year = {2019}, abstract = {Various feature selection algorithms have been proposed to identify cancer prognostic biomarkers. In recent years, however, their reproducibility is criticized. The performance of feature selection algorithms is shown to be affected by the datasets, underlying networks and evaluation metrics. One of the causes is the curse of dimensionality, which makes it hard to select the features that generalize well on independent data. Even the integration of biological networks does not mitigate this issue because the networks are large and many of their components are not relevant for the phenotype of interest. With the availability of multi-omics data, integrative approaches are being developed to build more robust predictive models. In this scenario, the higher data dimensions create greater challenges. We proposed a phenotype relevant network-based feature selection (PRNFS) framework and demonstrated its advantages in lung cancer prognosis prediction. We constructed cancer prognosis relevant networks based on epithelial mesenchymal transition (EMT) and integrated them with different types of omics data for feature selection. With less than 2.5\% of the total dimensionality, we obtained EMT prognostic signatures that achieved remarkable prediction performance (average AUC values above 0.8), very significant sample stratifications, and meaningful biological interpretations. In addition to finding EMT signatures from different omics data levels, we combined these single-omics signatures into multi-omics signatures, which improved sample stratifications significantly. Both single- and multi-omics EMT signatures were tested on independent multi-omics lung cancer datasets and significant sample stratifications were obtained.}, language = {en} } @article{GuptaPeterJungetal.2019, author = {Gupta, Pooja and Peter, Sarah and Jung, Markus and Lewin, Astrid and Hemmrich-Stanisak, Georg and Franke, Andre and von Kleist, Max and Sch{\"u}tte, Christof and Einspanier, Ralf and Sharbati, Soroush and zur Bruegge, Jennifer}, title = {Analysis of long non-coding RNA and mRNA expression in bovine macrophages brings up novel 2 aspects of Mycobacterium avium subspecies paratuberculosis infections}, volume = {9}, journal = {Scientific Reports in Nature}, doi = {10.1038/s41598-018-38141-x}, year = {2019}, abstract = {Paratuberculosis is a major disease in cattle that severely affects animal welfare and causes huge economic losses worldwide. Development of alternative diagnostic methods is of urgent need to control the disease. Recent studies suggest that long non-coding RNAs (lncRNAs) play a crucial role in regulating immune function and may confer valuable information about the disease. However, their role has not yet been investigated in cattle with respect to infection towards Paratuberculosis. Therefore, we investigated the alteration in genomic expression profiles of mRNA and lncRNA in bovine macrophages in response to Paratuberculosis infection using RNA-Seq. We identified 397 potentially novel lncRNA candidates in macrophages of which 38 were differentially regulated by the infection. A total of 820 coding genes were also significantly altered by the infection. Co-expression analysis of lncRNAs and their neighbouring coding genes suggest regulatory functions of lncRNAs in pathways related to immune response. For example, this included protein coding genes such as TNIP3, TNFAIP3 and NF-κB2 that play a role in NF-κB2 signalling, a pathway associated with immune response. This study advances our understanding of lncRNA roles during Paratuberculosis infection.}, language = {en} } @article{BernardSalamancaThunbergetal.2017, author = {Bernard, Florian and Salamanca, Luis and Thunberg, Johan and Tack, Alexander and Jentsch, Dennis and Lamecker, Hans and Zachow, Stefan and Hertel, Frank and Goncalves, Jorge and Gemmar, Peter}, title = {Shape-aware Surface Reconstruction from Sparse 3D Point-Clouds}, volume = {38}, journal = {Medical Image Analysis}, doi = {10.1016/j.media.2017.02.005}, pages = {77 -- 89}, year = {2017}, abstract = {The reconstruction of an object's shape or surface from a set of 3D points plays an important role in medical image analysis, e.g. in anatomy reconstruction from tomographic measurements or in the process of aligning intra-operative navigation and preoperative planning data. In such scenarios, one usually has to deal with sparse data, which significantly aggravates the problem of reconstruction. However, medical applications often provide contextual information about the 3D point data that allow to incorporate prior knowledge about the shape that is to be reconstructed. To this end, we propose the use of a statistical shape model (SSM) as a prior for surface reconstruction. The SSM is represented by a point distribution model (PDM), which is associated with a surface mesh. Using the shape distribution that is modelled by the PDM, we formulate the problem of surface reconstruction from a probabilistic perspective based on a Gaussian Mixture Model (GMM). In order to do so, the given points are interpreted as samples of the GMM. By using mixture components with anisotropic covariances that are "oriented" according to the surface normals at the PDM points, a surface-based fitting is accomplished. Estimating the parameters of the GMM in a maximum a posteriori manner yields the reconstruction of the surface from the given data points. We compare our method to the extensively used Iterative Closest Points method on several different anatomical datasets/SSMs (brain, femur, tibia, hip, liver) and demonstrate superior accuracy and robustness on sparse data.}, language = {en} } @article{WinkelmannSchuette2016, author = {Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {The Spatiotemporal Master Equation: Approximation of Reaction-Diffusion Dynamics via Markov State Modeling}, volume = {145}, journal = {Journal of Chemical Physics}, number = {21}, doi = {10.1063/1.4971163}, year = {2016}, abstract = {Accurate modeling and numerical simulation of reaction kinetics is a topic of steady interest.We consider the spatiotemporal chemical master equation (ST-CME) as a model for stochastic reaction-diffusion systems that exhibit properties of metastability. The space of motion is decomposed into metastable compartments and diffusive motion is approximated by jumps between these compartments. Treating these jumps as first-order reactions, simulation of the resulting stochastic system is possible by the Gillespie method. We present the theory of Markov state models (MSM) as a theoretical foundation of this intuitive approach. By means of Markov state modeling, both the number and shape of compartments and the transition rates between them can be determined. We consider the ST-CME for two reaction-diffusion systems and compare it to more detailed models. Moreover, a rigorous formal justification of the ST-CME by Galerkin projection methods is presented.}, language = {en} } @article{DuwalWinkelmannSchuetteetal.2015, author = {Duwal, Sulav and Winkelmann, Stefanie and Sch{\"u}tte, Christof and von Kleist, Max}, title = {Optimal Treatment Strategies in the Context of 'Treatment for Prevention' against HIV/1 in Resource-Poor Settings}, volume = {11}, journal = {PloS Computational Biology}, number = {4}, doi = {10.1371/journal.pcbi.1004200}, year = {2015}, abstract = {An estimated 2.7 million new HIV-1 infections occurred in 2010. `Treatment-for-prevention' may strongly prevent HIV-1 transmission. The basic idea is that immediate treatment initiation rapidly decreases virus burden, which reduces the number of transmittable viruses and thereby the probability of infection. However, HIV inevitably develops drug resistance, which leads to virus rebound and nullifies the effect of `treatment-for-prevention' for the time it remains unrecognized. While timely conducted treatment changes may avert periods of viral rebound, necessary treatment options and diagnostics may be lacking in resource-constrained settings. Within this work, we provide a mathematical platform for comparing different treatment paradigms that can be applied to many medical phenomena. We use this platform to optimize two distinct approaches for the treatment of HIV-1: (i) a diagnostic-guided treatment strategy, based on infrequent and patient-specific diagnostic schedules and (ii) a pro-active strategy that allows treatment adaptation prior to diagnostic ascertainment. Both strategies are compared to current clinical protocols (standard of care and the HPTN052 protocol) in terms of patient health, economic means and reduction in HIV-1 onward transmission exemplarily for South Africa. All therapeutic strategies are assessed using a coarse-grained stochastic model of within-host HIV dynamics and pseudo-codes for solving the respective optimal control problems are provided. Our mathematical model suggests that both optimal strategies (i)-(ii) perform better than the current clinical protocols and no treatment in terms of economic means, life prolongation and reduction of HIV-transmission. The optimal diagnostic-guided strategy suggests rare diagnostics and performs similar to the optimal pro-active strategy. Our results suggest that 'treatment-for-prevention' may be further improved using either of the two analyzed treatment paradigms.}, language = {en} } @article{SarichDjurdjevacConradBruckneretal.2014, author = {Sarich, Marco and Djurdjevac Conrad, Natasa and Bruckner, Sharon and Conrad, Tim and Sch{\"u}tte, Christof}, title = {Modularity revisited: A novel dynamics-based concept for decomposing complex networks}, volume = {1}, journal = {Journal of Computational Dynamics}, number = {1}, doi = {10.3934/jcd.2014.1.191}, pages = {191 -- 212}, year = {2014}, language = {en} } @article{AkbariShandizBoulosSavarssonetal.2018, author = {Akbari Shandiz, Mohsen and Boulos, Paul and S{\ae}varsson, Stefan and Ramm, Heiko and Fu, Chun Kit and Miller, Stephen and Zachow, Stefan and Anglin, Carolyn}, title = {Changes in Knee Shape and Geometry Resulting from Total Knee Arthroplasty}, volume = {232}, journal = {Journal of Engineering in Medicine}, number = {1}, doi = {10.1177/0954411917743274}, pages = {67 -- 79}, year = {2018}, abstract = {Changes in knee shape and geometry resulting from total knee arthroplasty can affect patients in numerous important ways: pain, function, stability, range of motion, and kinematics. Quantitative data concerning these changes have not been previously available, to our knowledge, yet are essential to understand individual experiences of total knee arthroplasty and thereby improve outcomes for all patients. The limiting factor has been the challenge of accurately measuring these changes. Our study objective was to develop a conceptual framework and analysis method to investigate changes in knee shape and geometry, and prospectively apply it to a sample total knee arthroplasty population. Using clinically available computed tomography and radiography imaging systems, the three-dimensional knee shape and geometry of nine patients (eight varus and one valgus) were compared before and after total knee arthroplasty. All patients had largely good outcomes after their total knee arthroplasty. Knee shape changed both visually and numerically. On average, the distal condyles were slightly higher medially and lower laterally (range: +4.5 mm to -4.4 mm), the posterior condyles extended farther out medially but not laterally (range: +1.8 to -6.4 mm), patellofemoral distance increased throughout flexion by 1.8-3.5 mm, and patellar thickness alone increased by 2.9 mm (range: 0.7-5.2 mm). External femoral rotation differed preop and postop. Joint line distance, taking cartilage into account, changed by +0.7 to -1.5 mm on average throughout flexion. Important differences in shape and geometry were seen between pre-total knee arthroplasty and post-total knee arthroplasty knees. While this is qualitatively known, this is the first study to report it quantitatively, an important precursor to identifying the reasons for the poor outcome of some patients. Using the developed protocol and visualization techniques to compare patients with good versus poor clinical outcomes could lead to changes in implant design, implant selection, component positioning, and surgical technique. Recommendations based on this sample population are provided. Intraoperative and postoperative feedback could ultimately improve patient satisfaction.}, language = {en} } @article{DibakJdelRazoDeSanchoetal.2018, author = {Dibak, Manuel and J. del Razo, Mauricio and De Sancho, David and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations}, volume = {148}, journal = {Journal of Chemical Physics}, number = {214107}, doi = {10.1063/1.5020294}, year = {2018}, abstract = {Molecular dynamics (MD) simulations can model the interactions between macromolecules with high spatiotemporal resolution but at a high computational cost. By combining high-throughput MD with Markov state models (MSMs), it is now possible to obtain long-timescale behavior of small to intermediate biomolecules and complexes. To model the interactions of many molecules at large lengthscales, particle-based reaction-diffusion (RD) simulations are more suitable but lack molecular detail. Thus, coupling MSMs and RD simulations (MSM/RD) would be highly desirable, as they could efficiently produce simulations at large time- and lengthscales, while still conserving the characteristic features of the interactions observed at atomic detail. While such a coupling seems straightforward, fundamental questions are still open: Which definition of MSM states is suitable? Which protocol to merge and split RD particles in an association/dissociation reaction will conserve the correct bimolecular kinetics and thermodynamics? In this paper, we make the first step towards MSM/RD by laying out a general theory of coupling and proposing a first implementation for association/dissociation of a protein with a small ligand (A + B <--> C). Applications on a toy model and CO diffusion into the heme cavity of myoglobin are reported.}, language = {en} } @article{TackMukhopadhyayZachow2018, author = {Tack, Alexander and Mukhopadhyay, Anirban and Zachow, Stefan}, title = {Knee Menisci Segmentation using Convolutional Neural Networks: Data from the Osteoarthritis Initiative}, volume = {26}, journal = {Osteoarthritis and Cartilage}, number = {5}, doi = {10.1016/j.joca.2018.02.907}, pages = {680 -- 688}, year = {2018}, abstract = {Abstract: Objective: To present a novel method for automated segmentation of knee menisci from MRIs. To evaluate quantitative meniscal biomarkers for osteoarthritis (OA) estimated thereof. Method: A segmentation method employing convolutional neural networks in combination with statistical shape models was developed. Accuracy was evaluated on 88 manual segmentations. Meniscal volume, tibial coverage, and meniscal extrusion were computed and tested for differences between groups of OA, joint space narrowing (JSN), and WOMAC pain. Correlation between computed meniscal extrusion and MOAKS experts' readings was evaluated for 600 subjects. Suitability of biomarkers for predicting incident radiographic OA from baseline to 24 months was tested on a group of 552 patients (184 incident OA, 386 controls) by performing conditional logistic regression. Results: Segmentation accuracy measured as Dice Similarity Coefficient was 83.8\% for medial menisci (MM) and 88.9\% for lateral menisci (LM) at baseline, and 83.1\% and 88.3\% at 12-month follow-up. Medial tibial coverage was significantly lower for arthritic cases compared to non-arthritic ones. Medial meniscal extrusion was significantly higher for arthritic knees. A moderate correlation between automatically computed medial meniscal extrusion and experts' readings was found (ρ=0.44). Mean medial meniscal extrusion was significantly greater for incident OA cases compared to controls (1.16±0.93 mm vs. 0.83±0.92 mm; p<0.05). Conclusion: Especially for medial menisci an excellent segmentation accuracy was achieved. Our meniscal biomarkers were validated by comparison to experts' readings as well as analysis of differences w.r.t groups of OA, JSN, and WOMAC pain. It was confirmed that medial meniscal extrusion is a predictor for incident OA.}, language = {en} } @inproceedings{AmbellanTackEhlkeetal.2018, author = {Ambellan, Felix and Tack, Alexander and Ehlke, Moritz and Zachow, Stefan}, title = {Automated Segmentation of Knee Bone and Cartilage combining Statistical Shape Knowledge and Convolutional Neural Networks: Data from the Osteoarthritis Initiative}, booktitle = {Medical Imaging with Deep Learning}, year = {2018}, abstract = {We present a method for the automated segmentation of knee bones and cartilage from magnetic resonance imaging, that combines a priori knowledge of anatomical shape with Convolutional Neural Networks (CNNs). The proposed approach incorporates 3D Statistical Shape Models (SSMs) as well as 2D and 3D CNNs to achieve a robust and accurate segmentation of even highly pathological knee structures. The method is evaluated on data of the MICCAI grand challenge "Segmentation of Knee Images 2010". For the first time an accuracy equivalent to the inter-observer variability of human readers has been achieved in this challenge. Moreover, the quality of the proposed method is thoroughly assessed using various measures for 507 manual segmentations of bone and cartilage, and 88 additional manual segmentations of cartilage. Our method yields sub-voxel accuracy. In conclusion, combining of anatomical knowledge using SSMs with localized classification via CNNs results in a state-of-the-art segmentation method.}, language = {en} } @article{KlusNueskePeitzetal.2020, author = {Klus, Stefan and N{\"u}ske, Feliks and Peitz, Sebastian and Niemann, Jan-Hendrik and Clementi, Cecilia and Sch{\"u}tte, Christof}, title = {Data-driven approximation of the Koopman generator: Model reduction, system identification, and control}, volume = {406}, journal = {Physica D: Nonlinear Phenomena}, doi = {10.1016/j.physd.2020.132416}, year = {2020}, language = {en} } @article{ChaumelSchotteBizzarroetal.2020, author = {Chaumel, J{\´u}lia and Schotte, Merlind and Bizzarro, Joseph J. and Zaslansky, Paul and Fratzl, Peter and Baum, Daniel and Dean, Mason N.}, title = {Co-aligned chondrocytes: Zonal morphological variation and structured arrangement of cell lacunae in tessellated cartilage}, volume = {134}, journal = {Bone}, doi = {10.1016/j.bone.2020.115264}, pages = {115264}, year = {2020}, abstract = {In most vertebrates the embryonic cartilaginous skeleton is replaced by bone during development. During this process, cartilage cells (chondrocytes) mineralize the extracellular matrix and undergo apoptosis, giving way to bone cells (osteocytes). In contrast, sharks and rays (elasmobranchs) have cartilaginous skeletons throughout life, where only the surface mineralizes, forming a layer of tiles (tesserae). Elasmobranch chondrocytes, unlike those of other vertebrates, survive cartilage mineralization and are maintained alive in spaces (lacunae) within tesserae. However, the function(s) of the chondrocytes in the mineralized tissue remain unknown. Applying a custom analysis workflow to high-resolution synchrotron microCT scans of tesserae, we characterize the morphologies and arrangements of stingray chondrocyte lacunae, using lacunar morphology as a proxy for chondrocyte morphology. We show that the cell density is comparable in unmineralized and mineralized tissue from our study species and that cells maintain the similar volume even when they have been incorporated into tesserae. This discovery supports previous hypotheses that elasmobranch chondrocytes, unlike those of other taxa, do not proliferate, hypertrophy or undergo apoptosis during mineralization. Tessera lacunae show zonal variation in their shapes—being flatter further from and more spherical closer to the unmineralized cartilage matrix and larger in the center of tesserae— and show pronounced organization into parallel layers and strong orientation toward neighboring tesserae. Tesserae also exhibit local variation in lacunar density, with the density considerably higher near pores passing through the tesseral layer, suggesting pores and cells interact (e.g. that pores contain a nutrient source). We hypothesize that the different lacunar types reflect the stages of the tesserae formation process, while also representing local variation in tissue architecture and cell function. Lacunae are linked by small passages (canaliculi) in the matrix to form elongate series at the tesseral periphery and tight clusters in the center of tesserae, creating a rich connectivity among cells. The network arrangement and the shape variation of chondrocytes in tesserae indicate that cells may interact within and between tesserae and manage mineralization differently from chondrocytes in other vertebrates, perhaps performing analogous roles to osteocytes in bone.}, language = {en} } @article{SchotteChaumelDeanetal.2020, author = {Schotte, Merlind and Chaumel, J{\´u}lia and Dean, Mason N. and Baum, Daniel}, title = {Image analysis pipeline for segmentation of a biological porosity network, the lacuno-canalicular system in stingray tesserae}, volume = {7}, journal = {MethodsX}, doi = {10.1016/j.mex.2020.100905}, pages = {100905}, year = {2020}, abstract = {A prerequisite for many analysis tasks in modern comparative biology is the segmentation of 3-dimensional (3D) images of the specimens being investigated (e.g. from microCT data). Depending on the specific imaging technique that was used to acquire the images and on the image resolution, different segmentation tools will be required. While some standard tools exist that can often be applied for specific subtasks, building whole processing pipelines solely from standard tools is often difficult. Some tasks may even necessitate the implementation of manual interaction tools to achieve a quality that is sufficient for the subsequent analysis. In this work, we present a pipeline of segmentation tools that can be used for the semi-automatic segmentation and quantitative analysis of voids in tissue (i.e. internal structural porosity). We use this pipeline to analyze lacuno-canalicular networks in stingray tesserae from 3D images acquired with synchrotron microCT. * The first step of this processing pipeline, the segmentation of the tesserae, was performed using standard marker-based watershed segmentation. The efficient processing of the next two steps, that is, the segmentation of all lacunae spaces belonging to a specific tessera and the separation of these spaces into individual lacunae required modern, recently developed tools. * For proofreading, we developed a graph-based interactive method that allowed us to quickly split lacunae that were accidentally merged, and to merge lacunae that were wrongly split. * Finally, the tesserae and their corresponding lacunae were subdivided into anatomical regions of interest (structural wedges) using a semi- manual approach.}, language = {en} } @misc{MollenhauerSchusterKlusetal.2020, author = {Mollenhauer, Mattes and Schuster, Ingmar and Klus, Stefan and Sch{\"u}tte, Christof}, title = {Singular Value Decomposition of Operators on Reproducing Kernel Hilbert Spaces}, volume = {304}, journal = {Advances om Dynamics, Optimization and Computation. Series: Studies in Systems, Decision and Control. A volume dedicated to Michael Dellnitz on his 60th birthday}, editor = {Junge, Oliver and Sch{\"u}tze, O. and Froyland, Gary and Ober-Blobaum, S. and Padberg-Gehle, K.}, publisher = {Springer International}, isbn = {978-3-030-51264-4}, doi = {10.1007/978-3-030-51264-4_5}, pages = {109 -- 131}, year = {2020}, language = {en} } @article{AlchikhConradMaetal.2019, author = {Alchikh, Maren and Conrad, Tim and Ma, Xiaolin and Broberg, Eeva K. and Penttinen, P. and Reiche, J. and Biere, Barbara and Schweiger, Brunhilde and Rath, Barbara and Hoppe, Christian}, title = {Are we missing respiratory viral infections in infants and children? Comparison of a hospital-based quality management system with standard of care}, volume = {25}, journal = {Clinical Microbiology and Infection}, number = {3}, issn = {1469-0691}, doi = {10.1016/j.cmi.2018.05.023}, pages = {380.e9 -- 380.e16}, year = {2019}, language = {en} } @article{HaaseSunkaraKohletal.2019, author = {Haase, Tobias and Sunkara, Vikram and Kohl, Benjamin and Meier, Carola and Bußmann, Patricia and Becker, Jessica and Jagielski, Michal and von Kleist, Max and Ertel, Wolfgang}, title = {Discerning the spatio-temporal disease patterns of surgically induced OA mouse models}, volume = {14}, journal = {PLOS One}, number = {4}, publisher = {PLOS One}, doi = {10.1371/journal.pone.0213734}, year = {2019}, abstract = {Osteoarthritis (OA) is the most common cause of disability in ageing societies, with no effective therapies available to date. Two preclinical models are widely used to validate novel OA interventions (MCL-MM and DMM). Our aim is to discern disease dynamics in these models to provide a clear timeline in which various pathological changes occur. OA was surgically induced in mice by destabilisation of the medial meniscus. Analysis of OA progression revealed that the intensity and duration of chondrocyte loss and cartilage lesion formation were significantly different in MCL-MM vs DMM. Firstly, apoptosis was seen prior to week two and was narrowly restricted to the weight bearing area. Four weeks post injury the magnitude of apoptosis led to a 40-60\% reduction of chondrocytes in the non-calcified zone. Secondly, the progression of cell loss preceded the structural changes of the cartilage spatio-temporally. Lastly, while proteoglycan loss was similar in both models, collagen type II degradation only occurred more prominently in MCL-MM. Dynamics of chondrocyte loss and lesion formation in preclinical models has important implications for validating new therapeutic strategies. Our work could be helpful in assessing the feasibility and expected response of the DMM- and the MCL-MM models to chondrocyte mediated therapies.}, language = {en} } @article{RamsConrad2020, author = {Rams, Mona and Conrad, Tim}, title = {Dictionary Learning for transcriptomics data reveals type-specific gene modules in a multi-class setting}, volume = {62}, journal = {it - Information Technology}, number = {3-4}, publisher = {De Gruyter}, address = {Oldenbourg}, issn = {2196-7032}, doi = {https://doi.org/10.1515/itit-2019-0048}, year = {2020}, language = {en} } @article{RaySunkaraSchuetteetal.2020, author = {Ray, Sourav and Sunkara, Vikram and Sch{\"u}tte, Christof and Weber, Marcus}, title = {How to calculate pH-dependent binding rates for receptor-ligand systems based on thermodynamic simulations with different binding motifs}, volume = {46}, journal = {Molecular Simulation}, number = {18}, publisher = {Taylor and Francis}, doi = {10.1080/08927022.2020.1839660}, pages = {1443 -- 1452}, year = {2020}, abstract = {Molecular simulations of ligand-receptor interactions are a computational challenge, especially when their association- ('on'-rate) and dissociation- ('off'-rate) mechanisms are working on vastly differing timescales. One way of tackling this multiscale problem is to compute the free-energy landscapes, where molecular dynamics (MD) trajectories are used to only produce certain statistical ensembles. The approach allows for deriving the transition rates between energy states as a function of the height of the activation-energy barriers. In this article, we derive the association rates of the opioids fentanyl and N-(3-fluoro-1-phenethylpiperidin-4-yl)-N-phenyl propionamide (NFEPP) in a μ-opioid receptor by combining the free-energy landscape approach with the square-root-approximation method (SQRA), which is a particularly robust version of Markov modelling. The novelty of this work is that we derive the association rates as a function of the pH level using only an ensemble of MD simulations. We also verify our MD-derived insights by reproducing the in vitro study performed by the Stein Lab.}, language = {en} } @article{LindowBruenigDercksenetal.2020, author = {Lindow, Norbert and Br{\"u}nig, Florian and Dercksen, Vincent J. and Fabig, Gunar and Kiewisz, Robert and Redemann, Stefanie and M{\"u}ller-Reichert, Thomas and Prohaska, Steffen and Baum, Daniel}, title = {Semi-automatic stitching of filamentous structures in image stacks from serial-section electron tomography}, journal = {bioRxiv}, doi = {10.1101/2020.05.28.120899}, year = {2020}, abstract = {We present a software-assisted workflow for the alignment and matching of filamentous structures across a 3D stack of serial images. This is achieved by combining automatic methods, visual validation, and interactive correction. After an initial alignment, the user can continuously improve the result by interactively correcting landmarks or matches of filaments. Supported by a visual quality assessment of regions that have been already inspected, this allows a trade-off between quality and manual labor. The software tool was developed to investigate cell division by quantitative 3D analysis of microtubules (MTs) in both mitotic and meiotic spindles. For this, each spindle is cut into a series of semi-thick physical sections, of which electron tomograms are acquired. The serial tomograms are then stitched and non-rigidly aligned to allow tracing and connecting of MTs across tomogram boundaries. In practice, automatic stitching alone provides only an incomplete solution, because large physical distortions and a low signal-to-noise ratio often cause experimental difficulties. To derive 3D models of spindles despite the problems related to sample preparation and subsequent data collection, semi-automatic validation and correction is required to remove stitching mistakes. However, due to the large number of MTs in spindles (up to 30k) and their resulting dense spatial arrangement, a naive inspection of each MT is too time consuming. Furthermore, an interactive visualization of the full image stack is hampered by the size of the data (up to 100 GB). Here, we present a specialized, interactive, semi-automatic solution that considers all requirements for large-scale stitching of filamentous structures in serial-section image stacks. The key to our solution is a careful design of the visualization and interaction tools for each processing step to guarantee real-time response, and an optimized workflow that efficiently guides the user through datasets.}, language = {en} } @article{PimentelSzengelEhlkeetal.2020, author = {Pimentel, Pedro and Szengel, Angelika and Ehlke, Moritz and Lamecker, Hans and Zachow, Stefan and Estacio, Laura and Doenitz, Christian and Ramm, Heiko}, title = {Automated Virtual Reconstruction of Large Skull Defects using Statistical Shape Models and Generative Adversarial Networks}, volume = {12439}, journal = {Towards the Automatization of Cranial Implant Design in Cranioplasty}, editor = {Li, Jianning and Egger, Jan}, edition = {1}, publisher = {Springer International Publishing}, doi = {10.1007/978-3-030-64327-0_3}, pages = {16 -- 27}, year = {2020}, abstract = {We present an automated method for extrapolating missing regions in label data of the skull in an anatomically plausible manner. The ultimate goal is to design patient-speci� c cranial implants for correcting large, arbitrarily shaped defects of the skull that can, for example, result from trauma of the head. Our approach utilizes a 3D statistical shape model (SSM) of the skull and a 2D generative adversarial network (GAN) that is trained in an unsupervised fashion from samples of healthy patients alone. By � tting the SSM to given input labels containing the skull defect, a First approximation of the healthy state of the patient is obtained. The GAN is then applied to further correct and smooth the output of the SSM in an anatomically plausible manner. Finally, the defect region is extracted using morphological operations and subtraction between the extrapolated healthy state of the patient and the defective input labels. The method is trained and evaluated based on data from the MICCAI 2020 AutoImplant challenge. It produces state-of-the art results on regularly shaped cut-outs that were present in the training and testing data of the challenge. Furthermore, due to unsupervised nature of the approach, the method generalizes well to previously unseen defects of varying shapes that were only present in the hidden test dataset.}, language = {en} } @article{OeltzeJaffraMeuschkeNeugebaueretal.2019, author = {Oeltze-Jaffra, Steffen and Meuschke, Monique and Neugebauer, Mathias and Saalfeld, Sylvia and Lawonn, Kai and Janiga, Gabor and Hege, Hans-Christian and Zachow, Stefan and Preim, Bernhard}, title = {Generation and Visual Exploration of Medical Flow Data: Survey, Research Trends, and Future Challenges}, volume = {38}, journal = {Computer Graphics Forum}, number = {1}, publisher = {Wiley}, doi = {10.1111/cgf.13394}, pages = {87 -- 125}, year = {2019}, abstract = {Simulations and measurements of blood and air flow inside the human circulatory and respiratory system play an increasingly important role in personalized medicine for prevention, diagnosis, and treatment of diseases. This survey focuses on three main application areas. (1) Computational Fluid Dynamics (CFD) simulations of blood flow in cerebral aneurysms assist in predicting the outcome of this pathologic process and of therapeutic interventions. (2) CFD simulations of nasal airflow allow for investigating the effects of obstructions and deformities and provide therapy decision support. (3) 4D Phase-Contrast (4D PC) Magnetic Resonance Imaging (MRI) of aortic hemodynamics supports the diagnosis of various vascular and valve pathologies as well as their treatment. An investigation of the complex and often dynamic simulation and measurement data requires the coupling of sophisticated visualization, interaction, and data analysis techniques. In this paper, we survey the large body of work that has been conducted within this realm. We extend previous surveys by incorporating nasal airflow, addressing the joint investigation of blood flow and vessel wall properties, and providing a more fine-granular taxonomy of the existing techniques. From the survey, we extract major research trends and identify open problems and future challenges. The survey is intended for researchers interested in medical flow but also more general, in the combined visualization of physiology and anatomy, the extraction of features from flow field data and feature-based visualization, the visual comparison of different simulation results, and the interactive visual analysis of the flow field and derived characteristics.}, language = {en} } @article{RathConradMylesetal.2017, author = {Rath, Barbara and Conrad, Tim and Myles, Puja and Alchikh, Maren and Ma, Xiaolin and Hoppe, Christian and Tief, Franziska and Chen, Xi and Obermeier, Patrick and Kisler, Bron and Schweiger, Brunhilde}, title = {Influenza and other respiratory viruses: standardizing disease severity in surveillance and clinical trials}, volume = {15}, journal = {Expert Review of Anti-infective Therapy}, number = {6}, doi = {10.1080/14787210.2017.1295847}, pages = {545 -- 568}, year = {2017}, abstract = {Introduction: Influenza-Like Illness is a leading cause of hospitalization in children. Disease burden due to influenza and other respiratory viral infections is reported on a population level, but clinical scores measuring individual changes in disease severity are urgently needed. Areas covered: We present a composite clinical score allowing individual patient data analyses of disease severity based on systematic literature review and WHO-criteria for uncomplicated and complicated disease. The 22-item ViVI Disease Severity Score showed a normal distribution in a pediatric cohort of 6073 children aged 0-18 years (mean age 3.13; S.D. 3.89; range: 0 to 18.79). Expert commentary: The ViVI Score was correlated with risk of antibiotic use as well as need for hospitalization and intensive care. The ViVI Score was used to track children with influenza, respiratory syncytial virus, human metapneumovirus, human rhinovirus, and adenovirus infections and is fully compliant with regulatory data standards. The ViVI Disease Severity Score mobile application allows physicians to measure disease severity at the point-of care thereby taking clinical trials to the next level.}, language = {en} } @article{MarzbanConradMarbanetal.2018, author = {Marzban, Forough and Conrad, Tim and Marban, Pouria and Sodoudi, Sahar}, title = {Estimation of the Near-Surface Air Temperature during the Day and Nighttime from MODIS in Berlin, Germany}, volume = {7}, journal = {International Journal of Advanced Remote Sensing and GIS}, number = {1}, doi = {10.23953/cloud.ijarsg.337}, pages = {2478 -- 2517}, year = {2018}, abstract = {Air temperature (Tair or T2m) is an important climatological variable for forest biosphere processes and climate change research. Due to the low density and the uneven distribution of weather stations, traditional ground-based observations cannot accurately capture the spatial distribution of Tair. In this study, Tair in Berlin is estimated during the day and night time over six land cover/land use (LC/LU) types by satellite remote sensing data over a large domain and a relatively long period (7 years). Aqua and Terra MODIS (Moderate Resolution Imaging Spectroradiometer) data and meteorological data for the period from 2007 to 2013 were collected to estimate Tair. Twelve environmental variables (land surface temperature (LST), normalized difference vegetation index (NDVI), Julian day, latitude, longitude, Emissivity31, Emissivity32, altitude, albedo, wind speed, wind direction and air pressure) were selected as predictors. Moreover, a comparison between LST from MODIS Terra and Aqua with daytime and night time air temperatures (Tday, Tnight) was done respectively and in addition, the spatial variability of LST and Tair relationship by applying a varying window size on the MODIS LST grid was examined. An analysis of the relationship between the observed Tair and the spatially averaged remotely sensed LST, indicated that 3 × 3 and 1 × 1 pixel size was the optimal window size for the statistical model estimating Tair from MODIS data during the day and night time, respectively. Three supervised learning methods (Adaptive Neuro Fuzzy Inference system (ANFIS), Artificial Neural Network (ANN) and Support vector machine (SVR)) were used to estimate Tair during the day and night time, and their performances were validated by cross-validation for each LC/LU. Moreover, tuning the hyper parameters of some models like SVR and ANN were investigated. For tuning the hyper parameters of SVR, Simulated Annealing (SA) was applied (SA-SVR model) and a multiple-layer feed-forward (MLF) neural networks with three layers and different nodes in hidden layers are used with Levenber-Marquardt back-propagation (LM-BP), in order to achieve higher accuracy in the estimation of Tair. Results indicated that the ANN model achieved better accuracy (RMSE= 2.16°C, MAE = 1.69°C, R2 = 0.95) than SA_SVR model (RMSE= 2.50°C, MAE = 1.92°C, R2 = 0.91) and ANFIS model (RMSE= 2.88°C, MAE= 2.2°C, R2 = 0.89) over six LC/LU during the day and night time. The Q-Q diagram of SA-SVR, ANFIS and NN show that all three models slightly tended to underestimate and overestimate the extreme and low temperatures for all LC/LU classes during the day and night time. The weak performance in the extreme and low temperatures are a consequence of the small number of data in these temperatures. These satisfactory results indicate that this approach is proper for estimating air temperature and spatial window size is an important factor that should be considered in the estimation of air temperature.}, language = {en} } @article{KoltaiWuNoeetal.2018, author = {Koltai, Peter and Wu, Hao and No{\´e}, Frank and Sch{\"u}tte, Christof}, title = {Optimal data-driven estimation of generalized Markov state models for non-equilibrium dynamics}, volume = {6}, journal = {Computation}, number = {1}, publisher = {MDPI}, address = {Basel, Switzerland}, doi = {10.3390/computation6010022}, pages = {22}, year = {2018}, language = {en} } @article{DibakdelRazodeSanchoetal.2018, author = {Dibak, Manuel and del Razo, Mauricio J. and de Sancho, David and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations}, volume = {148}, journal = {Journal of Chemical Physics}, number = {21}, doi = {10.1063/1.5020294}, year = {2018}, abstract = {Molecular dynamics (MD) simulations can model the interactions between macromolecules with high spatiotemporal resolution but at a high computational cost. By combining high-throughput MD with Markov state models (MSMs), it is now possible to obtain long time-scale behavior of small to intermediate biomolecules and complexes. To model the interactions of many molecules at large length scales, particle-based reaction-diffusion (RD) simulations are more suitable but lack molecular detail. Thus, coupling MSMs and RD simulations (MSM/RD) would be highly desirable, as they could efficiently produce simulations at large time and length scales, while still conserving the characteristic features of the interactions observed at atomic detail. While such a coupling seems straightforward, fundamental questions are still open: Which definition of MSM states is suitable? Which protocol to merge and split RD particles in an association/dissociation reaction will conserve the correct bimolecular kinetics and thermodynamics? In this paper, we make the first step toward MSM/RD by laying out a general theory of coupling and proposing a first implementation for association/dissociation of a protein with a small ligand (A + B ⇌ C). Applications on a toy model and CO diffusion into the heme cavity of myoglobin are reported.}, language = {en} } @inproceedings{EstacioEhlkeTacketal.2021, author = {Estacio, Laura and Ehlke, Moritz and Tack, Alexander and Castro-Gutierrez, Eveling and Lamecker, Hans and Mora, Rensso and Zachow, Stefan}, title = {Unsupervised Detection of Disturbances in 2D Radiographs}, booktitle = {2021 IEEE 18th International Symposium on Biomedical Imaging (ISBI)}, doi = {10.1109/ISBI48211.2021.9434091}, pages = {367 -- 370}, year = {2021}, abstract = {We present a method based on a generative model for detection of disturbances such as prosthesis, screws, zippers, and metals in 2D radiographs. The generative model is trained in an unsupervised fashion using clinical radiographs as well as simulated data, none of which contain disturbances. Our approach employs a latent space consistency loss which has the benefit of identifying similarities, and is enforced to reconstruct X-rays without disturbances. In order to detect images with disturbances, an anomaly score is computed also employing the Frechet distance between the input X-ray and the reconstructed one using our generative model. Validation was performed using clinical pelvis radiographs. We achieved an AUC of 0.77 and 0.83 with clinical and synthetic data, respectively. The results demonstrated a good accuracy of our method for detecting outliers as well as the advantage of utilizing synthetic data.}, language = {en} } @article{SekuboyinaBayatHusseinietal.2020, author = {Sekuboyina, Anjany and Bayat, Amirhossein and Husseini, Malek E. and L{\"o}ffler, Maximilian and Li, Hongwei and Tetteh, Giles and Kukačka, Jan and Payer, Christian and Štern, Darko and Urschler, Martin and Chen, Maodong and Cheng, Dalong and Lessmann, Nikolas and Hu, Yujin and Wang, Tianfu and Yang, Dong and Xu, Daguang and Ambellan, Felix and Amiranashvili, Tamaz and Ehlke, Moritz and Lamecker, Hans and Lehnert, Sebastian and Lirio, Marilia and de Olaguer, Nicol{\´a}s P{\´e}rez and Ramm, Heiko and Sahu, Manish and Tack, Alexander and Zachow, Stefan and Jiang, Tao and Ma, Xinjun and Angerman, Christoph and Wang, Xin and Wei, Qingyue and Brown, Kevin and Wolf, Matthias and Kirszenberg, Alexandre and Puybareau, {\´E}lodie and Valentinitsch, Alexander and Rempfler, Markus and Menze, Bj{\"o}rn H. and Kirschke, Jan S.}, title = {VerSe: A Vertebrae Labelling and Segmentation Benchmark for Multi-detector CT Images}, journal = {arXiv}, arxiv = {http://arxiv.org/abs/2001.09193}, year = {2020}, language = {en} } @article{WulkowConradDjurdjevacConradetal.2021, author = {Wulkow, Hanna and Conrad, Tim and Djurdjevac Conrad, Natasa and M{\"u}ller, Sebastian A. and Nagel, Kai and Sch{\"u}tte, Christof}, title = {Prediction of Covid-19 spreading and optimal coordination of counter-measures: From microscopic to macroscopic models to Pareto fronts}, volume = {16}, journal = {PLOS One}, number = {4}, publisher = {Public Library of Science}, doi = {10.1371/journal.pone.0249676}, year = {2021}, language = {en} } @article{BanischDjurdjevacConradSchuette2015, author = {Banisch, Ralf and Djurdjevac Conrad, Natasa and Sch{\"u}tte, Christof}, title = {Reactive flows and unproductive cycles for random walks on complex networks}, journal = {The European Physical Journal Special Topics, vol. 224, iss. 12 (2015) pp. 2369-2387}, doi = {10.1140/epjst/e2015-02417-8}, year = {2015}, language = {en} } @misc{BanischDjurdjevacConradSchuette2015, author = {Banisch, Ralf and Djurdjevac Conrad, Natasa and Sch{\"u}tte, Christof}, title = {Reactive flows and unproductive cycles for random walks on complex networks}, issn = {1438-0064}, doi = {10.1140/epjst/e2015-02417-8}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-54239}, year = {2015}, abstract = {We present a comprehensive theory for analysis and understanding of transition events between an initial set A and a target set B for general ergodic finite-state space Markov chains or jump processes, including random walks on networks as they occur, e.g., in Markov State Modelling in molecular dynamics. The theory allows us to decompose the probability flow generated by transition events between the sets A and B into the productive part that directly flows from A to B through reaction pathways and the unproductive part that runs in loops and is supported on cycles of the underlying network. It applies to random walks on directed networks and nonreversible Markov processes and can be seen as an extension of Transition Path Theory. Information on reaction pathways and unproductive cycles results from the stochastic cycle decomposition of the underlying network which also allows to compute their corresponding weight, thus characterizing completely which structure is used how often in transition events. The new theory is illustrated by an application to a Markov State Model resulting from weakly damped Langevin dynamics where the unproductive cycles are associated with periodic orbits of the underlying Hamiltonian dynamics.}, language = {en} } @article{WinkelmannSchuettevonKleist2014, author = {Winkelmann, Stefanie and Sch{\"u}tte, Christof and von Kleist, Max}, title = {Markov Control Processes with Rare State Observation: Theory and Application to Treatment Scheduling in HIV-1}, volume = {12}, journal = {Communications in Mathematical Sciences}, number = {5}, doi = {10.4310/CMS.2014.v12.n5.a4}, pages = {859 -- 877}, year = {2014}, abstract = {Markov Decision Processes (MDP) or Partially Observable MDPs (POMDP) are used for modelling situations in which the evolution of a process is partly random and partly controllable. These MDP theories allow for computing the optimal control policy for processes that can continuously or frequently be observed, even if only partially. However, they cannot be applied if state observation is very costly and therefore rare (in time). We present a novel MDP theory for rare, costly observations and derive the corresponding Bellman equation. In the new theory, state information can be derived for a particular cost after certain, rather long time intervals. The resulting information costs enter into the total cost and thus into the optimization criterion. This approach applies to many real world problems, particularly in the medical context, where the medical condition is examined rather rarely because examination costs are high. At the same time, the approach allows for efficient numerical realization. We demonstrate the usefulness of the novel theory by determining, from the national economic perspective, optimal therapeutic policies for the treatment of the human immunodeficiency virus (HIV) in resource-rich and resource-poor settings. Based on the developed theory and models, we discover that available drugs may not be utilized efficiently in resource-poor settings due to exorbitant diagnostic costs.}, language = {en} } @article{DjurdjevacConradBanischSchuette2015, author = {Djurdjevac Conrad, Natasa and Banisch, Ralf and Sch{\"u}tte, Christof}, title = {Modularity of Directed Networks: Cycle Decomposition Approach}, journal = {Journal of Computational Dynamics 2 (2015) pp. 1-24}, doi = {10.3934/jcd.2015.2.1}, year = {2015}, abstract = {The problem of decomposing networks into modules (or clusters) has gained much attention in recent years, as it can account for a coarsegrained description of complex systems, often revealing functional subunits of these systems. A variety of module detection algorithms have been proposed, mostly oriented towards finding hard partitionings of undirected networks. Despite the increasing number of fuzzy clustering methods for directed networks, many of these approaches tend to neglect important directional information. In this paper, we present a novel random walk based approach for finding fuzzy partitions of directed, weighted networks, where edge directions play a crucial role in defining how well nodes in a module are interconnected. We will show that cycle decomposition of a random walk process connects the notion of network modules and information transport in a network, leading to a new, symmetric measure of node communication. Finally, we will use this measure to introduce a communication graph, for which we will show that although being undirected it inherits all necessary information about modular structures from the original network.}, language = {en} } @misc{WangSchuette2014, author = {Wang, Han and Sch{\"u}tte, Christof}, title = {Building Markov State Models for Periodically Driven Non-Equilibrium Systems}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-53167}, year = {2014}, abstract = {Recent years have seen an increased interest in non-equilibrium molecular dynamics (NEMD) simulations, especially for molecular systems with periodic forcing by external fields, e.g., in the context of studying effects of electromagnetic radiation on the human body tissue. Lately, an NEMD methods with local thermostating has been proposed that allows for studying non-equilibrium processes in a statistically reliable and thermodynamically consistent way. In this article, we demonstrate how to construct Markov State Models (MSMs) for such NEMD simulations. MSM building has been well-established for systems in equilibrium where MSMs with just a few (macro-)states allow for accurate reproduction of the essential kinetics of the molecular system under consideration. Non-equilibrium MSMs have been lacking so far. The article presents how to construct such MSMs and illustrates their validity and usefulness for the case of conformation dynamics of alanine dipeptide in an external electric field.}, language = {en} } @misc{GulSchuetteBernhard2015, author = {Gul, Raheem and Sch{\"u}tte, Christof and Bernhard, Stefan}, title = {Mathematical modeling and sensitivity analysis of arterial anastomosis in arm arteries}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-54339}, year = {2015}, abstract = {Cardiovascular diseases are one of the major problems in medicine today and the number of patients increases worldwide. To find the most efficient treatment, prior knowledge about function and dysfunction of the cardiovas- cular system is required and methods need to be developed that identify the disease in an early stage. Mathematical modeling is a powerful tool for prediction and investigation of cardiovascular diseases. It has been shown that the Windkessel model, being based on an analogy between electrical circuits and fluid flow, is a simple but effective method to model the human cardiovascular system. In this paper, we have applied parametric local sensitivity analysis (LSA) to a linear elastic model of the arm arteries, to find and rank sensitive param- eters that may be helpful in clinical diagnosis. A computational model for end-to-side anastomosis (superior ulnar collateral anastomosis with posterior ulnar recurrent, SUC-PUR) is carried out to study the effects of some clinically relevant haemodynamic parameters like blood flow resistance and terminal re- sistance on pressure and flow at different locations of the arm artery. In this context, we also discuss the spatio-temporal dependency of local sensitivities. The sensitivities with respect to cardiovascular parameters reveal the flow resistance and diameter of the vessels as most sensitive parameters. These parameters play a key role in diagnosis of severe stenosis and aneurysms. In contrast, wall thickness and elastic modulus are found to be less sensitive.}, language = {en} } @inproceedings{Mukhopadhyay2016, author = {Mukhopadhyay, Anirban}, title = {Total Variation Random Forest: Fully automatic MRI segmentation in congenital heart disease}, volume = {LNCS 10129}, booktitle = {RAMBO 2016, HVSMR 2016: Reconstruction, Segmentation, and Analysis of Medical Images}, doi = {10.1007/978-3-319-52280-7_17}, pages = {165 -- 171}, year = {2016}, language = {en} } @inproceedings{MukhopadhyayMorilloZachowetal.2016, author = {Mukhopadhyay, Anirban and Morillo, Oscar and Zachow, Stefan and Lamecker, Hans}, title = {Robust and Accurate Appearance Models Based on Joint Dictionary Learning Data from the Osteoarthritis Initiative}, volume = {9993}, booktitle = {Lecture Notes in Computer Science, Patch-Based Techniques in Medical Imaging. Patch-MI 2016}, doi = {10.1007/978-3-319-47118-1_4}, pages = {25 -- 33}, year = {2016}, abstract = {Deformable model-based approaches to 3D image segmentation have been shown to be highly successful. Such methodology requires an appearance model that drives the deformation of a geometric model to the image data. Appearance models are usually either created heuristically or through supervised learning. Heuristic methods have been shown to work effectively in many applications but are hard to transfer from one application (imaging modality/anatomical structure) to another. On the contrary, supervised learning approaches can learn patterns from a collection of annotated training data. In this work, we show that the supervised joint dictionary learning technique is capable of overcoming the traditional drawbacks of the heuristic approaches. Our evaluation based on two different applications (liver/CT and knee/MR) reveals that our approach generates appearance models, which can be used effectively and efficiently in a deformable model-based segmentation framework.}, language = {en} } @article{SahuMukhopadhyaySzengeletal.2017, author = {Sahu, Manish and Mukhopadhyay, Anirban and Szengel, Angelika and Zachow, Stefan}, title = {Addressing multi-label imbalance problem of Surgical Tool Detection using CNN}, volume = {12}, journal = {International Journal of Computer Assisted Radiology and Surgery}, number = {6}, publisher = {Springer}, doi = {10.1007/s11548-017-1565-x}, pages = {1013 -- 1020}, year = {2017}, abstract = {Purpose: A fully automated surgical tool detection framework is proposed for endoscopic video streams. State-of-the-art surgical tool detection methods rely on supervised one-vs-all or multi-class classification techniques, completely ignoring the co-occurrence relationship of the tools and the associated class imbalance. Methods: In this paper, we formulate tool detection as a multi-label classification task where tool co-occurrences are treated as separate classes. In addition, imbalance on tool co-occurrences is analyzed and stratification techniques are employed to address the imbalance during Convolutional Neural Network (CNN) training. Moreover, temporal smoothing is introduced as an online post-processing step to enhance run time prediction. Results: Quantitative analysis is performed on the M2CAI16 tool detection dataset to highlight the importance of stratification, temporal smoothing and the overall framework for tool detection. Conclusion: The analysis on tool imbalance, backed by the empirical results indicates the need and superiority of the proposed framework over state-of-the-art techniques.}, language = {en} } @article{KlusSchuette2016, author = {Klus, Stefan and Sch{\"u}tte, Christof}, title = {Towards tensor-based methods for the numerical approximation of the Perron-Frobenius and Koopman operator}, volume = {3}, journal = {Journal of Computational Dynamics}, number = {2}, doi = {10.3934/jcd.2016007}, pages = {139 -- 161}, year = {2016}, abstract = {The global behavior of dynamical systems can be studied by analyzing the eigenvalues and corresponding eigenfunctions of linear operators associated with the system. Two important operators which are frequently used to gain insight into the system's behavior are the Perron-Frobenius operator and the Koopman operator. Due to the curse of dimensionality, computing the eigenfunctions of high-dimensional systems is in general infeasible. We will propose a tensor-based reformulation of two numerical methods for computing finite-dimensional approximations of the aforementioned infinite-dimensional operators, namely Ulam's method and Extended Dynamic Mode Decomposition (EDMD). The aim of the tensor formulation is to approximate the eigenfunctions by low-rank tensors, potentially resulting in a significant reduction of the time and memory required to solve the resulting eigenvalue problems, provided that such a low-rank tensor decomposition exists. Typically, not all variables of a high-dimensional dynamical system contribute equally to the system's behavior, often the dynamics can be decomposed into slow and fast processes, which is also reflected in the eigenfunctions. Thus, the weak coupling between different variables might be approximated by low-rank tensor cores. We will illustrate the efficiency of the tensor-based formulation of Ulam's method and EDMD using simple stochastic differential equations.}, language = {en} } @article{KlusKoltaiSchuette2016, author = {Klus, Stefan and Koltai, Peter and Sch{\"u}tte, Christof}, title = {On the numerical approximation of the Perron-Frobenius and Koopman operator}, volume = {3}, journal = {Journal of Computational Dynamics}, number = {1}, doi = {10.3934/jcd.2016003}, pages = {51 -- 77}, year = {2016}, abstract = {Information about the behavior of dynamical systems can often be obtained by analyzing the eigenvalues and corresponding eigenfunctions of linear operators associated with a dynamical system. Examples of such operators are the Perron-Frobenius and the Koopman operator. In this paper, we will review di� fferent methods that have been developed over the last decades to compute � infinite-dimensional approximations of these in� finite-dimensional operators - in particular Ulam's method and Extended Dynamic Mode Decomposition (EDMD) - and highlight the similarities and di� fferences between these approaches. The results will be illustrated using simple stochastic di� fferential equations and molecular dynamics examples.}, language = {en} } @article{HuttaryGoubergritsSchuetteetal.2017, author = {Huttary, Rudolf and Goubergrits, Leonid and Sch{\"u}tte, Christof and Bernhard, Stefan}, title = {Simulation, Identification and Statistical Variation in Cardiovascular Analysis (SISCA) - a Software Framework for Multi-compartment Lumped Modeling}, volume = {87}, journal = {Computers in Biology and Medicine}, doi = {10.1016/j.compbiomed.2017.05.021}, pages = {104 -- 123}, year = {2017}, language = {en} } @misc{GuptaGramatkeEinspanieretal.2017, author = {Gupta, Pooja and Gramatke, Annika and Einspanier, Ralf and Sch{\"u}tte, Christof and von Kleist, Max and Sharbati, Jutta}, title = {In silicio cytotoxicity assessment on cultured rat intestinal cells deduced from cellular impedance measurements}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-62666}, year = {2017}, abstract = {Early and reliable identification of chemical toxicity is of utmost importance. At the same time, reduction of animal testing is paramount. Therefore, methods that improve the interpretability and usability of in vitro assays are essential. xCELLigence's real-time cell analyzer (RTCA) provides a novel, fast and cost effective in vitro method to probe compound toxicity. We developed a simple mathematical framework for the qualitative and quantitative assessment of toxicity for RTCA measurements. Compound toxicity, in terms of its 50\% inhibitory concentration IC_{50} on cell growth, and parameters related to cell turnover were estimated on cultured IEC-6 cells exposed to 10 chemicals at varying concentrations. Our method estimated IC50 values of 113.05, 7.16, 28.69 and 725.15 μM for the apparently toxic compounds 2-acetylamino-fluorene, aflatoxin B1, benzo-[a]-pyrene and chloramphenicol in the tested cell line, in agreement with literature knowledge. IC_{50} values of all apparent in vivo non-toxic compounds were estimated to be non-toxic by our method. Corresponding estimates from RTCA's in-built model gave false positive (toxicity) predictions in 5/10 cases. Taken together, our proposed method reduces false positive predictions and reliably identifies chemical toxicity based on impedance measurements. The source code for the developed method including instructions is available at https://git.zib.de/bzfgupta/toxfit/tree/master.}, language = {en} } @article{KlusNueskeKoltaietal.2018, author = {Klus, Stefan and N{\"u}ske, Feliks and Koltai, Peter and Wu, Hao and Kevrekidis, Ioannis and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {Data-driven model reduction and transfer operator approximation}, volume = {28}, journal = {Journal of Nonlinear Science}, number = {3}, doi = {10.1007/s00332-017-9437-7}, pages = {985 -- 1010}, year = {2018}, language = {en} } @article{KoltaiCiccottiSchuette2016, author = {Koltai, Peter and Ciccotti, Giovanni and Sch{\"u}tte, Christof}, title = {On Markov state models for non-equilibrium molecular dynamics}, volume = {145}, journal = {The Journal of Chemical Physics}, number = {174103}, doi = {10.1063/1.4966157}, year = {2016}, language = {en} } @article{ConradKarschObermeieretal.2015, author = {Conrad, Tim and Karsch, K. and Obermeier, Patrick and Seeber, L. and Chen, X. and Tief, Franziska and Muehlhans, S. and Hoppe, Christian and Boettcher, Sindy and Diedrich, S. and Rath, Barbara}, title = {Human Parechovirus Infections Associated with Seizures and Rash: A Syndromic Surveillance Study in Children}, volume = {34}, journal = {The Pediatric Infectious Disease Journal}, number = {10}, year = {2015}, language = {en} } @article{ConradMireles2015, author = {Conrad, Tim and Mireles, Victor}, title = {Minimum-overlap clusterings and the sparsity of overcomplete decompositions of binary matrics}, volume = {51}, journal = {Procedia Computer Science}, doi = {10.1016/j.procs.2015.05.500}, pages = {2967 -- 2971}, year = {2015}, abstract = {Given a set of n binary data points, a widely used technique is to group its features into k clusters. In the case where n {\ensuremath{<}} k, the question of how overlapping are the clusters becomes of interest. In this paper we approach the question through matrix decomposition, and relate the degree of overlap with the sparsity of one of the resulting matrices. We present analytical results regarding bounds on this sparsity, and a heuristic to estimate the minimum amount of overlap that an exact grouping of features into k clusters must have. As shown below, adding new data will not alter this minimum amount of overlap.}, language = {en} } @article{ConradShao2015, author = {Conrad, Tim and Shao, Borong}, title = {Are NoSQL data stores useful for bioinformatics researchers?}, volume = {3}, journal = {International Journal on Recent and Innovation Trends in Computing and Communication (IJRITCC)}, number = {3}, doi = {10.17762/ijritcc2321-8169.1503176}, pages = {1704 -- 1708}, year = {2015}, language = {en} } @article{ConradLeichtleCeglareketal.2013, author = {Conrad, Tim and Leichtle, Alexander Benedikt and Ceglarek, Uta and Weinert, P. and Nakas, C.T. and Nuoffer, Jean-Marc and Kase, Julia and Witzigmann, Helmut and Thiery, Joachim and Fiedler, Georg Martin}, title = {Pancreatic carcinoma, pancreatitis, and healthy controls - metabolite models in a three-class diagnostic dilemma}, volume = {9}, journal = {Metabolomics}, number = {3}, doi = {10.1007/s11306-012-0476-7}, pages = {677 -- 687}, year = {2013}, abstract = {Background: Metabolomics as one of the most rapidly growing technologies in the ?-omics?field denotes the comprehensive analysis of low molecular-weight compounds and their pathways. Cancer-specific alterations of the metabolome can be detected by high-throughput massspectrometric metabolite profiling and serve as a considerable source of new markers for the early differentiation of malignant diseases as well as their distinction from benign states. However, a comprehensive framework for the statistical evaluation of marker panels in a multi-class setting has not yet been established. Methods: We collected serum samples of 40 pancreatic carcinoma patients, 40 controls, and 23 pancreatitis patients according to standard protocols and generated amino acid profiles by routine mass-spectrometry. In an intrinsic three-class bioinformatic approach we compared these profiles, evaluated their selectivity and computed multi-marker panels combined with the conventional tumor marker CA 19-9. Additionally, we tested for non-inferiority and superiority to determine the diagnostic surplus value of our multi-metabolite marker panels.  Results: Compared to CA 19-9 alone, the combined amino acid-based metabolite panel had a superior selectivity for the discrimination of healthy controls, pancreatitis, and pancreatic carcinoma patients [Volume under ROC surface (VUS) = 0.891 (95\\% CI 0.794 - 0.968)]. Conclusions: We combined highly standardized samples, a three-class study design, a highthroughput mass-spectrometric technique, and a comprehensive bioinformatic framework to identify metabolite panels selective for all three groups in a single approach. Our results suggest that metabolomic profiling necessitates appropriate evaluation strategies and ?despite all its current limitations? can deliver marker panels with high selectivity even in multi-class settings.}, language = {en} } @article{GuptaReinschSpoetteretal.2013, author = {Gupta, Pooja and Reinsch, Norbert and Sp{\"o}tter, Andreas and Conrad, Tim and Bienefeld, Kaspar}, title = {Accuracy of the unified approach in maternally influenced traits - illustrated by a simulation study in the honey bee (Apis mellifera)}, volume = {14}, journal = {BMC Genetics}, number = {36}, doi = {10.1186/1471-2156-14-36}, year = {2013}, language = {en} } @article{ConradYou2016, author = {Conrad, Tim and You, Xintian}, title = {Acfs: accurate circRNA identification and quantification from NGS data}, volume = {6}, journal = {Nature Scientific Reports}, doi = {10.1038/srep38820}, year = {2016}, abstract = {Circular RNAs (circRNAs) are a group of single-stranded RNAs in closed circular form. They are splicing-generated, widely expressed in various tissues and have functional implications in development and diseases. To facilitate genome-wide characterization of circRNAs using RNA-Seq data, we present a freely available software package named acfs. Acfs allows de novo, accurate and fast identification and abundance quantification of circRNAs from single- and paired-ended RNA-Seq data. On simulated datasets, acfs achieved the highest F1 accuracy and lowest false discovery rate among current state-of-the-art tools. On real-world datasets, acfs efficiently identified more bona fide circRNAs. Furthermore, we demonstrated the power of circRNA analysis on two leukemia datasets. We identified a set of circRNAs that are differentially expressed between AML and APL samples, which might shed light on the potential molecular classification of complex diseases using circRNA profiles. Moreover, chromosomal translocation, as manifested in numerous diseases, could produce not only fusion transcripts but also fusion circRNAs of clinical relevance. Featured with high accuracy, low FDR and the ability to identify fusion circRNAs, we believe that acfs is well suited for a wide spectrum of applications in characterizing the landscape of circRNAs from non-model organisms to cancer biology.}, language = {en} } @article{HoppeObermeierMehlhansetal.2016, author = {Hoppe, Christian and Obermeier, Patrick and Mehlhans, S. and Alchikh, Maren and Seeber, L. and Tief, Franziska and Karsch, K. and Chen, X. and Boettcher, Sindy and Diedrich, S. and Conrad, Tim}, title = {Innovative Digital Tools and Surveillance Systems for the Timely Detection of Adverse Events at the Point of Care: A Proof-of-Concept Study}, volume = {39}, journal = {Drug Safety}, number = {10}, doi = {10.1007/s40264-016-0437-6}, pages = {977 -- 988}, year = {2016}, abstract = {Regulatory authorities often receive poorly structured safety reports requiring considerable effort to investigate potential adverse events post hoc. Automated question-and-answer systems may help to improve the overall quality of safety information transmitted to pharmacovigilance agencies. This paper explores the use of the VACC-Tool (ViVI Automated Case Classification Tool) 2.0, a mobile application enabling physicians to classify clinical cases according to 14 pre-defined case definitions for neuroinflammatory adverse events (NIAE) and in full compliance with data standards issued by the Clinical Data Interchange Standards Consortium. METHODS: The validation of the VACC-Tool 2.0 (beta-version) was conducted in the context of a unique quality management program for children with suspected NIAE in collaboration with the Robert Koch Institute in Berlin, Germany. The VACC-Tool was used for instant case classification and for longitudinal follow-up throughout the course of hospitalization. Results were compared to International Classification of Diseases , Tenth Revision (ICD-10) codes assigned in the emergency department (ED). RESULTS: From 07/2013 to 10/2014, a total of 34,368 patients were seen in the ED, and 5243 patients were hospitalized; 243 of these were admitted for suspected NIAE (mean age: 8.5 years), thus participating in the quality management program. Using the VACC-Tool in the ED, 209 cases were classified successfully, 69 \\% of which had been missed or miscoded in the ED reports. Longitudinal follow-up with the VACC-Tool identified additional NIAE. CONCLUSION: Mobile applications are taking data standards to the point of care, enabling clinicians to ascertain potential adverse events in the ED setting and during inpatient follow-up. Compliance with Clinical Data Interchange Standards Consortium (CDISC) data standards facilitates data interoperability according to regulatory requirements.}, language = {en} } @inproceedings{ShaoConrad2016, author = {Shao, Borong and Conrad, Tim}, title = {Epithelial Mesenchymal Transition Regulatory Network-based Feature Selection in Lung Cancer Prognosis Prediction}, volume = {9656}, booktitle = {Lecture Notes in Computer Science (LNCS)}, doi = {10.1007/978-3-319-31744-1_13}, pages = {1235 -- 146}, year = {2016}, abstract = {Feature selection technique is often applied in identifying cancer prognosis biomarkers. However, many feature selection methods are prone to over-fitting or poor biological interpretation when applied on biological high-dimensional data. Network-based feature selection and data integration approaches are proposed to identify more robust biomarkers. We conducted experiments to investigate the advantages of the two approaches using epithelial mesenchymal transition regulatory network, which is demonstrated as highly relevant to cancer prognosis. We obtained data from The Cancer Genome Atlas. Prognosis prediction was made using Support Vector Machine. Under our experimental settings, the results showed that network-based features gave significantly more accurate predictions than individual molecular features, and features selected from integrated data (RNA-Seq and micro-RNA data) gave significantly more accurate predictions than features selected from single source data (RNA-Seq data). Our study indicated that biological network-based feature transformation and data integration are two useful approaches to identify robust cancer biomarkers.}, language = {en} } @article{TiefHoppeSeeberetal.2016, author = {Tief, Franziska and Hoppe, Christian and Seeber, L. and Obermeier, Patrick and Chen, X. and Karsch, K. and Muehlhans, S. and Adamou, E. and Conrad, Tim and Schweiger, Brunhilde and Adam, T. and Rath, Barbara}, title = {An inception cohort study assessing the role of bacterial co-infections in children with influenza and ILI and a clinical decision model for stringent antibiotic use}, volume = {21}, journal = {Antiviral Therapy}, doi = {10.3851/IMP3034}, pages = {413 -- 424}, year = {2016}, abstract = {BACKGROUND: Influenza-like illness (ILI) is a common reason for paediatric consultations. Viral causes predominate, but antibiotics are used frequently. With regard to influenza, pneumococcal coinfections are considered major contributors to morbidity/mortality. METHODS: In the context of a perennial quality management (QM) programme at the Charit{\'e} Departments of Paediatrics and Microbiology in collaboration with the Robert Koch Institute, children aged 0-18 years presenting with signs and symptoms of ILI were followed from the time of initial presentation until hospital discharge (Charit{\'e} Influenza-Like Disease = ChILD Cohort). An independent QM team performed highly standardized clinical assessments using a disease severity score based on World Health Organization criteria for uncomplicated and complicated/progressive disease. Nasopharyngeal and pharyngeal samples were collected for viral reverse transcription polymerase chain reaction and bacterial culture/sensitivity and MaldiTOF analyses. The term 'detection' was used to denote any evidence of viral or bacterial pathogens in the (naso)pharyngeal cavity. With the ChILD Cohort data collected, a standard operating procedure (SOP) was created as a model system to reduce the inappropriate use of antibiotics in children with ILI. Monte Carlo simulations were performed to assess cost-effectiveness. RESULTS: Among 2,569 ChILD Cohort patients enrolled from 12/2010 to 04/2013 (55\\% male, mean age 3.2 years, range 0-18, 19\\% {\ensuremath{>}}5 years), 411 patients showed laboratory-confirmed influenza, with bacterial co-detection in 35\\%. Influenza and pneumococcus were detected simultaneously in 12/2,569 patients, with disease severity clearly below average. Pneumococcal vaccination rates were close to 90\\%. Nonetheless, every fifth patient was already on antibiotics upon presentation; new antibiotic prescriptions were issued in an additional 20\\%. Simulation of the model SOP in the same dataset revealed that the proposed decision model could have reduced the inappropriate use of antibiotics significantly (P{\ensuremath{<}}0.01) with an incremental cost-effectiveness ratio of -99.55?. CONCLUSIONS: Physicians should be made aware that in times of pneumococcal vaccination the prevalence and severity of influenza infections complicated by pneumococci may decline. Microbiological testing in combination with standardized disease severity assessments and review of vaccination records could be cost-effective, as well as promoting stringent use of antibiotics and a personalized approach to managing children with ILI.}, language = {en} } @article{ObermeierMuehlhansHoppeetal.2016, author = {Obermeier, Patrick and Muehlhans, S. and Hoppe, Christian and Karsch, K. and Tief, Franziska and Seeber, L. and Chen, X. and Conrad, Tim and Boettcher, Sindy and Diedrich, S. and Rath, Barbara}, title = {Enabling Precision Medicine With Digital Case Classification at the Point-of-Care}, volume = {4}, journal = {EBioMedicine}, doi = {10.1016/j.ebiom.2016.01.008}, pages = {191 -- 196}, year = {2016}, abstract = {Infectious and inflammatory diseases of the central nervous system are difficult to identify early. Case definitions for aseptic meningitis, encephalitis, myelitis, and acute disseminated encephalomyelitis (ADEM) are available, but rarely put to use. The VACC-Tool (Vienna Vaccine Safety Initiative Automated Case Classification-Tool) is a mobile application enabling immediate case ascertainment based on consensus criteria at the point-of-care. The VACC-Tool was validated in a quality management program in collaboration with the Robert-Koch-Institute. Results were compared to ICD-10 coding and retrospective analysis of electronic health records using the same case criteria. Of 68,921 patients attending the emergency room in 10/2010-06/2013, 11,575 were hospitalized, with 521 eligible patients (mean age: 7.6 years) entering the quality management program. Using the VACC-Tool at the point-of-care, 180/521 cases were classified successfully and 194/521 ruled out with certainty. Of the 180 confirmed cases, 116 had been missed by ICD-10 coding, 38 misclassified. By retrospective application of the same case criteria, 33 cases were missed. Encephalitis and ADEM cases were most likely missed or misclassified. The VACC-Tool enables physicians to ask the right questions at the right time, thereby classifying cases consistently and accurately, facilitating translational research. Future applications will alert physicians when additional diagnostic procedures are required.}, language = {en} } @article{ConradBrucknerKayser2013, author = {Conrad, Tim and Bruckner, Sharon and Kayser, Bastian}, title = {Finding Modules in Networks with Non-modular Regions}, volume = {7933}, journal = {Lecture Notes in Computer Science (Proceedings of SEA 2013)}, doi = {10.1007/978-3-642-38527-8_18}, pages = {188 -- 199}, year = {2013}, abstract = {Most network clustering methods share the assumption that the network can be completely decomposed into modules, that is, every node belongs to (usually exactly one) module. Forcing this constraint can lead to misidentification of modules where none exist, while the true modules are drowned out in the noise, as has been observed e.g. for protein interaction networks. We thus propose a clustering model where networks contain both a modular region consisting of nodes that can be partitioned into modules, and a transition region containing nodes that lie between or outside modules. We propose two scores based on spectral properties to determine how well a network fits this model. We then evaluate three (partially adapted) clustering algorithms from the literature on random networks that fit our model, based on the scores and comparison to the ground truth. This allows to pinpoint the types of networks for which the different algorithms perform well.}, language = {en} } @article{ConradRathTiefetal.2013, author = {Conrad, Tim and Rath, Barbara and Tief, Franziska and Karsch, K. and Muehlhans, S. and Obermeier, Patrick and Adamou, E. and Chen, X. and Seeber, L. and Peiser, Ch. and Hoppe, Christian and von Kleist, Max and Schweiger, Brunhilde}, title = {Towards a personalized approach to managing of influenza infections in infants and children - food for thought and a note on oseltamivir}, volume = {13}, journal = {Infectious Disorders - Drug Targets}, number = {1}, pages = {25 -- 33}, year = {2013}, language = {en} } @article{ConradLeichtleNuofferetal.2012, author = {Conrad, Tim and Leichtle, Alexander Benedikt and Nuoffer, Jean-Marc and Ceglarek, Uta and Kase, Julia and Witzigmann, Helmut and Thiery, Joachim and Fiedler, Georg Martin}, title = {Serum amino acid profiles and their alterations in colorectal cancer}, journal = {Metabolomics}, doi = {10.1007/s11306-011-0357-5}, year = {2012}, abstract = {Mass spectrometry-based serum metabolic profiling is a promising tool to analyse complex cancer associated metabolic alterations, which may broaden our pathophysiological understanding of the disease and may function as a source of new cancer-associated biomarkers. Highly standardized serum samples of patients suffering from colon cancer (n = 59) and controls (n = 58) were collected at the University Hospital Leipzig. We based our investigations on amino acid screening profiles using electrospray tandem-mass spectrometry. Metabolic profiles were evaluated using the Analyst 1.4.2 software. General, comparative and equivalence statistics were performed by R 2.12.2. 11 out of 26 serum amino acid concentrations were significantly different between colorectal cancer patients and healthy controls. We found a model including CEA, glycine, and tyrosine as best discriminating and superior to CEA alone with an AUROC of 0.878 (95\\% CI 0.815?0.941). Our serum metabolic profiling in colon cancer revealed multiple significant disease-associated alterations in the amino acid profile with promising diagnostic power. Further large-scale studies are necessary to elucidate the potential of our model also to discriminate between cancer and potential differential diagnoses. In conclusion, serum glycine and tyrosine in combination with CEA are superior to CEA for the discrimination between colorectal cancer patients and controls.}, language = {en} } @article{GuptaConradSpoetteretal.2012, author = {Gupta, Pooja and Conrad, Tim and Sp{\"o}tter, Andreas and Reinsch, Norbert and Bienefeld, Kaspar}, title = {Simulating a base population in honey bee for molecular genetic studies}, journal = {Genetics Selection Evolution}, doi = {10.1186/1297-9686-44-14}, year = {2012}, abstract = {Over the past years, reports have indicated that honey bee populations are declining and that infestation by an ecto-parasitic mite (Varroa destructor) is one of the main causes. Selective breeding of resistant bees can help to prevent losses due to the parasite, but it requires that a robust breeding program and genetic evaluation are implemented. Genomic selection has emerged as an important tool in animal breeding programs and simulation studies have shown that it yields more accurate breeding values estimates, higher genetic gain and low rates of inbreeding. Since genomic selection relies on marker data, simulations conducted on a genomic dataset are a pre-requisite before selection can be implemented. Although genomic datasets have been simulated in other species undergoing genetic evaluation, simulation of a genomic dataset specific to the honey bee is required since this species has distinct genetic and reproductive biology characteristics. Our software program was aimed at constructing a base population by simulating a random mating honey bee population. A forward-time population simulation approach was applied since it allows modeling of genetic characteristics and reproductive behavior specific to the honey bee.  Results: Our software program yielded a genomic dataset for a base population in linkage disequilibrium. In addition, information was obtained on (1) the position of markers on each chromosome, (2) allele frequency, (3) ?2 statistics for Hardy- Weinberg equilibrium, (4) a sorted list of markers with a minor allele frequency less than or equal to the input value, (5) average r2 values of linkage disequilibrium between all simulated marker loci pair for all generations and (6) average r2 value of linkage disequilibrium in the last generation for selected markers with the highest minor allele frequency. Conclusion: We developed a software program that takes into account the genetic and reproductive biology characteristics specific to the honey bee and that can be used to constitute a genomic dataset compatible with the simulation studies necessary to optimize breeding programs. The source code together with an instruction file is freely accessible at http://msproteomics.org/Research/Misc/honeybeepopulationsimulator.html}, language = {en} } @article{Conrad2004, author = {Conrad, Tim}, title = {New Appraches for Visualizing and Analyzing Metabolic Pathways}, journal = {Proceedings of the Second Australian Undergraduate Students? Computing Conference}, year = {2004}, abstract = {Visualizing of metabolic pathways (or networks) has been done by many differentapproaches. In this work, we implemented and tested existing graph layout algorithms, and present a new approach to lay-out medium size metabolic pathways (500-20,000 vertices) by implementing and combining three well known graph lay-out algorithms (high dimension embedding, spring-embedder preprocessing, spring-embedder), through 3D space density analysis facilitated by the Octree technique. For the analysis of the results of metabolic pathways simulations we present two new techniques: rstly, a powerful technique to visualize pathways simulation data was created to unveil and understand concentration ows through metabolic pathways. This was achieved by mapping the color encoded concentration value of every substance from each time step of the simulation to its graphical representation in the layout. By combining all resulting images (from each time step) and displaying them as a movie, many characteristics such as subnetworks, alternative routes through the network, and differences between a modied pathway and its unmodied version can be revealed. Secondly, a new method to detect co-regulated substances in metabolic pathways and to recognize differences between two versions of a pathway, was established. To do this, we transformed the simulation data into a row-based representation, color-coded these rows, and reordered them with respect to similarity by using a Genetic Algorithm variant. From the arising discrete 2-dimensional matrix consisting of concentration values, a continuous 2-dimensional fourier row function was computed. This function can be used to measure properties, such as similarities in a pathway between time steps, or substances, or to detect and evaluate differences between modied versions of the same pathway.}, language = {en} } @article{GelssMateraSchuette2016, author = {Gelß, Patrick and Matera, Sebastian and Sch{\"u}tte, Christof}, title = {Solving the master equation without kinetic Monte Carlo: Tensor train approximations for a CO oxidation model}, volume = {314}, journal = {Journal of Computational Physics}, doi = {10.1016/j.jcp.2016.03.025}, pages = {489 -- 502}, year = {2016}, abstract = {In multiscale modeling of heterogeneous catalytic processes, one crucial point is the solution of a Markovian master equation describing the stochastic reaction kinetics. Usually, this is too high-dimensional to be solved with standard numerical techniques and one has to rely on sampling approaches based on the kinetic Monte Carlo method. In this study we break the curse of dimensionality for the direct solution of the Markovian master equation by exploiting the Tensor Train Format for this purpose. The performance of the approach is demonstrated on a first principles based, reduced model for the CO oxidation on the RuO2(110) surface. We investigate the complexity for increasing system size and for various reaction conditions. The advantage over the stochastic simulation approach is illustrated by a problem with increased}, language = {en} } @article{VegaSchuetteConrad2016, author = {Vega, Iliusi and Sch{\"u}tte, Christof and Conrad, Tim}, title = {Finding metastable states in real-world time series with recurrence networks}, volume = {445}, journal = {Physica A: Statistical Mechanics and its Applications}, doi = {10.1016/j.physa.2015.10.041}, pages = {1 -- 17}, year = {2016}, abstract = {In the framework of time series analysis with recurrence networks, we introduce a self-adaptive method that determines the elusive recurrence threshold and identifies metastable states in complex real-world time series. As initial step, we introduce a way to set the embedding parameters used to reconstruct the state space from the time series. We set them as the ones giving the maximum Shannon entropy of the diagonal line length distribution for the first simultaneous minima of recurrence rate and Shannon entropy. To identify metastable states, as well as the transitions between them, we use a soft partitioning algorithm for module finding which is specifically developed for the case in which a system shows metastability. We illustrate our method with a complex time series example. Finally, we show the robustness of our method for identifying metastable states. Our results suggest that our method is robust for identifying metastable states in complex time series, even when introducing considerable levels of noise and missing data points.}, language = {en} } @article{ConradGenzelCvetkovicetal.2017, author = {Conrad, Tim and Genzel, Martin and Cvetkovic, Nada and Wulkow, Niklas and Leichtle, Alexander Benedikt and Vybiral, Jan and Kytyniok, Gitta and Sch{\"u}tte, Christof}, title = {Sparse Proteomics Analysis - a compressed sensing-based approach for feature selection and classification of high-dimensional proteomics mass spectrometry data}, volume = {18}, journal = {BMC Bioinfomatics}, number = {160}, doi = {10.1186/s12859-017-1565-4}, year = {2017}, abstract = {Background: High-throughput proteomics techniques, such as mass spectrometry (MS)-based approaches, produce very high-dimensional data-sets. In a clinical setting one is often interested in how mass spectra differ between patients of different classes, for example spectra from healthy patients vs. spectra from patients having a particular disease. Machine learning algorithms are needed to (a) identify these discriminating features and (b) classify unknown spectra based on this feature set. Since the acquired data is usually noisy, the algorithms should be robust against noise and outliers, while the identified feature set should be as small as possible. Results: We present a new algorithm, Sparse Proteomics Analysis (SPA),based on thet heory of compressed sensing that allows us to identify a minimal discriminating set of features from mass spectrometry data-sets. We show (1) how our method performs on artificial and real-world data-sets, (2) that its performance is competitive with standard (and widely used) algorithms for analyzing proteomics data, and (3) that it is robust against random and systematic noise. We further demonstrate the applicability of our algorithm to two previously published clinical data-sets.}, language = {en} } @article{GuptaGramatkeEinspanieretal.2017, author = {Gupta, Pooja and Gramatke, Annika and Einspanier, Ralf and Sch{\"u}tte, Christof and von Kleist, Max and Sharbati, Jutta}, title = {In silico cytotoxicity assessment on cultured rat intestinal cells deduced from cellular impedance measurements}, volume = {41}, journal = {Toxicology in Vitro}, issn = {1438-0064}, pages = {179 -- 188}, year = {2017}, abstract = {Early and reliable identification of chemical toxicity is of utmost importance. At the same time, reduction of animal testing is paramount. Therefore, methods that improve the interpretability and usability of in vitro assays are essential. xCELLigence's real-time cell analyzer (RTCA) provides a novel, fast and cost effective in vitro method to probe compound toxicity. We developed a simple mathematical framework for the qualitative and quantitative assessment of toxicity for RTCA measurements. Compound toxicity, in terms of its 50\% inhibitory concentration IC50 on cell growth, and parameters related to cell turnover were estimated on cultured IEC-6 cells exposed to 10 chemicals at varying concentrations. Our method estimated IC50 values of 113.05, 7.16, 28.69 and 725.15 μM for the apparently toxic compounds 2-acetylamino-fluorene, aflatoxin B1, benzo-[a]-pyrene and chloramphenicol in the tested cell line, in agreement with literature knowledge. IC50 values of all apparent in vivo non-toxic compounds were estimated to be non-toxic by our method. Corresponding estimates from RTCA's in-built model gave false positive (toxicity) predictions in 5/10 cases. Taken together, our proposed method reduces false positive predictions and reliably identifies chemical toxicity based on impedance measurements. The source code for the developed method including instructions is available at https://git.zib.de/bzfgupta/toxfit/tree/master.}, language = {en} } @article{WilsonAnglinAmbellanetal.2017, author = {Wilson, David and Anglin, Carolyn and Ambellan, Felix and Grewe, Carl Martin and Tack, Alexander and Lamecker, Hans and Dunbar, Michael and Zachow, Stefan}, title = {Validation of three-dimensional models of the distal femur created from surgical navigation point cloud data for intraoperative and postoperative analysis of total knee arthroplasty}, volume = {12}, journal = {International Journal of Computer Assisted Radiology and Surgery}, number = {12}, publisher = {Springer}, doi = {10.1007/s11548-017-1630-5}, pages = {2097 -- 2105}, year = {2017}, abstract = {Purpose: Despite the success of total knee arthroplasty there continues to be a significant proportion of patients who are dissatisfied. One explanation may be a shape mismatch between pre and post-operative distal femurs. The purpose of this study was to investigate a method to match a statistical shape model (SSM) to intra-operatively acquired point cloud data from a surgical navigation system, and to validate it against the pre-operative magnetic resonance imaging (MRI) data from the same patients. Methods: A total of 10 patients who underwent navigated total knee arthroplasty also had an MRI scan less than 2 months pre-operatively. The standard surgical protocol was followed which included partial digitization of the distal femur. Two different methods were employed to fit the SSM to the digitized point cloud data, based on (1) Iterative Closest Points (ICP) and (2) Gaussian Mixture Models (GMM). The available MRI data were manually segmented and the reconstructed three-dimensional surfaces used as ground truth against which the statistical shape model fit was compared. Results: For both approaches, the difference between the statistical shape model-generated femur and the surface generated from MRI segmentation averaged less than 1.7 mm, with maximum errors occurring in less clinically important areas. Conclusion: The results demonstrated good correspondence with the distal femoral morphology even in cases of sparse data sets. Application of this technique will allow for measurement of mismatch between pre and post-operative femurs retrospectively on any case done using the surgical navigation system and could be integrated into the surgical navigation unit to provide real-time feedback.}, language = {en} } @article{vonTycowiczAmbellanMukhopadhyayetal.2018, author = {von Tycowicz, Christoph and Ambellan, Felix and Mukhopadhyay, Anirban and Zachow, Stefan}, title = {An Efficient Riemannian Statistical Shape Model using Differential Coordinates}, volume = {43}, journal = {Medical Image Analysis}, number = {1}, doi = {10.1016/j.media.2017.09.004}, pages = {1 -- 9}, year = {2018}, abstract = {We propose a novel Riemannian framework for statistical analysis of shapes that is able to account for the nonlinearity in shape variation. By adopting a physical perspective, we introduce a differential representation that puts the local geometric variability into focus. We model these differential coordinates as elements of a Lie group thereby endowing our shape space with a non-Euclidean structure. A key advantage of our framework is that statistics in a manifold shape space becomes numerically tractable improving performance by several orders of magnitude over state-of-the-art. We show that our Riemannian model is well suited for the identification of intra-population variability as well as inter-population differences. In particular, we demonstrate the superiority of the proposed model in experiments on specificity and generalization ability. We further derive a statistical shape descriptor that outperforms the standard Euclidean approach in terms of shape-based classification of morphological disorders.}, language = {en} } @article{BennHiepenOsterlandetal.2017, author = {Benn, Andreas and Hiepen, Christian and Osterland, Marc and Sch{\"u}tte, Christof and Zwijsen, An and Knaus, Petra}, title = {Role of bone morphogenetic proteins in sprouting angiogenesis: differential BMP receptor-dependent signaling pathways balance stalk vs. tip cell competence}, volume = {31}, journal = {FASEB Journal}, number = {11}, doi = {10.1096/fj.201700193RR}, pages = {4720 -- 4733}, year = {2017}, abstract = {Before the onset of sprouting angiogenesis, the endothelium is prepatterned for the positioning of tip and stalk cells. Both cell identities are not static, as endothelial cells (ECs) constantly compete for the tip cell position in a dynamic fashion. Here, we show that both bone morphogenetic protein (BMP) 2 and BMP6 are proangiogenic in vitro and ex vivo and that the BMP type I receptors, activin receptor-like kinase (ALK)3 and ALK2, play crucial and distinct roles in this process. BMP2 activates the expression of tip cell-associated genes, such as DLL4 (delta-like ligand 4) and KDR (kinase insert domain receptor), and p38-heat shock protein 27 (HSP27)-dependent cell migration, thereby generating tip cell competence. Whereas BMP6 also triggers collective cell migration via the p38-HSP27 signaling axis, BMP6 induces in addition SMAD1/5 signaling, thereby promoting the expression of stalk cell-associated genes, such as HES1 (hairy and enhancer of split 1) and FLT1 (fms-like tyrosine kinase 1). Specifically, ALK3 is required for sprouting from HUVEC spheroids, whereas ALK2 represses sprout formation. We demonstrate that expression levels and respective complex formation of BMP type I receptors in ECs determine stalk vs. tip cell identity, thus contributing to endothelial plasticity during sprouting angiogenesis. As antiangiogenic monotherapies that target the VEGF or ALK1 pathways have not fulfilled efficacy objectives in clinical trials, the selective targeting of the ALK2/3 pathways may be an attractive new approach.}, language = {en} } @article{WinkelmannSchuette2017, author = {Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {Hybrid models for chemical reaction networks: Multiscale theory and application to gene regulatory systems}, volume = {147}, journal = {The Journal of Chemical Physics}, number = {11}, doi = {10.1063/1.4986560}, pages = {114115-1 -- 114115-18}, year = {2017}, abstract = {Well-mixed stochastic chemical kinetics are properly modeled by the chemical master equation (CME) and associated Markov jump processes in molecule number space. If the reactants are present in large amounts, however, corresponding simulations of the stochastic dynamics become computationally expensive and model reductions are demanded. The classical model reduction approach uniformly rescales the overall dynamics to obtain deterministic systems characterized by ordinary differential equations, the well-known mass action reaction rate equations. For systems with multiple scales, there exist hybrid approaches that keep parts of the system discrete while another part is approximated either using Langevin dynamics or deterministically. This paper aims at giving a coherent overview of the different hybrid approaches, focusing on their basic concepts and the relation between them. We derive a novel general description of such hybrid models that allows expressing various forms by one type of equation. We also check in how far the approaches apply to model extensions of the CME for dynamics which do not comply with the central well-mixed condition and require some spatial resolution. A simple but meaningful gene expression system with negative self-regulation is analysed to illustrate the different approximation qualities of some of the hybrid approaches discussed. Especially, we reveal the cause of error in the case of small volume approximations.}, language = {en} } @misc{SchuetteConrad2014, author = {Sch{\"u}tte, Christof and Conrad, Tim}, title = {Showcase 3: Information-based medicine}, volume = {1}, journal = {MATHEON-Mathematics for Key Technologies}, editor = {Deuflhard, Peter and Gr{\"o}tschel, Martin and H{\"o}mberg, Dietmar and Horst, Ulrich and Kramer, J{\"u}rg and Mehrmann, Volker and Polthier, Konrad and Schmidt, Frank and Skutella, Martin and Sprekels, J{\"u}rgen}, publisher = {European Mathematical Society}, pages = {66 -- 67}, year = {2014}, language = {en} } @article{KryvenRoeblitzSchuette2015, author = {Kryven, Ivan and R{\"o}blitz, Susanna and Sch{\"u}tte, Christof}, title = {Solution of the chemical master equation by radial basis functions approximation with interface tracking}, volume = {9}, journal = {BMC Systems Biology}, number = {67}, doi = {10.1186/s12918-015-0210-y}, pages = {1 -- 12}, year = {2015}, abstract = {Background. The chemical master equation is the fundamental equation of stochastic chemical kinetics. This differential-difference equation describes temporal evolution of the probability density function for states of a chemical system. A state of the system, usually encoded as a vector, represents the number of entities or copy numbers of interacting species, which are changing according to a list of possible reactions. It is often the case, especially when the state vector is high-dimensional, that the number of possible states the system may occupy is too large to be handled computationally. One way to get around this problem is to consider only those states that are associated with probabilities that are greater than a certain threshold level. Results. We introduce an algorithm that significantly reduces computational resources and is especially powerful when dealing with multi-modal distributions. The algorithm is built according to two key principles. Firstly, when performing time integration, the algorithm keeps track of the subset of states with significant probabilities (essential support). Secondly, the probability distribution that solves the equation is parametrised with a small number of coefficients using collocation on Gaussian radial basis functions. The system of basis functions is chosen in such a way that the solution is approximated only on the essential support instead of the whole state space. Discussion. In order to demonstrate the effectiveness of the method, we consider four application examples: a) the self-regulating gene model, b) the 2-dimensional bistable toggle switch, c) a generalisation of the bistable switch to a 3-dimensional tristable problem, and d) a 3-dimensional cell differentiation model that, depending on parameter values, may operate in bistable or tristable modes. In all multidimensional examples the manifold containing the system states with significant probabilities undergoes drastic transformations over time. This fact makes the examples especially challenging for numerical methods. Conclusions. The proposed method is a new numerical approach permitting to approximately solve a wide range of problems that have been hard to tackle until now. A full representation of multi-dimensional distributions is recovered. The method is especially attractive when dealing with models that yield solutions of a complex structure, for instance, featuring multi-stability. Electronic version: http://www.biomedcentral.com/1752-0509/9/67}, language = {en} } @incollection{LameckerZachow2016, author = {Lamecker, Hans and Zachow, Stefan}, title = {Statistical Shape Modeling of Musculoskeletal Structures and Its Applications}, volume = {23}, booktitle = {Computational Radiology for Orthopaedic Interventions}, publisher = {Springer}, isbn = {978-3-319-23481-6}, doi = {10.1007/978-3-319-23482-3}, pages = {1 -- 23}, year = {2016}, abstract = {Statistical shape models (SSM) describe the shape variability contained in a given population. They are able to describe large populations of complex shapes with few degrees of freedom. This makes them a useful tool for a variety of tasks that arise in computer-aided madicine. In this chapter we are going to explain the basic methodology of SSMs and present a variety of examples, where SSMs have been successfully applied.}, language = {en} } @inproceedings{MukhopadhyayOksuzBevilacquaetal.2015, author = {Mukhopadhyay, Anirban and Oksuz, Ilkay and Bevilacqua, Marco and Dharmakumar, Rohan and Tsaftaris, Sotirios}, title = {Data-Driven Feature Learning for Myocardial Segmentation of CP-BOLD MRI}, volume = {9126}, booktitle = {Functional Imaging and Modeling of the Heart}, publisher = {Springer}, doi = {10.1007/978-3-319-20309-6_22}, pages = {189 -- 197}, year = {2015}, abstract = {Cardiac Phase-resolved Blood Oxygen-Level-Dependent (CP- BOLD) MR is capable of diagnosing an ongoing ischemia by detecting changes in myocardial intensity patterns at rest without any contrast and stress agents. Visualizing and detecting these changes require significant post-processing, including myocardial segmentation for isolating the myocardium. But, changes in myocardial intensity pattern and myocardial shape due to the heart's motion challenge automated standard CINE MR myocardial segmentation techniques resulting in a significant drop of segmentation accuracy. We hypothesize that the main reason behind this phenomenon is the lack of discernible features. In this paper, a multi scale discriminative dictionary learning approach is proposed for supervised learning and sparse representation of the myocardium, to improve the myocardial feature selection. The technique is validated on a challenging dataset of CP-BOLD MR and standard CINE MR acquired in baseline and ischemic condition across 10 canine subjects. The proposed method significantly outperforms standard cardiac segmentation techniques, including segmentation via registration, level sets and supervised methods for myocardial segmentation.}, language = {en} } @inproceedings{MukhopadhyayOksuzBevilacquaetal.2015, author = {Mukhopadhyay, Anirban and Oksuz, Ilkay and Bevilacqua, Marco and Dharmakumar, Rohan and Tsaftaris, Sotirios}, title = {Unsupervised myocardial segmentation for cardiac MRI}, volume = {LNCS 9351}, booktitle = {Medical Image Computing and Computer-Assisted Intervention -- MICCAI 2015}, doi = {10.1007/978-3-319-24574-4_2}, pages = {12 -- 20}, year = {2015}, abstract = {Though unsupervised segmentation was a de-facto standard for cardiac MRI segmentation early on, recently cardiac MRI segmentation literature has favored fully supervised techniques such as Dictionary Learning and Atlas-based techniques. But, the benefits of unsupervised techniques e.g., no need for large amount of training data and better potential of handling variability in anatomy and image contrast, is more evident with emerging cardiac MR modalities. For example, CP-BOLD is a new MRI technique that has been shown to detect ischemia without any contrast at stress but also at rest conditions. Although CP-BOLD looks similar to standard CINE, changes in myocardial intensity patterns and shape across cardiac phases, due to the heart's motion, BOLD effect and artifacts affect the underlying mechanisms of fully supervised segmentation techniques resulting in a significant drop in segmentation accuracy. In this paper, we present a fully unsupervised technique for segmenting myocardium from the background in both standard CINE MR and CP-BOLD MR. We combine appearance with motion information (obtained via Optical Flow) in a dictionary learning framework to sparsely represent important features in a low dimensional space and separate myocardium from background accordingly. Our fully automated method learns background-only models and one class classifier provides myocardial segmentation. The advantages of the proposed technique are demonstrated on a dataset containing CP-BOLD MR and standard CINE MR image sequences acquired in baseline and ischemic condition across 10 canine subjects, where our method outperforms state-of-the-art supervised segmentation techniques in CP-BOLD MR and performs at-par for standard CINE MR.}, language = {en} } @inproceedings{OksuzMukhopadhyayBevilacquaetal.2015, author = {Oksuz, Ilkay and Mukhopadhyay, Anirban and Bevilacqua, Marco and Dharmakumar, Rohan and Tsaftaris, Sotirios}, title = {Dictionary Learning Based Image Descriptor for Myocardial Registration of CP-BOLD MR}, volume = {9350}, booktitle = {Medical Image Computing and Computer-Assisted Intervention -- MICCAI 2015}, publisher = {Springer}, doi = {10.1007/978-3-319-24571-3_25}, pages = {205 -- 213}, year = {2015}, abstract = {Cardiac Phase-resolved Blood Oxygen-Level-Dependent (CP- BOLD) MRI is a new contrast agent- and stress-free imaging technique for the assessment of myocardial ischemia at rest. The precise registration among the cardiac phases in this cine type acquisition is essential for automating the analysis of images of this technique, since it can potentially lead to better specificity of ischemia detection. However, inconsistency in myocardial intensity patterns and the changes in myocardial shape due to the heart's motion lead to low registration performance for state- of-the-art methods. This low accuracy can be explained by the lack of distinguishable features in CP-BOLD and inappropriate metric defini- tions in current intensity-based registration frameworks. In this paper, the sparse representations, which are defined by a discriminative dictionary learning approach for source and target images, are used to improve myocardial registration. This method combines appearance with Gabor and HOG features in a dictionary learning framework to sparsely represent features in a low dimensional space. The sum of squared differences of these distinctive sparse representations are used to define a similarity term in the registration framework. The proposed descriptor is validated on a challenging dataset of CP-BOLD MR and standard CINE MR acquired in baseline and ischemic condition across 10 canines.}, language = {en} } @misc{OsterlandBennProhaskaetal.2015, author = {Osterland, Marc and Benn, Andreas and Prohaska, Steffen and Sch{\"u}tte, Christof}, title = {Single Cell Tracking in Phase-Contrast Microscopy}, journal = {EMBL Symposium 2015 - Seeing is Believing - Imaging the Processes of Life}, year = {2015}, abstract = {In this work, we developed an automatic algorithm to analyze cell migration in chemotaxis assays, based on phase-contrast time-lapse microscopy. While manual approaches are still widely used in recent publications, our algorithm is able to track hundreds of single cells per frame. The extracted paths are analysed with traditional geometrical approaches as well as diffusion-driven Markov state models (MSM). Based on these models, a detailed view on spatial and temporal effects is possible. Using our new approach on experimental data, we are able to distinguish between directed migration (e.g. towards a VEGF gradient) and random migration without favored direction. A calculation of the committor probabilities reveals that cells of the whole image area are more likely to migrate directly towards the VEGF than away from it during the first four hours. However, in absence of a chemoattractant, cells migrate more likely to their nearest image border. These conclusions are supported by the spatial mean directions. In a next step, the cell-cell interaction during migration and the migration of cell clusters will be analyzed. Furthermore, we want to observe phenotypical changes during migration based on fluorescence microscopy and machine learning. The algorithm is part of a collaborative platform which brings the experimental expertise of scientists from life sciences and the analytical knowledge of computer scientists together. This platform is built using web-based technologies with a responsive real-time user interface. All data, including raw and metadata as well as the accompanying results, will be stored in a secure and scalable compute cluster. The compute cluster provides sufficient space and computational power for modern image-based experiments and their analyses. Specific versions of data and results can be tagged to keep immutable records for archival.}, language = {en} } @misc{KoltaiCiccottiSchuette2016, author = {Koltai, Peter and Ciccotti, Giovanni and Sch{\"u}tte, Christof}, title = {On metastability and Markov state models for non-stationary molecular dynamics}, volume = {174103}, journal = {The Journal of Chemical Physics}, edition = {145}, issn = {1438-0064}, doi = {10.1063/1.4966157}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-57869}, year = {2016}, abstract = {We utilize the theory of coherent sets to build Markov state models for non- equilibrium molecular dynamical systems. Unlike for systems in equilibrium, "meta- stable" sets in the non-equilibrium case may move as time evolves. We formalize this concept by relying on the theory of coherent sets, based on this we derive finite-time non-stationary Markov state models, and illustrate the concept and its main differences to equilibrium Markov state modeling on simple, one-dimensional examples.}, language = {en} } @article{DjurdjevacConradWeberSchuette2016, author = {Djurdjevac Conrad, Natasa and Weber, Marcus and Sch{\"u}tte, Christof}, title = {Finding dominant structures of nonreversible Markov processes}, volume = {14}, journal = {Multiscale Modeling and Simulation}, number = {4}, doi = {10.1137/15M1032272}, pages = {1319 -- 1340}, year = {2016}, language = {en} } @misc{MukhopadhyayKumarBhandarkar2016, author = {Mukhopadhyay, Anirban and Kumar, Arun and Bhandarkar, Suchendra}, title = {Joint Geometric Graph Embedding for Partial Shape Matching in Images}, journal = {IEEE Winter Conference on Applications of Computer Vision}, edition = {IEEE Winter Conference on Applications of Computer Vision (WACV)}, publisher = {IEEE}, pages = {1 -- 9}, year = {2016}, abstract = {A novel multi-criteria optimization framework for matching of partially visible shapes in multiple images using joint geometric graph embedding is proposed. The proposed framework achieves matching of partial shapes in images that exhibit extreme variations in scale, orientation, viewpoint and illumination and also instances of occlusion; conditions which render impractical the use of global contour-based descriptors or local pixel-level features for shape matching. The proposed technique is based on optimization of the embedding distances of geometric features obtained from the eigenspectrum of the joint image graph, coupled with regularization over values of the mean pixel intensity or histogram of oriented gradients. It is shown to obtain successfully the correspondences denoting partial shape similarities as well as correspondences between feature points in the images. A new benchmark dataset is proposed which contains disparate image pairs with extremely challenging variations in viewing conditions when compared to an existing dataset [18]. The proposed technique is shown to significantly outperform several state-of-the-art partial shape matching techniques on both datasets.}, language = {en} } @inproceedings{GuptaKrauseRikeitetal.2014, author = {Gupta, Pooja and Krause, Carola and Rikeit, Paul and R{\"o}blitz, Susanna and Knaus, Petra and Sch{\"u}tte, Christof}, title = {Modeling of the BMP mediated co-regulation of the Smad and Non-Smad pathways in the context of cell density}, booktitle = {10th International BMP conference, 2014, Berlin, Germany}, year = {2014}, language = {en} } @article{PeppertvonKleistSchuetteetal.2022, author = {Peppert, Felix and von Kleist, Max and Sch{\"u}tte, Christof and Sunkara, Vikram}, title = {On the Sufficient Condition for Solving the Gap-Filling Problem Using Deep Convolutional Neural Networks}, volume = {33}, journal = {IEEE Transactions on Neural Networks and Learning Systems}, number = {11}, doi = {10.1109/TNNLS.2021.3072746}, pages = {6194 -- 6205}, year = {2022}, abstract = {Deep convolutional neural networks (DCNNs) are routinely used for image segmentation of biomedical data sets to obtain quantitative measurements of cellular structures like tissues. These cellular structures often contain gaps in their boundaries, leading to poor segmentation performance when using DCNNs like the U-Net. The gaps can usually be corrected by post-hoc computer vision (CV) steps, which are specific to the data set and require a disproportionate amount of work. As DCNNs are Universal Function Approximators, it is conceivable that the corrections should be obsolete by selecting the appropriate architecture for the DCNN. In this article, we present a novel theoretical framework for the gap-filling problem in DCNNs that allows the selection of architecture to circumvent the CV steps. Combining information-theoretic measures of the data set with a fundamental property of DCNNs, the size of their receptive field, allows us to formulate statements about the solvability of the gap-filling problem independent of the specifics of model training. In particular, we obtain mathematical proof showing that the maximum proficiency of filling a gap by a DCNN is achieved if its receptive field is larger than the gap length. We then demonstrate the consequence of this result using numerical experiments on a synthetic and real data set and compare the gap-filling ability of the ubiquitous U-Net architecture with variable depths. Our code is available at https://github.com/ai-biology/dcnn-gap-filling.}, language = {en} } @article{LiPimentelSzengeletal.2021, author = {Li, Jianning and Pimentel, Pedro and Szengel, Angelika and Ehlke, Moritz and Lamecker, Hans and Zachow, Stefan and Estacio, Laura and Doenitz, Christian and Ramm, Heiko and Shi, Haochen and Chen, Xiaojun and Matzkin, Franco and Newcombe, Virginia and Ferrante, Enzo and Jin, Yuan and Ellis, David G. and Aizenberg, Michele R. and Kodym, Oldrich and Spanel, Michal and Herout, Adam and Mainprize, James G. and Fishman, Zachary and Hardisty, Michael R. and Bayat, Amirhossein and Shit, Suprosanna and Wang, Bomin and Liu, Zhi and Eder, Matthias and Pepe, Antonio and Gsaxner, Christina and Alves, Victor and Zefferer, Ulrike and von Campe, Cord and Pistracher, Karin and Sch{\"a}fer, Ute and Schmalstieg, Dieter and Menze, Bjoern H. and Glocker, Ben and Egger, Jan}, title = {AutoImplant 2020 - First MICCAI Challenge on Automatic Cranial Implant Design}, volume = {40}, journal = {IEEE Transactions on Medical Imaging}, number = {9}, issn = {0278-0062}, doi = {10.1109/TMI.2021.3077047}, pages = {2329 -- 2342}, year = {2021}, abstract = {The aim of this paper is to provide a comprehensive overview of the MICCAI 2020 AutoImplant Challenge. The approaches and publications submitted and accepted within the challenge will be summarized and reported, highlighting common algorithmic trends and algorithmic diversity. Furthermore, the evaluation results will be presented, compared and discussed in regard to the challenge aim: seeking for low cost, fast and fully automated solutions for cranial implant design. Based on feedback from collaborating neurosurgeons, this paper concludes by stating open issues and post-challenge requirements for intra-operative use.}, language = {en} } @article{DaiFuellgrabePfeufferetal.2021, author = {Dai, Chengxin and F{\"u}llgrabe, Anja and Pfeuffer, Julianus and Solovyeva, Elizaveta M. and Deng, Jingwen and Moreno, Pablo and Kamatchinathan, Selvakumar and Kundu, Deepti Jaiswal and George, Nancy and Fexovy, Silvie and Gr{\"u}ning, Bj{\"o}rn and F{\"o}ll, Melanie Christine and Griss, Johannes and Vaudel, Marc and Audain, Enrique and Locard-Paulet, Marie and Turewicz, Michael and Eisenacher, Martin and Uszkoreit, Julian and Van Den Bossche, Tim and Schw{\"a}mmle, Veit and Webel, Henry and Schulze, Stefan and Bouyssi{\´e}, David and Jayaram, Savita and Duggineni, Vinay Kumar and Samaras, Patroklos and Wilhelm, Mathias and Choi, Meena and Wang, Mingxun and Kohlbacher, Oliver and Brazma, Alvis and Papatheodorou, Irene and Bandeira, Nuno and Deutsch, Eric W. and Vizca{\´i}no, Juan Antonio and Bai, Mingze and Sachsenberg, Timo and Levitsky, Lev I. and Perez-Riverol, Yasset}, title = {A proteomics sample metadata representation for multiomics integration and big data analysis}, volume = {12}, journal = {Nature Communications}, number = {5854}, doi = {https://doi.org/10.1038/s41467-021-26111-3}, year = {2021}, abstract = {The amount of public proteomics data is rapidly increasing but there is no standardized format to describe the sample metadata and their relationship with the dataset files in a way that fully supports their understanding or reanalysis. Here we propose to develop the transcriptomics data format MAGE-TAB into a standard representation for proteomics sample metadata. We implement MAGE-TAB-Proteomics in a crowdsourcing project to manually curate over 200 public datasets. We also describe tools and libraries to validate and submit sample metadata-related information to the PRIDE repository. We expect that these developments will improve the reproducibility and facilitate the reanalysis and integration of public proteomics datasets.}, language = {en} } @article{UmerZhuPfeufferetal.2021, author = {Umer, Husen M. and Zhu, Yafeng and Pfeuffer, Julianus and Sachsenberg, Timo and Lehti{\"o}, Janne and Branca, Rui and Perez-Riverol, Yasset}, title = {Generation of ENSEMBL-based proteogenomics databases boosts the identification of non-canonical peptides}, journal = {Bioinformatics}, number = {5}, edition = {38}, publisher = {Oxford Academic}, pages = {1470 -- 1472}, year = {2021}, abstract = {We have implemented the pypgatk package and the pgdb workflow to create proteogenomics databases based on ENSEMBL resources. The tools allow the generation of protein sequences from novel protein-coding transcripts by performing a three-frame translation of pseudogenes, lncRNAs, and other non-canonical transcripts, such as those produced by alternative splicing events. It also includes exonic out-of-frame translation from otherwise canonical protein-coding mRNAs. Moreover, the tool enables the generation of variant protein sequences from multiple sources of genomic variants including COSMIC, cBioportal, gnomAD, and mutations detected from sequencing of patient samples. pypgatk and pgdb provide multiple functionalities for database handling, notably optimized target/decoy generation by the algorithm DecoyPyrat. Finally, we perform a reanalysis of four public datasets in PRIDE by generating cell-type specific databases for 65 cell lines using the pypgatk and pgdb workflow, revealing a wealth of non-canonical or cryptic peptides amounting to more than 10\% of the total number of peptides identified (43,501 out of 402,512).}, language = {en} }