Evolutionary characters, phenotypes and ontologies: curating data from the systematic biology literature.

BACKGROUND: The wealth of phenotypic descriptions documented in the published articles, monographs, and dissertations of phylogenetic systematics is traditionally reported in a free-text format, and it is therefore largely inaccessible for linkage to biological databases for genetics, development, and phenotypes, and difficult to manage for large-scale integrative work. The Phenoscape project aims to represent these complex and detailed descriptions with rich and formal semantics that are amenable to computation and integration with phenotype data from other fields of biology. This entails reconceptualizing the traditional free-text characters into the computable Entity-Quality (EQ) formalism using ontologies. METHODOLOGY/PRINCIPAL FINDINGS: We used ontologies and the EQ formalism to curate a collection of 47 phylogenetic studies on ostariophysan fishes (including catfishes, characins, minnows, knifefishes) and their relatives with the goal of integrating these complex phenotype descriptions with information from an existing model organism database (zebrafish, http://zfin.org). We developed a curation workflow for the collection of character, taxonomic and specimen data from these publications. A total of 4,617 phenotypic characters (10,512 states) for 3,449 taxa, primarily species, were curated into EQ formalism (for a total of 12,861 EQ statements) using anatomical and taxonomic terms from teleost-specific ontologies (Teleost Anatomy Ontology and Teleost Taxonomy Ontology) in combination with terms from a quality ontology (Phenotype and Trait Ontology). Standards and guidelines for consistently and accurately representing phenotypes were developed in response to the challenges that were evident from two annotation experiments and from feedback from curators. CONCLUSIONS/SIGNIFICANCE: The challenges we encountered and many of the curation standards and methods for improving consistency that we developed are generally applicable to any effort to represent phenotypes using ontologies. This is because an ontological representation of the detailed variations in phenotype, whether between mutant or wildtype, among individual humans, or across the diversity of species, requires a process by which a precise combination of terms from domain ontologies are selected and organized according to logical relations. The efficiencies that we have developed in this process will be useful for any attempt to annotate complex phenotypic descriptions using ontologies. We also discuss some ramifications of EQ representation for the domain of systematics.

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.

Oldest known cranium of a juvenile New World monkey (Early Miocene, Patagonia, Argentina): implications for the taxonomy and the molar eruption pattern of early platyrrhines.

A juvenile cranium of Homunculus patagonicus Ameghino, 1891a from the late Early Miocene of Santa Cruz Province (Argentina) provides the first evidence of developing cranial anatomy for any fossil platyrrhine. The specimen preserves the rostral part of the cranium with deciduous and permanent alveoli and teeth. The dental eruption sequence in the new specimen and a reassessment of eruption patterns in living and fossil platyrrhines suggest that the ancestral platyrrhine pattern of tooth replacement was for the permanent incisors to erupt before M(1), not an accelerated molar eruption (before the incisors) as recently proposed. Two genera and species of Santacrucian monkeys are now generally recognized: H. patagonicus Ameghino, 1891a and Killikaike blakei Tejedor et al., 2006. Taxonomic allocation of Santacrucian monkeys to these species encounters two obstacles: 1) the (now lost) holotype and a recently proposed neotype of H. patagonicus are mandibles from different localities and different geologic members of the Santa Cruz Formation, separated by approximately 0.7 million years, whereas the holotype of K. blakei is a rostral part of a cranium without a mandible; 2) no Santacrucian monkey with associated cranium and mandible has ever been found. Bearing in mind these uncertainties, our examination of the new specimen as well as other cranial specimens of Santacrucian monkeys establishes the overall dental and cranial similarity between the holotype of Killikaike blakei, adult cranial material previously referred to H. patagonicus, and the new juvenile specimen. This leads us to conclude that Killikaike blakei is a junior subjective synonym of H. patagonicus.

Modeling and experiments on an isothermal fatigue test for solder joints

