Search

The global spread and continued evolution of SARS-CoV-2 has driven an unprecedented surge in viral genomic surveillance. Amplicon-based sequencing methods provide a sensitive, low-cost and rapid approach but suffer a high potential for contamination, which can undermine laboratory processes and results. This challenge will increase with the expanding global production of sequences across a variety of laboratories for epidemiological and clinical interpretation, as well as for genomic surveillance of emerging diseases in future outbreaks. We present SDSI + AmpSeq, an approach that uses 96 synthetic DNA spike-ins (SDSIs) to track samples and detect inter-sample contamination throughout the sequencing workflow. We apply SDSIs to the ARTIC Consortium's amplicon design, demonstrate their utility and efficiency in a real-time investigation of a suspected hospital cluster of SARS-CoV-2 cases and validate them across 6,676 diagnostic samples at multiple laboratories. We establish that SDSI + AmpSeq provides increased confidence in genomic data by detecting and correcting for relatively common, yet previously unobserved modes of error, including spillover and sample swaps, without impacting genome recovery.

Background Several Multilocus Sequence Typing (MLST) schemes have been developed for Chlamydia trachomatis. Bom’s MLST scheme for MLST is based on nested PCR amplification and sequencing of five hypervariable genes and ompA. In contrast to other Chlamydia MLST schemes, Bom’s MLST scheme gives higher resolution and phylogenetic trees that are comparable to those from whole genome sequencing. However, poor results have been obtained with Bom’s MLST scheme in clinical samples with low concentrations of Chlamydia DNA. Results In this work, we present an improved version of the scheme that is based on the same genes and MLST database as Bom’s MLST scheme, but with newly designed primers for nested-1 and nested-2 steps under stringent conditions. Furthermore, we introduce a third primer set for the sequencing step, which considerably improves the performance of the assay. The improved primers were tested in-silico using a dataset of 141 Whole Genome Sequences (WGS) and in a comparative analysis of 32 clinical samples. Based on cycle threshold and melting curve analysis values obtained during Real-Time PCR of nested-1 & 2 steps, we developed a simple scoring scheme and flow chart that allow identification of reaction inhibitors as well as to predict with high accuracy amplification success. The improved MLST version was used to obtain a genovars distribution in patients attending an STI clinic in Tel Aviv. Conclusions The newly developed MLST version showed great improvement of assay results for samples with very low concentrations of Chlamydia DNA. A similar concept could be applicable to other MLST schemes.

The aim of this study is the development and evaluation of a rapid and accurate quantitative PCR (qPCR)-based protocol for detection of zoonotic pathogen Streptococcus iniae in bacterial cultures and tissues of diseased fish. For this purpose, the lactate permease-encoding (lldY) gene was selected as a target for the design of S. iniae-specific primers based on comparative genomic analysis using 45 sequences retrieved from NCBI genome database. Specificity and applicability of these primers were tested using 115 bacterial strains and fish tissues infected with S. iniae. Sensitivity, reproducibility and efficiency of qPCR assay were also determined. The developed qPCR assay showed 100% specificity with pure bacterial cultures or DNA extracted from S. iniae or tissues of fish infected with the bacterium. The method has high sensitivity with a detection limit of 1.12 × 101 amplicon copies per assay (equivalent to 2 × 10-9 ng/µl) using bacterial DNA and of 1.44 × 101 gene copies in tissues of fish infected with S. iniae. In conclusion, this qPCR protocol provides an accurate and sensitive alternative for the identification of S. iniae and its detection on fish tissues that can be implemented as a routine tool in microbiological laboratories.

Background: Liriodendron chinense ranges widely in subtropical China and northern Vietnam; however, it inhabits several small, isolated populations and is now an endangered species due to its limited seed production. The objective of this study was to develop a set of nuclear SSR (simple sequence repeats) and multiple chloroplast genome markers for genetic studies in L. chinense and their characterization in diverse germplasm. Results: We performed low-coverage whole genome sequencing of the L. chinense from four genotypes, assembled the chloroplast genome and identified nuclear SSR loci by searching in contigs for SSR motifs. Comparative analysis of the four chloroplast genomes of L. chinense revealed 45 SNPs, 17 indels, 49 polymorphic SSR loci, and five small inversions. Most chloroplast intraspecific polymorphisms were located in the interspaces of single-copy regions. In total, 6147 SSR markers were isolated from low-coverage whole genome sequences. The most common SSR motifs were dinucleotide (70.09%), followed by trinucleotide motifs (23.10%). The motif AG/TC (33.51%) was the most abundant, followed by TC/AG (25.53%). A set of 13 SSR primer combinations were tested for amplification and their ability to detect polymorphisms in a set of 109 L. chinense individuals, representing distinct varieties or germplasm. The number of alleles per locus ranged from 8 to 28 with an average of 21 alleles. The expected heterozygosity (He) varied from 0.19 to 0.93 and the observed heterozygosity (Ho) ranged from 0.11 to 0.79. Conclusions: The genetic resources characterized and tested in this study provide a valuable tool to detect polymorphisms in L. chinense for future genetic studies and breeding programs.

