Frequent loss of heterozygosity in CRISPR-Cas9-edited early human embryos.

CRISPR-Cas9 genome editing is a promising technique for clinical applications, such as the correction of disease-associated alleles in somatic cells. The use of this approach has also been discussed in the context of heritable editing of the human germ line. However, studies assessing gene correction in early human embryos report low efficiency of mutation repair, high rates of mosaicism, and the possibility of unintended editing outcomes that may have pathologic consequences. We developed computational pipelines to assess single-cell genomics and transcriptomics datasets from OCT4 (POU5F1) CRISPR-Cas9-targeted and control human preimplantation embryos. This allowed us to evaluate on-target mutations that would be missed by more conventional genotyping techniques. We observed loss of heterozygosity in edited cells that spanned regions beyond the POU5F1 on-target locus, as well as segmental loss and gain of chromosome 6, on which the POU5F1 gene is located. Unintended genome editing outcomes were present in ∼16% of the human embryo cells analyzed and spanned 4-20 kb. Our observations are consistent with recent findings indicating complexity at on-target sites following CRISPR-Cas9 genome editing. Our work underscores the importance of further basic research to assess the safety of genome editing techniques in human embryos, which will inform debates about the potential clinical use of this technology.

Genetic profile of T-cell acute lymphoblastic leukemias with MYC translocations.

MYC translocations represent a genetic subtype of T-lineage acute lymphoblastic leukemia (T-ALL), which occurs at an incidence of ∼6%, assessed within a cohort of 196 T-ALL patients (64 adults and 132 children). The translocations were of 2 types; those rearranged with the T-cell receptor loci and those with other partners. MYC translocations were significantly associated with the TAL/LMO subtype of T-ALL (P = .018) and trisomies 6 (P < .001) and 7 (P < .001). Within the TAL/LMO subtype, gene expression profiling identified 148 differentially expressed genes between patients with and without MYC translocations; specifically, 77 were upregulated and 71 downregulated in those with MYC translocations. The poor prognostic marker, CD44, was among the upregulated genes. MYC translocations occurred as secondary abnormalities, present in subclones in one-half of the cases. Longitudinal studies indicated an association with induction failure and relapse.

The kinetics of relapse in DEK-NUP214-positive acute myeloid leukemia patients.

Preemptive treatment of relapse of acute myeloid leukemia (AML) holds the promise to improve the prognosis of this currently highly lethal condition. Proposed treatment modalities applicable in preemptive cytoreduction (e.g., demethylating agents or standard chemotherapy) differ substantially in interval from administration to antileukemic effect. The t(6;9) balanced translocation, producing the DEK-NUP214 fusion protein, is seen in only 1% of patients with AML. We hypothesized that in these patients, who relapse with a very high frequency, a more detailed knowledge of leukemic relapse growth kinetics would improve the personalized decision-making regarding re-administration of chemotherapy. Based on standardized quantitative PCR data, we therefore delineated the relapse kinetics in a cohort of 27 relapsing DEK-NUP214-positive patients treated in four different European countries. The prerelapse leukemic burden increased with a median doubling time of 13 d (range: 5-51 d, median: 0.71 logs/month, range: 0.18-1.91 logs/month), with FLT3-ITD-positive patients relapsing significantly faster than FLT3-ITD-negative ones (median: 0.9 vs. 0.6 logs/month, Wilcoxon rank sum test, P = 0.041). Peripheral blood and bone marrow were equally useful for minimal residual disease (MRD) detection, and thus, we found that with sampling intervals of 2 months, 94% of relapses would be detected with a median time from MRD detection to hematological relapse of 64 d. In conclusion, this data provide algorithms for handling the rare patients with DEK-NUP214-positive AML allowing for planning of both MRD follow-up and, upon molecular relapse, the timing of cytoreduction or possibly transplant procedures.

Hypomorphic homozygous mutations in phosphoglucomutase 3 (PGM3) impair immunity and increase serum IgE levels.

