The role of GC-biased gene conversion in shaping the fastest evolving regions of the human genome

GC-biased gene conversion (gBGC) is a recombination-associated evolutionary process that accelerates the fixation of guanine or cytosine alleles, regardless of their effects on fitness. gBGC can increase the overall rate of substitutions, a hallmark of positive selection. Many fast-evolving genes and noncoding sequences in the human genome have GC-biased substitution patterns, suggesting that gBGC-in contrast to adaptive processes-may have driven the human changes in these sequences. To investigate this hypothesis, we developed a substitution model for DNA sequence evolution that quantifies the nonlinear interacting effects of selection and gBGC on substitution rates and patterns. Based on this model, we used a series of lineage-specific likelihood ratio tests to evaluate sequence alignments for evidence of changes in mode of selection, action of gBGC, or both. With a false positive rate of less than 5% for individual tests, we found that the majority (76%) of previously identified human accelerated regions are best explained without gBGC, whereas a substantial minority (19%) are best explained by the action of gBGC alone. Further, more than half (55%) have substitution rates that significantly exceed local estimates of the neutral rate, suggesting that these regions may have been shaped by positive selection rather than by relaxation of constraint. By distinguishing the effects of gBGC, relaxation of constraint, and positive selection we provide an integrated analysis of the evolutionary forces that shaped the fastest evolving regions of the human genome, which facilitates the design of targeted functional studies of adaptation in humans.     

Single base-pair mutations in centromere element III cause aberrant chromosome segregation in Saccharomyces cerevisiae.

In this paper we show that a 211-base pair segment of CEN3 DNA is sufficient to confer wild-type centromere function in the yeast Saccharomyces cerevisiae. We used site-directed mutagenesis of the 211-base pair fragment to examine the sequence-specific functional requirements of a conserved 11-base pair segment of centromere DNA, element III (5'-TGATTTATCCGAA-3'). Element III is the most highly conserved of the centromeric DNA sequences, differing by only a single adenine X thymine base pair among the four centromere DNAs sequenced thus far. All of the element III sequences contain specific cytosine X guanine base pairs, including a 5'-CCG-3' arrangement, which we targeted for single cytosine-to-thymine mutations by using sodium bisulfite. The effects of element III mutations on plasmid and chromosome segregation were determined by mitotic stability assays. Conversion of CCG to CTG completely abolished centromere function both in plasmids and in chromosome III, whereas conversion of CCG to TCG decreased plasmid and chromosome stability moderately. The other two guanine X cytosine base pairs in element III could be independently converted to adenine X thymine base pairs without affecting plasmid or chromosome stability. We concluded that while some specific nucleotides within the conserved element III sequence are essential for proper centromere function, other conserved nucleotides can be changed.     

Evidence for strong selective constraint acting on the nucleotide composition of 16S ribosomal RNA genes

Previous studies have shown that the guanine plus cytosine (G+C) content of ribosomal RNAs (rRNAs) is highly correlated with bacterial growth temperatures. This correlation is strongest in the double-stranded stem regions of the rRNA, a fact that can be explained by selection for increased structural stability at high growth temperatures. In this study, we examined the single-stranded regions of 16S rRNAs. We reasoned that, since these regions of the molecule are subject to less structural constraint than the stem regions, their nucleotide content might simply reflect the overall nucleotide content of the genome. Contrary to this expectation, however, we found that all of the single-stranded regions are characterized by very high adenine (A) and relatively low cytosine (C) contents. Moreover, the nucleotide content of these single-stranded regions is surprisingly constant between species, despite dramatic differences in optimal growth temperatures, and despite large differences in the overall genomic G+C content. This provides compelling evidence for strong stabilizing selection acting on 16S rRNA single-stranded regions. We found that selection favors purines (A+G), and especially adenine (A), in the single-stranded regions of these rRNAs.     

Isolation and characterization of human DNA fragments with nucleotide sequence homologies with the simian virus 40 regulatory region.

