ThermoAlign: a genome-aware primer design tool for tiled amplicon resequencing

primerJinn: a tool for rationally designing multiplex PCR primer sets for amplicon sequencing and performing in silico PCR

BACKGROUND

Multiplex PCR amplifies numerous targets in a single tube reaction and is essential in molecular biology and clinical diagnostics. One of its most important applications is in the targeted sequencing of pathogens. Despite this importance, few tools are available for designing multiplex primers.

RESULTS

We developed primerJinn, a tool that designs a set of multiplex primers and allows for the in silico PCR evaluation of primer sets against numerous input genomes. We used primerJinn to create a multiplex PCR for the sequencing of drug resistance-conferring gene regions from Mycobacterium tuberculosis, which were then successfully sequenced.

CONCLUSIONS

primerJinn provides a user-friendly, efficient, and accurate method for designing multiplex PCR primers for targeted sequencing and performing in silico PCR. It can be used for various applications in molecular biology and bioinformatics research, including the design of assays for amplifying and sequencing drug-resistance-conferring regions in important pathogens.

primerJinn - a tool for rationally designing multiplex PCR primer sets and in silico PCR

Background

Multiplex PCR amplifies numerous targets in a single tube reaction and is essential in molecular biology and clinical diagnostics. One of its most important applications is in the targeted sequencing of pathogens. Despite this importance, few tools are available for designing multiplex primers.

Results

We developed primerJinn, a tool that designs a set of multiplex primers and allows for the in silico PCR evaluation of primer sets against numerous input genomes. We used primerJinn to create a multiplex PCR for the sequencing of drug resistance-conferring gene regions from Mycobacterium tuberculosis, which were then successfully sequenced.

Conclusions

primerJinn provides a user-friendly, efficient, and accurate method for designing multiplex PCR primers and performing in silico PCR. It can be used for various applications in molecular biology and bioinformatics research, including the design of assays for amplifying and sequencing drug-resistance-conferring regions in important pathogens.

Long amplicon HiFi sequencing for mitochondrial DNA genomes

