Lung cancer screening with CT: evaluation of radiologists and different computer assisted detection software (CAD) as first and second readers for lung nodule detection at different dose levels

OBJECTIVES To find the best pairing of first and second reader at highest sensitivity for detecting lung nodules with CT at various dose levels. MATERIALS AND METHODS An anthropomorphic lung phantom and artificial lung nodules were used to simulate screening CT-examination at standard dose (100 mAs, 120 kVp) and 8 different low dose levels, using 120, 100 and 80 kVp combined with 100, 50 and 25 mAs. At each dose level 40 phantoms were randomly filled with 75 solid and 25 ground glass nodules (5-12 mm). Two radiologists and 3 different computer aided detection softwares (CAD) were paired to find the highest sensitivity. RESULTS Sensitivities at standard dose were 92%, 90%, 84%, 79% and 73% for reader 1, 2, CAD1, CAD2, CAD3, respectively. Combined sensitivity for human readers 1 and 2 improved to 97%, (p1=0.063, p2=0.016). Highest sensitivities--between 97% and 99.0%--were achieved by combining any radiologist with any CAD at any dose level. Combining any two CADs, sensitivities between 85% and 88% were significantly lower than for radiologists combined with CAD (p<0.03). CONCLUSIONS Combination of a human observer with any of the tested CAD systems provide optimal sensitivity for lung nodule detection even at reduced dose at 25 mAs/80 kVp.

The applicability of using different energy levels in CT imaging for differentiation or identification of dental restorative materials.

PURPOSE The goal of this study was to investigate whether different computed tomography (CT) energy levels could supply additional information for the differentiation of dental materials for forensic investigations. METHODS Nine different commonly used restorative dental materials were investigated in this study. A total of 75 human third molars were filled with the restorative dental materials and then scanned using the forensic reference phantom in singlesource mode. The mean Hounsfield unit values and standard deviations (SDs) of each material were calculated at 120, 80 and 140 kVp. RESULTS Most of the dental materials could be differentiated at 120 kVp. We found that greater X-ray density of a material resulted in higher SDs and that the material volume could influence the measurements. CONCLUSION Differentiation of dental materials in CT was possible in many cases using single-energy CT scans at 120 kVp. Because of the number of dental restorative materials available and scanner and scan parameter dependence, as well as the CT imaging artifacts, the identification (in contrast to differentiation) was problematic.

Impact of different magnification levels on visual caries detection with ICDAS

OBJECTIVES The aim of this in vitro study was to examine the effect of different levels of magnification on the accuracy and reliability of visual caries detection using ICDAS criteria. METHODS Occlusal surfaces of 100 extracted molars were assessed by 14 examiners (3rd and the 4th year dental students and dentists) using no magnification aids, a 2.5× Galilean loupe, a 4.5× Keplerian loupe, or a surgical microscope with 10× magnification. The assessments were repeated on a different day. Sensitivity, specificity, AUC and reliabilities were calculated according to the gold standard of histology. RESULTS We found that with increasing magnification, the number of surfaces rated as "sound" (ICDAS code 0) decreased, while the number of surfaces with a localized enamel breakdown (ICDAS code 3) increased. While the sensitivities increased, the values of the specificities decreased to an unacceptably low level irrespective of the clinical experience of the examiners. CONCLUSIONS ICDAS seems to be optimized for natural vision up to 2.0× magnification and not for high magnifications. The use of powerful magnification in visual caries detection involves the risk of unnecessary and premature invasive treatment. CLINICAL SIGNIFICANCE This paper discusses when it does and does not make sense to use magnification devices for visual caries detection using ICDAS criteria. Strong magnifications should be refrained from for this purpose.

Operator variability using different polishing methods and surface geometry of a nanohybrid composite

This study evaluated the operator variability of different finishing and polishing techniques. After placing 120 composite restorations (Tetric EvoCeram) in plexiglassmolds, the surface of the specimens was roughened in a standardized manner. Twelve operators with different experience levels polished the specimens using the following finishing/polishing procedures: method 1 (40 ?m diamond [40D], 15 ?m diamond [15D], 42 ?m silicon carbide polisher [42S], 6 ?m silicon carbide polisher [6S] and Occlubrush [O]); method 2 (40D, 42S, 6S and O); method 3 (40D, 42S, 6S and PoGo); method 4 (40D, 42S and PoGo) and method 5 (40D, 42S and O). The mean surface roughness (Ra) was measured with a profilometer. Differences between the methods were analyzed with non-parametric ANOVA and pairwise Wilcoxon signed rank tests (?=0.05). All the restorations were qualitatively assessed using SEM. Methods 3 and 4 showed the best polishing results and method 5 demonstrated the poorest. Method 5 was also most dependent on the skills of the operator. Except for method 5, all of the tested procedures reached a clinically acceptable surface polish of Ra?0.2 ?m. Polishing procedures can be simplified without increasing variability between operators and without jeopardizing polishing results.