BACKGROUND: Recurrent bacterial and fungal infections, eczema, and increased serum IgE levels characterize patients with the hyper-IgE syndrome (HIES). Known genetic causes for HIES are mutations in signal transducer and activator of transcription 3 (STAT3) and dedicator of cytokinesis 8 (DOCK8), which are involved in signal transduction pathways. However, glycosylation defects have not been described in patients with HIES. One crucial enzyme in the glycosylation pathway is phosphoglucomutase 3 (PGM3), which catalyzes a key step in the synthesis of uridine diphosphate N-acetylglucosamine, which is required for the biosynthesis of N-glycans. OBJECTIVE: We sought to elucidate the genetic cause in patients with HIES who do not carry mutations in STAT3 or DOCK8. METHODS: After establishing a linkage interval by means of SNPchip genotyping and homozygosity mapping in 2 families with HIES from Tunisia, mutational analysis was performed with selector-based, high-throughput sequencing. Protein expression was analyzed by means of Western blotting, and glycosylation was profiled by using mass spectrometry. RESULTS: Mutational analysis of candidate genes in an 11.9-Mb linkage region on chromosome 6 shared by 2 multiplex families identified 2 homozygous mutations in PGM3 that segregated with disease status and followed recessive inheritance. The mutations predict amino acid changes in PGM3 (p.Glu340del and p.Leu83Ser). A third homozygous mutation (p.Asp502Tyr) and the p.Leu83Ser variant were identified in 2 other affected families, respectively. These hypomorphic mutations have an effect on the biosynthetic reactions involving uridine diphosphate N-acetylglucosamine. Glycomic analysis revealed an aberrant glycosylation pattern in leukocytes demonstrated by a reduced level of tri-antennary and tetra-antennary N-glycans. T-cell proliferation and differentiation were impaired in patients. Most patients had developmental delay, and many had psychomotor retardation. CONCLUSION: Impairment of PGM3 function leads to a novel primary (inborn) error of development and immunity because biallelic hypomorphic mutations are associated with impaired glycosylation and a hyper-IgE-like phenotype.

DNA methylation profiling of human chromosomes 6, 20 and 22.

DNA methylation is the most stable type of epigenetic modification modulating the transcriptional plasticity of mammalian genomes. Using bisulfite DNA sequencing, we report high-resolution methylation profiles of human chromosomes 6, 20 and 22, providing a resource of about 1.9 million CpG methylation values derived from 12 different tissues. Analysis of six annotation categories showed that evolutionarily conserved regions are the predominant sites for differential DNA methylation and that a core region surrounding the transcriptional start site is an informative surrogate for promoter methylation. We find that 17% of the 873 analyzed genes are differentially methylated in their 5' UTRs and that about one-third of the differentially methylated 5' UTRs are inversely correlated with transcription. Despite the fact that our study controlled for factors reported to affect DNA methylation such as sex and age, we did not find any significant attributable effects. Our data suggest DNA methylation to be ontogenetically more stable than previously thought.

Fine mapping of breast cancer genome-wide association studies loci in women of African ancestry identifies novel susceptibility markers.

