A common copy-number breakpoint of ERBB2 amplification in breast cancer colocalizes with a complex block of segmental duplications

Aneuploidy and complex genomic rearrangements in cancer evolution

Mutational processes that alter large genomic regions occur frequently in developing tumors. They range from simple copy number gains and losses to the shattering and reassembly of entire chromosomes. These catastrophic events, such as chromothripsis, chromoplexy and the formation of extrachromosomal DNA, affect the expression of many genes and therefore have a substantial effect on the fitness of the cells in which they arise. In this review, we cover large genomic alterations, the mechanisms that cause them and their effect on tumor development and evolution.

Double-strand breaks induce inverted duplication chromosome rearrangements by a DNA polymerase δ-dependent mechanism

Inverted duplications, also known as foldback inversions, are commonly observed in cancers and are the major class of chromosome rearrangement recovered from yeast cells lacking Mre11 nuclease activity. Foldback priming at DNA double-strand breaks (DSBs) is one mechanism proposed for the generation of inverted duplications. However, the other pathway steps have not been fully elucidated. Here, we show that a DSB induced near natural inverted repeats drives high frequency inverted duplication in Sae2 and Mre11-deficient cells. We find that DNA polymerase δ proof-reading activity, but not Rad1 nuclease, trims the heterologous flaps formed after foldback annealing. Additionally, Pol32 is required for the generation of inverted duplications, suggesting that Pol δ catalyzes fill-in synthesis primed from the foldback to create a hairpin-capped chromosome that is subsequently replicated to form a dicentric inversion chromosome. Finally, we show that stabilization of the dicentric chromosome after breakage involves telomere capture by non-reciprocal translocation mediated by repeat sequences or by deletion of one centromere.

Double-strand breaks induce inverted duplication chromosome rearrangements by a DNA polymerase δ and Rad51-dependent mechanism

Inverted duplications, also known as foldback inversions, are commonly observed in cancers and are the major class of chromosome rearrangement recovered from yeast cells lacking Mre11 nuclease. Foldback priming at naturally occurring inverted repeats is one mechanism proposed for the generation of inverted duplications. However, the initiating lesion for these events and the mechanism by which they form has not been fully elucidated. Here, we show that a DNA double-strand break (DSB) induced near natural short, inverted repeats drives high frequency inverted duplication in Sae2 and Mre11-deficient cells. We find that DNA polymerase δ proof-reading activity acts non-redundantly with Rad1 nuclease to remove heterologous tails formed during foldback annealing. Additionally, Pol32 is required for the generation of inverted duplications, suggesting that Pol δ catalyzes fill-in synthesis primed from the foldback to create a hairpin-capped chromosome that is subsequently replicated to form a dicentric isochromosome. Stabilization of the dicentric chromosome after breakage involves telomere capture by non-reciprocal translocation mediated by repeat sequences and requires Rad51.

Multi-omic cross-sectional cohort study of pre-malignant Barrett's esophagus reveals early structural variation and retrotransposon activity

Barrett's esophagus is a pre-malignant lesion that can progress to esophageal adenocarcinoma. We perform a multi-omic analysis of pre-cancer samples from 146 patients with a range of outcomes, comprising 642 person years of follow-up. Whole genome sequencing reveals complex structural variants and LINE-1 retrotransposons, as well as known copy number changes, occurring even prior to dysplasia. The structural variant burden captures the most variance across the cohort and genomic profiles do not always match consensus clinical pathology dysplasia grades. Increasing structural variant burden is associated with: high levels of chromothripsis and breakage-fusion-bridge events; increased expression of genes related to cell cycle checkpoint, DNA repair and chromosomal instability; and epigenetic silencing of Wnt signalling and cell cycle genes. Timing analysis reveals molecular events triggering genomic instability with more clonal expansion in dysplastic samples. Overall genomic complexity occurs early in the Barrett's natural history and may inform the potential for cancer beyond the clinically discernible phenotype.

Distribution of copy number variations and rearrangement endpoints in human cancers with a review of literature

Copy number variations (CNVs) which include deletions, duplications, inversions, translocations, and other forms of chromosomal re-arrangements are common to human cancers. In this report we investigated the pattern of these variations with the goal of understanding whether there exist specific cancer signatures. We used re-arrangement endpoint data deposited on the Catalogue of Somatic Mutations in Cancers (COSMIC) for our analysis. Indeed, we find that human cancers are characterized by specific patterns of chromosome rearrangements endpoints which in turn result in cancer specific CNVs. A review of the literature reveals tissue specific mutations which either drive these CNVs or appear as a consequence of CNVs because they confer an advantage to the cancer cell. We also identify several rearrangement endpoints hotspots that were not previously reported. Our analysis suggests that in addition to local chromosomal architecture, CNVs are driven by the internal cellular or nuclear physiology of each cancer tissue.

