Mutation spectrum of Drosophila CNVs revealed by breakpoint sequencing

Background

The detailed study of breakpoints associated with copy number variants (CNVs) can elucidate the mutational mechanisms that generate them and the comparison of breakpoints across species can highlight differences in genomic architecture that may lead to lineage-specific differences in patterns of CNVs. Here, we provide a detailed analysis of Drosophila CNV breakpoints and contrast it with similar analyses recently carried out for the human genome.

Results

By applying split-read methods to a total of 10x coverage of 454 shotgun sequence across nine lines of D. melanogaster and by re-examining a previously published dataset of CNVs detected using tiling arrays, we identified the precise breakpoints of more than 600 insertions, deletions, and duplications. Contrasting these CNVs with those found in humans showed that in both taxa CNV breakpoints fall into three classes: blunt breakpoints; simple breakpoints associated with microhomology; and breakpoints with additional nucleotides inserted/deleted and no microhomology. In both taxa CNV breakpoints are enriched with non-B DNA sequence structures, which may impair DNA replication and/or repair. However, in contrast to human genomes, non-allelic homologous-recombination (NAHR) plays a negligible role in CNV formation in Drosophila. In flies, non-homologous repair mechanisms are responsible for simple, recurrent, and complex CNVs, including insertions of de novo sequence as large as 60 bp.

Conclusions

Humans and Drosophila differ considerably in the importance of homology-based mechanisms for the formation of CNVs, likely as a consequence of the differences in the abundance and distribution of both segmental duplications and transposable elements between the two genomes.     

Integration of SNP and mRNA Arrays with MicroRNA Profiling Reveals That MiR-370 Is Upregulated and Targets NF1 in Acute Myeloid Leukemia

Background

Deregulated miRNA expression plays a crucial role in carcinogenesis. Recent studies show different mechanisms leading to miRNA deregulation in cancer; however, alterations affecting miRNAs by DNA copy number variations (CNV) remain poorly studied.

Results

Our integrative analysis including data from high resolution SNPs arrays, mRNA expression arrays, and miRNAs expression profiles in 16 myeloid cell lines highlights that CNV are alternative mechanisms to deregulate the expression of miRNAs in acute myeloid leukemia (AML), and represent a novel approach to identify novel candidate genes involved in AML. We found association between the expression levels of 19 miRNAs and CNVs affecting their loci. Functional analysis showed that NF1 is a direct target of miR-370, and that overexpression of miR-370 has similar effects that NF1 inactivation, increasing proliferation and colony formation in AML cells. Moreover, real time RT-PCR showed that NF1 downregulation is a recurrent event in AML (30.8%), and western blot analysis confirmed this result. MiR-370 overexpression and deletions affecting the NF1 locus were identified as alternative mechanisms to downregulate NF1.

Conclusions

NF1 downregulation is a common event in AML, and both deletions in the NF1 locus and overexpression of miR-370 are alternative mechanisms to downregulate NF1 in this disease. Our results suggest a leukemogenic role of miR-370 through NF1 downregulation in AML cells. Since NF1 deficiency leads to RAS activation, patients with AML and overexpression of miR-370 may potentially benefit from additional treatment with either RAS or mTOR inhibitors.     

DNA copy number evolution in Drosophila cell lines

Background

Structural rearrangements of the genome resulting in genic imbalance due to copy number change are often deleterious at the organismal level, but are common in immortalized cell lines and tumors, where they may be an advantage to cells. In order to explore the biological consequences of copy number changes in the Drosophila genome, we resequenced the genomes of 19 tissue-culture cell lines and generated RNA-Seq profiles.

Results

Our work revealed dramatic duplications and deletions in all cell lines. We found three lines of evidence indicating that copy number changes were due to selection during tissue culture. First, we found that copy numbers correlated to maintain stoichiometric balance in protein complexes and biochemical pathways, consistent with the gene balance hypothesis. Second, while most copy number changes were cell line-specific, we identified some copy number changes shared by many of the independent cell lines. These included dramatic recurrence of increased copy number of the PDGF/VEGF receptor, which is also over-expressed in many cancer cells, and of bantam, an anti-apoptosis miRNA. Third, even when copy number changes seemed distinct between lines, there was strong evidence that they supported a common phenotypic outcome. For example, we found that proto-oncogenes were over-represented in one cell line (S2-DRSC), whereas tumor suppressor genes were under-represented in another (Kc167).

