Deep Amplicon Sequencing Research Articles

Application of a new highly multiplexed amplicon sequencing tool to evaluate Plasmodium falciparum antimalarial resistance and relatedness in individual and pooled samples from Dschang, Cameroon.

Resistance to antimalarial drugs remains a major obstacle to malaria elimination. Multiplexed, targeted amplicon sequencing is being adopted for surveilling resistance and dissecting the genetics of complex malaria infections. Moreover, genotyping of parasites and detection of molecular markers drug resistance in resource-limited regions requires open-source protocols for processing samples, using accessible reagents, and rapid methods for processing numerous samples including pooled sequencing. P lasmodium falciparum S treamlined M ultiplex A ntimalarial R esistance and R elatedness T esting ( Pf -SMARRT) is a PCR-based amplicon panel consisting of 15 amplicons targeting antimalarial resistance mutations and 9 amplicons targeting hypervariable regions. This assay uses oligonucleotide primers in two pools and a non-proprietary library and barcoding approach. We evaluated Pf -SMARRT using control mocked dried blood spots (DBS) at varying levels of parasitemia and a mixture of 3D7 and Dd2 strains at known frequencies, showing the ability to genotype at low parasite density and recall within-sample allele frequencies. We then piloted Pf -SMARRT to genotype 100 parasite isolates collected from uncomplicated malaria cases at three health facilities in Dschang, Western Cameroon. Antimalarial resistance genotyping showed high levels of sulfadoxine-pyrimethamine resistance mutations, including 31% prevalence of the DHPS A613S mutation. No K13 candidate or validated artemisinin partial resistance mutations were detected, but one low-level non-synonymous change was observed. Pf -SMARRT's hypervariable targets, used to assess complexity of infections and parasite diversity and relatedness, showed similar levels and patterns compared to molecular inversion probe (MIP) sequencing. While there was strong concordance of antimalarial resistance mutations between individual samples and pools, low-frequency variants in the pooled samples were often missed. Overall, Pf -SMARRT is a robust tool for assessing parasite relatedness and antimalarial drug resistance markers from both individual and pooled samples. Control samples support that accurate genotyping as low as 1 parasite per microliter is routinely possible. Malaria remains a critical global public health problem. Antimalarial drug resistance has repeatedly undermined control and the emergence of artemisinin partial resistance in Africa is the latest major challenge. Malaria molecular surveillance (MMS) has emerged as a powerful tool to monitor molecular markers of resistance and changes in the parasite population. Streamlined methods are needed that can be readily adopted in endemic countries. We developed P lasmodium falciparum S treamlined M ultiplex A ntimalarial R esistance and R elatedness T esting ( Pf -SMARRT), a multiplex amplicon deep sequencing approach that uses easily accessible products without proprietary steps and can be sequenced on any Illumina sequencer. We validated this tool using controls, including mocked dried blood spots, and then implemented it to evaluate resistance and parasite relatedness among 100 samples from Cameroon. The assay was able to reliably assess the within-sample allele frequency of antimalarial resistance markers and discriminate strains within and between individuals. We also evaluated a more cost-effective surveillance approach for antimalarial resistance polymorphisms using pooled samples. While within-pool frequencies of mutations were accurate in pools with higher numbers of samples, this resulted in the loss of the ability to detect variants uncommon in the pool. Overall Pf- SMARRT provides a new protocol for conducting MMS that is easily implementable in Africa.

Genotyping methods to distinguish Plasmodium falciparum recrudescence from new infection for the assessment of antimalarial drug efficacy: an observational, single-centre, comparison study