Numerous single nucleotide polymorphisms (SNPs) associated with breast cancer susceptibility have been identified by genome-wide association studies (GWAS). However, these SNPs were primarily discovered and validated in women of European and Asian ancestry. Because linkage disequilibrium is ancestry-dependent and heterogeneous among racial/ethnic populations, we evaluated common genetic variants at 22 GWAS-identified breast cancer susceptibility loci in a pooled sample of 1502 breast cancer cases and 1378 controls of African ancestry. None of the 22 GWAS index SNPs could be validated, challenging the direct generalizability of breast cancer risk variants identified in Caucasians or Asians to other populations. Novel breast cancer risk variants for women of African ancestry were identified in regions including 5p12 (odds ratio [OR] = 1.40, 95% confidence interval [CI] = 1.11-1.76; P = 0.004), 5q11.2 (OR = 1.22, 95% CI = 1.09-1.36; P = 0.00053) and 10p15.1 (OR = 1.22, 95% CI = 1.08-1.38; P = 0.0015). We also found positive association signals in three regions (6q25.1, 10q26.13 and 16q12.1-q12.2) previously confirmed by fine mapping in women of African ancestry. In addition, polygenic model indicated that eight best markers in this study, compared with 22 GWAS-identified SNPs, could better predict breast cancer risk in women of African ancestry (per-allele OR = 1.21, 95% CI = 1.16-1.27; P = 9.7 × 10(-16)). Our results demonstrate that fine mapping is a powerful approach to better characterize the breast cancer risk alleles in diverse populations. Future studies and new GWAS in women of African ancestry hold promise to discover additional variants for breast cancer susceptibility with clinical implications throughout the African diaspora.

Clonal origins of relapse in ETV6-RUNX1 acute lymphoblastic leukemia.

B-cell precursor childhood acute lymphoblastic leukemia with ETV6-RUNX1 (TEL-AML1) fusion has an overall good prognosis, but relapses occur, usually after cessation of treatment and occasionally many years later. We have investigated the clonal origins of relapse by comparing the profiles of genomewide copy number alterations at presentation in 21 patients with those in matched relapse (12-119 months). We identified, in total, 159 copy number alterations at presentation and 231 at relapse (excluding Ig/TCR). Deletions of CDKN2A/B or CCNC (6q16.2-3) or both increased from 38% at presentation to 76% in relapse, suggesting that cell-cycle deregulation contributed to emergence of relapse. A novel observation was recurrent gain of chromosome 16 (2 patients at presentation, 4 at relapse) and deletion of plasmocytoma variant translocation 1 in 3 patients. The data indicate that, irrespective of time to relapse, the relapse clone was derived from either a major or minor clone at presentation. Backtracking analysis by FISH identified a minor subclone at diagnosis whose genotype matched that observed in relapse ∼ 10 years later. These data indicate subclonal diversity at diagnosis, providing a variable basis for intraclonal origins of relapse and extended periods (years) of dormancy, possibly by quiescence, for stem cells in ETV6-RUNX1(+) acute lymphoblastic leukemia.

Disruption to the FOXO-PRDM1 axis resulting from deletions of chromosome 6 in acute lymphoblastic leukaemia.

A common problem in the study of human malignancy is the elucidation of cancer driver mechanisms associated with recurrent deletion of regions containing multiple genes. Taking B-cell acute lymphoblastic leukaemia (B-ALL) and large deletions of 6q [del(6q)] as a model, we integrated analysis of functional cDNA clone tracking assays with patient genomic and transcriptomic data, to identify the transcription factors FOXO3 and PRDM1 as candidate tumour suppressor genes (TSG). Analysis of cell cycle and transcriptomic changes following overexpression of FOXO3 or PRDM1 indicated that they co-operate to promote cell cycle exit at the pre-B cell stage. FOXO1 abnormalities are absent in B-ALL, but like FOXO3, FOXO1 expression suppressed growth of TCF3::PBX1 and ETV6::RUNX1 B-ALL in-vitro. While both FOXOs induced PRDM1 and other genes contributing to late pre-B cell development, FOXO1 alone induced the key transcription factor, IRF4, and chemokine, CXCR4. CRISPR-Cas9 screening identified FOXO3 as a TSG, while FOXO1 emerged as essential for B-ALL growth. We relate this FOXO3-specific leukaemia-protective role to suppression of glycolysis based on integrated analysis of CRISPR-data and gene sets induced or suppressed by FOXO1 and FOXO3. Pan-FOXO agonist Selinexor induced the glycolysis inhibitor TXNIP and suppressed B-ALL growth at low dose (ID50 < 50 nM).

Epidermal expression of an Elovl4 transgene rescues neonatal lethality of homozygous Stargardt disease-3 mice.

