Micro-Anatomical Quantitative Imaging Towards Enabling Automated Diagnosis of Thick Tissues at the Point of Care

Histopathology is the clinical standard for tissue diagnosis. However, histopathology has several limitations including that it requires tissue processing, which can take 30 minutes or more, and requires a highly trained pathologist to diagnose the tissue. Additionally, the diagnosis is qualitative, and the lack of quantitation leads to possible observer-specific diagnosis. Taken together, it is difficult to diagnose tissue at the point of care using histopathology.

Several clinical situations could benefit from more rapid and automated histological processing, which could reduce the time and the number of steps required between obtaining a fresh tissue specimen and rendering a diagnosis. For example, there is need for rapid detection of residual cancer on the surface of tumor resection specimens during excisional surgeries, which is known as intraoperative tumor margin assessment. Additionally, rapid assessment of biopsy specimens at the point-of-care could enable clinicians to confirm that a suspicious lesion is successfully sampled, thus preventing an unnecessary repeat biopsy procedure. Rapid and low cost histological processing could also be potentially useful in settings lacking the human resources and equipment necessary to perform standard histologic assessment. Lastly, automated interpretation of tissue samples could potentially reduce inter-observer error, particularly in the diagnosis of borderline lesions.

To address these needs, high quality microscopic images of the tissue must be obtained in rapid timeframes, in order for a pathologic assessment to be useful for guiding the intervention. Optical microscopy is a powerful technique to obtain high-resolution images of tissue morphology in real-time at the point of care, without the need for tissue processing. In particular, a number of groups have combined fluorescence microscopy with vital fluorescent stains to visualize micro-anatomical features of thick (i.e. unsectioned or unprocessed) tissue. However, robust methods for segmentation and quantitative analysis of heterogeneous images are essential to enable automated diagnosis. Thus, the goal of this work was to obtain high resolution imaging of tissue morphology through employing fluorescence microscopy and vital fluorescent stains and to develop a quantitative strategy to segment and quantify tissue features in heterogeneous images, such as nuclei and the surrounding stroma, which will enable automated diagnosis of thick tissues.

To achieve these goals, three specific aims were proposed. The first aim was to develop an image processing method that can differentiate nuclei from background tissue heterogeneity and enable automated diagnosis of thick tissue at the point of care. A computational technique called sparse component analysis (SCA) was adapted to isolate features of interest, such as nuclei, from the background. SCA has been used previously in the image processing community for image compression, enhancement, and restoration, but has never been applied to separate distinct tissue types in a heterogeneous image. In combination with a high resolution fluorescence microendoscope (HRME) and a contrast agent acriflavine, the utility of this technique was demonstrated through imaging preclinical sarcoma tumor margins. Acriflavine localizes to the nuclei of cells where it reversibly associates with RNA and DNA. Additionally, acriflavine shows some affinity for collagen and muscle. SCA was adapted to isolate acriflavine positive features or APFs (which correspond to RNA and DNA) from background tissue heterogeneity. The circle transform (CT) was applied to the SCA output to quantify the size and density of overlapping APFs. The sensitivity of the SCA+CT approach to variations in APF size, density and background heterogeneity was demonstrated through simulations. Specifically, SCA+CT achieved the lowest errors for higher contrast ratios and larger APF sizes. When applied to tissue images of excised sarcoma margins, SCA+CT correctly isolated APFs and showed consistently increased density in tumor and tumor + muscle images compared to images containing muscle. Next, variables were quantified from images of resected primary sarcomas and used to optimize a multivariate model. The sensitivity and specificity for differentiating positive from negative ex vivo resected tumor margins was 82% and 75%. The utility of this approach was further tested by imaging the in vivo tumor cavities from 34 mice after resection of a sarcoma with local recurrence as a bench mark. When applied prospectively to images from the tumor cavity, the sensitivity and specificity for differentiating local recurrence was 78% and 82%. The results indicate that SCA+CT can accurately delineate APFs in heterogeneous tissue, which is essential to enable automated and rapid surveillance of tissue pathology.

