EST-SSRs as a resource for population genetic analyses

Transcriptome sequencing of the endangered land snail Karaftohelix adamsi from the Island Ulleung: De novo assembly, annotation, valuation of fitness genes and SSR markers

BACKGROUND

The Bradybaenidae snail Karaftohelix adamsi is endemic to Korea, with the species tracked from Island Ulleung in North Gyeongsang Province of South Korea. K. adamsi has been classified under the Endangered Wildlife Class II species of Korea and poses a severe risk of extinction following habitat disturbances. With no available information at the DNA (genome) or mRNA (transcriptome) level for the species, conservation by utilizing informed molecular resources seems difficult.

OBJECTIVE

In this study, we used the Illumina short-read sequencing and Trinity de novo assembly to draft the reference transcriptome of K. adamsi.

RESULTS

After assembly, 13,753 unigenes were obtained of which 10,511 were annotated to public databases (a maximum of 10,165 unigenes found homologs in PANM DB). A total of 6,351, 3,535, 358, and 3,407 unigenes were ascribed to the functional categories under KOG, GO, KEGG, and IPS, respectively. The transcripts such as the HSP 70, aquaporin, TLR, and MAPK, among others, were screened as putative functional resources for adaptation. DNA transposons were found to be thickly populated in comparison to retrotransposons in the assembled unigenes. Further, 2,164 SSRs were screened with the promiscuous presence of dinucleotide repeats such as AC/GT and AG/CT.

CONCLUSION

The transcriptome-guided discovery of molecular resources in K. adamsi will not only serve as a basis for functional genomics studies but also provide sustainable tools to be utilized for the protection of the species in the wild. Moreover, the development of polymorphic SSRs is valuable for the identification of species from newer habitats and cross-species genotyping.

Genetic Diversity of Rhipicephalus (Boophilus) microplus for a Global Scenario: A Comprehensive Review

Rhipicephalus microplus poses a substantial threat to livestock health and agricultural economies worldwide. Its remarkable adaptability to diverse environments and hosts is a testament to its extensive genetic diversity. This review delves into the genetic diversity of R. microplus, employing three pivotal genetic markers: the cytochrome c oxidase I (COX1) gene, ribosomal genes, and microsatellites. The COX1 gene, a crucial tool for genetic characterization and phylogenetic clustering, provides insights into the adaptability of ticks. Ribosomal genes, such as internal transcribed spacer regions (ITS-1 and2) as well as 18S and 28S, are routinely utilized for species differentiation. However, their use is limited due to indels (insertions and deletions). Microsatellites and minisatellites, known for their high polymorphism, have been successfully employed to study populations and genetic diversity across various tick species. Despite their effectiveness, challenges such as null alleles and marker variations warrant careful consideration. Bm86, a well-studied vaccine candidate, exhibits substantial genetic diversity. This diversity directly influences vaccine efficacy, posing challenges for developing a universally effective Bm86-based vaccine. Moreover, the review emphasizes the prevalence of genes associated with synthetic pyrethroid resistance. Identifying single nucleotide polymorphisms in the acaricide-resistant genes of R. microplus has facilitated the development of molecular markers for detecting and monitoring resistance against synthetic pyrethroids. However, mutations in sodium channels, the target site for synthetic pyrethroid, correlate well with the resistance status of R. microplus, which is not the case with other acaricide target genes. This study underscores the importance of understanding genetic diversity in developing effective tick management strategies. The choice of genetic marker should be tailored based on the level of taxonomic resolution and the group of ticks under investigation. A holistic approach combining multiple markers and integrating additional molecular and morphological data may offer a more comprehensive understanding of tick diversity and relationships. This research has far-reaching implications in formulating breeding programs and the development of vaccine against ticks and tick-borne diseases (TTBDs) as well as strategies for the management of resistant ticks.

Development and characterization of expressed sequence tag-simple sequence repeat markers for Anaphalis margaritacea var. yedoensis (Asteraceae)

RNA-sequencing was used to develop 16 microsatellite markers for the pearly everlasting, Anaphalis margaritacea var. yedoensis (Franch. et Sav.) Ohwi (Asteraceae), which inhabits gravel bars throughout the Japanese archipelago. The mean number of alleles for these 16 markers in two populations in the Hokkaido and Shizuoka Prefectures, was 3.5 and 4.0, respectively, while the mean expected heterozygosity was 0.525 and 0.560, respectively, with a significant genetic differentiation between the two populations. All markers could also be amplified in two conspecific taxa, A. margaritacea var. margaritacea and var. angustifolia, whereas 11 of the 16 markers were amplifiable in two congeneric species, A. sinica and A. alpicola. These newly developed microsatellite markers will support understanding of population genetics and mating systems in A. margaritacea var. yedoensis, and several will potentially be of use in similar studies in other Anaphalis species.

Validation of candidate gene-based EST-SSR markers for sugar yield in sugarcane

Sugarcane (Saccharum spp.) is a widely cultivated crop that fulfils approximately 75% of the sucrose demand worldwide. Owing to its polyploidy and complex genetic nature, it is difficult to identify and map genes related to complex traits, such as sucrose content. However, association mapping is one of the alternatives for identifying genes or markers for marker-assisted selection. In the present study, EST-SSR primers were obtained from in silico studies. The functionality of each primer was tested using Blast2Go software, and 30 EST-SSR primers related to sugar content were selected. These markers were validated using association analysis. A total of 70 F1 diverse genotypes for sugar content were phenotypes with two check lines. All parameters related to sugar content were recorded. The results showed a significant variation between the genotypes for sugar yield traits such as Brix value, purity, and sucrose content, etc. Correlation studies revealed that the Brix%, sucrose content, and sucrose recovery were significantly correlated. An association analysis was performed using mixed linear model to avoid false positive associations. The association analysis revealed that the SEM 407 marker was significantly associated with Brix% and sucrose content. The SEM 407 primers are putatively related to diphosphate-fructose-6-phosphate 1-phosphotransferase which is associated with Brix% and sucrose content. This functional marker can be used for marker-assisted selection for sugar yield traits in sugarcane that could accelerate the sugarcane breeding program.

Transcriptome analysis of Haemaphysalis flava female using Illumina HiSeq 4000 sequencing: de novo assembly, functional annotation and discovery of SSR markers

BACKGROUND

Ticks are ectoparasites capable of directly damaging their hosts and transmitting vector-borne diseases. The ixodid tick Haemaphysalis flava has a broad distribution that extends from East to South Asia. This tick is a reservoir of severe fever with thrombocytopenia syndrome virus (SFTSV) that causes severe hemorrhagic disease, with cases reported from China, Japan and South Korea. Recently, the distribution of H. flava in South Korea was found to overlap with the occurrence of SFTSV.

METHODS

This study was undertaken to discover the molecular resources of H. flava female ticks using the Illumina HiSeq 4000 system, the Trinity de novo sequence assembler and annotation against public databases. The locally curated Protostome database (PANM-DB) was used to screen the putative adaptation-related transcripts classified to gene families, such as angiotensin-converting enzyme, aquaporin, adenylate cyclase, AMP-activated protein kinase, glutamate receptors, heat shock proteins, molecular chaperones, insulin receptor, mitogen-activated protein kinase and solute carrier family proteins. Also, the repeats and simple sequence repeats (SSRs) were screened from the unigenes using RepeatMasker (v4.0.6) and MISA (v1.0) software tools, followed by the designing of SSRs flanking primers using BatchPrimer 3 (v1.0) software.

RESULTS

The transcriptome produced a total of 69,822 unigenes, of which 46,175 annotated to the homologous proteins in the PANM-DB. The unigenes were also mapped to the EuKaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) specializations. Promiscuous presence of protein kinase, zinc finger (C2H2-type), reverse transcriptase, and RNA recognition motif domains was observed in the unigenes. A total of 3480 SSRs were screened, of which 1907 and 1274 were found as tri- and dinucleotide repeats, respectively. A list of primer sequences flanking the SSR motifs was detailed for validation of polymorphism in H. flava and the related tick species.

CONCLUSIONS

The reference transcriptome information on H. flava female ticks will be useful for an enriched understanding of tick biology, its competency to act as a vector and the study of species diversity related to disease transmission.

Genome-wide survey and characterization of microsatellites in cashew and design of a web-based microsatellite database: CMDB

