Ferrets as a transmission model for influenza: sequence changes in HA1 of type A (H3N2) virus.

Ferrets were used as an animal model to study whether controlled transmission of type A influenza is similar to human transmission when sequence changes in HA1 are used as the outcome. Ferrets were infected initially with A/Sydney/5/97 (H3N2) or A/LA/1/87 (H3N2) intranasally, and transmission chains were established by housing infected ferrets with noninfected ferrets with no influenza antibody titer against the infecting virus. Ferrets infected with A/Sydney were seronegative for A/Sydney and A/LA; ferrets infected with A/LA were seronegative for A/LA but had hemagglutination inhibition titers against A/Sydney. Titers of naturally transmitted influenza were higher than those after direct intranasal infection, but lymphocyte counts from nasal washes diminished with transmission. Ferrets infected with A/LA had 2 amino acid differences in HA1 after transmission through 5 ferret cohorts, but those infected with A/Sydney had none. The results show the value of the ferret model. A/LA resembled the transmission of influenza in humans when under antibody pressure.

Evolution and characterization of tetraonine endogenous retrovirus: a new virus related to avian sarcoma and leukosis viruses.

In a previous study, we found avian sarcoma and leukosis virus (ASLV) gag genes in 19 species of birds in the order Galliformes including all grouse and ptarmigan (Tetraoninae) surveyed. Our data suggested that retroviruses had been transmitted horizontally among some host species. To further investigate these elements, we sequenced a replication-defective retrovirus, here named tetraonine endogenous retrovirus (TERV), from Bonasa umbellus (ruffed grouse). This is the first report of a complete, replication-defective ASLV provirus sequence from any bird other than the domestic chicken. We found a replication-defective proviral sequence consisting of putative Gag and Env proteins flanked by long terminal repeats. Reverse transcription-PCR analysis showed that retroviral gag sequences closely related to TERV are transcribed, supporting the hypothesis that TERV is an active endogenous retrovirus. Phylogenetic analyses suggest that TERV may have arisen via recombination between different retroviral lineages infecting birds. Southern blotting using gag probes showed that TERV occurs in tetraonines but not in chickens or ducks, suggesting that integration occurred after the earliest phasianid divergences but prior to the radiation of tetraonine birds.

Rooting a phylogeny with homologous genes on opposite sex chromosomes (gametologs): a case study using avian CHD.

We describe a previously unrecognized form of gene homology using the term "gametology," which we define as homology arising through lack of recombination and subsequent differentiation of sex chromosomes. We demonstrate use of gametologous genes to root each other in phylogenetic analyses of sex-specific avian Chromo-helicase-DNA binding gene (CHD) sequences. Phylogenetic analyses of a set of neognath bird sequences yield monophyletic groups for CHD-W and CHD-Z gametologs, as well as congruent relationships between these two clades and between them and current views of avian taxonomy. Phylogenetic analyses including paleognath bird CHD sequences and rooting with crocodilian CHD sequences, suggest an early divergence for paleognath CHD within the avian CHD clade. Based on our CHD analyses calibrated with avian fossil dates, we estimate the divergence between CHD-W and CHD-Z at 123 MYA, suggesting an early differentiation of sex chromosomes that predates most extant avian orders. In agreement with the notion of male-driven evolution, we find a faster rate of change in male-linked CHD-Z sequences.

Modeling evolution at the protein level using an adjustable amino acid fitness model.

An adjustable fitness model for amino acid site substitutions is investigated. This model, a generalization of previously developed evolutionary models, has several distinguishing characteristics: it separately accounts for the processes of mutation and substitution, allows for heterogeneity among substitution rates and among evolutionary constraints, and does not make any prior assumptions about which sites or characteristics of proteins are important to molecular evolution. While the model has fewer adjustable parameters than the general reversible mtREV model, when optimized it outperforms mtREV in likelihood analysis on protein-coding mitochondrial genes. In addition, the optimized fitness parameters of the model show correspondence to some biophysical characteristics of amino acids.