Conclusion

Our study illustrates how genome structure changes may contribute to selection of cell lines in vitro. This has implications for other cell-level natural selection progressions, including tumorigenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/gb-2014-15-8-r70) contains supplementary material, which is available to authorized users.     

Inferring mechanisms of copy number change from haplotype structures at the human DEFA1A3 locus

Background

The determination of structural haplotypes at copy number variable regions can indicate the mechanisms responsible for changes in copy number, as well as explain the relationship between gene copy number and expression. However, obtaining spatial information at regions displaying extensive copy number variation, such as the DEFA1A3 locus, is complex, because of the difficulty in the phasing and assembly of these regions. The DEFA1A3 locus is intriguing in that it falls within a region of high linkage disequilibrium, despite its high variability in copy number (n = 3–16); hence, the mechanisms responsible for changes in copy number at this locus are unclear.

Results

In this study, a region flanking the DEFA1A3 locus was sequenced across 120 independent haplotypes with European ancestry, identifying five common classes of DEFA1A3 haplotype. Assigning DEFA1A3 class to haplotypes within the 1000 Genomes project highlights a significant difference in DEFA1A3 class frequencies between populations with different ancestry. The features of each DEFA1A3 class, for example, the associated DEFA1A3 copy numbers, were initially assessed in a European cohort (n = 599) and replicated in the 1000 Genomes samples, showing within-class similarity, but between-class and between-population differences in the features of the DEFA1A3 locus. Emulsion haplotype fusion-PCR was used to generate 61 structural haplotypes at the DEFA1A3 locus, showing a high within-class similarity in structure.

Conclusions

Structural haplotypes across the DEFA1A3 locus indicate that intra-allelic rearrangement is the predominant mechanism responsible for changes in DEFA1A3 copy number, explaining the conservation of linkage disequilibrium across the locus. The identification of common structural haplotypes at the DEFA1A3 locus could aid studies into how DEFA1A3 copy number influences expression, which is currently unclear.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-614) contains supplementary material, which is available to authorized users.     

Lineage specific evolution of the VNTR composite retrotransposon central domain and its role in retrotransposition of gibbon LAVA elements

Background

VNTR (Variable Number of Tandem Repeats) composite retrotransposons - SVA (SINE-R-VNTR-Alu), LAVA (LINE-1-Alu-VNTR-Alu), PVA (PTGR2-VNTR-Alu) and FVA (FRAM-VNTR-Alu) - are specific to hominoid primates. Their assembly, the evolution of their 5’ and 3’ domains, and the functional significance of the shared 5’ Alu-like region are well understood. The central VNTR domain, by contrast, has long been assumed to represent a more or less random collection of 30-50 bp GC-rich repeats. It is only recently that it attracted attention in the context of regulation of SVA expression.

Results

Here we provide evidence that the organization of the VNTR is non-random, with conserved repeat unit (RU) arrays at both the 5’ and 3’ ends of the VNTRs of human, chimpanzee and orangutan SVA and gibbon LAVA. The younger SVA subfamilies harbour highly organized internal RU arrays. The composition of these arrays is specific to the human/chimpanzee and orangutan lineages, respectively. Tracing the development of the VNTR through evolution we show for the first time how tandem repeats evolve within the constraints set by a functional, non-autonomous non-LTR retrotransposon in two different families - LAVA and SVA - in different hominoid lineages. Our analysis revealed that a microhomology-driven mechanism mediates expansion/contraction of the VNTR domain at the DNA level.Elements of all four VNTR composite families have been shown to be mobilized by the autonomous LINE1 retrotransposon in trans. In case of SVA, key determinants of mobilization are found in the 5’ hexameric repeat/Alu-like region. We now demonstrate that in LAVA, by contrast, the VNTR domain determines mobilization efficiency in the context of domain swaps between active and inactive elements.

