The Post-Illumination Pupil Response (PIPR)

Purpose The post-illumination pupil response (PIPR) has been quantified using four metrics, but the spectral sensitivity of only one is known; here we determine the other three. To optimize the human PIPR measurement, we determine the protocol producing the largest PIPR, the duration of the PIPR, and the metric(s) with the lowest coefficient of variation. Methods The consensual pupil light reflex (PLR) was measured with a Maxwellian view pupillometer. - Experiment 1: Spectral sensitivity of four PIPR metrics [plateau, 6 s, area under curve (AUC) early and late recovery] was determined from a criterion PIPR to a 1s pulse and fitted with Vitamin A1 nomogram (λmax = 482nm). - Experiment 2: The PLR was measured as a function of three stimulus durations (1s, 10s, 30s), five irradiances spanning low to high melanopsin excitation levels (retinal irradiance: 9.8 to 14.8 log quanta.cm-2.s-1), and two wavelengths, one with high (465nm) and one with low (637nm) melanopsin excitation. Intra and inter-individual coefficients of variation (CV) were calculated. Results The melanopsin (opn4) photopigment nomogram adequately describes the spectral sensitivity of all four PIPR metrics. The PIPR amplitude was largest with 1s short wavelength pulses (≥ 12.8 log quanta.cm-2.s-1). The plateau and 6s PIPR showed the least intra and inter-individual CV (≤ 0.2). The maximum duration of the sustained PIPR was 83.0±48.0s (mean±SD) for 1s pulses and 180.1±106.2s for 30s pulses (465nm; 14.8 log quanta.cm-2.s-1). Conclusions All current PIPR metrics provide a direct measure of the intrinsic melanopsin photoresponse. To measure progressive changes in melanopsin function in disease, we recommend that the PIPR be measured using short duration pulses (e.g., ≤ 1s) with high melanopsin excitation and analyzed with plateau and/or 6s metrics. Our PIPR duration data provide a baseline for the selection of inter-stimulus intervals between consecutive pupil testing sequences.

The post-illumination pupil response of melanopsin expressing retinal ganglion cells

Purpose The post-illumination pupil response (PIPR) has been quantified in the literature by four metrics. The spectral sensitivity of only one metric is known and this study quantifies the other three. To optimize the measurement of the PIPR in humans, we also determine the stimulus protocol producing the largest PIPR, the duration of the PIPR, and the metric(s) with the lowest coefficient of variation. Methods The consensual pupil light reflex (PLR) was measured with a Maxwellian view pupillometer (35.6° diameter stimulus). - Experiment 1: Spectral sensitivity of four PIPR metrics [plateau, 6 s, area under curve (AUC) early and late recovery] was determined from a criterion PIPR (n = 2 participants) to a 1 s pulse at five wavelengths (409-592nm) and fitted with Vitamin A nomogram (ƛmax = 482 nm). - Experiment 2: The PLR was measured in five healthy participants [29 to 42 years (mean = 32.6 years)] as a function of three stimulus durations (1 s, 10 s, 30 s), five irradiances spanning low to high melanopsin excitation levels (retinal irradiance: 9.8 to 14.8 log quanta.cm-2.s-1), and two wavelengths, one with high (465 nm) and one with low (637 nm) melanopsin excitation. Intra and inter-individual coefficients of variation (CV) were calculated. Results The melanopsin (opn4) photopigment nomogram adequately described the spectral sensitivity derived from all four PIPR metrics. The largest PIPR amplitude was observed with 1 s short wavelength pulses (retinal irradiance ≥ 12.8 log quanta.cm-2.s-1). Of the 4 PIPR metrics, the plateau and 6 s PIPR showed the least intra and inter-individual CV (≤ 0.2). The maximum duration of the sustained PIPR was 83.4 ± 48.0 s (mean ± SD) for 1 s pulses and 180.1 ± 106.2 s for 30 s pulses (465 nm; 14.8 log quanta.cm-2.s-1). Conclusions All current PIPR metrics provide a direct measure of intrinsic melanopsin retinal ganglion cell function. To measure progressive changes in melanopsin function in disease, we recommend that the intrinsic melanopsin response should be measured using a 1 s pulse with high melanopsin excitation and the PIPR should be analyzed with the plateau and/or 6 s metrics. That the PIPR can have a sustained constriction for as long as 3 minutes, our PIPR duration data provide a baseline for the selection of inter-stimulus intervals between consecutive pupil testing sequences.

