Durable and efficient gene silencing in vivo by hit-and-run epigenome editing.

Permanent epigenetic silencing using programmable editors equipped with transcriptional repressors holds great promise for the treatment of human diseases1-3. However, to unlock its full therapeutic potential, an experimental confirmation of durable epigenetic silencing after the delivery of transient delivery of editors in vivo is needed. To this end, here we targeted Pcsk9, a gene expressed in hepatocytes that is involved in cholesterol homeostasis. In vitro screening of different editor designs indicated that zinc-finger proteins were the best-performing DNA-binding platform for efficient silencing of mouse Pcsk9. A single administration of lipid nanoparticles loaded with the editors' mRNAs almost halved the circulating levels of PCSK9 for nearly one year in mice. Notably, Pcsk9 silencing and accompanying epigenetic repressive marks also persisted after forced liver regeneration, further corroborating the heritability of the newly installed epigenetic state. Improvements in construct design resulted in the development of an all-in-one configuration that we term evolved engineered transcriptional repressor (EvoETR). This design, which is characterized by a high specificity profile, further reduced the circulating levels of PCSK9 in mice with an efficiency comparable with that obtained through conventional gene editing, but without causing DNA breaks. Our study lays the foundation for the development of in vivo therapeutics that are based on epigenetic silencing.

Interferon-inducible phospholipids govern IFITM3-dependent endosomal antiviral immunity.

The interferon-induced transmembrane proteins (IFITM) are implicated in several biological processes, including antiviral defense, but their modes of action remain debated. Here, taking advantage of pseudotyped viral entry assays and replicating viruses, we uncover the requirement of host co-factors for endosomal antiviral inhibition through high-throughput proteomics and lipidomics in cellular models of IFITM restriction. Unlike plasma membrane (PM)-localized IFITM restriction that targets infectious SARS-CoV2 and other PM-fusing viral envelopes, inhibition of endosomal viral entry depends on lysines within the conserved IFITM intracellular loop. These residues recruit Phosphatidylinositol 3,4,5-trisphosphate (PIP3) that we show here to be required for endosomal IFITM activity. We identify PIP3 as an interferon-inducible phospholipid that acts as a rheostat for endosomal antiviral immunity. PIP3 levels correlated with the potency of endosomal IFITM restriction and exogenous PIP3 enhanced inhibition of endocytic viruses, including the recent SARS-CoV2 Omicron variant. Together, our results identify PIP3 as a critical regulator of endosomal IFITM restriction linking it to the Pi3K/Akt/mTORC pathway and elucidate cell-compartment-specific antiviral mechanisms with potential relevance for the development of broadly acting antiviral strategies.

Circulating hematopoietic stem/progenitor cell subsets contribute to human hematopoietic homeostasis.

ABSTRACT: In physiological conditions, few circulating hematopoietic stem/progenitor cells (cHSPCs) are present in the peripheral blood, but their contribution to human hematopoiesis remain unsolved. By integrating advanced immunophenotyping, single-cell transcriptional and functional profiling, and integration site (IS) clonal tracking, we unveiled the biological properties and the transcriptional features of human cHSPC subpopulations in relationship to their bone marrow (BM) counterpart. We found that cHSPCs reduced in cell count over aging and are enriched for primitive, lymphoid, and erythroid subpopulations, showing preactivated transcriptional and functional state. Moreover, cHSPCs have low expression of multiple BM-retention molecules but maintain their homing potential after xenotransplantation. By generating a comprehensive human organ-resident HSPC data set based on single-cell RNA sequencing data, we detected organ-specific seeding properties of the distinct trafficking HSPC subpopulations. Notably, circulating multi-lymphoid progenitors are primed for seeding the thymus and actively contribute to T-cell production. Human clonal tracking data from patients receiving gene therapy (GT) also showed that cHSPCs connect distant BM niches and participate in steady-state hematopoietic production, with primitive cHSPCs having the highest recirculation capability to travel in and out of the BM. Finally, in case of hematopoietic impairment, cHSPCs composition reflects the BM-HSPC content and might represent a biomarker of the BM state for clinical and research purposes. Overall, our comprehensive work unveiled fundamental insights into the in vivo dynamics of human HSPC trafficking and its role in sustaining hematopoietic homeostasis. GT patients' clinical trials were registered at ClinicalTrials.gov (NCT01515462 and NCT03837483) and EudraCT (2009-017346-32 and 2018-003842-18).

