A new resource for the development of SSR markers: Millions of loci from a thousand plant transcriptomes

Genotyping-by-sequencing and weighted gene co-expression network analysis of genes responsive against Potato virus Y in commercial potato cultivars

Potato is considered a key component of the global food system and plays a vital role in strengthening world food security. A major constraint to potato production worldwide is the Potato Virus Y (PVY), belonging to the genus Potyvirus in the family of Potyviridae. Selective breeding of potato with resistance to PVY pathogens remains the best method to limit the impact of viral infections. Understanding the genetic diversity and population structure of potato germplasm is important for breeders to improve new cultivars for the sustainable use of genetic materials in potato breeding to PVY pathogens. While, genetic diversity improvement in modern potato breeding is facing increasingly narrow genetic basis and the decline of the genetic diversity. In this research, we performed genotyping-by-sequencing (GBS)-based diversity analysis on 10 commercial potato cultivars and weighted gene co-expression network analysis (WGCNA) to identify candidate genes related to PVY-resistance. WGCNA is a system biology technique that uses the WGCNA R software package to describe the correlation patterns between genes in multiple samples. In terms of consumption, these cultivars are a high rate among Iranian people. Using population structure analysis, the 10 cultivars were clustered into three groups based on the 118343 single nucleotide polymorphisms (SNPs) generated by GBS. Read depth ranged between 5 and 18. The average data size and Q30 of the reads were 145.98 Mb and 93.63%, respectively. Based on the WGCNA and gene expression analysis, the StDUF538, StGTF3C5, and StTMEM161A genes were associated with PVY resistance in the potato genome. Further, these three hub genes were significantly involved in defense mechanism where the StTMEM161A was involved in the regulation of alkalization apoplast, the StDUF538 was activated in the chloroplast degradation program, and the StGTF3C5 regulated the proteins increase related to defense in the PVY infected cells. In addition, in the genetic improvement programs, these hub genes can be used as genetic markers for screening commercial cultivars for PVY resistance. Our survey demonstrated that the combination of GBS-based genetic diversity germplasm analysis and WGCNA can assist breeders to select cultivars resistant to PVY as well as help design proper crossing schemes in potato breeding.

Efficient PCR Amplification Protocol of Nuclear Microsatellites for Exuviae-Derived DNA of Cicada, Yezoterpnosia nigricosta

Although insect exuviae-based genetics is challenging, it can be a valuable method for obtaining reliable DNA resources by non-invasive sampling. This approach is especially effective when the target species is endangered/endemic or when sampling the adult is difficult. One example is cicadas, which during molt leave their exoskeletons on tree trunks, making them easily collectable. While cicada exuviae-derived DNA has previously been employed for mitochondrial DNA sequencing, this study aimed to develop a reliable method for the PCR amplification of nuclear microsatellite loci from cicada exuviae derived DNA for application in molecular ecology, conservation and population genetics. Five different PCR amplification protocols were performed, and the fragment patterns compared with those obtained using DNA extracted from adult individuals. Moreover, the relationship between the freshness of the exuviae and genotyping success was evaluated. TaKaRa LA Taq provided the best performance in the PCR amplification of DNA isolated from cicada exuviae and the electropherogram showed a clear fragment pattern that was equivalent to that obtained from the DNA extracted from the adult individual. This result suggests that cicada exuviae-derived DNA can be amplified by PCR and that multiple independent loci of nuclear DNA microsatellite markers can be easily genotyped. This study demonstrates that fresh cicada exuviae provide high quality DNA, which can be used for microsatellite genotyping. The methods developed in this study are applicable not only for cicada but other insect species for which exuviae are available. Thus, this study can make a significant contribution to insect sciences.

Comparative analysis of genomic- and EST-SSRs in European plum (Prunus domestica L.): implications for the diversity analysis of polyploids

Simple sequence repeats (SSRs) are among the most useful DNA markers in plant science. The aim of this study was to compare the features and usefulness of genomic SSR (gSSR) and EST-SSRs in European plum (Prunus domestica L.), an economically important, hexaploid stone fruit crop globally cultivated to produce fleshy fruits and derived foodstuff. The analysis of an ample set of morphologically diverse varieties indicated that gSSRs and EST-SSRs provide different estimates of some of the locus-based indicators of diversity. Moreover, the two classes of SSRs gave different, weakly correlated, estimations of distance-based parameters with gSSRs being more powerful for discriminating purposes. The two SSR classes provide complementary information in European plum, making the contribution of EST-SSRs useful not only as non-neutral markers. The differences between SSR classes are discussed considering the neutral and non-neutral evolution, and the polyploidy and asexual propagation of the cultivated tree varieties.

