Intron self-complementarity enforces exon inclusion in a yeast pre-mRNA

Skipping of internal exons during removal of introns from pre-mRNA must be avoided for proper expression of most eukaryotic genes. Despite significant understanding of the mechanics of intron removal, mechanisms that ensure inclusion of internal exons in multi-intron pre-mRNAs remain mysterious. Using a natural two-intron yeast gene, we have identified distinct RNA–RNA complementarities within each intron that prevent exon skipping and ensure inclusion of internal exons. We show that these complementarities are positioned to act as intron identity elements, bringing together only the appropriate 5′ splice sites and branchpoints. Destroying either intron self-complementarity allows exon skipping to occur, and restoring the complementarity using compensatory mutations rescues exon inclusion, indicating that the elements act through formation of RNA secondary structure. Introducing new pairing potential between regions near the 5′ splice site of intron 1 and the branchpoint of intron 2 dramatically enhances exon skipping. Similar elements identified in single intron yeast genes contribute to splicing efficiency. Our results illustrate how intron secondary structure serves to coordinate splice site pairing and enforce exon inclusion. We suggest that similar elements in vertebrate genes could assist in the splicing of very large introns and in the evolution of alternative splicing.

Ancient mtDNA sequences in the human nuclear genome: A potential source of errors in identifying pathogenic mutations

Nuclear-localized mtDNA pseudogenes might explain a recent report describing a heteroplasmic mtDNA molecule containing five linked missense mutations dispersed over the contiguous mtDNA CO1 and CO2 genes in Alzheimer’s disease (AD) patients. To test this hypothesis, we have used the PCR primers utilized in the original report to amplify CO1 and CO2 sequences from two independent ρ° (mtDNA-less) cell lines. CO1 and CO2 sequences amplified from both of the ρ° cells, demonstrating that these sequences are also present in the human nuclear DNA. The nuclear pseudogene CO1 and CO2 sequences were then tested for each of the five “AD” missense mutations by restriction endonuclease site variant assays. All five mutations were found in the nuclear CO1 and CO2 PCR products from ρ° cells, but none were found in the PCR products obtained from cells with normal mtDNA. Moreover, when the overlapping nuclear CO1 and CO2 PCR products were cloned and sequenced, all five missense mutations were found, as well as a linked synonymous mutation. Unlike the findings in the original report, an additional 32 base substitutions were found, including two in adjacent tRNAs and a two base pair deletion in the CO2 gene. Phylogenetic analysis of the nuclear CO1 and CO2 sequences revealed that they diverged from modern human mtDNAs early in hominid evolution about 770,000 years before present. These data would be consistent with the interpretation that the missense mutations proposed to cause AD may be the product of ancient mtDNA variants preserved as nuclear pseudogenes.

A possible mitochondrial gene in the early-branching amitochondriate protist Trichomonas vaginalis

Trichomonads are anaerobic flagellated protists that, based on analyses of ribosomal RNA sequences, represent one of the earliest branching lineages among the eukaryotes. The absence of mitochondria in these organisms coupled with their deep phylogenetic position has prompted several authors to suggest that trichomonads, along with other deeply-branching amitochondriate protist groups, diverged from the main eukaryotic lineage prior to the endosymbiotic origin of mitochondria. In this report we describe the presence of a gene in Trichomonas vaginalis specifically related to mitochondrial chaperonin 60 (cpn60). A recent study indicates that a protein immunologically related to cpn60 is located in trichomonad hydrogenosomes. Together, these data provide evidence that ancestors of trichomonads perhaps harbored the endosymbiotic progenitors of mitochondria, but that these evolved into hydrogenosomes early in trichomonad evolution.

High relatedness and inbreeding at the origin of eusociality in gall-inducing thrips

Within the haplodiploid eusocial gall-inducing thrips, a species-level phylogeny combined with genetic data for five eusocial species enables an inference of levels of relatedness and inbreeding values for lineages at the origin of eusociality. Character optimization using data from five eusocial species indicates that the lineage or lineages where eusociality is inferred to have originated exhibit relatedness of 0.64–0.92, and FIS of 0.33–0.64. The high inbreeding coefficients found in these eusocial thrips have increased relatedness among and within both sexes and have reduced the haplodiploidy-induced relatedness asymmetries [Hamilton, W. D. (1964) J. Theor. Biol. 7, 1–52]. These results indicate that unusually high relatedness is associated with the origin of eusociality, and they suggest a role for inbreeding in the evolution of bisexual helping.

