Combining ganglion cell topology and data of patients with glaucoma to determine a structure-function map

PURPOSE: To introduce techniques for deriving a map that relates visual field locations to optic nerve head (ONH) sectors and to use the techniques to derive a map relating Medmont perimetric data to data from the Heidelberg Retinal Tomograph. METHODS: Spearman correlation coefficients were calculated relating each visual field location (Medmont M700) to rim area and volume measures for 10 degrees ONH sectors (HRT III software) for 57 participants: 34 with glaucoma, 18 with suspected glaucoma, and 5 with ocular hypertension. Correlations were constrained to be anatomically plausible with a computational model of the axon growth of retinal ganglion cells (Algorithm GROW). GROW generated a map relating field locations to sectors of the ONH. The sector with the maximum statistically significant (P < 0.05) correlation coefficient within 40 degrees of the angle predicted by GROW for each location was computed. Before correlation, both functional and structural data were normalized by either normative data or the fellow eye in each participant. RESULTS: The model of axon growth produced a 24-2 map that is qualitatively similar to existing maps derived from empiric data. When GROW was used in conjunction with normative data, 31% of field locations exhibited a statistically significant relationship. This significance increased to 67% (z-test, z = 4.84; P < 0.001) when both field and rim area data were normalized with the fellow eye. CONCLUSIONS: A computational model of axon growth and normalizing data by the fellow eye can assist in constructing an anatomically plausible map connecting visual field data and sectoral ONH data.

Intrinsically photosensitive melanopsin retinal ganglion cell contributions to the post-illumination pupil response and circadian rhythm

Intrinsically photosensitive retinal ganglion cells (ipRGCs) in the eye transmit the environmental light level, projecting to the suprachiasmatic nucleus (SCN) (Berson, Dunn & Takao, 2002; Hattar, Liao, Takao, Berson & Yau, 2002), the location of the circadian biological clock, and the olivary pretectal nucleus (OPN) of the pretectum, the start of the pupil reflex pathway (Hattar, Liao, Takao, Berson & Yau, 2002; Dacey, Liao, Peterson, Robinson, Smith, Pokorny, Yau & Gamlin, 2005). The SCN synchronizes the circadian rhythm, a cycle of biological processes coordinated to the solar day, and drives the sleep/wake cycle by controlling the release of melatonin from the pineal gland (Claustrat, Brun & Chazot, 2005). Encoded photic input from ipRGCs to the OPN also contributes to the pupil light reflex (PLR), the constriction and recovery of the pupil in response to light. IpRGCs control the post-illumination component of the PLR, the partial pupil constriction maintained for > 30 sec after a stimulus offset (Gamlin, McDougal, Pokorny, Smith, Yau & Dacey, 2007; Kankipati, Girkin & Gamlin, 2010; Markwell, Feigl & Zele, 2010). It is unknown if intrinsic ipRGC and cone-mediated inputs to ipRGCs show circadian variation in their photon-counting activity under constant illumination. If ipRGCs demonstrate circadian variation of the pupil response under constant illumination in vivo, when in vitro ipRGC activity does not (Weng, Wong & Berson, 2009), this would support central control of the ipRGC circadian activity. A preliminary experiment was conducted to determine the spectral sensitivity of the ipRGC post-illumination pupil response under the experimental conditions, confirming the successful isolation of the ipRGC response (Gamlin, et al., 2007) for the circadian experiment. In this main experiment, we demonstrate that ipRGC photon-counting activity has a circadian rhythm under constant experimental conditions, while direct rod and cone contributions to the PLR do not. Intrinsic ipRGC contributions to the post-illumination pupil response decreased 2:46 h prior to melatonin onset for our group model, with the peak ipRGC attenuation occurring 1:25 h after melatonin onset. Our results suggest a centrally controlled evening decrease in ipRGC activity, independent of environmental light, which is temporally synchronized (demonstrates a temporal phase-advanced relationship) to the SCN mediated release of melatonin. In the future the ipRGC post-illumination pupil response could be developed as a fast, non-invasive measure of circadian rhythm. This study establishes a basis for future investigation of cortical feedback mechanisms that modulate ipRGC activity.

Phenomenological modeling of cell-to-cell and beat-to-beat variability in isolated Guinea Pig ventricular myocytes

