Conservation and divergence in multigene families: alternatives to selection and drift

It is generally assumed that conservation and divergence of DNA signify function (selection) and no function (drift), respectively. This assumption is based on the view that a mutation is a unique event on a single chromosome, the fate of which depends on selection or drift. Knowledge of the rates, units and biases of widespread mechanisms of non-reciprocal DNA exchange, in particular within multigene families, provides alternative explanations for conservation and divergence, notwithstanding biological function. Such mechanisms of DNA turnover cause continual fluctuations in the copy-number of variant genes in an individual and, hence, promote the gradual and cohesive spread of a variant gene throughout a family (homogenization) and throughout a population (fixation). The dual processes (molecular drive) of homogenization and fixation are inextricably linked. Data are presented of the expected stages of transition in the spread of variant repeats by molecular drive in some non-genic families of DNA, seemingly not under the influence of selection. When a molecularly driven change in a given gene family is accompanied by the coevolution (mediated by selection) of other DNA, RNA or protein molecules that interact with the gene family then biological function is observed to be maintained despite sequence divergence. Conversely, the mechanics of DNA turnover and a turnover bias in favour of ancestral sequences can dramatically retard the rate of sequence change, in the absence of function. Examples of the maintenance of function by molecular coevolution and conservation of sequences in the absence of function, are drawn mainly from the rDNA multigene family.     

Rates of molecular evolution and the fraction of nucleotide positions free to vary

Summary Selective constraints on DNA sequence change were incorporated into a model of DNA divergence by restricting substitutions to a subset of nucleotide positions. A simple model showed that both mutation rate and the fraction of nucleotide positions free to vary are strong determinants of DNA divergence over time.When divergence between two species approaches the fraction of positions free to vary, standard methods that correct for multiple mutations yield severe underestimates of the number of substitutions per site. A modified method appropriate for use with DNA sequence, restriction site, or thermal renaturation data is derived taking this fraction into account. The model also showed that the ratio of divergence in two gene classes (e.g., nuclear and mitochondrial) may vary widely over time even if the ratio of mutation rates remains constant.DNA sequence divergence data are used increasingly to detect differences in rates of molecular evolution. Often, variation in divergence rate is assumed to represent variation in mutation rate. The present model suggests that differing divergence rates among comparisons (either among gene classes or taxa) should be interpreted cautiously. Differences in the fraction of nucleotide positions free to vary can serve as an important alternative hypothesis to explain differences in DNA divergence rates.     

Phylogenetic estimates of diversification rate are affected by molecular rate variation

Molecular phylogenies are increasingly being used to investigate the patterns and mechanisms of macroevolution. In particular, node heights in a phylogeny can be used to detect changes in rates of diversification over time. Such analyses rest on the assumption that node heights in a phylogeny represent the timing of diversification events, which in turn rests on the assumption that evolutionary time can be accurately predicted from DNA sequence divergence. But there are many influences on the rate of molecular evolution, which might also influence node heights in molecular phylogenies, and thus affect estimates of diversification rate. In particular, a growing number of studies have revealed an association between the net diversification rate estimated from phylogenies and the rate of molecular evolution. Such an association might, by influencing the relative position of node heights, systematically bias estimates of diversification time. We simulated the evolution of DNA sequences under several scenarios where rates of diversification and molecular evolution vary through time, including models where diversification and molecular evolutionary rates are linked. We show that commonly used methods, including metric‐based, likelihood and Bayesian approaches, can have a low power to identify changes in diversification rate when molecular substitution rates vary. Furthermore, the association between the rates of speciation and molecular evolution rate can cause the signature of a slowdown or speedup in speciation rates to be lost or misidentified. These results suggest that the multiple sources of variation in molecular evolutionary rates need to be considered when inferring macroevolutionary processes from phylogenies.     

Evolution of the genomic recombination rate in murid rodents

Although very closely related species can differ in their fine-scale patterns of recombination hotspots, variation in the average genomic recombination rate among recently diverged taxa has rarely been surveyed. We measured recombination rates in eight species that collectively represent several temporal scales of divergence within a single rodent family, Muridae. We used a cytological approach that enables in situ visualization of crossovers at meiosis to quantify recombination rates in multiple males from each rodent group. We uncovered large differences in genomic recombination rate between rodent species, which were independent of karyotypic variation. The divergence in genomic recombination rate that we document is not proportional to DNA sequence divergence, suggesting that recombination has evolved at variable rates along the murid phylogeny. Additionally, we document significant variation in genomic recombination rate both within and between subspecies of house mice. Recombination rates estimated in F(1) hybrids reveal evidence for sex-linked loci contributing to the evolution of recombination in house mice. Our results provide one of the first detailed portraits of genomic-scale recombination rate variation within a single mammalian family and demonstrate that the low recombination rates in laboratory mice and rats reflect a more general reduction in recombination rate across murid rodents.     

