Phenotypic Mutation 'candle' (pdf version)
Allelecandle
Mutation Type missense
Chromosome6
Coordinate124,705,382 bp (GRCm39)
Base Change T ⇒ C (forward strand)
Gene Ptpn6
Gene Name protein tyrosine phosphatase, non-receptor type 6
Synonym(s) Hcph, SHP-1, hcp, Ptp1C
Chromosomal Location 124,697,670-124,715,672 bp (-) (GRCm39)
MGI Phenotype FUNCTION: [Summary is not available for the mouse gene. This summary is for the human ortholog.] The protein encoded by this gene is a member of the protein tyrosine phosphatase (PTP) family. PTPs are known to be signaling molecules that regulate a variety of cellular processes including cell growth, differentiation, mitotic cycle, and oncogenic transformation. N-terminal part of this PTP contains two tandem Src homolog (SH2) domains, which act as protein phospho-tyrosine binding domains, and mediate the interaction of this PTP with its substrates. This PTP is expressed primarily in hematopoietic cells, and functions as an important regulator of multiple signaling pathways in hematopoietic cells. This PTP has been shown to interact with, and dephosphorylate a wide spectrum of phospho-proteins involved in hematopoietic cell signaling. Multiple alternatively spliced variants of this gene, which encode distinct isoforms, have been reported. [provided by RefSeq, Jul 2008]
PHENOTYPE: Homozygous mutants are immunodeficient and autoimmune and exhibit neutrophilic skin lesions that disrupt hair follicles and give the motheaten appearance. Alleles vary in severity, with death occurring at 6-9 weeks postnatally due to severe pneumonitis. [provided by MGI curators]
Accession Number

NCBI RefSeq: NM_013545 (variant 1), NM_001077705 (variant 2); MGI:96055

MappedYes 
Amino Acid Change Aspartic acid changed to Glycine
Institutional SourceBeutler Lab
Gene Model predicted gene model for protein(s): [ENSMUSP00000004377] [ENSMUSP00000108103] [ENSMUSP00000129124] [ENSMUSP00000133991] [ENSMUSP00000133747] [ENSMUSP00000133429]
AlphaFold P29351
SMART Domains Protein: ENSMUSP00000004377
Gene: ENSMUSG00000004266
AA Change: D224G

DomainStartEndE-ValueType
SH2 4 87 1.43e-28 SMART
SH2 110 202 1.45e-29 SMART
PTPc 245 519 7.51e-131 SMART
low complexity region 571 581 N/A INTRINSIC
Predicted Effect probably benign

PolyPhen 2 Score 0.028 (Sensitivity: 0.95; Specificity: 0.81)
(Using ENSMUST00000004377)
SMART Domains Protein: ENSMUSP00000108103
Gene: ENSMUSG00000004266
AA Change: D222G

DomainStartEndE-ValueType
SH2 2 85 4.05e-28 SMART
SH2 108 200 1.45e-29 SMART
PTPc 243 517 7.51e-131 SMART
low complexity region 569 579 N/A INTRINSIC
Predicted Effect probably benign

PolyPhen 2 Score 0.011 (Sensitivity: 0.96; Specificity: 0.78)
(Using ENSMUST00000112484)
SMART Domains Protein: ENSMUSP00000129124
Gene: ENSMUSG00000004266
AA Change: D224G

DomainStartEndE-ValueType
SH2 4 87 1.43e-28 SMART
SH2 110 202 1.45e-29 SMART
PTPc 245 519 7.51e-131 SMART
low complexity region 571 581 N/A INTRINSIC
Predicted Effect probably benign

PolyPhen 2 Score 0.028 (Sensitivity: 0.95; Specificity: 0.81)
(Using ENSMUST00000171549)
Predicted Effect probably benign
SMART Domains Protein: ENSMUSP00000133991
Gene: ENSMUSG00000004266
AA Change: D183G

DomainStartEndE-ValueType
Pfam:SH2 1 40 3.5e-6 PFAM
SH2 69 161 1.45e-29 SMART
PTPc 204 478 7.51e-131 SMART
low complexity region 530 540 N/A INTRINSIC
Predicted Effect probably damaging

PolyPhen 2 Score 1.000 (Sensitivity: 0.00; Specificity: 1.00)
(Using ENSMUST00000174265)
SMART Domains Protein: ENSMUSP00000133747
Gene: ENSMUSG00000004266

