A Faster Circular Binary Segmentation Algorithm for the Analysis of Array CGH Data

Motivation: Array CGH technologies enable the simultaneous measurement of DNA copy number for thousands of sites on a genome. We developed the circular binary segmentation (CBS) algorithm to divide the genome into regions of equal copy number (Olshen {\it et~al}, 2004). The algorithm tests for change-points using a maximal $t$-statistic with a permutation reference distribution to obtain the corresponding $p$-value. The number of computations required for the maximal test statistic is $O(N^2),$ where $N$ is the number of markers. This makes the full permutation approach computationally prohibitive for the newer arrays that contain tens of thousands markers and highlights the need for a faster. algorithm. Results: We present a hybrid approach to obtain the $p$-value of the test statistic in linear time. We also introduce a rule for stopping early when there is strong evidence for the presence of a change. We show through simulations that the hybrid approach provides a substantial gain in speed with only a negligible loss in accuracy and that the stopping rule further increases speed. We also present the analysis of array CGH data from a breast cancer cell line to show the impact of the new approaches on the analysis of real data. Availability: An R (R Development Core Team, 2006) version of the CBS algorithm has been implemented in the ``DNAcopy'' package of the Bioconductor project (Gentleman {\it et~al}, 2004). The proposed hybrid method for the $p$-value is available in version 1.2.1 or higher and the stopping rule for declaring a change early is available in version 1.5.1 or higher.

A Data-Analytic Strategy for Protein-Biomarker Discovery: Profiling of High-Dimensional Proteomic Data for Cancer Detection

With recent advances in mass spectrometry techniques, it is now possible to investigate proteins over a wide range of molecular weights in small biological specimens. This advance has generated data-analytic challenges in proteomics, similar to those created by microarray technologies in genetics, namely, discovery of "signature" protein profiles specific to each pathologic state (e.g., normal vs. cancer) or differential profiles between experimental conditions (e.g., treated by a drug of interest vs. untreated) from high-dimensional data. We propose a data analytic strategy for discovering protein biomarkers based on such high-dimensional mass-spectrometry data. A real biomarker-discovery project on prostate cancer is taken as a concrete example throughout the paper: the project aims to identify proteins in serum that distinguish cancer, benign hyperplasia, and normal states of prostate using the Surface Enhanced Laser Desorption/Ionization (SELDI) technology, a recently developed mass spectrometry technique. Our data analytic strategy takes properties of the SELDI mass-spectrometer into account: the SELDI output of a specimen contains about 48,000 (x, y) points where x is the protein mass divided by the number of charges introduced by ionization and y is the protein intensity of the corresponding mass per charge value, x, in that specimen. Given high coefficients of variation and other characteristics of protein intensity measures (y values), we reduce the measures of protein intensities to a set of binary variables that indicate peaks in the y-axis direction in the nearest neighborhoods of each mass per charge point in the x-axis direction. We then account for a shifting (measurement error) problem of the x-axis in SELDI output. After these pre-analysis processing of data, we combine the binary predictors to generate classification rules for cancer, benign hyperplasia, and normal states of prostate. Our approach is to apply the boosting algorithm to select binary predictors and construct a summary classifier. We empirically evaluate sensitivity and specificity of the resulting summary classifiers with a test dataset that is independent from the training dataset used to construct the summary classifiers. The proposed method performed nearly perfectly in distinguishing cancer and benign hyperplasia from normal. In the classification of cancer vs. benign hyperplasia, however, an appreciable proportion of the benign specimens were classified incorrectly as cancer. We discuss practical issues associated with our proposed approach to the analysis of SELDI output and its application in cancer biomarker discovery.

Development of a genotyping microarray for Usher syndrome