This paper investigates an isothermal fatigue test for solder joints developed at the NPL. The test specimen is a lap joint between two copper arms. During the test the displacement at the ends of the copper are controlled and the force measured. The modeling results in the paper show that the displacement across the solder joint is not equal to the displacement applied at the end of the specimen. This is due to deformation within the copper arms. A method is described to compensate for this difference. The strain distribution in the solder was determined by finite element analysis and compared to the distribution generated by a theoretical 'ideal' test which generates an almost pure shear mode in the solder. By using a damage-based constitutive law the shape of the crack generated in the specimen has been predicted for both the actual test and the ideal pure shear test. Results from the simulations are also compared with experimental data using SnAgCu solder.

The study of flow and temperature fields in conducting droplets suspended in a DC/AC combination field

Electromagnetic Levitation (EML) is a valuable method for measuring the thermo-physical properties of metals - surface tensions, viscosity, thermal/electrical conductivity, specific heat, hemispherical emissivity, etc. – beyond their melting temperature. In EML, a small amount of the test specimen is melted by Joule heating in a suspended AC coil. Once in liquid state, a small perturbation causes the liquid envelope to oscillate and the frequency of oscillation is then used to compute its surface tension by the well know Rayleigh formula. Similarly, the rate at which the oscillation is dampened relates to the viscosity. To measure thermal conductivity, a sinusoidally varying laser source may be used to heat the polar axis of the droplet and the temperature response measured at the polar opposite – the resulting phase shift yields thermal conductivity. All these theoretical methods assume that convective effects due to flow within the droplet are negligible compared to conduction, and similarly that the flow conditions are laminar; a situation that can only be realised under microgravity conditions. Hence the EML experiment is the method favoured for Spacelab experiments (viz. TEMPUS). Under terrestrial conditions, the full gravity force has to be countered by a much larger induced magnetic field. The magnetic field generates strong flow within the droplet, which for droplets of practical size becomes irrotational and turbulent. At the same time the droplet oscillation envelope is no longer ellipsoidal. Both these conditions invalidate simple theoretical models and prevent widespread EML use in terrestrial laboratories. The authors have shown in earlier publications that it is possible to suppress most of the turbulent convection generated in the droplet skin layer, through use of a static magnetic field. Using a pseudo-spectral discretisation method it is possible compute very accurately the dynamic variation in the suspended fluid envelope and simultaneously compute the time-varying electromagnetic, flow and thermal fields. The use of a DC field as a dampening agent was also demonstrated in cold crucible melting, where suppression of turbulence was achieved in a much larger liquid metal volume and led to increased superheat in the melt and reduction of heat losses to the water-cooled walls. In this paper, the authors describe the pseudo-spectral technique as applied to EML to compute the combined effects of AC and DC fields, accounting for all the flow-induced forces acting on the liquid volume (Lorentz, Maragoni, surface tension, gravity) and show example simulations.

The use of acoustic emission monitoring to rank paper materials with respect to their fracture toughness

In this study, a simplified Acoustic Emission (AE) equipment, in essence an AE signal conditioner and a USB (Universal Serial Bus) data acquisition system, is used to study what happens in paper structures during mechanical loading. By the use of such equipment, some parameters that can be extracted are e.g. the stress and strain at onset of AE, the stress and strain at the onset of rapid AE defined as some numerical factor (larger then one) times the initial emission rate, the emission rate at the first stage of loading and the stress and strain at final failure i.e. when the specimen loses its load carrying ability.In this study however, the interest is focused on one particular parameter i.e. the elastic strain energy density W c at onset of AE. This is a parameter with a clear physical meaning and in this study, the correlation between this parameter and a fracture toughness measure, is investigated.The conclusion is that when nine different paper materials (with a large span regarding properties) are considered, there is a correlation (however not linear) between these two parameters.

Comparisons between “sand blast” and “centripetal effect accelerator” type erosion testers

