Contrast-FEL—A Test for Differences in Selective Pressures at Individual Sites among Clades and Sets of Branches

A number of evolutionary hypotheses can be tested by comparing selective pressures among sets of branches in a phylogenetic tree. When the question of interest is to identify specific sites within genes that may be evolving differently, a common approach is to perform separate analyses on subsets of sequences and compare parameter estimates in a post hoc fashion. This approach is statistically suboptimal and not always applicable. Here, we develop a simple extension of a popular fixed effects likelihood method in the context of codon-based evolutionary phylogenetic maximum likelihood testing, Contrast-FEL. It is suitable for identifying individual alignment sites where any among the K≥2 sets of branches in a phylogenetic tree have detectably different ω ratios, indicative of different selective regimes. Using extensive simulations, we show that Contrast-FEL delivers good power, exceeding 90% for sufficiently large differences, while maintaining tight control over false positive rates, when the model is correctly specified. We conclude by applying Contrast-FEL to data from five previously published studies spanning a diverse range of organisms and focusing on different evolutionary questions.     

EPR properties of mixed-valent μ-oxo and μ-hydroxo dinuclear iron complexes produced by radiolytic reduction at 77 K

 Radiolytic reduction at 77 K of oxo-/hydroxo-bridged dinuclear iron(III) complexes in frozen solutions forms kinetically stabilized, mixed-valent species in high yields that model the mixed-valent sites of non-heme, diiron proteins. The mixed-valent species trapped at 77 K retain ligation geometry similar to the initial diferric clusters. The shapes of the mixed-valent EPR signals depend strongly on the bridging ligands. Spectra of the Fe(II)OFe(III) species reveal an S=1/2 ground state with small g-anisotropy as characterized by the uniaxial component (g _z –g _av /2<0.03) observable at temperatures as high as ∼100 K. In contrast, hydroxo-bridged mixed-valent species are characterized by large g-anisotropy (g _z –g _av /2>0.03) and are observable only below 30 K. Annealing at higher temperatures causes structural relaxation and changes in the EPR characteristics. EPR spectral properties allow the oxo- and hydroxo-bridged, mixed-valent diiron centers to be distinguished from each other and can help characterize the structure of mixed-valent centers in proteins. Received: 27 June 1998 / Accepted: 25 February 1999     

A unique,thermostable dimer linkage structure of RD114 retroviral RNA

Retroviruses package their genome as RNA dimers linked together primarily by base-pairing between palindromic stem–loop (psl) sequences at the 5′ end of genomic RNA. Retroviral RNA dimers usually melt in the range of 55°C–70°C. However, RNA dimers from virions of the feline endogenous gammaretrovirus RD114 were reported to melt only at 87°C. We here report that the high thermal stability of RD114 RNA dimers generated from in vitro synthesized RNA is an effect of multiple dimerization sites located in the 5′ region from the R region to sequences downstream from the splice donor (SD) site. By antisense oligonucleotide probing we were able to map at least five dimerization sites. Computational prediction revealed a possibility to form stems with autocomplementary loops for all of the mapped dimerization sites. Three of them were located upstream of the SD site. Mutant analysis supported a role of all five loop sequences in the formation and thermal stability of RNA dimers. Four of the five psls were also predicted in the RNA of two baboon endogenous retroviruses proposed to be ancestors of RD114. RNA fragments of the 5′ R region or prolonged further downstream could be efficiently dimerized in vitro. However, this was not the case for the 3′ R region linked to upstream U3 sequences, suggesting a specific mechanism of negative regulation of dimerization at the 3′ end of the genome, possibly explained by a long double-stranded RNA region at the U3-R border. Altogether, these data point to determinants of the high thermostability of the dimer linkage structure of the RD114 genome and reveal differences from other retroviruses.     

The lifetime of molecular (Davydov) solitons

The lifetime of Davydov solitons in a one-dimensional system is studied theoretically. The process of thermalization and the properties of solitons at finite temperature are investigated and the processes of soliton creation and disintegration are discussed.     

The Relationship between the Structure of the Tick-Borne Encephalitis Virus Strains and Their Pathogenic Properties

