Benchmarking computational variant effect predictors by their ability to infer human traits.

BACKGROUND: Computational variant effect predictors offer a scalable and increasingly reliable means of interpreting human genetic variation, but concerns of circularity and bias have limited previous methods for evaluating and comparing predictors. Population-level cohorts of genotyped and phenotyped participants that have not been used in predictor training can facilitate an unbiased benchmarking of available methods. Using a curated set of human gene-trait associations with a reported rare-variant burden association, we evaluate the correlations of 24 computational variant effect predictors with associated human traits in the UK Biobank and All of Us cohorts. RESULTS: AlphaMissense outperformed all other predictors in inferring human traits based on rare missense variants in UK Biobank and All of Us participants. The overall rankings of computational variant effect predictors in these two cohorts showed a significant positive correlation. CONCLUSION: We describe a method to assess computational variant effect predictors that sidesteps the limitations of previous evaluations. This approach is generalizable to future predictors and could continue to inform predictor choice for personal and clinical genetics.

A critical region of A20 unveiled by missense TNFAIP3 variations that lead to autoinflammation.

A20 haploinsufficiency (HA20) is an autoinflammatory disease caused by heterozygous loss-of-function variations in TNFAIP3, the gene encoding the A20 protein. Diagnosis of HA20 is challenging due to its heterogeneous clinical presentation and the lack of pathognomonic symptoms. While the pathogenic effect of TNFAIP3 truncating variations is clearly established, that of missense variations is difficult to determine. Herein, we identified a novel TNFAIP3 variation, p.(Leu236Pro), located in the A20 ovarian tumor (OTU) domain and demonstrated its pathogenicity. In the patients' primary cells, we observed reduced A20 levels. Protein destabilization was predicted in silico for A20_Leu236Pro and enhanced proteasomal degradation was confirmed in vitro through a flow cytometry-based functional assay. By applying this approach to the study of another missense variant, A20_Leu275Pro, for which no functional characterization has been performed to date, we showed that this variant also undergoes enhanced proteasomal degradation. Moreover, we showed a disrupted ability of A20_Leu236Pro to inhibit the NF-κB pathway and to deubiquitinate its substrate TRAF6. Structural modeling revealed that two residues involved in OTU pathogenic missense variations (i.e. Glu192Lys and Cys243Tyr) establish common interactions with Leu236. Interpretation of newly identified missense variations is challenging, requiring, as illustrated here, functional demonstration of their pathogenicity. Together with functional studies, in silico structure analysis is a valuable approach that allowed us (i) to provide a mechanistic explanation for the haploinsufficiency resulting from missense variations and (ii) to unveil a region within the OTU domain critical for A20 function.

Novel homozygous ADAMTS17 missense variant in Weill-Marchesani syndrome.

AIM: To explore the phenotype and genotype of Weill-Marchesani syndrome (WMS) in a Chinese family and review related literature. METHODS: Three WMS patients and other unaffected individuals in this family with a history of consanguineous marriage were included in this study. Medical history, comprehensive ophthalmic examinations, and systemic evaluation, as well as whole exome and Sanger sequencing of specific genomic regions, were performed. RESULTS: The three affected siblings presented with short stature, brachydactyly and ocular disorders, including very shallow anterior chamber, high myopia, microspherophakia lens subluxation with stretched zonules and glaucoma. Genetic analysis verified a homozygous missense mutation (c.2983C>T: p. Arg995Trp) in ADAMTS17, which was correlated with the diseases in this family, indicating an autosomal recessive inherited manner of WMS. This review aims to summarize the mutation sites of WMS genes, so as to prevent the disease and better guide clinical diagnosis and treatment. CONCLUSION: A novel homozygous missense variant of ADAMTS17 is identified in a WMS family with a history of consanguineous marriage. Our study expands the range of mutations associated with WMS and deepens our understanding of pathology in disease associated with ADAMTS17 variants.

Novel KCNC2 variant associated with developmental and epileptic encephalopathy.

