Multiplex Amplicon Sequencing Research Articles

Wastewater multiplex PCR amplicon sequencing revealed community transmission of SARS-CoV-2 lineages during the outbreak of infection in Chinese Mainland

During the COVID-19, wastewater-based epidemiology (WBE) has become a powerful epidemic surveillance tool widely used worldwide. However, the development and application of this technology in Chinese Mainland are relatively lagging. Herein, we for the first time monitored the community circulation of SARS-CoV-2 lineages using WBE methods in Chinese Mainland. During the peak period of infection outbreak at the end of 2022, six precious sewage samples were collected from the manhole in the student dormitory area on Wangjiang Campus of Sichuan University. RT-qPCR revealed that the six sewage samples were all positive for SARS-CoV-2 RNA. Multiplex PCR amplicon sequencing of the sewage samples reflected the local transmission of SARS-CoV-2 variants. The results of two deconvolution methods indicate that the main virus lineages have clear evolutionary genetic correlations. Furthermore, the sampling time is consistent with the timeline of concern for these virus lineages, as well as the timeline of uploading the nucleic acid sequences from the corresponding lineages in Sichuan to the database. These results demonstrate the reliability of the sewage sequencing results. Multiplex PCR amplicon sequencing is by far the most powerful analytical tool of WBE, enabling quantitative detection of virus lineages transmission and evolution at the community level.

Olivar: towards automated variant aware primer design for multiplex tiled amplicon sequencing of pathogens

Tiled amplicon sequencing has served as an essential tool for tracking the spread and evolution of pathogens. Over 15 million complete SARS-CoV-2 genomes are now publicly available, most sequenced and assembled via tiled amplicon sequencing. While computational tools for tiled amplicon design exist, they require downstream manual optimization both computationally and experimentally, which is slow and costly. Here we present Olivar, a first step towards a fully automated, variant-aware design of tiled amplicons for pathogen genomes. Olivar converts each nucleotide of the target genome into a numeric risk score, capturing undesired sequence features that should be avoided. In a direct comparison with PrimalScheme, we show that Olivar has fewer mismatches overlapping with primers and predicted PCR byproducts. We also compare Olivar head-to-head with ARTIC v4.1, the most widely used primer set for SARS-CoV-2 sequencing, and show Olivar yields similar read mapping rates (~90%) and better coverage to the manually designed ARTIC v4.1 amplicons. We also evaluate Olivar on real wastewater samples and found that Olivar has up to 3-fold higher mapping rates while retaining similar coverage. In summary, Olivar automates and accelerates the generation of tiled amplicons, even in situations of high mutation frequency and/or density. Olivar is available online as a web application at https://olivar.rice.edu and can be installed locally as a command line tool with Bioconda. Source code, installation guide, and usage are available at https://github.com/treangenlab/Olivar.

QTL discovery for agronomic and quality traits in diploid potato clones using PotatoMASH amplicon sequencing.

We genotyped a population of 618 diploid potato clones derived from six independent potato-breeding programmes from NW-Europe. The diploids were phenotyped for 23 traits, using standardised protocols and common check varieties, enabling us to derive whole population estimators for most traits. We subsequently performed a Genome-Wide Association Study (GWAS) to identify quantitative trait loci (QTL) for all traits with SNPs and short-read haplotypes derived from read-backed phasing. In this study, we used a marker platform called PotatoMASH (Potato Multi-Allele Scanning Haplotags); a pooled multiplex amplicon sequencing based approach. Through this method, neighbouring SNPs within an amplicon can be combined to generate multi-allelic short-read haplotypes (haplotags) that capture recombination history between the constituent SNPs, and reflect the allelic diversity of a given locus in a different way than single bi-allelic SNPs. We found a total of 37 unique QTL across both marker types. A core of 10 QTL were detected with SNPs as well as with haplotags. Haplotags allowed to detect an additional 14 QTL not found based on the SNP set. Conversely, the bi-allelic SNP set also found 13 QTL not detectable using the haplotag set. We conclude that both marker types should routinely be used in parallel to maximize the QTL detection power. We report 19 novel QTL for nine traits: Skin Smoothness, Sprout Dormancy, Total Tuber Number, Tuber Length, Yield, Chipping Colour, After-cooking Blackening, Cooking Type and Eye depth.

