Simultaneous Multiscale Polyaffine Registration by Incorporating Deformation Statistics

Locally affine (polyaffine) image registration methods capture intersubject non-linear deformations with a low number of parameters, while providing an intuitive interpretation for clinicians. Considering the mandible bone, anatomical shape differences can be found at different scales, e.g. left or right side, teeth, etc. Classically, sequential coarse to fine registration are used to handle multiscale deformations, instead we propose a simultaneous optimization of all scales. To avoid local minima we incorporate a prior on the polyaffine transformations. This kind of groupwise registration approach is natural in a polyaffine context, if we assume one configuration of regions that describes an entire group of images, with varying transformations for each region. In this paper, we reformulate polyaffine deformations in a generative statistical model, which enables us to incorporate deformation statistics as a prior in a Bayesian setting. We find optimal transformations by optimizing the maximum a posteriori probability. We assume that the polyaffine transformations follow a normal distribution with mean and concentration matrix. Parameters of the prior are estimated from an initial coarse to fine registration. Knowing the region structure, we develop a blockwise pseudoinverse to obtain the concentration matrix. To our knowledge, we are the first to introduce simultaneous multiscale optimization through groupwise polyaffine registration. We show results on 42 mandible CT images.

The simultaneous role of an alveolus as flow mixer and flow feeder for the deposition of inhaled submicron particles

In an effort to understand the fate of inhaled submicron particles in the small sacs, or alveoli, comprising the gas-exchange region of the lung, we calculated the flow in three-dimensional (3D) rhythmically expanding models of alveolated ducts. Since convection toward the alveolar walls is a precursor to particle deposition, it was the goal of this paper to investigate the streamline maps' dependence upon alveoli location along the acinar tree. On the alveolar midplane, the recirculating flow pattern exhibited closed streamlines with a stagnation saddle point. Off the midplane we found no closed streamlines but nested, funnel-like, spiral, structures (reminiscent of Russian nesting dolls) that were directed towards the expanding walls in inspiration, and away from the contracting walls in expiration. These nested, funnel-like, structures were surrounded by air that flowed into the cavity from the central channel over inspiration and flowed from the cavity to the central channel over expiration. We also found that fluid particle tracks exhibited similar nested funnel-like spiral structures. We conclude that these unique alveolar flow structures may be of importance in enhancing deposition. In addition, due to inertia, the nested, funnel-like, structures change shape and position slightly during a breathing cycle, resulting in flow mixing. Also, each inspiration feeds a fresh supply of particle-laden air from the central channel to the region surrounding the mixing region. Thus, this combination of flow mixer and flow feeder makes each individual alveolus an effective mixing unit, which is likely to play an important role in determining the overall efficiency of convective mixing in the acinus.

Epileptogenic developmental venous anomaly: insights from simultaneous EEG/fMRI

Developmental venous anomalies (DVAs) are associated with epileptic seizures; however, the role of DVA in the epileptogenesis is still not established. Simultaneous interictal electroencephalogram/functional magnetic resonance imaging (EEG/fMRI) recordings provide supplementary information to electroclinical data about the epileptic generators, and thus aid in the differentiation of clinically equivocal epilepsy syndromes. The main objective of our study was to characterize the epileptic network in a patient with DVA and epilepsy by simultaneous EEG/fMRI recordings. A 17-year-old woman with recently emerging generalized tonic-clonic seizures, and atypical generalized discharges, was investigated using simultaneous EEG/fMRI at the university hospital. Previous high-resolution MRI showed no structural abnormalities, except a DVA in the right frontal operculum. Interictal EEG recordings showed atypical generalized discharges, corresponding to positive focal blood oxygen level dependent (BOLD) correlates in the right frontal operculum, a region drained by the DVA. Additionally, widespread cortical bilateral negative BOLD correlates in the frontal and parietal lobes were delineated, resembling a generalized epileptic network. The EEG/fMRI recordings support a right frontal lobe epilepsy, originating in the vicinity of the DVA, propagating rapidly to both frontal and parietal lobes, as expressed on the scalp EEG by secondary bilateral synchrony. The DVA may be causative of focal epilepsies in cases where no concomitant epileptogenic lesions can be detected. Advanced imaging techniques, such as simultaneous EEG/fMRI, may thus aid in the differentiation of clinically equivocal epilepsy syndromes.

Determinants of prognostically relevant intracoronary electrocardiogram ST-segment shift during coronary balloon occlusion

The prognostic relevance of quantitative an intracoronary occlusive electrocardiographic (ECG) ST-segment shift and its determinants have not been investigated in humans. In 765 patients with chronic stable coronary artery disease, the following simultaneous quantitative measurements were obtained during a 1-minute coronary balloon occlusion: intracoronary ECG ST-segment shift (recorded by angioplasty guidewire), mean aortic pressure, mean distal coronary pressure, and mean central venous pressure (CVP). Collateral flow index (CFI) was calculated as follows: (mean distal coronary pressure minus CVP)/(mean aortic pressure minus CVP). During an average follow-up duration of 50 ± 34 months, the cumulative mortality rate from all causes was significantly lower in the group with an ST-segment shift <0.1 mV (n = 89) than in the group with an ST-segment shift ≥0.1 mV (n = 676, p = 0.0211). Factors independently related to intracoronary occlusive ECG ST-segment shift <0.1 mV (r(2) = 0.189, p <0.0001) were high CFI (p <0.0001), intracoronary occlusive RR interval (p = 0.0467), right coronary artery as the ischemic region (p <0.0001), and absence of arterial hypertension (p = 0.0132). "High" CFI according to receiver operating characteristics analysis was ≥0.217 (area under receiver operating characteristics curve 0.647, p <0.0001). In conclusion, absence of ECG ST-segment shift during brief coronary occlusion in patients with chronic coronary artery disease conveys a decreased mortality and is directly influenced by a well-developed collateral supply to the right versus left coronary ischemic region and by the absence of systemic hypertension in a patient's history.

