A bumpy ride on the diagnostic bench of massive parallel sequencing, the case of the mitochondrial genome

The advent of massive parallel sequencing (MPS) has revolutionized the field of human molecular genetics, including the diagnostic study of mitochondrial (mt) DNA dysfunction. The analysis of the complete mitochondrial genome using MPS platforms is now common and will soon outrun conventional sequencing. However, the development of a robust and reliable protocol is rather challenging. A previous pilot study for the re-sequencing of human mtDNA revealed an uneven coverage, affecting predominantly part of the plus strand. In an attempt to address this problem, we undertook a comparative study of standard and modified protocols for the Ion Torrent PGM system. We could not improve strand representation by altering the recommended shearing methodology of the standard workflow or omitting the DNA polymerase amplification step from the library construction process. However, we were able to associate coverage bias of the plus strand with a specific sequence motif. Additionally, we compared coverage and variant calling across technologies. The same samples were also sequenced on a MiSeq device which showed that coverage and heteroplasmic variant calling were much improved

AI is a viable alternative to high throughput screening: a 318-target study

: High throughput screening (HTS) is routinely used to identify bioactive small molecules. This requires physical compounds, which limits coverage of accessible chemical space. Computational approaches combined with vast on-demand chemical libraries can access far greater chemical space, provided that the predictive accuracy is sufficient to identify useful molecules. Through the largest and most diverse virtual HTS campaign reported to date, comprising 318 individual projects, we demonstrate that our AtomNet® convolutional neural network successfully finds novel hits across every major therapeutic area and protein class. We address historical limitations of computational screening by demonstrating success for target proteins without known binders, high-quality X-ray crystal structures, or manual cherry-picking of compounds. We show that the molecules selected by the AtomNet® model are novel drug-like scaffolds rather than minor modifications to known bioactive compounds. Our empirical results suggest that computational methods can substantially replace HTS as the first step of small-molecule drug discovery

Late\u2010stage diversification of pyrazoles as antileishmanial agents

Abstract: N \u2010Pyrazolylcarboxamides and N \u2010pyrazolylureas represent promising lead compounds for the development of novel antileishmanial drugs. Herein, we report the late\u2010stage diversification of 3\u2010bromopyrazoles 10\u2009A / B and 14\u2009A by Pd\u2010catalyzed Sonogashira and Suzuki\u2010Miyaura cross coupling reactions. The electron\u2010withdrawing properties of the cyano moiety in 4\u2010position of the pyrazole ring limited the acylation of the primary amino moiety in 5\u2010position. A large set of pyrazoles bearing diverse aryl and alkynyl substituents in 3\u2010position was prepared and the antileishmanial and antitrypanosomal activity was recorded. The urea 38 lacking the electron withdrawing cyano moiety in 4\u2010position and containing the large 4\u2010benzylpiperidinoo moiety exhibited a modest antileishmanial ( IC 50 =19\u2005\u3bcM) and antitrypanosomal activity ( IC 50 =7.9\u2005\u3bcM)). However, its considerable toxicity against the PMM and MRC\u20105 cells indicates low selectivity, i.\u2009e . a small gap between the desired antiparasitic activity and undesired cytotoxicity of <2\u2010 to 4\u2010fold

N-pyrazolyl- and N-triazolylamines and -ureas as antileishmanial and antitrypanosomal drugs

Abstract: Three types of modifications of antileishmanial pyrazole lead compounds 7 and 8 were conducted to expand understanding of the relationships between structural features and antileishmanial/antitrypanosomal activity: (1) the pyrazole core was retained or replaced by a 1,2,4-triazole ring; (2) various aryl moieties including 2-fluorophenyl, pyridin-3-yl and pyrazin-2-yl rings were attached at 3-position of the core azole; (3) either arylmethylamino or ureido substituents were introduced at 5-position of the azole core. The synthesis followed established routes starting with esters 9 or 15 and anhydride 21. The synthesized 3-arylpyrazoles and 3-aryl-1,2,4-triazoles had only very low antileishmanial activity. The 2-fluorophenyl-substituted pyrazole 18c revealed the highest antileishmanial activity of this series of compounds, but its IC50 value (20 mu M) still indicates low activity. However, low micromolar antitrypanosomal activity was detected for the pyridin-3-yl-substituted pyrazoles 12b (IC50=4.7 mu M) and 14a (IC50=2.1 mu M). Their IC50 values are comparable with the IC50 values of the reference compounds benznidazole and nifurtimox. Whereas only low unspecific cytotoxicity at the primary peritoneal mouse macrophages (PMM) was detected, considerable cytotoxicity at MRC-5 human fibroblast cells was found for both pyrazoles 12b an 14a. The activity of pyrazole 12b against T. cruzi is 4-fold higher than its unspecific MRC-5 cytotoxicity. Within the framework of "drugs for neglected disease initiative (DNDi)", novel relationships between the structure of 3-arylpyrazoles and 3-aryl-1,2,4-triazoles with different substituents in 5-position and their antileishmanial and antitrypanosomal activity are investigated Activity against L. infantum is not detected. However, low micromolar antitrypanosomal activity comparable with the antitrypanosomal activity of reference compounds benznidazole and nifurtimox.is observed for some pyrazoles. imag