Detection of genetically modified organisms using highly multiplexed amplicon sequencing

The circulation of products based on genetically modified (GM) organisms is highly regulated by some governments with strict implementation rules for the breeding, planting, marketing, labelling, and trading of such products. To ensure compliance, accurate detection methods for GM events are necessary, along with assurance that GM material falls within relevant threshold levels. The increasing complexity and potential of undocumented GM are a growing challenge for genetic screening. Here, we developed and assessed a highly multiplexed amplicon sequencing assay for the detection of GM events based on a microfluidics platform and next-generation sequencing (NGS). To probe GM events comprehensively, we designed a total of 230 new amplicons to cover flanking, promoter, junction, and coding sequences of GM sequences. In addition, we designed and implemented parallel amplification of ribosomal and chloroplast markers to define crop species identity from potentially mixed samples. Using reference GM material of 11 crop species and multiple amplicons, we successfully detected the presence of 10 known modifications per GM event. We also find that reported flanking sequences of GM events may not be all useful for diagnostic. We assessed the assay's potential to detect GM events in mixed samples as well as in highly diluted DNA. Finally, we performed a prospective search of potentially undocumented GM events in plant material. Our microfluidics-based amplicon GM detection approach fills important gaps in detecting potentially undocumented and complex GM events by recovering a wide range of specific amplicon sequences for evaluation. Integrating highly parallel amplicon assays in GM screening efforts should be an effective complement to aid post-market monitoring and regulatory compliance efforts.

ColiSeq: a multiplex amplicon assay that provides strain level resolution of Escherichia coli directly from clinical specimens.

Escherichia coli is a diverse pathogen, causing a range of disease in humans, from self-limiting diarrhea to urinary tract infections (UTIs). Uropathogenic E. coli (UPEC) is the most frequently observed uropathogen in UTIs, a common disease in high-income countries, incurring billions of dollars yearly in treatment costs. Although E. coli is easily grown and identified in the clinical laboratory, genotyping the pathogen is more complicated, yet critical for reducing the incidence of disease. These goals can be achieved through whole-genome sequencing of E. coli isolates, but this approach is relatively slow and typically requires culturing the pathogen in the laboratory. To genotype E. coli rapidly and inexpensively directly from clinical samples, including but not limited to urine, we developed and validated a multiplex amplicon sequencing assay, called ColiSeq. The assay consists of targets designed for E. coli species confirmation, high resolution genotyping, and mixture deconvolution. To demonstrate its utility, we screened the ColiSeq assay against 230 clinical urine samples collected from a hospital system in Flagstaff, Arizona, USA. A limit of detection analysis demonstrated the ability of ColiSeq to identify E. coli at a concentration of ~2 genomic equivalent (GEs)/mL and to generate high-resolution genotyping at a concentration of 1 × 105 GEs/mL. The results of this study suggest that ColiSeq could be a valuable method to understand the source of UPEC strains and guide infection mitigation efforts. As sequence-based diagnostics become accepted in the clinical laboratory, workflows such as ColiSeq will provide actionable information to improve patient outcomes.IMPORTANCEUrinary tract infections (UTIs), caused primarily by Escherichia coli, create an enormous health care burden in the United States and other high-income countries. The early detection of E. coli from clinical samples, including urine, is important to target therapy and prevent further patient complications. Additionally, understanding the source of E. coli exposure will help with future mitigation efforts. In this study, we developed, tested, and validated an amplicon sequencing assay focused on direct detection of E. coli from urine. The resulting sequence data were demonstrated to provide strain level resolution of the pathogen, not only confirming the presence of E. coli, which can focus treatment efforts, but also providing data needed for source attribution and contact tracing. This assay will generate inexpensive, rapid, and reproducible data that can be deployed by public health agencies to track, diagnose, and potentially mitigate future UTIs caused by E. coli.

A mutation associated with resistance to synthetic pyrethroids is widespread in US populations of the tropical lineage of Rhipicephalus sanguineus s.l