Two primary challenges were identified in the work in aim 1. First, while SCA can be used to isolate features, such as APFs, from heterogeneous images, its performance is limited by the contrast between APFs and the background. Second, while it is feasible to create mosaics by scanning a sarcoma tumor bed in a mouse, which is on the order of 3-7 mm in any one dimension, it is not feasible to evaluate an entire human surgical margin. Thus, improvements to the microscopic imaging system were made to (1) improve image contrast through rejecting out-of-focus background fluorescence and to (2) increase the field of view (FOV) while maintaining the sub-cellular resolution needed for delineation of nuclei. To address these challenges, a technique called structured illumination microscopy (SIM) was employed in which the entire FOV is illuminated with a defined spatial pattern rather than scanning a focal spot, such as in confocal microscopy.

Thus, the second aim was to improve image contrast and increase the FOV through employing wide-field, non-contact structured illumination microscopy and optimize the segmentation algorithm for new imaging modality. Both image contrast and FOV were increased through the development of a wide-field fluorescence SIM system. Clear improvement in image contrast was seen in structured illumination images compared to uniform illumination images. Additionally, the FOV is over 13X larger than the fluorescence microendoscope used in aim 1. Initial segmentation results of SIM images revealed that SCA is unable to segment large numbers of APFs in the tumor images. Because the FOV of the SIM system is over 13X larger than the FOV of the fluorescence microendoscope, dense collections of APFs commonly seen in tumor images could no longer be sparsely represented, and the fundamental sparsity assumption associated with SCA was no longer met. Thus, an algorithm called maximally stable extremal regions (MSER) was investigated as an alternative approach for APF segmentation in SIM images. MSER was able to accurately segment large numbers of APFs in SIM images of tumor tissue. In addition to optimizing MSER for SIM image segmentation, an optimal frequency of the illumination pattern used in SIM was carefully selected because the image signal to noise ratio (SNR) is dependent on the grid frequency. A grid frequency of 31.7 mm-1 led to the highest SNR and lowest percent error associated with MSER segmentation.

Once MSER was optimized for SIM image segmentation and the optimal grid frequency was selected, a quantitative model was developed to diagnose mouse sarcoma tumor margins that were imaged ex vivo with SIM. Tumor margins were stained with acridine orange (AO) in aim 2 because AO was found to stain the sarcoma tissue more brightly than acriflavine. Both acriflavine and AO are intravital dyes, which have been shown to stain nuclei, skeletal muscle, and collagenous stroma. A tissue-type classification model was developed to differentiate localized regions (75x75 µm) of tumor from skeletal muscle and adipose tissue based on the MSER segmentation output. Specifically, a logistic regression model was used to classify each localized region. The logistic regression model yielded an output in terms of probability (0-100%) that tumor was located within each 75x75 µm region. The model performance was tested using a receiver operator characteristic (ROC) curve analysis that revealed 77% sensitivity and 81% specificity. For margin classification, the whole margin image was divided into localized regions and this tissue-type classification model was applied. In a subset of 6 margins (3 negative, 3 positive), it was shown that with a tumor probability threshold of 50%, 8% of all regions from negative margins exceeded this threshold, while over 17% of all regions exceeded the threshold in the positive margins. Thus, 8% of regions in negative margins were considered false positives. These false positive regions are likely due to the high density of APFs present in normal tissues, which clearly demonstrates a challenge in implementing this automatic algorithm based on AO staining alone.

Thus, the third aim was to improve the specificity of the diagnostic model through leveraging other sources of contrast. Modifications were made to the SIM system to enable fluorescence imaging at a variety of wavelengths. Specifically, the SIM system was modified to enabling imaging of red fluorescent protein (RFP) expressing sarcomas, which were used to delineate the location of tumor cells within each image. Initial analysis of AO stained panels confirmed that there was room for improvement in tumor detection, particularly in regards to false positive regions that were negative for RFP. One approach for improving the specificity of the diagnostic model was to investigate using a fluorophore that was more specific to staining tumor. Specifically, tetracycline was selected because it appeared to specifically stain freshly excised tumor tissue in a matter of minutes, and was non-toxic and stable in solution. Results indicated that tetracycline staining has promise for increasing the specificity of tumor detection in SIM images of a preclinical sarcoma model and further investigation is warranted.

