Function-altering SNPs in the human multidrug transporter gene ABCB1 identified using a Saccharomyces-based assay.

The human ABCB1 (MDR1)-encoded multidrug transporter P-glycoprotein (P-gp) plays a major role in disposition and efficacy of a broad range of drugs including anticancer agents. ABCB1 polymorphisms could therefore determine interindividual variability in resistance to these drugs. To test this hypothesis we developed a Saccharomyces-based assay for evaluating the functional significance of ABCB1 polymorphisms. The P-gp reference and nine variants carrying amino-acid-altering single nucleotide polymorphisms (SNPs) were tested on medium containing daunorubicin, doxorubicin, valinomycin, or actinomycin D, revealing SNPs that increased (M89T, L662R, R669C, and S1141T) or decreased (W1108R) drug resistance. The R669C allele's highly elevated resistance was compromised when in combination with W1108R. Protein level or subcellular location of each variant did not account for the observed phenotypes. The relative resistance profile of the variants differed with drug substrates. This study established a robust new methodology for identification of function-altering polymorphisms in human multidrug transporter genes, identified polymorphisms affecting P-gp function, and provided a step toward genotype-determined dosing of chemotherapeutics.

Genome-wide identification of the regulatory targets of a transcription factor using biochemical characterization and computational genomic analysis.

BACKGROUND: A major challenge in computational genomics is the development of methodologies that allow accurate genome-wide prediction of the regulatory targets of a transcription factor. We present a method for target identification that combines experimental characterization of binding requirements with computational genomic analysis. RESULTS: Our method identified potential target genes of the transcription factor Ndt80, a key transcriptional regulator involved in yeast sporulation, using the combined information of binding affinity, positional distribution, and conservation of the binding sites across multiple species. We have also developed a mathematical approach to compute the false positive rate and the total number of targets in the genome based on the multiple selection criteria. CONCLUSION: We have shown that combining biochemical characterization and computational genomic analysis leads to accurate identification of the genome-wide targets of a transcription factor. The method can be extended to other transcription factors and can complement other genomic approaches to transcriptional regulation.

The Smk1p MAP kinase negatively regulates Gsc2p, a 1,3-beta-glucan synthase, during spore wall morphogenesis in Saccharomyces cerevisiae.

Spore formation in Saccharomyces cerevisiae involves the sequential deposition of multiple spore wall layers between the prospore membranes that surround each meiotic product. The Smk1p mitogen-activated protein (MAP) kinase plays a critical role in spore formation, but the proteins that interact with Smk1p to regulate spore morphogenesis have not been described. Using mass spectrometry, we identify Gsc2p as a Smk1p-associated protein. Gsc2p is a 1,3-beta-glucan synthase subunit involved in synthesizing an inner spore wall layer. We find that 1,3-beta-glucan synthase activity is elevated in smk1 mutants, suggesting that SMK1 negatively regulates GSC2. Although deposition of the two inner spore wall layers is normal in smk1 mutants, deposition of the outer layers is aberrant. However, eliminating GSC2 activity restores normal deposition of the third spore wall layer in smk1 mutants, indicating that negative regulation of GSC2 by SMK1 is important for spore wall deposition. Our findings suggest a model for the coordination of spore wall layer deposition in which Smk1p facilitates the transition between early and late phases of spore wall deposition by inhibiting a spore wall-synthesizing enzyme important for early phases of spore wall deposition.

Cell polarity protein Spa2P associates with proteins involved in actin function in Saccharomyces cerevisiae.