A recombinant library of human DNA sequences was screened with a segment of simian virus 40 (SV40) DNA that spans the viral origin of replication. One hundred and fifty phage were isolated that hybridized to this probe. Restriction enzyme and hybridization analyses indicated that these sequences were partially homologous to one another. Direct DNA sequencing of two such SV40-hybridizing segments indicated that this was not a highly conserved family of sequences, but rather a set of DNA fragments that contained repetitive regions of high guanine plus cytosine content. These sequences were not members of the previously described Alu family of repeats and hybridized to SV40 DNA more strongly than do Alu family members. Computer analyses showed that the human DNA segments contained multiple homologies with sequences throughout the SV40 origin region, although sequences on the late side of the viral origin contained the strongest cross-hybridizing sequences. Because of the number and complexity of the matches detected, we could not determine unambiguously which of the many possible heteroduplexes between these DNAs was thermodynamically most favored. No hybridization of these human DNA sequences to any other segment of the SV40 genome was detected. In contrast, the human DNA segments isolated cross-hybridized with many sequences within the human genome. We tested for the presence of several functional domains on two of these human DNA fragments. One SV40-hybridizing fragment, SVCR29, contained a sequence which enhanced the efficiency of thymidine kinase transformation in human cells by approximately 20-fold. This effect was seen in an orientation-independent manner when the sequence was present at the 3' end of the chicken thymidine kinase gene. We propose that this segment of DNA contains a sequence analogous to the 72-base-pair repeats of SV40. The existence of such an "activator" element in cellular DNA raises the possibility that families of these sequences may exist in the mammalian genome.     

Analysis of human accelerated DNA regions using archaic hominin genomes

Several previous comparisons of the human genome with other primate and vertebrate genomes identified genomic regions that are highly conserved in vertebrate evolution but fast-evolving on the human lineage. These human accelerated regions (HARs) may be regions of past adaptive evolution in humans. Alternatively, they may be the result of non-adaptive processes, such as biased gene conversion. We captured and sequenced DNA from a collection of previously published HARs using DNA from an Iberian Neandertal. Combining these new data with shotgun sequence from the Neandertal and Denisova draft genomes, we determine at least one archaic hominin allele for 84% of all positions within HARs. We find that 8% of HAR substitutions are not observed in the archaic hominins and are thus recent in the sense that the derived allele had not come to fixation in the common ancestor of modern humans and archaic hominins. Further, we find that recent substitutions in HARs tend to have come to fixation faster than substitutions elsewhere in the genome and that substitutions in HARs tend to cluster in time, consistent with an episodic rather than a clock-like process underlying HAR evolution. Our catalog of sequence changes in HARs will help prioritize them for functional studies of genomic elements potentially responsible for modern human adaptations.     

Exponential decay of GC content detected by strand-symmetric substitution rates influences the evolution of isochore structure

The distribution of guanine and cytosine nucleotides throughout a genome, or the GC content, is associated with numerous features in mammals; understanding the pattern and evolutionary history of GC content is crucial to our efforts to annotate the genome. The local GC content is decaying toward an equilibrium point, but the causes and rates of this decay, as well as the value of the equilibrium point, remain topics of debate. By comparing the results of 2 methods for estimating local substitution rates, we identify 620 Mb of the human genome in which the rates of the various types of nucleotide substitutions are the same on both strands. These strand-symmetric regions show an exponential decay of local GC content at a pace determined by local substitution rates. DNA segments subjected to higher rates experience disproportionately accelerated decay and are AT rich, whereas segments subjected to lower rates decay more slowly and are GC rich. Although we are unable to draw any conclusions about causal factors, the results support the hypothesis proposed by Khelifi A, Meunier J, Duret L, and Mouchiroud D (2006. GC content evolution of the human and mouse genomes: insights from the study of processed pseudogenes in regions of different recombination rates. J Mol Evol. 62:745-752.) that the isochore structure has been reshaped over time. If rate variation were a determining factor, then the current isochore structure of mammalian genomes could result from the local differences in substitution rates. We predict that under current conditions strand-symmetric portions of the human genome will stabilize at an average GC content of 30% (considerably less than the current 42%), thus confirming that the human genome has not yet reached equilibrium.     

