RAD51-mediated R-loop formation acts to repair transcription-associated DNA breaks driving antigenic variation in Trypanosoma brucei.

RNA-DNA hybrids are epigenetic features of all genomes that intersect with many processes, including transcription, telomere homeostasis, and centromere function. Increasing evidence suggests that RNA-DNA hybrids can provide two conflicting roles in the maintenance and transmission of genomes: They can be the triggers of DNA damage, leading to genome change, or can aid the DNA repair processes needed to respond to DNA lesions. Evasion of host immunity by African trypanosomes, such as Trypanosoma brucei, relies on targeted recombination of silent Variant Surface Glycoprotein (VSG) genes into a specialized telomeric locus that directs transcription of just one VSG from thousands. How such VSG recombination is targeted and initiated is unclear. Here, we show that a key enzyme of T. brucei homologous recombination, RAD51, interacts with RNA-DNA hybrids. In addition, we show that RNA-DNA hybrids display a genome-wide colocalization with DNA breaks and that this relationship is impaired by mutation of RAD51. Finally, we show that RAD51 acts to repair highly abundant, localised DNA breaks at the single transcribed VSG and that mutation of RAD51 alters RNA-DNA hybrid abundance at 70 bp repeats both around the transcribed VSG and across the silent VSG archive. This work reveals a widespread, generalised role for RNA-DNA hybrids in directing RAD51 activity during recombination and uncovers a specialised application of this interplay during targeted DNA break repair needed for the critical T. brucei immune evasion reaction of antigenic variation.

Immunoprecipitation of RNA-DNA hybrid interacting proteins in Trypanosoma brucei reveals conserved and novel activities, including in the control of surface antigen expression needed for immune evasion by antigenic variation.

RNA-DNA hybrids are epigenetic features of genomes that provide a diverse and growing range of activities. Understanding of these functions has been informed by characterising the proteins that interact with the hybrids, but all such analyses have so far focused on mammals, meaning it is unclear if a similar spectrum of RNA-DNA hybrid interactors is found in other eukaryotes. The African trypanosome is a single-cell eukaryotic parasite of the Discoba grouping and displays substantial divergence in several aspects of core biology from its mammalian host. Here, we show that DNA-RNA hybrid immunoprecipitation coupled with mass spectrometry recovers 602 putative interactors in T. brucei mammal- and insect-infective cells, some providing activities also found in mammals and some lineage-specific. We demonstrate that loss of three factors, two putative helicases and a RAD51 paralogue, alters T. brucei nuclear RNA-DNA hybrid and DNA damage levels. Moreover, loss of each factor affects the operation of the parasite immune survival mechanism of antigenic variation. Thus, our work reveals the broad range of activities contributed by RNA-DNA hybrids to T. brucei biology, including new functions in host immune evasion as well as activities likely fundamental to eukaryotic genome function.

Conditional knockout of RAD51-related genes in Leishmania major reveals a critical role for homologous recombination during genome replication.

Homologous recombination (HR) has an intimate relationship with genome replication, both during repair of DNA lesions that might prevent DNA synthesis and in tackling stalls to the replication fork. Recent studies led us to ask if HR might have a more central role in replicating the genome of Leishmania, a eukaryotic parasite. Conflicting evidence has emerged regarding whether or not HR genes are essential, and genome-wide mapping has provided evidence for an unorthodox organisation of DNA replication initiation sites, termed origins. To answer this question, we have employed a combined CRISPR/Cas9 and DiCre approach to rapidly generate and assess the effect of conditional ablation of RAD51 and three RAD51-related proteins in Leishmania major. Using this approach, we demonstrate that loss of any of these HR factors is not immediately lethal but in each case growth slows with time and leads to DNA damage and accumulation of cells with aberrant DNA content. Despite these similarities, we show that only loss of RAD51 or RAD51-3 impairs DNA synthesis and causes elevated levels of genome-wide mutation. Furthermore, we show that these two HR factors act in distinct ways, since ablation of RAD51, but not RAD51-3, has a profound effect on DNA replication, causing loss of initiation at the major origins and increased DNA synthesis at subtelomeres. Our work clarifies questions regarding the importance of HR to survival of Leishmania and reveals an unanticipated, central role for RAD51 in the programme of genome replication in a microbial eukaryote.

ASCIIGenome: a command line genome browser for console terminals.