Conclusions

The central domain of VNTR composites evolves in a lineage-specific manner which gives rise to distinct structures in gibbon LAVA, orangutan SVA, and human/chimpanzee SVA. The differences observed between the families and lineages are likely to have an influence on the expression and mobilization of the elements.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1543-z) contains supplementary material, which is available to authorized users.     

Guanine quadruplexes are formed by specific regions of human transposable elements

Background

Transposable elements form a significant proportion of eukaryotic genomes. Recently, Lexa et al. (Nucleic Acids Res 42:968-978, 2014) reported that plant long terminal repeat (LTR) retrotransposons often contain potential quadruplex sequences (PQSs) in their LTRs and experimentally confirmed their ability to adopt four-stranded DNA conformations.

Results

Here, we searched for PQSs in human retrotransposons and found that PQSs are specifically localized in the 3’-UTR of LINE-1 elements, in LTRs of HERV elements and are strongly accumulated in specific regions of SVA elements. Circular dichroism spectroscopy confirmed that most PQSs had adopted monomolecular or bimolecular guanine quadruplex structures. Evolutionarily young SVA elements contained more PQSs than older elements and their propensity to form quadruplex DNA was higher. Full-length L1 elements contained more PQSs than truncated elements; the highest proportion of PQSs was found inside transpositionally active L1 elements (PA2 and HS families).

Conclusions

Conservation of quadruplexes at specific positions of transposable elements implies their importance in their life cycle. The increasing quadruplex presence in evolutionarily young LINE-1 and SVA families makes these elements important contributors toward present genome-wide quadruplex distribution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1032) contains supplementary material, which is available to authorized users.     

Distinct Genetic Alterations in Colorectal Cancer

Background

Colon cancer (CRC) development often includes chromosomal instability (CIN) leading to amplifications and deletions of large DNA segments. Epidemiological, clinical, and cytogenetic studies showed that there are considerable differences between CRC tumors from African Americans (AAs) and Caucasian patients. In this study, we determined genomic copy number aberrations in sporadic CRC tumors from AAs, in order to investigate possible explanations for the observed disparities.

Methodology/Principal Findings

We applied genome-wide array comparative genome hybridization (aCGH) using a 105k chip to identify copy number aberrations in samples from 15 AAs. In addition, we did a population comparative analysis with aCGH data in Caucasians as well as with a widely publicized list of colon cancer genes (CAN genes). There was an average of 20 aberrations per patient with more amplifications than deletions. Analysis of DNA copy number of frequently altered chromosomes revealed that deletions occurred primarily in chromosomes 4, 8 and 18. Chromosomal duplications occurred in more than 50% of cases on chromosomes 7, 8, 13, 20 and X. The CIN profile showed some differences when compared to Caucasian alterations.

Conclusions/Significance

Chromosome X amplification in male patients and chromosomes 4, 8 and 18 deletions were prominent aberrations in AAs. Some CAN genes were altered at high frequencies in AAs with EXOC4, EPHB6, GNAS, MLL3 and TBX22 as the most frequently deleted genes and HAPLN1, ADAM29, SMAD2 and SMAD4 as the most frequently amplified genes. The observed CIN may play a distinctive role in CRC in AAs.     

seqCNA: an R package for DNA copy number analysis in cancer using high-throughput sequencing

Background

Deviations in the amount of genomic content that arise during tumorigenesis, called copy number alterations, are structural rearrangements that can critically affect gene expression patterns. Additionally, copy number alteration profiles allow insight into cancer discrimination, progression and complexity. On data obtained from high-throughput sequencing, improving quality through GC bias correction and keeping false positives to a minimum help build reliable copy number alteration profiles.

Results

We introduce seqCNA, a parallelized R package for an integral copy number analysis of high-throughput sequencing cancer data. The package includes novel methodology on (i) filtering, reducing false positives, and (ii) GC content correction, improving copy number profile quality, especially under great read coverage and high correlation between GC content and copy number. Adequate analysis steps are automatically chosen based on availability of paired-end mapping, matched normal samples and genome annotation.

