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.     

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 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.     

Evolutionary changes in the structural and regulatory regions of aminoglycoside-3'-phosphotransferase type I gene of E. coli

To investigate the evolutionary relationships between the aph(3') genes from different plasmids, the nucleotide sequence of the aph(3') gene from the E. coli R plasmid was determined and compared with the known aph(3') genes of Tn903 and Tn4352. Three point mutations in the structural part of the cloned aph(3') gene caused amino acid changes in the enzyme molecule at positions 19, 27 and 48 beginning from the start codon. The structural part of the gene was followed by two stop codons and a long DNA region containing no nucleotide sequences homologous to the sequences of Tn903 or Tn4352. Both the cloned aph(3') gene and Tn4352 were limited on the left by the spacer sequence and the insertion sequence IS176. Twenty one base pairs deletion abolished the -35 sequence of the promoter suggested for the aph(3') gene of Tn4352 and resulted in formation of a fusion promoter utilizing the -35 box of IS176 and the -10 box of the aph(3') gene. The distance between the -35 and -10 sequences changed from 18 to 17 bp. Changes in the cloned aph(3') gene and the flanking DNA regions resulted in formation of a new promoter and loss of the right IS176 element.     

Genes of resistance to aminoglycoside antibiotics in clinical strains of Serratia marcescens and the enzyme encoded by them]

The patterns of aminoglycoside inactivating enzymes were determined by AGRP in 31 clinical isolated of Serratia marcescens. The results were compared with the data on identification of the aminoglycoside resistance genes by the specific DNA probes. It was shown that all the isolates of Serratia marcescens contained the AAC(6')-Ic gene which was not expressed in some isolates. The other detected aminoglycoside inactivating enzymes were the following: AAC(3)-V in 17 isolates, ANT(2') in 7 isolates, AAC(3)-I in 4 isolates and APH(3')-I in 13 isolates. Reliability of the methods of AGRP and DNA-DNA hybridization was estimated in the assay of the aminoglycoside resistant clinical strains of Serratia marcescens.     

Aminoglycoside-inactivating enzymes in strains of Salmonella and genes encoding them

Aminoglycoside resistance patterns of 56 strains isolated from man, cattle and environment were determined. 34 out of 42 gentamicin-resistant strains were shown to produce AAC(3)-II and 7 strains produced ANT(2"). All the 48 kanamycin resistant strains produced APH(3')-I. Spot hybridization of the 42 gentamicin resistant strains with the inner fragment of the aacC2 gene revealed positive signals for all the strains. Hybridization of the 48 kanamycin-resistant strains with the aphA1 gene probe provided positive results in all the strains. The AAC(3)-IV encoding gene was not detected by DNA-DNA hybridization in the strains studied.     

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.