Interpretation of array comparative genome hybridization data: a major challenge.

The advent and application of high-resolution array-based comparative genome hybridization (array CGH) has led to the detection of large numbers of copy number variants (CNVs) in patients with developmental delay and/or multiple congenital anomalies as well as in healthy individuals. The notion that CNVs are also abundantly present in the normal population challenges the interpretation of the clinical significance of detected CNVs in patients. In this review we will illustrate a general clinical workflow based on our own experience that can be used in routine diagnostics for the interpretation of CNVs.

Novel MLPA procedure using self-designed probes allows comprehensive analysis for CNVs of the genes involved in Hirschsprung disease.

Background: Hirschsprung disease is characterized by the absence of intramural ganglion cells in the enteric plexuses, due to a fail during enteric nervous system formation. Hirschsprung has a complex genetic aetiology and mutations in several genes have been related to the disease. There is a clear predominance of missense/nonsense mutations in these genes whereas copy number variations (CNVs) have been seldom described, probably due to the limitations of conventional techniques usually employed for mutational analysis. In this study, we have looked for CNVs in some of the genes related to Hirschsprung (EDNRB, GFRA1, NRTN and PHOX2B) using the Multiple Ligation-dependent Probe Amplification (MLPA) approach. Methods: CNVs screening was performed in 208 HSCR patients using a self-designed set of MLPA probes, covering the coding region of those genes. Results: A deletion comprising the first 4 exons in GFRA1 gene was detected in 2 sporadic HSCR patients and in silico approaches have shown that the critical translation initiation signal in the mutant gene was abolished. In this study, we have been able to validate the reliability of this technique for CNVs screening in HSCR. Conclusions: The implemented MLPA based technique presented here allows CNV analysis of genes involved in HSCR that have not been not previously evaluated. Our results indicate that CNVs could be implicated in the pathogenesis of HSCR, although they seem to be an uncommon molecular cause of HSCR.

New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk.

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes.

Genome-wide arrays in routine diagnostics of hematological malignancies.

Over the last three decades, cytogenetic analysis of malignancies has become an integral part of disease evaluation and prediction of prognosis or responsiveness to therapy. In most diagnostic laboratories, conventional karyotyping, in conjunction with targeted fluorescence in situ hybridization analysis, is routinely performed to detect recurrent aberrations with prognostic implications. However, the genetic complexity of cancer cells requires a sensitive genome-wide analysis, enabling the detection of small genomic changes in a mixed cell population, as well as of regions of homozygosity. The advent of comprehensive high-resolution genomic tools, such as molecular karyotyping using comparative genomic hybridization or single-nucleotide polymorphism microarrays, has overcome many of the limitations of traditional cytogenetic techniques and has been used to study complex genomic lesions in, for example, leukemia. The clinical impact of the genomic copy-number and copy-neutral alterations identified by microarray technologies is growing rapidly and genome-wide array analysis is evolving into a diagnostic tool, to better identify high-risk patients and predict patients' outcomes from their genomic profiles. Here, we review the added clinical value of an array-based genome-wide screen in leukemia, and discuss the technical challenges and an interpretation workflow in applying arrays in the acquired cytogenetic diagnostic setting.

Structural variation and its effect on expression.

Structural variation, whether it is caused by copy number variants or present in a balanced form, such as reciprocal translocations and inversions, can have a profound and dramatic effect on the expression of genes mapping within and close to the rearrangement, as well as affecting others genome wide. These effects can be caused by altering the copy number of one or more genes or regulatory elements (dosage effect) or from physical disruption of links between regulatory elements and their associated gene or genes, resulting in perturbation of expression. Similarly, large-scale structural variants can result in genome-wide expression changes by altering the positions that chromosomes occupy within the nucleus, potentially disrupting not only local cis interactions, but also trans interactions that occur throughout the genome. Structural variation is, therefore, a significant factor in the study of gene expression and is discussed here in more detail.

Genetic and modifying factors that determine the risk of brain tumors