Dynamic expression of multiple scavenger receptor cysteine-rich genes in coelomocytes of the purple sea urchin

Coelomocytes, the heterogeneous population of sea urchin putative immune cells, were found to express a complex set of transcripts featuring scavenger receptor cysteine-rich (SRCR) repeats. SRCR domains define a metazoan superfamily of proteins, many of which are implicated in development and regulation of the immune system of vertebrates. Coelomocytes transcribe multiple SRCR genes from among a multigene family encoding an estimated number of 1,200 SRCR domains in specific patterns particular to each individual. Transcription levels for given SRCR genes may range from pronounced to undetectable, yet all tested animals harbor the genomic loci encoding these genes. Analysis of several SRCR genes revealed multiple loci corresponding to each type. In the case of one SRCR type, a cluster of at least three genes was detected within a 133-kb bacterial artificial chromosome insert, and conserved as well as unique regions were identified in sequences of three genomic clones derived from a single animal. Array hybridizations with repeated samples of coelomocyte messages revealed substantial alterations in levels of expression of many SRCR genes, with fluctuations of up to 10-fold in 1 week and up to 30-fold over a period of 3 months. This report is the first demonstration of genomic and transcriptional complexity in molecules expressed by invertebrate coelomocytes. The mechanisms controlling SRCR gene expression and the functional significance of this dynamic system await elucidation.

Isoprenoid biosynthesis: The evolution of two ancient and distinct pathways across genomes

Isopentenyl diphosphate (IPP) is the central intermediate in the biosynthesis of isoprenoids, the most ancient and diverse class of natural products. Two distinct routes of IPP biosynthesis occur in nature: the mevalonate pathway and the recently discovered deoxyxylulose 5-phosphate (DXP) pathway. The evolutionary history of the enzymes involved in both routes and the phylogenetic distribution of their genes across genomes suggest that the mevalonate pathway is germane to archaebacteria, that the DXP pathway is germane to eubacteria, and that eukaryotes have inherited their genes for IPP biosynthesis from prokaryotes. The occurrence of genes specific to the DXP pathway is restricted to plastid-bearing eukaryotes, indicating that these genes were acquired from the cyanobacterial ancestor of plastids. However, the individual phylogenies of these genes, with only one exception, do not provide evidence for a specific affinity between the plant genes and their cyanobacterial homologues. The results suggest that lateral gene transfer between eubacteria subsequent to the origin of plastids has played a major role in the evolution of this pathway.

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.

Characterization of MADS homeotic genes in the fern Ceratopteris richardii

The MADS genes encode a family of transcription factors, some of which control the identities of floral organs in flowering plants. To understand the role of MADS genes in the evolution of floral organs, five MADS genes (CMADS1, 2, 3, 4, and 6) were cloned from the fern Ceratopteris richardii, a nonflowering plant. A gene tree of partial amino acid sequences of seed plant and fern MADS genes showed that the fern genes form three subfamilies. All members of one of the fern MADS subfamilies have additional amino-terminal amino acids, which is a synapomorphic character of the AGAMOUS subfamily of the flowering plant MADS genes. Their structural similarity indicates a sister relationship between the two subfamilies. The temporal and spatial patterns of expression of the five fern MADS genes were assessed by Northern blot analyses and in situ hybridizations. CMADS1, 2, 3, and 4 are expressed similarly in the meristematic regions and primordia of sporophyte shoots and roots, as well as in reproductive structures, including sporophylls and sporangial initials, although the amount of expression in each tissue is different in each gene. CMADS6 is expressed in gametophytic tissues but not in sporophytic tissues. The lack of organ-specific expression of MADS genes in the reproductive structures of the fern sporophyte may indicate that the restriction of MADS gene expression to specific reproductive organs and the specialization of MADS gene functions as homeotic selector genes in the flowering plant lineage were important in floral organ evolution.