Aerobic spore-forming Bacillales are a highly diverse and ubiquitous group that includes organisms that cause foodborne illnesses and food spoilage. Classical microbiological and biochemical identification of members of the order Bacillales represents a challenge due to the diversity of organisms in this group as well as the fact that the phenotypic-based taxonomic assignment of some named species in this group is not consistent with their phylogenomic characteristics. DNA-sequencing-based tools, on the other hand, can be fast and cost-effective, and can provide for a more reliable identification and characterization of Bacillales isolates. In comparison to 16S rDNA , rpoB was shown to better discriminate between Bacillales isolates and to allow for improved taxonomic assignment to the species level. However, the lack of a publicly accessible rpoB database, as well as the lack of standardized protocols for rpoB-based typing and strain identification, is a major challenge. Here, we report (i) the curation of a DNA sequence database for rpoB-based subtype classification of Bacillales isolates; (ii) the development of standardized protocols for generating rpoB sequence data, and a scheme for rpoB- based initial taxonomic identification of Bacillales isolates at the species level; and (iii) the integration of the database in a publicly accessible online platform that allows for the analysis of rpoB sequence data from uncharacterized Bacillales isolates. Specifically, we curated a database of DNA sequences for a 632-nt internal variable region within the rpoB gene from representative Bacillales reference type strains and a large number of isolates that we have previously isolated and characterized through multiple projects. As of May 21, 2021, the rpoB database contained more than 8350 rpoB sequences representing 1902 distinct rpoB allelic types that can be classified into 160 different genera. The database also includes 1129 rpoB sequences for representative Bacillales reference type strains as available on May 21, 2021 in the NCBI database. The rpoB database is integrated into the online Food Microbe Tracker platform ( www.foodmicrobetracker.com ) and can be queried using the integrated BLAST tool to initially subtype and taxonomically identify aerobic and facultative anaerobic spore-formers. While whole-genome sequencing is increasingly used in bacterial taxonomy , the rpoB sequence-based identification scheme described here provides a valuable tool as it allows for rapid and cost-effective initial isolate characterization, which can help to identify and characterize foodborne pathogens and food spoilage bacteria. In addition, the database and primers described here can also be adopted for metagenomics approaches that include rpoB as a target, improving discriminatory power and identification over what can be achieved using 16S rDNA as a target.

Routinely used typing methods including MLST, rep-PCR and whole genome sequencing (WGS) are time-consuming, costly, and often low throughput. Here, we describe a novel mini-MLST scheme for Eschericha coli as an alternative method for rapid genotyping. Using the proposed mini-MLST scheme, 10,946 existing STs were converted into 1,038 Melting Types (MelTs). To validate the new mini-MLST scheme, in silico analysis was performed on 73,704 strains retrieved from EnteroBase resulting in discriminatory power D = 0.9465 (CI 95% 0.9726–0.9736) for mini-MLST and D = 0.9731 (CI 95% 0.9726–0.9736) for MLST. Moreover, validation on clinical isolates was conducted with a significant concordance between MLST, rep-PCR and WGS. To conclude, the great portability, efficient processing, cost-effectiveness, and high throughput of mini-MLST represents immense benefits, even when accompanied with a slightly lower discriminatory power than other typing methods. This study proved mini-MLST is an ideal method to screen and subgroup large sets of isolates and/or quick strain typing during outbreaks. In addition, our results clearly showed its suitability for prospective surveillance monitoring of emergent and high-risk E. coli clones’.

In this study, a TaqMan-based real-time polymerase chain reaction (PCR) method to detect canine astrovirus in clinical samples was developed. Primers and probes were designed to target conserved regions of the complete viral genome sequence. The results showed that the proposed method can detect a minimum of 101 copy numbers. No cross-reactivity with other canine and feline viruses was observed. The coefficient of variation was

Fragmented and primer ligated dsRNA sequencing (FLDS) is a sequencing method applicable to long double-stranded RNA (dsRNA) that enables the complete genome sequencing of both double- and single-stranded RNA viruses. However, the application of this method on a low amount of dsRNA has been hindered by adaptor dimer formation during cDNA amplification and sequence library preparation. We herein developed FLDS ver. 3 by optimizing the terminal modification of an oligonucleotide adaptor and the conditions of adaptor ligation. We also examined the concentration of Mg2+ in the PCR reaction for cDNA amplification and the purification method of amplified cDNA. Fine sequence reads were successfully obtained from metagenomic shotgun sequencing libraries constructed from 10 and 100 pg dsRNA, and these libraries exhibited weaker detection sensitivity for low-abundance dsRNAs (viral genomes and genome segments) than that constructed from 1‍ ‍ng of dsRNA. We also report the utility of capillary electrophoresis for dsRNA quantification. The FLDS ver. 3 package expands the frontiers of our knowledge in RNA virus diversity and evolution.

The long-lasting global COVID-19 pandemic demands timely genomic investigation of SARS-CoV-2 viruses. Here, we report a simple and efficient workflow for whole-genome sequencing utilizing one-step reverse transcription-PCR (RT-PCR) amplification on a microfluidic platform, followed by MiSeq amplicon sequencing., The long-lasting global COVID-19 pandemic demands timely genomic investigation of SARS-CoV-2 viruses. Here, we report a simple and efficient workflow for whole-genome sequencing utilizing one-step reverse transcription-PCR (RT-PCR) amplification on a microfluidic platform, followed by MiSeq amplicon sequencing. The method uses Fluidigm integrated fluidic circuit (IFC) and instruments to amplify 48 samples with 39 pairs of primers, including 35 custom-designed primer pairs and four additional primer pairs from the ARTIC network protocol v3. Application of this method on RNA samples from both viral isolates and clinical specimens demonstrates robustness and efficiency in obtaining the full genome sequence of SARS-CoV-2.