Lipid nanoparticles allow efficient and harmless ex vivo gene editing of human hematopoietic cells.

Ex vivo gene editing in T cells and hematopoietic stem/progenitor cells (HSPCs) holds promise for treating diseases. Gene editing encompasses the delivery of a programmable editor RNA or ribonucleoprotein, often achieved ex vivo via electroporation, and when aiming for homology-driven correction of a DNA template, often provided by viral vectors together with a nuclease editor. Although HSPCs activate a robust p53-dependent DNA damage response upon nuclease-based editing, the responses triggered in T cells remain poorly characterized. Here, we performed comprehensive multiomics analyses and found that electroporation is the main culprit of cytotoxicity in T cells, causing death and cell cycle delay, perturbing metabolism, and inducing an inflammatory response. Nuclease RNA delivery using lipid nanoparticles (LNPs) nearly abolished cell death and ameliorated cell growth, improving tolerance to the procedure and yielding a higher number of edited cells compared with using electroporation. Transient transcriptomic changes upon LNP treatment were mostly caused by cellular loading with exogenous cholesterol, whose potentially detrimental impact could be overcome by limiting exposure. Notably, LNP-based HSPC editing dampened p53 pathway induction and supported higher clonogenic activity and similar or higher reconstitution by long-term repopulating HSPCs compared with electroporation, reaching comparable editing efficiencies. Overall, LNPs may allow efficient and harmless ex vivo gene editing in hematopoietic cells for the treatment of human diseases.

Removal of innate immune barriers allows efficient transduction of quiescent human hematopoietic stem cells.

Quiescent human hematopoietic stem cells (HSC) are ideal targets for gene therapy applications due to their preserved stemness and repopulation capacities; however, they have not been exploited extensively because of their resistance to genetic manipulation. We report here the development of a lentiviral transduction protocol that overcomes this resistance in long-term repopulating quiescent HSC, allowing their efficient genetic manipulation. Mechanistically, lentiviral vector transduction of quiescent HSC was found to be restricted at the level of vector entry and by limited pyrimidine pools. These restrictions were overcome by the combined addition of cyclosporin H (CsH) and deoxynucleosides (dNs) during lentiviral vector transduction. Clinically relevant transduction levels were paired with higher polyclonal engraftment of long-term repopulating HSC as compared with standard ex vivo cultured controls. These findings identify the cell-intrinsic barriers that restrict the transduction of quiescent HSC and provide a means to overcome them, paving the way for the genetic engineering of unstimulated HSC.

Mesenchymal stromal cells improve the transplantation outcome of CRISPR-Cas9 gene-edited human HSPCs.

Mesenchymal stromal cells (MSCs) have been employed in vitro to support hematopoietic stem and progenitor cell (HSPC) expansion and in vivo to promote HSPC engraftment. Based on these studies, we developed an MSC-based co-culture system to optimize the transplantation outcome of clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 gene-edited (GE) human HSPCs. We show that bone marrow (BM)-MSCs produce several hematopoietic supportive and anti-inflammatory factors capable of alleviating the proliferation arrest and mitigating the apoptotic and inflammatory programs activated in GE-HSPCs, improving their expansion and clonogenic potential in vitro. The use of BM-MSCs resulted in superior human engraftment and increased clonal output of GE-HSPCs contributing to the early phase of hematological reconstitution in the peripheral blood of transplanted mice. In conclusion, our work poses the biological bases for a novel clinical use of BM-MSCs to promote engraftment of GE-HSPCs and improve their transplantation outcome.

Follicular helper T cell signature of replicative exhaustion, apoptosis, and senescence in common variable immunodeficiency.

Common variable immunodeficiency (CVID) is the most frequent primary antibody deficiency whereby follicular helper T (Tfh) cells fail to establish productive responses with B cells in germinal centers. Here, we analyzed the frequency, phenotype, transcriptome, and function of circulating Tfh (cTfh) cells in CVID patients displaying autoimmunity as an additional phenotype. A group of patients showed a high frequency of cTfh1 cells and a prominent expression of PD-1 and ICOS as well as a cTfh mRNA signature consistent with highly activated, but exhausted, senescent, and apoptotic cells. Plasmatic CXCL13 levels were elevated in this group and positively correlated with cTfh1 cell frequency and PD-1 levels. Monoallelic variants in RTEL1, a telomere length- and DNA repair-related gene, were identified in four patients belonging to this group. Their blood lymphocytes showed shortened telomeres, while their cTfh were more prone to apoptosis. These data point toward a novel pathogenetic mechanism in CVID, whereby alterations in DNA repair and telomere elongation might predispose to antibody deficiency. A Th1, highly activated but exhausted and apoptotic cTfh phenotype was associated with this form of CVID.

