CopywriteR: DNA copy number detection from off-target sequence data.

Current methods for detection of copy number variants (CNV) and aberrations (CNA) from targeted sequencing data are based on the depth of coverage of captured exons. Accurate CNA determination is complicated by uneven genomic distribution and non-uniform capture efficiency of targeted exons. Here we present CopywriteR, which eludes these problems by exploiting 'off-target' sequence reads. CopywriteR allows for extracting uniformly distributed copy number information, can be used without reference, and can be applied to sequencing data obtained from various techniques including chromatin immunoprecipitation and target enrichment on small gene panels. CopywriteR outperforms existing methods and constitutes a widely applicable alternative to available tools.

Clinical relevance of polymorphic DNA copy number variation (CNV) in African American women with stage I-II breast cancer

Objectives. The chief goal of this study was to analyze copy number variation (CNV) in breast cancer tumors from 25 African American women with early stage breast cancer (BC) using molecular inversion probes (MIP) in order to: (1) compare the degree of CNV in tumors compared to normal lymph nodes, and (2) determine whether gains and/or losses of genes in specific chromosomes differ between pathologic subtypes of breast cancer defined by known prognostic markers, (3) determine whether gains/losses in CN are associated with known oncogenes or tumor suppressor genes, and (4) determine whether increased gains/losses in CN for specific chromosomes were associated with differences in breast cancer recurrence. ^ Methods. Twenty to 37 nanograms of DNA extracted from 25 formalin-fixed paraffin embedded (FFPE) tumor samples and matched normal lymph nodes were added to individual tubes. Oligonucleotide probes with recognition sequences at each terminus were hybridized with a genomic target sequence to form a circular structure. Probes are released from genomic DNA obtained from FFPE samples, and those which have been correctly "circularized" in the proper allele/nucleotide reaction combination are amplified using polymerase chain reaction (PCR) primers. Amplicons were fluorescently labeled and the tag sequences released from the genome homology regions by treatment with uracil-N-glycosylase to cleave the probe at the site where uracils are present, and detected using a complementary tag array developed by Affymetrix. ^ Results. Analysis of CN gains and losses from tumors and normal tissues showed marked differences in tumors with numerous chromosomes affected. Similar changes were not observed in normal lymph nodes. When tumors were stratified into four groups based on expression or lack of expression of the estrogen receptor and HER2/neu, distinct patterns of CNV for different chromosomes were observed. Gains or losses in CN for specific chromosomes correlated with amplifications/deletions of particular oncogenes or tumor suppressor genes (i.e. such as found on chromosome 17) known to be associated with aggressive tumor phenotype and poor prognosis. There was a trend for increases in CN observed for chromosome 17 to correlate inversely with time to recurrence of BC (p=0.14 for trend). CNV was also observed for chromosomes 5, 8, 10, 11, and 16, which are known sites for several breast cancer susceptibility alleles. ^ Conclusions. This study is the first to validate the MIP technique, to correlate differences in gene expression with known prognostic tumor markers, and to correlate significant increases/decreases in CN with known tumor markers associated with prognosis. The results of this study may have far reaching public health implications towards identifying new high-risk groups based on genomic differences in CNP, both with respect to prognosis and response to therapy, and to eventually identify new therapeutic targets for prevention and treatment of this disease. ^

Quantitative analysis of chromosomal CGH in human breast tumors associates copy number abnormalities with p53 status and patient survival

We present a general method for rigorously identifying correlations between variations in large-scale molecular profiles and outcomes and apply it to chromosomal comparative genomic hybridization data from a set of 52 breast tumors. We identify two loci where copy number abnormalities are correlated with poor survival outcome (gain at 8q24 and loss at 9q13). We also identify a relationship between abnormalities at two loci and the mutational status of p53. Gain at 8q24 and loss at 5q15-5q21 are linked with mutant p53. The 9q and 5q losses suggest the possibility of gene products involved in breast cancer progression. The analytical techniques are general and also are applicable to the analysis of array-based expression data.

High copy number I-Ab transgenes induce production of IgE through an interluekin 4-dependent mechanism.

