The genomes of the yaws bacterium, Treponema pallidum subsp. pertenue, of nonhuman primate and human origin are not genomically distinct.

BACKGROUND: Treponema pallidum subsp. pertenue (TPE) is the causative agent of human yaws. Yaws is currently reported in 13 endemic countries in Africa, southern Asia, and the Pacific region. During the mid-20th century, a first yaws eradication effort resulted in a global 95% drop in yaws prevalence. The lack of continued surveillance has led to the resurgence of yaws. The disease was believed to have no animal reservoirs, which supported the development of a currently ongoing second yaws eradication campaign. Concomitantly, genetic evidence started to show that TPE strains naturally infect nonhuman primates (NHPs) in sub-Saharan Africa. In our current study we tested hypothesis that NHP- and human-infecting TPE strains differ in the previously unknown parts of the genomes. METHODOLOGY/PRINCIPAL FINDINGS: In this study, we determined complete (finished) genomes of ten TPE isolates that originated from NHPs and compared them to TPE whole-genome sequences from human yaws patients. We performed an in-depth analysis of TPE genomes to determine if any consistent genomic differences are present between TPE genomes of human and NHP origin. We were able to resolve previously undetermined TPE chromosomal regions (sequencing gaps) that prevented us from making a conclusion regarding the sequence identity of TPE genomes from NHPs and humans. The comparison among finished genome sequences revealed no consistent differences between human and NHP TPE genomes. CONCLUSION/SIGNIFICANCE: Our data show that NHPs are infected with strains that are not only similar to the strains infecting humans but are genomically indistinguishable from them. Although interspecies transmission in NHPs is a rare event and evidence for current spillover events is missing, the existence of the yaws bacterium in NHPs is demonstrated. While the low risk of spillover supports the current yaws treatment campaign, it is of importance to continue yaws surveillance in areas where NHPs are naturally infected with TPE even if yaws is successfully eliminated in humans.

Complete genome sequences of five Escherichia coli strains with probiotic attributes.

The complete genome sequences of five Escherichia coli strains with probiotic attributes were determined, including strain A0 34/86, a component of the probiotic product Colinfant New Born, and strains H22, 582, B771, and B1172 with published probiotic potential. The size of sequenced genomes ranged from 5,092 to 5,408 kb.

A workflow to study mechanistic indicators for driver gene prediction with Moonlight.

Prediction of driver genes (tumor suppressors and oncogenes) is an essential step in understanding cancer development and discovering potential novel treatments. We recently proposed Moonlight as a bioinformatics framework to predict driver genes and analyze them in a system-biology-oriented manner based on -omics integration. Moonlight uses gene expression as a primary data source and combines it with patterns related to cancer hallmarks and regulatory networks to identify oncogenic mediators. Once the oncogenic mediators are identified, it is important to include extra levels of evidence, called mechanistic indicators, to identify driver genes and to link the observed gene expression changes to the underlying alteration that promotes them. Such a mechanistic indicator could be for example a mutation in the regulatory regions for the candidate gene. Here, we developed new functionalities and released Moonlight2 to provide the user with a mutation-based mechanistic indicator as a second layer of evidence. These functionalities analyze mutations in a cancer cohort to classify them into driver and passenger mutations. Those oncogenic mediators with at least one driver mutation are retained as the final set of driver genes. We applied Moonlight2 to the basal-like breast cancer subtype, lung adenocarcinoma and thyroid carcinoma using data from The Cancer Genome Atlas. For example, in basal-like breast cancer, we found four oncogenes (COPZ2, SF3B4, KRTCAP2 and POLR2J) and nine tumor suppressor genes (KIR2DL4, KIF26B, ARL15, ARHGAP25, EMCN, GMFG, TPK1, NR5A2 and TEK) containing a driver mutation in their promoter region, possibly explaining their deregulation. Moonlight2R is available at https://github.com/ELELAB/Moonlight2R.

