Ecological factors rather than temporal factors dominate the evolution of vesicular stomatitis virus

Vesicular stomatitis New Jersey virus (VSV-NJ) is a rhabdovirus that causes economically important disease in cattle and other domestic animals in endemic areas from southeastern United States to northern South America. Its negatively stranded RNA genome is capable of undergoing rapid evolution, which allows phylogenetic analysis and molecular epidemiology studies to be performed. Previous epidemiological studies in Costa Rica showed the existence of at least two distinct ecological zones of high VSV-NJ activity, one located in the highlands (premontane tropical moist forest) and the other in the lowlands (tropical dry forest). We wanted to test the hypothesis that the viruses circulating in these ecological zones were genetically distinct. For this purpose, we sequenced the hypervariable region of the phosphoprotein gene for 50 VSV-NJ isolates from these areas. Phylogenetic analysis showed that viruses from each ecological zone had distinct genotypes. These genotypes were maintained in each area for periods of up to 8 years. This evolutionary pattern of VSV-NJ suggests an adaptation to ecological factors that could exert selective pressure on the virus. As previous data indicated an absence of virus adaptation to factors related to the bovine host (including immunological pressure), it appears that VSV genetic divergence represents positive selection to adapt to specific vectors and/or reservoirs at each ecological zone.

Australian lungfish neurohypophysial hormone genes encode vasotocin and [Phe2]mesotocin precursors homologous to tetrapod-type precursors

In view of the well-established role of neurohypophysial hormones in osmoregulation of terrestrial vertebrates, lungfishes are a key group for study of the molecular and functional evolution of the hypothalamo-neurohypophysial system. Here we report on the primary structure of the precursors encoding vasotocin (VT) and [Phe2]mesotocin ([Phe2]MT) of the Australian lungfish, Neoceratodus forsteri. Genomic sequence analysis and Northern blot analysis confirmed that [Phe2]MT is a native oxytocin family peptide in the Australian lungfish, although it has been reported that the lungfish neurohypophysis contains MT. The VT precursor consists of a signal peptide, VT, that is connected to a neurophysin by a Gly-Lys-Arg sequence, and a copeptin moiety that includes a Leu-rich core segment and a glycosylation site. In contrast, the [Phe2]MT precursor does not contain a copeptin moiety. These structural features of the lungfish precursors are consistent with those in tetrapods, but different from those in teleosts where both VT and isotocin precursors contain a copeptin-like moiety without a glycosylation site at the carboxyl terminals of their neurophysins. Comparison of the exon/intron organization also supports homology of the lungfish [Phe2]MT gene with tetrapod oxytocin/MT genes, rather than with teleost isotocin genes. Moreover, molecular phylogenetic analysis shows that neurohypophysial hormone genes of the lungfish are closely related to those of the toad. The present results along with previous morphological findings indicate that the hypothalamo-neurohypophysial system of the lungfish has evolved along the tetrapod lineage, whereas the teleosts form a separate lineage, both within the class Osteichthyes.

A functional homolog of a yeast tRNA splicing enzyme is conserved in higher eukaryotes and in Escherichia coli

tRNA splicing in the yeast Saccharomyces cerevisiae requires an endonuclease to excise the intron, tRNA ligase to join the tRNA half-molecules, and 2′-phosphotransferase to transfer the splice junction 2′-phosphate from ligated tRNA to NAD, producing ADP ribose 1′′–2′′ cyclic phosphate (Appr>p). We show here that functional 2′-phosphotransferases are found throughout eukaryotes, occurring in two widely divergent yeasts (Candida albicans and Schizosaccharomyces pombe), a plant (Arabidopsis thaliana), and mammals (Mus musculus); this finding is consistent with a role for the enzyme, acting in concert with ligase, to splice tRNA or other RNA molecules. Surprisingly, functional 2′-phosphotransferase is found also in the bacterium Escherichia coli, which does not have any known introns of this class, and does not appear to have a ligase that generates junctions with a 2′-phosphate. Analysis of the database shows that likely members of the 2′-phosphotransferase family are found also in one other bacterium (Pseudomonas aeruginosa) and two archaeal species (Archaeoglobus fulgidus and Pyrococcus horikoshii). Phylogenetic analysis reveals no evidence for recent horizontal transfer of the 2′-phosphotransferase into Eubacteria, suggesting that the 2′-phosphotransferase has been present there since close to the time that the three kingdoms diverged. Although 2′-phosphotransferase is not present in all Eubacteria, and a gene disruption experiment demonstrates that the protein is not essential in E. coli, the continued presence of 2′-phosphotransferase in Eubacteria over large evolutionary times argues for an important role for the protein.

Chloroplast small heat shock proteins: Evidence for atypical evolution of an organelle-localized protein

