Diversity, Distribution, and Ancient Taxonomic Relationships Within the TIR and Non-TIR NBS-LRR Resistance Gene Subfamilies

Phylogenetic relationships among the NBS-LRR (nucleotide binding site–leucine-rich repeat) resistance gene homologues (RGHs) from 30 genera and nine families were evaluated relative to phylogenies for these taxa. More than 800 NBS-LRR RGHs were analyzed, primarily from Fabaceae, Brassicaceae, Poaceae, and Solanaceae species, but also from representatives of other angiosperm and gymnosperm families. Parsimony, maximum likelihood, and distance methods were used to classify these RGHs relative to previously observed gene subfamilies as well as within more closely related sequence clades. Grouping sequences using a distance cutoff of 250 PAM units (point accepted mutations per 100 residues) identified at least five ancient sequence clades with representatives from several plant families: the previously observed TIR gene subfamily and a minimum of four deep splits within the non-TIR gene subfamily. The deep splits in the non-TIR subfamily are also reflected in comparisons of amino acid substitution rates in various species and in ratios of nonsynonymous-to-synonymous nucleotide substitution rates (K _A/K _S values) in Arabidopsis thaliana. Lower K _A/K _S values in the TIR than the non-TIR sequences suggest greater functional constraints in the TIR subfamily. At least three of the five identified ancient clades appear to predate the angiosperm–gymnosperm radiation. Monocot sequences are absent from the TIR subfamily, as observed in previous studies. In both subfamilies, clades with sequences separated by approximately 150 PAM units are family but not genus specific, providing a rough measure of minimum dates for the first diversification event within these clades. Within any one clade, particular taxa may be dramatically over- or underrepresented, suggesting preferential expansions or losses of certain RGH types within particular taxa and suggesting that no one species will provide models for all major sequence types in other taxa. Received: 13 June 2001 / Accepted: 22 October 2001     

Rapid Evolution of the Ribonuclease A Superfamily: Adaptive Expansion of Independent Gene Clusters in Rats and Mice

The two eosinophil ribonucleases, eosinophil-derived neurotoxin (EDN/RNase 2) and eosinophil cationic protein (ECP/RNase 3), are among the most rapidly evolving coding sequences known among primates. The eight mouse genes identified as orthologs of EDN and ECP form a highly divergent, species-limited cluster. We present here the rat ribonuclease cluster, a group of eight distinct ribonuclease A superfamily genes that are more closely related to one another than they are to their murine counterparts. The existence of independent gene clusters suggests that numerous duplications and diversification events have occurred at these loci recently, sometime after the divergence of these two rodent species (∼10–15 million years ago). Nonsynonymous substitutions per site (d _N) calculated for the 64 mouse/rat gene pairs indicate that these ribonucleases are incorporating nonsilent mutations at accelerated rates, and comparisons of nonsynonymous to synonymous substitution (d _N / d _S) suggest that diversity in the mouse ribonuclease cluster is promoted by positive (Darwinian) selection. Although the pressures promoting similar but clearly independent styles of rapid diversification among these primate and rodent genes remain uncertain, our recent findings regarding the function of human EDN suggest a role for these ribonucleases in antiviral host defense. Received: 8 April 1999 / Accepted: 22 June 1999     

Accelerated Evolution and Molecular Surface of Venom Phospholipase A2 Enzymes

Multiple phospholipase A₂ (PLA₂) isoenzymes found in a single snake venom induce a variety of pharmacological effects. These multiple forms are formed by gene duplication and accelerated evolution of exons. We examined the amino acid sequences of 127 snake venom PLA₂ enzymes and their homologues to study in which location most natural substitutions occur. Our data show that hot spots of amino acid substitutions in this group of proteins occur mostly on the surface. A logistic model correlating the substitution rates of each amino acid residue with their surface accessibility indicates that the probability of natural substitutions occurring in the fully exposed residue is 2.6–3.5 times greater than that of substitutions occurring in buried residues. These surface substitutions play a significant role in the evolution of new PLA₂ isoenzymes by altering the specificity of targeting to various tissues or cells, resulting in distinct pharmacological effects. Thus natural substitutions in PLA₂ enzymes, in contrast to popular belief, are not random substitutions but appear to be directed toward modifying the molecular surface. Received: 11 May 1998 / Accepted: 29 June 1998     

