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.

The human sex-reversing ATRX gene has a homologue on the marsupial Y chromosome, ATRY: Implications for the evolution of mammalian sex determination

Mutations in the ATRX gene on the human X chromosome cause X-linked α-thalassemia and mental retardation. XY patients with deletions or mutations in this gene display varying degrees of sex reversal, implicating ATRX in the development of the human testis. To explore further the role of ATRX in mammalian sex differentiation, the homologous gene was cloned and characterized in a marsupial. Surprisingly, active homologues of ATRX were detected on the marsupial Y as well as the X chromosome. The Y-borne copy (ATRY) displays testis-specific expression. This, as well as the sex reversal of ATRX patients, suggests that ATRY is involved in testis development in marsupials and may represent an ancestral testis-determining mechanism that predated the evolution of SRY as the primary mammalian male sex-determining gene. There is no evidence for a Y-borne ATRX homologue in mouse or human, implying that this gene has been lost in eutherians and its role supplanted by the evolution of SRY from SOX3 as the dominant determiner of male differentiation.

Directed evolution of a (βα)8-barrel enzyme to catalyze related reactions in two different metabolic pathways

Enzymes participating in different metabolic pathways often have similar catalytic mechanisms and structures, suggesting their evolution from a common ancestral precursor enzyme. We sought to create a precursor-like enzyme for N′-[(5′-phosphoribosyl)formimino]-5-aminoimidazole-4-carboxamide ribonucleotide (ProFAR) isomerase (HisA; EC 5.3.1.16) and phosphoribosylanthranilate (PRA) isomerase (TrpF; EC 5.3.1.24), which catalyze similar reactions in the biosynthesis of the amino acids histidine and tryptophan and have a similar (βα)8-barrel structure. Using random mutagenesis and selection, we generated several HisA variants that catalyze the TrpF reaction both in vivo and in vitro, and one of these variants retained significant HisA activity. A more detailed analysis revealed that a single amino acid exchange could establish TrpF activity on the HisA scaffold. These findings suggest that HisA and TrpF may have evolved from an ancestral enzyme of broader substrate specificity and underscore that (βα)8-barrel enzymes are very suitable for the design of new catalytic activities.

Cloning and functional expression of a cDNA encoding a pheromone gland-specific acyl-CoA Δ11-desaturase of the cabbage looper moth, Trichoplusia ni

Desaturation of coenzyme-A esters of saturated fatty acids is a common feature of sex pheromone biosynthetic pathways in the Lepidoptera. The enzymes that catalyze this step share several biochemical properties with the ubiquitous acyl-CoA Δ9-desaturases of animals and fungi, suggesting a common ancestral origin. Unlike metabolic acyl-CoA Δ9-desaturases, pheromone desaturases have evolved unusual regio- and stereoselective activities that contribute to the remarkable diversity of chemical structures used as pheromones in this large taxonomic group. In this report, we describe the isolation of a cDNA encoding a pheromone gland desaturase from the cabbage looper moth, Trichoplusia ni, a species in which all unsaturated pheromone products are produced via a Δ11Z-desaturation mechanism. The largest ORF of the ≈1,250-bp cDNA encodes a 349-aa apoprotein (PDesat-Tn Δ11Z) with a predicted molecular mass of 40,240 Da. Its hydrophobicity profile is similar overall to those of rat and yeast Δ9-desaturases, suggesting conserved transmembrane topology. A 182-aa core domain delimited by conserved histidine-rich motifs implicated in iron-binding and catalysis has 72 and 58% similarity (including conservative substitutions) to acyl-CoA Δ9Z-desaturases of rat and yeast, respectively. Northern blot analysis revealed an ≈1,250-nt PDesat-Tn Δ11Z mRNA that is consistent with the spatial and temporal distribution of Δ11-desaturase enzyme activity. Genetic transformation of a desaturase-deficient strain of the yeast Saccharomyces cerevisiae with an expression plasmid encoding PDesat-Tn Δ11Z resulted in complementation of the strain’s fatty acid auxotrophy and the production of Δ11Z-unsaturated fatty acids.

Using rare mutations to estimate population divergence times: A maximum likelihood approach

