Meta-analysis of erosive hand osteoarthritis identifies four common variants that associate with relatively large effect.

OBJECTIVES: Erosive hand osteoarthritis (EHOA) is a severe subset of hand osteoarthritis (OA). It is unclear if EHOA is genetically different from other forms of OA. Sequence variants at ten loci have been associated with hand OA but none with EHOA. METHODS: We performed meta-analysis of EHOA in 1484 cases and 550 680 controls, from 5 populations. To identify causal genes, we performed eQTL and plasma pQTL analyses, and developed one zebrafish mutant. We analysed associations of variants with other traits and estimated shared genetics between EHOA and other traits. RESULTS: Four common sequence variants associated with EHOA, all with relatively high effect. Rs17013495 (SPP1/MEPE, OR=1.40, p=8.4×10-14) and rs11243284 (6p24.3, OR=1.35, p=4.2×10-11) have not been associated with OA, whereas rs11631127 (ALDH1A2, OR=1.46, p=7.1×10-18), and rs1800801 (MGP, OR=1.37, p=3.6×10-13) have previously been associated with hand OA. The association of rs1800801 (MGP) was consistent with a recessive mode of inheritance in contrast to its additive association with hand OA (OR homozygotes vs non-carriers=2.01, 95% CI 1.71 to 2.37). All four variants associated nominally with finger OA, although with substantially lower effect. We found shared genetic components between EHOA and other OA measures, grip strength, urate levels and gout, but not rheumatoid arthritis. We identified ALDH1A2, MGP and BMP6 as causal genes for EHOA, with loss-of-function Bmp6 zebrafish mutants displaying EHOA-like phenotypes. CONCLUSIONS: We report on significant genetic associations with EHOA. The results support the view of EHOA as a form of severe hand OA and partly separate it from OA in larger joints.

NOD/RIPK2 signalling pathway contributes to osteoarthritis susceptibility.

OBJECTIVES: How inflammatory signalling contributes to osteoarthritis (OA) susceptibility is undetermined. An allele encoding a hyperactive form of the Receptor Interacting Protein Kinase 2 (RIPK2) proinflammatory signalling intermediate has been associated with familial OA. To test whether altered nucleotide-binding oligomerisation domain (NOD)/RIPK2 pathway activity causes heightened OA susceptibility, we investigated whether variants affecting additional pathway components are associated with familial OA. To determine whether the Ripk2104Asp disease allele is sufficient to account for the familial phenotype, we determined the effect of the allele on mice. METHODS: Genomic analysis of 150 independent families with dominant inheritance of OA affecting diverse joints was used to identify coding variants that segregated strictly with occurrence of OA. Genome editing was used to introduce the OA-associated RIPK2 (p.Asn104Asp) allele into the genome of inbred mice. The consequences of the Ripk2104Asp disease allele on physiology and OA susceptibility in mice were measured by histology, immunohistochemistry, serum cytokine levels and gene expression. RESULTS: We identified six novel variants affecting components of the NOD/RIPK2 inflammatory signalling pathway that are associated with familial OA affecting the hand, shoulder or foot. The Ripk2104Asp allele acts dominantly to alter basal physiology and response to trauma in the mouse knee. Whereas the knees of uninjured Ripk2Asp104 mice appear normal histologically, the joints exhibit a set of marked gene expression changes reminiscent of overt OA. Although the Ripk2104Asp mice lack evidence of chronically elevated systemic inflammation, they do exhibit significantly increased susceptibility to post-traumatic OA (PTOA). CONCLUSIONS: Two types of data support the hypothesis that altered NOD/RIPK2 signalling confers susceptibility to OA.

Shutdown corner, a large deletion mutant isolated from a haploid mutagenesis screen in zebrafish.

