Expressão e purificação da proteína E3 do vírus chikungunya (CHIKV)

RESUMO: O vírus chikungunya (CHIKV) é um vírus de RNA, com invólucro, da família Togaviridae, transmitido por mosquitos Aedes spp. Distribuído por largas regiões de África e Ásia, causa grandes epidemias de artrite grave. A semelhança de sintomas com outras doenças como a dengue e a malária e a persistência de IgM específicas, dificultam o diagnóstico da infeção por CHIKV. A deteção no sangue de E3, uma glicoproteína viral secretada, a incluir num ensaio imunoenzimático poderá melhorar o diagnóstico nos países onde as técnicas de biologia molecular são de difícil acesso. Para testar a utilidade de E3 num ensaio de diagnóstico, esta deverá ser expressa em quantidade, purificada e usada para produção de anticorpos específicos. Para expressar E3 numa forma solúvel, suscetível de ser purificada num único passo cromatográfico sem proteases, recorreu-se à estratégia da fusão com o domínio de ligação à quitina (CBD)-inteína (IMPACT™ System, NEB). A sequência codificadora de E3 foi amplificada a partir de RNA viral, clonada em pTYB21 e expressa em E. coli como uma proteína de fusão insolúvel de 64 kDa. A expressão a 12ºC induzida por IPTG 0,1 mM aumentou a solubilidade de CBD-inteína-E3. A aplicação de lisados celulares em colunas de quitina originou a retenção de CBD-inteína-E3 na matriz. Porém, a autoclivagem da inteína na coluna, induzida com reagentes tiol, foi pouco eficiente e mesmo a proteína E3 separada não eluiu da coluna. E3 foi ainda expressa em E. coli com uma cauda de seis histidinas (E3[His]6) por clonagem no vetor pET28b(+). Lisados celulares aplicados em colunas de níquel permitiram a eluição de uma proteína de 9 kDa, compatível com a massa molecular estimada para E3[His]6, ainda que com outros contaminantes proteicos. A identidade da proteína de 9 kDa será confirmada pela indução de anticorpos com esta preparação e reatividade daqueles com células infetadas com CHIKV.----------------ABSTRACT: Chikungunya virus (CHIKV) is an enveloped, positive strand RNA virus belonging to the family Togaviridae. Transmitted by Aedes spp mosquitoes, CHIKV causes large epidemics of severe arthritogenic disease in Africa and Asia and represents a serious threat in countries where vectors are present. Symptoms similarity with other diseases, e.g. dengue and malaria, along with CHIKV IgM persistence turns accurate CHIKV diagnosis a difficult task in low-income countries. Detection of E3, a small secreted viral glycoprotein, to be included in an immunoenzymatic test was envisaged as a possible improvement in CHIKV diagnosis. To test the diagnostic value of E3, recombinant E3 should be expressed and purified to generate antibodies. In order to express CHIKV E3 in a soluble form amenable to purification by a single step affinity chromatography, the chitin binding domain (CBD)-intein fusion strategy without proteases (IMPACT™ System, NEB) was employed. The E3 coding sequence was amplified from viral RNA, cloned in pTYB21 and expressed in E. coli ER2566 as an insoluble 64 kDa CBD-intein-E3 fusion protein. Solubility was partially achieved by lowering the expression temperature to 12ºC and the inducer (IPTG) concentration to 0.1 mM. Clarified cell lysate loaded onto a chitin column allowed ligation of the fusion protein but the intein-mediated cleavage efficiency was low and E3 failed to elute from the column as demonstrated by SDS-PAGE. E3 was further expressed with a six histidine tag, E3[His]6, employing the pET System (Novagen). E3[His]6 was expressed in E. coli Rosetta (30ºC, 0.4 mM IPTG) as a 9 kDa protein. Soluble cell extracts in 20-40 mM imidazole, applied onto a nickel column and eluted with 500 mM imidazole yielded a protein preparation enriched in the 9kDa protein. The 9 kDa will be used as antigen to generate antibodies that upon reaction with CHIKV infected cells will confirm its identity.

