The adaptation ability of neuronal models subject to a current step stimulus

Three neuronal models of the spike initiating process were investigated with respect to their ability to show adaptation to a current step: (i) the perfect integrator model (PIM), (ii) the leaky integrator model (LIM), and (iii) the Hodgkin-Huxley (HH-) model. It was found that although each neuronal model will generate different response spike trains to a given stimulus, all responses fulfilled the criteria of a deterministic neural response (Awiszus 1988). The results show that both PIM and LIM are unable to show adaptation regardless of the choice of model parameters whereas the HH-model shows a clear rate of discharge adaptation. The reason for this adaptation lies in the fact that there are conditions for the HH-model where a step stimulus is highly effective. These conditions are investigated by means of a phase plane analysis. Consequences of these results for the explanation of neuronal adaptation and the validity of the neuronal models investigated are discussed.     

Experimental and theoretical studies of the optical and electrical characteristics of a high-pressure pulsed discharge in argon

The optical and electrical characteristics of pulsed discharges in pure Ar at pressures of up to 7 atm, at which the discharge becomes unstable, are studied in a simple experimental device with automatic preionization. The gas temperature in the discharge is estimated from the width of the recorded emission spectrum. An analytical model of the vibrational relaxation of Ar ^*₂ (v) is used to better determine the constants of the vibrational-translational relaxation of Ar ^*₂ (v) molecules in their collisions with Ar atoms. The zerodimensional numerical model of a pulsed discharge in Ar is modified. The experimental and calculated results are compared in detail. Good agreement is achieved between the measured and calculated time dependences of the electrode voltage and the intensity of spontaneous emission in the pressure range of 1–6 atm, as well as between the measured and calculated values of the gas temperature at pressures of 3–6 atm. Preliminary results from numerical studies of the possibility of achieving generation are discussed.

     

A stochastic afterhyperpolarization model of repetitive activity in vestibular afferents

A stochastic version of Kernell's (1968, 1972) model with cumulative afterhyperpolarization (AHP) was simulated. A characteristic of the model is that the AHP is the result of an increased potassium conductance (g _K) that is time-dependent but not voltage-dependent. Quantal synaptic inputs are assumed to be the only source of interspike interval variability. The model reproduces many features of the steady-state discharge of peripheral vestibular afferents, provided that firing rates are higher than 40 spikes/s. Among the results accounted for are the interspike interval statistics occurring during natural stimulation, their alteration by externally applied galvanic currents and the increase in the interspike interval following an interposed shock. Empirical studies show that some vestibular afferents have a regular spacing of action potentials, others an irregular spacing (Goldberg and Fernández 1971b; Fernández and Goldberg 1976). Irregularly discharging afferents have a higher sensitivity to externally applied galvanic currents than do regular afferents (Goldberg et al. 1984). To explain the relation between galvanic sensitivity and discharge regularity requires the assumption that neurons differ in both their synaptic noise (_v) and the slopes of their postspike voltage trajectories (d _v/dt). The more irregular the neuron's discharge at a given firing frequency, the greater is _v and the smaller is d _v/dt. Of the two factors, d _v/dt is estimated to be four times more influential in determining discharge regularity across the afferent population. The shortcomings of the model are considered, as are possible remedies. Our conclusions are compared to previous discussions of mechanisms responsible for differences in the discharge regularity of vestibular afferents.     

The initial reaction velocities of lactate dehydrogenase in various cell types

Summary The initial reaction velocities (v _v ) of lactate dehydrogenase in hepatocytes, cardiac muscle fibres, skeletal (gastrocnemius) muscle fibres, gastric parietal cells, ductal epithelial and acinar cells of the parotid gland, and oocytes were determined, by computer-assisted image analysis, in unfixed sections of these tissues incubated at 37°C on substrate-containing agarose gel films. They were found to fit the equations v _i = a₁A (equation 1) and v _i – v = a₂A (equation 2) reported previously for mouse hepatocytes (Nakae & Stoward, 1993a, b), where v and A are, respectively, the gradients (or steady-state velocities) and the intercepts on the absorbance axis of the linear regression lines of the absorbance (A) of the finalreaction product on incubation times between 1 and 3 min, and a ₁ and a ₂ are constants. Both equations 1 and 2 fitted the observed v _i closely for mouse (a ₁ = 2.7, a ₂ = 2.2) and human (a ₁ = 3.0, a ₂ = 1.9) hepatocytes. However, equation 2 fitted the observed v _i better than equation 1 for mouse cardiac muscle fibres (a ₂ = 1.5), skeletal muscle fibres (a ₂ = 1.2), gastric parietal cells (a ₂ = 1.7), acinar (a ₂ = 1.4) and striated ductal (a ₂ = 2.2) epithelial cells of the parotid gland, and oocytes (a ₂ = 1.6). The values of v _i calculated from the two equations agreed with the observed v _i to within about 11%. They ranged from 105 mole hydrogen equivalents/cm³ cell/min units in hepatocytes to 24 units in parotid acinar cells, but for other cell types they were between 46 and 61 units. These are all considerably higher than values reported previously.     