Distinguishing Plasmodium falciparum recrudescence from new infections is crucial for the assessment of antimalarial drug efficacy against Pfalciparum. We aimed to compare the efficacy of different genotyping methods to assess their effect on drug efficacy estimates, particularly in patients from high-transmission settings with polyclonal infections. In this head-to-head comparison study, we compared five different genotyping methods currently used: fast capillary electrophoresis (F-CE) using msp1, msp2, and glurp; high-resolution capillary electrophoresis (H-CE) using msp1, msp2, and glurp; H-CE using microsatellites; targeted amplicon deep sequencing (TADS) using single nucleotide polymorphism (SNP)-rich markers; and high-resolution melting (HRM) analysis using msp1 and msp2. We assessed their sensitivity in detecting minority clones in polyclonal infections, their reproducibility, and the genetic diversity of the markers used. Our study used four well characterised Pfalciparum laboratory strains mixed in varying ratios, and 20paired samples collected from an in-vivo clinical trial. The experiments were performed at the Swiss Tropical and Public Health Institute in Basel, Switzerland between May 5, 2020, and Aug 23,2021. H-CE using msp1 and msp2 and TADS revealed the highest sensitivity in detecting minority clones (up to ratios of 1:100 for H-CE and 50:1:1:1 for TADS in the FCB1:HB3 and 3D7:K1:HB3:FCB1 laboratory strain mixtures, respectively), highest reproducibility (intra-assay: 99% and 91% for H-CE and TADS, respectively; inter-assay: 98% and 92% for H-CE and TADS, respectively), and highest genetic diversity in the used markers (up to 36and32unique genotypes in 20paired samples for H-CE using msp2 and TADS using cpmp, respectively). Microsatellites assessed by H-CE had a lower genetic diversity compared with msp1, msp2, and glurp assessed by H-CE and the SNP-rich markers assessed by TADS, with a maximum of 13unique genotypes, and some genotypes having allelic frequencies larger than 30%.Markers used by TADS gave the most consistent results in distinguishing recrudescence from new infection across all methods (in 18of 20pairs of samples vs 15of 20pairs for H-CE). WHO currently recommends replacing glurp with microsatellites. However, in this study, the replacement of glurp with microsatellites did not change the genotyping outcome, probably due to the lower genetic diversity of microsatellites. More studies with large sample sizes are required to identify the most suitablemicrosatellites that could replace glurp. Our study indicates that TADS should be considered the gold standardfor genotyping to distinguish recrudescence from new infection, and that it should be used to validate other methods. Swiss Tropical and Public Health Institute.

Application of multiple mosaic callers improves post-zygotic mutation detection from exome sequencing data

PurposeThe gold standard for identification of post-zygotic variants (PZVs) is droplet digital polymerase chain reaction or high-depth sequencing across multiple tissues types. These approaches are yet to be systematically implemented for monogenic disorders. We developed PZV detection pipelines for correct classification of de novo variants. MethodOur pipelines detect PZV in parents (gonosomal mosaicism [pGoM]) and children (somatic mosaicism, “M3”). We applied them to research exome sequencing (ES) data from the Australian Cerebral Palsy Biobank (n = 145 trios) and Simons Simplex Collection (n = 405 families). Candidate mosaic variants were validated using deep amplicon sequencing or droplet digital polymerase chain reaction. Results69.2% (M3trio), 63.9% (M3single), and 92.7% (pGoM) of detected variants were validated, with 48.6%, 56.7%, and 26.2% of variants, respectively, meeting strict criteria for mosaicism. In the Australian Cerebral Palsy Biobank, 16.6% of probands and 20.7% of parents had at least 1 true-positive somatic or pGoM variant, respectively. A large proportion of PZVs detected in Simons Simplex Collection parents (79.8%) and child (94.5%) were not previously reported. We reclassified 3.7% to 8.0% of germline de novo variants as mosaic. ConclusionMany PZVs were incorrectly classified as germline variants or missed by previous approaches. Systematic application of our pipelines could increase genetic diagnostic rate, improve estimates of recurrence risk in families, and benefit novel disease gene identification.

Allele-specific CRISPR-Cas9 editing inactivates a single nucleotide variant associated with collagen VI muscular dystrophy

The application of allele-specific gene editing tools can expand the therapeutic options for dominant genetic conditions, either via gene correction or via allelic gene inactivation in situations where haploinsufficiency is tolerated. Here, we used allele-targeted CRISPR/Cas9 guide RNAs (gRNAs) to introduce inactivating frameshifting indels at a single nucleotide variant in the COL6A1 gene (c.868G>A; G290R), a variant that acts as dominant negative and that is associated with a severe form of congenital muscular dystrophy. We expressed spCas9 along with allele-targeted gRNAs, without providing a repair template, in primary fibroblasts derived from four patients and one control subject. Amplicon deep-sequencing for two gRNAs tested showed that single nucleotide deletions accounted for the majority of indels introduced. While activity of the two gRNAs was greater at the G290R allele, both gRNAs were also active at the wild-type allele. To enhance allele-selectivity, we introduced deliberate additional mismatches to one gRNA. One of these optimized gRNAs showed minimal activity at the WT allele, while generating productive edits and improving collagen VI matrix in cultured patient fibroblasts. This study strengthens the potential of gene editing to treat dominant-negative disorders, but also underscores the challenges in achieving allele selectivity with gRNAs.

Origins, diversity, and adaptive evolution of DWV in the honey bees of the Azores: the impact of the invasive mite Varroa destructor.

