Introduction to molecular combing: genomics, DNA replication, and cancer

Sequencing and Optical Genome Mapping for the Adventurous Chemist

This review provides a comprehensive overview of the chemistries and workflows of the sequencing methods that have been or are currently commercially available, providing a very brief historical introduction to each method. The main optical genome mapping approaches are introduced in the same manner, although only a subset of these are or have ever been commercially available. The review comes with a deck of slides containing all of the figures for ease of access and consultation.

Molecular combing and its application in clinical settings

Molecular combing technology (MCT) is an effective means for stretching DNA molecules and making them thus accessible for in situ studies. MCT uses the force exerted in the process of liquid flow via surface tension to stretch DNA molecules and spread them on solid surfaces, i.e. glass cover slips. Many DNA molecules can be stretched at the same time in parallel and neatly arranged side-by-side, making the approach convenient for statistical analysis. Accordingly, DNA replication and transcription can be studied at the single molecule level. In this paper, the principle, experimental methods, important applications, advantages and shortcuts of MCT in medical field are presented and discussed.

Analysis of Omics Data Reveals Nucleotide Excision Repair-Related Genes Signature in Highly-Grade Serous Ovarian Cancer to Predict Prognosis

Most of the high-grade serous ovarian cancers (HGSOC) are accompanied by P53 mutations, which are related to the nucleotide excision repair (NER) pathway. This study aims to construct a risk signature based on NER-related genes that could effectively predict the prognosis for advanced patients with HGSOC. In our study, we found that two clusters of HGSOC with significantly different overall survival (OS) were identified by consensus clustering and principal component analysis (PCA). Then, a 7-gene risk signature (DDB2, POLR2D, CCNH, XPC, ERCC2, ERCC4, and RPA2) for OS prediction was developed subsequently based on TCGA cohort, and the risk score-based signature was identified as an independent prognostic indicator for HGSOC. According to the risk score, HGSOC patients were divided into high-risk group and low-risk group, in which the distinct OS and the predictive power were also successfully verified in the GEO validation sets. Then we constructed a nomogram, including the risk signature and clinical-related risk factors (age and treatment response) that predicted an individual’s risk of OS, which can be validated by assessing calibration curves. Furthermore, GSEA showed that the genes in the high-risk group were significantly enriched in cancer-related pathways, such as “MAPK signaling pathway”, “mTOR signaling pathway”, “VEGF signaling pathway” and so on. In conclusion, our study has developed a robust NER-related genes-based molecular signature for prognosis prediction, and the nomogram could be used as a convenient tool for OS evaluation and guidance of therapeutic strategies in advanced patients with HGSOC.

CtIP is essential for telomere replication

The maintenance of telomere length is critical to longevity and survival. Specifically, the failure to properly replicate, resect, and/or form appropriate telomeric structures drives telomere shortening and, in turn, genomic instability. The endonuclease CtIP is a DNA repair protein that is well-known to promote genome stability through the resection of endogenous DNA double-stranded breaks. Here, we describe a novel role for CtIP. We show that in the absence of CtIP, human telomeres shorten rapidly to non-viable lengths. This telomere dysfunction results in an accumulation of fusions, breaks, and frank telomere loss. Additionally, CtIP suppresses the generation of circular, extrachromosomal telomeric DNA. These latter structures appear to arise from arrested DNA replication forks that accumulate in the absence of CtIP. Hence, CtIP is required for faithful replication through telomeres via its roles at stalled replication tracts. Our findings demonstrate a new role for CtIP as a protector of human telomere integrity.

Mild replication stress causes aneuploidy by deregulating microtubule dynamics in mitosis

Chromosomal instability (CIN) causes structural and numerical chromosome aberrations and represents a hallmark of cancer. Replication stress (RS) has emerged as a driver for structural chromosome aberrations while mitotic defects can cause whole chromosome missegregation and aneuploidy. Recently, first evidence indicated that RS can also influence chromosome segregation in cancer cells exhibiting CIN, but the underlying mechanisms remain unknown. Here, we show that chromosomally unstable cancer cells suffer from very mild RS, which allows efficient proliferation and which can be mimicked by treatment with very low concentrations of aphidicolin. Both, endogenous RS and aphidicolin-induced very mild RS cause chromosome missegregation during mitosis leading to the induction of aneuploidy. Moreover, RS triggers an increase in microtubule plus end growth rates in mitosis, an abnormality previously identified to cause chromosome missegregation in cancer cells. In fact, RS-induced chromosome missegregation is mediated by increased mitotic microtubule growth rates and is suppressed after restoration of proper microtubule growth rates and upon rescue of replication stress. Hence, very mild and cancer-relevant RS triggers aneuploidy by deregulating microtubule dynamics in mitosis.