DomainStartEndE-ValueType
SH2 2 64 2.35e-6 SMART
Predicted Effect probably benign
Predicted Effect probably benign
Meta Mutation Damage Score 0.1372 question?
Is this an essential gene? Possibly nonessential (E-score: 0.481) question?
Phenotypic Category Autosomal Recessive
Candidate Explorer Status loading ...
Single pedigree
Linkage Analysis Data
Penetrance  
Alleles Listed at MGI

All Mutations and Alleles(16) : Chemically induced (ENU)(4) Gene trapped(3) Spontaneous (4) Targeted(5)

Lab Alleles
AlleleSourceChrCoordTypePredicted EffectPPH Score
IGL00710:Ptpn6 APN 6 124709319 splice site probably null
IGL01490:Ptpn6 APN 6 124705307 missense probably damaging 1.00
IGL01865:Ptpn6 APN 6 124709428 missense probably damaging 1.00
IGL02017:Ptpn6 APN 6 124709449 missense probably damaging 0.98
IGL02272:Ptpn6 APN 6 124698171 missense probably damaging 0.99
IGL02276:Ptpn6 APN 6 124705828 missense probably null 1.00
IGL02556:Ptpn6 APN 6 124705623 missense probably benign 0.00
caterpillar UTSW 6 124701947 missense probably benign
farfalla_notturna UTSW 6 124709398 missense probably damaging 1.00
Flutterby UTSW 6 124698821 missense possibly damaging 0.89
Hawk UTSW 6 124705748 missense probably damaging 1.00
Lepidopteran UTSW 6 124705135 missense probably damaging 1.00
Malachite UTSW 6 124705614 missense possibly damaging 0.84
Moth UTSW 6 124705135 missense possibly damaging 0.89
Naphthalene UTSW 6 124698752 missense probably benign 0.42
spin UTSW 6 124705522 missense probably damaging 1.00
spin2 UTSW 6 124709332 missense probably damaging 1.00
Vermeil UTSW 6 124709913 missense probably benign 0.10
R0183:Ptpn6 UTSW 6 124705914 missense probably damaging 1.00
R0254:Ptpn6 UTSW 6 124705113 missense probably damaging 1.00
R0636:Ptpn6 UTSW 6 124702242 missense probably benign
R0835:Ptpn6 UTSW 6 124704499 critical splice acceptor site probably null
R1383:Ptpn6 UTSW 6 124698856 missense probably damaging 1.00
R1638:Ptpn6 UTSW 6 124698148 missense probably benign
R1900:Ptpn6 UTSW 6 124705896 missense probably benign 0.15
R2047:Ptpn6 UTSW 6 124698752 missense probably benign 0.42
R2143:Ptpn6 UTSW 6 124701947 missense probably benign 0.01
R3907:Ptpn6 UTSW 6 124702239 missense possibly damaging 0.86
R4082:Ptpn6 UTSW 6 124705382 missense probably damaging 1.00
R4382:Ptpn6 UTSW 6 124704361 missense possibly damaging 0.86
R5319:Ptpn6 UTSW 6 124709913 missense probably benign 0.10
R5807:Ptpn6 UTSW 6 124701947 missense probably benign
R5878:Ptpn6 UTSW 6 124705748 missense probably damaging 1.00
R6056:Ptpn6 UTSW 6 124709398 missense probably damaging 1.00
R6374:Ptpn6 UTSW 6 124709532 splice site probably null
R7238:Ptpn6 UTSW 6 124698821 missense possibly damaging 0.89
R7381:Ptpn6 UTSW 6 124705135 missense probably damaging 1.00
R7935:Ptpn6 UTSW 6 124709425 missense possibly damaging 0.93
R8297:Ptpn6 UTSW 6 124705614 missense possibly damaging 0.84
R8863:Ptpn6 UTSW 6 124709309 missense probably damaging 1.00
R9160:Ptpn6 UTSW 6 124705135 missense possibly damaging 0.89
R9176:Ptpn6 UTSW 6 124702249 missense probably benign
R9448:Ptpn6 UTSW 6 124709771 missense probably damaging 1.00
R9594:Ptpn6 UTSW 6 124704728 missense probably benign 0.04
R9756:Ptpn6 UTSW 6 124705592 missense probably damaging 1.00
Z1176:Ptpn6 UTSW 6 124702039 nonsense probably null
Mode of Inheritance Autosomal Recessive
Local Stock
Repository
Last Updated 2019-09-04 9:43 PM by Diantha La Vine
Record Created 2016-03-04 6:05 AM by Bruce Beutler
Record Posted 2016-06-02
Phenotypic Description