A method for directed evolution and functional cloning of enzymes

A general scheme is described for the in vitro evolution of protein catalysts in a biologically amplifiable system. Substrate is covalently and site specifically attached by a flexible tether to the pIII coat protein of a filamentous phage that also displays the catalyst. Intramolecular conversion of substrate to product provides a basis for selecting active catalysts from a library of mutants, either by release from or attachment to a solid support. This methodology has been developed with the enzyme staphylococcal nuclease as a model. An analysis of factors influencing the selection efficiency is presented, and it is shown that phage displaying staphylococcal nuclease can be enriched 100-fold in a single step from a library-like ensemble of phage displaying noncatalytic proteins. Additionally, this approach should allow one to functionally clone natural enzymes, based on their ability to catalyze specific reactions (e.g., glycosyl transfer, sequence-specific proteolysis or phosphorylation, polymerization, etc.) rather than their sequence- or structural homology to known enzymes.

Molecular evolution of two vertebrate aryl hydrocarbon (dioxin) receptors (AHR1 and AHR2) and the PAS family

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor through which halogenated aromatic hydrocarbons such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) cause altered gene expression and toxicity. The AHR belongs to the basic helix–loop–helix/Per-ARNT-Sim (bHLH-PAS) family of transcriptional regulatory proteins, whose members play key roles in development, circadian rhythmicity, and environmental homeostasis; however, the normal cellular function of the AHR is not yet known. As part of a phylogenetic approach to understanding the function and evolutionary origin of the AHR, we sequenced the PAS homology domain of AHRs from several species of early vertebrates and performed phylogenetic analyses of these AHR amino acid sequences in relation to mammalian AHRs and 24 other members of the PAS family. AHR sequences were identified in a teleost (the killifish Fundulus heteroclitus), two elasmobranch species (the skate Raja erinacea and the dogfish Mustelus canis), and a jawless fish (the lamprey Petromyzon marinus). Two putative AHR genes, designated AHR1 and AHR2, were found both in Fundulus and Mustelus. Phylogenetic analyses indicate that the AHR2 genes in these two species are orthologous, suggesting that an AHR gene duplication occurred early in vertebrate evolution and that multiple AHR genes may be present in other vertebrates. Database searches and phylogenetic analyses identified four putative PAS proteins in the nematode Caenorhabditis elegans, including possible AHR and ARNT homologs. Phylogenetic analysis of the PAS gene family reveals distinct clades containing both invertebrate and vertebrate PAS family members; the latter include paralogous sequences that we propose have arisen by gene duplication early in vertebrate evolution. Overall, our analyses indicate that the AHR is a phylogenetically ancient protein present in all living vertebrate groups (with a possible invertebrate homolog), thus providing an evolutionary perspective to the study of dioxin toxicity and AHR function.

Forced unfolding modulated by disulfide bonds in the Ig domains of a cell adhesion molecule

Cell adhesion molecules (CAMs) mediate cell attachment and stress transfer through extracellular domains. Here we forcibly unfold the Ig domains of a prototypical Ig superfamily CAM that contains intradomain disulfide bonds. The Ig domains of all such CAMs have conformations homologous to cadherin extracellular domains, titin Ig-type domains, and fibronectin type-III (FNIII) domains. Atomic force microscopy has been used to extend the five Ig domains of Mel-CAM (melanoma CAM)—a protein that is overexpressed in metastatic melanomas—under conditions where the disulfide bonds were either left intact or disrupted through reduction. Under physiological conditions where intradomain disulfide bonds are intact, partial unfolding was observed at forces far smaller than those reported previously for either titin's Ig-type domains or tenascin's FNIII domains. This partial unfolding under low force may be an important mechanism for imparting elasticity to cell–cell contacts, as well as a regulatory mechanism for adhesive interactions. Under reducing conditions, Mel-CAM's Ig domains were found to fully unfold through a partially folded state and at slightly higher forces. The results suggest that, in divergent evolution of all such domains, stabilization imparted by disulfide bonds relaxes requirements for strong, noncovalent, folded-state interactions.

Computational method to reduce the search space for directed protein evolution

