Integrating copy number polymorphisms into array CGH analysis using a robust HMM

MOTIVATION: Array comparative genomic hybridization (aCGH) is a pervasive technique used to identify chromosomal aberrations in human diseases, including cancer. Aberrations are defined as regions of increased or decreased DNA copy number, relative to a normal sample. Accurately identifying the locations of these aberrations has many important medical applications. Unfortunately, the observed copy number changes are often corrupted by various sources of noise, making the boundaries hard to detect. One popular current technique uses hidden Markov models (HMMs) to divide the signal into regions of constant copy number called segments; a subsequent classification phase labels each segment as a gain, a loss or neutral. Unfortunately, standard HMMs are sensitive to outliers, causing over-segmentation, where segments erroneously span very short regions. RESULTS: We propose a simple modification that makes the HMM robust to such outliers. More importantly, this modification allows us to exploit prior knowledge about the likely location of "outliers", which are often due to copy number polymorphisms (CNPs). By "explaining away" these outliers with prior knowledge about the locations of CNPs, we can focus attention on the more clinically relevant aberrated regions. We show significant improvements over the current state of the art technique (DNAcopy with MergeLevels) on previously published data from mantle cell lymphoma cell lines, and on published benchmark synthetic data augmented with outliers. AVAILABILITY: Source code written in Matlab is available from http://www.cs.ubc.ca/~sshah/acgh.     

Biosurfactants and oil bioremediation

Oil pollution is an environmental problem of increasing importance. Hydrocarbon-degrading microorganisms, adapted to grow and thrive in oil-containing environments, have an important role in the biological treatment of this pollution. One of the limiting factors in this process is the bioavailability of many fractions of the oil. The hydrocarbon-degrading microorganisms produce biosurfactants of diverse chemical nature and molecular size. These surface-active materials increase the surface area of hydrophobic water-insoluble substrates and increase their bioavailability, thereby enhancing the growth of bacteria and the rate of bioremediation.     

Effect of an extract of nereis virens on mammalian spermatozoa

Epididymal sperm of the mouse, rat, and guinea pig and ejaculated sperm of rabbits are cleaved at the head-tail junction by an extract of Nereis virens. Annelids are extracted with water and the extract is purified by ion exchange chromatography. Electron microscopy shows that the extract acts on the filaments connecting the capitulum of the tail with the basal plate lining the nuclear envelope. Following detachment, the basal plate remains with the head. The extract contains proteases as indicated by hydrolysis of tosyl arginine methyl ester (TAME), benzoyl arginine ethyl ester (BAEE), and Azocoll, a general protease substrate. The hydrolysis of TAME is inhibited by tosyl lysine chloromethyl ketone (TLCK), a trypsin inhibitor, but TLCK does not prevent head-tail separation by the Nereis extract. Similarly tosyl phenylalanine chloromethyl ketone (TPCK), a chymotrypsin inhibitor, and phosphoramidon and leucyltryptophan, both thermolysin and acrolysin inhibitors — singly or in combination — do not prevent hydrolysis of Azocoll. Head-tail separation activity of the extract was inhibited by dithiothreitol, which reduces disulfide bonds, and phenylmethyl sulfonyl fluoride, an inhibitor of serine proteases. These results indicate that the extract is a mixture of proteases, one being a serine protease similar to trypsin. Digestion of the connecting filaments with the pure proteases, trypsin and Staphylococcus aureus V8 protease, has yielded the following information on the proteins of the filaments. The accessibility of arginine and/or lysine peptide bonds to enzyme action is highest in rat sperm filaments, whereas those in the filaments of mouse, rabbit, and guinea pig sperm are less accessible than in the rat. Another possibility is that the total content of arginine and/or lysine varies between the species. The most dramatic difference is the enzymatic action on glutamyl peptide bonds of the filaments, the order being: mouse 〉 rat 〉 rabbit, with guinea pig sperm filaments completely resistant over the time course of the experiment.     

Cysteine residues in the nucleotide binding domains regulate the conductance state of CFTR channels

