RhythmicDB: A Database of Predicted Multi-Frequency Rhythmic Transcripts

The physiology and behavior of living organisms are featured by time-related variations driven by molecular clockworks that arose during evolution stochastically and heterogeneously. Over the years, several high-throughput experiments were performed to evaluate time-dependent gene expression in different cell types across several species and experimental conditions. Here, these were retrieved, manually curated, and analyzed by two software packages, BioCycle and MetaCycle, to infer circadian or ultradian transcripts across different species. These transcripts were stored in RhythmicDB and made publically available

Clinical and molecular description of the first Italian cohort of 33 subjects with hypophosphatasia

IntroductionHypophosphatasia (HPP) is a rare genetic disease caused by inactivating variants of the ALPL gene. Few data are available on the clinical presentation in Italy and/or on Italian HPP surveys.MethodsThere were 30 suspected HPP patients recruited from different Italian tertiary cares. Biological samples and related clinical, biochemical, and anamnestic data were collected and the ALPL gene sequenced. Search for large genomic deletions at the ALPL locus (1p36) was done. Phylogenetic conservation and modeling were applied to infer the effect of the variants on the protein structure.ResultsThere were 21 ALPL variants and one large genomic deletion found in 20 out of 30 patients. Unexpectedly, NGS-driven differential diagnosis allowed uncovering three hidden additional HPP cases, for a total of 33 HPP subjects. Eight out of 24 coding variants were novel and classified as “pathogenic”, “likely pathogenic”, and “variants of uncertain significance”. Bioinformatic analysis confirmed that all the variants strongly destabilize the homodimer structure. There were 10 cases with low ALP and high VitB6 that resulted negative to genetic testing, whereas two positive cases have an unexpected normal ALP value. No association was evident with other biochemical/clinical parameters.DiscussionWe present the survey of HPP Italian patients with the highest ALPL mutation rate so far reported and confirm the complexity of a prompt recognition of the syndrome, mostly for HPP in adults. Low ALP and high VitB6 values are mandatory for the genetic screening, this latter remaining the gold standard not only to confirm the clinical diagnosis but also to make differential diagnosis, to identify carriers, to avoid likely dangerous therapy in unrecognized cases

Phosphorylation within Intrinsic Disordered Region Discriminates Histone Variant macroH2A1 Splicing Isoforms-macroH2A1.1 and macroH2A1.2

Background: Gene expression in eukaryotic cells can be governed by histone variants, which replace replication-coupled histones, conferring unique chromatin properties. MacroH2A1 is a histone H2A variant containing a domain highly similar to H2A and a large non-histone (macro) domain. MacroH2A1, in turn, is present in two alternatively exon-spliced isoforms: macroH2A1.1 and macroH2A1.2, which regulate cell plasticity and proliferation in a remarkably distinct manner. The N-terminal and the C-terminal tails of H2A histones stem from the nucleosome core structure and can be target sites for several post-translational modifications (PTMs). MacroH2A1.1 and macroH2A1.2 isoforms differ only in a few amino acids and their ability to bind NAD-derived metabolites, a property allegedly conferring their different functions in vivo. Some of the modifications on the macroH2A1 variant have been identified, such as phosphorylation (T129, S138) and methylation (K18, K123, K239). However, no study to our knowledge has analyzed extensively, and in parallel, the PTM pattern of macroH2A1.1 and macroH2A1.2 in the same experimental setting, which could facilitate the understanding of their distinct biological functions in health and disease. Methods: We used a mass spectrometry-based approach to identify the sites for phosphorylation, acetylation, and methylation in green fluorescent protein (GFP)-tagged macroH2A1.1 and macroH2A1.2 expressed in human hepatoma cells. The impact of selected PTMs on macroH2A1.1 and macroH2A1.2 structure and function are demonstrated using computational analyses. Results: We identified K7 as a new acetylation site in both macroH2A1 isoforms. Quantitative comparison of histone marks between the two isoforms revealed significant differences in the levels of phosphorylated T129 and S170. Our computational analysis provided evidence that the phosphorylation status in the intrinsically disordered linker region in macroH2A1 isoforms might represent a key regulatory element contributing to their distinct biological responses. Conclusions: Taken together, our results report different PTMs on the two macroH2A1 splicing isoforms as responsible for their distinct features and distribution in the cell

Are Gaming-Enabled Graphic Processing Unit Cards Convenient for Molecular Dynamics Simulation?