Genomic and transcriptomic analyses reveal a tandem amplification unit of 11 genes and mutations in mismatch repair genes in methotrexate-resistant HT-29 cells

DHFR gene amplification is commonly present in methotrexate (MTX)-resistant colon cancer cells and acute lymphoblastic leukemia. In this study, we proposed an integrative framework to characterize the amplified region by using a combination of single-molecule real-time sequencing, next-generation optical mapping, and chromosome conformation capture (Hi-C). We identified an amplification unit spanning 11 genes, from the DHFR gene to the ATP6AP1L gene position, with high adjusted interaction frequencies on chromosome 5 (~2.2 Mbp) and a twenty-fold tandemly amplified region, and novel inversions at the start and end positions of the amplified region as well as frameshift insertions in most of the MSH and MLH genes were detected. These mutations might stimulate chromosomal breakage and cause the dysregulation of mismatch repair. Characterizing the tandem gene-amplified unit may be critical for identifying the mechanisms that trigger genomic rearrangements. These findings may provide new insight into the mechanisms underlying the amplification process and the evolution of drug resistance.

Sequencing a large region of DNA containing many surplus copies of genes linked to drug resistance in colon cancer cells may illuminate how these genomic rearrangements arise. Such regions of gene amplification are highly repetitive, making them impossible to sequence using ordinary methods, and little is known about how they are generated. Using advanced methods, Jeong-Sun Seo at Seoul National University Bundang Hospital in South Korea and co-workers sequenced a region of gene amplification in colon cancer cells. The amplified region was approximately 20 times the length of that in healthy cells and contained many copies of an eleven-gene segment, including a gene implicated in drug resistance. The region also contained mutations in chromosomal repair genes which would disrupt repair pathways. These results illuminate the genetic changes that lead to gene amplification and drug resistance in cancer cells.

Next Generation Sequencing Technology in the Clinic and Its Challenges

Simple Summary

Precise identification and annotation of mutations are of utmost importance in clinical oncology. Insights of the DNA sequence can provide meaningful knowledge to unravel the underlying genetics of disease. Hence, tailoring of personalized medicine often relies on specific genomic alteration for treatment efficacy. The aim of this review is to highlight that sequencing harbors much more than just four nucleotides. Moreover, the gradual transition from first to second generation sequencing technologies has led to awareness for choosing the most appropriate bioinformatic analytic tools based on the aim, quality and demand for a specific purpose. Thus, the same raw data can lead to various results reflecting the intrinsic features of different datamining pipelines.

Abstract

Data analysis has become a crucial aspect in clinical oncology to interpret output from next-generation sequencing-based testing. NGS being able to resolve billions of sequencing reactions in a few days has consequently increased the demand for tools to handle and analyze such large data sets. Many tools have been developed since the advent of NGS, featuring their own peculiarities. Increased awareness when interpreting alterations in the genome is therefore of utmost importance, as the same data using different tools can provide diverse outcomes. Hence, it is crucial to evaluate and validate bioinformatic pipelines in clinical settings. Moreover, personalized medicine implies treatment targeting efficacy of biological drugs for specific genomic alterations. Here, we focused on different sequencing technologies, features underlying the genome complexity, and bioinformatic tools that can impact the final annotation. Additionally, we discuss the clinical demand and design for implementing NGS.

Functional Common and Rare ERBB2 Germline Variants Cooperate in Familial and Sporadic Cancer Susceptibility

