Framework for morphometric classification of cells in imaging flow cytometry

Imaging flow cytometry is an emerging technology that combines the statistical power of flow cytometry with spatial and quantitative morphology of digital microscopy. It allows high-throughput imaging of cells with good spatial resolution, while they are in flow. This paper proposes a general framework for the processing/classification of cells imaged using imaging flow cytometer. Each cell is localized by finding an accurate cell contour. Then, features reflecting cell size, circularity and complexity are extracted for the classification using SVM. Unlike the conventional iterative, semi-automatic segmentation algorithms such as active contour, we propose a noniterative, fully automatic graph-based cell localization. In order to evaluate the performance of the proposed framework, we have successfully classified unstained label-free leukaemia cell-lines MOLT, K562 and HL60 from video streams captured using custom fabricated cost-effective microfluidics-based imaging flow cytometer. The proposed system is a significant development in the direction of building a cost-effective cell analysis platform that would facilitate affordable mass screening camps looking cellular morphology for disease diagnosis. Lay description In this article, we propose a novel framework for processing the raw data generated using microfluidics based imaging flow cytometers. Microfluidics microscopy or microfluidics based imaging flow cytometry (mIFC) is a recent microscopy paradigm, that combines the statistical power of flow cytometry with spatial and quantitative morphology of digital microscopy, which allows us imaging cells while they are in flow. In comparison to the conventional slide-based imaging systems, mIFC is a nascent technology enabling high throughput imaging of cells and is yet to take the form of a clinical diagnostic tool. The proposed framework process the raw data generated by the mIFC systems. The framework incorporates several steps: beginning from pre-processing of the raw video frames to enhance the contents of the cell, localising the cell by a novel, fully automatic, non-iterative graph based algorithm, extraction of different quantitative morphological parameters and subsequent classification of cells. In order to evaluate the performance of the proposed framework, we have successfully classified unstained label-free leukaemia cell-lines MOLT, K562 and HL60 from video streams captured using cost-effective microfluidics based imaging flow cytometer. The cell lines of HL60, K562 and MOLT were obtained from ATCC (American Type Culture Collection) and are separately cultured in the lab. Thus, each culture contains cells from its own category alone and thereby provides the ground truth. Each cell is localised by finding a closed cell contour by defining a directed, weighted graph from the Canny edge images of the cell such that the closed contour lies along the shortest weighted path surrounding the centroid of the cell from a starting point on a good curve segment to an immediate endpoint. Once the cell is localised, morphological features reflecting size, shape and complexity of the cells are extracted and used to develop a support vector machine based classification system. We could classify the cell-lines with good accuracy and the results were quite consistent across different cross validation experiments. We hope that imaging flow cytometers equipped with the proposed framework for image processing would enable cost-effective, automated and reliable disease screening in over-loaded facilities, which cannot afford to hire skilled personnel in large numbers. Such platforms would potentially facilitate screening camps in low income group countries; thereby transforming the current health care paradigms by enabling rapid, automated diagnosis for diseases like cancer.

Einstein Relation in n-Channel Inversion Layers of Nonlinear Optical, Optoelectronic and Related Materials: Simplified Theory, Relative Comparison and Suggestion for an Experimental Determination

In this paper, we study the Einstein relation for the diffusivity to mobility ratio (DMR) in n-channel inversion layers of non-linear optical materials on the basis of a newly formulated electron dispersion relation by considering their special properties within the frame work of k.p formalism. The results for the n-channel inversion layers of III-V, ternary and quaternary materials form a special case of our generalized analysis. The DMR for n-channel inversion layers of II-VI, IV-VI and stressed materials has been investigated by formulating the respective 2D electron dispersion laws. It has been found, taking n-channel inversion layers of CdGeAs2, Cd(3)AS(2), InAs, InSb, Hg1-xCdxTe, In1-xGaxAsyP1-y lattice matched to InP, CdS, PbTe, PbSnTe, Pb1-xSnxSe and stressed InSb as examples, that the DMR increases with the increasing surface electric field with different numerical values and the nature of the variations are totally band structure dependent. The well-known expression of the DMR for wide gap materials has been obtained as a special case under certain limiting conditions and this compatibility is an indirect test for our generalized formalism. Besides, an experimental method of determining the 2D DMR for n-channel inversion layers having arbitrary dispersion laws has been suggested.

Species-area and species-individual relationships for tropical trees: A comparison of three 50-ha plots