In this paper we propose a method to estimate by maximum likelihood the divergence time between two populations, specifically designed for the analysis of nonrecurrent rare mutations. Given the rapidly growing amount of data, rare disease mutations affecting humans seem the most suitable candidates for this method. The estimator RD, and its conditional version RDc, were derived, assuming that the population dynamics of rare alleles can be described by using a birth–death process approximation and that each mutation arose before the split of a common ancestral population into the two diverging populations. The RD estimator seems more suitable for large sample sizes and few alleles, whose age can be approximated, whereas the RDc estimator appears preferable when this is not the case. When applied to three cystic fibrosis mutations, the estimator RD could not exclude a very recent time of divergence among three Mediterranean populations. On the other hand, the divergence time between these populations and the Danish population was estimated to be, on the average, 4,500 or 15,000 years, assuming or not a selective advantage for cystic fibrosis carriers, respectively. Confidence intervals are large, however, and can probably be reduced only by analyzing more alleles or loci.

Independent deletions of a pathogen-resistance gene in Brassica and Arabidopsis

Plant disease resistance (R) genes confer race-specific resistance to pathogens and are genetically defined on the basis of intra-specific functional polymorphism. Little is known about the evolutionary mechanisms that generate this polymorphism. Most R loci examined to date contain alternate alleles and/or linked homologs even in disease-susceptible plant genotypes. In contrast, the resistance to Pseudomonas syringae pathovar maculicola (RPM1) bacterial resistance gene is completely absent (rpm1-null) in 5/5 Arabidopsis thaliana accessions that lack RPM1 function. The rpm1-null locus contains a 98-bp segment of unknown origin in place of the RPM1 gene. We undertook comparative mapping of RPM1 and flanking genes in Brassica napus to determine the ancestral state of the RPM1 locus. We cloned two B. napus RPM1 homologs encoding hypothetical proteins with ≈81% amino acid identity to Arabidopsis RPM1. Collinearity of genes flanking RPM1 is conserved between B. napus and Arabidopsis. Surprisingly, we found four additional B. napus loci in which the flanking marker synteny is maintained but RPM1 is absent. These B. napus rpm1-null loci have no detectable nucleotide similarity to the Arabidopsis rpm1-null allele. We conclude that RPM1 evolved before the divergence of the Brassicaceae and has been deleted independently in the Brassica and Arabidopsis lineages. These results suggest that functional polymorphism at R gene loci can arise from gene deletions.

Regulation of number and size of digits by posterior Hox genes: A dose-dependent mechanism with potential evolutionary implications

The proper development of digits, in tetrapods, requires the activity of several genes of the HoxA and HoxD homeobox gene complexes. By using a variety of loss-of-function alleles involving the five Hox genes that have been described to affect digit patterning, we report here that the group 11, 12, and 13 genes control both the size and number of murine digits in a dose-dependent fashion, rather than through a Hox code involving differential qualitative functions. A similar dose–response is observed in the morphogenesis of the penian bone, the baculum, which further suggests that digits and external genitalia share this genetic control mechanism. A progressive reduction in the dose of Hox gene products led first to ectrodactyly, then to olygodactyly and adactyly. Interestingly, this transition between the pentadactyl to the adactyl formula went through a step of polydactyly. We propose that in the distal appendage of polydactylous short-digited ancestral tetrapods, such as Acanthostega, the HoxA complex was predominantly active. Subsequent recruitment of the HoxD complex contributed to both reductions in digit number and increase in digit length. Thus, transition through a polydactylous limb before reaching and stabilizing the pentadactyl pattern may have relied, at least in part, on asynchronous and independent changes in the regulation of HoxA and HoxD gene complexes.

A shark antibody heavy chain encoded by a nonsomatically rearranged VDJ is preferentially expressed in early development and is convergent with mammalian IgG

In most vertebrate embryos and neonates studied to date unique antigen receptors (antibodies and T cell receptors) are expressed that possess a limited immune repertoire. We have isolated a subclass of IgM, IgM1gj, from the nurse shark Ginglymostoma cirratum that is preferentially expressed in neonates. The variable (V) region gene encoding the heavy (H) chain underwent V-D-J rearrangement in germ cells (“germline-joined”). Such H chain V genes were discovered over 10 years ago in sharks but until now were not shown to be expressed at appreciable levels; we find expression of H1gj in primary and secondary lymphoid tissues early in life, but in adults only in primary lymphoid tissue, which is identified in this work as the epigonal organ. H1gj chain associates covalently with light (L) chains and is most similar in sequence to IgM H chains, but like mammalian IgG has three rather than the four IgM constant domains; deletion of the ancestral IgM C2 domain thus defines both IgG and IgM1gj. Because sharks are the members of the oldest vertebrate class known to possess antibodies, unique or specialized antibodies expressed early in ontogeny in sharks and other vertebrates were likely present at the inception of the adaptive immune system.