Deformed wing virus (DWV) is a honey bee virus, whose emergence from relative obscurity is driven by the recent host-switch, adaptation, and global dispersal of the ectoparasitic mite Varroa destructor (a highly efficient vector of DWV) to reproduction on honey bees (Apis mellifera). Our study examines how varroa affects the continuing evolution of DWV, using the Azores archipelago, where varroa is present on only three out of the eight Islands, as a natural experimental system for comparing different evolutionary conditions and trajectories. We combined qPCR of 494 honey bee colonies sampled across the archipelago with amplicon deep sequencing to reveal how the DWV genetic landscape is altered by varroa. Two of the varroa-free Islands were also free of DWV, while a further two Islands were intriguingly dominated by the rare DWV-C major variant. The other four Islands, including the three varroa-infested Islands, were dominated by the common DWV-A and DWV-B variants. The varroa-infested Islands had, as expected, an elevated DWV prevalence relative to the uninfested Islands, but not elevated DWV loads, due the relatively high prevalence and loads of DWV-C on the varroa-free Islands. This establishes the Azores as a stable refuge for DWV-C and provides the most convincing evidence to date that at least some major strains of DWV may be capable of not just surviving, but actually thriving in honey bees in the absence of varroa-mediated transmission. We did not detect any change in DWV genetic diversity associated with island varroa status but did find a positive association of DWV diversity with virus load, irrespective of island varroa status.

Bidirectional transfer of human cytomegalovirus strains in donor and recipient seropositive lung transplant patients

AbstractDonor and recipient human cytomegalovirus (HCMV) seropositive (D+R+) lung transplant recipients (LTRs) often harbor multiple strains of HCMV, likely due to transmitted donor (D) strains and reactivated recipient (R) strains. To date, the extent and timely occurrence of each likely source in shaping the post‐transplantation (post‐Tx) strain population is unknown. Here, we deciphered the D and R origin of the post‐Tx HCMV strain composition in blood, bronchoalveolar lavage (BAL), and CD45+ BAL cell subsets. We investigated either D and/or R formalin‐fixed paraffin‐embedded blocks or fresh D lung tissue from four D+R+ LTRs obtained before transplantation. HCMV strains were characterized by short amplicon deep sequencing. In two LTRs, we show that the transplanted lung is reseeded by R strains within the first 6 months after transplantation, likely by infiltrating CD14+ CD163+/− alveolar macrophages. In three LTRs, we demonstrate both rapid D‐strain dissemination and persistence in the transplanted lung for >1 year post‐Tx. Broad inter‐host diversity contrasts with intra‐host genotype sequence stability upon transmission, during follow‐up and across compartments. In D+R+ LTRs, HCMV strains of both, D and R origin can emerge first and dominate long‐term in subsequent episodes of infection, indicating replication of both sources despite pre‐existing immunity.

Deep Amplicon Sequencing Reveals Culture Selection of Mycobacterium Tuberculosis by Clinical Samples.

The commonly-used drug susceptibility testing (DST) relies on bacterial culture and faces shortcomings such as long turnaround time and clone/subclone selection. We developed a targeted deep amplification sequencing (DAS) method directly applied to clinical specimens. In this DAS panel, we examined 941 drug-resistant mutations associated with 20 anti-tuberculosis drugs with an initial amount of 4 pg DNA and reduced clinical testing time from 20 days to two days. A prospective study was conducted using 115 clinical specimens mainly with Xpert® Mycobacterium tuberculosis/rifampicin (Xpert MTB/RIF) assay positive to evaluate drug-resistant mutation detection. DAS was performed on culture-free specimens, while culture-dependent isolates were used for phenotypic DST, DAS, and whole-genome sequencing (WGS). For in silico molecular DST, our result based on DAS panel revealed the similar accuracy to three published reports based on WGS. For 82 isolates, application of DAS showed better sensitivity (93.03% vs. 92.16%), specificity (96.10% vs. 95.02%), and accuracy (91.33% vs. 90.62%) than Mykrobe software using WGS. Compared to culture-dependent WGS, culture-free DAS provides a full picture of sequence variation at population level, exhibiting in detail the gain-and-loss variants caused by bacterial culture. Our study performs a systematic verification of the advantages of DAS in clinical applications and comprehensively illustrates the discrepancy in Mycobacterium tuberculosis before and after culture.

Use of Deep-Amplicon Sequencing (DAS), Real-Time PCR and In Situ Hybridization to Detect H. pylori and Other Pathogenic Helicobacter Species in Feces from Children.