Deep oncopanel sequencing reveals within block position-dependent quality degradation in FFPE processed samples.

BACKGROUND: Clinical laboratories routinely use formalin-fixed paraffin-embedded (FFPE) tissue or cell block cytology samples in oncology panel sequencing to identify mutations that can predict patient response to targeted therapy. To understand the technical error due to FFPE processing, a robustly characterized diploid cell line was used to create FFPE samples with four different pre-tissue processing formalin fixation times. A total of 96 FFPE sections were then distributed to different laboratories for targeted sequencing analysis by four oncopanels, and variants resulting from technical error were identified. RESULTS: Tissue sections that fail more frequently show low cellularity, lower than recommended library preparation DNA input, or target sequencing depth. Importantly, sections from block surfaces are more likely to show FFPE-specific errors, akin to "edge effects" seen in histology, while the inner samples display no quality degradation related to fixation time. CONCLUSIONS: To assure reliable results, we recommend avoiding the block surface portion and restricting mutation detection to genomic regions of high confidence.

Evolution of TP53 abnormalities during CLL disease course is associated with telomere length changes.

BACKGROUND: Telomeres are protective structures at chromosome ends which shorten gradually with increasing age. In chronic lymphocytic leukemia (CLL), short telomeres have been associated with unfavorable disease outcome, but the link between clonal evolution and telomere shortening remains unresolved. METHODS: We investigated relative telomere length (RTL) in a well-characterized cohort of 198 CLL patients by qPCR and focused in detail on a subgroup 26 patients who underwent clonal evolution of TP53 mutations (evolTP53). In the evolTP53 subgroup we explored factors influencing clonal evolution and corresponding changes in telomere length through measurements of telomerase expression, lymphocyte doubling time, and BCR signaling activity. RESULTS: At baseline, RTL of the evolTP53 patients was scattered across the entire RTL spectrum observed in our CLL cohort. RTL changed in the follow-up samples of 16/26 (62%) evolTP53 cases, inclining to reach intermediate RTL values, i.e., longer telomeres shortened compared to baseline while shorter ones prolonged. For the first time we show that TP53 clonal shifts are linked to RTL change, including unexpected RTL prolongation. We further investigated parameters associated with RTL changes. Unstable telomeres were significantly more frequent among younger patients (P = 0.032). Shorter telomeres were associated with decreased activity of the B-cell receptor signaling components p-ERK1/2, p-ZAP-70/SYK, and p-NFκB (P = 0.04, P = 0.01, and P = 0.02, respectively). CONCLUSIONS: Our study revealed that changes of telomere length reflect evolution in leukemic subclone proportion, and are associated with specific clinico-biological features of the explored cohort.

Hierarchical distribution of somatic variants in newly diagnosed chronic myeloid leukaemia at diagnosis and early follow-up.

There is an emerging body of evidence that patients with chronic myeloid leukaemia (CML) may carry not only breakpoint cluster region-Abelson murine leukaemia viral oncogene homologue 1 (BCR-ABL1) kinase domain mutations (BCR-ABL1 KD mutations), but also mutations in other genes. Their occurrence is highest during progression or at failure, but their impact at diagnosis is unclear. In the present study, we prospectively screened for mutations in 18 myeloid neoplasm-associated genes and BCR-ABL1 KD in the following populations: bulk leucocytes, CD34+ CD38+ progenitors and CD34+ CD38- stem cells, at diagnosis and early follow-up. In our cohort of chronic phase CML patients, nine of 49 patients harboured somatic mutations in the following genes: six ASXL1 mutations, one SETBP1, one TP53, one JAK2, but no BCR-ABL1 KD mutations. In seven of the nine patients, mutations were detected in multiple hierarchical populations including bulk leucocytes at diagnosis. The mutation dynamics reflected the BCR-ABL1 transcript decline induced by treatment in eight of the nine cases, suggesting that mutations were acquired in the Philadelphia chromosome (Ph)-positive clone. In one patient, the JAK2 V617F mutation correlated with a concomitant Ph-negative myeloproliferative neoplasm and persisted despite a 5-log reduction of the BCR-ABL1 transcript. Only two of the nine patients with mutations failed first-line therapy. No correlation was found between the mutation status and survival or response outcomes.