Photosynthetic performance of cold-sensitive mutants of maize at low temperature

The virescent character is a genetic variant in pigmentation characterized by a delay in greening. Seedlings of the virescent mutants v1, v2, v3, v4, v13, v16, v18, v19 and v26 of maize exhibit chlorosis when grown at low temperature. Chlorotic leaves contain plastids that appear to have been arrested at an early stage of development. The results indicated that V16, V2, V3 and V4 loci control early stages of chloroplast development while V1, V13 and V19 may play a role at the end of development. The mutations in the V18 and V26 loci may control an intermediate step. At the pigment level, the virescent mutants of maize differ widely from analogous mutations existing in other plants. The mutations were characterized by a reduced amount of chlorophyll a and b (up to 100 times in v16) and chlorophyll a/b ratio above normal (up to 13.7 in v16). Lutein content was reduced in all mutants (less than 3% in v16 compared to wild type) but v13, while pigments of the xanthophyll cycle were found at higher levels in v1 and v13 (more than 10 and 90%, respectively). The v2, v3, v4, v16 and v18 mutants that are most depleted in β-carotene (36 times less in average than wild type) are also deprived in D1 and D2 polypeptides. Moreover, the v2, v3, v4, v16 and v18 mutants characterized by a lower accumulation in lutein are most depleted of light-harvesting complex II. All mutants possess a functioning, fully reversible, non-photochemical quenching mechanism. This is most developed in the v13 and v19 mutants (φ_n = 0.48 and 0.44, respectively). These two mutants also have a relatively high primary photochemical yield for photosystem II and a functioning photosystem I (φ_p = 0.23 and 0.39, respectively). The most interesting mutant is v13 that shows severe chlorosis and possesses the most effective non-photochemical quenching mechanism(s), which is thought to provide protection against excess photon absorption by photosystem II.     

A mathematical model for wing dimorphism in maleCardiocondyla ants

Several species ofCardiocondyla ants have dimorphic males: wingless (ergatoid) and winged (alate) males, while otherCardiocondyla species includingC. nuda have only wingless males. We made an evolutionarily stable strategy model for explaning the male polymorphism and the ratios of wingless males in the genusCardiocondyla. Wingless males emerge earlier than winged males in each reproductive season. Females (F ₁) which have emerged before winged males copulate only with wingless males, and females (F ₂) which emerge after winged males copulate with both wingless and winged males. Wingless males have a lower copulation ability (b _n) than winged males (b _w). The reproductive success of females which copulate at the early stage (v ₁) is assumed to be larger than that of females which copulate at the late stage (v ₂). The model predicts that there are 3 different evolutionarily stable states: 2 monomorphic states of wingless and winged males, and a dimorphism of the 2 types of males. In the dimorphic state, the rate of wingless males increases as the survival rate of wingless males (s) increases,v ₁/v ₂ increases,F ₁/F ₂ increases andb _n/b _w increases. For dimorphism to exist,s b _n/b _w<1 must be satisfied, and this condition corresponds to the value of observed data. The value ofv ₁/v ₂ would be difficult to be obtained by actual data, but we can estimate this value with the model.     

Cell growth and P(3HB) accumulation from CO<Subscript>2</Subscript> of a carbon monoxide-tolerant hydrogen-oxidizing bacterium, <Emphasis Type="Italic">Ideonella</Emphasis> sp. O-1