Orchestration of base excision repair by controlling the rates of enzymatic activities

Base excision repair (BER) is one of the major pathways for repair of simple DNA base lesions and is carried out through a series of coordinated reactions relying on several different enzymatic activities and accessory proteins. Imbalance of BER activities has been reported to be linked to genetic instability and cancer. To experimentally address the mechanisms orchestrating BER, we monitored both the overall rate and the rate-limiting steps in the repair in cell-free extracts of five different endogenously occurring DNA lesions (abasic site, uracil, 8-oxoguanine, hypoxanthine and 5,6-dihydrouracil) and the effect of addition of rate-limiting BER components on the rate and co-ordination of BER reactions. We find that several mechanisms including regulation of DNA glycosylase turnover and involvement of poly(ADP-ribose) polymerase participate in synchronization of the repair events. We also find that repair of different DNA lesions involves different mechanisms for optimizing repair rates without accumulation of intermediates. Repair of some lesions such as 8-oxoguanine is regulated by glycosylase turnover and progress without substantial accumulation of repair intermediates. However, during repair of the apurinic/apyrimidinic (AP) sites or 5,6-dihydrouracil, poly(ADP-ribose) polymerase plays an important role in the coordination of the rates of repair reactions.     

Studies on DNA sequences in the Osmundaceae

Summary Phylogenetic relationships ofOsmunda cinnamomea, O. claytoniana, andO. regalis were explored by means of DNA sequence comparisons. Hydroxyapatite thermal elution profiles of self-reassociated repetitive DNA fragments were very similar, indicating the absence of gross differences in the amount of recent amplification or addition of repetitive DNA in any of these three genomes. Interspecific DNA sequence comparisons showed, in contrast to our earlier interpretation, that repeated DNA sequences ofO. claytoniana are nearly equally diverged from those ofO. cinnamomea andO. regalis. Differences between repetitive sequences of the three species can be interpreted as reflecting amplification events which occurred subsequent to speciation. The data obtained suggest that the threeOsmunda species most likely arose more or less simultaneously from a common ancestor. These findings were verified in experiments with tracer DNA preparations enriched for single copy sequences. On the basis of the hybridization data presented here and of the fossil record, the rate of single copy sequence divergence in the ferns is comparable to that in the primates, although slower than that observed in other animal taxa. From this first evaluation of rates of DNA evolution in plants it would seem that the rates for plants and animals are roughly comparable. The evidence suggests that species divergence is accompanied by further reiteration of preexisting repeat sequences. The rate of addition of repetitive sequences probably is slower in ferns than in angiosperms. This difference might be attributable to the much larger effective generation time in ferns.     

Interhomologue sequence variation of alpha satellite DNA from human chromosome 17: Evidence for concerted evolution along haplotypic lineages

Alpha satellite DNA is a family of tandemly repeated DNA found at the centromeres of all primate chromosomes. Different human chromosomes 17 in the population are characterized by distinct alpha satellite haplotypes, distinguished by the presence of variant repeat forms that have precise monomeric deletions. Pairwise comparisons of sequence diversity between variant repeat units from each haplotype show that they are closely related in sequence. Direct sequencing of PCR-amplified alpha satellite reveals heterogeneous positions between the repeat units on a chromosome as two bands at the same position on a sequencing ladder. No variation was detected in the sequence and location of these heterogeneous positions between chromosomes 17 from the same haplotype, but distinct patterns of variation were detected between chromosomes from different haplotypes. Subsequent sequence analysis of individual repeats from each haplotype confirmed the presence of extensive haplotype-specific sequence variation. Phylogenetic inference yielded a tree that suggests these chromosome 17 repeat units evolve principally along haplotypic lineages. These studies allow insight into the relative rates and/or timing of genetic turnover processes that lead to the homogenization of tandem DNA families. Correspondence to: H.F. Willard     

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

