Simulate_PCR for amplicon prediction and annotation from multiplex, degenerate primers and probes

Development of a rapid and sensitive real-time diagnostic assay to detect and quantify Aphanomyces invadans, the causative agent of epizootic ulcerative syndrome

The oomycete Aphanomyces invadans causes epizootic ulcerative syndrome (EUS), a World Organization for Animal Health (WOAH)-listed disease that has seriously impacted a wide range of fish worldwide. Currently, only three conventional polymerase chain reaction (PCR) assays are recommended for the detection of A. invadans. The robust quantitative PCR (qPCR) assay has recently become more important due to its highly accurate nature and the applicability of qPCR-based environmental DNA (eDNA) detection in the monitoring of pathogens in aquatic environments. Therefore, in this study, we developed a novel TaqMan probe-based qPCR method to sensitively and quantitatively detect A. invadans. The assay limit of detection was determined using 10-fold serial dilutions of linearized A. invadans plasmid. Assay sensitivity was assessed in the presence of interfering substances and compared to three WOAH-listed primers using the mycelia and zoospores of A. invadans with and without fish muscle tissue. The assay specificity was also theoretically and experimentally assessed against other oomycetes, fish muscle tissue, and water samples. The assay's repeatability and reproducibility were determined. In this study, the limit of detection of the developed assay was 7.24 copies of A. invadans genomic DNA per reaction (95% confidence interval (CI): 2.75 to 19.05 copies/reaction). The assay showed the same sensitivity in the presence of other substances. Compared to the WOAH-recommended PCR assays, this assay had 10-times higher sensitivity for all tested samples. There were no cross-reactions with other closely related oomycetes, fish muscle, or water samples, indicating that the assay was highly specific for A. invadans. The repeatability and reproducibility tests showed little variation, ranging from 0.1-0.9% and 0.04-1.1%, respectively, indicating the high consistency, repeatability, and reliability of the developed assay. This highly rapid, sensitive, specific, and consistent EUS qPCR assay would be of importance in transboundary disease management and the monitoring of pathogens in aquatic environments.

Development, testing and validation of a SARS-CoV-2 multiplex panel for detection of the five major variants of concern on a portable PCR platform

Many SARS-CoV-2 variants have emerged during the course of the COVID-19 pandemic. These variants have acquired mutations conferring phenotypes such as increased transmissibility or virulence, or causing diagnostic, therapeutic, or immune escape. Detection of Alpha and the majority of Omicron sublineages by PCR relied on the so-called S gene target failure due to the deletion of six nucleotides coding for amino acids 69-70 in the spike (S) protein. Detection of hallmark mutations in other variants present in samples relied on whole genome sequencing. However, whole genome sequencing as a diagnostic tool is still in its infancy due to geographic inequities in sequencing capabilities, higher cost compared to other molecular assays, longer turnaround time from sample to result, and technical challenges associated with producing complete genome sequences from samples that have low viral load and/or high background. Hence, there is a need for rapid genotyping assays. In order to rapidly generate information on the presence of a variant in a given sample, we have created a panel of four triplex RT-qPCR assays targeting 12 mutations to detect and differentiate all five variants of concern: Alpha, Beta, Gamma, Delta, and Omicron. We also developed an expanded pentaplex assay that can reliably distinguish among the major sublineages (BA.1-BA.5) of Omicron. In silico, analytical and clinical testing of the variant panel indicate that the assays exhibit high sensitivity and specificity. This panel can help fulfill the need for rapid identification of variants in samples, leading to quick decision making with respect to public health measures, as well as treatment options for individuals. Compared to sequencing, these genotyping PCR assays allow much faster turn-around time from sample to results-just a couple hours instead of days or weeks.

Relative virulence of Staphylococcus aureus bovine mastitis strains representing the main Canadian spa types and clonal complexes as determined using in vitro and in vivo mastitis models

