Estimation of evolutionary distances under stationary and nonstationary models of nucleotide substitution

Estimation of evolutionary distances has always been a major issue in the study of molecular evolution because evolutionary distances are required for estimating the rate of evolution in a gene, the divergence dates between genes or organisms, and the relationships among genes or organisms. Other closely related issues are the estimation of the pattern of nucleotide substitution, the estimation of the degree of rate variation among sites in a DNA sequence, and statistical testing of the molecular clock hypothesis. Mathematical treatments of these problems are considerably simplified by the assumption of a stationary process in which the nucleotide compositions of the sequences under study have remained approximately constant over time, and there now exist fairly extensive studies of stationary models of nucleotide substitution, although some problems remain to be solved. Nonstationary models are much more complex, but significant progress has been recently made by the development of the paralinear and LogDet distances. This paper reviews recent studies on the above issues and reports results on correcting the estimation bias of evolutionary distances, the estimation of the pattern of nucleotide substitution, and the estimation of rate variation among the sites in a sequence.

Quantitative ER ↔ Golgi Transport Kinetics and Protein Separation upon Golgi Exit Revealed by Vesicular Integral Membrane Protein 36 Dynamics in Live CellsV⃞

To quantitatively investigate the trafficking of the transmembrane lectin VIP36 and its relation to cargo-containing transport carriers (TCs), we analyzed a C-terminal fluorescent-protein (FP) fusion, VIP36-SP-FP. When expressed at moderate levels, VIP36-SP-FP localized to the endoplasmic reticulum, Golgi apparatus, and intermediate transport structures, and colocalized with epitope-tagged VIP36. Temperature shift and pharmacological experiments indicated VIP36-SP-FP recycled in the early secretory pathway, exhibiting trafficking representative of a class of transmembrane cargo receptors, including the closely related lectin ERGIC53. VIP36-SP-FP trafficking structures comprised tubules and globular elements, which translocated in a saltatory manner. Simultaneous visualization of anterograde secretory cargo and VIP36-SP-FP indicated that the globular structures were pre-Golgi carriers, and that VIP36-SP-FP segregated from cargo within the Golgi and was not included in post-Golgi TCs. Organelle-specific bleach experiments directly measured the exchange of VIP36-SP-FP between the Golgi and endoplasmic reticulum (ER). Fitting a two-compartment model to the recovery data predicted first order rate constants of 1.22 ± 0.44%/min for ER → Golgi, and 7.68 ± 1.94%/min for Golgi → ER transport, revealing a half-time of 113 ± 70 min for leaving the ER and 1.67 ± 0.45 min for leaving the Golgi, and accounting for the measured steady-state distribution of VIP36-SP-FP (13% Golgi/87% ER). Perturbing transport with AlF4− treatment altered VIP36-SP-GFP distribution and changed the rate constants. The parameters of the model suggest that relatively small differences in the first order rate constants, perhaps manifested in subtle differences in the tendency to enter distinct TCs, result in large differences in the steady-state localization of secretory components.

Silver ion high pressure liquid chromatography provides unprecedented separation of sterols: application to the enzymatic formation of cholesta-5,8-dien-3 beta-ol.

We report that silver ion HPLC provides remarkable separations of C27 sterols differing only in the number or location of olefinic double bonds. This technique has been extended to LC-MS, analysis of purified components by GC, GC-MS, and 1H NMR, and to its use on a semipreparative scale. The application of this methodology for the demonstration of the catalysis, by rat liver microsomes, of the conversion of 7-dehydrocholesterol to cholesta-5,8-dien-3 beta-ol is also presented.

Highly plastic chromosomal organization in Salmonella typhi.