Elongase of very long chain fatty acids-4 (ELOVL4) is the only mammalian enzyme known to synthesize C28-C36 fatty acids. In humans, ELOVL4 mutations cause Stargardt disease-3 (STGD3), a juvenile dominant macular degeneration. Heterozygous Stgd3 mice that carry a pathogenic mutation in the mouse Elovl4 gene demonstrate reduced levels of retinal C28-C36 acyl phosphatidylcholines (PC) and epidermal C28-C36 acylceramides. Homozygous Stgd3 mice die shortly after birth with signs of disrupted skin barrier function. In this study, we report generation of transgenic (Tg) mice with targeted Elovl4 expression driven by an epidermal-specific involucrin promoter. In homozygous Stgd3 mice, this transgene reinstates both epidermal Elovl4 expression and synthesis of two missing epidermal lipid groups: C28-C36 acylceramides and (O-linoleoyl)-omega-hydroxy C28-C36 fatty acids. Transgene expression also restores skin barrier function and rescues the neonatal lethality of homozygous Stgd3 mice. These studies establish the critical requirement for epidermal C28-C36 fatty acid synthesis for animal viability. In addition to the skin, Elovl4 is also expressed in other tissues, including the retina, brain, and testes. Thus, these mice will facilitate future studies to define the roles of C28-C36 fatty acids in the Elovl4-expressing tissues.

The HBS1L-MYB intergenic interval associated with elevated HbF levels shows characteristics of a distal regulatory region in erythroid cells.

HBS1L-MYB intergenic polymorphism (HMIP) on chromosome 6q23 is associated with elevated fetal hemoglobin levels and has pleiotropic effects on several hematologic parameters. To investigate potential regulatory activity in the region, we have measured sensitivity of the sequences to DNase I cleavage that identified 3 tissue-specific DNase I hypersensitive sites in the core intergenic interval. Chromatin immunoprecipitation with microarray (ChIP-chip) analysis showed strong histone acetylation in a defined interval of 65 kb corresponding to the core HBS1L-MYB intergenic region in primary human erythroid cells but not in non-MYB-expressing HeLa cells. ChIP-chip analysis also identified several potential cis-regulatory elements as strong GATA-1 signals that coincided with the DNase I hypersensitive sites present in MYB-expressing erythroid cells. We suggest that HMIP contains regulatory sequences that could be important in hematopoiesis by controlling MYB expression. This study provides the functional link between genetic association of HMIP with control of fetal hemoglobin and other hematologic parameters. We also present a large-scale analysis of histone acetylation as well as RNA polymerase II and GATA-1 interactions on chromosome 6q, and alpha and beta globin gene loci. The data suggest that GATA-1 regulates numerous genes of various functions on chromosome 6q.

Hetero-oligomeric interactions of an ELOVL4 mutant protein: implications in the molecular mechanism of Stargardt-3 macular dystrophy.

PURPOSE: Stargardt disease 3 (STGD3) is a juvenile macular dystrophy caused by mutations in the elongase of very long-chain fatty acids-like 4 (ELOVL4) gene, which encodes an elongase involved in the production of extremely long-chain fatty acids. The STGD3-related mutations cause production of C-terminally truncated proteins (ELOVL4ΔC). STGD3 is transmitted in an autosomal dominant manner. To date, molecular mechanisms of this pathology have been proposed based solely on the interaction between wild-type ELOVL4 and ELOVL4ΔC. However, analyses of Elovl4ΔC knockin mice revealed reduced levels of not only ELOVL4 substrates, but also of fatty acids with a broad spectrum of chain lengths. Therefore, we investigated the molecular mechanisms responsible for ELOVL4ΔC affecting the entire very long-chain fatty acid (VLCFA) elongation pathway. METHODS: The ELOVL4ΔC protein was expressed in HEK 293T cells, and its effect on elongase activities toward several acyl-CoAs were examined. We also investigated the homo- and hetero-oligomerization of ELOVL4ΔC with other elongases (ELOVL1-7) or with other enzymes involved in VLCFA elongation using coimmunoprecipitation experiments. RESULTS: We found that ELOVL4ΔC forms a homo-oligomer more strongly than wild-type ELOVL4. ELOVL4ΔC also interacts strongly with other elongases, although similar interactions for wild-type ELOVL4 were observed as only weak. In addition, ELOVL4ΔC is able to form an elongase complex by interacting with other components of the VLCFA elongation machinery, similar to wild-type ELOVL4. CONCLUSIONS: We propose that not only the ELOVL4-ELOVL4ΔC homo-oligomeric interaction, but also several hetero-oligomeric interactions, may contribute to the pathology of STGD3.