The cashew is an edible tree nut crop having a wide range of food and industrial applications. Despite great economic importance, the genome-wide characterization of microsatellites [simple sequence repeats (SSRs)] in cashew is lacking. In this study, we carried out the first comprehensive genome-wide microsatellites/SSRs characterization in cashew and developed polymorphic markers and a web-based microsatellite database. A total of 54526 SSRs were discovered in the cashew genome, with a mean frequency of 153 SSRs/Mb. Among the mined genome-wide SSRs (2-6 bp size motifs), the dinucleotide repeat motifs were dominant (68.98%) followed by the trinucleotides (24.56%). The Class I type of SSRs (≥20 bp) were 45.10%, while Class II repeat motifs (≥12-<20 bp) were 54.89% of the total genomic SSRs discovered here. Further, the AT-rich SSRs occurred more frequently in the cashew genome (84%) compared to the GC-rich SSRs. The validation of the in silico-mined genome-wide SSRs by PCR screening in cashew genotypes resulted in the development of 59 polymorphic SSR markers, and the polymorphism information content (PIC) of the polymorphic SSR markers ranged from 0.19 to 0.84. Further, a web-based database, "Cashew Microsatellite Database (CMDB)," was constructed to provide access to the genome-wide SSRs mined in this study as well as transcriptome-based SSRs from our previous study to the research community through a user-friendly searchable interface. Besides, CMDB provides information on experimentally validated SSRs. CMDB permits the retrieval of SSR markers information with the customized search options. Altogether, the genome-wide SSRs characterization, the polymorphic markers and CMDB database developed in this study would serve as valuable marker resources for DNA fingerprinting, germplasm characterization, genetic studies, and molecular breeding in cashew and related Anacardium species.

Applying molecular and genetic methods to trees and their fungal communities

Forests provide invaluable economic, ecological, and social services. At the same time, they are exposed to several threats, such as fragmentation, changing climatic conditions, or increasingly destructive pests and pathogens. Trees, the inherent species of forests, cannot be viewed as isolated organisms. Manifold (micro)organisms are associated with trees playing a pivotal role in forest ecosystems. Of these organisms, fungi may have the greatest impact on the life of trees. A multitude of molecular and genetic methods are now available to investigate tree species and their associated organisms. Due to their smaller genome sizes compared to tree species, whole genomes of different fungi are routinely compared. Such studies have only recently started in forest tree species. Here, we summarize the application of molecular and genetic methods in forest conservation genetics, tree breeding, and association genetics as well as for the investigation of fungal communities and their interrelated ecological functions. These techniques provide valuable insights into the molecular basis of adaptive traits, the impacts of forest management, and changing environmental conditions on tree species and fungal communities and can enhance tree-breeding cycles due to reduced time for field testing. It becomes clear that there are multifaceted interactions among microbial species as well as between these organisms and trees. We demonstrate the versatility of the different approaches based on case studies on trees and fungi. KEY POINTS: • Current knowledge of genetic methods applied to forest trees and associated fungi. • Genomic methods are essential in conservation, breeding, management, and research. • Important role of phytobiomes for trees and their ecosystems.

Genome survey sequencing-based SSR marker development and their validation in Dendrocalamus longispathus

Bamboo is an important genetic resource in India, supporting rural livelihood and industries. Unfortunately, most Indian bamboo taxa are devoid of basic genomic or marker information required to comprehend the genetic processes for further conservation and management. In this study, we perform genome survey sequencing for development of de novo genomic SSRs in Dendrocalamus longispathus, a socioeconomically important bamboo species of northeast India. Using Illumina platform, 69.49 million raw reads were generated and assembled into 1,145,321 contig with GC content 43% and N50 1228 bp. In total, 46,984 microsatellite repeats were mined-out wherein di-nucleotide repeats were most abundant (54.71%) followed by mono- (31.91%) and tri-repeats (9.85%). Overall, AT-rich repeats were predominant in the genome, but GC-rich motifs were more frequent in tri-repeats. Afterwards, 21,596 SSR loci were successfully tagged with the primer pairs, and a subset of 50 were validated through polymerase chain reaction amplification. Of these, 36 SSR loci were successfully amplified, and 16 demonstrated polymorphism. Using 13 polymorphic SSRs, a moderate level of gene diversity (He = 0.480; Ar = 3.52) was recorded in the analysed populations of D. longispathus. Despite the high gene flow (Nm = 4.928) and low genetic differentiation (F_ST = 0.119), severe inbreeding (F_IS = 0.407) was detected. Further, genetic clustering and STRUCTURE analysis revealed that the entire genetic variability is captured under two major gene pools. Conclusively, we present a comprehensive set of novel SSR markers in D. longispathus as well as other taxa of tropical woody bamboos.

Maximizing genetic gain through unlocking genetic variation in different ecotypes of kalmegh (Andrographis paniculata (Burm. f.) Nee)

Andrographis paniculata, commonly known as kalmegh is among the most popular medicinal herbs in Southeast Asia. It is widely cultivated for medicinal purposes. The bioactive molecule, Andrographolide accumulated in herb leaves has immense therapeutic and economic potential. However, comprehensive information regarding genetic diversity is very limited in this species. The present study assessed genetic diversity between and within the six populations (ecotypes) of twenty-four kalmegh accessions using multiple datasets (agro-morphological traits, phytochemical traits, and genic markers). This is the established report where EST-SSR (Expressed sequence tags-Simple Sequence Repeat) markers have been used to unlock genetic variation in kalmegh. Here, we identified and developed ninety-one metabolic pathway-specific EST-SSR markers. Finally, 32 random EST-SSR primer pairs were selected for genetic diversity assessment. Multivariate analysis to unveil the agro-morphological, phytochemical and genotypic variability was helpful in discriminating various germplasms studied in the present study. Among all the morphological discriptors used in present study, days to fifty percent flowering and dry herb yield were found as potential selection index for AP genetic improvement. Hierarchical cluster analysis built with agro-morphological data identified three major groups. However, corresponding analysis with phytochemical and molecular data generated two clear-cut groups among the studied individuals. Moreover, the grouping of individuals into different clusters using multiple datasets was geographically independent, and also showed inconsistency in grouping among agromorphological, phytochemical and molecular dataset based clusters. However, joint analysis using agro-morphological, phytochemical and genotypic information generated two genetic groups, which could be a valuable resource for identifying complementary crossing panels in the kalmegh breeding program. The accessions AP7, AP13, AP5, AP3 belong to cluster I and accessions AP17, AP18 belong to cluster II could be utilized as potential donors for high dry herb yield and andrographolide content, respectively in different selective breeding programs of AP. Thus, our results provided useful information about the overall genetic diversity and variation in economic traits useful for initiating selective breeding programs for contrasting traits of interest and maximizing genetic gain in kalmegh.

De novo transcriptome assembly using Illumina sequencing and development of EST-SSR markers in a monoecious herb Sagittaria trifolia Linn

Background

Sagittaria trifolia Linn. is a widespread macrophyte in Asia and southeast Europe and cultivated in parts of Asia. Although a few genomic studies have been conducted for S. trifolia var. sinensis, a crop breed, there is limited genomic information on the wild species of S. trifolia. Effective microsatellite markers are also lacking.

Objective

To assemble transcriptome sequence and develop effective EST-SSR markers for S. trifolia.

Methods

Here we developed microsatellite markers based on tri-, tetra-, penta-, and hexa-nucleotide repeat sequences by comparatively screening multiple transcriptome sequences of eleven individuals from ten natural populations of S. trifolia.

Results

A total of 107,022 unigenes were de novo assembled, with a mean length of 730 bp and an N50 length of 1,378 bp. The main repeat types were mononucleotide, trinucleotide, and dinucleotide, accounting for 55.83%, 23.51%, and 17.56% of the total repeats, respectively. A total of 86 microsatellite loci were identified with repeats of tri-, tetra-, penta-, and hexa-nucleotide. For SSR verification, 28 polymorphic loci from 41 randomly picked markers were found to produce stable and polymorphic bands, with the number of alleles per locus ranging from 2 to 11 and a mean of 5.2. The range of polymorphic information content (PIC) of each SSR locus varied from 0.25 to 0.80, with an average of 0.58. The expected heterozygosity ranged from 0.29 to 0.82, whereas the observed heterozygosity ranged from 0.25 to 0.90.

Conclusion

The assembled transcriptome and annotated unigenes of S. trifolia provide a basis for future studies on gene functions, pathways, and molecular mechanisms associated with this species and other related. The newly developed EST-SSR markers could be effective in examining population genetic structure, differentiation, and parentage analyses in ecological and evolutionary studies of S. trifolia.

Genetic Diversity Study on Geographical Populations of the Multipurpose Species Elsholtzia stauntonii Using Transferable Microsatellite Markers

