IgG4-related disease and B-cell malignancy due to an IKZF1 gain-of-function variant.

BACKGROUND: Monoallelic loss-of-function IKZF1 (IKAROS) variants cause B-cell deficiency or combined immunodeficiency, whereas monoallelic gain-of-function (GOF) IKZF1 variants have recently been reported to cause hypergammaglobulinemia, abnormal plasma cell differentiation, autoimmune and allergic manifestations, and infections. OBJECTIVE: We studied 7 relatives with autoimmune/inflammatory and lymphoproliferative manifestations to identify the immunologic disturbances and the genetic cause of their disease. METHODS: We analyzed biopsy results and performed whole-exome sequencing and immunologic studies. RESULTS: Disease onset occurred at a mean age of 25.2 years (range, 10-64, years). Six patients suffered from autoimmune/inflammatory diseases, 4 had confirmed IG4-related disease (IgG4-RD), and 5 developed B-cell malignancies: lymphoma in 4 and multiple myeloma in the remaining patient. Patients without immunosuppression were not particularly prone to infectious diseases. Three patients suffered from life-threatening coronavirus disease 2019 pneumonia, of whom 1 had autoantibodies neutralizing IFN-α. The recently described IKZF1 GOF p.R183H variant was found in the 5 affected relatives tested and in a 6-year-old asymptomatic girl. Immunologic analysis revealed hypergammaglobulinemia and high frequencies of certain lymphocyte subsets (exhausted B cells, effector memory CD4 T cells, effector memory CD4 T cells that have regained surface expression of CD45RA and CD28-CD57+ CD4+ and CD8+ T cells, TH2, and Tfh2 cells) attesting to immune dysregulation. Partial clinical responses to rituximab and corticosteroids were observed, and treatment with lenalidomide, which promotes IKAROS degradation, was initiated in 3 patients. CONCLUSIONS: Heterozygosity for GOF IKZF1 variants underlies autoimmunity/inflammatory diseases, IgG4-RD, and B-cell malignancies, the onset of which may occur in adulthood. Clinical and immunologic data are similar to those for patients with unexplained IgG4-RD. Patients may therefore benefit from treatments inhibiting pathways displaying IKAROS-mediated overactivity.

Novel genes and sex differences in COVID-19 severity.

Here, we describe the results of a genome-wide study conducted in 11 939 coronavirus disease 2019 (COVID-19) positive cases with an extensive clinical information that were recruited from 34 hospitals across Spain (SCOURGE consortium). In sex-disaggregated genome-wide association studies for COVID-19 hospitalization, genome-wide significance (P < 5 × 10-8) was crossed for variants in 3p21.31 and 21q22.11 loci only among males (P = 1.3 × 10-22 and P = 8.1 × 10-12, respectively), and for variants in 9q21.32 near TLE1 only among females (P = 4.4 × 10-8). In a second phase, results were combined with an independent Spanish cohort (1598 COVID-19 cases and 1068 population controls), revealing in the overall analysis two novel risk loci in 9p13.3 and 19q13.12, with fine-mapping prioritized variants functionally associated with AQP3 (P = 2.7 × 10-8) and ARHGAP33 (P = 1.3 × 10-8), respectively. The meta-analysis of both phases with four European studies stratified by sex from the Host Genetics Initiative (HGI) confirmed the association of the 3p21.31 and 21q22.11 loci predominantly in males and replicated a recently reported variant in 11p13 (ELF5, P = 4.1 × 10-8). Six of the COVID-19 HGI discovered loci were replicated and an HGI-based genetic risk score predicted the severity strata in SCOURGE. We also found more SNP-heritability and larger heritability differences by age (<60 or ≥60 years) among males than among females. Parallel genome-wide screening of inbreeding depression in SCOURGE also showed an effect of homozygosity in COVID-19 hospitalization and severity and this effect was stronger among older males. In summary, new candidate genes for COVID-19 severity and evidence supporting genetic disparities among sexes are provided.

A genome-wide association study of adults with community-acquired pneumonia.