Human genome organization: Alu, lines, and the molecular structure of metaphase chromosome bands

Combining high resolution in situ hybridization with quantitative solid state imaging, we show that human metaphase chromosome Giemsa/Quinacrine and Reverse bands are each characterized by distinct families of interspersed repeated sequences: the SINES, Alu family dominates in Reverse bands, and the LINES, L1 family dominates in Giemsa/Quinacrine positive bands. Alu is 56% guanine plus cytosine, and L1 is 58% adenine plus thymine, and each may comprise 13%-18% of the total DNA in a chromosome band. Therefore, the distribution of these sequences alone may account for a large part of human chromosome banding seen with fluorescent dyes. With the exception of some telomeric regions, and the chromosomal regions of simple sequence DNA, Alu and L1 are precisely inversely distributed, suggesting an inverse functional relationship. This finding links genome organization with chromosome structure and function.     

Compilation and analysis of Bacillus subtilis sigma A-dependent promoter sequences: evidence for extended contact between RNA polymerase and upstream promoter DNA.

Sequence analysis of 236 promoters recognized by the Bacillus subtilis sigma A-RNA polymerase reveals an extended promoter structure. The most highly conserved bases include the -35 and -10 hexanucleotide core elements and a TG dinucleotide at position -15, -14. In addition, several weakly conserved A and T residues are present upstream of the -35 region. Analysis of dinucleotide composition reveals A2- and T2-rich sequences in the upstream promoter region (-36 to -70) which are phased with the DNA helix: An tracts are common near -43, -54 and -65; Tn tracts predominate at the intervening positions. When compared with larger regions of the genome, upstream promoter regions have an excess of An and Tn sequences for n > 4. These data indicate that an RNA polymerase binding site affects DNA sequence as far upstream as -70. This sequence conservation is discussed in light of recent evidence that the alpha subunits of the polymerase core bind DNA and that the promoter may wrap around RNA polymerase.     

Evidence for widespread degradation of gene control regions in hominid genomes

Although sequences containing regulatory elements located close to protein-coding genes are often only weakly conserved during evolution, comparisons of rodent genomes have implied that these sequences are subject to some selective constraints. Evolutionary conservation is particularly apparent upstream of coding sequences and in first introns, regions that are enriched for regulatory elements. By comparing the human and chimpanzee genomes, we show here that there is almost no evidence for conservation in these regions in hominids. Furthermore, we show that gene expression is diverging more rapidly in hominids than in murids per unit of neutral sequence divergence. By combining data on polymorphism levels in human noncoding DNA and the corresponding human–chimpanzee divergence, we show that the proportion of adaptive substitutions in these regions in hominids is very low. It therefore seems likely that the lack of conservation and increased rate of gene expression divergence are caused by a reduction in the effectiveness of natural selection against deleterious mutations because of the low effective population sizes of hominids. This has resulted in the accumulation of a large number of deleterious mutations in sequences containing gene control elements and hence a widespread degradation of the genome during the evolution of humans and chimpanzees.     

Nucleotide sequences of a feline leukemia virus subgroup A envelope gene and long terminal repeat and evidence for the recombinational origin of subgroup B viruses.