The KCNC2 gene encodes Kv3.2, which is a member of the voltage-gated potassium channel subfamily. It is crucial for the generation of fast-spiking properties in cortical GABAergic interneurons. Recently, KCNC2 variations were found to be associated with epileptic encephalopathy in unrelated individuals. Here, we report a Chinese patient with developmental and epileptic encephalopathy (DEE) and motor development delay. Whole-exome sequencing (WES) revealed a novel heterozygous variant in the KCNC2 gene NM_139137.4:c.1163T>C (p.Phe388Ser), and subsequent Sanger sequencing showed that it was a de novo mutation. We identified the KCNC2 likely pathogenic variant in a DEE patient by reanalysis of WES data in a Chinese family. Our study enriched the variation spectrum of the KCNC2 gene and promoted the application of WES technology and data reanalysis in the diagnosis of epilepsy.

Han family with essential tremor caused by the P421L variant of the TENM4 gene in China.

BACKGROUND: Essential tremor (ET) is an autosomal dominant inheritance disorder. Mutations in fusion sarcoma (FUS), mitochondrial serine peptidase 2 (HTRA2), teneurin transmembrane protein 4 (TENM4), sortilin1 (SORT1), SCN11A, and notch2N-terminal-like (NOTCH2NLC) genes are associated with familial ET. METHODS: A proband with ET was tested using whole-exome sequencing and repeat-primed polymerase chain reaction. Subsequently, the family members were screened for the suspected mutation, and the results were verified using Sanger sequencing. The relationship between pedigree and phenotype was also analyzed, and structural and functional changes in the variants were predicted using bioinformatics analysis. RESULTS: In a family with ET, the proband (III4) and the proband's father (II1), grandfather (I1), uncle (II2), and cousin (III5) all presented with involuntary tremors of both upper limbs. The responsible mutation was identified as TENM4 c.1262C > T (p.P421L), which showed genetic co-segregation in the family survey. AlphaFold predicted a change in the spatial position of TENM4 after the P421L mutation, which may have affected its stability. AlphaFold also predicted P421L to be a deleterious variation, which would lead to lower degrees of freedom of the TENM4 protein, thereby affecting the protein's structure and stability. According to the bioinformatics analysis, TENM4 (p.P421L) may reduce the signal reaching the nucleus by affecting the expression of TENM4 messenger RNA (mRNA), thereby impairing the normal oligodendrocyte differentiation process and leading to impaired myelination. CONCLUSION: This study revealed that the TENM4 (p.P421L) pathogenic missense variation was responsible for ET in the proband.

Case Report of Neonatal Sotos Syndrome with a New Missense Mutation in the NSD1 Gene and Literature Analysis in the Chinese Han Population.

Currently, no consensus exists regarding Sotos syndrome in the Chinese population. Here, we present a case of neonatal Sotos syndrome, followed by a retrospective analysis of five cases of neonatal Sotos syndrome, reported in China. The study subject was a twin premature infant, heavier than gestational age, with characteristic facial features, limb shaking, and hypertonia. Transient hypoglycemia, abnormal cranial magnetic resonance imaging, multiple nodules in polycystic kidneys and liver, abnormal hearing, patent ductus arteriosus, and an atrial septal defect were also noted. The subject showed overgrowth and developmental retardation at 3 months of age. Sequencing revealed a novel missense mutation, c.5000C>A, in the nuclear receptor binding the SET domain protein 1 gene, resulting in an alanine-to-glutamate substitution. The bioinformatics analysis suggested high pathogenicity at this site. This study provides insights into diagnosis of neonatal Sotos syndrome based on specific phenotypes. Subsequent treatment and follow-up should focus on developmental retardation, epilepsy, and scoliosis.

Mutational and bioinformatics analysis of the NKX2.1 gene in a cohort of Iranian pediatric patients with congenital hypothyroidism (CH).

Congenital hypothyroidism (CH) occurs with a relatively alarming prevalence in infants, and if not diagnosed and treated in time, it can have devastating consequences for the development of the nervous system. CH is associated with genetic changes in several genes that encode transcription factors responsible for thyroid development, including mutations in the NK2 homeobox 1 (NKX2.1) gene, which encodes the thyroid transcription factor-1 (TTF-1). Although CH is frequently observed in pediatric populations, there is still a limited understanding of the genetic factors and molecular mechanisms contributing to this disease. The sequence of the NKX2.1 gene was investigated in 75 pediatric patients with CH by polymerase chain reaction (PCR), single-stranded conformation polymorphism (SSCP), and direct DNA sequencing. Four missense heterozygous variations were identified in exon 3 of the NKX2.1 gene, including three novel missense variations, namely c.708A>G, p.Gln202Arg; c.713T>G, p.Tyr204Asp; c.833T>G, p.Tyr244Asp, and a previously reported variant rs781133468 (c.772C>G, p.His223Gln). Importantly, these variations occur in highly conserved residues of the TTF-1 DNA-binding domain and were predicted by bioinformatics analysis to alter the protein structure, with a probable alteration in the protein function. These results indicate that nucleotide changes in the NKX2.1 gene may contribute to CH pathogenesis.

