AICAR confers prophylactic cardioprotection in doxorubicin-induced heart failure in rats.

Doxorubicin (DOX) is a widely used chemotherapeutic agent that can cause serious cardiotoxic side effects, leading to heart failure (HF). Impaired mitochondrial function is thought to be key factor driving progression into HF. We have previously shown in a rat model of DOX-HF that heart failure with reduced ejection fraction correlates with mitochondrial loss and dysfunction. Adenosine monophosphate-dependent kinase (AMPK) is a cellular energy sensor, regulating mitochondrial biogenesis and energy metabolism, including fatty acid oxidation. We hypothesised that AMPK activation could restore mitochondrial function and therefore be a novel cardioprotective strategy for the prevention of DOX-HF. Consequently, we set out to assess whether 5-aminoimidazole-4-carboxamide 1-ß-D-ribofuranoside (AICAR), an activator of AMPK, could prevent cardiac functional decline in this chronic intravenous rat model of DOX-HF. In line with our hypothesis, AICAR improved cardiac systolic function. AICAR furthermore improved cardiac mitochondrial fatty acid oxidation, independent of mitochondrial number, and in the absence of observable AMPK-activation. In addition, we found that AICAR prevented loss of myocardial mass. RNAseq analysis showed that this may be driven by normalisation of pathways associated with ribosome function and protein synthesis, which are impaired in DOX-treated rat hearts. AICAR furthermore prevented dyslipidemia and excessive body-weight loss in DOX-treated rats, which may contribute to preservation of myocardial mass. Though it is unclear whether AICAR exerted its cardioprotective effect through cardiac or extra-cardiac AMPK-activation or via an AMPK-independent effect, these results show promise for the use of AICAR as a cardioprotective agent in DOX-HF to both preserve cardiac function and mass.

Decreased ATM Function Causes Delayed DNA Repair and Apoptosis in Common Variable Immunodeficiency Disorders.

PURPOSE: Common variable immunodeficiency disorders (CVID) is characterized by low/absent serum immunoglobulins and susceptibility to bacterial infection. Patients can develop an infections-only phenotype or a complex disease course with inflammatory, autoimmune, and/or malignant complications. We hypothesized that deficient DNA repair mechanisms may be responsible for the antibody deficiency and susceptibility to inflammation and cancer in some patients. METHODS: Germline variants were identified following targeted sequencing of n = 252 genes related to DNA repair in n = 38 patients. NanoString nCounter PlexSet assay measured gene expression in n = 20 CVID patients and n = 7 controls. DNA damage and apoptosis were assessed by flow cytometry in n = 34 CVID patients and n = 11 controls. RESULTS: Targeted sequencing supported enrichment of rare genetic variants in genes related to DNA repair pathways with novel and rare likely pathogenic variants identified and an altered gene expression signature that distinguished patients from controls and complex patients from those with an infections-only phenotype. Consistent with this, flow cytometric analyses of lymphocytes following DNA damage revealed a subset of CVID patients whose immune cells have downregulated ATM, impairing the recruitment of other repair factors, delaying repair and promoting apoptosis. CONCLUSION: These data suggest that germline genetics and altered gene expression predispose a subset of CVID patients to increased sensitivity to DNA damage and reduced DNA repair capacity.

Impact of spliceosome mutations on RNA splicing in myelodysplasia: dysregulated genes/pathways and clinical associations.

SF3B1, SRSF2, and U2AF1 are the most frequently mutated splicing factor genes in the myelodysplastic syndromes (MDS). We have performed a comprehensive and systematic analysis to determine the effect of these commonly mutated splicing factors on pre-mRNA splicing in the bone marrow stem/progenitor cells and in the erythroid and myeloid precursors in splicing factor mutant MDS. Using RNA-seq, we determined the aberrantly spliced genes and dysregulated pathways in CD34+ cells of 84 patients with MDS. Splicing factor mutations result in different alterations in splicing and largely affect different genes, but these converge in common dysregulated pathways and cellular processes, focused on RNA splicing, protein synthesis, and mitochondrial dysfunction, suggesting common mechanisms of action in MDS. Many of these dysregulated pathways and cellular processes can be linked to the known disease pathophysiology associated with splicing factor mutations in MDS, whereas several others have not been previously associated with MDS, such as sirtuin signaling. We identified aberrantly spliced events associated with clinical variables, and isoforms that independently predict survival in MDS and implicate dysregulation of focal adhesion and extracellular exosomes as drivers of poor survival. Aberrantly spliced genes and dysregulated pathways were identified in the MDS-affected lineages in splicing factor mutant MDS. Functional studies demonstrated that knockdown of the mitosis regulators SEPT2 and AKAP8, aberrantly spliced target genes of SF3B1 and SRSF2 mutations, respectively, led to impaired erythroid cell growth and differentiation. This study illuminates the effect of the common spliceosome mutations on the MDS phenotype and provides novel insights into disease pathophysiology.