A simple and novel DNA combing methodology for Fiber-FISH and optical mapping

Single molecule analysis can help us study genomics efficiently. It involves studying single DNA molecules for genomic studies. DNA combing is one of such techniques which allowed us to study single DNA molecules for multiple uses. DNA combing technology can be used to perform Fiber-FISH and optical mapping. Physical mapping of genomes can be studied by restriction digestion of combed DNA on glass slides. Restriction fragments can be arranged into optical maps by gathering fluorescent intensity data by CCD camera and image analysis by softwares. Physical mapping and DNA segment rearrangements can be studied by Fiber-FISH which involves application of probes on genomic DNA combed over glass slides. We developed a novel methodology involving combing solution optimization, denatured combed DNA and performed restriction digestion of combed DNA. Thus we provided an efficient and robust combing platform for its application in Fiber-FISH and optical mapping.

Influence of concentration on distribution properties of stretched-DNA in the MEC studied with fluorescence imaging and drop shape analyzing

Stretching and manipulating DNA efficiently is significant for exploring the properties and applications of single DNA molecules. Here, the influence of concentrations of buffer and DNA on properties of stretched DNA molecules in the molecular evaporation combing (MEC) is investigated systematically with the single molecule fluorescence imaging microscopy and the high-precision drop shape analyzing technology. The stretched degree and uniformity of combed DNA molecules decrease as the buffer concentration are increased from 7 to 20mM. When the buffer concentration changes from 12 to 15mM, the stretched DNA molecules are apt to form a ringlike pattern. During the MEC process, there exist two kinds of evaporation modes, i.e., the constant contact angle mode and the constant contact radius mode. The former only takes effect in the lower concentration of buffer and DNA, enabling the uniform stretching. While the latter plays the leading role in the higher concentration, promoting the formation of the ringlike pattern of DNA molecules.

Molecular Combing of λ-DNA using Self-Propelled Water Droplets on Wettability Gradient Surfaces

Surface wettability gradients were used to elongate and align double stranded λ-DNA. Gradients were prepared by vapor phase deposition of octyltrichlorosilane (C8-silane) and fluorinated octyltrichlorosilane (F-silane) precursors. Gradient formation was confirmed by water contact angle and ellipsometric film thickness measurements. Placement of a droplet of aqueous DNA solution on the hydrophobic end of each gradient led to spontaneous motion of the droplet toward the hydrophilic end and deposition of the DNA. Fluorescence imaging of surface-adsorbed YOYO-1 labeled DNA molecules revealed that they are elongated and aligned perpendicular to the droplet-surface contact line at all positions along the gradient, consistent with a dominant role played by surface tension forces in elongating the DNA. The density of adsorbed DNA was found to be greatest on the C8-silane gradient at its hydrophobic end. DNA density decreased toward the hydrophilic end, while the length of the elongated DNA was less dependent on position. The elongation of DNA molecules by spontaneous droplet motion on chemical gradient surfaces has possible applications in DNA barcoding and studies of DNA-protein interactions.

Folate levels modulate oncogene-induced replication stress and tumorigenicity

Chromosomal instability in early cancer stages is caused by replication stress. One mechanism by which oncogene expression induces replication stress is to drive cell proliferation with insufficient nucleotide levels. Cancer development is driven by alterations in both genetic and environmental factors. Here, we investigated whether replication stress can be modulated by both genetic and non-genetic factors and whether the extent of replication stress affects the probability of neoplastic transformation. To do so, we studied the effect of folate, a micronutrient that is essential for nucleotide biosynthesis, on oncogene-induced tumorigenicity. We show that folate deficiency by itself leads to replication stress in a concentration-dependent manner. Folate deficiency significantly enhances oncogene-induced replication stress, leading to increased DNA damage and tumorigenicity in vitro. Importantly, oncogene-expressing cells, when grown under folate deficiency, exhibit a significantly increased frequency of tumor development in mice. These findings suggest that replication stress is a quantitative trait affected by both genetic and non-genetic factors and that the extent of replication stress plays an important role in cancer development.