Conclusions

seqCNA, available through Bioconductor, provides accurate copy number predictions in tumoural data, thanks to the extensive filtering and better GC bias correction, while providing an integrated and parallelized workflow.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-178) contains supplementary material, which is available to authorized users.     

Association of G-quadruplex forming sequences with human mtDNA deletion breakpoints

Background

Mitochondrial DNA (mtDNA) deletions cause disease and accumulate during aging, yet our understanding of the molecular mechanisms underlying their formation remains rudimentary. Guanine-quadruplex (GQ) DNA structures are associated with nuclear DNA instability in cancer; recent evidence indicates they can also form in mitochondrial nucleic acids, suggesting that these non-B DNA structures could be associated with mtDNA deletions. Currently, the multiple types of GQ sequences and their association with human mtDNA stability are unknown.

Results

Here, we show an association between human mtDNA deletion breakpoint locations (sites where DNA ends rejoin after deletion of a section) and sequences with G-quadruplex forming potential (QFP), and establish the ability of selected sequences to form GQ in vitro. QFP contain four runs of either two or three consecutive guanines (2G and 3G, respectively), and we identified four types of QFP for subsequent analysis: intrastrand 2G, intrastrand 3G, duplex derived interstrand (ddi) 2G, and ddi 3G QFP sequences. We analyzed the position of each motif set relative to either 5'' or 3'' unique mtDNA deletion breakpoints, and found that intrastrand QFP sequences, but not ddi QFP sequences, showed significant association with mtDNA deletion breakpoint locations. Moreover, a large proportion of these QFP sequences occur at smaller distances to breakpoints relative to distribution-matched controls. The positive association of 2G QFP sequences persisted when breakpoints were divided into clinical subgroups. We tested in vitro GQ formation of representative mtDNA sequences containing these 2G QFP sequences and detected robust GQ structures by UV–VIS and CD spectroscopy. Notably, the most frequent deletion breakpoints, including those of the "common deletion", are bounded by 2G QFP sequence motifs.

Conclusions

The potential for GQ to influence mitochondrial genome stability supports a high-priority investigation of these structures and their regulation in normal and pathological mitochondrial biology. These findings emphasize the potential importance of helicases that subsequently resolve GQ to maintain the stability of the mitochondrial genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-677) contains supplementary material, which is available to authorized users.     

Individual and combined effects of DNA methylation and copy number alterations on miRNA expression in breast tumors

Background

The global effect of copy number and epigenetic alterations on miRNA expression in cancer is poorly understood. In the present study, we integrate genome-wide DNA methylation, copy number and miRNA expression and identify genetic mechanisms underlying miRNA dysregulation in breast cancer.

Results

We identify 70 miRNAs whose expression was associated with alterations in copy number or methylation, or both. Among these, five miRNA families are represented. Interestingly, the members of these families are encoded on different chromosomes and are complementarily altered by gain or hypomethylation across the patients. In an independent breast cancer cohort of 123 patients, 41 of the 70 miRNAs were confirmed with respect to aberration pattern and association to expression. In vitro functional experiments were performed in breast cancer cell lines with miRNA mimics to evaluate the phenotype of the replicated miRNAs. let-7e-3p, which in tumors is found associated with hypermethylation, is shown to induce apoptosis and reduce cell viability, and low let-7e-3p expression is associated with poorer prognosis. The overexpression of three other miRNAs associated with copy number gain, miR-21-3p, miR-148b-3p and miR-151a-5p, increases proliferation of breast cancer cell lines. In addition, miR-151a-5p enhances the levels of phosphorylated AKT protein.

Conclusions

Our data provide novel evidence of the mechanisms behind miRNA dysregulation in breast cancer. The study contributes to the understanding of how methylation and copy number alterations influence miRNA expression, emphasizing miRNA functionality through redundant encoding, and suggests novel miRNAs important in breast cancer.     

Reference genome of wild goat (capra aegagrus) and sequencing of goat breeds provide insight into genic basis of goat domestication

Background