Staphylococcus aureus is one of the main pathogens leading to both clinical and subclinical bovine mastitis in dairy cattle. Prediction of disease evolution based on the characteristics of Staph. aureus isolates that cause intramammary infections and understanding the host-pathogen interactions may improve management of mastitis in dairy herds. For this study, several strains were selected from each of the 6 major Canadian spa types associated with mastitis (t267, t359, t529, t605, t2445, and t13401). Adherence to host cells and intracellular persistence of these strains were studied using a bovine mammary gland epithelial cell line (MAC-T). Additionally, relative virulence and host response (cytokines production) were also studied in vivo using a mouse model of mastitis. Whole-genome sequencing was performed on all strains and associations between clonal complex, sequence type, and presence of certain virulence factors were also investigated. Results show that spa type t2445 was correlated with persistence in MAC-T cells. Strains from spa t359 and t529 showed better ability to colonize mouse mammary glands. The exception was strain sa3154 (spa t529), which showed less colonization of glands compared with other t359 and t529 strains but possessed the highest number of superantigen genes including tst. All strains possessed hemolysins, but spa types t529 and t2445 showed the largest diameter of β-hemolysis on blood agar plates. Although several spa types possessed 2 or 3 serine-aspartate rich proteins (Sdr) believed to be involved in many pathogenic processes, most t529 strains expressed only an allelic variant of sdrE. The spa types t605 (positive for the biofilm associated protein gene; bap+) and t13401 (bap-), that produced the largest amounts of biofilm in vitro, were the least virulent in vivo. Finally, strains from spa type t529 (ST151) elicited a cytokine expression profile (TNF-α, IL-1β and IL-12) that suggests a potential for severe inflammation. This study suggests that determination of the spa type may help predict the severity of the disease and the ability of the immune system to eliminate intramammary infections caused by Staph. aureus.

Chromosome-Level Genome Assemblies Expand Capabilities of Genomics for Conservation Biology

Genome assemblies are in the process of becoming an increasingly important tool for understanding genetic diversity in threatened species. Unfortunately, due to limited budgets typical for the area of conservation biology, genome assemblies of threatened species, when available, tend to be highly fragmented, represented by tens of thousands of scaffolds not assigned to chromosomal locations. The recent advent of high-throughput chromosome conformation capture (Hi-C) enables more contiguous assemblies containing scaffolds spanning the length of entire chromosomes for little additional cost. These inexpensive contiguous assemblies can be generated using Hi-C scaffolding of existing short-read draft assemblies, where N50 of the draft contigs is larger than 0.1% of the estimated genome size and can greatly improve analyses and facilitate visualization of genome-wide features including distribution of genetic diversity in markers along chromosomes or chromosome-length scaffolds. We compared distribution of genetic diversity along chromosomes of eight mammalian species, including six listed as threatened by IUCN, where both draft genome assemblies and newer chromosome-level assemblies were available. The chromosome-level assemblies showed marked improvement in localization and visualization of genetic diversity, especially where the distribution of low heterozygosity across the genomes of threatened species was not uniform.

Microbiota and Metabolite Modifications after Dietary Exclusion of Dairy Products and Reduced Consumption of Fermented Food in Young and Older Men

The gut microbiota adapts to age-related changes in host physiology but is also affected by environmental stimuli, like diet. As a source of both pre- and probiotics, dairy and fermented foods modulate the gut microbiota composition, which makes them interesting food groups to use for the investigation of interactions between diet and ageing. Here we present the effects of excluding dairy products and limiting fermented food consumption for 19 days on gut microbiota composition and circulating metabolites of 28 healthy, young (YA) and older (OA) adult men. The intervention affected gut microbial composition in both groups, with significant increases in Akkermansia muciniphila and decreases in bacteria of the Clostridiales order. Lower fasting levels of glucose and insulin, as well as dairy-associated metabolites like lactose and pentadecanoic acid, were observed after the intervention, with no effect of age. The intervention also decreased HDL and LDL cholesterol levels. Dairy fat intake was positively associated with the HDL cholesterol changes but not with the LDL/HDL ratio. In conclusion, restricting the intake of dairy and fermented foods in men modified their gut microbiota and blood metabolites, while the impact of the dietary restrictions on these outcomes was more marked than the effect of age.

Mycoplasma bovis in Nordic European Countries: Emergence and Dominance of a New Clone

Mycoplasma (M.) bovis is an important pathogen of cattle implicated in a broad range of clinical manifestations that adversely impacts livestock production worldwide. In the absence of a safe, effective commercial vaccine in Europe, the reported reduced susceptibility to antimicrobials for this organism has contributed to difficulties in controlling infection. Despite global presence, some countries have only recently experienced outbreaks of this pathogen. In the present study, M. bovis isolates collected in Denmark between 1981 and 2016 were characterized to determine (i) genetic diversity and phylogenetic relationships using whole genome sequencing and various sequence-based typing methods and (ii) patterns of antimicrobial resistance compared to other European isolates. The M. bovis population in Denmark was found to be highly homogeneous genomically and with respect to the antimicrobial resistance profile. Previously dominated by an old genotype shared by many other countries (ST17 in the PubMLST legacy scheme), a new predominant type represented by ST94-adh1 has emerged. The same clone is also found in Sweden and Finland, where M. bovis introduction is more recent. Although retrieved from the Netherlands, it appears absent from France, two countries with a long history of M. bovis infection where the M. bovis population is more diverse.