There are two major types of erosion testing devices that are used throughout the world for quantifying particle impact erosion against a solid surface. The first of these uses pressurised air to accelerate abrasive particles through a nozzle so that they impinge upon a target specimen. The second adopts a rotating disc to accelerate abrasive particles using the centripetal effect so that they impinge upon a series of targets arranged around the periphery of the disc. This paper reports the findings of a collaborative project that was designed to compare the performance and results obtained from a rig of each of the two types mentioned above. The sand blast type rig was provided by The Department of Powder Science Technology (POSTEC) at The Telemark Technological Research and Development Centre (TEL-TEK), Porsgrunn, Norway while the centripetal effect accelerator was provided by The Wolfson Centre for Bulk Solids Handling Technology, University of Greenwich, London, UK. The test programme included tests against a wide range of materials that are commonly used in pneumatic handling facilities. (Pneumatic handling is a means of conveying and transporting powders and granular solid materials in bulk in industrial process plant, through pipelines using a gas as the carrier medium.) Olivine sand was used as the abrasive and it was projected against the test specimens at velocities and concentrations commensurate with those seen in pneumatic conveyors. In all instances the materials used in the test programme were taken from the same batch so that scatter of experimental results due to specimen variation was minimised. The paper contains a series of recommendations for erosion testing equipment. A discussion based on the results and their applicability to the prediction of wear in pneumatic conveyors concludes the paper.

Notes and memoranda

On a specimen of Leptocephalus morrisii. New record of Cuthona? aurantiaca off Plymouth (UK)

A unusual pinnotherid zoea attributed to Afropinnotheres monodi Manning, 1993 (Decapoda: Brachyura: Pinnotheridea) from the Selvagens Islands (Eastern Atlantic Ocean)

A pinnotherid zoea taken in a plankton sample from the Selvagens Islands during the TFMCBM/SELVAGENS’2000 Cruise, organized by the Natural Sciences Museum of Tenerife (Canary Islands), differs from previously described zoeas of the family. The specimen has dorsal and lateral spines but no rostrum, a combination of characters not previously described from pinnotherid zoeas. The lateral spines originate behind the upper part of the eye, unlike all previously described Pinnotheridae zoeas where these spines, when present, originate near the ventrolateral margin of the carapace. The specimen is attributed to Afropinnotheres monodi, the only pinnotherid species known from the area with undescribed larvae.

Granular gland transcriptomes in stimulated amphibian skin secretions.

Amphibian defensive skin secretions are complex, species-specific cocktails of biologically active molecules, including many uncharacterized peptides. The study of such secretions for novel peptide discovery is time-limited, as amphibians are in rapid global decline. While secretion proteome analysis is non-lethal, transcriptome analysis has until now required killing of specimens prior to skin dissection for cDNA library construction. Here we present the discovery that polyadenylated mRNAs encoding dermal granular gland peptides are present in defensive skin secretions, stabilized by endogenous nucleic acid-binding amphipathic peptides. Thus parallel secretory proteome and transcriptome analyses can be performed without killing the specimen in this model amphibian system--a finding that has important implications in conservation of biodiversity within this threatened vertebrate taxon and whose mechanistics may have broader implications in biomolecular science.

Partial structure of the phylloxin gene from the giant monkey frog, Phyllomedusa bicolor: parallel cloning of precursor cDNA and genomic DNA from lyophilized skin secretion

Phylloxin is a novel prototype antimicrobial peptide from the skin of Phyllomedusa bicolor. Here, we describe parallel identification and sequencing of phylloxin precursor transcript (mRNA) and partial gene structure (genomic DNA) from the same sample of lyophilized skin secretion using our recently-described cloning technique. The open-reading frame of the phylloxin precursor was identical in nucleotide sequence to that previously reported and alignment with the nucleotide sequence derived from genomic DNA indicated the presence of a 175 bp intron located in a near identical position to that found in the dermaseptins. The highly-conserved structural organization of skin secretion peptide genes in P. bicolor can thus be extended to include that encoding phylloxin (plx). These data further reinforce our assertion that application of the described methodology can provide robust genomic/transcriptomic/peptidomic data without the need for specimen sacrifice.

Dermatoxin and phylloxin from the waxy monkey frog, Phyllomedusa sauvagei: cloning of precursor cDNAs and structural characterization from lyophilized skin secretion