Knowledge of the origin and evolution of gene families is critical to our understanding of the evolution of protein function. To gain a detailed understanding of the evolution of the small heat shock proteins (sHSPs) in plants, we have examined the evolutionary history of the chloroplast (CP)-localized sHSPs. Previously, these nuclear-encoded CP proteins had been identified only from angiosperms. This study reveals the presence of the CP sHSPs in a moss, Funaria hygrometrica. Two clones for CP sHSPs were isolated from a F. hygrometrica heat shock cDNA library that represent two distinct CP sHSP genes. Our analysis of the CP sHSPs reveals unexpected evolutionary relationships and patterns of sequence conservation. Phylogenetic analysis of the CP sHSPs with other plant CP sHSPs and eukaryotic, archaeal, and bacterial sHSPs shows that the CP sHSPs are not closely related to the cyanobacterial sHSPs. Thus, they most likely evolved via gene duplication from a nuclear-encoded cytosolic sHSP and not via gene transfer from the CP endosymbiont. Previous sequence analysis had shown that all angiosperm CP sHSPs possess a methionine-rich region in the N-terminal domain. The primary sequence of this region is not highly conserved in the F. hygrometrica CP sHSPs. This lack of sequence conservation indicates that sometime in land plant evolution, after the divergence of mosses from the common ancestor of angiosperms but before the monocot–dicot divergence, there was a change in the selective constraints acting on the CP sHSPs.

Duplicated genes evolve independently after polyploid formation in cotton

Of the many processes that generate gene duplications, polyploidy is unique in that entire genomes are duplicated. This process has been important in the evolution of many eukaryotic groups, and it occurs with high frequency in plants. Recent evidence suggests that polyploidization may be accompanied by rapid genomic changes, but the evolutionary fate of discrete loci recently doubled by polyploidy (homoeologues) has not been studied. Here we use locus-specific isolation techniques with comparative mapping to characterize the evolution of homoeologous loci in allopolyploid cotton (Gossypium hirsutum) and in species representing its diploid progenitors. We isolated and sequenced 16 loci from both genomes of the allopolyploid, from both progenitor diploid genomes and appropriate outgroups. Phylogenetic analysis of the resulting 73.5 kb of sequence data demonstrated that for all 16 loci (14.7 kb/genome), the topology expected from organismal history was recovered. In contrast to observations involving repetitive DNAs in cotton, there was no evidence of interaction among duplicated genes in the allopolyploid. Polyploidy was not accompanied by an obvious increase in mutations indicative of pseudogene formation. Additionally, differences in rates of divergence among homoeologues in polyploids and orthologues in diploids were indistinguishable across loci, with significant rate deviation restricted to two putative pseudogenes. Our results indicate that most duplicated genes in allopolyploid cotton evolve independently of each other and at the same rate as those of their diploid progenitors. These indications of genic stasis accompanying polyploidization provide a sharp contrast to recent examples of rapid genomic evolution in allopolyploids.

Tracking molecular evolution of photosynthesis by characterization of a major photosynthesis gene cluster from Heliobacillus mobilis

A DNA sequence has been obtained for a 35.6-kb genomic segment from Heliobacillus mobilis that contains a major cluster of photosynthesis genes. A total of 30 ORFs were identified, 20 of which encode enzymes for bacteriochlorophyll and carotenoid biosynthesis, reaction-center (RC) apoprotein, and cytochromes for cyclic electron transport. Donor side electron-transfer components to the RC include a putative RC-associated cytochrome c553 and a unique four-large-subunit cytochrome bc complex consisting of Rieske Fe-S protein (encoded by petC), cytochrome b6 (petB), subunit IV (petD), and a diheme cytochrome c (petX). Phylogenetic analysis of various photosynthesis gene products indicates a consistent grouping of oxygenic lineages that are distinct and descendent from anoxygenic lineages. In addition, H. mobilis was placed as the closest relative to cyanobacteria, which form a monophyletic origin to chloroplast-based photosynthetic lineages. The consensus of the photosynthesis gene trees also indicates that purple bacteria are the earliest emerging photosynthetic lineage. Our analysis also indicates that an ancient gene-duplication event giving rise to the paralogous bchI and bchD genes predates the divergence of all photosynthetic groups. In addition, our analysis of gene duplication of the photosystem I and photosystem II core polypeptides supports a “heterologous fusion model” for the origin and evolution of oxygenic photosynthesis.

Evolution of paired domains: Isolation and sequencing of jellyfish and hydra Pax genes related to Pax-5 and Pax-6