Quantification of nosZ genes and transcripts in activated sludge microbiomes with novel group-specific qPCR methods validated with metagenomic analyses

Substantial N₂O emission results from activated sludge nitrogen removal processes. N₂O-reducing organisms possessing NosZ-type N₂O reductases have been recognized to play crucial roles in suppressing emission of N₂O produced in anoxic activated sludge via denitrification; however, which of the diverse nosZ-possessing organisms function as the major N₂O sink in situ remains largely unknown. Here, nosZ genes and transcripts in wastewater microbiomes were analyzed with the group-specific qPCR assays designed de novo combining culture-based and computational approaches. A sewage sample was enriched in a batch reactor fed continuous stream of N₂ containing 20-10,000 ppmv N₂O with excess amount (10 mM) of acetate as the source of carbon and electrons, where 14 genera of potential N₂O-reducers were identified. All available amino acid sequences of NosZ affiliated to these taxa were grouped into five subgroups (two clade I and three clade II groups), and primers/probe sets exclusively and comprehensively targeting the subgroups were designed and validated with in silico PCR. Four distinct activated sludge samples from three different wastewater treatment plants in Korea were analyzed with the qPCR assays and the results were validated with the shotgun metagenome analysis results. With these group-specific qPCR assays, the nosZ genes and transcripts of six additional activated sludge samples were analyzed and the results of the analyses clearly indicated the dominance of two clade II nosZ subgroups (Flavobacterium-like and Dechloromonas-like) among both nosZ gene and transcript pools.

SNP Genotyping with Target Amplicon Sequencing Using a Multiplexed Primer Panel and Its Application to Genomic Prediction in Japanese Cedar, Cryptomeria japonica (L.f.) D.Don

Along with progress in sequencing technology and accumulating knowledge of genome and gene sequences, molecular breeding techniques have been developed for predicting the genetic potential of individual genotypes and for selecting superior individuals. For Japanese cedar (Cryptomeria japonica (L.f.) D.Don), which is the most common coniferous species in Japanese forestry, we constructed a custom primer panel for target amplicon sequencing in order to simultaneously determine 3034 informative single nucleotide polymorphisms (SNPs). We performed primary evaluation of the custom primer panel with actual sequencing and in silico PCR. Genotyped SNPs had a distribution over almost the entire region of the C. japonica linkage map and verified the high reproducibility of genotype calls compared to SNPs obtained by genotyping arrays. Genotyping was performed for 576 individuals of the F1 population, and genomic prediction models were constructed for growth and wood property-related traits using the genotypes. Amplicon sequencing with the custom primer panel enables efficient obtaining genotype data in order to perform genomic prediction, manage clones, and advance forest tree breeding.

Large scale genome skimming from herbarium material for accurate plant identification and phylogenomics

BACKGROUND

Herbaria are valuable sources of extensive curated plant material that are now accessible to genetic studies because of advances in high-throughput, next-generation sequencing methods. As an applied assessment of large-scale recovery of plastid and ribosomal genome sequences from herbarium material for plant identification and phylogenomics, we sequenced 672 samples covering 21 families, 142 genera and 530 named and proposed named species. We explored the impact of parameters such as sample age, DNA concentration and quality, read depth and fragment length on plastid assembly error. We also tested the efficacy of DNA sequence information for identifying plant samples using 45 specimens recently collected in the Pilbara.

RESULTS

Genome skimming was effective at producing genomic information at large scale. Substantial sequence information on the chloroplast genome was obtained from 96.1% of samples, and complete or near-complete sequences of the nuclear ribosomal RNA gene repeat were obtained from 93.3% of samples. We were able to extract sequences for the core DNA barcode regions rbcL and matK from 96 to 93.3% of samples, respectively. Read quality and DNA fragment length had significant effects on sequencing outcomes and error correction of reads proved essential. Assembly problems were specific to certain taxa with low GC and high repeat content (Goodenia, Scaevola, Cyperus, Bulbostylis, Fimbristylis) suggesting biological rather than technical explanations. The structure of related genomes was needed to guide the assembly of repeats that exceeded the read length. DNA-based matching proved highly effective and showed that the efficacy for species identification declined in the order cpDNA >> rDNA > matK >> rbcL.