Domestic goats (Capra hircus) have been selected to play an essential role in agricultural production systems, since being domesticated from their wild progenitor, bezoar (Capra aegagrus). A detailed understanding of the genetic consequences imparted by the domestication process remains a key goal of evolutionary genomics.

Results

We constructed the reference genome of bezoar and sequenced representative breeds of domestic goats to search for genomic changes that likely have accompanied goat domestication and breed formation. Thirteen copy number variation genes associated with coat color were identified in domestic goats, among which ASIP gene duplication contributes to the generation of light coat-color phenotype in domestic goats. Analysis of rapidly evolving genes identified genic changes underlying behavior-related traits, immune response and production-related traits.

Conclusion

Based on the comparison studies of copy number variation genes and rapidly evolving genes between wild and domestic goat, our findings and methodology shed light on the genetic mechanism of animal domestication and will facilitate future goat breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1606-1) contains supplementary material, which is available to authorized users.     

Short Copy Number Variations Potentially Associated with Tonic Immobility Responses in Newly Hatched Chicks

Introduction

Tonic immobility (TI) is fear-induced freezing that animals may undergo when confronted by a threat. It is principally observed in prey species as defence mechanisms. In our preliminary research, we detected large inter-individual variations in the frequency and duration of freezing behavior among newly hatched domestic chicks (Gallus gallus). In this study we aim to identify the copy number variations (CNVs) in the genome of chicks as genetic candidates that underlie the behavioral plasticity to fearful stimuli.

Methods

A total of 110 domestic chicks were used for an association study between TI responses and copy number polymorphisms. Array comparative genomic hybridization (aCGH) was conducted between chicks with high and low TI scores using an Agilent 4×180 custom microarray. We specifically focused on 3 genomic regions (>60 Mb) of chromosome 1 where previous quantitative trait loci (QTL) analysis showed significant F-values for fearful responses.

Results

ACGH successfully detected short CNVs within the regions overlapping 3 QTL peaks. Eleven of these identified loci were validated by real-time quantitative polymerase chain reaction (qPCR) as copy number polymorphisms. Although there wkas no significant p value in the correlation analysis between TI scores and the relative copy number within each breed, several CNV loci showed significant differences in the relative copy number between 2 breeds of chicken (White Leghorn and Nagoya) which had different quantitative characteristics of fear-induced responses.

Conclusion

Our data shows the potential CNVs that may be responsible for innate fear response in domestic chicks.     

PacBio-LITS: a large-insert targeted sequencing method for characterization of human disease-associated chromosomal structural variations

Background

Generation of long (>5 Kb) DNA sequencing reads provides an approach for interrogation of complex regions in the human genome. Currently, large-insert whole genome sequencing (WGS) technologies from Pacific Biosciences (PacBio) enable analysis of chromosomal structural variations (SVs), but the cost to achieve the required sequence coverage across the entire human genome is high.

Results

We developed a method (termed PacBio-LITS) that combines oligonucleotide-based DNA target-capture enrichment technologies with PacBio large-insert library preparation to facilitate SV studies at specific chromosomal regions. PacBio-LITS provides deep sequence coverage at the specified sites at substantially reduced cost compared with PacBio WGS. The efficacy of PacBio-LITS is illustrated by delineating the breakpoint junctions of low copy repeat (LCR)-associated complex structural rearrangements on chr17p11.2 in patients diagnosed with Potocki–Lupski syndrome (PTLS; MIM#610883). We successfully identified previously determined breakpoint junctions in three PTLS cases, and also were able to discover novel junctions in repetitive sequences, including LCR-mediated breakpoints. The new information has enabled us to propose mechanisms for formation of these structural variants.

Conclusions

The new method leverages the cost efficiency of targeted capture-sequencing as well as the mappability and scaffolding capabilities of long sequencing reads generated by the PacBio platform. It is therefore suitable for studying complex SVs, especially those involving LCRs, inversions, and the generation of chimeric Alu elements at the breakpoints. Other genomic research applications, such as haplotype phasing and small insertion and deletion validation could also benefit from this technology.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1370-2) contains supplementary material, which is available to authorized users.     