Spa2p is a nonessential protein that regulates yeast cell polarity. It localizes early to the presumptive bud site and remains at sites of growth throughout the cell cycle. To understand how Spa2p localization is regulated and to gain insight into its molecular function in cell polarity, we used a coimmunoprecipitation strategy followed by tandem mass spectrometry analysis to identify proteins that associate with Spa2p in vivo. We identified Myo1p, Myo2p, Pan1p, and the protein encoded by YFR016c as proteins that interact with Spa2p. Strikingly, all of these proteins are involved in cell polarity and/or actin function. Here we focus on the functional significance of the interactions of Spa2p with Myo2p and Myo1p. We find that localization of Spa2GFP to sites of polarized growth depends on functional Myo2p but not on Myo1p. We also find that Spa2p, like Myo2p, cosediments with F-actin in an ATP-sensitive manner. We hypothesize that Spa2p associates with actin via a direct or indirect interaction with Myo2p and that Spa2p may be involved in mediating polarized localization of polarity proteins via Myo2p. In addition, we observe an enhanced cell-separation defect in a myo1spa2 strain at 37 degrees C. This provides further evidence that Spa2p is involved in cytokinesis and cell wall morphogenesis.

Genomic dissection of the cell-type-specification circuit in Saccharomyces cerevisiae.

The budding yeast Saccharomyces cerevisiae has three cell types (a cells, alpha cells, and a/alpha cells), each of which is specified by a unique combination of transcriptional regulators. This transcriptional circuit has served as an important model for understanding basic features of the combinatorial control of transcription and the specification of cell type. Here, using genome-wide chromatin immunoprecipitation, transcriptional profiling, and phylogenetic comparisons, we describe the complete cell-type-specification circuit for S. cerevisiae. We believe this work represents a complete description of cell-type specification in a eukaryote.

Anatomy and dynamics of DNA replication fork movement in yeast telomeric regions.

Replication initiation and replication fork movement in the subtelomeric and telomeric DNA of native Y' telomeres of yeast were analyzed using two-dimensional gel electrophoresis techniques. Replication origins (ARSs) at internal Y' elements were found to fire in early-mid-S phase, while ARSs at the terminal Y' elements were confirmed to fire late. An unfired Y' ARS, an inserted foreign (bacterial) sequence, and, as previously reported, telomeric DNA each were shown to impose a replication fork pause, and pausing is relieved by the Rrm3p helicase. The pause at telomeric sequence TG(1-3) repeats was stronger at the terminal tract than at the internal TG(1-3) sequences located between tandem Y' elements. We show that the telomeric replication fork pause associated with the terminal TG(1-3) tracts begins approximately 100 bp upstream of the telomeric repeat tract sequence. Telomeric pause strength was dependent upon telomere length per se and did not require the presence of a variety of factors implicated in telomere metabolism and/or known to cause telomere shortening. The telomeric replication fork pause was specific to yeast telomeric sequence and was independent of the Sir and Rif proteins, major known components of yeast telomeric heterochromatin.

Unique and redundant roles for HOG MAPK pathway components as revealed by whole-genome expression analysis.

The Saccharomyces cerevisiae high osmolarity glycerol (HOG) mitogen-activated protein kinase pathway is required for osmoadaptation and contains two branches that activate a mitogen-activated protein kinase (Hog1) via a mitogen-activated protein kinase kinase (Pbs2). We have characterized the roles of common pathway components (Hog1 and Pbs2) and components in the two upstream branches (Ste11, Sho1, and Ssk1) in response to elevated osmolarity by using whole-genome expression profiling. Several new features of the HOG pathway were revealed. First, Hog1 functions during gene induction and repression, cross talk inhibition, and in governing the regulatory period. Second, the phenotypes of pbs2 and hog1 mutants are identical, indicating that the sole role of Pbs2 is to activate Hog1. Third, the existence of genes whose induction is dependent on Hog1 and Pbs2 but not on Ste11 and Ssk1 suggests that there are additional inputs into Pbs2 under our inducing conditions. Fourth, the two upstream pathway branches are not redundant: the Sln1-Ssk1 branch has a much more prominent role than the Sho1-Ste11 branch for activation of Pbs2 by modest osmolarity. Finally, the general stress response pathway and both branches of the HOG pathway all function at high osmolarity. These studies demonstrate that cells respond to increased osmolarity by using different signal transduction machinery under different conditions.