A Survey of the Molecular Evolutionary Dynamics of Twenty-Five Multigene Families from Four Grass Taxa

We surveyed the molecular evolutionary characteristics of 25 plant gene families, with the goal of better understanding general processes in plant gene family evolution. The survey was based on 247 GenBank sequences representing four grass species (maize, rice, wheat, and barley). For each gene family, orthology and paralogy relationships were uncertain. Recognizing this uncertainty, we characterized the molecular evolution of each gene family in four ways. First, we calculated the ratio of nonsynonymous to synonymous substitutions (d _N/d _S) both on branches of gene phylogenies and across codons. Our results indicated that the d _N/d _S ratio was statistically heterogeneous across branches in 17 of 25 (68%) gene families. The vast majority of d _N/d _S estimates were <<1.0, suggestive of selective constraint on amino acid replacements, and no estimates were >1.0, either across phylogenetic lineages or across codons. Second, we tested separately for nonsynonymous and synonymous molecular clocks. Sixty-eight percent of gene families rejected a nonsynonymous molecular clock, and 52% of gene families rejected a synonymous molecular clock. Thus, most gene families in this study deviated from clock-like evolution at either synonymous or nonsynonymous sites. Third, we calculated the effective number of codons and the proportion of G+C synonymous sites for each sequence in each gene family. One or both quantities vary significantly within 18 of 25 gene families. Finally, we tested for gene conversion, and only six gene families provided evidence of gene conversion events. Altogether, evolution for these 25 gene families is marked by selective constraint that varies among gene family members, a lack of molecular clock at both synonymous and nonsynonymous sites, and substantial variation in codon usage. Received: 25 May 2000 / Accepted: 16 October 2000     

Expected Relationship Between the Silent Substitution Rate and the GC Content: Implications for the Evolution of Isochores

The relationship between the silent substitution rate (K _s) and the GC content along the genome is a focal point of the debate about the origin of the isochore structure in vertebrates. Recent estimation of the silent substitution rate showed a positive correlation between K _s and GC content, in contradiction with the predictions of both the regional mutation bias model and the selection or biased gene conversion model. The aim of this paper is to help resolve this contradiction between theoretical studies and data. We analyzed the relationship between K _s and GC content under (1) uniform mutation bias, (2) a regional mutation bias, and (3) mutation bias and selection. We report that an increase in K _s with GC content is expected under mutation bias because of either nonequilibrium of the isochore structure or an increasing mutation rate from AT toward GC nucleotides in GC-richer isochores. We show by simulations that CpG deamination tends to increase the mutation rate with GC content in a regional mutation bias model. We also demonstrate that the relationship between K _s and GC under the selectionist or biased gene conversion model is positive under weak selection if the mutation selection equilibrium GC frequency is less than 0.5. Received: 28 March 2001 / Accepted: 16 May 2001     

Evidence for horizontal transfer from streptococcus to escherichia coli of the kfid gene encoding the K5-specific UDP-glucose dehydrogenase

Capsular polysaccharides are important virulence factors both in Gram-positive and Gram-negative bacteria. A similar cluster organization of the genes involved in the synthesis of bacterial exopolysaccharides has been postulated in both cases, suggesting that these clusters evolved by module assembly. Horizontal gene transfer has been postulated to explain the polymorphism found in these cellular polymers. The cap1K and cap3A genes coding for the pneumococcal type 1 and type 3 UDP-glucose dehydrogenases, respectively, have been compared with other UDP-sugar dehydrogenases. We have observed that the evolutionary distance between Cap1K and Cap3A is approximately equal to that found between Cap1K (or Cap3A) and other UDP-GlcDH of families evolutionarily distant like KfiD, the dehydrogenase from Escherichia coli K5. On the basis of comparisons of G + C content, patterns of synonymous and nonsynonymous substitutions, dinucleotide frequencies, and codon usage bias, we conclude that the kfiD gene has been introduced into E. coli from an exogenous source, probably from a streptococcal species. Received: 26 May 1997 / Accepted: 30 July 1997     

Evolutionary History of 4.5SI RNA and Indication That It Is Functional