Experimental action potential (AP) recordings in isolated ventricular myoctes display significant temporal beat-to-beat variability in morphology and duration. Furthermore, significant cell-to-cell differences in AP also exist even for isolated cells originating from the same region of the same heart. However, current mathematical models of ventricular AP fail to replicate the temporal and cell-to-cell variability in AP observed experimentally. In this study, we propose a novel mathematical framework for the development of phenomenological AP models capable of capturing cell-to-cell and temporal variabilty in cardiac APs. A novel stochastic phenomenological model of the AP is developed, based on the deterministic Bueno-Orovio/Fentonmodel. Experimental recordings of AP are fit to the model to produce AP models of individual cells from the apex and the base of the guinea-pig ventricles. Our results show that the phenomenological model is able to capture the considerable differences in AP recorded from isolated cells originating from the location. We demonstrate the closeness of fit to the available experimental data which may be achieved using a phenomenological model, and also demonstrate the ability of the stochastic form of the model to capture the observed beat-to-beat variablity in action potential duration.

Sensitivity of edge detection methods for quantifying cell migration assays

Quantitative imaging methods to analyze cell migration assays are not standardized. Here we present a suite of two–dimensional barrier assays describing the collective spreading of an initially–confined population of 3T3 fibroblast cells. To quantify the motility rate we apply two different automatic image detection methods to locate the position of the leading edge of the spreading population after 24, 48 and 72 hours. These results are compared with a manual edge detection method where we systematically vary the detection threshold. Our results indicate that the observed spreading rates are very sensitive to the choice of image analysis tools and we show that a standard measure of cell migration can vary by as much as 25% for the same experimental images depending on the details of the image analysis tools. Our results imply that it is very difficult, if not impossible, to meaningfully compare previously published measures of cell migration since previous results have been obtained using different image analysis techniques and the details of these techniques are not always reported. Using a mathematical model, we provide a physical interpretation of our edge detection results. The physical interpretation is important since edge detection algorithms alone do not specify any physical measure, or physical definition, of the leading edge of the spreading population. Our modeling indicates that variations in the image threshold parameter correspond to a consistent variation in the local cell density. This means that varying the threshold parameter is equivalent to varying the location of the leading edge in the range of approximately 1–5% of the maximum cell density.

Procaspase 8 overexpression in non-small-cell lung cancer promotes apoptosis induced by FLIP silencing

We found that procaspase 8 was overexpressed in non-small-cell lung cancers (NSCLCs) compared with matched normal tissues. The caspase 8 inhibitor FLICE-inhibitory protein (FLIP) was also overexpressed in the majority of NSCLCs. Silencing FLIP induced caspase 8 activation and apoptosis in NSCLC cell lines, but not in normal lung cell lines. Apoptosis induced by FLIP silencing was mediated by the TRAIL death receptors DR4 and DR5, but was not dependent on ligation of the receptors by TRAIL. Furthermore, the apoptosis induced by FLIP silencing was dependent on the overexpression of procaspase 8 in NSCLC cells. Moreover, in NSCLC cells, but not in normal cells, FLIP silencing induced co-localization of DR5 and ceramide, and disruption of this co-localization abrogated apoptosis. FLIP silencing supra-additively increased TRAIL-induced apoptosis of NSCLC cells; however, normal lung cells were resistant to TRAIL, even when FLIP was silenced. Importantly, FLIP silencing sensitized NSCLC cells but not normal cells to chemotherapy in vitro, and silencing FLIP in vivo retarded NSCLC xenograft growth and enhanced the anti-tumour effects of cisplatin. Collectively, our results suggest that due to frequent procaspase 8 overexpression, NSCLCs may be particularly sensitive to FLIP-targeted therapies.

Are in vitro estimates of cell diffusivity and cell proliferation rate sensitive to assay geometry?

Cells respond to various biochemical and physical cues during wound–healing and tumour progression. In vitro assays used to study these processes are typically conducted in one particular geometry and it is unclear how the assay geometry affects the capacity of cell populations to spread, or whether the relevant mechanisms, such as cell motility and cell proliferation, are somehow sensitive to the geometry of the assay. In this work we use a circular barrier assay to characterise the spreading of cell populations in two different geometries. Assay 1 describes a tumour–like geometry where a cell population spreads outwards into an open space. Assay 2 describes a wound–like geometry where a cell population spreads inwards to close a void. We use a combination of discrete and continuum mathematical models and automated image processing methods to obtain independent estimates of the effective cell diffusivity, D, and the effective cell proliferation rate, λ. Using our parameterised mathematical model we confirm that our estimates of D and λ accurately predict the time–evolution of the location of the leading edge and the cell density profiles for both assay 1 and assay 2. Our work suggests that the effective cell diffusivity is up to 50% lower for assay 2 compared to assay 1, whereas the effective cell proliferation rate is up to 30% lower for assay 2 compared to assay 1.