Coevolution of restriction factors and retroviruses

Les virus exploitent la machinerie cellulaire de l'hôte pour se répliquer. Ils doivent s'adapter pour infecter la cellule hôte de manière optimale tout en échappant à la vigilance du système de défense de l'hôte. Ainsi l'hôte et les virus se livrent à de constantes batailles évolutives. Mon travail de thèse a porté sur l'étude des signatures évolutives de facteurs de l'hôte agissant comme des 'facteurs de restriction' en bloquant la réplication rétrovirale chez les primates. Plus spécifiquement, mon travail a visé à utiliser des données évolutives pour renseigner les analyses fonctionnelles et la biologie. Nous avons étudié le facteur anti-VIH-1 nommé TRIM5a (i) chez les prosimiens pour mieux comprendre son rôle dans le contrôle d'un lentivirus endogène, (ii) dans son activité contre d'autres anciennes infections représentées par des rétrovirus endogènes humains et (iii) en tant que protéine capable de générer des mutants de la capside. Premièrement nous nous sommes intéressés à TRIM5a chez deux espèces de lémuriens dont Microcebus murinus qui porte le lentivirus endogène PSIV dans son génome depuis plusieurs millions d'années,. Nous avons observé que TRIM5a chez M. murinus a un spectre d'activité antivirale réduit à l'opposé de TRIM5a chez le Lemur catta - non porteur du PSIV endogène - qui bloque une large variété de rétrovirus dont le PSIV. De ce fait TRIM5a aurait pu contribuer à protéger certaines espèces de lémuriens vis-à-vis d'anciennes infections par le PSIV. A l'inverse du PSIV, des virus dérivés des rétrovirus endogènes humains HERV-K and HERV-H se sont révélés largement résistants à l'inhibition par TRIM5a. Ces données illustrent une absence de protection par TRIM5a face à d'autres anciennes infections rétrovirales. Puis, pour évaluer l'impact de la protéine TRIM5a humaine sur le VIH-1, nous avons testé l'effet de mutations des résidues sous sélection positive dans la capside du VIH-1 sur l'inhibition par TRIM5a. Nos résultats montrent que TRIM5a ne jouerait pas un rôle significatif dans l'évolution de la capside du VIH-1. Enfin notre travail a porté sur le facteur anti-VIH-1 SAMHD1 récemment découvert, que nous avons séquencé chez 25 espèces de primates. L'analyse évolutive des sites sous sélection positive et des expériences fonctionnelles ont permis d'identifier le domaine de SAMHD1 interagissant avec la protéine lentivirale Vpx. De même que d'autres protéines virales contrecarrent les facteurs de restriction en les menant à la dégradation, nous avons observé que Vpx induit la dégradation de SAMHD1 de manière spécifique à l'espèce. Ces découvertes contribuent à comprendre comment les facteurs de restriction et les virus co-évoluent pour se neutraliser l'un l'autre. - Viruses hijack the host cellular machinery to replicate. They adapt to infect optimally host cells while escaping host defense systems. Viruses and the host coevolve in an evolutionary struggle. My thesis work has been devoted to study the evolutionary signatures of host factors acting as restriction factors that block retroviral replication in primates. Specifically, my work aimed at using evolutionary data to inform functional analyses and biology. We studied the anti-HIV-1 factor TRIM5a (i) in prosimians to better understand its possible role in the control of an endogenous lentivirus, (ii) in its activity against other ancient infections - as represented by HERVs, and (iii) as a protein capable of generating escape mutants in the viral capsid. First, my work focused on two lemur species, one of which was the gray mouse lemur that carries the endogenous lentivirus PSIV integrated in its genome for several million years. TRIM5a from gray mouse lemur exhibited a limited antiviral spectrum as opposed to TRIM5a from ring-tailed lemur - not a host of PSIV - that is able to block diverse retroviruses notably PSIV. These results support the possible contribution of TRIM5a in protecting lemur species from ancient infection by PSIV. In contrast, chimeric viruses derived from two human endogenous retroviruses were broadly resistant to TRIM5a-mediated restriction, suggesting TRIM5a lack of activity against other types of ancient infections. To evaluate the recent impact of human TRIM5a on HIV-1 evolution, we tested whether variants at positively selected sites in the HIV-1 capsid affected the ability of human TRIM5a alleles to restrict HIV-1. Our results indicate that TRIM5a does not play a significant role in the evolution of HIV1 capsid. At last, our work concentrated on the newly discovered anti-HIV-1 restriction factor SAMHD1. We determined its coding sequence in a panel of 25 species of primates. Evolutionary analyses of positively selected sites in SAMHD1 domains and functional assays identified the domain of SAMHD1 interacting with the lentiviral protein Vpx. Similar to other viral countermeasures targeting cellular restriction factors, Vpx was responsible of the degradation of SAMHD1 orthologs in a species-specific manner. These findings contributed to understanding how restriction factors and viruses evolve to counteract each other.