Cell growth and accumulation of polyhydroxybutyric acid, P(3HB), from CO₂ in autotrophic condition of a newly isolated hydrogen-oxidizing bacterium, the strain O-1, was investigated. The bacterium, which was deposited in the Japan Collection of Microorganisms as JCM17105, autotrophically grows by assimilating H₂, O₂, and CO₂ as substrate. 16S rRNA gene sequence of the bacterium was the closest to Ideonella dechloratans (99%). Specific growth rate of the strain O-1 was faster than a hydrogen-oxidizing bacterium, Ralstonia eutropha, which is well-known P(3HB)-producing microorganism. The strain O-1 is tolerant to high O₂ concentration and it can grow above 30% (v/v) O₂, while the growth of R. eutropha and Alcaligenes latus was seriously inhibited. In culture medium containing 1 g/L (NH₄)₂SO₄, cell concentration of the strain O-1 and P(3HB) increased to 6.75 and 5.26 g/L, respectively. The content of P(3HB) in the cells was 77.9% (w/w). The strain O-1 was very tolerant to carbon monoxide (CO) and it grew even at 70% (v/v) CO, while the growth of R. eutropha and A. latus were seriously inhibited at 5% (v/v) CO. From these results, it is expected that the strain O-1 will be useful in the manufacture of P(3HB) because the industrial exhaust gas containing CO₂, H₂, and CO can be directly used as the substrate in the fermentation process.     

Analytical reconstruction of the neuronal input current from spike train data

The time course of the current driving action potential generation at a neuron investigated experimentally is in general not measurable directly. In this paper an indirect method is introduced that allows estimation of this unknown current time course using only spike train data. Assuming the leaky integrator model as valid for the action potential encoding site of the investigated neuron, the unknown input current is obtained by determining (analytically) a current time course that upon injection into the leaky integrator model evokes action potential sequences identical to those observed experimentally. Applications of this current-reconstruction procedure to neuronal output data obtained from a leaky integrator model showed that the procedure allows a good estimation of the underlying input current even if the membrane time constant of the investigated neuron is not known exactly. Additionally, an application of current reconstruction to experimental data obtained from a cat muscle spindle primary afferent subject to repeated -stimuli is demonstrated.     

Sodium/calcium selectivity of cloned calcium T-type channels

We studied the peculiarities of permeability with respect to the main extracellular cations, Na⁺ and Ca²⁺, of cloned low-threshold calcium channels (LTCCs) of three subtypes, Ca_v3.1 (α1G), Ca_v3.2 (α 1H), and Ca_v3.3 (α1I), functionally expressed in Xenopus oocytes. In a calcium-free solution containing 100 mM Na⁺ and 5 mM calcium-chelating EGTA buffer (to eliminate residual concentrations of Ca²⁺) we observed considerable integral currents possessing the kinetics of inactivation typical of LTCCs and characterized by reversion potentials of −10 ± 1, −12 ± 1, and −18 ± 2 mV, respectively, for Ca_v3.1, Ca_v3.2, and Ca_v3.3 channels. The presence of Ca²⁺ in the extracellular solution exerted an ambiguous effect on the examined currents. On the one hand, Ca²⁺ effectively blocked the current of monovalent cations through cloned LTCCs (K _d = 2, 10, and 18 μM for currents through channels Ca_v3.1, Ca_v3.2, and Ca_v3.3, respectively). On the other hand, at the concentration of 1 to 100 mM, Ca²⁺ itself functioned as a carrier of the inward current. Despite the fact that the calcium current reached the level of saturation in the presence of 5 mM Ca²⁺ in the external solution, extracellular Na⁺ influenced the permeability of these channels even in the presence of 10 mM Ca²⁺. The Ca_v3.3 channels were more permeable with respect to Na⁺ (P _Ca/P _Na ∼ 21) than Ca_v3.1 and Ca_v3.2 (P _Ca/P _Na ∼ 66). As a whole, our data indicate that cloned LTCCs form multi-ion Ca²⁺-selective pores, as these ions possess a high affinity for certain binding sites. Monovalent cations present together with Ca²⁺ in the external solution modulate the calcium permeability of these channels. Among the above-mentioned subtypes, Ca_v3.3 channels show the minimum selectivity with respect to Ca²⁺ and are most permeable for monovalent cations. Neirofiziologiya/Neurophysiology, Vol. 38, No. 3, pp. 183–192, May–June, 2006.     

An unstable gene controlling developmental variegation in okra (Abelmoschus esculentus (L.) Moench)