Differential rates of nucleotide substitution among different gene segments and between distinct evolutionary lineages is well documented among mitochondrial genes and is likely a consequence of locus-specific selective constraints that delimit mutational divergence over evolutionary time. We compared sequence variation of 18 homologous loci (15 coding genes and 3 parts of the control region) among 10 mammalian mitochondrial DNA genomes which allowed us to describe different mitochondrial evolutionary patterns and to produce an estimation of the relative order of gene divergence. The relative rates of divergence of mitochondrial DNA genes in the family Felidae were estimated by comparing their divergence from homologous counterpart genes included in nuclear mitochondrial DNA (Numt, pronounced "new might"), a genomic fossil that represents an ancient transfer of 7.9 kb of mitochondrial DNA to the nuclear genome of an ancestral species of the domestic cat (Felis catus). Phylogenetic analyses of mitochondrial (mtDNA) sequences with multiple outgroup species were conducted to date the ancestral node common to the Numt and the cytoplasmic (Cymt) mtDNA genes and to calibrate the rate of sequence divergence of mitochondrial genes relative to nuclear homologous counterparts. By setting the fastest substitution rate as strictly mutational, an empirical "selective retardation index" is computed to quantify the sum of all constraints, selective and otherwise, that limit sequence divergence of mitochondrial gene sequences over time.      

Nucleotide sequence of a highly repetitive component of rat DNA.

A highly repetitive component of rat DNA which could not yet be enriched by density gradient centrifugation was isolated with the help of the restriction nuclease Sau3AI. This nuclease converted the bulk of the DNA to small fragments and left a repetitive DNA component as large fragments which were subsequently purified by gel filtration and electrophoresis. This DNA component which was termed rat satellite DNA I is composed of tandemly repeated 370 bp blocks. According to sequence analysis the 370 bp repeats consist of alternating 92 and 93 bp units with homologous but not identical sequences. Methylation of CpG residues was correlated to the rate of cleavage by restriction nucleases. Significant homologies exist between the sequences of rat satellite DNA I and satellite DNAs of several other organisms. The divergence of the sequence of rat satellite DNA I was discussed with respect to evolutionary considerations.     

Evolution of R1 and R2 in the rDNA units of the genus Drosophila

R1 and R2 are non-long terminal repeat (non-LTR) retrotransposable elements that specifically insert in the 28S ribosomal RNA (rRNA) genes of insects. Using the Drosophila genus, which includes some of the best characterized insect taxa, we have conducted a number of studies on the evolution of these elements. We find that R1 and R2 are subject to the same recombinational forces that give rise to the concerted evolution of the rDNA units. The turnover of R1 and R2 elements can be readily documented in different strains of D. melanogaster using 5′ truncated elements as restriction-length polymorphisms. This turnover leads to uniform populations of elements with nucleotide sequence divergence of different copies averaging only 0.23% for the R2 and 0.47% for the R1 elements. Molecular phylogenetic analysis of elements from 16 different species of Drosophila suggests that these elements have been stable components of the rDNA locus for the 50–70 million year history of the Drosophila genus. Using changes at synonymous positions within the protein-encoding regions as estimates of the baseline substitution rate, it could be shown that R1 and R2 are evolving at rates similar to that of typical protein encoding genes provided corrections are made for the low codon bias of the elements. R1 and R2 are clearly well-adapted for their existence in the rDNA units of their host. This revised version was published online in August 2006 with corrections to the Cover Date.     

Intercontinental and intracontinental biogeography patterns and methods

The study of biogeography has benefited from the exponential increase of DNA sequence data from recent molecular systematic studies, the development of analytical methods in the last decade concerning divergence time estimation and geographic area analyses, and the availability of large-scale distributiofi data of species in many groups of organisms. The underlying principle of divergence time estimation from DNA and protein data is that sequence divergence depends on the product of evolutionary rate and time. With their molecular clock hypothesis, Zuckerkandl and Pauling （1965） separated rates of molecular evolution from time by incorporating fossil evidence. Originally,     

Insights into the history of a bacterial group II intron remnant from the genomes of the nitrogen-fixing symbionts Sinorhizobium meliloti and Sinorhizobium medicae

