Hyperactive HRAS dysregulates energetic metabolism in fibroblasts from patients with Costello syndrome via enhanced production of reactive oxidizing species.

Germline-activating mutations in HRAS cause Costello syndrome (CS), a cancer prone multisystem disorder characterized by reduced postnatal growth. In CS, poor weight gain and growth are not caused by low caloric intake. Here, we show that constitutive plasma membrane translocation and activation of the GLUT4 glucose transporter, via reactive oxygen species-dependent AMP-activated protein kinase α and p38 hyperactivation, occurs in primary fibroblasts of CS patients, resulting in accelerated glycolysis and increased fatty acid synthesis and storage as lipid droplets. An accelerated autophagic flux was also identified as contributing to the increased energetic expenditure in CS. Concomitant inhibition of p38 and PI3K signaling by wortmannin was able to rescue both the dysregulated glucose intake and accelerated autophagic flux. Our findings provide a mechanistic link between upregulated HRAS function, defective growth and increased resting energetic expenditure in CS, and document that targeting p38 and PI3K signaling is able to revert this metabolic dysfunction.

Expanding the molecular spectrum of pathogenic SHOC2 variants underlying Mazzanti syndrome.

We previously molecularly and clinically characterized Mazzanti syndrome, a RASopathy related to Noonan syndrome that is mostly caused by a single recurrent missense variant (c.4A > G, p.Ser2Gly) in SHOC2, which encodes a leucine-rich repeat-containing protein facilitating signal flow through the RAS-mitogen-associated protein kinase (MAPK) pathway. We also documented that the pathogenic p.Ser2Gly substitution causes upregulation of MAPK signaling and constitutive targeting of SHOC2 to the plasma membrane due to the introduction of an N-myristoylation recognition motif. The almost invariant occurrence of the pathogenic c.4A > G missense change in SHOC2 is mirrored by a relatively homogeneous clinical phenotype of Mazzanti syndrome. Here, we provide new data on the clinical spectrum and molecular diversity of this disorder and functionally characterize new pathogenic variants. The clinical phenotype of six unrelated individuals carrying novel disease-causing SHOC2 variants is delineated, and public and newly collected clinical data are utilized to profile the disorder. In silico, in vitro and in vivo characterization of the newly identified variants provides evidence that the consequences of these missense changes on SHOC2 functional behavior differ from what had been observed for the canonical p.Ser2Gly change but converge toward an enhanced activation of the RAS-MAPK pathway. Our findings expand the molecular spectrum of pathogenic SHOC2 variants, provide a more accurate picture of the phenotypic expression associated with variants in this gene and definitively establish a gain-of-function behavior as the mechanism of disease.

SPRED2 loss-of-function causes a recessive Noonan syndrome-like phenotype.

Upregulated signal flow through RAS and the mitogen-associated protein kinase (MAPK) cascade is the unifying mechanistic theme of the RASopathies, a family of disorders affecting development and growth. Pathogenic variants in more than 20 genes have been causally linked to RASopathies, the majority having a dominant role in promoting enhanced signaling. Here, we report that SPRED2 loss of function is causally linked to a recessive phenotype evocative of Noonan syndrome. Homozygosity for three different variants-c.187C>T (p.Arg63∗), c.299T>C (p.Leu100Pro), and c.1142_1143delTT (p.Leu381Hisfs∗95)-were identified in four subjects from three families. All variants severely affected protein stability, causing accelerated degradation, and variably perturbed SPRED2 functional behavior. When overexpressed in cells, all variants were unable to negatively modulate EGF-promoted RAF1, MEK, and ERK phosphorylation, and time-course experiments in primary fibroblasts (p.Leu100Pro and p.Leu381Hisfs∗95) documented an increased and prolonged activation of the MAPK cascade in response to EGF stimulation. Morpholino-mediated knockdown of spred2a and spred2b in zebrafish induced defects in convergence and extension cell movements indicating upregulated RAS-MAPK signaling, which were rescued by expressing wild-type SPRED2 but not the SPRED2Leu381Hisfs∗95 protein. The clinical phenotype of the four affected individuals included developmental delay, intellectual disability, cardiac defects, short stature, skeletal anomalies, and a typical facial gestalt as major features, without the occurrence of the distinctive skin signs characterizing Legius syndrome. These features, in part, characterize the phenotype of Spred2-/- mice. Our findings identify the second recessive form of Noonan syndrome and document pleiotropic consequences of SPRED2 loss of function in development.

