APE1 controls DICER1 expression in NSCLC through miR-33a and miR-130b.

Increasing evidence suggests different, not completely understood roles of microRNA biogenesis in the development and progression of lung cancer. The overexpression of the DNA repair protein apurinic/apyrimidinic endodeoxyribonuclease 1 (APE1) is an important cause of poor chemotherapeutic response in lung cancer and its involvement in onco-miRNAs biogenesis has been recently described. Whether APE1 regulates miRNAs acting as prognostic biomarkers of lung cancer has not been investigated, yet. In this study, we analyzed miRNAs differential expression upon APE1 depletion in the A549 lung cancer cell line using high-throughput methods. We defined a signature of 13 miRNAs that strongly correlate with APE1 expression in human lung cancer: miR-1246, miR-4488, miR-24, miR-183, miR-660, miR-130b, miR-543, miR-200c, miR-376c, miR-218, miR-146a, miR-92b and miR-33a. Functional enrichment analysis of this signature revealed its biological relevance in cancer cell proliferation and survival. We validated DICER1 as a direct functional target of the APE1-regulated miRNA-33a-5p and miR-130b-3p. Importantly, IHC analyses of different human tumors confirmed a negative correlation existing between APE1 and Dicer1 protein levels. DICER1 downregulation represents a prognostic marker of cancer development but the mechanisms at the basis of this phenomenon are still completely unknown. Our findings, suggesting that APE1 modulates DICER1 expression via miR-33a and miR-130b, reveal new mechanistic insights on DICER1 regulation, which are of relevance in lung cancer chemoresistance and cancer invasiveness.

HIV-1 Nef promotes infection by excluding SERINC5 from virion incorporation.

HIV-1 Nef, a protein important for the development of AIDS, has well-characterized effects on host membrane trafficking and receptor downregulation. By an unidentified mechanism, Nef increases the intrinsic infectivity of HIV-1 virions in a host-cell-dependent manner. Here we identify the host transmembrane protein SERINC5, and to a lesser extent SERINC3, as a potent inhibitor of HIV-1 particle infectivity that is counteracted by Nef. SERINC5 localizes to the plasma membrane, where it is efficiently incorporated into budding HIV-1 virions and impairs subsequent virion penetration of susceptible target cells. Nef redirects SERINC5 to a Rab7-positive endosomal compartment and thereby excludes it from HIV-1 particles. The ability to counteract SERINC5 was conserved in Nef encoded by diverse primate immunodeficiency viruses, as well as in the structurally unrelated glycosylated Gag from murine leukaemia virus. These examples of functional conservation and convergent evolution emphasize the fundamental importance of SERINC5 as a potent anti-retroviral factor.

ddSeeker: a tool for processing Bio-Rad ddSEQ single cell RNA-seq data.

BACKGROUND: New single-cell isolation technologies are facilitating studies on the transcriptomics of individual cells. Bio-Rad ddSEQ is a droplet-based microfluidic system that, when coupled with downstream Illumina library preparation and sequencing, enables the monitoring of thousands of genes per cell. Sequenced reads show unique features that do not permit the use of freely available tools to perform single cell demultiplexing. RESULTS: We present ddSeeker, a tool to perform initial processing and quality metrics of reads generated through Bio-Rad ddSEQ/Illumina experiments. Its application to the Illumina test dataset demonstrates that ddSeeker performs better than Illumina BaseSpace software, enabling a higher recovery of valid reads. We also show its utility in the analysis of an in-house dataset including two read sets characterized by low and high sequencing quality. ddSeeker and its source code are available at https://github.com/cgplab/ddSeeker . CONCLUSIONS: ddSeeker is a freely available tool to perform initial processing and quality metrics of reads generated through Bio-Rad ddSEQ/Illumina single cell transcriptomic experiments.

Characterization of 17 strains belonging to the Mycobacterium simiae complex and description of Mycobacterium paraense sp. nov.

Fourteen mycobacterial strains isolated from pulmonary samples of independent patients in the state of Pará (Brazil), and three strains isolated in Italy, were characterized using a polyphasic approach. Thorough genetic investigation, including whole-genome sequencing, demonstrated that the strains belong to the M. simiae complex, being most closely related to Mycobacterium interjectum. For 14 of the strains, evidence emerged supporting their inclusion in a previously unreported species of the genus Mycobacterium, for which the name Mycobacterium paraense sp. nov. is proposed (type strain, IEC26(T) = DSM 46749(T) = CCUG 66121(T)). The novel species is characterized by slow growth, unpigmented or pale yellow scotochromogenic colonies, and a HPLC mycolic acid profile different from other known mycobacteria. In different genetic regions, high sequence microheterogeneity was detected.