The brown dog tick, Rhipicephalus sanguineus sensu lato (s.l.), is an important vector for Rickettsia rickettsii, causative agent of Rocky Mountain spotted fever. Current public health prevention and control efforts to protect people involve preventing tick infestations on domestic animals and in and around houses. Primary prevention tools rely on acaricides, often synthetic pyrethroids (SPs); resistance to this chemical class is widespread in ticks and other arthropods. Rhipicephalus sanguineus s.l. is a complex that likely contains multiple unique species and although the distribution of this complex is global, there are differences in morphology, ecology, and perhaps vector competence among these major lineages. Two major lineages within Rh. sanguineus s.l., commonly referred to as temperate and tropical, have been documented from multiple locations in North America, but are thought to occupy different ecological niches. To evaluate potential acaricide resistance and better define the distributions of the tropical and temperate lineages throughout the US and in northern Mexico, we employed a highly multiplexed amplicon sequencing approach to characterize sequence diversity at: 1) three loci within the voltage-gated sodium channel (VGSC) gene, which contains numerous genetic mutations associated with resistance to SPs; 2) a region of the gamma-aminobutyric acid-gated chloride channel gene (GABA-Cl) containing several mutations associated with dieldrin/fipronil resistance in other species; and 3) three mitochondrial genes (COI, 12S, and 16S). We utilized a geographically diverse set of Rh sanguineus s.l. collected from domestic pets in the US in 2013 and a smaller set of ticks collected from canines in Baja California, Mexico in 2021. We determined that a single nucleotide polymorphism (T2134C) in domain III segment 6 of the VGSC, which has previously been associated with SP resistance in Rh. sanguineus s.l., was widespread and abundant in tropical lineage ticks (>50 %) but absent from the temperate lineage, suggesting that resistance to SPs may be common in the tropical lineage. We found evidence of multiple copies of GABA-Cl in ticks from both lineages, with some copies containing mutations associated with fipronil resistance in other species, but the effects of these patterns on fipronil resistance in Rh. sanguineus s.l. are currently unknown. The tropical lineage was abundant and geographically widespread, accounting for 79 % of analyzed ticks and present at 13/14 collection sites. The temperate and tropical lineages co-occurred in four US states, and as far north as New York. None of the ticks we examined were positive for Rickettsia rickettsii or Rickettsia massiliae.

Haitian coffee agroforestry systems harbor complex arabica variety mixtures and under-recognized genetic diversity.

Though facing significant challenges, coffee (Coffea arabica) grown in Haitian agroforestry systems are important contributors to rural livelihoods and provide several ecosystem services. However, little is known about their genetic diversity and the variety mixtures used. In light of this, there is a need to characterize Haitian coffee diversity to help inform revitalization of this sector. We sampled 28 diverse farms in historically important coffee growing regions of northern and southern Haiti. We performed KASP-genotyping of SNP markers and HiPlex multiplex amplicon sequencing for haplotype calling on our samples, as well as several Ethiopian and commercial accessions from international collections. This allowed us to assign Haitian samples to varietal groups. Our analyses revealed considerable genetic diversity in Haitian farms, higher in fact than many farmers realized. Notably, genetic structure analyses revealed the presence of clusters related to Typica, Bourbon, and Catimor groups, another group that was not represented in our reference accession panel, and several admixed individuals. Across the study areas, we found both mixed-variety farms and monovarietal farms with the historical and traditional Typica variety. This study is, to our knowledge, the first to genetically characterize Haitian C. arabica variety mixtures, and report the limited cultivation of C. canephora (Robusta coffee) in the study area. Our results show that some coffee farms are repositories of historical, widely-abandoned varieties while others are generators of new diversity through genetic mixing.

Genomic surveillance of Anopheles mosquitoes on the Bijagós Archipelago using custom targeted amplicon sequencing identifies mutations associated with insecticide resistance

