Germline variants in tumor suppressor FBXW7 lead to impaired ubiquitination and a neurodevelopmental syndrome

Copy number variants at 4q31.3 affecting the regulatory region of FBXW7 associated with neurodevelopmental delay

Emerging research has demonstrated that genomic alterations disrupting topologically associated domains (TADs) and chromatin interactions underlie the pathogenic mechanisms of specific copy number variants (CNVs) in neurodevelopmental disorders. We report two patients with a de novo deletion and a duplication in chromosome 4q31, potentially causing FBX-related neurodevelopmental syndrome by affecting the regulatory region of FBXW7. High-throughput chromosome conformation capture (Hi-C) analysis using available capture data in neural progenitor cells revealed the rewiring of the TAD boundary close to FBXW7. Both patients exhibited facial dysmorphisms, cardiac and limb abnormalities, and neurodevelopmental delays, showing significant clinical overlap with previously reported FBXW7-related features. We also included an additional 10 patients with CNVs in the 4q31 region from the literature and the DECIPHER database for Hi-C analysis, which confirmed that disruption of the regulatory region of FBXW7 likely contributes to the developmental defects observed in these patients.

Reciprocal chromosome translocation t(3;4)(q27;q31.2) with deletion of 3q27 and reduced FBXW7 expression in a patient with developmental delay, hypotonia, and seizures

Reciprocal chromosomal translocation is one of genomic variations. When cytogenetically de novo reciprocal translocations are identified in patients with some clinical manifestations, the genes in the breakpoints are considered to be related to the clinical features. In this study, we encountered a patient with severe developmental delay, intractable epilepsy, growth failure, distinctive features, and skeletal manifestations. Conventional karyotyping revealed a de novo translocation described as 46,XY,t(3;4)(q27;q31.2). Chromosomal microarray testing detected a 1.25-Mb microdeletion at 3q27.3q28. Although the skeletal manifestations may have been affected by this deletion, the neurological features of this patient were severe and could not be fully explained by this deletion. Since no genomic copy number aberration was detected on chromosome 4, long-read whole-genome sequencing analysis was performed and a precise breakpoint was confirmed. A 460-bp deletion was detected between the two breakpoints; however, no gene was disrupted. FBXW7, the gene responsible for developmental delay, hypotonia, and impaired language, is in the 0.5-Mb telomeric region. Most of the patient's clinical features were considered consistent with symptoms of FBXW7-related disorders, but were more severe. FBXW7 expression in the immortalized lymphoblasts of the patient was reduced compared to that in controls. Based on these findings, we suspect that FBXW7 is affected by downstream position effects of chromosomal translocations. The severe neurological features of the patient may have been affected not only by the 3q27-q28 deletion but also by impaired expression of FBXW7 derived from the breakage of chromosome 4.

Further evidence that the neurodevelopmental gene FBXW7 predisposes to Wilms tumor

Somatic variants in the NOTCH pathway regulator FBXW7 are frequently seen in a variety of malignancies. Heterozygous loss-of-function germline variants in FBXW7 have recently been described as causative for a neurodevelopmental syndrome. Independently, FBXW7 was also considered as a susceptibility gene for Wilms tumor due to a few observations of heterozygous germline variants in patients with Wilms tumor. Whether the same FBXW7 variants are implicated in both, neurodevelopmental delay and Wilms tumor formation, remained unclear. By clinical testing, we now observed a patient with neurodevelopmental delay due to a de novo constitutional mosaic FBXW7 splice site pathogenic variant who developed Wilms tumor. In the tumor, we identified a second hit frameshift variant in FBXW7. Immunohistochemical staining was consistent with mosaic loss of FBXW7 protein expression in the tumor. Our data support the role of constitutional FBXW7 pathogenic variants in both, neurodevelopmental disorder and the etiology of Wilms tumor. Therefore, Wilms tumor screening should be considered in individuals with constitutional or germline pathogenic variants in FBXW7 and associated neurodevelopmental syndrome.

MYCN in human development and diseases