Gating of cystic fibrosis transmembrane conductance regulator (CFTR) channels requires intermolecular or interdomain interactions, but the exact nature and physiological significance of those interactions remains uncertain. Subconductance states of the channel may result from alterations in interactions among domains, and studying mutant channels enriched for a single conductance type may elucidate those interactions. Analysis of CFTR channels in inside-out patches revealed that mutation of cysteine residues in NBD1 and NBD2 affects the frequency of channel opening to the full-size versus a 3-pS subconductance. Mutating cysteines in NBD1 resulted in channels that open almost exclusively to the 3-pS subconductance, while mutations of cysteines in NBD2 decreased the frequency of subconductance openings. Wild-type channels open to both size conductances and make fast transitions between them within a single open burst. Full-size and subconductance openings of both mutant and wild-type channels are similarly activated by ATP and phosphorylation. However, the different size conductances open very differently in the presence of a nonhydrolyzable ATP analog, with subconductance openings significantly shortened by ATPgammaS, while full-size channels are locked open. In wild-type channels, reducing conditions increase the frequency and decrease the open time of subconductance channels, while oxidizing conditions decrease the frequency of subconductance openings. In contrast, in the cysteine mutants studied, altering redox potential has little effect on gating of the subconductance.     

Identification of the sex‐determining region in flathead grey mullet (Mugil cephalus)

Elucidation of the sex‐determination mechanism in flathead grey mullet (Mugil cephalus) is required to exploit its economic potential by production of genetically determined monosex populations and application of hormonal treatment to parents rather than to the marketed progeny. Our objective was to construct a first‐generation linkage map of the M. cephalus in order to identify the sex‐determining region and sex‐determination system. Deep‐sequencing data of a single male was assembled and aligned to the genome of Nile tilapia (Oreochromis niloticus). A total 245 M. cephalus microsatellite markers were designed, spanning the syntenic tilapia genome assembly at intervals of 10 Mb. In the mapping family of full‐sib progeny, 156 segregating markers were used to construct a first‐generation linkage map of 24 linkage groups (LGs), corresponding to the number of chromosomes. The linkage map spanned approximately 1200 cM with an average inter‐marker distance of 10.6 cM. Markers segregating on LG9 in two independent mapping families showed nearly complete concordance with gender (R² = 0.95). The sex determining locus was fine mapped within an interval of 8.6 cM on LG9. The sex of offspring was determined only by the alleles transmitted from the father, thus indicating an XY sex‐determination system.     

ExPASy: The proteomics server for in-depth protein knowledge and analysis