We investigated a Spanish and Catalan family in which multiple cancer types tracked across three generations, but for which no genetic etiology had been identified. Whole-exome sequencing of germline DNA from multiple affected family members was performed to identify candidate variants to explain this occurrence of familial cancer. We discovered in all cancer-affected family members a single rare heterozygous germline variant (I654V, rs1801201) in ERBB2/HER2, which is located in a transmembrane glycine zipper motif critical for ERBB2-mediated signaling and in complete linkage disequilibrium (D' = 1) with a common polymorphism (I655V, rs1136201) previously reported in some populations as associated with cancer risk. Because multiple cancer types occurred in this family, we tested both the I654V and the I655V variants for association with cancer across multiple tumor types in 6,371 cases of Northern European ancestry drawn from The Cancer Genome Atlas and 6,647 controls, and found that the rare variant (I654V) was significantly associated with an increased risk for cancer (OR = 1.40; P = 0.021; 95% confidence interval (CI), 1.05-1.89). Functional assays performed in HEK 293T cells revealed that both the I655V single mutant (SM) and the I654V;I655V double mutant (DM) stabilized ERBB2 protein and activated ERBB2 signaling, with the DM activating ERBB2 significantly more than the SM alone. Thus, our results suggest a model whereby heritable genetic variation in the transmembrane domain activating ERBB2 signaling is associated with both sporadic and familial cancer risk, with increased ERBB2 stabilization and activation associated with increased cancer risk. PREVENTION RELEVANCE: By performing whole-exome sequencing on germline DNA from multiple cancer-affected individuals belonging to a family in which multiple cancer types track across three generations, we identified and then characterized functional common and rare variation in ERBB2 associated with both sporadic and familial cancer. Our results suggest that heritable variation activating ERBB2 signaling is associated with risk for multiple cancer types, with increases in signaling correlated with increases in risk, and modified by ancestry or family history.

Mechanisms Generating Cancer Genome Complexity: Back to the Future

Understanding the mechanisms underlying cancer genome evolution has been a major goal for decades. A recent study combining live cell imaging and single-cell genome sequencing suggested that interwoven chromosome breakage-fusion-bridge cycles, micronucleation events and chromothripsis episodes drive cancer genome evolution. Here, I discuss the "interphase breakage model," suggested from prior fluorescent in situ hybridization data that led to a similar conclusion. In this model, the rapid genome evolution observed at early stages of gene amplification was proposed to result from the interweaving of an amplification mechanism (breakage-fusion-bridge cycles) and of a deletion mechanism (micronucleation and stitching of DNA fragments retained in the nucleus).

Expandable and reversible copy number amplification drives rapid adaptation to antifungal drugs

Previously, we identified long repeat sequences that are frequently associated with genome rearrangements, including copy number variation (CNV), in many diverse isolates of the human fungal pathogen Candida albicans (Todd et al., 2019). Here, we describe the rapid acquisition of novel, high copy number CNVs during adaptation to azole antifungal drugs. Single-cell karyotype analysis indicates that these CNVs appear to arise via a dicentric chromosome intermediate and breakage-fusion-bridge cycles that are repaired using multiple distinct long inverted repeat sequences. Subsequent removal of the antifungal drug can lead to a dramatic loss of the CNV and reversion to the progenitor genotype and drug susceptibility phenotype. These findings support a novel mechanism for the rapid acquisition of antifungal drug resistance and provide genomic evidence for the heterogeneity frequently observed in clinical settings.

Comparing copy-number profiles under multi-copy amplifications and deletions

BACKGROUND

During cancer progression, malignant cells accumulate somatic mutations that can lead to genetic aberrations. In particular, evolutionary events akin to segmental duplications or deletions can alter the copy-number profile (CNP) of a set of genes in a genome. Our aim is to compute the evolutionary distance between two cells for which only CNPs are known. This asks for the minimum number of segmental amplifications and deletions to turn one CNP into another. This was recently formalized into a model where each event is assumed to alter a copy-number by 1 or -1, even though these events can affect large portions of a chromosome.

RESULTS

We propose a general cost framework where an event can modify the copy-number of a gene by larger amounts. We show that any cost scheme that allows segmental deletions of arbitrary length makes computing the distance strongly NP-hard. We then devise a factor 2 approximation algorithm for the problem when copy-numbers are non-zero and provide an implementation called cnp2cnp. We evaluate our approach experimentally by reconstructing simulated cancer phylogenies from the pairwise distances inferred by cnp2cnp and compare it against two other alternatives, namely the MEDICC distance and the Euclidean distance.

CONCLUSIONS

The experimental results show that our distance yields more accurate phylogenies on average than these alternatives if the given CNPs are error-free, but that the MEDICC distance is slightly more robust against error in the data. In all cases, our experiments show that either our approach or the MEDICC approach should preferred over the Euclidean distance.

Tetrasomy 21 pter→q21.3 due to an extra +dic(21;21)mat in a severely psychomotor-retarded female patient without Down syndrome phenotype

Complete or partial tetrasomy 21 has been reported only in rare cases. We report a Japanese female patient with tetrasomy 21 due to an extra chromosome derived from chromosome 21 (Chr21). The patient had severe psychomotor retardation without Down syndrome (DS) phenotype; she showed short stature, microcephaly, round face, cleft lip and palate, and other dysmorphic features. The chromosome analyses for the patient detected an extra dicentric Chr21 consisting of two partial Chr21 copies fused together within their long arms. Her karyotype was revealed to be 47,XX,+dic(21;21). Allelic ratios of heterozygous SNPs observed in the patient indicated the maternal origin of the extra Chr21. Copy number and structural variant analyses using whole genome sequencing data indicated that the distal breakpoint of the dicentric Chr21 was located within 21q21.3 and that the extra Chr21 did not simply consist of inverted duplications of the pter→q21.3 region, but likely contained multiple partial deletions, duplications, and inversions within it. Fluorescence in situ hybridization results were consistent with the karyotype and genomic analyses. The patient's lack of DS phenotype turned out to be due to the normal copy number of the DS critical region (21q22.13-22.3). A possible molecular mechanism leading to the complex genomic rearrangements in the tetrasomic region consists mainly of breakage-fusion-bridge cycles with an unequal crossing-over event.

Progression inference for somatic mutations in cancer

Computational methods were employed to determine progression inference of genomic alterations in commonly occurring cancers. Using cross-sectional TCGA data, we computed evolutionary trajectories involving selectivity relationships among pairs of gene-specific genomic alterations such as somatic mutations, deletions, amplifications, downregulation, and upregulation among the top 20 driver genes associated with each cancer. Results indicate that the majority of hierarchies involved TP53, PIK3CA, ERBB2, APC, KRAS, EGFR, IDH1, VHL, etc. Research into the order and accumulation of genomic alterations among cancer driver genes will ever-increase as the costs of nextgen sequencing subside, and personalized/precision medicine incorporates whole-genome scans into the diagnosis and treatment of cancer.

Palindromic amplification of the ERBB2 oncogene in primary HER2-positive breast tumors

Oncogene amplification confers a growth advantage to tumor cells for clonal expansion. There are several, recurrently amplified oncogenes throughout the human genome. However, it remains unclear whether this recurrent amplification is solely a manifestation of increased fitness resulting from random amplification mechanisms, or if a genomic locus-specific amplification mechanism plays a role. Here we show that the ERBB2 oncogene at 17q12 is susceptible to palindromic gene amplification, a mechanism characterized by the inverted (palindromic) duplication of genomic segments, in HER2-positive breast tumors. We applied two genomic approaches to investigate amplification mechanisms: sequencing of DNA libraries enriched with tumor-derived palindromic DNA (Genome-wide Analysis of Palindrome Formation) and whole genome sequencing (WGS). We observed significant enrichment of palindromic DNA within amplified ERBB2 genomic segments. Palindromic DNA was particularly enriched at amplification peaks and at boundaries between amplified and normal copy-number regions. Thus, palindromic gene amplification shaped the amplified ERBB2 locus. The enrichment of palindromic DNA throughout the amplified segments leads us to propose that the ERBB2 locus is amplified through the mechanism that repeatedly generates palindromic DNA, such as Breakage-Fusion-Bridge cycles. The genomic architecture surrounding ERBB2 in the normal genome, such as segmental duplications, could promote the locus-specific mechanism.

A whole-genome sequence and transcriptome perspective on HER2-positive breast cancers

HER2-positive breast cancer has long proven to be a clinically distinct class of breast cancers for which several targeted therapies are now available. However, resistance to the treatment associated with specific gene expressions or mutations has been observed, revealing the underlying diversity of these cancers. Therefore, understanding the full extent of the HER2-positive disease heterogeneity still remains challenging. Here we carry out an in-depth genomic characterization of 64 HER2-positive breast tumour genomes that exhibit four subgroups, based on the expression data, with distinctive genomic features in terms of somatic mutations, copy-number changes or structural variations. The results suggest that, despite being clinically defined by a specific gene amplification, HER2-positive tumours melt into the whole luminal-basal breast cancer spectrum rather than standing apart. The results also lead to a refined ERBB2 amplicon of 106 kb and show that several cases of amplifications are compatible with a breakage-fusion-bridge mechanism.

Mre11-Sae2 and RPA Collaborate to Prevent Palindromic Gene Amplification

Foldback priming at DNA double-stranded breaks is one mechanism proposed to initiate palindromic gene amplification, a common feature of cancer cells. Here, we show that small (5-9 bp) inverted repeats drive the formation of large palindromic duplications, the major class of chromosomal rearrangements recovered from yeast cells lacking Sae2 or the Mre11 nuclease. RPA dysfunction increased the frequency of palindromic duplications in Sae2 or Mre11 nuclease-deficient cells by ∼ 1,000-fold, consistent with intra-strand annealing to create a hairpin-capped chromosome that is subsequently replicated to form a dicentric isochromosome. The palindromic duplications were frequently associated with duplication of a second chromosome region bounded by a repeated sequence and a telomere, suggesting the dicentric chromosome breaks and repairs by recombination between dispersed repeats to acquire a telomere. We propose secondary structures within single-stranded DNA are potent instigators of genome instability, and RPA and Mre11-Sae2 play important roles in preventing their formation and propagation, respectively.

Mammalian keratin associated proteins (KRTAPs) subgenomes: disentangling hair diversity and adaptation to terrestrial and aquatic environments

Background

Adaptation of mammals to terrestrial life was facilitated by the unique vertebrate trait of body hair, which occurs in a range of morphological patterns. Keratin associated proteins (KRTAPs), the major structural hair shaft proteins, are largely responsible for hair variation.

Results

We exhaustively characterized the KRTAP gene family in 22 mammalian genomes, confirming the existence of 30 KRTAP subfamilies evolving at different rates with varying degrees of diversification and homogenization. Within the two major classes of KRTAPs, the high cysteine (HS) subfamily experienced strong concerted evolution, high rates of gene conversion/recombination and high GC content. In contrast, high glycine-tyrosine (HGT) KRTAPs showed evidence of positive selection and low rates of gene conversion/recombination. Species with more hair and of higher complexity tended to have more KRATP genes (gene expansion). The sloth, with long and coarse hair, had the most KRTAP genes (175 with 141 being intact). By contrast, the “hairless” dolphin had 35 KRTAPs and the highest pseudogenization rate (74% relative to the 19% mammalian average). Unique hair-related phenotypes, such as scales (armadillo) and spines (hedgehog), were correlated with changes in KRTAPs. Gene expression variation probably also influences hair diversification patterns, for example human have an identical KRTAP repertoire as apes, but much less hair.

Conclusions

We hypothesize that differences in KRTAP gene repertoire and gene expression, together with distinct rates of gene conversion/recombination, pseudogenization and positive selection, are likely responsible for micro and macro-phenotypic hair diversification among mammals in response to adaptations to ecological pressures.

Electronic supplementary material

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

Targeted next-generation sequencing at copy-number breakpoints for personalized analysis of rearranged ends in solid tumors

BACKGROUND

The concept of the utilization of rearranged ends for development of personalized biomarkers has attracted much attention owing to its clinical applicability. Although targeted next-generation sequencing (NGS) for recurrent rearrangements has been successful in hematologic malignancies, its application to solid tumors is problematic due to the paucity of recurrent translocations. However, copy-number breakpoints (CNBs), which are abundant in solid tumors, can be utilized for identification of rearranged ends.

METHOD

As a proof of concept, we performed targeted next-generation sequencing at copy-number breakpoints (TNGS-CNB) in nine colon cancer cases including seven primary cancers and two cell lines, COLO205 and SW620. For deduction of CNBs, we developed a novel competitive single-nucleotide polymorphism (cSNP) microarray method entailing CNB-region refinement by competitor DNA.

RESULT

Using TNGS-CNB, 19 specific rearrangements out of 91 CNBs (20.9%) were identified, and two polymerase chain reaction (PCR)-amplifiable rearrangements were obtained in six cases (66.7%). And significantly, TNGS-CNB, with its high positive identification rate (82.6%) of PCR-amplifiable rearrangements at candidate sites (19/23), just from filtering of aligned sequences, requires little effort for validation.

CONCLUSION

Our results indicate that TNGS-CNB, with its utility for identification of rearrangements in solid tumors, can be successfully applied in the clinical laboratory for cancer-relapse and therapy-response monitoring.

Homology-mediated end-capping as a primary step of sister chromatid fusion in the breakage-fusion-bridge cycles

Breakage-fusion-bridge (BFB) cycle is a series of chromosome breaks and duplications that could lead to the increased copy number of a genomic segment (gene amplification). A critical step of BFB cycles leading to gene amplification is a palindromic fusion of sister chromatids following the rupture of a dicentric chromosome during mitosis. It is currently unknown how sister chromatid fusion is produced from a mitotic break. To delineate the process, we took an integrated genomic, cytogenetic and molecular approach for the recurrent MCL1 amplicon at chromosome 1 in human tumor cells. A newly developed next-generation sequencing-based approach identified a cluster of palindromic fusions within the amplicon at ∼50-kb intervals, indicating a series of breaks and fusions by BFB cycles. The physical location of the amplicon (at the end of a broken chromosome) further indicated BFB cycles as underlying processes. Three palindromic fusions were mediated by the homologies between two nearby inverted Alu repeats, whereas the other two fusions exhibited microhomology-mediated events. Such breakpoint sequences indicate that homology-mediated fold-back capping of broken ends followed by DNA replication is an underlying mechanism of sister chromatid fusion. Our results elucidate nucleotide-level events during BFB cycles and end processing for naturally occurring mitotic breaks.