Kobuvirus infections are common among humans, rodents, carnivores, pigs, and ruminants. We report herein the complete genome sequence of a novel caprine kobuvirus (MN604700) from diarrheic kids in Minnesota. Whole-genome sequencing revealed a kobuvirus genome of 8,139 nt with a single ORF region encoding a polyprotein of 2,480 amino acids. Further analysis revealed nt substitutions along the genome compared with that of the caprine kobuvirus reference strain, with 93% identity. Phylogenetic analysis indicated that the clade of the caprine kobuvirus was most closely related to porcine kobuviruses rather than bovine or ovine kobuviruses. Using primers designed from this genome, caprine kobuvirus was identified in the stools of other goats. Sanger sequencing of PCR products indicated 3D and VP1 gene nucleotides of this latter strain were 95% and 91% identical with those of MN604700, respectively. There were 35 and 101 nt substitutions in 3D and VP1 genes, respectively. Findings of kobuvirus over a 2-y period may indicate an endemic state, which needs further research. In addition, screening for kobuviruses over large geographic areas is needed to identify the evolutionary connections among different strains.

Deep knowledge of the genetic features of SARS-CoV-2 is essential to track the ongoing pandemic through different geographical areas and to design and develop early diagnostic procedures, therapeutic strategies, public health interventions, and vaccines. We describe protocols and first results of the Ion AmpliSeq&trade, SARS-CoV-2 Research Panel by a massively parallel sequencing (MPS) assay. The panel allows for targeted sequencing by overlapping amplicons, thereby providing specific, accurate, and high throughput analysis. A modified reverse transcription reaction, which consists of the use of a SARS-CoV-2 specific primers pool from the Ion AmpliSeq SARS-CoV-2 Research Panel, was assessed in order to promote viral RNA specific reverse transcription. The aim of this study was to evaluate the effectiveness of the Ion AmpliSeq&trade, SARS-CoV-2 Research Panel in sequencing the entire viral genome in different samples. SARS-CoV-2 sequence data were obtained from ten viral isolates and one nasopharyngeal swab from different patients. The ten isolate samples amplified with 12 PCR cycles displayed high mean depth values compared to those of the two isolates amplified with 20 PCR cycles. High mean depth values were also obtained for the nasopharyngeal swab processed by use of a target-specific reverse transcription. The relative depth of coverage (rDoC) analysis showed that when 12 PCR cycles were used, all target regions were amplified with high sequencing coverage, while in libraries amplified at 20 cycles, a poor uniformity of amplification, with absent or low coverage of many target regions, was observed. Our results show that the Ion AmpliSeq SARS-CoV-2 Research Panel can achieve rapid and high throughput SARS-CoV-2 whole genome sequencing from 10 ng of DNA-free viral RNA from isolates and from 1 ng of DNA-free viral RNA from a nasopharyngeal swab using 12 PCR cycles for library amplification. The modified RT-PCR protocol yielded superior results on the nasopharyngeal swab compared to the reverse transcription reaction set up according to the manufacturer&rsquo, s instructions.

Background Polyploid organisms such as wheat complicate even the simplest of procedures in molecular biology. Whilst knowledge of genomic sequences in crops is increasing rapidly, the scientific community is still a long way from producing a full pan-genome for every species. Polymerase chain reaction and Sanger sequencing therefore remain widely used as methods for characterizing gene sequences in many varieties of crops. High sequence similarity between genomes in polyploids means that if primers are not homeologue-specific via the incorporation of a SNP at the 3’ tail, sequences other than the target sequence will also be amplified. Current consensus for gene cloning in wheat is to manually perform many steps in a long bioinformatics pipeline. Results Here we present AutoCloner (www.autocloner.com), a fully automated pipeline for crop gene cloning that includes a free-to-use web interface for users. AutoCloner takes a sequence of interest from the user and performs a basic local alignment search tool (BLAST) search against the genome assembly for their particular polyploid crop. Homologous sequences are then compiled with the input sequence into a multiple sequence alignment which is mined for single-nucleotide polymorphisms (SNPs). Various combinations of potential primers that cover the entire gene of interest are then created and evaluated by Primer3; the set of primers with the highest score, as well as all possible primers at every SNP location, are then returned to the user for polymerase chain reaction (PCR). We have successfully used AutoCloner to clone various genes of interest in the Apogee wheat variety, which has no current genome sequence. In addition, we have successfully run the pipeline on ~ 80,000 high-confidence gene models from a wheat genome assembly. Conclusion AutoCloner is the first tool to fully-automate primer design for gene cloning in polyploids, where previously the consensus within the wheat community was to perform this process manually. The web interface for AutoCloner provides a simple and effective polyploid primer-design method for gene cloning, with no need for researchers to download software or input any other details other than their sequence of interest.