The ExPASy (the Expert Protein Analysis System) World Wide Web server (http://www.expasy.org), is provided as a service to the life science community by a multidisciplinary team at the Swiss Institute of Bioinformatics (SIB). It provides access to a variety of databases and analytical tools dedicated to proteins and proteomics. ExPASy databases include SWISS-PROT and TrEMBL, SWISS-2DPAGE, PROSITE, ENZYME and the SWISS-MODEL repository. Analysis tools are available for specific tasks relevant to proteomics, similarity searches, pattern and profile searches, post-translational modification prediction, topology prediction, primary, secondary and tertiary structure analysis and sequence alignment. These databases and tools are tightly interlinked: a special emphasis is placed on integration of database entries with related resources developed at the SIB and elsewhere, and the proteomics tools have been designed to read the annotations in SWISS-PROT in order to enhance their predictions. ExPASy started to operate in 1993, as the first WWW server in the field of life sciences. In addition to the main site in Switzerland, seven mirror sites in different continents currently serve the user community.     

Keratinocyte growth factor

Keratinocyte growth factor (KGF) is a member of the heparin-binding fibroblast growth factor family (FGF-7) with a distinctive pattern of target-cell specificity. Studies performed in cell culture suggested that KGF was mitogenically active only on epithelial cells, albeit from a variety of tissues. In contrast, KGF was produced solely by cells of mesenchymal origin, leading to the hypothesis that it might function as a paracrine mediator of mesenchymal-epithelial communication. Biochemical analysis and molecular cloning established that the KGF receptor (KGFR) was a tyrosine kinase isoform encoded by the fgfr-2 gene. Many detailed investigations of KGF and KGFR expression in whole tissue and cell lines largely substantiated the pattern initially perceived in vitro of mesenchymal and epithelial distribution, respectively. Moreover, functional assays in organ culture and in vivo and studies of KGF regulation by sex sterorid hormones reinforced the idea that KGF acts predominantly on epithelial cells to elicit a variety of responses including proliferation, migration and morphogenesis.     

Ventilatory response to helium-oxygen breathing during exercise: effect of airway anesthesia

Krishnan, Bharath S., Ron E. Clemens, Trevor A. Zintel,Martin J. Stockwell, and Charles G. Gallagher. Ventilatory response to helium-oxygen breathing during exercise: effect of airwayanesthesia. J. Appl. Physiol. 83(1):82-88, 1997.The substitution of a normoxic helium mixture(HeO₂) for room air (Air) during exercise results in a sustained hyperventilation, which is present evenin the first breath. We hypothesized that this response is dependent onintact airway afferents; if so, airway anesthesia (Anesthesia) shouldaffect this response. Anesthesia was administered to the upper airwaysby topical application and to lower central airways by aerosolinhalation and was confirmed to be effective for over 15 min. Subjectsperformed constant work-rate exercise (CWE) at 69 ± 2 (SE) % maximal work rate on a cycle ergometer on three separate days: twiceafter saline inhalation (days 1 and3) and once after Anesthesia(day 2). CWE commenced after a briefwarm-up, with subjects breathing Air for the first 5 min (Air-1),HeO₂ for the next 3 min, and Airagain until the end of CWE (Air-2). The resistance of the breathingcircuit was matched for Air andHeO₂. BreathingHeO₂ resulted in a small butsignificant increase in minute ventilation(I) anddecrease in alveolar PCO₂ in both theSaline (average of 2 saline tests; not significant) and Anesthesiatests. Although Anesthesia had no effect on the sustainedhyperventilatory response to HeO₂breathing, theI transientswithin the first six breaths ofHeO₂ were significantly attenuatedwith Anesthesia. We conclude that theI response to HeO₂ is not simply due to areduction in external tubing resistance and that, in humans, airwayafferents mediate the transient but not the sustained hyperventilatoryresponse to HeO₂ breathing duringexercise.

     

Application of a predictive approach to estimate exposure to non‐smoking urban sub‐populations to background levels of Benzene in Ontario

An indirect approach was used to estimate exposure to background levels of atmospheric benzene for a selection of Ontario non‐smoking urban sub‐populations. Activity codes obtained from nationally representative time‐budget surveys were allocated to five general microenvironments (home, work (or school), outdoors, commuting, and other indoors) and further combined with inhalation rates corresponding to specific levels of physical activity in order to develop average activity patterns for six sub‐populations believed to be differently exposed to atmospheric benzene in urban areas (children, teenagers, office workers, outdoor workers, transit workers, and seniors). These activity patterns were then combined with representative concentrations of benzene measured in the selected microenvironments in Ontario in order to evaluate exposure. Two metrics were used for this purpose, integrated exposure and potential average daily dose (intake). Potential lifetime average daily doses were also estimated for three composite sub‐groups representing average office, outdoor, and transit workers. A probabilistic approach using a Monte‐Carlo sampling procedure was used in order to estimate possible ranges of exposures experienced by the various sub‐populations. Results obtained suggested that the highest daily integrated exposure (mean: 131 μg‐hrs/m³) was associated with the average transit worker while comparable levels were estimated for the other sub‐populations investigated (mean: 77–86 μg‐hrs/m³). These levels corresponded to 24‐hours time‐weighted average (TWA)‐equivalent concentrations of 5.5 μg/m³ and 3.5 μg/m³, respectively. Statistical distributions of integrated exposures obtained with the probabilistic approach indicated levels as high as 343 μg‐hrs/m³ (97.5th percentile) in the case of the average transit worker, corresponding to TWA‐equivalents in excess of 15 μg/m³. When levels of physical activities and inhalation rates were integrated in the calculation of exposure, the highest potential average daily dose was found to be associated with the average child (mean: 3.1 μg/kg‐day; 97.5th percentile: 6.0 μg/kg‐day) whereas comparable amounts were estimated for teenager and transit workers (mean: 2.1 μg/kg‐day; 97.5th: 4.1 and 6.9 μg/kg‐day, respectively). Indoor microenvironments (home, office/school, other indoors) were identified as the major contributors to total exposure and intake of benzene (≥50%) although their relative importance varied depending on the exposure metric utilized. Potential lifetime average daily doses estimated for transit workers varied from 2.1 (mean) to 5.4 (97.5th) μg/kg‐day. This was slightly higher than those estimated for the average office and outdoor workers (mean: 1.5–1.7 μg/kg‐day). These projections suggest that average non‐smoking children and teenagers are the most exposed sub‐populations among those investigated to background levels of atmospheric benzene as a result of their daily activities. However, these projections must be interpreted with caution in view of uncertainties associated with some of the assumptions adopted, the limited database used to document benzene levels, and as a result of unaccounted sources of emissions which, under certain circumstances, can significantly modify these conclusions.