Pax proteins are a family of transcription factors with a highly conserved paired domain; many members also contain a paired-type homeodomain and/or an octapeptide. Nine mammalian Pax genes are known and classified into four subgroups: Pax-1/9, Pax-2/5/8, Pax-3/7, and Pax-4/6. Most of these genes are involved in nervous system development. In particular, Pax-6 is a key regulator that controls eye development in vertebrates and Drosophila. Although the Pax-4/6 subgroup seems to be more closely related to Pax-2/5/8 than to Pax-3/7 or Pax-1/9, its evolutionary origin is unknown. We therefore searched for a Pax-6 homolog and related genes in Cnidaria, which is the lowest phylum of animals that possess a nervous system and eyes. A sea nettle (a jellyfish) genomic library was constructed and two pax genes (Pax-A and -B) were isolated and partially sequenced. Surprisingly, unlike most known Pax genes, the paired box in these two genes contains no intron. In addition, the complete cDNA sequences of hydra Pax-A and -B were obtained. Hydra Pax-B contains both the homeodomain and the octapeptide, whereas hydra Pax-A contains neither. DNA binding assays showed that sea nettle Pax-A and -B and hydra Pax-A paired domains bound to a Pax-5/6 site and a Pax-5 site, although hydra Pax-B paired domain bound neither. An alignment of all available paired domain sequences revealed two highly conserved regions, which cover the DNA binding contact positions. Phylogenetic analysis showed that Pax-A and especially Pax-B were more closely related to Pax-2/5/8 and Pax-4/6 than to Pax-1/9 or Pax-3/7 and that the Pax genes can be classified into two supergroups: Pax-A/Pax-B/Pax-2/5/8/4/6 and Pax-1/9/3/7. From this analysis and the gene structure, we propose that modern Pax-4/6 and Pax-2/5/8 genes evolved from an ancestral gene similar to cnidarian Pax-B, having both the homeodomain and the octapeptide.

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.

Identification of a Second Myosin-II in Schizosaccharomyces pombe:: Myp2p Is Conditionally Required for Cytokinesis

As in many eukaryotic cells, fission yeast cytokinesis depends on the assembly of an actin ring. We cloned myp2+, a myosin-II in Schizosaccharomyces pombe, conditionally required for cytokinesis. myp2+, the second myosin-II identified in S. pombe, does not completely overlap in function with myo2+. The catalytic domain of Myp2p is highly homologous to known myosin-IIs, and phylogenetic analysis places Myp2p in the myosin-II family. The Myp2p sequence contains well-conserved ATP- and actin-binding motifs, as well as two IQ motifs. However, the tail sequence is unusual, since it is predicted to form two long coiled-coils separated by a stretch of sequence containing 19 prolines. Disruption of myp2+ is not lethal but under nutrient limiting conditions cells lacking myp2+ function are multiseptated, elongated, and branched, indicative of a defect in cytokinesis. The presence of salt enhances these morphological defects. Additionally, Δmyp2 cells are cold sensitive in high salt, failing to form colonies at 17°C. Thus, myp2+ is required under conditions of stress, possibly linking extracellular growth conditions to efficient cytokinesis and cell growth. GFP-Myp2p localizes to a ring in the middle of late mitotic cells, consistent with a role in cytokinesis. Additionally, we constructed double mutants of Δmyp2 with temperature-sensitive mutant strains defective in cytokinesis. We observed synthetic lethal interactions between Δmyp2 and three alleles of cdc11ts, as well as more modest synthetic interactions with cdc14ts and cdc16ts, implicating myp2+ function for efficient cytokinesis under normal conditions.

Multiple independent transpositions of mitochondrial DNA control region sequences to the nucleus

Transpositions of mtDNA sequences to the nuclear genome have been documented in a wide variety of individual taxa, but little is known about their taxonomic frequency or patterns of variation. We provide evidence of nuclear sequences homologous to the mtDNA control region in seven species of diving ducks (tribe Aythyini). Phylogenetic analysis places each nuclear sequence as a close relative of the mtDNA haplotypes of the specie(s) in which it occurs, indicating that they derive from six independent transposition events, all occurring within the last ≈1.5 million years. Relative-rate tests and comparison of intraspecific variation in nuclear and mtDNA sequences confirm the expectation of a greatly reduced rate of evolution in the nuclear copies. By representing mtDNA haplotypes from ancestral populations, nuclear insertions may be valuable in some phylogenetic analyses, but they also confound the accurate determination of mtDNA sequences. In particular, our data suggest that the presumably nonfunctional but more slowly evolving nuclear sequences often will not be identifiable by changes incompatible with function and may be preferentially amplified by PCR primers based on mtDNA sequences from related taxa.

Cryptocyanin, a crustacean molting protein: Evolutionary link with arthropod hemocyanins and insect hexamerins