CONCLUSIONS

We showed that a large-scale approach to genome sequencing using herbarium specimens produces high-quality complete cpDNA and rDNA sequences as a source of data for DNA barcoding and phylogenomics.

Whole-genome comparative analysis of Campylobacter jejuni strains isolated from patients with diarrhea in northeastern Poland

BACKGROUND

Campylobacter jejuni is the leading cause of bacterial gastroenteritis (campylobacteriosis) in humans worldwide, and the most frequent pathogen associated with Guillain-Barré syndrome (GBS) and Miller-Fisher syndrome (MFS). The study was designed in order to assess similarities between genomes of Campylobacter jejuni strains, isolated from children suffering from acute diarrhea in northeastern Poland, in comparison to C. jejuni genomes stored in public databases. The analysis involved phylogeny, resistome and virulome. In addition, the Campylobacter PubMLST database was used to estimate the prevalence of the analyzed C. jejuni sequence type (STs) in other countries.

RESULTS

Campylobacter jejuni ST50, ST257 and ST51 represented 5.3%, 4.5% and 2.2% of the PubMLST records, respectively. Overall, strains representing the STs showed common resistance to tetracyclines (51.3%) and fluoroquinolones (31.8%), mediated through the tetO gene (98.2%) and point mutation (T86I) in the gyrA gene (100%). However, the latter was present in all our isolates. The major differences in virulence patterns concerned serotypes, lipooligosaccharide (LOS) classes and certain clinically relevant genes.

CONCLUSIONS

Campylobacter jejuni ST50, ST51 and ST257 are among the top ten of STs isolated in Europe. WGS revealed diversity of serotypes and LOS classes in ST50 strains, that deserves further clinical and epidemiological investigations as it might be related to a risk of post-infectious neurological sequels such as Guillain-Barré syndrome. Additionally, the results implicate lower pathogenic potential and distinct transmission chains or reservoirs for C. jejuni ST51 isolates responsible for campylobacteriosis in northeastern Poland.

A single test approach for accurate and sensitive detection and taxonomic characterization of Trypanosomes by comprehensive analysis of internal transcribed spacer 1 amplicons

To improve our knowledge on the epidemiological status of African trypanosomiasis, better tools are required to monitor Trypanosome genotypes circulating in both mammalian hosts and tsetse fly vectors. This is important in determining the diversity of Trypanosomes and understanding how environmental factors and control efforts affect Trypanosome evolution. We present a single test approach for molecular detection of different Trypanosome species and subspecies using newly designed primers to amplify the Internal Transcribed Spacer 1 region of ribosomal RNA genes, coupled to Illumina sequencing of the amplicons. The protocol is based on Illumina’s widely used 16s bacterial metagenomic analysis procedure that makes use of multiplex PCR and dual indexing. Results from analysis of wild tsetse flies collected from Zambia and Zimbabwe show that conventional methods for Trypanosome species detection based on band size comparisons on gels is not always able to accurately distinguish between T. vivax and T. godfreyi. Additionally, this approach shows increased sensitivity in the detection of Trypanosomes at species level with the exception of the Trypanozoon subgenus. We identified subspecies of T. congolense, T. simiae, T. vivax, and T. godfreyi without the need for additional tests. Results show T. congolense Kilifi subspecies is more closely related to T. simiae than to other T. congolense subspecies. This agrees with previous studies using satellite DNA and 18s RNA analysis. While current classification does not list any subspecies for T. godfreyi, we observed two distinct clusters for these species. Interestingly, sequences matching T. congolense Tsavo (now classified as T. simiae Tsavo) clusters distinctly from other T. simiae Tsavo sequences suggesting the Nannomonas group is more divergent than currently thought thus the need for better classification criteria. This method presents a simple but comprehensive way of identification of Trypanosome species and subspecies-specific using one PCR assay for molecular epidemiology of trypanosomes.