Elsholtzia stauntonii Benth. (Lamiaceae) is an economically important ornamental, medicinal and aromatic plant species. To meet the increasing market demand for E. stauntonii, it is necessary to assess genetic diversity within the species to accelerate the process of genetic improvement. Analysis of the transferability of simple sequence repeat (SSR) markers from related species or genera is a fast and economical method to evaluate diversity, and can ensure the availability of molecular markers in crops with limited genomic resources. In this study, the cross-genera transferability of 497 SSR markers selected from other members of the Lamiaceae (Salvia L., Perilla L., Mentha L., Hyptis Jacq., Leonurus L., Pogostemon Desf., Rosmarinus L., and Scutella L.) to E. stauntonii was 9.05% (45 primers). Among the 45 transferable markers, 10 markers revealed relatively high polymorphism in E. stauntonii. The genetic variation among 825 individuals from 18 natural populations of E. stauntonii in Hebei Province of China was analyzed using the 10 polymorphic SSR markers. On the basis of the SSR data, the average number of alleles (N _A), expected heterozygosity (H _E), and Shannon's information index (I) of the 10 primers pairs were 7.000, 0.478, and 0.688, respectively. Lower gene flow (N _m = 1.252) and high genetic differentiation (F _st = 0.181) were detected in the populations. Analysis of molecular variance (AMOVA) revealed that most of the variation (81.47%) was within the populations. Integrating the results of STRUCTURE, UPGMA (Unweighted Pair Group Method with Arithmetic Mean) clustering, and principal coordinate analysis, the 825 samples were grouped into two clusters associated with geographical provenance (southwestern and northeastern regions), which was consistent with the results of a Mantel test (r = 0.56, p < 0.001). Overall, SSR markers developed in related genera were effective to study the genetic structure and genetic diversity in geographical populations of E. stauntonii. The results provide a theoretical basis for conservation of genetic resources, genetic improvement, and construction of a core collection for E. stauntonii.

Basic Characteristics of Flower Transcriptome Data and Derived Novel EST-SSR Markers of Luculia yunnanensis, an Endangered Species Endemic to Yunnan, Southwestern China

Luculia yunnanensis (Rubiaceae), an evergreen shrub or small tree, is endemic to China and confined to Nujiang Prefecture, Yunnan Province. This plant is of high ornamental value owing to its attractive pink flowers, sweet fragrance, and long flowering period. Due to the influence of climate change and human factors, the distribution range of L. yunnanensis has exhibited a significant shrinking trend, and it has become a vulnerable species that is in urgent need of conservation and rational utilization research. In this study, the flower transcriptome sequencing of L. yunnanensis was conducted using an Illumina HiSeq platform. We designed and developed a series of EST-SSR primers based on the flower transcriptome data of L. yunnanensis. The results showed that 98,389 unigenes were obtained from the L. yunnanensis flower transcriptome, all of which were aligned with sequences in public databases. Nr, Nt, Pfam, KOG/COG, Swiss-Prot, KEGG, and GO annotated 31,859, 13,853, 22,684, 10,947, 21,416, 9722, and 23,390 unigenes, respectively. The MISA (Microsatellite) tool was used to identify SSR loci from all unigenes, and a total of 15,384 SSRs were identified. Repeat motifs were given priority with mononucleotides, dinucleotides, and trinucleotides. The 81 primer pairs were synthesized randomly, of which 44 pairs showed effective amplification. A total of 17 primers showed stable amplification, and rich polymorphism was observed in 6 populations. We concluded via genetic diversity analysis that the average effective number of alleles (Ne), Shannon’s information index (I), and polymorphism information content (PIC) were 1.925, 0.837, and 0.403, respectively. In conclusion, 17 EST-SSR primers can be used for subsequent population genetic diversity analysis and molecular-marker-assisted breeding, which is of great significance for formulating resource conservation and utilization strategies for L. yunnanensis.

Plastome structure, phylogenomics and evolution of plastid genes in Swertia (Gentianaceae) in the Qing-Tibetan Plateau

BACKGROUND

The genus Swertia is of great medicinal importance and one of the most taxonomically challenging taxa within Gentianaceae, largely due to the morphological similarities of species within this genus and with its closely related genera. Previous molecular studies confirmed its polyphyly but suffered from low phylogenetic resolutions because only limited sequence loci were used. Thus, we conducted the structural, gene evolutionary, and phylogenetic analyses of 11 newly obtained plastomes of Swertia. Our result greatly improved the phylogenetic resolutions in Swertia, shed new light on the plastome evolution and phylogenetic relationships of this genus.

RESULTS

The 11 Swertia plastomes together with the published seven species proved highly similar in overall size, structure, gene order, and content, but revealed some structural variations caused by the expansion and contraction of the IRb region into the LSC region, due to the heterogeneous length of the ψycf1. The gene rps16 was found to be in a state flux with pseudogenes or completely lost. Similar situation was also documented in other genera of Gentianaceae. This might imply loss of the gene in the common ancestor of Gentianaceae. The distribution plot of ENC vs. GC3 showed all these plastomes arranging very close in the Wright line with an expected ENC value (49-52%), suggesting the codon usage of Swertia was mainly constrained by a GC mutation bias. Most of the genes remained under the purifying selection, however, the cemA was identified under positive selection, possibly reflecting an adaptive response to low CO₂ atmospheric conditions during the Late Miocene. Our phylogenomic analyses, based on 74 protein-coding genes (CDS), supported the polyphyly of Swertia with its close allies in the subtribe Swertiinae, presumably due to recent rapid radiation. The topology inferred from our phylogenetic analyses partly supported the current taxonomic treatment. Finally, several highly variable loci were identified, which can be used in future phylogenetic studies and accurate identification of medicinal genuineness of Swertia.

CONCLUSIONS

Our study confirmed the polyphyly of Swertia and demonstrated the power of plastome phylogenomics in improvement of phylogenetic resolution, thus contributing to a better understanding of the evolutionary history of Swertia.

Comparative plastid genomics of Mazaceae: focusing on a new recognized genus, Puchiumazus

Six Mazaceae plastomes were assembled in this study, including the newly recognized genus, Puchiumazus. Comparative plastid genomic analysis provided new insights into Mazaceae. The phylogenetic categorization of Mazus lanceifolius (Mazaceae) has long been uncertain. In 2021, the scholars Bo Li, D. G. Zhang, and C. L. Xiang republished M. lanceifolius as a new species Puchiumazus lanceifolius, within a new genus Puchiumazus. However, there is little plastome information on Mazaceae. Following the publishing of the new genus Puchiumazus, it is now necessary to study the Mazaceae plastome features to comprehensively understand this young family. The Mazaceae plastomes all have a typical quartile structure. The plastomes have a size ranging from 152,388 to 154,252 bp, and each plastome contains 112 unique genes, including 78 protein-coding genes, 4 rRNA genes, and 30 tRNA genes. A comparative analysis showed that these plastome sequences are highly conserved. Furthermore, we identified four relatively hypervariable regions (trnQ-UUC-psbK, trnS-GCU- trnS-CGA, trnT-UGU-trnL-UAA and ycf1) that can be used as potential DNA barcodes for the identification of this clade. Phylogenetic relationships based on the whole plastome sequences of 25 samples of 14 genera of Lamiales placed M. lanceifolius in the basal clade of the family Mazaceae, with 100% bootstrap support. In summary, the M. lanceifolius results indicate that a new monotype genus (Puchiumazus) should be established at the whole-plastome level. This study provides plastid genomic resources for exploring the phylogeny of Mazaceae.

Development of a new set of genic SSR markers in the genus Gentiana: in silico mining, characterization and validation

Gentiana is an important genus of around 360 medicinally important species, majority of which are not well characterized. Despite its importance, very few genomic resources are available for Gentiana L. Till date, the number of informative and robust simple sequence repeat (SSR)-based markers is limited and very few efforts have been made for their development. A set of robust, freely accessible and informative SSR markers for Gentiana is a pre-requisite for any molecular systematic as well as improvement studies in this group of pharmacologically valuable plants. In view of the importance of these plants, Expressed Sequence Tag (EST) sequences of 18 Gentiana species were surveyed for the development of a large set of non-redundant SSR markers. A total of 5808 perfect SSR with an average length of 17 bp and relative abundance of 214 loci/Mb were identified in the analysed 47,487 EST sequences using Krait software. Mapping of the ESTs resulted in gene ontology annotations of 49.14% of the sequences. Based on these perfect SSRs, 2902 primer pairs were designed, and 60 markers were randomly selected and validated on a set of Gentiana kurroo Royle accessions. Among the screened markers, 39 (65%) were found to be cross-species transferable. This is the first report of the largest set of functional, novel genic SSR markers in Gentiana, which will be a valuable resource for future characterization, genotype identification, conservation and genomic studies in the various species of this group of important medicinal plants.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02969-4.