Two different neurodegenerative diseases caused by proteins with similar structures

The downstream prion-like protein (doppel, or Dpl) is a paralog of the cellular prion protein, PrPC. The two proteins have ≈25% sequence identity, but seem to have distinct physiologic roles. Unlike PrPC, Dpl does not support prion replication; instead, overexpression of Dpl in the brain seems to cause a completely different neurodegenerative disease. We report the solution structure of a fragment of recombinant mouse Dpl (residues 26–157) containing a globular domain with three helices and a small amount of β-structure. Overall, the topology of Dpl is very similar to that of PrPC. Significant differences include a marked kink in one of the helices in Dpl, and a different orientation of the two short β-strands. Although the two proteins most likely arose through duplication of a single ancestral gene, the relationship is now so distant that only the structures retain similarity; the functions have diversified along with the sequence.

Chlorophyll and carotenoid binding in a simple red algal light-harvesting complex crosses phylogenetic lines

The membrane proteins of peripheral light-harvesting complexes (LHCs) bind chlorophylls and carotenoids and transfer energy to the reaction centers for photosynthesis. LHCs of chlorophytes, chromophytes, dinophytes, and rhodophytes are similar in that they have three transmembrane regions and several highly conserved Chl-binding residues. All LHCs bind Chl a, but in specific taxa certain characteristic pigments accompany Chl a: Chl b and lutein in chlorophytes, Chl c and fucoxanthin in chromophytes, Chl c and peridinin in dinophytes, and zeaxanthin in rhodophytes. The specificity of pigment binding was examined by in vitro reconstitution of various pigments with a simple light-harvesting protein (LHCaR1), from a red alga (Porphyridium cruentum), that normally has eight Chl a and four zeaxanthin molecules. The pigments typical of a chlorophyte (Spinacea oleracea), a chromophyte (Thallasiosira fluviatilis), and a dinophyte (Prorocentrum micans) were found to functionally bind to this protein as evidenced by their participation in energy transfer to Chl a, the terminal pigment. This is a demonstration of a functional relatedness of rhodophyte and higher plant LHCs. The results suggest that eight Chl-binding sites per polypeptide are an ancestral trait, and that the flexibility to bind various Chl and carotenoid pigments may have been retained throughout the evolution of LHCs.

KH domain: one motif, two folds

The K homology (KH) module is a widespread RNA-binding motif that has been detected by sequence similarity searches in such proteins as heterogeneous nuclear ribonucleoprotein K (hnRNP K) and ribosomal protein S3. Analysis of spatial structures of KH domains in hnRNP K and S3 reveals that they are topologically dissimilar and thus belong to different protein folds. Thus KH motif proteins provide a rare example of protein domains that share significant sequence similarity in the motif regions but possess globally distinct structures. The two distinct topologies might have arisen from an ancestral KH motif protein by N- and C-terminal extensions, or one of the existing topologies may have evolved from the other by extension, displacement and deletion. C-terminal extension (deletion) requires β-sheet rearrangement through the insertion (removal) of a β-strand in a manner similar to that observed in serine protease inhibitors serpins. Current analysis offers a new look on how proteins can change fold in the course of evolution.

Feeding specialization and host-derived chemical defense in Chrysomeline leaf beetles did not lead to an evolutionary dead end

Combination of molecular phylogenetic analyses of Chrysomelina beetles and chemical data of their defensive secretions indicate that two lineages independently developed, from an ancestral autogenous metabolism, an energetically efficient strategy that made the insect tightly dependent on the chemistry of the host plant. However, a lineage (the interrupta group) escaped this subordination through the development of a yet more derived mixed metabolism potentially compatible with a large number of new host-plant associations. Hence, these analyses on leaf beetles document a mechanism that can explain why high levels of specialization do not necessarily lead to “evolutionary dead ends.”

The nucleotide changes governing cuticular hydrocarbon variation and their evolution in Drosophila melanogaster