Tsetse flies are central actors in the transmission of Trypanosomes to vertebrate hosts. Therefore, detection of Trypanosomes in the tsetse flies is important for understanding the epidemiology of African trypanosomiasis as a component of new control or surveillance strategies. We have developed a method that combines multiplex PCR and next-generation sequencing for the detection of different Trypanosome species and subspecies. Similar to the widely used bacterial metagenomic analysis protocol, this method uses a modular, two-step PCR process followed by sequencing of all amplicons in a single run, making sequencing of amplicons more efficient and cost-effective when dealing with large sample sizes. As part of this approach, we designed novel Internal Transcribed Spacer 1 primers optimized for short read sequencing and have slightly better sensitivity than conventional primers. Taxonomic identification of amplicons is based on BLAST searches against the constantly updated NCBI’s nt database. Our approach is more accurate than traditional gel-based analyses which are prone to misidentification of species. It is also able to discriminate between subspecies of T. congolense, T. simiae, T. vivax, and T. godfreyi species. This method has the potential to provide new insights into the epidemiology of different Trypanosome genotypes and the discovery of new ones.

Rapid multiplex gene expression assays for monitoring metabolic resistance in the major malaria vector Anopheles gambiae

BACKGROUND

Metabolic resistance of the major malaria vector Anopheles gambiae (s.l.) to insecticides is operationally significant, particularly in combination with target site resistance. However, detection of metabolic resistance is not trivial and relies on laborious bioassays, unspecific biochemical methods, or sophisticated and expensive molecular approaches using transcriptomics.

METHODS

Rapid one-step multiplex TaqMan-probe based RT-qPCR assays were developed and optimised to measure the expression levels of genes associated with metabolic insecticide resistance in An. gambiae (s.l.). Primers and probes were designed to target the mRNA of cytochrome P450-dependent monooxygenases CYP6P3, CYP6M2, CYP9K1, CYP6P4 and CYP6Z1, and the glutathione-S-transferase GSTE2. The novel assays were validated versus gold standard methods with a range of phenotyped mosquito specimens. The assays were also tested directly on lysates of RNAlater®-preserved mosquitoes without an RNA extraction step.

RESULTS

The novel assays are efficient (reaction efficiencies = 95-109%), sensitive (covering a > 10.0 Ct range with R² values > 0.99), specific (TaqMan chemistry), reproducible (%CV = 4.46-12.07%), as well as readily expandable to capture additional loci as they evolve or to cover additional species. The assays were successfully validated in terms of expression levels against standard two-step singleplex qPCR assays (overall % difference = -17.6%, 95% CI = -38.7-3.43%) and microarrays, using laboratory strains and field-caught samples. The assays can also be applied directly on lysates of mosquito specimens, without RNA extraction or DNase treatment.

CONCLUSIONS

The novel multiplex assays for monitoring the levels of major detoxification genes and metabolic resistance in An. gambiae (s.l.) are simple to perform, robust and rapid. They may complement current diagnostic assays to provide evidence-based and operationally relevant information for insecticide resistance management.

FastPCR: An in silico tool for fast primer and probe design and advanced sequence analysis

Polymerase chain reaction (PCR) is one of the most important laboratory techniques used in molecular biology, genetics and molecular diagnostics. The success of a PCR-based method largely depends on the correct nucleic acid sequence analysis in silico prior to a wet-bench experiment. Here, we report the development of an online Java-based software for virtual PCR on linear or circular DNA templates and multiple primer or probe search from large or small databases. Primer or probe sensitivity and specificity are predicted by searching a database to find sequences with an optimal number of mismatches, similarity and stability. The software determines primer location, orientation, efficiency of binding and calculates primer melting temperatures for standard and degenerate oligonucleotides. The software is suitable for batch file processing, which is essential for automation when working with large amounts of data. The online Java software is available for download at http://primerdigital.com/tools/pcr.html. Accession numbers for the sequences resulting from this study: EU140956 EU177767 EU867815 EU882730 FJ975775-FJ975780 HM481419 HM481420 KC686837-KC686839 KM262797.

Ribosomal background of the Bacillus cereus group thermotypes

In this study we reconstructed the architecture of Bacillus cereus sensu lato population based on ribosomal proteins, and identified a link between the ribosomal proteins’ variants and thermal groups (thermotypes) of the bacilli. The in silico phyloproteomic analysis of 55 ribosomal proteins (34 large and 21 small subunit r-proteins) of 421 strains, representing 14 well-established or plausible B. cereus sensu lato species, revealed several ribosomal clusters (r-clusters), which in general were well correlated with the strains’ affiliation to phylogenetic/thermal groups I–VII. However, a conformity and possibly a thermal characteristic of certain phylogenetic groups, e.g. the group IV, were not supported by a distribution of the corresponding r-clusters, and consequently neither by the analysis of cold-shock proteins (CSPs) nor by a content of heat shock proteins (HSPs). Furthermore, a preference for isoleucine and serine over valine and alanine in r-proteins along with a lack of HSP16.4 were recognized in non-mesophilic thermotypes. In conclusion, we suggest that the observed divergence in ribosomal proteins may be connected with an adaptation of B. cereus sensu lato members to various thermal niches.

