Evolution of chlorophyll and bacteriochlorophyll: the problem of invariant sites in sequence analysis.

Competing hypotheses seek to explain the evolution of oxygenic and anoxygenic processes of photosynthesis. Since chlorophyll is less reduced and precedes bacteriochlorophyll on the modern biosynthetic pathway, it has been proposed that chlorophyll preceded bacteriochlorophyll in its evolution. However, recent analyses of nucleotide sequences that encode chlorophyll and bacteriochlorophyll biosynthetic enzymes appear to provide support for an alternative hypothesis. This is that the evolution of bacteriochlorophyll occurred earlier than the evolution of chlorophyll. Here we demonstrate that the presence of invariant sites in sequence datasets leads to inconsistency in tree building (including maximum-likelihood methods). Homologous sequences with different biological functions often share invariant sites at the same nucleotide positions. However, different constraints can also result in additional invariant sites unique to the genes, which have specific and different biological functions. Consequently, the distribution of these sites can be uneven between the different types of homologous genes. The presence of invariant sites, shared by related biosynthetic genes as well as those unique to only some of these genes, has misled the recent evolutionary analysis of oxygenic and anoxygenic photosynthetic pigments. We evaluate an alternative scheme for the evolution of chlorophyll and bacteriochlorophyll.

The exaptive excellence of spandrels as a term and prototype

In 1979, Lewontin and I borrowed the architectural term “spandrel” (using the pendentives of San Marco in Venice as an example) to designate the class of forms and spaces that arise as necessary byproducts of another decision in design, and not as adaptations for direct utility in themselves. This proposal has generated a large literature featuring two critiques: (i) the terminological claim that the spandrels of San Marco are not true spandrels at all and (ii) the conceptual claim that they are adaptations and not byproducts. The features of the San Marco pendentives that we explicitly defined as spandrel-properties—their necessary number (four) and shape (roughly triangular)—are inevitable architectural byproducts, whatever the structural attributes of the pendentives themselves. The term spandrel may be extended from its particular architectural use for two-dimensional byproducts to the generality of “spaces left over,” a definition that properly includes the San Marco pendentives. Evolutionary biology needs such an explicit term for features arising as byproducts, rather than adaptations, whatever their subsequent exaptive utility. The concept of biological spandrels—including the examples here given of masculinized genitalia in female hyenas, exaptive use of an umbilicus as a brooding chamber by snails, the shoulder hump of the giant Irish deer, and several key features of human mentality—anchors the critique of overreliance upon adaptive scenarios in evolutionary explanation. Causes of historical origin must always be separated from current utilities; their conflation has seriously hampered the evolutionary analysis of form in the history of life.

A mitochondrial DNA clone is associated with increased risk for Alzheimer disease.

Severe mitochondrial genetic mutations lead to early degeneration of specific human tissues; milder mitochondrial mutations may cause degeneration at a later point in life. A mutation at position 4336 was reported to occur at increased frequency in individuals with Alzheimer disease (AD) and Parkinson disease [Shoffner, J. M., Brown, M. D., Torroni, A., Lott, M. T., Cabell, M. F., Mirra, S. S., Beal, M. F., Yang, C.-C., Gearing, M., Salvo, R., Watts, R. L., Juncos, J. L., Hansen, L. A., Crain, B. J., Fayad, M., Reckord, C. L. & Wallace, D. C. (1993) Genomics 17, 171-184]. We have investigated the notion that this mutation leads to excess risk of AD by using a case-control study design of 72 AD autopsies and 296 race- and age-matched controls. The 4336G mutation occurred at higher frequency in AD autopsies than age-matched controls, a statistically significant difference. Evolutionary analysis of mtDNAs bearing the 4336G mutation indicated they were more closely related to each other than to other mtDNAs, consistent with the model of a single origin for this mutation. The tight evolutionary relatedness and homoplasmy of mtDNAs that confer elevated risk for a late-onset disease contrast strikingly with the distant relatedness and heteroplasmy of mitochondrial genomes that cause early-onset disease. The dichotomy can be explained by a lack of selection against mutations that confer a phenotype at advanced age during most of the evolution of humans. We estimate that approximately 1.5 million Caucasians in the United States bear the 4336G mutation and are at significantly increased risk of developing mitochondrial AD in their lifetime. A mechanism for 4336G-mediated cell death is proposed.

