De novo variants in ATP2B1 lead to neurodevelopmental delay

Heterozygous loss-of-function variants in DOCK4 cause neurodevelopmental delay and microcephaly

Neurons form the basic anatomical and functional structure of the nervous system, and defects in neuronal differentiation or formation of neurites are associated with various psychiatric and neurodevelopmental disorders. Dynamic changes in the cytoskeleton are essential for this process, which is, inter alia, controlled by the dedicator of cytokinesis 4 (DOCK4) through the activation of RAC1. Here, we clinically describe 7 individuals (6 males and one female) with variants in DOCK4 and overlapping phenotype of mild to severe global developmental delay. Additional symptoms include coordination or gait abnormalities, microcephaly, nonspecific brain malformations, hypotonia and seizures. Four individuals carry missense variants (three of them detected de novo) and three individuals carry null variants (two of them maternally inherited). Molecular modeling of the heterozygous missense variants suggests that the majority of them affect the globular structure of DOCK4. In vitro functional expression studies in transfected Neuro-2A cells showed that all missense variants impaired neurite outgrowth. Furthermore, Dock4 knockout Neuro-2A cells also exhibited defects in promoting neurite outgrowth. Our results, including clinical, molecular and functional data, suggest that loss-of-function variants in DOCK4 probable cause a variable spectrum of a novel neurodevelopmental disorder with microcephaly.

Protein-protein interaction network-based integration of GWAS and functional data for blood pressure regulation analysis

BACKGROUND

It is valuable to analyze the genome-wide association studies (GWAS) data for a complex disease phenotype in the context of the protein-protein interaction (PPI) network, as the related pathophysiology results from the function of interacting polyprotein pathways. The analysis may include the design and curation of a phenotype-specific GWAS meta-database incorporating genotypic and eQTL data linking to PPI and other biological datasets, and the development of systematic workflows for PPI network-based data integration toward protein and pathway prioritization. Here, we pursued this analysis for blood pressure (BP) regulation.

METHODS

The relational scheme of the implemented in Microsoft SQL Server BP-GWAS meta-database enabled the combined storage of: GWAS data and attributes mined from GWAS Catalog and the literature, Ensembl-defined SNP-transcript associations, and GTEx eQTL data. The BP-protein interactome was reconstructed from the PICKLE PPI meta-database, extending the GWAS-deduced network with the shortest paths connecting all GWAS-proteins into one component. The shortest-path intermediates were considered as BP-related. For protein prioritization, we combined a new integrated GWAS-based scoring scheme with two network-based criteria: one considering the protein role in the reconstructed by shortest-path (RbSP) interactome and one novel promoting the common neighbors of GWAS-prioritized proteins. Prioritized proteins were ranked by the number of satisfied criteria.

RESULTS

The meta-database includes 6687 variants linked with 1167 BP-associated protein-coding genes. The GWAS-deduced PPI network includes 1065 proteins, with 672 forming a connected component. The RbSP interactome contains 1443 additional, network-deduced proteins and indicated that essentially all BP-GWAS proteins are at most second neighbors. The prioritized BP-protein set was derived from the union of the most BP-significant by any of the GWAS-based or the network-based criteria. It included 335 proteins, with ~ 2/3 deduced from the BP PPI network extension and 126 prioritized by at least two criteria. ESR1 was the only protein satisfying all three criteria, followed in the top-10 by INSR, PTN11, CDK6, CSK, NOS3, SH2B3, ATP2B1, FES and FINC, satisfying two. Pathway analysis of the RbSP interactome revealed numerous bioprocesses, which are indeed functionally supported as BP-associated, extending our understanding about BP regulation.

CONCLUSIONS

The implemented workflow could be used for other multifactorial diseases.

Biallelic ATP2B1 variants as a likely cause of a novel neurodevelopmental malformation syndrome with primary hypoparathyroidism