Quantification of atelectatic lung volumes in two different porcine models of ARDS

BACKGROUND: Cyclic recruitment during mechanical ventilation contributes to ventilator associated lung injury. Two different pathomechanisms in acute respiratory distress syndrome (ARDS) are currently discussed: alveolar collapse vs persistent flooding of small airways and alveoli. We compare two different ARDS animal models by computed tomography (CT) to describe different recruitment and derecruitment mechanisms at different airway pressures: (i) lavage-ARDS, favouring alveolar collapse by surfactant depletion; and (ii) oleic acid ARDS, favouring alveolar flooding by capillary leakage. METHODS: In 12 pigs [25 (1) kg], ARDS was randomly induced, either by saline lung lavage or oleic acid (OA) injection, and 3 animals served as controls. A respiratory breathhold manoeuvre without spontaneous breathing at different continuous positive airway pressure (CPAP) was applied in random order (CPAP levels of 5, 10, 15, 30, 35 and 50 cm H(2)O) and spiral-CT scans of the total lung were acquired at each CPAP level (slice thickness=1 mm). In each spiral-CT the volume of total lung parenchyma, tissue, gas, non-aerated, well-aerated, poorly aerated, and over-aerated lung was calculated. RESULTS: In both ARDS models non-aerated lung volume decreased significantly from CPAP 5 to CPAP 50 [oleic acid lung injury (OAI): 346.9 (80.1) to 96.4 (48.8) ml, P<0.001; lavage-ARDS: 245 17.6) to 42.7 (4.8) ml, P<0.001]. In lavage-ARDS poorly aerated lung volume decreased at higher CPAP levels [232 (45.2) at CPAP 10 to 84 (19.4) ml at CPAP 50, P<0.001] whereas in OAI poorly aerated lung volume did not vary at different airway pressures. CONCLUSIONS: In both ARDS models well-aerated and non-aerated lung volume respond to different CPAP levels in a comparable fashion: Thus, a cyclical alveolar collapse seems to be part of the derecruitment process also in the OA-ARDS. In OA-ARDS, the increase in poorly aerated lung volume reflects the specific initial lesion, that is capillary leakage with interstitial and alveolar oedema.

Stochastic rank correlation: a robust merit function for 2D/3D registration of image data obtained at different energies

In this article, the authors evaluate a merit function for 2D/3D registration called stochastic rank correlation (SRC). SRC is characterized by the fact that differences in image intensity do not influence the registration result; it therefore combines the numerical advantages of cross correlation (CC)-type merit functions with the flexibility of mutual-information-type merit functions. The basic idea is that registration is achieved on a random subset of the image, which allows for an efficient computation of Spearman's rank correlation coefficient. This measure is, by nature, invariant to monotonic intensity transforms in the images under comparison, which renders it an ideal solution for intramodal images acquired at different energy levels as encountered in intrafractional kV imaging in image-guided radiotherapy. Initial evaluation was undertaken using a 2D/3D registration reference image dataset of a cadaver spine. Even with no radiometric calibration, SRC shows a significant improvement in robustness and stability compared to CC. Pattern intensity, another merit function that was evaluated for comparison, gave rather poor results due to its limited convergence range. The time required for SRC with 5% image content compares well to the other merit functions; increasing the image content does not significantly influence the algorithm accuracy. The authors conclude that SRC is a promising measure for 2D/3D registration in IGRT and image-guided therapy in general.

Threshold-dependent effects on peripheral nerve in vivo excitability properties in the rat

Various factors, including maturity, have been shown to influence peripheral nerve excitability measures, but little is known about differences in these properties between axons with different stimulation thresholds. Multiple nerve excitability tests were performed on the caudal motor axons of immature and mature female rats, recording from tail muscles at three target compound muscle action potential (CMAP) levels: 10%, 40% ("standard" level), and 60% of the maximum CMAP amplitude. Compared to lower target levels, axons at high target levels have the following characteristics: lower strength-duration time constant, less threshold reduction during depolarizing currents and greater threshold increase to hyperpolarizing currents, most notably to long hyperpolarizing currents in mature rats. Threshold-dependent effects on peripheral nerve excitability properties depend on the maturation stage, especially inward rectification (Ih), which becomes inversely related to threshold level. Performing nerve excitability tests at different target levels is useful in understanding the variation in membrane properties between different axons within a nerve. Because of the threshold effects on nerve excitability and the possibility of increased variability between axons and altered electric recruitment order in disease conditions, excitability parameters measured only at the "standard" target level should be interpreted with caution, especially the responses to hyperpolarizing currents.

