Disease spreading modeling and analysis: a survey.

MOTIVATION: The control of the diffusion of diseases is a critical subject of a broad research area, which involves both clinical and political aspects. It makes wide use of computational tools, such as ordinary differential equations, stochastic simulation frameworks and graph theory, and interaction data, from molecular to social granularity levels, to model the ways diseases arise and spread. The coronavirus disease 2019 (COVID-19) is a perfect testbench example to show how these models may help avoid severe lockdown by suggesting, for instance, the best strategies of vaccine prioritization. RESULTS: Here, we focus on and discuss some graph-based epidemiological models and show how their use may significantly improve the disease spreading control. We offer some examples related to the recent COVID-19 pandemic and discuss how to generalize them to other diseases.

Usmani-Riazuddin syndrome can have a recognizable phenotype: Report of a novel AP1G1 variant.

Usmani-Riazuddin syndrome (USRISR, MIM# 619548; USRISD, MIM#619467) is a very rare genetic condition. recently associated with deleterious variants in AP1G1 (MIM* 603533). It is characterized by multisystemic involvement including intellectual disability, speech and developmental delay, behavioral anomalies, muscular tone disorders, seizures, limb defects, and unspecified facial gestalt. In this report, we describe this syndrome for the second time, in association to a novel AP1G1 variant identified in a toddler with multisystemic involvement including intellectual disability, speech and developmental delay, behavioral anomalies, arrhythmias, hearing loss, skin changes, and limb defects. Next generation sequencing (NGS) analysis through clinical exome disclosed AP1G1: c.1969C>G (p.Leu657Val), de novo, likely pathogenic variant, according to ACMG classification criteria. Proband's facial features resembled the spectrum of chromatinopathies. Clinical pictures were analyzed and a clinical overlap was supported by DeepGestalt analysis (www.face2gene.com). The system identified 6 chromatin disorders out of 30 possible diagnoses. The remaining 24 included 9 miscellaneous cryptic chromosomal abnormalities (excluded due to normal microarray study). To the best of our knowledge, this is the first description of likely distinctive facial features in a patient with Usmani-Riazuddin syndrome. Further multicentric analyses are needed for a better definition of this aspect.

A comparative benchmark of classic DNA motif discovery tools on synthetic data.

Hundreds of human proteins were found to establish transient interactions with rather degenerated consensus DNA sequences or motifs. Identifying these motifs and the genomic sites where interactions occur represent one of the most challenging research goals in modern molecular biology and bioinformatics. The last twenty years witnessed an explosion of computational tools designed to perform this task, whose performance has been last compared fifteen years ago. Here, we survey sixteen of them, benchmark their ability to identify known motifs nested in twenty-nine simulated sequence datasets, and finally report their strengths, weaknesses, and complementarity.

MitImpact 3: modeling the residue interaction network of the Respiratory Chain subunits.

Numerous lines of evidence have shown that the interaction between the nuclear and mitochondrial genomes ensures the efficient functioning of the OXPHOS complexes, with substantial implications in bioenergetics, adaptation, and disease. Their interaction is a fascinating and complex trait of the eukaryotic cell that MitImpact explores with its third major release. MitImpact expands its collection of genomic, clinical, and functional annotations of all non-synonymous substitutions of the human mitochondrial genome with new information on putative Compensated Pathogenic Deviations and co-varying amino acid sites of the Respiratory Chain subunits. It further provides evidence of energetic and structural residue compensation by techniques of molecular dynamics simulation. MitImpact is freely accessible at http://mitimpact.css-mendel.it.

Clinical and Molecular Traits of a Novel SPECC1L-ALK Fusion in a Patient with Advanced Non-Small Cell Lung Cancer.

Anaplastic lymphoma kinase (ALK) fusions account for 5-7% of non-small cell lung cancer (NSCLC) patients, the therapeutic approaches for which have significantly evolved in the last few years. However, the response to target therapies remains heterogeneous, partially due to the many different ALK fusion variants reported to date. Rare fusion variants have also been discovered, but their role in influencing responses to ALK inhibitors (ALKis) remains poorly elucidated. Laboratory investigation at both the tissue and protein levels, and a molecular profile by next-generation sequencing (NGS) were performed on a lung biopsy sample from one patient with poorly differentiated adenocarcinoma. An in silico prediction model using ColabFold software v1.5.5 was used to model and predict the entire structure of the chimeric aberrant ALK protein. Here, we report a case of a former smoker, a 60-year-old man, diagnosed with NSCLC and undergoing ALK translocation. He received first-, second- and third-generation ALK protein inhibitors (ALKis), showing a clinical benefit for about 4 years. IHC analysis and the molecular examination of the tissue biopsy indicated a positive staining for ALK and a novel ALK gene fusion variant, involving the sperm antigen with calponin homology and coiled-coil domain 1-like (SPECC1L) gene with an unreported breakpoint in exon 7. The novel SPECCL1::ALK fusion was identified using Anchored Multiplex PCR (AMP)-NGS technology and was predicted to retain the Pkinase_Tyr domain at the carboxy-terminal region of the resulting chimeric protein. To the best of our knowledge, this is the first case of an ALK-positive NSCLC patient carrying the SPECC1L exon 7 fusion breakpoint and one of the few reports about clinical outcomes related to SPECC1L::ALK fusion. The in silico hypothesized biological role of this newly identified fusion variant corroborates the observed clinical response to multiple ALKis. The molecular findings also reinforce the utility of AMP-based NGS technology as a valuable tool for the identification of rare chromosomal events that may be related to the variability of patient outcomes to different ALKis treatments.

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

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.

Exploiting in silico structural analysis to introduce emerging genotype-phenotype correlations in DHCR24-related sterol biosynthesis disorder: a case study.

Desmosterolosis is a rare sterol biosynthesis disorder characterized by multiple congenital anomalies, failure to thrive, severe developmental delay, progressive epileptic encephalopathy, and elevated levels of desmosterol caused by biallelic mutations of DHCR24 encoding 3-ß-hydroxysterol Δ-24-reductase. DHCR24 is regarded as the key enzyme of cholesterol synthesis in the metabolism of brain cholesterol as it catalyzes the reduction of the Δ-24 double bond of sterol intermediates during cholesterol biosynthesis. To date, 15 DHCR24 variants, detected in 2 related and 14 unrelated patients, have been associated with the desmosterolosis disorder. Here, we describe a proband harboring the never-described DHCR24 homozygous missense variant NM_014762.4:c.506T>C, NP_055577.1:p.M169T, whose functional validation was confirmed through biochemical assay. By using molecular dynamics simulation techniques, we investigated the impact of this variant on the protein stability and interaction network with the flavin adenine dinucleotide cofactor, thereby providing a preliminary assessment of its mechanistic role in comparison to all known pathogenic variants, the wild-type protein, and a known benign DHCR24 variant. This report expands the clinical and molecular spectra of the DHCR24-related disorder, reports on a novel DHCR24 deleterious variant associated with desmosterolosis, and gives new insights into genotype-phenotype correlations.

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.

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

Introduction: Hypophosphatasia (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. Methods: There 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. Results: There 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. Discussion: We 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.

Beyond COVID-19 pandemic: Topology-aware optimization of vaccination strategy for minimizing virus spreading.

The mitigation of an infectious disease spreading has recently gained considerable attention from the research community. It may be obtained by adopting sanitary measurements (e.g., vaccination, wearing masks), social rules (e.g., social distancing), together with an extensive vaccination campaign. Vaccination is currently the primary way for mitigating the Coronavirus Disease (COVID-19) outbreak without severe lockdown. Its effectiveness also depends on the number and timeliness of administrations and thus demands strict prioritization criteria. Almost all countries have prioritized similar classes of exposed workers: healthcare professionals and the elderly, obtaining to maximize the survival of patients and years of life saved. Nevertheless, the virus is currently spreading at high rates, and any prioritization criterion so far adopted did not account for the structural organization of the contact networks. We reckon that a network where nodes are people while the edges represent their social contacts may efficiently model the virus's spreading. It is known that tailored interventions (e.g., vaccination) on central nodes may efficiently stop the propagation, thereby eliminating the "bridge edges." We then introduce such a model and consider both synthetic and real datasets. We present the benefits of a topology-aware versus an age-based vaccination strategy to mitigate the spreading of the virus. The code is available at https://github.com/mazzalab/playgrounds.

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.

Connecting the dots: A practical evaluation of web-tools for describing protein dynamics as networks.

Protein Structure Networks (PSNs) are a well-known mathematical model for estimation and analysis of the three-dimensional protein structure. Investigating the topological architecture of PSNs may help identify the crucial amino acid residues for protein stability and protein-protein interactions, as well as deduce any possible mutational effects. But because proteins go through conformational changes to give rise to essential biological functions, this has to be done dynamically over time. The most effective method to describe protein dynamics is molecular dynamics simulation, with the most popular software programs for manipulating simulations to infer interaction networks being RING, MD-TASK, and NAPS. Here, we compare the computational approaches used by these three tools-all of which are accessible as web servers-to understand the pathogenicity of missense mutations and talk about their potential applications as well as their advantages and disadvantages.

KDM6A missense variants hamper H3 histone demethylation in lung squamous cell carcinoma.

KDM6A is the disease causative gene of type 2 Kabuki Syndrome, a rare multisystem disease; it is also a known cancer driver gene, with multiple somatic mutations found in a few cancer types. In this study, we looked at eleven missense variants in lung squamous cell carcinoma, one of the most common lung cancer subtypes, to see how they affect the KDM6A catalytic mechanisms. We found that they influence the interaction with histone H3 and the exposure of the trimethylated Lys27, which is critical for wild-type physiological function to varying degrees, by altering the conformational transition.

Potassium Channel KCNH1 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-à-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.

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.

GDF5 mutation case report and a systematic review of molecular and clinical spectrum: Expanding current knowledge on genotype-phenotype correlations.

INTRODUCTION: Brachydactyly is a bone development abnormality presenting with variable phenotypes and different transmission patterns. Mutations in GDF5 (Growth and Differentiation Factor 5, MIM *601146) account for a significant amount of cases. Here, we report on a three-generation family, where the proband and the grandfather have an isolated brachydactyly with features of both type A1 (MIM #112500) and type C (MIM #113100), while the mother shows only subtle hand phenotype signs. MATERIALS AND METHODS: Whole Exome Sequencing (WES) was performed on the two affected individuals. An in-depth analysis of GDF5 genotype-phenotype correlations was performed through literature reviewing and retrieving information from several databases to elucidate GDF5-related molecular pathogenic mechanisms. RESULTS: WES analysis disclosed a pathogenic variant in GDF5 (NM_000557.5:c.157dup; NP_000548.2:p.Leu53Profs*41; rs778834209), segregating with the phenotype. The frameshift variant was previously associated with Brachydactyly type C (MIM #113100), in heterozygosity, and with the severe Grebe type chondrodysplasia (MIM #200700), in homozygosity. In-depth analysis of literature and databases allowed to retrieve GDF5 mutations and correlations to phenotypes. We disclosed the association of 49 GDF5 pathogenic mutations with eight phenotypes, with both autosomal dominant and recessive transmission patterns. Clinical presentations ranged from severe defects of limb morphogenesis to mild redundant ossification. We suggest that such clinical gradient can be linked to a continuum of GDF5-activity variation, with loss of GDF5 activity underlying bone development defects, and gain of function causing disorders with excessive bone formation. CONCLUSIONS: Our analysis of GDF5 pathogenicity mechanisms furtherly supports that mutation and zygosity backgrounds resulting in the same level of GDF5 activity may lead to similar phenotypes. This information can aid in interpreting the potential pathogenic effect of new variants and in supporting an appropriate genetic counseling.

Mechanisms of pathogenesis of missense mutations on the KDM6A-H3 interaction in type 2 Kabuki Syndrome.

Mutations in genes encoding for histone methylation proteins are associated with several developmental disorders. Among them, KDM6A is the disease causative gene of type 2 Kabuki Syndrome, a rare multisystem disease. While nonsense mutations and short insertions/deletions are known to trigger pathogenic mechanisms, the functional effects of missense mutations are still uncharacterized. In this study, we demonstrate that a selected set of missense mutations significantly hamper the interaction between KDM6A and the histone H3, by modifying the dynamics of the linker domain, and then causing a loss of function effect.

Pyntacle: a parallel computing-enabled framework for large-scale network biology analysis.

BACKGROUND: Some natural systems are big in size, complex, and often characterized by convoluted mechanisms of interaction, such as epistasis, pleiotropy, and trophism, which cannot be immediately ascribed to individual natural events or biological entities but that are often derived from group effects. However, the determination of important groups of entities, such as genes or proteins, in complex systems is considered a computationally hard task. RESULTS: We present Pyntacle, a high-performance framework designed to exploit parallel computing and graph theory to efficiently identify critical groups in big networks and in scenarios that cannot be tackled with traditional network analysis approaches. CONCLUSIONS: We showcase potential applications of Pyntacle with transcriptomics and structural biology data, thereby highlighting the outstanding improvement in terms of computational resources over existing tools.

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>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.