HECT-Type E3 Ubiquitin Ligases in Cancer.

Ubiquitination, a post-translational modification that involves a covalent attachment of ubiquitin to a protein substrate, is essential for cellular homeostatic maintenance. At the end of a three-enzyme cascade, E3 ubiquitin ligases (E3s) recruit substrates and promote or directly catalyze ubiquitin transfer to targets. These enzymes largely determine the specificity of the ubiquitination reaction. Genetic alteration, abnormal expression, or dysfunction of E3s account for the occurrence and progression of human cancers. Indeed, excessive degradation of relevant tumor-suppressor molecules and impaired disposal of oncogenic proteins have been linked to tumorigenesis. This review focuses on the emerging roles of HECT-type E3s in tumorigenesis, and emphasizes how perturbations of these enzymes contribute to cancer pathogenesis.

REDIportal: millions of novel A-to-I RNA editing events from thousands of RNAseq experiments.

RNA editing is a relevant epitranscriptome phenomenon able to increase the transcriptome and proteome diversity of eukaryotic organisms. ADAR mediated RNA editing is widespread in humans in which millions of A-to-I changes modify thousands of primary transcripts. RNA editing has pivotal roles in the regulation of gene expression or modulation of the innate immune response or functioning of several neurotransmitter receptors. Massive transcriptome sequencing has fostered the research in this field. Nonetheless, different aspects of the RNA editing biology are still unknown and need to be elucidated. To support the study of A-to-I RNA editing we have updated our REDIportal catalogue raising its content to about 16 millions of events detected in 9642 human RNAseq samples from the GTEx project by using a dedicated pipeline based on the HPC version of the REDItools software. REDIportal now allows searches at sample level, provides overviews of RNA editing profiles per each RNAseq experiment, implements a Gene View module to look at individual events in their genic context and hosts the CLAIRE database. Starting from this novel version, REDIportal will start collecting non-human RNA editing changes for comparative genomics investigations. The database is freely available at http://srv00.recas.ba.infn.it/atlas/index.html.

Posterior fossa ependymoma in neurodevelopmental syndrome caused by a de novo germline pathogenic POLR2A variant.

Ependymoma is the third most common pediatric brain tumor. Predisposition to develop ependymomas has been reported in different hereditary diseases, but the pathogenic variants related to the familial syndromes have rarely been detected in sporadic ependymomas. De novo variants in POLR2A, the gene encoding the largest subunit of RNA polymerase II, cause a neurodevelopmental disorder with a wide range of clinical manifestations, characterized by severe infantile-onset hypotonia, developmental delay, feeding difficulties, palatal anomalies, and facial dysmorphisms. As somatic events, POLR2A mutations represent a recurrent somatic lesion in benign meningiomas. Here we describe a case of ependymoma in a 2-year-old male with a de novo pathogenic variant in POLR2A predicted to impair proper interaction of the subunit with transcription-elongation factor TFIIS, whose function is required for back-tracking of the enzyme due to elongation blocks or nucleotide misincorporation, and expected to result in an increased error and reduced elongation rates. To date, ependymoma has never been reported in patients harboring pathogenic POLR2A variants. Further information is required to explore the possibility of a differential clinical and functional impact of the pathogenic POLR2A variants and the eventual inclusion of the POLR2A neurodevelopmental disorder among the cancer predisposition syndromes with the possible development of ependymomas.

The Spike Mutants Website: A Worldwide Used Resource against SARS-CoV-2.

A large number of SARS-CoV-2 mutations in a short period of time has driven scientific research related to vaccines, new drugs, and antibodies to combat the new variants of the virus. Herein, we present a web portal containing the structural information, the tridimensional coordinates, and the molecular dynamics trajectories of the SARS-CoV-2 spike protein and its main variants. The Spike Mutants website can serve as a rapid online tool for investigating the impact of novel mutations on virus fitness. Taking into account the high variability of SARS-CoV-2, this application can help the scientific community when prioritizing molecules for experimental assays, thus, accelerating the identification of promising drug candidates for COVID-19 treatment. Below we describe the main features of the platform and illustrate the possible applications for speeding up the drug discovery process and hypothesize new effective strategies to overcome the recurrent mutations in SARS-CoV-2 genome.