Some modifying factors may determine the risk of brain tumors. Until now, it could not be attempted to identify people at risk and also to improve significantly disease progression. Current therapy consists of surgical resection, followed by radiation therapy and chemotherapy. Despite of these treatments, the prognosis for patients is poor. In this review, we highlight general aspects concerning genetic alterations in brain tumors, namely astrocytomas, glioblastomas, oligodendrogliomas, medulloblastomas and ependymomas. The influence of these genetic alterations in patients' prognosis is discussed. Mutagen sensivity is associated with cancer risk. The convincing studies that linked DNA damages and DNA repair alterations with brain tumors are also described. Another important modifying factor is immunity. General immune response against cancer, tumor microenvironment and immune response, mechanisms of tumor escape, CNS tumor immunology, immune defects that impair anti-tumor systemic immunity in brain tumor patients and local immunosuppressive factors within CNS are also reviewed. New hope to treatment perspectives, as dendritic-cell-based vaccines is summarized too. Concluding, it seems well established that there is association between brain tumor risk and mutagen sensivity, which is highly heritable. Primary brain tumors cause depression in systemic host immunity; local immunosuppressive factors and immunological characteristics of tumor cells may explain the poor prognosis and DNA damages responses can alert immune system. However, it is necessary to clarify if individuals with both constitutional defects in immune functions and genetic instability have higher risk of developing brain tumors. Cytogenetic prospective studies and gene copy number variations analysis also must be performed in peripheral lymphocytes from brain tumor patients. © 2011 Bentham Science Publishers Ltd.

Detection of classical 17p11.2 deletions, an atypical deletion and RAI1 alterations in patients with features suggestive of Smith-Magenis syndrome

Smith-Magenis syndrome (SMS) is a complex disorder whose clinical features include mild to severe intellectual disability with speech delay, growth failure, brachycephaly, flat midface, short broad hands, and behavioral problems. SMS is typically caused by a large deletion on 17p11.2 that encompasses multiple genes including the retinoic acid induced 1, RAI1, gene or a mutation in the RAI1 gene. Here we have evaluated 30 patients with suspected SMS and identified SMS-associated classical 17p11.2 deletions in six patients, an atypical deletion of ∼139 kb that partially deletes the RAI1 gene in one patient, and RAI1 gene nonsynonymous alterations of unknown significance in two unrelated patients. The RAI1 mutant proteins showed no significant alterations in molecular weight, subcellular localization and transcriptional activity. Clinical features of patients with or without 17p11.2 deletions and mutations involving the RAI1 gene were compared to identify phenotypes that may be useful in diagnosing patients with SMS. © 2012 Macmillan Publishers Limited All rights reserved.

An Integrative Genomic and Transcriptomic Analysis Reveals Potential Targets Associated with Cell Proliferation in Uterine Leiomyomas

Background: Uterine Leiomyomas (ULs) are the most common benign tumours affecting women of reproductive age. ULs represent a major problem in public health, as they are the main indication for hysterectomy. Approximately 40-50% of ULs have non-random cytogenetic abnormalities, and half of ULs may have copy number alterations (CNAs). Gene expression microarrays studies have demonstrated that cell proliferation genes act in response to growth factors and steroids. However, only a few genes mapping to CNAs regions were found to be associated with ULs. Methodology: We applied an integrative analysis using genomic and transcriptomic data to identify the pathways and molecular markers associated with ULs. Fifty-one fresh frozen specimens were evaluated by array CGH (JISTIC) and gene expression microarrays (SAM). The CONEXIC algorithm was applied to integrate the data. Principal Findings: The integrated analysis identified the top 30 significant genes (P<0.01), which comprised genes associated with cancer, whereas the protein-protein interaction analysis indicated a strong association between FANCA and BRCA1. Functional in silico analysis revealed target molecules for drugs involved in cell proliferation, including FGFR1 and IGFBP5. Transcriptional and protein analyses showed that FGFR1 (P = 0.006 and P<0.01, respectively) and IGFBP5 (P = 0.0002 and P = 0.006, respectively) were up-regulated in the tumours when compared with the adjacent normal myometrium. Conclusions: The integrative genomic and transcriptomic approach indicated that FGFR1 and IGFBP5 amplification, as well as the consequent up-regulation of the protein products, plays an important role in the aetiology of ULs and thus provides data for potential drug therapies development to target genes associated with cellular proliferation in ULs. © 2013 Cirilo et al.