The nature and nurture of cell heterogeneity: accounting for macrophage gene-environment interactions with single-cell RNA-Seq.

BACKGROUND: Single-cell RNA-Seq can be a valuable and unbiased tool to dissect cellular heterogeneity, despite the transcriptome's limitations in describing higher functional phenotypes and protein events. Perhaps the most important shortfall with transcriptomic 'snapshots' of cell populations is that they risk being descriptive, only cataloging heterogeneity at one point in time, and without microenvironmental context. Studying the genetic ('nature') and environmental ('nurture') modifiers of heterogeneity, and how cell population dynamics unfold over time in response to these modifiers is key when studying highly plastic cells such as macrophages. RESULTS: We introduce the programmable Polaris™ microfluidic lab-on-chip for single-cell sequencing, which performs live-cell imaging while controlling for the culture microenvironment of each cell. Using gene-edited macrophages we demonstrate how previously unappreciated knockout effects of SAMHD1, such as an altered oxidative stress response, have a large paracrine signaling component. Furthermore, we demonstrate single-cell pathway enrichments for cell cycle arrest and APOBEC3G degradation, both associated with the oxidative stress response and altered proteostasis. Interestingly, SAMHD1 and APOBEC3G are both HIV-1 inhibitors ('restriction factors'), with no known co-regulation. CONCLUSION: As single-cell methods continue to mature, so will the ability to move beyond simple 'snapshots' of cell populations towards studying the determinants of population dynamics. By combining single-cell culture, live-cell imaging, and single-cell sequencing, we have demonstrated the ability to study cell phenotypes and microenvironmental influences. It's these microenvironmental components - ignored by standard single-cell workflows - that likely determine how macrophages, for example, react to inflammation and form treatment resistant HIV reservoirs.

Clinical whole-genome sequencing in severe early-onset epilepsy reveals new genes and improves molecular diagnosis.

In severe early-onset epilepsy, precise clinical and molecular genetic diagnosis is complex, as many metabolic and electro-physiological processes have been implicated in disease causation. The clinical phenotypes share many features such as complex seizure types and developmental delay. Molecular diagnosis has historically been confined to sequential testing of candidate genes known to be associated with specific sub-phenotypes, but the diagnostic yield of this approach can be low. We conducted whole-genome sequencing (WGS) on six patients with severe early-onset epilepsy who had previously been refractory to molecular diagnosis, and their parents. Four of these patients had a clinical diagnosis of Ohtahara Syndrome (OS) and two patients had severe non-syndromic early-onset epilepsy (NSEOE). In two OS cases, we found de novo non-synonymous mutations in the genes KCNQ2 and SCN2A. In a third OS case, WGS revealed paternal isodisomy for chromosome 9, leading to identification of the causal homozygous missense variant in KCNT1, which produced a substantial increase in potassium channel current. The fourth OS patient had a recessive mutation in PIGQ that led to exon skipping and defective glycophosphatidyl inositol biosynthesis. The two patients with NSEOE had likely pathogenic de novo mutations in CBL and CSNK1G1, respectively. Mutations in these genes were not found among 500 additional individuals with epilepsy. This work reveals two novel genes for OS, KCNT1 and PIGQ. It also uncovers unexpected genetic mechanisms and emphasizes the power of WGS as a clinical tool for making molecular diagnoses, particularly for highly heterogeneous disorders.

Coexpression network analysis in abdominal and gluteal adipose tissue reveals regulatory genetic loci for metabolic syndrome and related phenotypes.