BACKGROUND: Community-acquired pneumonia (CAP) is associated with high morbidity and hospitalization rate. In infectious diseases, host genetics plays a critical role in susceptibility and immune response, and the immune pathways involved are highly dependent on the microorganism and its route of infection. Here we aimed to identify genetic risk loci for CAP using a case-control genome-wide association study (GWAS). METHODS: We performed a GWAS on 3,765 Spanish individuals, including 257 adult patients hospitalized with CAP and 3,508 population controls. Pneumococcal CAP was documented in 30% of patients; the remaining 70% were selected among patients with unidentified microbiological etiology. We tested 7,6 million imputed genotypes using logistic regressions. UK Biobank GWAS of bacterial pneumonia were used for results validation. Subsequently, we prioritized genes and likely causal variants based on Bayesian fine mapping and functional evidence. Imputation and association of classical HLA alleles and amino acids were also conducted. RESULTS: Six independent sentinel variants reached the genome-wide significance (p < 5 × 10-8), three on chromosome 6p21.32, and one for each of the chromosomes 4q28.2, 11p12, and 20q11.22. Only one variant at 6p21.32 was validated in independent GWAS of bacterial and pneumococcal pneumonia. Our analyses prioritized C4orf33 on 4q28.2, TAPBP on 6p21.32, and ZNF341 on 20q11.22. Interestingly, genetic defects of TAPBP and ZNF341 are previously known inborn errors of immunity predisposing to bacterial pneumonia, including pneumococcus and Haemophilus influenzae. Associations were all non-significant for the classical HLA alleles. CONCLUSIONS: We completed a GWAS of CAP and identified four novel risk loci involved in CAP susceptibility.

PCSK6 and Survival in Idiopathic Pulmonary Fibrosis.

Rationale: Idiopathic pulmonary fibrosis (IPF) is a devastating disease characterized by limited treatment options and high mortality. A better understanding of the molecular drivers of IPF progression is needed. Objectives: To identify and validate molecular determinants of IPF survival. Methods: A staged genome-wide association study was performed using paired genomic and survival data. Stage I cases were drawn from centers across the United States and Europe and stage II cases from Vanderbilt University. Cox proportional hazards regression was used to identify gene variants associated with differential transplantation-free survival (TFS). Stage I variants with nominal significance (P < 5 × 10-5) were advanced for stage II testing and meta-analyzed to identify those reaching genome-wide significance (P < 5 × 10-8). Downstream analyses were performed for genes and proteins associated with variants reaching genome-wide significance. Measurements and Main Results: After quality controls, 1,481 stage I cases and 397 stage II cases were included in the analysis. After filtering, 9,075,629 variants were tested in stage I, with 158 meeting advancement criteria. Four variants associated with TFS with consistent effect direction were identified in stage II, including one in an intron of PCSK6 (proprotein convertase subtilisin/kexin type 6) reaching genome-wide significance (hazard ratio, 4.11 [95% confidence interval, 2.54-6.67]; P = 9.45 × 10-9). PCSK6 protein was highly expressed in IPF lung parenchyma. PCSK6 lung staining intensity, peripheral blood gene expression, and plasma concentration were associated with reduced TFS. Conclusions: We identified four novel variants associated with IPF survival, including one in PCSK6 that reached genome-wide significance. Downstream analyses suggested that PCSK6 protein plays a potentially important role in IPF progression.

Admixture mapping analysis reveals differential genetic ancestry associated with Chagas disease susceptibility in the Colombian population.

Chagas disease is an infection caused by the parasite Trypanosoma cruzi, endemic in Latino America. Leveraging the three-way admixture between Native American (AMR), European (EUR) and African (AFR) populations in Latin Americans, we aimed to better understand the genetic basis of Chagas disease by performing an admixture mapping study in a Colombian population. A two-stage study was conducted, and subjects were classified as seropositive and seronegative for T. cruzi. In stage 1, global and local ancestries were estimated using reference data from the 1000 Genomes Project (1KGP), and local ancestry associations were performed by logistic regression models. The AMR ancestry showed a protective association with Chagas disease within the major histocompatibility complex region [Odds ratio (OR) = 0.74, 95% confidence interval (CI) = 0.66-0.83, lowest P-value = 4.53 × 10-8]. The fine mapping assessment on imputed genotypes combining data from stage 1 and 2 from an independent Colombian cohort, revealed nominally associated variants in high linkage disequilibrium with the top signal (rs2032134, OR = 0.93, 95% CI = 0.90-0.97, P-value = 3.54 × 10-4) in the previously associated locus. To assess ancestry-specific adaptive signals, a selective sweep scan in an AMR reference population from 1KGP together with an in silico functional analysis highlighted the Tripartite Motif family and the human leukocyte antigen genes, with crucial role in the immune response against pathogens. Furthermore, these analyses emphasized the macrophages, neutrophils and eosinophils, as key players in the defense against T. cruzi. This first admixture mapping study in Chagas disease provided novel insights underlying the host immune response in the pathogenesis of this neglected disease.