Elusive Copy Number Variation in the Mouse Genome

Background

Array comparative genomic hybridization (aCGH) to detect copy number variants (CNVs) in mammalian genomes has led to a growing awareness of the potential importance of this category of sequence variation as a cause of phenotypic variation. Yet there are large discrepancies between studies, so that the extent of the genome affected by CNVs is unknown. We combined molecular and aCGH analyses of CNVs in inbred mouse strains to investigate this question.

Principal Findings

Using a 2.1 million probe array we identified 1,477 deletions and 499 gains in 7 inbred mouse strains. Molecular characterization indicated that approximately one third of the CNVs detected by the array were false positives and we estimate the false negative rate to be more than 50%. We show that low concordance between studies is largely due to the molecular nature of CNVs, many of which consist of a series of smaller deletions and gains interspersed by regions where the DNA copy number is normal.

Conclusions

Our results indicate that CNVs detected by arrays may be the coincidental co-localization of smaller CNVs, whose presence is more likely to perturb an aCGH hybridization profile than the effect of an isolated, small, copy number alteration. Our findings help explain the hitherto unexplored discrepancies between array-based studies of copy number variation in the mouse genome.     

Genome Rearrangements Detected by SNP Microarrays in Individuals with Intellectual Disability Referred with Possible Williams Syndrome

Background

Intellectual disability (ID) affects 2–3% of the population and may occur with or without multiple congenital anomalies (MCA) or other medical conditions. Established genetic syndromes and visible chromosome abnormalities account for a substantial percentage of ID diagnoses, although for ∼50% the molecular etiology is unknown. Individuals with features suggestive of various syndromes but lacking their associated genetic anomalies pose a formidable clinical challenge. With the advent of microarray techniques, submicroscopic genome alterations not associated with known syndromes are emerging as a significant cause of ID and MCA.

Methodology/Principal Findings

High-density SNP microarrays were used to determine genome wide copy number in 42 individuals: 7 with confirmed alterations in the WS region but atypical clinical phenotypes, 31 with ID and/or MCA, and 4 controls. One individual from the first group had the most telomeric gene in the WS critical region deleted along with 2 Mb of flanking sequence. A second person had the classic WS deletion and a rearrangement on chromosome 5p within the Cri du Chat syndrome (OMIM:123450) region. Six individuals from the ID/MCA group had large rearrangements (3 deletions, 3 duplications), one of whom had a large inversion associated with a deletion that was not detected by the SNP arrays.

Conclusions/Significance

Combining SNP microarray analyses and qPCR allowed us to clone and sequence 21 deletion breakpoints in individuals with atypical deletions in the WS region and/or ID or MCA. Comparison of these breakpoints to databases of genomic variation revealed that 52% occurred in regions harboring structural variants in the general population. For two probands the genomic alterations were flanked by segmental duplications, which frequently mediate recurrent genome rearrangements; these may represent new genomic disorders. While SNP arrays and related technologies can identify potentially pathogenic deletions and duplications, obtaining sequence information from the breakpoints frequently provides additional information.     

Amplified and Homozygously Deleted Genes in Glioblastoma: Impact on Gene Expression Levels

Background

Glioblastoma multiforme (GBM) displays multiple amplicons and homozygous deletions that involve relevant pathogenic genes and other genes whose role remains unknown.

Methodology

Single-nucleotide polymorphism (SNP)-arrays were used to determine the frequency of recurrent amplicons and homozygous deletions in GBM (n = 46), and to evaluate the impact of copy number alterations (CNA) on mRNA levels of the genes involved.

Principal Findings

Recurrent amplicons were detected for chromosomes 7 (50%), 12 (22%), 1 (11%), 4 (9%), 11 (4%), and 17 (4%), whereas homozygous deletions involved chromosomes 9p21 (52%) and 10q (22%). Most genes that displayed a high correlation between DNA CNA and mRNA levels were coded in the amplified chromosomes. For some amplicons the impact of DNA CNA on mRNA expression was restricted to a single gene (e.g., EGFR at 7p11.2), while for others it involved multiple genes (e.g., 11 and 5 genes at 12q14.1–q15 and 4q12, respectively). Despite homozygous del(9p21) and del(10q23.31) included multiple genes, association between these DNA CNA and RNA expression was restricted to the MTAP gene.