The leech and the physician: biology, etymology, and medical practice with Hirudinea medicinalis.

The history of the word "leech" and the practice of leeching reveal interconnected social histories. We give the linguistic and medical histories of the word, and explore its biology and clinical history. Our historical account extends from the earliest known record of leeching to current research. Despite historical variation in its reputation as a therapeutic technique, leeching remains useful today in a number of applications. Further investigation may well disclose even more uses for the leech, particularly for its enzymes with anesthetic, anticoagulant, and antimetastatic properties.

Phylogenetic position of turtles among amniotes: evidence from mitochondrial and nuclear genes.

Maximum likelihood analysis, accounting for site-heterogeneity in evolutionary rate with the Gamma-distribution model, was carried out with amino acid sequences of 12 mitochondrial proteins and nucleotide sequences of mitochondrial 12S and 16S rRNAs from three turtles, one squamate, one crocodile, and eight birds. The analysis strongly suggests that turtles are closely related to archosaurs (birds+crocodilians), and it supports both Tree-2: (((birds, crocodilians), turtles), squamates) and Tree-3: ((birds, (crocodilians, turtles)), squamates). A more traditional Tree-1: (((birds, crocodilians), squamates), turtles) and a tree in which turtles are basal to other amniotes were rejected with high statistical significance. Tree-3 has recently been proposed by Hedges and Poling [Science 283 (1999) 998-1001] based mainly on nuclear genes. Therefore, we re-analyzed their data using the maximum likelihood method, and evaluated the total evidence of the analyses of mitochondrial and nuclear data sets. Tree-1 was again rejected strongly. The most likely hypothesis was Tree-3, though Tree-2 remained a plausible candidate.

Cospeciation and horizontal transmission of avian sarcoma and leukosis virus gag genes in galliform birds.

In a study of the evolution and distribution of avian retroviruses, we found avian sarcoma and leukosis virus (ASLV) gag genes in 26 species of galliform birds from North America, Central America, eastern Europe, Asia, and Africa. Nineteen of the 26 host species from whom ASLVs were sequenced were not previously known to contain ASLVs. We assessed congruence between ASLV phylogenies based on a total of 110 gag gene sequences and ASLV-host phylogenies based on mitochondrial 12S ribosomal DNA and ND2 sequences to infer coevolutionary history for ASLVs and their hosts. Widespread distribution of ASLVs among diverse, endemic galliform host species suggests an ancient association. Congruent ASLV and host phylogenies for two species of Perdix, two species of Gallus, and Lagopus lagopus and L. mutus also indicate an old association with vertical transmission and cospeciation for these ASLVs and hosts. An inference of horizontal transmission of ASLVs among some members of the Tetraoninae subfamily (grouse and ptarmigan) is supported by ASLV monophyletic groups reflecting geographic distribution and proximity of hosts rather than host species phylogeny. We provide a preliminary phylogenetic taxonomy for the new ASLVs, in which named taxa denote monophyletic groups.

Primers for a PCR-based approach to mitochondrial genome sequencing in birds and other vertebrates.

A PCR-based approach to sequencing complete mitochondrial genomes is described along with a set of 86 primers designed primarily for avian mitochondrial DNA (mtDNA). This PCR-based approach allows an accurate determination of complete mtDNA sequences that is faster than sequencing cloned mtDNA. The primers are spaced at about 500-base intervals along both DNA strands. Many of the primers incorporate degenerate positions to accommodate variation in mtDNA sequence among avian taxa and to reduce the potential for preferential amplification of nuclear pseudogenes. Comparison with published vertebrate mtDNA sequences suggests that many of the primers will have broad taxonomic utility. In addition, these primers should make available a wider variety of mitochondrial genes for studies based on smaller data sets.

Using physical-chemistry-based substitution models in phylogenetic analyses of HIV-1 subtypes.

