IDSSR: An Efficient Pipeline for Identifying Polymorphic Microsatellites from a Single Genome Sequence

Comparative analysis of the medicinal plant Polygonatum kingianum (Asparagaceae) with related verticillate leaf types of the Polygonatum species based on chloroplast genomes

Background

Polygonatum kingianum has been widely used as a traditional Chinese medicine as well as a healthy food. Because of its highly variable morphology, this medicinal plant is often difficult to distinguish from other related verticillate leaf types of the Polygonatum species. The contaminants in P. kingianum products not only decrease the products' quality but also threaten consumer safety, seriously inhibiting the industrial application of P. kingianum.

Methods

Nine complete chloroplast (cp) genomes of six verticillate leaf types of the Polygonatum species were de novo assembled and systematically analyzed.

Results

The total lengths of newly sequenced cp genomes ranged from 155,437 to 155,977 bp, including 86/87 protein-coding, 38 tRNA, and 8 rRNA genes, which all exhibited well-conserved genomic structures and gene orders. The differences in the IR/SC (inverted repeats/single-copy) boundary regions and simple sequence repeats were detected among the verticillate leaf types of the Polygonatum cp genomes. Comparative cp genomes analyses revealed that a higher similarity was conserved in the IR regions than in the SC regions. In addition, 11 divergent hotspot regions were selected, providing potential molecular markers for the identification of the Polygonatum species with verticillate leaf types. Phylogenetic analysis indicated that, as a super barcode, plastids realized a fast and efficient identification that clearly characterized the relationships within the verticillate leaf types of the Polygonatum species. In brief, our results not only enrich the data on the cp genomes of the genus Polygonatum but also provide references for the P. kingianum germplasm resource protection, herbal cultivation, and drug production.

Conclusion

This study not only accurately identifies P. kingianum species, but also provides valuable information for the development of molecular markers and phylogenetic analyses of the Polygonatum species with verticillate leaf types.

Pipeline for developing polymorphic microsatellites in species without reference genomes

Microsatellites, also known as simple sequence repeats (SSRs), are the preferred type of marker for many genetic applications. In conjunction with the ongoing development of next-generation sequencing, several bioinformatic tools have been developed for identifying SSRs from genomic or transcriptomic sequences. Although these tools are handy for generating polymorphic SSRs, their application almost always depends on an existing reference genome or self-assembly of the reference genome. With this in mind, we propose a pipeline for developing polymorphic SSRs that may be applied to species without reference genomes. Using a species without a reference genome (black Amur bream; Megalobrama terminalis Richardson, 1846) as a model, our pipeline was able to effectively discover polymorphic SSRs. Under different R parameters of a reference-free single nucleotide polymorphisms (SNPs) caller (ebwt2InDel), a total of 258, 208, 102, and 11 polymorphic SSRs were mined. To quantify the accuracy of the polymorphic SSRs detected using our pipeline, we analyzed 25 SSRs with PCR experiments. All primers were successfully amplified, and most SSRs (23 SSRs, 92%) were polymorphic. From the 36 individual black Amur bream, we acquired an average of 3.36 alleles per locus, ranging from one to 11. This demonstrates the effectiveness of our pipeline in identifying polymorphic SSRs and designing primers for SSR genotyping. Ultimately, our pipeline can effectively mine polymorphic SSRs for species without reference genomes, complementing SSR mining approaches based on reference genomes and helping to resolve biological issues that accompany these methods.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03313-0.

Comparative analysis of the plastid and mitochondrial genomes of Artemisia giraldii Pamp

Artemisia giraldii Pamp. is an herbaceous plant distributed only in some areas in China. To understand the evolutionary relationship between plastid and mitochondria in A. giraldii, we sequenced and analysed the plastome and mitogenome of A. giraldii on the basis of Illumina and Nanopore DNA sequencing data. The mitogenome was 194,298 bp long, and the plastome was 151,072 bp long. The mitogenome encoded 56 genes, and the overall GC content was 45.66%. Phylogenetic analysis of the two organelle genomes revealed that A. giraldii is located in the same branching position. We found 13 pairs of homologous sequences between the plastome and mitogenome, and only one of them might have transferred from the plastid to the mitochondria. Gene selection pressure analysis in the mitogenome showed that ccmFc, nad1, nad6, atp9, atp1 and rps12 may undergo positive selection. According to the 18 available plastome sequences, we found 17 variant sites in two hypervariable regions that can be used in completely distinguishing 18 Artemisia species. The most interesting discovery was that the mitogenome of A. giraldii was only 43,226 bp larger than the plastome. To the best of our knowledge, this study represented one of the smallest differences between all sequenced mitogenomes and plastomes from vascular plants. The above results can provide a reference for future taxonomic and molecular evolution studies of Asteraceae species.