In several fields of research, molecular dynamics simulation techniques are exploited to evaluate the temporal motion of particles constituting water, ions, small molecules, macromolecules, or more complex systems over time. These techniques are considered difficult to setup, computationally demanding and require high specialization and scientific skills. Moreover, they need specialized computing infrastructures to run faster and make the simulation of big systems feasible. Here, we have simulated 3 systems of increasing sizes on scientific- and gaming-enabled graphic processing unit (GPU) cards with Amber, GROMACS, and NAMD and measured their performance accounting also for the market prices of the GPU cards where they were run on

Investigating Mitochondrial Gene Expression Patterns in <i>Drosophila melanogaster</i> Using Network Analysis to Understand Aging Mechanisms

The process of aging is a complex phenomenon that involves a progressive decline in physiological functions required for survival and fertility. To better understand the mechanisms underlying this process, the scientific community has utilized several tools. Among them, mitochondrial DNA has emerged as a crucial factor in biological aging, given that mitochondrial dysfunction is thought to significantly contribute to this phenomenon. Additionally, Drosophila melanogaster has proven to be a valuable model organism for studying aging due to its low cost, capacity to generate large populations, and ease of genetic manipulation and tissue dissection. Moreover, graph theory has been employed to understand the dynamic changes in gene expression patterns associated with aging and to investigate the interactions between aging and aging-related diseases. In this study, we have integrated these approaches to examine the patterns of gene co-expression in Drosophila melanogaster at various stages of development. By applying graph-theory techniques, we have identified modules of co-expressing genes, highlighting those that contain a significantly high number of mitochondrial genes. We found important mitochondrial genes involved in aging and age-related diseases in Drosophila melanogaster, including UQCR-C1, ND-B17.2, ND-20, and Pdhb. Our findings shed light on the role of mitochondrial genes in the aging process and demonstrate the utility of Drosophila melanogaster as a model organism and graph theory in aging research

Potassium Channel {KCNH}1 Activating Variants Cause Altered Functional and Morphological Ciliogenesis