Metabolic Syndrome (MetS) is highly prevalent and has considerable public health impact, but its underlying genetic factors remain elusive. To identify gene networks involved in MetS, we conducted whole-genome expression and genotype profiling on abdominal (ABD) and gluteal (GLU) adipose tissue, and whole blood (WB), from 29 MetS cases and 44 controls. Co-expression network analysis for each tissue independently identified nine, six, and zero MetS-associated modules of coexpressed genes in ABD, GLU, and WB, respectively. Of 8,992 probesets expressed in ABD or GLU, 685 (7.6%) were expressed in ABD and 51 (0.6%) in GLU only. Differential eigengene network analysis of 8,256 shared probesets detected 22 shared modules with high preservation across adipose depots (D(ABD-GLU) = 0.89), seven of which were associated with MetS (FDR P<0.01). The strongest associated module, significantly enriched for immune response-related processes, contained 94/620 (15%) genes with inter-depot differences. In an independent cohort of 145/141 twins with ABD and WB longitudinal expression data, median variability in ABD due to familiality was greater for MetS-associated versus un-associated modules (ABD: 0.48 versus 0.18, P = 0.08; GLU: 0.54 versus 0.20, P = 7.8×10(-4)). Cis-eQTL analysis of probesets associated with MetS (FDR P<0.01) and/or inter-depot differences (FDR P<0.01) provided evidence for 32 eQTLs. Corresponding eSNPs were tested for association with MetS-related phenotypes in two GWAS of >100,000 individuals; rs10282458, affecting expression of RARRES2 (encoding chemerin), was associated with body mass index (BMI) (P = 6.0×10(-4)); and rs2395185, affecting inter-depot differences of HLA-DRB1 expression, was associated with high-density lipoprotein (P = 8.7×10(-4)) and BMI-adjusted waist-to-hip ratio (P = 2.4×10(-4)). Since many genes and their interactions influence complex traits such as MetS, integrated analysis of genotypes and coexpression networks across multiple tissues relevant to clinical traits is an efficient strategy to identify novel associations.

Pervasive haplotypic variation in the spliceo-transcriptome of the human major histocompatibility complex.

The human major histocompatibility complex (MHC) on chromosome 6p21 is a paradigm for genomics, showing remarkable polymorphism and striking association with immune and non-immune diseases. The complex genomic landscape of the MHC, notably strong linkage disequilibrium, has made resolving causal variants very challenging. A promising approach is to investigate gene expression levels considered as tractable intermediate phenotypes in mapping complex diseases. However, how transcription varies across the MHC, notably relative to specific haplotypes, remains unknown. Here, using an original hybrid tiling and splice junction microarray that includes alternate allele probes, we draw the first high-resolution strand-specific transcription map for three common MHC haplotypes (HLA-A1-B8-Cw7-DR3, HLA-A3-B7-Cw7-DR15, and HLA-A26-B18-Cw5-DR3-DQ2) strongly associated with autoimmune diseases including type 1 diabetes, systemic lupus erythematosus, and multiple sclerosis. We find that haplotype-specific differences in gene expression are common across the MHC, affecting 96 genes (46.4%), most significantly the zing finger protein gene ZFP57. Differentially expressed probes are correlated with polymorphisms between haplotypes, consistent with cis effects that we directly demonstrate for ZFP57 in a cohort of healthy volunteers (P = 1.2 × 10(-14)). We establish that alternative splicing is significantly more frequent in the MHC than genome-wide (72.5% vs. 62.1% of genes, P ≤ 1 × 10(-4)) and shows marked haplotypic differences. We also unmask novel and abundant intergenic transcription involving 31% of transcribed blocks identified. Our study reveals that the renowned MHC polymorphism also manifests as transcript diversity, and our novel haplotype-based approach marks a new step toward identification of regulatory variants involved in the control of MHC-associated phenotypes and diseases.

GREVE: Genomic Recurrent Event ViEwer to assist the identification of patterns across individual cancer samples.