AutomAted RepeaT Identifier (AARTI): A tool to identify common, polymorphic, and unique microsatellites

Here we are presenting an automated computational pipeline used to mine 5976 mitochondrial genomes to identify common, polymorphic, and unique microsatellites also known as simple sequence repeats (SSRs). Microsatellites are repetitive motifs of 1-6 bases in a DNA sequence. Due to their abundance and highly polymorphic nature, microsatellites have become one of the widely used molecular/genetic markers valuable for many studies including gene tagging, genetic diversity, and species identification. Several computational tools dedicated to mine and categorize microsatellites in nucleotide sequences were developed; however, there is no tool which can identify unique, common and polymorphic microsatellites between each pair of nucleotide sequences. To explore such microsatellites, we have developed a fully automated computational pipeline named AutomAted RepeaT Identifier (AARTI). The AARTI is the only tool till date, which identifies common, polymorphic, and unique microsatellites between each pair of nucleotide sequences. The computational pipeline was constructed using the PERL programming language and the web server for the pipeline was developed with the help of PHP, HTML, CSS, and JavaScript. It was successfully tested to reproduce the results of previous study on 7 mitochondrial genomes of genus Orthotrichum. Moreover, the pipeline was also applied on 5846 (Metazoa) and 130 (Viridiplantae) mitochondrial genomes. The AARTI is freely available at https://lms.snu.edu.in/aarti/ and will certainly accelerate the studies of length variation in microsatellites between species. Additionally, it will be useful in comparative genomic studies, for the computational characterization of microsatellites, and has the potential to be a routine genome analysis pipeline for mitochondrial genomes.

BigFiRSt: A Software Program Using Big Data Technique for Mining Simple Sequence Repeats From Large-Scale Sequencing Data

Background

Simple Sequence Repeats (SSRs) are short tandem repeats of nucleotide sequences. It has been shown that SSRs are associated with human diseases and are of medical relevance. Accordingly, a variety of computational methods have been proposed to mine SSRs from genomes. Conventional methods rely on a high-quality complete genome to identify SSRs. However, the sequenced genome often misses several highly repetitive regions. Moreover, many non-model species have no entire genomes. With the recent advances of next-generation sequencing (NGS) techniques, large-scale sequence reads for any species can be rapidly generated using NGS. In this context, a number of methods have been proposed to identify thousands of SSR loci within large amounts of reads for non-model species. While the most commonly used NGS platforms (e.g., Illumina platform) on the market generally provide short paired-end reads, merging overlapping paired-end reads has become a common way prior to the identification of SSR loci. This has posed a big data analysis challenge for traditional stand-alone tools to merge short read pairs and identify SSRs from large-scale data.

Results

In this study, we present a new Hadoop-based software program, termed BigFiRSt, to address this problem using cutting-edge big data technology. BigFiRSt consists of two major modules, BigFLASH and BigPERF, implemented based on two state-of-the-art stand-alone tools, FLASH and PERF, respectively. BigFLASH and BigPERF address the problem of merging short read pairs and mining SSRs in the big data manner, respectively. Comprehensive benchmarking experiments show that BigFiRSt can dramatically reduce the execution times of fast read pairs merging and SSRs mining from very large-scale DNA sequence data.

Conclusions

The excellent performance of BigFiRSt mainly resorts to the Big Data Hadoop technology to merge read pairs and mine SSRs in parallel and distributed computing on clusters. We anticipate BigFiRSt will be a valuable tool in the coming biological Big Data era.

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.

The complete plastomes of red fleshed pitaya (Selenicereus monacanthus) and three related Selenicereus species: insights into gene losses, inverted repeat expansions and phylogenomic implications

Selenicereus is a genus of perennial shrub from the family Cactaceae, and some of them play an important role in the food industry, pharmaceuticals, cosmetics and medicine. To date, there are few reports on Selenicereus plastomes, which limits our understanding of this genus. Here, we have reported the complete plastomes of four Selenicereus species (S. monacanthus, S. annthonyanus, S. grandifloras, and S. validus) and carried out a comprehensive comparative analysis. All four Selenicereus plastomes have a typical quartile structure. The plastome size ranged from 133,146 to 134,450 bp, and contained 104 unique genes, including 30 tRNA genes, 4 rRNA genes and 70 protein-coding genes. Comparative analysis showed that there were massive losses of ndh genes in Selenicereus. Besides, we observed the inverted repeat regions had undergone a dramatic expansion and formed a previously unreported small single copy/inverted repeat border in the intron region of the atpF gene. Furthermore, we identified 6 hypervariable regions (trnF-GAA-rbcL, ycf1, accD, clpP-trnS-GCU, clpP-trnT-CGU and rpl22-rps19) that could be used as potential DNA barcodes for the identification of Selenicereus species. Our study enriches the plastome in the family Cactaceae, and provides the basis for the reconstruction of phylogenetic relationships.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01121-z.