Detecting Helicobacter pylori in fecal samples is easier and more comfortable than invasive techniques, especially in children. Thus, the objective of the present work was to detect H. pylori in feces from children by molecular methods as an alternative for diagnostic and epidemiological studies. Forty-five fecal samples were taken from pediatric patients who presented symptoms compatible with H. pylori infection. HpSA test, culture, real-time quantitative PCR (qPCR), fluorescence in situ hybridization (FISH), direct viable count associated with FISH (DVC-FISH), and Illumina-based deep-amplicon sequencing (DAS) were applied. No H. pylori colonies were isolated from the samples. qPCR analysis detected H. pylori in the feces of 24.4% of the patients. In comparison, DVC-FISH analysis showed the presence of viable H. pylori cells in 53.3% of the samples, 37% of which carried 23S rRNA mutations that confer resistance to clarithromycin. After DAS, H. pylori-specific 16S rDNA sequences were detected in 26 samples. In addition, DNA from H. hepaticus was identified in 10 samples, and H. pullorum DNA was detected in one sample. The results of this study show the presence of H. pylori, H. hepaticus, and H. pullorum in children's stools, demonstrating the coexistence of more than one Helicobacter species in the same patient. The DVC-FISH method showed the presence of viable, potentially infective H. pylori cells in a high percentage of the children's stools. These results support the idea that fecal-oral transmission is probably a common route for H. pylori and suggest possible fecal-oral transmission of other pathogenic Helicobacter species.

Worldwide absence of canonical benzimidazole resistance-associated mutations within β-tubulin genes from Ascaris.

The giant roundworm Ascaris is an intestinal nematode, causing ascariasis by infecting humans and pigs worldwide. Recent estimates suggest that Ascaris infects over half a billion people, with chronic infections leading to reduced growth and cognitive ability. Ascariasis affects innumerable pigs worldwide and is known to reduce production yields via decreased growth and condemnation of livers. The predominant anthelminthic drugs used to treat ascariasis are the benzimidazoles. Benzimidazoles interact with β-tubulins and block their function, and several benzimidazole resistance-associated mutations have been described in the β-tubulins of ruminant nematodes. Recent research on ascarids has shown that these canonical benzimidazole resistance-associated mutations are likely not present in the β-tubulins of Ascaris, Ascaridia or Parascaris, even in phenotypically resistant populations. To further determine the putative absence of key β-tubulin polymorphisms, we screened two β-tubulin isotypes of Ascaris, highly expressed in adult worms. Using adult and egg samples of Ascaris obtained from pigs and humans worldwide, we performed deep amplicon sequencing to look for canonical resistance-associated mutations in Ascaris β-tubulins. Subsequently, we examined these data in closer detail to study the population dynamics of Ascaris and genetic diversity within the two isotypes and tested whether genotypes appeared to partition across human and pig hosts. In the 187 isolates, 69 genotypes were found, made up of eight haplotypes of β-tubulin isotype A and 20 haplotypes of isotype B. Single nucleotide polymorphisms were seen at 14 and 37 positions for β-tubulin isotype A and isotype B, respectively. No evidence of any canonical benzimidazole resistance-associated mutations was found in either human- or pig-derived Ascaris isolates. There was, however, a difference in the genetic diversity of each isotype and distribution of β-tubulin genotypes between human- and pig-derived Ascaris. Statistical tests of population differentiation show significant differences (p < 0.001) between pig- and human-derived worms; however, more diversity was seen between worms from different populations than worms from different hosts. Our work suggests an absence of canonical β-tubulin mutations within Ascaris, but alternative modes of anthelminthic resistance may emerge necessitating continued genetic scrutiny alongside monitoring of drug efficacy.

A Targeted Deep Sequencing Method to Quantify Endogenous Retrovirus Gag Sequence Variants and Open Reading Frames Expressed in Nonobese Diabetic Mice.

Endogenous retroviruses (ERVs) are involved in autoimmune diseases such as type 1 diabetes (T1D). ERV gene products homologous to murine leukemia retroviruses are expressed in the pancreatic islets of NOD mice, a model of T1D. One ERV gene, Gag, with partial or complete open reading frames (ORFs), is detected in the islets, and it contains many sequence variants. An amplicon deep sequencing analysis was established by targeting a conserved region within the Gag gene to compare NOD with T1D-resistant mice or different ages of prediabetic NOD mice. We observed that the numbers of different Gag variants and ORFs are linked to T1D susceptibility. More importantly, these numbers change during the course of diabetes development and can be quantified to calculate the levels of disease progression. Sequence alignment analysis led to identification of additional markers, including nucleotide mismatching and amino acid consensus at specific positions that can distinguish the early and late stages, before diabetes onset. Therefore, the expression of sequence variants and ORFs of ERV genes, particularly Gag, can be quantified as biomarkers to estimate T1D susceptibility and disease progression.

