Mining and survey of simple sequence repeats in expressed sequence tags of dicotyledonous species

Development and Validation of Gene-Based SSR Markers in the Genus Mesembryanthemum

Bioinformatics tools have been employed for the direct development of gene-based simple sequence repeat (SSR) markers. Through the analysis of 28,056 Mesembryanthemum expressed sequence tag (EST) sequences, a total of 5,851 ESTs containing SSRs were identified, amounting to approximately 17.07 Mb. Among these, 938 EST sequences harbored more than one SSR marker, and 788 EST-SSR sequences were found in compound form. The most prevalent types of SSR motifs were mononucleotide repeats (MNRs), accounting for 44%, followed by di-nucleotide repeats (DNRs) at 37%, and trinucleotide repeats (TNRs) at 16%. Notably, TNR or longer SSR motifs primarily consisted of shorter repeat lengths, with only 51 motifs containing 10 or more repeats. The BLASTX analysis successfully assigned functions to 4,623 (79%) of the EST sequences. Among the developed primer sets, 21 primers amplified a total of 65 alleles, with primer PMA79 EST-SSR exhibiting the maximum of six alleles. The polymorphic information content (PIC) values ranged from 0 to 0.76, with a mean of 0.47. The marker index (MI) and discriminating power (D) values reached 0.66 (primer PMA63) and 0.95 (primer PMA20), respectively. Utilizing the unweighted pair group method with arithmetic mean (UPGMA), a dendrogram was constructed, successfully segregating the 24 Mesembryanthemum genotypes into three distinct clusters, with a similarity coefficient ranging from 0.96 to 0.38. In this study, we have developed a total of 83 EST-SSR primer pairs specific to the Mesembryanthemum genus. These newly developed EST-SSRs will serve as valuable tools for researchers, particularly molecular breeders, enabling gene-based identification and trait selection through marker-assisted breeding approaches.

Genome-wide in silico characterization, validation, and cross-species transferability of microsatellite markers in Mallard and Muscovy ducks

BACKGROUND

Microsatellites are important markers for livestock including ducks. The development of microsatellites is expensive and labor-intensive. Meanwhile, the in silico approach for mining for microsatellites became a practicable alternative. Therefore, the current study aimed at comparing whole-genome and chromosome-wise microsatellite mining approaches in Muscovy and Mallard ducks and testing the transferability of markers between them. The GMATA software was used for the in silico study, and validation was performed using 26 primers.

RESULTS

The total number of the detected microsatellites using chromosome-wise was 250,053 and 226,417 loci compared to 260,059 and 238,462 loci using whole genome in Mallards and Muscovies. The frequencies of different motifs had similar patterns using the two approaches. Dinucleotide motifs were predominant (> 50%) in both Mallards and Muscovies. The amplification of the genomes revealed an average number of alleles of 5.08 and 4.96 in Mallards and Muscovies. One locus was monographic in Mallards, and two were monomorphic in Muscovies. The average expected heterozygosity was higher in Muscovy than in Mallards (0.45 vs. 0.43) with no significant difference between the two primer sets, which indicated the usefulness of cross-species amplification of different primers.

CONCLUSION

The current study developed a whole-genome SSR panel for ducks for the first time, and the results could prove that using chromosome-wise mining did not generate different results compared to the whole-genome approach.

Disentangling potential genotypes for macro and micro nutrients and polymorphic markers in Chickpea

The present investigation was conducted to assess the nutritional diverseness and identify novel genetic resources to be utilized in chickpea breeding for macro and micro nutrients. The plants were grown in randomized block design. Nutritional and phytochemical properties of nine chickpea genotypes were estimated. The EST sequences from NCBI database were downloaded in FASTA format, clustered into contigs using CAP3, mined for novel SSRs using TROLL analysis and primer pairs were designed using Primer 3 software. Jaccard's similarity coefficients were used to compare the nutritional and molecular indexes followed by dendrograms construction employing UPGMA approach. The genotypes PUSA-1103, K-850, PUSA-1108, PUSA-1053 and the EST-SSR markers including the 5 newly designed namely ICCeM0012, ICCeM0049, ICCeM0067, ICCeM0070, ICCeM0078, SVP55, SVP95, SVP96, SVP146, and SVP217 were found as potential donor/marker resources for the macro-micro nutrients. The genotypes differed (p < 0.05) for nutritional properties. Amongst newly designed primers, 6 were found polymorphic with median PIC (0.46). The alleles per primer ranged 1 to 8. Cluster analysis based on nutritional and molecular diversities partially matched to each other in principle. The identified novel genetic resources may be used to widen the germplasm base, prepare maintainable catalogue and identify systematic blueprints for future chickpea breeding strategies targeting macro-micro nutrients.

Development of unigene-derived SSR markers from RNA-seq data of Uraria lagopodioides (Fabaceae) and their application in the genus Uraria Desv. (Fabaceae)

BACKGROUND

Uraria Desv. belongs to the tribe Desmodieae (Fabaceae), a group of legume plants, some of which have medicinal properties. However, due to a lack of genomic information, the interspecific relationships, genetic diversity, population genetics, and identification of functional genes within Uraria species are still unclear.

RESULTS

Using RNA-Seq, a total of 66,026 Uraria lagopodioides unigenes with a total sequence content of 52,171,904 bp were obtained via de novo assembly and annotated using GO, KEGG, and KOG databases. 17,740 SSRs were identified from a set of 66,026 unigenes. Cross-species amplification showed that 54 out of 150 potential unigene-derived SSRs were transferable in Uraria, of which 19 polymorphic SSRs were developed. Cluster analysis based on polymorphisms successfully distinguished seven Uraria species and revealed their interspecific relationships. Seventeen samples of seven Uraria species were clustered into two monophyletic clades, and phylogenetic relationships of Uraria species based on unigene-derived SSRs were consistent with classifications based on morphological characteristics.

CONCLUSIONS

Unigenes annotated in the present study will provide new insights into the functional genomics of Uraria species. Meanwhile, the unigene-derived SSR markers developed here will be invaluable for assessing the genetic diversity and evolutionary history of Uraria and relatives.

De Novo Transcriptome Assembly and EST-SSR Marker Development and Application in Chrysosplenium macrophyllum

Chrysosplenium macrophyllum Oliv., belonging to the family Saxifragaceae, is a traditional and unique Chinese herbal medicine. However, the lack of adequate molecular markers has hampered the progress regarding population genetics and evolution within this species. In this research, we used the DNBSEQ-T7 Sequencer (MGI) sequencing assay to analyze the transcriptome profiles of C. macrophyllum. SSR markers were developed on the basis of transcriptomic sequences and further validated on C. macrophyllum and other Chrysosplenium species. The genetic diversity and structure of the 12 populations were analyzed by using polymorphic expressed sequence tag simple sequence repeat (EST-SSR) markers. A potential pool of 3127 non-redundant EST-SSR markers were identified for C. macrophyllum in this study. The developed EST-SSR markers had high amplification rates and cross-species transferability in Chrysosplenium. Our results also showed that the natural populations of C. macrophyllum had a high level of genetic diversity. Genetic distance, principal component analysis, and popular structure analysis revealed that all 60 samples clustered into two major groups that were consistent with their geographical origins. This study provided a batch of highly polymorphic EST-SSR molecular markers that were developed via transcriptome sequencing. These markers will be of great significance for the study of the genetic diversity and evolutionary history of C. macrophyllum and other Chrysosplenium species.

Transcriptome analysis and development of EST-SSR markers in Anoectochilus emeiensis (Orchidaceae)

Anoectochilus emeiensis K. Y. Lang, together with other Anoectochilus species, has long been used as the main source of many traditional Chinese medicines. Owing to the shortcomings of molecular markers, the study of the genetic diversity and medicinal component synthesis mechanism of the endemic Anoectochilus species has been delayed. In this study, I carried out a transcriptome analysis of A. emeiensis. A total of 78,381 unigenes were assembled from 64.2 million reads, and 47,541 (60.65%) unigenes were matched to known proteins in the public databases. Then, 9284 expressed sequence tag-derived simple sequence repeats (EST-SSRs) were identified, and the frequency of SSRs in the A. emeiensis transcriptome was 9.88%. Trinucleotide repeats (3699, 39.84%) were the most common type, followed by dinucleotide (3251, 35.02%) and mononucleotide (1750, 18.85%) repeats. Based on the SSR sequence, 6683 primer pairs were successfully designed, 40 primer pairs were randomly selected, and 10 primer pairs were identified as polymorphic loci from 186 individuals of A. emeiensis. The EST-SSR markers examined in this study will be informative for future population genetic studies of A. emeiensis.