Morphogenesis, the formation of three-dimensional organ structures, requires precise coupling of genetic regulation and complex cell behaviors. The genetic networks governing many morphogenetic systems, including that of the embryonic eye, are poorly understood. In zebrafish, several forward genetic screens have sought to identify factors regulating eye development. These screens often look for eye defects at stages after the optic cup is formed and when retinal neurogenesis is under way. This approach can make it difficult to identify mutants specific for morphogenesis, as opposed to neurogenesis. To this end, we carried out a forward genetic, small-scale haploid mutagenesis screen in zebrafish (Danio rerio) to identify factors that govern optic cup morphogenesis. We screened ∼100 genomes and isolated shutdown corner (sco), a mutant that exhibits multiple tissue defects and harbors a ∼10-Mb deletion that encompasses 89 annotated genes. Using a combination of live imaging and antibody staining, we found cell proliferation, cell death, and tissue patterning defects in the sco optic cup. We also observed other phenotypes, including paralysis, neuromuscular defects, and ocular vasculature defects. To date, the largest deletion mutants reported in zebrafish are engineered using CRISPR-Cas9 and are less than 300 kb. Because of the number of genes within the deletion interval, shutdown corner [Df(Chr05:sco)z207] could be a useful resource to the zebrafish community, as it may be helpful for gene mapping, understanding genetic interactions, or studying many genes lost in the mutant.

The Paf1 complex and P-TEFb have reciprocal and antagonist roles in maintaining multipotent neural crest progenitors.

Multipotent progenitor populations are necessary for generating diverse tissue types during embryogenesis. We show the RNA polymerase-associated factor 1 complex (Paf1C) is required to maintain multipotent progenitors of the neural crest (NC) lineage in zebrafish. Mutations affecting each Paf1C component result in near-identical NC phenotypes; alyron mutant embryos carrying a null mutation in paf1 were analyzed in detail. In the absence of zygotic paf1 function, definitive premigratory NC progenitors arise but fail to maintain expression of the sox10 specification gene. The mutant NC progenitors migrate aberrantly and fail to differentiate appropriately. Blood and germ cell progenitor development is affected similarly. Development of mutant NC could be rescued by additional loss of positive transcription elongation factor b (P-TEFb) activity, a key factor in promoting transcription elongation. Consistent with the interpretation that inhibiting/delaying expression of some genes is essential for maintaining progenitors, mutant embryos lacking the CDK9 kinase component of P-TEFb exhibit a surfeit of NC progenitors and their derivatives. We propose Paf1C and P-TEFb act antagonistically to regulate the timing of the expression of genes needed for NC development.

Highly Efficient CRISPR-Cas9-Based Methods for Generating Deletion Mutations and F0 Embryos that Lack Gene Function in Zebrafish.

Inconsistent activity limits the use of CRISPR-Cas9 in zebrafish. We show supernumerary guanine nucleotides at the 5' ends of single guide RNAs (sgRNAs) account for diminished CRISPR-Cas9 activity in zebrafish embryos. Genomic sequences can be targeted consistently with extremely high efficiency using Cas9 ribonucleoproteins (RNPs) containing either a sgRNA molecule or a synthetic crRNA:tracrRNA duplex that perfectly matches the protospacer target site. Following injection of zebrafish eggs with such RNPs, virtually every copy of a targeted locus harbors an induced indel mutation. Loss of gene function is often complete, as F0 embryos closely resemble true null mutants without detectable non-specific effects. Mosaicism is sufficiently low in F0 embryos that cell non-autonomous gene functions can be probed effectively and redundant activities of genes can be uncovered when two genes are targeted simultaneously. Finally, heritable deletion mutations of at least 50 kbp can be readily induced using pairs of duplex guide RNPs targeted to a single chromosome.

Interactions among ryanodine receptor isotypes contribute to muscle fiber type development and function.