SNP rs17079281 decreases lung cancer risk through creating an YY1-binding site to suppress DCBLD1 expression.

Genome-wide association studies (GWAS) have identified numerous genetic variants that are associated with lung cancer risk, but the biological mechanisms underlying these associations remain largely unknown. Here we investigated the functional relevance of a genetic region in 6q22.2 which was identified to be associated with lung cancer risk in our previous GWAS. We performed linkage disequilibrium (LD) analysis and bioinformatic prediction to screen functional SNPs linked to a tagSNP in 6q22.2 loci, followed by two case-control studies and a meta-analysis with 4403 cases and 5336 controls to identify if these functional SNPs were associated with lung cancer risk. A novel SNP rs17079281 in the DCBLD1 promoter was identified to be associated with lung cancer risk in Chinese populations. Compared with those with C allele, patients with T allele had lower risk of adenocarcinoma (adjusted OR = 0.86; 95% CI: 0.80-0.92), but not squamous cell carcinoma (adjusted OR = 0.99; 95% CI: 0.91-1.10), and patients with the C/T or T/T genotype had lower levels of DCBLD1 expression than those with C/C genotype in lung adenocarcinoma tissues. We performed functional assays to characterize its biological relevance. The results showed that the T allele of rs17079281 had higher binding affinity to transcription factor YY1 than the C allele, which suppressed DCBLD1 expression. DCBLD1 behaved like an oncogene, promoting tumor growth by influencing cell cycle progression. These findings suggest that the functional variant rs17079281C>T decreased lung adenocarcinoma risk by creating an YY1-binding site to suppress DCBLD1 expression, which may serve as a biomarker for assessing lung cancer susceptibility.

"Reversible" myelodysplastic syndrome or ineffectual clonal haematopoiesis? - add(6p) myeloid neoplasm with a spontaneous cytogenetic remission.

Cytotoxic chemotherapy has inherent mutagenic potential and alters the bone marrow microenvironment after therapy. In some cases, this potentiates expansion of an aberrant clone and may lead to a therapy-related myeloid neoplasm if the clone overcomes selective pressure. We present the case of a 43-year-old woman diagnosed with an indolent, therapy-related myeloid neoplasm with an isolated chromosome 6p abnormality following treatment for de novo Acute Myeloid Leukaemia (AML), who manifest a sustained spontaneous cytogenetic remission two years later, possibly due to an ineffectual or non-dominant founding clone. This case reminds us to be mindful of the possibility that clonal haematopoiesis may not always equate to clinically relevant disease, even in the setting of an abnormal clonal karyotype.

Genetic mutations associated with neonatal diabetes mellitus in Omani patients.