The DNA-repair protein APE1 participates with hnRNPA2B1 to motif-enriched and prognostic miRNA secretion.

The base excision repair (BER) Apurinic/apyrimidinic endonuclease 1 (APE1) enzyme is endowed with several non-repair activities including miRNAs processing. APE1 is overexpressed in many cancers but its causal role in the tumorigenic processes is largely unknown. We recently described that APE1 can be actively secreted by mammalian cells through exosomes. However, APE1 role in EVs or exosomes is still unknown, especially regarding a putative regulatory function on vesicular small non-coding RNAs. Through dedicated transcriptomic analysis on cellular and vesicular small RNAs of different APE1-depleted cancer cell lines, we found that miRNAs loading into EVs is a regulated process, dependent on APE1, distinctly conveying RNA subsets into vesicles. We identified APE1-dependent secreted miRNAs characterized by enriched sequence motifs and possible binding sites for APE1. In 33 out of 34 APE1-dependent-miRNA precursors, we surprisingly found EXO-motifs and proved that APE1 cooperates with hnRNPA2B1 for the EV-sorting of a subset of miRNAs, including miR-1246, through direct binding to GGAG stretches. Using TCGA-datasets, we showed that these miRNAs identify a signature with high prognostic significance in cancer. In summary, we provided evidence that the ubiquitous DNA-repair enzyme APE1 is part of the EV protein cargo with a novel post-transcriptional role for this ubiquitous DNA-repair enzyme that could explain its role in cancer progression. These findings could open new translational perspectives in cancer biology.

Identification of new p53 target microRNAs by bioinformatics and functional analysis.

BACKGROUND: The tumor suppressor p53 is a sequence-specific transcription factor that regulates an extensive network of coding genes, long non-coding RNAs and microRNAs, that establish intricate gene regulatory circuits influencing many cellular responses beyond the prototypical control of cell cycle, apoptosis and DNA repair. METHODS: Using bioinformatic approaches, we identified an additional group of candidate microRNAs (miRs) under direct p53 transcriptional control. To validate p53 family-mediated responsiveness of the newly predicted target miRs we first evaluated the potential for wild type p53, p63ß and p73ß to transactivate from p53 response elements (REs) mapped in the miR promoters, using an established yeast-based assay. RESULTS: The REs found in miR-10b, -23b, -106a, -151a, -191, -198, -202, -221, -320, -1204, -1206 promoters were responsive to p53 and 8 of them were also responsive to p63ß or p73ß. The potential for germline p53 mutations to drive transactivation at selected miR-associated REs was also examined. Chromatin Immuno-Precipitation (ChIP) assays conducted in doxorubicin-treated MCF7 cells and HCT116 p53+/+ revealed moderate induction of p53 occupancy at the miR-202, -1204, -1206, -10b RE-containing sites, while weak occupancy was observed for the miR-23b-associated RE only in MCF7 cells. RT-qPCR analyses cells showed modest doxorubicin- and/or Nutlin-dependent induction of the levels of mature miR-10b, -23b, -151a in HCT116 p53+/+ and MCF7 cells. The long noncoding RNA PVT1 comprising miR-1204 and -1206 was weakly induced only in HCT116 p53+/+ cells, but the mature miRs were not detected. miR-202 expression was not influenced by p53-activating stimuli in our cell systems. CONCLUSIONS: Our study reveals additional miRs, particularly miR-10b and miR-151a, that could be directly regulated by the p53-family of transcription factors and contribute to the tuning of p53-induced responses.

Alterations of ribosomal RNA pseudouridylation in human breast cancer.

