Umeå University's logo

The genomes of immortalized cell lines (and cancer cells) are characterized by multiple types of aberrations, ranging from single nucleotide polymorphisms (SNPs) to structural rearrangements that have accumulated over time. Consequently, it is difficult to estimate the relative impact of different aberrations, the order of events, and which gene functions were under selective pressure at the early stage towards cellular immortalization. Here, we have established novel cell cultures derived from Drosophila melanogaster embryos that were sampled at multiple time points over a one-year period. Using short-read DNA sequencing, we show that copy-number gain in preferentially stress-related genes were acquired in a dominant fraction of cells in 300-days old cultures. Furthermore, transposable elements were active in cells of all cultures. Only a few (<1%) SNPs could be followed over time, and these showed no trend to increase or decrease. We conclude that the early cellular responses of a novel culture comprise sequence duplication and transposable element activity. During immortalization, positive selection first occurs on genes that are related to stress response before shifting to genes that are related to growth.

Cells are the vehicles that allows genetic code to proliferate in the world, taking on various forms – as illustrated by the tree of life. The cell features are determined by the manufacturing of proteins, a process that has its blueprints encoded as genes in the genome. It is crucial for all cells to have the right amount of protein, regardless of context (part of a multicellular organism or self-sustained). The protein landscape (amount and type) vary depending on the environment. Cells of the multicellular organism should maintain the protein balance to provide its’ intended function in the organism tissue. The cells of multicellular organisms are faced with an imbalance due to sex-related chromosomal imbalances and other genome effects that change the number of gene copies. Restoration from the imbalance is done by dosage compensation systems. Cells that are isolated from the organism and grown inside the lab are common in research, known as cell lines. Cancer cells are similar to cell lines and have lost their original function in the organism in favor of a self-sustained lifestyle. The new environment (context) for these isolated cells impose a challenge; the cells must reorganize their genomes (holding the blueprints for proteins) to obtain autonomy.

In this thesis, the genome evolution of isolated cells, cell lines, has been studied using Drosophila melanogaster (the fruit fly). Compared to normal cells of the host organism, cell line genomes are highly mutated and rearranged. With the emergence of novel sequencing technologies that can read long fragments of the genome, this complexity of cell line genomes can be captured. On the topic of novel sequencing technologies, new software implementations are presented and the future of software for long reads and complex genomes is discussed. The main focus of this thesis is to describe how an established and commonly used cell line has reorganized its’ genome to sustain a culture environment. Important information about the genome structure is provided to the research community. The thesis also describes the genome reorganization in new cell lines, covering the early adaptations to cell autonomy. Together, these investigations are of high relevance to human cancer research. This thesis has also studied the fundamentals for regulation of protein balance in organismal cells. Specifically, a recognition sequence to the X chromosome of the protein Painting of Fourth. This protein is related to dosage compensation and primarily enhance transcription from the 4 thchromosome in Drosophila melanogaster, but has been observed tooccasionally bind to the X chromosome.

Background: Immortalized cell lines are widely used model systems whose genomes are often highly rearranged and polyploid. However, their genome structure is seldom deciphered and is thus not accounted for during analyses. We therefore used linked short- and long-read sequencing to perform haplotype-level reconstruction of the genome of a Drosophila melanogaster cell line (S2-DRSC) with a complex genome structure.

Results: Using a custom implementation (that is designed to use ultra-long reads in complex genomes with nested rearrangements) to call structural variants (SVs), we found that the most common SV was repetitive sequence insertion or deletion (> 80% of SVs), with Gypsy retrotransposon insertions dominating. The second most common SV was local sequence duplication. SNPs and other SVs were rarer, but several large chromosomal translocations and mitochondrial genome insertions were observed. Haplotypes were highly similar at the nucleotide level but structurally very different. Insertion SVs existed at various haplotype frequencies and were unlinked on chromosomes, demonstrating that haplotypes have different structures and suggesting the existence of a mechanism that allows SVs to propagate across haplotypes. Finally, using public short-read data, we found that transposable element insertions and local duplications are common in other D. melanogaster cell lines.