Mutation screening of the glutamate cysteine ligase modifier (GCLM) gene in patients with schizophrenia.

BACKGROUND: Experimental evidences show that glutathione and its rate-limiting synthesizing enzyme, the glutamate-cysteine ligase (GCL), are involved in the pathogenesis of schizophrenia. Furthermore, genetic association has been previously reported between two single nucleotide polymorphisms lying in noncoding regions of glutamate cysteine ligase modifier (GCLM) gene, which specifies for the modifier subunit of GCL and schizophrenia. OBJECTIVE: We wanted to investigate the presence of GCLM true functional mutations, likely in linkage disequilibrium with the previously identified single nucleotide polymorphism alleles, in the same set of cases that allowed the detection of the original association signal. METHODS: We screened all the coding regions of GCLM and their intronic flanking vicinities in 353 patients with schizophrenia by direct DNA sequencing. RESULTS: Ten sequence variations were identified, five of which were not previously described. None of these DNA changes was within the GCLM coding sequence and in-silico analysis failed to indicate functional impairment induced by these variations. Furthermore, screening of normal controls and downstream statistical analyses revealed no significant relationship of any of these DNA variants with schizophrenia. CONCLUSION: It is unlikely that functional mutations in the GCLM gene could play a major role in genetic predisposition to schizophrenia and further studies will be required to assess its etiological function in the disease.

Survey of the irp2 gene among Yersinia pestis strains isolated during several plague outbreaks in northeast Brazil

The irp2 gene codes for a 190 kDa protein (HMWP2) synthesidez when highly pathogenic Yersinia are grown under conditions of iron starvation. In this work, the presence of irp2 in strains of Y. pestis isolated from different hosts during several plague outbreaks in the foci of Northeast Brazil wasstudied. For this purpose, 53 strains were spotted onto nylon filters and their DNA was hybridized with the A13 probe which is a 1 kb fragment of the irp2 coding sequence. All strains except two hybridized with the probe. However, when the initial stock culture of these two strains were analyzed, they both proved to bepositive with the A13 probe, indicating that the locus was lost after subculturein vitro but was always present in vivo. To examine the degree of conservation of the chromosomal fragment carrying irp2 among Brazilian strains, the hybridization profiles of 15 strains from different outbreaks, different hosts and different foci were compared. The hybridization profiles of these strains were all identical when their DNA was digested with either EcoRI, EcorRV or AvaII, indicatingthat the restriction sites surrounding the irp2 locus are very well conserved among Northeast Brazilian strains of Y. pestis. Altogether, these results suggest that the irp2 chromosomal region should be of prime importance for the bacteria during their multiplication in the host.