Detection and Gene Mutation Analysis of Three Variations in Two Unrelated Chinese Hereditary Coagulation Factor XI Deficiency Families.

INTRODUCTION: Three variations including a novel F11 gene variation were detected in two unrelated Chinese families with coagulation factor XI deficiency, and their possible pathogenesis was elucidated. METHODS: The genomic DNA of the probands' pedigrees was extracted, and all exons and flanking sequences of F11 gene were subjected to PCR amplification and Sanger sequencing. ClustalX-2.1-win, Mutation Taster, and Swiss-Pdb Viewer software were used to analyze the conservation and impact of the variations on protein function and structure. RESULTS: DNA sequencing showed that the proband one had p.Gly350Glu and p.Trp501stop complex heterozygous variations, while the proband two took p.Pro338Leu and p.Trp501stop compound heterozygous variations. Conservation, structural, and functional analysis of variant amino acids indicated that these three variations were harmful and probably affected the structure and function of the variable protein. CONCLUSIONS: Three variations including p.Pro338Leu, p.Gly350Glu, and p.Trp501stop responsible for the reduction of the FXI activities were herein detected. Notably, the p.Pro338Leu variation was discovered for the first time in the world. Furthermore, the p.Gly350Glu was first reported in China.

Understanding the impact of ZBTB18 missense variation on transcription factor function in neurodevelopment and disease.

Mutations to genes that encode DNA-binding transcription factors (TFs) underlie a broad spectrum of human neurodevelopmental disorders. Here, we highlight the pathological mechanisms arising from mutations to TF genes that influence the development of mammalian cerebral cortex neurons. Drawing on recent findings for TF genes including ZBTB18, we discuss how functional missense mutations to such genes confer non-native gene regulatory actions in developing neurons, leading to cell-morphological defects, neuroanatomical abnormalities during foetal brain development and functional impairment. Further, we discuss how missense variation to human TF genes documented in the general population endow quantifiable changes to transcriptional regulation, with potential cell biological effects on the temporal progression of cerebral cortex neuron development and homeostasis. We offer a systematic approach to investigate the functional impact of missense variation in brain TFs and define their direct molecular and cellular actions in foetal neurodevelopment, tissue homeostasis and disease states.

Novel TBK1 variant associated with Frontotemporal Dementia overlap syndrome.

BACKGROUND: Recently, TANK binding kinase 1 (TBK1) mutation has been reported as a causative gene for overlap frontotemporal dementia (FTD)-amyotrophic lateral sclerosis (ALS) syndrome. However, there are no reports from families of South Asian ethnicity. OBJECTIVE: To report a case study of a family with the proband having overlap FTD-ALS syndrome caused by a novel TBK1 variant. MATERIALS AND METHODS: Clinical, brain imaging, genetic analysis and laboratory data of the patient with FTD-ALS were performed. In addition, family-based segregation analysis of identified novel variants was also done. RESULTS: This study pertains to genetic analysis in 11 members in a family with only one member affected with overlap FTD-ALS syndrome. The whole-exome sequencing analysis in the symptomatic member showed a novel loss-of-function (LoF) variant c.1810G>T(p.E604X) in the TBK1 gene. Neuroimaging showed a pattern of asymmetric frontotemporal atrophy and hypometabolism. Segregation analysis of the variation demonstrated its presence in several family members, although none of the other members was symptomatic. Further, we observed another missense variation in the NEFH gene (p.Pro683Leu) which was seen in the symptomatic and two asymptomatic family members, the pathogenicity of which is unclear. CONCLUSION: This is the first study of a rare novel TBK1 variant associated with FTD-ALS from India. Asymptomatic family members with the variant have important clinical implications and necessitate the genetic evaluation and long-term follow-up of family members of patients detected with TBK1 mutations. Therefore, although infrequent, genetic screening for the TBK1 gene should be considered when encountering overlap FTD syndromes.