To date, the small nuclear 4.5S_I RNA has only been studied in the rat (Rattus norvegicus). Combining PCR and hybridization analyses, we have revealed 4.5S_I RNA homologues sequences in the genomes of four myomorph rodent families (Muridae, Cricetidae, Spalicidae, and Rhizomyidae), and not in other myomorph families (Dipodidae, Zapodidae, Geomyidae, and Heteromyidae) or sciuromorph and caviomorph rodents. By Northern-hybridization, 4.5S_I RNA has been detected in the common rat (R. norvegicus, Muridae), golden hamster (Mesocricetus auratus, Cricetidae), and Russian mole rat (Spalax microphthalmus, Spalacidae), but not in the related great jerboa (Allactaga jaculus, Dipodidae) or in four non-myomorph rodent species tested. cDNA derived from 4.5S_I RNA of M. auratus and S. microphthalmus has been cloned and sequenced. The hamster RNA is found to differ from rat 4.5S_I RNA by only one nucleotide substitution. For the mole rat, two variants of 4.5S_I RNA are detected: short (S) and long (L) with length 101 and 108 nt, respectively. The L variant differs from the S variant as well as from murid and cricetid 4.5S_I RNAs by both a 7 nt insertion and a varying number of nucleotide substitutions. The sequence similarity between the spalacid S-variant and murid/crecitid variants of 4.5S_I RNA is 90%. Judging from species distribution, 4.5S_I RNA genes emerged during the same period of time as the related short interspersed element B2 arose. This occurred after the divergence of Dipodidae lineage but before the branching of Spalicidae/Rhizomyidae lineage from a common myomorph rodent stem. S variant genes seemed to emerge in a common ancestor of spalacids and rhizomyds whereas L variant genes formed in spalacids following the divergence of these two families. The low rate of evolutionary changes of 4.5S_I RNA, at least, in murids and cricetids (6 × 10⁻⁴ substitutions per site per million years), suggests that this RNA is under selection constraint and have a function. This is a remarkable fact if the recent origin and narrow species distribution range of 4.5S_I RNA genes is taken into account. Genes with narrow species distribution are proposed to be referred to as stenogenes. Received: 11 December 2000 / Accepted: 27 August 2001     

Blood-Brain Barrier Permeation: Molecular Parameters Governing Passive Diffusion

53 compounds with clinically established ability to cross or not to cross the blood-brain barrier by passive diffusion were characterized by means of surface activity measurements in terms of three parameters, i.e., the air-water partition coefficient, K _aw , the critical micelle concentration, CMC _D , and the cross-sectional area, A _D . A three-dimensional plot in which the surface area, A _D , is plotted as a function of K ⁻¹ _aw and CMC _D shows essentially three groups of compounds: (i) very hydrophobic compounds with large air-water partition coefficients and large cross-sectional areas, A _D > 80 ?² which do not cross the blood-brain barrier, (ii) compounds with lower air-water partition coefficients and an average cross-sectional area, A _D ≅ 50 ?² which easily cross the blood-brain barrier, and (iii) hydrophilic compounds with low air-water partition coefficients (A _D < 50 ?²) which cross the blood-brain barrier only if applied at high concentrations. It was shown that the lipid membrane-water partition coefficient, K _lw , measured previously, can be correlated with the air-water partition coefficient if the additional work against the internal lateral bilayer pressure, π _bi = 34 ± 4 mN/m is taken into account. The partitioning into anisotropic lipid membranes decreases exponentially with increasing cross-sectional areas, A _D , according to K _lw =const. K _aw exp(−A _D π _bi /kT) where kT is the thermal energy. The cross-sectional area of the molecule oriented at a hydrophilic-hydrophobic interface is thus the main determinant for membrane permeation provided the molecule is surface active and has a pK _a > 4 for acids and a pK _a < 10 for bases. Received: 7 April 1998/Revised: 25 June 1998     

Nucleotide divergence of the rp49 gene region between Drosophila melanogaster and two species of the Obscura group of Drosophila

Summary A 2.1-kb SStI fragment including the rp49 gene and the 3 end of the -serendipity gene has been cloned and sequenced in Drosophila pseudoobscura. rp49 maps at region 62 on the tip of chromosome II of this species. Both the coding and flanking regions have been aligned and compared with those of D. subobscura. There is no evidence for heterogeneity in the rate of silent substitution between the rp49 coding region and the rate of substitutions in flanking regions, the overall silent divergence per site being 0.19. Noncoding regions also differ between both species by different insertions/deletions, some of which are related to repeated sequences. The rp49 region of D. pseudoobscura shows a strong codon bias similar to those of D. subobscura and D. melanogaster. Comparison of the rates of silent (K _S ) and nonsilent (K _a ) substitutions of the rp49 gene and other genes completely sequenced in D. pseudoobscura and D. melanogaster confirms previous results indicating that rp49 is evolving slowly both at silent and nonsilent sites. According to the data for the rp49 region, D. pseudoobscura and D. subobscura lineages would have diverged some 9 Myr ago, if one assumes a divergence time of 30 Myr for the melanogaster and obscura groups.Offprint requests to: C. Segarra     

