Alu-Alu fusion sequences identified at junction sites of copy number amplified regions in cancer cell lines

Transition from Transient DNA Rereplication to Inherited Gene Amplification Following Prolonged Environmental Stress

Cells require the ability to adapt to changing environmental conditions, however, it is unclear how these changes elicit stable permanent changes in genomes. We demonstrate that, in response to environmental metal exposure, the metallothionein (MT) locus undergoes DNA rereplication generating transient site-specific gene amplifications (TSSGs). Chronic metal exposure allows transition from MT TSSG to inherited MT gene amplification through homologous recombination within and outside of the MT locus. DNA rereplication of the MT locus is suppressed by H3K27me3 and EZH2. Long-term ablation of EZH2 activity eventually leads to integration and inheritance of MT gene amplifications without the selective pressure of metal exposure. The rereplication and inheritance of MT gene amplification is an evolutionarily conserved response to environmental metal from yeast to human. Our results describe a new paradigm for adaptation to environmental stress where targeted, transient DNA rereplication precedes stable inherited gene amplification.

Chimeric chromosome landscapes of human somatic cell cultures show dependence on stress and regulation of genomic repeats by CGGBP1

Genomes of somatic cells in culture are prone to spontaneous mutations due to errors in replication and DNA repair. Some of these errors, such as chromosomal fusions, are not rectifiable and subject to selection or elimination in growing cultures. Somatic cell cultures are thus expected to generate background levels of potentially stable chromosomal chimeras. A description of the landscape of such spontaneously generated chromosomal chimeras in cultured cells will help understand the factors affecting somatic mosaicism. Here we show that short homology-associated non-homologous chromosomal chimeras occur in normal human fibroblasts and HEK293T cells at genomic repeats. The occurrence of chromosomal chimeras is enhanced by heat stress and depletion of a repeat regulatory protein CGGBP1. We also present evidence of homologous chromosomal chimeras between allelic copies in repeat-rich DNA obtained by methylcytosine immunoprecipitation. The formation of homologous chromosomal chimeras at Alu and L1 repeats increases upon depletion of CGGBP1. Our data are derived from de novo sequencing from three different cell lines under different experimental conditions and our chromosomal chimera detection pipeline is applicable to long as well as short read sequencing platforms. These findings present significant information about the generation, sensitivity and regulation of somatic mosaicism in human cell cultures.

Structural variation and its potential impact on genome instability: Novel discoveries in the EGFR landscape by long-read sequencing

Structural variation (SV) is typically defined as variation within the human genome that exceeds 50 base pairs (bp). SV may be copy number neutral or it may involve duplications, deletions, and complex rearrangements. Recent studies have shown SV to be associated with many human diseases. However, studies of SV have been challenging due to technological constraints. With the advent of third generation (long-read) sequencing technology, exploration of longer stretches of DNA not easily examined previously has been made possible. In the present study, we utilized third generation (long-read) sequencing techniques to examine SV in the EGFR landscape of four haplotypes derived from two human samples. We analyzed the EGFR gene and its landscape (+/- 500,000 base pairs) using this approach and were able to identify a region of non-coding DNA with over 90% similarity to the most common activating EGFR mutation in non-small cell lung cancer. Based on previously published Alu-element genome instability algorithms, we propose a molecular mechanism to explain how this non-coding region of DNA may be interacting with and impacting the stability of the EGFR gene and potentially generating this cancer-driver gene. By these techniques, we were also able to identify previously hidden structural variation in the four haplotypes and in the human reference genome (hg38). We applied previously published algorithms to compare the relative stabilities of these five different EGFR gene landscape haplotypes to estimate their relative potentials to generate the EGFR exon 19, 15 bp canonical deletion. To our knowledge, the present study is the first to use the differences in genomic architecture between targeted cancer-linked phased haplotypes to estimate their relative potentials to form a common cancer-linked driver mutation.

Contrasting duplication patterns reflect functional diversities of ubiquitin and ubiquitin-like protein modifiers in plants