Low-burden TP53 mutations in CLL: clinical impact and clonal evolution within the context of different treatment options.

Patients with chronic lymphocytic leukemia (CLL) bearing TP53 mutations experience chemorefractory disease and are therefore candidates for targeted therapy. However, the significance of low-burden TP53 mutations with <10% variant allele frequency (VAF) remains a matter for debate. Herein, we describe clonal evolution scenarios of low-burden TP53 mutations, the clinical impact of which we analyzed in a "real-world" CLL cohort. TP53 status was assessed by targeted next-generation sequencing (NGS) in 511 patients entering first-line treatment with chemo- and/or immunotherapy and 159 patients in relapse before treatment with targeted agents. Within the pretherapy cohort, 16% of patients carried low-burden TP53 mutations (0.1% to 10% VAF). Although their presence did not significantly shorten event-free survival after first-line therapy, it affected overall survival (OS). In a subgroup with TP53 mutations of 1% to 10% VAF, the impact on OS was observed only in patients with unmutated IGHV who had not received targeted therapy, as patients benefited from switching to targeted agents, regardless of initial TP53 mutational status. Analysis of the clonal evolution of low-burden TP53 mutations showed that the highest expansion rates were associated with fludarabine, cyclophosphamide, and rituximab regimen in both first- and second-line treatments (median VAF increase, 14.8× and 11.8×, respectively) in contrast to treatment with less intense treatment regimens (1.6×) and no treatment (0.8×). In the relapse cohort, 33% of patients carried low-burden TP53 mutations, which did not expand significantly upon targeted treatment (median VAF change, 1×). Sporadic cases of TP53 mutations' clonal shifts were connected with the development of resistance-associated mutations. Altogether, our data support the incorporation of low-burden TP53 variants in clinical decision making.

Novel Illumina-based next generation sequencing approach with one-round amplification provides early and reliable detection of BCR-ABL1 kinase domain mutations in chronic myeloid leukemia.

The occurrence of mutations in the BCR-ABL1 kinase domain (KD) can lead to treatment resistance in chronic myeloid leukaemia patients. Nowadays, next-generation sequencing (NGS) is an alternative method for the detection of kinase domain mutations, compared to routinely used Sanger sequencing, providing a higher sensitivity of mutation detection. However, in the protocols established so far multiple rounds of amplification limit reliable mutation detection to approximately 5% variant allele frequency. Here, we present a simplified, one-round amplification NGS protocol for the Illumina platform, which offers a robust early detection of BCR-ABL1 KD mutations with a reliable detection limit of 3% variant allele frequency.

Clonal hierarchy of main molecular lesions in acute myeloid leukaemia.

Genetic mutations in acute myeloid leukaemia (AML) are assumed to occur in a sequential order; however, the predominant hierarchical roles of specific mutated genes have not been fully described. In this study, we aimed to determine the clonal involvement of the most frequent AML-associated mutations. Using a targeted sequencing panel for 18 genes, we traced changes and relative clonal contribution of mutations in 52 patients. We analysed 35 pairs of diagnosis and relapse samples, 27 pairs of primary samples and corresponding patient-derived xenografts, and 34 pairs of total leukocytes and corresponding isolated primitive cells or blast populations. In both relapse and xenografts, we observed conservation of main leukaemic clones and variability was limited to subclones with late-acquired mutations. AML evolution thus mainly involved modification of subclones while the clonal background remained unchanged. NPM1 mutations were identified as the most probable leukaemia-transformation lesion, remaining conserved in contrast to high variation of accompanying subclonal FLT3 and NRAS mutations. DNMT3A mutations represented the most stable mutations forming a preleukaemic background in most samples. Mutations in genes IDH1/2, TET2, RUNX1, ASXL1 and U2AF1 were detected both as preleukaemic and as subclonal lesions, suggesting a non-specific order of acquisition.

