THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

Signal relay between two integral membrane proteins: The sensory rhodopsin I/HtrI transducer molecular complex

The molecular complex of sensory rhodopsin I (SRI) and its transducer HtrI mediate color-sensitive phototaxis in the archaeon Halobacterium salinarum. Orange light causes an attractant response by a one-photon reaction and white light causes a repellent response by a two-photon reaction. Three aspects of this molecular complex were explored: (i) We determined the stoichiometry of SRI and HtrI to be 2:2 by gene fusion analysis. A SRI-HtrI fusion protein was expressed in H. salinarum and shown to mediate 1-photon and 2-photon phototaxis responses comparable to wild-type complex. Disulfide crosslinking demonstrated that the fusion protein is a homodimer in the membrane. Measurement of photochemical reaction kinetics and pH titration of absorption spectra established that both SRI domains are complexed to HtrI in the fusion protein, and therefore the stoichiometry is 2:2. (ii) Cytoplasmic channel closure of SRI by HtrI, an important aspect of their interaction, was investigated by incremental HtrI truncation. We found that binding of the membrane-embedded portion of HtrI is insufficient for channel closure, whereas cytoplasmic extension of the second HtrI transmembrane helix by 13 residues blocks proton conduction through the channel as well as full-length HtrI. The closure activity is localized to 5 specific residues, each of which incrementally contributes to reduction of proton conductivity. Moreover, these same residues in the dark incrementally and proportionally increase the pKa of the Asp76 counterion to the protonated Schiff base chromophore. We conclude that this critical region of HtrI alters the dark conformation of SRI as well as light-induced channel opening. (iii) We developed a procedure for reconstituting HtrI-free SRI and the SRI/HtrI complex into liposomes, which exhibit photocycles with opened and closed cytoplasmic channels, respectively, as in the membrane. This opens the way for study of the light-induced conformational change and the interaction in vitro by fluorescence and spin-labeling. Single-cysteine mutations were introduced into helix F of SRI, labeled with a nitroxide spin probe and a fluorescence probe, reconstituted into proteoliposomes, and light-induced conformational changes detected in the complex. The probe signals can now be used as the readout of signaling to analyze mutants and the kinetics of signal relay. ^

Molecular mechanisms of response in gastrointestinal stromal tumors treated with imatinib mesylate

Gastrointestinal stromal tumors (GIST) represent 80% of sarcoma arising from the GI tract. The inciting event in tumor progression is mutation of the kit or, rarely, platelet derived growth factor receptor-α (PDGFR) gene. These mutations encode ligand independent, constitutively active proteins: Kit or PDGFR. ^ These tumors are notoriously chemo and radio resistant. Historically, patients with advanced disease realized a median overall survival of 9 months. However, with modern management of GIST with imatinib mesylate (Novartis), a small molecule inhibitor of the Kit, PDGFR, and Abl tyrosine kinases, patients now realize a median overall survival greater than 30 months. However, almost half of patients present with surgically resectable GIST and the utility of imatinib in this context has not been prospectively studied. Also, therapeutic benefit of imatinib is variable from patient to patient and alternative targeted therapy is emerging as potential alternatives to imatinib. Thus, elucidating prognostic factors for patients with GIST in the imatinib-era is crucial to providing optimal care to each particular patient. Moreover, the exact mechanism of action of imatinib in GIST is not fully understood. Therefore, physicians find difficulty in accurately predicting which patient will benefit from imatinib, how to assess response to therapy, and the time at which to assess response. ^ I have hypothesized that imatinib is tolerable and clinically beneficial in the context of surgery, VEGF expression and kit non-exon 11 genotypes portend poor survival on imatinib therapy, and imatinib's mechanism of action is in part due to anti-vascular effects and inhibition of the Kit/SCF signaling axis of tumor-associated endothelial cells. ^ Results herein demonstrate that imatinib is safe and increases the duration of disease-free survival when combined with surgery. Radiographic and molecular (namely, apoptosis) changes occur within 3 days of imatinib initiation. I illustrate that non-exon 11 mutant genotypes and VEGF are poor prognostic factors for patients treated with imatinib. These findings may allow for patient stratification to emerging therapies rather than imatinib. I show that imatinib has anti-vascular effects via inducing tumor endothelial cell apoptosis perhaps by abrogation of the Kit/SCF signaling axis. ^