Ubiquitin (Ub) and Ub-like proteins, collectively forming the ubiquiton family, regulate nearly all aspects of cellular processes via post-translational modifications. Studies devoted to specific members suggested a large expansion of this family in plants; however, a lack of systematic analysis hinders the comparison of individual members at both evolutionary history and functional divergence levels, which may provide new insight into biological functions. In this work, we first retrieved a total of 5856 members of 17 known ubiquiton subfamilies in 50 plant genomes by searching both prior annotations and missing loci in each genome. We then applied this list to analyze the duplication history of major ubiquiton subfamilies in plants. We show that autophagy-related protein 8 (ATG8), membrane-anchored Ub-fold (MUB), small Ub-like modifier (SUMO) and Ub loci encode 88% of the plant ubiquiton family. Although whole genome duplications (WGDs) significantly expanded the family, we discovered contrasting duplication patterns both in species and in subfamilies. Within the family, the ATG8 and MUB members were primarily duplicated through WGDs, whereas a significant number of Ub and SUMO loci were generated through retroposition and tandem duplications, respectively. Although Ub coding regions are highly conserved in plants, promoter activity analysis demonstrated lineage-specific expression patterns of polyUb genes in Oryza sativa (rice) and Arabidopsis, confirming their retroposition origin. Based on the theory of dosage balance constraints, our study suggests that ubiquiton members duplicated through WGDs play crucial roles in plants, and that the regulatory pathways involving ATG8 and MUB are more conserved than those controlled by Ub and SUMO.

Alteration of the exDNA profile in blood serum of LLC-bearing mice under the decrease of tumour invasion potential by bovine pancreatic DNase I treatment

Taking into account recently obtained data indicating the participation of circulating extracellular DNA (exDNA) in tumorigenesis, enzymes with deoxyribonucleic activity have again been considered as potential antitumour and antimetastatic drugs. Previously, using murine Lewis lung carcinoma and hepatocellular carcinoma A1 tumour models, we have shown the antimetastatic activity of bovine DNase I, which correlates with an increase of DNase activity and a decrease of exDNA concentration in the blood serum of tumour-bearing mice. In this work, using next-generation sequencing on the ABS SOLiD^™ 5.500 platform, we performed a search for molecular targets of DNase I by comparing the exDNA profiles of healthy animals, untreated animals with Lewis lung carcinoma (LLC) and those with LLC treated with DNase I. We found that upon DNase I treatment of LLC-bearing mice, together with inhibition of metastasis, a number of strong alterations in the patterns of exDNA were observed. The major differences in exDNA profiles between groups were: i) the level of GC-poor sequences increased during tumour development was reduced to that of healthy mice; ii) levels of sequences corresponding to tumour-associated genes Hmga2, Myc and Jun were reduced in the DNase I-treated group in comparison with non-treated mice; iii) 224 types of tandem repeat over-presented in untreated LLC-bearing mice were significantly reduced after DNase I treatment. The most important result obtained in the work is that DNase I decreased the level of B-subfamily repeats having homology to human ALU repeats, known as markers of carcinogenesis, to the level of healthy animals. Thus, the obtained data lead us to suppose that circulating exDNA plays a role in tumour dissemination, and alteration of multiple molecular targets in the bloodstream by DNase I reduces the invasive potential of tumours.

Heavy Metal Exposure Influences Double Strand Break DNA Repair Outcomes

Heavy metals such as cadmium, arsenic and nickel are classified as carcinogens. Although the precise mechanism of carcinogenesis is undefined, heavy metal exposure can contribute to genetic damage by inducing double strand breaks (DSBs) as well as inhibiting critical proteins from different DNA repair pathways. Here we take advantage of two previously published culture assay systems developed to address mechanistic aspects of DNA repair to evaluate the effects of heavy metal exposures on competing DNA repair outcomes. Our results demonstrate that exposure to heavy metals significantly alters how cells repair double strand breaks. The effects observed are both specific to the particular metal and dose dependent. Low doses of NiCl2 favored resolution of DSBs through homologous recombination (HR) and single strand annealing (SSA), which were inhibited by higher NiCl2 doses. In contrast, cells exposed to arsenic trioxide preferentially repaired using the "error prone" non-homologous end joining (alt-NHEJ) while inhibiting repair by HR. In addition, we determined that low doses of nickel and cadmium contributed to an increase in mutagenic recombination-mediated by Alu elements, the most numerous family of repetitive elements in humans. Sequence verification confirmed that the majority of the genetic deletions were the result of Alu-mediated non-allelic recombination events that predominantly arose from repair by SSA. All heavy metals showed a shift in the outcomes of alt-NHEJ repair with a significant increase of non-templated sequence insertions at the DSB repair site. Our data suggest that exposure to heavy metals will alter the choice of DNA repair pathway changing the genetic outcome of DSBs repair.

Heavy Metal Exposure Influences Double Strand Break DNA Repair Outcomes