Clinical and genetic findings in a large cohort of patients with congenital myopathies due to mutations in the skeletal muscle ryanodine receptor (RYR1) gene

RYR1 mutations are the most common cause of structural congenital myopathies and may exhibit both dominant and recessive inheritance. Histopathological findings are variable and include central cores, multi-minicores, type 1 predominance/ uniformity, fibre type disproportion, increased internal nucleation and fatty and connective tissue. Until recently, diagnostic RYR1 sequencing was limited to mutational hotspots due to the large size of the gene. Since the introduction of full RYR1 sequencing in 2007 we have detected pathogenic mutations in 77 families: 39 had dominant inheritance and 38 recessive inheritance. In some cases with presumably recessive inheritance, only one heterozygous mutation inherited from an asymptomatic parent was identified. Of 28 dominant mutations, 6 were novel; 37 of the 59 recessive mutations were also novel. Dominant mutations were more frequently in recognized hotspot regions, while recessive mutations were distributed throughout the coding sequence. Dominant mutations were predominantly missense, whereas recessive mutations included many nonsense and splice mutations expected to result in reduced RyR1 protein. There was wide clinical variability in patients with both dominant and recessive inheritance. As a group, those with dominant mutations were generally more mildly affected than those with recessive inheritance, who had earlier onset and were weaker with more functional limitations. Extraocular muscle involvement was almost exclusively observed in the recessive group. Bulbar involvement was also more prominent in this group, resulting in a larger number requiring gastrostomy insertion. In conclusion, genomic sequencing of the entire RYR1 leads to the detection of many novel mutations, but may miss large genetic rearrangements in some cases. Assigning pathogenicity to novel mutations is often difficult and interpretation of genetic results in the context of clinical, histological and, increasingly, muscle MRI findings is essential.

Use of the chicken lysozyme 5' matrix attachment region to generate high producer CHO cell lines.

Scaffold or matrix attachment region (S/MAR) genetic elements have previously been proposed to insulate transgenes from repressive effects linked to their site of integration within the host cell genome. We have evaluated their use in various stable transfection settings to increase the production of recombinant proteins such as monoclonal antibodies from Chinese hamster ovary (CHO) cell lines. Using the green fluorescent protein coding sequence, we show that S/MAR elements mediate a dual effect on the population of transfected cells. First, S/MAR elements almost fully abolish the occurrence of cell clones that express little transgene that may result from transgene integration in an unfavorable chromosomal environment. Second, they increase the overall expression of the transgene over the whole range of expression levels, allowing the detection of cells with significantly higher levels of transgene expression. An optimal setting was identified as the addition of a S/MAR element both in cis (on the transgene expression vector) and in trans (co-transfected on a separate plasmid). When used to express immunoglobulins, the S/MAR element enabled cell clones with high and stable levels of expression to be isolated following the analysis of a few cell lines generated without transgene amplification procedures.

Comprehensive analysis of NRG1 common and rare variants in Hirschsprung patients.

Hirschsprung disease (HSCR, OMIM 142623) is a developmental disorder characterized by the absence of ganglion cells along variable lengths of the distal gastrointestinal tract, which results in tonic contraction of the aganglionic gut segment and functional intestinal obstruction. The RET proto-oncogene is the major gene for HSCR with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. Many other genes have been described to be associated with the pathology, as NRG1 gene (8p12), encoding neuregulin 1, which is implicated in the development of the enteric nervous system (ENS), and seems to contribute by both common and rare variants. Here we present the results of a comprehensive analysis of the NRG1 gene in the context of the disease in a series of 207 Spanish HSCR patients, by both mutational screening of its coding sequence and evaluation of 3 common tag SNPs as low penetrance susceptibility factors, finding some potentially damaging variants which we have functionally characterized. All of them were found to be associated with a significant reduction of the normal NRG1 protein levels. The fact that those mutations analyzed alter NRG1 protein would suggest that they would be related with HSCR disease not only in Chinese but also in a Caucasian population, which reinforces the implication of NRG1 gene in this pathology.