Gene order in the chromosomes of Escherichia coli K-12 and Salmonella typhimurium LT2, and in many other species of Salmonella, is strongly conserved, even though the genera diverged about 160 million years ago. However, partial digestion of chromosomal DNA of Salmonella typhi, the causal organism of typhoid fever, with the endonuclease I-CeuI followed by separation of the DNA fragments by pulsed-field gel electrophoresis showed that the chromosomes of independent wild-type isolates of S. typhi are rearranged due to homologous recombination between the seven rrn genes that code for ribosomal RNA. The order of genes within the I-CeuI fragments is largely conserved, but the order of the fragments on the chromosome is rearranged. Twenty-one different orders of the I-CeuI fragments were detected among the 127 wild-type strains we examined. Duplications and deletions were not found, but transpositions and inversions were common. Transpositions of I-CeuI fragments into sites that do not change their distance from the origin of replication (oriC) are frequently detected among the wild-type strains, but transpositions that move the fragments much further from oriC were rare. This supports the gene dosage hypothesis that genes at different distances from oriC have different gene dosages and, hence, different gene expression, and that during evolution genes become adapted to their specific location; thus, cells with changes in gene location due to transpositions may be less fit. Therefore, gene dosage may be one of the forces that conserves gene order, although its effects seem less strong in S. typhi than in other enteric bacteria. However, both the gene dosage and the genomic balance hypotheses, the latter of which states that the origin (oriC) and terminus (TER) of replication must be separated by 180 degrees C, need further investigation.

Genetic separation of tumor growth and hemorrhagic phenotypes in an estrogen-induced tumor.

Chronic administration of estrogen to the Fischer 344 (F344) rat induces growth of large, hemorrhagic pituitary tumors. Ten weeks of diethylstilbestrol (DES) treatment caused female F344 rat pituitaries to grow to an average of 109.2 +/- 6.3 mg (mean +/- SE) versus 11.3 +/- 1.4 mg for untreated rats, and to become highly hemorrhagic. The same DES treatment produced no significant growth (8.9 +/- 0.5 mg for treated females versus 8.7 +/- 1.1 for untreated females) or morphological changes in Brown Norway (BN) rat pituitaries. An F1 hybrid of F344 and BN exhibited significant pituitary growth after 10 weeks of DES treatment with an average mass of 26.3 +/- 0.7 mg compared with 8.6 +/- 0.9 mg for untreated rats. Surprisingly, the F1 hybrid tumors were not hemorrhagic and had hemoglobin content and outward appearance identical to that of BN. Expression of both growth and morphological changes is due to multiple genes. However, while DES-induced pituitary growth exhibited quantitative, additive inheritance, the hemorrhagic phenotype exhibited recessive, epistatic inheritance. Only 5 of the 160 F2 pituitaries exhibited the hemorrhagic phenotype; 36 of the 160 F2 pituitaries were in the F344 range of mass, but 31 of these were not hemorrhagic, indicating that the hemorrhagic phenotype is not merely a consequence of extensive growth. The hemorrhagic F2 pituitaries were all among the most massive, indicating that some of the genes regulate both phenotypes.

Separation of the transcriptional coactivator and antirepression functions of transcription factor IIA.

Human transcription factor IIA (TFIIA) is composed of three subunits (alpha, beta, and gamma). TFIIA interacts with the TATA-box binding protein and can overcome repression of transcription. TFIIA was found to be necessary for VP16-mediated transcriptional activation through a coactivator function. We have separated the coactivator and antirepression activities of TFIIA. A TFIIA lacking the alpha subunit was isolated from HeLa cells. This "mini-TFIIA" interacts with the TATA-box binding protein and can overcome repression of transcription, but it is defective in transcriptional coactivator function.

Long-time quantum simulation of the primary charge separation in bacterial photosynthesis.

Accurate quantum mechanical simulations of the primary charge transfer in photosynthetic reaction centers are reported. The process is modeled by three coupled electronic states corresponding to the photoexcited chlorophyll special pair (donor), the reduced bacteriopheophytin (acceptor), and the reduced accessory chlorophyll (bridge) that interact with a dissipative medium of protein and solvent degrees of freedom. The time evolution of the excited special pair is followed over 17 ps by using a fully quantum mechanical path integral scheme. We find that a free energy of the reduced accessory chlorophyll state approximately equal to 400 cm(-1) lower than that of the excited special pair state yields state populations in agreement with experimental results on wild-type and modified reaction centers. For this energetic configuration electron transfer is a two-step process.

Phase separation in aqueous solutions of lens gamma-crystallins: special role of gamma s.