Analysis of genetic diversity and relationships of Perilla frutescens using novel EST-SSR markers derived from transcriptome between wild-type and mutant Perilla

BACKGROUND

Perilla frutescens (Lamiaceae) is distributed in East Asia and is classified into var. frutescens and crispa. P. frutescens is multipurpose crop for human health because of a variety of secondary metabolites such as phenolic compound and essential oil. However, a lack of genetic information has hindered the development and utilization of Perilla genotypes.

METHODS AND RESULTS

This study was performed to develop expressed sequence tag-simple sequence repeat (EST-SSR) markers from P. frutescens var. crispa (wild type) and Antisperill (a mutant cultivar) and used them to assess the genetic diversity of, and relationships among, 94 P. frutescens genotypes. We obtained 65 Gb of sequence data comprising 632,970 transcripts by de novo RNA-sequencing. Of the 14,780 common SSRs, 102 polymorphic EST-SSRs were selected using in silico polymerase chain reaction (PCR). Overall, successful amplification from 58 EST-SSRs markers revealed remarkable genetic diversity and relationships among 94 P. frutescens genotypes. In total, 268 alleles were identified, with an average of 4.62 alleles per locus (range 2-11 alleles/locus). The average polymorphism information content (PIC) value was 0.50 (range 0.04-0.86). In phylogenetic and population structure analyses, the genotypes formed two major groups: Group I (var. crispa) and Group II (var. frutescens).

CONCLUSION

This results suggest that 58 novel EST-SSR markers derived from wild-type cultivar (var. crispa) and its mutant cultivar (Antisperill) have potential uses for population genetics and recombinant inbred line mapping analyses, which will provide comprehensive insights into the genetic diversity and relationship of P. frutescens.

Genetic Diversity and Structure of Rear Edge Populations of Sorbus aucuparia (Rosaceae) in the Hyrcanian Forest

Sorbus aucuparia (Rosaceae) is a small tree species widely distributed in Eurasia. The Hyrcanian forest is the southernmost distribution limit of this species. Severe habitat degradation and inadequate human interventions have endangered the long-term survival of this species in this region, and it is necessary to develop and apply appropriate management methods to prevent the loss of its genetic diversity. In this study, we used 10 SSR markers in order to evaluate the genetic diversity of this taxon. Leaf samples were collected from five known populations of S. aucuparia throughout its distribution area in the Hyrcanian forest. Expected heterozygosity ranged from 0.61 (ASH) to 0.73, and according to the M-ratio, all populations showed a significant reduction in effective population size, indicating a genetic bottleneck. Global FST was not statistically significant and attained the same values with and without excluding null alleles (ENA) correction (FST = 0.12). Bayesian analysis performed with STRUCTURE defined two genetic clusters among the five known populations, while the results of discriminant analysis of principal components (DAPC) identified three distinct groups. The average proportion of migrants was 22. In general, the gene flow was asymmetrical, with the biggest differences between immigration and emigration in Barzekoh and Asbehriseh. The Mantel test showed that there was no significant correlation between genetic distance (FST) and geographic distance in S. aucuparia. The best pathway for theoretical gene flow is located across the coast of the Caspian Sea and significant spatial autocorrelation was observed in only one population. In order to reduce the extinction risk of very small and scattered populations of S. aucuparia in the Hyrcanian forest, it is very important to establish and/or enhance the connectivity through habitat restoration or genetic exchange.

Protocol development for somatic embryogenesis, SSR markers and genetic modification of Stipagrostis pennata (Trin.) De Winter

BACKGROUND

Stipagrostis pennata (Trin.) De Winter is an important species for fixing sand in shifting and semi-fixed sandy lands, for grazing, and potentially as a source of lignocellulose fibres for pulp and paper industry. The seeds have low viability, which limits uses for revegetation. Somatic embryogenesis offers an alternative method for obtaining large numbers of plants from limited seed sources.

RESULTS

A protocol for plant regeneration from somatic embryos of S. pennata was developed. Somatic embryogenesis was induced on Murashige & Skoog (MS) medium supplemented with 3 mg·L-1 2,4-D subsequently shoots were induced on MS medium and supplemented with 5 mg·L-1 zeatin riboside. The highest shoots induction was obtained when embryogenic callus derived from mature embryos (96%) in combination with MS filter-sterilized medium was used from Khuzestan location. The genetic stability of regenerated plants was analysed using ten simple sequence repeats (SSR) markers from S. pennata which showed no somaclonal variation in regenerated plants from somatic embryos of S. pennata. The regenerated plants of S. pennata showed genetic stability without any somaclonal variation for the four pairs of primers that gave the expected amplicon sizes. This data seems very reliable as three of the PCR products belonged to the coding region of the genome. Furthermore, stable expression of GUS was obtained after Agrobacterium-mediated transformation using a super binary vector carried by a bacterial strain LBA4404.

CONCLUSION

To our knowledge, the current work is the first attempt to develop an in vitro protocol for somatic embryogenesis including the SSR marker analyses of regenerated plants, and Agrobacterium-mediated transformation of S. pennata that can be used for its large-scale production for commercial purposes.

Development of polymorphic microsatellites for genetic studies of white scar oyster (Crassostrea belcheri) using paired-end shotgun sequencing

White scar oyster Crassostrea belcheri is a commercially important bivalve species in Thailand. Appropriate genetic markers are needed for effective management to elevate its production efficiency. Type II microsatellites of C. belcheri were identified and characterized using an Illumina paired-end shotgun sequencing. A total of 14,743,710 reads were generated for which 198,849 reads containing microsatellites and 217,998 microsatellite loci were found. Twenty out of 60 microsatellite loci (33.33%) were polymorphic and these microsatellites were further tested against DNA bulks (N = 10 each) originating from 7 different geographic locations in Thai waters. Results indicated that newly developed microsatellites can be used for genetic diversity analysis of C. belcheri. Genotyping of C. belcheri collected from Surat Thani (Gulf of Thailand; N = 50) were performed. The number of alleles per locus ranged from 2 to 12 (average = 4.95). Observed and expected heterozygosities ranged from 0.0000 to 0.9400 (average = 0.3419) and 0.1139 to 0.8190 (average = 0.5844), respectively. Genome information and 20 newly isolated microsatellites will facilitate further studies in population genetics, stock management, and genetic improvement of C. belcheri in Thailand.

Genetic Distinctiveness but Low Diversity Characterizes Rear-Edge Thuja standishii (Gordon) Carr. (Cupressaceae) Populations in Southwest Japan

Rear-edge populations are of significant scientific interest because they can contain allelic variation not found in core-range populations. However, such populations can differ in their level of genetic diversity and divergence reflecting variation in life-history traits, demographic histories and human impacts. Using 13 EST-microsatellites, we investigated the genetic diversity and differentiation of rear-edge populations of the Japanese endemic conifer Thuja standishii (Gordon) Carr. in southwest Japan from the core-range in northeast Japan. Range-wide genetic differentiation was moderate (Fst = 0.087), with northeast populations weakly differentiated (Fst = 0.047), but harboring high genetic diversity (average population-level Ar = 4.76 and Ho = 0.59). In contrast, rear-edge populations were genetically diverged (Fst = 0.168), but contained few unique alleles with lower genetic diversity (Ar = 3.73, Ho = 0.49). The divergence between rear-edge populations exceeding levels observed in the core-range and results from ABC analysis and species distribution modelling suggest that these populations are most likely relicts of the Last Glacial Maximum. However, despite long term persistence, low effective population size, low migration between populations and genetic drift have worked to promote the genetic differentiation of southwest Japan populations of T. standishii without the accumulation of unique alleles.

Development and Characterization of 15 Novel Genomic SSRs for Viburnum farreri

The Viburnum genus is of particular interest to horticulturalists, phylogeneticists, and biogeographers. Despite its popularity, there are few existing molecular markers to investigate genetic diversity in this large genus, which includes over 160 species. There are also few polymorphic molecular tools that can delineate closely related species within the genus. Viburnum farreri, a member of the Solenotinus subclade and one of the centers of diversity for Viburnum, was selected for DNA sequencing and development of genomic simple sequence repeats (gSSRs). In this study, 15 polymorphic gSSRs were developed and characterized for a collection of 19 V. farreri samples. Number of alleles per locus ranged from two- to- eight and nine loci had four or more alleles. Observed heterozygosity ranged from 0 to 0.84 and expected heterozygosity ranged from 0.10 to 0.80 for the 15 loci. Shannon diversity index values across these loci ranged from 0.21 to 1.62. The markers developed in this study add to the existing molecular toolkit for the genus and will be used in future studies investigating cross-transferability, genetic variation, and species and cultivar delimitation in the Viburnum genus and closely allied genera in the Adoxaceae and Caprifoliaceae.