Somatic mutations in MYCN have been identified across various tumors, playing pivotal roles in tumorigenesis, tumor progression, and unfavorable prognoses. Despite its established notoriety as an oncogenic driver, there is a growing interest in exploring the involvement of MYCN in human development. While MYCN variants have traditionally been associated with Feingold syndrome type 1, recent discoveries highlight gain-of-function variants, specifically p.(Thr58Met) and p.(Pro60Leu), as the cause for megalencephaly-polydactyly syndrome. The elucidation of cellular and murine analytical data from both loss-of-function (Feingold syndrome model) and gain-of-function models (megalencephaly-polydactyly syndrome model) is significantly contributing to a comprehensive understanding of the physiological role of MYCN in human development and pathogenesis. This review discusses the MYCN's functional implications for human development by reviewing the clinical characteristics of these distinct syndromes, Feingold syndrome, and megalencephaly-polydactyly syndrome, providing valuable insights into the understanding of pathophysiological backgrounds of other syndromes associated with the MYCN pathway and the overall comprehension of MYCN's role in human development.

Craniosynostosis Associated With Novel TUBG1 Mutation (NM_001070.4:c.821C>T) (p.Thr274Ile)

TUBG1, a tubulin gene, plays an important role in neurodevelopment. Here we describe a case of a novel TUGB1 mutation (NM_001070.4:c.821C>T) (p.Thr274Ile). This patient presented similarly to previous cases with features including microcephaly, epilepsy, and speech and motor delay. Unique characteristics were also present such as trigonocephaly, tethered frenulum, scoliosis, nystagmus, and a concurrent FBXW7 mutation. This case expands our breadth of knowledge on TUBG1 genotypic and phenotypic variation. However, further work is needed to fully understand this rare mutation and the associations between TUBG1 and FBXW7 mutations.

Diagnostic yield and clinical impact of germline sequencing in children with CNS and extracranial solid tumors-a nationwide, prospective Swedish study

Background

Childhood cancer predisposition (ChiCaP) syndromes are increasingly recognized as contributing factors to childhood cancer development. Yet, due to variable availability of germline testing, many children with ChiCaP might go undetected today. We report results from the nationwide and prospective ChiCaP study that investigated diagnostic yield and clinical impact of integrating germline whole-genome sequencing (gWGS) with tumor sequencing and systematic phenotyping in children with solid tumors.

Methods

gWGS was performed in 309 children at diagnosis of CNS (n = 123, 40%) or extracranial (n = 186, 60%) solid tumors and analyzed for disease-causing variants in 189 known cancer predisposing genes. Tumor sequencing data were available for 74% (227/309) of patients. In addition, a standardized clinical assessment for underlying predisposition was performed in 95% (293/309) of patients.

Findings

The prevalence of ChiCaP diagnoses was 11% (35/309), of which 69% (24/35) were unknown at inclusion (diagnostic yield 8%, 24/298). A second-hit and/or relevant mutational signature was observed in 19/21 (90%) tumors with informative data. ChiCaP diagnoses were more prevalent among patients with retinoblastomas (50%, 6/12) and high-grade astrocytomas (37%, 6/16), and in those with non-cancer related features (23%, 20/88), and ≥2 positive ChiCaP criteria (28%, 22/79). ChiCaP diagnoses were autosomal dominant in 80% (28/35) of patients, yet confirmed de novo in 64% (18/28). The 35 ChiCaP findings resulted in tailored surveillance (86%, 30/35) and treatment recommendations (31%, 11/35).

Interpretation

Overall, our results demonstrate that systematic phenotyping, combined with genomics-based diagnostics of ChiCaP in children with solid tumors is feasible in large-scale clinical practice and critically guides personalized care in a sizable proportion of patients.

Funding

The study was supported by the Swedish Childhood Cancer Fund and the Ministry of Health and Social Affairs.

Characterizing and predicting ccRCC-causing missense mutations in Von Hippel-Lindau disease

BACKGROUND