Mutations affecting ryanodine receptor (RyR) calcium release channels commonly underlie congenital myopathies. Although these channels are known principally for their essential roles in muscle contractility, mutations in the human RYR1 gene result in a broad spectrum of phenotypes, including muscle weakness, altered proportions of fiber types, anomalous muscle fibers with cores or centrally placed nuclei, and dysmorphic craniofacial features. Currently, it is unknown which phenotypes directly reflect requirements for RyRs and which result secondarily to aberrant muscle function. To identify biological processes requiring RyR function, skeletal muscle development was analyzed in zebrafish embryos harboring protein-null mutations. RyR channels contribute to both muscle fiber development and function. Loss of some RyRs had modest effects, altering muscle fiber-type specification in the embryo without compromising viability. In addition, each RyR-encoding gene contributed to normal swimming behavior and muscle function. The RyR channels do not function in a simple additive manner. For example, although isoform RyR1a is sufficient for muscle contraction in the absence of RyR1b, RyR1a normally attenuates the activity of the co-expressed RyR1b channel in slow muscle. RyR3 also acts to modify the functions of other RyR channels. Furthermore, diminished RyR-dependent contractility affects both muscle fiber maturation and craniofacial development. These findings help to explain some of the heterogeneity of phenotypes that accompany RyR1 mutations in humans.

tbx6l and tbx16 are redundantly required for posterior paraxial mesoderm formation during zebrafish embryogenesis.

BACKGROUND: T-box genes encode a large transcription factor family implicated in many aspects of development. We are focusing on two related zebrafish T-box genes, tbx6l and tbx16, that are expressed in highly overlapping patterns in embryonic paraxial mesoderm. tbx16 mutants are deficient in trunk, but not tail, somites; we explored whether presence of tail somites in tbx16 mutants was due to compensatory function provided by the tbx6l gene. RESULTS: We generated two zebrafish tbx6l mutant alleles. Loss of tbx6l has no apparent effect on embryonic development, nor does tbx6l loss enhance the phenotype of two other T-box gene mutants, ta and tbx6, or of the mesp family gene mutant msgn1. In contrast, loss of tbx6l function dramatically enhances the paraxial mesoderm deficiency of tbx16 mutants. CONCLUSIONS: These data demonstrate that tbx6l and tbx16 genes function redundantly to direct tail somite development. tbx6l single mutants develop normally because tbx16 fully compensates for loss of tbx6l function. However, tbx6l only partially compensates for loss of tbx16 function. These results resolve the question of why loss of function of tbx16 gene, which is expressed throughout the ventral and paraxial mesoderm, profoundly affects somite development in the trunk but not the tail. Developmental Dynamics 246:759-769, 2017. © 2017 Wiley Periodicals, Inc.

Quadruple zebrafish mutant reveals different roles of Mesp genes in somite segmentation between mouse and zebrafish.

The segmental pattern of somites is generated by sequential conversion of the temporal periodicity provided by the molecular clock. Whereas the basic structure of this clock is conserved among different species, diversity also exists, especially in terms of the molecular network. The temporal periodicity is subsequently converted into the spatial pattern of somites, and Mesp2 plays crucial roles in this conversion in the mouse. However, it remains unclear whether Mesp genes play similar roles in other vertebrates. In this study, we generated zebrafish mutants lacking all four zebrafish Mesp genes by using TALEN-mediated genome editing. Contrary to the situation in the mouse Mesp2 mutant, in the zebrafish Mesp quadruple mutant embryos the positions of somite boundaries were clearly determined and morphological boundaries were formed, although their formation was not completely normal. However, each somite was caudalized in a similar manner to the mouse Mesp2 mutant, and the superficial horizontal myoseptum and lateral line primordia were not properly formed in the quadruple mutants. These results clarify the conserved and species-specific roles of Mesp in the link between the molecular clock and somite morphogenesis.

Precise Editing of the Zebrafish Genome Made Simple and Efficient.