Sequencing and Analyzing the Transcriptomes of a Thousand Species Across the Tree of Life for Green Plants

The 1,000 Plants (1KP) initiative was the first large-scale effort to collect next-generation sequencing (NGS) data across a phylogenetically representative sampling of species for a major clade of life, in this case theViridiplantae, or green plants. As an international multidisciplinary consortium, we focused on plant evolution and its practical implications. Among the major outcomes were the inference of a reference species tree for green plants by phylotranscriptomic analysis of low-copy genes, a survey of paleopolyploidy (whole-genome duplications) across the Viridiplantae, the inferred evolutionary histories for many gene families and biological processes, the discovery of novel light-sensitive proteins for optogenetic studies in mammalian neuroscience, and elucidation of the genetic network for a complex trait (C₄ photosynthesis). Altogether, 1KP demonstrated how value can be extracted from a phylodiverse sequencing data set, providing a template for future projects that aim to generate even more data, including complete de novo genomes, across the tree of life.

Plant virus interaction mechanism and associated pathways in mosaic disease of small cardamom (Elettaria cardamomum Maton) by RNA-Seq approach

Small cardamom (Elettaria cardamomum), grown in limited coastal tropical countries is one of the costliest and widely exported agri-produce having global turnover of >10 billion USD. Mosaic/marble disease is one of the major impediments that requires understanding of disease at molecular level. Neither whole genome sequence nor any genomic resources are available, thus RNA seq approach can be a rapid and economical alternative. De novo transcriptome assembly was done with Illumina Hiseq data. A total of 5317 DEGs, 2267 TFs, 114 pathways and 175,952 genic region putative markers were obtained. Gene regulatory network analysis deciphered molecular events involved in marble disease. This is the first transcriptomic report revealing disease mechanism mediated by perturbation in auxin homeostasis and ethylene signalling leading to senescence. The web-genomic resource (SCMVTDb) catalogues putative molecular markers, candidate genes and transcript information. SCMVTDb can be used in germplasm improvement against mosaic disease in endeavour of small cardamom productivity. Availability of genomic resource, SCMVTDb: http://webtom.cabgrid.res.in/scmvtdb/.

Optimizing depth and type of high-throughput sequencing data for microsatellite discovery

PREMISE

Simple sequence repeat (SSR) markers (microsatellites) are a mainstay of many labs, especially when working on a limited budget, carrying out preliminary analyses, and in teaching. Whether SSRs mined from plant genomes or transcriptomes are preferred for certain applications, and the depth of sequencing needed to allow efficient SSR discovery, has not been tested.

METHODS

I used genome and transcriptome high-throughput sequencing data at a range of sequencing depths to compare efficacy of SSR identification. I then tested primers from tomato for amplification, polymorphism, and transferability to related species.

RESULTS

Small assemblies (two million read pairs) identified ca. 200-2000 potential markers from the genome assemblies and ca. 600-3650 from the transcriptome assemblies. Genome-derived contigs were often short, potentially precluding primer design. Genomic SSR primers were less transferable across species but exhibited greater variation (partially explained by being composed of more repeat units) than transcriptome-derived primers.

DISCUSSION

Small high-throughput sequencing resources may be sufficient for identification of hundreds of SSRs. Genomic data may be preferable in species with low polymorphism, but transcriptome data may result in longer loci (more amenable to primer design) and primers may be more transferable to related species.

An RNA Sequencing Transcriptome Analysis and Development of EST-SSR Markers in Chinese Hawthorn through Illumina Sequencing