Arsenic trioxide promotes mitochondrial DNA mutation and cell apoptosis in primary APL cells and NB4 cell line

This study aimed to investigate the effects of arsenic trioxide (As2O3) on the mitochondrial DNA (mtDNA) of acute promyelocytic leukemia (APL) cells. The NB4 cell line was treated with 2.0 μmol/L As2O3in vitro, and the primary APL cells were treated with 2.0 μmol/L As2O3in vitro and 0.16 mg kg-1 d-1 As2O3in vivo. The mitochondrial DNA of all the cells above was amplified by PCR, directly sequenced and analyzed by Sequence Navigatore and Factura software. The apoptosis rates were assayed by flow cytometry. Mitochondrial DNA mutation in the D-loop region was found in NB4 and APL cells before As2O3 use, but the mutation spots were remarkably increased after As2O3 treatment, which was positively correlated to the rates of cellular apoptosis, the correlation coefficient: rNB4-As2O3=0.973818, and rAPL-As2O3=0.934703. The mutation types include transition, transversion, codon insertion or deletion, and the mutation spots in all samples were not constant and regular. It is revealed that As2O3 aggravates mtDNA mutation in the D-loop region of acute promyelocytic leukemia cells both in vitro and in vivo. Mitochondrial DNA might be one of the targets of As2O3 in APL treatment.

Correlation of carotid atheromatous plaque inflammation with biomechanical stress: Utility of USPIO enhanced MR imaging and finite element analysis

Objective: The aim of this study was to explore whether there is a relationship between the degree of MR-defined inflammation using ultra small super-paramagnetic iron oxide (USPIO) particles, and biomechanical stress using finite element analysis (FEA) techniques, in carotid atheromatous plaques. Methods and Results: 18 patients with angiographically proven carotid stenoses underwent multi-sequence MR imaging before and 36 h after USPIO infusion. T2 * weighted images were manually segmented into quadrants and the signal change in each quadrant normalised to adjacent muscle was calculated after USPIO administration. Plaque geometry was obtained from the rest of the multi-sequence dataset and used within a FEA model to predict maximal stress concentration within each slice. Subsequently, a new statistical model was developed to explicitly investigate the form of the relationship between biomechanical stress and signal change. The Spearman's rank correlation coefficient for USPIO enhanced signal change and maximal biomechanical stress was -0.60 (p = 0.009). Conclusions: There is an association between biomechanical stress and USPIO enhanced MR-defined inflammation within carotid atheroma, both known risk factors for plaque vulnerability. This underlines the complex interaction between physiological processes and biomechanical mechanisms in the development of carotid atheroma. However, this is preliminary data that will need validation in a larger cohort of patients.

Assessment of carotid plaque vulnerability using structural and geometrical determinants

Background Because many acute cerebral ischemic events are caused by rupture of vulnerable carotid atheroma and subsequent thrombosis, the present study used both idealized and patient-specific carotid atheromatous plaque models to evaluate the effect of structural determinants on stress distributions within plaque. Methods and Results Using a finite element method, structural analysis was performed using models derived from in vivo high-resolution magnetic resonance imaging (MRI) of carotid atheroma in 40 non-consecutive patients (20 symptomatic, 20 asymptomatic). Plaque components were modeled as hyper-elastic materials. The effects of varying fibrous cap thickness, lipid core size and lumen curvature on plaque stress distributions were examined. Lumen curvature and fibrous cap thickness were found to be major determinants of plaque stress. The size of the lipid core did not alter plaque stress significantly when the fibrous cap was relatively thick. The correlation between plaque stress and lumen curvature was significant for both symptomatic (p = 0.01; correlation coefficient: 0.689) and asymptomatic patients (p = 0.01; correlation coefficient: 0.862). Lumen curvature in plaques of symptomatic patients was significantly larger than those of asymptomatic patients (1.50±1.0mm-1 vs 1.25±0.75 mm-1; p = 0.01). Conclusion Specific plaque morphology (large lumen curvature and thin fibrous cap) is closely related to plaque vulnerability. Structural analysis using high-resolution MRI of carotid atheroma may help in detecting vulnerable atheromatous plaque and aid the risk stratification of patients with carotid disease.