1 Species-accumulation curves for woody plants were calculated in three tropical forests, based on fully mapped 50-ha plots in wet, old-growth forest in Peninsular Malaysia, in moist, old-growth forest in central Panama, and in dry, previously logged forest in southern India. A total of 610 000 stems were identified to species and mapped to < Im accuracy. Mean species number and stem number were calculated in quadrats as small as 5 m x 5 m to as large as 1000 m x 500 m, for a variety of stem sizes above 10 mm in diameter. Species-area curves were generated by plotting species number as a function of quadrat size; species-individual curves were generated from the same data, but using stem number as the independent variable rather than area. 2 Species-area curves had different forms for stems of different diameters, but species-individual curves were nearly independent of diameter class. With < 10(4) stems, species-individual curves were concave downward on log-log plots, with curves from different forests diverging, but beyond about 104 stems, the log-log curves became nearly linear, with all three sites having a similar slope. This indicates an asymptotic difference in richness between forests: the Malaysian site had 2.7 times as many species as Panama, which in turn was 3.3 times as rich as India. 3 Other details of the species-accumulation relationship were remarkably similar between the three sites. Rectangular quadrats had 5-27% more species than square quadrats of the same area, with longer and narrower quadrats increasingly diverse. Random samples of stems drawn from the entire 50 ha had 10-30% more species than square quadrats with the same number of stems. At both Pasoh and BCI, but not Mudumalai. species richness was slightly higher among intermediate-sized stems (50-100mm in diameter) than in either smaller or larger sizes, These patterns reflect aggregated distributions of individual species, plus weak density-dependent forces that tend to smooth the species abundance distribution and 'loosen' aggregations as stems grow. 4 The results provide support for the view that within each tree community, many species have their abundance and distribution guided more by random drift than deterministic interactions. The drift model predicts that the species-accumulation curve will have a declining slope on a log-log plot, reaching a slope of O.1 in about 50 ha. No other model of community structure can make such a precise prediction. 5 The results demonstrate that diversity studies based on different stem diameters can be compared by sampling identical numbers of stems. Moreover, they indicate that stem counts < 1000 in tropical forests will underestimate the percentage difference in species richness between two diverse sites. Fortunately, standard diversity indices (Fisher's sc, Shannon-Wiener) captured diversity differences in small stem samples more effectively than raw species richness, but both were sample size dependent. Two nonparametric richness estimators (Chao. jackknife) performed poorly, greatly underestimating true species richness.

Purification and functional characterization of type II DNA topoisomerase from rat testis and comparison with topoisomerase II from liver

A number of studies in yeast have shown that DNA topoisomerase TI is essential for chromosome condensation and disjunction during mitosis at the metaphase/anaphase transition and meiosis I. Accordingly, kinetic and mechanistic studies have implied a role for topoisomerase rr in chromosome disjunction. As a step toward understanding the nature and role of topoisomerase II in a mammalian germline in vivo, we have purified topoisomerase II from rat testis to homogeneity and ascertained several of its catalytic activities in conjunction with that of the purified enzyme from liver. The purified enzymes appeared to be monomers under denaturing conditions; however, they differed in their relative molecular mass. Topoisomerase II from testis and liver have apparent molecular masses of 150 +/- 10 kDa and 160 +/- 10 kDa, respectively. The native molecular mass of testis topoisomerase II as assayed by immunoblot analysis of cell-foe extracts, prepared in the presence of SDS and a number of protease inhibitors, corroborated with the size of the purified enzyme. Both enzymes are able to promote decatenation and relax supercoiled DNA substrates in an ATP and Mg2+-dependent manner. However, quantitative comparison of catalytic properties of topoisomerase II from testis with that of the enzyme from liver displayed significant differences in their efficiencies. Optimal pH values for testis enzyme are 6.5 to 8.5 while they are 6 to 7.5 for the liver enzyme. Intriguingly, the relaxation activity of liver topoisomerase II was inhibited by potassium glutamate at 1 M, whereas testis enzyme required about half its concentration. These findings argue that topoisomerase II from rat testis is structurally distinct from that of its somatic form and the functional differences between the two enzymes parallels with the physiological environment that is unique to these two tissues.

Charge transport in transparent conductors: A comparison

A comparative investigation of charge transport properties is presented, for polymeric [poly(3,4-ethylenedioxythiophene):poly(styrenesulfonate) (PEDOT:PSS)], single-wall carbon nanotube (SWNT) and inorganic (indium tin oxide, ITO), transparent conducting electrodes. The polymeric and nanotube systems show hopping transport at low temperatures, in contrast with the disordered-metal transport in ITO. The low temperature magnetotransport (up to 11 T) and high electric-field transport (up to 500 V/cm) indicate the significant role of nanoscopic scale disorder for charge transport in polymer and nanotube based systems. The results show that characteristic length scales like localization length correlates with the nanomorphology in these systems. Further, the high frequency conductivity measurements (up to 30 MHz) in PEDOT:PSS and SWNT follow the extended pair approximation model [σ(ω)=σ(0)[1+(ω/ω0)s].

Comparison of Interproton Distances in DNA Models with Nuclear Overhauser Enhancement Data

The conformational flexibility inherent in the polynucleotide chain plays an important role in deciding its three-dimensonal structure and enables it to undergo structural transitions in order to fulfil all its functions. Following certain stereochemical guidelines, both right and left handed double-helical models have been built in our laboratory and they are in reasonably good agreement with the fibre patterns for various polymorphous forms of DNA. Recently, nuclear magnetic resonance spectroscopy has become an important technique for studying the solution conformation and polymorphism of nucleic acids. Several workers have used 1H nuclear magnetic resonance nuclear Overhauser enhancement measurements to estimate the interproton distances for the various DNA oligomers and compared them with the interproton distances for particular models of A and Β form DNA. In some cases the solution conformation does not seem to fit either of these models. We have been studying various models for DNA with a view to exploring the full conformational space allowed for nucleic acid polymers. In this paper, the interproton distances calculated for the different stereochemically feasible models of DNA are presented and they are compared and correlated against those obtained from 1Η nuclear magnetic resonance nuclear Overhauser enhancement measurements of various nucleic acid oligomers.

A comparison of the molecular orbital (ab initio, CNDO and EHT) methods applied to the conformational study of thiocarbonohydrazide in its basic and diprotonated forms

Conformational preferences of thiocarbonohydrazide (H2NNHCSNHNH2) in its basic and N,N′-diprotonated forms are examined by calculating the barrier to internal rotation around the C---N bonds, using the theoretical LCAO—MO (ab initio and semiempirical CNDO and EHT) methods. The calculated and experimental results are compared with each other and also with values for N,N′-dimethylthiourea which is isoelectronic with thiocarbonohydrazide. The suitability of these methods for studying rotational isomerism seems suspect when lone pair interactions are present.

Comparison of Rectangular and Polar Quantizers in Digital Holography

Complex amplitude encoded in any digital hologram must undergo quantization, usually in either polar or rectangular format . In this paper these two schemes are compared under the constraints and conditions inherent in digital holography . For Fourier transform holograms when the spectrum is levelled through phase coding, the rectangular format is shown to be optimal . In the absence of phase coding, and also if the amplitude spectrum has a large dynamic range, the polar format may be preferable .

Comparison Of If Estimation Methods For Noise Robustness

Non-stationary signal modeling is a well addressed problem in the literature. Many methods have been proposed to model non-stationary signals such as time varying linear prediction and AM-FM modeling, the later being more popular. Estimation techniques to determine the AM-FM components of narrow-band signal, such as Hilbert transform, DESA1, DESA2, auditory processing approach, ZC approach, etc., are prevalent but their robustness to noise is not clearly addressed in the literature. This is critical for most practical applications, such as in communications. We explore the robustness of different AM-FM estimators in the presence of white Gaussian noise. Also, we have proposed three new methods for IF estimation based on non-uniform samples of the signal and multi-resolution analysis. Experimental results show that ZC based methods give better results than the popular methods such as DESA in clean condition as well as noisy condition.

Estrogen modulation of retinol-binding protein in immature chicks: Comparison with riboflavin carrier protein

The kinetics of estrogen (E) modulation of retinol-binding protein (RBP) production in the liver of immature chicks were compared with those governing de novo induction of riboflavin carrier protein (RCP) in the same tissue. A single dose of E markedly enhanced the plasma levels of RBP without any detectable lag period to reach peak value by 24 h and this was followed by a decline to attain the baseline by 4 days. There was no amplification of the response during secondary stimulation unlike the case with RCP induction. With multiple E administration, the 4-fold increased plasma RBP concentrations were sustained at a steady state during both primary and secondary stimulations, whereas concomitant RCP concentration progressively increased with each hormone administration and this response was further amplified during secondary stimulation. Unlike RCP induction, enhanced RBP accumulation was not strictly E dose dependent although a minimal threshold level of the steroid was required to elicit measurable response. Progesterone (P) could neither modulate nor substitute for E in enhancing plasma levels of either of the 2 proteins while the anti-estrogens, en- and zuclomifene citrate severely suppressed the production of both the proteins. RCP induction was completely inhibited by both α-amanitin and cycloheximide for prolonged periods while E-stimulated RBP production was affected only partially by α-amanitin. Likewise, cycloheximide inhibition of RBP accumulation followed a pattern similar to that of hepatic general protein synthesis.

Comparison of Elasticity, Shell Core, and Sandwich Shell Theories

The problem of an infinite circular sandwich shell subjected to an a\isymmetric radial line load is investigated using three-dimensional elasticity theory, shell core method, and sandwich shell theory due to Fulton and Schmidt. A comparison of the stresses and displacements with an exact elasticity solution is carried out for various shell parameters in order to clearly bring out the limitations of sandwich shell theories of Fulton and Schmidt as well as the shell core solution.

How Good is the Assumption About Visible Surface Reflectance in MODIS Aerosol Retrieval Over Land? A Comparison With Aircraft Measurements Over an Urban Site in India

In this paper, we discuss the measurements of spectral surface reflectance (rho(s)(lambda)) in the wavelength range 350-2500 nm measured using a spectroradiometer onboard a low-flying aircraft over Bangalore (12.95 degrees N, 77.65 degrees E), an urban site in southern India. The large discrepancies in the retrieval of aerosol propertiesover land by the Moderate-Resolution Imaging Spectroradiometer (MODIS), which could be attributed to the inaccurate estimation of surface reflectance at many sites in India and elsewhere, provided motivation for this paper. The aim of this paper was to verify the surface reflectance relationships assumed by the MODIS aerosol algorithm for the estimation of surface reflectance in the visible channels (470 and 660 nm) from the surface reflectance at 2100 nm for aerosol retrieval over land. The variety of surfaces observed in this paper includes green and dry vegetations, bare land, and urban surfaces. The measuredreflectance data were first corrected for the radiative effects of atmosphere lying between the ground and aircraft using the Second Simulation of Satellite Signal in the Solar Spectrum (6S) radiative transfer code. The corrected surface reflectance in the MODIS's blue (rho(s)(470)), red (rho(s)(660)), and shortwave-infrared (SWIR) channel (rho(s)(2100)) was linearly correlated. We found that the slope of reflectance relationship between 660 and 2100 nm derived from the forward scattering data was 0.53 with an intercept of 0.07, whereas the slope for the relationship between the reflectance at 470 and 660 nm was 0.85. These values are much higher than the slope (similar to 0.49) for either wavelengths assumed by the MODIS aerosol algorithm over this region. The reflectance relationship for the backward scattering data has a slope of 0.39, with an intercept of 0.08 for 660 nm, and 0.65, with an intercept of 0.08 for 470 nm. The large values of the intercept (which is very small in the MODIS reflectance relationships) result in larger values of absolute surface reflectance in the visible channels. The discrepancy between the measured and assumed surface reflectances could lead to error in the aerosol retrieval. The reflectance ratio (rho(s)(660)/rho(s)(2100)) showed a clear dependence on the N D V I-SWIR where the ratio increased from 0.5 to 1 with an increase in N V I-SWIR from 0 to 0.5. The high correlation between the reflectance at SWIR wavelengths (2100, 1640, and 1240 nm) indicated an opportunity to derive the surface reflectance and, possibly, aerosol properties at these wavelengths. We need more experiments to characterize the surface reflectance and associated inhomogeneity of land surfaces, which play a critical role in the remote sensing of aerosols over land.

Hormonal induction of riboflavin carrier protein in the chicken oviduct and liver: a comparison of kinetics and modulation

Estrogen (E) induction of riboflavin carrier protein (RCP) in the chicken oviduct and liver was investigated to compare and contrast the kinetics, hormonal specificity and modulation of its elaboration in the 2 steroid-responsive tissues. During primary stimulation, continued daily E administration to immature female chicks elicited, after an initial lag, rapid growth and RCP content of the oviduct; neither progesterone (P) nor testosterone (T) could substitute for E in this respect. Furthermore, P given along with E curtailed tissue growth and its RCP content, whereas E + T had a synergistic effect on tissue growth only. During secondary stimulation, E administration steeply enhanced both tissue weight and RCP content without any lag. Interestingly, P (but not T) could substitute for E in augmenting magnum RCP concentration to a comparable extent while a concomitant effect on tissue growth was less marked. In contrast, hepatic induction of RCP was absolutely E-specific during both primary and secondary stimulations. Secondary stimulation with either E or P of E-primed birds enhanced the rates of RCP synthesis in the oviduct relative to that of total protein, whereas in the liver only E was effective in this regard. The absolute rate of E-induced RCP synthesis in both the steroid-stimulated tissues was significantly higher than that of general protein elaboration.

Short-window spectral analysis using AMVAR and multitaper methods: a comparison

We compare two popular methods for estimating the power spectrum from short data windows, namely the adaptive multivariate autoregressive (AMVAR) method and the multitaper method. By analyzing a simulated signal (embedded in a background Ornstein-Uhlenbeck noise process) we demonstrate that the AMVAR method performs better at detecting short bursts of oscillations compared to the multitaper method. However, both methods are immune to jitter in the temporal location of the signal. We also show that coherence can still be detected in noisy bivariate time series data by the AMVAR method even if the individual power spectra fail to show any peaks. Finally, using data from two monkeys performing a visuomotor pattern discrimination task, we demonstrate that the AMVAR method is better able to determine the termination of the beta oscillations when compared to the multitaper method.