The primary cilium is a non-motile sensory organelle that extends from the surface of most vertebrate cells and transduces signals regulating proliferation, differentiation, and migration. Primary cilia dysfunctions have been observed in cancer and in a group of heterogeneous disorders called ciliopathies, characterized by renal and liver cysts, skeleton and limb abnormalities, retinal degeneration, intellectual disability, ataxia, and heart disease and, recently, in autism spectrum disorder, schizophrenia, and epilepsy. The potassium voltage-gated channel subfamily H member 1 (KCNH1) gene encodes a member of the EAG (ether-a-go-go) family, which controls potassium flux regulating resting membrane potential in both excitable and non-excitable cells and is involved in intracellular signaling, cell proliferation, and tumorigenesis. KCNH1 missense variants have been associated with syndromic neurodevelopmental disorders, including Zimmermann-Laband syndrome 1 (ZLS1, MIM #135500), Temple-Baraitser syndrome (TMBTS, MIM #611816), and, recently, with milder phenotypes as epilepsy. In this work, we provide evidence that KCNH1 localizes at the base of the cilium in pre-ciliary vesicles and ciliary pocket of human dermal fibroblasts and retinal pigment epithelial (hTERT RPE1) cells and that the pathogenic missense variants (L352V and R330Q; NP_002229.1) perturb cilia morphology, assembly/disassembly, and Sonic Hedgehog signaling, disclosing a multifaceted role of the protein. The study of KCNH1 localization, its functions related to primary cilia, and the alterations introduced by mutations in ciliogenesis, cell cycle coordination, cilium morphology, and cilia signaling pathways could help elucidate the molecular mechanisms underlying neurological phenotypes and neurodevelopmental disorders not considered as classical ciliopathies but for which a significant role of primary cilia is emerging

Heterozygous nonsense ARX mutation in a family highlights the complexity of clinical and molecular diagnosis in case of chromosomal and single gene disorder co-inheritance

Background: Corpus callosum agenesis (ACC) is one of the most frequent Central Nervous System (CNS) malformations. However, genetics underlying isolated forms is still poorly recognized. Here, we report on two female familial cases with partial ACC. The proband shows isolated partial ACC and a mild neurodevelopmental phenotype. A fetus from a previous interrupted pregnancy exhibited a complex phenotype including partial ACC and the occurrence of a de novo 17q12 microduplication, which was interpreted as probably disease-causing. Methods: A trio-based clinical exome sequencing (CES) was performed. Results: Clinical exome sequencing data analysis led to identifying a heterozygous nonsense variant (NM_139058.3:c.922G&gt;T; NP_620689.1:p.Glu308Ter) in the aristaless related homeobox gene (ARX) in the proband, with a putative de novo occurrence, producing a hypothetical protein lacking two essential domains. Sanger analysis confirmed the wild-type status of both parents in different tissues, and disclosed the occurrence of the nonsense variant in the fetus of the interrupted pregnancy, suggesting a formerly unrecognized contribution of the ARX mutation to the fetus' phenotype and gonadal or gonadosomatic mosaicism in one of the parents. Conclusion: This study describes the phenotype associated with a heterozygous loss of function variant in ARX. Moreover, it highlights the importance of investigating both chromosomal and genetic contributions in cases of complex syndromic phenotypes involving CNS

APOGEE 2: multi-layer machine-learning model for the interpretable prediction of mitochondrial missense variants

Abstract Mitochondrial dysfunction has pleiotropic effects and is frequently caused by mitochondrial DNA mutations. However, factors such as significant variability in clinical manifestations make interpreting the pathogenicity of variants in the mitochondrial genome challenging. Here, we present APOGEE 2, a mitochondrially-centered ensemble method designed to improve the accuracy of pathogenicity predictions for interpreting missense mitochondrial variants. Built on the joint consensus recommendations by the American College of Medical Genetics and Genomics/Association for Molecular Pathology, APOGEE 2 features an improved machine learning method and a curated training set for enhanced performance metrics. It offers region-wise assessments of genome fragility and mechanistic analyses of specific amino acids that cause perceptible long-range effects on protein structure. With clinical and research use in mind, APOGEE 2 scores and pathogenicity probabilities are precompiled and available in MitImpact. APOGEE 2’s ability to address challenges in interpreting mitochondrial missense variants makes it an essential tool in the field of mitochondrial genetics

Molecular dynamics analysis of superoxide dismutase 1 mutations suggests decoupling between mechanisms underlying ALS onset and progression

Mutations in the superoxide dismutase 1 (SOD1) gene are the second most common known cause of ALS. SOD1 variants express high phenotypic variability and over 200 have been reported in people with ALS. It was previously proposed that variants can be broadly classified in two groups, ‘wild-type like’ (WTL) and ‘metal binding region’ (MBR) variants, based on their structural location and biophysical properties. MBR variants, but not WTL variants, were associated with a reduction of SOD1 enzymatic activity. In this study we used molecular dynamics and large clinical datasets to characterise the differences in the structural and dynamic behaviour of WTL and MBR variants with respect to the wild-type SOD1, and how such differences influence the ALS clinical phenotype. Our study identified marked structural differences, some of which are observed in both variant groups, while others are group specific. Moreover, collecting clinical data of approximately 500 SOD1 ALS patients carrying variants, we showed that the survival time of patients carrying an MBR variant is generally longer (∼6 years median difference, p < 0.001) with respect to patients with a WTL variant. In conclusion, our study highlighted key differences in the dynamic behaviour between WTL and MBR SOD1 variants, and between variants and wild-type SOD1 at an atomic and molecular level, that could be further investigated to explain the associated phenotypic variability. Our results support the hypothesis of a decoupling between mechanisms of onset and progression of SOD1 ALS, and an involvement of loss-of-function of SOD1 with the disease progression

Clinical variability in DYNC2H1-related skeletal ciliopathies includes Ellis-van Creveld syndrome

: Deleterious variants of DYNC2H1 gene are associated with a wide spectrum of skeletal ciliopathies (SC). We used targeted parallel sequencing to analyze 25 molecularly unsolved families with different SCs. Deleterious DYNC2H1 variants were found in six sporadic patients and two monozygotic (MZ) twins. Clinical diagnoses included short rib-polydactyly type 3 in two cases, and asphyxiating thoracic dystrophy (ATD) in one case. Remarkably, clinical diagnosis fitted with EvC, mixed ATD/EvC and short rib-polydactyly/EvC phenotypes in three sporadic patients and the MZ twins. EvC/EvC-like features always occurred in compound heterozygotes sharing a previously unreported splice site change (c.6140-5A>G) or compound heterozygotes for two missense variants. These results expand the DYNC2H1 mutational repertoire and its clinical spectrum, suggesting that EvC may be occasionally caused by DYNC2H1 variants presumably acting as hypomorphic alleles