Case Report: Corneal Leucoma as a Novel Clinical Presentation of Nail-Patella Syndrome in a 5-Year-Old Girl.

Nail-patella syndrome (NPS) is a rare autosomal-dominant disorder characterized by the classic tetrad of absent or hypoplastic finger and toe nails, absent or hypoplastic patella, skeletal deformities involving the elbow joints, and iliac horns. This disease is caused by heterozygous pathogenic variations in the LMX1B gene, which encodes the LIM homeodomain transcription factor protein (LMX1B). We report a case of corneal leucoma and dysplasia prior to overt steroid-resistant nephrotic syndrome (SRNS) in a patient with NPS. At presentation, the parents of a 5-year-old female patient reported their daughter had corneal leucoma, psychomotor delay and speech defect. We also noted the presence of bilateral edema of the lower extremities, hypertension, nail dystrophy, and the bilateral absence of patella. She developed steroid-resistant nephrotic syndrome. Lowe oculocerebrorenal syndrome and NPS were the conditions considered in differential diagnosis. Trio-based whole genome sequencing indicated a heterozygous de novo likely pathogenic variation in the LMX1B gene (c.805A>C [p.Asn269His]). Patients with NPS often develop nail, ocular, or orthopedic symptoms prior to nephrotic syndrome. Corneal leucoma may be a novel clinical presentation of NPS.

Analysis of Missense Variants in the Human Histamine Receptor Family Reveals Increased Constitutive Activity of E4106.30×30K Variant in the Histamine H1 Receptor.

The Exome Aggregation Consortium has collected the protein-encoding DNA sequences of almost 61,000 unrelated humans. Analysis of this dataset for G protein-coupled receptor (GPCR) proteins (available at GPCRdb) revealed a total of 463 naturally occurring genetic missense variations in the histamine receptor family. In this research, we have analyzed the distribution of these missense variations in the four histamine receptor subtypes concerning structural segments and sites important for GPCR function. Four missense variants R1273.52×52H, R13934.57×57H, R4096.29×29H, and E4106.30×30K, were selected for the histamine H1 receptor (H1R) that were hypothesized to affect receptor activity by interfering with the interaction pattern of the highly conserved D(E)RY motif, the so-called ionic lock. The E4106.30×30K missense variant displays higher constitutive activity in G protein signaling as compared to wild-type H1R, whereas the opposite was observed for R1273.52×52H, R13934.57×57H, and R4096.29×29H. The E4106.30×30K missense variant displays a higher affinity for the endogenous agonist histamine than wild-type H1R, whereas antagonist affinity was not affected. These data support the hypothesis that the E4106.30×30K mutation shifts the equilibrium towards active conformations. The study of these selected missense variants gives additional insight into the structural basis of H1R activation and, moreover, highlights that missense variants can result in pharmacologically different behavior as compared to wild-type receptors and should consequently be considered in the drug discovery process.

Molecular signatures of local adaptation to light in Norway spruce.

MAIN CONCLUSION: Transcriptomic and exome capture analysis reveal an adaptive cline for shade tolerance in Norway spruce. Genes involved in the lignin pathway and immunity seem to play a potential role in contributing towards local adaptation to light. The study of natural variation is an efficient method to elucidate how plants adapt to local climatic conditions, a key process for the evolution of a species. Norway spruce is a shade-tolerant conifer in which the requirement of far-red light for growth increases latitudinally northwards. The objective of the study is to characterize the genetic control of local adaptation to light enriched in far-red in Norway spruce, motivated by a latitudinal gradient for the Red:Far-red (R:FR) ratio to which Norway spruce has been proven to be genetically adapted. We have established the genomic signatures of local adaptation by conducting transcriptomic (total RNA-sequencing) and genomic analyses (exome capture), for the identification of genes differentially regulated along the cline. RNA-sequencing revealed 274 differentially expressed genes in response to SHADE (low R:FR light), between the southern and northern natural populations in Sweden. Exome capture included analysis of a uniquely large data set (1654 trees) that revealed missense variations in coding regions of nine differentially expressed candidate genes, which followed a latitudinal cline in allele and genotype frequencies. These genes included five transcription factors involved in vital processes like bud-set/bud-flush, lignin pathway, and cold acclimation and other genes that take part in cell-wall remodeling, secondary cell-wall thickening, response to starvation, and immunity. Based on these results, we suggest that the northern populations might not only be able to adjust their growing season in response to low R:FR light, but they may also be better adapted towards disease resistance by up-regulation of the lignin pathway that is linked to immunity. This forms a concrete basis for local adaptation to light quality in Norway spruce, one of the most economically important conifer tree species in Sweden.