Novel (60%) and recurrent (40%) androgen receptor gene mutations in a series of 59 patients with a 46,XY disorder of sex development.

BACKGROUND Androgen receptor (AR) gene mutations are the most frequent cause of 46,XY disorders of sex development (DSD) and are associated with a variety of phenotypes, ranging from phenotypic women [complete androgen insensitivity syndrome (CAIS)] to milder degrees of undervirilization (partial form or PAIS) or men with only infertility (mild form or MAIS). OBJECTIVE The aim of the study was to characterize the contribution of the AR gene to the molecular cause of 46,XY DSD in a series of Spanish patients. SETTING We studied a series of 133 index patients with 46,XY DSD in whom gonads were differentiated as testes, with phenotypes including varying degrees of undervirilization, and in whom the AR gene was the first candidate for a molecular analysis. METHODS The AR gene was sequenced (exons 1 to 8 with intronic flanking regions) in all patients and in family members of 61% of AR-mutated gene patients. RESULTS AR gene mutations were found in 59 individuals (44.4% of index patients), of whom 46 (78%) were CAIS and 13 (22%) PAIS. Fifty-seven different mutations were found: 21.0% located in exon 1, 15.8% in exons 2 and 3, 57.9% in exons 4-8, and 5.3% intronic. Twenty-three mutations (40.4%) had been previously described and 34 (59.6%) were novel. CONCLUSIONS AR gene mutation is the most frequent cause of 46,XY DSD, with a clearly higher frequency in the complete phenotype. Mutations spread along the whole coding sequence, including exon 1. This series shows that 60% of mutations detected during the period 2002-2009 were novel.

ROQUIN/RC3H1 alterations are not found in angioimmunoblastic T-cell lymphoma.

Angioimmunoblastic T-cell Lymphoma (AITL) is one of the most frequent T-cell lymphoma entities. Follicular helper T lymphocytes (TFH) are recognized as the normal cellular counterpart of the neoplastic component. Despite a clonal T-cell feature and few described recurrent cytogenetic abnormalities, a driving oncogenic event has not been identified so far. It has been recently reported that in mice, heterozygous inactivation of Roquin/Rc3h1, a RING type E3 ubiquitine ligase, recapitulates many of the clinical, histological, and cellular features associated with human AITL. In this study we explored whether ROQUIN alterations could be an initial event in the human AITL oncogenic process. Using microarray and RT-PCR analyses, we investigated the levels of ROQUIN transcripts in TFH tumor cells purified from AITL (n = 8) and reactive tonsils (n = 12) and found similar levels of ROQUIN expression in both. Moreover, we also demonstrated that ROQUIN protein was expressed by AITL TFH (PD1+) cells. We then analysed ROQUIN coding sequence in 12 tumor cell-rich AITL samples and found no mutation in any of the samples. Finally, we analysed the expression of MiR101, a putative partner of ROQUIN involved in the modulation of ICOS expression and found similar levels of expression in tumor and reactive TFH. Altogether, this study shows that neither alteration of ROQUIN gene nor deregulation of miR101 expression is likely to be a frequent recurrent event in AITL.

Clinical and genetic findings in a large cohort of patients with ryanodine receptor 1 gene-associated myopathies.

