Despite intense investigation, human replication origins and termini remain elusive. Existing data have shown strong discrepancies. Here we sequenced highly purified Okazaki fragments from two cell types and, for the first time, quantitated replication fork directionality and delineated initiation and termination zones genome-wide. Replication initiates stochastically, primarily within non-transcribed, broad (up to 150 kb) zones that often abut transcribed genes, and terminates dispersively between them. Replication fork progression is significantly co-oriented with the transcription. Initiation and termination zones are frequently contiguous, sometimes separated by regions of unidirectional replication. Initiation zones are enriched in open chromatin and enhancer marks, even when not flanked by genes, and often border ‘topologically associating domains’ (TADs). Initiation zones are enriched in origin recognition complex (ORC)-binding sites and better align to origins previously mapped using bubble-trap than λ-exonuclease. This novel panorama of replication reveals how chromatin and transcription modulate the initiation process to create cell-type-specific replication programs.

An esterase gene from the mosquito Culex quinquefasciatus that is responsible for resistance to a variety of organophosphorus (OP) insecticides was cloned in λgt11 phage. This gene was used to investigate the genetic mechanism of the high production of the esterase B1 it encodes in OP-resistant Culex quinquefasciatus Say (Tem-R strain) from California. Adults of the Tem-R strain were found to possess at least 250 times more copies of the gene than adults of a susceptible strain (S-Lab). The finding that selection by pesticides may result in the amplification of genes encoding detoxifying enzymes in whole, normally developed, reproducing insects emphasizes the biological importance of this mechanism and opens new areas of investigation in pesticide resistance management.

In early Xenopus embryos, in which ribosomal RNA genes (rDNA) are not transcribed, rDNA replication initiates and terminates at 9- to 12-kilobase pair intervals, with no detectable dependence on specific DNA sequences. Resumption of ribosomal RNA (rRNA) synthesis at late blastula and early gastrula is accompanied by a specific repression of replication initiation within transcription units; the frequency of initiation within intergenic spacers remains as high as in early blastula. These results demonstrate that for rRNA genes, circumscribed zones of replication initiation emerge in intergenic DNA during the time in metazoan development when the chromatin is remodeled to allow gene transcription.

Neutral nucleotide substitutions occur at varying rates along genomes, and it remains a major issue to unravel the mechanisms that cause these variations and to analyze their evolutionary consequences. Here, we study the role of replication in the neutral substitution pattern. We obtained a high-resolution replication timing profile of the whole human genome by massively parallel sequencing of nascent BrdU-labeled replicating DNA. These data were compared to the neutral substitution rates along the human genome, obtained by aligning human and chimpanzee genomes using macaque and orangutan as outgroups. All substitution rates increase monotonously with replication timing even after controlling for local or regional nucleotide composition, crossover rate, distance to telomeres, and chromatin compaction. The increase in non-CpG substitution rates might result from several mechanisms including the increase in mutation-prone activities or the decrease in efficiency of DNA repair during the S phase. In contrast, the rate of C → T transitions in CpG dinucleotides increases in later-replicating regions due to increasing DNA methylation level that reflects a negative correlation between timing and gene expression. Similar results are observed in the mouse, which indicates that replication timing is a main factor affecting nucleotide substitution dynamics at non-CpG sites and constitutes a major neutral process driving mammalian genome evolution.