Coupling biomechanics to a cellular level model: An approach to patient-specific image driven multi-scale and multi-physics tumor simulation

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20\% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

ABCtoolbox: a versatile toolkit for approximate Bayesian computations

Background The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. Results Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. Conclusion ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

ARIS-Campaign: Intercomparison of three ground based 22 GHz radiometers for middle atmospheric water vapor at the Zugspitze in winter 2009

This paper presents the Alpine Radiometer Intercomparison at the Schneefernerhaus (ARIS), which took place in winter 2009 at the high altitude station at the Zugspitze, Germany (47.42° N, 10.98° E, 2650 m). This campaign was the first direct intercomparison between three new ground based 22 GHz water vapor radiometers for middle atmospheric profiling with the following instruments participating: MIRA 5 (Karlsruhe Institute of Technology), cWASPAM3 (Max Planck Institute for Solar System Research, Katlenburg-Lindau) and MIAWARA-C (Institute of Applied Physics, University of Bern). Even though the three radiometers all measure middle atmospheric water vapor using the same rotational transition line and similar fundamental set-ups, there are major differences between the front ends, the back ends, the calibration concepts and the profile retrieval. The spectrum comparison shows that all three radiometers measure spectra without severe baseline artifacts and that the measurements are in good general agreement. The measurement noise shows good agreement to the values theoretically expected from the radiometer noise formula. At the same time the comparison of the noise levels shows that there is room for instrumental and calibration improvement, emphasizing the importance of low elevation angles for the observation, a low receiver noise temperature and an efficient calibration scheme. The comparisons of the retrieved profiles show that the agreement between the profiles of MIAWARA-C and cWASPAM3 with the ones of MLS is better than 0.3 ppmv (6%) at all altitudes. MIRA 5 has a dry bias of approximately 0.5 ppm (8%) below 0.1 hPa with respect to all other instruments. The profiles of cWASPAM3 and MIAWARA-C could not be directly compared because the vertical region of overlap was too small. The comparison of the time series at different altitude levels show a similar evolution of the H2O volume mixing ratio (VMR) for the ground based instruments as well as the space borne sensor MLS.

Effects of nutrient restriction on mammary cell turnover and mammary gland remodeling in lactating dairy cows

The aim of this study was to investigate the effects of a severe nutrient restriction on mammary tissue morphology and remodeling, mammary epithelial cell (MEC) turnover and activity, and hormonal status in lactating dairy cows. We used 16 Holstein x Normande crossbred dairy cows, divided into 2 groups submitted to different feeding levels (basal and restricted) from 2 wk before calving to wk 11 postpartum. Restricted-diet cows had lower 11-wk average daily milk yield from calving to slaughter than did basal-diet cows (20.5 vs. 33.5 kg/d). Feed restriction decreased milk fat, protein, and lactose yields. Restriction also led to lower plasma insulin-like growth factor 1 and higher growth hormone concentrations. Restricted-diet cows had lighter mammary glands than did basal-diet cows. The total amount of DNA in the mammary gland and the size of the mammary acini were smaller in the restricted-diet group. Feed restriction had no significant effect on MEC proliferation at the time of slaughter but led to a higher level of apoptosis in the mammary gland. Gelatin zymography highlighted remodeling of the mammary extracellular matrix in restricted-diet cows. Udders from restricted-diet cows showed lower transcript expression of alpha-lactalbumin and kappa-casein. In conclusion, nutrient restriction resulted in lower milk yield in lactating dairy cows, partly due to modulation of MEC activity and a lower number of mammary cells. An association was found between feed restriction-induced changes in the growth hormone-insulin-like growth factor-1 axis and mammary epithelial cell dynamics.

High interstitial fluid pressure is associated with low tumour penetration of diagnostic monoclonal antibodies applied for molecular imaging purposes

The human epithelial cell adhesion molecule (EpCAM) is highly expressed in a variety of clinical tumour entities. Although an antibody against EpCAM has successfully been used as an adjuvant therapy in colon cancer, this therapy has never gained wide-spread use. We have therefore investigated the possibilities and limitations for EpCAM as possible molecular imaging target using a panel of preclinical cancer models. Twelve human cancer cell lines representing six tumour entities were tested for their EpCAM expression by qPCR, flow cytometry analysis and immunocytochemistry. In addition, EpCAM expression was analyzed in vivo in xenograft models for tumours derived from these cells. Except for melanoma, all cell lines expressed EpCAM mRNA and protein when grown in vitro. Although they exhibited different mRNA levels, all cell lines showed similar EpCAM protein levels upon detection with monoclonal antibodies. When grown in vivo, the EpCAM expression was unaffected compared to in vitro except for the pancreatic carcinoma cell line 5072 which lost its EpCAM expression in vivo. Intravenously applied radio-labelled anti EpCAM MOC31 antibody was enriched in HT29 primary tumour xenografts indicating that EpCAM binding sites are accessible in vivo. However, bound antibody could only be immunohistochemically detected in the vicinity of perfused blood vessels. Investigation of the fine structure of the HT29 tumour blood vessels showed that they were immature and prone for higher fluid flux into the interstitial space. Consistent with this hypothesis, a higher interstitial fluid pressure of about 12 mbar was measured in the HT29 primary tumour via "wick-in-needle" technique which could explain the limited diffusion of the antibody into the tumour observed by immunohistochemistry.