A Kinetic Model with Ordered Cytoplasmic Dissociation for SUC1, an Arabidopsis H+/Sucrose Cotransporter Expressed in Xenopus Oocytes

To elucidate the kinetic properties of the Arabidopsis H⁺/sucrose cotransporter, SUC1, with respect to transmembrane voltage and ligand concentrations, the transport system was heterologously expressed in Xenopus laevis oocytes. Steady-state plasma membrane currents associated with transport of sucrose were measured with two-electrode voltage clamp over the voltage range −180 to +40 mV as a function of extracellular pH and sugar concentrations. At any given voltage, currents exhibited hyperbolic kinetics with respect to extracellular H⁺ and sugar concentrations, and this enabled determination of values for the maximum currents in the presence of each ligand (i ^H _max , i ^S _max for H⁺ and sucrose) and of the ligand concentrations eliciting half-maximal currents (K ^H _m , K ^S _m ). The i ^H _max and i ^S _max exhibited marked and statistically significant increases as a function of increasingly negative membrane potential. However, the K ^H _m and K ^S _m decreased with increasingly negative membrane potential. Furthermore, at any given voltage, i ^S _max increased and K ^S _m decreased as a function of the external H⁺ concentration. Eight six-state carrier models—which comprised the four possible permutations of intracellular and extracellular ligand binding order, each with charge translocation on the sugar-loaded or -unloaded forms of the carrier—were analyzed algebraically with respect to their competence to account for the ensemble of kinetic observations. Of these, two models (first-on, first-off and last-on, first-off with respect to sucrose binding as it passes from outside to inside the cell and with charge translocation on the loaded form of the carrier) exhibit sufficient kinetic flexibility to describe the observations. Combining these two, a single model emerges in which the binding on the external side can be random, but it can only be ordered on the inside, with the sugar dissociating before the proton. Received: 23 January 1996/Revised: 16 April 1997     

Mutations in the M4 Domain of the Torpedo californica Nicotinic Acetylcholine Receptor Alter Channel Opening and Closing

We studied the functional effects of single amino acid substitutions in the postulated M4 transmembrane domains of Torpedo californica nicotinic acetylcholine receptors (nAChRs) expressed in Xenopus oocytes at the single-channel level. At low ACh concentrations and cold temperatures, the replacement of wild-type α418Cys residues with the large, hydrophobic amino acids tryptophan or phenylalanine increased mean open times 26-fold and 3-fold, respectively. The mutation of a homologous cysteine in the β subunit (β447Trp) had similar but smaller effects on mean open time. Coexpression of α418Trp and β447Trp had the largest effect on channel open time, increasing mean open time 58-fold. No changes in conductance or ion selectivity were detected for any of the single subunit amino acid substitutions tested. However, the coexpression of the α418Trp and β447Trp mutated subunits also produced channels with at least two additional conductance levels. Block by acetylcholine was apparent in the current records from α418Trp mutants. Burst analysis of the α418Trp mutations showed an increase in the channel open probability, due to a decrease in the apparent channel closing rate and a probable increase in the effective opening rate. Our results show that modifications in the primary structure of the α- and β subunit M4 domain, which are postulated to be at the lipid-protein interface, can significantly alter channel gating, and that mutations in multiple subunits act additively to increase channel open time. Received: 27 September 1996/Revised: 28 January 1997     

GABA Concentration Sets the Conductance of Delayed GABAA Channels in Outside-out Patches from Rat Hippocampal Neurons