Chinese hawthorn (Crataegus pinnatifida) is an important ornamental and economic horticultural plant. However, the lack of molecular markers has limited the development and utilization of hawthorn germplasm resources. Simple sequence repeats (SSRs) derived from expressed sequence tags (ESTs) allow precise and effective cultivar characterization and are routinely used for genetic diversity analysis. Thus, we first reported the development of polymorphic EST-SSR markers in C. pinnatifida with perfect repeats using Illumina RNA-Seq technique. In total, we investigated 14,364 unigenes, from which 5091 EST-SSR loci were mined. Di-nucleotides (2012, 39.52%) were the most abundant SSRs, followed by mono- (1989, 39.07%), and tri-nucleotides (1024, 20.11%). On the basis of these EST-SSRs, a total of 300 primer pairs were designed and used for polymorphism analysis in 70 accessions collected from different geographical regions of China. Of 239 (79.67%) pairs of primer-generated amplification products, 163 (54.33%) pairs of primers showed polymorphism. Finally, 33 primers with high polymorphism were selected for genetic diversity analysis and tested on 70 individuals with low-cost fluorescence-labeled M13 primers using capillary electrophoresis genotyping platform. A total of 108 alleles were amplified by 33 SSR markers, with the number of alleles (Na) ranging from 2 to 14 per locus (mean: 4.939), and the effective number of alleles (Ne) ranging from 1.258 to 3.214 (mean: 2.221). The mean values of gene diversity (He), observed heterozygosity (Ho), and polymorphism information content (PIC) were 0.524 (range 0.205–0.689), 0.709 (range 0.132–1.000), and 0.450 (range 0.184–0.642), respectively. Furthermore, the dendrogram constructed based on the EST-SSR separated the cultivars into two main clusters. In sum, our study was the first comprehensive study on the development and analysis of a large set of SSR markers in hawthorn. The results suggested that the use of NGS techniques for SSR development represented a powerful tool for genetic studies. Additionally, fluorescence-labeled M13 markers proved to be a valuable method for genotyping. All of these EST-SSR markers have agronomic potential and constitute a scientific basis for future studies on the identification, classification, and innovation of hawthorn germplasms.

Transcriptome-wide mining, characterization, and development of microsatellite markers in Lychnis kiusiana (Caryophyllaceae)

BACKGROUND

Lychnis kiusiana Makino is an endangered perennial herb native to wetland areas in Korea and Japan. Despite its conservational and evolutionary significance, population genetic resources are lacking for this species. Next-generation sequencing has been accepted as a rapid and cost-effective solution for the identification of microsatellite markers in nonmodel plants.

RESULTS

Using Illumina HiSeq 2000 sequencing technology, we assembled 67,498,600 reads into 91,900 contigs and identified 11,403 microsatellite repeat motifs in 9563 contigs. A total of 4510 microsatellite-containing transcripts had Gene Ontology (GO) annotations, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified 124 pathways with significant scores. Many microsatellites in the L. kiusiana leaf transcriptome were linked to genes involved in the plant response to light intensity, salt stress, temperature stimulus, and nutrient and water deprivation. A total of 12,486 single-nucleotide polymorphisms (SNPs) were identified on transcripts harboring microsatellites. The analysis of nucleotide substitution rates for 2389 unigenes indicated that 39 genes were under strong positive selection. The primers of 6911 microsatellites were designed, and 40 of 50 selected primer pairs were consistently and successfully amplified from 51 individuals. Twenty-five of these were polymorphic, and the average number of alleles per SSR locus was 6.96, with a range from 2 to 15. The observed and expected heterozygosities ranged from 0.137 to 0.902 and 0.131 to 0.827, respectively, and locus-specific F_IS estimates ranged from - 0.116 to 0.290. Eleven of the 25 primer pairs were successfully amplified in three additional species of Lychnis: 56% in L. wilfordii, 64% in L. cognata and 80% in L. fulgens.

CONCLUSIONS

The transcriptomic SSR markers of Lychnis kiusiana provide a valuable resource for understanding the population genetics, evolutionary history, and effective conservation management of this species. Furthermore, the identified microsatellite loci linked to the annotated genes should be useful for developing functional markers of L. kiusiana. The developed markers represent a potentially valuable source of transcriptomic SSR markers for population genetic analyses with moderate levels of cross-taxon portability.

High-throughput SSR marker development and its application in a centipedegrass (Eremochloa ophiuroides (Munro) Hack.) genetic diversity analysis

Centipedegrass (Eremochloa ophiuroides (Munro) Hack.) is a perennial, warm-season C4 grass species that shows great potential for use as a low-maintenance turfgrass species in tropical and subtropical regions. However, limited genetic and genomic information is available for this species, which has impeded systematic studies on the enhancement of its turf quality and resistance against biotic and abiotic stress. In this study, Illumina HiSeq high-throughput sequencing technology was performed to generate centipedegrass transcriptome sequences. A total of 352,513 assembled sequences were used to search for simple sequence repeat (SSR) loci, and 64,470 SSR loci were detected in 47,638 SSR containing sequences. The tri-nucleotides were the most frequent repeat motif, followed by di-nucleotides, tetra-nucleotides hexnucleotides, and pentanucleotides. A total of 48,061 primer pairs were successfully designed in the flanking sequences of the SSRs, and 100 sets of primers were randomly selected for the initial validation in four centipedegrass accessions. In total, 56 (56.0%) of the 100 primer pairs tested successfully amplified alleles from all four centipedegrass accessions, while 50 were identified as polymorphic markers and were then used to assess the level of genetic diversity among 43 centipedegrass core collections. The genetic diversity analysis exhibited that the number of alleles (Na) per locus ranged from 3 to 13, and the observed heterozygosity (Ho) ranged from 0.17 to 0.83. The polymorphism information content (PIC) value of the markers ranged from 0.15 to 0.78, and the genetic distances (coefficient Nei72) between the accessions varied from 0.07 to 0.48. The UPGMA-based dendrogram clustered all 43 core collections into two main groups and six subgroups, which further validated the effectiveness of these newly developed SSR markers. Hence, these newly developed SSR markers will be valuable and potentially useful for future genetic and genomic studies of E. ophiuroides.