SUMMARY: Current genome browsers are designed to work via graphical user interfaces (GUIs), which, however intuitive, are not amenable to operate within console terminals and therefore are difficult to streamline or integrate in scripts. To circumvent these limitations, ASCIIGenome runs exclusively via command line interface to display genomic data directly in a terminal window. By following the same philosophy of UNIX tools, ASCIIGenome aims to be easily integrated with the command line, including batch processing of data, and therefore enables an effective exploration of the data. IMPLEMENTATION: ASCIIGenome is written in Java. Consequently, it is a cross-platform tool and requires minimal or no installation. Some of the common genomic data types are supported and data access on remote ftp servers is possible. Speed and memory footprint are comparable to or better than those of common genome browsers. AVAILABILITY AND IMPLEMENTATION: Software and source code (MIT License) are available at https://github.com/dariober/ASCIIGenome with detailed documentation at http://asciigenome.readthedocs.io . CONTACT: Dario.beraldi@cruk.cam.ac.uk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Sequencing 5-Hydroxymethyluracil at Single-Base Resolution.

5-hydroxymethyluracil (5hmU) is formed through oxidation of thymine both enzymatically and non-enzymatically in various biological systems. Although 5hmU has been reported to affect biological processes such as protein-DNA interactions, the consequences of 5hmU formation in genomes have not been yet fully explored. Herein, we report a method to sequence 5hmU at single-base resolution. We employ chemical oxidation to transform 5hmU to 5-formyluracil (5fU), followed by the polymerase extension to induce T-to-C base changes owing to the inherent ability of 5fU to form 5fU:G base pairing. In combination with the Illumina next generation sequencing technology, we developed polymerase chain reaction (PCR) conditions to amplify the T-to-C base changes and demonstrate the method in three different synthetic oligonucleotide models as well as part of the genome of a 5hmU-rich eukaryotic pathogen. Our method has the potential capability to map 5hmU in genomic DNA and thus will contribute to promote the understanding of this modified base.

Identification and annotation of conserved promoters and macrophage-expressed genes in the pig genome.

BACKGROUND: The FANTOM5 consortium used Cap Analysis of Gene Expression (CAGE) tag sequencing to produce a comprehensive atlas of promoters and enhancers within the human and mouse genomes. We reasoned that the mapping of these regulatory elements to the pig genome could provide useful annotation and evidence to support assignment of orthology. RESULTS: For human transcription start sites (TSS) associated with annotated human-mouse orthologs, 17% mapped to the pig genome but not to the mouse, 10% mapped only to the mouse, and 55% mapped to both pig and mouse. Around 17% did not map to either species. The mapping percentages were lower where there was not clear orthology relationship, but in every case, mapping to pig was greater than to mouse, and the degree of homology was also greater. Combined mapping of mouse and human CAGE-defined promoters identified at least one putative conserved TSS for >16,000 protein-coding genes. About 54% of the predicted locations of regulatory elements in the pig genome were supported by CAGE and/or RNA-Seq analysis from pig macrophages. CONCLUSIONS: Comparative mapping of promoters and enhancers from humans and mice can provide useful preliminary annotation of other animal genomes. The data also confirm extensive gain and loss of regulatory elements between species, and the likelihood that pigs provide a better model than mice for human gene regulation and function.

Comparative analysis of monocyte subsets in the pig.

Human and mouse monocyte can be divided into two different subpopulations based on surface marker expression: CD14/16 and Ly6C/CX3CR1, respectively. Monocyte subpopulations in the pig were identified based on reciprocal expression of CD14 and the scavenger receptor CD163. The two populations, CD14(hi)-CD163(low) and CD14(low)-CD163(hi), show approximately equal abundance in the steady-state. Culture of pig PBMCs in CSF1 indicates that the two populations are a maturation series controlled by this growth factor. Gene expression in pig monocyte subpopulations was profiled using the newly developed and annotated pig whole genome snowball microarray. Previous studies have suggested a functional equivalence between human and mouse subsets, but certain genes such as CD36, CLEC4E, or TREM-1 showed human-specific expression. The same genes were expressed selectively in pig monocyte subsets. However, the profiles suggest that the pig CD14(low)-CD163(high) cells are actually equivalent to intermediate human monocytes, and there is no CD14(-) CD16(+) "nonclassical" population. The results are discussed in terms of the relevance of the pig as a model for understanding human monocyte function.

Pig bone marrow-derived macrophages resemble human macrophages in their response to bacterial lipopolysaccharide.