Validation of deep amplicon sequencing of Dicrocoelium in small ruminants from Northern regions of Pakistan.

Dicrocoelium lancet flukes cause significant production loss in ruminant livestock. Although co-infection with multiple Dicrocoelium species within a host is common, techniques for studying the composition of these complex parasite communities are lacking. The pathogenicity, epidemiology, and therapeutic susceptibility of different helminth species vary, and little is known about the interactions that take place between co-infecting species and their hosts. Here, we describe the first applicationof metabarcoding deep amplicon sequencing method to studythe Dicrocoelium species in sheep and goats. First, rDNA ITS-2 sequences of four Dicrocoelium species (Dicrocoelium dendriticum, Dicrocoelium hospes, Dicrocoelium orientalis, and Dicrocoelium chinensis) were extracted from the NCBI public database. Phylogenetic analysis revealed separate clades of Dicrocoelium species; hence, molecular differentiation between each species is possible in co-infections. Second, 202 flukes belonging to seventeen host populations (morphologically verified as belonging to the Dicrocoelium genus) were evaluated to determine the deep amplicon sequencing read threshold of an individual fluke for each of the four species. The accuracy of the method in proportional quantification of samples collected from single hosts was further assessed. Overall, 198 (98.01%) flukes were confirmed as D. dendriticum and 1.98% produced no reads. The comparison of genetic distances between rDNA ITS-2 revealed 86% to 98% identity between the Dicrocoelium species. Phylogenetic analysis demonstrated a distinct clustering of species, apart from D. orientalis and D. chinensis, which sit very close to each other in a single large clade whereas D. hospes and D. dendriticum are separated into their own clade. In conclusion each sample was identified as D. dendriticum based on the proportion of MiSeq reads and validated the presence of this group of parasites in the Gilgit Baltistan and Khyber Pakhtunkhwa provinces of Pakistan. The metabarcoding deep amplicon sequencing technology and bioinformatics pathway have several potential applications, including species interactions during co-infections, identifying the host and geographical distribution of Dicrocoelium in livestock, drug therapy response evaluation and understanding of the emergence and spread of drug resistance.

Genotyping Plasmodium falciparum gametocytes using amplicon deep sequencing

BackgroundUnderstanding the dynamics of gametocyte production in polyclonal Plasmodium falciparum infections requires a genotyping method that detects distinct gametocyte clones and estimates their relative frequencies. Here, a marker was identified and evaluated to genotype P. falciparum mature gametocytes using amplicon deep sequencing.MethodsA data set of polymorphic regions of the P. falciparum genome was mined to identify a gametocyte genotyping marker. To assess marker resolution, the number of unique haplotypes in the marker region was estimated from 95 Malawian P. falciparum whole genome sequences. Specificity of the marker for detection of mature gametocytes was evaluated using reverse transcription-polymerase chain reaction of RNA extracted from NF54 mature gametocytes and rings from a non-gametocyte-producing strain of P. falciparum. Amplicon deep sequencing was performed on experimental mixtures of mature gametocytes from two distinct parasite clones, as well as gametocyte-positive P. falciparum field isolates to evaluate the quantitative ability and determine the limit of detection of the genotyping approach.ResultsA 400 bp region of the pfs230 gene was identified as a gametocyte genotyping marker. A larger number of unique haplotypes was observed at the pfs230 marker (34) compared to the sera-2 (18) and ama-1 (14) markers in field isolates from Malawi. RNA and DNA genotyping accurately estimated gametocyte and total parasite clone frequencies when evaluating agreement between expected and observed haplotype frequencies in gametocyte mixtures, with concordance correlation coefficients of 0.97 [95% CI: 0.92–0.99] and 0.92 [95% CI: 0.83–0.97], respectively. The detection limit of the genotyping method for male gametocytes was 0.41 pfmget transcripts/µl [95% CI: 0.28–0.72] and for female gametocytes was 1.98 ccp4 transcripts/µl [95% CI: 1.35–3.68].ConclusionsA region of the pfs230 gene was identified as a marker to genotype P. falciparum gametocytes. Amplicon deep sequencing of this marker can be used to estimate the number and relative frequency of parasite clones among mature gametocytes within P. falciparum infections. This gametocyte genotyping marker will be an important tool for studies aimed at understanding dynamics of gametocyte production in polyclonal P. falciparum infections.