Chloroplast genome analyses and genomic resource development for epilithic sister genera Oresitrophe and Mukdenia (Saxifragaceae), using genome skimming data

Background

Epilithic sister genera Oresitrophe and Mukdenia (Saxifragaceae) have an epilithic habitat (rocky slopes) and a parapatric distribution in East Asia, which makes them an ideal model for a more comprehensive understanding of the demographic and divergence history and the influence of climate changes in East Asia. However, the genetic background and resources for these two genera are scarce.

Results

The complete chloroplast (cp) genomes of two Oresitrophe rupifraga and one Mukdenia rossii individuals were reconstructed and comparative analyses were conducted to examine the evolutionary pattern of chloroplast genomes in Saxifragaceae. The cp genomes ranged from 156,738 bp to 156,960 bp in length and had a typical quadripartite structure with a conserved genome arrangement. Comparative analysis revealed the intron of rpl2 has been lost in Heuchera parviflora, Tiarella polyphylla, M. rossii and O. rupifraga but presents in the reference genome of Penthorum chinense. Seven cp hotspot regions (trnH-psbA, trnR-atpA, atpI-rps2, rps2-rpoC2, petN-psbM, rps4-trnT and rpl33-rps18) were identified between Oresitrophe and Mukdenia, while four hotspots (trnQ-psbK, trnR-atpA, trnS-psbZ and rpl33-rps18) were identified within Oresitrophe. In addition, 24 polymorphic cpSSR loci were found between Oresitrophe and Mukdenia. Most importantly, we successfully developed 126 intergeneric polymorphic gSSR markers between Oresitrophe and Mukdenia, as well as 452 intrageneric ones within Oresitrophe. Twelve randomly selected intergeneric gSSRs have shown that these two genera exhibit a significant genetic structure.

Conclusions

In this study, we conducted genome skimming for Oresitrophe rupifraga and Mukdenia rossii. Using these data, we were able to not only assemble their complete chloroplast genomes, but also develop abundant genetic resources (cp hotspots, cpSSRs, polymorphic gSSRs). The genomic patterns and genetic resources presented here will contribute to further studies on population genetics, phylogeny and conservation biology in Saxifragaceae.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4633-x) contains supplementary material, which is available to authorized users.

Adding loci improves phylogeographic resolution in red mangroves despite increased missing data: comparing microsatellites and RAD-Seq and investigating loci filtering

The widespread adoption of RAD-Seq data in phylogeography means genealogical relationships previously evaluated using relatively few genetic markers can now be addressed with thousands of loci. One challenge, however, is that RAD-Seq generates complete genotypes for only a small subset of loci or individuals. Simulations indicate that loci with missing data can produce biased estimates of key population genetic parameters, although the influence of such biases in empirical studies is not well understood. Here we compare microsatellite data (8 loci) and RAD-Seq data (six datasets ranging from 239 to 25,198 loci) from red mangroves (Rhizophora mangle) in Florida to evaluate how different levels of data filtering influence phylogeographic inferences. For all datasets, we calculated population genetic statistics and evaluated population structure, and for RAD-Seq datasets, we additionally examined population structure using coalescence. We found higher F_ST using microsatellites, but that RAD-Seq-based estimates approached those based on microsatellites as more loci with more missing data were included. Analyses of RAD-Seq datasets resolved the classic Gulf-Atlantic coastal phylogeographic break, which was not significant in the microsatellite analyses. Applying multiple levels of filtering to RAD-Seq datasets can provide a more complete picture of potential biases in the data and elucidate subtle phylogeographic patterns.

Transcriptomic analysis, genic SSR development, and genetic diversity of proso millet (Panicum miliaceum; Poaceae)

PREMISE OF THE STUDY

Proso millet (Panicum miliaceum; Poaceae) is a minor crop with good nutritional qualities and strong tolerance to drought stress and soil infertility. However, studies on genetic diversity have been limited due to a lack of efficient genetic markers.

METHODS