Microsatellite Variation in the Most Devastating Beetle Pests (Coleoptera: Curculionidae) of Agricultural and Forest Crops

Weevils, classified in the family Curculionidae (true weevils), constitute a group of phytophagous insects of which many species are considered significant pests of crops. Within this family, the red palm weevil (RPW), Rhynchophorus ferrugineus, has an integral role in destroying crops and has invaded all countries of the Middle East and many in North Africa, Southern Europe, Southeast Asia, Oceania, and the Caribbean Islands. Simple sequence repeats (SSRs), also termed microsatellites, have become the DNA marker technology most applied to study population structure, evolution, and genetic diversity. Although these markers have been widely examined in many mammalian and plant species, and draft genome assemblies are available for many species of true weevils, very little is yet known about SSRs in weevil genomes. Here we carried out a comparative analysis examining and comparing the relative abundance, relative density, and GC content of SSRs in previously sequenced draft genomes of nine true weevils, with an emphasis on R. ferrugineus. We also used Illumina paired-end sequencing to generate draft sequence for adult female RPW and characterized it in terms of perfect SSRs with 1–6 bp nucleotide motifs. Among weevil genomes, mono- to trinucleotide SSRs were the most frequent, and mono-, di-, and hexanucleotide SSRs exhibited the highest GC content. In these draft genomes, SSR number and genome size were significantly correlated. This work will aid our understanding of the genome architecture and evolution of Curculionidae weevils and facilitate exploring SSR molecular marker development in these species.

Expression Profiles of Microsatellites in Fruit Tissues of Akebia trifoliata and Development of Efficient EST-SSR Markers

Akebia trifoliata, a member of the family Lardizabalaceae, has high exploitation potential for multiple economic purposes, so genetic improvements to meet requirements for commercial demand are needed. However, this progress is largely impeded by a lack of effective selection markers. In this study, we obtained 271.49 Gb of clean transcriptomic data from 12 samples (three tissues at four developmental stages) of A. trifoliata fruit. We identified 175,604, 194,370, and 207,906 SSRs from the de novo assembled 416,363, 463,756, and 491,680 unigene sequences obtained from the flesh, seed, and rind tissues, respectively. The profile and proportion of SSR motifs expressed in each fruit tissue and developmental stage were remarkably similar, but many trinucleotide repeats had differential expression levels among different tissues or at different developmental stages. In addition, we successfully designed 16,869 functional EST-SSR primers according to the annotated unigenes. Finally, 94 and 72 primer pairs out of 100 randomly selected primer pairs produced clear bands and polymorphic bands, respectively. These results were also used to elucidate the expression profiles of different tissues at various stages. Additionally, we provided a set of effective, polymorphic, and reliable EST-SSR markers sufficient for accelerating the discovery of metabolic and pathway-specific functional genes for genetic improvement and increased commercial productivity.

Transcriptome Analysis and Identification of a Female-Specific SSR Marker in Pistacia chinensis Based on Illumina Paired-End RNA Sequencing

Pistacia chinensis Bunge (P. chinensis), a dioecious plant species, has been widely found in China. The female P. chinensis plants are more important than male plants in agricultural production, as their seeds can serve as an ideal feedstock for biodiesel. However, the sex of P. chinensis plants is hard to distinguish during the seedling stage due to the scarcity of available transcriptomic and genomic information. In this work, Illumina paired-end RNA sequencing assay was conducted to unravel the transcriptomic profiles of female and male P. chinensis flower buds. In total, 50,925,088 and 51,470,578 clean reads were obtained from the female and male cDNA libraries, respectively. After quality checks and de novo assembly, a total of 83,370 unigenes with a mean length of 1.3 kb were screened. Overall, 64,539 unigenes (77.48%) could be matched in at least one of the NR, NT, Swiss-Prot, COG, KEGG, and GO databases, 71 of which were putatively related to the floral development of P. chinensis. Additionally, 21,662 simple sequence repeat (SSR) motifs were identified in 17,028 unigenes of P. chinensis, and the mononucleotide motif was the most dominant type of repeats (52.59%) in P. chinensis, followed by dinucleotide (22.29%), trinucleotide (20.15%). The most abundant repeats were AG/CT (13.97%), followed by AAC/GTT (6.75%) and AT/TA (6.10%). Based on these SSR, 983 EST-SSR primers were designed, 151 of which were randomly chosen for validation. Of these validated EST-SSR markers, 25 SSR markers were found to be polymorphic between male and female plants. One SSR marker, namelyPCSSR55, displayed excellent specificity in female plants, which could clearly distinguish between male and female P. chinensis. Altogether, our findings not only reveal that the EST-SSR marker is extremely effective in distinguishing between male and female P. chinensis but also provide a solid framework for sex determination of plant seedlings.

Rapid Genome-Wide Location-Specific Polymorphic SSR Marker Discovery in Black Pepper by GBS Approach