Molecular epidemiological studies on genetic predispositions to squamous cell carcinoma of the head and neck

Two molecular epidemiological studies were conducted to examine associations between genetic variation and risk of squamous cell carcinoma of the head and neck (SCCHN). In the first study, we hypothesized that genetic variation in p53 response elements (REs) may play roles in the etiology of SCCHN. We selected and genotyped five polymorphic p53 REs as well as a most frequently studied p53 codon 72 (Arg72Pro, rs1042522) polymorphism in 1,100 non-Hispanic White SCCHN patients and 1,122 age-and sex-matched cancer-free controls recruited at The University of Texas M. D. Anderson Cancer Center. In multivariate logistic regression analysis with adjustment for age, sex, smoking and drinking status, marital status and education level, we observed that the EOMES rs3806624 CC genotype had a significant effect of protection against SCCHN risk (adjusted odds ratio= 0.79, 95% confidence interval =0.64–0.98), compared with the -838TT+CT genotypes. Moreover, a significantly increased risk associated with the combined genotypes of p53 codon 72CC and EOMES -838TT+CT was observed, especially in the subgroup of non-oropharyneal cancer patients. The values of false-positive report probability were also calculated for significant findings. In the second study, we assessed the association between SCCHN risk and four potential regulatory single nucleotide polymorphisms (SNPs) of DEC1 (deleted in esophageal cancer 1) gene, a candidate tumor suppressor gene for esophageal cancer. After adjustment for age, sex, and smoking and drinking status, the variant -606CC (i.e., -249CC) homozygotes had a significantly reduced SCCHN risk (adjusted odds ratio = 0.71, 95% confidence interval = 0.52–0.99), compared with the -606TT homozygotes. Stratification analyses showed that a reduced risk associated with the -606CC genotype was more pronounced in subgroups of non-smokers, non-drinkers, younger subjects (defined as ≤ 57 years), carriers of TP53 Arg/Arg (rs1042522) genotype, patients with oropharyngeal cancer or late-stage SCCHN. Further in silico analysis revealed that the -249 T-to-C change led to a gain of a transcription factor binding site. Additional functional analysis showed that the -249T-to-C change significantly enhanced transcriptional activity of the DEC1 promoter and the DNA-protein binding activity. We conclude that the DEC1 promoter -249 T>C (rs2012775) polymorphism is functional, modulating susceptibility to SCCHN among non-Hispanic Whites. Additional large-scale, preferably population-based studies are needed to validate our findings.^

Molecular and evolutionary genetics of the X -linked visual pigment genes in humans and New World monkeys

Normal humans have one red and at least one green visual pigment genes. These genes are tightly linked as tandem repeats on the X chromosome and each of them has six exons. There is only one X-linked visual pigment gene in New World monkeys (NWMs) but the locus has three polymorphic alleles encoding red, yellow and green visual pigments, respectively. The spectral properties of the squirrel monkey and the marmoset (both NWMs) have been studied and partial sequences of the three alleles are available. To study the evolutionary history of these X-linked opsin genes in humans and NWMs, coding and intron sequences of the three squirrel monkey alleles and the three marmoset alleles were amplified by PCR followed by subcloning and sequencing. Introns 2 and 4 of the human red and green pigment genes were also sequenced. The results obtained are as follows: (1) The sequences of introns 2 and 4 of the human red and green opsin genes are significantly more similar between the two genes than are coding sequences, contrary to the usual situation where coding regions are better conserved in evolution than are introns. The high similarities in the two introns are probably due to recent gene conversion events during evolution of the human lineage. (2) Phylogenetic analysis of both intron and exon sequences indicates that the phylogenetic tree of the available primate opsin genes is the same as the species tree. The two human genes were derived from a gene duplication event after the divergence of the human and NWM lineages. The three alleles in each of the two NWM species diverged after the split of the two NWMs but have persisted in the population for at least 5 million years. (3) Allelic gene conversion might have occurred between the three squirrel monkey alleles. (4) A model of additive effect of hydroxyl-bearing amino acids on spectral tuning is proposed by treating some unknown variables as groups. Under the assumption that some residues have no effect, it is found that at least five amino acid residues, at positions 178 (3 nm), 180 (5 nm), 230 ($-$4 nm), 277 (9 nm) and 285 (13 nm), have linear spectral tuning effects. (5) Adaptive evolution of the opsin genes to different spectral peaks was observed at four residues that are important for spectral tuning. ^