Lactone-ring-cleaving enzyme: Genetic analysis, novel RNA editing, and evolutionary implications

A lactonohydrolase from Fusarium oxysporum AKU 3702 is an enzyme catalyzing the hydrolysis of aldonate lactones to the corresponding aldonic acids. The amino acid sequences of the NH2 terminus and internal peptide fragments of the enzyme were determined to prepare synthetic oligonucleotides as primers for the PCR. An approximate 1,000-base genomic DNA fragment thus amplified was used as the probe to clone both genomic DNA and cDNA for the enzyme. The lactonohydrolase genomic gene consists of six exons separated by five short introns. A novel type of RNA editing, in which lactonohydrolase mRNA included the insertion of guanosine and cytidine residues, was observed. The predicted amino acid sequence of the cloned lactonohydrolase cDNA showed significant similarity to those of the gluconolactonase from Zymomonas mobilis, and paraoxonases from human and rabbit, forming a unique superfamily consisting of C-O cleaving enzymes and P-O cleaving enzymes. Lactonohydrolase was expressed under the control of the lac promoter in Escherichia coli.

Molecular phylogenetic analysis of evolutionary trends in stonefly wing structure and locomotor behavior

Insects in the order Plecoptera (stoneflies) use a form of two-dimensional aerodynamic locomotion called surface skimming to move across water surfaces. Because their weight is supported by water, skimmers can achieve effective aerodynamic locomotion even with small wings and weak flight muscles. These mechanical features stimulated the hypothesis that surface skimming may have been an intermediate stage in the evolution of insect flight, which has perhaps been retained in certain modern stoneflies. Here we present a phylogeny of Plecoptera based on nucleotide sequence data from the small subunit rRNA (18S) gene. By mapping locomotor behavior and wing structural data onto the phylogeny, we distinguish between the competing hypotheses that skimming is a retained ancestral trait or, alternatively, a relatively recent loss of flight. Our results show that basal stoneflies are surface skimmers, and that various forms of surface skimming are distributed widely across the plecopteran phylogeny. Stonefly wings show evolutionary trends in the number of cross veins and the thickness of the cuticle of the longitudinal veins that are consistent with elaboration and diversification of flight-related traits. These data support the hypothesis that the first stoneflies were surface skimmers, and that wing structures important for aerial flight have become elaborated and more diverse during the radiation of modern stoneflies.

Evolutionary EST analysis identifies rapidly evolving male reproductive proteins in Drosophila

Sequence comparisons of genomes or expressed sequence tags (ESTs) from related organisms provide insight into functional conservation and diversification. We compare the sequences of ESTs from the male accessory gland of Drosophila simulans to their orthologs in its close relative Drosophila melanogaster, and demonstrate rapid divergence of many of these reproductive genes. Nineteen (∼11%) of 176 independent genes identified in the EST screen contain protein-coding regions with an excess of nonsynonymous over synonymous changes, suggesting that their divergence has been accelerated by positive Darwinian selection. Genes that encode putative accessory gland-specific seminal fluid proteins had a significantly elevated level of nonsynonymous substitution relative to nonaccessory gland-specific genes. With the 57 new accessory gland genes reported here, we predict that ∼90% of the male accessory gland genes have been identified. The evolutionary EST approach applied here to identify putative targets of adaptive evolution is readily applicable to other tissues and organisms.

Sequence similarity analysis of Escherichia coli proteins: functional and evolutionary implications.

A computer analysis of 2328 protein sequences comprising about 60% of the Escherichia coli gene products was performed using methods for database screening with individual sequences and alignment blocks. A high fraction of E. coli proteins--86%--shows significant sequence similarity to other proteins in current databases; about 70% show conservation at least at the level of distantly related bacteria, and about 40% contain ancient conserved regions (ACRs) shared with eukaryotic or Archaeal proteins. For > 90% of the E. coli proteins, either functional information or sequence similarity, or both, are available. Forty-six percent of the E. coli proteins belong to 299 clusters of paralogs (intraspecies homologs) defined on the basis of pairwise similarity. Another 10% could be included in 70 superclusters using motif detection methods. The majority of the clusters contain only two to four members. In contrast, nearly 25% of all E. coli proteins belong to the four largest superclusters--namely, permeases, ATPases and GTPases with the conserved "Walker-type" motif, helix-turn-helix regulatory proteins, and NAD(FAD)-binding proteins. We conclude that bacterial protein sequences generally are highly conserved in evolution, with about 50% of all ACR-containing protein families represented among the E. coli gene products. With the current sequence databases and methods of their screening, computer analysis yields useful information on the functions and evolutionary relationships of the vast majority of genes in a bacterial genome. Sequence similarity with E. coli proteins allows the prediction of functions for a number of important eukaryotic genes, including several whose products are implicated in human diseases.