Synergistic approach of PCR-based fragment length analysis and amplicon deep sequencing reveals rich diversity of S-alleles in sweet cherries from the Caucasian region of origin.

The self-incompatibility system in sweet cherry (Prunus avium L.) prevents fertilization with own or genetically related pollen, and is genetically determined by the multi-allelic S-locus. Therefore, determining S-alleles is crucial for plant breeding and fruit production, as it enables the selection of compatible combinations of S-genotypes for successful pollination. In this study, S-alleles were identified in a total of 260 genotypes from the Caucasian region, the species' center of origin. S-allele genotyping was conducted using PCR fragment length analysis with the standard marker PaConsI-F/R2 and reference genotypes, complemented by sequence analysis through amplicon deep sequencing. The genotypes collected from Azerbaijan and Turkey exhibit a high allelic richness at the S-locus, particularly compared to modern sweet cherry cultivars worldwide. Nine previously undescribed S-alleles were identified and designated as S45, S46, S47, S48, S49, S50, S51, S52 and S53. Given the expected high diversity for other traits, this plant material represents a valuable resource for further breeding research and introgression of new traits in future breeding programs. Furthermore, our results underscore that fragment length alone may not be sufficient for unambiguous assignment of S-alleles due to minimal length differences between different alleles. To address this issue, an S-allele reference ladder was developed using the rich diversity for precise assignment of the S-alleles. This tool can be applied in future experiments as a robust and cost-effective method for accurate S-genotyping across different runs and laboratories. Additionally, several selected S-genotypes were planted in a trial field and will be maintained as an S-allele reference collection.

Abstract 987: Novel ALK mutations in EML4::ALK+ NSCLC resistant to TKIs identified by liquid biopsy

Abstract The identification of Anaplastic Lymphoma Kinase (ALK) chromosomal translocations led to the development of effective targeted therapies in a subset of Non-Small Cell Lung Cancer (NSCLC) patients. Several drugs are currently in clinical practice for advanced ALK+ NSCLC. However, resistance remains a challenge in these patients. The possibility to analyze circulating tumor DNA (ctDNA) with non-invasive methods can greatly improve prognostic capabilities, by allowing early detection of relapse and of drug-resistant mutations from peripheral blood. We report a series of 15 consecutive relapsed NSCLC patients (12 ALK+, 3 ROS1+) treated at the Fondazione IRCCS-San Gerardo dei Tintori (Monza, Italy) that were investigated by liquid biopsy at tyrosine kinase inhibitor (TKI) failure, through amplicon deep sequencing of the ALK/ROS1 kinase domains. Median DNA yield was 4.8 ng per ml of plasma and mean fragment size was 176 bp. Mutational hotspot coding regions of ALK and ROS1 genes were amplified by high-fidelity PCR and sequenced. Mean coverage of called variants was 24,687x. The threshold for mutation calling was set at 0.6%, after correcting for mapping quality, minimum coverage, and strand bias. Analysis of background signal on healthy control DNA identified 0.0024 mismatches per mapped base. Overall, plasma ctDNA analysis detected an ALK mutation in 7/18 relapse samples from 7/15 patients. No ROS1 mutations were found in ROS1 fusion-positive patients. Among 12 patients carrying an EML4::ALK fusion, 10 variants potentially driving resistance were found. Along with mutations that were previously identified (L1196M/G1202R, E1154K, F1174L), novel ALK variants were detected such as L1198R, T1211I, C1235R, C1237Y, L1196P. Mutations were characterized for sensitivity to ALK inhibitors in a BaF3-EML4::ALK cell model, using cell proliferation assays and western blotting to evaluate the effects on ALK phosphorylation. C1237Y and L1196P mutants appeared to resist all TKIs; hence, we further tested their sensitivity to novel investigational drugs, zotizalkib and repotrectinib. Both mutants demonstrated resistance, suggesting that these two mutations confer broad resistance to TKIs. Molecular modelling of C1237Y and L1196P mutants was run to elucidate the mechanisms of resistance. The ADP-C1237Y complex showed a much smaller root-mean-square deviation than ADP-WT, indicating a higher affinity/activity of the mutant. L1196 interacts hydrophobically with a methyl group of lorlatinib, P1196 is predicted to lose such interaction, while ATP binding is not affected. The presence of C1237Y and L1196P variants, is a relevant finding due to their resistance, even after the exposure to the fourth-generation inhibitor zotizalkib. The clonal selection of these mutations, made by second/third generation ALK TKIs may lead to a complete refractory status, which would then render TKIs useless in a sequential therapeutic strategy. Citation Format: Federica Malighetti, Matteo Villa, Elisa Sala, Geeta Sharma, Giulia Arosio, Maria Gemelli, Chiara Manfroni, Diletta Fontana, Nicoletta Cordani, Raffaella Meneveri, Alfonso Zambon, Rocco Piazza, Fabio Pagni, Diego Cortinovis, Luca Mologni. Novel ALK mutations in EML4::ALK+ NSCLC resistant to TKIs identified by liquid biopsy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 987.