A genome-wide association study of survival in patients with sepsis.

BACKGROUND: Sepsis is a severe systemic inflammatory response to infections that is accompanied by organ dysfunction and has a high mortality rate in adult intensive care units. Most genetic studies have identified gene variants associated with development and outcomes of sepsis focusing on biological candidates. We conducted the first genome-wide association study (GWAS) of 28-day survival in adult patients with sepsis. METHODS: This study was conducted in two stages. The first stage was performed on 687 European sepsis patients from the GEN-SEP network and 7.5 million imputed variants. Association testing was conducted with Cox regression models, adjusting by sex, age, and the main principal components of genetic variation. A second stage focusing on the prioritized genetic variants was performed on 2,063 ICU sepsis patients (1362 European Americans and 701 African-Americans) from the MESSI study. A meta-analysis of results from the two stages was conducted and significance was established at p < 5.0 × 10-8. Whole-blood transcriptomic, functional annotations, and sensitivity analyses were evaluated on the identified genes and variants. FINDINGS: We identified three independent low-frequency variants associated with reduced 28-day sepsis survival, including a missense variant in SAMD9 (hazard ratio [95% confidence interval] = 1.64 [1.37-6.78], p = 4.92 × 10-8). SAMD9 encodes a possible mediator of the inflammatory response to tissue injury. INTERPRETATION: We performed the first GWAS of 28-day sepsis survival and identified novel variants associated with reduced survival. Larger sample size studies are needed to better assess the genetic effects in sepsis survival and to validate the findings.

Transactive Response DNA-Binding Protein (TARDBP/TDP-43) Regulates Cell Permissivity to HIV-1 Infection by Acting on HDAC6.

The transactive response DNA-binding protein (TARDBP/TDP-43) influences the processing of diverse transcripts, including that of histone deacetylase 6 (HDAC6). Here, we assessed TDP-43 activity in terms of regulating CD4+ T-cell permissivity to HIV-1 infection. We observed that overexpression of wt-TDP-43 increased both mRNA and protein levels of HDAC6, resulting in impaired HIV-1 infection independently of the viral envelope glycoprotein complex (Env) tropism. Consistently, using an HIV-1 Env-mediated cell-to-cell fusion model, the overexpression of TDP-43 levels negatively affected viral Env fusion capacity. Silencing of endogenous TDP-43 significantly decreased HDAC6 levels and increased the fusogenic and infection activities of the HIV-1 Env. Using pseudovirus bearing primary viral Envs from HIV-1 individuals, overexpression of wt-TDP-43 strongly reduced the infection activity of Envs from viremic non-progressors (VNP) and rapid progressors (RP) patients down to the levels of the inefficient HIV-1 Envs observed in long-term non-progressor elite controllers (LTNP-EC). On the contrary, silencing endogenous TDP-43 significantly favored the infectivity of primary Envs from VNP and RP individuals, and notably increased the infection of those from LTNP-EC. Taken together, our results indicate that TDP-43 shapes cell permissivity to HIV-1 infection, affecting viral Env fusion and infection capacities by altering the HDAC6 levels and associated tubulin-deacetylase anti-HIV-1 activity.

Admixture mapping of asthma in southwestern Europeans with North African ancestry influences.

The prevalence of asthma symptoms in Canary Islanders, a southwestern European population from Spain, is almost three times higher than the country average. Because the genetic risks identified so far explain <5% of asthma heritability, here we aimed to discover new asthma loci by completing the first admixture mapping study in Canary Islanders leveraging their distinctive genetic makeup, where significant northwest African influences coexist in the European genetic diversity landscape. A 2-stage study was conducted in 1,491 unrelated individuals self-declaring having a Canary Islands origin for the 4 grandparents. Local ancestry estimates were obtained for the shared positions with reference data from putative ancestral populations from Europe, North Africa, and sub-Saharan Africa. Case-control deviations in local ancestry were tested for each ancestry separately using logistic regressions adjusted for principal components, followed by fine-mapping analyses based on imputed genotypes and analyses of the likely deleterious exonic variants. The admixture mapping analysis of asthma detected that local North African ancestry in a locus spanning 365 kb of chromosome 16q23.3 was associated with asthma risk at study-wide significance [lowest P = 1.12 × 10-4; odds ratio (OR) = 2.05; 95% confidence interval (CI) = 1.41-3.00]. Fine-mapping studies identified a variant associated with asthma, and results were replicated in independent samples (rs3852738, OR = 1.34; 95% CI = 1.13-1.59, P = 7.58 × 10-4). Whole exome sequencing data from a subset of individuals revealed an enrichment of likely deleterious variants among asthma cases in 16q23.3, particularly in the phospholipase Cγ2 (PLCG2) gene (P = 3.67 × 10-4). By completing the first mapping study of asthma in admixed populations from Europe, our results revealed a new plausible asthma locus.