Figure 1. Candle mice exhibit decreased frequencies of peripheral B cells. Flow cytometric analysis of peripheral blood was utilized to determine B cell frequency. Normalized data are shown. Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.

Figure 2. Candle mice exhibit decreased frequencies of peripheral B2 cells. Flow cytometric analysis of peripheral blood was utilized to determine B2 cell frequency. Normalized data are shown. Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 3. Candle mice exhibit a decreased percentage of peripheral IgD+ B cells. Flow cytometric analysis of peripheral blood was utilized to determine IgD+ B cell percentage. Normalized data are shown. Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 4. Candle mice exhibit decreased frequencies of peripheral IgM+ B cells. Flow cytometric analysis of peripheral blood was utilized to determine IgM+ B cell frequency. Normalized data are shown. Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 5. Candle mice exhibit increased frequencies of peripheral B1 cells. Flow cytometric analysis of peripheral blood was utilized to determine B1 cell frequency. Normalized data are shown. Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 6. Candle mice exhibit increased frequencies of peripheral B1a cells. Flow cytometric analysis of peripheral blood was utilized to determine B1a cell frequency. Normalized data are shown. Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 7. Candle mice exhibit increased frequencies of peripheral B1a cells in B cells. Flow cytometric analysis of peripheral blood was utilized to determine B1a cell frequency. Normalized data are shown. Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 8. Candle mice exhibit increased frequencies of peripheral B1b cells. Flow cytometric analysis of peripheral blood was utilized to determine B1b cell frequency. Normalized data are shown. Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 9. Candle mice exhibit increased frequencies of peripheral CD11c+ dendritic cells. Flow cytometric analysis of peripheral blood was utilized to determine CD11c+ dendritic cell frequency. Normalized data are shown. Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 10. Candle mice exhibit increased frequencies of peripheral macrophages. Flow cytometric analysis of peripheral blood was utilized to determine macrophage frequency. Normalized data are shown. Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 11. Candle mice exhibit increased CD44 mean fluorescence intensity on peripheral CD8+ T cells. Flow cytometric analysis of peripheral blood was utilized to determine CD8+ T MFI. Normalized data are shown. Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.

The candle phenotype was identified among G3 mice of the pedigree R4082, some of which showed a decrease in the frequency of B cells (Figure 1), B2 cells (Figure 2), the percentage of IgD+ B cells (Figure 3), and the frequency of IgM+ B cells (Figure 4) with a concomitant increase in the frequency of B1 cells (Figure 5), B1a cells (Figure 6), B1a cells in B cells (Figure 7), B1b cells (Figure 8), CD11c+ dendritic cells (Figure 9), and macrophages (Figure 10) as well as an increase in the CD44 mean fluorescence intensity on CD8+ T cells (Figure 11), all in the peripheral blood.

Nature of Mutation

Figure 12. Linkage mapping of the increased frequency of peripheral blood IgM+ B cells using a recessive model of inheritance. Manhattan plot shows -log10 P values (Y-axis) plotted against the chromosome positions of 68 mutations (X-axis) identified in the G1 male of pedigree R4082. Normalized phenotype data are shown for single locus linkage analysis with consideration of G2 dam identity.  Horizontal pink and red lines represent thresholds of P = 0.05, and the threshold for P = 0.05 after applying Bonferroni correction, respectively.

Whole exome HiSeq sequencing of the G1 grandsire identified 68 mutations. All of the above anomalies were linked by continuous variable mapping to a mutation in Ptpn6:  an A to G transition at base pair 124,728,419 (v38) on chromosome 6, or base pair 10,291 in the GenBank genomic region NC_000072 encoding Ptpn6.  The strongest association was found with a recessive model of linkage to the normalized frequency of peripheral blood IgM+ B cells, wherein seven variant homozygotes departed phenotypically from 11 homozygous reference mice and 17 heterozygous mice with a P value of 1.603 x 10-8 (Figure 12).  A substantial semidominant effect was observed in most of the assays but the mutation is preponderantly recessive, and in no assay was a purely dominant effect observed. 

The mutation corresponds to residue 866 in the mRNA sequence NM_013545 within exon 6 of 16 total exons and residue 868 in the mRNA sequence NM_001077705 within exon 6 of 16 total exons.

 

10275 CGGGTAAACGCAGCTGACATTGAGAATCGGGTC