Ryanodine receptor 1 (RYR1) mutations are a common cause of congenital myopathies associated with both dominant and recessive inheritance. Histopathological findings frequently feature central cores or multi-minicores, more rarely, type 1 predominance/uniformity, fiber-type disproportion, increased internal nucleation, and fatty and connective tissue. We describe 71 families, 35 associated with dominant RYR1 mutations and 36 with recessive inheritance. Five of the dominant mutations and 35 of the 55 recessive mutations have not been previously reported. Dominant mutations, typically missense, were frequently located in recognized mutational hotspot regions, while recessive mutations were distributed throughout the entire coding sequence. Recessive mutations included nonsense and splice mutations expected to result in reduced RyR1 protein. There was wide clinical variability. As a group, dominant mutations were associated with milder phenotypes; patients with recessive inheritance had earlier onset, more weakness, and functional limitations. Extraocular and bulbar muscle involvement was almost exclusively observed in the recessive group. In conclusion, our study reports a large number of novel RYR1 mutations and indicates that recessive variants are at least as frequent as the dominant ones. Assigning pathogenicity to novel mutations is often difficult, and interpretation of genetic results in the context of clinical, histological, and muscle magnetic resonance imaging findings is essential.

Patterns of evolution of host proteins involved in retroviral pathogenesis.

BACKGROUND: Evolutionary analysis may serve as a useful approach to identify and characterize host defense and viral proteins involved in genetic conflicts. We analyzed patterns of coding sequence evolution of genes with known (TRIM5alpha and APOBEC3G) or suspected (TRIM19/PML) roles in virus restriction, or in viral pathogenesis (PPIA, encoding Cyclophilin A), in the same set of human and non-human primate species. RESULTS AND CONCLUSION: This analysis revealed previously unidentified clusters of positively selected sites in APOBEC3G and TRIM5alpha that may delineate new virus-interaction domains. In contrast, our evolutionary analyses suggest that PPIA is not under diversifying selection in primates, consistent with the interaction of Cyclophilin A being limited to the HIV-1M/SIVcpz lineage. The strong sequence conservation of the TRIM19/PML sequences among primates suggests that this gene does not play a role in antiretroviral defense.

PRPF31 alternative splicing and expression in human retina

PURPOSE: To provide a mechanistic link between mutations in PRPF31, and essential and ubiquitously expressed gene, and retinitis pigmentosa, a disorder restricted to the eye. METHODS: We investigated the existence of retina-specific PRPF31 isoforms and the expression of this gene in human retina and other tissues, as well as in cultured human cell lines. PRPF31 transcripts were examined by RT-PCR, quantitative PCR, cloning and sequencing. RESULTS: Database searching revealed the presence of a retina-specific PRPF31 isoform in mouse. However, this isoform could not be experimentally identified in transcripts from human retina or from a human whole eye. Nevertheless, four different PRPF31 isoforms, that were common to all analyzed tissues and cell lines, were isolated. Three of these harbored the full-length PRPF31 coding sequence, whereas the fourth was very short and probably non-coding. The amount of PRPF31 mRNA was previously found to be lower in patients with mutations in this gene than in healthy individuals, making it likely that retinal cells are more sensitive to variation in PRPF31 expression. However, quantitative PCR experiments revealed that PRPF31 mRNA levels in human retina were comparable to those detected in other tissues. CONCLUSIONS: Our results show that the retina-restricted phenotype caused by PRPF31 mutations cannot be explained by the presence of tissue-specific isoforms, or by differential expression of PRPF31 in the retina. As a consequence, the etiology of PRPF31-associated retinitis pigmentosa likely relies on other, probably more subtle molecular mechanisms.

GENCODE: producing a reference annotation for ENCODE

