chr2-190910320-A-G

Variant summary

Our verdict is Uncertain significance. The variant received 0 ACMG points: 1P and 1B. PP3BS1_Supporting

The NM_014905.5(GLS):​c.1037A>G​(p.Lys346Arg) variant causes a missense, splice region change involving the alteration of a conserved nucleotide. The variant allele was found at a frequency of 0.0000457 in 1,554,886 control chromosomes in the GnomAD database, with no homozygous occurrence. In-silico tool predicts a benign outcome for this variant. 3/3 splice prediction tools predicting alterations to normal splicing. Variant has been reported in ClinVar as Uncertain significance (no stars). Another variant affecting the same amino acid position, but resulting in a different missense (i.e. K346Q) has been classified as Uncertain significance.

Frequency

Genomes: 𝑓 0.000020 ( 0 hom., cov: 32)
Exomes 𝑓: 0.000048 ( 0 hom. )

Consequence

GLS
NM_014905.5 missense, splice_region

Scores

1
7
11
Splicing: ADA: 0.9994
2

Clinical Significance

Uncertain significance no assertion criteria provided U:1

Conservation

PhyloP100: 8.64

Publications

1 publications found
Variant links:
Genes affected
GLS (HGNC:4331): (glutaminase) This gene encodes the K-type mitochondrial glutaminase. The encoded protein is an phosphate-activated amidohydrolase that catalyzes the hydrolysis of glutamine to glutamate and ammonia. This protein is primarily expressed in the brain and kidney plays an essential role in generating energy for metabolism, synthesizing the brain neurotransmitter glutamate and maintaining acid-base balance in the kidney. Alternate splicing results in multiple transcript variants. [provided by RefSeq, Jan 2012]
STAT1 (HGNC:11362): (signal transducer and activator of transcription 1) The protein encoded by this gene is a member of the STAT protein family. In response to cytokines and growth factors, STAT family members are phosphorylated by the receptor associated kinases, and then form homo- or heterodimers that translocate to the cell nucleus where they act as transcription activators. The protein encoded by this gene can be activated by various ligands including interferon-alpha, interferon-gamma, EGF, PDGF and IL6. This protein mediates the expression of a variety of genes, which is thought to be important for cell viability in response to different cell stimuli and pathogens. The protein plays an important role in immune responses to viral, fungal and mycobacterial pathogens. Mutations in this gene are associated with Immunodeficiency 31B, 31A, and 31C. [provided by RefSeq, Jun 2020]
STAT1 Gene-Disease associations (from GenCC):
  • autoimmune enteropathy and endocrinopathy - susceptibility to chronic infections syndrome
    Inheritance: AD Classification: DEFINITIVE, STRONG, SUPPORTIVE Submitted by: Ambry Genetics, Labcorp Genetics (formerly Invitae), Illumina, Orphanet
  • immunodeficiency 31B
    Inheritance: AR Classification: DEFINITIVE, STRONG, SUPPORTIVE, LIMITED Submitted by: G2P, Ambry Genetics, Labcorp Genetics (formerly Invitae), Orphanet
  • Mendelian susceptibility to mycobacterial diseases due to partial STAT1 deficiency
    Inheritance: AD Classification: STRONG, MODERATE, SUPPORTIVE Submitted by: Ambry Genetics, Orphanet, Labcorp Genetics (formerly Invitae)

Genome browser will be placed here

ACMG classification

Classification was made for transcript

Our verdict: Uncertain_significance. The variant received 0 ACMG points.

PP3
Splicing predictors support a deleterious effect. Scorers claiming Pathogenic: dbscSNV1_ADA, dbscSNV1_RF, max_spliceai. No scorers claiming Uncertain. No scorers claiming Benign.
BS1
Variant frequency is greater than expected in population nfe. GnomAd4 allele frequency = 0.0000197 (3/152144) while in subpopulation NFE AF = 0.0000441 (3/67996). AF 95% confidence interval is 0.0000117. There are 0 homozygotes in GnomAd4. There are 3 alleles in the male GnomAd4 subpopulation. Median coverage is 32. This position passed quality control check. Existence of Clinvar submissions makes me limit the strength of this signal to Supporting

Transcripts

RefSeq

Gene Transcript HGVSc HGVSp Effect Exon rank MANE Protein UniProt
GLSNM_014905.5 linkc.1037A>G p.Lys346Arg missense_variant, splice_region_variant Exon 7 of 18 ENST00000320717.8 NP_055720.3 O94925-1

Ensembl

Gene Transcript HGVSc HGVSp Effect Exon rank TSL MANE Protein Appris UniProt
GLSENST00000320717.8 linkc.1037A>G p.Lys346Arg missense_variant, splice_region_variant Exon 7 of 18 1 NM_014905.5 ENSP00000317379.3 O94925-1

Frequencies

GnomAD3 genomes
AF:
0.0000197
AC:
3
AN:
152144
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.0000441
Gnomad OTH
AF:
0.00
GnomAD2 exomes
AF:
0.0000207
AC:
5
AN:
241608
AF XY:
0.0000382
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 FIN exome
AF:
0.00
Gnomad NFE exome
AF:
0.0000455
Gnomad OTH exome
AF:
0.00
GnomAD4 exome
AF:
0.0000485
AC:
68
AN:
1402742
Hom.:
0
Cov.:
21
AF XY:
0.0000514
AC XY:
36
AN XY:
700286
show subpopulations
African (AFR)
AF:
0.00
AC:
0
AN:
31988
American (AMR)
AF:
0.00
AC:
0
AN:
42906
Ashkenazi Jewish (ASJ)
AF:
0.00
AC:
0
AN:
25422
East Asian (EAS)
AF:
0.00
AC:
0
AN:
39024
South Asian (SAS)
AF:
0.00
AC:
0
AN:
82206
European-Finnish (FIN)
AF:
0.00
AC:
0
AN:
52926
Middle Eastern (MID)
AF:
0.00
AC:
0
AN:
5658
European-Non Finnish (NFE)
AF:
0.0000611
AC:
65
AN:
1064296
Other (OTH)
AF:
0.0000514
AC:
3
AN:
58316
Allele Balance Distribution
Red line indicates average allele balance
Average allele balance: 0.433
Heterozygous variant carriers
0
4
7
11
14
18
0.00
0.20
0.40
0.60
0.80
0.95
Allele balance