Impact of gene mutations and chromosomal aberrations on progression-free survival in chronic lymphocytic leukemia patients treated with front-line chemoimmunotherapy: Clinical practice experience.

The impact of genetic aberrations on rituximab-based therapeutic regimens has been intensely studied in chronic lymphocytic leukemia (CLL). According to the current consensus chemoimmunotherapy consisting of rituximab and DNA-damaging drugs is not suitable for patients with TP53 defects. In our study, we focused on CLL patients with an intact TP53 gene and investigated four recurrently mutated genes in CLL, genomic aberrations by FISH, and IGHV status with the aim of analyzing their impact on progression-free survival (PFS) after front-line therapy with FCR (fludarabine, cyclophosphamide, rituximab) or BR (bendamustine, rituximab) regimens. Using next-generation sequencing, we analyzed 120 patients treated with FCR and 57 patients treated with BR at a university hospital. We used a 10% cut-off for variant allele frequency and recorded the following mutation frequencies in the pre-therapy samples: ATM 23%, SF3B1 20%, NOTCH1 19% and BIRC3 11%. The data on cytogenetic aberrations (11q22, 13q14, trisomy 12) and IGHV mutation status were also considered in PFS analyses. In univariate analyses, we observed a negative impact of BIRC3 mutations and 11q22 deletion in both regimens, while the unmutated IGHV status was associated with a significantly shorter PFS only in the FCR-treated cohort. In a multivariate analysis, only deletion 11q22 in both regimens, and the unmutated IGHV in the FCR cohort maintained an independent association with the reduced PFS. Notably, sole 11q22 deletion, without an ATM mutation on the other allele, manifested the shortest PFS of all analyzed markers. Deletion 11q22 and IGHV status predict PFS in previously untreated CLL patients.

Prognostic significance of mutation profile at diagnosis and mutation persistence during disease remission in adult acute myeloid leukaemia patients.

In this multi-centre study, we analysed the prognostic impact of mutations in 19 genes associated with myeloid malignancies in 258 newly diagnosed acute myeloid leukaemia patients (aged 19-70 years) undergoing intensive therapy. We identified five patient groups with different prognostic risks and different benefits from allogeneic hematopoietic stem cell transplantation (alloHSCT) within the intermediate cytogenetic risk group patients (n = 184). The most adverse prognosis was observed in patients with DNMT3A and FLT3-ITD co-mutation, whose survival could be significantly improved with alloHSCT. In contrast, the most favourable prognosis without any further benefit from alloHSCT was identified in patients with mutations in NPM1 or CEBPA, after exclusion of the unfavourable prognostic groups defined by mutations in DNMT3A, RUNX1 or genes from chromatin/spliceosome group. An additional analysis of 113 diagnosis-remission paired samples revealed that persistence of non-DNMT3A mutations (above 2% VAF) represented a further negative prognostic factor. The proposed model offers a possible molecular stratification and treatment guidance for intermediate cytogenetic risk group patients.

ATM and TP53 mutations show mutual exclusivity but distinct clinical impact in mantle cell lymphoma patients.

Mantle cell lymphoma (MCL) is characterized by the hallmark t(11;14)(q13;q32) translocation, leading to cyclin D1 over-expression. Additionally, disrupting the DNA damage response pathway through ATM or TP53 defects plays an important role in MCL pathogenesis. Using deep next-generation sequencing we analyzed the mutual composition of ATM and TP53 mutations in 72 MCL patients, and assessed their impact on progression-free survival (PFS) and overall survival (OS). Mutated ATM and TP53 alleles were found in 51% (37/72) and 22% (16/72) of the cases examined, respectively, with only three patients harboring mutations in both genes. Only a mutated TP53 gene was associated with the significantly reduced PFS and OS and the same output was observed when ATM and TP53 defective groups included also sole deletions 11q and 17p, respectively. Determining the exact ATM/p53 pathway dysfunction may influence the selection of MCL patients for innovative therapies based on the targeted inhibition of selected proteins.