Insights from the cDNA and EST analysis of Antrodia cinnamomea

It is of interest to document the insights gleaned from the cDNA and EST analysis of Antrodia cinnamomea (a fungal species). Hence a library of sequences was constructed and analysed using standard procedures to gain new insights. Therefore, 65 ESTs, with size ranging from 300-2000 bp, were constructed. This included 46 ESTs with definite annotation, 18 ESTs were hypothetical and 1 new protein derived from BLAST analysis. We assigned 227 Gene Ontology terms linked to cell composition, transport, catalytic activity, and regulation functions in these sequences. Moreover, 56 matching genes were found in 8 Kyoto Encyclopedia of Genes and Genomes pathways. Data also showed 271 SSRs from Antrodia cinnamomea ESTs with an occurrence frequency of 96.82%. The STRING data analysis showed 29 genes encoded enzymes highly involved in protein-to-protein interactions linked to expression of regulation function. Thus, we documented some insights from the cDNA and EST analysis of Antrodia cinnamomea for further data mining.

Genetic structure of coast redwood (Sequoia sempervirens [D. Don] Endl.) populations in and outside of the natural distribution range based on nuclear and chloroplast microsatellite markers

Coast redwood (Sequoia sempervirens) naturally growing in southern Oregon and northern California is one of the few conifer tree species that are polyploid. Despite its unique ecological and economic importance, its population genetic structure is still insufficiently studied. To obtain additional data on its population genetic structure we genotyped 317 samples collected from populations in California (data set C) and 144 trees growing in a provenance trial in France (data set F) using 12 nuclear (five random nuclear genomic nSSRs and seven expressed sequence tag EST-SSRs) and six chloroplast (cpSSRs) microsatellite or simple sequence repeat (SSR) markers, respectively. These data sets were also used as reference to infer the origin of 147 coast redwood trees growing in Germany (data set G). Coast redwood was introduced to Europe, including Germany as an ornamental species, decades ago. Due to its fast growth and high timber quality, it could be considered as a potential commercial timber species, especially in perspective to climate warming that makes more regions in Germany suitable for its growing. The well performing trees in colder Germany could be potential frost resistant genotypes, but their genetic properties and origin are mostly unknown. Within the natural range in southern Oregon and northern California, only two relatively weak clusters were identified, one northern and one southern, separated by the San Francisco Bay. High genetic diversity, but low differentiation was found based on the 12 nuclear SSR markers for all three data sets F, C and G. We found that investigated 147 German trees represented only 37 different genotypes. They showed genetic diversity at the level less than diversity observed within the natural range in the northern or southern cluster, but more similar to the diversity observed in the southern cluster. It was difficult to assign German trees to the original single native populations using the six cpSSR markers, but rather to either the northern or southern cluster. The high number of haplotypes found in the data sets based on six cpSSR markers and low genetic differentiation based on 12 nuclear SSRs found in this study helps us study and better understand population genetic structure of this complex polyploid tree and supports the selection of potential genotypes for German forestry.

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.

Transcriptome Sequencing and Development of Novel Genic SSR Markers From Pistacia vera L

In this study, we aimed to develop novel genic simple sequence repeat (eSSR) markers and to study phylogenetic relationship among Pistacia species. Transcriptome sequencing was performed in different tissues of Siirt and Atl cultivars of pistachio (Pistacia vera). A total of 37.5-Gb data were used in the assembly. The number of total contigs and unigenes was calculated as 98,831, and the length of N50 was 1,333 bp after assembly. A total of 14,308 dinucleotide, trinucleotide, tetranucleotide, pentanucleotide, and hexanucleotide SSR motifs (4–17) were detected, and the most abundant SSR repeat types were trinucleotide (29.54%), dinucleotide (24.06%), hexanucleotide (20.67%), pentanucleotide (18.88%), and tetranucleotide (6.85%), respectively. Overall 250 primer pairs were designed randomly and tested in eight Pistacia species for amplification. Of them, 233 were generated polymerase chain reaction products in at least one of the Pistacia species. A total of 55 primer pairs that had amplifications in all tested Pistacia species were used to characterize 11 P. vera cultivars and 78 wild Pistacia genotypes belonging to nine Pistacia species (P. khinjuk, P. eurycarpa, P. atlantica, P. mutica, P. integerrima, P. chinensis, P. terebinthus, P. palaestina, and P. lentiscus). A total of 434 alleles were generated from 55 polymorphic eSSR loci with an average of 7.89 alleles per locus. The mean number of effective allele was 3.40 per locus. Polymorphism information content was 0.61, whereas observed (Ho) and expected heterozygosity (He) values were 0.39 and 0.65, respectively. UPGMA (unweighted pair-group method with arithmetic averages) and STRUCTURE analysis divided 89 Pistacia genotypes into seven populations. The closest species to P. vera was P. khinjuk. P. eurycarpa was closer P. atlantica than P. khinjuk. P. atlantica–P. mutica and P. terebinthus–P. palaestina pairs of species were not clearly separated from each other, and they were suggested as the same species. The present study demonstrated that eSSR markers can be used in the characterization and phylogenetic analysis of Pistacia species and cultivars, as well as genetic linkage mapping and QTL (quantitative trait locus) analysis.

Species diversity and phylogeography of Cornus kousa (Asian dogwood) captured by genomic and genic microsatellites

Cornus kousa (Asian dogwood), an East Asia native tree, is the most economically important species of the dogwood genus, owing to its desirable horticultural traits and ability to hybridize with North America-native dogwoods. To assess the species genetic diversity and to better inform the ongoing and future breeding efforts, we assembled an herbarium and arboretum collection of 131 noncultivated C. kousa specimens. Genotyping and capillary electrophoresis analyses of our C. kousa collection with the newly developed genic and published nuclear genomic microsatellites permitted assessment of genetic diversity and evolutionary history of the species. Regardless of the microsatellite type used, the study yielded generally similar insights into the C. kousa diversity with subtle differences deriving from and underlining the marker used. The accrued evidence pointed to the species distinct genetic pools related to the plant country of origin. This can be helpful in the development of the commercial cultivars for this important ornamental crop with increased pyramided utility traits. Analyses of the C. kousa evolutionary history using the accrued genotyping datasets pointed to an unsampled ancestor population, possibly now extinct, as per the phylogeography of the region. To our knowledge, there are few studies utilizing the same gDNA collection to compare performance of genomic and genic microsatellites. This is the first detailed report on C. kousa species diversity and evolutionary history inference.

Development of polymorphic EST-SSR markers and their applicability in genetic diversity evaluation in Rhododendron arboreum

The genus Rhododendron, known for large impressive flowers is widely distributed throughout the world. Rhododendrons have limited genetic information, despite of comprising high species diversity, morphological overlap and weak genetic barrier. In present study, expressed sequence tag (EST) data from Rhododendron catawbiense Michx (Subgenus Hymenanthes, Section Ponticum) and Rhododendron mucronatum var. ripense (Makino) E.H. Wilson (Subgenus Tsutsusi, Section Tsutsusi) were utilized for mining and identification of the SSRs for genetic diversity analysis of R. arboreum Smith (Subgenus Tsutsusi, Section Tsutsusi). A total of 249 SSRs were developed from 1767 contigs. Di-nucleotide was found to be most abundant repeat followed by tri- and tetra-nucleotide repeats. The motif AG/CT was most common di-nucleotide motif (31.73%), whereas, AAC/GTT (8.43%), ACG/CGT (8.03%), AAG/CTT (7.23%) and AGG/CCT (6.43%) were most abundant tri-nucleotide repeat motif. Among these SSRs, 168 sequences were only fit into the criteria to design flanking primer pairs. A total of 30 randomly selected primer pairs were utilized for validation and genetic diversity study in 36 genotypes of R. arboreum collected from western Himalayan region. In aggregate, 26 SSR markers (86.66%) produced good and repeatable amplifications. Expected heterozygosity (H_E) ranged from 0.322 to 0.841 and observed heterozygosity (H_O) ranged from 0.327 to 1.000 and PIC value ranged from 0.008 to 0.786. These primers were able to distinguish the geographic differences of occurrence based on cluster analysis. These developed EST-SSRs can be useful in future population genetics analysis and micro-evolutionary studies in Rhododendron species.

Cross-amplification of microsatellite markers across agarwood-producing species of the Aquilarieae tribe (Thymelaeaceae)