Measurement of stenotic carotid arterial compliance with MRI

Arterial compliance has been shown to correlate well with overall cardiovascular outcome and it may also be a potential risk factor for the development of atheromatous disease. This study assesses the utility of 2-D phase contrast Magnetic Resonance (MR) imaging with intra-sequence blood pressure measurement to determine carotid compliance and distensibility. 20 patients underwent 2-D phase contrast MR imaging and also ultrasound-based wall tracking measurements. Values for carotid compliance and distensibility were derived from the two different modalities and compared. Linear regression analysis was utilised to determine the extent of correlation between MR and ultrasound derived parameters. In those variables that could be directly compared, an agreement analysis was undertaken. MR measures of compliance showed a good correlation with measures based on ultrasound wall-tracking (r=0.61, 95% CI 0.34 to 0.81 p=0.0003). Vessels that had undergone carotid endarterectomy previously were significantly less compliant than either diseased or normal contralateral vessels (p=0.04). Agreement studies showed a relatively poor intra-class correlation coefficient (ICC) between diameter-based measures of compliance through either MR or ultrasound (ICC=0.14). MRI based assessment of local carotid compliance appears to be both robust and technically feasible in most subjects. Measures of compliance correlate well with ultrasound-based values and correlate best when cross-sectional area change is used rather than derived diameter changes. If validated by further larger studies, 2-D phase contrast imaging with intra-sequence blood pressure monitoring and off-line radial artery tonometry may provide a useful tool in further assessment of patients with carotid atheroma.

Sediment concentration prediction and statistical evaluation for annual load estimation

We consider the development of statistical models for prediction of constituent concentration of riverine pollutants, which is a key step in load estimation from frequent flow rate data and less frequently collected concentration data. We consider how to capture the impacts of past flow patterns via the average discounted flow (ADF) which discounts the past flux based on the time lapsed - more recent fluxes are given more weight. However, the effectiveness of ADF depends critically on the choice of the discount factor which reflects the unknown environmental cumulating process of the concentration compounds. We propose to choose the discount factor by maximizing the adjusted R-2 values or the Nash-Sutcliffe model efficiency coefficient. The R2 values are also adjusted to take account of the number of parameters in the model fit. The resulting optimal discount factor can be interpreted as a measure of constituent exhaustion rate during flood events. To evaluate the performance of the proposed regression estimators, we examine two different sampling scenarios by resampling fortnightly and opportunistically from two real daily datasets, which come from two United States Geological Survey (USGS) gaging stations located in Des Plaines River and Illinois River basin. The generalized rating-curve approach produces biased estimates of the total sediment loads by -30% to 83%, whereas the new approaches produce relatively much lower biases, ranging from -24% to 35%. This substantial improvement in the estimates of the total load is due to the fact that predictability of concentration is greatly improved by the additional predictors.

Effects of fish density distribution and effort distribution on catchability

The effects of fish density distribution and effort distribution on the overall catchability coefficient are examined. Emphasis is also on how aggregation and effort distribution interact to affect overall catch rate [catch per unit effort (cpue)]. In particular, it is proposed to evaluate three indices, the catchability index, the knowledge parameter, and the aggregation index, to describe the effectiveness of targeting and the effects on overall catchability in the stock area. Analytical expressions are provided so that these indices can easily be calculated. The average of the cpue calculated from small units where fishing is random is a better index for measuring the stock abundance. The overall cpue, the ratio of lumped catch and effort, together with the average cpue, can be used to assess the effectiveness of targeting. The proposed methods are applied to the commercial catch and effort data from the Australian northern prawn fishery. The indices are obtained assuming a power law for the effort distribution as an approximation of targeting during the fishing operation. Targeting increased catchability in some areas by 10%, which may have important implications on management advice.