HIV-1 subtype phylogeny is investigated using a previously developed computational model of natural amino acid site substitutions. This model, based on Boltzmann statistics and Metropolis kinetics, involves an order of magnitude fewer adjustable parameters than traditional substitution matrices and deals more effectively with the issue of protein site heterogeneity. When optimized for sequences of HIV-1 envelope (env) proteins from a few specific subtypes, our model is more likely to describe the evolutionary record for other subtypes than are methods using a single substitution matrix, even a matrix optimized over the same data. Pairwise distances are calculated between various probabilistic ancestral subtype sequences, and a distance matrix approach is used to find the optimal phylogenetic tree. Our results indicate that the relationships between subtypes B, C, and D and those between subtypes A and H may be closer than previously thought.

Assessing the Cretaceous superordinal divergence times within birds and placental mammals by using whole mitochondrial protein sequences and an extended statistical framework.

Using the set of all vertebrate mtDNA protein sequences published as of May 1998, plus unpublished examples for elephant and birds, we examined divergence times in Placentalia and Aves. Using a parsimony-based test, we identified a subset of slower evolutionary rate placental sequences that do not appear to violate the clock assumption. Analyzing just these sequences decreases support for Marsupionta and the carnivore + perissodactyl group but increases support for armadillo diverging earlier than rabbit (which may represent the whole Glires group). A major theme of the paper is to use more comprehensive estimates of divergence time standard error (SE). From the well-studied horse/rhino split, estimated to be 55 million years before present (mybp), the splitting time within carnivores is confidently shown to be older than 50 million years. Some of our estimates of divergence times within placentals are relatively old, at up to 169 million years, but are within 2 SE of other published estimates. The whale/cow split at 65 mybp may be older than commonly assumed. All the sampled splits between the main groups of fereuungulates (the clade of carnivores, cetartiodactyls, perissodactyls, and pholidotes) seem to be distinctly before the Cretaceous/Tertiary boundary. Analyses suggest a close relationship between elephants (representing Afrotheria) and armadillos (Xenarthra), and our timing of this splitting is coincident with the opening of the South Atlantic, a major vicariant event. Recalibrating with this event (at 100 mybp), we obtain younger estimates for the earliest splits among placentals. Divergence times within birds are also assessed by using previously unpublished sequences. We fail to reject a clock for all bird taxa available. Unfortunately, available deep calibration points for birds are questionable, so a new calibration based on the age of the Anseriform stem lineage is estimated. The divergence time of rhea and ostrich may be much more recent than commonly assumed, while that of passerines may be older. Our major concern is the rooting point of the bird subtree, as the nearest outgroup (alligator) is very distant.

Interordinal relationships of birds and other reptiles based on whole mitochondrial genomes.

Several different groups of birds have been proposed as being the oldest or earliest diverging extant lineage within the avian phylogenetic tree, particularly ratites (Struthioniformes), waterfowl (Anseriformes), and shorebirds (Charadriiformes). Difficulty in resolving this issue stems from several factors, including the relatively rapid radiation of primary (ordinal) bird lineages and the lack of characters from an extant outgroup for birds that is closely related to them by measure of time. To help resolve this question, we have sequenced entire mitochondrial genomes for five birds (a rhea, a duck, a falcon, and two perching birds), one crocodilian, and one turtle. Maximum parsimony and maximum likelihood analyses of these new sequences together with published sequences (18 taxa total) yield the same optimal tree topology, in which a perching bird (Passeriformes) is sister to all the other bird taxa. A basal position for waterfowl among the bird study taxa is rejected by maximum likelihood analyses. However, neither the conventional view, in which ratites (including rhea) are basal to other birds, nor tree topologies with falcon or chicken basal among birds could be rejected in the same manner. In likelihood analyses of a subset of seven birds, alligator, and turtle (9 taxa total), we find that increasing the number of parameters in the model shifts the optimal topology from one with a perching bird basal among birds to the conventional view with ratites diverging basally; moreover, likelihood scores for the two trees are not significantly different. Thus, although our largest set of taxa and characters supports a tree with perching birds diverging basally among birds, the position of this earliest divergence among birds appears unstable. Our analyses indicate a sister relationship between a waterfowl/chicken clade and ratites, relative to perching birds and falcon. We find support for a sister relationship between turtles and a bird/crocodilian clade, and for rejecting both the Haemothermia hypothesis (birds and mammals as sister taxa) and the placement of turtles as basal within the phylogenetic tree for amniote animals.