Background: The GENCODE consortium was formed to identify and map all protein-coding genes within the ENCODE regions. This was achieved by a combination of initial manualannotation by the HAVANA team, experimental validation by the GENCODE consortium and a refinement of the annotation based on these experimental results.Results: The GENCODE gene features are divided into eight different categories of which onlythe first two (known and novel coding sequence) are confidently predicted to be protein-codinggenes. 5’ rapid amplification of cDNA ends (RACE) and RT-PCR were used to experimentallyverify the initial annotation. Of the 420 coding loci tested, 229 RACE products have beensequenced. They supported 5’ extensions of 30 loci and new splice variants in 50 loci. In addition,46 loci without evidence for a coding sequence were validated, consisting of 31 novel and 15putative transcripts. We assessed the comprehensiveness of the GENCODE annotation byattempting to validate all the predicted exon boundaries outside the GENCODE annotation. Outof 1,215 tested in a subset of the ENCODE regions, 14 novel exon pairs were validated, only twoof them in intergenic regions.Conclusions: In total, 487 loci, of which 434 are coding, have been annotated as part of theGENCODE reference set available from the UCSC browser. Comparison of GENCODEannotation with RefSeq and ENSEMBL show only 40% of GENCODE exons are contained withinthe two sets, which is a reflection of the high number of alternative splice forms with uniqueexons annotated. Over 50% of coding loci have been experimentally verified by 5’ RACE forEGASP and the GENCODE collaboration is continuing to refine its annotation of 1% humangenome with the aid of experimental validation.

SelenoDB 1.0 : a database of selenoprotein genes, proteins and SECIS elements

Selenoproteins are a diverse group of proteinsusually misidentified and misannotated in sequencedatabases. The presence of an in-frame UGA (stop)codon in the coding sequence of selenoproteingenes precludes their identification and correctannotation. The in-frame UGA codons are recodedto cotranslationally incorporate selenocysteine,a rare selenium-containing amino acid. The developmentof ad hoc experimental and, more recently,computational approaches have allowed the efficientidentification and characterization of theselenoproteomes of a growing number of species.Today, dozens of selenoprotein families have beendescribed and more are being discovered in recentlysequenced species, but the correct genomic annotationis not available for the majority of thesegenes. SelenoDB is a long-term project that aims toprovide, through the collaborative effort of experimentaland computational researchers, automaticand manually curated annotations of selenoproteingenes, proteins and SECIS elements. Version 1.0 ofthe database includes an initial set of eukaryoticgenomic annotations, with special emphasis on thehuman selenoproteome, for immediate inspectionby selenium researchers or incorporation into moregeneral databases. SelenoDB is freely available athttp://www.selenodb.org.

Loss-of-Function Mutations in PTPN11 Cause Metachondromatosis, but Not Ollier Disease or Maffucci Syndrome.

Metachondromatosis (MC) is a rare, autosomal dominant, incompletely penetrant combined exostosis and enchondromatosis tumor syndrome. MC is clinically distinct from other multiple exostosis or multiple enchondromatosis syndromes and is unlinked to EXT1 and EXT2, the genes responsible for autosomal dominant multiple osteochondromas (MO). To identify a gene for MC, we performed linkage analysis with high-density SNP arrays in a single family, used a targeted array to capture exons and promoter sequences from the linked interval in 16 participants from 11 MC families, and sequenced the captured DNA using high-throughput parallel sequencing technologies. DNA capture and parallel sequencing identified heterozygous putative loss-of-function mutations in PTPN11 in 4 of the 11 families. Sanger sequence analysis of PTPN11 coding regions in a total of 17 MC families identified mutations in 10 of them (5 frameshift, 2 nonsense, and 3 splice-site mutations). Copy number analysis of sequencing reads from a second targeted capture that included the entire PTPN11 gene identified an additional family with a 15 kb deletion spanning exon 7 of PTPN11. Microdissected MC lesions from two patients with PTPN11 mutations demonstrated loss-of-heterozygosity for the wild-type allele. We next sequenced PTPN11 in DNA samples from 54 patients with the multiple enchondromatosis disorders Ollier disease or Maffucci syndrome, but found no coding sequence PTPN11 mutations. We conclude that heterozygous loss-of-function mutations in PTPN11 are a frequent cause of MC, that lesions in patients with MC appear to arise following a "second hit," that MC may be locus heterogeneous since 1 familial and 5 sporadically occurring cases lacked obvious disease-causing PTPN11 mutations, and that PTPN11 mutations are not a common cause of Ollier disease or Maffucci syndrome.