Structural Evolution of Delta (B.1.617.2) and Omicron (BA.1) Spike Glycoproteins.

The vast amount of epidemiologic and genomic data that were gathered as a global response to the COVID-19 pandemic that was caused by SARS-CoV-2 offer a unique opportunity to shed light on the structural evolution of coronaviruses and in particular on the spike (S) glycoprotein, which mediates virus entry into the host cell by binding to the human ACE2 receptor. Herein, we carry out an investigation into the dynamic properties of the S glycoprotein, focusing on the much more transmissible Delta and Omicron variants. Notwithstanding the great number of mutations that have accumulated, particularly in the Omicron S glycoprotein, our data clearly showed the conservation of some structural and dynamic elements, such as the global motion of the receptor binding domain (RBD). However, our studies also revealed structural and dynamic alterations that were concentrated in the aa 627-635 region, on a small region of the receptor binding motif (aa 483-485), and the so-called "fusion-peptide proximal region". In particular, these last two S regions are known to be involved in the human receptor ACE2 recognition and membrane fusion. Our structural evidence, therefore, is likely involved in the observed different transmissibility of these S mutants. Finally, we highlighted the role of glycans in the increased RBD flexibility of the monomer in the up conformation of Omicron.

Broadening the phenotypic spectrum of Beta3GalT6-associated phenotypes.

Biallelic mutations in B3GALT6, coding for a galactosyltransferase involved in the synthesis of glycosaminoglycans (GAGs), have been associated with various clinical conditions, causing spondyloepimetaphyseal dysplasia with joint laxity type 1 (SEMDJL1 or SEMDJL Beighton type), Al-Gazali syndrome (ALGAZ), and a severe progeroid form of Ehlers-Danlos syndrome (EDSSPD2). In the 2017 Ehlers-Danlos syndrome (EDS) classification, Beta3GalT6-related disorders were grouped in the spondylodysplastic EDSs together with spondylodysplastic EDSs due to B4GALT7 and SLC39A13 mutations. Herein, we describe a patient with a previously unreported homozygous pathogenic B3GALT6 variant resulting in a complex phenotype more severe than spondyloepimetaphyseal dysplasia with joint laxity type 1, and having dural ectasia and aortic dilation as additionally associated features, further broadening the phenotypic spectrum of the Beta3GalT6-related syndromes. We also document the utility of repeating sequencing in patients with uninformative exomes, particularly when performed by using "first generations" enrichment capture methods.

Altered Local Interactions and Long-Range Communications in UK Variant (B.1.1.7) Spike Glycoprotein.

The COVID-19 pandemic is caused by SARS-CoV-2. Currently, most of the research efforts towards the development of vaccines and antibodies against SARS-CoV-2 were mainly focused on the spike (S) protein, which mediates virus entry into the host cell by binding to ACE2. As the virus SARS-CoV-2 continues to spread globally, variants have emerged, characterized by multiple mutations of the S glycoprotein. Herein, we employed microsecond-long molecular dynamics simulations to study the impact of the mutations of the S glycoprotein in SARS-CoV-2 Variant of Concern 202012/01 (B.1.1.7), termed the "UK variant", in comparison with the wild type, with the aim to decipher the structural basis of the reported increased infectivity and virulence. The simulations provided insights on the different dynamics of UK and wild-type S glycoprotein, regarding in particular the Receptor Binding Domain (RBD). In addition, we investigated the role of glycans in modulating the conformational transitions of the RBD. The overall results showed that the UK mutant experiences higher flexibility in the RBD with respect to wild type; this behavior might be correlated with the increased transmission reported for this variant. Our work also adds useful structural information on antigenic "hotspots" and epitopes targeted by neutralizing antibodies.

Transcriptomic Characterization of Cow, Donkey and Goat Milk Extracellular Vesicles Reveals Their Anti-Inflammatory and Immunomodulatory Potential.