BACKGROUND: Usher syndrome, a combination of retinitis pigmentosa (RP) and sensorineural hearing loss with or without vestibular dysfunction, displays a high degree of clinical and genetic heterogeneity. Three clinical subtypes can be distinguished, based on the age of onset and severity of the hearing impairment, and the presence or absence of vestibular abnormalities. Thus far, eight genes have been implicated in the syndrome, together comprising 347 protein-coding exons. METHODS: To improve DNA diagnostics for patients with Usher syndrome, we developed a genotyping microarray based on the arrayed primer extension (APEX) method. Allele-specific oligonucleotides corresponding to all 298 Usher syndrome-associated sequence variants known to date, 76 of which are novel, were arrayed. RESULTS: Approximately half of these variants were validated using original patient DNAs, which yielded an accuracy of >98%. The efficiency of the Usher genotyping microarray was tested using DNAs from 370 unrelated European and American patients with Usher syndrome. Sequence variants were identified in 64/140 (46%) patients with Usher syndrome type I, 45/189 (24%) patients with Usher syndrome type II, 6/21 (29%) patients with Usher syndrome type III and 6/20 (30%) patients with atypical Usher syndrome. The chip also identified two novel sequence variants, c.400C>T (p.R134X) in PCDH15 and c.1606T>C (p.C536S) in USH2A. CONCLUSION: The Usher genotyping microarray is a versatile and affordable screening tool for Usher syndrome. Its efficiency will improve with the addition of novel sequence variants with minimal extra costs, making it a very useful first-pass screening tool.

Error models for microarray intensities

We derive the additive-multiplicative error model for microarray intensities, and describe two applications. For the detection of differentially expressed genes, we obtain a statistic whose variance is approximately independent of the mean intensity. For the post hoc calibration (normalization) of data with respect to experimental factors, we describe a method for parameter estimation.

Bayesian Hidden Markov Modeling of Array CGH Data

Genomic alterations have been linked to the development and progression of cancer. The technique of Comparative Genomic Hybridization (CGH) yields data consisting of fluorescence intensity ratios of test and reference DNA samples. The intensity ratios provide information about the number of copies in DNA. Practical issues such as the contamination of tumor cells in tissue specimens and normalization errors necessitate the use of statistics for learning about the genomic alterations from array-CGH data. As increasing amounts of array CGH data become available, there is a growing need for automated algorithms for characterizing genomic profiles. Specifically, there is a need for algorithms that can identify gains and losses in the number of copies based on statistical considerations, rather than merely detect trends in the data. We adopt a Bayesian approach, relying on the hidden Markov model to account for the inherent dependence in the intensity ratios. Posterior inferences are made about gains and losses in copy number. Localized amplifications (associated with oncogene mutations) and deletions (associated with mutations of tumor suppressors) are identified using posterior probabilities. Global trends such as extended regions of altered copy number are detected. Since the posterior distribution is analytically intractable, we implement a Metropolis-within-Gibbs algorithm for efficient simulation-based inference. Publicly available data on pancreatic adenocarcinoma, glioblastoma multiforme and breast cancer are analyzed, and comparisons are made with some widely-used algorithms to illustrate the reliability and success of the technique.

Pathotyping Escherichia coli by using miniaturized DNA microarrays

The detection of virulence determinants harbored by pathogenic Escherichia coli is important for establishing the pathotype responsible for infection. A sensitive and specific miniaturized virulence microarray containing 60 oligonucleotide probes was developed. It detected six E. coli pathotypes and will be suitable in the future for high-throughput use.

Clinical significance of cell cycle- and apoptosis-related markers in biliary tract cancer: a tissue microarray-based approach revealing a distinctive immunophenotype for intrahepatic and extrahepatic cholangiocarcinomas

Cholangiocarcinoma is the second most common malignant tumor of the liver. We analyzed, immunohistochemically, the significance of cell cycle- and apoptosis-related markers in 128 cholangiocarcinomas (42 intrahepatic, 70 extrahepatic, and 16 gallbladder carcinomas) combined in a tissue microarray. Follow-up was available for 57 patients (44.5%). In comparison with normal tissue (29 specimens), cholangiocarcinomas expressed significantly more frequently p53, bcl-2, bax, and COX-2 (P.05 <). Intrahepatic tumors were significantly more frequently bcl-2+ and p16+, whereas extrahepatic tumors were more often p53+ (P < .05). Loss of p16 expression was associated with reduced survival of patients. Our data show that p53, bcl-2, bax, and COX-2 have an important role in the pathogenesis of cholangiocarcinomas. The differential expression of p16, bcl-2, and p53 between intrahepatic and extrahepatic tumors demonstrates that there are location-related differences in the phenotype and the genetic profiles of these tumors. Moreover, p16 was identified as an important prognostic marker in cholangiocarcinomas.

Catapult-like release of mitochondrial DNA by eosinophils contributes to antibacterial defense