Funder: We thank the Human Research Tissue Bank, which is supported by the UK National Institute for Health Research (NIHR) Cambridge Biomedical Research Centre, from Addenbrooke’s Hospital. Additional infrastructure support was provided from the Cancer Research UK–funded Experimental Cancer Medicine Centre, Mutational processes acting on cancer genomes can be traced by investigating mutational signatures. Because high sequencing costs limit current studies to small numbers of good-quality samples, we propose a robust, cost- and time-effective method, called mutREAD, to detect mutational signatures from small quantities of DNA, including degraded samples. We show that mutREAD recapitulates mutational signatures identified by whole genome sequencing, and will ultimately allow the study of mutational signatures in larger cohorts and, by compatibility with formalin-fixed paraffin-embedded samples, in clinical settings.

Since December 2019, the coronavirus disease 2019 (COVID-19) caused by a novel coronavirus SARS-CoV-2 has rapidly spread to almost every nation in the world. Soon after the pandemic was recognized by epidemiologists, a group of biologists comprising the ARTIC Network, has devised a multiplexed polymerase chain reaction (PCR) protocol and primer set for targeted whole-genome amplification of SARS-CoV-2. The ARTIC primer set amplifies 98 amplicons, which are separated only in two PCRs, across a nearly entire viral genome. The original primer set and protocol showed a fairly small amplification bias when clinical samples with relatively high viral loads were used. However, as sample’s viral load become low, rapid decrease in abundances of several amplicons were seen. In this report, we will show that dimer formations between some primers are the major cause of coverage bias in the multiplex PCR. Based on this, we propose 12 alternative primers in total in the ARTIC primer set that were predicted to be involved in 14 primer interactions. The resulting primer set, version N1 (NIID-1), exhibits improved overall coverage compared to the ARTIC Network’s original (V1) and modified (V3) primer set.

Whole genome sequencing can be used to characterize and to trace viral outbreaks. Nanopore-based whole genome sequencing protocols have been described for several different viruses. These approaches utilize an overlapping amplicon-based approach which can be used to target a specific virus or group of genetically related viruses. In addition to confirmation of the virus presence, sequencing can be used for genomic epidemiology studies, to track viruses and unravel origins, reservoirs and modes of transmission. For such applications, it is crucial to understand possible effects of the error rate associated with the platform used. Routine application in clinical and public health settings require that this is documented with every important change in the protocol. Previously, a protocol for whole genome Usutu virus sequencing on the nanopore sequencing platform was validated (R9.4 flowcell) by direct comparison to Illumina sequencing. Here, we describe the method used to determine the required read coverage, using the comparison between the R10 flow cell and Illumina sequencing as an example.

Intoduction: Inherited cataract, opacification of the lens, is the most common worldwide cause of blindness in children. We aimed to identify the genetic cause of autosomal dominant (AD) posterior nuclear cataract in a four generation British family.Whole genome sequence (WGS) was performed on two affected and one unaffected individual of the family and further validated by direct sequencing. Haplotype analysis was performed via genotying.A splice-site mutation c.2826-9GA in the gene EPHA2, encoding EPH receptor A2 was identified and found to co-segregate with disease.We have identified a recurrent splice-site mutation c.2826-9GA in EPHA2 causing isolated posterior nuclear cataract, providing evidence of further phenotypic heterogeneity associated with this variant.

Objective Legionella pneumophila serogroup 1 (Lp1) sequence type (ST) 1 is globally widespread in the environment and accounts for a significant proportion of Legionella infections, including nosocomial Legionnaires' disease (LD). This study aimed to design a sensitive and specific detection method for Lp ST1 that will underpin epidemiological investigations and risk assessment. Methods A total of 628 Lp genomes (126 ST1s) were analyzed by comparative genomics. Interrogation of more than 900 accessory genes revealed seven candidate targets for specific ST1 detection and specific primers and hydrolysis probes were designed and evaluated. The analytical sensitivity and specificity of the seven primer and probe sets were evaluated on serially diluted DNA extracted from the reference strain CIP107629 and via qPCR applied on 200 characterized isolates. The diagnostic performance of the assay was evaluated on 142 culture-proven clinical samples from LD cases and a real-life investigation of a case cluster. Results Of seven qPCR assays that underwent analytical validation, one PCR target (lpp1868) showed higher sensitivity and specificity for ST1 and ST1-like strains. The diagnostic performance of the assay using respiratory samples corresponded to a sensitivity of 95% (19/20) (95% CI (75.1–99.9)) and specificity of 100% (122/122) (95% CI (97–100)). The ST1 PCR assay could link two out of three culture-negative hospitalized LD cases to ST1 during a real-time investigation. Conclusion Using whole genome sequencing (WGS) data, we developed and validated a sensitive and specific qPCR assay for the detection of Lp1 belonging to the ST1 clonal complex by amplification of the lpp1868 gene. The ST1 qPCR is expected to deliver an added value for Lp control and prevention, in conjunction with other recently developed molecular assays.

In this study, we evaluated extended-spectrum β-lactamase (ESBL)-producing bacteria with the newly developed primer and probe sets to detect bla