Assembly of the complete mitochondrial genome of an endemic plant, Scutellaria tsinyunensis, revealed the existence of two conformations generated by a repeat-mediated recombination

We assembled the complete mitochondrial genome of Scutellaria tsinyunensis in this study. Repeat-mediated recombination resulted in the formation of two conformations of the mitochondrial genome in S. tsinyunensis. Scutellaria tsinyunensis belongs to the family Lamiaceae, distributed only in the Jinyun Mountain, Chongqing, China. As a valuable endemic and small population species, it is regarded as a natural resource potentially with significant economic and ecological importance. In this study, we assembled a complete and gap-free mitochondrial genome of S. tsinyunensis. This genome had a length of 354,073 bp and the base composition of the genome was A (27.44%), T (27.30%), C (22.58%), and G (22.68%). This genome encodes 59 genes, including 32 protein-coding genes, 24 tRNA genes, and 3 rRNA genes. The Sanger sequencing and Oxford Nanopore sequencing confirmed a pair of direct repeats had mediated genome recombination, resulting in the formation of two conformations. The gene conversation between plastome and mitochondrial genome was also observed in S. tsinyunensis by detecting gene migration, including six tRNA genes (namely, trnW-CCA, trnI-CAU, trnH-UUU, trnD-GUC, trnN-GUU, and trnM-CAU), five protein-coding gene fragments, and the fragments from 2 rRNA genes. Moreover, the dN/dS analysis revealed the atp9 gene had undergone strong negative selection, and four genes (atp4, mttB, ccmFc, and ccmB) probably had undergone positive selection during evolution in Lamiales. This work reported the first mitochondrial genome of S. tsinyunensis, which could be used as a reference genome for the important medicinal plants of the genus Scutellaria, and also provide much-desired information for molecular breeding.

Comparative plastid genomics of four Pilea (Urticaceae) species: insight into interspecific plastid genome diversity in Pilea

BACKGROUND

Pilea is a genus of perennial herbs from the family Urticaceae, and some species are used as courtyard ornamentals or for medicinal purposes. At present, there is no information about the plastid genome of Pilea, which limits our understanding of this genus. Here, we report 4 plastid genomes of Pilea taxa (Pilea mollis, Pilea glauca 'Greizy', Pilea peperomioides and Pilea serpyllacea 'Globosa') and performed comprehensive comparative analysis.

RESULTS

The four plastid genomes all have a typical quartile structure. The lengths of the plastid genomes ranged from 150,398 bp to 152,327 bp, and each genome contained 113 unique genes, including 79 protein-coding genes, 4 rRNA genes, and 30 tRNA genes. Comparative analysis showed a rather high level of sequence divergence in the four genomes. Moreover, eight hypervariable regions were identified (petN-psbM, psbZ-trnG-GCC, trnT-UGU-trnL-UAA, accD-psbI, ndhF-rpl32, rpl32-trnL-UAG, ndhA-intron and ycf1), which are proposed for use as DNA barcode regions. Phylogenetic relationships based on the plastid genomes of 23 species of 14 genera of Urticaceae resulted in the placement of Pilea in the middle and lower part of the phylogenetic tree, with 100% bootstrap support within Urticaceae.

CONCLUSION

Our results enrich the resources concerning plastid genomes. Comparative plastome analysis provides insight into the interspecific diversity of the plastid genome of Pilea. The identified hypervariable regions could be used for developing molecular markers applicable in various research areas.

SSRMMD: A Rapid and Accurate Algorithm for Mining SSR Feature Loci and Candidate Polymorphic SSRs Based on Assembled Sequences

Microsatellites or simple sequence repeats (SSRs) are short tandem repeats of DNA widespread in genomes and transcriptomes of diverse organisms and are used in various genetic studies. Few software programs that mine SSRs can be further used to mine polymorphic SSRs, and these programs have poor portability, have slow computational speed, are highly dependent on other programs, and have low marker development rates. In this study, we develop an algorithm named Simple Sequence Repeat Molecular Marker Developer (SSRMMD), which uses improved regular expressions to rapidly and exhaustively mine perfect SSR loci from any size of assembled sequence. To mine polymorphic SSRs, SSRMMD uses a novel three-stage method to assess the conservativeness of SSR flanking sequences and then uses the sliding window method to fragment each assembled sequence to assess its uniqueness. Furthermore, molecular biology assays support the polymorphic SSRs identified by SSRMMD. SSRMMD is implemented using the Perl programming language and can be downloaded from https://github.com/GouXiangJian/SSRMMD.