ATP2B1 encodes plasma membrane calcium-transporting-ATPase1 and plays an essential role in maintaining intracellular calcium homeostasis that regulates diverse signaling pathways. Heterozygous de novo missense and truncating ATP2B1 variants are associated with a neurodevelopmental phenotype of variable expressivity. We describe a proband with distinctive craniofacial gestalt, Pierre-Robin sequence, neurodevelopmental and growth deficit, periventricular heterotopia, brachymesophalangy, cutaneous syndactyly, and persistent hypocalcemia from primary hypoparathyroidism. Proband-parent trio exome sequencing identified compound heterozygous ATP2B1 variants: a maternally inherited splice-site (c.3060+2 T > G) and paternally inherited missense c.2938 G > T; p.(Val980Leu). Reverse-transcription-PCR on the proband's fibroblast-derived mRNA showed aberrantly spliced ATP2B1 transcripts targeted for nonsense-mediated decay. All correctly-spliced ATP2B1 mRNA encoding p.(Val980Leu) functionally causes decreased cellular Ca2+ extrusion. Immunoblotting showed reduced fibroblast ATP2B1. We conclude that biallelic ATP2B1 variants are the likely cause of the proband's phenotype, strengthening the association of ATP2B1 as a neurodevelopmental gene and expanding the phenotypic characterization of a biallelic loss-of-function genotype.

Evaluation of input data modality choices on functional gene embeddings

Functional gene embeddings, numerical vectors capturing gene function, provide a promising way to integrate functional gene information into machine learning models. These embeddings are learnt by applying self-supervised machine-learning algorithms on various data types including quantitative omics measurements, protein-protein interaction networks and literature. However, downstream evaluations comparing alternative data modalities used to construct functional gene embeddings have been lacking. Here we benchmarked functional gene embeddings obtained from various data modalities for predicting disease-gene lists, cancer drivers, phenotype-gene associations and scores from genome-wide association studies. Off-the-shelf predictors trained on precomputed embeddings matched or outperformed dedicated state-of-the-art predictors, demonstrating their high utility. Embeddings based on literature and protein-protein interactions inferred from low-throughput experiments outperformed embeddings derived from genome-wide experimental data (transcriptomics, deletion screens and protein sequence) when predicting curated gene lists. In contrast, they did not perform better when predicting genome-wide association signals and were biased towards highly-studied genes. These results indicate that embeddings derived from literature and low-throughput experiments appear favourable in many existing benchmarks because they are biased towards well-studied genes and should therefore be considered with caution. Altogether, our study and precomputed embeddings will facilitate the development of machine-learning models in genetics and related fields.

ATP2B2 de novo variants as a cause of variable neurodevelopmental disorders that feature dystonia, ataxia, intellectual disability, behavioral symptoms, and seizures

PURPOSE

ATP2B2 encodes the variant-constrained plasma-membrane calcium-transporting ATPase-2, expressed in sensory ear cells and specialized neurons. ATP2B2/Atp2b2 variants were previously linked to isolated hearing loss in patients and neurodevelopmental deficits with ataxia in mice. We aimed to establish the association between ATP2B2 and human neurological disorders.

METHODS

Multinational case recruitment, scrutiny of trio-based genomics data, in silico analyses, and functional variant characterization were performed.

RESULTS

We assembled 7 individuals harboring rare, predicted deleterious heterozygous ATP2B2 variants. The alleles comprised 5 missense substitutions that affected evolutionarily conserved sites and 2 frameshift variants in the penultimate exon. For 6 variants, a de novo status was confirmed. Unlike described patients with hearing loss, the individuals displayed a spectrum of neurological abnormalities, ranging from ataxia with dystonic features to complex neurodevelopmental manifestations with intellectual disability, autism, and seizures. Two cases with recurrent amino-acid variation showed distinctive overlap with cerebellar atrophy-associated ataxia and epilepsy. In cell-based studies, all variants caused significant alterations in cytosolic calcium handling with both loss- and gain-of-function effects.

CONCLUSION