Tick-borne encephalitis virus (TBEV) is transmitted to vertebrates by taiga or forest ticks through bites, inducing disease of variable severity. The reasons underlying these differences in the severity of the disease are unknown. In order to identify genetic factors affecting the pathogenicity of virus strains, we have sequenced and compared the complete genomes of 34 Far-Eastern subtype (FE) TBEV strains isolated from patients with different disease severity (Primorye, the Russian Far East). We analyzed the complete genomes of 11 human pathogenic strains isolated from the brains of dead patients with the encephalitic form of the disease (Efd), 4 strains from the blood of patients with the febrile form of TBE (Ffd), and 19 strains from patients with the subclinical form of TBE (Sfd). On the phylogenetic tree, pathogenic Efd strains formed two clusters containing the prototype strains, Senzhang and Sofjin, respectively. Sfd strains formed a third separate cluster, including the Oshima strain. The strains that caused the febrile form of the disease did not form a separate cluster. In the viral proteins, we found 198 positions with at least one amino acid residue substitution, of which only 17 amino acid residue substitutions were correlated with the variable pathogenicity of these strains in humans and they authentically differed between the groups. We considered the role of each amino acid substitution and assumed that the deletion of 111 amino acids in the capsid protein in combination with the amino acid substitutions R16K and S45F in the NS3 protease may affect the budding process of viral particles. These changes may be the major reason for the diminished pathogenicity of TBEV strains. We recommend Sfd strains for testing as attenuation vaccine candidates.     

Changes in synapses after gamma-irradiation of the rat head

High reactivity and, at the same time, flexibility of interneuronal contacts were observed after exposure of rat head to 2-100 Gy radiation. At high doses (200-400 Gy) radiation-induced changes played a major role in the development of the cerebral form of radiation sickness. A complete asynapsis is probably one of the causes of the animals death "under the ray" (irradiation of the head with a dose of 1000 Gy).     

Volume-dependent regulation of ion transport and membrane phosphorylation in human and rat erythrocytes

Summary Osmotic swelling of human and rat erythrocytes does not induce regulatory volume decrease. Regulatory volume increase was observed in shrunken erythrocytes of rats only. This reaction was blocked by the inhibitors of Na⁺/H⁺ exchange. Cytoplasmic acidification in erythrocytes of both species increases the amiloride-inhibited component of²²Na influx by five- to eight-fold. Both the osmotic and isosmotic shrinkage of rat erythrocytes results in the 10- to 30-fold increase of amiloride-inhibited²²Na influx and a two-fold increase of furosemide-inhibited⁸⁶Rb influx. We failed to indicate any significant changes of these ion transport systems in shrunken human erythrocytes. The shrinking of quin 2-loaded human and rat erythrocytes results in the two- to threefold increase of the rate of⁴⁵Ca influx, which is completely blocked by amiloride. The dependence of volume-induced²²Na influx in rat erythrocytes and⁴⁵Ca influx in human erythrocytes on amiloride concentration does not differ. The rate of⁴⁵Ca influx in resealed ghosts was reduced by one order of magnitude when intravesicular potassium and sodium were replaced by choline. It is assumed that the erythrocyte shrinkage increases the rate of a nonselective Ca _o ²⁺ (Na _i ⁺ , K _i ⁺ ) exchange. Erythrocyte shrinking does not induce significant phosphorylation of membrane protein but increases the³²P incorporation in diphosphoinositides. The effect of shrinkage on the³²P labeling of phosphoinositides is diminished after addition of amiloride. It is assumed that volume-induced phosphoinositide response plays an essential role in the mechanism of the activation of transmembrane ion movements.     

Mathematical modelling of runoff and material transport from drainage areas into recipient water bodies

This paper describes a set of mathematical models to calculate runoff, soil erosion and chemical substance yields for drainage areas. The models represent different degrees of complexity and details of describing hydrological and hydrochemical processes. The availability of initial information and the requirements of the task determine the choice of a model for each particular case. These models have been written as PC/AT programs, to enable their use by a wide range of specialists. The models have been tested in different physical and geographical conditions, e.g. in some small watersheds on the Karelian Isthmus, west of Lake Ladoga, and in experimental catchment areas located in the Valdai Hills area in the southern part of Lake Ladoga drainage area. In further applications the drainage area models have been combined with models of the recipient water body, to obtain a single model of a watershed-lake system. The results of tests have proved that the models can be used for evaluation, simulation and forecasting of the processes of rainfall and snowmelt runoff, infiltration, evaporation, water erosion and pollutant wash-off from drainage areas.