rs267608240

Variant summary

Our verdict is Uncertain significance. Variant got 5 ACMG points: 5P and 0B. PM2PP3_ModeratePP5

The NM_000287.4(PEX6):​c.2362G>A​(p.Val788Met) variant causes a missense, splice region change involving the alteration of a conserved nucleotide. The variant allele was found at a frequency of 0.00000744 in 1,613,844 control chromosomes in the GnomAD database, with no homozygous occurrence. In-silico tool predicts a pathogenic outcome for this variant. 3/3 splice prediction tools predicting alterations to normal splicing. Variant has been reported in ClinVar as Conflicting classifications of pathogenicity (no stars). Synonymous variant affecting the same amino acid position (i.e. V788V) has been classified as Benign.

Frequency

Genomes: 𝑓 0.0000066 ( 0 hom., cov: 32)
Exomes 𝑓: 0.0000075 ( 0 hom. )

Consequence

PEX6
NM_000287.4 missense, splice_region

Scores

11
6
2
Splicing: ADA: 0.9999
2

Clinical Significance

Conflicting classifications of pathogenicity criteria provided, conflicting classifications P:6U:1

Conservation

PhyloP100: 7.68
Variant links:
Genes affected
PEX6 (HGNC:8859): (peroxisomal biogenesis factor 6) This gene encodes a member of the AAA (ATPases associated with diverse cellular activities) family of ATPases. This member is a predominantly cytoplasmic protein, which plays a direct role in peroxisomal protein import and is required for PTS1 (peroxisomal targeting signal 1, a C-terminal tripeptide of the sequence ser-lys-leu) receptor activity. Mutations in this gene cause peroxisome biogenesis disorders of complementation group 4 and complementation group 6. Several transcript variants encoding different isoforms have been found for this gene. [provided by RefSeq, Oct 2015]

Genome browser will be placed here

ACMG classification

Classification made for transcript

Verdict is Uncertain_significance. Variant got 5 ACMG points.

PM2
Very rare variant in population databases, with high coverage;
PP3
MetaRNN computational evidence supports a deleterious effect, 0.933
PP5
Variant 6-42966044-C-T is Pathogenic according to our data. Variant chr6-42966044-C-T is described in ClinVar as [Conflicting_classifications_of_pathogenicity]. Clinvar id is 556244.We mark this variant Likely_pathogenic, oryginal submissions are: {Uncertain_significance=1, Pathogenic=5, Likely_pathogenic=1}. Variant chr6-42966044-C-T is described in Lovd as [Pathogenic]. Variant chr6-42966044-C-T is described in Lovd as [Pathogenic].

Transcripts

RefSeq

Gene Transcript HGVSc HGVSp Effect Exon rank MANE Protein UniProt
PEX6NM_000287.4 linkc.2362G>A p.Val788Met missense_variant, splice_region_variant Exon 12 of 17 ENST00000304611.13 NP_000278.3 Q13608-1A0A024RD09

Ensembl

Gene Transcript HGVSc HGVSp Effect Exon rank TSL MANE Protein Appris UniProt
PEX6ENST00000304611.13 linkc.2362G>A p.Val788Met missense_variant, splice_region_variant Exon 12 of 17 1 NM_000287.4 ENSP00000303511.8 Q13608-1
PEX6ENST00000244546.4 linkc.2116-255G>A intron_variant Intron 10 of 14 1 ENSP00000244546.4 Q13608-2

Frequencies

GnomAD3 genomes
AF:
0.00000657
AC:
1
AN:
152192
Hom.:
0
Cov.:
32
show subpopulations
Gnomad AFR
AF:
0.00
Gnomad AMI
AF:
0.00
Gnomad AMR
AF:
0.00
Gnomad ASJ
AF:
0.00
Gnomad EAS
AF:
0.00
Gnomad SAS
AF:
0.00
Gnomad FIN
AF:
0.00
Gnomad MID
AF:
0.00
Gnomad NFE
AF:
0.0000147
Gnomad OTH
AF:
0.00
GnomAD3 exomes
AF:
0.0000159
AC:
4
AN:
251456
Hom.:
0
AF XY:
0.0000221
AC XY:
3
AN XY:
135902
show subpopulations
Gnomad AFR exome
AF:
0.00
Gnomad AMR exome
AF:
0.00
Gnomad ASJ exome
AF:
0.00
Gnomad EAS exome
AF:
0.00
Gnomad SAS exome
AF:
0.00
Gnomad FIN exome
AF:
0.00
Gnomad NFE exome
AF:
0.0000352
Gnomad OTH exome
AF:
0.00
GnomAD4 exome
AF:
0.00000753
AC:
11
AN:
1461652
Hom.:
0
Cov.:
34
AF XY:
0.0000124
AC XY:
9
AN XY:
727130
show subpopulations
Gnomad4 AFR exome
AF:
0.00
Gnomad4 AMR exome
AF:
0.00
Gnomad4 ASJ exome
AF:
0.00
Gnomad4 EAS exome
AF:
0.00
Gnomad4 SAS exome
AF:
0.00
Gnomad4 FIN exome
AF:
0.00
Gnomad4 NFE exome
AF:
0.00000989
Gnomad4 OTH exome
AF:
0.00
GnomAD4 genome
AF:
0.00000657
AC:
1
AN:
152192
Hom.:
0
Cov.:
32
AF XY:
0.0000135
AC XY:
1
AN XY:
74334
show subpopulations
Gnomad4 AFR
AF:
0.00
Gnomad4 AMR
AF:
0.00
Gnomad4 ASJ
AF:
0.00
Gnomad4 EAS
AF:
0.00
Gnomad4 SAS
AF:
0.00
Gnomad4 FIN
AF:
0.00
Gnomad4 NFE
AF:
0.0000147
Gnomad4 OTH
AF:
0.00
Alfa
AF:
0.0000868
Hom.:
0
Bravo
AF:
0.0000113
EpiCase
AF:
0.000109
EpiControl
AF:
0.00