Tree species in the Aquilarieae tribe of the Thymelaeaceae family produce agarwood, a natural product highly valued for its fragrance, but the species are under threat due to indiscriminate harvesting. For conservation of these species, molecular techniques such as DNA profiling have been used. In this study, we assessed cross-amplification of microsatellite markers, initially developed for three Aquilaria species (A. crassna, A. malaccensis, and A. sinensis), on ten other agarwood-producing species, including members of Aquilaria (A. beccariana, A. hirta, A. microcarpa, A. rostrata, A. rugosa, A. subintegra, and A. yunnanensis) and Gyrinops (G. caudata, G. versteegii, and G. walla), both from the Aquilarieae tribe. Primers for 18 out of the 30 microsatellite markers successfully amplified bands of expected sizes in 1 sample each of at least 10 species. These were further used to genotype 74 individuals representing all the 13 studied species, yielding 13 cross-amplifiable markers, of which only 1 being polymorphic across all species. At each locus, the number of alleles ranged from 7 to 23, indicating a rather high variability. Four markers had relatively high species discrimination power. Our results demonstrated that genetic fingerprinting can be an effective tool in helping to manage agarwood genetic resources by potentially supporting the chain-of-custody of agarwood and its products in the market.

Genetic Divergence and Relationship Among Opisthopappus Species Identified by Development of EST-SSR Markers

Opisthopappus Shih is an endemic and endangered genus restricted to the Taihang Mountains that has important ornamental and economic value. According to the Flora Reipublicae Popularis Sinicae (FRPS, Chinese version), this genus contains two species (Opisthopappus longilobus and Opisthopappus taihangensis), whereas in the Flora of China (English version) only one species O. taihangensis is present. The interspecific phylogenetic relationship remains unclear and undefined, which might primarily be due to the lack of specific molecular markers for phylogenetic analysis. For this study, 2644 expressed sequence tag-simple sequence repeats (EST-SSRs) from 33,974 unigenes using a de novo transcript assembly of Opisthopappus were identified with a distribution frequency of 7.78% total unigenes. Thereinto, mononucleotides (1200, 45.39%) were the dominant repeat motif, followed by trinucleotides (992, 37.52%), and dinucleotides (410, 15.51%). The most dominant trinucleotide repeat motif was ACC/GGT (207, 20.87%). Based on the identified EST-SSRs, 245 among 1444 designed EST-SSR primers were selected for the development of potential molecular markers. Among these markers, 63 pairs of primers (25.71%) generated clear and reproducible bands with expected sizes. Eventually, 11 primer pairs successfully amplified all individuals from the studied populations. Through the EST-SSR markers, a high level of genetic diversity was detected between Opisthopappus populations. A significant genetic differentiation between the O. longilobus and O. taihangensis populations was found. All studied populations were divided into two clusters by UPGMA, NJ, STRUCTURE, and PCoA. These results fully supported the view of the FRPS, namely, that O. longilobus and O. taihangensis should be regarded as two distinct species. Our study demonstrated that transcriptome sequences, as a valuable tool for the quick and cost-effective development of molecular markers, was helpful toward obtaining comprehensive EST-SSR markers that could contribute to an in-depth assessment of the genetic and phylogenetic relationships between Opisthopappus species.

Development and characterization of 20 novel EST-SSR markers for Pteroceltis tatarinowii, a relict tree in China

PREMISE

Pteroceltis tatarinowii (Ulmaceae), the only species of the genus Pteroceltis, is an endangered tree in China. Here, novel expressed sequence tag-simple sequence repeat (EST-SSR) markers were developed to illuminate its genetic diversity for conservation and assisted breeding.

METHODS AND RESULTS

Based on Illumina transcriptome data from P. tatarinowii, a total of 70 EST-SSR markers were initially designed and tested. Forty-eight of 70 loci (68.6%) were successfully amplified, of which 20 were polymorphic. The number of alleles per locus ranged from two to six, and the levels of observed and expected heterozygosity ranged from 0.018 to 0.781 and from 0.023 to 0.702, respectively. Additionally, cross-amplification was successful for 17 loci in two related species, Ulmus gaussenii and U. chenmoui.

CONCLUSIONS

These new EST-SSR markers are valuable transcriptomic resources for P. tatarinowii and will facilitate population genetics and molecular breeding of this species and its relatives in Ulmaceae.

Genetic Structure of Norway Spruce Ecotypes Studied by SSR Markers

Norway spruce is a widespread and economically highly important tree species in Central Europe which occurs there in different morphotypic forms (also known as ecotypes). Previously established common garden experiments indicated that the morphological differentiation is most likely genetically determined. The genetic structure of Norway spruce morphological variants might be an indicator (marker) of specific sustainability in forest ecosystems. In this study, we investigated 436 individuals from autochthonous populations belonging to three different ecotypes. The main aim was to evaluate a level of genetic intra and interpopulation diversity among the low, medium and high-elevation ecotypes using both expressed sequence tag simple sequence repeats (EST – SSR) and genomic SSR markers. Sixteen highly polymorphic microsatellite loci folded in two newly designed multiplexes were used to depicture the genetic structure of targeted trees. Important allele frequency parameters, such as the mean expected (0.722, SE = 0.061) and observed (0.585, SE = 0.062) heterozygosity and mean effective number of alleles (Ne = 5.943, SE = 1.279), were estimated. The low genetic differentiation detected among different ecotypes (Fst = 0.008) was further discussed and clarified.

Genetic Diversity and Population Genetic Structure of Cinnamomum camphora in South China Revealed by EST-SSR Markers

Cinnamomum camphora is a valuable broad-leaf tree indigenous to South China and East Asia and has been widely cultivated and utilized by humans since ancient times. However, owing to its overutilization for essential oil extraction, the Transplanting Big Trees into Cities Program, and over deforestation to make furniture, its wild populations have been detrimentally affected and are declining rapidly. In the present study, the genetic diversity and population structure of 180 trees sampled from 41 populations in South China were investigated with 22 expressed sequence tag-simple sequence repeat (EST-SSR) markers. In total, 61 alleles were harbored across 180 individuals, and medium genetic diversity level was inferred from the observed heterozygosity (Ho), expected heterozygosity (He), and Nei’ gene diversity (GD), which were 0.45, 0.44, and 0.44, respectively. Among the 41 wild populations, C. camphora had an average of 44 alleles, 2.02 effective alleles, and He ranging from 0.30 (SC) to 0.61 (HK). Analysis of molecular variance (AMOVA) showed that 17% of the variation among populations and the average pairwise genetic differentiation coefficient (FST) between populations was 0.162, indicating relatively low genetic population differentiations. Structure analysis suggested two groups for the 180 individuals, which was consistent with the principal coordinate analysis (PCoA) and unweighted pair-group method with arithmetic means (UPGMA). Populations grouped to cluster I were nearly all distributed in Jiangxi Province (except population XS in Zhejiang Province), and cluster II mainly comprised populations from other regions, indicating a significant geographical distribution. Moreover, the Mantel test showed that this geographical distance was significantly correlated with genetic distance. The findings of this research will assist in future C. camphora conservation management and breeding programs.

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.

Transcriptomic Analysis of Marine Gastropod Hemifusus tuba Provides Novel Insights into Conotoxin Genes

The marine gastropod Hemifusus tuba is served as a luxury food in Asian countries and used in traditional Chinese medicine to treat lumbago and deafness. The lack of genomic data on H. tuba is a barrier to aquaculture development and functional characteristics of potential bioactive molecules are poorly understood. In the present study, we used high-throughput sequencing technologies to generate the first transcriptomic database of H. tuba. A total of 41 unique conopeptides were retrieved from 44 unigenes, containing 6-cysteine frameworks belonging to four superfamilies. Duplication of mature regions and alternative splicing were also found in some of the conopeptides, and the de novo assembly identified a total of 76,306 transcripts with an average length of 824.6 nt, of which including 75,620 (99.1%) were annotated. In addition, simple sequence repeats (SSRs) detection identified 14,000 unigenes containing 20,735 SSRs, among which, 23 polymorphic SSRs were screened. Thirteen of these markers could be amplified in Hemifusus ternatanus and seven in Rapana venosa. This study provides reports of conopeptide genes in Buccinidae for the first time as well as genomic resources for further drug development, gene discovery and population resource studies of this species.

Transcriptome analysis and codominant markers development in caper, a drought tolerant orphan crop with medicinal value