Mouse bone marrow-derived macrophages (BMDM) grown in M-CSF (CSF-1) have been used widely in studies of macrophage biology and the response to TLR agonists. We investigated whether similar cells could be derived from the domestic pig using human rCSF-1 and whether porcine macrophages might represent a better model of human macrophage biology. Cultivation of pig bone marrow cells for 5-7 d in presence of human rCSF-1 generated a pure population of BMDM that expressed the usual macrophage markers (CD14, CD16, and CD172a), were potent phagocytic cells, and produced TNF in response to LPS. Pig BMDM could be generated from bone marrow cells that had been stored frozen and thawed so that multiple experiments can be performed on samples from a single animal. Gene expression in pig BMDM from outbred animals responding to LPS was profiled using Affymetrix microarrays. The temporal cascade of inducible and repressible genes more closely resembled the known responses of human than mouse macrophages, sharing with humans the regulation of genes involved in tryptophan metabolism (IDO, KYN), lymphoattractant chemokines (CCL20, CXCL9, CXCL11, CXCL13), and the vitamin D3-converting enzyme, Cyp27B1. Conversely, in common with published studies of human macrophages, pig BMDM did not strongly induce genes involved in arginine metabolism, nor did they produce NO. These results establish pig BMDM as an alternative tractable model for the study of macrophage transcriptional control.

The impact of breed and tissue compartment on the response of pig macrophages to lipopolysaccharide.

BACKGROUND: The draft genome of the domestic pig (Sus scrofa) has recently been published permitting refined analysis of the transcriptome. Pig breeds have been reported to differ in their resistance to infectious disease. In this study we examine whether there are corresponding differences in gene expression in innate immune cells RESULTS: We demonstrate that macrophages can be harvested from three different compartments of the pig (lungs, blood and bone-marrow), cryopreserved and subsequently recovered and differentiated in CSF-1. We have performed surface marker analysis and gene expression profiling on macrophages from these compartments, comparing twenty-five animals from five different breeds and their response to lipopolysaccharide. The results provide a clear distinction between alveolar macrophages (AM) and monocyte-derived (MDM) and bone-marrow-derived macrophages (BMDM). In particular, the lung macrophages express the growth factor, FLT1 and its ligand, VEGFA at high levels, suggesting a distinct pathway of growth regulation. Relatively few genes showed breed-specific differential expression, notably CXCR2 and CD302 in alveolar macrophages. In contrast, there was substantial inter-individual variation between pigs within breeds, mostly affecting genes annotated as being involved in immune responses. CONCLUSIONS: Pig macrophages more closely resemble human, than mouse, in their set of macrophage-expressed and LPS-inducible genes. Future research will address whether inter-individual variation in macrophage gene expression is heritable, and might form the basis for selective breeding for disease resistance.

Structural and functional annotation of the porcine immunome.

BACKGROUND: The domestic pig is known as an excellent model for human immunology and the two species share many pathogens. Susceptibility to infectious disease is one of the major constraints on swine performance, yet the structure and function of genes comprising the pig immunome are not well-characterized. The completion of the pig genome provides the opportunity to annotate the pig immunome, and compare and contrast pig and human immune systems. RESULTS: The Immune Response Annotation Group (IRAG) used computational curation and manual annotation of the swine genome assembly 10.2 (Sscrofa10.2) to refine the currently available automated annotation of 1,369 immunity-related genes through sequence-based comparison to genes in other species. Within these genes, we annotated 3,472 transcripts. Annotation provided evidence for gene expansions in several immune response families, and identified artiodactyl-specific expansions in the cathelicidin and type 1 Interferon families. We found gene duplications for 18 genes, including 13 immune response genes and five non-immune response genes discovered in the annotation process. Manual annotation provided evidence for many new alternative splice variants and 8 gene duplications. Over 1,100 transcripts without porcine sequence evidence were detected using cross-species annotation. We used a functional approach to discover and accurately annotate porcine immune response genes. A co-expression clustering analysis of transcriptomic data from selected experimental infections or immune stimulations of blood, macrophages or lymph nodes identified a large cluster of genes that exhibited a correlated positive response upon infection across multiple pathogens or immune stimuli. Interestingly, this gene cluster (cluster 4) is enriched for known general human immune response genes, yet contains many un-annotated porcine genes. A phylogenetic analysis of the encoded proteins of cluster 4 genes showed that 15% exhibited an accelerated evolution as compared to 4.1% across the entire genome. CONCLUSIONS: This extensive annotation dramatically extends the genome-based knowledge of the molecular genetics and structure of a major portion of the porcine immunome. Our complementary functional approach using co-expression during immune response has provided new putative immune response annotation for over 500 porcine genes. Our phylogenetic analysis of this core immunome cluster confirms rapid evolutionary change in this set of genes, and that, as in other species, such genes are important components of the pig's adaptation to pathogen challenge over evolutionary time. These comprehensive and integrated analyses increase the value of the porcine genome sequence and provide important tools for global analyses and data-mining of the porcine immune response.