Multiple independent origins of mitochondrial gene order in birds.

Mitochondrial genomes of all vertebrate animals analyzed to date have the same 37 genes, whose arrangement in the circular DNA molecule varies only in the relative position of a few genes. This relative conservation suggests that mitochondrial gene order characters have potential utility as phylogenetic markers for higher-level vertebrate taxa. We report discovery of a mitochondrial gene order that has had multiple independent originations within birds, based on sampling of 137 species representing 13 traditionally recognized orders. This provides evidence of parallel evolution in mitochondrial gene order for animals. Our results indicate operation of physical constraints on mitochondrial gene order changes and support models for gene order change based on replication error. Bird mitochondria have a displaced OL (origin of light-strand replication site) as do various other Reptilia taxa prone to gene order changes. Our findings point to the need for broad taxonomic sampling in using mitochondrial gene order for phylogenetic analyses. We found, however, that the alternative mitochondrial gene orders distinguish the two primary groups of songbirds (order Passeriformes), oscines and suboscines, in agreement with other molecular as well as morphological data sets. Thus, although mitochondrial gene order characters appear susceptible to some parallel evolution because of mechanistic constraints, they do hold promise for phylogenetic studies.

Dihydropteroate synthase polymorphisms in Pneumocystis carinii.

Sulfa drugs are widely used in the treatment and prophylaxis of Pneumocystis carinii pneumonia. The nucleotide sequences of the sulfa target enzyme, dihydropteroate synthase (DHPS), differed substantially in human-, rat-, and mouse-derived P. carinii. Sequence variation also existed in the DHPSs from human-derived isolates. Six nucleotide changes were found in 6 human isolates; each was nonsynonymous and resulted in an amino acid change. Several of these changes were in highly conserved regions and are similar to those that cause sulfa resistance in other organisms. These data suggest that the human-derived P. carinii DHPS may be evolving under positive selective pressure from sulfa drugs.

Mitochondrial gene order adjacent to the control region in crocodile, turtle, and tuatara.

We used the polymerase chain reaction and sequencing of mitochondrial gene junctions adjacent to the 5' end of the control region (light strand orientation) and the 3' end of ND6 to assess whether a representative crocodilian, turtle, or tuatara shares a unique mitochondrial gene order that is found in birds but not in mammals or amphibians. Turtles and crocodiles have the same gene order as mammals, except that crocodile has a tRNAPhe gene inserted between tRNAPro and the 5' end of the control region. Two different arrangements were detected at the 5' end of the control region in the tuatara, one resembling the mammalian (but with tRNAThr absent) and one resembling the avian gene order. Our data are consistent with the hypothesis that some tuatara mtDNAs within a single individual have undergone a deletion that removed the genes coding for cytochrome b and tRNAPro as well as 87 bp of the control region.

Positive selection and rates of evolution in immunodeficiency viruses from humans and chimpanzees.