Mechanism of molecular regulation of nuclear -encoded mitochondrial gene expression by electrical stimulation in the cultured neonatal rat cardiac myocytes

To answer the question whether increased energy demand resulting from myocyte hypertrophy and enhanced $\beta$-myosin heavy chain mRNA, contractile protein synthesis and assembly leads to mitochondrial proliferation and differentiation, we set up an electrical stimulation model of cultured neonatal rat cardiac myocytes. We describe, as a result of increased contractile activity, increased mitochondrial profiles, cytochrome oxidase mRNA, and activity, as well as a switch in mitochondrial carnitine palmitoyltransferase-I (CPT-I) from the liver to muscle isoform. We investigate physiological pathways that lead to accumulation of gene transcripts for nuclear encoded mitochondrial proteins in the heart. Cardiomyocytes were stimulated for varying times up to 72 hr in serum-free culture. The mRNA contents for genes associated with transcriptional activation (c-fos, c-jun, junB, nuclear respiratory factor 1 (Nrf-1)), mitochondrial proliferation (cytochrome c (Cyt c), cytochrome oxidase), and mitochondrial differentiation (carnitine palmitonyltransferase I (CPT-I) isoforms) were measured. The results establish a temporal pattern of mRNA induction beginning with c-fos (0.25-3 hr) and followed by c-jun (0.5-3 hr), junB (0.5-6 hr), NRF-1 (1-12 hr), Cyt c (12-72 hr), cytochrome c oxidase (12-72 hr). Induction of the latter was accompanied by a marked decrease in the liver-specific CPT-I mRNA. Electrical stimulation increased c-fos, $\beta$-myosin heavy chain, and Cyt c promoter activities. These increases coincided with a rise in their respective endogenous gene transcripts. NRF-1, cAMP response element (CRE), and Sp-1 site mutations within the Cyt c promoter reduced luciferase expression in both stimulated and nonstimulated myocytes. Mutations in the Nrf-1 and CRE sites inhibited the induction by electrical stimulation or by transfection of c-jun into non-paced cardiac myocytes whereas mutation of the Sp-1 site maintained or increased the fold induction. This is consistent with the appearance of NRF-1 and fos/jun mRNAs prior to that of Cyt c. Overexpression of c-jun by transfection also activates the Nrf-1 and Cyt c mRNA sequentially. Electrical stimulation of cardiac myocytes activates the c-Jun-N-terminal kinase so that the fold-activation of the cyt c promoter is increased by pacing when either c-jun or c-fos/c-jun are cotransfected. We have identified physical association of Nrf-1 protein with the Nrf-1 enhancer element and of c-Jun with the CRE binding sites on the Cyt c promoter. This is the first demonstration that induction of Nrf-1 and c-Jun by pacing of cardiac myocytes directly mediates Cyt c gene expression and mitochondrial proliferation in response to hypertrophic stimuli in the heart.^ Subsequent to gene activation pathways that lead to mitochondrial proliferation, we observed an isoform switch in CPT-I from the liver to muscle mRNA. We have found that the half-life for the muscle CPT-I is not affected by electrical stimulation, but electrical decrease the T1/2 in the liver CPT-I by greater than 50%. This suggests that the liver CPT-I switch to muscle isoform is due to (1) a decrease in T1/2 of liver CPT-I and (2) activation of muscle CPT-Itranscripts by electrical stimulation. (Abstract shortened by UMI.) ^