New antiviral drugs are needed as a first line of defense in the event of a novel influenza pandemic. Favipiravir is a broad-spectrum antiviral which is effective against influenza. The exact mechanism of how favipiravir works to inhibit influenza is still unclear. We used next-generation sequencing (NGS) to demonstrate that favipiravir causes mutations in influenza RNA. The greater depth of NGS sequence information over traditional sequencing methods allowed us to precisely determine the bias of particular mutations caused by favipiravir. NGS can also be used in a standard diagnostic pipeline to show that favipiravir is acting on the virus by revealing the mutation bias pattern typical to the drug. Our work will aid in testing whether viruses are resistant to favipiravir and may help demonstrate the effect of favipiravir on viruses in a clinical setting. This will be important if favipiravir is used during a future influenza pandemic., Favipiravir is a broad-spectrum antiviral drug that may be used to treat influenza. Previous research has identified that favipiravir likely acts as a mutagen, but the precise mutation bias that favipiravir induces in influenza virus RNAs has not been described. Here, we use next-generation sequencing (NGS) with barcoding of individual RNA molecules to accurately and quantitatively detect favipiravir-induced mutations and to sample orders of magnitude more mutations than would be possible through Sanger sequencing. We demonstrate that favipiravir causes mutations and show that favipiravir primarily acts as a guanine analogue and secondarily as an adenine analogue resulting in the accumulation of transition mutations. We also use a standard NGS pipeline to show that the mutagenic effect of favipiravir can be measured by whole-genome sequencing of virus. IMPORTANCE New antiviral drugs are needed as a first line of defense in the event of a novel influenza pandemic. Favipiravir is a broad-spectrum antiviral which is effective against influenza. The exact mechanism of how favipiravir works to inhibit influenza is still unclear. We used next-generation sequencing (NGS) to demonstrate that favipiravir causes mutations in influenza RNA. The greater depth of NGS sequence information over traditional sequencing methods allowed us to precisely determine the bias of particular mutations caused by favipiravir. NGS can also be used in a standard diagnostic pipeline to show that favipiravir is acting on the virus by revealing the mutation bias pattern typical to the drug. Our work will aid in testing whether viruses are resistant to favipiravir and may help demonstrate the effect of favipiravir on viruses in a clinical setting. This will be important if favipiravir is used during a future influenza pandemic.

African wild rice Oryza brachyantha (FF), a distant relative of cultivated rice Oryza sativa (AA), carries genes for pests and disease resistance. Molecular marker assisted alien gene introgression from this wild species to its domesticated counterpart is largely impeded due to the scarce availability of cross-transferable and polymorphic molecular markers that can clearly distinguish these two species. Availability of the whole genome sequence (WGS) of both the species provides a unique opportunity to develop markers, which are cross-transferable. We observed poor cross-transferability (~0.75 %) of O. sativa specific sequence tagged microsatellite (STMS) markers to O. brachyantha. By utilizing the genome sequence information, we developed a set of 45 low cost PCR based co-dominant polymorphic markers (STS and CAPS). These markers were found cross-transferrable (84.78 %) between the two species and could distinguish them from each other and thus allowed tracing alien genome introgression. Finally, we validated a Monosomic Alien Addition Line (MAAL) carrying chromosome 1 of O. brachyantha in O. sativa background using these markers, as a proof of concept. Hence, in this study, we have identified a set molecular marker (comprising of STMS, STS and CAPS) that are capable of detecting alien genome introgression from O. brachyantha to O. sativa.

Highlights • The whole genome sequence of feline morbillivirus (FeMV) Piuma/2015 has been obtained by the combination of SISPA and NGS starting from the infected urine sample of a cat suffering of chronic kidney disease. • The highest sequence identity was detected with early FeMVs from Hong Kong. • Sequence heterogeneity exists within European FeMVs as for the existence of FeMVs in Germany and Turkey divergent from Piuma/2015., Feline morbillivirus (FeMV) has been recently identified by RT-PCR in the urine sample of a nephropathic cat in Italy. In this report, we describe the whole genome sequence of strain Piuma/2015 obtained by combination of sequence independent single primer amplification method (SISPA) and next generation sequencing (NGS) starting from RNA purified from the infected urine sample. The existence in Germany and Turkey of FeMVs from cats divergent from Piuma/2015, suggests the presence of FeMV heterogeneity in Europe as it has been described previously in Japan and China.

Background Viral infection by dengue virus is a major public health problem in tropical countries. Early diagnosis and detection are increasingly based on quantitative reverse transcriptase real-time polymerase chain reaction (RT-qPCR) directed against genomic regions conserved between different isolates. Genetic variation can however result in mismatches of primers and probes with their targeted nucleic acid regions. Whole genome sequencing allows to characterize and track such changes, which in turn enables to evaluate, optimize, and (re-)design novel and existing RT-qPCR methods. The immense amount of available sequence data renders this however a labour-intensive and complex task. Results We present a bioinformatics approach that enables in silico evaluation of primers and probes intended for routinely employed RT-qPCR methods. This approach is based on analysing large amounts of publically available whole genome data, by first employing BLASTN to mine the genomic regions targeted by the RT-qPCR method(s), and afterwards using BLASTN-SHORT to evaluate whether primers and probes will anneal based on a set of simple in silico criteria. Using dengue virus as a case study, we evaluated 18 published RT-qPCR methods using more than 3000 publically available genomes in the NCBI Virus Variation Resource, and provide a systematic overview of method performance based on in silico sensitivity and specificity. Conclusions We provide a comprehensive overview of dengue virus RT-qPCR method performance that will aid appropriate method selection allowing to take specific measures that aim to contain and prevent viral spread in afflicted regions. Notably, we find that primer-template mismatches at their 3′ end may represent a general issue for dengue virus RT-qPCR detection methods that merits more attention in their development process. Our approach is also available as a public tool, and demonstrates how utilizing genomic data can provide meaningful insights in an applied public health setting such as the detection of viral species in human diagnostics. Electronic supplementary material The online version of this article (10.1186/s12859-018-2313-0) contains supplementary material, which is available to authorized users.