The hunt for origins of DNA replication in multicellular eukaryotes

Origins of DNA replication (ORIs) occur at defined regions in the genome. Although DNA sequence defines the position of ORIs in budding yeast, the factors for ORI specification remain elusive in metazoa. Several methods have been used recently to map ORIs in metazoan genomes with the hope that features for ORI specification might emerge. These methods are reviewed here with analysis of their advantages and shortcomings. The various factors that may influence ORI selection for initiation of DNA replication are discussed.

Cellular response to the genotoxic insult: the question of threshold for genotoxic carcinogens

Maintenance of cellular integrity is crucial for its physiological function, which is constantly threatened by DNA damage arising from numerous intrinsic and environmental sources.

Best practices for mapping replication origins in eukaryotic chromosomes

Understanding the regulatory principles ensuring complete DNA replication in each cell division is critical for deciphering the mechanisms that maintain genomic stability. Recent advances in genome sequencing technology facilitated complete mapping of DNA replication sites and helped move the field from observing replication patterns at a handful of single loci to analyzing replication patterns genome-wide. These advances address issues, such as the relationship between replication initiation events, transcription, and chromatin modifications, and identify potential replication origin consensus sequences. This unit summarizes the technological and fundamental aspects of replication profiling and briefly discusses novel insights emerging from mining large datasets, published in the last 3 years, and also describes DNA replication dynamics on a whole-genome scale.

Transcriptional dynamics in colorectal carcinogenesis: new insights into the role of c-Myc and miR17 in benign to cancer transformation

Colorectal cancer develops in a sequential, evolutionary process, leading to a heterogenic tumor. Comprehensive molecular studies of colorectal cancer have been previously performed; still, the process of carcinogenesis is not fully understood. We utilized gene expression patterns from 94 samples including normal, adenoma, and adenocarcinoma colon biopsies and performed a coexpression network analysis to determine gene expression trajectories of 8,000 genes across carcinogenesis. We found that the majority of gene expression changes occur in the transition from normal tissue to adenoma. The upregulated genes, known to be involved in cellular proliferation, included c-Myc along with its targets. In a cellular model system, we show that physiologic upregulation of c-Myc can lead to cellular proliferation without DNA replication stress. Our analysis also found that carcinogenesis involves a progressive downregulation of genes that are markers of colonic tissue and propose that this reflects a perturbed differentiation of colon cells during carcinogenesis. The analysis of miRNAs targets pointed toward the involvement of miR17 in the regulation of colon cell differentiation. Finally, we found that copy-number variations (CNV) enriched in colon adenocarcinoma tend to occur in genes whose expression changes already in adenoma, with deletions occurring in genes downregulated and duplications in genes upregulated in adenomas. We suggest that the CNVs are selected to reinforce changes in gene expression, rather than initiate them. Together, these findings shed new light into the molecular processes that underlie the transformation of colon tissue from normal to cancer and add a temporal context that has been hitherto lacking.

Toward single-molecule optical mapping of the epigenome

The past decade has seen an explosive growth in the utilization of single-molecule techniques for the study of complex systems. The ability to resolve phenomena otherwise masked by ensemble averaging has made these approaches especially attractive for the study of biological systems, where stochastic events lead to inherent inhomogeneity at the population level. The complex composition of the genome has made it an ideal system to study at the single-molecule level, and methods aimed at resolving genetic information from long, individual, genomic DNA molecules have been in use for the last 30 years. These methods, and particularly optical-based mapping of DNA, have been instrumental in highlighting genomic variation and contributed significantly to the assembly of many genomes including the human genome. Nanotechnology and nanoscopy have been a strong driving force for advancing genomic mapping approaches, allowing both better manipulation of DNA on the nanoscale and enhanced optical resolving power for analysis of genomic information. During the past few years, these developments have been adopted also for epigenetic studies. The common principle for these studies is the use of advanced optical microscopy for the detection of fluorescently labeled epigenetic marks on long, extended DNA molecules. Here we will discuss recent single-molecule studies for the mapping of chromatin composition and epigenetic DNA modifications, such as DNA methylation.