We present simple and efficient methods for creating heritable modifications of the zebrafish genome. Precisely modified alleles are generated by homologous recombination between the host genome and dsDNA donor molecules, stimulated by the induction of chromosomally targeted double-strand breaks. Several kilobase-long tracts of genome sequence can be replaced. Tagging donor sequences with reporter genes that can be subsequently excised improves recovery of edited alleles by an order of magnitude and facilitates recovery of recessive and phenotypically silent conditional mutations. We generate and demonstrate functionality of (1) alleles with a single codon change, (2) an allele encoding an epitope-tagged version of an endogenous protein, (3) alleles expressing reporter proteins, and (4) a conditional allele in which an exon is flanked by recombinogenic loxP sites. Our methods make recovery of a broad range of genome editing events very practicable, significantly advancing applicability of the zebrafish for studying normal biological processes and modeling diseases.

Close Association of Carbonic Anhydrase (CA2a and CA15a), Na(+)/H(+) Exchanger (Nhe3b), and Ammonia Transporter Rhcg1 in Zebrafish Ionocytes Responsible for Na(+) Uptake.

Freshwater (FW) fishes actively absorb salt from their environment to tolerate low salinities. We previously reported that vacuolar-type H(+)-ATPase/mitochondrion-rich cells (H-MRCs) on the skin epithelium of zebrafish larvae (Danio rerio) are primary sites for Na(+) uptake. In this study, in an attempt to clarify the mechanism for the Na(+) uptake, we performed a systematic analysis of gene expression patterns of zebrafish carbonic anhydrase (CA) isoforms and found that, of 12 CA isoforms, CA2a and CA15a are highly expressed in H-MRCs at larval stages. The ca2a and ca15a mRNA expression were salinity-dependent; they were upregulated in 0.03 mM Na(+) water whereas ca15a but not ca2a was down-regulated in 70 mM Na(+) water. Immunohistochemistry demonstrated cytoplasmic distribution of CA2a and apical membrane localization of CA15a. Furthermore, cell surface immunofluorescence staining revealed external surface localization of CA15a. Depletion of either CA2a or CA15a expression by Morpholino antisense oligonucleotides resulted in a significant decrease in Na(+) accumulation in H-MRCs. An in situ proximity ligation assay demonstrated a very close association of CA2a, CA15a, Na(+)/H(+) exchanger 3b (Nhe3b), and Rhcg1 ammonia transporter in H-MRC. Our findings suggest that CA2a, CA15a, and Rhcg1 play a key role in Na(+)uptake under FW conditions by forming a transport metabolon with Nhe3b.

Efficient identification of TALEN-mediated genome modifications using heteroduplex mobility assays.

The heteroduplex mobility assay (HMA) is widely used to characterize strain variants of human viruses. To determine whether it can detect small sequence differences in homologous templates, we constructed a series of deletion constructs (1-10 bp deletions) in the multiple cloning site (MCS) of pBluescript II. After PCR amplification of the MCS using a mixture of wild-type and one of the deletion constructs, the resulting PCR amplicons were electrophoresed using 15% polyacrylamide gels. Two types of heteroduplexes exhibited retarded electrophoretic migration compared with individual homoduplexes. Therefore, we applied this HMA to detect transcription activator-like effector nucleases (TALEN)-induced insertion and/or deletion (indel) mutations at an endogenous locus. We found that TALEN in vivo activity was easily estimated by the degree of multiple HMA profiles derived from TALEN-injected F0 embryos. Furthermore, TALEN-injected F0 founder fish produced several unique HMA profiles in F1 embryos. Sequence analysis confirmed that the different HMA profiles contained distinct indel mutations. Thus, HMA is a rapid and sensitive analytical method for the detection of the TALEN-mediated genome modifications.

Zebrafish foxP2 zinc finger nuclease mutant has normal axon pathfinding.