The genus Escherichia is composed of Escherichia albertii, E. fergusonii, five cryptic Escherichia clades and E. coli sensu stricto. Furthermore, the E. coli species can be divided into seven main phylogroups termed A, B1, B2, C, D, E and F. As specific lifestyles and/or hosts can be attributed to these species/phylogroups, their identification is meaningful for epidemiological studies. Classical phenotypic tests fail to identify non-sensu stricto E. coli as well as phylogroups. Clermont and colleagues have developed PCR assays that allow the identification of most of these species/phylogroups, the triplex/quadruplex PCR for E. coli phylogroup determination being the most popular. With the growing availability of whole genome sequences, we have developed the ClermonTyping method and its associated web-interface, the ClermonTyper, that allows a given strain sequence to be assigned to E. albertii, E. fergusonii, Escherichia clades I–V, E. coli sensu stricto as well as to the seven main E. coli phylogroups. The ClermonTyping is based on the concept of in vitro PCR assays and maintains the principles of ease of use and speed that prevailed during the development of the in vitro assays. This in silico approach shows 99.4 % concordance with the in vitro PCR assays and 98.8 % with the Mash genome-clustering tool. The very few discrepancies result from various errors occurring mainly from horizontal gene transfers or SNPs in the primers. We propose the ClermonTyper as a freely available resource to the scientific community at: http://clermontyping.iame-research.center/.

Highlights • Identification of two new serovars of Actinobacillus pleuropneumoniae. • Serological confirmation of specific reactivity with homologous antisera. • Characterization of the capsule loci of serovars 17 and 18. • Development of PCRs for molecular diagnostics., The aim of this study was to investigate isolates of Actinobacillus pleuropneumoniae previously designated serologically either as non-typable (NT) or as ‘K2:07’, which did not produce serovar-specific amplicons in PCR assays. We used whole genome sequencing to identify the capsule (CPS) loci of six previously designated biovar 1 NT and two biovar 1 ‘K2:O7’ isolates of A. pleuropneumoniae from Denmark, as well as a recent biovar 2 NT isolate from Canada. All of the NT isolates have the same six-gene type I CPS locus, sharing common cpsABC genes with serovars 2, 3, 6, 7, 8, 9, 11 and 13. The two ‘K2:O7’ isolates contain a unique three-gene type II CPS locus, having a cpsA gene similar to that of serovars 1, 4, 12, 14 and 15. The previously NT isolates share the same O-antigen genes, found between erpA and rpsU, as serovars 3, 6, 8, and 15. Whereas the ‘K2:O7’ isolates, have the same O-antigen genes as serovar 7, which likely contributed to their previous mis-identification. All of the NT and ‘K2:O7’ isolates have only the genes required for production of ApxII (apxIICA structural genes, and apxIBD export genes). Rabbit polyclonal antisera raised against representative isolates with these new CPS loci demonstrated distinct reactivity compared to the 16 known serovars. The serological and genomic results indicate that the isolates constitute new serovars 17 (previously NT) and 18 (previously ‘K2:O7’). Primers designed for amplification of specific serovar 17 and 18 sequences for molecular diagnostics will facilitate epidemiological tracking of these two new serovars of A. pleuropneumoniae.

Genomics has greatly transformed our understanding of phage biology; however, traditional methods of DNA isolation for whole genome sequencing have required phages to be grown to high titers in large-scale preparations, potentially selecting for only those phages that can grow efficiently under laboratory conditions. This may also select for mutations or deletions that enable more efficient growth in culture. The ability to sequence a bacteriophage genome from a single isolated plaque reduces these risks while decreasing the time and complexity of bacteriophage genome sequencing. A method of amplification and library preparation is described, utilizing Sequence Independent Single Primer Amplification (SISPA), that can be used for whole genome shotgun sequencing of bacteriophages from a single isolated plaque.

Coronaviruses are of major importance for both animal and human health. With the emergence of novel coronaviruses such as SARS and MERS, the need for fast genome characterisation is ever so important. Further, in order to understand the influence of quasispecies of these viruses in relation to biology, techniques for deep-sequence and full-length viral genome analysis are needed. In the present study, we compared the efficiency of two sequence-independent approaches [sequence-independent single primer amplification (SISPA) and single primer isothermal amplification (SPIA, represented by the Ovation kit)] coupled with high-throughput sequencing to generate the full-length genome of bovine coronavirus (BCoV) from a nasal swab. Both methods achieved high genome coverage (100% for SPIA and 99% for SISPA), however, there was a clear difference in the percentage of reads that mapped to BCoV. While approximately 45% of the Ovation reads mapped to BCoV (sequence depth of 169-284 944), only 0.07% of the SISPA reads (sequence depth of 0-249) mapped to the reference genome. Although BCoV was the focus of the study we also identified a bovine rhinitis B virus (BRBV) in the data sets. The trend for this virus was similar to that observed for BCoV regarding Ovation vs. SISPA, but with fewer sequences mapping to BRBV due to a lower amount of this virus. In summary, the SPIA approach used in this study produced coverage of the entire BCoV (high copy number) and BRBV (low copy number) and a high sequence/genome depth compared to SISPA. Although this is a limited study, the results indicate that the Ovation method could be a preferred approach for full genome sequencing if a low copy number of viral RNA is expected and if high sequence depth is desired.

