Origin of Genes

We discuss two tests of the hypothesis that the first genes were assembled from exons. The hypothesis of exon shuffling in the progenote predicts that intron phases will be correlated so that exons will be an integer number of codons and predicts that the exons will be correlated with compact regions of polypeptide chain. These predictions have been tested on ancient conserved proteins (proteins without introns in prokaryotes but with introns in eukaryotes) and hold with high statistical significance. We conclude that introns are correlated with compact features of proteins 15-, 22-, or 30-amino acid residues long, as was predicted by “The Exon Theory of Genes.”

Transposable elements as sources of variation in animals and plants

A tremendous wealth of data is accumulating on the variety and distribution of transposable elements (TEs) in natural populations. There is little doubt that TEs provide new genetic variation on a scale, and with a degree of sophistication, previously unimagined. There are many examples of mutations and other types of genetic variation associated with the activity of mobile elements. Mutant phenotypes range from subtle changes in tissue specificity to dramatic alterations in the development and organization of tissues and organs. Such changes can occur because of insertions in coding regions, but the more sophisticated TE-mediated changes are more often the result of insertions into 5′ flanking regions and introns. Here, TE-induced variation is viewed from three evolutionary perspectives that are not mutually exclusive. First, variation resulting from the intrinsic parasitic nature of TE activity is examined. Second, we describe possible coadaptations between elements and their hosts that appear to have evolved because of selection to reduce the deleterious effects of new insertions on host fitness. Finally, some possible cases are explored in which the capacity of TEs to generate variation has been exploited by their hosts. The number of well documented cases in which element sequences appear to confer useful traits on the host, although small, is growing rapidly.

Long term trends in the evolution of H(3) HA1 human influenza type A

We have studied the HA1 domain of 254 human influenza A(H3N2) virus genes for clues that might help identify characteristics of hemagglutinins (HAs) of circulating strains that are predictive of that strain’s epidemic potential. Our preliminary findings include the following. (i) The most parsimonious tree found requires 1,260 substitutions of which 712 are silent and 548 are replacement substitutions. (ii) The HA1 portion of the HA gene is evolving at a rate of 5.7 nucleotide substitutions/year or 5.7 × 10−3 substitutions/site per year. (iii) The replacement substitutions are distributed randomly across the three positions of the codon when allowance is made for the number of ways each codon can change the encoded amino acid. (iv) The replacement substitutions are not distributed randomly over the branches of the tree, there being 2.2 times more changes per tip branch than for non-tip branches. This result is independent of how the virus was amplified (egg grown or kidney cell grown) prior to sequencing or if sequencing was carried out directly on the original clinical specimen by PCR. (v) These excess changes on the tip branches are probably the result of a bias in the choice of strains to sequence and the detection of deleterious mutations that had not yet been removed by negative selection. (vi) There are six hypervariable codons accumulating replacement substitutions at an average rate that is 7.2 times that of the other varied codons. (vii) The number of variable codons in the trunk branches (the winners of the competitive race against the immune system) is 47 ± 5, significantly fewer than in the twigs (90 ± 7), which in turn is significantly fewer variable codons than in tip branches (175 ± 8). (viii) A minimum of one of every 12 branches has nodes at opposite ends representing viruses that reside on different continents. This is, however, no more than would be expected if one were to randomly reassign the continent of origin of the isolates. (ix) Of 99 codons with at least four mutations, 31 have ratios of non-silent to silent changes with probabilities less than 0.05 of occurring by chance, and 14 of those have probabilities <0.005. These observations strongly support positive Darwinian selection. We suggest that the small number of variable positions along the successful trunk lineage, together with knowledge of the codons that have shown positive selection, may provide clues that permit an improved prediction of which strains will cause epidemics and therefore should be used for vaccine production.

Genes, peoples, and languages

