Automating microsatellite screening and primer design from multi-individual libraries using Micro-Primers

Analysis of intra- and inter-population diversity has become important for defining the genetic status and distribution patterns of a species and a powerful tool for conservation programs, as high levels of inbreeding could lead into whole population extinction in few generations. Microsatellites (SSR) are commonly used in population studies but discovering highly variable regions across species’ genomes requires demanding computation and laboratorial optimization. In this work, we combine next generation sequencing (NGS) with automatic computing to develop a genomic-oriented tool for characterizing SSRs at the population level. Herein, we describe a new Python pipeline, named Micro-Primers, designed to identify, and design PCR primers for amplification of SSR loci from a multi-individual microsatellite library. By combining commonly used programs for data cleaning and microsatellite mining, this pipeline easily generates, from a fastq file produced by high-throughput sequencing, standard information about the selected microsatellite loci, including the number of alleles in the population subset, and the melting temperature and respective PCR product of each primer set. Additionally, potential polymorphic loci can be identified based on the allele ranges observed in the population, to easily guide the selection of optimal markers for the species. Experimental results show that Micro-Primers significantly reduces processing time in comparison to manual analysis while keeping the same quality of the results. The elapsed times at each step can be longer depending on the number of sequences to analyze and, if not assisted, the selection of polymorphic loci from multiple individuals can represent a major bottleneck in population studies.

LinearTurboFold: Linear-time global prediction of conserved structures for RNA homologs with applications to SARS-CoV-2

The constant emergence of COVID-19 variants reduces the effectiveness of existing vaccines and test kits. Therefore, it is critical to identify conserved structures in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes as potential targets for variant-proof diagnostics and therapeutics. However, the algorithms to predict these conserved structures, which simultaneously fold and align multiple RNA homologs, scale at best cubically with sequence length and are thus infeasible for coronaviruses, which possess the longest genomes (∼30,000 nt) among RNA viruses. As a result, existing efforts on modeling SARS-CoV-2 structures resort to single-sequence folding as well as local folding methods with short window sizes, which inevitably neglect long-range interactions that are crucial in RNA functions. Here we present LinearTurboFold, an efficient algorithm for folding RNA homologs that scales linearly with sequence length, enabling unprecedented global structural analysis on SARS-CoV-2. Surprisingly, on a group of SARS-CoV-2 and SARS-related genomes, LinearTurboFold’s purely in silico prediction not only is close to experimentally guided models for local structures, but also goes far beyond them by capturing the end-to-end pairs between 5′ and 3′ untranslated regions (UTRs) (∼29,800 nt apart) that match perfectly with a purely experimental work. Furthermore, LinearTurboFold identifies undiscovered conserved structures and conserved accessible regions as potential targets for designing efficient and mutation-insensitive small-molecule drugs, antisense oligonucleotides, small interfering RNAs (siRNAs), CRISPR-Cas13 guide RNAs, and RT-PCR primers. LinearTurboFold is a general technique that can also be applied to other RNA viruses and full-length genome studies and will be a useful tool in fighting the current and future pandemics.

Ecological surveillance of bat coronaviruses in Sarawak, Malaysian Borneo

Objective Coronaviruses (CoVs) are natural commensals of bats. Two subgenera, namely Sarbecoviruses and Merbecoviruses have a high zoonotic potential and have been associated with three separate spillover events in the past 2 decades, making surveillance of bat-CoVs crucial for the prevention of the next epidemic. The study was aimed to elucidate the presence of coronavirus in fresh bat guano sampled from Wind Cave Nature Reserve (WCNR) in Sarawak, Malaysian Borneo. Samples collected were placed into viral transport medium, transported on ice within the collection day, and preserved at − 80 °C. Nucleic acid was extracted using the column method and screened using consensus PCR primers targeting the RNA-dependent RNA polymerase (RdRp) gene. Amplicons were sequenced bidirectionally using the Sanger method. Phylogenetic tree with maximum-likelihood bootstrap and Bayesian posterior probability were constructed. Results CoV-RNA was detected in ten specimens (47.6%, n = 21). Six alphacoronavirus and four betacoronaviruses were identified. The bat-CoVs can be phylogenetically grouped into four novel clades which are closely related to Decacovirus-1 and Decacovirus-2, Sarbecovirus, and an unclassified CoV. CoVs lineages unique to the Island of Borneo were discovered in Sarawak, Malaysia, with one of them closely related to Sarbecovirus. All of them are distant from currently known human coronaviruses.

More DNA and RNA of HBV SP1 splice variants are detected in genotypes B and C at low viral replication