In conclusion, this work presents the development of a combination of tools that is capable of automated segmentation and quantification of micro-anatomical images of thick tissue. When compared to the fluorescence microendoscope, wide-field multispectral fluorescence SIM imaging provided improved image contrast, a larger FOV with comparable resolution, and the ability to image a variety of fluorophores. MSER was an appropriate and rapid approach to segment dense collections of APFs from wide-field SIM images. Variables that reflect the morphology of the tissue, such as the density, size, and shape of nuclei and nucleoli, can be used to automatically diagnose SIM images. The clinical utility of SIM imaging and MSER segmentation to detect microscopic residual disease has been demonstrated by imaging excised preclinical sarcoma margins. Ultimately, this work demonstrates that fluorescence imaging of tissue micro-anatomy combined with a specialized algorithm for delineation and quantification of features is a means for rapid, non-destructive and automated detection of microscopic disease, which could improve cancer management in a variety of clinical scenarios.

Predictive Modeling of Loss-of-Heterozygosity Events in Yeast

During mitotic cell cycles, DNA experiences many types of endogenous and exogenous damaging agents that could potentially cause double strand breaks (DSB). In S. cerevisiae, DSBs are primarily repaired by mitotic recombination and as a result, could lead to loss-of-heterozygosity (LOH). Genetic recombination can happen in both meiosis and mitosis. While genome-wide distribution of meiotic recombination events has been intensively studied, mitotic recombination events have not been mapped unbiasedly throughout the genome until recently. Methods for selecting mitotic crossovers and mapping the positions of crossovers have recently been developed in our lab. Our current approach uses a diploid yeast strain that is heterozygous for about 55,000 SNPs, and employs SNP-Microarrays to map LOH events throughout the genome. These methods allow us to examine selected crossovers and unselected mitotic recombination events (crossover, noncrossover and BIR) at about 1 kb resolution across the genome. Using this method, we generated maps of spontaneous and UV-induced LOH events. In this study, we explore machine learning and variable selection techniques to build a predictive model for where the LOH events occur in the genome.

Randomly from the yeast genome, we simulated control tracts resembling the LOH tracts in terms of tract lengths and locations with respect to single-nucleotide-polymorphism positions. We then extracted roughly 1,100 features such as base compositions, histone modifications, presence of tandem repeats etc. and train classifiers to distinguish control tracts and LOH tracts. We found interesting features of good predictive values. We also found that with the current repertoire of features, the prediction is generally better for spontaneous LOH events than UV-induced LOH events.

Modulation of bacterial outer membrane vesicle production by envelope structure and content.

BACKGROUND: Vesiculation is a ubiquitous secretion process of Gram-negative bacteria, where outer membrane vesicles (OMVs) are small spherical particles on the order of 50 to 250 nm composed of outer membrane (OM) and lumenal periplasmic content. Vesicle functions have been elucidated in some detail, showing their importance in virulence factor secretion, bacterial survival, and biofilm formation in pathogenesis. Furthermore, OMVs serve as an envelope stress response, protecting the secreting bacteria from internal protein misfolding stress, as well as external envelope stressors. Despite their important functional roles very little is known about the regulation and mechanism of vesicle production. Based on the envelope architecture and prior characterization of the hypervesiculation phenotypes for mutants lacking the lipoprotein, Lpp, which is involved in the covalent OM-peptidoglycan (PG) crosslinks, it is expected that an inverse relationship exists between OMV production and PG-crosslinked Lpp. RESULTS: In this study, we found that subtle modifications of PG remodeling and crosslinking modulate OMV production, inversely correlating with bound Lpp levels. However, this inverse relationship was not found in strains in which OMV production is driven by an increase in "periplasmic pressure" resulting from the accumulation of protein, PG fragments, or lipopolysaccharide. In addition, the characterization of an nlpA deletion in backgrounds lacking either Lpp- or OmpA-mediated envelope crosslinks demonstrated a novel role for NlpA in envelope architecture. CONCLUSIONS: From this work, we conclude that OMV production can be driven by distinct Lpp concentration-dependent and Lpp concentration-independent pathways.