GABA_A channels were activated by GABA in outside-out patches from rat cultured hippocampal neurons. They were blocked by bicuculline and potentiated by diazepam. In 109 of 190 outside-out patches, no channels were active before exposure to GABA (silent patches). The other 81 patches showed spontaneous channel activity. In patches containing spontaneous channel activity, rapid application of GABA rapidly activated channels. In 93 of the silent patches, channels could be activated by GABA but only after a delay that was sometimes as long as 10 minutes. The maximum channel conductance of the channels activated after a delay increased with GABA concentration from less than 10 pS (0.5 μm GABA) to more than 100 pS (10 mm GABA). Fitting the data with a Hill-type equation gave an EC ₅₀ value of 33 μm and a Hill coefficient of 0.6. The channels showed outward rectification and were chloride selective. In the presence of 1 μm diazepam, the GABA EC ₅₀ decreased to 0.2 μm but the maximum conductance was unchanged. Diazepam decreased the average latency for channel opening. Bicuculline, a GABA antagonist, caused a concentration-dependent decrease in channel conductance. In channels activated with 100 μm GABA the bicuculline IC ₅₀ was 19 μm. The effect of GABA on channel conductance shows that the role of the ligand in GABA_A receptor channel function is more complex than previously thought. Received: 23 October 2000/Revised: 27 February 2001     

Identification of a Region Critically Involved in the Interaction of Phlorizin with the Rabbit Sodium-D-Glucose Cotransporter SGLT1

In order to define potential interaction sites of SGLT1 with the transport inhibitor phlorizin, mutagenesis studies were performed in a hydrophobic region of loop 13 (aa 604–610), located extracellularly, close to the C-terminus. COS 7 cells were transiently transfected with the mutants and the kinetic parameters of α-methyl-d-glucopyranoside (AMG) uptake into the cells were investigated. Replacement of the respective amino acids with lysine reduced the maximal uptake rate: Y604K showed 2.2%, L606K 48.4%, F607K 15.1%, C608K 13.1%, G609K 14.1%, and L610K 17.2% of control. In all mutants the apparent K _i for phlorizin increased at least by a factor of 5 compared to the wild-type K _i of 4.6 ± 0.7 μmol/l; most striking changes were observed for Y604K (K _i = 75.3 ± 19.0 μmol/l) and C608K (K _i = 83.6 ± 13.9 μmol/l). Replacement of these amino acids with a nonpolar amino acid instead of lysine such as in Y604F, Y604G and C608A showed markedly higher affinities for phlorizin. In cells expressing the mutants the apparent affinity of AMG uptake for the sugar was not statistically different from that of the wild type (K _m = 0.8 ± 0.2 mmol/l). These studies suggest that the region between amino acids 604 and 610 is involved in the interaction between SGLT1 and phlorizin, probably by providing a hydrophobic pocket for one of the aromatic rings of the aglucone moiety of the glycoside. Received: 29 March 2001/Revised: 15 June 2001     

Modulation of I A Potassium Current in Adrenal Cortical Cells by a Series of Ten Lanthanide Elements

The modulation of I _A K⁺ current by ten trivalent lanthanide (Ln³⁺) cations spanning the series with ionic radii ranging from 0.99 ? to 1.14 ? was characterized by the whole-cell patch clamp technique in bovine adrenal zona fasciculata (AZF) cells. Each of the ten Ln³⁺s reduced I _A amplitude measured at +20 mV in a concentration-dependent manner. Smaller Ln³⁺s were the most potent and half-maximally effective concentrations (EC₅₀s) varied inversely with ionic radius for the larger elements. Estimation of EC₅₀s yielded the following potency sequence: Lu³⁺ (EC₅₀= 3.0 μm) ≈ Yb³⁺ (EC₅₀= 2.7 μm) > Er³⁺ (EC₅₀= 3.7 μm) ≥ Dy³⁺ (EC₅₀= 4.7 μm) > Gd³⁺ (EC₅₀= 6.7 μm) ≈ Sm³⁺ (EC₅₀= 6.9 μm) > Nd³⁺ (EC₅₀= 11.2 μm) > Pr³⁺ (EC₅₀= 22.3 μm) > Ce³⁺ (EC₅₀= 28.0 μm) > La³⁺ (EC₅₀= 33.7 μm). Ln³⁺s altered selected voltage-dependent gating and kinetic parameters of I _A with a potency and order of effectiveness that paralleled the reduction of I _A amplitude. Ln³⁺s markedly slowed activation kinetics and shifted the voltage-dependence of I _A gating such that activation and steady-state inactivation occurred at more depolarized potentials. In contrast, Ln³⁺s did not measurably alter inactivation or deactivation kinetics and only slightly slowed kinetics of inactivated channels returning to the closed state. Replacement of external Ca²⁺ with Mg²⁺ had no effect on the concentration-dependent inhibition of I _A by Ln³⁺s. In contrast to their action on I _A K⁺ current, Ln³⁺s inhibited T-type Ca²⁺ currents in AZF cells without slowing activation kinetics. These results indicate that Ln³⁺ modulate I _A K⁺ channels through binding to a site on I _A channels located within the electric field but which is not specific for Ca²⁺. They are consistent with a model where Ln³⁺ binding to negative charges on the gating apparatus alters the voltage-dependence and kinetics of channel opening. Ln³⁺s modulate transient K⁺ and Ca²⁺ currents by two fundamentally different mechanisms. Received: 21 January 1997/Revised: 3 April 1998     