The genetic history of a group of populations is usually analyzed by reconstructing a tree of their origins. Reliability of the reconstruction depends on the validity of the hypothesis that genetic differentiation of the populations is mostly due to population fissions followed by independent evolution. If necessary, adjustment for major population admixtures can be made. Dating the fissions requires comparisons with paleoanthropological and paleontological dates, which are few and uncertain. A method of absolute genetic dating recently introduced uses mutation rates as molecular clocks; it was applied to human evolution using microsatellites, which have a sufficiently high mutation rate. Results are comparable with those of other methods and agree with a recent expansion of modern humans from Africa. An alternative method of analysis, useful when there is adequate geographic coverage of regions, is the geographic study of frequencies of alleles or haplotypes. As in the case of trees, it is necessary to summarize data from many loci for conclusions to be acceptable. Results must be independent from the loci used. Multivariate analyses like principal components or multidimensional scaling reveal a number of hidden patterns and evaluate their relative importance. Most patterns found in the analysis of human living populations are likely to be consequences of demographic expansions, determined by technological developments affecting food availability, transportation, or military power. During such expansions, both genes and languages are spread to potentially vast areas. In principle, this tends to create a correlation between the respective evolutionary trees. The correlation is usually positive and often remarkably high. It can be decreased or hidden by phenomena of language replacement and also of gene replacement, usually partial, due to gene flow.

DNA variation at the Sod locus of Drosophila melanogaster: An unfolding story of natural selection

Patterns of variation at the Sod locus of Drosophila melanogaster suggest that the protein polymorphism at this locus has very recently arisen. In addition, it appears that a previously rare DNA variant has been recently and rapidly driven to intermediate frequency. From the size of the region (>20 kb) that has been swept along with this rare variant, and patterns of linkage disequilibrium in the region, it is inferred that strength of selection was large (s > 0.01) and that the sweep occurred more than 25,000 generations ago. In addition, there are striking similarities to patterns of variation observed at the Est6 and Est-P loci, which are located approximately 1,000 kb from Sod.

Neutral behavior of shared polymorphism

Several cases have been described in the literature where genetic polymorphism appears to be shared between a pair of species. Here we examine the distribution of times to random loss of shared polymorphism in the context of the neutral Wright–Fisher model. Order statistics are used to obtain the distribution of times to loss of a shared polymorphism based on Kimura’s solution to the diffusion approximation of the Wright–Fisher model. In a single species, the expected absorption time for a neutral allele having an initial allele frequency of ½ is 2.77 N generations. If two species initially share a polymorphism, that shared polymorphism is lost as soon as either of two species undergoes fixation. The loss of a shared polymorphism thus occurs sooner than loss of polymorphism in a single species and has an expected time of 1.7 N generations. Molecular sequences of genes with shared polymorphism may be characterized by the count of the number of sites that segregate in both species for the same nucleotides (or amino acids). The distribution of the expected numbers of these shared polymorphic sites also is obtained. Shared polymorphism appears to be more likely at genetic loci that have an unusually large number of segregating alleles, and the neutral coalescent proves to be very useful in determining the probability of shared allelic lineages expected by chance. These results are related to examples of shared polymorphism in the literature.

Divergent and conserved features in the spatial expression of the Drosophila pseudoobscura esterase-5B gene and the esterase-6 gene of Drosophila melanogaster

The regulatory regions of homologous genes encoding esterase 6 (Est-6) of Drosophila melanogaster and esterase 5B (Est-5B) of Drosophila pseudoobscura show very little similarity. We have undertaken a comparative study of the pattern of expression directed by the Est-5B and Est-6 5′-flanking DNA to attempt to reveal conserved elements regulating tissue-specific expression in adults. Esterase regulatory sequences were linked to a lacZ reporter gene and transformed into D. melanogaster embryos. Est-5B, 5′ upstream elements, give rise to a β-galactosidase expression pattern that coincides with the wild-type expression of Est-5B in D. pseudoobscura. The expression patterns of the Est-5B/lacZ construct are different from those of a fusion gene containing the upstream region of Est-6. Common sites of expression for both kinds of constructs are the third segment of antenna, the maxillary palps, and salivary glands. In vitro deletion mutagenesis has shown that the two genes have a different organization of regulatory elements controlling expression in both the third segment of antenna and maxillary palps. The results suggest that the conservation of the expression pattern in genes that evolved from a common ancestor may not be accompanied by preservation of the corresponding cis-regulatory elements.

A demographic approach to selection

The concepts of demography provide a means of combining the ecological approach to population growth with the genetical approach to natural selection. We have utilized the demographic theory of natural selection developed by Norton and Charlesworth to analyze life history schedules of births and deaths for populations of genotypes in Drosophila pseudoobscura. Our populations illustrate a stable genetic equilibrium, an unstable genetic equilibrium, and a case of no equilibrium. We have estimated population growth rates and Darwinian fitnesses of the genotypes and have explored the role of population growth in determining natural selection. The age-specific rates of births and deaths provide insights into components of selection. Both viability and fertility are important components in our populations.