Interactions of chromatin context, binding site sequence content, and sequence evolution in stress-induced p53 occupancy and transactivation.

Cellular stresses activate the tumor suppressor p53 protein leading to selective binding to DNA response elements (REs) and gene transactivation from a large pool of potential p53 REs (p53REs). To elucidate how p53RE sequences and local chromatin context interact to affect p53 binding and gene transactivation, we mapped genome-wide binding localizations of p53 and H3K4me3 in untreated and doxorubicin (DXR)-treated human lymphoblastoid cells. We examined the relationships among p53 occupancy, gene expression, H3K4me3, chromatin accessibility (DNase 1 hypersensitivity, DHS), ENCODE chromatin states, p53RE sequence, and evolutionary conservation. We observed that the inducible expression of p53-regulated genes was associated with the steady-state chromatin status of the cell. Most highly inducible p53-regulated genes were suppressed at baseline and marked by repressive histone modifications or displayed CTCF binding. Comparison of p53RE sequences residing in different chromatin contexts demonstrated that weaker p53REs resided in open promoters, while stronger p53REs were located within enhancers and repressed chromatin. p53 occupancy was strongly correlated with similarity of the target DNA sequences to the p53RE consensus, but surprisingly, inversely correlated with pre-existing nucleosome accessibility (DHS) and evolutionary conservation at the p53RE. Occupancy by p53 of REs that overlapped transposable element (TE) repeats was significantly higher (p<10-7) and correlated with stronger p53RE sequences (p<10-110) relative to nonTE-associated p53REs, particularly for MLT1H, LTR10B, and Mer61 TEs. However, binding at these elements was generally not associated with transactivation of adjacent genes. Occupied p53REs located in L2-like TEs were unique in displaying highly negative PhyloP scores (predicted fast-evolving) and being associated with altered H3K4me3 and DHS levels. These results underscore the systematic interaction between chromatin status and p53RE context in the induced transactivation response. This p53 regulated response appears to have been tuned via evolutionary processes that may have led to repression and/or utilization of p53REs originating from primate-specific transposon elements.

The Dietary Approaches to Stop Hypertension (DASH) eating pattern in special populations.

The Dietary Approaches to Stop Hypertension (DASH) trial showed that a diet rich in fruits, vegetables, low-fat dairy products with reduced total and saturated fat, cholesterol, and sugar-sweetened products effectively lowers blood pressure in individuals with prehypertension and stage I hypertension. Limited evidence is available on the safety and efficacy of the DASH eating pattern in special patient populations that were excluded from the trial. Caution should be exercised before initiating the DASH diet in patients with chronic kidney disease, chronic liver disease, and those who are prescribed renin-angiotensin-aldosterone system antagonist, but these conditions are not strict contraindications to DASH. Modifications to the DASH diet may be necessary to facilitate its use in patients with chronic heart failure, uncontrolled diabetes mellitus type II, lactose intolerance, and celiac disease. In general, the DASH diet can be adopted by most patient populations and initiated simultaneously with medication therapy and other lifestyle interventions.

Brain evolution by brain pathway duplication.

Understanding the mechanisms of evolution of brain pathways for complex behaviours is still in its infancy. Making further advances requires a deeper understanding of brain homologies, novelties and analogies. It also requires an understanding of how adaptive genetic modifications lead to restructuring of the brain. Recent advances in genomic and molecular biology techniques applied to brain research have provided exciting insights into how complex behaviours are shaped by selection of novel brain pathways and functions of the nervous system. Here, we review and further develop some insights to a new hypothesis on one mechanism that may contribute to nervous system evolution, in particular by brain pathway duplication. Like gene duplication, we propose that whole brain pathways can duplicate and the duplicated pathway diverge to take on new functions. We suggest that one mechanism of brain pathway duplication could be through gene duplication, although other mechanisms are possible. We focus on brain pathways for vocal learning and spoken language in song-learning birds and humans as example systems. This view presents a new framework for future research in our understanding of brain evolution and novel behavioural traits.