Caper (Capparis spinosa L.) is a xerophytic shrub cultivated for its flower buds and fruits, used as food and for their medicinal properties. Breeding programs and even proper taxonomic classification of the genus Capparis has been hampered so far by the lack of reliable genetic information and molecular markers. Here, we present the first genomic resource for C. spinosa, generated by transcriptomic approach and de novo assembly. The sequencing effort produced nearly 80 million clean reads assembled into 124,723 unitranscripts. Careful annotation and comparison with public databases revealed homologs to genes with a key role in important metabolic pathways linked to abiotic stress tolerance and bio-compounds production, such purine, thiamine and phenylpropanoid biosynthesis, α-linolenic acid and lipid metabolism. Additionally, a panel of genes involved in stomatal development/distribution and encoding for Stress Associated Proteins (SAPs) was also identified. We also used the transcriptomic data to uncover novel molecular markers for caper. Out of 50 SSRs tested, 14 proved polymorphic and represent the first set of SSR markers for the genus Capparis. This transcriptome will be an important contribution to future studies and breeding programs for this orphan crop, aiding to the development of improved varieties to sustain agriculture in arid conditions.

Discovery of novel genic-SSR markers from transcriptome dataset of an important non-human primate, Macaca fascicularis

Macaca fascicularis, also known as the cynomolgus macaque, is an important non-human primate animal model used in biomedical research. It is an Old-World primate widely distributed in Southeast Asia and is one of the most abundant macaque species in Malaysia. However, the genetic structure of wild cynomolgus macaque populations in Malaysia has not been thoroughly elucidated. In this study, we developed genic-simple sequence repeat (genic-SSR) markers from an in-house transcriptome dataset generated from the Malaysian cynomolgus macaque via RNA sequencing, and applied these markers on 26 cynomolgus macaque individuals. A collection of 14,751 genic-SSRs were identified, where 13,709 were perfect SSRs. Dinucleotide repeats were the most common repeat motifs with a frequency of 65.05%, followed by trinucleotide repeats (20.55%). Subsequently, we designed 300 pairs of primers based on perfect di- and trinucleotide SSRs, in which 105 SSRs were associated with functional genes. A subset of 30 SSR markers were randomly selected and validated, yielding 19 polymorphic markers with an average polymorphism information content value of 0.431. The development of genic-SSR markers in this study is indeed timely to provide useful markers for functional and population genetic studies of the cynomolgus macaque and other related non-human primate species.

Development of EST-SSR markers in Saxifraga sinomontana (Saxifragaceae) and cross-amplification in three related species

PREMISE

Saxifraga sinomontana (Saxifragaceae) is a widespread alpine species in the Qinghai-Tibetan Plateau and its flanking mountains. We developed a set of expressed sequence tag-simple sequence repeat (EST-SSR) markers to investigate the genetic diversity and evolutionary history of the species.

METHODS AND RESULTS

We initially designed 50 EST-SSR markers based on transcriptome data of S. sinomontana. Nineteen of 50 loci (38%) were successfully amplified, 13 of which were polymorphic. These were tested on 71 individuals from four populations. Three to 18 alleles per locus were detected, and the levels of observed and expected heterozygosity ranged from 0.2817 to 0.9155 and 0.2585 to 0.8495, respectively. In addition, cross-amplification was successful for all 13 loci in three congeneric species, S. tangutica, S. heleonastes, and S. congestiflora.

CONCLUSIONS

These EST-SSR markers will be useful for studying the genetic diversity of S. sinomontana and disentangling the phylogenetic relationships of related species.

Comparative transcriptomic analysis reveals genetic divergence and domestication genes in Diospyros

BACKGROUND

Persimmon (Diospyros kaki) is the most economically cultivated species belonging to the genus Diospyros. However, little is known about the interspecific diversity and mechanism of domestication, partly due to the lack of genomic information that is available for closely related species of D. kaki (DK). Here, we performed transcriptome sequencing on nine samples, including DK, a variety of DK and seven closely related species, to evaluate the interspecific genetic divergence and to identify candidate genes involved in persimmon domestication.

RESULTS

We obtained a total of 483,421 unigenes with N50 at 1490 bp in the nine Diospyros samples and identified 2603 orthogroups that were shared among all the samples using OrthoMCL analysis. A phylogenetic tree was established based on the tandem 2603 one-to-one single copy gene alignments, showing that DK was closely related to D. kaki var. silvestris (DKV) and that it clustered with the clade of D. deyangnsis (DD) and was farthest from the D. cathayensis (DC) species. The nonsynonymous substitutions (Ka), via synonymous substitution (Ks) ratios, was directly proportional to the genetic relationship of the different species. The higher the Ka/Ks ratios, the longer the distance was. Moreover, 31 positively selected genes (PSGs) involved in carbohydrate metabolism and phenolic metabolism were identified and isolated, and nearly all PSGs except the MATE gene had a high expression in the DK or DKV species. It was hypothesized that these genes might contribute to the domestication of the DK species. Finally, we developed the expressed sequence tag-simple sequence repeat (EST-SSR) and identified 2 unique amplicons DKSSR10 and DKSSR39: the former was absent in the DC species but was present in the other species, the latter had a long amplification product in the DJ species.

CONCLUSION

This study presents the first transcriptome resources for the closely related species of persimmon and reveals interspecific genetic divergence. It is speculated that DK is derived from the hybridization of DD and DO species. Furthermore, our analysis suggests candidate PSGs that may be crucial for the adaptation, domestication, and speciation of persimmon relatives and suggests that DKSSR10 and DKSSSR39 could potentially serve as species-specific molecular markers.

Transcriptome Sequencing of Different Avocado Ecotypes: de novo Transcriptome Assembly, Annotation, Identification and Validation of EST-SSR Markers

Avocado (Persea americana Mill.) could be considered as an important tropical and subtropical woody oil crop with high economic and nutritional value. Despite the importance of this species, genomic information is currently unavailable for avocado and closely related congeners. In this study, we generated more than 216 million clean reads from different avocado ecotypes using Illumina HiSeq high-throughput sequencing technology. The high-quality reads were assembled into 154,310 unigenes with an average length of 922 bp. A total of 55,558 simple sequence repeat (SSR) loci detected among the 43,270 SSR-containing unigene sequences were used to develop 74,580 expressed sequence tag (EST)-SSR markers. From these markers, a subset of 100 EST-SSR markers was randomly chosen to identify polymorphic EST-SSR markers in 28 avocado accessions. Sixteen EST-SSR markers with moderate to high polymorphism levels were detected, with polymorphism information contents ranging from 0.33 to 0.84 and averaging 0.63. These 16 polymorphic EST-SSRs could clearly and effectively distinguish the 28 avocado accessions. In summary, our study is the first presentation of transcriptome data of different avocado ecotypes and comprehensive study on the development and analysis of a set of EST-SSR markers in avocado. The application of next-generation sequencing techniques for SSR development is a potentially powerful tool for genetic studies.

Mining and characterization of novel EST-SSR markers of Parrotia subaequalis (Hamamelidaceae) from the first Illumina-based transcriptome datasets

Parrotia subaequalis is an endangered Tertiary relict tree from eastern China. Despite its important ecological and horticultural value, no transcriptomic data and limited molecular markers are currently available in this species. In this study, we first performed high-throughput transcriptome sequencing of two individuals representing the northernmost (TX) and southernmost (SJD) population of P. subaequalis on the Illumina HiSeq 2500 platform. We gathered a total of 69,135 unigenes for P. subaequalis (TX) and 84,009 unigenes for P. subaequalis (SJD). From two unigenes datasets, 497 candidate polymorphic novel expressed sequence tag-simple sequence repeats (EST-SSRs) were identified using CandiSSR. Among these repeats, di-nucleotide repeats were the most abundant repeat type (62.78%) followed by tri-, tetra- and hexa-nucleotide repeats. We then randomly selected 54 primer pairs for polymorphism validation, of which 27 (50%) were successfully amplified and showed polymorphisms in 96 individuals from six natural populations of P. subaequalis. The average number of alleles per locus and the polymorphism information content values were 3.70 and 0.343; the average observed and expected heterozygosity were 0.378 and 0.394. A relatively high level of genetic diversity (H_T = 0.393) and genetic differentiation level (F_ST = 0.171) were surveyed, indicating P. subaequalis maintained high levels of species diversity in the long-term evolutionary history. Additionally, a high level of cross-transferability (92.59%) was displayed in five congeneric Hamamelidaceae species. Therefore, these new transcriptomic data and novel polymorphic EST-SSR markers will pinpoint genetic resources and facilitate future studies on population genetics and molecular breeding of P. subaequalis and other Hamamelidaceae species.

Development of EST-derived microsatellite markers to investigate the population structure of sparganum - the causative agent of zoonotic sparganosis