HBV produces unspliced and spliced RNAs during replication. Encapsidated spliced RNA is converted into DNA generating defective virions that are detected in plasma and associated with HCC development. Herein we describe a quantitative real-time PCR detection of splice variant SP1 DNA/RNA in HBV plasma. Three PCR primers/probe sets were designed detecting the SP1 variants, unspliced core, or X gene. Plasmids carrying the three regions were constructed for the nine HBV genotypes to evaluate the three sets, which were also tested on DNA/RNA extracted from 193 HBV plasma with unknown HCC status. The assay had an LOD of 80 copies/ml and was equally efficient for detecting all nine genotypes and three targets. In testing 84 specimens for both SP1 DNA (77.4%) and RNA (82.1%), higher viral loads resulted in increased SP1 levels. Most samples yielded < 1% of SP1 DNA, while the average SP1 RNA was 3.29%. At viral load of ≤ 5 log copies/ml, the detectable SP1 DNA varied by genotype, with 70% for B, 33.3% for C, 10.5% for E, 4% for D and 0% for A, suggesting higher levels of splicing in B and C during low replication. At > 5 log, all samples regardless of genotype had detectable SP1 DNA.

Development of Multiplex PCR Coupled DNA Chip Technology for Assessment of Endogenous and Exogenous Allergens in GM Soybean

Allergenicity assessment of transgenic plants and foods is important for food safety, labeling regulations, and health protection. The aim of this study was to develop an effective multi-allergen diagnostic approach for transgenic soybean assessment. For this purpose, multiplex polymerase chain reaction (PCR) coupled with DNA chip technology was employed. The study was focused on the herbicide-resistant Roundup Ready soya (RRS) using a set of certified reference materials consisting of 0, 0.1%, 0.5%, and 10% RRS. Technically, the procedure included design of PCR primers and probes; genomic DNA extraction; development of uniplex and multiplex PCR systems; DNA analysis by agarose gel electrophoresis; microarray development, hybridization, and scanning. The use of the asymmetric multiplex PCR method is shown to be very efficient for DNA hybridization with biochip probes. We demonstrate that newly developed fourplex PCR methods coupled with DNA-biochips enable simultaneous identification of three major endogenous allergens, namely, Gly m Bd 28K, Gly m Bd 30K, and lectin, as well as exogenous 5-enolppyruvyl shikimate-phosphate synthase (epsps) expressed in herbicide-resistant roundup ready GMOs. The approach developed in this study can be used for accurate, cheap, and fast testing of food allergens.

Use of universal primers for the 18S ribosomal RNA gene and whole soil DNAs to reveal the taxonomic structures of soil nematodes by high-throughput amplicon sequencing

Nematodes are abundant metazoans that play crucial roles in nutrient recycle in the pedosphere. Although high-throughput amplicon sequencing is a powerful tool for the taxonomic profiling of soil nematodes, polymerase chain reaction (PCR) primers for amplification of the 18S ribosomal RNA (SSU) gene and preparation of template DNAs have not been sufficiently evaluated. We investigated nematode community structure in copse soil using four nematode-specific (regions 1–4) and two universal (regions U1 and U2) primer sets for the SSU gene regions with two DNAs prepared from copse-derived mixed nematodes and whole soil. The major nematode-derived sequence variants (SVs) identified in each region was detected in both template DNAs. Order level taxonomy and feeding type of identified nematode-derived SVs were distantly related between the two DNA preparations, and the region U2 was closely related to region 4 in the non-metric multidimensional scaling (NMDS) based on Bray-Curtis dissimilarity. Thus, the universal primers for region U2 could be used to analyze soil nematode communities. We further applied this method to analyze the nematodes living in two sampling sites of a sweet potato-cultivated field, where the plants were differently growing. The structure of nematode-derived SVs from the two sites was distantly related in the principal coordinate analysis (PCoA) with weighted unifrac distances, suggesting their distinct soil environments. The resultant ecophysiological status of the nematode communities in the copse and field on the basis of feeding behavior and maturity indices was fairly consistent with those of the copse- and the cultivated house garden-derived nematodes in prior studies. These findings will be useful for the DNA metabarcoding of soil eukaryotes, including nematodes, using soil DNAs.

Shine: To explore specific, sensitive and conserved biomarkers from massive microbial genomic data within intrapopulations.