Amphibian skin is a morphologically, biochemically and physiologically complex organ that performs the wide range of functions necessary for amphibian survival. Here we describe the primary structures of representatives of two novel classes of amphibian skin antimicrobials, dermatoxin and phylloxin, from the skin secretion of Phyllomedusa sauvagei, deduced from their respective precursor encoding cDNAs cloned from a lyophilized skin secretion library. A degenerate primer, designed to a highly conserved domain in the 5'-untranslated region of analogous peptide precursor cDNAs from Phyllomedusa bicolor, was employed in a 3'-RACE reaction. Peptides with molecular masses coincident with precursor-deduced mature toxin peptides were identified in LC/MS fractions of skin secretion and primary structures were confirmed by MS/MS fragmentation. This integrated experimental approach can thus rapidly expedite the primary structural characterization of amphibian skin peptides in a manner that circumvents specimen sacrifice whilst preserving robustness of scientific data.

Isolation of scorpion (Androctonus amoreuxi) alpha-neurotoxins and parallel cloning of their respective cDNAs from a single sample of venom

The venoms of buthid scorpions are known to contain basic, single-chain protein toxins (alpha toxins) consisting of 60–70 amino acid residues that are tightly folded by four disulfide bridges. Here we describe isolation and sequencing of three novel putative alpha toxins (AamH1-3) from the venom of the North African scorpion, Androctonus amoreuxi, and subsequent cloning of their precursor cDNAs from the same sample of venom. This experimental approach can expedite functional genomic analyses of the protein toxins from this group of venomous animals and does not require specimen sacrifice for cloning of protein toxin precursor cDNAs.

Unmasking venom gland transcriptomes in reptile venoms

While structural studies of reptile venom toxins can be achieved using lyophilized venom samples, until now the cloning of precursor cDNAs required sacrifice of the specimen for dissection of the venom glands. Here we describe a simple and rapid technique that unmasks venom protein mRNAs present in lyophilized venom samples. To illustrate the technique we have RT-PCR-amplified a range of venom protein transcripts from cDNA libraries derived from the venoms of a hemotoxic snake, the Chinese copperhead (Deinagkistrodon acutus), a neurotoxic snake, the black mamba (Dendroaspis polylepis), and a venomous lizard, the Gila monster (Heloderma suspectum). These include a metalloproteinase and phospholipase A2 from D. acutus, a potassium channel blocker, dendrotoxin K, from D. polylepis, and exendin-4 from H. suspectum. These findings imply that the apparent absence and/or lability of mRNA in complex biological matrices is not always real and paves the way for accelerated acquisition of molecular genetic data on venom toxins for scientific and potential therapeutic purposes without sacrifice of endangered herpetofauna.

Detection of Staphylococcus aureus by 16S rRNA directed in situ hybridisation in a patient with a brain abscess caused by small colony variants.

Kipp F, Ziebuhr W, Becker K, Krimmer V, Höbeta N, Peters G, Von Eiff C. Institute of Medical Microbiology, Hospital and Clinics, University of Münster, Germany. A 45 year old man was admitted to hospital with a right sided facial paralysis and three month history of seizures. Computed tomography showed a left temporal mass including both intracerebral and extracerebral structures. Ten years earlier the patient had undergone a neurosurgical intervention in the same anatomical region to treat a subarachnoid haemorrhage. In tissue samples and pus obtained during neurosurgery, Staphylococcus aureus was detected by a 16S rRNA-directed in situ hybridisation technique. Following long term cultivation, small colony variants (SCV) of methicillin resistant S aureus were identified. The patient was treated successfully with a combination of vancomycin and rifampin followed by prolonged treatment with teicoplanin, with no sign of infection on follow up nine months after discharge. This is the first report in which S aureus SCV have been identified as causative organisms in a patient with brain abscess and in which in situ hybridisation has been used to detect S aureus in a clinical specimen containing SCV. Antimicrobial agents such as rifampin which have intracellular activity should be included in treatment of infections caused by S aureus SCV.