The optimization principle in phylogenetic analysis tends to give incorrect topologies when the number of nucleotides or amino acids used is small

In the maximum parsimony (MP) and minimum evolution (ME) methods of phylogenetic inference, evolutionary trees are constructed by searching for the topology that shows the minimum number of mutational changes required (M) and the smallest sum of branch lengths (S), respectively, whereas in the maximum likelihood (ML) method the topology showing the highest maximum likelihood (A) of observing a given data set is chosen. However, the theoretical basis of the optimization principle remains unclear. We therefore examined the relationships of M, S, and A for the MP, ME, and ML trees with those for the true tree by using computer simulation. The results show that M and S are generally greater for the true tree than for the MP and ME trees when the number of nucleotides examined (n) is relatively small, whereas A is generally lower for the true tree than for the ML tree. This finding indicates that the optimization principle tends to give incorrect topologies when n is small. To deal with this disturbing property of the optimization principle, we suggest that more attention should be given to testing the statistical reliability of an estimated tree rather than to finding the optimal tree with excessive efforts. When a reliability test is conducted, simplified MP, ME, and ML algorithms such as the neighbor-joining method generally give conclusions about phylogenetic inference very similar to those obtained by the more extensive tree search algorithms.

A structural view of evolutionary divergence

Two directed evolution experiments on p-nitrobenzyl esterase yielded one enzyme with a 100-fold increased activity in aqueous-organic solvents and another with a 17°C increase in thermostability. Structures of the wild type and its organophilic and thermophilic counterparts are presented at resolutions of 1.5 Å, 1.6 Å, and 2.0 Å, respectively. These structures identify groups of interacting mutations and demonstrate how directed evolution can traverse complex fitness landscapes. Early-generation mutations stabilize flexible loops not visible in the wild-type structure and set the stage for further beneficial mutations in later generations. The mutations exert their influence on the esterase structure over large distances, in a manner that would be difficult to predict. The loops with the largest structural changes generally are not the sites of mutations. Similarly, none of the seven amino acid substitutions in the organophile are in the active site, even though the enzyme experiences significant changes in the organization of this site. In addition to reduction of surface loop flexibility, thermostability in the evolved esterase results from altered core packing, helix stabilization, and the acquisition of surface salt bridges, in agreement with other comparative studies of mesophilic and thermophilic enzymes. Crystallographic analysis of the wild type and its evolved counterparts reveals networks of mutations that collectively reorganize the active site. Interestingly, the changes that led to diversity within the α/β hydrolase enzyme family and the reorganization seen in this study result from main-chain movements.

Character displacement in some Cnemidophorus lizards revisited: A phylogenetic analysis

Ecological studies have demonstrated the role of competition in structuring communities; however, the importance of competition as a vehicle for evolution by natural selection and speciation remains unresolved. Study systems of insular faunas have provided several well known cases where ecological character displacement, coevolution of competitors leading to increased morphological separation, is thought to have occurred (e.g., anoline lizards and geospizine finches). Whiptail lizards (genus Cnemidophorus) from the islands of the Sea of Cortez and the surrounding mainland demonstrate a biogeographic pattern of morphological variation suggestive of character displacement. Two species of Cnemidophorus occur on the Baja peninsula, one relatively large (Cnemidophorus tigris) and one smaller (Cnemidophorus hyperythrus). Oceanic islands in the Sea of Cortez contain only single species, five of six having sizes intermediate to both species found on the Baja peninsula. On mainland Mexico C. hyperythrus is absent, whereas C. tigris is the smaller species in whiptail guilds. Here we construct a phylogeny using nucleotide sequences of the cytochrome b gene to infer the evolutionary history of body size change and historical patterns of colonization in the Cnemidophorus system. The phylogenetic analysis indicates that (i) oceanic islands have been founded at least five times from mainland sources by relatives of either C. tigris or C. hyperythrus, (ii) there have been two separate instances of character relaxation on oceanic islands for C. tigris, and (iii) there has been colonization of the oceanic island Cerralvo with retention of ancestral size for Cnemidophorus ceralbensis, a relative of C. hyperythrus. Finally, the phylogenetic analysis reveals potential cryptic species within mainland populations of C. tigris.