Milk extracellular vesicles (mEVs) seem to be one of the main maternal messages delivery systems. Extracellular vesicles (EVs) are micro/nano-sized membrane-bound structures enclosing signaling molecules and thus acting as signal mediators between distant cells and/or tissues, exerting biological effects such as immune modulation and pro-regenerative activity. Milk is also a unique, scalable, and reliable source of EVs. Our aim was to characterize the RNA content of cow, donkey, and goat mEVs through transcriptomic analysis of mRNA and small RNA libraries. Over 10,000 transcripts and 2000 small RNAs were expressed in mEVs of each species. Among the most represented transcripts, 110 mRNAs were common between the species with cow acting as the most divergent. The most represented small RNA class was miRNA in all the species, with 10 shared miRNAs having high impact on the immune regulatory function. Functional analysis for the most abundant mRNAs shows epigenetic functions such as histone modification, telomere maintenance, and chromatin remodeling for cow; lipid catabolism, oxidative stress, and vitamin metabolism for donkey; and terms related to chemokine receptor interaction, leukocytes migration, and transcriptional regulation in response to stress for goat. For miRNA targets, shared terms emerged as the main functions for all the species: immunity modulation, protein synthesis, cellular cycle regulation, transmembrane exchanges, and ion channels. Moreover, donkey and goat showed additional terms related to epigenetic modification and DNA maintenance. Our results showed a potential mEVs immune regulatory purpose through their RNA cargo, although in vivo validation studies are necessary.

HPC-REDItools: a novel HPC-aware tool for improved large scale RNA-editing analysis.

BACKGROUND: RNA editing is a widespread co-/post-transcriptional mechanism that alters primary RNA sequences through the modification of specific nucleotides and it can increase both the transcriptome and proteome diversity. The automatic detection of RNA-editing from RNA-seq data is computational intensive and limited to small data sets, thus preventing a reliable genome-wide characterisation of such process. RESULTS: In this work we introduce HPC-REDItools, an upgraded tool for accurate RNA-editing events discovery from large dataset repositories. AVAILABILITY: https://github.com/BioinfoUNIBA/REDItools2 . CONCLUSIONS: HPC-REDItools is dramatically faster than the previous version, REDItools, enabling big-data analysis by means of a MPI-based implementation and scaling almost linearly with the number of available cores.

ELIXIR-IT HPC@CINECA: high performance computing resources for the bioinformatics community.

BACKGROUND: The advent of Next Generation Sequencing (NGS) technologies and the concomitant reduction in sequencing costs allows unprecedented high throughput profiling of biological systems in a cost-efficient manner. Modern biological experiments are increasingly becoming both data and computationally intensive and the wealth of publicly available biological data is introducing bioinformatics into the "Big Data" era. For these reasons, the effective application of High Performance Computing (HPC) architectures is becoming progressively more recognized also by bioinformaticians. Here we describe HPC resources provisioning pilot programs dedicated to bioinformaticians, run by the Italian Node of ELIXIR (ELIXIR-IT) in collaboration with CINECA, the main Italian supercomputing center. RESULTS: Starting from April 2016, CINECA and ELIXIR-IT launched the pilot Call "ELIXIR-IT HPC@CINECA", offering streamlined access to HPC resources for bioinformatics. Resources are made available either through web front-ends to dedicated workflows developed at CINECA or by providing direct access to the High Performance Computing systems through a standard command-line interface tailored for bioinformatics data analysis. This allows to offer to the biomedical research community a production scale environment, continuously updated with the latest available versions of publicly available reference datasets and bioinformatic tools. Currently, 63 research projects have gained access to the HPC@CINECA program, for a total handout of ~ 8 Millions of CPU/hours and, for data storage, ~ 100 TB of permanent and ~ 300 TB of temporary space. CONCLUSIONS: Three years after the beginning of the ELIXIR-IT HPC@CINECA program, we can appreciate its impact over the Italian bioinformatics community and draw some considerations. Several Italian researchers who applied to the program have gained access to one of the top-ranking public scientific supercomputing facilities in Europe. Those investigators had the opportunity to sensibly reduce computational turnaround times in their research projects and to process massive amounts of data, pursuing research approaches that would have been otherwise difficult or impossible to undertake. Moreover, by taking advantage of the wealth of documentation and training material provided by CINECA, participants had the opportunity to improve their skills in the usage of HPC systems and be better positioned to apply to similar EU programs of greater scale, such as PRACE. To illustrate the effective usage and impact of the resources awarded by the program - in different research applications - we report five successful use cases, which have already published their findings in peer-reviewed journals.

Molecular dynamics recipes for genome research.