foxP2, a forkhead-domain transcription factor, is critical for speech and language development in humans, but its role in the establishment of CNS connectivity is unclear. While in vitro studies have identified axon guidance molecules as targets of foxP2 regulation, and cell culture assays suggest a role for foxP2 in neurite outgrowth, in vivo studies have been lacking regarding a role for foxP2 in axon pathfinding. We used a modified zinc finger nuclease methodology to generate mutations in the zebrafish foxP2 gene. Using PCR-based high resolution melt curve analysis (HRMA) of G0 founder animals, we screened and identified three mutants carrying nonsense mutations in the 2(nd) coding exon: a 17 base-pair (bp) deletion, an 8bp deletion, and a 4bp insertion. Sequence analysis of cDNA confirmed that these were frameshift mutations with predicted early protein truncations. Homozygous mutant fish were viable and fertile, with unchanged body morphology, and no apparent differences in CNS apoptosis, proliferation, or patterning at embryonic stages. There was a reduction in expression of the known foxP2 target gene cntnap2 that was rescued by injection of wild-type foxP2 transcript. When we examined axon pathfinding using a pan-axonal marker or transgenic lines, including a foxP2-neuron-specific enhancer, we did not observe any axon guidance errors. Our findings suggest that foxP2 is not necessary for axon pathfinding during development.

Wnt signaling regulates postembryonic hypothalamic progenitor differentiation.

Previous studies have raised the possibility that Wnt signaling may regulate both neural progenitor maintenance and neuronal differentiation within a single population. Here we investigate the role of Wnt/ß-catenin activity in the zebrafish hypothalamus and find that the pathway is first required for the proliferation of unspecified hypothalamic progenitors in the embryo. At later stages, including adulthood, sequential activation and inhibition of Wnt activity is required for the differentiation of neural progenitors and negatively regulates radial glia differentiation. The presence of Wnt activity is conserved in hypothalamic progenitors of the adult mouse, where it plays a conserved role in inhibiting the differentiation of radial glia. This study establishes the vertebrate hypothalamus as a model for Wnt-regulated postembryonic neural progenitor differentiation and defines specific roles for Wnt signaling in neurogenesis.

Simple methods for generating and detecting locus-specific mutations induced with TALENs in the zebrafish genome.

The zebrafish is a powerful experimental system for uncovering gene function in vertebrate organisms. Nevertheless, studies in the zebrafish have been limited by the approaches available for eliminating gene function. Here we present simple and efficient methods for inducing, detecting, and recovering mutations at virtually any locus in the zebrafish. Briefly, double-strand DNA breaks are induced at a locus of interest by synthetic nucleases, called TALENs. Subsequent host repair of the DNA lesions leads to the generation of insertion and deletion mutations at the targeted locus. To detect the induced DNA sequence alterations at targeted loci, genomes are examined using High Resolution Melt Analysis, an efficient and sensitive method for detecting the presence of newly arising sequence polymorphisms. As the DNA binding specificity of a TALEN is determined by a custom designed array of DNA recognition modules, each of which interacts with a single target nucleotide, TALENs with very high target sequence specificities can be easily generated. Using freely accessible reagents and Web-based software, and a very simple cloning strategy, a TALEN that uniquely recognizes a specific pre-determined locus in the zebrafish genome can be generated within days. Here we develop and test the activity of four TALENs directed at different target genes. Using the experimental approach described here, every embryo injected with RNA encoding a TALEN will acquire targeted mutations. Multiple independently arising mutations are produced in each growing embryo, and up to 50% of the host genomes may acquire a targeted mutation. Upon reaching adulthood, approximately 90% of these animals transmit targeted mutations to their progeny. Results presented here indicate the TALENs are highly sequence-specific and produce minimal off-target effects. In all, it takes about two weeks to create a target-specific TALEN and generate growing embryos that harbor an array of germ line mutations at a pre-specified locus.

Identification of zebrafish Fxyd11a protein that is highly expressed in ion-transporting epithelium of the gill and skin and its possible role in ion homeostasis.