Activation-induced deaminase and its splice variants associate with trisomy 12 in chronic lymphocytic leukemia.

Activation-induced cytidine deaminase (AID) is a mutator enzyme essential for somatic hypermutation (SHM) and class switch recombination (CSR) during effective adaptive immune responses. Its aberrant expression and activity have been detected in lymphomas, leukemias, and solid tumors. In chronic lymphocytic leukemia (CLL) increased expression of alternatively spliced AID variants has been documented. We used real-time RT-PCR to quantify the expression of AID and its alternatively spliced transcripts (AIDΔE4a, AIDΔE4, AIDivs3, and AIDΔE3E4) in 149 CLL patients and correlated this expression to prognostic markers including recurrent chromosomal aberrations, the presence of complex karyotype, mutation status of the immunoglobulin heavy chain variable gene, and recurrent mutations. We report a previously unappreciated association between higher AID transcript levels and trisomy of chromosome 12. Functional analysis of AID splice variants revealed loss of their activity with respect to SHM, CSR, and induction of double-strand DNA breaks. In silico modeling provided insight into the molecular interactions and structural dynamics of wild-type AID and a shortened AID variant closely resembling AIDΔE4, confirming its loss-of-function phenotype.

Detection and Functional Analysis of TP53 Mutations in CLL.

Chronic lymphocytic leukemia (CLL) represents a prototype disease in which TP53 gene defects lead to inferior prognosis. Here, we present two distinct methodologies which can be used to identify TP53 mutations in CLL patients; both protocols are primarily intended for research purposes. The functional analysis of separated alleles in yeast (FASAY) can be flexibly adapted to a variable number of samples and provides an immediate functional readout of identified mutations. Amplicon-based next-generation sequencing then allows for a high throughput and accurately detects subclonal TP53 variants (sensitivity <1% of mutated cells).

ToTem: a tool for variant calling pipeline optimization.

BACKGROUND: High-throughput bioinformatics analyses of next generation sequencing (NGS) data often require challenging pipeline optimization. The key problem is choosing appropriate tools and selecting the best parameters for optimal precision and recall. RESULTS: Here we introduce ToTem, a tool for automated pipeline optimization. ToTem is a stand-alone web application with a comprehensive graphical user interface (GUI). ToTem is written in Java and PHP with an underlying connection to a MySQL database. Its primary role is to automatically generate, execute and benchmark different variant calling pipeline settings. Our tool allows an analysis to be started from any level of the process and with the possibility of plugging almost any tool or code. To prevent an over-fitting of pipeline parameters, ToTem ensures the reproducibility of these by using cross validation techniques that penalize the final precision, recall and F-measure. The results are interpreted as interactive graphs and tables allowing an optimal pipeline to be selected, based on the user's priorities. Using ToTem, we were able to optimize somatic variant calling from ultra-deep targeted gene sequencing (TGS) data and germline variant detection in whole genome sequencing (WGS) data. CONCLUSIONS: ToTem is a tool for automated pipeline optimization which is freely available as a web application at https://totem.software .

The influence of mutational status and biological characteristics of acute myeloid leukemia on xenotransplantation outcomes in NOD SCID gamma mice.