Phylogenetics and the origin of species

A recent criticism that the biological species concept (BSC) unduly neglects phylogeny is examined under a novel modification of coalescent theory that considers multiple, sex-defined genealogical pathways through sexual organismal pedigrees. A competing phylogenetic species concept (PSC) also is evaluated from this vantage. Two analytical approaches are employed to capture the composite phylogenetic information contained within the braided assemblages of hereditary pathways of a pedigree: (i) consensus phylogenetic trees across allelic transmission routes and (ii) composite phenograms from quantitative values of organismal coancestry. Outcomes from both approaches demonstrate that the supposed sharp distinction between biological and phylogenetic species concepts is illusory. Historical descent and reproductive ties are related aspects of phylogeny and jointly illuminate biotic discontinuity.

Assortative fertilization in Drosophila

The concept of gametic isolation has its origins in the 1937 edition of T. Dobzhansky’s Genetics and the Origin of Species. Involving either positive assortative fertilization (as opposed to self-incompatibility) or negative assortative fertilization, it occurs after mating but prior to fertilization. Gametic isolation is generally subsumed under either prezygotic or postmating isolation and thus has not been the subject of extensive investigation. Examples of assortative fertilization in Drosophila are reviewed and compared with those of other organisms. Potential mechanisms leading to assortative fertilization are discussed, as are their evolutionary implications.

Incipient species formation in salamanders of the Ensatina complex

The Ensatina eschscholtzii complex of plethodontid salamanders, a well-known “ring species,” is thought to illustrate stages in the speciation process. Early research, based on morphology and coloration, has been extended by the incorporation of studies of protein variation and mitochondrial DNA sequences. The new data show that the complex includes a number of geographically and genetically distinct components that are at or near the species level. The complex is old and apparently has undergone instances of range contraction, isolation, differentiation, and then expansion and secondary contact. While the hypothesis that speciation is retarded by gene flow around the ring is not supported by molecular data, the general biogeographical hypothesis is supported. There is evidence of a north to south range expansion along two axes, with secondary contact and completion of the ring in southern California. Current research targets regions once thought to show primary intergradation, but which molecular markers reveal to be zones of secondary contact. Here emphasis is on the subspecies E. e. xanthoptica, which is involved in four distinct secondary contacts in central California. There is evidence of renewed genetic interactions upon recontact, with greater genetic differentiation within xanthoptica than between it and some of the interacting populations. The complex presents a full array of intermediate conditions between well-marked species and geographically variable populations. Geographically differentiated segments represent a diversity of depths of time of isolation and admixture, reflecting the complicated geomorphological history of California. Ensatina illustrates the continuing difficulty in making taxonomic assignments in complexes studied during species formation.

Genetics and the origin of bird species

External (environmental) factors affecting the speciation of birds are better known than the internal (genetic) factors. The opposite is true for several groups of invertebrates, Drosophila being the outstanding example. Ideas about the genetics of speciation in general trace back to Dobzhansky who worked with Drosophila. These ideas are an insufficient guide for reconstructing speciation in birds for two main reasons. First, speciation in birds proceeds with the evolution of behavioral barriers to interbreeding; postmating isolation usually evolves much later, perhaps after gene exchange has all but ceased. As a consequence of the slow evolution of postmating isolating factors the scope for reinforcement of premating isolation is small, whereas the opportunity for introgressive hybridization to influence the evolution of diverging species is large. Second, premating isolation may arise from nongenetic, cultural causes; isolation may be affected partly by song, a trait that is culturally inherited through an imprinting-like process in many, but not all, groups of birds. Thus the genetic basis to the origin of bird species is to be sought in the inheritance of adult traits that are subject to natural and sexual selection. Some of the factors involved in premating isolation (plumage, morphology, and behavior) are under single-gene control, most are under polygenic control. The genetic basis of the origin of postmating isolating factors affecting the early development of embryos (viability) and reproductive physiology (sterility) is almost completely unknown. Bird speciation is facilitated by small population size, involves few genetic changes, and occurs relatively rapidly.

Vagaries of the molecular clock