Mutations within the Von Hippel-Lindau (VHL) tumor suppressor gene are known to cause VHL disease, which is characterized by the formation of cysts and tumors in multiple organs of the body, particularly clear cell renal cell carcinoma (ccRCC). A major challenge in clinical practice is determining tumor risk from a given mutation in the VHL gene. Previous efforts have been hindered by limited available clinical data and technological constraints.

METHODS

To overcome this, we initially manually curated the largest set of clinically validated VHL mutations to date, enabling a robust assessment of existing predictive tools on an independent test set. Additionally, we comprehensively characterized the effects of mutations within VHL using in silico biophysical tools describing changes in protein stability, dynamics and affinity to binding partners to provide insights into the structure-phenotype relationship. These descriptive properties were used as molecular features for the construction of a machine learning model, designed to predict the risk of ccRCC development as a result of a VHL missense mutation.

RESULTS

Analysis of our model showed an accuracy of 0.81 in the identification of ccRCC-causing missense mutations, and a Matthew's Correlation Coefficient of 0.44 on a non-redundant blind test, a significant improvement in comparison to the previous available approaches.

CONCLUSION

This work highlights the power of using protein 3D structure to fully explore the range of molecular and functional consequences of genomic variants. We believe this optimized model will better enable its clinical implementation and assist guiding patient risk stratification and management.

Exploring the effects of missense mutations on protein thermodynamics through structure-based approaches: findings from the CAGI6 challenges

Missense mutations are known contributors to diverse genetic disorders, due to their subtle, single amino acid changes imparted on the resultant protein. Because of this, understanding the impact of these mutations on protein stability and function is crucial for unravelling disease mechanisms and developing targeted therapies. The Critical Assessment of Genome Interpretation (CAGI) provides a valuable platform for benchmarking state-of-the-art computational methods in predicting the impact of disease-related mutations on protein thermodynamics. Here we report the performance of our comprehensive platform of structure-based computational approaches to evaluate mutations impacting protein structure and function on 3 challenges from CAGI6: Calmodulin, MAPK1 and MAPK3. Our stability predictors have achieved correlations of up to 0.74 and AUCs of 1 when predicting changes in ΔΔG for MAPK1 and MAPK3, respectively, and AUC of up to 0.75 in the Calmodulin challenge. Overall, our study highlights the importance of structure-based approaches in understanding the effects of missense mutations on protein thermodynamics. The results obtained from the CAGI6 challenges contribute to the ongoing efforts to enhance our understanding of disease mechanisms and facilitate the development of personalised medicine approaches.

AI-Driven Enhancements in Drug Screening and Optimization

The greatest challenge in drug discovery remains the high rate of attrition across the different phases of the process, which cost the industry billions of dollars every year. While all phases remain crucial to ensure pharmaceutical-level safety, quality, and efficacy of the end product, streamlining these efforts toward compounds with success potential is pivotal for a more efficient and cost-effective process. The use of artificial intelligence (AI) within the pharmaceutical industry aims at just this, and has applications in preclinical screening for biological activity, optimization of pharmacokinetic properties for improved drug formulation, early toxicity prediction which reduces attrition, and pre-emptively screening for genetic changes in the biological target to improve therapeutic longevity. Here, we present a series of in silico tools that address these applications in small molecule development and describe how they can be embedded within the current pharmaceutical development pipeline.

Simultaneous CRISPR screening and spatial transcriptomics reveals intracellular, intercellular, and functional transcriptional circuits

Pooled optical screens have enabled the study of cellular interactions, morphology, or dynamics at massive scale, but have not yet leveraged the power of highly-plexed single-cell resolved transcriptomic readouts to inform molecular pathways. Here, we present Perturb-FISH, which bridges these approaches by combining imaging spatial transcriptomics with parallel optical detection of in situ amplified guide RNAs. We show that Perturb-FISH recovers intracellular effects that are consistent with Perturb-seq results in a screen of lipopolysaccharide response in cultured monocytes, and uncover new intercellular and density-dependent regulation of the innate immune response. We further pair Perturb-FISH with a functional readout in a screen of autism spectrum disorder risk genes, showing common calcium activity phenotypes in induced pluripotent stem cell derived astrocytes and their associated genetic interactions and dysregulated molecular pathways. Perturb-FISH is thus a generally applicable method for studying the genetic and molecular associations of spatial and functional biology at single-cell resolution.