Although eosinophils are considered useful in defense mechanisms against parasites, their exact function in innate immunity remains unclear. The aim of this study is to better understand the role of eosinophils within the gastrointestinal immune system. We show here that lipopolysaccharide from Gram-negative bacteria activates interleukin-5 (IL-5)- or interferon-gamma-primed eosinophils to release mitochondrial DNA in a reactive oxygen species-dependent manner, but independent of eosinophil death. Notably, the process of DNA release occurs rapidly in a catapult-like manner--in less than one second. In the extracellular space, the mitochondrial DNA and the granule proteins form extracellular structures able to bind and kill bacteria both in vitro and under inflammatory conditions in vivo. Moreover, after cecal ligation and puncture, Il5-transgenic but not wild-type mice show intestinal eosinophil infiltration and extracellular DNA deposition in association with protection against microbial sepsis. These data suggest a previously undescribed mechanism of eosinophil-mediated innate immune responses that might be crucial for maintaining the intestinal barrier function after inflammation-associated epithelial cell damage, preventing the host from uncontrolled invasion of bacteria.

Ethanol-induced cytochrome P4502E1 causes carcinogenic etheno-DNA lesions in alcoholic liver disease

Oxidative stress is thought to play a major role in the pathogenesis of hepatocellular cancer (HCC), a frequent complication of alcoholic liver disease (ALD). However, the underlying mechanisms are poorly understood. In hepatocytes of ALD patients, we recently detected by immunohistochemistry significantly increased levels of carcinogenic etheno-DNA adducts that are formed by the reaction of the major lipid peroxidation product, 4-hydroxynonenal (4-HNE) with nucleobases. In the current study, we show that protein-bound 4-HNE and etheno-DNA adducts both strongly correlate with cytochrome P450 2E1 (CYP2E1) expression in patients with ALD (r = 0.9, P < 0.01). Increased levels of etheno-DNA adducts were also detected in the liver of alcohol-fed lean (Fa/?) and obese (fa/fa) Zucker rats. The number of nuclei in hepatocytes stained positively for etheno-DNA adducts correlated significantly with CYP2E1 expression (r = 0.6, P = 0.03). To further assess the role of CYP2E1 in the formation of etheno-DNA adducts, HepG2 cells stably transfected with human CYP2E1 were exposed to ethanol with or without chlormethiazole (CMZ), a specific CYP2E1 inhibitor. Ethanol increased etheno-DNA adducts in the nuclei of CYP2E1-transfected HepG2 cells in a concentration-dependent and time-dependent manner, but not in vector mock-transfected control cells. CMZ blocked the generation of etheno-DNA adducts by 70%-90% (P < 0.01). CONCLUSION: Our data support the assumption that ethanol-mediated induction of hepatic CYP2E1 leading inter alia to highly miscoding lipid peroxidation-derived DNA lesions may play a central role in hepatocarcinogenesis in patients with ALD.

Local gene transfer and expression following intramuscular administration of FGF-1 plasmid DNA in patients with critical limb ischemia

NV1FGF is an expression plasmid encoding sp.FGF-1(21-154) currently under investigation for therapeutic angiogenesis in clinical trials. NV1FGF plasmid distribution and transgene expression following intramuscular (IM) injection in patients is unknown. The study involved six patients with chronic critical limb ischemia (CLI) planned to undergo amputation. A total dose of 0.5, 2, or 4 mg NV1FGF was administered as eight IM injections (0.006, 0.25, or 0.5 mg per injection) 3-5 days before amputation. Injected sites (30 cm(3)) were divided into equally sized smaller pieces to assess spatial distribution of NV1FGF sequences (PCR), NV1FGF mRNA (reverse transcriptase-PCR), and fibroblast growth factor-1 (FGF-1)-expressing cells (immunohistochemistry). Data indicated gene expression at all doses. The distribution area was within 5-12 cm for NV1FGF sequences containing the expression cassette, up to 5 cm for NV1FGF mRNA, and up to 3 cm for FGF-1-expressing myofibers. All FGF receptors were detected indicating robust potential for bioactivity after NV1FGF gene transfer. Circulating levels of NV1FGF sequences were shown to decrease within days after injection. Data support demonstration of plasmid-mediated gene transfer and expression in muscles from patients with CLI. FGF-1 expression was shown to be limited to injection sites, which supports the concept of multiple-site injection for therapeutic use.