Molecular clones of the subgroup A feline leukemia virus FeLV-A/Glasgow-1 have been obtained. Nucleotide sequence analysis of the 3' end of the proviral genome and comparison with the published sequence of FeLV-B/Gardner-Arnstein showed that the most extensive differences are located within the 5' domain of the env gene. Within this domain, several divergent regions of env are separated by more conserved segments. The 3' end of env is highly conserved, with only a single amino acid coding difference in p15env. The proviral long terminal repeats are also highly conserved, differing by only eight base substitutions and one base insertion. Specific probes constructed from the FeLV-A or FeLV-B env genes were used to compare the env genes of various exogenous FeLV isolates and the endogenous FeLV-related proviruses of normal cat DNA. An FeLV-A-derived env probe showed no hybridization to normal cat DNA but detected all FeLV-A and FeLV-C isolates tested. In contrast, an FeLV-B env probe detected independent FeLV-B isolates and a family of endogenous FeLV-related proviruses. Our observations provide strong evidence to support the hypothesis that FeLV-B viruses have arisen by recombination between FeLV-A and endogenous proviral elements in cat DNA.     

Low fidelity DNA synthesis by a y family DNA polymerase due to misalignment in the active site

Sulfolobus solfataricus DNA polymerase IV (Dpo4) is a member of the Y family of DNA polymerases whose crystal structure has recently been solved. As a model for other evolutionarily conserved Y family members that perform translesion DNA synthesis and have low fidelity, we describe here the base substitution and frameshift fidelity of DNA synthesis by Dpo4. Dpo4 generates all 12 base-base mismatches at high rates, 11 of which are similar to those of its human homolog, DNA polymerase kappa. This result is consistent with the Dpo4 structure, implying lower geometric selection for correct base pairs. Surprisingly, Dpo4 generates C.dCMP mismatches at an unusually high average rate and preferentially at cytosine flanked by 5'-template guanine. Dpo4 also has very low frameshift fidelity and frequently generates deletions of even noniterated nucleotides, especially cytosine flanked by a 5'-template guanine. Both unusual features of error specificity suggest that Dpo4 can incorporate dNTP precursors when two template nucleotides are present in the active site binding pocket. These results have implications for mutagenesis resulting from DNA synthesis by Y family polymerases.     

The human genome: a view from above.

The genome of vertebrates (and of eukaryotes in general) is not simply formed by genes that are randomly scattered over vast expanses of "junk DNA", but is organized in a system which obeys precise rules, that amount to a genomic code. Moreover, genes are concentrated in the chromosomal regions which are the richest in G (guanine) and C (cytosine) and seem to correspond to the telomeric regions of certain chromosome arms (T-bands). The study of the genome organization in different vertebrate classes allowed us to approach in a novel way a number of fundamental problems of genome evolution.     

Score functions for determining regional conservation in two-species local alignments.

We construct several score functions for use in locating unusually conserved regions in a genomewide search of aligned DNA from two species. We test these functions on regions of the human genome aligned to the mouse genome. These score functions are derived from properties of neutrally evolving sites on the mouse and human genome and can be adjusted to the local background rate of conservation. The aim of these functions is to try to identify regions of the human genome that are conserved by evolutionary selection because they have an important function, rather than by chance. We use them to get a very rough estimate of the amount of DNA in the human genome that is under selection.     

Repetitive sequences in complete and defective genomes of Herpesvirus saimiri.

Two types of Herpesvirus saimiri genomes can be isolated from purified virions: (i) the M genome is a double-stranded, liniear DNA molecule with a mean contour length corresponding to 89 times 10-6 daltons. The M genome contains about 70% of unique sequences (light DNA, 36% guanine plus cytosine) and 30% reiterated sequences (heavy DNA, 71% guanine plus cytosine). (ii) the H genome is composed of heavy DNA only and is more heterogeneous in size. The sequences in the H genome are up to 40-fold reiterated, indicating defectiveness of this type of genome. The repetitions in the H genome and the M genome cross-hybridize almost completely and have identical kinetic complexity (2.8 times 10-6 daltons). DNA infectivity studies by using the calcium phosphate and the DEAE-dextran method gave further evidence that H genomes are defective: no infectious virus was recovered from permissive cells treated with heavy DNA, whereas M genome-infected cells developed cytopathic changes after 11 to 56 days. Defective H genomes were present in the progeny virus two passages after transfection.