ClinVar

Significance: Conflicting classifications of pathogenicity
Submissions summary: Pathogenic:6Uncertain:1
Revision: criteria provided, conflicting classifications
LINK: link

Submissions by phenotype

not provided Pathogenic:2
Pathogenic, criteria provided, single submitterclinical testingGeneDxOct 12, 2021Published functional studies demonstrate an impact on splicing by causing skipping of exon 11 and 12 and an insertion of 21 nucleotides, therefore demonstrating a damaging effect (Matsumoto et al., 2001); Not observed at significant frequency in large population cohorts (gnomAD); In silico analysis supports that this missense variant has a deleterious effect on protein structure/function; This variant is associated with the following publications: (PMID: 19877282, 15542397, 26700162, 19142205, 11355018) -
Pathogenic, criteria provided, single submitterclinical testingCeGaT Center for Human Genetics TuebingenSep 01, 2020- -
Peroxisome biogenesis disorder Pathogenic:1Uncertain:1
Pathogenic, criteria provided, single submitterclinical testingWomen's Health and Genetics/Laboratory Corporation of America, LabCorpMay 31, 2019Variant summary: PEX6 c.2362G>A (p.Val788Met) results in a conservative amino acid change located in the ATPase, AAA-type, core domain of the encoded protein sequence. Four of five in-silico tools predict a damaging effect of the variant on protein function. The variant is the last exonic nucleotide at an exon-intron junction, suggesting it may affect splicing. At least one publication, Matsumoto_2001, reports experimental evidence that this variant affects mRNA splicing. The variant allele was found at a frequency of 1.6e-05 in 251656 control chromosomes (gnomAD). The variant, c.2362G>A, has been reported in the literature in multiple individuals affected with Zellweger Syndrome, both as a homozygous and compound heterozygous allele. (Steinberg_2004, Krause_2009, Ebberink_2010, Lusebrink_2016). These data indicate that the variant is very likely to be associated with disease. One publication reports experimental evidence evaluating an impact on protein function, however, does not allow convincing conclusions about the variant effect (Lusebrink_2016). One ClinVar submission from other clinical diagnostic laboratories (evaluation after 2014) cites the variant as likely pathogenic. Based on the evidence outlined above, the variant was classified as pathogenic. -
Uncertain significance, criteria provided, single submitterclinical testingLabcorp Genetics (formerly Invitae), LabcorpSep 02, 2021This sequence change replaces valine with methionine at codon 788 of the PEX6 protein (p.Val788Met). The valine residue is highly conserved and there is a small physicochemical difference between valine and methionine. RNA analysis indicates that this missense change induces altered splicing and may result in an absent or disrupted protein product. This variant is not present in population databases (ExAC no frequency). This missense change has been observed in individual(s) with clinical features of neonatal adrenoleukodystrophy or Zellweger syndrome (PMID: 3515938, 19142205, 19877282, 26700162). This variant is also known as PEX6del2095–2362/21ins. ClinVar contains an entry for this variant (Variation ID: 556244). Algorithms developed to predict the effect of missense changes on protein structure and function are either unavailable or do not agree on the potential impact of this missense change (SIFT: "Deleterious"; PolyPhen-2: "Probably Damaging"; Align-GVGD: "Class C0"). Variants that disrupt the consensus splice site are a relatively common cause of aberrant splicing (PMID: 17576681, 9536098). Studies have shown that this missense change results in skipping of exons 11 and 12 and an insertion of 21 nucleotides, which introduces a frameshift and introduces a premature termination codon (PMID: 11355018). The resulting mRNA is expected to undergo nonsense-mediated decay. In summary, the available evidence is currently insufficient to determine the role of this variant in disease. Therefore, it has been classified as a Variant of Uncertain Significance. -
Peroxisome biogenesis disorder 4A (Zellweger);C3553937:Peroxisome biogenesis disorder 4B Pathogenic:1
Likely pathogenic, criteria provided, single submitterclinical testingCounsylMar 30, 2018- -
Peroxisome biogenesis disorder 4A (Zellweger) Pathogenic:1
Pathogenic, criteria provided, single submitterclinical testingVictorian Clinical Genetics Services, Murdoch Childrens Research InstituteJul 17, 2023Based on the classification scheme VCGS_Germline_v1.3.4, this variant is classified as Pathogenic. Following criteria are met: 0102 - Loss of function is a known mechanism of disease in this gene and is associated with peroxisome biogenesis disorder 4B (MIM#614863), Heimler syndrome 2 (MIM#616617) and peroxisome biogenesis disorder 4A (Zellweger) (MIM#614862). (I) 0108 - This gene is associated with both recessive and dominant disease. A recurrent variant p.(Arg860Trp) is associated with autosomal dominant peroxisome biogenesis disorder 4B (PMID: 29220678). In addition, Heimler syndrome 2 is milder due to one allele encoding for a protein with residual function (OMIM). (I) 0200 - Variant is predicted to result in a missense amino acid change from valine to methionine. It was noted that ClinVar entries cited Matsumoto 2001 (PMID: 11355018), describing a splicing defect caused by this variant. However, upon asssessing this manuscript, this assay may have amplified a naturally occurring transcript (NR_133009.2, UCSC), and therefore, the information from this publication was treated with caution during variant classification. (I) 0251 - This variant is heterozygous. (I) 0304 - Variant is present in gnomAD <0.01 for a recessive condition (v2: 5 heterozygotes, 0 homozygotes). (SP) 0502 - Missense variant with conflicting in silico predictions and very high conservation. In addition, this variant affects the last base pair of an exon. Abnormal splicing is predicted by in silico tools and the affected nucleotide is highly conserved. (SP) 0600 - Variant is located in the annotated AAA ATPase domain (DECIPHER). (I) 0705 - No comparable missense variants have previous evidence for pathogenicity. (I) 0801 - This variant has strong previous evidence of pathogenicity in unrelated individuals. It has been described as compound heterozygous with p.(Arg601Trp) in an individual with Zellweger syndrome spectrum (ZSS) (PMID: 26700162). Two other heterozygotes with ZSS were reported; however, their second allele was unspecified (PMID: 19142205, 19877282). It has also been classified mutliple times as pathogenic/likely pathogenic and once as a VUS by diagnostic laboratories in ClinVar. (SP) 1007 - No published functional evidence has been identified for this variant. (I) 1201 - Heterozygous variant detected in trans with a pathogenic heterozygous variant (NM_000287.3(PEX6):c.1801C>T; p.(Arg601Trp)) in a recessive disease. (SP) 1206 - This variant has been shown to be paternally inherited (by trio analysis). (I) Legend: (SP) - Supporting pathogenic, (I) - Information, (SB) - Supporting benign -
Heimler syndrome 2 Pathogenic:1
Pathogenic, criteria provided, single submitterclinical testingBaylor GeneticsMar 12, 2024- -

Computational scores

Source: dbNSFP v4.3

Name
Calibrated prediction
Score
Prediction
AlphaMissense
Pathogenic
0.62
BayesDel_addAF
Pathogenic
0.42
D
BayesDel_noAF
Pathogenic
0.48
CADD
Pathogenic
35
DANN
Uncertain
1.0
DEOGEN2
Uncertain
0.54
D
Eigen
Pathogenic
0.82
Eigen_PC
Pathogenic
0.80
FATHMM_MKL
Pathogenic
0.99
D
LIST_S2
Pathogenic
0.99
D
M_CAP
Pathogenic
0.39
D
MetaRNN
Pathogenic
0.93
D
MetaSVM
Pathogenic
0.87
D
MutationAssessor
Benign
1.5
L
PrimateAI
Uncertain
0.65
T
PROVEAN
Uncertain
-2.7
D
REVEL
Pathogenic
0.92
Sift
Uncertain
0.0010
D
Sift4G
Uncertain
0.0060
D
Polyphen
1.0
D
Vest4
0.84
MutPred
0.69
Gain of disorder (P = 0.107);
MVP
0.94
MPC
0.64
ClinPred
0.86
D
GERP RS
5.5
Varity_R
0.72
gMVP
0.80

Splicing

Name
Calibrated prediction
Score
Prediction
dbscSNV1_ADA
Pathogenic
1.0
dbscSNV1_RF
Pathogenic
0.98
SpliceAI score (max)
0.76
Details are displayed if max score is > 0.2
DS_DL_spliceai
0.76
Position offset: 0

Find out detailed SpliceAI scores and Pangolin per-transcript scores at spliceailookup.broadinstitute.org

Publications

LitVar

Below is the list of publications found by LitVar. It may be empty.

Other links and lift over

dbSNP: rs267608240; hg19: chr6-42933782; API