NanoDJ: a Dockerized Jupyter notebook for interactive Oxford Nanopore MinION sequence manipulation and genome assembly.

BACKGROUND: The Oxford Nanopore Technologies (ONT) MinION portable sequencer makes it possible to use cutting-edge genomic technologies in the field and the academic classroom. RESULTS: We present NanoDJ, a Jupyter notebook integration of tools for simplified manipulation and assembly of DNA sequences produced by ONT devices. It integrates basecalling, read trimming and quality control, simulation and plotting routines with a variety of widely used aligners and assemblers, including procedures for hybrid assembly. CONCLUSIONS: With the use of Jupyter-facilitated access to self-explanatory contents of applications and the interactive visualization of results, as well as by its distribution into a Docker software container, NanoDJ is aimed to simplify and make more reproducible ONT DNA sequence analysis. The NanoDJ package code, documentation and installation instructions are freely available at https://github.com/genomicsITER/NanoDJ .

Genomic Analyses of Human European Diversity at the Southwestern Edge: Isolation, African Influence and Disease Associations in the Canary Islands.

Despite the genetic resemblance of Canary Islanders to other southern European populations, their geographical isolation and the historical admixture of aborigines (from North Africa) with sub-Saharan Africans and Europeans have shaped a distinctive genetic makeup that likely affects disease susceptibility and health disparities. Based on single nucleotide polymorphism array data and whole genome sequencing (30×), we inferred that the last African admixture took place ∼14 generations ago and estimated that up to 34% of the Canary Islander genome is of recent African descent. The length of regions in homozygosis and the ancestry-related mosaic organization of the Canary Islander genome support the view that isolation has been strongest on the two smallest islands. Furthermore, several genomic regions showed significant and large deviations in African or European ancestry and were significantly enriched in genes involved in prevalent diseases in this community, such as diabetes, asthma, and allergy. The most prominent of these regions were located near LCT and the HLA, two well-known targets of selection, at which 40‒50% of the Canarian genome is of recent African descent according to our estimates. Putative selective signals were also identified in these regions near the SLC6A11-SLC6A1, KCNMB2, and PCDH20-PCDH9 genes. Taken together, our findings provide solid evidence of a significant recent African admixture, population isolation, and adaptation in this part of Europe, with the favoring of African alleles in some chromosome regions. These findings may have medical implications for populations of recent African ancestry.

Modelling the demographic history of human North African genomes points to a recent soft split divergence between populations.

BACKGROUND: North African human populations present a complex demographic scenario due to the presence of an autochthonous genetic component and population substructure, plus extensive gene flow from the Middle East, Europe, and sub-Saharan Africa. RESULTS: We conducted a comprehensive analysis of 364 genomes to construct detailed demographic models for the North African region, encompassing its two primary ethnic groups, the Arab and Amazigh populations. This was achieved through an Approximate Bayesian Computation with Deep Learning (ABC-DL) framework and a novel algorithm called Genetic Programming for Population Genetics (GP4PG). This innovative approach enabled us to effectively model intricate demographic scenarios, utilizing a subset of 16 whole genomes at > 30X coverage. The demographic model suggested by GP4PG exhibited a closer alignment with the observed data compared to the ABC-DL model. Both point to a back-to-Africa origin of North African individuals and a close relationship with Eurasian populations. Results support different origins for Amazigh and Arab populations, with Amazigh populations originating back in Epipaleolithic times, while GP4PG supports Arabization as the main source of Middle Eastern ancestry. The GP4PG model includes population substructure in surrounding populations (sub-Saharan Africa and Middle East) with continuous decaying gene flow after population split. Contrary to ABC-DL, the best GP4PG model does not require pulses of admixture from surrounding populations into North Africa pointing to soft splits as drivers of divergence in North Africa. CONCLUSIONS: We have built a demographic model on North Africa that points to a back-to-Africa expansion and a differential origin between Arab and Amazigh populations.