Methyl-combing: single-molecule analysis of DNA methylation on stretched DNA fibers

The methyl-combing technique combines the dynamic molecular combing method with the detection of DNA modifications. The assay allows the single-molecule analysis of epigenetic marks on regularly stretched DNA fibers, at the megabase scale with kilobase resolution. The protocol presented in this chapter is based on proof-of-principle experiments where the single-molecule detection of DNA methylation has been performed on unmethylated and in vitro methylated λ phage DNA.

Capillary electrophoretic separation-based approach to determine the labeling kinetics of oligodeoxynucleotides

With the recent advances in electron microscopy (EM), computation, and nanofabrication, the original idea of reading DNA sequence directly from an image can now be tested. One approach is to develop heavy atom labels that can provide the contrast required for EM imaging. While evaluating tentative labels for the respective nucleobases in synthetic oligodeoxynucleotides (oligos), we developed a streamlined CE protocol to assess the label stability, reactivity, and selectivity. We report our protocol using osmium tetroxide 2,2'-bipyridine (Osbipy) as a thymidine (T) specific label. The observed rates show that the labeling process is kinetically independent of both the oligo length, and the base composition. The conditions, i.e. temperature, optimal Osbipy concentration, and molar ratio of reagents, to promote 100% conversion of the starting oligo to labeled product were established. Hence, the optimized conditions developed with the oligos could be leveraged to allow osmylation of effectively all Ts in ssDNA, while achieving minimal mislabeling. In addition, the approach and methods employed here may be adapted to the evaluation of other prospective contrasting agents/labels to facilitate next-generation DNA sequencing by EM.

A general approach to break the concentration barrier in single-molecule imaging

Single-molecule fluorescence imaging is often incompatible with physiological protein concentrations, as fluorescence background overwhelms an individual molecule's signal. We solve this problem with a new imaging approach called PhADE (PhotoActivation, Diffusion and Excitation). A protein of interest is fused to a photoactivatable protein (mKikGR) and introduced to its surface-immobilized substrate. After photoactivation of mKikGR near the surface, rapid diffusion of the unbound mKikGR fusion out of the detection volume eliminates background fluorescence, whereupon the bound molecules are imaged. We labeled the eukaryotic DNA replication protein flap endonuclease 1 with mKikGR and added it to replication-competent Xenopus laevis egg extracts. PhADE imaging of high concentrations of the fusion construct revealed its dynamics and micrometer-scale movements on individual, replicating DNA molecules. Because PhADE imaging is in principle compatible with any photoactivatable fluorophore, it should have broad applicability in revealing single-molecule dynamics and stoichiometry of macromolecular protein complexes at previously inaccessible fluorophore concentrations.

A diagnostic genetic test for the physical mapping of germline rearrangements in the susceptibility breast cancer genes BRCA1 and BRCA2

The BRCA1 and BRCA2 genes are involved in breast and ovarian cancer susceptibility. About 2 to 4% of breast cancer patients with positive family history, negative for point mutations, can be expected to carry large rearrangements in one of these two genes. We developed a novel diagnostic genetic test for the physical mapping of large rearrangements, based on molecular combing (MC), a FISH-based technique for direct visualization of single DNA molecules at high resolution. We designed specific Genomic Morse Codes (GMCs), covering the exons, the noncoding regions, and large genomic portions flanking both genes. We validated our approach by testing 10 index cases with positive family history of breast cancer and 50 negative controls. Large rearrangements, corresponding to deletions and duplications with sizes ranging from 3 to 40 kb, were detected and characterized on both genes, including four novel mutations. The nature of all the identified mutations was confirmed by high-resolution array comparative genomic hybridization (aCGH) and breakpoints characterized by sequencing. The developed GMCs allowed to localize several tandem repeat duplications on both genes. We propose the developed genetic test as a valuable tool to screen large rearrangements in BRCA1 and BRCA2 to be combined in clinical settings with an assay capable of detecting small mutations.

The MMS22L-TONSL complex mediates recovery from replication stress and homologous recombination