Primaquine-induced differential gene expression analysis in mice liver using DNA microarrays

Primaquine (PQ). a clinically important derivative of 8-aminoquinoline used against the hepatic stages (hypnozoites) of Plasmodium vivax and Plasmodium ova Ie. was studied to evaluate and compare between mRNA expression. and biochemical and histological parameters of hepatic stress in adult Swiss mice (Mus musculus). Following single oral dose of PQ (40 mglkg. bw). alanine aminotransferase (ALT) and aspartate aminotransferase (AST) along with hematoxylin and eosin stained liver sections did not show any signs of hepatic stress at 6. 12 and 24 h except for ALT activity at 6 h. However. analysis at RNA transcript level revealed consistent and significant deregulation (p<0.01 and twofold) of 16 probes corresponding to important cellular processes such as protein transportation. transcription regulation. intracellular signaling. protein synthesis, hematopoiesis, cell adhesion and cell proliferation. Pathway analysis identified large number of affected genes corresponding to 40 Gene Ontology terms having a z score greaibr than 2. These results indicate that PQ at high doses may affect gene expression in liver and may produce undesirable outcomes if consumed for longer durations.

Virulence typing of Escherichia coli using microarrays

We describe a microarray based broad-range screening technique for Escherichia coli virulence typing. Gene probes were amplified by PCR from a plasmid bank of characterised E. coli virulence genes and were spotted onto a glass slide to form an array of capture probes. Genomic DNA from E. coli strains which were to be tested for the presence of these virulence gene sequences was labelled with fluorescent cyanine dyes by random amplification and then hybridised against the array of probes. The hybridisation, washing and data analysis conditions were optimised for glass slides, and the applicability of the method for identifying the presence of the virulence genes was determined using reference strains and clinical isolates. It was found to be a sensitive screening method for detecting virulence genes, and a powerful tool for determining the pathotype of E. coli. It will be possible to expand and automate this microarray technique to make it suitable for rapid and reliable diagnostic screening of bacterial isolates.

ATG5 is induced by DNA-damaging agents and promotes mitotic catastrophe independent of autophagy

Anticancer drug therapy activates both molecular cell death and autophagy pathways. Here we show that even sublethal concentrations of DNA-damaging drugs, such as etoposide and cisplatin, induce the expression of autophagy-related protein 5 (ATG5), which is both necessary and sufficient for the subsequent induction of mitotic catastrophe. We demonstrate that ATG5 translocates to the nucleus, where it physically interacts with survivin in response to DNA-damaging agents both in vitro and in carcinoma tissues obtained from patients who had undergone radiotherapy and/or chemotherapy. As a consequence, elements of the chromosomal passenger complex are displaced during mitosis, resulting in chromosome misalignment and segregation defects. Pharmacological inhibition of autophagy does not prevent ATG5-dependent mitotic catastrophe, but shifts the balance to an early caspase-dependent cell death. Our data suggest a dual role for ATG5 in response to drug-induced DNA damage, where it acts in two signalling pathways in two distinct cellular compartments, the cytosol and the nucleus.

The Molecular Crosstalk between the MET Receptor Tyrosine Kinase and the DNA Damage Response - Biological and Clinical Aspects

Radiation therapy remains an imperative treatment modality for numerous malignancies. Enduring significant technical achievements both on the levels of treatment planning and radiation delivery have led to improvements in local control of tumor growth and reduction in healthy tissue toxicity. Nevertheless, resistance mechanisms, which presumably also involve activation of DNA damage response signaling pathways that eventually may account for loco-regional relapse and consequent tumor progression, still remain a critical problem. Accumulating data suggest that signaling via growth factor receptor tyrosine kinases, which are aberrantly expressed in many tumors, may interfere with the cytotoxic impact of ionizing radiation via the direct activation of the DNA damage response, leading eventually to so-called tumor radioresistance. The aim of this review is to overview the current known data that support a molecular crosstalk between the hepatocyte growth factor receptor tyrosine kinase MET and the DNA damage response. Apart of extending well established concepts over MET biology beyond its function as a growth factor receptor, these observations directly relate to the role of its aberrant activity in resistance to DNA damaging agents, such as ionizing radiation, which are routinely used in cancer therapy and advocate tumor sensitization towards DNA damaging agents in combination with MET targeting.