BackgroundInsecticide resistance is reducing the efficacy of vector control interventions, consequently threatening efforts to control vector-borne diseases, including malaria. Investigating the prevalence of molecular markers of resistance is a useful tool for monitoring the spread of insecticide resistance in disease vectors. The Bijagós Archipelago (Bijagós) in Guinea-Bissau is a region of stable malaria transmission where insecticide-treated nets are the mainstay for malaria control. However, the prevalence of molecular markers of insecticide resistance in malaria vectors is not well understood.MethodsA total of 214 Anopheles mosquitoes were analysed from 13 islands across the Bijagós. These mosquitoes were collected using CDC light traps in November 2019, during the peak malaria transmission season. High-throughput multiplex amplicon sequencing was used to investigate the prevalence of 17 different molecular markers associated with insecticide resistance in four genes: vgsc, rdl, ace1 and gste2.ResultsOf the 17 screened mutations, four were identified in mosquitoes from the Bijagós: vgsc L995F (12.2%), N1570Y (6.2%) and A1746S (0.7%) and rdl A269G (1.1%). This study is the first to report the L995F knock-down resistance (kdr)-west allele in Anopheles melas on the Archipelago. An additional eight non-synonymous single-nucleotide polymorphisms were identified across the four genes which have not been described previously. The prevalences of the vgsc L995F and N1570Y mutations were higher on Bubaque Island than on the other islands in this study; Bubaque is the most populous island in the archipelago, with the greatest population mobility and connection to continental Guinea-Bissau.ConclusionsThis study provides the first surveillance data for genetic markers present in malaria vectors from islands across the Bijagós Archipelago. Overall prevalence of insecticide resistance mutations was found to be low. However, the identification of the vgsc L995F and N1570Y mutations associated with pyrethroid resistance warrants further monitoring. This is particularly important as the mainstay of malaria control on the islands is the use of pyrethroid insecticide-treated nets.Graphical

Genetic diversity and phylogenetic analyses of Asian lineage Zika virus whole genome sequences derived from Culex quinquefasciatus mosquitoes and urine of patients during the 2020 epidemic in Thailand

Zika virus (ZIKV), a mosquito-borne flavivirus, has been continually emerging and re-emerging since 2010, with sporadic cases reported annually in Thailand, peaking at over 1000 confirmed positive cases in 2016. Leveraging high-throughput sequencing technologies, specifically whole genome sequencing (WGS), has facilitated rapid pathogen genome sequencing. In this study, we used multiplex amplicon sequencing on the Illumina Miseq instrument to describe ZIKV WGS. Six ZIKV WGS were derived from three samples of field-caught Culex quinquefasciatus mosquitoes (two males and one female) and three urine samples collected from patients in three different provinces of Thailand. Additionally, successful isolation of a ZIKV isolate occurred from a female Cx. quinquefasciatus. The WGS analysis revealed a correlation between the 2020 outbreak and the acquisition of five amino acid changes in the Asian lineage ZIKV strains from Thailand (2006), Cambodia (2010 and 2019), and the Philippines (2012). These changes, including C-T106A, prM-V1A, E-V473M, NS1-A188V, and NS5-M872V, were identified in all seven WGS, previously linked to significantly higher mortality rates. Furthermore, phylogenetic analysis indicated that the seven ZIKV sequences belonged to the Asian lineage. Notably, the genomic region of the E gene showed the highest nucleotide diversity (0.7–1.3%). This data holds significance in informing the development of molecular tools that enhance our understanding of virus patterns and evolution. Moreover, it may identify targets for improved methods to prevent and control future ZIKV outbreaks.

Assessing shared respiratory pathogens between domestic (Ovis aries) and bighorn (Ovis canadensis) sheep; methods for multiplex PCR, amplicon sequencing, and bioinformatics to characterize respiratory flora.

Respiratory disease is responsible for dramatic population declines in bighorn sheep (Ovis canadensis), and respiratory pathogen diagnostics contribute to the management of bighorn populations. To create a comprehensive and consistent approach to bighorn sheep respiratory diagnostics, we created a culture-independent assay to detect and strain type Mannheimia haemolytica, Bibersteinia trehalosi, Pasteurella multocida, and Mycoplasma ovipneumoniae. The assay also detects and characterizes the Pasteurellaceae leukotoxin A gene, and broadly assesses the bacterial composition of each sample based on 16S rRNA sequences. The assay is based on a three-step approach: 1) Multiplex PCR to amplify targets including eight loci for each bacterial species, the Pasteurellaceae lktA gene, and the 16S rRNA gene 2) Library preparation, barcoding, and short-read Illumina sequencing to determine the genetic sequences of each target, and 3) Bioinformatics in the form of automated software to analyze genetic sequences. The assay was designed to assess shared pathogens between domestic and bighorn sheep, but could be useful for many applications in bighorn sheep respiratory disease research and management.

Detection of Borrelia burgdorferi sensu lato species in host-seeking Ixodes species ticks in the United States