Porcine kobuvirus has been detected from pig fecal samples in the USA, but there is still no information on the full-length genomes. In this study, we characterized the first complete genomic sequence of a US porcine kobuvirus strain OH/RV50/2011. The viral genome is 8123 nucleotides (nt) long, including a 576-nt 5'-untranslated region (UTR), a 7380-nt polyprotein encoding sequence, and a 167-nt 3'-UTR. A complete genome sequence alignment suggested that two types of porcine kobuviruses were found based on whether a 30-aa deletion existed in the 2B encoding region. Furthermore, several conserved motifs that can be used for the design of universal kobuvirus or porcine kobuvirus-specific primers were verified in non-structural protein genes. Phylogenetic analysis based on the complete genome sequence showed that RV50 was grouped with other porcine kobuviruses and more closely related to Chinese strains. Secondary structure analysis of the 5'-UTR showed that RV50 has three stem-loop domains in the first 108 nt and has a potential hepacivirus-/pestivirus-like type IV group-B-like internal ribosomal entry site, like the porcine kobuvirus prototype strain S-1. Codon usage analysis showed that the most preferred usage tends to be C or U at the end of a codon in a porcine kobuvirus genome. These results will be useful in understanding the evolution of porcine kobuviruses .

We have developed and evaluated a real-time reverse transcriptase PCR (RT-PCR) assay for the detection of human enterovirus D68 (EV-D68) in clinical specimens. This assay was developed in response to the unprecedented 2014 nationwide EV-D68 outbreak in the United States associated with severe respiratory illness. As part of our evaluation of the outbreak, we sequenced and published the genome sequence of the EV-D68 virus circulating in St. Louis, MO. This sequence, along with other GenBank sequences from past EV-D68 occurrences, was used to computationally select a region of EV-D68 appropriate for targeting in a strain-specific RT-PCR assay. The RT-PCR assay amplifies a segment of the VP1 gene, with an analytic limit of detection of 4 copies per reaction, and it was more sensitive than commercially available assays that detect enteroviruses and rhinoviruses without distinguishing between the two, including three multiplex respiratory panels approved for clinical use by the FDA. The assay did not detect any other enteroviruses or rhinoviruses tested and did detect divergent strains of EV-D68, including the first EV-D68 strain (Fermon) identified in California in 1962. This assay should be useful for identifying and studying current and future outbreaks of EV-D68 viruses.

Diverse members of the genus Clostridium produce botulinum neurotoxins (BoNTs), which cause a flaccid paralysis known as botulism. While multiple species of clostridia produce BoNTs, the majority of human botulism cases have been attributed to Clostridium botulinum groups I and II. Recent comparative genomic studies have demonstrated the genomic diversity within these BoNT-producing species. This report introduces a multiplex PCR assay for differentiating members of C. botulinum group I, C. sporogenes , and two major subgroups within C. botulinum group II. Coding region sequences unique to each of the four species/subgroups were identified by in silico analyses of thousands of genome assemblies, and PCR primers were designed to amplify each marker. The resulting multiplex PCR assay correctly assigned 41 tested isolates to the appropriate species or subgroup. A separate PCR assay to determine the presence of the ntnh gene (a gene associated with the botulinum neurotoxin gene cluster) was developed and validated. The ntnh gene PCR assay provides information about the presence or absence of the botulinum neurotoxin gene cluster and the type of gene cluster present ( ha positive [ ha + ] or orfX + ). The increased availability of whole-genome sequence data and comparative genomic tools enabled the design of these assays, which provide valuable information for characterizing BoNT-producing clostridia. The PCR assays are rapid, inexpensive tests that can be applied to a variety of sample types to assign isolates to species/subgroups and to detect clostridia with botulinum neurotoxin gene ( bont ) clusters. IMPORTANCE Diverse clostridia produce the botulinum neurotoxin, one of the most potent known neurotoxins. In this study, a multiplex PCR assay was developed to differentiate clostridia that are most commonly isolated in connection with human botulism cases: C. botulinum group I, C. sporogenes , and two major subgroups within C. botulinum group II. Since BoNT-producing and nontoxigenic isolates can be found in each species, a PCR assay to determine the presence of the ntnh gene, which is a universally present component of bont gene clusters, and to provide information about the type ( ha + or orfX + ) of bont gene cluster present in a sample was also developed. The PCR assays provide simple, rapid, and inexpensive tools for screening uncharacterized isolates from clinical or environmental samples. The information provided by these assays can inform epidemiological studies, aid with identifying mixtures of isolates and unknown isolates in culture collections, and confirm the presence of bacteria of interest.