Algae-Based Purification of Landfill Biogas Using a CO2 Removal System and Helical Photobioreactor in Series

Biogas is a mixture of methane and other gases. In its crude state, it contains carbon dioxide (CO2) that reduces its energy efficiency and hydrogen sulfide (H2S) that is toxic and highly corrosive. Because chemical methods of removal are expensive and environmentally hazardous, this project investigated an algal-based system to remove CO2 from biogas. An anaerobic digester was used to mimic landfill biogas. Iron oxide and an alkaline spray were used to remove H2S and CO2 respectively. The CO2-laden alkali solution was added to a helical photobioreactor where the algae metabolized the dissolved CO2 to generate algal biomass. Although technical issues prevented testing of the complete system for functionality, cost analysis was completed and showed that the system, in its current state, is not economically feasible. However, modifications may reduce operation costs.

The Ets transcription factors of the PEA3 group: Transcriptional regulators in metastasis

The PEA3 group is composed of three highly conserved Ets transcription factors: Erm, Er81, and Pea3. These proteins regulate transcription of multiple genes, and their transactivating potential is affected by post-translational modifications. Among their target genes are several matrix metalloproteases (MMPs), which are enzymes degrading the extracellular matrix during normal remodelling events and cancer metastasis. In fact, PEA3-group genes are often over-expressed in different types of cancers that also over-express these MMPs and display a disseminating phenotype. Experimental regulation of the synthesis of PEA3 group members influences the metastatic process. This suggests that these factors play a key role in metastasis. © 2006 Elsevier B.V. All rights reserved.

Methylation-dependent T cell immunity to Mycobacterium tuberculosis heparin-binding hemagglutinin.

Although post-translational modifications of protein antigens may be important componenets of some B cell epitopes, the determinants of T cell immunity are generally nonmodified peptides. Here we show that methylation of the Mycobacterium tuberculosis heparin-binding hemagglutinin (HBHA) by the bacterium is essential for effective T cell immunity to this antigen in infected healthy humans and in mice. Methylated HBHA provides high levels of protection against M. tuberculosis challenge in mice, whereas nonmethylated HBHA does not. Protective immunity induced by methylated HBHA is comparable to that afforded by vaccination with bacille Calmette et Guérin, the only available anti-tuberculosis vaccine. Thus, post-translational modifications of proteins may be crucial for their ability to induce protective T cell-mediated immunity against infectious diseases such as tuberculosis.

Pathogenicity of Salmonella: SopE-mediated membrane ruffling is independent of inositol phosphate signals

Studies [Zhou, D. Chen, L.-M. Hernandez, L. Shears, S.B. and Galán, J.E. (2001) A Salmonella inositol polyphosphatase acts in conjunction with other bacterial effectors to promote host-cell actin cytoskeleton rearrangements and bacterial internalization. Mol. Microbiol. 39, 248-259] with engineered Salmonella mutants showed that deletion of SopE attenuated the pathogen's ability to deplete host-cell InsP5 and remodel the cytoskeleton. We pursued these observations: In SopE-transfected host-cells, membrane ruffling was induced, but SopE did not dephosphorylate InsP5, nor did it recruit PTEN (a cytosolic InsP5 phosphatase) for this task. However, PTEN strengthened SopE-mediated membrane ruffling. We conclude SopE promotes host-cell InsP5 hydrolysis only with the assistance of other Salmonella proteins. Our demonstration that Salmonella-mediated cytoskeletal modifications are independent of inositolphosphates will focus future studies on elucidating alternate pathogenic consequences of InsP5 metabolism, including ion channel conductance and apoptosis.

The HITRAN molecular spectroscopic database: Edition of 2000 including updates through 2001