Lyme disease is the most commonly reported vector-borne disease in the United States and is transmitted by Ixodes scapularis in the eastern US and I. pacificus in the west. The causative agents, Borrelia burgdorferi sensu stricto (Bbss) and B. mayonii belong to the B. burgdorferi sensu lato (Bbsl) species complex. An additional eight species of Bbsl have been identified in Ixodes species ticks in the US, but their geographic distribution, vector associations, human encounter rates and pathogenicity in humans are poorly defined. To better understand the geographic distribution and vector associations of Bbsl spirochetes in frequent and infrequent human-biting Ixodes species ticks in the US, we previously screened 29,517 host-seeking I. scapularis or I. pacificus ticks and 692 ticks belonging to eight other Ixodes species for Borrelia spirochetes using a previously described tick testing algorithm that utilizes a combination of real-time PCR and Sanger sequencing for Borrelia species identification. The assay was designed to detect known human pathogens spread by Ixodes species ticks, but it was not optimized to detect Bbsl co-infections. To determine if such co-infections were overlooked particularly in ticks infected with Bbss, we retested and analyzed a subsample of 845 Borrelia infected ticks using a next generation sequencing multiplex PCR amplicon sequencing (MPAS) assay that can identify Borrelia species and Bbsl co-infections. The assay also includes targets that can molecularly confirm identifications of Ixodes species ticks to better inform pathogen-vector associations. We show that Bbss is the most prevalent species in I. scapularis and I. pacificus; other Bbsl species were rarely detected in I. scapularis and the only Bbsl co-infections identified in I. scapularis were with Bbss and B. mayonii. We detected B. andersonii in I. dentatus in the Mid-Atlantic and Upper Midwest regions, B. kurtenbachii in I. scapularis in the Upper Midwest, B. bissettiae in I. pacificus and I. spinipalpis in the Northwest, and B. carolinensis in I. affinis in the Mid-Atlantic and Southeast, and B. lanei in I. spinipalpis in the Northwest. Twelve of 62 (19.4%) Borrelia-infected I. affinis from the Mid-Atlantic region were co-infected with Bbss and B. carolinensis. Our data support the notion that Bbsl species are maintained in largely independent enzootic cycles, with occasional spill-over resulting in multiple Bbsl species detected in Ixodes species ticks.

Olivar: automated variant aware primer design for multiplex tiled amplicon sequencing of pathogens.

Tiled amplicon sequencing has served as an essential tool for tracking the spread and evolution of pathogens. Over 2 million complete SARS-CoV-2 genomes are now publicly available, most sequenced and assembled via tiled amplicon sequencing. While computational tools for tiled amplicon design exist, they require downstream manual optimization both computationally and experimentally, which is slow and costly. Here we present Olivar, a first step towards a fully automated, variant-aware design of tiled amplicons for pathogen genomes. Olivar converts each nucleotide of the target genome into a numeric risk score, capturing undesired sequence features that should be avoided. In a direct comparison with PrimalScheme, we show that Olivar has fewer SNPs overlapping with primers and predicted PCR byproducts. We also compared Olivar head-to-head with ARTIC v4.1, the most widely used primer set for SARS-CoV-2 sequencing, and show Olivar yields similar read mapping rates (~90%) and better coverage to the manually designed ARTIC v4.1 amplicons. We also evaluated Olivar on real wastewater samples and found that Olivar had up to 3-fold higher mapping rates while retaining similar coverage. In summary, Olivar automates and accelerates the generation of tiled amplicons, even in situations of high mutation frequency and/or density. Olivar is available as a web application at https://olivar.rice.edu. Olivar can also be installed locally as a command line tool with Bioconda. Source code, installation guide and usage are available at https://github.com/treangenlab/Olivar.

Molecular surveillance of Plasmodium falciparum drug-resistance markers in Vietnam using multiplex amplicon sequencing (2000–2016)