211   -R--V--N--A--A--D--I--E--N--R--V- (variant 1)

219   -R--V--N--A--A--D--I--E--N--R--V- (variant 2)

Genomic numbering corresponds to NC_000072. The mutated nucleotide is indicated in red.  The mutation results in an aspartic acid (D) to glycine (G) substitution at position 219 (D219G) in variant 1 and D224G in variant 2 of the SHP1 protein.

Illustration of Mutations in
Gene & Protein
Protein Prediction

Figure 13. Domain structure of SHP1. The position of the candle mutation is indicated. Other mutations found in SHP1 are noted in red. Click on the mutations for more specific information.

Ptpn6 encodes SHP1, a Src-homology 2 (SH2) domain-containing cytoplasmic protein tyrosine phosphatase. SHP1 has 4 isoforms. Three isoforms of SHP1 contain variations in their N-termini; the fourth isoform is a longer form with an extended C-terminus. SHP1 contains two tandem N-terminal SH2 domains (residues 1-108 and 116-208), a central catalytic domain (residues 270-532), and a C-terminal tail (Figure 13) [(1), discussed in (2)]. The C-terminal tail contains multiple sites for tyrosine and serine phosphorylation. The candle mutation results in an aspartic acid to glycine amino acid change at amino acid 219 within the region between the C-terminal SH2 domain and the phosphatase domain. The candle mutation is predicted to affect all of the SHP1 isoforms.

For more information on Ptpn6, please see the record for spin.

Putative Mechanism

The phenotype of the candle mice suggests decreased SHP1 function, and optimal signaling through SHP1 requires an intact connecting region between the SH2 domain and the catalytic domain. The phenotype of candle animals is less severe than the spin (3) and spin2 phenotypes as foot lesions were not observed in the candle mice. The foot lesions observed in the spin mice were associated with the development of splenomegaly and an increased number of erythroid and myeloid cells in the spleen, as well as a reduction in mature B cell numbers in the peripheral blood, spleen and bone marrow. Spontaneously occurring motheaten (me) and viable motheaten (me-v) mutants are immunodeficient and exhibit multiple defects stemming from increased inflammation, including alopecia, glomerulonephritis, dermatitis, inflammation of the paws, and interstitial pneumonitis which ultimately causes death by 3 and 9 weeks of age in Ptpn6me/me and Ptpn6me-v/me-v mice, respectively. The phenotype of the candle animals is significantly less severe than the Ptpn6me-v allele, which encodes a SHP1 protein with approximately 20% catalytic activity (4).  It is likely that the SHP1 protein encoded by the candle allele retains catalytic activity that is greater than 20% of normal, and greater than the SHP1 protein encoded by the spin and spin2 alleles.

Primers PCR Primer
candle_pcr_F: CTCCTGCTTTTGTAGACTCTGAAATG
candle_pcr_R: ACCTGGTGGAGCACTTCAAG

Sequencing Primer
candle_seq_F: TCAGAGAACTGAGTCCCCCG
candle_seq_R: CTGGTGGAGCACTTCAAGAAGAC
Genotyping

PCR program

1) 94°C 2:00
2) 94°C 0:30
3) 55°C 0:30
4) 72°C 1:00
5) repeat steps (2-4) 40x
6) 72°C 10:00
7) 4°C hold


The following sequence of 425 nucleotides is amplified (chromosome 6, - strand):


1   acctggtgga gcacttcaag aagacaggga ttgaggaggc ctcgggtgcc tttgtctacc
61  tgcggcaggt gaggggccga catgcttgcc tcttcccctg agccatcgga gatgtggctt
121 tctgaggtct ctcggccgct gacttctcgt cctccctcac ctcagcctta ctacgctact
181 cgggtaaacg cagctgacat tgagaatcgg gtcttggaac tgaacaagaa gcaggagtcg
241 gaggacacag ccaaggctgg cttctgggag gagtttgagg tgtgtggtgg ggacagcagg
301 gctgggatgg atgaggtcac cttgccttct ttatatgcct ctgacctatg gcctctgcga
361 ccgtggatgc aaccttgggc gggggactca gttctctgac atttcagagt ctacaaaagc
421 aggag


Primer binding sites are underlined and the sequencing primers are highlighted; the mutated nucleotide is shown in red.

References
Science Writers Anne Murray
Illustrators Peter Jurek
AuthorsXue Zhong, Ming Zeng, and Bruce Beutler