To better understand the role of class II major histocompatibility complex molecules in both normal and autoimmune responses, we have produced a series of I-Ab transgenic mice. One of these transgenic constructs, designated NOD.PD, has the sequence of the NOD beta chain (Abeta(g7)) except at positions 56 and 57, where Pro-Asp replaces His-Ser. Several NOD.PD transgenic lines have been produced. One line of these mice carried a very high number of copies (>50) of the NOD.PD transgene. As has been described in other mice carrying high copy numbers of I-Ab transgenes, B-cell development was abnormal. The steady state numbers of mature B cells (IgM+/IgD(hi)) in the periphery were greatly reduced in transgenic mice compared to nontransgenic littermates. Surprisingly, rather than being accompanied by a generalized hypogammaglobulinemia, this B-cell deficiency was accompanied by elevated concentrations of IgG1 and IgE in the serum. Conversely, the levels of IgG2a were reduced in transgenic mice compared to nontransgenic littermates. Because this isotype pattern was characteristic of interleukin (IL)-4-induced class-switching, we then investigated the role of IL-4 in causing the observed phenotype. We crossed the high copy number transgenic mice with an IL-4-deficient strain of mice. As expected, the elevated levels of IgE in high copy number transgenic mice were eliminated when the IL-4 gene was inactivated. However, the reduction in the number of B cells was not ameliorated. These data indicate that the primary defect caused by the transgene was to reduce the number of B cells in these mice. This reduction was accompanied by a secondary increase in IL-4 production, which drove the remaining B cells toward the production of IgGl and IgE.

A compendium of myeloma-associated chromosomal copy number abnormalities and their prognostic value.

To obtain a comprehensive genomic profile of presenting multiple myeloma cases we performed high-resolution single nucleotide polymorphism mapping array analysis in 114 samples alongside 258 samples analyzed by U133 Plus 2.0 expression array (Affymetrix). We examined DNA copy number alterations and loss of heterozygosity (LOH) to define the spectrum of minimally deleted regions in which relevant genes of interest can be found. The most frequent deletions are located at 1p (30%), 6q (33%), 8p (25%), 12p (15%), 13q (59%), 14q (39%), 16q (35%), 17p (7%), 20 (12%), and 22 (18%). In addition, copy number-neutral LOH, or uniparental disomy, was also prevalent on 1q (8%), 16q (9%), and X (20%), and was associated with regions of gain and loss. Based on fluorescence in situ hybridization and expression quartile analysis, genes of prognostic importance were found to be located at 1p (FAF1, CDKN2C), 1q (ANP32E), and 17p (TP53). In addition, we identified common homozygously deleted genes that have functions relevant to myeloma biology. Taken together, these analyses indicate that the crucial pathways in myeloma pathogenesis include the nuclear factor-κB pathway, apoptosis, cell-cycle regulation, Wnt signaling, and histone modifications. This study was registered at http://isrctn.org as ISRCTN68454111.

Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits.

Background: Copy number variations (CNVs) have been shown to account for substantial portions of observed genomic variation and have been associated with qualitative and quantitative traits and the onset of disease in a number of species. Information from high-resolution studies to detect, characterize and estimate population-specific variant frequencies will facilitate the incorporation of CNVs in genomic studies to identify genes affecting traits of importance. Results: Genome-wide CNVs were detected in high-density single nucleotide polymorphism (SNP) genotyping data from 1,717 Nelore (Bos indicus) cattle, and in NGS data from eight key ancestral bulls. A total of 68,007 and 12,786 distinct CNVs were observed, respectively. Cross-comparisons of results obtained for the eight resequenced animals revealed that 92 % of the CNVs were observed in both datasets, while 62 % of all detected CNVs were observed to overlap with previously validated cattle copy number variant regions (CNVRs). Observed CNVs were used for obtaining breed-specific CNV frequencies and identification of CNVRs, which were subsequently used for gene annotation. A total of 688 of the detected CNVRs were observed to overlap with 286 non-redundant QTLs associated with important production traits in cattle. All of 34 CNVs previously reported to be associated with milk production traits in Holsteins were also observed in Nelore cattle. Comparisons of estimated frequencies of these CNVs in the two breeds revealed 14, 13, 6 and 14 regions in high (>20 %), low (<20 %) and divergent (NEL > HOL, NEL < HOL) frequencies, respectively. Conclusions: Obtained results significantly enriched the bovine CNV map and enabled the identification of variants that are potentially associated with traits under selection in Nelore cattle, particularly in genome regions harboring QTLs affecting production traits.

Genomic Copy Number Alterations Of Primary And Secondary Metastasizing Pleomorphic Adenomas.