Emergence and spread of Plasmodium falciparum resistance to artemisinin-based combination therapies (ACT) is a major challenge for Greater Mekong Subregion countries in their goal to eliminate malaria by 2030. Tools to efficiently monitor drug resistance beyond resource-demanding therapeutic efficacy studies are necessary. A custom multiplex amplicon sequencing assay based on Illumina technology was designed to target the marker of partial resistance to artemisinin (K13), five candidate modulators of artemisinin resistance, the marker of resistance to chloroquine (crt), and four neutral microsatellite loci. The assay was used to genotype 635 P. falciparum-positive blood samples collected across seven provinces of Vietnam and one of Cambodia between 2000 and 2016. Markers of resistance to artemisinin partner-drugs piperaquine (copy number of plasmepsin-2) and mefloquine (copy number of multidrug-resistance 1) were determined by qPCR. Parasite population structure was further assessed using a 101-SNP barcode. Validated mutations of artemisinin partial resistance in K13 were found in 48.1% of samples, first detection was in 2000, and by 2015 prevalence overcame > 50% in Central Highlands and Binh Phuoc province. K13-C580Y variant became predominant country-wide, quickly replacing an outbreak of K13-I543T in Central Highlands. Mutations in candidate artemisinin resistance modulator genes paralleled the trends of K13 mutants, whereas resistance to piperaquine and mefloquine remained low (≈ 10%) by 2015–2016. Genomic tools applied to malaria surveillance generate comprehensive information on dynamics of drug resistance and population structure and reflect drug efficacy profiles from in vivo studies.

Characterization of the genetic composition and establishment of a core collection for the INERA Robusta coffee (Coffea canephora) field genebank from the Democratic Republic of Congo

Cultivation of Robusta coffee is likely to gain importance because of its high disease resistance and climate envelope. Robusta coffee genetic resources conserved in field genebanks can play an important role to further improve its cupping quality and other agronomic traits, but such Coffea canephora collections are limited and still poorly characterized. In this study, we characterized the genetic composition of the historically important but until recently neglected INERA Coffee Collection in Yangambi (the Democratic Republic of Congo). We used GBS to discover genome-wide genetic diversity, created and validated a novel multiplex amplicon sequencing (HiPlex) screening assay to genetically screen 730 coffee shrubs of the Yangambi Coffee Collection, grouped clonal material and delineated 263 accessions with unique genetic fingerprints. Comparison to reference material of three genetic origins revealed that the majority of the Yangambi accessions were assigned a ‘Lula’ cultivar origin, four accessions were assigned to Congolese subgroup A and nine accessions were most closely related to local wild accessions. About one-quarter of the accessions was likely derived from hybridization between these groups, which could result from seed-based propagation of the collection, breeding efforts, or natural cross-pollination. Parental analyses discovered eight preferentially used accessions, which may correspond to historically selected founders, or direct descendants thereof, whose seed material was once widely used to establish coffee plantations. Finally, two core collections were proposed using the maximization strategy (CC-I; 100 accessions) and genetic distance method (CC-X; 10 accessions). Our study demonstrates a method for the genetic characterization of Robusta coffee collections in general and contributes to the re-evaluation and exploration of the Robusta coffee genetic resources in the Democratic Republic of the Congo in particular.

Quantitative Amplicon Sequencing Is Necessary to Identify Differential Taxa and Correlated Taxa Where Population Sizes Differ.

High-throughput, multiplexed-amplicon sequencing has become a core tool for understanding environmental microbiomes. As researchers have widely adopted sequencing, many open-source analysis pipelines have been developed to compare microbiomes using compositional analysis frameworks. However, there is increasing evidence that compositional analyses do not provide the information necessary to accurately interpret many community assembly processes. This is especially true when there are large gradients that drive distinct community assembly processes. Recently, sequencing has been combined with Q-PCR (among other sources of total quantitation) to generate "Quantitative Sequencing" (QSeq) data. QSeq more accurately estimates the true abundance of taxa, is a more reliable basis for inferring correlation, and, ultimately, can be more reliably related to environmental data to infer community assembly processes. In this paper, we use a combination of published data sets, synthesis, and empirical modeling to offer guidance for which contexts QSeq is advantageous. As little as 5% variation in total abundance among experimental groups resulted in more accurate inference by QSeq than compositional methods. Compositional methods for differential abundance and correlation unreliably detected patterns in abundance and covariance when there was greater than 20% variation in total abundance among experimental groups. Whether QSeq performs better for beta diversity analysis depends on the question being asked, and the analytic strategy (e.g., what distance metric is being used); for many questions and methods, QSeq and compositional analysis are equivalent for beta diversity analysis. QSeq is especially useful for taxon-specific analysis; QSeq transformation and analysis should be the default for answering taxon-specific questions of amplicon sequence data. Publicly available bioinformatics pipelines should incorporate support for QSeq transformation and analysis.