Black pepper (Piper nigrum), the "King of Spices," is an economically important spice in India and is known for its medicinal and cultural values. SSRs, the tandem repeats of small DNA sequences, are often polymorphic in nature with diverse applications. For population structure, QTL/gene discovery, MAS, and diversity analysis, it is imperative to have their location specificity. The existing PinigSSRdb catalogs ~70K putative SSR markers but these are anonymous (unknown chromosomal location), based on 916 scaffolds rather than 26 chromosomes. Under this study, we generated ddRAD sequence data of 29 black pepper genotypes from all over India, being low-cost and most efficient technique for the identification of polymorphic markers. The major limitation of ddRAD with compromised/non-uniform coverage has been successfully overcome by taking advantage of chromosome-wise data availability. The latest black pepper genome assembly was used to extract genome-wide SSRs. A total of 276,230 genomic SSRs were mined distributed over 26 chromosomes, with relative density of 362.88 SSRs/Mb and average distance of 2.76 Kb between two SSRs. This assembly was also used to find the polymorphic SSRs in the generated GBS data of 29 black pepper genotypes utilizing rapid and cost-effective method giving 3,176 polymorphic SSRs, out of which 2015 were found to be hypervariable. The developed web-genomic resource, BlackP2MSATdb (http://webtom.cabgrid.res.in/blackp2msatdb/), is the largest and first reported web resource for genomic and polymorphic SSRs of black pepper, which is useful to develop varietal signature, coreset, physical map, QTL/gene identification, and MAS in endeavor of black pepper production.

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.

Full-length transcriptome analysis of Adiantum flabellulatum gametophyte

Ferns are important components of plant communities on earth, but their genomes are generally very large, with many redundant genes, making whole genome sequencing of ferns prohibitively expensive and time-consuming. This means there is a significant lack of fern reference genomes, making molecular biology research difficult. The gametophytes of ferns can survive independently, are responsible for sexual reproduction and the feeding of young sporophytes, and play an important role in the alternation of generations. For this study, we selected Adiantum flabellulatum as it has both ornamental and medicinal value and is also an indicator plant of acidic soil. The full-length transcriptome sequencing of its gametophytes was carried out using PacBio three-generation sequencing technology. A total of 354,228 transcripts were obtained, and 231,705 coding sequences (CDSs) were predicted, including 5,749 transcription factors (TFs), 2,214 transcription regulators (TRs) and 4,950 protein kinases (PKs). The transcripts annotated by non-redundant protein sequence database (NR), Kyoto encyclopedia of genes and genomes (KEGG), eukaryotic ortholog groups (KOG), Swissprot, protein family (Pfma), nucleotide sequence database (NT) and gene ontology (GO) were 251,501, 197,474, 193,630, 194,639, 195,956, 113,069 and 197,883, respectively. In addition, 138,995 simple sequence repeats (SSRs) and 111,793 long non-coding RNAs (lncRNAs) were obtained. We selected nine chlorophyll synthase genes for qRT-PCR, and the results showed that the full-length transcript sequences and the annotation information were reliable. This study can provide a reference gene set for subsequent gene expression quantification.

SSR Loci Analysis in Transcriptome and Molecular Marker Development in Polygonatum sibiricum

To study the SSR loci information and develop molecular markers, a total of 435,858 unigenes in transcriptome of Polygonatum sibiricum were used to explore SSR. The distribution frequency of SSR and the basic characteristics of repeat motifs were analyzed using MISA software, and SSR primers were designed by Primer 3.0 software and then validated by PCR. Moreover, the gene function analysis of SSR Unigene was obtained by Blast. The results showed that 112,728 SSR loci were found in the transcriptome of Polygonatum sibiricum, which distributed in 435,858 unigenes with a distribution frequency of 25.86%. Mo-nucleotide and Di-nucleotide repeat were the main types, accounted for 83.83% of all SSRs. The repeat motifs of A/T and AC/GT were the predominant repeat types of Mo-nucleotide and Di-nucleotide, respectively. A total of 113,305 pairs of SSR primers with the potential to produce polymorphism were designed for maker development. One hundred and fifty-four of the 500 randomly selected primers not only produced fragments with expected molecular size but also had high polymorphism, which could accurately separate the tested varieties. The gene function of unigenes containing SSR was mostly related to the molecular function of Polygonatum sibiricum. The SSR markers in transcriptome of Polygonatum sibiricum show rich type, strong specificity, and high potential of polymorphism, which will benefit the candidate gene mining and marker-assisted breeding. The developed markers can also provide technical methods for molecular identification of intraspecific species of Polygonatum Mill. and maker-assisted breeding of superior varieties of Polygonatum Mill.

Development of EST-SSR markers based on transcriptome and its validation in ginger (Zingiber officinale Rosc.)

Ginger (Zingiber officinale Rosc.) is an economically important and valuable spice crop around the world. It is used as food, spice, condiment, and medicine. A considerable extent of genetic diversity in ginger occurs in the Western Ghats and North-Eastern India. However, genetic diversity studies at the molecular level in ginger is limited due to limited availability of genetic and genomic information. In the present study, for the first time, we have identified and validated expressed sequence tag (EST)-simple sequence repeat (SSR) markers from ginger. We obtained 16,790 EST-SSR loci from 78987 unigenes, and 4597 SSR loci in the predicted 76929 coding sequences from RNA-Seq assembled contigs of ginger through Illumina paired-end sequencing. Gene ontology results indicate that the unigenes with SSR loci participate in various biological processes such as metabolism, growth, and development in ginger. One hundred and twenty-five primer pairs were designed from unigenes and coding sequences. These primers were tested for PCR optimization, characterization, and amplification and identified 12 novel EST-SSR markers. Twelve flanking polymorphic EST-SSR primers were validated using 48 ginger genotypes representing North-Eastern India and different eco-geographical adaptations by PCR amplification and allele sizing through capillary electrophoresis. Twelve EST-SSR primers generated a total of 111 alleles with an average of 9.25 alleles per locus and allele sizes ranging between 115-189bp. This study implies that the SSR markers designed from transcriptome sequences provides ample EST-SSR resources, which are helpful for genetic diversity analysis of Zingiberaceae species and molecular verification of ginger genotypes.

Transcriptome sequencing and microsatellite marker discovery in Ailanthus altissima (Mill.) Swingle (Simaroubaceae)

Ailanthus altissima Swingle, is a tree species native to East Asia and has a great potential in decorative, bioenergy and industrial applications in many countries. To date, despite its commercial importance, the genomic and genetic resources available for this species are still insufficient. In this study, we characterized the transcriptome of A. altissima and developed thirteen EST-SSRs (expressed sequence tag-simple sequence repeats) based on Illumina paired-end RNA sequencing (RNA-seq). Besides, we developed ten polymorphic chloroplast microsatellite (cpSSR) markers using the available chloroplast genome of A. altissima. The transcriptome data produced 87,797 unigenes, of which 64,891 (73.91%) unigenes were successfully annotated in at least one protein database. For cpSSR markers the number of detected alleles (N) per marker varied from three at cpSSR12 to twelve at cpSSR8, the unbiased haploid diversity indices (uh) varied from 0.111 to 0.485, and haploid diversity indices (h) ranged from 0.101 to 0.444 with an average unbiased haploid diversity index (uh) of 0.274. Overall, a total of 65 different cpSSR alleles were identified at the ten loci among 165 individuals of A. altissima. The allele number per locus for EST-SSRs varied from 2.143 to 9.357, and the values of observed and expected heterozygosity ranged from 0.312 to 1.000 and 0.505 to 0.826, respectively. The molecular markers developed in this study will facilitate future genetic diversity, population structure, long distance-gene transfer and pollen-based gene flow analyses of A. altissima populations from its known distribution ranges in China focusing on planted and natural forest stands.

Genome survey sequencing of Atractylodes lancea and identification of its SSR markers

Atractylodes lancea (Thunb.) DC. is a traditional Chinese medicine rich in sesquiterpenes that has been widely used in China and Japan for the treatment of viral infections. Despite its important pharmacological value, genomic information regarding A. lancea is currently unavailable. In the present study, the whole genome sequence of A. lancea was obtained using an Illumina sequencing platform. The results revealed an estimated genome size for A. lancea of 4,159.24 Mb, with 2.28% heterozygosity, and a repeat rate of 89.2%, all of which indicate a highly heterozygous genome. Based on the genomic data of A. lancea, 27,582 simple sequence repeat (SSR) markers were identified. The differences in representation among nucleotide repeat types were large, e.g., the mononucleotide repeat type was the most abundant (54.74%) while the pentanucleotide repeats were the least abundant (0.10%), and sequence motifs GA/TC (31.17%) and TTC/GAA (7.23%) were the most abundant among the dinucleotide and trinucleotide repeat motifs, respectively. A total of 93,434 genes matched known genes in common databases including 48,493 genes in the Gene Ontology (GO) database and 34,929 genes in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. This is the first report to sequence and characterize the whole genome of A. lancea and will provide a theoretical basis and reference for further genome-wide deep sequencing and SSR molecular marker development of A. lancea.

Transcriptome wide SSR discovery cross-taxa transferability and development of marker database for studying genetic diversity population structure of Lilium species

Lily belongs to family liliaceae, which mainly propagates vegetatively. Therefore, sufficient number of polymorphic, informative, and functional molecular markers are essential for studying a wide range of genetic parameters in Lilium species. We attempted to develop, characterize and design SSR (simple sequence repeat) markers using online genetic resources for analyzing genetic diversity and population structure of Lilium species. We found di-nucleotide repeat motif were more frequent (4684) within 0.14 gb (giga bases) transcriptome than other repeats, of which was two times higher than tetra-repeat motifs. Frequency of di-(AG/CT), tri-(AGG/CTT), tetra-(AAAT), penta-(AGAGG), and hexa-(AGAGGG) repeats was 34.9%, 7.0%, 0.4%, 0.3%, and 0.2%, respectively. A total of 3607 non-redundant SSR primer pairs was designed based on the sequences of CDS, 5'-UTR and 3'-UTR region covering 34%, 14%, 23%, respectively. Among them, a sub set of primers (245 SSR) was validated using polymerase chain reaction (PCR) amplification, of which 167 primers gave expected PCR amplicon and 101 primers showed polymorphism. Each locus contained 2 to 12 alleles on average 0.82 PIC (polymorphic information content) value. A total of 87 lily accessions was subjected to genetic diversity analysis using polymorphic SSRs and found to separate into seven groups with 0.73 to 0.79 heterozygosity. Our data on large scale SSR based genetic diversity and population structure analysis may help to accelerate the breeding programs of lily through utilizing different genomes, understanding genetics and characterizing germplasm with efficient manner.

Development of EST-SSR markers and association mapping with floral traits in Syringa oblata

BACKGROUND

Lilac (Syringa oblata) is an important woody plant with high ornamental value. However, very limited genetic marker resources are currently available, and little is known about the genetic architecture of important ornamental traits for S. oblata, which is hindering its genetic studies. Therefore, it is of great significance to develop effective molecular markers and understand the genetic architecture of complex floral traits for the genetic research of S. oblata.

RESULTS

In this study, a total of 10,988 SSRs were obtained from 9864 unigene sequences with an average of one SSR per 8.13 kb, of which di-nucleotide repeats were the dominant type (32.86%, 3611). A set of 2042 primer pairs were validated, out of which 932 (45.7%) exhibited successful amplifications, and 248 (12.1%) were polymorphic in eight S. oblata individuals. In addition, 30 polymorphic EST-SSR markers were further used to assess the genetic diversity and the population structure of 192 cultivated S. oblata individuals. Two hundred thirty-four alleles were detected, and the PIC values ranged from 0.23 to 0.88 with an average of 0.51, indicating a high level of genetic diversity within this cultivated population. The analysis of population structure showed two major subgroups in the association population. Finally, 20 significant associations were identified involving 17 markers with nine floral traits using the mixed linear model. Moreover, marker SO104, SO695 and SO790 had significant relationship with more than one trait.

CONCLUSION

The results showed newly developed markers were valuable resource and provided powerful tools for genetic breeding of lilac. Beyond that, our study could serve an efficient foundation for further facilitate genetic improvement of floral traits for lilac.

De novo transcriptome assembly and population genetic analyses of an important coastal shrub, Apocynum venetum L

BACKGROUND

Apocynum venetum L. is an important medicinal plant that is mainly distributed in the coastal areas and northwest of China. In addition to its high medical and economic value, its adaptation to saline-alkali and coastal saline lands makes A. venetum an ideal candidate for use in vegetation restoration. To date, the study of A. venetum has been limited in the northwest region of China, little attention has been paid to the genetic diversity and population structure of A. venetum populations in the coastal region. Here, we performed transcriptome sequencing of total RNA from A. venetum leaves and developed efficient expressed sequence tag-simple sequence repeat (EST-SSR) markers for analyzing the genetic diversity and population structure of A. venetum in the coastal region.

RESULTS

A total of 86,890 unigenes were generated after de novo assembly, and 68,751 of which were successfully annotated by searching against seven protein databases. Furthermore, 14,072 EST-SSR loci were detected and 10,243 primer pairs were successfully designed from these loci. One hundred primer pairs were randomly selected and synthesized, twelve primer pairs were identified as highly polymorphic and further used for population genetic analysis. Population genetic analyses showed that A. venetum exhibited low level of genetic diversity (mean alleles per locus, N_A = 3.3; mean expected heterozygosity, H_E = 0.342) and moderate level of genetic differentiation among the populations (genetic differentiation index, F_ST = 0.032-0.220) in the coastal region. Although the contemporary (mean m_c = 0.056) and historical (mean m_h = 0.106) migration rates among the six A. venetum populations were moderate, a decreasing trend over the last few generations was detected. Bayesian structure analysis clustered six populations into two major groups, and genetic bottlenecks were found to have occurred in two populations (QG, BH).

CONCLUSIONS

Using novel EST-SSR markers, we evaluated the genetic variation of A. venetum in the coastal region and determined conservation priorities based on these findings. The large dataset of unigenes and SSRs identified in our study, combining samples from a broader range, will support further research on the conservation and evolution of this important coastal plant and its related species.

De novo assembly and Transcriptome characterization of an endemic species of Vietnam, Panax vietnamensis Ha et Grushv., including the development of EST-SSR markers for population genetics

BACKGROUND

Understanding the genetic diversity in endangered species that occur inforest remnants is necessary to establish efficient strategies for the species conservation, restoration and management. Panax vietnamensis Ha et Grushv. is medicinally important, endemic and endangered species of Vietnam. However, genetic diversity and structure of population are unknown due to lack of efficient molecular markers.

RESULTS

In this study, we employed Illumina HiSeq™ 4000 sequencing to analyze the transcriptomes of P. vietnamensis (roots, leaves and stems). Raw reads total of 23,741,783 was obtained and then assembled, from which the generated unigenes were 89,271 (average length = 598.3191 nt). The 31,686 unigenes were annotated in different databases i.e. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Nucleotide Collection (NR/NT) and Swiss-Prot for functional annotation. Further, 11,343 EST-SSRs were detected. From 7774 primer pairs, 101 were selected for polymorphism validation, in which; 20 primer pairs were successfully amplified to DNA fragments and significant amounts of polymorphism was observed within population. The nine polymorphic microsatellite loci were used for population structure and diversity analyses. The obtained results revealed high levels of genetic diversity in populations, the average observed and expected heterozygosity were HO = 0.422 and HE = 0.479, respectively. During the Bottleneck analysis using TPM and SMM models (p < 0.01) shows that targeted population is significantly heterozygote deficient. This suggests sign of the bottleneck in all populations. Genetic differentiation between populations was moderate (FST = 0.133) and indicating slightly high level of gene flow (Nm = 1.63). Analysis of molecular variance (AMOVA) showed 63.17% of variation within individuals and 12.45% among populations. Our results shows two genetic clusters related to geographical distances.

CONCLUSION

Our study will assist conservators in future conservation management, breeding, production and habitats restoration of the species.

Development and characterization of bZIP transcription factor based SSRs in wheat

Assessment of existing diversity is the key for germplasm conservation and crop improvement. Wheat (Triticum aestivum L.) is among the most important cereal crop and consumed by two billion world's populations. DNA-based markers are predominantly used for diversity characterization because they are easy to develop and not influenced by environment. Among them microsatellites (simple sequence repeats, SSRs) are most suitable due to their genome-wide distribution, hypervariability and reproducibility for their applications in diversity, genetic improvement, and molecular breeding. bZIP transcription factors play major roles in plants in light and stress signalling, seed development, and defence. A total of 846 SSRs were identified from 370 wheat cDNA sequences and a sub-set of 35 polymorphic TabZIPMS (TriticumaestivumbZIP MicroSatellites) was used for diversity and genetic structure analysis of 92 Indian wheat varieties and related species. 114 SSR variants ranging from 2 to 5 per SSR locus were detected for 35 SSRs in the varieties. Average polymorphic information content (PIC) and observed heterozygosity was found to be 0.135 and 0.838, respectively. Thirty-four SSRs showed cross-transferability into different related species. Combined Bayesian model and Jaccard's similarity based genetic clustering analysis revealed two clusters of 80 bread wheat varieties and one separate cluster of related species. In this study, a total 35 novel bZIP-derived SSRs were identified in a set 370 bZIP genes and shown high polymorphism and cross-species transferability in wheat. The findings provide resources for future utilization in genetic resource conservation, trait introgression, breeding and varietal development.

De novo Transcriptome Assembly, Gene Annotation and SSR Marker Development in the Moon Seed Genus Menispermum (Menispermaceae)

The moonseed genus Menispermum L. (Menispermaceae) is disjunctly distributed in East Asia and eastern North America. Although Menispermum has important medicinal value, genetic and genomic information is scarce, with very few available molecular markers. In the current study, we used Illumina transcriptome sequencing and de novo assembly of the two Menispermum species to obtain in-depth genetic knowledge. From de novo assembly, 53,712 and 78,921 unigenes were generated for M. canadense and M. dauricum, with 37,527 (69.87%) and 55,211 (69.96%) showing significant similarities against the six functional databases, respectively. Moreover, 521 polymorphic EST-SSRs were identified. Of them, 23 polymorphic EST-SSR markers were selected to investigate the population genetic diversity within the genus. The newly developed EST-SSR markers also revealed high transferability among the three examined Menispermaceae species. Overall, we provide the very first transcriptomic analyses of this important medicinal genus. In addition, the novel microsatellite markers developed here will aid future studies on the population genetics and phylogeographic patterns of Menispermum at the intercontinental geographical scale.

Genome-wide characterization of simple sequence repeats in Palmae genomes

Background

Microsatellites or simple sequence repeats (SSRs) have become the most significant DNA marker technology used in genetic research. The availability of complete draft genomes for a number of Palmae species has made it possible to perform genome-wide analysis of SSRs in these species. Palm trees are tropical and subtropical plants with agricultural and economic importance due to the nutritional value of their fruit cultivars.

Objective

This is the first comprehensive study examining and comparing microsatellites in completely-sequenced draft genomes of Palmae species.

Methods

We identified and compared perfect SSRs with 1–6 bp nucleotide motifs to characterize microsatellites in Palmae species using PERF v0.2.5. We analyzed their relative abundance, relative density, and GC content in five palm species: Phoenix dactylifera, Cocos nucifera, Calamus simplicifolius, Elaeis oleifera, and Elaeis guineensis.

Results

A total of 118241, 328189, 450753, 176608, and 70694 SSRs were identified, respectively. The six repeat types were not evenly distributed across the five genomes. Mono- and dinucleotide SSRs were the most abundant, and GC content was highest in tri- and hexanucleotide SSRs.

Conclusion

We envisage that this analysis would further substantiate more in-depth computational, biochemical, and molecular studies on the roles SSRs may play in the genome organization of the palm species. The current study contributes a detailed characterization of simple sequence repeats in palm genomes.

Electronic supplementary material

The online version of this article (10.1007/s13258-020-00924-w) contains supplementary material, which is available to authorized users.

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.

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 of genome-wide polymorphic microsatellite markers for Trichinella spiralis

Background

Trichinella nematodes are globally distributed food-borne pathogens, in which Trichinella spiralis is the most common species in China. Microsatellites are a powerful tool in population genetics and phylogeographic analysis. However, only a few microsatellite markers were reported in T. spiralis. Thus, there is a need to develop and validate genome-wide microsatellite markers for T. spiralis.

Methods

Microsatellites were selected from shotgun genomic sequences using MIcroSAtellite identification tool (MISA). The identified markers were validated in 12 isolates of T. spiralis in China.

Results

A total of 93,140 microsatellites were identified by MISA from 9267 contigs in T. spiralis genome sequences, in which 16 polymorphic loci were selected for validation by PCR with single larvae from 12 isolates of T. spiralis in China. There were 7–19 alleles per locus (average 11.25 alleles per locus). The observed heterozygosity (H_O) and expected heterozygosity (H_E) ranged from 0.325 to 0.750 and 0.737 to 0.918, respectively. The polymorphism information content (PIC) ranged from 0.719 to 0.978 (average 0.826). Among the 16 loci, markers for 10 loci could be amplified from all 12 international standard strains of Trichinella spp.

Conclusions

Sixteen highly polymorphic markers were selected and validated for T. spiralis. Primary phylogenetic analysis showed that these markers might serve as a useful tool for genetic studies of Trichinella parasites.

Genome-wide characterization and analysis of microsatellite sequences in camelid species

Microsatellites or simple sequence repeats (SSRs) are among the genetic markers most widely utilized in research. This includes applications in numerous fields such as genetic conservation, paternity testing, and molecular breeding. Though ordered draft genome assemblies of camels have been announced, including for the Arabian camel, systemic analysis of camel SSRs is still limited. The identification and development of informative and robust molecular SSR markers are essential for marker assisted breeding programs and paternity testing. Here we searched and compared perfect SSRs with 1–6 bp nucleotide motifs to characterize microsatellites for draft genome sequences of the Camelidae. We analyzed and compared the occurrence, relative abundance, relative density, and guanine-cytosine (GC) content in four taxonomically different camelid species: Camelus dromedarius, C. bactrianus, C. ferus, and Vicugna pacos. A total of 546762, 544494, 547974, and 437815 SSRs were mined, respectively. Mononucleotide SSRs were the most frequent in the four genomes, followed in descending order by di-, tetra-, tri-, penta-, and hexanucleotide SSRs. GC content was highest in dinucleotide SSRs and lowest in mononucleotide SSRs. Our results provide further evidence that SSRs are more abundant in noncoding regions than in coding regions. Similar distributions of microsatellites were found in all four species, which indicates that the pattern of microsatellites is conserved in family Camelidae.

Development and Application of EST-SSR Markers for DNA Fingerprinting and Genetic Diversity Analysis of the Main Cultivars of Black Locust (Robinia pseudoacacia L.) in China

Black locust (Robinia pseudoacacia L.) is an economically and ecologically important tree species which is used for pillar construction, honey production and soil improvement. More EST-SSR (Expressed sequence tag simple sequence repeat) markers of black locust can be used as a complement and improvement of Genomic-SSR markers for the identification of the function of gene and the construction of genetic map. Additionally, currently there is no simple method for identifying black locust cultivars. In this study, we obtained 2702 unigenes from 3095 expressed sequence tags (ESTs) from the National Center of Biotechnology Information (NCBI) database to identify simple sequence repeats (SSRs) in R. pseudoacacia samples. A total of 170 SSR loci were found to be distributed in 162 non-redundant sequences with a frequency of 6.29%. Dinucleotide repeats were the most predominant types among microsatellites (62.35%), followed by tri-nucleotide repeats (25.88%); the remaining SSRs accounted for less than 12%. The repeat motifs AG/TC (29.25%) and CT/GA (29.25%) were the most abundant among dinucleotides, and AAT/TTA (15.91%) was the most common among tri-nucleotides. A total of 62 primer pairs were designed to screen polymorphic and stable SSR loci. The resulting 25 EST-SSR markers capable of amplifying polymorphic, stable, and repeatable products. Eight newly developed EST-SSR markers and four published SSR markers were selected for DNA fingerprinting and genetic diversity analysis of the 123 main R. pseudoacacia cultivars in China. The 12 SSR loci amplified 102 alleles, with an average number of alleles per locus of 8.5 (range 4–15). The average polymorphism information content at the 12 SSR loci for the 123 cultivars was 0.670 (range 0.427–0.881). The 123 cultivars clustered into six main groups based on similarity coefficients, with most cultivars in one subgroup. Fingerprinting was performed using eight SSR markers; 110 black locust cultivars were distinguished. The results of this study increase the availability of EST-SSR markers in black locust and make it a simple method for checking the collection, the certification, and the correct attribution of clones and cultivars.

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.

De novo assembly of transcriptomes, mining, and development of novel EST-SSR markers in Curcuma alismatifolia (Zingiberaceae family) through Illumina sequencing

Curcuma alismatifolia widely used as an ornamental plant in Thailand and Cambodia. This species of herbaceous perennial from the Zingiberaceae family, includes cultivars with a wide range of colours and long postharvest life, and is used as an ornamental cut flower, as a potted plant, and in exterior landscapes. For further genetic improvement, however, little genomic information and no specific molecular markers are available. The present study used Illumina sequencing and de novo transcriptome assembly of two C. alismatifolia cvs, 'Chiang Mai Pink' and 'UB Snow 701', to develop simple sequence repeat markers for genetic diversity studies. After de novo assembly, 62,105 unigenes were generated and 48,813 (78.60%) showed significant similarities versus six functional protein databases. In addition, 9,351 expressed sequence tag-simple sequence repeats (EST-SSRs) were identified with a distribution frequency of 12.5% total unigenes. Out of 8,955 designed EST-SSR primers, 150 primers were selected for the development of potential molecular markers. Among these markers, 17 EST-SSR markers presented a moderate level of genetic diversity among three C. alismatifolia cultivars, one hybrid, three Curcuma, and two Zingiber species. Three different genetic groups within these species were revealed using EST-SSR markers, indicating that the markers developed in this study can be effectively applied to the population genetic analysis of Curcuma and Zingiber species. This report describes the first analysis of transcriptome data of an important ornamental ginger cultivars, also provides a valuable resource for gene discovery and marker development in the genus Curcuma.

Identification of genes involved in fruit development/ripening in Capsicum and development of functional markers

The molecular mechanism of the underlying genes involved in the process of fruit ripening in Capsicum (family Solanaceae) is not clearly known. In the present study, we identified orthologs of 32 fruit development/ripening genes of tomato in Capsicum, and validated their expression in fruit development stages in C. annuum, C. frutescens, and C. chinense. In silico expression analysis using transcriptome data identified a total of 12 out of 32 genes showing differential expression during different stages of fruit development in Capsicum. Real time expression identified gene LOC107847473 (ortholog of MADS-RIN) had substantially higher expression (>500 folds) in breaker and mature fruits, which suggested the non-climacteric ripening behaviour of Capsicum. However, differential expression of Ehtylene receptor 2-like (LOC107873245) gene during fruit maturity supported the climacteric behaviour of only C. frutescens (hot pepper). Furthermore, development of 49 gene based simple sequence repeat (SSR) markers would help in selection of identified genes in Capsicum breeding.

Generation and classification of transcriptomes in two Croomia species and molecular evolution of CYC/TB1 genes in Stemonaceae

The genus Croomia (Stemonaceae) is an excellent model for studying the evolution of the Eastern Asia (EA)-Eastern North America (ENA) floristic disjunction and the genetic mechanisms of floral zygomorphy formation. In addition to the presence of both actinomorphic and zygomorphic flowers within the genus, species are disjunctively distributed between EA and ENA. However, due to the limited availability of genomic resources, few studies of Croomia have examined these questions. In this study, we sequenced the floral and leaf transcriptomes of the zygomorphic flowered C roomia heterosepala and the actinomorphic flowered Croomia japonica, and used comparative genomic approaches to investigate the transcriptome evolution of the two closely related species. The sequencing and de novo assembly of transcriptomes from flowers of C. heterosepala (ChFlower), flowers of C. japonica (CjFlower), and leaves of C. japonica (CjLeaf) yielded 57,193, 62,131 and 64,448 unigenes, respectively. In addition, estimation of Ka/Ks ratios for 11,566 potential orthologous groups between ChFlower and CjFlower revealed that only six pairs had Ka/Ks ratios significantly greater than 1 and are likely under positive selection. A total of 429 single copy nuclear genes (SCNGs) and 21,460 expression sequence tags-simple sequence repeats (EST-SSRs) were identified in this study. Specifically, we identified seven CYC/TB1-like genes from Stemonaceae. Phylogenetic and molecular evolution analyses indicated that these CYC/TB1-like genes formed a monophyletic clade (SteTBL1) and were subject to strong purifying selection. The shifts of floral symmetry in Stemonaceae do not appear to be correlated with TBL copy number.

Transcriptome Analysis of Cinnamomum chago: A Revelation of Candidate Genes for Abiotic Stress Response and Terpenoid and Fatty Acid Biosyntheses

Cinnamomum chago, an endangered species endemic to Yunnan province, possesses large economic and phylogenetic values in Lauraceae. However, the genomic information of this species remains relatively unexplored. In this study, we used RNAseq technology to characterize and annotate the C. chago transcriptome and identify candidate genes involved in special metabolic pathways and gene-associated simple sequence repeats (SSRs) and single-nucleotide polymorphism (SNP). A total of 129,097 unigenes, with a mean length of 667 bp and an N50 length of 1,062 bp, were assembled. Among these genes, 56,887 (44.07%) unigenes were successfully annotated using at least one database. Furthermore, 47 and 46 candidate genes were identified in terpenoid biosynthesis and fatty acid biosynthesis, respectively. A total of 22 candidate genes participated in at least one abiotic stress response of C. chago. Additionally, a total of 25,654 SSRs and 640 SNPs were also identified. Based on these potential loci, 55 novel expressed sequence tag (EST)-SSR primers were successfully developed. This work provides comprehensive transcriptomic data that can be used to establish a valuable information platform for gene prediction, signaling pathway investigation, and molecular marker development for C. chago and other related species. Such a platform can facilitate further studies on germplasm conservation and utilization of Lauraceae species.

De Novo assembly, characterization and development of EST-SSRs from Bletilla striata transcriptomes profiled throughout the whole growing period

Bletilla striata is an endangered orchid that has been used for millennia as a medicinal herb, in cosmetics and as a horticultural plant. To construct the first nucleotide database for this species and to develop abundant EST-SSR markers for facilitating further studies, various tissues and organs of plants in the main developmental stages were harvested for mRNA isolation and subsequent RNA sequencing. A total of 106,054,784 clean reads were generated by using Illumina paired-end sequencing technology. The reads were assembled into 127,261 unigenes by the Trinity package; the unigenes had an average length of 612 bp and an N50 of 957 bp. Of these unigenes, 67,494 (51.86%) were annotated in a series of databases. Of these annotated unigenes, 41,818 and 24,615 were assigned to gene ontology categories and clusters of orthologous groups, respectively. Additionally, 20,764 (15.96%) unigenes were mapped onto 275 pathways using the KEGG database. In addition, 25,935 high-quality EST-SSR primer pairs were developed from the 15,433 unigenes by MISA mining. To validate the accuracy of the newly designed markers, 87 of 100 randomly selected primers were effectively amplified; 63 of those yielded PCR products of the expected size, and 25 yielded products with significant amounts of polymorphism among the 4 landraces. Furthermore, the transferability test of the 25 polymorphic markers was performed in 6 individuals of two closely related genus Phalaenopsis and dendrobium. Which results showed a total of 5 markers can successfully amplified among these populations. This research provides a comprehensive nucleotide database and lays a solid foundation for functional gene mining and genomic research in B. striata. The developed EST-SSR primers could facilitate phylogenetic studies and breeding.

GinMicrosatDb: a genome-wide microsatellite markers database for sesame (Sesamum indicum L.)

Molecular breeding in sesame is still at infancy due to limited number of microsatellite markers available and the low level of polymorphism exhibited by them. Therefore, whole genome sequencing was used for development of microsatellite markers so as to ensure availability of substantial number of polymorphic markers for use in marker assisted breeding programs. Whole genome sequencing of sesame variety 'Swetha' was done using Illumina paired-end sequencing and Roche 454 shotgun sequencing technologies (GCA_000975565.1 in GenBank). 'GinMicrosatDb', a genome-wide microsatellite marker database has been developed using the whole genome sequence data of sesame variety 'Swetha'. The database consists of microsatellites localized on both linkage groups and scaffolds with their genomic co-ordinates. It provides five sets of forward and reverse primers for each of the microsatellite loci along with the flanking sequences, primer GC content, product size and melting temperature etc. The distribution of microsatellites can be viewed and selected through a genome browser as well as through a physical map. The newly identified microsatellite markers are expected to help sesame breeders in developing marker tags for traits of economic importance thereby bringing about greater efficiency in marker-assisted selection programs.

De Novo Transcriptome Assembly and Population Genetic Analyses for an Endangered Chinese Endemic Acer miaotaiense (Aceraceae)

Acer miaotaiense (P. C. Tsoong) is a rare and highly endangered plant in China. Because of the lack of genomic information and the limited number of available molecular markers, there are insufficient tools to determine the genetic diversity of this species. Here, 93,305 unigenes were obtained by multiple assembled contigs with a transcriptome sequencing program. Furthermore, 12,819 expressed sequence tag-derived simple sequence repeat (EST-SSR) markers were generated, 300 were randomly selected and synthesized, 19 primer pairs were identified as highly polymorphic (average number of alleles (N_a) = 8, expected heterozygosity (H_e) = 0.635, polymorphism information content (PIC) = 0.604) and were further used for population genetic analysis. All 261 samples were grouped into two genetic clusters by UPGMA, a principal component analyses and a STRUCTURE analyses. A moderate level of genetic differentiation (genetic differentiation index (F_st) = 0.059–0.116, gene flow = 1.904–3.993) among the populations and the major genetic variance (81.01%) within populations were revealed by the AMOVA. Based on the results, scientific conservation strategies should be established using in situ and ex situ conservation strategies. The study provides useful genetic information for the protection of precious wild resources and for further research on the origin and evolution of this endangered plant and its related species.

Characterization of leaf transcriptome, development and utilization of unigenes-derived microsatellite markers in sugarcane (Saccharum sp. hybrid)

Sugarcane (Saccharum species hybrid) is the major source of sugar (> 80% sugar) in the world and is cultivated in more than 115 countries. It has recently gained attention as a source of biofuel (ethanol). Due to genomic complexity, the development of new genomic resources is imperative in understanding the gene regulation and function, and to fine tune the genetic improvement of sugarcane. In this study, a cDNA library was constructed from mature leaves so as to develop ESTs resources which were further compared with nucleotide and protein databases to explore the functional identity of sugarcane genes. The non-redundant ESTs (unigenes) were categorized into 18 metabolic functions. The major categories were bioenergetics and photosynthesis (4%), cell metabolism (5%), development related protein (3%), membrane-related, mobile genetic elements (5%), signal transduction (2%), DNA (1%), RNA (1%) and protein (2%) metabolism, other metabolic processes (3%), transcription factors (1%), transport (4%) and proteins related to stress/defense (4%). From 540 unique ESTs, 212 simple sequence repeats were identified, of which 206 were from 463 singlets and six were mined from 77 contig sequences. A total of 540 unique EST sequences were used for SSR search of which 97 (17.9%) contained specified SSR motifs, generating 212 unique SSRs. The genes characterized in this study and the EST-derived microsatellite markers identified from the cDNA library will enrich genomic resources for association- and linkage-mapping studies in sugarcane.

Genetic Evaluation of Natural Populations of the Endangered Conifer Thuja koraiensis Using Microsatellite Markers by Restriction-Associated DNA Sequencing

Thuja koraiensis Nakai is an endangered conifer of high economic and ecological value in Jilin Province, China. However, studies on its population structure and conservation genetics have been limited by the lack of genomic data. Here, 37,761 microsatellites (simple sequence repeat, SSR) were detected based on 875,792 de novo-assembled contigs using a restriction-associated DNA (RAD) approach. Among these SSRs, 300 were randomly selected to test for polymorphisms and 96 obtained loci were able to amplify a fragment of expected size. Twelve polymorphic SSR markers were developed to analyze the genetic diversity and population structure of three natural populations. High genetic diversity (mean NA = 5.481, HE = 0.548) and moderate population differentiation (pairwise Fst = 0.048–0.078, Nm = 2.940–4.958) were found in this species. Molecular variance analysis suggested that most of the variation (83%) existed within populations. Combining the results of STRUCTURE, principal coordinate, and neighbor-joining analysis, the 232 individuals were divided into three genetic clusters that generally correlated with their geographical distributions. Finally, appropriate conservation strategies were proposed to protect this species. This study provides genetic information for the natural resource conservation and utilization of T. koraiensis and will facilitate further studies of the evolution and phylogeography of the species.

Mining and comparative survey of EST-SSR markers among members of Euphorbiaceae family

Euphorbiaceae represents flowering plants family of tropical and sub-tropical region rich in secondary metabolites of economic importance. To understand and assess the genetic makeup among the members, this study was undertaken to characterize and compare SSR markers from publicly available ESTs and GSSs of nine selected species of the family. Mining of SSRs was performed by MISA, primer designing by Primer3, while functional annotation, gene ontology (GO) and enrichment analysis were performed by Blast2GO. A total 12,878 number of SSRs were detected from 101,701 number of EST sequences. SSR density ranged from 1 SSR/3.22 kb to 1 SSR/15.65 kb. A total of 1873 primer pairs were designed for the annotated SSR-Contigs. About 77.07% SSR-ESTs could be assigned a significant match to the protein database. 3037 unique SSR-FDM were assigned and IPR003657 (WRKY Domain) was found to be the most dominant FDM among the members. 1810 unique GO terms obtained were further subjected to enrichment analysis to obtain 513 statistically significant GO terms mapped to the SSR containing ESTs. Most frequent enriched GO terms were, GO:0003824 for molecular function, GO:0006350 for biological process and GO:0005886 for cellular component, justifying the richness of defensive secondary metabolites and phytomedicine within the family. The results from this study provides tangible insight to genetic make-up and distribution of SSRs. Functional annotation corresponded many genes of unknown functions which may be considered as novel genes or genes responsible for stress specific secondary metabolites. Further studies are required to understand stress specific genes accountable for leveraging the synthesis of secondary metabolites.

Transcriptome sequence analysis and mining of SSRs in Jhar Ber (Ziziphus nummularia (Burm.f.) Wight & Arn) under drought stress

Ziziphus nummularia (Burm.f.) Wight & Arn., a perennial shrub that thrives in the arid regions, is naturally tolerant to drought. However, there are limited studies on the genomics of drought tolerance in Ziziphus sp. In this study, RNA-sequencing of one month old seedlings treated with PEG 6000 was performed using Roche GS-FLX454 Titanium pyrosequencing. A total of 367,176 raw sequence reads were generated, and upon adapter trimming and quality filtration 351,872 reads were assembled de novo into 32,739 unigenes. Further characterization of the unigenes indicated that 73.25% had significant hits in the protein database. Kyoto encyclopedia of genes and genomes database (KEGG) identified 113 metabolic pathways from the obtained unigenes. A large number of drought-responsive genes were obtained and among them differential gene expression of 16 highly induced genes was validated by qRT-PCR analysis. To develop genic-markers, 3,425 simple sequence repeats (SSRs) were identified in 2,813 unigene sequences. The data generated shall serve as an important reservoir for the identification and characterization of drought stress responsive genes for development of drought tolerant crops.

Transcriptome analysis and development of simple sequence repeat (SSR) markers in Zingiber striolatum Diels

Illumina-based paired-end sequencing technology was used for the high-throughput transcriptome sequencing of combined Zingiber striolatum Diels tissues (i.e., root, stem, leaf, flower, and fruit tissues). More than 130 million sequencing reads were generated, and a de novo assembly yielded 287,959 contigs and 112,107 unigenes with an average length of 1029 and 28,891 bp, respectively. Similarity searches with known sequences led to the identification of 51,804 (46.21%) genes. Of the annotated unigenes, 6867 and 51,987 were assigned to Gene Ontology and Clusters of Orthologous Groups categories, respectively. Additionally, 8384 simple sequence repeats (SSRs) were identified as potential molecular markers in the unigenes. Thirty pairs of polymerase chain reaction primers were designed and used to validate the unigenes and assess the associated genomic polymorphism. The PCR amplification products for 25 primer pairs were of the expected size. These primers may represent usable molecular markers. The thousands of SSR markers identified in the present study may be useful for analyses of genetic diversity, genetic linkage mapping, and the identification and improvement of varieties during the breeding of Z. striolatum Diels. The unigene sequences and SSR markers described herein may serve as valuable resources for future investigations of Z. striolatum Diels.

Identification of Putative Precursor Genes for the Biosynthesis of Cannabinoid-Like Compound in Radula marginata

The liverwort Radula marginata belongs to the bryophyte division of land plants and is a prospective alternate source of cannabinoid-like compounds. However, mechanistic insights into the molecular pathways directing the synthesis of these cannabinoid-like compounds have been hindered due to the lack of genetic information. This prompted us to do deep sequencing, de novo assembly and annotation of R. marginata transcriptome, which resulted in the identification and validation of the genes for cannabinoid biosynthetic pathway. In total, we have identified 11,421 putative genes encoding 1,554 enzymes from 145 biosynthetic pathways. Interestingly, we have identified all the upstream genes of the central precursor of cannabinoid biosynthesis, cannabigerolic acid (CBGA), including its two first intermediates, stilbene acid (SA) and geranyl diphosphate (GPP). Expression of all these genes was validated using quantitative real-time PCR. We have characterized the protein structure of stilbene synthase (STS), which is considered as a homolog of olivetolic acid in R. marginata. Moreover, the metabolomics approach enabled us to identify CBGA-analogous compounds using electrospray ionization mass spectrometry (ESI-MS/MS) and gas chromatography mass spectrometry (GC-MS). Transcriptomic analysis revealed 1085 transcription factors (TF) from 39 families. Comparative analysis showed that six TF families have been uniquely predicted in R. marginata. In addition, the bioinformatics analysis predicted a large number of simple sequence repeats (SSRs) and non-coding RNAs (ncRNAs). Our results collectively provide mechanistic insights into the putative precursor genes for the biosynthesis of cannabinoid-like compounds and a novel transcriptomic resource for R. marginata. The large-scale transcriptomic resource generated in this study would further serve as a reference transcriptome to explore the Radulaceae family.

Genome-wide mining of perfect microsatellites and tetranucleotide orthologous microsatellites estimates in six primate species

Advancement in genome sequencing and in silico mining tools have provided new opportunities for comparative primate genomics of microsatellites. The SSRs (simple sequence repeats) numbers were not correlated with the genome size (Pearson, r=0.310, p=0.550), and were positively correlated with the total length of SSRs (Pearson, r=0.992, p=0.00). A total of 224,289 tetranucleotide orthologous microsatellites families and 367 single-copy orthologous SSRs loci were found in six primate species by homologous alignment. The inner mutation types of single-copy orthologous SSRs loci included the copy number variance, point mutation, and chromosomal translocation. The accumulated repeat times and average length of tetranucleotide orthologous microsatellites in Rhinopithecus roxellana, Papio anubis and Macaca mulatta were longer than Homo sapiens and Pan troglodytes, which showed the tetranucleotide orthologous SSRs loci had more repeat times and longer average length on the branches with earlier divergence time, one exception may be Microcebus murinus as a primitive monkey with a smallest morphology in Malagasy. Our conclusion indicated that single-copy tetranucleotide orthologous SSRs sequences accumulated individual mutation more slowly through time in H. sapiens and P. troglodytes than in R. roxellanae, P. anubis and M. mulatta. However, such divergence wouldn't arise uniformly in all branches of the primate tree. A comparison of genomic sequence assemblages would offer remarkable insights about comparisons and contrasts, and the evolutionary processes of the microsatellites involved in human and nonhuman primate species.

Mining and characterization of EST-SSR markers for Zingiber officinale Roscoe with transferability to other species of Zingiberaceae

Zingiber officinale is a model spice herb, well known for its medicinal value. It is primarily a vegetatively propagated commercial crop. However, considerable diversity in its morphology, fiber content and chemoprofiles has been reported. The present study explores the utility of EST-derived markers in studying genetic diversity in different accessions of Z. officinale and their cross transferability within the Zingiberaceae family. A total of 38,115 ESTs sequences were assembled to generate 7850 contigs and 10,762 singletons. SSRs were searched in the unigenes and 515 SSR-containing ESTs were identified with a frequency of 1 SSR per 25.21 kb of the genome. These ESTs were also annotated using BLAST2GO. Primers were designed for 349 EST-SSRs and 25 primer pairs were randomly picked for EST SSR study. Out of these, 16 primer pairs could be optimized for amplification in different accessions of Z. officinale as well as other species belonging to Zingiberaceae. GES454, GES466, GES480 and GES486 markers were found to exhibit 100% cross-transferability among different members of Zingiberaceae.

Genome-wide mining and comparative analysis of microsatellites in three macaque species

Microsatellites are found in taxonomically different organisms, and such repeats are related with genomic structure, function and certain diseases. To characterize microsatellites for macaques, we searched and compared SSRs with 1-6 bp nucleotide motifs in rhesus, cynomolgus and pigtailed macaque. A total of 1395671, 1284929 and 1266348 perfect SSRs were mined, respectively. The most frequent perfect SSRs were mononucleotide SSRs. The most GC-content was in dinucleotide SSRs and the least was in the mononucleotide SSRs. Chromosome size was positively correlated with SSR number and negatively correlated with the relative frequency and density of SSRs. The GC content of chromosome SSRs were negatively correlated with relative frequency of SSRs and GC content of chromosome sequences. The features of microsatellite distribution in assembled genomes of the three species were greatly similar, which revealed that the distributional pattern of microsatellites is probably conservative in genus Macaca. The degenerated number of repeat motifs was found to be different in pentanucleotide and hexanucleotide repeats. Species-specific motifs for each macaque were significantly underrepresented. Overall, SSR frequencies of each chromosome in rhesus macaque were higher than in cynomolgus macaque. The maximum repeat times of mono- to pentanucleotide repeats in cynomolgus macaque was more than other two macaques. These results emphasize the genetic diversity and phylogenetic relationship of genus Macaca species. Our data will be beneficial for comparative genome mapping, understanding the distribution of SSRs and genome structure between these animal models, and provide a foundation for further development and identification of more macaque-specific SSRs.

Development of Novel SSR Markers for Flax (Linum usitatissimum L.) Using Reduced-Representation Genome Sequencing

Flax (Linum usitatissimum L.) is a major fiber and oil yielding crop grown in northeastern China. Identification of flax molecular markers is a key step toward improving flax yield and quality via marker-assisted breeding. Simple sequence repeat (SSR) markers, which are based on genomic structural variation, are considered the most valuable type of genetic marker for this purpose. In this study, we screened 1574 microsatellites from Linum usitatissimum L. obtained using reduced representation genome sequencing (RRGS) to systematically identify SSR markers. The resulting set of microsatellites consisted mainly of trinucleotide (56.10%) and dinucleotide (35.23%) repeats, with each motif consisting of 5-8 repeats. We then evaluated marker sensitivity and specificity based on samples of 48 flax isolates obtained from northeastern China. Using the new SSR panel, the results demonstrated that fiber flax and oilseed flax varieties clustered into two well separated groups. The novel SSR markers developed in this study show potential value for selection of varieties for use in flax breeding programs.

Enriching Genomic Resources and Marker Development from Transcript Sequences of Jatropha curcas for Microgravity Studies

Jatropha (Jatropha curcas L.) is an economically important species with a great potential for biodiesel production. To enrich the jatropha genomic databases and resources for microgravity studies, we sequenced and annotated the transcriptome of jatropha and developed SSR and SNP markers from the transcriptome sequences. In total 1,714,433 raw reads with an average length of 441.2 nucleotides were generated. De novo assembling and clustering resulted in 115,611 uniquely assembled sequences (UASs) including 21,418 full-length cDNAs and 23,264 new jatropha transcript sequences. The whole set of UASs were fully annotated, out of which 59,903 (51.81%) were assigned with gene ontology (GO) term, 12,584 (10.88%) had orthologs in Eukaryotic Orthologous Groups (KOG), and 8,822 (7.63%) were mapped to 317 pathways in six different categories in Kyoto Encyclopedia of Genes and Genome (KEGG) database, and it contained 3,588 putative transcription factors. From the UASs, 9,798 SSRs were discovered with AG/CT as the most frequent (45.8%) SSR motif type. Further 38,693 SNPs were detected and 7,584 remained after filtering. This UAS set has enriched the current jatropha genomic databases and provided a large number of genetic markers, which can facilitate jatropha genetic improvement and many other genetic and biological studies.

Development and Application of Transcriptome-Derived Microsatellites in Actinidia eriantha (Actinidiaceae)

Actinidia eriantha Benth. is a diploid perennial woody vine native to China and is recognized as a valuable species for commercial kiwifruit improvement with high levels of ascorbic acid as well as having been used in traditional Chinese medicine. Due to the lack of genomic resources for the species, microsatellite markers for population genetics studies are scarce. In this study, RNASeq was conducted on fruit tissue of A. eriantha, yielding 5,678,129 reads with a total output of 3.41 Gb. De novo assembly yielded 69,783 non-redundant unigenes (41.3 Mb), of which 21,730 were annotated using protein databases. A total of 8,658 EST-SSR loci were identified in 7,495 unigene sequences, for which primer pairs were successfully designed for 3,842 loci (44.4%). Among these, 183 primer pairs were assayed for PCR amplification, yielding 69 with detectable polymorphism in A. eriantha. Additionally, 61 of the 69 polymorphic loci could be successfully amplified in at least one other Actinidia species. Of these, 14 polymorphic loci (mean N_A = 6.07 ± 2.30) were randomly selected for assessing levels of genetic diversity and population structure within A. eriantha. Finally, a neighbor-joining tree and Bayesian clustering analysis showed distinct clustering into two groups (K = 2), agreeing with the geographical distributions of these populations. Overall, our results will facilitate further studies of genetic diversity within A. eriantha and will aid in discriminating outlier loci involved in local adaptation.

De Novo Assembly of Transcriptome and Development of Novel EST-SSR Markers in Rhododendron rex Lévl. through Illumina Sequencing

Transcriptome sequences generated by next-generation sequencing (NGS) technologies can be utilized to rapidly detect and characterize a large number of gene-based microsatellites from different plants. Rhododendron rex Lévl. is a perennial woody species from the family Ericaceae and an endangered plant with high ornamental value endemic to Southwestern China. Nevertheless, the genetic and genomic information of R. rex remain unknown. In this study, we performed transcriptome sequencing for R. rex leaf samples, and generated large transcript sequences for functional characterization and development gene-associated SSR markers. A total of 164,242 unigenes were assembled and 115,089 (70.07%) unigenes were successfully annotated in public databases. In addition, a total of 15,314 potential EST-SSRs were identified, and the frequency of SSRs in the R. rex unigenes was 9.32%, with an average of one EST-SSR per 5.65 kb. The most abundant type was repeated di-nucleotide (54.63%), followed by mono- (26.03%) and tri-nucleotide (18.51%) repeats. Based on the SSR-containing sequence, 100 primer pairs were randomly selected and synthesized and used for assessment of the polymorphism. Thirty-six primer pairs were polymorphic and revealed polymorphism among 20 individuals from four R. rex populations. A total of 197 alleles were identified, with an average of 5.472 alleles per locus. The Polymorphism Information Content ranged from 0.154 to 0.870, with a mean of 0.482. The newly developed EST-SSR markers exhibited high transferability (58.33-83.33%) among the six subgenera. Thus, these novel EST-SSR markers developed would provide valuable sequence resources for population structure, genetic diversity analysis, and genetic resource assessments of R. rex and its related species.