Michiko Kawamoto

BayGenomics: a resource of insertional mutations in mouse embryonic stem cells

The BayGenomics gene-trap resource (http://baygenomics.ucsf.edu) provides researchers with access to thousands of mouse embryonic stem (ES) cell lines harboring characterized insertional mutations in both known and novel genes. Each cell line contains an insertional mutation in a specific gene. The identity of the gene that has been interrupted can be determined from a DNA sequence tag. Approximately 75% of our cell lines contain insertional mutations in known mouse genes or genes that share strong sequence similarities with genes that have been identified in other organisms. These cell lines readily transmit the mutation to the germline of mice and many mutant lines of mice have already been generated from this resource. BayGenomics provides facile access to our entire database, including sequence tags for each mutant ES cell line, through the World Wide Web. Investigators can browse our resource, search for specific entries, download any portion of our database and BLAST sequences of interest against our entire set of cell line sequence tags. They can then obtain the mutant ES cell line for the purpose of generating knockout mice.

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.

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 (π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 conserved/evolutionarily 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 conserved/evolutionarily 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.

Polymorphisms in a human kidney xenobiotic transporter, OCT2, exhibit altered function

The completion of the Human Genome Project and the development of high-throughput polymorphism identification methods have allowed researchers to carry out full genetic analyses of many clinically relevant genes. However, few studies have combined genetic analysis with in vitro phenotyping to better understand the relationship between genetic variation and protein function. Many transporters in the kidney are thought to play key roles in defense against a variety of foreign substances. The goal of this study was to understand the relationship between variation in a gene encoding a major renal xenobiotic transporter, OCT2, and transporter function. We report a comprehensive genetic analysis and functional characterization of variants of OCT2. Twenty-eight variable sites in the OCT2 gene were identified in a collection of 247 ethnically diverse DNA samples. Eight caused non-synonymous amino acid changes, of which four were present at >/= 1% in an ethnic population. All four of these altered transporter function assayed in Xenopus laevis oocytes. Analysis of nucleotide diversity (pi) revealed a higher prevalence of synonymous (pi = 22.4 x 10-4) versus non-synonymous (pi = 2.1 x 10-4) changes in OCT2 than in other genes. In addition, the non-synonymous sites had a significant tendency to exhibit more skewed allele frequencies (more negative Tajimas D-values) compared to synonymous sites. The population-genetic analysis, together with the functional characterization, suggests that selection has acted against amino acid changes in OCT2. This selection may be due to a necessary role of OCT2 in the renal elimination of endogenous amines or xenobiotics, including environmental toxins, neurotoxic amines and therapeutic drugs.

Interaction of Methotrexate with Organic-Anion Transporting Polypeptide 1A2 and Its Genetic Variants

Methotrexate (MTX) is used in patients with malignant and autoimmune diseases. This drug is primarily excreted unchanged in the urine, and its net excretion occurs via active secretory and reabsorptive processes. We characterized the interaction of MTX with human organic-anion transporting polypeptide transporter (OATP) 1A2, which is expressed in tissues important for MTX disposition and toxicity, such as the intestine, kidney, liver, and endothelial cells of the blood-brain barrier. In Xenopus laevis oocytes expressing OATP1A2, the uptake of the model substrate, estrone-3-sulfate (ES), was enhanced 30-fold compared with uninjected oocytes. MTX uptake in oocytes expressing OATP1A2 was saturable (Km = 457 ± 118 μM; Vmax = 17.5 ± 4.9 pmol/oocyte/60 min) and sensitive to extracellular pH. That is, acidic pHs stimulated MTX uptake by as much as 7-fold. Seven novel protein-altering variants were identified in 270 ethnically diverse DNA samples. Four protein-altering variants in OATP1A2 exhibited altered transport of ES and/or MTX. The common variant, protein reference sequence (p.) Ile13Thr, was hyperfunctional for ES and MTX and showed a 2-fold increase in the Vmax for ES. The common variant, p. Glu172Asp, exhibited reduced maximal transport capacity for ES and MTX. p. Arg168Cys was hypofunctional, and p. Asn277DEL was nonfunctional. Because of its expression on the apical membrane of the distal tubule and in tissues relevant to MTX disposition and toxicity, these findings suggest that OATP1A2 may play a role in active tubular reabsorption of MTX and in MTX-induced toxicities. Furthermore, genetic variation in OATP1A2 may contribute to variation in MTX disposition and response.

The International Gene Trap Consortium Website: a portal to all publicly available gene trap cell lines in mouse

Gene trapping is a method of generating murine embryonic stem (ES) cell lines containing insertional mutations in known and novel genes. A number of international groups have used this approach to create sizeable public cell line repositories available to the scientific community for the generation of mutant mouse strains. The major gene trapping groups worldwide have recently joined together to centralize access to all publicly available gene trap lines by developing a user-oriented Website for the International Gene Trap Consortium (IGTC). This collaboration provides an impressive public informatics resource comprising ∼45 000 well-characterized ES cell lines which currently represent ∼40% of known mouse genes, all freely available for the creation of knockout mice on a non-collaborative basis. To standardize annotation and provide high confidence data for gene trap lines, a rigorous identification and annotation pipeline has been developed combining genomic localization and transcript alignment of gene trap sequence tags to identify trapped loci. This information is stored in a new bioinformatics database accessible through the IGTC Website interface. The IGTC Website () allows users to browse and search the database for trapped genes, BLAST sequences against gene trap sequence tags, and view trapped genes within biological pathways. In addition, IGTC data have been integrated into major genome browsers and bioinformatics sites to provide users with outside portals for viewing this data. The development of the IGTC Website marks a major advance by providing the research community with the data and tools necessary to effectively use public gene trap resources for the large-scale characterization of mammalian gene function.

Functional analysis of polymorphisms in the organic anion transporter, SLC22A6 (OAT1)

Objectives The organic anion transporter, OAT1 (SLC22A6), plays a role in the renal elimination of many drugs and environmental toxins. The goal of this study was to identify and functionally characterize OAT1 variants as a first step towards understanding whether genetic variation in OAT1 may contribute to interindividual differences in renal elimination of xenobiotics. Methods As part of a larger study, 276 DNA samples from an ethnically diverse population were screened and 12 coding region variants of OAT1 were identified. The non-synonymous variants were then constructed and characterized in Xenopus laevis oocytes. A small family-based clinical study was conducted to determine the renal elimination of a model OAT1 substrate, adefovir (an antiviral agent) in human subjects who possessed a non-functional variant, OAT1-R454Q. Results Six non-synonymous variants were identified; two (OAT1-R50 H and OAT1-R293W) were present at ≥1% in at least one ethnic population. These two variants exhibited normal uptake of p-aminohippurate, ochratoxin A and methotrexate assayed in X. laevis oocytes. One variant, OAT1-R454Q, was non-functional with respect to the above substrates. In the clinical study, there was no significant decrease in the renal secretory clearance of adefovir in family members heterozygous for OAT1-454Q in comparison to those with the reference transporter, OAT1-454R. Conclusions These data indicate that the coding region of OAT1 has low genetic and functional diversity and suggest that coding region variants of OAT1 may not contribute substantially to interindividual differences in renal elimination of xenobiotics.

The Human Multidrug Resistance Protein 4 (MRP4, ABCC4) : Functional Analysis of a Highly Polymorphic Gene

ABCC4 encodes multidrug resistance protein 4 (MRP4), a member of the ATP-binding cassette family of membrane transporters involved in the efflux of endogenous and xenobiotic molecules. The aims of this study were to identify single nucleotide polymorphisms of ABCC4 and to functionally characterize selected nonsynonymous variants. Resequencing was performed in a large ethnically diverse population. Ten nonsynonymous variants were selected for analysis of transport function based on allele frequencies and evolutionary conservation. The reference and variant MRP4 cDNAs were constructed by site-directed mutagenesis and transiently transfected into human embryonic kidney cells (HEK 293T). The function of MRP4 variants was compared by measuring the intracellular accumulation of two antiviral agents, azidothymidine (AZT) and adefovir (PMEA). A total of 98 variants were identified in the coding and flanking intronic regions of ABCC4. Of these, 43 variants are in the coding region, and 22 are nonsynonymous. In a functional screen of ten variants, there was no evidence for a complete loss of function allele. However, two variants (G187W and G487E) showed a significantly reduced function compared to reference with both substrates, as evidenced by higher intracellular accumulation of AZT and PMEA compared to the reference MRP4 (43 and 69% increase in accumulation for G187W compared with the reference MRP4, with AZT and PMEA, respectively). The G187W variant also showed decreased expression following transient transfection of HEK 293T cells. Further studies are required to assess the clinical significance of this altered function and expression and to evaluate substrate specificity of this functional change.

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.

Genetic analysis and functional characterization of polymorphisms in the human concentrative nucleoside transporter, CNT2.

The concentrative nucleoside transporter CNT2 (SPNT1; SLC28A2) plays a role in the absorption and disposition of naturally occurring nucleosides, as well as nucleoside analog drugs. The aim of the present study was to characterize genetic variation in SLC28A2, the gene encoding CNT2, and to functionally analyse non-synonymous variants of CNT2, as a first step towards understanding whether genetic variation in this nucleoside transporter contributes to variation in response to nucleoside analogs. As part of a larger study, DNA samples from an ethnically diverse population (100 African-Americans, 100 European-Americans, 30 Asians, 10 Mexicans and seven Pacific Islanders) were screened and 10 coding region variants of CNT2 were identified. The non-synonymous variants were then constructed and characterized in Xenopus laevis oocytes. Six non-synonymous variants were identified, and all were able to transport guanosine. The four common variants (>1% in the sample population) were further characterized with the anti-viral nucleoside analog drug ribavirin. No differences were observed among the four common variants in the uptake kinetics of 3H-ribavirin (Km in μM: 35.6±9.27 for CNT2-reference, 40.7±6.47 for CNT2-P22 l, 31.2±15.8 for CNT2-S75R, 26.7±6.13 for CNT2-S245T and 49.9±14.6 for CNT2-F355S). The variant CNT2-F355S exhibited a change in specificity for the naturally occurring nucleosides, inosine and uridine. All non-synonymous variants of CNT2 took up guanosine, and the four variants examined showed no significant difference in ribavirin kinetics. However, CNT2-F355S (3% allele frequency in the African-American sample) was found to alter specificity for naturally occurring nucleosides, which may have implications for nucleoside homeostasis.