Introduction of a multiplex amplicon sequencing assay to quantify DNA methylation in target cytosine markers underlying four selected epigenetic clocks

BackgroundDNA methylation analysis has proven to be a powerful tool for age assessment. However, the implementation of epigenetic age prediction in diagnostics or routine forensic casework requires appropriate laboratory methods. In this study, we aimed to compare the performance of large-scale DNA methylation analysis protocols that show promise in terms of accuracy, throughput, multiplexing capacity, and high sensitivity.ResultsThe protocols were designed to target a predefined panel of 161 genomic CG/CA sites from four known estimators of epigenetic age-related parameters, optimized and validated using artificially methylated controls or blood samples. We successfully targeted 96% of these loci using two enrichment protocols: Ion AmpliSeq™, an amplicon-based method integrated with Ion Torrent S5, and SureSelectXT Methyl-Seq, a hybridization-based method followed by MiSeq FGx sequencing. Both protocols demonstrated high accuracy and robustness. Although hybridization assays have greater multiplexing capabilities, the best overall performance was observed for the amplicon-based protocol with the lowest variability in DNA methylation at 25 ng of starting DNA, mean observed marker coverage of ~ 6.7 k reads, and accuracy of methylation quantification with a mean absolute difference between observed and expected methylation beta value of 0.054. The Ion AmpliSeq method correlated strongly with genome-scale EPIC microarray data (R = 0.91) and showed superiority in terms of methylation measurement accuracy. Method-to-method bias was accounted for by the use of linear transformation, which provided a highly accurate prediction of calendar age with a mean absolute error of less than 5 years for the VISAGE and Hannum age clocks used. The pace of aging (PoAm) and the mortality risk score (MRS) estimators included in our panel represent next-generation clocks, were found to have low to moderate correlations with the VISAGE and Hannum models (R < 0.75), and thus may capture different aspects of epigenetic aging.ConclusionsWe propose a laboratory tool that allows the quantification of DNA methylation in cytosines underlying four different clocks, thus providing broad information on epigenetic aging while maintaining a reasonable number of CpG markers, opening the way to a wide range of applications in forensics, medicine, and healthcare.

Ticks and tick-borne microbes identified through passive and active surveillance in Alaska.

Rapid environmental change in Alaska and other regions of the Arctic and sub-Arctic has raised concerns about increasing human exposure to ticks and the pathogens they carry. We tested a sample of ticks collected through a combination of passive and active surveillance from humans, domestic animals, and wildlife hosts in Alaska for a panel of the most common tick-borne pathogens in the contiguous United States to characterize the diversity of microbes present in this region. We tested 189 pooled tick samples collected in 2019-2020 for Borrelia spp., Anaplasma spp., Ehrlichia spp., and Babesia spp. using a multiplex PCR amplicon sequencing assay. We found established populations of Ixodes angustus Neumann (Acari: Ixodidae), Ixodes uriae White (Acari: Ixodidae), and Haemaphysalis leporispalustris Packard (Acari: Ixodidae) in Alaska, with I. angustus found on a variety of hosts including domestic companion animals (dogs and cats), small wild mammals, and humans. Ixodes angustus were active from April through October with peaks in adult and nymphal activity observed in summer months (mainly July). Although no known human pathogens were detected, Babesia microti-like parasites and candidatus Ehrlichia khabarensis were identified in ticks and small mammals. The only human pathogen detected (B. burgdorferi s.s.) was found in a tick associated with a dog that had recently traveled to New York, where Lyme disease is endemic. This study highlights the value of a combined passive and active tick surveillance system to detect introduced tick species and pathogens and to assess which tick species and microbes are locally established.

A bioinformatics pipeline for a tick pathogen surveillance multiplex amplicon sequencing assay