RNA modifications are key regulatory factors for several biological and pathological processes. They are abundantly represented on ribosomal RNA (rRNA), where they contribute to regulate ribosomal function in mRNA translation. Altered RNA modification pathways have been linked to tumorigenesis as well as to other human diseases. In this study we quantitatively evaluated the site-specific pseudouridylation pattern in rRNA in breast cancer samples exploiting the RBS-Seq technique involving RNA bisulfite treatment coupled with a new NGS approach. We found a wide variability among patients at different sites. The most dysregulated positions in tumors turned out to be hypermodified with respect to a reference RNA. As for 2'O-methylation level of rRNA modification, we detected variable and stable pseudouridine sites, with the most stable sites being the most evolutionary conserved. We also observed that pseudouridylation levels at specific sites are related to some clinical and bio-pathological tumor features and they are able to distinguish different patient clusters. This study is the first example of the contribution that newly available high-throughput approaches for site specific pseudouridine detection can provide to the understanding of the intrinsic ribosomal changes occurring in human tumors.

Clinical Significance of Apurinic/Apyrimidinic Endodeoxyribonuclease 1 and MicroRNA Axis in Hepatocellular Carcinoma.

Background and Aims: Identification of prognostic factors for hepatocellular carcinoma (HCC) opens new perspectives for therapy. Circulating and cellular onco-miRNAs are noncoding RNAs which can control the expression of genes involved in oncogenesis through post-transcriptional mechanisms. These microRNAs (miRNAs) are considered novel prognostic and predictive factors in HCC. The apurinic/apyrimidinic endodeoxyribonuclease 1 (APE1) contributes to the quality control and processing of specific onco-miRNAs and is a negative prognostic factor in several tumors. The present work aims to: a) define APE1 prognostic value in HCC; b) identify miRNAs regulated by APE1 and their relative target genes and c) study their prognostic value. Methods: We used The Cancer Genome Atlas (commonly known as TCGA) data analysis to evaluate the expression of APE1 in HCC. To identify differentially-expressed miRNAs (DEmiRNAs) upon APE1 depletion through specific small interfering RNA, we used NGS and nanostring approaches in the JHH-6 HCC tumor cell line. Bioinformatics analyses were performed to identify signaling pathways involving APE1-regulated miRNAs. Microarray analysis was performed to identify miRNAs correlating with serum APE1 expression. Results: APE1 is considerably overexpressed in HCC tissues compared to normal liver, according to the TCGA-liver HCC (known as LIHC) dataset. Enrichment analyses showed that APE1-regulated miRNAs are implicated in signaling and metabolic pathways linked to cell proliferation, transformation, and angiogenesis, identifying Cyclin Dependent Kinase 6 and Lysosomal Associated Membrane Protein 2 as targets. miR-33a-5p, miR-769, and miR-877 are related to lower overall survival in HCC patients. Through array profiling, we identified eight circulating DE-miRNAs associated with APE1 overexpression. A training phase identified positive association between sAPE1 and miR-3180-3p and miR-769. Conclusions: APE1 regulates specific miRNAs having prognostic value in HCC.

Alpha-1 Adrenergic Antagonists Sensitize Neuroblastoma to Therapeutic Differentiation.

Neuroblastoma (NB) is an aggressive childhood tumor, with high-risk cases having a 5-year overall survival probability of approximately 50%. The multimodal therapeutic approach for NB includes treatment with the retinoid isotretinoin (13-cis retinoic acid; 13cRA), which is used in the post-consolidation phase as an antiproliferation and prodifferentiation agent to minimize residual disease and prevent relapse. Through small-molecule screening, we identified isorhamnetin (ISR) as a synergistic compound with 13cRA in inhibiting up to 80% of NB cell viability. The synergistic effect was accompanied by a marked increase in the expression of the adrenergic receptor α1B (ADRA1B) gene. Genetic knockout of ADRA1B or its specific blockade using α1/α1B adrenergic antagonists led to selective sensitization of MYCN-amplified NB cells to cell viability reduction and neural differentiation induced by 13cRA, thus mimicking ISR activity. Administration of doxazosin, a safe α1-antagonist used in pediatric patients, in combination with 13cRA in NB xenografted mice exerted marked control of tumor growth, whereas each drug alone was ineffective. Overall, this study identified the α1B adrenergic receptor as a pharmacologic target in NB, supporting the evaluation of adding α1-antagonists to the post-consolidation therapy of NB to more efficiently control residual disease. SIGNIFICANCE: Targeting α-adrenergic receptors synergizes with isotretinoin to suppress growth and to promote differentiation of neuroblastoma, revealing a combinatorial approach for more effective management of the disease and prevention of relapse.

Shallow Whole-Genome Sequencing of Aedes japonicus and Aedes koreicus from Italy and an Updated Picture of Their Evolution Based on Mitogenomics and Barcoding.