The Cdk-activating kinase Cak1p promotes meiotic S phase through Ime2p.

CAK1 encodes an essential protein kinase in Saccharomyces cerevisiae that is required for activation of the Cdc28p Cdk. CAK1 also has several CDC28-independent functions that are unique to meiosis. The earliest of these functions is to induce S phase, which is regulated differently in meiosis than in mitosis. In mitosis, Cdc28p controls its own S-phase-promoting activity by signaling the destruction of its inhibitor, Sic1p. In meiosis, Sic1p destruction is signaled by the meiosis-specific Ime2p protein kinase. Our data show that Cak1p is required to activate Ime2p through a mechanism that requires threonine 242 and tyrosine 244 in Ime2p's activation loop. This activation promotes autophosphorylation and accumulation of multiply phosphorylated forms of Ime2p during meiotic development. Consistent with Cak1p's role in activating Ime2p, cells lacking Cak1p are deficient in degrading Sic1p. Deletion of SIC1 or overexpression of IME2 can partially suppress the S-phase defect in cak1 mutant cells, suggesting that Ime2p is a key target of Cak1p regulation. These data show that Cak1p is required for the destruction of Sic1p in meiosis, as in mitosis, but in meiosis, it functions through a sporulation-specific kinase.

Sequence diversity and haplotype structure in the human ABCB1 (MDR1, multidrug resistance transporter) gene.

OBJECTIVES: There is increasing evidence that polymorphism of the ABCB1 (MDR1) gene contributes to interindividual variability in bioavailability and tissue distribution of P-glycoprotein substrates. The aim of the present study was to (1) identify and describe novel variants in the ABCB1 gene, (2) understand the extent of variation in ABCB1 at the population level, (3) analyze how variation in ABCB1 is structured in haplotypes, and (4) functionally characterize the effect of the most common amino acid change in P-glycoprotein. METHODS AND RESULTS: Forty-eight variant sites, including 30 novel variants and 13 coding for amino acid changes, were identified in a collection of 247 ethnically diverse DNA samples. These variants comprised 64 statistically inferred haplotypes, 33 of which accounted for 92% of chromosomes analyzed. The two most common haplotypes, ABCB1*1 and ABCB1*13, differed at six sites (three intronic, two synonymous, and one non-synonymous) and were present in 36% of all chromosomes. Significant population substructure was detected at both the nucleotide and haplotype level. Linkage disequilibrium was significant across the entire ABCB1 gene, especially between the variant sites found in ABCB1*13, and recombination was inferred. The Ala893Ser change found in the common ABCB1*13 haplotype did not affect P-glycoprotein function. CONCLUSION: This study represents a comprehensive analysis of ABCB1 nucleotide diversity and haplotype structure in different populations and illustrates the importance of haplotype considerations in characterizing the functional consequences of ABCB1 polymorphisms.

Control of landmark events in meiosis by the CDK Cdc28 and the meiosis-specific kinase Ime2.

Meiosis is thought to require the protein kinase Ime2 early for DNA replication and the cyclin-dependent kinase Cdc28 late for chromosome segregation. To elucidate the roles of these kinases, we inhibited their activities early and late using conditional mutants that are sensitive to chemical inhibitors. Our studies reveal that both Cdc28 and Ime2 have critical roles in meiotic S phase and M phase. Early inhibition of analog-sensitive cdc28-as1 blocked DNA replication, revealing a previously undetected role for Cdc28. Yet Cdc28 was dispensable for one of its functions in the mitotic cell cycle, degradation of Sic1. Late addition of inhibitor to ime2-as1 revealed unexpected roles of Ime2 in the initiation and execution of chromosome segregation. The requirement of Ime2 for M phase is partially explained by its stimulation of the key meiotic transcription factor Ndt80, which is needed in turn for high Cdc28 activity. In accordance with a late role for Ime2, we observed an increase in its activity during M phase that depended on Cdc28 and Ndt80. We speculate that several unique features of the meiotic cell division reflect a division of labor and regulatory coordination between Ime2 and Cdc28.