Metastasizing pleomorphic adenoma (MPA) is a rare tumour, and its mechanism of metastasis still is unknown. To date, there has been no study on MPA genomics. We analysed primary and secondary MPAs with array comparative genomic hybridization to identify somatic copy number alterations and affected genes. Tumour DNA samples from primary (parotid salivary gland) and secondary (scalp skin) MPAs were subjected to array comparative genomic hybridization investigation, and the data were analysed with NEXUS COPY NUMBER DISCOVERY. The primary MPA showed copy number losses affecting 3p22.2p14.3 and 19p13.3p123, and a complex pattern of four different deletions at chromosome 6. The 3p deletion encompassed several genes: CTNNB1, SETD2, BAP1, and PBRM1, among others. The secondary MPA showed a genomic profile similar to that of the primary MPA, with acquisition of additional copy number changes affecting 9p24.3p13.1 (loss), 19q11q13.43 (gain), and 22q11.1q13.33 (gain). Our findings indicated a clonal origin of the secondary MPA, as both tumours shared a common profile of genomic copy number alterations. Furthermore, we were able to detect in the primary tumour a specific pattern of copy number alterations that could explain the metastasizing characteristic, whereas the secondary MPA showed a more unbalanced genome.

EDD, the human orthologue of the hyperplastic discs tumour suppressor gene, is amplified and overexpressed in cancer

EDD (E3 isolated by differential display), located at chromosome 8q22.3, is the human orthologue of the Drosophila melanogaster tumour suppressor gene 'hyperplastic discs' and encodes a HECT domain E3 ubiquitin protein-ligase. To investigate the possible involvement of EDD in human cancer, several cancers from diverse tissue sites were analysed for allelic gain or loss (allelic imbalance, AI) at the EDD locus using an EDD-specific microsatellite, CEDD, and other polymorphic microsatellites mapped in the vicinity of the 8q22.3 locus. Of 143 cancers studied, 38 had AI at CEDD (42% of 90 informative cases). In 14 of these cases, discrete regions of imbalance encompassing 8q22.3 were present, while the remainder had more extensive 8q aberrations. AI of CEDD was most frequent in ovarian cancer (22/47 informative cases, 47%), particularly in the serous subtype (16/22, 73%), but was rare in benign and borderline ovarian tumours. AI was also common in breast cancer (31%), hepatocellular carcinoma (46%), squamous cell carcinoma of the tongue (50%) and metastatic melanoma (18%). AI is likely to represent amplification of the EDD gene locus rather than loss of heterozygosity, as quantitative RT-PCR and immunohistochemistry showed that EDD mRNA and protein are frequently overexpressed in breast and ovarian cancers, while among breast cancer cell lines EDD overexpression and increased gene copy number were correlated. These results demonstrate that AI at the EDD locus is common in a diversity of carcinomas and that the EDD gene is frequently overexpressed in breast and ovarian cancer, implying a potential role in cancer progression.

De novo copy number variants identify new genes and loci in isolated sporadic tetralogy of Fallot

Tetralogy of Fallot (TOF), the most common severe congenital heart malformation, occurs sporadically, without other anomaly, and from unknown cause in 70% of cases. Through a genome-wide survey of 114 subjects with TOF and their unaffected parents, we identified 11 de novo copy number variants (CNVs) that were absent or extremely rare (<0.1%) in 2,265 controls. We then examined a second, independent TOF cohort (n = 398) for additional CNVs at these loci. We identified CNVs at chromosome 1q21.1 in 1% (5/512, P = 0.0002, OR = 22.3) of nonsyndromic sporadic TOF cases. We also identified recurrent CNVs at 3p25.1, 7p21.3 and 22q11.2. CNVs in a single subject with TOF occurred at six loci, two that encode known (NOTCH1, JAG1) disease-associated genes. Our findings predict that at least 10% (4.5-15.5%, 95% confidence interval) of sporadic nonsyndromic TOF cases result from de novo CNVs and suggest that mutations within these loci might be etiologic in other cases of TOF.

Copy number variation and chromatin structure