Myelomonocytic cells in giant cell arteritis activate trained immunity programs sustaining inflammation and cytokine production.

OBJECTIVE: Trained immunity (TI) is a de facto memory program of innate immune cells, characterized by immunometabolic and epigenetic changes sustaining enhanced production of cytokines. TI evolved as a protective mechanism against infections; however, inappropriate activation can cause detrimental inflammation and might be implicated in the pathogenesis of chronic inflammatory diseases. In this study, we investigated the role of TI in the pathogenesis of giant cell arteritis (GCA), a large-vessel vasculitis characterized by aberrant macrophage activation and excess cytokine production. METHODS: Monocytes from GCA patients and from age- and sex-matched healthy donors were subjected to polyfunctional studies, including cytokine production assays at baseline and following stimulation, intracellular metabolomics, chromatin immunoprecipitation-qPCR, and combined ATAC/RNA sequencing. Immunometabolic activation (i.e. glycolysis) was assessed in inflamed vessels of GCA patients with FDG-PET and immunohistochemistry (IHC), and the role of this pathway in sustaining cytokine production was confirmed with selective pharmacologic inhibition in GCA monocytes. RESULTS: GCA monocytes exhibited hallmark molecular features of TI. Specifically, these included enhanced IL-6 production upon stimulation, typical immunometabolic changes (e.g. increased glycolysis and glutaminolysis) and epigenetic changes promoting enhanced transcription of genes governing pro-inflammatory activation. Immunometabolic changes of TI (i.e. glycolysis) were a feature of myelomonocytic cells in GCA lesions and were required for enhanced cytokine production. CONCLUSIONS: Myelomonocytic cells in GCA activate TI programs sustaining enhanced inflammatory activation with excess cytokine production.

Oncogene-induced maladaptive activation of trained immunity in the pathogenesis and treatment of Erdheim-Chester disease.

Trained immunity (TI) is a proinflammatory program induced in monocyte/macrophages upon sensing of specific pathogens and is characterized by immunometabolic and epigenetic changes that enhance cytokine production. Maladaptive activation of TI (ie, in the absence of infection) may result in detrimental inflammation and development of disease; however, the exact role and extent of inappropriate activation of TI in the pathogenesis of human diseases is undetermined. In this study, we uncovered the oncogene-induced, maladaptive induction of TI in the pathogenesis of a human inflammatory myeloid neoplasm (Erdheim-Chester disease, [ECD]), characterized by the BRAFV600E oncogenic mutation in monocyte/macrophages and excess cytokine production. Mechanistically, myeloid cells expressing BRAFV600E exhibit all molecular features of TI: activation of the AKT/mammalian target of rapamycin signaling axis; increased glycolysis, glutaminolysis, and cholesterol synthesis; epigenetic changes on promoters of genes encoding cytokines; and enhanced cytokine production leading to hyperinflammatory responses. In patients with ECD, effective therapeutic strategies combat this maladaptive TI phenotype; in addition, pharmacologic inhibition of immunometabolic changes underlying TI (ie, glycolysis) effectively dampens cytokine production by myeloid cells. This study revealed the deleterious potential of inappropriate activation of TI in the pathogenesis of human inflammatory myeloid neoplasms and the opportunity for inhibition of TI in conditions characterized by maladaptive myeloid-driven inflammation.

Hematopoietic stem cell function in ß-thalassemia is impaired and is rescued by targeting the bone marrow niche.

Hematopoietic stem cells (HSCs) are regulated by signals from the bone marrow (BM) niche that tune hematopoiesis at steady state and in hematologic disorders. To understand HSC-niche interactions in altered nonmalignant homeostasis, we selected ß-thalassemia, a hemoglobin disorder, as a paradigm. In this severe congenital anemia, alterations secondary to the primary hemoglobin defect have a potential impact on HSC-niche cross talk. We report that HSCs in thalassemic mice (th3) have an impaired function, caused by the interaction with an altered BM niche. The HSC self-renewal defect is rescued after cell transplantation into a normal microenvironment, thus proving the active role of the BM stroma. Consistent with the common finding of osteoporosis in patients, we found reduced bone deposition with decreased levels of parathyroid hormone (PTH), which is a key regulator of bone metabolism but also of HSC activity. In vivo activation of PTH signaling through the reestablished Jagged1 and osteopontin levels correlated with the rescue of the functional pool of th3 HSCs by correcting HSC-niche cross talk. Reduced HSC quiescence was confirmed in thalassemic patients, along with altered features of the BM stromal niche. Our findings reveal a defect in HSCs in ß-thalassemia induced by an altered BM microenvironment and provide novel and relevant insight for improving transplantation and gene therapy approaches.