Genomic epidemiology of the primary methicillin-resistant Staphylococcus aureus clones causing invasive infections in Paraguayan children.

Methicillin-resistant Staphylococcus aureus (MRSA) is one of the major human pathogens. It could carry numerous resistance genes and virulence factors in its genome, some of which are related to the severity of the infection. An observational, descriptive, cross-sectional study was designed to molecularly analyze MRSA isolates that cause invasive infections in Paraguayan children from 2009 to 2013. Ten representative MRSA isolates of the main clonal complex identified were analyzed with short-read paired-end sequencing and assessed for the virulome, resistome, and phylogenetic relationships. All the genetically linked MRSA isolates were recovered from diverse clinical sources, patients, and hospitals at broad gap periods. The pan-genomic analysis of these clones revealed three major and different clonal complexes (CC30, CC5, and CC8), each composed of clones closely related to each other. The CC30 genomes prove to be a successful clone, strongly installed and disseminated throughout our country, and closely related to other CC30 public genomes from the region and the world. The CC5 shows the highest genetic variability, and the CC8 carried the complete arginine catabolic mobile element (ACME), closely related to the USA300-NAE-ACME+, identified as the major cause of CA-MRSA infections in North America. Multiple virulence and resistance genes were identified for the first time in this study, highlighting the complex virulence profiles of MRSA circulating in the country. This study opens a wide range of new possibilities for future projects and trials to improve the existing knowledge on the epidemiology of MRSA circulating in Paraguay. IMPORTANCE: The increasing prevalence of methicillin-resistant Staphylococcus aureus (MRSA) is a public health problem worldwide. The most frequent MRSA clones identified in Paraguay in previous studies (including community and hospital acquired) were the Pediatric (CC5-ST5-IV), the Cordobes-Chilean (CC5-ST5-I), the SouthWest Pacific (CC30-ST30-IV), and the Brazilian (CC8-ST239-III) clones. In this study, the pan-genomic analysis of the most representative MRSA clones circulating in invasive infection in Paraguayan children over the years 2009-2013, such as the CC30-ST30-IV, CC5-ST5-IV, and CC8-ST8-IV, was carried out to evaluate their genetic diversity, their repertoire of virulence factors, and antimicrobial resistance determinants. This revealed multiple virulence and resistance genes, highlighting the complex virulence profiles of MRSA circulating in Paraguay. Our work is the first genomic study of MRSA in Paraguay and will contribute to the development of genomic surveillance in the region and our understanding of the global epidemiology of this pathogen.

Novel risk loci for COVID-19 hospitalization among admixed American populations.

The genetic basis of severe COVID-19 has been thoroughly studied, and many genetic risk factors shared between populations have been identified. However, reduced sample sizes from non-European groups have limited the discovery of population-specific common risk loci. In this second study nested in the SCOURGE consortium, we conducted a genome-wide association study (GWAS) for COVID-19 hospitalization in admixed Americans, comprising a total of 4702 hospitalized cases recruited by SCOURGE and seven other participating studies in the COVID-19 Host Genetic Initiative. We identified four genome-wide significant associations, two of which constitute novel loci and were first discovered in Latin American populations (BAZ2B and DDIAS). A trans-ethnic meta-analysis revealed another novel cross-population risk locus in CREBBP. Finally, we assessed the performance of a cross-ancestry polygenic risk score in the SCOURGE admixed American cohort. This study constitutes the largest GWAS for COVID-19 hospitalization in admixed Latin Americans conducted to date. This allowed to reveal novel risk loci and emphasize the need of considering the diversity of populations in genomic research.

Benchmarking of human Y-chromosomal haplogroup classifiers with whole-genome and whole-exome sequence data.

In anthropological, medical, and forensic studies, the nonrecombinant region of the human Y chromosome (NRY) enables accurate reconstruction of pedigree relationships and retrieval of ancestral information. Using high-throughput sequencing (HTS) data, we present a benchmarking analysis of command-line tools for NRY haplogroup classification. The evaluation was performed using paired Illumina data from whole-genome sequencing (WGS) and whole-exome sequencing (WES) experiments from 50 unrelated donors. Additionally, as a validation, we also used paired WGS/WES datasets of 54 individuals from the 1000 Genomes Project. Finally, we evaluated the tools on data from third-generation HTS obtained from a subset of donors and one reference sample. Our results show that WES, despite typically offering less genealogical resolution than WGS, is an effective method for determining the NRY haplogroup. Y-LineageTracker and Yleaf showed the highest accuracy for WGS data, classifying precisely 98% and 96% of the samples, respectively. Yleaf outperforms all benchmarked tools in the WES data, classifying approximately 90% of the samples. Yleaf, Y-LineageTracker, and pathPhynder can correctly classify most samples (88%) sequenced with third-generation HTS. As a result, Yleaf provides the best performance for applications that use WGS and WES. Overall, our study offers researchers with a guide that allows them to select the most appropriate tool to analyze the NRY region using both second- and third-generation HTS data.