Long-term hematopoietic stem cells trigger quiescence in Leishmania parasites

Abstract: Addressing the challenges of quiescence and post-treatment relapse is of utmost importance in the microbiology field. This study shows that Leishmania infantum and L. donovani parasites rapidly enter into quiescence after an estimated 2\u20133 divisions in both human and mouse bone marrow stem cells. Interestingly, this behavior is not observed in macrophages, which are the primary host cells of the Leishmania parasite. Transcriptional comparison of the quiescent and non-quiescent metabolic states confirmed the overall decrease of gene expression as a hallmark of quiescence. Quiescent amastigotes display a reduced size and signs of a rapid evolutionary adaptation response with genetic alterations. Our study provides further evidence that this quiescent state significantly enhances resistance to treatment. Moreover, transitioning through quiescence is highly compatible with sand fly transmission and increases the potential of parasites to infect cells. Collectively, this work identified stem cells in the bone marrow as a niche where Leishmania quiescence occurs, with important implications for antiparasitic treatment and acquisition of virulence traits

Accurate and comprehensive analysis of single nucleotide variants and large deletions of the human mitochondrial genome in DNA and single cells

Massive parallel sequencing (MPS) can accurately quantify mitochondrial DNA (mtDNA) single nucleotide variants (SNVs), but no MPS methods are currently validated to simultaneously and accurately establish the breakpoints and frequency of large deletions at low heteroplasmic loads. Here we present the thorough validation of an MPS protocol to quantify the load of very low frequency, large mtDNA deletions in bulk DNA and single cells, along with SNV calling by standard methods. We used a set of well-characterized DNA samples, DNA mixes and single cells to thoroughly control the study. We developed a custom script for the detection of mtDNA rearrangements that proved to be more accurate in detecting and quantifying deletions than pre-existing tools. We also show that PCR conditions and primersets must be carefully chosen to avoid biases in the retrieved variants and an increase in background noise, and established a lower detection limit of 0.5% heteroplasmic load for large deletions, and 1.5 and 2% for SNVs, for bulk DNA and single cells, respectively. Finally, the analysis of different single cells provided novel insights into mtDNA cellular mosaicism.European Journal of Human Genetics advance online publication, 23 August 2017; doi:10.1038/ejhg.2017.129.status: publishe

Liquid Biopsy-Derived DNA Sources as Tools for Comprehensive Mutation Profiling in Multiple Myeloma: A Comparative Study

The analysis of bone marrow (BM) samples in multiple myeloma (MM) patients can lead to the underestimation of the genetic heterogeneity within the tumor. Blood-derived liquid biopsies may provide a more comprehensive approach to genetic characterization. However, no thorough comparison between the currently available circulating biomarkers as tools for mutation profiling in MM has been published yet and the use of extracellular vesicle-derived DNA for this purpose in MM has not yet been investigated. Therefore, we collected BM aspirates and blood samples in 30 patients with active MM to isolate five different DNA types, i.e., cfDNA, EV-DNA, BM-DNA and DNA isolated from peripheral blood mononucleated cells (PBMNCs-DNA) and circulating tumor cells (CTC-DNA). DNA was analyzed for genetic variants with targeted gene sequencing using a 165-gene panel. After data filtering, 87 somatic and 39 germline variants were detected among the 149 DNA samples used for sequencing. cfDNA showed the highest concordance with the mutation profile observed in BM-DNA and outperformed EV-DNA, CTC-DNA and PBMNCs-DNA. Of note, 16% of all the somatic variants were only detectable in circulating biomarkers. Based on our analysis, cfDNA is the preferable circulating biomarker for genetic characterization in MM and its combined use with BM-DNA allows for comprehensive mutation profiling in MM