Group II introns are self-splicing catalytic RNAs that act as mobile retroelements. In bacteria, they are thought to be tolerated to some extent because they self-splice and home preferentially to sites outside of functional genes, generally within intergenic regions or in other mobile genetic elements, by mechanisms including the divergence of DNA target specificity to prevent target site saturation. RmInt1 is a mobile group II intron that is widespread in natural populations of Sinorhizobium meliloti and was first described in the GR4 strain. Like other bacterial group II introns, RmInt1 tends to evolve toward an inactive form by fragmentation, with loss of the 3′ terminus. We identified genomic evidence of a fragmented intron closely related to RmInt1 buried in the genome of the extant S. meliloti/S. medicae species. By studying this intron, we obtained evidence for the occurrence of intron insertion before the divergence of ancient rhizobial species. This fragmented group II intron has thus existed for a long time and has provided sequence variation, on which selection can act, contributing to diverse genetic rearrangements, and to generate pan-genome divergence after strain differentiation. The data presented here suggest that fragmented group II introns within intergenic regions closed to functionally important neighboring genes may have been microevolutionary forces driving adaptive evolution of these rhizobial species.     

New insights into the evolutionary rate of HIV-1 at the within-host and epidemiological levels

Over calendar time, HIV-1 evolves considerably faster within individuals than it does at the epidemic level. This is a surprising observation since, from basic population genetic theory, we would expect the genetic substitution rate to be similar across different levels of biological organization. Three different mechanisms could potentially cause the observed mismatch in phylogenetic rates of divergence: temporal changes in selection pressure during the course of infection; frequent reversion of adaptive mutations after transmission; and the storage of the virus in the body followed by the preferential transmission of stored ancestral virus. We evaluate each of these mechanisms to determine whether they are likely to make a major contribution to the mismatch in phylogenetic rates. We conclude that the cycling of the virus through very long-lived memory CD4(+) T cells, a process that we call 'store and retrieve', is probably the major contributing factor to the rate mismatch. The preferential transmission of ancestral virus needs to be integrated into evolutionary models if we are to accurately predict the evolution of immune escape, drug resistance and virulence in HIV-1 at the population level. Moreover, early infection viruses should be the major target for vaccine design, because these are the viral strains primarily involved in transmission.     

Diversity in a chorion multigene family created by tandem duplications and a putative gene-conversion event

Summary Two families of high-cysteine chorion proteins inBombyx mori are encoded in 15 tandemly arranged nonidentical gene pairs. It is assumed that this locus arose by duplication with subsequent sequence divergence. We have compared DNA sequences from two such neighboring pairs of genes in an attempt to understand the manner in which diversity has been generated and/or removed. A high level of sequence identity (91%–99%) was found between the repeats throughout the transcribed and flanking regions, with two significant exceptions. First, in the DNA segment encoding a conserved region of the chorion proteins, ten substitutions were detected in a 39-base-pair region. This localized region of high variability would suggest an intergene conversion-like event. Second, a length difference of 141 base pairs was detected in a region encoding the carboxy-terminal arm of the protein. This difference can be explained by three separate reiterations of single codons (3 base pairs) separated in time by duplication or triplication events.     

Related base sequences in the DNA of simple and complex organisms. VI. The extent of base sequence divergence among the DNAs of various rodents

The magnitude of differences which occur in the base sequence of DNAs from related animals was investigated through studies of four rodent DNAs. Fractionation on preparative CsCl density gradients and detailed examination of thermal denaturation profiles revealed the existence of several distinguishable components comprising the main band in addition to satellite bands. No difference was apparent in the distribution of unique and partially redundant sequences among these main band components, or in the apparent rate at which base substitution has occurred in evolution. Comparisons of sequence were also made by annealing reactions in solution at high concentration for extended periods with either labeled total mouse or rat DNA, or labeled unique sequences of mouse DNA. These experiments permitted assessment of the absolute extent of cross-reaction, as well as the extent of sequence divergence which had occurred between the two groups of interacting DNA molecules. Some implications of the apparent rapidity of nucleotide substitution in DNA are discussed.This research was supported by a grant from the National Science Foundation (GB 6099).     

Exploring the correlations between sequence evolution rate and phenotypic divergence across the Mammalian tree provides insights into adaptive evolution