Molecular dynamics (MD) simulation allows one to predict the time evolution of a system of interacting particles. It is widely used in physics, chemistry and biology to address specific questions about the structural properties and dynamical mechanisms of model systems. MD earned a great success in genome research, as it proved to be beneficial in sorting pathogenic from neutral genomic mutations. Considering their computational requirements, simulations are commonly performed on HPC computing devices, which are generally expensive and hard to administer. However, variables like the software tool used for modeling and simulation or the size of the molecule under investigation might make one hardware type or configuration more advantageous than another or even make the commodity hardware definitely suitable for MD studies. This work aims to shed lights on this aspect.

New Insights into the Effect of Residue Mutations on the Rotavirus VP1 Function Using Molecular Dynamic Simulations.

Rotavirus group A remains a major cause of diarrhea in infants and young children worldwide. The permanent emergence of new genotypes puts the potential effectiveness of vaccines under serious questions. Thirteen VP1 structures with mutations mapping to the RNA entry site were analyzed using molecular dynamics simulations, and the results were combined with the experimental findings reported previously. The results revealed structural fluctuations in the protein-protein recognition sites and in the bottleneck of the RNA entry site that may affect the interaction of different proteins and delay the initiation of the viral replication, respectively. Altogether, the structural analysis of VP1 in the region crucial for the initiation of the viral replication, mainly the bottleneck site, may boost efforts to develop antivirals, as they might complement the available vaccines.

Multiple Recombination Events and Strong Purifying Selection at the Origin of SARS-CoV-2 Spike Glycoprotein Increased Correlated Dynamic Movements.

Our evolutionary and structural analyses revealed that the severe acute respiratory syndrome (SARS) coronavirus 2 (SARS-CoV-2) spike gene is a complex mosaic resulting from several recombination events. Additionally, the fixation of variants has mainly been driven by purifying selection, suggesting the presence of conserved structural features. Our dynamic simulations identified two main long-range covariant dynamic movements of the novel glycoprotein, and showed that, as a result of the evolutionary duality, they are preserved. The first movement involves the receptor binding domain with the N-terminal domain and the C-terminal domain 2 and is maintained across human, bat and pangolin coronaviruses. The second is a complex network of long-range dynamics specific to SARS-CoV-2 involving the novel PRRA and the conserved KR*SF cleavage sites, as well as conserved segments in C-terminal domain 3. These movements, essential for host cell binding, are maintained by hinges conserved across human, bat, and pangolin coronaviruses glycoproteins. The hinges, located around Threonine 333 and Proline 527 within the N-terminal domain and C-terminal domain 2, represent candidate targets for the future development of novel pan-coronavirus inhibitors. In summary, we show that while recombination created a new configuration that increased the covariant dynamic movements of the SARS-CoV-2 glycoprotein, negative selection preserved its inter-domain structure throughout evolution in different hosts and inter-species transmissions.

DROPA: DRIP-seq optimized peak annotator.

BACKGROUND: R-loops are three-stranded nucleic acid structures that usually form during transcription and that may lead to gene regulation or genome instability. DRIP (DNA:RNA Immunoprecipitation)-seq techniques are widely used to map R-loops genome-wide providing insights into R-loop biology. However, annotation of DRIP-seq peaks to genes can be a tricky step, due to the lack of strand information when using the common basic DRIP technique. RESULTS: Here, we introduce DRIP-seq Optimized Peak Annotator (DROPA), a new tool for gene annotation of R-loop peaks based on gene expression information. DROPA allows a full customization of annotation options, ranging from the choice of reference datasets to gene feature definitions. DROPA allows to assign R-loop peaks to the DNA template strand in gene body with a false positive rate of less than 7%. A comparison of DROPA performance with three widely used annotation tools show that it identifies less false positive annotations than the others. CONCLUSIONS: DROPA is a fully customizable peak-annotation tool optimized for co-transcriptional DRIP-seq peaks, which allows a finest gene annotation based on gene expression information. Its output can easily be integrated into pipelines to perform downstream analyses, while useful and informative summary plots and statistical enrichment tests can be produced.

Clinical and functional characterization of a novel RASopathy-causing SHOC2 mutation associated with prenatal-onset hypertrophic cardiomyopathy.