Reinfection rate and disease severity of the BA.5 Omicron SARS-CoV-2 lineage compared to previously circulating variants of concern in the Canary Islands (Spain).

ABSTRACTThe emergence of the Omicron SARS-CoV-2 variant of concern has changed the COVID-19 scenario as this variant is characterized by high transmissibility and immune evasion ability. To evaluate the impact of this variant on the Canary Islands (Spain) population, we determined the reinfection rates and disease severity associated with the Omicron sublineages and the previously circulating variants of concern. We performed a retrospective observational study on 21,745 SARS-CoV-2 viral genomes collected from December 2020 to July 2022 in the Canary Islands (Spain). We compared the reinfection rates between lineages using pairwise proportion and Fisher's exact tests. To assess disease severity, we studied the association of Alpha, Delta, BA.1, BA.2, BA.5, and other risk factors on 28-day hospital mortality using logistic regression and Cox proportional hazard models. We observed 127 bona fide reinfection cases throughout the study period. We found that BA.5 had the highest reinfection rate compared to other lineages (vs. Delta p = 2.89 × 10-25; vs. BA.1 p = 5.17 × 10-11; vs. BA.2 p = 0.002). Among the 1,094 hospitalized patients, multivariate logistic regression showed that Alpha (Odds Ratio [OR] = 0.45, 95% Confidence Interval [CI] = 0.23-0.87, p = 0.02), BA.2 (OR = 0.38, 95% CI = 0.22-0.63, p = 1.91 × 10-4), and BA.5 (OR = 0.30, 95% CI = 0.16-0.55, p = 1.05 × 10-4) had lower 28-day hospital mortality compared to Delta. These results were confirmed by using Cox proportional hazard models. Omicron lineages, and in particular BA.5, were associated with higher reinfection rates and lower disease severity (28-day hospital mortality) than previously circulating variants of concern.

Digging into the admixture strata of current-day Canary Islanders based on mitogenomes.

The conquest of the Canary Islands by Europeans began at the beginning of the 15th century and culminated in 1496 with the surrender of the aborigines. The collapse of the aboriginal population during the conquest and the arrival of settlers caused a drastic change in the demographic composition of the archipelago. To shed light on this historical process, we analyzed 896 mitogenomes of current inhabitants from the seven main islands. Our findings confirm the continuity of aboriginal maternal contributions and the persistence of their genetic footprints in the current population, even at higher levels (>60% on average) than previously evidenced. Moreover, the age estimates for most autochthonous founder lineages support a first aboriginal arrival to the islands at the beginning of the first millennium. We also revealed for the first time that the main recognizable genetic influences from Europe are from Portuguese and Galicians.

Bioinformatic approaches to draft the viral genome sequence of Canary Islands cases related to the multicountry mpox virus 2022-outbreak.

On July 23, 2022, monkeypox disease (mpox) was declared a Public Emergency of International Concern (PHEIC) by the World Health Organization (WHO) due to a multicountry outbreak. In Europe, several cases of mpox virus (MPXV) infection related to this outbreak were detected in the Canary Islands (Spain). Here we describe the combination of viral DNA sequencing and bioinformatic approaches, including methods for de novo genome assembly and short- and long-read technologies, used to reconstruct the first MPXV genome isolated in the Canary Islands on the 31st of May 2022 from a male adult patient with mild symptoms. The same sequencing and bioinformatic approaches were then validated with three other positive cases of MPXV infection from the same mpox outbreak. We obtained the best results using a reference-based approach with short reads, evidencing 46-79 nucleotide variants against viral sequences from the 2018-2019 mpox outbreak and placing the viral sequences in the new B.1 sublineage of clade IIb of the MPXV classification. This study of MPXV demonstrates the potential of metagenomics sequencing for rapid and precise pathogen identification.