Aedes japonicus and Aedes koreicus are two invasive mosquitoes native to East Asia that are quickly establishing in temperate regions of Europe. Both species are vectors of arboviruses, but we currently lack a clear understanding of their evolution. Here, we present new short-read, shallow genome sequencing of A. japonicus and A. koreicus individuals from northern Italy, which we used for downstream phylogenetic and barcode analyses. We explored associated microbial DNA and found high occurrences of Delftia bacteria in both samples, but neither Asaia nor Wolbachia. We then assembled complete mitogenomes and used these data to infer divergence times estimating the split of A. japonicus from A. koreicus in the Oligocene, which was more recent than that previously reported using mitochondrial markers. We recover a younger age for most other nodes within Aedini and other Culicidae. COI barcoding and phylogenetic analyses indicate that A. japonicus yaeyamensis, A. japonicus amamiensis, and the two A. koreicus sampled from Europe should be considered as separate species within a monophyletic species complex. Our studies further clarify the evolution of A. japonicus and A. koreicus, and indicate the need to obtain whole-genome data from putative species in order to disentangle their complex patterns of evolution.

Suppressing gain-of-function proteins via CRISPR/Cas9 system in SCA1 cells.

SCAs are autosomal dominant neurodegenerative disorders caused by a gain-of-function protein with toxic activities, containing an expanded polyQ tract in the coding region. There are no treatments available to delay the onset, stop or slow down the progression of these pathologies. In this work we focus our attention on SCA1 which is one of the most common genotypes circulating in Italy. Here, we develop a CRISPR/Cas9-based approach to reduce both forms of the ATXN1 protein, normal and mutated with expanded polyQ. We started with the screening of 10 different sgRNAs able to target Exon 8 of the ATXN1 gene. The two most promising sgRNAs were validated in fibroblasts isolated from SCA1 patients, following the identification of the best transfection method for this type of cell. Our silencing approach significantly downregulated the expression of ataxin1, due to large deletions and the introduction of small changes in the ATXN1 gene, evidenced by NGS analysis, without major effects on cell viability. Furthermore, very few significant guide RNA-dependent off-target effects were observed. These preliminary results not only allowed us to identify the best transfection method for SCA1 fibroblasts, but strongly support CRISPR/Cas9 as a promising approach for the treatment of expanded polyQ diseases. Further investigations will be needed to verify the efficacy of our silencing system in SCA1 neurons and animal models.

Human dyskerin binds to cytoplasmic H/ACA-box-containing transcripts affecting nuclear hormone receptor dependence.

BACKGROUND: Dyskerin is a nuclear protein involved in H/ACA box snoRNA-guided uridine modification of RNA. In humans, its defective function is associated with cancer development and induces specific post-transcriptional alterations of gene expression. In this study, we seek to unbiasedly identify mRNAs regulated by dyskerin in human breast cancer-derived cells. RESULTS: We find that dyskerin depletion affects the expression and the association with polysomes of selected mRNA isoforms characterized by the retention of H/ACA box snoRNA-containing introns. These snoRNA retaining transcripts (snoRTs) are bound by dyskerin in the cytoplasm in the form of shorter 3' snoRT fragments. We then characterize the whole cytoplasmic dyskerin RNA interactome and find both H/ACA box snoRTs and protein-coding transcripts which may be targeted by the snoRTs' guide properties. Since a fraction of these protein-coding transcripts is involved in the nuclear hormone receptor binding, we test to see if this specific activity is affected by dyskerin. Obtained results indicate that dyskerin dysregulation may alter the dependence on nuclear hormone receptor ligands in breast cancer cells. These results are paralleled by consistent observations on the outcome of primary breast cancer patients stratified according to their tumor hormonal status. Accordingly, experiments in nude mice show that the reduction of dyskerin levels in estrogen-dependent cells favors xenograft development in the absence of estrogen supplementation. CONCLUSIONS: Our work suggests a cytoplasmic function for dyskerin which could affect mRNA post-transcriptional networks relevant for nuclear hormone receptor functions.

Next Generation Sequencing Analysis of MODY-X Patients: A Case Report Series.