FXYD proteins, small single-transmembrane proteins, have been proposed to be auxiliary regulatory subunits of Na(+)-K(+)-ATPase and have recently been implied in ion osmoregulation of teleost fish. In freshwater (FW) fish, numerous ions are actively taken up through mitochondrion-rich cells (MRCs) of the gill and skin epithelia, using the Na(+) electrochemical gradient generated by Na(+)-K(+)-ATPase. In the present study, to understand the molecular mechanism for the regulation of Na(+)-K(+)-ATPase in MRCs of FW fish, we sought to identify FXYD proteins expressed in MRCs of zebrafish. Reverse-transcriptase PCR studies of adult zebrafish tissues revealed that, out of eight fxyd genes found in zebrafish database, only zebrafish fxyd11 (zfxyd11) mRNA exhibited a gill-specific expression. Double immunofluorescence staining showed that zFxyd11 is abundantly expressed in MRCs rich in Na(+)-K(+)-ATPase (NaK-MRCs) but not in those rich in vacuolar-type H(+)-transporting ATPase. An in situ proximity ligation assay demonstrated its close association with Na(+)-K(+)-ATPase in NaK-MRCs. The zfxyd11 mRNA expression was detectable at 1 day postfertilization, and its expression levels in the whole larvae and adult gills were regulated in response to changes in environmental ionic concentrations. Furthermore, knockdown of zFxyd11 resulted in a significant increase in the number of Na(+)-K(+)-ATPase-positive cells in the larval skin. These results suggest that zFxyd11 may regulate the transport ability of NaK-MRCs by modulating Na(+)-K(+)-ATPase activity, and may be involved in the regulation of body fluid and electrolyte homeostasis.

Chromosomal position mediates spinal cord expression of a dbx1a enhancer.

Dbx homeodomain proteins are important for the production of multiple spinal cord cell types. To examine the regulation of Dbx genes in more detail, we have generated transgenic zebrafish in which fluorescent protein expression is driven by predicted dbx1a enhancers. We identified three areas of sequence conservation upstream of the dbx1a coding sequence and generated fluorescent reporter constructs driven by these predicted enhancer elements and the endogenous dbx1a promoter. In multiple stable insertions of a 3.5-kb enhancer fragment, we observed that there was additional reporter expression in the dorsal spinal cord not normally observed by dbx1a in situ hybridization. In addition, these lines exhibited only transient reporter expression, unlike the endogenous gene. Surprisingly, a single insertion line expressed the reporter in the endogenous pattern, indicating that other local regulatory elements modulate gene expression through the 3.5-kb enhancer.

Mechanism of development of ionocytes rich in vacuolar-type H(+)-ATPase in the skin of zebrafish larvae.

Mitochondrion-rich cells (MRCs), or ionocytes, play a central role in aquatic species, maintaining body fluid ionic homeostasis by actively taking up or excreting ions. Since their first description in 1932 in eel gills, extensive morphological and physiological analyses have yielded important insights into ionocyte structure and function, but understanding the developmental pathway specifying these cells remains an ongoing challenge. We previously succeeded in identifying a key transcription factor, Foxi3a, in zebrafish larvae by database mining. In the present study, we analyzed a zebrafish mutant, quadro (quo), deficient in foxi1 gene expression and found that foxi1 is essential for development of an MRC subpopulation rich in vacuolar-type H(+)-ATPase (vH-MRC). foxi1 acts upstream of Delta-Notch signaling that determines sporadic distribution of vH-MRC and regulates foxi3a expression. Through gain- and loss-of-function assays and cell transplantation experiments, we further clarified that (1) the expression level of foxi3a is maintained by a positive feedback loop between foxi3a and its downstream gene gcm2 and (2) Foxi3a functions cell-autonomously in the specification of vH-MRC. These observations provide a better understanding of the differentiation and distribution of the vH-MRC subtype.

Drug-induced phototoxicity evoked by inhibition of human ABC transporter ABCG2: development of in vitro high-speed screening systems.