Is there a role for simultaneous hepatic and colorectal resections? A contemporary view from NSQIP

The optimal timing of primary and metastatic tumor management in patients with synchronous hepatic colorectal metastases remains controversial. We aimed to compare perioperative outcomes of simultaneous colorectal/liver resection (SCLR) with isolated resections utilizing a national clinical database.

Validated capillary electrophoresis assay for the simultaneous enantioselective determination of propafenone and its major metabolites in biological samples

A robust, inexpensive, and fully validated CE method for the simultaneous determination of the enantiomers of propafenone (PPF), 5-hydroxy-propafenone (5OH-PPF) and N-despropyl-propafenone (NOR-PPF) in serum and in in vitro media is described. It is based upon liquid-liquid extraction at alkaline pH followed by analysis of the reconstituted extract by CE in presence of a pH 2.0 running buffer composed of 100 mM sodium phosphate, 19% methanol, and 0.6% highly sulfated beta-CD. For each compound, the S-enantiomers are shown to migrate ahead of their antipodes, and the overall run time is about 30 min. Enantiomer levels between 25 and 1000 ng/mL provide linear calibration graphs, and the LOD for all enantiomers is between 10 and 12 ng/mL. The assay is shown to be suitable for the determination of the enantiomers of PPF and its metabolites in in vitro incubations comprising human liver microsomes or single CYP450 enzymes (SUPERSOMES). Incubations with CYP2D6 SUPERSOMES revealed, for the first time, the simultaneous formation of the enantiomers of 5OH-PPF and NOR-PPF with that enzyme. CE data can be used for the evaluation of the enzymatic N-dealkylation and hydroxylation rates.

[Simultaneous defect of visual fields and loss of libido--a coincidence?]

We present the case of a 60 year old male patient with incidentally detected visual abnormalities. Detailed personal history revealed a hypogonadism that had been present for several years. Further investigations established the diagnosis of an infiltrative macroadenoma. Medical treatment with cabergoline led to a rapid regression of ophthalmologic symptoms and, subsequently, of tumor size. In male subjects symptoms of hypogonadism are often reported only late in the course of the disease, thereby leading to a generally larger tumor size at the point of diagnosis. In contrast to other pituitary tumors that are mainly treated by surgery, medical treatment with dopamine agonists is the principal therapeutic option in prolactinomas.

QTc interval and resting heart rate as long-term predictors of mortality in type 1 and type 2 diabetes mellitus: a 23-year follow-up

We evaluated the association of QT interval corrected for heart rate (QT(c)) and resting heart rate (rHR) with mortality (all-causes, cardiovascular, cardiac, and ischaemic heart disease) in subjects with type 1 and type 2 diabetes.

Estimation with Interval Censored Data in Longitudinal Studies

In biostatistical applications interest often focuses on the estimation of the distribution of a time-until-event variable T. If one observes whether or not T exceeds an observed monitoring time at a random number of monitoring times, then the data structure is called interval censored data. We extend this data structure by allowing the presence of a possibly time-dependent covariate process that is observed until end of follow up. If one only assumes that the censoring mechanism satisfies coarsening at random, then, by the curve of dimensionality, typically no regular estimators will exist. To fight the curse of dimensionality we follow the approach of Robins and Rotnitzky (1992) by modeling parameters of the censoring mechanism. We model the right-censoring mechanism by modeling the hazard of the follow up time, conditional on T and the covariate process. For the monitoring mechanism we avoid modeling the joint distribution of the monitoring times by only modeling a univariate hazard of the pooled monitoring times, conditional on the follow up time, T, and the covariates process, which can be estimated by treating the pooled sample of monitoring times as i.i.d. In particular, it is assumed that the monitoring times and the right-censoring times only depend on T through the observed covariate process. We introduce inverse probability of censoring weighted (IPCW) estimator of the distribution of T and of smooth functionals thereof which are guaranteed to be consistent and asymptotically normal if we have available correctly specified semiparametric models for the two hazards of the censoring process. Furthermore, given such correctly specified models for these hazards of the censoring process, we propose a one-step estimator which will improve on the IPCW estimator if we correctly specify a lower-dimensional working model for the conditional distribution of T, given the covariate process, that remains consistent and asymptotically normal if this latter working model is misspecified. It is shown that the one-step estimator is efficient if each subject is at most monitored once and the working model contains the truth. In general, it is shown that the one-step estimator optimally uses the surrogate information if the working model contains the truth. It is not optimal in using the interval information provided by the current status indicators at the monitoring times, but simulations in Peterson, van der Laan (1997) show that the efficiency loss is small.

Smooth Hazard Function Estimation for Interval Censored Data with Time Varying Covariates

In this paper we propose methods for smooth hazard estimation of a time variable where that variable is interval censored. These methods allow one to model the transformed hazard in terms of either smooth (smoothing splines) or linear functions of time and other relevant time varying predictor variables. We illustrate the use of this method on a dataset of hemophiliacs where the outcome, time to seroconversion for HIV, is interval censored and left-truncated.

The Two-Interval Line-Segment Problem

In this paper, the NPMLE in the one-dimensional line segment problem is defined and studied, where line segments on the real line through two non-overlapping intervals are observed. The self-consistency equations for the NPMLE are defined and a quick algorithm for solving them is provided. Supnorm weak convergence to a Gaussian process and efficiency of the NPMLE is proved. The problem has a strong geological application in the study of the lifespan of species.