Functional characterization in yeast of genetic variants in the human equilibrative nucleoside transporter, ENT1.

The human equilibrative nucleoside transporter, ENT1, appears to play a critical role in the disposition of nucleosides and nucleoside analogs used clinically as anti-viral and anti-cancer drugs. Recently, we identified variants of ENT1 in an ethnically diverse DNA sample set from 247 individuals, focusing primarily on the coding region. The goal of the present study was to analyse the haplotype structure and functionally characterize the variants of ENT1. We observed that a single haplotype, ENT1*1, accounted for 91.3% of the 494 chromosomes. Functional analysis in Saccharomyces cerevisiae revealed no differences in the kinetics of uptake of nucleosides and nucleoside analogs by the two non-synonymous variant transporters, ENT1-I216T and ENT1-E391K, and the reference ENT1. These results, together with the observation that there are few haplotypes of ENT1, indicate that coding region variants of ENT1 do not contribute to inter-individual differences in response to nucleoside analog drugs.

Natural variation in human membrane transporter genes reveals evolutionary and functional constraints.

Membrane transporters maintain cellular and organismal homeostasis by importing nutrients and exporting toxic compounds. Transporters also play a crucial role in drug response, serving as drug targets and setting drug levels. As part of a pharmacogenetics project, we screened exons and flanking intronic regions for variation in a set of 24 membrane transporter genes (96 kb; 57% coding) in 247 DNA samples from ethnically diverse populations. We identified 680 single nucleotide polymorphisms (SNPs), of which 175 were synonymous and 155 caused amino acid changes, and 29 small insertions and deletions. Amino acid diversity (pi(NS)) in transmembrane domains (TMDs) was significantly lower than in loop domains, suggesting that TMDs have special functional constraints. This difference was especially striking in the ATP-binding cassette superfamily and did not parallel evolutionary conservation: there was little variation in the TMDs, even in evolutionarily unconserved residues. We used allele frequency distribution to evaluate different scoring systems (Grantham, blosum62, SIFT, and evolutionarily conservedevolutionarily unconserved) for their ability to predict which SNPs affect function. Our underlying assumption was that alleles that are functionally deleterious will be selected against and thus under represented at high frequencies and over represented at low frequencies. We found that evolutionary conservation of orthologous sequences, as assessed by evolutionarily conservedevolutionarily unconserved and SIFT, was the best predictor of allele frequency distribution and hence of function. European Americans had an excess of high frequency alleles in comparison to African Americans, consistent with a historic bottleneck. In addition, African Americans exhibited a much higher frequency of population specific medium-frequency alleles than did European Americans.

Evolutionary conservation predicts function of variants of the human organic cation transporter, OCT1.

The organic cation transporter, OCT1, is a major hepatic transporter that mediates the uptake of many organic cations from the blood into the liver where the compounds may be metabolized or secreted into the bile. Because OCT1 interacts with a variety of structurally diverse organic cations, including clinically used drugs as well as toxic substances (e.g., N-methylpyridinium, MPP(+)), it is an important determinant of systemic exposure to many xenobiotics. To understand the genetic basis of extensive interindividual differences in xenobiotic disposition, we functionally characterized 15 protein-altering variants of the human liver organic cation transporter, OCT1, in Xenopus oocytes. All variants that reduced or eliminated function (OCT1-R61C, OCT1-P341L, OCT1-G220V, OCT1-G401S, and OCT1-G465R) altered evolutionarily conserved amino acid residues. In general, variants with decreased function had amino acid substitutions that resulted in more radical chemical changes (higher Grantham values) and were less evolutionarily favorable (lower blosum62 values) than variants that maintained function. A variant with increased function (OCT1-S14F) changed an amino acid residue such that the human protein matched the consensus of the OCT1 mammalian orthologs. Our results indicate that changes at evolutionarily conserved positions of OCT1 are strong predictors of decreased function and suggest that a combination of evolutionary conservation and chemical change might be a stronger predictor of function.