Age Distribution

Exome Het
Variant carriers
0
4
8
12
16
20
<30
30-35
35-40
40-45
45-50
50-55
55-60
60-65
65-70
70-75
75-80
>80
Age
GnomAD4 genome
AF:
0.0000197
AC:
3
AN:
152144
Hom.:
0
Cov.:
32
AF XY:
0.0000404
AC XY:
3
AN XY:
74330
show subpopulations
African (AFR)
AF:
0.00
AC:
0
AN:
41460
American (AMR)
AF:
0.00
AC:
0
AN:
15284
Ashkenazi Jewish (ASJ)
AF:
0.00
AC:
0
AN:
3470
East Asian (EAS)
AF:
0.00
AC:
0
AN:
5200
South Asian (SAS)
AF:
0.00
AC:
0
AN:
4828
European-Finnish (FIN)
AF:
0.00
AC:
0
AN:
10584
Middle Eastern (MID)
AF:
0.00
AC:
0
AN:
316
European-Non Finnish (NFE)
AF:
0.0000441
AC:
3
AN:
67996
Other (OTH)
AF:
0.00
AC:
0
AN:
2094
Allele Balance Distribution
Red line indicates average allele balance
Average allele balance: 0.442
Heterozygous variant carriers
0
0
1
1
2
2
0.00
0.20
0.40
0.60
0.80
0.95
Allele balance

Age Distribution

Genome Het
Variant carriers
0
2
4
6
8
10
<30
30-35
35-40
40-45
45-50
50-55
55-60
60-65
65-70
70-75
75-80
>80
Age
Alfa
AF:
0.00
Hom.:
0
Bravo
AF:
0.0000113
ExAC
AF:
0.0000165
AC:
2

ClinVar

Significance: Uncertain significance
Submissions summary: Uncertain:1
Revision: no assertion criteria provided
LINK: link

Submissions by phenotype

not provided Uncertain:1
-
Department of Pathology and Laboratory Medicine, Sinai Health System
Significance:Uncertain significance
Review Status:no assertion criteria provided
Collection Method:clinical testing

The GLS p.Lys346Arg variant was not identified in the literature nor was it identified in ClinVar or LOVD 3.0. The variant was identified in dbSNP (ID: rs770047438) and in control databases in 5 of 241608 chromosomes at a frequency of 0.00002069 (Genome Aggregation Database March 6, 2019, v2.1.1). The variant was observed in the European (non-Finnish) population in 5 of 109908 chromosomes (freq: 0.000045), but was not observed in the African, Latino, Ashkenazi Jewish, East Asian, European (Finnish), Other, or South Asian populations. The p.Lys346 residue is not conserved in mammals and four out of five computational analyses (PolyPhen-2, SIFT, AlignGVGD, BLOSUM, MutationTaster) do not suggest a high likelihood of impact to the protein; however, this information is not predictive enough to rule out pathogenicity. The p.Lys346Arg variant occurs in the second last base of the exon; this position has been shown to be part of the splicing consensus sequence and variants involving this position sometimes affect splicing. In addition, 2 of 4 in silico or computational prediction software programs (SpliceSiteFinder, MaxEntScan, NNSPLICE, GeneSplicer) predict a greater than 10% difference in splicing; however this is not predictive enough to assume pathogenicity. In summary, based on the above information the clinical significance of this variant cannot be determined with certainty at this time. This variant is classified as a variant of uncertain significance. -

Computational scores

Source: dbNSFP v4.3

Name
Calibrated prediction
Score
Prediction
AlphaMissense
Benign
0.10
BayesDel_addAF
Benign
-0.045
T
BayesDel_noAF
Benign
-0.17
CADD
Pathogenic
28
DANN
Uncertain
1.0
DEOGEN2
Benign
0.33
T;.
Eigen
Uncertain
0.36
Eigen_PC
Uncertain
0.46
FATHMM_MKL
Pathogenic
0.99
D
LIST_S2
Uncertain
0.95
D;D
M_CAP
Benign
0.016
T
MetaRNN
Benign
0.34
T;T
MetaSVM
Benign
-0.74
T
MutationAssessor
Uncertain
2.1
M;M
PhyloP100
8.6
PrimateAI
Uncertain
0.78
T
PROVEAN
Uncertain
-2.7
D;D
REVEL
Benign
0.21
Sift
Benign
0.091
T;T
Sift4G
Benign
0.096
T;T
Polyphen
0.31
B;B
Vest4
0.40
MutPred
0.49
Loss of ubiquitination at K346 (P = 0.0412);Loss of ubiquitination at K346 (P = 0.0412);
MVP
0.75
MPC
1.3
ClinPred
0.64
D
GERP RS
5.3
Varity_R
0.47
gMVP
0.74
Mutation Taster
=21/79
disease causing

Splicing

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

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

Publications

Other links and lift over

dbSNP: rs770047438; hg19: chr2-191775046; API