Illumina sequencing technology was used to generate short read sequences of proso millet, and de novo transcriptome assemblies were used to develop a de novo assembly of proso millet. Genic simple sequence repeat (SSR) markers were identified and used to detect polymorphism among 56 accessions. Population structure and genetic similarity coefficient were estimated.

RESULTS

In total, 25,341 unique gene sequences and 4724 SSR loci were obtained from the transcriptome, of which 229 pairs of SSR primers were validated, which resulted in 14 polymorphic genic SSR primers exhibiting 43 total alleles. According to the ratio of polymorphic markers (6.1%, 14/229), there are potentially 288 polymorphic genic SSR markers available for genetic assay development in the future. Bayesian population analyses showed that the 56 accessions comprised two distinct groups.

DISCUSSION

A genetic structure and cluster assay indicated that the accessions from the Loess Plateau of China shared a high genetic similarity coefficient with those from other regions and that there was no correlation between genetic diversity and geographic origin. The transcriptome sequencing data and millet-specific SSR markers developed in this study establish an excellent resource for gene discovery and may improve the development of breeding programs in proso millet in the future.

Development and characterization of EST-SSR markers in Stipa breviflora (Poaceae)

PREMISE OF THE STUDY

Stipa breviflora (Poaceae) is one of the dominant species of the desert steppe in the eastern Eurasian grasslands. Simple sequence repeat (SSR) markers were developed for use in genetic diversity studies of this species.

METHODS AND RESULTS

A total of 1954 potentially polymorphic loci were obtained by comparing transcriptome data of eight different S. breviflora individuals. We selected 81 loci to verify polymorphism and 63 loci amplified, of which 21 loci exhibited polymorphism. The number of alleles per locus varied from two to 24, the observed heterozygosity ranged from 0.083 to 0.958, and the expected heterozygosity ranged from 0.396 to 0.738.

CONCLUSIONS

These newly identified SSR loci can be used for population genetic and landscape genetic studies of S. breviflora. In addition, 14 loci also amplified in six related Stipa species (S. grandis, S. krylovii, S. bungeana, S. aliena, S. gobica, and S. purpurea).

Transcriptomic resources and marker validation for diploid and polyploid Veronica (Plantaginaceae) from New Zealand and Europe

PREMISE OF THE STUDY

Polyploidy may generate novel variation, leading to adaptation and species diversification. An excellent natural system to study polyploid evolution in a comparative framework is Veronica (Plantaginaceae), which comprises several parallel, recently evolved polyploid series.

METHODS

Over 105 million Illumina paired-end sequence reads were generated from cDNA libraries of leaf tissue from eight individuals representing three European and four New Zealand species. Forty-eight simple sequence repeat (SSR) and 48 low-copy nuclear (LCN) markers were developed and validated with Fluidigm microfluidic PCR and Illumina MiSeq amplicon sequencing on 48 different individuals each.

RESULTS

Individual Trinity assemblies were similar regarding annotated transcripts (13,009-14,271), mean contig length (635-742 bp), N50 value (916-1133 bp), E90N50 value (1099-1308 bp), contigs with positive BLAST hits (42-63%), and gene ontology terms. Analyses of 29,738 single-nucleotide polymorphisms (8746 phylogenetically informative) mined from these transcriptomes plus two outgroups (Picrorhiza kurrooa and Plantago ovata) showed moderate to high bootstrap support for all branches and reticulation among sampled European Veronica.

DISCUSSION

The transcriptome sequences themselves, as well as the validated SSR (40/48) and LCN (11/48) markers derived from them, show inter- and intraspecific genetic variation. These resources will be invaluable for future population genetic, phylogenetic, and functional genetic investigations in polyploid Veronica.

The report of my death was an exaggeration: A review for researchers using microsatellites in the 21st century

Microsatellites, or simple sequence repeats (SSRs), have long played a major role in genetic studies due to their typically high polymorphism. They have diverse applications, including genome mapping, forensics, ascertaining parentage, population and conservation genetics, identification of the parentage of polyploids, and phylogeography. We compare SSRs and newer methods, such as genotyping by sequencing (GBS) and restriction site associated DNA sequencing (RAD-Seq), and offer recommendations for researchers considering which genetic markers to use. We also review the variety of techniques currently used for identifying microsatellite loci and developing primers, with a particular focus on those that make use of next-generation sequencing (NGS). Additionally, we review software for microsatellite development and report on an experiment to assess the utility of currently available software for SSR development. Finally, we discuss the future of microsatellites and make recommendations for researchers preparing to use microsatellites. We argue that microsatellites still have an important place in the genomic age as they remain effective and cost-efficient markers.