Achilles tendon loading patterns during barefoot walking and slow running on a treadmill: An ultrasonic propagation study

Measurement of tendon loading patterns during gait is important for understanding the pathogenesis of tendon "overuse" injury. Given that the speed of propagation of ultrasound in tendon is proportional to the applied load, this study used a noninvasive ultrasonic transmission technique to measure axial ultrasonic velocity in the right Achilles tendon of 27 healthy adults (11 females and 16 males; age, 26 ± 9 years; height, 1.73 ± 0.07 m; weight, 70.6 ± 21.2 kg), walking at self-selected speed (1.1 ± 0.1 m/s), and running at fixed slow speed (2 m/s) on a treadmill. Synchronous measures of ankle kinematics, spatiotemporal gait parameters, and vertical ground reaction forces were simultaneously measured. Slow running was associated with significantly higher cadence, shorter step length, but greater range of ankle movement, higher magnitude and rate of vertical ground reaction force, and higher ultrasonic velocity in the tendon than walking (P < 0.05). Ultrasonic velocity in the Achilles tendon was highly reproducible during walking and slow running (mean within-subject coefficient of variation < 2%). Ultrasonic maxima (P1, P2) and minima (M1, M2) were significantly higher and occurred earlier in the gait cycle (P1, M1, and M2) during running than walking (P < 0.05). Slow running was associated with higher and earlier peaks in loading of the Achilles tendon than walking.

Morphology-based interslice interpolation on manual segmentations of joint bones and muscles in MRI

This paper presents a validation study on the application of a novel interslice interpolation technique for musculoskeletal structure segmentation of articulated joints and muscles on human magnetic resonance imaging data. The interpolation technique is based on morphological shape-based interpolation combined with intensity based voxel classification. Shape-based interpolation in the absence of the original intensity image has been investigated intensively. However, in some applications of medical image analysis, the intensity image of the slice to be interpolated is available. For example, when manual segmentation is conducted on selected slices, the segmentation on those unselected slices can be obtained by interpolation. We proposed a two- step interpolation method to utilize both the shape information in the manual segmentation and local intensity information in the image. The method was tested on segmentations of knee, hip and shoulder joint bones and hamstring muscles. The results were compared with two existing interpolation methods. Based on the calculated Dice similarity coefficient and normalized error rate, the proposed method outperformed the other two methods.

Atmospheric visibility and PM10 as indicators of new particle formation in an urban environment

It is well-known that new particle formation (NPF) in the atmosphere is inhibited by pre-existing particles in the air that act as condensation sinks to decrease the concentration and, thus, the supersaturation of precursor gases. In this study, we investigate the effects of two parameters - atmospheric visibility, expressed as the particle back-scatter coefficient (BSP), and PM10 particulate mass concentration, on the occurrences of NPF events in an urban environment where the majority of precursor gases originate from motor vehicle and industrial sources. This is the first attempt to derive direct relationships between each of these two parameters and the occurrence of NPF. NPF events were identified from data obtained with a neutral cluster and air ion spectrometer over 245 days within a calendar year. Bayesian logistic regression was used to determine the probability of observing NPF as functions of BSP and PM10. We show that the BSP at 08 h on a given day is a reliable indicator of an NPF event later that day. The posterior median probability of observing an NPF event was greater than 0.5 (95%) when the BSP at 08 h was less than 6.8 Mm-1.

Calculating the Fickian diffusivity for a lattice-based random walk with agents and obstacles of different shapes and sizes

Random walk models are often used to interpret experimental observations of the motion of biological cells and molecules. A key aim in applying a random walk model to mimic an in vitro experiment is to estimate the Fickian diffusivity (or Fickian diffusion coefficient),D. However, many in vivo experiments are complicated by the fact that the motion of cells and molecules is hindered by the presence of obstacles. Crowded transport processes have been modeled using repeated stochastic simulations in which a motile agent undergoes a random walk on a lattice that is populated by immobile obstacles. Early studies considered the most straightforward case in which the motile agent and the obstacles are the same size. More recent studies considered stochastic random walk simulations describing the motion of an agent through an environment populated by obstacles of different shapes and sizes. Here, we build on previous simulation studies by analyzing a general class of lattice-based random walk models with agents and obstacles of various shapes and sizes. Our analysis provides exact calculations of the Fickian diffusivity, allowing us to draw conclusions about the role of the size, shape and density of the obstacles, as well as examining the role of the size and shape of the motile agent. Since our analysis is exact, we calculateDdirectly without the need for random walk simulations. In summary, we find that the shape, size and density of obstacles has a major influence on the exact Fickian diffusivity. Furthermore, our results indicate that the difference in diffusivity for symmetric and asymmetric obstacles is significant.