Long-read sequencing technology is a powerful approach with application in various genetic and genomic research. Herein, we developed the pipeline for long amplicon high-fidelity (HiFi) sequencing and then applied it for sequencing mitochondrial DNA (mtDNA) genomes from pools of 79 Tibetan Mastiffs. We amplified the mtDNA genome with long-range PCR using two pairs of primers. Two rounds of circular consensus sequencing (CCS) were conducted and their accuracy was evaluated. The results indicate that the second round of CCS can improve the accuracy of HiFi reads. In addition, the analysis of 79 high-quality mtDNA genomes shows the Tibetan Mastiffs from outside of the Tibetan Plateau experienced hybridization with other dogs. The high quality reads generator (HQGR) software is provided to facilitate data analyses, which is publicly accessible on GitHub (https://github.com/Caizf-script/HQGR). Our long amplicon HiFi sequencing pipeline can also be applied in various target enrichment strategies for small genomes and candidate genes.

The web-based multiplex PCR primer design software Ultiplex and the associated experimental workflow: up to 100- plex multiplicity

Background

A large number of variants have been employed in various medical applications, such as providing medication instructions, disease susceptibility testing, paternity testing, and tumour diagnosis. A high multiplicity PCR will outperform other technologies because of its lower cost, reaction time and sample consumption. To conduct a multiplex PCR with higher than 100 plex multiplicity, primers need to be carefully designed to avoid the formation of secondary structures and nonspecific amplification between primers, templates and products. Thus, a user-friendly, highly automated and highly user-defined web-based multiplex PCR primer design software is needed to minimize the work of primer design and experimental verification.

Results

Ultiplex was developed as a free online multiplex primer design tool with a user-friendly web-based interface (http://ultiplex.igenebook.cn). To evaluate the performance of Ultiplex, 294 out of 295 (99.7%) target primers were successfully designed. A total of 275 targets produced qualified primers after primer filtration, and 271 of those targets were successfully clustered into one compatible PCR group and could be covered by 108 primers. The designed primer group stably detected the rs28934573(C > T) mutation at lower than a 0.25% mutation rate in a series of samples with different ratios of HCT-15 and HaCaT cell line DNA.

Conclusion

Ultiplex is a web-based multiplex PCR primer tool that has several functions, including batch design and compatibility checking for the exclusion of mutual secondary structures and mutual false alignments across the whole genome. It offers flexible arguments for users to define their own references, primer Tm values, product lengths, plex numbers and tag oligos. With its user-friendly reports and web-based interface, Ultiplex will provide assistance for biological applications and research involving genomic variants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08149-1.

Comparative analysis reveals within-population genome size variation in a rotifer is driven by large genomic elements with highly abundant satellite DNA repeat elements

BACKGROUND

Eukaryotic genomes are known to display an enormous variation in size, but the evolutionary causes of this phenomenon are still poorly understood. To obtain mechanistic insights into such variation, previous studies have often employed comparative genomics approaches involving closely related species or geographically isolated populations within a species. Genome comparisons among individuals of the same population remained so far understudied-despite their great potential in providing a microevolutionary perspective to genome size evolution. The rotifer Brachionus asplanchnoidis represents one of the most extreme cases of within-population genome size variation among eukaryotes, displaying almost twofold variation within a geographic population.

RESULTS

Here, we used a whole-genome sequencing approach to identify the underlying DNA sequence differences by assembling a high-quality reference genome draft for one individual of the population and aligning short reads of 15 individuals from the same geographic population including the reference individual. We identified several large, contiguous copy number variable regions (CNVs), up to megabases in size, which exhibited striking coverage differences among individuals, and whose coverage overall scaled with genome size. CNVs were of remarkably low complexity, being mainly composed of tandemly repeated satellite DNA with only a few interspersed genes or other sequences, and were characterized by a significantly elevated GC-content. CNV patterns in offspring of two parents with divergent genome size and CNV patterns in several individuals from an inbred line differing in genome size demonstrated inheritance and accumulation of CNVs across generations.

CONCLUSIONS

By identifying the exact genomic elements that cause within-population genome size variation, our study paves the way for studying genome size evolution in contemporary populations rather than inferring patterns and processes a posteriori from species comparisons.

AnthOligo: automating the design of oligonucleotides for capture/enrichment technologies

SUMMARY

A number of methods have been devised to address the need for targeted genomic resequencing. One of these methods, region-specific extraction (RSE) is characterized by the capture of long DNA fragments (15-20 kb) by magnetic beads, after enzymatic extension of oligonucleotides hybridized to selected genomic regions. Facilitating the selection of the most appropriate capture oligos for targeting a region of interest, satisfying the properties of temperature (Tm) and entropy (ΔG), while minimizing the formation of primer-dimers in a pooled experiment, is therefore necessary. Manual design and selection of oligos becomes very challenging, complicated by factors such as length of the target region and number of targeted regions. Here we describe, AnthOligo, a web-based application developed to optimally automate the process of generation of oligo sequences used to target and capture the continuum of large and complex genomic regions. Apart from generating oligos for RSE, this program may have wider applications in the design of customizable internal oligos to be used as baits for gene panel analysis or even probes for large-scale comparative genomic hybridization array processes. AnthOligo was tested by capturing the Major Histocompatibility Complex (MHC) of a random sample.The application provides users with a simple interface to upload an input file in BED format and customize parameters for each task. The task of probe design in AnthOligo commences when a user uploads an input file and concludes with the generation of a result-set containing an optimal set of region-specific oligos. AnthOligo is currently available as a public web application with URL: http://antholigo.chop.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Defining Blood Group Gene Reference Alleles by Long-Read Sequencing: Proof of Concept in the ACKR1 Gene Encoding the Duffy Antigens

Background

In the novel era of blood group genomics, (re-)defining reference gene/allele sequences of blood group genes has become an important goal to achieve, both for diagnostic and research purposes. As novel potent sequencing technologies are available, we thought to investigate the variability encountered in the three most common alleles of ACKR1, the gene encoding the clinically relevant Duffy antigens, at the haplotype level by a long-read sequencing approach.

Materials and Methods

After long-range PCR amplification spanning the whole ACKR1 gene locus (∼2.5 kilobases), amplicons generated from 81 samples with known genotypes were sequenced in a single read by using the Pacific Biosciences (PacBio) single molecule, real-time (SMRT) sequencing technology.

Results

High-quality sequencing reads were obtained for the 162 alleles (accuracy >0.999). Twenty-two nucleotide variations reported in databases were identified, defining 19 haplotypes: four, eight, and seven haplotypes in 46 ACKR1*01, 63 ACKR1*02, and 53 ACKR1*02N.01 alleles, respectively.

Discussion

Overall, we have defined a subset of reference alleles by third-generation (long-read) sequencing. This technology, which provides a "longitudinal" overview of the loci of interest (several thousand base pairs) and is complementary to the second-generation (short-read) next-generation sequencing technology, is of critical interest for resolving novel, rare, and null alleles.

Weighted gene co-expression network analysis unveils gene networks associated with the Fusarium head blight resistance in tetraploid wheat

BACKGROUND

Fusarium head blight (FHB) resistance in the durum wheat breeding gene pool is rarely reported. Triticum turgidum ssp. carthlicum line Blackbird is a tetraploid relative of durum wheat that offers partial FHB resistance. Resistance QTL were identified for the durum wheat cv. Strongfield × Blackbird population on chromosomes 1A, 2A, 2B, 3A, 6A, 6B and 7B in a previous study. The objective of this study was to identify the defense mechanisms underlying the resistance of Blackbird and report candidate regulator defense genes and single nucleotide polymorphism (SNP) markers within these genes for high-resolution mapping of resistance QTL reported for the durum wheat cv. Strongfield/Blackbird population.

RESULTS

Gene network analysis identified five networks significantly (P < 0.05) associated with the resistance to FHB spread (Type II FHB resistance) one of which showed significant correlation with both plant height and relative maturity traits. Two gene networks showed subtle differences between Fusarium graminearum-inoculated and mock-inoculated plants, supporting their involvement in constitutive defense. The candidate regulator genes have been implicated in various layers of plant defense including pathogen recognition (mainly Nucleotide-binding Leucine-rich Repeat proteins), signaling pathways including the abscisic acid and mitogen activated protein (MAP) kinase, and downstream defense genes activation including transcription factors (mostly with dual roles in defense and development), and cell death regulator and cell wall reinforcement genes. The expression of five candidate genes measured by quantitative real-time PCR was correlated with that of RNA-seq, corroborating the technical and analytical accuracy of RNA-sequencing.

CONCLUSIONS

Gene network analysis allowed identification of candidate regulator genes and genes associated with constitutive resistance, those that will not be detected using traditional differential expression analysis. This study also shed light on the association of developmental traits with FHB resistance and partially explained the co-localization of FHB resistance with plant height and maturity QTL reported in several previous studies. It also allowed the identification of candidate hub genes within the interval of three previously reported FHB resistance QTL for the Strongfield/Blackbird population and associated SNPs for future high resolution mapping studies.

In-silico Design of DNA Oligonucleotides: Challenges and Approaches

DNA oligonucleotides are essential components of a high number of technologies in molecular biology. The key event of each oligonucleotide-based assay is the specific binding between oligonucleotides and their target DNA. However, single-stranded DNA molecules also tend to bind to unintended targets or themselves. The probability of such unspecific binding increases with the complexity of an assay. Therefore, accurate data management and design workflows are necessary to optimize the in-silico design of primers and probes. Important considerations concerning computational infrastructure and run time need to be made for both data management and the design process. Data retrieval, data updates, storage, filtering and analysis are the main parts of a sequence data management system. Each part needs to be well-implemented as the resulting sequences form the basis for the oligonucleotide design. Important key features, such as the oligonucleotide length, melting temperature, secondary structures and primer dimer formation, as well as the specificity, should be considered for the in-silico selection of oligonucleotides. The development of an efficient oligonucleotide design workflow demands the right balance between the precision of the applied computer models, the general expenditure of time, and computational workload. This paper gives an overview of important parameters during the design process, starting from the data retrieval, up to the design parameters for optimized oligonucleotide design.

Clustering of circular consensus sequences: accurate error correction and assembly of single molecule real-time reads from multiplexed amplicon libraries

BACKGROUND

Targeted resequencing with high-throughput sequencing (HTS) platforms can be used to efficiently interrogate the genomes of large numbers of individuals. A critical issue for research and applications using HTS data, especially from long-read platforms, is error in base calling arising from technological limits and bioinformatic algorithms. We found that the community standard long amplicon analysis (LAA) module from Pacific Biosciences is prone to substantial bioinformatic errors that raise concerns about findings based on this pipeline, prompting the need for a new method.

RESULTS

A single molecule real-time (SMRT) sequencing-error correction and assembly pipeline, C3S-LAA, was developed for libraries of pooled amplicons. By uniquely leveraging the structure of SMRT sequence data (comprised of multiple low quality subreads from which higher quality circular consensus sequences are formed) to cluster raw reads, C3S-LAA produced accurate consensus sequences and assemblies of overlapping amplicons from single sample and multiplexed libraries. In contrast, despite read depths in excess of 100X per amplicon, the standard long amplicon analysis module from Pacific Biosciences generated unexpected numbers of amplicon sequences with substantial inaccuracies in the consensus sequences. A bootstrap analysis showed that the C3S-LAA pipeline per se was effective at removing bioinformatic sources of error, but in rare cases a read depth of nearly 400X was not sufficient to overcome minor but systematic errors inherent to amplification or sequencing.

CONCLUSIONS

C3S-LAA uses a divide and conquer processing algorithm for SMRT amplicon-sequence data that generates accurate consensus sequences and local sequence assemblies. Solving the confounding bioinformatic source of error in LAA allowed for the identification of limited instances of errors due to DNA amplification or sequencing of homopolymeric nucleotide tracts. For research and development in genomics, C3S-LAA allows meaningful conclusions and biological inferences to be made from accurately polished sequence output.