De novo missense variants disrupting protein-protein interactions affect risk for autism through gene co-expression and protein networks in neuronal cell types.

BACKGROUND: Whole-exome sequencing studies have been useful for identifying genes that, when mutated, affect risk for autism spectrum disorder (ASD). Nonetheless, the association signal primarily arises from de novo protein-truncating variants, as opposed to the more common missense variants. Despite their commonness in humans, determining which missense variants affect phenotypes and how remains a challenge. We investigate the functional relevance of de novo missense variants, specifically whether they are likely to disrupt protein interactions, and nominate novel genes in risk for ASD through integrated genomic, transcriptomic, and proteomic analyses. METHODS: Utilizing our previous interactome perturbation predictor, we identify a set of missense variants that are likely disruptive to protein-protein interactions. For genes encoding the disrupted interactions, we evaluate their expression patterns across developing brains and within specific cell types, using both bulk and inferred cell-type-specific brain transcriptomes. Connecting all disrupted pairs of proteins, we construct an "ASD disrupted network." Finally, we integrate protein interactions and cell-type-specific co-expression networks together with published association data to implicate novel genes in ASD risk in a cell-type-specific manner. RESULTS: Extending earlier work, we show that de novo missense variants that disrupt protein interactions are enriched in individuals with ASD, often affecting hub proteins and disrupting hub interactions. Genes encoding disrupted complementary interactors tend to be risk genes, and an interaction network built from these proteins is enriched for ASD proteins. Consistent with other studies, genes identified by disrupted protein interactions are expressed early in development and in excitatory and inhibitory neuronal lineages. Using inferred gene co-expression for three neuronal cell types-excitatory, inhibitory, and neural progenitor-we implicate several hundred genes in risk (FDR [Formula: see text]0.05), ~ 60% novel, with characteristics of genuine ASD genes. Across cell types, these genes affect neuronal morphogenesis and neuronal communication, while neural progenitor cells show strong enrichment for development of the limbic system. LIMITATIONS: Some analyses use the imperfect guilt-by-association principle; results are statistical, not functional. CONCLUSIONS: Disrupted protein interactions identify gene sets involved in risk for ASD. Their gene expression during brain development and within cell types highlights how they relate to ASD.

In Silico Study of Correlation between Missense Variations of F8 Gene and Inhibitor Formation in Severe Hemophilia A

Objective: Deleterious substitutions of the F8 gene are responsible for causing hemophilia A, which is an inherited bleeding disorder resulting from reduced or absent activity of the coagulant protein factor VIII (FVIII). The most important complication in treatment is inhibitor development toward therapeutic factor VIII. In this study, we aimed to analyze the effects of deleterious substitutions in the F8 gene upon protein structure and function. Materials and Methods: All tests were conducted by computational methods from the CHAMP (CDC Hemophilia A Mutation Project) database. We performed an in silico analysis of deleterious variations using five software programs, Sift, PolyPhen-2, Align-GVGD, KD4v, and MutationTaster, in order to analyze the correlation between variation and the disease. We also studied the correlation between these variations and inhibitor formation. Results: Our analysis showed that these in silico tools are coherent and that there are more variations in the A than the C domains. Moreover, we noticed that there are more deleterious variations than neutral variations in each of the A and C domains. We also found that 13.51% of the patients suffered from a severe form of hemophilia A and that carriers of missense variations developed inhibitors. Also, for the first time, we determined that variation nature is not associated with inhibitor formation. Furthermore, this analysis showed that the risk of developing inhibitors increases when the variation causes a change of amino acid class. Conclusion: This study will help to correctly associate variations with inhibitor development and aid in early characterization of novel variants.

Genetic lesions in the UGT1A1 genes among Gilbert's syndrome patients from India.