Summary Genetic studies on radiation-induced chlorina and variegated mutants of okra (Abelmoschus esculentus (L.) Moench) revealed the existence of an unstable gene. The normal green color of the leaves is controlled by duplicate genes C₁ and C₂, either of which produces the green colour. The chlorina plants are C ₁ C ₁ C ₂ C ₂. The allele c ₁ ^v is dominant to both C ₁ and C ₂ but is unstable. The homozygote c ₁ ^v c ₁ ^v c ₂ c ₂ is a normal green while the heterozygote c _i ^v c ₁ c ₂ c ₂ has a variegated phenotype as a result of the mutation of c ₁ ^v to c ₁ during development. In green plants with a c ₁ ^v c{sh1/v}c ₂ c ₂ genotype, the autonomous mutation of one of the c ₁ ^v alleles to c ₁ may take place at the pre-meiotic stage. In the variegated genotype (c ₁ ^v c ₁ c ₂ c ₂), the mutation of c ₁ to c ₁ ^v may take place in early ontogeny, thus producing green plants. The allele C ₁, when associated with c ₁ ^v in a heterozygous condition, mutates to c ₁ at the pre-meiotic stage even in the presence of the allele C ₂.     

The Frequency Response, Coherence, and Information Capacity of Two Neuronal Models

Two neuronal models are analyzed in which subthreshold inputs are integrated either without loss (perfect integrator) or with a decay which follows an exponential time course (leaky integrator). Linear frequency response functions for these models are compared using sinusoids, Poisson-distributed impulses, or gaussian white noise as inputs. The responses of both models show the nonlinear behavior characteristic of a rectifier for sinusoidal inputs of sufficient amplitude. The leaky integrator shows another nonlinearity in which responses become phase locked to cyclic stimuli. Addition of white noise reduces the distortions due to phase locking. Both models also show selective attenuation of high-frequency components with white noise inputs. Input, output, and cross-spectra are computed using inputs having a broad frequency spectrum. Measures of the coherence and information transmission between the input and output of the models are also derived. Steady inputs, which produce a constant “carrier” rate, and intrinsic sources, which produce variability in the discharge of neurons, may either increase or decrease coherence; however, information transmission using inputs with a broad spectrum is generally increased by steady inputs and reduced by intrinsic variability.     

Chromosome polymorphism in Avena ventricosa

The same karyotype which was described previously in A. ventricosa plants from Oran, Algeria was found in plants from Cyprus and in most of the plants from the Apsheron Peninsula, U.S.S.R. This karyotype, which consists of subterminal chromosomes, was designated C _v ¹ . In the karyotype of about one-third of the Apsheron plants a single pair of submedian chromosomes replaced a pair of long subterminals. This karyotype was designated C _v ² . The chromosomes formed 7^II at meiosis in both C _v ¹ C _v ¹ and C _v ¹ C _v ² hybrids regardless of geographic origin. A heteromorphic bivalent was apparent in PMCs in the C _v ¹ C _v ² hybrids. A pericentric inversion was postulated as the origin of the submedian pair, the C _v ² C _v ² karyotype being the inversion homozygote. The absence of duplication-deficiency gametes was inferred from the good fertility of the heterozygote. Members of the heteromorphic pair were transmitted in equal proportion and the segregates C _v ¹ C _v ¹ , C _v ¹ C _v ² and C _v ² C _v ² were equally fertile; thus the inversion did not handicap reproduction. The small contemporary colonies of A. ventricosa, with peripheral locations and vast distances between them, because of the full homology and interfertility, are considered remnants of a formerly large central population.Contribution No 282 from the Research Station, Central Experimental Farm, Ottawa, Ontario.This paper is dedicated to the memory of Dr. Barna Györffy, whose untimely death on August 5, 1970 ended a distinguished career that was an example of integrity and scholarship to a troubled generation of Hungarian geneticists.     

Discharge dynamics and the production of active particles in a cathode-directed streamer

The emission spectroscopy technique is used to analyze a cathode-directed streamer discharge in air at atmospheric pressure in point-plane geometry at interelectrode distances of up to 100 mm and a high-voltage pulse amplitude of 18 kV. The densities of molecules in the N₂(C ³Π_u, v=0), N ₂ ⁺ (B ²Σ _u ⁺ , v=0) and NO(A ²Σ⁺, v=0) states are determined, and the reduced electric field in the streamer head is estimated. It is shown that the increase in the average electric field in the discharge gap substantially intensifies the production of active particles in the discharge plasma and makes the plasma more homogeneous. This effect is only related to the increase in the fraction of regions with a high electric field in the discharge gap and, as a result, the reduction of the discharge energy losses via rapidly thermalized degrees of freedom. The active particles are only produced in the streamer head, including the case in which the interelectrode gap is bridged by the streamer channel.     