The plerocercoid (sparganum) of Spirometra erinaceieuropaei is the main aetiological agent of human sparganosis. To improve the current knowledge on S. erinaceieuropaei evolution, we performed multi-locus microsatellite typing of sparganum isolates from China for the first time. All available expressed sequence tag (EST) sequences for the Spirometra were downloaded from the GenBank. The identification and localization of microsatellites in ESTs was accomplished by MISA. Based on the selected microsatellites, the genetic structure of 64 sparganum isolates collected from 11 geographical locations in southwest China were investigated through principal component analysis, STRUCTURE analysis and neighbour-joining clustering. A total of 522 non-redundant ESTs containing 915 simple sequence repeats were identified from 12 481 ESTs screened. Five primer pairs were finally selected. Using these loci, a total of 12 alleles were detected in 64 sparganum isolates. Little variability was observed within each of geographical population, especially among isolates derived from Kunming of Yunnan (YN-KM) province. Both STRUCTURE analysis and the clustering analysis supported that two genotypes existed among the sparganum isolates from southwest China. In conclusion, five microsatellite markers were successfully developed, and sparganum population was observed to harbour low genetic variation, further investigation with deeper sampling was needed to elucidate the population structure.

Taraxacum kok-saghyz (rubber dandelion) genomic microsatellite loci reveal modest genetic diversity and cross-amplify broadly to related species

Taraxacum kok-saghyz (TKS) carries great potential as alternative natural rubber source. To better inform future breeding efforts with TKS and gain a deeper understanding of its genetic diversity, we utilized de novo sequencing to generate novel genomic simple sequence repeats markers (gSSRs). We utilized 25 gSSRs on a collection of genomic DNA (gDNA) samples from germplasm bank, and two gDNA samples from historical herbarium specimens. PCR coupled with capillary electrophoresis and an array of population genetics tools were employed to analyze the dataset of our study as well as a dataset of the recently published genic SSRs (eSSRs) generated on the same germplasm. Our results using both gSSRs and eSSRs revealed that TKS has low- to- moderate genetic diversity with most of it partitioned to the individuals and individuals within populations, whereas the species lacked population structure. Nineteen of the 25 gSSR markers cross-amplified to other Taraxacum spp. collected from Southeastern United States and identified as T. officinale by ITS sequencing. We used a subset of 14 gSSRs to estimate the genetic diversity of the T. officinale gDNA collection. In contrast to the obligatory outcrossing TKS, T. officinale presented evidence for population structure and clonal reproduction, which agreed with the species biology. We mapped the molecular markers sequences from this study and several others to the well-annotated sunflower genome. Our gSSRs present a functional tool for the biodiversity analyses in Taraxacum, but also in the related genera, as well as in the closely related tribes of the Asteraceae.

The Promise of Molecular and Genomic Techniques for Biodiversity Research and DNA Barcoding of the Arabian Peninsula Flora

The Arabian Peninsula is known to have a comprehensive and rich endowment of unique and genetically diverse plant genetic resources. Analysis and conservation of biological diversity is a crucial issue to the whole Arabian Peninsula. The rapid and accurate delimitation and identification of a species is crucial to genetic diversity analysis and the first critical step in the assessment of distribution, population abundance and threats related to a particular target species. During the last two decades, classical strategies of evaluating genetic variability, such as morphology and physiology, have been greatly complemented by phylogenetic, taxonomic, genetic diversity and breeding research molecular studies. At present, initiatives are taking place around the world to generate DNA barcode libraries for vascular plant flora and to make these data available in order to better understand, conserve and utilize biodiversity. The number of herbarium collection-based plant evolutionary genetics and genomics studies being conducted has been increasing worldwide. The herbaria provide a rich resource of already preserved and identified material, and these as well as freshly collected samples from the wild can be used for creating a reference DNA barcode library for the vascular plant flora of a region. This review discusses the main molecular and genomic techniques used in plant identification and biodiversity analysis. Hence, we highlight studies emphasizing various molecular techniques undertaken during the last 10 years to study the plant biodiversity of the Arabian Peninsula. Special emphasis on the role of DNA barcoding as a powerful tool for plant biodiversity analysis is provided, along with the crucial role of herbaria in creating a DNA barcode library.

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.

Microsatellite development from genome skimming and transcriptome sequencing: comparison of strategies and lessons from frog species

Background

Even though microsatellite loci frequently have been isolated using recently developed next-generation sequencing (NGS) techniques, this task is still difficult because of the subsequent polymorphism screening requires a substantial amount of time. Selecting appropriate polymorphic microsatellites is a critical issue for ecological and evolutionary studies. However, the extent to which assembly strategy, read length, sequencing depth, and library layout produce a measurable effect on microsatellite marker development remains unclear. Here, we use six frog species for genome skimming and two frog species for transcriptome sequencing to develop microsatellite markers, and investigate the effect of different isolation strategies on the yield of microsatellites.

Results

The results revealed that the number of isolated microsatellites increases with increased data quantity and read length. Assembly strategy could influence the yield and the polymorphism of microsatellite development. Larger k-mer sizes produced fewer total number of microsatellite loci, but these loci had a longer repeat length, suggesting greater polymorphism. However, the proportion of each type of nucleotide repeats was not affected; dinucleotide repeats were always the dominant type. Finally, the transcriptomic microsatellites displayed lower levels of polymorphisms and were less abundant than genomic microsatellites, but more likely to be functionally linked loci.

Conclusions

These observations provide deep insight into the evolution and distribution of microsatellites and how different isolation strategies affect microsatellite development using NGS.

Electronic supplementary material

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

High Morphological Differentiation in Crown Architecture Contrasts with Low Population Genetic Structure of German Norway Spruce Stands

High elevation sites in the low mountain ranges in Germany are naturally covered by Norway spruce (Picea abies (Karst.) L.) stands. Historically, large scale anthropogenic range expansion starting in the mid to late 18th century had a huge impact on the forest composition throughout Germany. Utilisation and exploitation often led to artificial regeneration, mostly carried out using seeds from allochthonous provenances. Usually, autochthonous (natural) high elevation Norway spruce trees have narrow crown phenotypes, whereas lowland trees have broader crowns. Narrow crown phenotypes are likely the result of adaptation to heavy snow loads combined with high wind speeds. In the present study, neighbouring stand pairs of putative autochthonous and allochthonous origin with contrasting phenotypes in high elevation sites were investigated with 200 samples each. These stands are located in the Ore Mountains, the Thuringian Forest, and the Harz Mountains. Additionally, a relict population with the typical narrow high elevation phenotypes was sampled in Thuringia, known as “Schlossbergfichte”. The objective of the study was to quantify supposedly adaptive phenotypic differences in crown architecture and the genetic differentiation of 11 putatively neutral nuclear microsatellite markers (i.e., simple sequence repeats (nSSRs)). The high differentiation of morphological traits (PST = 0.952–0.989) between the neighbouring autochthonous and allochthonous stands of similar age contrasts with the very low neutral genetic differentiation (FST = 0.002–0.007; G″ST = 0.002–0.030), suggesting that directional selection at adaptive gene loci was involved in phenotypic differentiation. Comparing the regions, a small isolation by distance effect for the Harz Mountains was detected, suggesting landscape resistance restricting gene flow. Finally, the differentiation of the very old autochthonous (up to 250 years) stand “Schlossbergfichte” with typical high elevation phenotypes could cohere with the sampling of a relict genepool.

Analysis of Environment-Marker Associations in American Chestnut

American chestnut (Castanea dentata Borkh.) was a dominant tree species in its native range in eastern North America until the accidentally introduced fungus Cryphonectria parasitica (Murr.) Barr, that causes chestnut blight, led to a collapse of the species. Different approaches (e.g., genetic engineering or conventional breeding) are being used to fight against chestnut blight and to reintroduce the species with resistant planting stock. Because of large climatic differences within the distribution area of American chestnut, successful reintroduction of the species requires knowledge and consideration of local adaptation to the prevailing environmental conditions. Previous studies revealed clear patterns of genetic diversity along the northeast-southwest axis of the Appalachian Mountains, but less is known about the distribution of potentially adaptive genetic variation within the distribution area of this species. In this study, we investigated neutral and potentially adaptive genetic variation in nine American chestnut populations collected from sites with different environmental conditions. In total, 272 individuals were genotyped with 24 microsatellite (i.e., simple sequence repeat (SSR)) markers (seven genomic SSRs and 17 EST-SSRs). An FST-outlier analysis revealed five outlier loci. The same loci, as well as five additional ones, were significantly associated with environmental variables of the population sites in an environmental association analysis. Four of these loci are of particular interest, since they were significant in both methods, and they were associated with environmental variation, but not with geographic variation. Hence, these loci might be involved in (temperature-related) adaptive processes in American chestnut. This work aims to help understanding the genetic basis of adaptation in C. dentata, and therefore the selection of suitable provenances for further breeding efforts.