BACKGROUND: Classic criteria for a maturity-onset diabetes of the young (MODY) diagnosis are often unable to identify all subjects, and traditional Sanger sequencing, using a candidate gene approach, leads to a high prevalence of missed genetic diagnosis, classified as MODY-X. Next generation sequencing (NGS) panels provide a highly sensitive method even for rare forms. METHODS: We investigated 28 pediatric subjects suspected for MODY-X, utilizing a 15-gene NGS panel for monogenic diabetes (MD). RESULTS: NGS detected variants of uncertain significance (VUS), likely pathogenic or pathogenic for rarer subtypes of MODY, in six patients. We found variants in the wolframin gene (WFS1), traditionally not considered in MD genetic screening panels, in three patients; KCNJ11 gene mutation, typically responsible for neonatal diabetes and rarely causing isolated diabetes in adolescents; INS gene mutation; a variant in the HNF1B gene in a young male with diabetes on sulfonylurea treatment. CONCLUSION: In our cohort, the availability of an NGS panel for MD was determined for the correct identification of MD subtypes in six patients with MODY-X. Our study underlines how a precise diagnosis utilizing NGS may have an impact on the management of different forms of MODY and, thus, lead to a tailored treatment and enable genetic counselling of other family members.

YES1 and MYC Amplifications as Synergistic Resistance Mechanisms to Different Generation ALK Tyrosine Kinase Inhibitors in Advanced NSCLC: Brief Report of Clinical and Preclinical Proofs.

INTRODUCTION: ALK tyrosine kinase inhibitors (TKIs) are the standard treatment for advanced ALK-positive NSCLC. Nevertheless, drug resistance inevitably occurs. Here, we report a case of a patient with metastatic ALK-positive lung adenocarcinoma with an impressive resistance to sequential treatment with ALK TKIs mediated by YES1 and MYC amplification in a contest of epithelial-to-mesenchymal transition and high progressive chromosomal instability. METHODS: The patient received, after chemotherapy and 7 months of crizotinib, brigatinib and lorlatinib with no clinical benefit to both treatments. A study of resistance mechanisms was performed with whole exome sequencing on different biological samples; primary cell lines were established from pleural effusion after lorlatinib progression. RESULTS: At whole exome sequencing analysis, YES1 and MYC amplifications were observed both in the pericardial biopsy and the pleural effusion samples collected at brigatinib and lorlatinib progression, respectively. Increasing chromosomal instability from diagnostic biopsy to pleural effusion was also observed. The addition of dasatinib to brigatinib or lorlatinib restored the sensitivity in primary cell lines; data were confirmed also in H3122_ALK-positive model overexpressing both YES1 and MYC. CONCLUSIONS: In conclusion, YES1 and MYC amplifications are candidates to justify a rapid acquired resistance to crizotinib entailing primary brigatinib and lorlatinib resistance. In this context, a combination strategy of ALK TKI with dasatinib could be effective to overcome a rapid resistance.

A Community Study of SARS-CoV-2 Detection by RT-PCR in Saliva: A Reliable and Effective Method.

Efficient, wide-scale testing for SARS-CoV-2 is crucial for monitoring the incidence of the infection in the community. The gold standard for COVID-19 diagnosis is the molecular analysis of epithelial secretions from the upper respiratory system captured by nasopharyngeal (NP) or oropharyngeal swabs. Given the ease of collection, saliva has been proposed as a possible substitute to support testing at the population level. Here, we used a novel saliva collection device designed to favour the safe and correct acquisition of the sample, as well as the processivity of the downstream molecular analysis. We tested 1003 nasopharyngeal swabs and paired saliva samples self-collected by individuals recruited at a public drive-through testing facility. An overall moderate concordance (68%) between the two tests was found, with evidence that neither system can diagnose the infection in 100% of the cases. While the two methods performed equally well in symptomatic individuals, their discordance was mainly restricted to samples from convalescent subjects. The saliva test was at least as effective as NP swabs in asymptomatic individuals recruited for contact tracing. Our study describes a testing strategy of self-collected saliva samples, which is reliable for wide-scale COVID-19 screening in the community and is particularly effective for contact tracing.

The altered transcriptome of pediatric myelodysplastic syndrome revealed by RNA sequencing.

Pediatric myelodysplastic syndrome (PMDS) is a very rare and still poorly characterized disorder. In this work, we identified novel potential targets of PMDS by determining genes with aberrant expression, which can be correlated with PMDS pathogenesis. We identified 291 differentially expressed genes (DEGs) in PMDS patients, comprising genes involved in the regulation of apoptosis and the cell cycle, ribosome biogenesis, inflammation and adaptive immunity. Ten selected DEGs were then validated, confirming the sequencing data. These DEGs will potentially represent new molecular biomarkers and therapeutic targets for PMDS.