Presentations in our series recapitulate key phenotypic aspects of Atp2b2-mouse models and underline the importance of precise calcium regulation for neurodevelopment and cerebellar function. Our study documents a role for ATP2B2 variants in causing heterogeneous neurodevelopmental and movement-disorder syndromes.

DNA Methylation of Birthweight-Blood Pressure Genes and Changes of Blood Pressure in Response to Weight-Loss Diets in the POUNDS Lost Trial

BACKGROUND

DNA methylation (DNAm) may play a critical role in bridging prenatal adverse events and cardiometabolic disorders including hypertension in later life.

METHODS

We included 672 adult participants with overweight or obesity, who participated in a 2-year randomized weight-loss dietary intervention study. We defined the regional DNAm levels as the average methylation level of 5'-cytosine-phosphate-guanine-3' within 500 bp of LINC00319 (cg01820192), ATP2B1 (cg00508575), and LMNA (cg12593793), respectively. Generalized linear regression models were used to assess the association between the regional DNAm and 2-year blood pressure changes. Trajectory analysis was used to identify subgroups that shared a similar underlying trajectory of 2-year blood pressure changes.

RESULTS

The regional DNAm at LINC00319, showed significantly different associations with 2-year changes in systolic blood pressure and diastolic blood pressure among participants assigned to low- or high-fat diets (P for interaction<0.05 for all). In response to the low-fat diet, per SD higher regional DNAm at LINC00319 was associated with greater reductions in both 2-year changes in systolic blood pressure (β, -1.481; P=0.020) and diastolic blood pressure (β, -1.096; P=0.009). Three trajectories of changes in systolic blood pressure or diastolic blood pressure were identified, and participants with higher regional DNAm at LINC00319 were more likely to experience and maintain decreased systolic blood pressure and diastolic blood pressure (odds ratio of being in decrease-stable versus stable [95% CI], 1.542 [1.146-2.076] and 1.463 [1.125-1.902]).

CONCLUSIONS

Our findings suggest that DNAm could be a metabolic memory bridging early and later life, and an indicator of more benefits from eating a low-fat weight-loss diet.

Hypoxic Preconditioned Neural Stem Cell-Derived Extracellular Vesicles Contain Distinct Protein Cargo from Their Normal Counterparts

Hypoxic preconditioning has been demonstrated to increase the resistance of neural stem cells (NSCs) to hypoxic conditions, as well as to improve their capacity for differentiation and neurogenesis. Extracellular vesicles (EVs) have recently emerged as critical mediators of cell–cell communication, but their role in this hypoxic conditioning is presently unknown. Here, we demonstrated that three hours of hypoxic preconditioning triggers significant neural stem cell EV release. Proteomic profiling of EVs from normal and hypoxic preconditioned neural stem cells identified 20 proteins that were upregulated and 22 proteins that were downregulated after hypoxic preconditioning. We also found an upregulation of some of these proteins by qPCR, thus indicating differences also at the transcript level within the EVs. Among the upregulated proteins are CNP, Cyfip1, CASK, and TUBB5, which are well known to exhibit significant beneficial effects on neural stem cells. Thus, our results not only show a significant difference of protein cargo in EVs consequent to hypoxic exposure, but identify several candidate proteins that might play a pivotal role in the cell-to-cell mediated communication underlying neuronal differentiation, protection, maturation, and survival following exposure to hypoxic conditions.

Calcium signaling and genetic rare diseases: An auditory perspective

Deafness is a highly heterogeneous disorder which stems, for 50%, from genetic origins. Sensory transduction relies mainly on sensory hair cells of the cochlea, in the inner ear. Calcium is key for the function of these cells and acts as a fundamental signal transduction. Its homeostasis depends on three factors: the calcium influx, through the mechanotransduction channel at the apical pole of the hair cell as well as the voltage-gated calcium channel at the base of the cells; the calcium buffering via Ca²⁺-binding proteins in the cytoplasm, but also in organelles such as mitochondria and the reticulum endoplasmic mitochondria-associated membranes with specialized proteins; and the calcium extrusion through the Ca-ATPase pump, located all over the plasma membrane. In addition, the synaptic transmission to the central nervous system is also controlled by calcium. Genetic studies of inherited deafness have tremendously helped understand the underlying molecular pathways of calcium signaling. In this review, we discuss these different factors in light of the associated genetic diseases (syndromic and non-syndromic deafness) and the causative genes.