Advantages of using graph databases to explore chromatin conformation capture experiments.

BACKGROUND: High-throughput sequencing Chromosome Conformation Capture (Hi-C) allows the study of DNA interactions and 3D chromosome folding at the genome-wide scale. Usually, these data are represented as matrices describing the binary contacts among the different chromosome regions. On the other hand, a graph-based representation can be advantageous to describe the complex topology achieved by the DNA in the nucleus of eukaryotic cells. METHODS: Here we discuss the use of a graph database for storing and analysing data achieved by performing Hi-C experiments. The main issue is the size of the produced data and, working with a graph-based representation, the consequent necessity of adequately managing a large number of edges (contacts) connecting nodes (genes), which represents the sources of information. For this, currently available graph visualisation tools and libraries fall short with Hi-C data. The use of graph databases, instead, supports both the analysis and the visualisation of the spatial pattern present in Hi-C data, in particular for comparing different experiments or for re-mapping omics data in a space-aware context efficiently. In particular, the possibility of describing graphs through statistical indicators and, even more, the capability of correlating them through statistical distributions allows highlighting similarities and differences among different Hi-C experiments, in different cell conditions or different cell types. RESULTS: These concepts have been implemented in NeoHiC, an open-source and user-friendly web application for the progressive visualisation and analysis of Hi-C networks based on the use of the Neo4j graph database (version 3.5). CONCLUSION: With the accumulation of more experiments, the tool will provide invaluable support to compare neighbours of genes across experiments and conditions, helping in highlighting changes in functional domains and identifying new co-organised genomic compartments.

γ-TRIS: a graph-algorithm for comprehensive identification of vector genomic insertion sites.

SUMMARY: Retroviruses and their vector derivatives integrate semi-randomly in the genome of host cells and are inherited by their progeny as stable genetic marks. The retrieval and mapping of the sequences flanking the virus-host DNA junctions allows the identification of insertion sites in gene therapy or virally infected patients, essential for monitoring the evolution of genetically modified cells in vivo. However, since ∼30% of insertions land in low complexity or repetitive regions of the host cell genome, they cannot be correctly assigned and are currently discarded, limiting the accuracy and predictive power of clonal tracking studies. Here, we present γ-TRIS, a new graph-based genome-free alignment tool for identifying insertion sites even if embedded in low complexity regions. By using γ-TRIS to reanalyze clinical studies, we observed improvements in clonal quantification and tracking. AVAILABILITY AND IMPLEMENTATION: Source code at https://bitbucket.org/bereste/g-tris. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Expanded circulating hematopoietic stem/progenitor cells as novel cell source for the treatment of TCIRG1 osteopetrosis.

Allogeneic hematopoietic stem cell transplantation is the treatment of choice for autosomal recessive osteopetrosis caused by defects in the TCIRG1 gene. Despite recent progress in conditioning, a relevant number of patients are not eligible for allogeneic stem cell transplantation because of the severity of the disease and significant transplant-related morbidity. We exploited peripheral CD34+ cells, known to circulate at high frequency in the peripheral blood of TCIRG1-deficient patients, as a novel cell source for autologous transplantation of gene corrected cells. Detailed phenotypical analysis showed that circulating CD34+ cells have a cellular composition that resembles bone marrow, supporting their use in gene therapy protocols. Transcriptomic profile revealed enrichment in genes expressed by hematopoietic stem and progenitor cells (HSPCs). To overcome the limit of bone marrow harvest/ HSPC mobilization and serial blood drawings in TCIRG1 patients, we applied UM171-based ex-vivo expansion of HSPCs coupled with lentiviral gene transfer. Circulating CD34+ cells from TCIRG1-defective patients were transduced with a clinically-optimized lentiviral vector (LV) expressing TCIRG1 under the control of phosphoglycerate promoter and expanded ex vivo. Expanded cells maintained long-term engraftment capacity and multi-lineage repopulating potential when transplanted in vivo both in primary and secondary NSG recipients. Moreover, when CD34+ cells were differentiated in vitro, genetically corrected osteoclasts resorbed the bone efficiently. Overall, we provide evidence that expansion of circulating HSPCs coupled to gene therapy can overcome the limit of stem cell harvest in osteopetrotic patients, thus opening the way to future gene-based treatment of skeletal diseases caused by bone marrow fibrosis.