The functional consequences of structural variation in the human genome range from adaptation, to phenotypic variation, to predisposition to diseases. Copy number variation (CNV) was shown to influence the phenotype by modifying, in a somewhat dose-dependent manner, the expression of genes that map within them, as well as that of genes located on their flanks. To assess the possible mechanism(s) behind this neighboring effect, we compared histone modification status of cell lines from patients affected by Williams-Beuren, Williams-Beuren region duplication, Smith-Magenis or DiGeorge Syndrome and control individuals using a high-throughput version of chromatin immuno-precipitation method (ChIP), called ChlP-seq. We monitored monomethylation of lysine K20 on histone H4 and trimethylation of lysine K27 on histone H3, as proxies for open and condensed chromatin, respectively. Consistent with the changes in expression levels observed for multiple genes mapping on the entire length of chromosomes affected by structural variants, we also detected regions with modified histone status between samples, up- and downstream from the critical regions, up to the end of the rearranged chromosome. We also gauged the intrachromosomal interactions of these cell lines utilizing chromosome conformation capture (4C-seq) technique. We observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together, possibly forming an interacting cluster with each other and the WBSCR. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We conclude, that large genomic rearrangements can lead to changes in the state of the chromatin spreading far away from the critical region, thus possibly affecting expression globally and as a result modifying the phenotype of the patients. - Les conséquences fonctionnelles des variations structurelles dans le génome humain sont vastes, allant de l'adaptation, en passant par les variations phénotypiques, aux prédispositions à certaines maladies. Il a été démontré que les variations du nombre de copies (CNV) influencent le phénotype en modifiant, d'une manière plus ou moins dose-dépendante, l'expression des gènes se situant à l'intérieur de ces régions, mais également celle des gènes se trouvant dans les régions flanquantes. Afin d'étudier les mécanismes possibles sous-jacents à cet effet de voisinage, nous avons comparé les états de modification des histones dans des lignées cellulaires dérivées de patients atteints du syndrome de Williams-Beuren, de la duplication de la région Williams-Beuren, du syndrome de Smith-Magenis ou du syndrome de Di- George et d'individus contrôles en utilisant une version haut-débit de la méthode d'immunoprécipitation de la chromatine (ChIP), appelée ChIP-seq. Nous avons suivi la mono-méthylation de la lysine K20 sur l'histone H4 et la tri-méthylation de la lysine K27 sur l'histone H3, marqueurs respectifs de la chromatine ouverte et fermée. En accord avec les changements de niveaux d'expression observés pour de multiples gènes tout le long des chromosomes affectés par les CNVs, nous avons aussi détecté des régions présentant des modifications d'histones entre les échantillons, situées de part et d'autre des régions critiques, jusqu'aux extrémités du chromosome réarrangé. Nous avons aussi évalué les interactions intra-chromosomiques ayant lieu dans ces cellules par l'utilisation de la technique de capture de conformation des chromosomes (4C-seq). Nous avons observé qu'un groupe de gènes flanquants la région critique du syndrome de Williams-Beuren (WBSCR) forment souvent une boucle, constituant un groupe d'interactions privilégiées entre ces gènes et la WBSCR. La délétion de la WBSCR perturbe l'expression de ce groupe de gènes flanquants, mais également les interactions à grande échelle entre eux et la région réarrangée. Nous en concluons que les larges réarrangements génomiques peuvent aboutir à des changements de l'état de la chromatine pouvant s'étendre bien plus loin que la région critique, affectant donc potentiellement l'expression de manière globale et ainsi modifiant le phénotype des patients.

The effect of translocation-induced nuclear reorganization on gene expression.

Translocations are known to affect the expression of genes at the breakpoints and, in the case of unbalanced translocations, alter the gene copy number. However, a comprehensive understanding of the functional impact of this class of variation is lacking. Here, we have studied the effect of balanced chromosomal rearrangements on gene expression by comparing the transcriptomes of cell lines from controls and individuals with the t(11;22)(q23;q11) translocation. The number of differentially expressed transcripts between translocation-carrying and control cohorts is significantly higher than that observed between control samples alone, suggesting that balanced rearrangements have a greater effect on gene expression than normal variation. Many of the affected genes are located along the length of the derived chromosome 11. We show that this chromosome is concomitantly altered in its spatial organization, occupying a more central position in the nucleus than its nonrearranged counterpart. Derivative 22-mapping chromosome 22 genes, on the other hand, remain in their usual environment. Our results are consistent with recent studies that experimentally altered nuclear organization, and indicated that nuclear position plays a functional role in regulating the expression of some genes in mammalian cells. Our study suggests that chromosomal translocations can result in hitherto unforeseen, large-scale changes in gene expression that are the consequence of alterations in normal chromosome territory positioning. This has consequences for the patterns of gene expression change seen during tumorigenesis-associated genome instability and during the karyotype changes that lead to speciation.

Prospective analysis of KRAS, BRAF, and PIK3CA mutations and EGFR copy number in patients (pts) with locally advanced rectal cancer: A translational substudy of a clinical trial (SAKK 41/07) evaluating the effect of neoadjuvant chemoradiation (CRT) with or without panitumumab