Fluorescence in situ hybridization analysis of keratinocyte growth factor gene amplification and dispersion in evolution of great apes and humans

Keratinocyte growth factor (KGF) is a member of the fibroblast growth factor family. Portions of the gene encoding KGF were amplified during primate evolution and are present in multiple nonprocessed copies in the human genome. Nucleotide analysis of a representative sampling of these KGF-like sequences indicated that they were at least 95% identical to corresponding regions of the KGF gene. To localize these sequences to specific chromosomal sites in human and higher primates, we used fluorescence in situ hybridization. In human, using a cosmid probe encoding KGF exon 1, we assigned the location of the KGF gene to chromosome 15q15–21.1. In addition, copies of KGF-like sequences hybridizing only with a cosmid probe encoding exons 2 and 3 were localized to dispersed sites on chromosome 2q21, 9p11, 9q12–13, 18p11, 18q11, 21q11, and 21q21.1. The distribution of KGF-like sequences suggests a role for alphoid DNA in their amplification and dispersion. In chimpanzee, KGF-like sequences were observed at five chromosomal sites, which were each homologous to sites in human, while in gorilla, a subset of four of these homologous sites was identified; in orangutan two sites were identified, while gibbon exhibited only a single site. The chromosomal localization of KGF sequences in human and great ape genomes indicates that amplification and dispersion occurred in multiple discrete steps, with initial KGF gene duplication and dispersion taking place in gibbon and involving loci corresponding to human chromosomes 15 and 21. These findings support the concept of a closer evolutionary relationship of human and chimpanzee and a possible selective pressure for such dispersion during the evolution of higher primates.

Molecular dynamics of MHC genesis unraveled by sequence analysis of the 1,796,938-bp HLA class I region

The intensely studied MHC has become the paradigm for understanding the architectural evolution of vertebrate multigene families. The 4-Mb human MHC (also known as the HLA complex) encodes genes critically involved in the immune response, graft rejection, and disease susceptibility. Here we report the continuous 1,796,938-bp genomic sequence of the HLA class I region, linking genes between MICB and HLA-F. A total of 127 genes or potentially coding sequences were recognized within the analyzed sequence, establishing a high gene density of one per every 14.1 kb. The identification of 758 microsatellite provides tools for high-resolution mapping of HLA class I-associated disease genes. Most importantly, we establish that the repeated duplication and subsequent diversification of a minimal building block, MIC-HCGIX-3.8–1-P5-HCGIV-HLA class I-HCGII, engendered the present-day MHC. That the currently nonessential HLA-F and MICE genes have acted as progenitors to today’s immune-competent HLA-ABC and MICA/B genes provides experimental evidence for evolution by “birth and death,” which has general relevance to our understanding of the evolutionary forces driving vertebrate multigene families.

Phylogenetic analysis of “Volvocacae” for comparative genetic studies

Sequence analysis based on multiple isolates representing essentially all genera and species of the classic family Volvocaeae has clarified their phylogenetic relationships. Cloned internal transcribed spacer sequences (ITS-1 and ITS-2, flanking the 5.8S gene of the nuclear ribosomal gene cistrons) were aligned, guided by ITS transcript secondary structural features, and subjected to parsimony and neighbor joining distance analysis. Results confirm the notion of a single common ancestor, and Chlamydomonas reinharditii alone among all sequenced green unicells is most similar. Interbreeding isolates were nearest neighbors on the evolutionary tree in all cases. Some taxa, at whatever level, prove to be clades by sequence comparisons, but others provide striking exceptions. The morphological species Pandorina morum, known to be widespread and diverse in mating pairs, was found to encompass all of the isolates of the four species of Volvulina. Platydorina appears to have originated early and not to fall within the genus Eudorina, with which it can sometimes be confused by morphology. The four species of Pleodorina appear variously associated with Eudorina examples. Although the species of Volvox are each clades, the genus Volvox is not. The conclusions confirm and extend prior, more limited, studies on nuclear SSU and LSU rDNA genes and plastid-encoded rbcL and atpB. The phylogenetic tree suggests which classical taxonomic characters are most misleading and provides a framework for molecular studies of the cell cycle-related and other alterations that have engendered diversity in both vegetative and sexual colony patterns in this classical family.