BACKGROUND: Neonatal diabetes mellitus (NDM) is a rare disorder worldwide where diabetes is diagnosed in the first 6 months of life. However, Oman has a relatively high incidence of NDM. METHODS: In this study, we investigated the genetic etiologies underlying NDM and their prevalence in Oman. We collected a cohort of 24 NDM patients, with and without genetic diagnosis, referred to our center from 2007 to 2015. All patients without a genetic diagnosis were tested for mutations in 23 NDM-associated genes using a custom-targeted next-generation sequencing (NGS) panel and methylation analysis of the 6q24 locus. RESULTS: A genetic abnormality was detected in 15/24 (62.5%) of our Omani NDM patients. We report the detection of 6q24 methylation abnormalities and KCNJ11 mutations for the first time in Omani NDM patients. Unlike Western populations where NDM is predominantly due to mutations in the KCNJ11, ABCC8 and INS genes, NDM due to homozygous GCK gene mutations were most prevalent in Oman, having been observed in seven out of 15 NDM patients in whom we established the genetic etiology. This reflects the high degree of consanguinity which makes recessive conditions more likely. CONCLUSIONS: The results of this study are likely to impact any future strategy to introduce genetic testing for NDM disorders within the national healthcare system in Oman.

Generic features of tertiary chromatin structure as detected in natural chromosomes.

Knowledge of tertiary chromatin structure in mammalian interphase chromosomes is largely derived from artificial tandem arrays. In these model systems, light microscope images reveal fibers or beaded fibers after high-density targeting of transactivators to insertional domains spanning several megabases. These images of fibers have lent support to chromonema fiber models of tertiary structure. To assess the relevance of these studies to natural mammalian chromatin, we identified two different approximately 400-kb regions on human chromosomes 6 and 22 and then examined light microscope images of interphase tertiary chromatin structure when the regions were transcriptionally active and inactive. When transcriptionally active, these natural chromosomal regions elongated, yielding images characterized by a series of adjacent puncta or "beads", referred to hereafter as beaded images. These elongated structures required transcription for their maintenance. Thus, despite marked differences in the density and the mode of transactivation, the natural and artificial systems showed similarities, suggesting that beaded images are generic features of transcriptionally active tertiary chromatin. We show here, however, that these images do not necessarily favor chromonema fiber models but can also be explained by a radial-loop model or even a simple nucleosome affinity, random-chain model. Thus, light microscope images of tertiary structure cannot distinguish among competing models, although they do impose key constraints: chromatin must be clustered to yield beaded images and then packaged within each cluster to enable decondensation into adjacent clusters.

An integrated model for detecting significant chromatin interactions from high-resolution Hi-C data.

Here we present HiC-DC, a principled method to estimate the statistical significance (P values) of chromatin interactions from Hi-C experiments. HiC-DC uses hurdle negative binomial regression account for systematic sources of variation in Hi-C read counts-for example, distance-dependent random polymer ligation and GC content and mappability bias-and model zero inflation and overdispersion. Applied to high-resolution Hi-C data in a lymphoblastoid cell line, HiC-DC detects significant interactions at the sub-topologically associating domain level, identifying potential structural and regulatory interactions supported by CTCF binding sites, DNase accessibility, and/or active histone marks. CTCF-associated interactions are most strongly enriched in the middle genomic distance range (∼700 kb-1.5 Mb), while interactions involving actively marked DNase accessible elements are enriched both at short (<500 kb) and longer (>1.5 Mb) genomic distances. There is a striking enrichment of longer-range interactions connecting replication-dependent histone genes on chromosome 6, potentially representing the chromatin architecture at the histone locus body.

FOXO3 longevity interactome on chromosome 6.