Evolutionary theory predicts the recent spread of primate immunodeficiency viruses (PIVs) to new human populations to be accompanied by positive selection in response to new host environments and/or by random genetic drift. I assess evidence for positive selection in human and chimpanzee PIVs type I (PIV1s), using ratios of synonymous to nonsynonymous nucleotide change based on branch lengths and outgroup rooting. Ratios are smaller for PIV1s from humans than for PIV1 from a chimpanzee for the pol, gag, and env glycoprotein 120 (gp120) regions, indicating greater effects of positive selection in PIV1s from humans. Parsimony-based relative rate tests for amino acid changes showed significant differences between PIV1s from humans and chimpanzees in 18 of 48 pairwise comparisons, with all 18 showing faster rates of change in PIV1s from humans. This study indicates that in some instances, the recent evolution of human PIV1s follows a speciational pattern, in which increased diversification of taxa is correlated with greater amounts of character change appearing and being maintained through time. This extends the generality of the speciational pattern to a group of organisms (viruses) having the fastest known rates of anagenetic change for nucleotide characters and indicates that comprehensive understanding of PIV1 evolution requires consideration of both anagenetic change within viral lineages and the relative historical success of different viral clades. Phylogenetic analyses show that neither PIV1s infecting humans nor those infecting chimpanzees represent monophyletic groups and suggest multiple host-species shifts for PIV1s.

Phylogenetic comparison of two polycyclic aromatic hydrocarbon-degrading mycobacteria.

Two mycobacterial strains previously isolated from fossil-fuel-contaminated environments and shown to degrade four- and/or five-ring polycyclic aromatic hydrocarbons were further characterized. The two strains, PYR-I and RJGII-135, had similar growth characteristics, colony morphologies, and scotochromogenic pigmentations. DNA amplification fingerprints obtained with total genomic DNA indicated some strain similarities but with several distinctly different bands. Moreover, phylogenetic analysis based upon essentially full-length 16S rRNA gene sequences separates the two strains as distinct species within the fast-growing group of mycobacteria. Although both strains are thermosensitive, strain PYR-I has the bulged U between positions 184 and 193 characteristic of thermotolerant mycobacteria. Both strains are of potential use for reintroduction into and bioremediation of polycyclic aromatic hydrocarbon-contaminated soils.

Evolution of the WANCY region in amniote mitochondrial DNA.

In most vertebrate mitochondrial genomes, the site for initiation of light-strand replication, OL, is found within a cluster of five transfer RNA (tRNA) genes (tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), and tRNA(Tyr)). This region and part of the adjacent cytochrome c oxydase subunit I (COI) gene were sequenced for two crocodilian, two turtle, and one snake species and for Sphenodon punctatus; part of the adjacent nicotinamide adenine dinucleotide dehydrogenase subunit 2 (ND2) gene was also sequenced for the crocodilian and turtle species. All had the typical vertebrate gene order. The turtles and the snake have a lengthy noncoding sequence between the tRNA(Asn) and tRNA(Cys) genes that we assumed to be homologous to the mammalian OL. The crocodilians and Sphenodon lack such a sequence, a condition they share with birds. Most proposed phylogenies for the amniotes require that OL at this position was lost at least twice during their diversification or was evolved independently more than once. Within the five tRNA genes, frequencies of substitutions are much higher in loops than in stems. Many loops vary dramatically in size among the species; in the most extreme case, the D-arm of the Sphenodon tRNA(Cys) is a "D-arm replacement" loop of seven nucleotides. Frequency of transitions in stems is relatively uniform across tRNAs, but frequency of transversions varies greatly. Mismatches in stems are infrequent, and their relative frequency in a specific tRNA is unrelated to the frequency of substitution in the corresponding gene. Several features of mammalian mitochondrial tRNAs are conserved in WANCY tRNAs throughout amniotes. The inferred initiation codon for COI is GTG in crocodilians, turtles, and the snake, a condition they share with fishes, certain amphibians, and birds. TTG appears to be the initiation codon for COI in Sphenodon; if correct, this would be a novel initiation codon for vertebrate mitochondrial DNA. Phylogenetic analyses of the inferred amino acid sequences of ND2 and COI support the sister-group relationship of birds and crocodilians and suggest that mammals are an early derived lineage within the amniotes.