Background Although the sequencing landscape is rapidly evolving and sequencing costs are continuously decreasing, whole genome sequencing is still too expensive for use on a routine basis. Targeted resequencing of only the regions of interest decreases both costs and the complexity of the downstream data-analysis. Various target enrichment strategies are available, but none of them obtain the degree of coverage uniformity, flexibility and specificity of PCR-based enrichment. On the other hand, the biggest limitation of target enrichment by PCR is the need to design large numbers of partially overlapping assays to cover the target. Results To overcome the aforementioned hurdles, we have developed primerXL, a state-of-the-art PCR primer design pipeline for targeted resequencing. It uses an optimized design criteria relaxation cascade and a thorough downstream in silico evaluation process to generate high quality singleplex PCR assays, reducing the need for amplicon normalization, and outperforming other target enrichment strategies and similar primer design tools when considering assay quality, coverage uniformity and target coverage. Results of four different sequencing projects with 2348 amplicons in total covering 470 kb are presented. PrimerXL can be accessed at www.primerxl.org. Conclusion PrimerXL is an state-of-the-art, easy to use primer design webtool capable of generating high-quality targeted resequencing assays. The workflow is fully customizable to suit every researchers’ needs, while an innovative relaxation cascade ensures maximal target coverage. Electronic supplementary material The online version of this article doi:(10.1186/s12859-017-1809-3) contains supplementary material, which is available to authorized users.

Hepatitis B virus genotype F (HBV/F) is endemic in Central and South America with a minor proportion in Mexico and North America. HBV/F is divided into subgenotypes and subtypes with particular geographic circulation patterns. Here, we report the complete genome sequence and molecular characterization of HBV/F from three isolates. Phylogenetic analysis with all available HBV/F sequences showed that our sequences belonged to the F1b subtype and, in addition, the absence of the previously reported F1a subtype in Mexican isolates. Our findings suggest the circulation of HBV/F1b, the first phylogenomic study of HBV/F in Mexico.

We have developed methods for full-genome sequencing of Zika viruses (ZIKVs) based on a targeted amplification approach. We used alignments of publicly available complete genome data to design a primer set that selectively amplifies ZIKVs. The approach includes amplification strategies for templates present at both high- and low-copy number, and PCR cycling conditions that have been normalized across genome fragments in order to streamline laboratory handling. Abundant templates can be amplified using a strategy that uses 6 overlapping amplicons to cover the complete viral genome, whereas scarce templates can be amplified using a strategy that uses 11 overlapping amplicons of smaller size. The workflow is sequencing platform agnostic, and thus, can be used in low resource settings where access to traditional Sanger sequencing is the only option available. Given the scarcity of tools for ZIKV, this approach should facilitate epidemiological surveillance and other studies that require the generation of complete viral genomic information quickly and cost-effectively.

Mimosa scabrella Benth., popularly known as ''bracatinga'', is a pioneer and endemic species of Brazil, occurring in Mixed Ombrophilous Forest associated with Brazilian Atlantic Rainforest biomes. It is a fast-growing tree of the Fabaceae family that facilitates the dynamics of ecological succession. SSR development, when there is no genome sequence, is time and labor intensive and there are no molecular markers for M. scabrella. We developed and validated the first microsatellite markers for this tetraploid species, evaluating mother trees and progenies. Using Illumina sequencing, we identified 290 SSR loci and 211 primer pairs. After 31 SSR loci PCR/agarose electrophoresis selection, a subset of 11 primer pairs was synthetized with fluorescence in the forward primer for PCR and capillary electrophoresis validation with leaf DNA of 33 adult and 411 progeny individuals. Polymorphic locus percentage was 36, 4 in 11 loci, 3 chloroplast SSRs, and 1 nuclear SSR. Allele number of polymorphic loci ranged from 2 to 11 alleles considering all sampling. All 11 primer pairs were also tested for cross-species amplification for five Fabaceae-Mimosoideae species, ranging from 2 loci transferred to Calliandra tweedii Benth. and all 11 loci transferred to Mimosa taimbensis Burkart. The assessed and validated SSR markers for M. scabrella are suitable and useful for analysis and population genetic studies.

Whole genome sequencing of six diagnostic brucellaphages, Tbilisi (Tb), Firenze (Fz), Weybridge (Wb), S708, Berkeley (Bk) and R/C, was followed with genomic comparisons including recently described genomes of the Tb phage from Mexico (TbM) and Pr phage to elucidate genomic diversity and candidate host range determinants. Comparative whole genome analysis revealed high sequence homogeneity among these brucellaphage genomes and resolved three genetic groups consistent with defined host range phenotypes. Group I was composed of Tb and Fz phages that are predominantly lytic for Brucella abortus and Brucella neotomae; Group II included Bk, R/C, and Pr phages that are lytic mainly for B. abortus, Brucella melitensis and Brucella suis; Group III was composed of Wb and S708 phages that are lytic for B. suis, B. abortus and B. neotomae. We found that the putative phage collar protein is a variable locus with features that may be contributing to the host specificities exhibited by different brucellaphage groups. The presence of several candidate host range determinants is illustrated herein for future dissection of the differential host specificity observed among these phages.