Polo-like kinase Cdc5 promotes chiasmata formation and cosegregation of sister centromeres at meiosis I.

During meiosis, two rounds of chromosome segregation occur after a single round of DNA replication, producing haploid progeny from diploid progenitors. Three innovations in chromosome behaviour during meiosis I accomplish this unique division. First, crossovers between maternal and paternal sister chromatids (detected cytologically as chiasmata) bind replicated maternal and paternal chromosomes together. Second, sister kinetochores attach to microtubules from the same pole (mono-polar orientation), causing maternal and paternal centromere pairs (and not sister chromatids) to be separated. Third, sister chromatid cohesion near centromeres is preserved at anaphase I when cohesion along chromosome arms is destroyed. The finding that destruction of mitotic cohesion is regulated by Polo-like kinases prompted us to investigate the meiotic role of the yeast Polo-like kinase Cdc5. We show here that cells lacking Cdc5 synapse homologues and initiate recombination normally, but fail to efficiently resolve recombination intermediates as crossovers. They also fail to properly localize the Lrs4 (ref. 3) and Mam1 (ref. 4) monopolin proteins, resulting in bipolar orientation of sister kinetochores. Cdc5 is thus required both for the formation of chiasmata and for cosegregation of sister centromeres at meiosis I.

SNP analysis and presentation in the Pharmacogenetics of Membrane Transporters Project.

The multidisciplinary UCSF Pharmacogenetics of Membrane Transporters project seeks to systematically identify sequence variants in transporters and to determine the functional significance of these variants through evaluation of relevant cellular and clinical phenotypes. The project is structured around four interacting cores: genomics, cellular phenotyping, clinical phenotyping, and bioinformatics. The bioinformatics core is responsible for collecting, storing, and analyzing the information obtained by the other cores and for presenting the results, in particular, for the genomic data. Most of this process is automated using locally developed software written in Python, an open source language well suited for rapid, modular development that meets requirements that are themselves constantly evolving. Here we present the details of transforming ABI trace file data into useful information for project investigators and a description of the types of data analysis and display that we have developed.

Bud morphogenesis and the actin and microtubule cytoskeletons during budding in the corn smut fungus, Ustilago maydis.

Ustilago maydis is a dimorphic Basidiomycete fungus with a yeast-like form and a hyphal form. Here we present a comprehensive analysis of bud formation and the actin and microtubule cytoskeletons of the yeast-like form during the cell cycle. We show that bud morphogenesis entails a series of shape changes, initially a tubular or conical structure, culminating in a cigar-shaped cell connected to the mother cell by a narrow neck. Labelling of cells with concanavalin A demonstrated that growth occurs at bud tip. Indirect immunofluorescence studies revealed that the actin cytoskeleton consists of patches and cables that polarize to the presumptive bud site and the bud tip and an actin ring that forms at the neck region. Because the bud tip corresponds to the site of active cell wall growth, we hypothesize that actin is involved in secretion of cell wall components. The microtubule cytoskeleton has recently been shown to consist of a cytoplasmic network during interphase that disassembles at mitosis when a spindle and astral microtubules are formed. We have carried out studies of U. maydis cells synchronized by the microtubule-depolymerizing drug thiabendazole which allow us to construct a temporal sequence of steps in spindle formation and spindle elongation during the cell cycle. These studies suggest that astral microtubules may be involved in early stages of spindle orientation and migration of the nucleus into the bud and that the spindle pole bodies may be involved in reestablishment of the cytoplasmic microtubule network.