Unbiased estimation of the rates of synonymous and nonsynonymous substitution

Summary The current convention in estimating the number of substitutions per synonymous site (K _S ) and per nonsynonymous site (K _A ) between two protein-coding genes is to count each twofold degenerate site as one-third synonymous and two-thirds nonsynonymous because one of the three possible changes at such a site is synonymous and the other two are nonsynonymous. This counting rule can considerably overestimate theK _S value because transitional mutations tend to occur more often than transversional mutations and because most transitional mutations at twofold degenerate sites are synonymous. A new method that gives unbiased estimates is proposed. An application of the new and the old method to 14 pairs of mouse and rat genes shows that the new method gives aK _S value very close to the number of substitutions per fourfold degenerate site whereas the old method gives a value 30% higher. Both methods give aK _A value close to the number of substitutions per nondegenerate site.     

III. Ion Channel Expression in PMA-differentiated Human THP-1 Macrophages

Ion channel expression was studied in THP-1 human monocytic leukemia cells induced to differentiate into macrophage-like cells by exposure to the phorbol ester, phorbol 12-myristate 13-acetate (PMA). Inactivating delayed rectifier K⁺ currents, I _DR, present in almost all undifferentiated THP-1 monocytes, were absent from PMA-differentiated macrophages. Two K⁺ channels were observed in THP-1 cells only after differentiation into macrophages, an inwardly rectifying K⁺ channel (I _IR) and a Ca²⁺-activated maxi-K channel (I _BK). I _IR was a classical inward rectifier, conducting large inward currents negative to E _K and very small outward currents. I _IR was blocked in a voltage-dependent manner by Cs⁺, Na⁺, and Ba²⁺, block increasing with hyperpolarization. Block by Na⁺ and Ba²⁺ was time-dependent, whereas Cs⁺ block was too fast to resolve. Rb⁺ was sparingly permeant. In cell-attached patches with high [K⁺] in the pipette, the single I _IR channel conductance was ∼30 pS and no outward current could be detected. I _BK channels were observed in cell-attached or inside-out patches and in whole-cell configuration. In cell-attached patches the conductance was ∼200–250 pS and at potentials positive to ∼100 mV a negative slope conductance of the unitary current was observed, suggesting block by intracellular Na⁺. I _BK was activated at large positive potentials in cell-attached patches; in inside-out patches the voltage-activation relationship was shifted to more negative potentials by increased [Ca²⁺]. Macroscopic I _BK was blocked by external TEA⁺ with half block at 0.35 mm. THP-1 cells were found to contain mRNA for Kv1.3 and IRK1. Levels of mRNA coding for these K⁺ channels were studied by competitive PCR (polymerase chain reaction), and were found to change upon differentiation in the same direction as did channel expression: IRK1 mRNA increased at least 5-fold, and Kv1.3 mRNA decreased on average 7-fold. Possible functional correlates of the changes in ion channel expression during differentiation of THP-1 cells are discussed. Received: 19 September 1995/Revised: 14 March 1996     

Synonymous and Nonsynonymous Substitutions in Genes from Gramineae: Intragenic Correlations