FOXO3 has been implicated in longevity in multiple populations. By DNA sequencing in long-lived individuals, we identified all single nucleotide polymorphisms (SNPs) in FOXO3 and showed 41 were associated with longevity. Thirteen of these had predicted alterations in transcription factor binding sites. Those SNPs appeared to be in physical contact, via RNA polymerase II binding chromatin looping, with sites in the FOXO3 promoter, and likely function together as a cis-regulatory unit. The SNPs exhibited a high degree of LD in the Asian population, in which they define a specific longevity haplotype that is relatively common. The haplotype was less frequent in whites and virtually nonexistent in Africans. We identified distant contact points between FOXO3 and 46 neighboring genes, through long-range physical contacts via CCCTC-binding factor zinc finger protein (CTCF) binding sites, over a 7.3 Mb distance on chromosome 6q21. When activated by cellular stress, we visualized movement of FOXO3 toward neighboring genes. FOXO3 resides at the center of this early-replicating and highly conserved syntenic region of chromosome 6. Thus, in addition to its role as a transcription factor regulating gene expression genomewide, FOXO3 may function at the genomic level to help regulate neighboring genes by virtue of its central location in chromatin conformation via topologically associated domains. We believe that the FOXO3 'interactome' on chromosome 6 is a chromatin domain that defines an aging hub. A more thorough understanding of the functions of these neighboring genes may help elucidate the mechanisms through which FOXO3 variants promote longevity and healthy aging.

Analysis of Plasminogen Genetic Variants in Multiple Sclerosis Patients.

Multiple sclerosis (MS) is a prevalent neurological disease of complex etiology. Here, we describe the characterization of a multi-incident MS family that nominated a rare missense variant (p.G420D) in plasminogen (PLG) as a putative genetic risk factor for MS. Genotyping of PLG p.G420D (rs139071351) in 2160 MS patients, and 886 controls from Canada, identified 10 additional probands, two sporadic patients and one control with the variant. Segregation in families harboring the rs139071351 variant, identified p.G420D in 26 out of 30 family members diagnosed with MS, 14 unaffected parents, and 12 out of 30 family members not diagnosed with disease. Despite considerably reduced penetrance, linkage analysis supports cosegregation of PLG p.G420D and disease. Genotyping of PLG p.G420D in 14446 patients, and 8797 controls from Canada, France, Spain, Germany, Belgium, and Austria failed to identify significant association with disease (P = 0.117), despite an overall higher prevalence in patients (OR = 1.32; 95% CI = 0.93-1.87). To assess whether additional rare variants have an effect on MS risk, we sequenced PLG in 293 probands, and genotyped all rare variants in cases and controls. This analysis identified nine rare missense variants, and although three of them were exclusively observed in MS patients, segregation does not support pathogenicity. PLG is a plausible biological candidate for MS owing to its involvement in immune system response, blood-brain barrier permeability, and myelin degradation. Moreover, components of its activation cascade have been shown to present increased activity or expression in MS patients compared to controls; further studies are needed to clarify whether PLG is involved in MS susceptibility.

Human protein tyrosine phosphatase-sigma: alternative splicing and inhibition by bisphosphonates.

Two forms of the transmembrane human protein tyrosine phosphatase (PTP sigma), generated by alternative splicing, were identified by cDNA cloning and Northern hybridization with selective cDNA probes. The larger form of PTP sigma is expressed in various human tissues, human osteosarcoma, and rat tibia. The hPTP sigma cDNA codes for a protein of 1911 amino acid residues and is composed of a cytoplasmic region with two PTP domains and an extracellular region that can be organized into three tandem repeats of immunoglobulin-like domains and eight tandem repeats of fibronectin type III-like domains. In the brain, the major transcript of PTP sigma is an alternatively spliced mRNA, in which the coding region for the fibronectin type III-like domains number four to seven are spliced out, thus coding for a protein of 1502 amino acid residues similar to the rat PTP sigma and rat PTP-NE3. Using in situ hybridization, we assigned hPTP sigma to chromosome 6, arm 6q and band 6q15. The bacterial-expressed hPTP sigma exhibits PTPase activity that was inhibited by orthovanadate (IC50 = 0.02 microM) and by two bisphosphonates used for the treatment of bone diseases, alendronate (ALN) (IC50 = 0.5 microM) and etidronate (IC50 = 0.2 microM). In quiescent calvaria osteoblasts, micromolar concentrations of vanadate, ALN and etidronate stimulate cellular proliferation. These findings show tissue-specific alternative splicing of PTP sigma and suggest that PTPs are putative targets of bisphosphonate action.