GenHap: a novel computational method based on genetic algorithms for haplotype assembly.

BACKGROUND: In order to fully characterize the genome of an individual, the reconstruction of the two distinct copies of each chromosome, called haplotypes, is essential. The computational problem of inferring the full haplotype of a cell starting from read sequencing data is known as haplotype assembly, and consists in assigning all heterozygous Single Nucleotide Polymorphisms (SNPs) to exactly one of the two chromosomes. Indeed, the knowledge of complete haplotypes is generally more informative than analyzing single SNPs and plays a fundamental role in many medical applications. RESULTS: To reconstruct the two haplotypes, we addressed the weighted Minimum Error Correction (wMEC) problem, which is a successful approach for haplotype assembly. This NP-hard problem consists in computing the two haplotypes that partition the sequencing reads into two disjoint sub-sets, with the least number of corrections to the SNP values. To this aim, we propose here GenHap, a novel computational method for haplotype assembly based on Genetic Algorithms, yielding optimal solutions by means of a global search process. In order to evaluate the effectiveness of our approach, we run GenHap on two synthetic (yet realistic) datasets, based on the Roche/454 and PacBio RS II sequencing technologies. We compared the performance of GenHap against HapCol, an efficient state-of-the-art algorithm for haplotype phasing. Our results show that GenHap always obtains high accuracy solutions (in terms of haplotype error rate), and is up to 4× faster than HapCol in the case of Roche/454 instances and up to 20× faster when compared on the PacBio RS II dataset. Finally, we assessed the performance of GenHap on two different real datasets. CONCLUSIONS: Future-generation sequencing technologies, producing longer reads with higher coverage, can highly benefit from GenHap, thanks to its capability of efficiently solving large instances of the haplotype assembly problem. Moreover, the optimization approach proposed in GenHap can be extended to the study of allele-specific genomic features, such as expression, methylation and chromatin conformation, by exploiting multi-objective optimization techniques. The source code and the full documentation are available at the following GitHub repository: https://github.com/andrea-tango/GenHap .

How computer science can help in understanding the 3D genome architecture.

Chromosome conformation capture techniques are producing a huge amount of data about the architecture of our genome. These data can provide us with a better understanding of the events that induce critical regulations of the cellular function from small changes in the three-dimensional genome architecture. Generating a unified view of spatial, temporal, genetic and epigenetic properties poses various challenges of data analysis, visualization, integration and mining, as well as of high performance computing and big data management. Here, we describe the critical issues of this new branch of bioinformatics, oriented at the comprehension of the three-dimensional genome architecture, which we call 'Nucleome Bioinformatics', looking beyond the currently available tools and methods, and highlight yet unaddressed challenges and the potential approaches that could be applied for tackling them. Our review provides a map for researchers interested in using computer science for studying 'Nucleome Bioinformatics', to achieve a better understanding of the biological processes that occur inside the nucleus.

VISPA2: a scalable pipeline for high-throughput identification and annotation of vector integration sites.