A gene expression atlas of the domestic pig.

BACKGROUND: This work describes the first genome-wide analysis of the transcriptional landscape of the pig. A new porcine Affymetrix expression array was designed in order to provide comprehensive coverage of the known pig transcriptome. The new array was used to generate a genome-wide expression atlas of pig tissues derived from 62 tissue/cell types. These data were subjected to network correlation analysis and clustering. RESULTS: The analysis presented here provides a detailed functional clustering of the pig transcriptome where transcripts are grouped according to their expression pattern, so one can infer the function of an uncharacterized gene from the company it keeps and the locations in which it is expressed. We describe the overall transcriptional signatures present in the tissue atlas, where possible assigning those signatures to specific cell populations or pathways. In particular, we discuss the expression signatures associated with the gastrointestinal tract, an organ that was sampled at 15 sites along its length and whose biology in the pig is similar to human. We identify sets of genes that define specialized cellular compartments and region-specific digestive functions. Finally, we performed a network analysis of the transcription factors expressed in the gastrointestinal tract and demonstrate how they sub-divide into functional groups that may control cellular gastrointestinal development. CONCLUSIONS: As an important livestock animal with a physiology that is more similar than mouse to man, we provide a major new resource for understanding gene expression with respect to the known physiology of mammalian tissues and cells. The data and analyses are available on the websites http://biogps.org and http://www.macrophages.com/pig-atlas.

Host cell maturation modulates parasite invasion and sexual differentiation in Plasmodium berghei.

Malaria remains a global health problem causing more than 400,000 deaths annually. Plasmodium parasites, the causative agents of malaria, replicate asexually in red blood cells (RBCs) of their vertebrate host, while a subset differentiates into sexual stages (gametocytes) for mosquito transmission. Parasite replication and gametocyte maturation in the erythropoietic niches of the bone marrow and spleen contribute to pathogenesis and drive transmission, but the mechanisms underlying this organ enrichment remain unknown. Here, we performed a comprehensive analysis of rodent P. berghei infection by flow cytometry and single-cell RNA sequencing. We identified CD71 as a host receptor for reticulocyte invasion and found that parasites metabolically adapt to the host cell environment. Transcriptional analysis and functional assays further revealed a nutrient-dependent tropism for gametocyte formation in reticulocytes. Together, we provide a thorough characterization of host-parasite interactions in erythropoietic niches and define host cell maturation state as the key driver of parasite adaptation.

Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population.

Understanding the genetic architecture of phenotypic variation in natural populations is a fundamental goal of evolutionary genetics. Wild Soay sheep (Ovis aries) have an inherited polymorphism for horn morphology in both sexes, controlled by a single autosomal locus, Horns. The majority of males have large normal horns, but a small number have vestigial, deformed horns, known as scurs; females have either normal horns, scurs or no horns (polled). Given that scurred males and polled females have reduced fitness within each sex, it is counterintuitive that the polymorphism persists within the population. Therefore, identifying the genetic basis of horn type will provide a vital foundation for understanding why the different morphs are maintained in the face of natural selection. We conducted a genome-wide association study using ∼36000 single nucleotide polymorphisms (SNPs) and determined the main candidate for Horns as RXFP2, an autosomal gene with a known involvement in determining primary sex characters in humans and mice. Evidence from additional SNPs in and around RXFP2 supports a new model of horn-type inheritance in Soay sheep, and for the first time, sheep with the same horn phenotype but different underlying genotypes can be identified. In addition, RXFP2 was shown to be an additive quantitative trait locus (QTL) for horn size in normal-horned males, accounting for up to 76% of additive genetic variation in this trait. This finding contrasts markedly from genome-wide association studies of quantitative traits in humans and some model species, where it is often observed that mapped loci only explain a modest proportion of the overall genetic variation.