Conclusions: The S2-DRSC cell line evolved through retrotransposon activity and vast local sequence duplications, that we hypothesize were the products of DNA re-replication events. Additionally, mutations can propagate across haplotypes (possibly explained by mitotic recombination), which enables fine-tuning of mutational impact and prevents accumulation of deleterious events, an inherent problem of clonal reproduction. We conclude that traditional linear homozygous genome representation conceals the complexity when dealing with rearranged and heterozygous clonal cells.

Chromosome-specific regulatory mechanisms provide a model to understand the coordinated regulation of genes on entire chromosomes or on larger genomic regions. In fruit flies, two chromosome-wide systems have been characterized: The male-specific lethal (MSL) complex, which mediates dosage compensation and primarily acts on the male X-chromosome, and Painting of fourth (POF), which governs chromosome-specific regulation of genes located on the 4th chromosome. How targeting of one specific chromosome evolves is still not understood; but repeated sequences, in forms of satellites and transposable elements, are thought to facilitate the evolution of chromosome-specific targeting. The highly repetitive 1.688 satellite has been functionally connected to both these systems. Considering the rapid evolution and the necessarily constant adaptation of regulatory mechanisms, such as dosage compensation, we hypothesised that POF and/or 1.688 may still show traces of dosage-compensation functions. Here, we test this hypothesis by transcriptome analysis. We show that loss of Pof decreases not only chromosome 4 expression but also reduces the X-chromosome expression in males. The 1.688 repeat deletion, Zhr¹(Zygotic hybrid rescue), does not affect male dosage compensation detectably; however, Zhr¹ in females causes a stimulatory effect on X-linked genes with a strong binding affinity to the MSL complex (genes close to high-affinity sites). Lack of pericentromeric 1.688 also affected 1.688 expression in trans and was linked to the differential expression of genes involved in eggshell formation. We discuss our results with reference to the connections between POF, the 1.688 satellite and dosage compensation, and the role of the 1.688 satellite in hybrid lethality.

In specific cases, chromatin clearly forms long-range loops that place distant regulatory elements in close proximity to transcription start sites, but we have limited understanding of many loops identified by Chromosome Conformation Capture (such as Hi-C) analyses. In efforts to elucidate their characteristics and functions, we have identified highly interacting regions (HIRs) using intra-chromosomal Hi-C datasets with a new computational method based on looking at the eigenvector that corresponds to the smallest eigenvalue (here unity). Analysis of these regions using ENCODE data shows that they are in general enriched in bound factors involved in DNA damage repair and have actively transcribed genes. However, both highly transcribed regions as well as transcriptionally inactive regions can form HIRs. The results also indicate that enhancers and super-enhancers in particular form long-range interactions within the same chromosome. The accumulation of DNA repair factors in most identified HIRs suggests that protection from DNA damage in these regions is essential for avoidance of detrimental rearrangements.

Repetitive DNA, represented by transposons and satellite DNA, constitutes a large portion of eukaryotic genomes, being the major component of constitutive heterochromatin. There is a growing body of evidence that it regulates several nuclear functions including chromatin state and the proper functioning of centromeres and telomeres. The 1.688 satellite is one of the most abundant repetitive sequences in Drosophila melanogaster, with the longest array being located in the pericentromeric region of the X-chromosome. Short arrays of 1.688 repeats are widespread within the euchromatic part of the X-chromosome, and these arrays were recently suggested to assist in recognition of the X-chromosome by the dosage compensation male-specific lethal complex. We discovered that a short array of 1.688 satellite repeats is essential for recruitment of the protein POF to a previously described site on the X-chromosome (PoX2) and to various transgenic constructs. On an isolated target, i.e., an autosomic transgene consisting of a gene upstream of 1.688 satellite repeats, POF is recruited to the transgene in both males and females. The sequence of the satellite, as well as its length and position within the recruitment element, are the major determinants of targeting. Moreover, the 1.688 array promotes POF targeting to the roX1-proximal PoX1 site in trans Finally, binding of POF to the 1.688-related satellite-enriched sequences is conserved in evolution. We hypothesize that the 1.688 satellite functioned in an ancient dosage compensation system involving POF targeting to the X-chromosome.