SUMMARY: GREVE has been developed to assist with the identification of recurrent genomic aberrations across cancer samples. The exact characterization of such aberrations remains a challenge despite the availability of increasing amount of data, from SNParray to next-generation sequencing. Furthermore, genomic aberrations in cancer are especially difficult to handle because they are, by nature, unique to the patients. However, their recurrence in specific regions of the genome has been shown to reflect their relevance in the development of tumors. GREVE makes use of previously characterized events to identify such regions and focus any further analysis. AVAILABILITY: GREVE is available through a web interface and open-source application (http://www.well.ox.ac.uk/GREVE).

Positional cloning of a quantitative trait locus on chromosome 13q14 that influences immunoglobulin E levels and asthma.

Atopic or immunoglobulin E (IgE)-mediated diseases include the common disorders of asthma, atopic dermatitis and allergic rhinitis. Chromosome 13q14 shows consistent linkage to atopy and the total serum IgE concentration. We previously identified association between total serum IgE levels and a novel 13q14 microsatellite (USAT24G1; ref. 7) and have now localized the underlying quantitative-trait locus (QTL) in a comprehensive single-nucleotide polymorphism (SNP) map. We found replicated association to IgE levels that was attributed to several alleles in a single gene, PHF11. We also found association with these variants to severe clinical asthma. The gene product (PHF11) contains two PHD zinc fingers and probably regulates transcription. Distinctive splice variants were expressed in immune tissues and cells.

Positional cloning of a novel gene influencing asthma from chromosome 2q14.

Asthma is a common disease in children and young adults. Four separate reports have linked asthma and related phenotypes to an ill-defined interval between 2q14 and 2q32 (refs. 1-4), and two mouse genome screens have linked bronchial hyper-responsiveness to the region homologous to 2q14 (refs. 5,6). We found and replicated association between asthma and the D2S308 microsatellite, 800 kb distal to the IL1 cluster on 2q14. We sequenced the surrounding region and constructed a comprehensive, high-density, single-nucleotide polymorphism (SNP) linkage disequilibrium (LD) map. SNP association was limited to the initial exons of a solitary gene of 3.6 kb (DPP10), which extends over 1 Mb of genomic DNA. DPP10 encodes a homolog of dipeptidyl peptidases (DPPs) that cleave terminal dipeptides from cytokines and chemokines, and it presents a potential new target for asthma therapy.

Genomic regulation of Krüppel-like-factor family members by corticosteroid receptors in the rat brain.

Hippocampal mineralocorticoid receptors (MRs) and glucocorticoid receptors (GRs) mediate glucocorticoid hormone (GC) action in the hippocampus. These receptors bind to glucocorticoid responsive elements (GREs) within target genes, eliciting transcriptional effects in response to stress and circadian variation. Until recently, little was known about the genome-wide targets of hippocampal MRs and GRs under physiological conditions. Following on from our genome-wide MR and GR ChIP-seq and Ribo-Zero RNA-seq studies on rat hippocampus, we investigated the Krüppel-like factors (KLFs) as targets of MRs and GRs throughout the brain under circadian variation and after acute stress. In particular, Klf2, Klf9 and Klf15 are known to be stress and/or GC responsive and play a role in neurobiological processes including synaptic plasticity and neuronal differentiation. We found increased binding of MR and GR to GREs within Klf2, Klf9 and Klf15 in the hippocampus, amygdala, prefrontal cortex, and neocortex after acute stress and resulting from circadian variation, which was accompanied by upregulation of corresponding hnRNA and mRNA levels. Adrenalectomy abolished transcriptional upregulation of specific Klf genes. These results show that MRs and GRs regulate Klf gene expression throughout the brain following exposure to acute stress or in response to circadian variation, likely alongside other transcription factors.

miR-3156-5p is downregulated in serum of MEN1 patients and regulates expression of MORF4L2.

Multiple endocrine neoplasia type 1 (MEN1), caused by mutations in the MEN1 gene encoding menin, is an autosomal dominant disorder characterised by the combined occurrence of parathyroid, pituitary and pancreatic neuroendocrine tumours (NETs). Development of these tumours is associated with wide variations in their severity, order and ages (from <5 to >80 years), requiring life-long screening. To improve tumour surveillance and quality of life, better circulating biomarkers, particularly for pancreatic NETs that are associated with higher mortality, are required. We, therefore, examined the expression of circulating miRNA in the serum of MEN1 patients. Initial profiling analysis followed by qRT-PCR validation studies identified miR-3156-5p to be significantly downregulated (-1.3 to 5.8-fold, P < 0.05-0.0005) in nine MEN1 patients, compared to matched unaffected relatives. MEN1 knock-down experiments in BON-1 human pancreatic NET cells resulted in reduced MEN1 (49%, P < 0.05), menin (54%, P < 0.05) and miR-3156-5p expression (20%, P < 0.005), compared to control-treated cells, suggesting that miR-3156-5p downregulation is a consequence of loss of MEN1 expression. In silico analysis identified mortality factor 4-like 2 (MOR4FL2) as a potential target of miR-3156-5p, and in vitro functional studies in BON-1 cells transfected with either miR-3156-5p mimic or inhibitors showed that the miR-3156-5p mimic significantly reduced MORF4L2 protein expression (46%, P < 0.005), while miR-3156-5p inhibitor significantly increased MORF4L2 expression (1.5-fold, P < 0.05), compared to control-treated cells, thereby confirming that miR-3156-5p regulates MORF4L2 expression. Thus, the inverse relationship between miR-3156-5p and MORF4L2 expression represents a potential serum biomarker that could facilitate the detection of NET occurrence in MEN1 patients.

Distinct regulation of hippocampal neuroplasticity and ciliary genes by corticosteroid receptors.

Glucocorticoid hormones (GCs) - acting through hippocampal mineralocorticoid receptors (MRs) and glucocorticoid receptors (GRs) - are critical to physiological regulation and behavioural adaptation. We conducted genome-wide MR and GR ChIP-seq and Ribo-Zero RNA-seq studies on rat hippocampus to elucidate MR- and GR-regulated genes under circadian variation or acute stress. In a subset of genes, these physiological conditions resulted in enhanced MR and/or GR binding to DNA sequences and associated transcriptional changes. Binding of MR at a substantial number of sites however remained unchanged. MR and GR binding occur at overlapping as well as distinct loci. Moreover, although the GC response element (GRE) was the predominant motif, the transcription factor recognition site composition within MR and GR binding peaks show marked differences. Pathway analysis uncovered that MR and GR regulate a substantial number of genes involved in synaptic/neuro-plasticity, cell morphology and development, behavior, and neuropsychiatric disorders. We find that MR, not GR, is the predominant receptor binding to >50 ciliary genes; and that MR function is linked to neuronal differentiation and ciliogenesis in human fetal neuronal progenitor cells. These results show that hippocampal MRs and GRs constitutively and dynamically regulate genomic activities underpinning neuronal plasticity and behavioral adaptation to changing environments.

Early detection of doxorubicin-induced cardiotoxicity in rats by its cardiac metabolic signature assessed with hyperpolarized MRI.

Doxorubicin (DOX) is a widely used chemotherapeutic agent that can cause serious cardiotoxic side effects culminating in congestive heart failure (HF). There are currently no clinical imaging techniques or biomarkers available to detect DOX-cardiotoxicity before functional decline. Mitochondrial dysfunction is thought to be a key factor driving functional decline, though real-time metabolic fluxes have never been assessed in DOX-cardiotoxicity. Hyperpolarized magnetic resonance imaging (MRI) can assess real-time metabolic fluxes in vivo. Here we show that cardiac functional decline in a clinically relevant rat-model of DOX-HF is preceded by a change in oxidative mitochondrial carbohydrate metabolism, measured by hyperpolarized MRI. The decreased metabolic fluxes were predominantly due to mitochondrial loss and additional mitochondrial dysfunction, and not, as widely assumed hitherto, to oxidative stress. Since hyperpolarized MRI has been successfully translated into clinical trials this opens up the potential to test cancer patients receiving DOX for early signs of cardiotoxicity.

Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans.

Several strands of evidence question the dogma that human mitochondrial DNA (mtDNA) is inherited exclusively down the maternal line, most recently in three families where several individuals harbored a 'heteroplasmic haplotype' consistent with biparental transmission. Here we report a similar genetic signature in 7 of 11,035 trios, with allelic fractions of 5-25%, implying biparental inheritance of mtDNA in 0.06% of offspring. However, analysing the nuclear whole genome sequence, we observe likely large rare or unique nuclear-mitochondrial DNA segments (mega-NUMTs) transmitted from the father in all 7 families. Independently detecting mega-NUMTs in 0.13% of fathers, we see autosomal transmission of the haplotype. Finally, we show the haplotype allele fraction can be explained by complex concatenated mtDNA-derived sequences rearranged within the nuclear genome. We conclude that rare cryptic mega-NUMTs can resemble paternally mtDNA heteroplasmy, but find no evidence of paternal transmission of mtDNA in humans.

Author Correction: Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

GLIDERS--a web-based search engine for genome-wide linkage disequilibrium between HapMap SNPs.

BACKGROUND: A number of tools for the examination of linkage disequilibrium (LD) patterns between nearby alleles exist, but none are available for quickly and easily investigating LD at longer ranges (>500 kb). We have developed a web-based query tool (GLIDERS: Genome-wide LInkage DisEquilibrium Repository and Search engine) that enables the retrieval of pairwise associations with r2 >or= 0.3 across the human genome for any SNP genotyped within HapMap phase 2 and 3, regardless of distance between the markers. DESCRIPTION: GLIDERS is an easy to use web tool that only requires the user to enter rs numbers of SNPs they want to retrieve genome-wide LD for (both nearby and long-range). The intuitive web interface handles both manual entry of SNP IDs as well as allowing users to upload files of SNP IDs. The user can limit the resulting inter SNP associations with easy to use menu options. These include MAF limit (5-45%), distance limits between SNPs (minimum and maximum), r2 (0.3 to 1), HapMap population sample (CEU, YRI and JPT+CHB combined) and HapMap build/release. All resulting genome-wide inter-SNP associations are displayed on a single output page, which has a link to a downloadable tab delimited text file. CONCLUSION: GLIDERS is a quick and easy way to retrieve genome-wide inter-SNP associations and to explore LD patterns for any number of SNPs of interest. GLIDERS can be useful in identifying SNPs with long-range LD. This can highlight mis-mapping or other potential association signal localisation problems.

Phenotypic and transcriptomic characterization of canine myeloid-derived suppressor cells.

Myeloid-derived suppressor cells (MDSCs) are key players in immune evasion, tumor progression and metastasis. MDSCs accumulate under various pathological states and fall into two functionally and phenotypically distinct subsets that have been identified in humans and mice: polymorphonuclear (PMN)-MDSCs and monocytic (M)-MDSCs. As dogs are an excellent model for human tumor development and progression, we set out to identify PMN-MDSCs and M-MDSCs in clinical canine oncology patients. Canine hypodense MHC class II-CD5-CD21-CD11b+ cells can be subdivided into polymorphonuclear (CADO48A+CD14-) and monocytic (CADO48A-CD14+) MDSC subsets. The transcriptomic signatures of PMN-MDSCs and M-MDSCs are distinct, and moreover reveal a statistically significant similarity between canine and previously published human PMN-MDSC gene expression patterns. As in humans, peripheral blood frequencies of canine PMN-MDSCs and M-MDSCs are significantly higher in dogs with cancer compared to healthy control dogs (PMN-MDSCs: p < 0.001; M-MDSCs: p < 0.01). By leveraging the power of evolution, we also identified additional conserved genes in PMN-MDSCs of multiple species that may play a role in MDSC function. Our findings therefore validate the dog as a model for studying MDSCs in the context of cancer.

A practical solution for preserving single cells for RNA sequencing.

The design and implementation of single-cell experiments is often limited by their requirement for fresh starting material. We have adapted a method for histological tissue fixation using dithio-bis(succinimidyl propionate) (DSP), or Lomant's Reagent, to stabilise cell samples for single-cell transcriptomic applications. DSP is a reversible cross-linker of free amine groups that has previously been shown to preserve tissue integrity for histology while maintaining RNA integrity and yield in bulk RNA extractions. Although RNA-seq data from DSP-fixed single cells appears to be prone to characteristic artefacts, such as slightly reduced yield of cDNA and a detectable 3' bias in comparison with fresh cells, cell preservation using DSP does not appear to substantially reduce RNA complexity at the gene level. In addition, there is evidence that instantaneous fixation of cells can reduce inter-cell technical variability. The ability of DSP-fixed cells to retain commonly used dyes, such as propidium iodide, enables the tracking of experimental sub-populations and the recording of cell viability at the point of fixation. Preserving cells using DSP will remove several barriers in the staging of single-cell experiments, including the transport of samples and the scheduling of shared equipment for downstream single-cell isolation and processing.

SW-ARRAY: a dynamic programming solution for the identification of copy-number changes in genomic DNA using array comparative genome hybridization data.

Comparative genome hybridization (CGH) to DNA microarrays (array CGH) is a technique capable of detecting deletions and duplications in genomes at high resolution. However, array CGH studies of the human genome noting false negative and false positive results using large insert clones as probes have raised important concerns regarding the suitability of this approach for clinical diagnostic applications. Here, we adapt the Smith-Waterman dynamic-programming algorithm to provide a sensitive and robust analytic approach (SW-ARRAY) for detecting copy-number changes in array CGH data. In a blind series of hybridizations to arrays consisting of the entire tiling path for the terminal 2 Mb of human chromosome 16p, the method identified all monosomies between 267 and 1567 kb with a high degree of statistical significance and accurately located the boundaries of deletions in the range 267-1052 kb. The approach is unique in offering both a nonparametric segmentation procedure and a nonparametric test of significance. It is scalable and well-suited to high resolution whole genome array CGH studies that use array probes derived from large insert clones as well as PCR products and oligonucleotides.