Total parasite biomass but not peripheral parasitaemia is associated with endothelial and haematological perturbations in Plasmodium vivax patients.

Plasmodium vivax is the major cause of human malaria in the Americas. How P. vivax infection can lead to poor clinical outcomes, despite low peripheral parasitaemia, remains a matter of intense debate. Estimation of total P. vivax biomass based on circulating markers indicates existence of a predominant parasite population outside of circulation. In this study, we investigate associations between both peripheral and total parasite biomass and host response in vivax malaria. We analysed parasite and host signatures in a cohort of uncomplicated vivax malaria patients from Manaus, Brazil, combining clinical and parasite parameters, multiplexed analysis of host responses, and ex vivo assays. Patterns of clinical features, parasite burden, and host signatures measured in plasma across the patient cohort were highly heterogenous. Further data deconvolution revealed two patient clusters, here termed Vivaxlow and Vivaxhigh. These patient subgroups were defined based on differences in total parasite biomass but not peripheral parasitaemia. Overall Vivaxlow patients clustered with healthy donors and Vivaxhigh patients showed more profound alterations in haematological parameters, endothelial cell (EC) activation, and glycocalyx breakdown and levels of cytokines regulating different haematopoiesis pathways compared to Vivaxlow. Vivaxhigh patients presented more severe thrombocytopenia and lymphopenia, along with enrichment of neutrophils in the peripheral blood and increased neutrophil-to-lymphocyte ratio (NLCR). When patients' signatures were combined, high association of total parasite biomass with a subset of markers of EC activation, thrombocytopenia, and lymphopenia severity was observed. Finally, machine learning models defined a combination of host parameters measured in the circulation that could predict the extent of parasite infection outside of circulation. Altogether, our data show that total parasite biomass is a better predictor of perturbations in host homeostasis in P. vivax patients than peripheral parasitaemia. This supports the emerging paradigm of a P. vivax tissue reservoir, particularly in the haematopoietic niche of bone marrow and spleen.

Experimental evolution, genetic analysis and genome re-sequencing reveal the mutation conferring artemisinin resistance in an isogenic lineage of malaria parasites.

BACKGROUND: Classical and quantitative linkage analyses of genetic crosses have traditionally been used to map genes of interest, such as those conferring chloroquine or quinine resistance in malaria parasites. Next-generation sequencing technologies now present the possibility of determining genome-wide genetic variation at single base-pair resolution. Here, we combine in vivo experimental evolution, a rapid genetic strategy and whole genome re-sequencing to identify the precise genetic basis of artemisinin resistance in a lineage of the rodent malaria parasite, Plasmodium chabaudi. Such genetic markers will further the investigation of resistance and its control in natural infections of the human malaria, P. falciparum. RESULTS: A lineage of isogenic in vivo drug-selected mutant P. chabaudi parasites was investigated. By measuring the artemisinin responses of these clones, the appearance of an in vivo artemisinin resistance phenotype within the lineage was defined. The underlying genetic locus was mapped to a region of chromosome 2 by Linkage Group Selection in two different genetic crosses. Whole-genome deep coverage short-read re-sequencing (Illumina Solexa) defined the point mutations, insertions, deletions and copy-number variations arising in the lineage. Eight point mutations arise within the mutant lineage, only one of which appears on chromosome 2. This missense mutation arises contemporaneously with artemisinin resistance and maps to a gene encoding a de-ubiquitinating enzyme. CONCLUSIONS: This integrated approach facilitates the rapid identification of mutations conferring selectable phenotypes, without prior knowledge of biological and molecular mechanisms. For malaria, this model can identify candidate genes before resistant parasites are commonly observed in natural human malaria populations.

A strong quantitative trait locus for wing length on chromosome 2 in a wild population of great reed warblers.

Wing length is a key character for essential behaviours related to bird flight such as migration and foraging. In the present study, we initiate the search for the genes underlying wing length in birds by studying a long-distance migrant, the great reed warbler (Acrocephalus arundinaceus). In this species wing length is an evolutionary interesting trait with pronounced latitudinal gradient and sex-specific selection regimes in local populations. We performed a quantitative trait locus (QTL) scan for wing length in great reed warblers using phenotypic, genotypic, pedigree and linkage map data from our long-term study population in Sweden. We applied the linkage analysis mapping method implemented in GridQTL (a new web-based software) and detected a genome-wide significant QTL for wing length on chromosome 2, to our knowledge, the first detected QTL in wild birds. The QTL extended over 25 cM and accounted for a substantial part (37%) of the phenotypic variance of the trait. A genome scan for tarsus length (a body-size-related trait) did not show any signal, implying that the wing-length QTL on chromosome 2 was not associated with body size. Our results provide a first important step into understanding the genetic architecture of avian wing length, and give opportunities to study the evolutionary dynamics of wing length at the locus level.