BACKGROUND: Photosensitivity depends on both genetic and environmental factors. Pheophorbide a, present in various plant-derived foods and food supplements, can be absorbed by the small intestine. Accumulation of pheophorbide a and porphyrins in the systemic blood circulation can result in phototoxic lesions on light-exposed skin. OBJECTIVE: As the human ATP-binding cassette (ABC) transporter ABCG2 has been suggested to be critically involved in porphyrin-mediated photosensitivity, we aimed to develop in vitro screening systems for drug-induced phototoxicity. CONCLUSION: Functional impairment owing to inhibition of ABCG2 by drugs or its genetic polymorphisms can lead to the disruption of porphyrin homeostasis. This review article provides an overview on drug-induced photosensitivity, as well as our hypothesis on a potential role of ABCG2 in phototoxicity.

Zebrafish early cardiac connexin, Cx36.7/Ecx, regulates myofibril orientation and heart morphogenesis by establishing Nkx2.5 expression.

Heart development is a precisely coordinated process of cellular proliferation, migration, differentiation, and integrated morphogenetic interactions, and therefore it is highly susceptible to developmental anomalies such as the congenital heart disease (CHD). One of the major causes of CHD has been shown to be the mutations in key cardiac transcription factors, including nkx2.5. Here, we report the analysis of zebrafish mutant ftk that showed a progressive heart malformation in the later stages of heart morphogenesis. Our analyses revealed that the cardiac muscle maturation and heart morphogenesis in ftk mutants were impaired because of the disorganization of myofibrils. Notably, we found that the expression of nkx2.5 was down-regulated in the ftk heart despite the normal expression of gata4 and tbx5, suggesting a common mechanism for the occurrence of ftk phenotype and CHD. We identified ftk to be a loss-of-function mutation in a connexin gene, cx36.7/early cardiac connexin (ecx), expressed during early heart development. We further showed by a rescue experiment that Nkx2.5 is the downstream mediator of Ecx-mediated signaling. From these results, we propose that the cardiac connexin Ecx and its downstream signaling are crucial for establishing nkx2.5 expression, which in turn promotes unidirectional, parallel alignment of myofibrils and the subsequent proper heart morphogenesis.

Ubiquitin-mediated proteasomal degradation of non-synonymous SNP variants of human ABC transporter ABCG2.

Clinical relevance is implicated between the genetic polymorphisms of the ABC (ATP-binding cassette) transporter ABCG2 (ABC subfamily G, member 2) and the individual differences in drug response. We expressed a total of seven non-synonymous SNP (single nucleotide polymorphism) variants in Flp-In-293 cells by using the Flp (flippase) recombinase system. Of these, ABCG2 F208S and S441N variants were found to be expressed at markedly low levels, whereas their mRNA levels were equal to those of the other SNP variants and ABCG2 WT (wild-type). Interestingly, protein expression levels of the ABCG2 F208S and S441N variants increased 6- to 12-fold when Flp-In-293 cells were treated with MG132, a proteasome inhibitor. Immunoprecipitation followed by immunoblot analysis showed that the ABCG2 F208S and S441N variant proteins were endogenously ubiquitinated in Flp-In-293 cells, and treatment with MG132 significantly enhanced the level of these ubiquitinated variants. Immunofluorescence microscopy demonstrated that MG132 greatly affected the ABCG2 F208S and S441N variants in terms of both protein levels and intracellular distribution. Immunoblot analysis revealed that those variants were N-glycosylated; however, their oligosaccharides were immature compared with those present on ABCG2 WT. The ABCG2 F208S and S441N variant proteins do not appear to be processed in the Golgi apparatus, but undergo ubiquitin-mediated protein degradation in proteasomes, whereas ABCG2 WT is sorted to the plasma membrane and then degraded via the lysosomal pathway. The present study provides the first evidence that certain genetic polymorphisms can affect the protein stability of ABCG2. Control of proteasomal degradation of ABCG2 would provide a novel approach in cancer chemotherapy to circumvent multidrug resistance of human cancers.