SHOC2 is a scaffold protein mediating RAS-promoted activation of mitogen-activated protein kinase (MAPK) signaling in response to extracellular stimuli. A recurrent activating mutation in SHOC2 (p.Ser2Gly) causes Mazzanti syndrome, a RASopathy characterized by features resembling Noonan syndrome and distinctive ectodermal abnormalities. A second mutation (p.Met173Ile) supposed to cause loss-of-function was more recently identified in two individuals with milder phenotypes. Here, we report on the third RASopathy-causing SHOC2 mutation (c.807_808delinsTT, p.Gln269_His270delinsHisTyr), which was found associated with prenatal-onset hypertrophic cardiomyopathy. Structural analyses indicated a possible impact of the mutation on the relative orientation of the two SHOC2's leucine-rich repeat domains. Functional studies provided evidence of its activating role, revealing enhanced binding of the mutant protein to MRAS and PPP1CB, and increased signaling through the MAPK cascade. Differing from SHOC2 S2G , SHOC2 Q269_H270delinsHY is not constitutively targeted to the plasma membrane. These data document that diverse mechanisms in SHOC2 functional dysregulation converge toward MAPK signaling upregulation.

Biallelic Mutations in TBCD, Encoding the Tubulin Folding Cofactor D, Perturb Microtubule Dynamics and Cause Early-Onset Encephalopathy.

Microtubules are dynamic cytoskeletal elements coordinating and supporting a variety of neuronal processes, including cell division, migration, polarity, intracellular trafficking, and signal transduction. Mutations in genes encoding tubulins and microtubule-associated proteins are known to cause neurodevelopmental and neurodegenerative disorders. Growing evidence suggests that altered microtubule dynamics may also underlie or contribute to neurodevelopmental disorders and neurodegeneration. We report that biallelic mutations in TBCD, encoding one of the five co-chaperones required for assembly and disassembly of the αß-tubulin heterodimer, the structural unit of microtubules, cause a disease with neurodevelopmental and neurodegenerative features characterized by early-onset cortical atrophy, secondary hypomyelination, microcephaly, thin corpus callosum, developmental delay, intellectual disability, seizures, optic atrophy, and spastic quadriplegia. Molecular dynamics simulations predicted long-range and/or local structural perturbations associated with the disease-causing mutations. Biochemical analyses documented variably reduced levels of TBCD, indicating relative instability of mutant proteins, and defective ß-tubulin binding in a subset of the tested mutants. Reduced or defective TBCD function resulted in decreased soluble α/ß-tubulin levels and accelerated microtubule polymerization in fibroblasts from affected subjects, demonstrating an overall shift toward a more rapidly growing and stable microtubule population. These cells displayed an aberrant mitotic spindle with disorganized, tangle-shaped microtubules and reduced aster formation, which however did not alter appreciably the rate of cell proliferation. Our findings establish that defective TBCD function underlies a recognizable encephalopathy and drives accelerated microtubule polymerization and enhanced microtubule stability, underscoring an additional cause of altered microtubule dynamics with impact on neuronal function and survival in the developing brain.

CoVaCS: a consensus variant calling system.

BACKGROUND: The advent and ongoing development of next generation sequencing technologies (NGS) has led to a rapid increase in the rate of human genome re-sequencing data, paving the way for personalized genomics and precision medicine. The body of genome resequencing data is progressively increasing underlining the need for accurate and time-effective bioinformatics systems for genotyping - a crucial prerequisite for identification of candidate causal mutations in diagnostic screens. RESULTS: Here we present CoVaCS, a fully automated, highly accurate system with a web based graphical interface for genotyping and variant annotation. Extensive tests on a gold standard benchmark data-set -the NA12878 Illumina platinum genome- confirm that call-sets based on our consensus strategy are completely in line with those attained by similar command line based approaches, and far more accurate than call-sets from any individual tool. Importantly our system exhibits better sensitivity and higher specificity than equivalent commercial software. CONCLUSIONS: CoVaCS offers optimized pipelines integrating state of the art tools for variant calling and annotation for whole genome sequencing (WGS), whole-exome sequencing (WES) and target-gene sequencing (TGS) data. The system is currently hosted at Cineca, and offers the speed of a HPC computing facility, a crucial consideration when large numbers of samples must be analysed. Importantly, all the analyses are performed automatically allowing high reproducibility of the results. As such, we believe that CoVaCS can be a valuable tool for the analysis of human genome resequencing studies. CoVaCS is available at: https://bioinformatics.cineca.it/covacs .