Genome duplication in Leishmania major relies on persistent subtelomeric DNA replication.

DNA replication is needed to duplicate a cell's genome in S phase and segregate it during cell division. Previous work in Leishmania detected DNA replication initiation at just a single region in each chromosome, an organisation predicted to be insufficient for complete genome duplication within S phase. Here, we show that acetylated histone H3 (AcH3), base J and a kinetochore factor co-localise in each chromosome at only a single locus, which corresponds with previously mapped DNA replication initiation regions and is demarcated by localised G/T skew and G4 patterns. In addition, we describe previously undetected subtelomeric DNA replication in G2/M and G1-phase-enriched cells. Finally, we show that subtelomeric DNA replication, unlike chromosome-internal DNA replication, is sensitive to hydroxyurea and dependent on 9-1-1 activity. These findings indicate that Leishmania's genome duplication programme employs subtelomeric DNA replication initiation, possibly extending beyond S phase, to support predominantly chromosome-internal DNA replication initiation within S phase.

Gene mapping in the wild with SNPs: guidelines and future directions.

One of the biggest challenges facing evolutionary biologists is to identify and understand loci that explain fitness variation in natural populations. This review describes how genetic (linkage) mapping with single nucleotide polymorphism (SNP) markers can lead to great progress in this area. Strategies for SNP discovery and SNP genotyping are described and an overview of how to model SNP genotype information in mapping studies is presented. Finally, the opportunity afforded by new generation sequencing and typing technologies to map fitness genes by genome-wide association studies is discussed.

Feasibility and clinical utility of endoscopic ultrasound guided biopsy of pancreatic cancer for next-generation molecular profiling.

Next-generation sequencing is enabling molecularly guided therapy for many cancer types, yet failure rates remain relatively high in pancreatic cancer (PC). The aim of this study is to investigate the feasibility of genomic profiling using endoscopic ultrasound (EUS) biopsy samples to facilitate personalised therapy for PC. Ninty-five patients underwent additional research biopsies at the time of diagnostic EUS. Diagnostic formalin-fixed (FFPE) and fresh frozen EUS samples underwent DNA extraction, quantification and targeted gene sequencing. Whole genome (WGS) and RNA sequencing was performed as proof of concept. Only 2 patients (2%) with a diagnosis of PC had insufficient material for targeted sequencing in both FFPE and frozen specimens. Targeted panel sequencing (n=54) revealed mutations in PC genes (KRAS, GNAS, TP53, CDKN2A, SMAD4) in patients with histological evidence of PC, including potentially actionable mutations (BRCA1, BRCA2, ATM, BRAF). WGS (n=5) of EUS samples revealed mutational signatures that are potential biomarkers of therapeutic responsiveness. RNA sequencing (n=35) segregated patients into clinically relevant molecular subtypes based on transcriptome. Integrated multi-omic analysis of PC using standard EUS guided biopsies offers clinical utility to guide personalized therapy and study the molecular pathology in all patients with PC.

DNA G-quadruplex structures mold the DNA methylome.

Control of DNA methylation level is critical for gene regulation, and the factors that govern hypomethylation at CpG islands (CGIs) are still being uncovered. Here, we provide evidence that G-quadruplex (G4) DNA secondary structures are genomic features that influence methylation at CGIs. We show that the presence of G4 structure is tightly associated with CGI hypomethylation in the human genome. Surprisingly, we find that these G4 sites are enriched for DNA methyltransferase 1 (DNMT1) occupancy, which is consistent with our biophysical observations that DNMT1 exhibits higher binding affinity for G4s as compared to duplex, hemi-methylated, or single-stranded DNA. The biochemical assays also show that the G4 structure itself, rather than sequence, inhibits DNMT1 enzymatic activity. Based on these data, we propose that G4 formation sequesters DNMT1 thereby protecting certain CGIs from methylation and inhibiting local methylation.