Serial explant culture provides novel insights into the potential location and phenotype of corneal endothelial progenitor cells

The routine cultivation of human corneal endothelial cells, with the view to treating patients with endothelial dysfunction, remains a challenging task. While progress in this field has been buoyed by the proposed existence of progenitor cells for the corneal endothelium at the corneal limbus, strategies for exploiting this concept remain unclear. In the course of evaluating methods for growing corneal endothelial cells, we have noted a case where remarkable growth was achieved using a serial explant culture technique. Over the course of 7 months, a single explant of corneal endothelium, acquired from cadaveric human tissue, was sequentially seeded into 7 culture plates and on each occasion produced a confluent cell monolayer. Sample cultures were confirmed as endothelial in origin by positive staining for glypican-4. On each occasion, small cells, closest to the tissue explant, developed into a highly compact layer with an almost homogenous structure. This layer was resistant to removal with trypsin and produced continuous cell outgrowth during multiple culture periods. The small cells gave rise to larger cells with phase-bright cell boundaries and prominent immunostaining for both nestin and telomerase. Nestin and telomerase were also strongly expressed in small cells immediately adjacent to the wound site, following transfer of the explant to another culture plate. These findings are consistent with the theory that progenitor cells for the corneal endothelium reside within the limbus and provide new insights into expected expression patterns for nestin and telomerase within the differentiation pathway.

CD1 expression and CD1-restricted T cell activity in normal and tumour-bearing human liver

CD1d-restricted natural killer T (NKT) cells expressing invariant Valpha14Jalpha18 T cell receptor alpha-chains are abundant in murine liver and are implicated in the control of malignancy, infection and autoimmunity. Invariant NKT cells have potent anti-metastatic effects in mice and phase I clinical trials involving their homologues in humans are ongoing. However, invariant NKT cells are less abundant in human liver ( approximately 0.5% of hepatic T cells) than in murine liver (up to 50%) and it is not known if other hepatic T cells are CD1-restricted. We have examined expression of CD1a, CD1b, CD1c and CD1d mRNA and protein in human liver and evaluated the reactivity of mononuclear cells (MNC) from histologically normal and tumour-bearing human liver specimens against these CD1 isoforms. Messenger RNA for all CD1 isotypes was detectable in all liver samples. CD1c and CD1d were expressed at the protein level by hepatic MNC. CD1d, only, was detectable at the cell surface, but CD1c and CD1d were found at an intracellular location in significant numbers of liver MNC. CD1b was not expressed by MNC from healthy livers but was detectable within MNC in all tumour samples tested. Hepatic T cells exhibited reactivity against C1R cells expressing transfected CD1c and CD1d, but neither CD1a nor CD1b. These cells secreted interferon-gamma (IFN-gamma) but not interleukin-4 (IL-4) upon stimulation. In contrast, similar numbers of peripheral T cells released 13- and 16-fold less IFN-gamma in response to CD1c and CD1d, respectively. CD1c and CD1d expression and T cell reactivity were not altered in tumour-bearing liver specimens compared to histologically normal livers. These data suggest that, in addition to invariant CD1d-restricted NKT cells, autoreactive T cells that recognise CD1c and CD1d and release inflammatory cytokines are abundant in human liver.

Effective Cell-Centred Time-Domain Maxwell's Equations Numerical Solvers

This research work analyses techniques for implementing a cell-centred finite-volume time-domain (ccFV-TD) computational methodology for the purpose of studying microwave heating. Various state-of-the-art spatial and temporal discretisation methods employed to solve Maxwell's equations on multidimensional structured grid networks are investigated, and the dispersive and dissipative errors inherent in those techniques examined. Both staggered and unstaggered grid approaches are considered. Upwind schemes using a Riemann solver and intensity vector splitting are studied and evaluated. Staggered and unstaggered Leapfrog and Runge-Kutta time integration methods are analysed in terms of phase and amplitude error to identify which method is the most accurate and efficient for simulating microwave heating processes. The implementation and migration of typical electromagnetic boundary conditions. from staggered in space to cell-centred approaches also is deliberated. In particular, an existing perfectly matched layer absorbing boundary methodology is adapted to formulate a new cell-centred boundary implementation for the ccFV-TD solvers. Finally for microwave heating purposes, a comparison of analytical and numerical results for standard case studies in rectangular waveguides allows the accuracy of the developed methods to be assessed.