The hypothesis of the molecular evolutionary clock asserts that informational macromolecules (i.e., proteins and nucleic acids) evolve at rates that are constant through time and for different lineages. The clock hypothesis has been extremely powerful for determining evolutionary events of the remote past for which the fossil and other evidence is lacking or insufficient. I review the evolution of two genes, Gpdh and Sod. In fruit flies, the encoded glycerol-3-phosphate dehydrogenase (GPDH) protein evolves at a rate of 1.1 × 10−10 amino acid replacements per site per year when Drosophila species are compared that diverged within the last 55 million years (My), but a much faster rate of ≈4.5 × 10−10 replacements per site per year when comparisons are made between mammals (≈70 My) or Dipteran families (≈100 My), animal phyla (≈650 My), or multicellular kingdoms (≈1100 My). The rate of superoxide dismutase (SOD) evolution is very fast between Drosophila species (16.2 × 10−10 replacements per site per year) and remains the same between mammals (17.2) or Dipteran families (15.9), but it becomes much slower between animal phyla (5.3) and still slower between the three kingdoms (3.3). If we assume a molecular clock and use the Drosophila rate for estimating the divergence of remote organisms, GPDH yields estimates of 2,500 My for the divergence between the animal phyla (occurred ≈650 My) and 3,990 My for the divergence of the kingdoms (occurred ≈1,100 My). At the other extreme, SOD yields divergence times of 211 My and 224 My for the animal phyla and the kingdoms, respectively. It remains unsettled how often proteins evolve in such erratic fashion as GPDH and SOD.

Evolution of codon usage bias in Drosophila

We first review what is known about patterns of codon usage bias in Drosophila and make the following points: (i) Drosophila genes are as biased or more biased than those in microorganisms. (ii) The level of bias of genes and even the particular pattern of codon bias can remain phylogenetically invariant for very long periods of evolution. (iii) However, some genes, even very tightly linked genes, can change very greatly in codon bias across species. (iv) Generally G and especially C are favored at synonymous sites in biased genes. (v) With the exception of aspartic acid, all amino acids contribute significantly and about equally to the codon usage bias of a gene. (vi) While most individual amino acids that can use G or C at synonymous sites display a preference for C, there are exceptions: valine and leucine, which prefer G. (vii) Finally, smaller genes tend to be more biased than longer genes. We then examine possible causes of these patterns and discount mutation bias on three bases: there is little evidence of regional mutation bias in Drosophila, mutation bias is likely toward A+T (the opposite of codon usage bias), and not all amino acids display the preference for the same nucleotide in the wobble position. Two lines of evidence support a selection hypothesis based on tRNA pools: highly biased genes tend to be highly and/or rapidly expressed, and the preferred codons in highly biased genes optimally bind the most abundant isoaccepting tRNAs. Finally, we examine the effect of bias on DNA evolution and confirm that genes with high codon usage bias have lower rates of synonymous substitution between species than do genes with low codon usage bias. Surprisingly, we find that genes with higher codon usage bias display higher levels of intraspecific synonymous polymorphism. This may be due to opposing effects of recombination.

The evolution of plant nuclear genes

We analyze the evolutionary dynamics of three of the best-studied plant nuclear multigene families. The data analyzed derive from the genes that encode the small subunit of ribulose-1,5-bisphosphate carboxylase (rbcS), the gene family that encodes the enzyme chalcone synthase (Chs), and the gene family that encodes alcohol dehydrogenases (Adh). In addition, we consider the limited evolutionary data available on plant transposable elements. New Chs and rbcS genes appear to be recruited at about 10 times the rate estimated for Adh genes, and this is correlated with a much smaller average gene family size for Adh genes. In addition, duplication and divergence in function appears to be relatively common for Chs genes in flowering plant evolution. Analyses of synonymous nucleotide substitution rates for Adh genes in monocots reject a linear relationship with clock time. Replacement substitution rates vary with time in a complex fashion, which suggests that adaptive evolution has played an important role in driving divergence following gene duplication events. Molecular population genetic studies of Adh and Chs genes reveal high levels of molecular diversity within species. These studies also reveal that inter- and intralocus recombination are important forces in the generation allelic novelties. Moreover, illegitimate recombination events appear to be an important factor in transposable element loss in plants. When we consider the recruitment and loss of new gene copies, the generation of allelic diversity within plant species, and ectopic exchange among transposable elements, we conclude that recombination is a pervasive force at all levels of plant evolution.