This paper describes the status circa 2001, of the HITRAN compilation that comprises the public edition available through 2001. The HITRAN compilation consists of several components useful for radiative transfer calculation codes: high-resolution spectroscopic parameters of molecules in the gas phase, absorption cross-sections for molecules with very dense spectral features, aerosol refractive indices, ultraviolet line-by-line parameters and absorption cross-sections, and associated database management software. The line-by-line portion of the database contains spectroscopic parameters for 38 molecules and their isotopologues and isotopomers suitable for calculating atmospheric transmission and radiance properties. Many more molecular species are presented in the infrared cross-section data than in the previous edition, especially the chlorofluorocarbons and their replacement gases. There is now sufficient representation so that quasi-quantitative simulations can be obtained with the standard radiance codes. In addition to the description and justification of new or modified data that have been incorporated since the last edition of HITRAN (1996), future modifications are indicated for cases considered to have a significant impact on remote-sensing experiments. © 2003 Elsevier Ltd. All rights reserved.

A 3/2 algorithm for two-machine open shop with route-dependent processing times

This paper considers the problem of minimizing the schedule length of a two-machine shop in which not only can a job be assigned any of the two possible routes, but also the processing times depend on the chosen route. This problem is known to be NP-hard. We describe a simple approximation algorithm that guarantees a worst-case performance ratio of 2. We also present some modifications to this algorithm that improve its performance and guarantee a worst-case performance ratio of 3=2.

Analysing the evacuation procedures employed on a Thames passenger boat using the maritimeEXODUS evacuation model

On the 19 June 2001, a Thames passenger/tour boat underwent several evacuation trials. This work was conducted in order to collect data for the validation of marine-based computer models. The trials involved 111 participants who were distributed throughout the vessel. The boat had two decks and two points of exit from the lower deck placed on either side of the craft, forward and aft. The boat had a twin set of staircases towards the rear of the craft, just forward of the rear exits. maritimeEXODUS was used to simulate the full-scale evacuation trials conducted. The simulation times generated were compared against the original results and categorised according to the exit point availability. The predictions closely approximate the original results, differing by an average of 6.6% across the comparisons, with numerous qualitative similarities between the predictions and experimental results. The maritimeEXODUS evacuation model was then used to examine the evacuation procedure currently employed on the vessel. This was found to have potential to produce long evacuation times. maritimeEXODUS was used to suggest modifications to the mustering procedures. These theoretical results suggest that it is possible to significantly reduce evacuation times.

Implementation of a proactive load sharing scheme

This paper presents a proactive approach to load sharing and describes the architecture of a scheme, Concert, based on this approach. A proactive approach is characterized by a shift of emphasis from reacting to load imbalance to avoiding its occurrence. In contrast, in a reactive load sharing scheme, activity is triggered when a processing node is either overloaded or underloaded. The main drawback of this approach is that a load imbalance is allowed to develop before costly corrective action is taken. Concert is a load sharing scheme for loosely-coupled distributed systems. Under this scheme, load and task behaviour information is collected and cached in advance of when it is needed. Concert uses Linux as a platform for development. Implemented partially in kernel space and partially in user space, it achieves transparency to users and applications whilst keeping the extent of kernel modifications to a minimum. Non-preemptive task transfers are used exclusively, motivated by lower complexity, lower overheads and faster transfers. The goal is to minimize the average response-time of tasks. Concert is compared with other schemes by considering the level of transparency it provides with respect to users, tasks and the underlying operating system.

Predicting toxic species concentrations in a full-scale vitiated fire

A toxicity model on dividing the computational domain into two parts, a control region (CR) and a transport region (TR), for species calculation was recently developed. The model can be incorporated with either the heat source approach or the eddy dissipation model (EDM). The work described in this paper is a further application of the toxicity model with modifications of the EDM for vitiated fires. In the modified EDM, chemical reaction only occurs within the CR. This is consistent with the approach used in the species concentration calculations within the toxicity model in which yields of combustion products only change within the CR. A vitiated large room-corridor fire, in which the carbon monoxide (CM) concentrations are very high and the temperatures are relatively low at locations distant from the original fire source, is simulated using the modified EDM coupled with the toxicity model. Compared with the EDM, the modified EDM provide significant improvements in the predictions of temperatures at remote locations. Predictions of species concentrations at various locations follow the measured trends. Good agreements between the measured and predicted species concentrations are obtained at the vitiated fire stage.