BACKGROUND: Bioinformatics tools designed to identify lentiviral or retroviral vector insertion sites in the genome of host cells are used to address the safety and long-term efficacy of hematopoietic stem cell gene therapy applications and to study the clonal dynamics of hematopoietic reconstitution. The increasing number of gene therapy clinical trials combined with the increasing amount of Next Generation Sequencing data, aimed at identifying integration sites, require both highly accurate and efficient computational software able to correctly process "big data" in a reasonable computational time. RESULTS: Here we present VISPA2 (Vector Integration Site Parallel Analysis, version 2), the latest optimized computational pipeline for integration site identification and analysis with the following features: (1) the sequence analysis for the integration site processing is fully compliant with paired-end reads and includes a sequence quality filter before and after the alignment on the target genome; (2) an heuristic algorithm to reduce false positive integration sites at nucleotide level to reduce the impact of Polymerase Chain Reaction or trimming/alignment artifacts; (3) a classification and annotation module for integration sites; (4) a user friendly web interface as researcher front-end to perform integration site analyses without computational skills; (5) the time speedup of all steps through parallelization (Hadoop free). CONCLUSIONS: We tested VISPA2 performances using simulated and real datasets of lentiviral vector integration sites, previously obtained from patients enrolled in a hematopoietic stem cell gene therapy clinical trial and compared the results with other preexisting tools for integration site analysis. On the computational side, VISPA2 showed a > 6-fold speedup and improved precision and recall metrics (1 and 0.97 respectively) compared to previously developed computational pipelines. These performances indicate that VISPA2 is a fast, reliable and user-friendly tool for integration site analysis, which allows gene therapy integration data to be handled in a cost and time effective fashion. Moreover, the web access of VISPA2 ( http://openserver.itb.cnr.it/vispa/ ) ensures accessibility and ease of usage to researches of a complex analytical tool. We released the source code of VISPA2 in a public repository ( https://bitbucket.org/andreacalabria/vispa2 ).

A novel molecular dynamics approach to evaluate the effect of phosphorylation on multimeric protein interface: the αB-Crystallin case study.

BACKGROUND: Phosphorylation is one of the most important post-translational modifications (PTM) employed by cells to regulate several cellular processes. Studying the effects of phosphorylations on protein structures allows to investigate the modulation mechanisms of several proteins including chaperones, like the small HSPs, which display different multimeric structures according to the phosphorylation of a few serine residues. In this context, the proposed study is aimed at finding a method to correlate different PTM patterns (in particular phosphorylations at the monomers interface of multimeric complexes) with the dynamic behaviour of the complex, using physicochemical parameters derived from molecular dynamics simulations in the timescale of nanoseconds. RESULTS: We have developed a methodology relying on computing nine physicochemical parameters, derived from the analysis of short MD simulations, and combined with N identifiers that characterize the PTMs of the analysed protein. The nine general parameters were validated on three proteins, with known post-translational modified conformation and unmodified conformation. Then, we applied this approach to the case study of αB-Crystallin, a chaperone which multimeric state (up to 40 units) is supposed to be controlled by phosphorylation of Ser45 and Ser59. Phosphorylation of serines at the dimer interface induces the release of hexamers, the active state of αB-Crystallin. 30 ns of MD simulation were obtained for each possible combination of dimer phosphorylation state and average values of structural, dynamic, energetic and functional features were calculated on the equilibrated portion of the trajectories. Principal Component Analysis was applied to the parameters and the first five Principal Components, which summed up to 84 % of the total variance, were finally considered. CONCLUSIONS: The validation of this approach on multimeric proteins, which structures were known both modified and unmodified, allowed us to propose a new approach that can be used to predict the impact of PTM patterns in multi-modified proteins using data collected from short molecular dynamics simulations. Analysis on the αB-Crystallin case study clusters together all-P dimers with all-P hexamers and no-P dimer with no-P hexamer and results suggest a great influence of Ser59 phosphorylation on chain B.

PWHATSHAP: efficient haplotyping for future generation sequencing.

BACKGROUND: Haplotype phasing is an important problem in the analysis of genomics information. Given a set of DNA fragments of an individual, it consists of determining which one of the possible alleles (alternative forms of a gene) each fragment comes from. Haplotype information is relevant to gene regulation, epigenetics, genome-wide association studies, evolutionary and population studies, and the study of mutations. Haplotyping is currently addressed as an optimisation problem aiming at solutions that minimise, for instance, error correction costs, where costs are a measure of the confidence in the accuracy of the information acquired from DNA sequencing. Solutions have typically an exponential computational complexity. WHATSHAP is a recent optimal approach which moves computational complexity from DNA fragment length to fragment overlap, i.e., coverage, and is hence of particular interest when considering sequencing technology's current trends that are producing longer fragments. RESULTS: Given the potential relevance of efficient haplotyping in several analysis pipelines, we have designed and engineered PWHATSHAP, a parallel, high-performance version of WHATSHAP. PWHATSHAP is embedded in a toolkit developed in Python and supports genomics datasets in standard file formats. Building on WHATSHAP, PWHATSHAP exhibits the same complexity exploring a number of possible solutions which is exponential in the coverage of the dataset. The parallel implementation on multi-core architectures allows for a relevant reduction of the execution time for haplotyping, while the provided results enjoy the same high accuracy as that provided by WHATSHAP, which increases with coverage. CONCLUSIONS: Due to its structure and management of the large datasets, the parallelisation of WHATSHAP posed demanding technical challenges, which have been addressed exploiting a high-level parallel programming framework. The result, PWHATSHAP, is a freely available toolkit that improves the efficiency of the analysis of genomics information.