In this work, we have investigated the relationships between synonymous and nonsynonymous rates and base composition in coding sequences from Gramineae to analyze the factors underlying the variation in substitutional rates. We have shown that in these genes the rates of nucleotide divergence, both synonymous and nonsynonymous, are, to some extent, dependent on each other and on the base composition. In the first place, the variation in nonsynonymous rate is related to the GC level at the second codon position (the higher the GC₂ level, the higher the amino acid replacement rate). The correlation is especially strong with T₂, the coefficients being significant in the three data sets analyzed. This correlation between nonsynonymous rate and base composition at the second codon position is also detectable at the intragenic level, which implies that the factors that tend to increase the intergenic variance in nonsynonymous rates also affect the intragenic variance. On the other hand, we have shown that the synonymous rate is strongly correlated with the GC₃ level. This correlation is observed both across genes and at the intragenic level. Similarly, the nonsynonymous rate is also affected at the intragenic level by GC₃ level, like the silent rate. In fact, synonymous and nonsynonymous rates exhibit a parallel behavior in relation to GC₃ level, indicating that the intragenic patterns of both silent and amino acid divergence rates are influenced in a similar way by the intragenic variation of GC₃. This result, taken together with the fact that the number of genes displaying intragenic correlation coefficients between synonymous and nonsynonymous rates is not very high, but higher than random expectation (in the three data sets analyzed), strongly suggests that the processes of silent and amino acid replacement divergence are, at least in part, driven by common evolutionary forces in genes from Gramineae. Received: 2 July 1998 / Accepted: 18 April 1999     

Multimeric hemoglobin of the Australian brine shrimp Parartemia

The hemoglobin molecule of the commercially important brine shrimp Artemia sp. has been used extensively as a model for the study of molecular evolution. It consists of nine globin domains joined by short linker sequences, and these domains are believed to have originated through a series of duplications from an original globin gene. In addition, in Artemia, two different polymers of hemoglobin, called C and T, are found which differ by 11.7% at the amino acid level and are believed to have diverged about 60 MYA. This provides a set of data of 18 globin domain sequences that have evolved in the same organism. The pattern of amino acid substitution between these two polymers is unusual, with pairs of equivalent domains displaying differences of up to 2.7-fold in total amino acid substitution. Such differences would reflect a similar range of molecular-clock rates in what appear to be duplicate, structurally equivalent domains. In order to provide a reference outgroup, we sequenced the cDNA for a nine-domain hemoglobin (P) from another genus of brine shrimp, Parartemia zietziana, which differs morphologically and ecologically from Artemia and is endemic to Australia. Parartemia produces only one hundredth the amount of hemoglobin that Artemia produces and does not upregulate production in response to low oxygen partial pressure. Comparison of the globin domains at the amino acid and DNA levels suggests that the Artemia globin T gene has accumulated substitutions differently from the Parartemia P and Artemia C globin genes. We discuss the questions of accelerated evolution after duplication and possible functions for the Parartemia globin.     

Phylogenetic Analysis of Reptilian Hemoglobins: Trees, Rates, and Divergences

Phylogenetic relationships among reptiles were examined using previously published and newly determined hemoglobin sequences. Trees reconstructed from these sequences using maximum-parsimony, neighbor-joining, and maximum-likelihood algorithms were compared with a phylogenetic tree of Amniota, which was assembled on the basis of published morphological data. All analyses differentiated α chains into α^A and α^D types, which are present in all reptiles except crocodiles, where only α^A chains are expressed. The occurrence of the α^D chain in squamates (lizards and snakes only in this study) appears to be a general characteristic of these species. Lizards and snakes also express two types of β chains (βI and βII), while only one type of β chain is present in birds and crocodiles. Reconstructed hemoglobin trees for both α and β sequences did not yield the monophyletic Archosauria (i.e., crocodilians + birds) and Lepidosauria (i.e., Sphenodon+ squamates) groups defined by the morphology tree. This discrepancy, as well as some other poorly resolved nodes, might be due to substantial heterogeneity in evolutionary rates among single hemoglobin lineages. Estimation of branch lengths based on uncorrected amino acid substitutions and on distances corrected for multiple substitutions (PAM distances) revealed that relative rates for squamate α^A and α^D chains and crocodilian β chains are at least twice as high as those of the rest of the chains considered. In contrast to these rate inequalities between reptilian orders, little variation was found within squamates, which allowed determination of absolute evolutionary rates for this subset of hemoglobins. Rate estimates for hemoglobins of lizards and snakes yielded 1.7 (α^A) and 3.3 (β) million years/PAM when calibrated with published divergence time vs. PAM distance correlates for several speciation events within snakes and for the squamate ↔ sphenodontid split. This suggests that hemoglobin chains of squamate reptiles evolved ∼3.5 (α^A) or ∼1.7 times (β) faster than their mammalian equivalents. These data also were used to obtain a first estimate of some intrasquamate divergence times. Received: 15 September 1997 / Accepted: 4 February 1998