Sequence evolution behaves in a relatively consistent manner, leading to one of the fundamental paradigms in biology, the existence of a ??molecular clock??. The molecular clock can be distilled to the concept of accumulation of substitutions, through time yielding a stable rate from which we can estimate lineage divergence. Over the last 50?years, evolutionary biologists have obtained an in-depth understanding of this clock??s nuances. It has been fine-tuned by taking into account the vast heterogeneity in rates across lineages and genes, leading to ??relaxed?? molecular clock methods for timetree reconstruction. Sequence rate varies with life history traits including body size, generation time and metabolic rate, and we review recent studies on this topic. However, few studies have explicitly examined correlates between molecular evolution and morphological evolution. The patterns observed across diverse lineages suggest that rates of molecular and morphological evolution are largely decoupled. We discuss how identifying the molecular mechanisms behind rapid functional radiations are central to understanding evolution. The vast functional divergence within mammalian lineages that have relatively ??slow?? sequence evolution refutes the hypotheses that pulses in diversification yielding major phenotypic change are the result of steady accumulation of substitutions. Patterns rather suggest phenotypic divergence is likely caused by regulatory alterations mediated through mechanisms such as insertions/deletions in functional regions. These can rapidly arise and sweep to fixation faster than predicted from a lineage??s sequence neutral substitution rate, enabling species to leapfrog between phenotypic ??islands??. We suggest research directions that could illuminate mechanisms behind the functional diversity we see today.     

Rates of genomic divergence in humans,chimpanzees and their lice

The rate of DNA mutation and divergence is highly variable across the tree of life. However, the reasons underlying this variation are not well understood. Comparing the rates of genetic changes between hosts and parasite lineages that diverged at the same time is one way to begin to understand differences in genetic mutation and substitution rates. Such studies have indicated that the rate of genetic divergence in parasites is often faster than that of their hosts when comparing single genes. However, the variation in this relative rate of molecular evolution across different genes in the genome is unknown. We compared the rate of DNA sequence divergence between humans, chimpanzees and their ectoparasitic lice for 1534 protein-coding genes across their genomes. The rate of DNA substitution in these orthologous genes was on average 14 times faster for lice than for humans and chimpanzees. In addition, these rates were positively correlated across genes. Because this correlation only occurred for substitutions that changed the amino acid, this pattern is probably produced by similar functional constraints across the same genes in humans, chimpanzees and their ectoparasites.     

Evolution of the cetacean mitochondrial D-loop region.

We sequenced the mitochondrial DNA D-loop regions from two cetacean species and compared these with the published D-loop sequences of several other mammalian species, including one other cetacean. Nucleotide substitution rates, DNA sequence simplicity, possible open reading frames (ORFs), and potential RNA secondary structure were investigated. The substitution rate is an order of magnitude lower than would be expected on the basis of reports on human sequence variation in this region but are consistent with interspecific primate and rodent D-loop sequence variation and with estimates of substitution rates from whole mitochondrial genomes. Deletions/insertions are less common in the cetacean D-loop than in other vertebrate species. Areas of high sequence simplicity (clusters of short repetitive motifs) across the region correspond to areas of high sequence divergence. Three regions predicted to form secondary structures are homologous to such putative structures in other species; however, the presumptive structures most conserved in cetaceans are different from those reported for other taxa. While all three species have possible long ORFs, only a short sequence of seven amino acids is shared with other mammalian species, and those changes that had occurred within it are all nonsynonymous. We conclude that DNA slippage, in addition to point mutation, contributes to the evolution of the D-loop and that regions of conserved secondary structure in cetaceans and an ORF are unlikely to contribute significantly to the conservation of the central region.     

Transition stages of molecular drive in multiple-copy DNA families in Drosophila

Multigene and non-genic DNA families are in a state of turnover and hence are continually being replaced throughout a population by new variant repeats. To quantify such molecular processes, in the absence of selection, it is necessary to find and compare stages of transistion during the homogenization of at least two non-genic families evolving in parallel in a closely related group of species. Detailed sequence analysis of patterns of variation, at each nucleotide position considered independently, amongst repeats of two tandem DNA families from seven related Drosophila species, reveals all stages of transition during the spread of randomly produced variant repeats. Variant repeats are found at different stages of homogenization and fixation in a population, irrespective of the loci, chromosomes or individuals from which they were cloned. Differences between the families in the relatively small number of variants at each transition stage and the greater number of fully homogenized and fixed variants between species of greater divergence indicate that the process of spread (molecular drive) is rapid relative to the mutation rate and occurs at seemingly different constant rates for each family. Occasional gene conversions, in addition to unequal exchanges, have contributed to family turnover. The significance of these results to the evolution of functional multigene families and divergence and conservation of sequences is discussed.