Heavy metals such as cadmium, arsenic and nickel are classified as carcinogens. Although the precise mechanism of carcinogenesis is undefined, heavy metal exposure can contribute to genetic damage by inducing double strand breaks (DSBs) as well as inhibiting critical proteins from different DNA repair pathways. Here we take advantage of two previously published culture assay systems developed to address mechanistic aspects of DNA repair to evaluate the effects of heavy metal exposures on competing DNA repair outcomes. Our results demonstrate that exposure to heavy metals significantly alters how cells repair double strand breaks. The effects observed are both specific to the particular metal and dose dependent. Low doses of NiCl₂ favored resolution of DSBs through homologous recombination (HR) and single strand annealing (SSA), which were inhibited by higher NiCl₂ doses. In contrast, cells exposed to arsenic trioxide preferentially repaired using the “error prone” non-homologous end joining (alt-NHEJ) while inhibiting repair by HR. In addition, we determined that low doses of nickel and cadmium contributed to an increase in mutagenic recombination-mediated by Alu elements, the most numerous family of repetitive elements in humans. Sequence verification confirmed that the majority of the genetic deletions were the result of Alu-mediated non-allelic recombination events that predominantly arose from repair by SSA. All heavy metals showed a shift in the outcomes of alt-NHEJ repair with a significant increase of non-templated sequence insertions at the DSB repair site. Our data suggest that exposure to heavy metals will alter the choice of DNA repair pathway changing the genetic outcome of DSBs repair.

Lack of major genome instability in tumors of p53 null rats

Tumorigenesis is often associated with loss of tumor suppressor genes (such as TP53), genomic instability and telomere lengthening. Previously, we generated and characterized a rat p53 knockout model in which the homozygous rats predominantly develop hemangiosarcomas whereas the heterozygous rats mainly develop osteosarcomas. Using genome-wide analyses, we find that the tumors that arise in the heterozygous and homozygous Tp53^C273X mutant animals are also different in their genomic instability profiles. While p53 was fully inactivated in both heterozygous and homozygous knockout rats, tumors from homozygous animals show very limited aneuploidy and low degrees of somatic copy number variation as compared to the tumors from heterozygous animals. In addition, complex structural rearrangements such as chromothripsis and breakage-fusion-bridge cycles were never found in tumors from homozygous animals, while these were readily detectable in tumors from heterozygous animals. Finally, we measured telomere length and telomere lengthening pathway activity and found that tumors of homozygous animals have longer telomeres but do not show clear telomerase or alternative lengthening of telomeres (ALT) activity differences as compared to the tumors from heterozygous animals. Taken together, our results demonstrate that host p53 status in this rat p53 knockout model has a large effect on both tumor type and genomic instability characteristics, where full loss of functional p53 is not the main driver of large-scale structural variations. Our results also suggest that chromothripsis primarily occurs under p53 heterozygous rather than p53 null conditions.

A comparison of 100 human genes using an alu element-based instability model

The human retrotransposon with the highest copy number is the Alu element. The human genome contains over one million Alu elements that collectively account for over ten percent of our DNA. Full-length Alu elements are randomly distributed throughout the genome in both forward and reverse orientations. However, full-length widely spaced Alu pairs having two Alus in the same (direct) orientation are statistically more prevalent than Alu pairs having two Alus in the opposite (inverted) orientation. The cause of this phenomenon is unknown. It has been hypothesized that this imbalance is the consequence of anomalous inverted Alu pair interactions. One proposed mechanism suggests that inverted Alu pairs can ectopically interact, exposing both ends of each Alu element making up the pair to a potential double-strand break, or "hit". This hypothesized "two-hit" (two double-strand breaks) potential per Alu element was used to develop a model for comparing the relative instabilities of human genes. The model incorporates both 1) the two-hit double-strand break potential of Alu elements and 2) the probability of exon-damaging deletions extending from these double-strand breaks. This model was used to compare the relative instabilities of 50 deletion-prone cancer genes and 50 randomly selected genes from the human genome. The output of the Alu element-based genomic instability model developed here is shown to coincide with the observed instability of deletion-prone cancer genes. The 50 cancer genes are collectively estimated to be 58% more unstable than the randomly chosen genes using this model. Seven of the deletion-prone cancer genes, ATM, BRCA1, FANCA, FANCD2, MSH2, NCOR1 and PBRM1, were among the most unstable 10% of the 100 genes analyzed. This algorithm may lay the foundation for comparing genetic risks posed by structural variations that are unique to specific individuals, families and people groups.