Background: Prognostic and predictive markers are of great importance for future study designs and essential for the interpretation of clinical trials incorporating an EGFR-inhibitor. The current study prospectively assessed and validated KRAS, BRAF and PIK3CA mutations in rectal cancer patients screened for the trial SAKK41/07 of concomitant preoperative radio-chemotherapy with or without panitumumab.Methods: Macrodissection was performed on pretreatment formalin fixed paraffin embedded biopsy tissue sections to arrive at a minimum of 50% of tumor cells. DNA was extracted with the Maxwell 16 FFPE Tissue LEV DNA purification kit. After PCR amplification, mutations were identified by pyrosequencing. We prospectively analysed pretreatment biopsy material from 149 rectal cancer pts biopsies for KRAS (exon 2 codon 12 [2-12] and 13 [2-13], exon 3 codon 59 [3-59]) and 61 [3-61], exon 4 codon 117 [4-117] and 146 [4-146]). Sixty-eight pts (KRASwt exon 2, 3 only) were further analysed for BRAF (exon 15 codon 600) and PIK3CA (exon 9 codon 542, 545 and 546, exon 20 codon 1043 [20-1043] and 1047 [20-1047]) mutations, and EGFR copy number by qPCR. For the calculation of the EGFR copy number, we used KRAS copy number as internal reference standard. The calculation was done on the basis of the two standard curves relative quantification method.Results: In 149 screened pts with rectal cancer, the prevalence of KRAS mutations was 36%. Among the 68 pts enrolled in SAKK 41/07 based on initially presumed KRASwt status (exon 2/codons 12+13), 18 pts (26%) had a total of 23 mutations in the RAS/PIK3CA-pathways upon validation analysis. Twelve pts had a KRAS mutation, 7 pts had a PIK3CA mutation, 3 pts had a NRAS mutation, 1 patient a BRAF mutation. Surprisingly, five of these pts had double- mutations, including 4 pts with KRAS plus PIK3CA mutations, and 1 pt with NRAS plus PIK3CA mutations. The median normalized EGFR copy number was 1. Neither mutations of KRAS, BRAF, and PIK3CA, nor EGFR copy number were statistically associated with the primary study endpoint pCR (pathological complete regression).Conclusions: The prevalence of KRAS mutations in rectal and in colon cancer appears to be similar. BRAF mutations are rare; PIK3CA mutations are more common (10%). EGFR copy number is not increased in rectal cancer. A considerable number or KRAS exon 2 wt tumors harbored KRAS exon 3+4 mutations. Further study is needed to determine if KRAS testing should include exons 2-4.

X-chromosome tiling path array detection of copy number variants in patients with chromosome X-linked mental retardation

Background: Aproximately 5–10% of cases of mental retardation in males are due to copy number variations (CNV) on the X chromosome. Novel technologies, such as array comparative genomic hybridization (aCGH), may help to uncover cryptic rearrangements in X-linked mental retardation (XLMR) patients. We have constructed an X-chromosome tiling path array using bacterial artificial chromosomes (BACs) and validated it using samples with cytogenetically defined copy number changes. We have studied 54 patients with idiopathic mental retardation and 20 controls subjects. Results: Known genomic aberrations were reliably detected on the array and eight novel submicroscopic imbalances, likely causative for the mental retardation (MR) phenotype, were detected. Putatively pathogenic rearrangements included three deletions and five duplications (ranging between 82 kb to one Mb), all but two affecting genes previously known to be responsible for XLMR. Additionally, we describe different CNV regions with significant different frequencies in XLMR and control subjects (44% vs. 20%). Conclusion:This tiling path array of the human X chromosome has proven successful for the detection and characterization of known rearrangements and novel CNVs in XLMR patients.

Genome-wide association analysis of copy number variation in recurrent depressive disorder.

Large, rare copy number variants (CNVs) have been implicated in a variety of psychiatric disorders, but the role of CNVs in recurrent depression is unclear. We performed a genome-wide analysis of large, rare CNVs in 3106 cases of recurrent depression, 459 controls screened for lifetime-absence of psychiatric disorder and 5619 unscreened controls from phase 2 of the Wellcome Trust Case Control Consortium (WTCCC2). We compared the frequency of cases with CNVs against the frequency observed in each control group, analysing CNVs over the whole genome, genic, intergenic, intronic and exonic regions. We found that deletion CNVs were associated with recurrent depression, whereas duplications were not. The effect was significant when comparing cases with WTCCC2 controls (P=7.7 × 10(-6), odds ratio (OR) =1.25 (95% confidence interval (CI) 1.13-1.37)) and to screened controls (P=5.6 × 10(-4), OR=1.52 (95% CI 1.20-1.93). Further analysis showed that CNVs deleting protein coding regions were largely responsible for the association. Within an analysis of regions previously implicated in schizophrenia, we found an overall enrichment of CNVs in our cases when compared with screened controls (P=0.019). We observe an ordered increase of samples with deletion CNVs, with the lowest proportion seen in screened controls, the next highest in unscreened controls and the highest in cases. This may suggest that the absence of deletion CNVs, especially in genes, is associated with resilience to recurrent depression.