Genomic gains in prostatic carcinoma and proliferative inflammatory atrophy in dogs

The dog can spontaneously develop prostate cancer and consequently can be used as an experimental model for prostatic diseases associated with aging, including benign prostate hyperplasia (BPH) and prostate carcinoma (PCa). DNA copy number variations (CNVs) have been used to identify genes associated with cancer development and progression. DNA microarray based comparative genomic hybridization (aCGH) is a technique that allows to identify copy number of thousands of genes throughout the genome. aCGH was used to identify genomic regions with significantly different DNA copy number in three benign prostatic hyperplasia (BPH), four proliferative inflammatory atrophy (PIA), and 14 canine prostate carcinoma (PCa). Five histologically normal prostate tissue were used as reference. Genomic DNA was extracted from formalin fixed and paraffin embedded samples and CNVs data was evaluated in Canine Genome CGH Microarray 4x44K (G2519F, Design ID021193, Agilent). Data analysis was performed using Genomic Workbench Standard Edition 5.0.14 (Agilent). PCa showed higher number of altered genes related to canonical diseases process, cellular functions and molecular pathways as well as greater inter-relationship between genes, compared with PIA and BPH. In conclusion, PCa showed a more complex genotype, being losses the most frequent genomic changes. Some discrepancies between genomic alterations in human and canine carcinomas may indicate the different clinical behavior of these tumors in these two species. In addition, it was observed was an ascending pattern of genomic complexity in BPH, PIA and CA consistent with a model of multistep tumor progression.

KEAP1 E3 ligase-mediated downregulation of NF-kappaB signaling by targeting IKKbeta.

IkappaB kinase beta (IKKbeta) is involved in tumor development and progression through activation of the nuclear factor (NF)-kappaB pathway. However, the molecular mechanism that regulates IKKbeta degradation remains largely unknown. Here, we show that a Cullin 3 (CUL3)-based ubiquitin ligase, Kelch-like ECH-associated protein 1 (KEAP1), is responsible for IKKbeta ubiquitination. Depletion of KEAP1 led to the accumulation and stabilization of IKKbeta and to upregulation of NF-kappaB-derived tumor angiogenic factors. A systematic analysis of the CUL3, KEAP1, and RBX1 genomic loci revealed a high percentage of genome loss and missense mutations in human cancers that failed to facilitate IKKbeta degradation. Our results suggest that the dysregulation of KEAP1-mediated IKKbeta ubiquitination may contribute to tumorigenesis.

Genome-wide algorithm for detecting CNV associations with diseases

SNP genotyping arrays have been developed to characterize single-nucleotide polymorphisms (SNPs) and DNA copy number variations (CNVs). The quality of the inferences about copy number can be affected by many factors including batch effects, DNA sample preparation, signal processing, and analytical approach. Nonparametric and model-based statistical algorithms have been developed to detect CNVs from SNP genotyping data. However, these algorithms lack specificity to detect small CNVs due to the high false positive rate when calling CNVs based on the intensity values. Association tests based on detected CNVs therefore lack power even if the CNVs affecting disease risk are common. In this research, by combining an existing Hidden Markov Model (HMM) and the logistic regression model, a new genome-wide logistic regression algorithm was developed to detect CNV associations with diseases. We showed that the new algorithm is more sensitive and can be more powerful in detecting CNV associations with diseases than an existing popular algorithm, especially when the CNV association signal is weak and a limited number of SNPs are located in the CNV.^

8p23 beta-defensin copy number determination by single-locus pseudogene-based paralog ratio tests risk bias due to low-frequency sequence variations.