Targeted amplicon deep sequencing of ama1 and mdr1 to track within-host P. falciparum diversity throughout treatment in a clinical drug trial

Introduction Antimalarial therapeutic efficacy studies are routinely conducted in malaria-endemic countries to assess the effectiveness of antimalarial treatment strategies. Targeted amplicon sequencing (AmpSeq) uniquely identifies and quantifies genetically distinct parasites within an infection. In this study, AmpSeq of Plasmodium falciparum apical membrane antigen 1 (ama1), and multidrug resistance gene 1 (mdr1), were used to characterise the complexity of infection (COI) and drug-resistance genotypes, respectively. Methods P. falciparum-positive samples were obtained from a triple artemisinin combination therapy clinical trial conducted in 30 children under 13 years of age between 2018 and 2019 in Kilifi, Kenya. Nine of the 30 participants presented with recurrent parasitemia from day 26 (624h) onwards. The ama1 and mdr1 genes were amplified and sequenced, while msp1, msp2 and glurp data were obtained from the original clinical study. Results The COI was comparable between ama1 and msp1, msp2 and glurp; overall, ama1 detected more microhaplotypes. Based on ama1, a stable number of microhaplotypes were detected throughout treatment until day 3. Additionally, a recrudescent infection was identified with an ama1 microhaplotype initially observed at 30h and later in an unscheduled follow-up visit. Using the relative frequencies of ama1 microhaplotypes and parasitemia, we identified a fast (&lt;1h) and slow (&gt;5h) clearing microhaplotype. As expected, only two mdr1 microhaplotypes (NF and NY) were identified based on the combination of amino acid polymorphisms at codons 86 and 184. Conclusions This study highlights AmpSeq as a tool for highly-resolution tracking of parasite microhaplotypes throughout treatment and can detect variation in microhaplotype clearance estimates. AmpSeq can also identify slow-clearing microhaplotypes, a potential early sign of selection during treatment. Consequently, AmpSeq has the capability of improving the discriminatory power to distinguish recrudescences from reinfections accurately.

Identification and molecular characterization of Mycobacterium bovis DNA in GeneXpert® MTB/RIF ultra-positive, culture-negative sputum from a rural community in South Africa

This study investigated the presence of Mycobacterium bovis (M. bovis) DNA in archived human sputum samples previously collected from residents who reside adjacent to the M. bovis-endemic Hluhluwe-iMfolozi wildlife park, South Africa (SA). Sixty-eight sputum samples were GeneXpert MTB/RIF Ultra-positive for M. tuberculosis complex (MTBC) DNA but culture negative for M. tuberculosis. Amplification and Sanger sequencing of hsp65 and rpoB genes from DNA extracted from stored heat-inactivated sputum samples confirmed the presence of detectable amounts of MTBC from 20 out of the 68 sputum samples. Region of difference PCR, spoligotyping and gyrB long-read amplicon deep sequencing identified M. bovis (n = 10) and M. tuberculosis (n = 7). Notably, M. bovis spoligotypes SB0130 and SB1474 were identified in 4 samples, with SB0130 previously identified in local cattle and wildlife and SB1474 exclusively in African buffaloes in the adjacent park. M. bovis DNA in sputum, from people living near the park, underscores zoonotic transmission potential in SA. Identification of spoligotypes specifically associated with wildlife only and spoligotypes found in livestock as well as wildlife, highlights the complexity of TB epidemiology at wildlife-livestock-human interfaces. These findings support the need for integrated surveillance and control strategies to curb potential spillover and for the consideration of human M. bovis infection in SA patients with positive Ultra results.

Metabarcoding study to reveal the structural community of strongylid nematodes in domesticated horses in Thailand