The Centers for Disease Control and Prevention's national tick and tick-borne pathogen surveillance program collects information to better understand the regional distribution, prevalence, and exposure risk of host-seeking medically important ticks in the United States. A recently developed next generation sequencing (NGS) targeted multiplex PCR amplicon sequencing (MPAS) assay has enhanced the detection capabilities for Ixodes-associated human pathogens found in Ixodes scapularis and Ixodes pacificus ticks compared to the routinely used real-time PCR assay. To operationalize the MPAS assay for the large number of tick surveillance submissions processed each year, a reproducible high throughput bioinformatics pipeline is needed. We describe the development and validation of the MPAS pipeline, a bioinformatics pipeline that identifies and summarizes amplicon sequences produced by the MPAS assay. This pipeline is portable and reproducible across different computing environments, and flexible by allowing modifications to input parameters, assay primer and reference sequences. The automation of the summary report, BLAST report, and phylogenetic analysis reduces the amount of time needed for downstream analysis. To validate this pipeline, we compared the analysis of a MPAS assay dataset consisting of 175 I. scapularis nymphs with the MPAS pipeline and previously published results analyzed with a CLC Genomic Workbench workflow. The MPAS pipeline identified the same number of positive ticks for Anaplasma phagocytophilum and Babesia species as the original analysis, but the MPAS pipeline provided enhanced sequencing resolution of Borrelia burgdorferi sensu lato co-infected samples. The reproducibility, flexibility, analysis automation, and improved sequence resolution of the MPAS pipeline make it well suited for a high throughput tick pathogen surveillance program.

Characterization of 15 CYP2J2 variants identified in the Chinese Han population on the metabolism of ebastine and terfenadine in vitro.

Genetic polymorphism of the cytochrome P450 (CYP) gene can significantly influence the metabolism of endogenous and xenobiotic compounds. However, few studies have focused on the polymorphism of CYP2J2 and its impact on drug catalytic activity, especially in the Chinese Han population. In this study, we sequenced the promoter and exon regions of CYP2J2 in 1,163 unrelated healthy Chinese Han individuals using the multiplex PCR amplicon sequencing method. Then, the catalytic activities of the detected CYP2J2 variants were evaluated after recombinant expression in S. cerevisiae microsomes. As a result, CYP2J2*7, CYP2J2*8, 13 variations in the promoter region and 15 CYP2J2 nonsynonymous variants were detected, of which V15A, G24R, V68A, L166F and A391T were novel missense variations. Immunoblotting results showed that 11 of 15 CYP2J2 variants exhibited lower protein expression than wild-type CYP2J2.1. In vitro functional analysis results revealed that the amino acid changes of 14 variants could significantly influence the drug metabolic activity of CYP2J2 toward ebastine or terfenadine. Specifically, 4 variants with relatively higher allele frequencies, CYP2J2.8, 173_173del, K267fs and R446W, exhibited extremely low protein expression and defective catalytic activities for both substrates. Our results indicated that a high genetic polymorphism of CYP2J2 could be detected in the Chinese Han population, and most genetic variations in CYP2J2 could influence the expression and catalytic activity of CYP2J2. Our data significantly enrich the knowledge of genetic polymorphisms in CYP2J2 and provide new theoretical information for corresponding individualized medication in Chinese and other Asian populations.

Antibiotic resistance gene sequencing is necessary to reveal the complex dynamics of immigration from sewers to activated sludge.

Microbial community composition has increasingly emerged as a key determinant of antibiotic resistance gene (ARG) content. However, in activated sludge wastewater treatment plants (AS-WWTPs), a comprehensive understanding of the microbial community assembly process and its impact on the persistence of antimicrobial resistance (AMR) remains elusive. An important part of this process is the immigration dynamics (or community coalescence) between the influent and activated sludge. While the influent wastewater contains a plethora of ARGs, the persistence of a given ARG depends initially on the immigration success of the carrying population, and the possible horizontal transfer to indigenously resident populations of the WWTP. The current study utilized controlled manipulative experiments that decoupled the influent wastewater composition from the influent microbial populations to reveal the fundamental mechanisms involved in ARG immigration between sewers and AS-WWTP. A novel multiplexed amplicon sequencing approach was used to track different ARG sequence variants across the immigration interface, and droplet digital PCR was used to quantify the impact of immigration on the abundance of the targeted ARGs. Immigration caused an increase in the abundance of over 70 % of the quantified ARGs. However, monitoring of ARG amplicon sequence variants (ARG-ASVs) at the immigration interface revealed various immigration patterns such as (i) suppression of the indigenous mixed liquor ARG-ASV by the immigrant, or conversely (ii) complete immigration failure of the influent ARG-ASV. These immigration profiles are reported for the first time here and highlight the crucial information that can be gained using our novel multiplex amplicon sequencing techniques. Future studies aiming to reduce AMR in WWTPs should consider the impact of influent immigration in process optimisation and design.