We introduce a computational method to optimize the in vitro evolution of proteins. Simulating evolution with a simple model that statistically describes the fitness landscape, we find that beneficial mutations tend to occur at amino acid positions that are tolerant to substitutions, in the limit of small libraries and low mutation rates. We transform this observation into a design strategy by applying mean-field theory to a structure-based computational model to calculate each residue's structural tolerance. Thermostabilizing and activity-increasing mutations accumulated during the experimental directed evolution of subtilisin E and T4 lysozyme are strongly directed to sites identified by using this computational approach. This method can be used to predict positions where mutations are likely to lead to improvement of specific protein properties.

Directed evolution of polymerase function by compartmentalized self-replication

We describe compartmentalized self-replication (CSR), a strategy for the directed evolution of enzymes, especially polymerases. CSR is based on a simple feedback loop consisting of a polymerase that replicates only its own encoding gene. Compartmentalization serves to isolate individual self-replication reactions from each other. In such a system, adaptive gains directly (and proportionally) translate into genetic amplification of the encoding gene. CSR has applications in the evolution of polymerases with novel and useful properties. By using three cycles of CSR, we obtained variants of Taq DNA polymerase with 11-fold higher thermostability than the wild-type enzyme or with a >130-fold increased resistance to the potent inhibitor heparin. Insertion of an extra stage into the CSR cycle before the polymerase reaction allows its application to enzymes other than polymerases. We show that nucleoside diphosphate kinase and Taq polymerase can form such a cooperative CSR cycle based on reciprocal catalysis, whereby nucleoside diphosphate kinase produces the substrates required for the replication of its own gene. We also find that in CSR the polymerase genes themselves evolve toward more efficient replication. Thus, polymerase genes and their encoded polypeptides cooperate to maximize postselection copy number. CSR should prove useful for the directed evolution of enzymes, particularly DNA or RNA polymerases, as well as for the design and study of in vitro self-replicating systems mimicking prebiotic evolution and viral replication.

The evolution of RNA viruses: A population genetics view

RNA viruses are excellent experimental models for studying evolution under the theoretical framework of population genetics. For a proper justification of this thesis we have introduced some properties of RNA viruses that are relevant for studying evolution. On the other hand, population genetics is a reductionistic theory of evolution. It does not consider or make simplistic assumptions on the transformation laws within and between genotypic and phenotypic spaces. However, such laws are minimized in the case of RNA viruses because the phenotypic space maps onto the genotypic space in a much more linear way than on higher DNA-based organisms. Under experimental conditions, we have tested the role of deleterious and beneficial mutations in the degree of adaptation of vesicular stomatitis virus (VSV), a nonsegmented virus of negative strand. We also have studied how effective population size, initial genetic variability in populations, and environmental heterogeneity shapes the impact of mutations in the evolution of vesicular stomatitis virus. Finally, in an integrative attempt, we discuss pros and cons of the quasispecies theory compared with classic population genetics models for haploid organisms to explain the evolution of RNA viruses.

Gene genealogies and population variation in plants

Early in the development of plant evolutionary biology, genetic drift, fluctuations in population size, and isolation were identified as critical processes that affect the course of evolution in plant species. Attempts to assess these processes in natural populations became possible only with the development of neutral genetic markers in the 1960s. More recently, the application of historically ordered neutral molecular variation (within the conceptual framework of coalescent theory) has allowed a reevaluation of these microevolutionary processes. Gene genealogies trace the evolutionary relationships among haplotypes (alleles) with populations. Processes such as selection, fluctuation in population size, and population substructuring affect the geographical and genealogical relationships among these alleles. Therefore, examination of these genealogical data can provide insights into the evolutionary history of a species. For example, studies of Arabidopsis thaliana have suggested that this species underwent rapid expansion, with populations showing little genetic differentiation. The new discipline of phylogeography examines the distribution of allele genealogies in an explicit geographical context. Phylogeographic studies of plants have documented the recolonization of European tree species from refugia subsequent to Pleistocene glaciation, and such studies have been instructive in understanding the origin and domestication of the crop cassava. Currently, several technical limitations hinder the widespread application of a genealogical approach to plant evolutionary studies. However, as these technical issues are solved, a genealogical approach holds great promise for understanding these previously elusive processes in plant evolution.