E3-ubiquitin ligase, FBXW7 regulates mitotic progression by targeting BubR1 for ubiquitin-mediated degradation

Faithful chromosome segregation requires correct attachment of kinetochores with the spindle microtubules. Erroneously-attached kinetochores recruit proteins to activate Spindle assembly checkpoint (SAC), which senses the errors and signals cells to delay anaphase progression for error correction. Temporal control of the levels of SAC activating-proteins is critical for checkpoint activation and silencing, but its mechanism is not fully understood. Here, we show that E3 ubiquitin ligase, SCF-FBXW7 targets BubR1 for ubiquitin-mediated degradation and thereby controls SAC in human cells. Depletion of FBXW7 results in prolonged metaphase arrest with increased stabilization of BubR1 at kinetochores. Similar kinetochore stabilization is also observed for BubR1-interacting protein, CENP-E. FBXW7 induced ubiquitination of both BubR1 and the BubR1-interacting kinetochore-targeting domain of CENP-E, but CENP-E domain degradation is dependent on BubR1. Interestingly, Cdk1 inhibition disrupts FBXW7-mediated BubR1 targeting and further, phospho-resistant mutation of Cdk1-targeted phosphorylation site, Thr 620 impairs BubR1-FBXW7 interaction and FBXW7-mediated BubR1 ubiquitination, supporting its role as a phosphodegron for FBXW7. The results demonstrate SCF-FBXW7 as a key regulator of spindle assembly checkpoint that controls stability of BubR1 and its associated CENP-E at kinetochores. They also support that upstream Cdk1 specific BubR1 phosphorylation signals the ligase to activate the process.

Differential gene expression profiling implicates altered network development in rat postnatal day 4 cortex following 4-Methylimidazole (4-MeI) induced maternal seizures

Compromised maternal health leading to maternal seizures can have adverse effects on the healthy development of offspring. This may be the result of inflammation, hypoxia-ischemia, and altered GABA signaling. The current study examined cortical tissue from F2b (2nd litter of the 2nd generation) postnatal day 4 (PND4) offspring of female Harlan SD rats chronically exposed to the seizuregenic compound, 4-Methylimidazole (0, 750, or 2500 ppm 4-MeI). Maternal seizures were evident only at 2500 ppm 4-MeI. GABA related gene expression as examined by qRT-PCR and whole genome microarray showed no indication of disrupted GABA or glutamatergic signaling. Canonical pathway hierarchical clustering and multi-omics combinatory genomic (CNet) plots of differentially expressed genes (DEG) showed alterations in genes associated with regulatory processes of cell development including neuronal differentiation and synaptogenesis. Functional enrichment analysis showed a similarity of cellular processes across the two exposure groups however, the genes comprising each cluster were primarily unique rather than shared and often showed different directionality. A dose-related induction of cytokine signaling was indicated however, pathways associated with individual cytokine signaling were not elevated, suggesting an alternative involvement of cytokine signaling. Pathways related to growth process and cell signaling showed a negative activation supporting an interpretation of disruption or delay in developmental processes at the 2500 ppm 4-MeI exposure level with maternal seizures. Thus, while GABA signaling was not altered as has been observed with maternal seizures, the pattern of DEG suggested a potential for alteration in neuronal network formation.

Mono-phosphorylation at Ser4 of barrier-to-autointegration factor (Banf1) significantly reduces its DNA binding capability by inducing critical changes in its local conformation and DNA binding surface