The cuticular hydrocarbon (CH) pheromones in Drosophila melanogaster exhibit strong geographic variation. African and Caribbean populations have a high ratio of 5,9 heptacosadiene/7,11 heptacosadiene (the “High” CH type), whereas populations from all other areas have a low ratio (“Low” CH type). Based on previous genetic mapping, DNA markers were developed that localized the genetic basis of this CH polymorphism to within a 13-kb region. We then carried out a hierarchical search for diagnostic nucleotide sites starting with four lines, and increasing to 24 and 43 lines from a worldwide collection. Within the 13-kb region, only one variable site shows a complete concordance with the CH phenotype. This is a 16-bp deletion in the 5′ region of a desaturase gene (desat2) that was recently suggested to be responsible for the CH polymorphism on the basis of its expression [Dallerac, R., Labeur, C., Jallon, J.-M., Knipple, D. C., Roelofs, W. L. & Wicker-Thomas, C. (2000) Proc. Natl. Acad. Sci. 97, 9449–9454]. The cosmopolitan Low type is derived from the ancestral High type, and DNA sequence variations suggest that the former spread worldwide with the aid of positive selection. Whether this CH variation could be a component of the sexual isolation between Zimbabwe and other cosmopolitan populations remains an interesting and unresolved question.

Molecular Analysis of (R)-(+)-Mandelonitrile Lyase Microheterogeneity in Black Cherry1

The flavoprotein (R)-(+)-mandelonitrile lyase (MDL; EC 4.1.2.10), which plays a key role in cyanogenesis in rosaceous stone fruits, occurs in black cherry (Prunus serotina Ehrh.) homogenates as several closely related isoforms. Biochemical and molecular biological methods were used to investigate MDL microheterogeneity and function in this species. Three novel MDL cDNAs of high sequence identity (designated MDL2, MDL4, and MDL5) were isolated. Like MDL1 and MDL3 cDNAs (Z. Hu, J.E. Poulton [1997] Plant Physiol 115: 1359–1369), they had open reading frames that predicted a flavin adenine dinucleotide-binding site, multiple N-glycosylation sites, and an N-terminal signal sequence. The N terminus of an MDL isoform purified from seedlings matched the derived amino acid sequence of the MDL4 cDNA. Genomic sequences corresponding to the MDL1, MDL2, and MDL4 cDNAs were obtained by polymerase chain reaction amplification of genomic DNA. Like the previously reported mdl3 gene, these genes are interrupted at identical positions by three short, conserved introns. Given their overall similarity, we conclude that the genes mdl1, mdl2, mdl3, mdl4, and mdl5 are derived from a common ancestral gene and constitute members of a gene family. Genomic Southern-blot analysis showed that this family has approximately eight members. Northern-blot analysis using gene-specific probes revealed differential expression of the genes mdl1, mdl2, mdl3, mdl4, and mdl5.

Global survey of genetic variation in CCR5, RANTES, and MIP-1α: Impact on the epidemiology of the HIV-1 pandemic

Expression of CC chemokine receptor 5 (CCR5), the major coreceptor for HIV-1 cell entry, and its ligands (e.g., RANTES and MIP-1α) is widely regarded as central to the pathogenesis of HIV-1 infection. By surveying nearly 3,000 HIV+ and HIV− individuals from worldwide populations for polymorphisms in the genes encoding RANTES, MIP-1α, and CCR5, we show that the evolutionary histories of human populations have had a significant impact on the distribution of variation in these genes, and that this may be responsible, in part, for the heterogeneous nature of the epidemiology of the HIV-1 pandemic. The varied distribution of RANTES haplotypes (AC, GC, and AG) associated with population-specific HIV-1 transmission- and disease-modifying effects is a striking example. Homozygosity for the AC haplotype was associated with an increased risk of acquiring HIV-1 as well as accelerated disease progression in European Americans, but not in African Americans. Yet, the prevalence of the ancestral AC haplotype is high in individuals of African origin, but substantially lower in non-Africans. In a Japanese cohort, AG-containing RANTES haplotype pairs were associated with a delay in disease progression; however, we now show that their contribution to HIV-1 pathogenesis and epidemiology in other parts of the world is negligible because the AG haplotype is infrequent in non-Far East Asians. Thus, the varied distribution of RANTES, MIP-1α, and CCR5 haplotype pairs and their population-specific phenotypic effects on HIV-1 susceptibility and disease progression results in a complex pattern of biological determinants of HIV-1 epidemiology. These findings have important implications for the design, assessment, and implementation of effective HIV-1 intervention and prevention strategies.