BackgroundMixed strongylid infections significantly impact equine health and performance. Traditional microscopy-based methods exhibit limitations in accurately identifying strongylid species. Nemabiome deep amplicon sequencing approach previously succeeded in describing the strongylid communities in livestock including equids. However, there are no available studies that describe the structural communities of strongylid parasites in horses in Thailand. Therefore, this study was undertaken encompassing the ITS-2 rDNA metabarcoding assay to characterize strongylid species within horse fecal samples collected from a cohort of yearlings at the largest domesticated stud farm in Thailand. In addition, to investigate the capability of ITS-2 rDNA in assessing the phylogenetic relationships among the identified strongylid species.ResultsThe study identified 14 strongylid species in the examined equine populations, each with varying prevalence. Notably, Cylicocyclus nassatus and Cylicostephanus longibursatus were identified as the predominant species, with Strongylus spp. conspicuously absent. The phylogenetic analysis of 207 amplicon sequence variants (ASVs) displayed a complex relationship among the investigated cyathostomin species, with some species are positioned across multiple clades, demonstrating close associations with various species and genera.ConclusionThe ITS-2 nemabiome sequencing technique provided a detailed picture of horse strongylid parasite species in the studied population. This establishes a foundation for future investigations into the resistance status of these parasites and enables efforts to mitigate their impact.

Natural History of Germline BRCA1 Mutated and BRCA Wild-type Triple-negative Breast Cancer.

In germline BRCA1 mutated and BRCA wild-type patients, TNBC shows a branching evolutionary pattern of mutations with a single founding clone, are polyclonal throughout their disease course, and have widespread copy-number aberrations. This evolutionary pattern may be associated with treatment resistance or sensitivity and could be therapeutically exploited.

Lizard host abundances and climatic factors explain phylogenetic diversity and prevalence of blood parasites on anoceanic island.

Host abundance might favour the maintenance of a high phylogenetic diversity of some parasites via rapid transmission rates. Blood parasites of insular lizards represent a good model to test this hypothesis because these parasites can be particularly prevalent in islands and host lizards highly abundant. We applied deep amplicon sequencing and analysed environmental predictors of blood parasite prevalence and phylogenetic diversity in the endemic lizard Gallotia galloti across 24 localities on Tenerife, an island in the Canary archipelago that has experienced increasing warming and drought in recent years. Parasite prevalence assessed by microscopy was over 94%, and a higher proportion of infected lizards was found in warmer and drier locations. A total of 33 different 18s rRNA parasite haplotypes were identified, and the phylogenetic analyses indicated that they belong to two genera of Adeleorina (Apicomplexa: Coccidia), with Karyolysus as the dominant genus. The most important predictor of between-locality variation in parasite phylogenetic diversity was the abundance of lizard hosts. We conclude that a combination of climatic and host demographic factors associated with an insular syndrome may be favouring a rapid transmission of blood parasites among lizards on Tenerife, which may favour the maintenance of a high phylogenetic diversity of parasites.

Analytic optimization of Plasmodium falciparum marker gene haplotype recovery from amplicon deep sequencing of complex mixtures.

Molecular epidemiologic studies of malaria parasites and other pathogens commonly employ amplicon deep sequencing (AmpSeq) of marker genes derived from dried blood spots (DBS) to answer public health questions related to topics such as transmission and drug resistance. As these methods are increasingly employed to inform direct public health action, it is important to rigorously evaluate the risk of false positive and false negative haplotypes derived from clinically-relevant sample types. We performed a control experiment evaluating haplotype recovery from AmpSeq of 5 marker genes (ama1, csp, msp7, sera2, and trap) from DBS containing mixtures of DNA from 1 to 10 known P. falciparum reference strains across 3 parasite densities in triplicate (n = 270 samples). While false positive haplotypes were present across all parasite densities and mixtures, we optimized censoring criteria to remove 83% (148/179) of false positives while removing only 8% (67/859) of true positives. Post-censoring, the median pairwise Jaccard distance between replicates was 0.83. We failed to recover 35% (477/1365) of haplotypes expected to be present in the sample. Haplotypes were more likely to be missed in low-density samples with <1.5 genomes/μL (OR: 3.88, CI: 1.82-8.27, vs. high-density samples with ≥75 genomes/μL) and in samples with lower read depth (OR per 10,000 reads: 0.61, CI: 0.54-0.69). Furthermore, minority haplotypes within a sample were more likely to be missed than dominant haplotypes (OR per 0.01 increase in proportion: 0.96, CI: 0.96-0.97). Finally, in clinical samples the percent concordance across markers for multiplicity of infection ranged from 40%-80%. Taken together, our observations indicate that, with sufficient read depth, the majority of haplotypes can be successfully recovered from DBS while limiting the false positive rate.