Barrier-to-autointegration factor (Banf1) is a small DNA-bridging protein. The binding status of Banf1 to DNA is regulated by its N-terminal phosphorylation and dephosphorylation, which plays a critical role in cell proliferation. Banf1 can be phosphorylated at Ser4 into mono-phosphorylated Banf1, which is further phosphorylated at Thr3 to form di-phosphorylated Banf1. It was observed decades ago that mono-phosphorylated Banf1 cannot bind to DNA. However, the underlying molecular- and atomic-level mechanisms remain unclear. A clear understanding of these mechanisms will aid in interfering with the cell proliferation process for better global health. Herein, we explored the detailed atomic bases of unphosphorylated Banf1-DNA binding and how mono- and di-phosphorylation of Banf1 impair these atomic bases to eliminate its DNA-binding capability, followed by exploring the DNA-binding capability of mono- and di-phosphorylation Banf1, using comprehensive and systematic molecular modelling and molecular dynamics simulations. This work presented in detail the residue-level binding energies, hydrogen bonds and water bridges between Banf1 and DNA, some of which have not been reported. Moreover, we revealed that mono-phosphorylation of Banf1 causes its N-terminal secondary structure changes, which in turn induce significant changes in Banf1's DNA binding surface, thus eliminating its DNA-binding capability. At the atomic level, we also uncovered the alterations in interactions due to the induction of mono-phosphorylation that result in the N-terminal secondary structure changes of Banf1. Additionally, our modelling showed that phosphorylated Banf1 with their dominant N-terminal secondary structures bind to DNA with a significantly lower affinity and the docked binding pose are not stable in MD simulations. These findings help future studies in predicting effect of mutations in Banf1 on its DNA-binding capability and open a novel avenue for the development of therapeutics such as cancer drugs, targeting cell proliferation by inducing conformational changes in Banf1's N-terminal domain.

Redefining cerebral palsies as a diverse group of neurodevelopmental disorders with genetic aetiology

Cerebral palsy is a clinical descriptor covering a diverse group of permanent, non-degenerative disorders of motor function. Around one-third of cases have now been shown to have an underlying genetic aetiology, with the genetic landscape overlapping with those of neurodevelopmental disorders including intellectual disability, epilepsy, speech and language disorders and autism. Here we review the current state of genomic testing in cerebral palsy, highlighting the benefits for personalized medicine and the imperative to consider aetiology during clinical diagnosis. With earlier clinical diagnosis now possible, we emphasize the opportunity for comprehensive and early genomic testing as a crucial component of the routine diagnostic work-up in people with cerebral palsy.

A widely used pan-isoform-FBXW7 antibody used in cell cycle studies does not detect FBXW7

FBXW7 is the substrate recognition component of the E3 ubiquitin ligase SCF^FBW7 complex which controls the levels of CYCLINE, c-MYC and HIF1α proteins crucial for cell growth and differentiation. Mutations in FBXW7 are frequently associated with tumourigenesis. While examining FBXW7 regulation we were compelled to reevaluate a commonly used anti-FBXW7 antibody. Retinal microvascular endothelial cells (RMEC) were exposed to normoxia (21% oxygen) or hypoxia (1% oxygen) for 24 h or treated with MG132 and protein extracted for western blotting. Flag-tagged FBXW7-α, β or γ isoforms were transfected into HEK293A cells and processed using denaturing and native extraction protocols for western blotting or immunoprecipitation analysis. Two anti-FBXW7 antibodies were used, one raised to the unique FBXW7α N-terminus and the other to the C-terminus region common to all isoforms. Initial studies showed that the pan-isoform C-terminus antibody detected a single 64kDa band in RMEC rather than any of the predicted sizes for FBXW7. In contrast, expression of the isoform-specific constructs, detected with an anti-Flag antibody, confirmed the expected migratory distance of 110kDa, 68kDa and 65kDa for α, β and γ respectfully. Similarly, the N-terminus FBXW7α antibody also detected the 110kDa product. Notably, the C-terminus antibody did not recognize any of the isoforms but continued to detect a 64kDa band in all samples, including the non-transfected controls. Immunoprecipitation confirmed this lack of specificity and the inability to detect overexpressed or endogenous FBXW7α in HEK293A cells and RMEC. A commonly used C-terminus FBXW7 antibody does not detect FBXW7 under standard western blotting conditions.

Neurodevelopmental disorders, like cancer, are connected to impaired chromatin remodelers, PI3K/mTOR, and PAK1-regulated MAPK