Genome integrity is jeopardized each time DNA replication forks stall or collapse. Here we report the identification of a complex composed of MMS22L (C6ORF167) and TONSL (NFKBIL2) that participates in the recovery from replication stress. MMS22L and TONSL are homologous to yeast Mms22 and plant Tonsoku/Brushy1, respectively. MMS22L-TONSL accumulates at regions of ssDNA associated with distressed replication forks or at processed DNA breaks, and its depletion results in high levels of endogenous DNA double-strand breaks caused by an inability to complete DNA synthesis after replication fork collapse. Moreover, cells depleted of MMS22L are highly sensitive to camptothecin, a topoisomerase I poison that impairs DNA replication progression. Finally, MMS22L and TONSL are necessary for the efficient formation of RAD51 foci after DNA damage, and their depletion impairs homologous recombination. These results indicate that MMS22L and TONSL are genome caretakers that stimulate the recombination-dependent repair of stalled or collapsed replication forks.

Replication dynamics at common fragile site FRA6E

The replication dynamics at common fragile site FRA6E has been evaluated by molecular combing and interphase fluorescent in situ hybridisation (FISH) in primary human lymphocytes cultured under normal or aphidicolin-induced stress conditions. FRA6E is one of the most frequently expressed common fragile sites of the human genome. It harbours several genes, PARK2 being regarded as the most relevant one. According to the results obtained from interphase FISH analysis, FRA6E can be considered a mid-late-replicating sequence characterised by heterogeneous replication timing. Molecular combing did not reveal specific replication parameters at the fragile site: fork rates were highly comparable to those detected at an early replicating locus (LMNB2) used as control and in very good agreement with the whole-genome data obtained in parallel. The same indication applied to the density of initiation zones, the inter-origin distances from adjacent ongoing forks, the frequencies of unidirectional forks, fork arrest events and asynchronous forks. Interestingly, PARK2 appeared embedded in an early/late replication transition zone, corresponding to intron 8 (162 kb) and to the fragility core of FRA6E. In cells exposed to aphidicolin, few forks progressing at a rather slow rate were observed, the majority of them being unidirectional, but again a specific response of the fragile site was not observed. In summary, at FRA6E the replication process is not impaired per se, but chromosome breakages occur preferentially at an early/late replication transition zone. Aphidicolin might increase the occurrence of breakage events at FRA6E by prolonging the time interval separating the replication of early and late replication domains. These results may be of general significance to address the problem of fragile site instability.

A winding road to origin discovery

Studies in our laboratory over the last three decades have shown that the Chinese hamster dihydrofolate reductase (DHFR) origin of replication corresponds to a broad zone of inefficient initiation sites distributed throughout the spacer between the convergently transcribed DHFR and 2BE2121 genes. It is clear from mutational analysis that none of these sites is genetically required for controlling origin activity. However, the integrity of the promoter of the DHFR gene is needed to activate the downstream origin, while the 3' processing signals prevent invasion and inactivation of the downstream origin by transcription forks. Several other origins in metazoans have been shown to correspond to zones of inefficient sites, while a different subset appears to be similar to the fixed replicators that characterize origins in S. cerevisiae and lower organisms. These observations have led us to suggest a model in which the mammalian genome is dotted with a hierarchy of degenerate, redundant, and inefficient replicators at intervals of a kilobase or less, some of which may have evolved to be highly circumscribed and efficient. The activities of initiation sites are proposed to be largely regulated by local transcription and chromatin architecture. Recently, we and others have devised strategies for identifying active origins on a genome-wide scale in order to define their distributions between fixed and dispersive origin types and to detect relationships among origins, genes, and epigenetic markers. The global pictures emerging are suggestive but far from complete and appear to be plagued by some of the same uncertainties that have led to conflicting views of individual origins in the past (particularly DHFR). In this paper, we will trace the history of origin discovery in mammalian genomes, primarily using the well-studied DHFR origin as a model, because it has been analyzed by nearly every available origin mapping technique in several different laboratories, while many origins have been identified by only one. We will address the strengths and shortcomings of the various methods utilized to identify and characterize origins in complex genomes and will point out how we and others were sometimes led astray by false assumptions and biases, as well as insufficient information. The goal is to help guide future experiments that will provide a truly comprehensive and accurate portrait of origins and their regulation. After all, in the words of George Santayana, "Those who do not learn from history are doomed to repeat it."