Flight of the honey bee VII: metabolic power versus flight speed relation

The existing experimental data on metabolic power P _m of honey bees are critically discussed, partly corrected for real flight conditions and plotted as a function of flight speed v. New wind tunnel measurements of tethered flight under near-natural conditions are added in the range 3.3<v<5.1 m·s^-1, derived from exhaustion flight measurements. Within this small sector the latter measurements can be characterised by a linear correlation: P _m(mW)=6.72v (m·s^-1)+13.83, the slope of which is significantly different from zero. The over-all P _m(v) curve is significantly not a straight line of zero slope but a U-shaped minimum curve and may be approximated by a second-order polynom: P _m=49.2-8.9v+1.5v ². The same is true for relative metabolic power, P _{m rel (e)} related to empty body mass of 76.5 mg: P _{m rel(e)}=630.0-114.0v+19.2v ² (P _m in mW: P _{m rel} in mW·g^-1; v in m·s^-1). The data support the existence of a U-shaped power-versus-speed curve in bees.Abbreviations bm body mass (mg) - f full - e empty - mu muscles - P _m (mJ·s^-1=mW) metabolic power (input) - P _{m rel} (mW·g^-1) relative metabolic power - P _mec (mW) mechanical power (output) - efficiency (of the flight musculature) - t(s) flight time - v (m·s^-1) relative speed between bee and air     

Theoretical study of BN<Subscript>4</Subscript>: potential precursors of high energy density materials (HEDMs)

Ab initio (MP2) and density functional theory (DFT) methods were used to examine nine isomers of the doublet BN₄ species with the 6-311 + G(d) basis set. To our knowledge, these nine structures are all first reported here. Energy analysis indicates that the C _2v branched structure is the global minimum of potential energy surface. Research results show that the C _2v branched, the cis-linear, the C _4v pyramidal, and the C _S five-membered ring structures are likely to be stable and to be observed experimentally. Among these four kinetically stable species, the last three are suitable to be used as potential precursors of HEDMs due to their high dissociation energies. However, the C _2v bent, the trans-linear, the D ₂ bicyclic, the C _2v four-membered ring, and the C _2v cage structures are kinetically unstable due to their low dissociation or isomerization barriers. Two synthesis pathways of the C _2v branched isomer were located. It seems more feasible to synthesize this species by linear NBN and N₂.     

Measurements of the gas temperature in an oxygen plasma by spectroscopy of the O2( b 1Σ g+ ) → O2( X 3Σ g? ) transition

A method for measuring the gas temperature in an oxygen plasma by spectroscopy of the electronic transition from the O₂(b ¹Σ _g ⁺ , v = 0) metastable state of molecular oxygen into the O₂(X ³Σ _g ⁻ , v = 0) ground state is considered in detail. The method is verified experimentally for the plasma of dc glow discharge in pure oxygen. It is shown that the gas temperature can be determined by analyzing high-resolution spectra of the P branch of this transition, no matter whether its fine structure (^P P and ^P Q branches) is resolved or masked, provided that the rotational structure of the spectrum is resolved. The feasibility of the method proposed in 1999 by P. Maco and P. Veis for determining the gas temperature from the ratio between the intensity maxima of the R and P branches of the O₂(b ¹Σ _g ⁺ , v = 0) → O₂(X ³Σ _g ⁻ , v = 0) transition in a poorly resolved spectrum was studied experimentally. It is shown that, in order to use this method, it is necessary to know the spectrograph instrumental function. The effect of the spatial inhomogeneity of the temperature and concentration of O₂(b ¹Σ _g ⁺ ) molecules on the accuracy of integral (over the plasma volume) measurements of the gas temperature is investigated using spatially resolved spectroscopy of the O₂(b ¹Σ _g ⁺ , v = 0) → O₂(X ³Σ _g ⁻ , v = 0) transition. It is shown that precise measurements of the temperature require that the optical measurement system be thoroughly adjusted in order for the temperature and concentration of the emitting particles to vary insignificantly over the optically selected volume. Original Russian Text ? S.M. Zyryanov, D.V. Lopaev, 2007, published in Fizika Plazmy, 2007, Vol. 33, No. 6, pp. 563–574.     

Dysregulation of Cav1.4 channels disrupts the maturation of photoreceptor synaptic ribbons in congenital stationary night blindness type 2