Neurodevelopmental disorders (NDDs) and cancer share proteins, pathways, and mutations. Their clinical symptoms are different. However, individuals with NDDs have higher probabilities of eventually developing cancer. Here, we review the literature and ask how the shared features can lead to different medical conditions and why having an NDD first can increase the chances of malignancy. To explore these vital questions, we focus on dysregulated PI3K/mTOR, a major brain cell growth pathway in differentiation, and MAPK, a critical pathway in proliferation, a hallmark of cancer. Differentiation is governed by chromatin organization, making aberrant chromatin remodelers highly likely agents in NDDs. Dysregulated chromatin organization and accessibility influence the lineage of specific cell brain types at specific embryonic development stages. PAK1, with pivotal roles in brain development and in cancer, also regulates MAPK. We review, clarify, and connect dysregulated pathways with dysregulated proliferation and differentiation in cancer and NDDs and highlight PAK1 role in brain development and MAPK regulation. Exactly how PAK1 activation controls brain development, and why specific chromatin remodeler components, e.g., BAF170 encoded by SMARCC2 in autism, await clarification.

DNA methylation signature of psychological resilience in young adults: Constructing a methylation risk score using a machine learning method

Resilience is a process associated with the ability to recover from stress and adversity. We aimed to explore the resilience-associated DNA methylation signatures and evaluate the abilities of methylation risk scores to discriminate low resilience (LR) individuals. The study recruited 78 young adults and used Connor-Davidson Resilience Scale (CD-RISC) to divide them into low and high resilience groups. We randomly allocated all participants of two groups to the discovery and validation sets. We used the blood DNA of the subjects to conduct a genome-wide methylation scan and identify the significant methylation differences of CpG Sites in the discovery set. Moreover, the classification accuracy of the DNA methylation probes was confirmed in the validation set by real-time quantitative methylation-specific polymerase chain reaction. In the genome-wide methylation profiling between LR and HR individuals, seventeen significantly differentially methylated probes were detected. In the validation set, nine DNA methylation signatures within gene coding regions were selected for verification. Finally, three methylation probes [cg18565204 (AARS), cg17682313 (FBXW7), and cg07167608 (LINC01107)] were included in the final model of the methylation risk score for LR versus HR. These methylation risk score models of low resilience demonstrated satisfactory discrimination by logistic regression and support vector machine, with an AUC of 0.81 and 0.93, accuracy of 72.3% and 87.1%, sensitivity of 75%, and 87.5%, and specificity of 70% and 80%. Our findings suggest that methylation signatures can be utilized to identify individuals with LR and establish risk score models that may contribute to the field of psychology.

Novel compound heterozygous variants in EMC1 associated with global developmental delay: a lesson from a non-silent synonymous exonic mutation

Background

The endoplasmic reticulum-membrane protein complex (EMC) as a molecular chaperone is required for the proper synthesis, folding and traffic of several transmembrane proteins. Variants in the subunit 1 of EMC (EMC1) have been implicated in neurodevelopmental disorders.

Methods

Whole exome sequencing (WES) with Sanger sequencing validation was performed for a Chinese family, including the proband (a 4-year-old girl who displayed global developmental delay, severe hypotonia and visual impairment), her affected younger sister and her non-consanguineous parents. RT-PCR assay and Sanger sequencing were used to detect abnormal RNA splicing.

Results

Novel compound heterozygous variants in EMC1, including the maternally inherited chr1: 19566812_1956800delinsATTCTACTT[hg19];NM_015047.3:c.765_777delins ATTCTACTT;p.(Leu256fsTer10) and the paternally inherited chr1:19549890G> A[hg19];NM_015047.3:c.2376G>A;p.(Val792=) are identified in the proband and her affected sister. RT-PCR assay followed by Sanger sequencing reveals that the c.2376G>A variant leads to aberrant splicing, with retention of intron 19 (561bp) in the mature mRNA, which is presumed to introduce a premature translational termination codon (p.(Val792fsTer31)).

Conclusion

Novel compound heterozygous variants in EMC1 have been identified in individuals with global developmental delay. Non-silent synonymous mutations should be kept in mind in genetic analysis.