PURPOSE: This study aimed at analyzing the association of gene mutations and other acute myeloid leukemia (AML) characteristics with engraftment outcomes in immunodeficient mice and to select the engraftment outcomes that best reflect patient survival. METHODS: Mutations in 19 genes as well as leukemia- and patient-related characteristics were analyzed for a group of 47 de novo AML samples with respect to three engraftment outcomes: engraftment ability, engraftment intensity (percentage of hCD45+ cells) and engraftment latency. Leukemia-related characteristics were additionally analyzed in an extended group of 68 samples that included the 47 de novo samples, and additional 21 samples from refractory and relapsed cases. Engraftment outcomes were compared with overall and event-free survival of the patients. RESULTS: For the 47 de novo samples, no single mutation influenced engraftment, whereas the NPM1 mut /DNMT3A mut co-mutation was associated with higher engraftment ability. NPM1 mut /FLT3-ITD neg had lower engraftment intensity. Among leukemia-related characteristics, a complex karyotype was associated with higher engraftment intensity. Among patient-related characteristics, higher cytogenetic risk was associated with higher engraftment intensity, and failure to achieve clinical remission was associated with shorter engraftment latency. In the extended group of 68 samples, white blood count was associated with higher engraftment ability, and the presence of a complex karyotype was associated with higher engraftment intensity. Association with patient overall survival was seen only for engraftment intensity. CONCLUSIONS: The engraftment of AML was influenced by mutation-interactions and other AML characteristics, rather than by single mutated genes, and engraftment intensity best reflected clinical penetrance of AML.

Transcription factor YY1 can control AID-mediated mutagenesis in mice.

Activation-induced cytidine deminase (AID) is crucial for controlling the immunoglobulin (Ig) diversification processes of somatic hypermutation (SHM) and class switch recombination (CSR). AID initiates these processes by deamination of cytosine, ultimately resulting in mutations or double strand DNA breaks needed for SHM and CSR. Levels of AID control mutation rates, and off-target non-Ig gene mutations can contribute to lymphomagenesis. Therefore, factors that control AID levels in the nucleus can regulate SHM and CSR, and may contribute to disease. We previously showed that transcription factor YY1 can regulate the level of AID in the nucleus and Ig CSR. Therefore, we hypothesized that conditional knock-out of YY1 would lead to reduction in AID localization at the Ig locus, and reduced AID-mediated mutations. Using mice that overexpress AID (IgκAID yy1f/f ) or that express normal AID levels (yy1f/f ), we found that conditional knock-out of YY1 results in reduced AID nuclear levels, reduced localization of AID to the Sµ switch region, and reduced AID-mediated mutations. We find that the mechanism of YY1 control of AID nuclear accumulation is likely due to YY1-AID physical interaction which blocks AID ubiquitination.

TP53 mutation analysis in chronic lymphocytic leukemia: comparison of different detection methods.

TP53 gene defects represent a strong adverse prognostic factor for patient survival and treatment resistance in chronic lymphocytic leukemia (CLL). Although various methods for TP53 mutation analysis have been reported, none of them allow the identification of all occurring sequence variants, and the most suitable methodology is still being discussed. The aim of this study was to determine the limitations of commonly used methods for TP53 mutation examination in CLL and propose an optimal approach for their detection. We examined 182 CLL patients enriched for high-risk cases using denaturing high-performance liquid chromatography (DHPLC), functional analysis of separated alleles in yeast (FASAY), and the AmpliChip p53 Research Test in parallel. The presence of T53 gene mutations was also evaluated using ultra-deep next generation sequencing (NGS) in 69 patients. In total, 79 TP53 mutations in 57 (31 %) patients were found; among them, missense substitutions predominated (68 % of detected mutations). Comparing the efficacy of the methods used, DHPLC and FASAY both combined with direct Sanger sequencing achieved the best results, identifying 95 % and 93 % of TP53-mutated patients. Nevertheless, we showed that in CLL patients carrying low-proportion TP53 mutation, the more sensitive approach, e.g., ultra-deep NGS, might be more appropriate. TP53 gene analysis using DHPLC or FASAY is a suitable approach for mutation detection. Ultra-deep NGS has the potential to overcome shortcomings of methods currently used, allows the detection of minor proportion mutations, and represents thus a promising methodology for near future.