To improve the quality of nucleic acid detection reagents, we provided a new strategy, Shine, to explore specific, sensitive and conserved biomarkers from massive microbial genomic data within intrapopulations in order to improve detection sensitivity and accuracy. It is obvious that the more comprehensive genomic data are, the more effective the detection biomarkers. Here, we demonstrated that our method could detect undiscovered multicopy conserved species-specific or even subspecies-specific target fragments, according to several clinical projects. In particular, this approach was effective for any pathogenic microorganism even in incompletely assembled motifs. Based on our strategy, the detection device designed with quantitative PCR primers and probes for systematic and automated detection of pathogenic microorganisms in biological samples may cover all pathogenic microorganisms without limits based on genome annotation. On the website https://bioinfo.liferiver.com.cn, users may select different configuration parameters depending on the purpose of the project to realize routine clinical detection practices. Therefore, it is recommended that our strategy is suitable to identify shared universal phylogenetic markers with few false positive or false negative errors and to automate the design of minimal primers and probes to detect pathogenic communities with cost-effective predictive power.

The invasive Korean bush mosquito Aedes koreicus (Diptera: Culicidae) in Germany as of 2020

Background The Korean bush mosquito Aedes koreicus was recently reported to have established a population in western Germany (Wiesbaden) in 2016. The species is difficult to distinguish morphologically from its close relative, the invasive Japanese bush mosquito Ae. japonicus, which is already widely distributed in many parts of Germany, including the area colonised by Ae. koreicus. Genetic confirmation of morphologically identified “Ae. japonicus” collection material, however, had only been done exceptionally before the German Ae. koreicus population became known. Methods Dried archived “Ae. japonicus” specimens both from the municipality of Wiesbaden and from deliberately and randomly selected distribution sites all over Germany were re-examined morphologically and genetically for admixture by Ae. koreicus. Moreover, cemeteries in the greater Wiesbaden area were sampled in 2019 and 2020 to check for Ae. koreicus spread. Korean and Japanese bush mosquitoes submitted to the German citizen science mosquito monitoring scheme “Mueckenatlas” in 2019 and 2020 were also subjected to particularly thorough species identification. The ND4 DNA sequences generated in this study in the context of species identification were phylogenetically compared to respective GenBank entries of Ae. koreicus. As a by-product, several genetic markers were evaluated for their suitability to identify Ae. koreicus. Results Aedes koreicus specimens could be identified in mosquito collection material and submissions from Wiesbaden from 2015 onwards, suggesting establishment to have happened in the same year as Ae. japonicus establishment. Detections of Ae. koreicus from 2019 and 2020 in Wiesbaden indicate a negligible enlargement of the populated area as described for 2018. Two Ae. koreicus specimens were also submitted from the city of Munich, southern Germany, in 2019 but further specimens could not be identified during immediate local inspections. Comparison of ND4 sequences generated in this and other studies demonstrate a high degree of homology, suggesting that this DNA region is not informative enough for clarification of origins and relationships of Ae. koreicus populations. For genetic identification of Ae. koreicus, PCR primers used for classical CO1 barcoding were found to lead to mismatches and produce no or incorrect amplicons. Alternative CO1 primers or a validated ND4 marker should be used instead. Conclusions Aedes koreicus is probably introduced into Germany every now and then but rarely succeeds in becoming established. As with most European populations, the German population is characterised by a limited expansion tendency. Since Ae. koreicus is a potential vector, however, Asian bush mosquitoes found at new places should be examined quite carefully and known distribution areas of Ae. japonicus regularly checked for the presence of Ae. koreicus. Graphical Abstract

Development of novel PCR primer sets for DNA metabarcoding of aquatic insects, and the discovery of some cryptic species

DNA barcoding is a powerful tool that provides rapid, accurate, and automatable species identification by using standardized genetic region(s). It can be a powerful tool in various fields of biology such as for revealing the existence of cryptic species and/or rare species and in environmental science such as when monitoring river biota. Biodiversity reduction in recent times has become one of the most serious environmental issues on a worldwide scale. DNA barcoding techniques require the development of sets of universal PCR primers for DNA metabarcoding. We tried to develop universal primer sets for the DNA barcoding of all insect groups. In this study, we succeeded in designing not only universal primer sets for DNA barcoding regions of almost all insects, which were designed to include a hypervariable site between highly conserved sites, but also primer sets for longer fragment sequences for registration in a database. We confirmed successful amplification for 14 orders, 43 families, and 68 species with DNA barcoding in the mtDNA 16S rRNA region, and for 13 orders, 42 families, and 66 species with DNA barcoding in the mtDNA 12S rRNA region. A key feature is that the DNA fragments of the DNA barcoding regions amplified by these primer sets are both short at about 200-bp, and longer fragment sequences will increase the level of data registration in the DNA database. Such resulting database enhancements will serve as a powerful tool for increasingly accurate assessment of biodiversity and genetic diversity.