The short-term impact of probiotic consumption on the oral cavity microbiome.

The dysbiosis of the oral microbiome is associated with both localized and systemic diseases. Modulating the resident microbial communities by the dietary consumption of probiotics has become an appealing means to promote host health by either restoring host-microbe balance or preventing dysbiosis. Most probiotics strategies target the intestinal microbiome, but little is known about their impact on the oral microbiome. We analyzed here the saliva microbiome from 21 volunteers, longitudinally collected before, during, and after consumption of a commercial probiotic and a standard yoghurt using 16S amplicon sequencing. The alpha diversity of the saliva microbiome had a statistically significant increase (P-value = 0.0011) in one of the groups that consumed the probiotic. The overall structure of the microbiome was however not significantly impacted by the probiotic, although oligotyping analysis revealed that both Streptococci and Lactobacilli present in the probiotic product persisted in the saliva microbiome. In contrast, non-probiotic yoghurt consumption had a lesser impact on the overall diversity and Lactobacillus and Streptococcus persistence. Our results suggest that consumption of commercial probiotics in healthy subjects increase the overall diversity of the oral cavity microbiome in the short term, but such dietary interventions are not able to substantially modify the structure of the microbiome.

Mother-to-Infant Microbial Transmission from Different Body Sites Shapes the Developing Infant Gut Microbiome.

The acquisition and development of the infant microbiome are key to establishing a healthy host-microbiome symbiosis. The maternal microbial reservoir is thought to play a crucial role in this process. However, the source and transmission routes of the infant pioneering microbes are poorly understood. To address this, we longitudinally sampled the microbiome of 25 mother-infant pairs across multiple body sites from birth up to 4 months postpartum. Strain-level metagenomic profiling showed a rapid influx of microbes at birth followed by strong selection during the first few days of life. Maternal skin and vaginal strains colonize only transiently, and the infant continues to acquire microbes from distinct maternal sources after birth. Maternal gut strains proved more persistent in the infant gut and ecologically better adapted than those acquired from other sources. Together, these data describe the mother-to-infant microbiome transmission routes that are integral in the development of the infant microbiome.

Unexplored diversity and strain-level structure of the skin microbiome associated with psoriasis.

Psoriasis is an immune-mediated inflammatory skin disease that has been associated with cutaneous microbial dysbiosis by culture-dependent investigations and rRNA community profiling. We applied, for the first time, high-resolution shotgun metagenomics to characterise the microbiome of psoriatic and unaffected skin from 28 individuals. We demonstrate psoriatic ear sites have a decreased diversity and psoriasis is associated with an increase in Staphylococcus, but overall the microbiomes of psoriatic and unaffected sites display few discriminative features at the species level. Finer strain-level analysis reveals strain heterogeneity colonisation and functional variability providing the intriguing hypothesis of psoriatic niche-specific strain adaptation or selection. Furthermore, we accessed the poorly characterised, but abundant, clades with limited sequence information in public databases, including uncharacterised Malassezia spp. These results highlight the skins hidden diversity and suggests strain-level variations could be key determinants of the psoriatic microbiome. This illustrates the need for high-resolution analyses, particularly when identifying therapeutic targets. This work provides a baseline for microbiome studies in relation to the pathogenesis of psoriasis.

Genomic characterization of Nontuberculous Mycobacteria.

Mycobacterium tuberculosis and Mycobacterium leprae have remained, for many years, the primary species of the genus Mycobacterium of clinical and microbiological interest. The other members of the genus, referred to as nontuberculous mycobacteria (NTM), have long been underinvestigated. In the last decades, however, the number of reports linking various NTM species with human diseases has steadily increased and treatment difficulties have emerged. Despite the availability of whole genome sequencing technologies, limited effort has been devoted to the genetic characterization of NTM species. As a consequence, the taxonomic and phylogenetic structure of the genus remains unsettled and genomic information is lacking to support the identification of these organisms in a clinical setting. In this work, we widen the knowledge of NTMs by reconstructing and analyzing the genomes of 41 previously uncharacterized NTM species. We provide the first comprehensive characterization of the genomic diversity of NTMs and open new venues for the clinical identification of opportunistic pathogens from this genus.