We have studied liquid-liquid phase separation in aqueous ternary solutions of calf lens gamma-crystallin proteins. Specifically, we have examined two ternary systems containing gamma s--namely, gamma IVa with gamma s in water and gamma II with gamma s in water. For each system, the phase-separation temperatures (Tph (phi)) alpha as a function of the overall protein volume fraction phi at various fixed compositions alpha (the "cloud-point curves") were measured. For the gamma IVa, gamma s, and water ternary solution, a binodal curve composed of pairs of coexisting points, (phi I, alpha 1) and (phi II, alpha II), at a fixed temperature (20 degrees C) was also determined. We observe that on the cloud-point curve the critical point is at a higher volume fraction than the maximum phase-separation temperature point. We also find that typically the difference in composition between the coexisting phases is at least as significant as the difference in volume fraction. We show that the asymmetric shape of the cloud-point curve is a consequence of this significant composition difference. Our observation that the phase-separation temperature of the mixtures in the high volume fraction region is strongly suppressed suggests that gamma s-crystallin may play an important role in maintaining the transparency of the lens.

Microsatellite variability and genetic distances.

We analyze the within- and between-population dynamics of the distribution of the number of repeats at multiple microsatellite DNA loci subject to stepwise mutation. Analytical expressions for moments up to the fourth order within a locus and the variance of between-locus variance at mutation-drift equilibrium have been obtained. These statistics may be used to test the appropriateness of the one-step mutation model and to detect between-locus variation in the mutation rate. Published data are compatible with the one-step mutation model, although they do not reject the two-step model. Using both multinomial sampling and diffusion approximations for the analysis of the genetic distance introduced by Goldstein et al. [Goldstein, D. B., Linares, A. R., Cavalli-Sforza, L. L. & Feldman, M. W. (1995) Proc. Natl. Acad. Sci. USA 92, 6723-6727], we show that this distance follows a chi 2 distribution with degrees of freedom equal to the number of loci when there is no variation in mutation rates among the loci. In the presence of such variation, the variance of the distance is obtained. We conclude that the number of microsatellite loci required for the construction of phylogenetic trees with reliable branch lengths may be several hundred. Also, mutations that change repeat scores by several units, even though extremely rare, may dramatically influence estimates of population parameters.

Kinetics of spindle pole body separation in budding yeast.

In the budding yeast Saccharomyces cerevisiae, the spindle pole body (SPB) serves as the microtubule-organizing center and is the functional analog of the centrosome of higher organisms. By expressing a fusion of a yeast SPB-associated protein to the Aequorea victoria green fluorescent protein, the movement of the SPBs in living yeast cells undergoing mitosis was observed by fluorescence microscopy. The ability to visualize SPBs in vivo has revealed previously unidentified mitotic events. During anaphase, the mitotic spindle has four sequential activities: alignment at the mother-daughter junction, fast elongation, translocation into the bud, and slow elongation. These results indicate that distinct forces act upon the spindle at different times during anaphase.

Identification of receptor ligands and receptor subtypes using antagonists in a capillary electrophoresis single-cell biosensor separation system.

A capillary electrophoresis system with single-cell biosensors as a detector has been used to separate and identify ligands in complex biological samples. The power of this procedure was significantly increased by introducing antagonists that inhibited the cellular response from selected ligand-receptor interactions. The single-cell biosensor was based on the ligand-receptor binding and G-protein-mediated signal transduction pathways in PC12 and NG108-15 cell lines. Receptor activation was measured as increases in cytosolic free calcium ion concentration by using fluorescence microscopy with the intracellular calcium ion indicator fluo-3-acetoxymethyl ester. Specifically, a mixture of bradykinin (BK) and acetylcholine (ACh) was fractionated and the components were identified by inhibiting the cellular response with icatibant (HOE 140), a selective antagonist to the BK B2 receptor subtype (B2BK), and atropine, an antagonist to muscarinic ACh receptor subtypes. Structurally related forms of BK were also identified based on inhibiting B2BK receptors. Applications of this technique include identification of endogenous BK in a lysate of human hepatocellular carcinoma cells (Hep G2) and screening for bioactivity of BK degradation products in human blood plasma. The data demonstrate that the use of antagonists with a single-cell biosensor separation system aids identification of separated components and receptor subtypes.