Ca2+ Homeostasis by Plasma Membrane Ca2+ ATPase (PMCA) 1 Is Essential for the Development of DP Thymocytes

The strength of Ca²⁺ signaling is a hallmark of T cell activation, yet the role of Ca²⁺ homeostasis in developing T cells before expressing a mature T cell receptor is poorly understood. We aimed to unveil specific functions of the two plasma membrane Ca²⁺ ATPases expressed in T cells, PMCA1 and PMCA4. On a transcriptional and protein level we found that PMCA4 was expressed at low levels in CD4^-CD8^- double negative (DN) thymocytes and was even downregulated in subsequent stages while PMCA1 was present throughout development and upregulated in CD4⁺CD8⁺ double positive (DP) thymocytes. Mice with a targeted deletion of Pmca1 in DN3 thymocytes had an almost complete block of DP thymocyte development with an accumulation of DN4 thymocytes but severely reduced numbers of CD8⁺ immature single positive (ISP) thymocytes. The DN4 thymocytes of these mice showed strongly elevated basal cytosolic Ca²⁺ levels and a pre-mature CD5 expression, but in contrast to the DP thymocytes they were only mildly prone to apoptosis. Surprisingly, mice with a germline deletion of Pmca4 did not show any signs of altered progression through the developmental thymocyte stages, nor altered Ca²⁺ homeostasis throughout this process. PMCA1 is, therefore, non-redundant in keeping cellular Ca²⁺ levels low in the early thymocyte development required for the DN to DP transition.

The ataxia-linked E1081Q mutation affects the sub-plasma membrane Ca2+-microdomains by tuning PMCA3 activity

Calcium concentration must be finely tuned in all eukaryotic cells to ensure the correct performance of its signalling function. Neuronal activity is exquisitely dependent on the control of Ca²⁺ homeostasis: its alterations ultimately play a pivotal role in the origin and progression of many neurodegenerative processes. A complex toolkit of Ca²⁺ pumps and exchangers maintains the fluctuation of cytosolic Ca²⁺ concentration within the appropriate threshold. Two ubiquitous (isoforms 1 and 4) and two neuronally enriched (isoforms 2 and 3) of the plasma membrane Ca²⁺ATPase (PMCA pump) selectively regulate cytosolic Ca²⁺ transients by shaping the sub-plasma membrane (PM) microdomains. In humans, genetic mutations in ATP2B1, ATP2B2 and ATP2B3 gene have been linked with hearing loss, cerebellar ataxia and global neurodevelopmental delay: all of them were found to impair pump activity. Here we report three additional mutations in ATP2B3 gene corresponding to E1081Q, R1133Q and R696H amino acids substitution, respectively. Among them, the novel missense mutation (E1081Q) immediately upstream the C-terminal calmodulin-binding domain (CaM-BD) of the PMCA3 protein was present in two patients originating from two distinct families. Our biochemical and molecular studies on PMCA3 E1081Q mutant have revealed a splicing variant-dependent effect of the mutation in shaping the sub-PM [Ca²⁺]. The E1081Q substitution in the full-length b variant abolished the capacity of the pump to reduce [Ca²⁺] in the sub-PM microdomain (in line with the previously described ataxia-related PMCA mutations negatively affecting Ca²⁺ pumping activity), while, surprisingly, its introduction in the truncated a variant selectively increased Ca²⁺ extrusion activity in the sub-PM Ca²⁺ microdomains. These results highlight the importance to set a precise threshold of [Ca²⁺] by fine-tuning the sub-PM microdomains and the different contribution of the PMCA splice variants in this regulation.