Mutations in the gene encoding Ca_v1.4, CACNA1F_, are associated with visual disorders including X-linked incomplete congenital stationary night blindness type 2 (CSNB2). In mice lacking Ca_v1.4 channels, there are defects in the development of “ribbon” synapses formed between photoreceptors (PRs) and second-order neurons. However, many CSNB2 mutations disrupt the function rather than expression of Ca_v1.4 channels. Whether defects in PR synapse development due to altered Ca_v1.4 function are common features contributing to the pathogenesis of CSNB2 is unknown. To resolve this issue, we profiled changes in the subcellular distribution of Ca_v1.4 channels and synapse morphology during development in wild-type (WT) mice and mouse models of CSNB2. Using Ca_v1.4-selective antibodies, we found that Ca_v1.4 channels associate with ribbon precursors early in development and are concentrated at both rod and cone PR synapses in the mature retina. In mouse models of CSNB2 in which the voltage-dependence of Ca_v1.4 activation is either enhanced (Ca_v1.4_I756T) or inhibited (CaBP4 KO), the initial stages of PR synaptic ribbon formation are largely unaffected. However, after postnatal day 13, many PR ribbons retain the immature morphology. This synaptic abnormality corresponds in severity to the defect in synaptic transmission in the adult mutant mice, suggesting that lack of sufficient mature synapses contributes to vision impairment in Ca_v1.4_I756T and CaBP4 KO mice. Our results demonstrate the importance of proper Ca_v1.4 function for efficient PR synapse maturation, and that dysregulation of Ca_v1.4 channels in CSNB2 may have synaptopathic consequences.     

Crystals of an antibody FV fragment against an integral membrane protein diffracting to 1.28 Å resolution

The F_v fragment of a monoclonal antibody, 7E2 (IgG₁, κ, murine), which is directed against the integral membrane protein cytochrome c oxidase (EC 1.9.3.1) from Paracoccus denitrificans, was cloned and produced in Escherichia coli. Crystals suitable for highresolution X-ray analysis were obtained by microdialysis under low salt conditions. The crystals belong to the orthorhombic space group P2₁2₁2₁ with unit cell dimensions of a = 51.51 Å, b = 56.15 Å, c = 99.86 Å (1 Å = 0.1 nm) and contain one F _v fragment per asymmetric unit. Using synchrotron radiation diffraction data were collected up to 1.28 Å resolution. This high resolution is very unusual for a heterodimeric protein. The crystals should open the way for refining not only the atomic positions, but also for obtaining information about internal dynamics. © 1995 Wiley-Liss, Inc.     

On transition bias in mitochondrial genes of pocket gophers

The relative contribution of mutation and purifying selection to transition bias has not been quantitatively assessed in mitochondrial protein genes. The observed transition/transversion (s/v) ratio is (μ _s P _s)/(μ _v P _v), where μ _s and μ _v denote mutation rate of transitions and transversions, respectively, andP _s andP _v denote fixation probabilities of transitions and transversions, respectively. Because selection against synonymous transitions can be assumed to be roughly equal to that against synonymous transversions,P _s/P_v ≈ 1 at fourfold degenerate sites, so that thes/v ratio at fourfold degenerate sites is approximately μ _s /μ _v , which is a measure of mutational contribution to transition bias. Similarly, thes/v ratio at nondegenerate sites is also an estimate of μ _s /μ _v if we assume that selection against nonsynonymous transitions is roughly equal to that against nonsynonymous transversions. In two mitochondrial genes, cytochrome oxidase subunit I (COI) and cytochromeb (cyt-b) in pocket gophers, thes/v ratio is about two at nondegenerate and fourfold degenerate sites for both the COI and the cyt-b genes. This implies that mutation contribution to transition bias is relatively small. In contrast, thes/v ratio is much greater at twofold degenerate sites, being 48 for COI and 40 for cyt-b. Given that the μ _s /μ _v ratio is about 2, theP _s/P_v ratio at twofold degenerate sites must be on the order of 20 or greater. This suggests a great effect of purifying selection on transition bias in mitochondrial protein genes because transitions are synonymous and transversions are nonsynonymous at twofold degenerate sites in mammalian mitochondrial genes. We also found that nonsynonymous mutations at twofold degenerate sites are more neutral than nonsynonymous mutations at nondegenerate sites, and that the COI gene is subject to stronger purifying selection than is the cyt-b gene. A model is presented to integrate the effect of purifying selection, codon bias, DNA repair and GC content ons/v ratio of protein-coding genes. Correspondence to: X. Xia