Conclusions

Overall, our results showed a high frequency of amplicons and homozygous deletions in GBM with variable impact on the expression of the genes involved, and they contributed to the identification of other potentially relevant genes.     

Micro-Scale Genomic DNA Copy Number Aberrations as Another Means of Mutagenesis in Breast Cancer

Introduction

In breast cancer, the basal-like subtype has high levels of genomic instability relative to other breast cancer subtypes with many basal-like-specific regions of aberration. There is evidence that this genomic instability extends to smaller scale genomic aberrations, as shown by a previously described micro-deletion event in the PTEN gene in the Basal-like SUM149 breast cancer cell line.

Methods

We sought to identify if small regions of genomic DNA copy number changes exist by using a high density, gene-centric Comparative Genomic Hybridizations (CGH) array on cell lines and primary tumors. A custom tiling array for CGH (244,000 probes, 200 bp tiling resolution) was created to identify small regions of genomic change, which was focused on previously identified basal-like-specific, and general cancer genes. Tumor genomic DNA from 94 patients and 2 breast cancer cell lines was labeled and hybridized to these arrays. Aberrations were called using SWITCHdna and the smallest 25% of SWITCHdna-defined genomic segments were called micro-aberrations (<64 contiguous probes, ∼ 15 kb).

Results

Our data showed that primary tumor breast cancer genomes frequently contained many small-scale copy number gains and losses, termed micro-aberrations, most of which are undetectable using typical-density genome-wide aCGH arrays. The basal-like subtype exhibited the highest incidence of these events. These micro-aberrations sometimes altered expression of the involved gene. We confirmed the presence of the PTEN micro-amplification in SUM149 and by mRNA-seq showed that this resulted in loss of expression of all exons downstream of this event. Micro-aberrations disproportionately affected the 5′ regions of the affected genes, including the promoter region, and high frequency of micro-aberrations was associated with poor survival.

Conclusion

Using a high-probe-density, gene-centric aCGH microarray, we present evidence of small-scale genomic aberrations that can contribute to gene inactivation. These events may contribute to tumor formation through mechanisms not detected using conventional DNA copy number analyses.     

Array-comparative genomic hybridization reveals loss of SOCS6 is associated with poor prognosis in primary lung squamous cell carcinoma

Background

Primary tumor recurrence commonly occurs after surgical resection of lung squamous cell carcinoma (SCC). Little is known about the genes driving SCC recurrence.

Methods

We used array comparative genomic hybridization (aCGH) to identify genes affected by copy number alterations that may be involved in SCC recurrence. Training and test sets of resected primary lung SCC were assembled. aCGH was used to determine genomic copy number in a training set of 62 primary lung SCCs (28 with recurrence and 34 with no evidence of recurrence) and the altered copy number of candidate genes was confirmed by quantitative PCR (qPCR). An independent test set of 72 primary lung SCCs (20 with recurrence and 52 with no evidence of recurrence) was used for biological validation. mRNA expression of candidate genes was studied using qRT-PCR. Candidate gene promoter methylation was evaluated using methylation microarrays and Sequenom EpiTYPER analysis.

Results

18q22.3 loss was identified by aCGH as being significantly associated with recurrence (p = 0.038). Seven genes within 18q22.3 had aCGH copy number loss associated with recurrence but only SOCS6 copy number was both technically replicated by qPCR and biologically validated in the test set. SOCS6 copy number loss correlated with reduced mRNA expression in the study samples and in the samples with copy number loss, there was a trend for increased methylation, albeit non-significant. Overall survival was significantly poorer in patients with SOCS6 loss compared to patients without SOCS6 loss in both the training (30 vs. 43 months, p = 0.023) and test set (27 vs. 43 months, p = 0.010).

Conclusion

Reduced copy number and mRNA expression of SOCS6 are associated with disease recurrence in primary lung SCC and may be useful prognostic biomarkers.