Survival probability for a diffusive process on a growing domain

We consider the motion of a diffusive population on a growing domain, 0 < x < L(t ), which is motivated by various applications in developmental biology. Individuals in the diffusing population, which could represent molecules or cells in a developmental scenario, undergo two different kinds of motion: (i) undirected movement, characterized by a diffusion coefficient, D, and (ii) directed movement, associated with the underlying domain growth. For a general class of problems with a reflecting boundary at x = 0, and an absorbing boundary at x = L(t ), we provide an exact solution to the partial differential equation describing the evolution of the population density function, C(x,t ). Using this solution, we derive an exact expression for the survival probability, S(t ), and an accurate approximation for the long-time limit, S = limt→∞ S(t ). Unlike traditional analyses on a nongrowing domain, where S ≡ 0, we show that domain growth leads to a very different situation where S can be positive. The theoretical tools developed and validated in this study allow us to distinguish between situations where the diffusive population reaches the moving boundary at x = L(t ) from other situations where the diffusive population never reaches the moving boundary at x = L(t ). Making this distinction is relevant to certain applications in developmental biology, such as the development of the enteric nervous system (ENS). All theoretical predictions are verified by implementing a discrete stochastic model.

Magnetic resonance microimaging of cancer cell spheroid constructs

BACKGROUND Hydrogel-based cell cultures are excellent tools for studying physiological events occurring in the growth and proliferation of cells, including cancer cells. Diffusion magnetic resonance is a physical technique that has been widely used for the characterisation of biological systems as well as hydrogels. In this work, we applied diffusion magnetic resonance imaging (MRI) to hydrogel-based cultures of human ovarian cancer cells. METHODS Diffusion-weighted spin-echo MRI measurements were used to obtain spatially-resolved maps of apparent diffusivities for hydrogel samples with different compositions, cell loads and drug (Taxol) treatment regimes. The samples were then characterised using their diffusivity histograms, mean diffusivities and the respective standard deviations, and pairwise Mann-Whitney tests. The elastic moduli of the samples were determined using mechanical compression testing. RESULTS The mean apparent diffusivity of the hydrogels was sensitive to the polymer content, cell load and Taxol treatment. For a given sample composition, the mean apparent diffusivity and the elastic modulus of the hydrogels exhibited a negative correlation. CONCLUSIONS Diffusivity of hydrogel-based cancer cell culture constructs is sensitive to both cell proliferation and Taxol treatment. This suggests that diffusion-weighted imaging is a promising technique for non-invasive monitoring of cancer cell proliferation in hydrogel-based, cellularly-sparse 3D cell cultures. The negative correlation between mean apparent diffusivity and elastic modulus suggests that the diffusion coefficient is indicative of the average density of the physical microenvironment within the hydrogel construct.

GWAS of 126,559 individuals identifies genetic variants associated with educational attainment

A genome-wide association study (GWAS) of educational attainment was conducted in a discovery sample of 101,069 individuals and a replication sample of 25,490. Three independent single-nucleotide polymorphisms (SNPs) are genome-wide significant (rs9320913, rs11584700, rs4851266), and all three replicate. Estimated effects sizes are small (coefficient of determination R(2) approximately 0.02%), approximately 1 month of schooling per allele. A linear polygenic score from all measured SNPs accounts for approximately 2% of the variance in both educational attainment and cognitive function. Genes in the region of the loci have previously been associated with health, cognitive, and central nervous system phenotypes, and bioinformatics analyses suggest the involvement of the anterior caudate nucleus. These findings provide promising candidate SNPs for follow-up work, and our effect size estimates can anchor power analyses in social-science genetics.