Cryptocyanin, a copper-free hexameric protein in crab (Cancer magister) hemolymph, has been characterized and the amino acid sequence has been deduced from its cDNA. It is markedly similar in sequence, size, and structure to hemocyanin, the copper-containing oxygen-transport protein found in many arthropods. Cryptocyanin does not bind oxygen, however, and lacks three of the six highly conserved copper-binding histidine residues of hemocyanin. Cryptocyanin has no phenoloxidase activity, although a phenoloxidase is present in the hemolymph. The concentration of cryptocyanin in the hemolymph is closely coordinated with the molt cycle and reaches levels higher than hemocyanin during premolt. Cryptocyanin resembles insect hexamerins in the lack of copper, molt cycle patterns of biosynthesis, and potential contributions to the new exoskeleton. Phylogenetic analysis of sequence similarities between cryptocyanin and other members of the hemocyanin gene family shows that cryptocyanin is closely associated with crustacean hemocyanins and suggests that cryptocyanin arose as a result of a hemocyanin gene duplication. The presence of both hemocyanin and cryptocyanin in one animal provides an example of how insect hexamerins might have evolved from hemocyanin. Our results suggest that multiple members of the hemocyanin gene family—hemocyanin, cryptocyanin, phenoloxidase, and hexamerins—may participate in two vital functions of molting animals, oxygen binding and molting. Cryptocyanin may provide important molecular data to further investigate evolutionary relationships among all molting animals.

Molecular studies suggest that cartilaginous fishes have a terminal position in the piscine tree

The Chondrichthyes (cartilaginous fishes) are commonly accepted as being sister group to the other extant Gnathostomata (jawed vertebrates). To clarify gnathostome relationships and to aid in resolving and dating the major piscine divergences, we have sequenced the complete mtDNA of the starry skate and have included it in phylogenetic analysis along with three squalomorph chondrichthyans—the common dogfish, the spiny dogfish, and the star spotted dogfish—and a number of bony fishes and amniotes. The direction of evolution within the gnathostome tree was established by rooting it with the most closely related non-gnathostome outgroup, the sea lamprey, as well as with some more distantly related taxa. The analyses placed the chondrichthyans in a terminal position in the piscine tree. These findings, which also suggest that the origin of the amniote lineage is older than the age of the oldest extant bony fishes (the lungfishes), challenge the evolutionary direction of several morphological characters that have been used in reconstructing gnathostome relationships. Applying as a calibration point the age of the oldest lungfish fossils, 400 million years, the molecular estimate placed the squalomorph/batomorph divergence at ≈190 million years before present. This dating is consistent with the occurrence of the earliest batomorph (skates and rays) fossils in the paleontological record. The split between gnathostome fishes and the amniote lineage was dated at ≈420 million years before present.

Identification of a third distinct estrogen receptor and reclassification of estrogen receptors in teleosts

This paper describes three distinct estrogen receptor (ER) subtypes: ERα, ERβ, and a unique type, ERγ, cloned from a teleost fish, the Atlantic croaker Micropogonias undulatus; the first identification of a third type of classical ER in vertebrate species. Phylogenetic analysis shows that ERγ arose through gene duplication from ERβ early in the teleost lineage and indicates that ERγ is present in other teleosts, although it has not been recognized as such. The Atlantic croaker ERγ shows amino acid differences in regions important for ligand binding and receptor activation that are conserved in all other ERγs. The three ER subtypes are genetically distinct and have different distribution patterns in Atlantic croaker tissues. In addition, ERβ and ERγ fusion proteins can each bind estradiol-17β with high affinity. The presence of three functional ERs in one species expands the role of ER multiplicity in estrogen signaling systems and provides a unique opportunity to investigate the dynamics and mechanisms of ER evolution.

Evolutionary specialization of the nuclear targeting apparatus

The α- and β-karyopherins (Kaps), also called importins, mediate the nuclear transport of proteins. All α-Kaps contain a central domain composed of eight approximately 40 amino acid, tandemly arranged, armadillo-like (Arm) repeats. The number and order of these repeats have not changed since the common origin of fungi, plants, and mammals. Phylogenetic analysis suggests that the various α-Kaps fall into two groups, α1 and α2. Whereas animals encode both types, the yeast genome encodes only an α1-Kap. The β-Kaps are characterized by 14–15 tandemly arranged HEAT motifs. We show that the Arm repeats of α-Kaps and the HEAT motifs of β-Kaps are similar, suggesting that the α-Kaps and β-Kaps (and for that matter, all Arm and HEAT repeat-containing proteins) are members of the same protein superfamily. Phylogenetic analysis indicates that there are at least three major groups of β-Kaps, consistent with their proposed cargo specificities. We present a model in which an α-independent β-Kap progenitor gave rise to the α-dependent β-Kaps and the α-Kaps.

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.