Liquid Biopsy-Derived DNA Sources as Tools for Comprehensive Mutation Profiling in Multiple Myeloma: A Comparative Study

The analysis of bone marrow (BM) samples in multiple myeloma (MM) patients can lead to the underestimation of the genetic heterogeneity within the tumor. Blood-derived liquid biopsies may provide a more comprehensive approach to genetic characterization. However, no thorough comparison between the currently available circulating biomarkers as tools for mutation profiling in MM has been published yet and the use of extracellular vesicle-derived DNA for this purpose in MM has not yet been investigated. Therefore, we collected BM aspirates and blood samples in 30 patients with active MM to isolate five different DNA types, i.e. cfDNA, EV-DNA, BM-DNA and DNA isolated from peripheral blood mononucleated cells (PBMNCs-DNA) and circulating tumor cells (CTC-DNA). DNA was analyzed for genetic variants with targeted gene sequencing using a 165-gene panel. After data filtering, 87 somatic and 39 germline variants were detected among the 149 DNA samples used for sequencing. cfDNA showed the highest concordance with the mutation profile observed in BM-DNA and outperformed EV-DNA, CTC-DNA and PBMNCs-DNA. Of note, 16% of all the somatic variants were only detectable in circulating biomarkers. Based on our analysis, cfDNA is the preferable circulating biomarker for genetic characterization in MM and its combined use with BM-DNA allows for comprehensive mutation profiling in MM.SCOPUS: ar.jinfo:eu-repo/semantics/publishe

Identification of candidate cancer predisposing variants by performing whole-exome sequencing on index patients from BRCA1 and BRCA2-negative breast cancer families

Abstract Background In the majority of familial breast cancer (BC) families, the etiology of the disease remains unresolved. To identify missing BC heritability resulting from relatively rare variants (minor allele frequency ≤ 1%), we have performed whole exome sequencing followed by variant analysis in a virtual panel of 492 cancer-associated genes on BC patients from BRCA1 and BRCA2 negative families with elevated BC risk. Methods BC patients from 54 BRCA1 and BRCA2-negative families with elevated BC risk and 120 matched controls were considered for germline DNA whole exome sequencing. Rare variants identified in the exome and in a virtual panel of cancer-associated genes [492 genes associated with different types of (hereditary) cancer] were compared between BC patients and controls. Nonsense, frame-shift indels and splice-site variants (strong protein-damaging variants, called PDAVs later on) observed in BC patients within the genes of the panel, which we estimated to possess the highest probability to predispose to BC, were further validated using an alternative sequencing procedure. Results Exome- and cancer-associated gene panel-wide variant analysis show that there is no significant difference in the average number of rare variants found in BC patients compared to controls. However, the genes in the cancer-associated gene panel with nonsense variants were more than two-fold over-represented in women with BC and commonly involved in the DNA double-strand break repair process. Approximately 44% (24 of 54) of BC patients harbored 31 PDAVs, of which 11 were novel. These variants were found in genes associated with known or suspected BC predisposition (PALB2, BARD1, CHEK2, RAD51C and FANCA) or in predisposing genes linked to other cancer types but not well-studied in the context of familial BC (EXO1, RECQL4, CCNH, MUS81, TDP1, DCLRE1A, DCLRE1C, PDE11A and RINT1) and genes associated with different hereditary syndromes but not yet clearly associated with familial cancer syndromes (ABCC11, BBS10, CD96, CYP1A1, DHCR7, DNAH11, ESCO2, FLT4, HPS6, MYH8, NME8 and TTC8). Exome-wide, only a few genes appeared to be enriched for PDAVs in the familial BC patients compared to controls. Conclusions We have identified a series of novel candidate BC predisposition variants/genes. These variants/genes should be further investigated in larger cohorts/case-control studies. Other studies including co-segregation analyses in affected families, locus-specific loss of heterozygosity and functional studies should shed further light on their relevance for BC risk

RNAseq data from quiescent amastigotes recovered from mouse LT-HSC.

(a) Principal component analysis (PCA) of the RNAseq data revealing distant clustering of the independent DsRedhi and quiescent samples. (b) Euclidean distance matrix between the samples illustrating the Poisson Distance. (c) Sorted amastigotes (DsRedhi and quiescent) of infected LT-HSC were RNA extracted and subjected to RT-dPCR for 18 target genes (S2 Table). (d) GO term analysis of 167 genes that are found to be expressed in the three independent quiescent Leishmania amastigote samples. Visual representation of GO terms enriched in biological processes, molecular function and cellular components. (TIF)</p