BACKGROUND The copy number variation (CNV) in beta-defensin genes (DEFB) on human chromosome 8p23 has been proposed to contribute to the phenotypic differences in inflammatory diseases. However, determination of exact DEFB CN is a major challenge in association studies. Quantitative real-time PCR (qPCR), paralog ratio tests (PRT) and multiplex ligation-dependent probe amplification (MLPA) have been extensively used to determine DEFB CN in different laboratories, but inter-method inconsistencies were observed frequently. In this study we asked which one is superior among the three methods for DEFB CN determination. RESULTS We developed a clustering approach for MLPA and PRT to statistically correlate data from a single experiment. Then we compared qPCR, a newly designed PRT and MLPA for DEFB CN determination in 285 DNA samples. We found MLPA had the best convergence and clustering results of the raw data and the highest call rate. In addition, the concordance rates between MLPA or PRT and qPCR (32.12% and 37.99%, respectively) were unacceptably low with underestimated CN by qPCR. Concordance rate between MLPA and PRT (90.52%) was high but PRT systematically underestimated CN by one in a subset of samples. In these samples a sequence variant which caused complete PCR dropout of the respective DEFB cluster copies was found in one primer binding site of one of the targeted paralogous pseudogenes. CONCLUSION MLPA is superior to PRT and even more to qPCR for DEFB CN determination. Although the applied PRT provides in most cases reliable results, such a test is particularly sensitive to low-frequency sequence variations preferably accumulating in loci like pseudogenes which are most likely not under selective pressure. In the light of the superior performance of multiplex assays, the drawbacks of such single PRTs could be overcome by combining more test markers.

DNA Topoisomerase I Gene Copy Number and mRNA Expression Assessed as Predictive Biomarkers for Adjuvant Irinotecan in Stage II/III Colon Cancer.

PURPOSE: Prospective-retrospective assessment of theTOP1gene copy number andTOP1mRNA expression as predictive biomarkers for adjuvant irinotecan in stage II/III colon cancer. EXPERIMENTAL DESIGN: Formalin-fixed, paraffin-embedded tissue microarrays were obtained from an adjuvant colon cancer trial (PETACC3) where patients were randomized to 5-fluorouracil/folinic acid with or without additional irinotecan.TOP1copy number status was analyzed by fluorescencein situhybridization (FISH) using aTOP1/CEN20 dual-probe combination.TOP1mRNA data were available from previous analyses. RESULTS: TOP1FISH and follow-up data were obtained from 534 patients.TOP1gain was identified in 27% using a single-probe enumeration strategy (≥4TOP1signals per cell) and in 31% when defined by aTOP1/CEN20 ratio ≥ 1.5. The effect of additional irinotecan was not dependent onTOP1FISH status.TOP1mRNA data were available from 580 patients with stage III disease. Benefit of irinotecan was restricted to patients characterized byTOP1mRNA expression ≥ third quartile (RFS: HRadjusted, 0.59;P= 0.09; OS: HRadjusted, 0.44;P= 0.03). The treatment byTOP1mRNA interaction was not statistically significant, but in exploratory multivariable fractional polynomial interaction analysis, increasingTOP1mRNA values appeared to be associated with increasing benefit of irinotecan. CONCLUSIONS: In contrast to theTOP1copy number, a trend was demonstrated for a predictive property ofTOP1mRNA expression. On the basis ofTOP1mRNA, it might be possible to identify a subgroup of patients where an irinotecan doublet is a clinically relevant option in the adjuvant setting of colon cancer.Clin Cancer Res; 22(7); 1621-31. ©2015 AACR.

Quantifying single gene copy number by measuring fluorescent probe lengths on combed genomic DNA

An approach was developed for the quantification of subtle gains and losses of genomic DNA. The approach relies on a process called molecular combing. Molecular combing consists of the extension and alignment of purified molecules of genomic DNA on a glass coverslip. It has the advantage that a large number of genomes can be combed per coverslip, which allows for a statistically adequate number of measurements to be made on the combed DNA. Consequently, a high-resolution approach to mapping and quantifying genomic alterations is possible. The approach consists of applying fluorescence hybridization to the combed DNA by using probes to identify the amplified region. Measurements then are made on the linear hybridization signals to ascertain the region's exact size. The reliability of the approach first was tested for low copy number amplifications by determining the copy number of chromosome 21 in a normal and trisomy 21 cell line. It then was tested for high copy number amplifications by quantifying the copy number of an oncogene amplified in the tumor cell line GTL-16. These results demonstrate that a wide range of amplifications can be accurately and reliably quantified. The sensitivity and resolution of the approach likewise was assessed by determining the copy number of a single allele (160 kb) alteration.