The present study was undertaken to investigate genetic variations present in the coding regions of the UGT1A1 gene among the Gilbert's syndrome patients. Analysis of genetic variations was performed by direct DNA sequencing among the patients that do not have any polymorphic variations in the promoter regions of the UGT1A1 gene. We identified seven different sequence variations among Gilbert's Syndrome patients, of which four were novel. Out of seven variants, six missense and one silent single nucleotide substitutions were present in the UGT1A1 gene. In addition, molecular modeling of UGT1A1 (H55R, P152S and N212H) variants suggested a reduced activity of the enzyme. This study demonstrates that different variations present in the UGT1A1 gene and specifically, the H55R variation had a significant effect on bilirubin levels and could be genetic risk factors for hyperbilirubinemia.

GP6 Haplotype of Missense Variants is Associated with Sticky Platelet Syndrome Manifested by Fetal Loss.

Disequilibrium of hemostasis is central to the pathogenesis of all thromboses, and platelets are essential for primary hemostasis. The platelet membrane glycoprotein receptor is involved in the clot formation in blood; therefore, the changes in related genes could impair platelet aggregation in patients with sticky platelet syndrome (SPS). Patients with SPS who experienced fetal loss were shown to harbor a risk haplotype at GP6 locus. The aim of the study was to examine the genetic linkage of this selected risk haplotype with single nucleotide variations (SNVs) in the coding sequence of the GP6 gene in order to identify possible functional SNVs in association with SPS and fetal loss. A total of 37 patients with SPS manifested fetal loss, and 42 healthy controls were enrolled in the study. The SPS was diagnosed with platelet aggregometry. The SNVs were determined by dideoxy sequencing and high-resolution melting analysis. The missense variations were detected in patients with risk haplotype only. The association analysis showed association of the minor alleles with the SPS manifested by fetal loss as follows-rs1671152 (odds ratio [OR]: 4.667, 95% confidence interval [CI]: 1.462-14.89, P = .006), rs2304167 (OR: 5.085, 95% CI: 1.605-16.10, P = .003), and rs1654416 (OR: 5.085, 95% CI: 1.605-16.10, P = .003). Using the Expectation-Maximization (EM) algorithm, the estimated minor haplotype with predicted protein residue PEAN was significantly associated with the given phenotype (OR: 4.746, 95% CI: 1.486-15.15, P = .005). We have shown that haplotype PEAN associated with SPS and manifested by fetal loss and suggest that the mechanism involved in the action of GPVI has significant effect on GPVI-mediated signal transduction through Syk-phosphorylation.

An Alpha-kinase 2 Gene Variant Disrupts Filamentous Actin Localization in the Surface Cells of Colorectal Cancer Spheroids.

BACKGROUND/AIM: Alpha-kinase 2 (ALPK2), suggested to be a novel tumour-suppressor gene down-regulated by oncogenic KRAS, plays a pivotal role in luminal apoptosis in normal colonic crypts. The aim of this study was to determine the association between ALPK2 germline variants and colorectal cancer. MATERIALS AND METHODS: Missense single nucleotide variants in the exons of the ALPK2 gene in 2,343 consecutive autopsy cases (1,446 cases with cancer and 897 cases without cancer) were screened using HumanExome BeadChip arrays. To address the functional effect of a missense ALPK2 variant, a 3D floating cell culture was performed using HCT116-derived human colorectal cancer cells stably expressing wild-type (wt) ALPK2 (HCT116-wtALPK2) or amino acid-substituted (sub) ALPK2 (HCT116-subALPK2). RESULTS: We identified that one of the ALPK2 germline variants, rs55674018 (p.Q1853E), was significantly associated with the presence of cancer (adjusted odds ratio(OR)=4.39; 95% confidence interval(CI)=1.31-14.78, p=0.001). The p.Q1853E variant was present in the East Asian population and located in the immunoglobulin-like domain. Notably, the basolateral polarity of actin in the surface of HCT116-wtALPK2 spheroids was more attenuated compared to that of HCT116-subALPK2 spheroids. Furthermore, luminal apoptosis and cell aggregation were promoted by wtALPK2, but not by subALPK2 in 3D culture. CONCLUSION: The p.Q1853E variant of ALPK2, which had been accumulating in the Japanese population, induced a metastatic phenotype by disrupting ALPK2 function.