In Silico PCR Tools for a Fast Primer, Probe, and Advanced Searching

The polymerase chain reaction (PCR) is fundamental to molecular biology and is the most important practical molecular technique for the research laboratory. The principle of this technique has been further used and applied in plenty of other simple or complex nucleic acid amplification technologies (NAAT). In parallel to laboratory "wet bench" experiments for nucleic acid amplification technologies, in silico or virtual (bioinformatics) approaches have been developed, among which in silico PCR analysis. In silico NAAT analysis is a useful and efficient complementary method to ensure the specificity of primers or probes for an extensive range of PCR applications from homology gene discovery, molecular diagnosis, DNA fingerprinting, and repeat searching. Predicting sensitivity and specificity of primers and probes requires a search to determine whether they match a database with an optimal number of mismatches, similarity, and stability. In the development of in silico bioinformatics tools for nucleic acid amplification technologies, the prospects for the development of new NAAT or similar approaches should be taken into account, including forward-looking and comprehensive analysis that is not limited to only one PCR technique variant. The software FastPCR and the online Java web tool are integrated tools for in silico PCR of linear and circular DNA, multiple primer or probe searches in large or small databases and for advanced search. These tools are suitable for processing of batch files that are essential for automation when working with large amounts of data. The FastPCR software is available for download at http://primerdigital.com/fastpcr.html and the online Java version at http://primerdigital.com/tools/pcr.html .

Bioinformatic approach in the identification of arabidopsis gene homologous in amaranthus

Bioinfomatics offers an efficient tool for molecular genetics applications and sequence homology search algorithms became an inevitable part for many different research strategies. Appropriate managing of known data that are stored in public available databases can be used in many ways in the research. Here, we report the identification of RmlC-like cupins superfamily protein DNA sequence than is known in Arabidopsis genome for the Amaranthus - plant specie where this sequence was still not sequenced. A BLAST based approach was used to identify the homologous sequences in the nucleotide database and to find suitable parts of the Arabidopsis sequence were primers can be designed. In total, 64 hits were found in nucleotide database for Arabidopsis RmlC-like cupins sequence. A query cover ranged from 10% up to the 100% among RmlC-like cupins nucleotides and its homologues that are actually stored in public nucleotide databases. The most conserved region was identified for matches that posses nucleotides in the range of 1506 up to the 1925 bp of RmlC-like cupins DNA sequence stored in the database. The in silico approach was subsequently used in PCR analysis where the specifity of designed primers was approved. A unique, 250 bp long fragment was obtained for Amaranthus cruentus and a hybride Amaranthus hypochondriacus x hybridus in our analysis. Bioinformatic based analysis of unknown parts of the plant genomes as showed in this study is a very good additional tool in PCR based analysis of plant variability. This approach is suitable in the case for plants, where concrete genomic data are still missing for the appropriate genes, as was demonstrated for Amaranthus. 

Multiplex degenerate primer design for targeted whole genome amplification of many viral genomes

Background. Targeted enrichment improves coverage of highly mutable viruses at low concentration in complex samples. Degenerate primers that anneal to conserved regions can facilitate amplification of divergent, low concentration variants, even when the strain present is unknown. Results. A tool for designing multiplex sets of degenerate sequencing primers to tile overlapping amplicons across multiple whole genomes is described. The new script, run_tiled_primers, is part of the PriMux software. Primers were designed for each segment of South American hemorrhagic fever viruses, tick-borne encephalitis, Henipaviruses, Arenaviruses, Filoviruses, Crimean-Congo hemorrhagic fever virus, Rift Valley fever virus, and Japanese encephalitis virus. Each group is highly diverse with as little as 5% genome consensus. Primer sets were computationally checked for nontarget cross reactions against the NCBI nucleotide sequence database. Primers for murine hepatitis virus were demonstrated in the lab to specifically amplify selected genes from a laboratory cultured strain that had undergone extensive passage in vitro and in vivo. Conclusions. This software should help researchers design multiplex sets of primers for targeted whole genome enrichment prior to sequencing to obtain better coverage of low titer, divergent viruses. Applications include viral discovery from a complex background and improved sensitivity and coverage of rapidly evolving strains or variants in a gene family.