A global ”imaging” view on systems approaches in immunology

The immune system exhibits an enormous complexity. High throughput methods such as the "-omic'' technologies generate vast amounts of data that facilitate dissection of immunological processes at ever finer resolution. Using high-resolution data-driven systems analysis, causal relationships between complex molecular processes and particular immunological phenotypes can be constructed. However, processes in tissues, organs, and the organism itself (so-called higher level processes) also control and regulate the molecular (lower level) processes. Reverse systems engineering approaches, which focus on the examination of the structure, dynamics and control of the immune system, can help to understand the construction principles of the immune system. Such integrative mechanistic models can properly describe, explain, and predict the behavior of the immune system in health and disease by combining both higher and lower level processes. Moving from molecular and cellular levels to a multiscale systems understanding requires the development of methodologies that integrate data from different biological levels into multiscale mechanistic models. In particular, 3D imaging techniques and 4D modeling of the spatiotemporal dynamics of immune processes within lymphoid tissues are central for such integrative approaches. Both dynamic and global organ imaging technologies will be instrumental in facilitating comprehensive multiscale systems immunology analyses as discussed in this review.

Speciation leads to divergent methylmercury accumulation in sympatric whitefish

Central European lake whitefish (Coregonus spp.) colonized Swiss lakes following the last glacial retreat and have undergone rapid speciation and adaptive radiation. Up to six species have been shown to coexist in some lakes, and individual species occupy specific ecological niches and have distinct feeding and reproductive ecologies. We studied methylmercury (MeHg) accumulation in sympatric whitefish species from seven Swiss lakes to determine if ecological divergence has led to different rates of MeHg bioaccumulation. In four of seven lakes, sympatric species had distinctly different MeHg levels, which varied by up to a factor of two between species. Generally, species with greater MeHg levels were smaller in body size and planktivorous, and species with lower MeHg were larger and benthivorous. While modest disparities in trophic position between species might be expected a priori to explain the divergence in MeHg, δ15N of bulk tissue did not correlate with fish MeHg in five of seven lakes. Results of a nested ANCOVA analysis across all lakes indicated that only two factors (species, lake) explained substantial portions of the variance, with species accounting for more variance (52 %) than inter-lake differences (32 %). We suggest that differences in MeHg accumulation were likely caused by diverging metabolic traits between species, such as differences in energy partitioning between anabolism and catabolism, potentially interacting with species-specific prey resource utilization. These results indicate substantial variability in MeHg accumulation between closely related fish species, illustrating that ecological speciation in fish can lead to divergent MeHg accumulation patterns.

Dissociated lateralization of transient and sustained blood oxygen level-dependent signal components in human primary auditory cortex

Among other auditory operations, the analysis of different sound levels received at both ears is fundamental for the localization of a sound source. These so-called interaural level differences, in animals, are coded by excitatory-inhibitory neurons yielding asymmetric hemispheric activity patterns with acoustic stimuli having maximal interaural level differences. In human auditory cortex, the temporal blood oxygen level-dependent (BOLD) response to auditory inputs, as measured by functional magnetic resonance imaging (fMRI), consists of at least two independent components: an initial transient and a subsequent sustained signal, which, on a different time scale, are consistent with electrophysiological human and animal response patterns. However, their specific functional role remains unclear. Animal studies suggest these temporal components being based on different neural networks and having specific roles in representing the external acoustic environment. Here we hypothesized that the transient and sustained response constituents are differentially involved in coding interaural level differences and therefore play different roles in spatial information processing. Healthy subjects underwent monaural and binaural acoustic stimulation and BOLD responses were measured using high signal-to-noise-ratio fMRI. In the anatomically segmented Heschl's gyrus the transient response was bilaterally balanced, independent of the side of stimulation, while in opposite the sustained response was contralateralized. This dissociation suggests a differential role at these two independent temporal response components, with an initial bilateral transient signal subserving rapid sound detection and a subsequent lateralized sustained signal subserving detailed sound characterization.