Characterization and development of EST-derived SSR markers in cultivated sweetpotato (Ipomoea batatas)

Genetic fingerprint construction and genetic diversity analysis of sweet potato (Ipomoea batatas) germplasm resources

BACKGROUND

China is the largest producer of sweet potato in the world, accounting for 57.0% of the global output. Germplasm resources are the basis for promoting innovations in the seed industry and ensuring food security. Individual and accurate identification of sweet potato germplasm is an important part of conservation and efficient utilization.

RESULTS

In this study, nine pairs of simple sequence repeat molecular markers and 16 morphological markers were used to construct genetic fingerprints for sweet potato individual identification. Combined with basic information, typical phenotypic photographs, genotype peak graphs, and a two-dimensional code for detection and identification were generated. Finally, a genetic fingerprint database containing 1021 sweet potato germplasm resources in the "National Germplasm Guangzhou Sweet Potato Nursery Genebank in China" was constructed. Genetic diversity analysis of the 1021 sweet potato genotypes using the nine pairs of simple sequence repeat markers revealed a narrow genetic variation range of Chinese native sweet potato germplasm resources, and Chinese germplasm was close to that from Japan and the United States, far from that from the Philippines and Thailand, and the furthest from that from Peru. Sweet potato germplasm resources from Peru had the richest genetic diversity, supporting the view that Peru is the center of origin and domestication of sweet potato varieties.

CONCLUSIONS

Overall, this study provides scientific guidance for the conservation, identification, and utilization of sweet potato germplasm resources and offers a reference to facilitate the discovery of important genes to boost sweet potato breeding.

Identification and functional characterization of a flavonol synthase gene from sweet potato [Ipomoea batatas (L.) Lam.]

Flavonol synthase (FLS) is a key enzyme of the flavonoid biosynthetic pathway, which catalyzes the conversion of dihydroflavonols into flavonols. In this study, the FLS gene IbFLS1 was cloned and characterized from sweet potato. The resulting IbFLS1 protein showed a high similarity with other plant FLSs. The conserved amino acids (HxDxnH motifs) binding ferrous iron and residues (RxS motifs) binding 2-oxoglutarate were found in IbFLS1 at conserved positions, as in other FLSs, suggesting that IbFLS1 belongs to the 2-oxoglutarate-dependent dioxygenases (2-ODD) superfamily. qRT-PCR analysis showed an organ-specific pattern of expression of the IbFLS1 gene, which was predominantly expressed in young leaves. The recombinant IbFLS1 protein could catalyze the conversion of dihydrokaempferol and dihydroquercetin to kaempferol and quercetin, respectively. The results of subcellular localization studies indicated that IbFLS1 was found mainly in the nucleus and cytomembrane. Furthermore, silencing the IbFLS gene in sweet potato changed the color of the leaves to purple, substantially inhibiting the expression of IbFLS1 and upregulating the expression of genes involved in the downstream pathway of anthocyanin biosynthesis (i.e., DFR, ANS, and UFGT). The total anthocyanin content in the leaves of the transgenic plants was dramatically increased, whereas the total flavonol content was significantly reduced. Thus, we conclude that IbFLS1 is involved in the flavonol biosynthetic pathway and is a potential candidate gene of color modification in sweet potato.

Evaluation of Yield and Nutraceutical Traits of Orange-Fleshed Sweet Potato Storage Roots in Two Agro-Climatic Zones of Northern Ethiopia

This study evaluated the genotype by environment interactions in the yield and nutraceutical traits of the orange-fleshed sweet potato (OFSP) storage root in different agro-climatic zones of northern Ethiopia. Five OFSP genotypes were cultivated at three different locations following a randomized complete block design, and the yield, dry matter, beta-carotene, flavonoids, polyphenols, soluble sugars, starch, soluble proteins, and free radical scavenging activity were measured in the storage root. The results showed consistent variations in the nutritional traits of the OFSP storage root depending on both the genotype and the location, as well as on their interaction. Ininda, Gloria, and Amelia were the genotypes that provided the higher yield and dry matter, as well as the higher content of starch and beta-carotene; they also showed a high antioxidant power. These findings suggest that the studied genotypes have the potential to alleviate vitamin A deficiency. This study demonstrated a high possibility of sweet potato production for storage root yield in arid agro-climate regions with limited production inputs. Moreover, the results suggest that it is possible to enhance the yield, dry matter content, beta-carotene, starch, and polyphenols of the OFSP storage root through genotype selection.

First Report on Development of Genome-Wide Microsatellite Markers for Stock (Matthiola incana L.)

Stock (Matthiola incana (L.) R. Br.) is a famous annual ornamental plant with important ornamental and economic value. The lack of DNA molecular markers has limited genetic analysis, genome evolution, and marker-assisted selective breeding studies of M. incana. Therefore, more DNA markers are needed to support the further elucidation of the biology and genetics of M. incana. In this study, a high-quality genome of M. incana was initially assembled and a set of effective SSR primers was developed at the whole-genome level using genome data. A total of 45,612 loci of SSRs were identified; the di-nucleotide motifs were the most abundant (77.35%). In total, 43,540 primer pairs were designed, of which 300 were randomly selected for PCR validation, and as the success rate for amplification. In addition, 22 polymorphic SSR markers were used to analyze the genetic diversity of 40 stock varieties. Clustering analysis showed that all varieties could be divided into two clusters with a genetic distance of 0.68, which were highly consistent with their flower shape (potted or cut type). Moreover, we have verified that these SSR markers are effective and transferable within the Brassicaceae family. In this study, potential SSR molecular markers were successfully developed for 40 M. incana varieties using whole genome analysis, providing an important genetic tool for theoretical and applied research on M. incana.

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.

Metabolomic and Transcriptomic Analyses of the Flavonoid Biosynthetic Pathway for the Accumulation of Anthocyanins and Other Flavonoids in Sweetpotato Root Skin and Leaf Vein Base

Sweetpotato [Ipomoea batatas (L.) Lam.] is a major tuberous root crop that is rich in flavonoids. Here, we discovered a spontaneous mutation in the color of the leaf vein base (LVB) and root skin (RS) in the Zheshu 81 cultivar. The flavonoid and anthocyanin metabolites and molecular mechanism were analyzed using metabolome and transcriptome data. Compared to the wild type, 13 differentially accumulated metabolites (DAMs) in the LVB and 59 DAMs in the RS were all significantly downregulated. Moreover, all the anthocyanin metabolites decreased significantly. The differentially expressed genes (DEGs) encoding the key enzymes in the later enzymatic reaction of anthocyanin and flavonoid were significantly downregulated in the mutant. The expression trends of the transcription factor MYB were evidently related to the anthocyanin content. These results offer insights into the coloration in the LVB and RS and a theoretical basis for determining the regulation of flavonoid and anthocyanin synthesis in sweetpotato.

Assessment of genetic diversity and variety identification based on developed retrotransposon-based insertion polymorphism (RBIP) markers in sweet potato (Ipomoea batatas (L.) Lam.)

Sweet potato, a dicotyledonous and perennial plant, is the third tuber/root crop species behind potato and cassava in terms of production. Long terminal repeat (LTR) retrotransposons are highly abundant in sweet potato, contributing to genetic diversity. Retrotransposon-based insertion polymorphism (RBIP) is a high-throughput marker system to study the genetic diversity of plant species. To date, there have been no transposon marker-based genetic diversity analyses of sweet potato. Here, we reported a structure-based analysis of the sweet potato genome, a total of 21555 LTR retrotransposons, which belonged to the main LTR-retrotransposon subfamilies Ty3-gypsy and Ty1-copia were identified. After searching and selecting using Hidden Markov Models (HMMs), 1616 LTR retrotransposon sequences containing at least two models were screened. A total of 48 RBIP primers were synthesized based on the high copy numbers of conserved LTR sequences. Fifty-six amplicons with an average polymorphism of 91.07% were generated in 105 sweet potato germplasm resources based on RBIP markers. A Unweighted Pair Group Method with Arithmatic Mean (UPGMA) dendrogram, a model-based genetic structure and principal component analysis divided the sweet potato germplasms into 3 groups containing 8, 53, and 44 germplasms. All the three analyses produced significant groupwise consensus. However, almost all the germplasms contained only one primary locus. The analysis of molecular variance (AMOVA) among the groups indicated higher intergroup genetic variation (53%) than intrapopulation genetic variation. In addition, long-term self-retention may cause some germplasm resources to exhibit variable segregation. These results suggest that these sweet potato germplasms are not well evolutionarily diversified, although geographic speciation could have occurred at a limited level. This study highlights the utility of RBIP markers for determining the intraspecies variability of sweet potato and have the potential to be used as core primer pairs for variety identification, genetic diversity assessment and linkage map construction. The results could provide a good theoretical reference and guidance for germplasm research and breeding.

Kibr G. Effect of blending ratio of wheat, orange fleshed sweet potato and haricot bean flour on proximate compositions, β-carotene, physicochemical properties and sensory acceptability of biscuits'

Background: Protein-energy malnutrition and vitamin A deficiency (VAD) are the most important public health issues, and a food-based strategy is crucial to combat those health problems among the vulnerable group of people. Methods: Composite biscuits were made with 100:0:0, 90:5:5, 80:10:10, 70:15:15, 60:20:20, and 50:25:25 percent wheat, haricot bean, and orange-fleshed sweet potato (OFSP) flours.Standard methods were used to evaluate the proximate compositions, β-carotene, physical properties, functional properties, and sensory acceptability. A one-way analysis of variance model was used to statistically evaluate the data using the statistical analysis system software package, version 9.0 standard methods. Results: The results showed that partially replacing wheat with haricot beans and OFSP increased the β-carotene and proximate compositions significantly. When wheat was replaced with haricot beans and OFSP, the physical characteristics of the biscuits did not vary significantly from those of biscuits made entirely of wheat flour. Sensory acceptability (appearance, color, flavor, taste and overall acceptability) was higher in the composite biscuits with up to 40% wheat substitution than in the 100% wheat flour biscuits. Conclusion: Based on the findings of this report, replacing wheat with OFSP and haricot beans in biscuit formulation appears to be promising in improving nutritional quality, sensory acceptability, and beta carotene. It is proposed that these products can mitigate food insecurity and deficiency of vitamin A.

Kibr G. Effect of blending ratio of wheat, orange fleshed sweet potato (Ipomoea batatas L.) powder and haricot bean ( Phaseolus vulgarisL.) flour on proximate composition, physical properties and sensory acceptability of biscuits

Background Protein-energy deficiency (PEM) is still a major health issue in developing countries, and it is the leading cause of disease and death in children under the age of 5 years. Methods 100:0:0; 90:5:5; 80:10:10; 70:15:15; 60:20:20; and 50:25:25 per cent wheat:haricot bean: orange-fleshed sweet potato (OFSP) flours were used to make composite cookies. Standard methods were used to evaluate the proximate structure, physical properties, and sensory assessment. A one-way analysis of variance (ANOVA) model was used to statistically evaluate the data using the statistical analysis system (SAS) software package, version 9.0 standard methods. Results The results showed that partially replacing wheat with haricot bean and OFSP increased the proximate composition significantly. When wheat was replaced with haricot bean and OFSP, the physical characteristics of the biscuits did not vary significantly from those of biscuits made entirely of wheat flour. Sensory acceptability (appearance, color, flavor, taste, and overall acceptability) was higher in the composite biscuits with up to 40% wheat substitution than in the 100% wheat flour biscuits. Conclusions Based on the findings of this report, using OFSP and a haricot bean to wheat flour blend in biscuit formulation appears to be promising in terms of nutritional quality, acceptability, and cost. It is proposed that these products be marketed to vitamin A deficiency (VAD) customers as a newly established product to help mitigate food insecurity.

Development of transcriptome-wide SSR markers for genetic diversity and structure analysis in Macrotyloma uniflorum (Lam.) Verdc

Horsegram is an important drought resistant pulse crop from Fabaceae and can be easily grown in dry lands with no irrigation facilities. However, it remained neglected since long and has been considered as orphan legume which requires immediate attention for its improvement and for the development of new promising varieties in future. In the present study, 7352 simple sequence repeat (SSR) markers were developed from the transcriptome data and 150 SSR were randomly synthesized for validation and diversity analysis in a panel of 58 horsegram genotypes. The synthesized primers included all types of repeats spanning direpeats to hexarepeats. Of the validated SSR markers, 33 markers were polymorphic and produced 40 loci which were used to analyze the genetic diversity and structure of horsegram. In total, 130 alleles were produced in a range of 2-9 alleles with maximum alleles produced by primer HTSSR 155. Expected heterozygosity (He) ranged from 0.03 to 1.00 and observed heterozygosity (Ho) ranged from 0.13 to 0.81. Polymorphism information content value ranged from 0.065 to 0.78. Dendrogram based on UPGMA and principal component analysis showed four groups of the 58 genotypes of horsegram. Structure analysis showed three genetic stocks for the analyzed germplasm. Thus, the developed SSRs can be useful in future population genetics analysis, molecular breeding studies and mapping works in horsegram germplasm as well as in related legume species.

Identification of candidate miRNAs related in storage root development of sweet potato by high throughput sequencing

Sweet potato (Ipomoea batatas L.) is a food consumed worldwide, an industrial raw material and new energy crop. The storage root is the most economical part of the crop. However, the mechanism of storage root initiation and development is still unclear. In this study, conserved and novel miRNAs during storage root development were identified by high-throughput sequencing technology by constructing small RNA libraries from sweet potato fibrous roots (F) and storage roots at four different developmental stages (storage roots with different diameters: 1 cm, D1; 3 cm, D3; 5 cm, D5 and 10 cm, D10). A total of 61 known miRNAs and 471 novel miRNAs were identified. In addition, 145 differentially expressed miRNAs were identified in the F library compared with the four storage root libraries, with 30 known miRNAs and 115 novel miRNAs. Moreover, the targets of the differentially expressed miRNAs were predicted and their network was further investigated by GO analysis using our previous transcriptome data. The GO analysis revealed that antioxidant activity and binding process were the most enriched terms of the target genes. The secondary structure and expression of six candidate miRNAs including three conserved miRNAs and three novel miRNAs were investigated and their predicted targets were validated by qRT-PCR. The results showed that the expression levels of the miRNAs were all consistent with the sequencing data. Most of the miRNAs and their corresponding targets had obvious negative correlations. This study contributed to elucidating the potential miRNA mediated regulatory mechanism of storage root development in sweet potato. The specific differentially expressed miRNAs in sweet potato storage roots can be used to breed high-yield sweet potatoes and other tuberous root crops.

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 genetic diversity detection and population structure analysis in sweetpotato (Ipomoea batatas) using RAD-seq

Sweetpotato (Ipomoea batatas L.) is one of the most important food and grain-forage crops globally. It has been planted in >100 countries. Due to the complexity of the sweetpotato genome, its research is far behind other major food crops. At present, limited information about the sweetpotato genome is available. Thus, it is central to find an efficient approach for the investigation of sweetpotato genome. In this study, RAD-seq (Restriction site-associated DNA sequencing) was used to evaluate sweetpotato genetic structure diversity and to develop relevant SSR markers. The study yielded >128 Gb reliable sequence data from 81 sweetpotato accessions. By analyzing polymorphic tags from each accession, a total of 55,622 restriction-site associated DNA sequencing tags (RAD-seq) were found, containing 907,010 SNP. Genetic analysis divided 81 accessions into five major clusters based on their SNP genotype, which matches the results of genetic analysis and the genetic family tree. In addition, 18,320 SSRs loci were detected and 9336 SSR primer pairs were developed. Eighty-three primer pairs were amplified in different sweetpotato genotypes, 76 of which successfully amplified polymorphism bands. These results provide significant information about sweetpotato genome, which can be used to identify novel gene and to further develop the gene chip. And more significant, clustering results based on the SNP genotype provide an essential reference for breeders to match parent plants in breeding program. Additionally, SSR markers developed in this study will supply a wealth of markers for marker-assisted selection in sweetpotato breeding.

Diversity Analysis of Sweet Potato Genetic Resources Using Morphological and Qualitative Traits and Molecular Markers

The European Union (EU) market for sweet potatoes has increased by 100% over the last five years, and sweet potato cultivation in southern European countries is a new opportunity for the EU to exploit and introduce new genotypes. In view of this demand, the origins of the principal Italian sweet potato clones, compared with a core collection of genotypes from Central and Southern America, were investigated for the first time. This was accomplished by combining a genetic analysis, exploiting 14 hypervariable microsatellite markers, with morphological and chemical measurements based on 16 parameters. From the molecular analyses, Italian accessions were determined to be genetically very similar to the South American germplasm, but they were sub-clustered into two groups. This finding was subsequently confirmed by the morphological and chemical measurements. Moreover, the analysis of the genetic structure of the population suggested that one of the two groups of Italian genotypes may have descended from one of the South American accessions, as predicted on the basis of the shared morphological characteristics and molecular fingerprints. Overall, the combination of two different characterization methods, genetic markers and agronomic traits, was effective in differentiating or clustering the sweet potato genotypes, in agreement with their geographical origin or phenotypic descriptors. This information could be exploited by both breeders and farmers to detect and protect commercial varieties, and hence for traceability purposes.

Elucidation of Galactomannan Biosynthesis Pathway Genes through Transcriptome Sequencing of Seeds Collected at Different Developmental Stages of Commercially Important Indian Varieties of Cluster Bean (Cyamopsis tetragonoloba L.)

Cyamopsis tetragonoloba (L) endosperm predominantly contains guar gum a polysaccharide, which has tremendous industrial applications in food, textile, paper, oil drilling and water treatment. In order to understand the genes controlling galactomannan biosynthesis, mRNA was isolated from seeds collected at different developmental stages; young pods, mature pods and young leaf from two guar varieties, HG365 and HG870 and subjected to Illumina sequencing. De novo assembly of fourteen individual read files from two varieties of guar representing seven developmental stages gave a total of 1,13,607 contigs with an N50 of 1,244 bases. Annotation of assemblies with GO mapping revealed three levels of distribution, namely, Biological Processes, Molecular Functions and Cellular Components. GO studies identified major genes involved in galactomannan biosynthesis: Cellulose synthase D1 (CS D1) and GAUT-like gene families. Among the polysaccharide biosynthetic process (GO:0000271) genes the transcript abundance for CS was found to be predominantly more in leaf samples, whereas, the transcript abundance for GAUT-like steadily increased from 65% to 90% and above from stage1 to stage5 indicating accumulation of galactomannan in developing seeds; and validated by qRT-PCR analysis. Galactomannan quantification by HPLC showed HG365 (12.98–20.66%) and HG870 (7.035–41.2%) gradually increasing from stage1 to stage 5 (10–50 DAA) and highest accumulation occurred in mature and dry seeds with 3.8 to 7.1 fold increase, respectively. This is the first report of transcriptome sequencing and complete profiling of guar seeds at different developmental stages, young pods, mature pods and young leaf material from two commercially important Indian varieties and elucidation of galactomannan biosynthesis pathway. It is envisaged that the data presented herein will be very useful for improvement of guar through biotechnological interventions in future.

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.

RNA-Seq and iTRAQ reveal multiple pathways involved in storage root formation and development in sweet potato (Ipomoea batatas L.)

BACKGROUND

Sweet potato (Ipomoea batatas L.) is the sixth most important food crop in the world. The formation and development of storage roots in sweet potato is a highly complicated and genetically programmed process. However, the underlying mechanisms of storage root development have not yet been elucidated.

RESULTS

To better understand the molecular mechanisms involved in storage root development, a combined analysis of the transcriptome and proteome of sweet potato fibrous roots (F) and storage roots at four different stages (D1, D3, D5 and D10) was performed in the present study. A total of 26,273 differentially expressed genes were identified in a comparison between the fibrous root library and four storage root libraries, while 2558 proteins showed a 1.0-fold or greater expression difference as indicated by isobaric tags for relative and absolute quantitation (iTRAQ) analysis. The combination of the transcriptome and proteome analyses and morphological and physiological data revealed several critical pathways involved in storage root formation and development. First, genes/proteins involved in the development of meristems/cambia and starch biosynthesis were all significantly upregulated in storage roots compared with fibrous roots. Second, multiple phytohormones and the genes related to their biosynthesis showed differential expression between fibrous roots and storage roots. Third, a large number of transcription factors were differentially expressed during storage root initiation and development, which suggests the importance of transcription factor regulation in the development of storage roots. Fourth, inconsistent gene expression was found between the transcriptome and proteome data, which indicated posttranscriptional regulatory activity during the development of storage roots.

CONCLUSION

Overall, these results reveal multiple events associated with storage root development and provide new insights into the molecular mechanisms underlying the regulatory networks involved in storage root development.

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.

In silico prediction and segregation analysis of putative virus defense genes based on SSR markers in sweet potato F1 progenies of cultivars 'New Kawogo' and 'Resisto'

In sweet potato, an anti-virus defense mechanism termed reversion has been postulated to lead to virus freedom from once infected plants. The objectives of this study were to identify anti-virus defense genes and evaluate their segregation in progenies. Reference genes from different plant species were used to assemble transcript sequences of each sweet potato defense gene in silico. Sequences were used for evaluate phylogenetic relationships with similar genes from different plant species, mining respective defense genes and thereafter developing simple sequence repeats (SSRs) for segregation analysis. Eight potential defense genes were identified: RNA dependent RNA polymerases 1, 2, 5, and 6; Argonaute 1, and Dicer-like 1, 2, and 4. Identified genes were differentially related to those of other plants and were observed on different chromosomes. The defense genes contained mono-, di-, tri-, tetra, penta-, and hexa-nucleotide repeat motifs. The SSR markers within progenies were segregated in disomic, co-segregation, nullisomic, monosomic, and trisomic modes. These findings indicate the possibility of deriving and utilizing SSRs using published genomic information. Furthermore, and given that the SSR markers were derived from known genes on defined chromosomes, this work will contribute to future molecular breeding and development of resistance gene analogs in this economically important crop.

Identification and characterization of SSR, SNP and InDel molecular markers from RNA-Seq data of guar (Cyamopsis tetragonoloba, L. Taub.) roots

Background

Guar [Cyamopsis tetragonoloba, L. Taub.] is an important industrial crop because of the commercial applications of the galactomannan gum contained in its seeds. Plant breeding programmes based on marker-assisted selection require a rich resource of molecular markers. As limited numbers of such markers are available for guar, molecular breeding programmes have not been undertaken for the genetic improvement of this important crop. Hence, the present work was done to enrich the molecular markers resource of guar by identifying high quality SSR, SNP and InDel markers from the RNA-Seq data of the roots of two guar varieties.

Results

We carried out RNA-Seq analysis of the roots of two guar varieties, namely, RGC-1066 and M-83. A total of 102,479 unigenes with an average length of 1016 bp were assembled from about 30 million high quality pair-end reads generated by an Illumina HiSeq 2500 platform. The assembled unigenes had 86.55% complete and 97.71% partially conserved eukaryotic genes (CEGs). The functional annotation of assembled unigenes using BLASTX against six databases showed that the guar unigenes were most similar to Glycine max. We could assign GO terms to 45,200 unigenes using the UniProt database. The screening of 102,479 unigenes with MISA and SAMtools version 1.4 softwares resulted in the identification of 25,040 high-confidence molecular markers which consisted of 18,792 SSRs, 5999 SNPs and 249 InDels. These markers tagged most of the genes involved in root development, stress tolerance and other general metabolic activities. Each of the 25,040 molecular markers was characterized, particularly with respect to its position in the unigene. For 71% of the molecular markers, we could determine the names, products and functions of the unigenes. About 80% of the markers, from a random sample of molecular markers, showed PCR amplification.

Conclusions

We have identified and characterized 25,040 high confidence SSR, SNP and InDel molecular markers in guar. It is expected that these markers will be useful in molecular breeding programmes and will also be helpful in studying molecular mechanisms of root development, stress tolerance and gum synthesis in guar.

Electronic supplementary material

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

Comparative analysis of microsatellite, SNP, and InDel markers in four Rhododendron species based on RNA-seq

Rhododendron possesses valuable horticultural and medicinal properties. However, the genetic studies have been hindered due to the lack of genetic markers. Based on RNA-seq, large-scale molecular markers were developed from four Rhododendron species endemic to Dabie Mountains (central China): R. fortunei (5.25 Gb; SSRs, 12,756, one/2.37 kb, 147 types; SNPs, 38,313; InDels, 3,174), R. simsii (5.80 Gb; SSRs, 13,294, one/2.58 kb, 167 types; SNPs, 136,590; InDels, 6,258), R. mariesii (6.53 Gb; SSRs, 15,724, one/2.51 kb, 170 types; SNPs, 44,942; InDels, 4,126), and R. molle (4.35 Gb; SSRs, 10,214, one/2.49 kb, 110 types; SNPs, 77,829; InDels, 3,416). Di-nucleotide repeats were the main type (59.126%-64.314%), and AG/CT repeat (55.18%-61.22%) was the most. In particular, 89 species-specific types had been found. Furthermore, C:G→T:A mutation was the main SNP type (30.475%-34.99%). However, C:G→G:C mutation was the least type in R. fortunei, while T:A→G:C mutation was the least in the other three species. InDels with length of 3 nt was most in R. fortunei, but 1 nt InDels were the main type in the other three species. Twelve microsatellite markers developed from R. simsii reveled high genetic diversity in the four populations, and heterozygote excess was observed. This research would benefit the genetic study, molecular marker-assisted selection, and breeding studies in Rhododendron species.

Genetic diversity and population structure analyses of Plectranthus edulis (Vatke) Agnew collections from diverse agro-ecologies in Ethiopia using newly developed EST-SSRs marker system

BACKGROUND

Plectranthus edulis (Vatke) Agnew (locally known as Ethiopian dinich or Ethiopian potato) is one of the most economically important edible tuber crops indigenous to Ethiopia. Evaluating the extent of genetic diversity within and among populations is one of the first and most important steps in breeding and conservation measures. Hence, this study was aimed at evaluating the genetic diversity and population structure of this crop using collections from diverse agro-ecologies in Ethiopia.

RESULTS

Twenty polymorphic expressed sequence tag based simple sequence repeat (EST-SSRs) markers were developed for P. edulis based on EST sequences of P. barbatus deposited in the GenBank. These markers were used for genetic diversity analyses of 287 individual plants representing 12 populations, and a total of 128 alleles were identified across the entire loci and populations. Different parameters were used to estimate the genetic diversity within populations; and gene diversity index (GD) ranged from 0.31 to 0.39 with overall mean of 0.35. Hierarchical analysis of molecular variance (AMOVA) showed significant but low population differentiation with only 3% of the total variation accounted for variation among populations. Likewise, cluster and STRUCTURE analyses did not group the populations into sharply distinct clusters, which could be attributed to historical and contemporary gene flow and the reproductive biology of the crop.

CONCLUSIONS

These newly developed EST-SSR markers are highly polymorphic within P. edulis and hence are valuable genetic tools that can be used to evaluate the extent of genetic diversity and population structure of not only P. edulis but also various other species within the Lamiaceae family. Among the 12 populations studied, populations collected from Wenbera, Awi and Wolaita showed a higher genetic diversity as compared to other populations, and hence these areas can be considered as hot spots for in-situ conservation as well as for identification of genotypes that can be used in breeding programs.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

Transcriptome profiling and digital gene expression analysis of sweet potato for the identification of putative genes involved in the defense response against Fusarium oxysporum f. sp. batatas

Sweet potato production is constrained by Fusarium wilt, which is caused by Fusarium oxysporum f. sp. batatas (Fob). The identification of genes related to disease resistance and the underlying mechanisms will contribute to improving disease resistance via sweet potato breeding programs. In the present study, we performed de novo transcriptome assembly and digital gene expression (DGE) profiling of sweet potato challenged with Fob using Illumina HiSeq technology. In total, 89,944,188 clean reads were generated from 12 samples and assembled into 101,988 unigenes with an average length of 666 bp; of these unigenes, 62,605 (61.38%) were functionally annotated in the NCBI non-redundant protein database by BLASTX with a cutoff E-value of 10-5. Clusters of Orthologous Groups (COG), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations were examined to explore the unigenes' functions. We constructed four DGE libraries for the sweet potato cultivars JinShan57 (JS57, highly resistant) and XinZhongHua (XZH, highly susceptible), which were challenged with pathogenic Fob. Genes that were differentially expressed in the four libraries were identified by comparing the transcriptomes. Various genes that were differentially expressed during defense, including chitin elicitor receptor kinase 1 (CERK), mitogen-activated protein kinase (MAPK), WRKY, NAC, MYB, and ethylene-responsive transcription factor (ERF), as well as resistance genes, pathogenesis-related genes, and genes involved in salicylic acid (SA) and jasmonic acid (JA) signaling pathways, were identified. These data represent a sequence resource for genetic and genomic studies of sweet potato that will enhance the understanding of the mechanism of disease resistance.

Construction of a genetic map using EST-SSR markers and QTL analysis of major agronomic characters in hexaploid sweet potato (Ipomoea batatas (L.) Lam)

The Sweet potato, Ipomoea batatas (L.) Lam, is difficult to study in genetics and genomics because it is a hexaploid. The sweet potato study not have been performed domestically or internationally. In this study was performed to construct genetic map and quantitative trait loci (QTL) analysis. A total of 245 EST-SSR markers were developed, and the map was constructed by using 210 of those markers. The total map length was 1508.1 cM, and the mean distance between markers was 7.2 cM. Fifteen characteristics were investigated for QTLs analysis. According to those, the Four QTLs were identified, and The LOD score was 3.0. Further studies need to develop molecular markers in terms of EST-SSR markers for doing to be capable of efficient breeding. The genetic map created here using EST-SSR markers will facilitate planned breeding of sweet potato cultivars with various desirable traits.

Development of SSR Markers Based on Transcriptome Sequencing and Association Analysis with Drought Tolerance in Perennial Grass Miscanthus from China

Drought has become a critical environmental stress affecting on plant in temperate area. As one of the promising bio-energy crops to sustainable biomass production, the genus Miscanthus has been widely studied around the world. However, the most widely used hybrid cultivar among this genus, Miscanthus × giganteus is proved poor drought tolerance compared to some parental species. Here we mainly focused on Miscanthus sinensis, which is one of the progenitors of M. × giganteus providing a comparable yield and well abiotic stress tolerance in some places. The main objectives were to characterize the physiological and photosynthetic respond to drought stress and to develop simple sequence repeats (SSRs) markers associated with drought tolerance by transcriptome sequencing within an originally collection of 44 Miscanthus genotypes from southwest China. Significant phenotypic differences were observed among genotypes, and the average of leaf relative water content (RWC) were severely affected by drought stress decreasing from 88.27 to 43.21%, which could well contribute to separating the drought resistant and drought sensitive genotype of Miscanthus. Furthermore, a total of 16,566 gene-associated SSRs markers were identified based on Illumina RNA sequencing under drought conditions, and 93 of them were randomly selected to validate. In total, 70 (75.3%) SSRs were successfully amplified and the generated loci from 30 polymorphic SSRs were used to estimate the genetic differentiation and population structure. Finally, two optimum subgroups of the population were determined by structure analysis and based on association analysis, seven significant associations were identified including two markers with leaf RWC and five markers with photosynthetic traits. With the rich sequencing resources annotation, such associations would serve an efficient tool for Miscanthus drought response mechanism study and facilitate genetic improvement of drought resistant for this species.

Genetic analysis and association of simple sequence repeat markers with storage root yield, dry matter, starch and β-carotene content in sweetpotato

Molecular markers are needed for enhancing the development of elite sweetpotato (Ipomoea batatas (L.) Lam) cultivars with a wide range of commercially important traits in sub-Saharan Africa. This study was conducted to estimate the heritability and determine trait correlations of storage root yield, dry matter, starch and β-carotene content in a cross between 'New Kawogo' × 'Beauregard'. The study was also conducted to identify simple sequence repeat (SSR) markers associated with these traits. A total of 287 progeny and the parents were evaluated for two seasons at three sites in Uganda and genotyped with 250 SSR markers. Broad sense heritability (H2) for storage root yield, dry matter, starch and β-carotene content were 0.24, 0.68, 0.70 and 0.90, respectively. Storage root β-carotene content was negatively correlated with dry matter (r = -0.59, P < 0.001) and starch (r = -0.93, P < 0.001) content, while storage root yield was positively correlated with dry matter (r = 0.57, P = 0.029) and starch (r = 0.41, P = 0.008) content. Through logistic regression, a total of 12, 4, 6 and 8 SSR markers were associated with storage root yield, dry matter, starch and β-carotene content, respectively. The SSR markers used in this study may be useful for quantitative trait loci analysis and selection for these traits in future.

RNA-Seq of Guar (Cyamopsis tetragonoloba, L. Taub.) Leaves: De novo Transcriptome Assembly, Functional Annotation and Development of Genomic Resources

Genetic improvement in industrially important guar (Cyamopsis tetragonoloba, L. Taub.) crop has been hindered due to the lack of sufficient genomic or transcriptomic resources. In this study, RNA-Seq technology was employed to characterize the transcriptome of leaf tissues from two guar varieties, namely, M-83 and RGC-1066. Approximately 30 million high-quality pair-end reads of each variety generated by Illumina HiSeq platform were used for de novo assembly by Trinity program. A total of 62,146 non-redundant unigenes with an average length of 679 bp were obtained. The quality assessment of assembled unigenes revealed 87.50% of complete and 97.18% partial core eukaryotic genes (CEGs). Sequence similarity analyses and annotation of the unigenes against non-redundant protein (Nr) and Gene Ontology (GO) databases identified 175,882 GO annotations. A total of 11,308 guar unigenes were annotated with various enzyme codes (EC) and categorized in six categories with 55 subclasses. The annotation of biochemical pathways resulted in a total of 11,971 unigenes assigned with 145 KEGG maps and 1759 enzyme codes. The species distribution analysis of the unigenes showed highest similarity with Glycine max genes. A total of 5773 potential simple sequence repeats (SSRs) and 3594 high-quality single nucleotide polymorphisms (SNPs) were identified. Out of 20 randomly selected SSRs for wet laboratory validation, 13 showed consistent PCR amplification in both guar varieties. In silico studies identified 145 polymorphic SSR markers in two varieties. To the best of our knowledge, this is the first report on transcriptome analysis and SNPs identification in guar till date.

Comparative Transcriptome Analysis Reveals Critical Function of Sucrose Metabolism Related-Enzymes in Starch Accumulation in the Storage Root of Sweet Potato

The starch properties of the storage root (SR) affect the quality of sweet potato (Ipomoea batatas (L.) Lam.). Although numerous studies have analyzed the accumulation and properties of starch in sweet potato SRs, the transcriptomic variation associated with starch properties in SR has not been quantified. In this study, we measured the starch and sugar contents and analyzed the transcriptome profiles of SRs harvested from sweet potatoes with high, medium, and extremely low starch contents, at five developmental stages [65, 80, 95, 110, and 125 days after transplanting (DAP)]. We found that differences in both water content and starch accumulation in the dry matter affect the starch content of SRs in different sweet potato genotypes. Based on transcriptome sequencing data, we assembled 112336 unigenes, and identified several differentially expressed genes (DEGs) involved in starch and sucrose metabolism, and revealed the transcriptional regulatory network controlling starch and sucrose metabolism in sweet potato SRs. Correlation analysis between expression patterns and starch and sugar contents suggested that the sugar-starch conversion steps catalyzed by sucrose synthase (SuSy) and UDP-glucose pyrophosphorylase (UGPase) may be essential for starch accumulation in the dry matter of SRs, and IbβFRUCT2, a vacuolar acid invertase, might also be a key regulator of starch content in the SRs. Our results provide valuable resources for future investigations aimed at deciphering the molecular mechanisms determining the starch properties of sweet potato SRs.

Characterization and development of EST-SSR markers in sweet potato (Ipomoea batatas (L.) Lam)

In this study, a cDNA library was constructed from the total RNA of sweet potato leaves. A total of 789 copies of the cDNA were cloned in Escherichia coli by employing the pGEM-T Easy vector. Sequencing was carried out by Solgent Co. (Korea). As many as 579 expressed sequence tag-simple sequence repeat (EST-SSR) markers were designed (73.38%) from the known cDNA nucleotide base sequences. The lengths of the developed EST-SSR markers ranged from 100 to 499 bp (average length 238 bp). Their motif sequence types were varied, with most being dinucleotides and pentanucleotides, and the most commonly found motifs were CAGAAT (29.0%) and TCT (2.8%). Based on these SSR-containing sequences, 619 pairs of high-quality SSR primers were designed using WebSat and Primer3web. The total number of primers designed was 144. Polymorphism was evident in 82 EST-SSR markers among 20 Korean sweet potato cultivars tested and in 90 EST-SSR markers in the two parents of a mapping population, Yeseumi and Annobeny. In this study, the hexaploid sweet potato (2n = 6x = 90) EST-SSR markers were developed in the absence of full-sequence data. Moreover, by acting as a molecular tag for particular traits, the EST-SSR marker can also simultaneously identify information about the corresponding gene. These EST-SSR markers will allow the molecular analysis of sweet potato to be done more efficiently. Thus, we can develop high-quality sweet potato while overcoming the challenges from climate change and other unfavorable conditions.

A genome-wide BAC-end sequence survey provides first insights into sweetpotato (Ipomoea batatas (L.) Lam.) genome composition

Background

Sweetpotato, Ipomoea batatas (L.) Lam., is an important food crop widely grown in the world. However, little is known about the genome of this species because it is a highly heterozygous hexaploid. Gaining a more in-depth knowledge of sweetpotato genome is therefore necessary and imperative. In this study, the first bacterial artificial chromosome (BAC) library of sweetpotato was constructed. Clones from the BAC library were end-sequenced and analyzed to provide genome-wide information about this species.

Results

The BAC library contained 240,384 clones with an average insert size of 101 kb and had a 7.93–10.82 × coverage of the genome, and the probability of isolating any single-copy DNA sequence from the library was more than 99%. Both ends of 8310 BAC clones randomly selected from the library were sequenced to generate 11,542 high-quality BAC-end sequences (BESs), with an accumulative length of 7,595,261 bp and an average length of 658 bp. Analysis of the BESs revealed that 12.17% of the sweetpotato genome were known repetitive DNA, including 7.37% long terminal repeat (LTR) retrotransposons, 1.15% Non-LTR retrotransposons and 1.42% Class II DNA transposons etc., 18.31% of the genome were identified as sweetpotato-unique repetitive DNA and 10.00% of the genome were predicted to be coding regions. In total, 3,846 simple sequences repeats (SSRs) were identified, with a density of one SSR per 1.93 kb, from which 288 SSRs primers were designed and tested for length polymorphism using 20 sweetpotato accessions, 173 (60.07%) of them produced polymorphic bands. Sweetpotato BESs had significant hits to the genome sequences of I. trifida and more matches to the whole-genome sequences of Solanum lycopersicum than those of Vitis vinifera, Theobroma cacao and Arabidopsis thaliana.

Conclusions

The first BAC library for sweetpotato has been successfully constructed. The high quality BESs provide first insights into sweetpotato genome composition, and have significant hits to the genome sequences of I. trifida and more matches to the whole-genome sequences of Solanum lycopersicum. These resources as a robust platform will be used in high-resolution mapping, gene cloning, assembly of genome sequences, comparative genomics and evolution for sweetpotato.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3302-1) contains supplementary material, which is available to authorized users.

De novo Assembly and Characterization of the Transcriptome of Broomcorn Millet (Panicum miliaceum L.) for Gene Discovery and Marker Development

Broomcorn millet (Panicum miliaceum L.) is one of the world's oldest cultivated cereals, which is well-adapted to extreme environments such as drought, heat, and salinity with an efficient C4 carbon fixation. Discovery and identification of genes involved in these processes will provide valuable information to improve the crop for meeting the challenge of global climate change. However, the lack of genetic resources and genomic information make gene discovery and molecular mechanism studies very difficult. Here, we sequenced and assembled the transcriptome of broomcorn millet using Illumina sequencing technology. After sequencing, a total of 45,406,730 and 51,160,820 clean paired-end reads were obtained for two genotypes Yumi No. 2 and Yumi No. 3. These reads were mixed and then assembled into 113,643 unigenes, with the length ranging from 351 to 15,691 bp, of which 62,543 contings could be assigned to 315 gene ontology (GO) categories. Cluster of orthologous groups and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses assigned could map 15,514 unigenes into 202 KEGG pathways and 51,020 unigenes to 25 COG categories, respectively. Furthermore, 35,216 simple sequence repeats (SSRs) were identified in 27,055 unigene sequences, of which trinucleotides were the most abundant repeat unit, accounting for 66.72% of SSRs. In addition, 292 differentially expressed genes were identified between the two genotypes, which were significantly enriched in 88 GO terms and 12 KEGG pathways. Finally, the expression patterns of four selected transcripts were validated through quantitative reverse transcription polymerase chain reaction analysis. Our study for the first time sequenced and assembled the transcriptome of broomcorn millet, which not only provided a rich sequence resource for gene discovery and marker development in this important crop, but will also facilitate the further investigation of the molecular mechanism of its favored agronomic traits and beyond.

Development and Identification of SSR Markers Associated with Starch Properties and β-Carotene Content in the Storage Root of Sweet Potato (Ipomoea batatas L.)

Sweet potato (Ipomoea batatas L.) is a nutritious food crop and, based on the high starch content of its storage root, a potential bioethanol feedstock. Enhancing the nutritional value and starch quantity of storage roots are important goals of sweet potato breeding programs aimed at developing improved varieties for direct consumption, processing, and industrial uses. However, developing improved lines of sweet potato is challenging due to the genetic complexity of this plant and the lack of genome information. Short sequence repeat (SSR) markers are powerful molecular tools for tracking important loci in crops and for molecular-based breeding strategies; however, few SSR markers and marker-trait associations have hitherto been identified in sweet potato. In this study, we identified 1824 SSRs by using a de novo assembly of publicly available ESTs and mRNAs in sweet potato, and designed 1476 primer pairs based on SSR-containing sequences. We mapped 214 pairs of primers in a natural population comprised of 239 germplasms, and identified 1278 alleles with an average of 5.972 alleles per locus and a major allele frequency of 0.7702. Population structure analysis revealed two subpopulations in this panel of germplasms, and phenotypic characterization demonstrated that this panel is suitable for association mapping of starch-related traits. We identified 32, 16, and 17 SSR markers associated with starch content, β-carotene content, and starch composition in the storage root, respectively, using association analysis and further evaluation of a subset of sweet potato genotypes with various characteristics. The SSR markers identified here can be used to select varieties with desired traits and to investigate the genetic mechanism underlying starch and carotenoid formation in the starchy roots of sweet potato.

Characterization and Development of EST-SSRs by Deep Transcriptome Sequencing in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Simple sequence repeats (SSRs) are among the most important markers for population analysis and have been widely used in plant genetic mapping and molecular breeding. Expressed sequence tag-SSR (EST-SSR) markers, located in the coding regions, are potentially more efficient for QTL mapping, gene targeting, and marker-assisted breeding. In this study, we investigated 51,694 nonredundant unigenes, assembled from clean reads from deep transcriptome sequencing with a Solexa/Illumina platform, for identification and development of EST-SSRs in Chinese cabbage. In total, 10,420 EST-SSRs with over 12 bp were identified and characterized, among which 2744 EST-SSRs are new and 2317 are known ones showing polymorphism with previously reported SSRs. A total of 7877 PCR primer pairs for 1561 EST-SSR loci were designed, and primer pairs for twenty-four EST-SSRs were selected for primer evaluation. In nineteen EST-SSR loci (79.2%), amplicons were successfully generated with high quality. Seventeen (89.5%) showed polymorphism in twenty-four cultivars of Chinese cabbage. The polymorphic alleles of each polymorphic locus were sequenced, and the results showed that most polymorphisms were due to variations of SSR repeat motifs. The EST-SSRs identified and characterized in this study have important implications for developing new tools for genetics and molecular breeding in Chinese cabbage.

Temporal patterns of gene expression associated with tuberous root formation and development in sweetpotato (Ipomoea batatas)

BACKGROUND

The tuberous root of sweetpotato is undisputedly an important organ from agronomic and biological perspectives. Little is known regarding the regulatory networks programming tuberous root formation and development.

RESULTS

Here, as a first step toward understanding these networks, we analyzed and characterized the genome-wide transcriptional profiling and dynamics of sweetpotato root in seven distinct developmental stages using a customized microarray containing 39,724 genes. Analysis of these genes identified temporal programs of gene expression, including hundreds of transcription factor (TF) genes. We found that most genes active in roots were shared across all developmental stages, although significant quantitative changes in gene abundance were observed for 5,368 (including 435 TFs) genes. Clustering analysis of these differentially expressed genes pointed out six distinct expression patterns during root development. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that genes involved in different processes were enriched at specific stages of root development. In contrast with the large number of shared expressed genes in root development, each stage or period of root development has only a small number of specific genes. In total, 712 (including 27 TFs) and 1,840 (including 115 TFs) genes were identified as root-stage and root-period specific, respectively at the level of microarray. Several of the specific TF genes are known regulators of root development, including DA1-related protein, SHORT-ROOT and BEL1-like. The remaining TFs with unknown roles would also play critical regulatory roles during sweetpotato tuberous root formation and development.

CONCLUSIONS

The results generated in this study provided spatiotemporal patterns of root gene expression in support of future efforts for understanding the underlying molecular mechanism that control sweetpotato yield and quality.

De novo Assembly, Characterization of Immature Seed Transcriptome and Development of Genic-SSR Markers in Black Gram [Vigna mungo (L.) Hepper]

Black gram [V. mungo (L.) Hepper] is an important legume crop extensively grown in south and south-east Asia, where it is a major source of dietary protein for its predominantly vegetarian population. However, lack of genomic information and markers has become a limitation for genetic improvement of this crop. Here, we report the transcriptome sequencing of the immature seeds of black gram cv. TU94-2, by Illumina paired end sequencing technology to generate transcriptome sequences for gene discovery and genic-SSR marker development. A total of 17.2 million paired-end reads were generated and 48,291 transcript contigs (TCS) were assembled with an average length of 443 bp. Based on sequence similarity search, 33,766 TCS showed significant similarity to known proteins. Among these, only 29,564 TCS were annotated with gene ontology (GO) functional categories. A total number of 138 unique KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were identified, of which majority of TCS are grouped into purine metabolism (678) followed by pyrimidine metabolism (263). A total of 48,291 TCS were searched for SSRs and 1,840 SSRs were identified in 1,572 TCS with an average frequency of one SSR per 11.9 kb. The tri-nucleotide repeats were most abundant (35%) followed by di-nucleotide repeats (32%). PCR primer pairs were successfully designed for 933 SSR loci. Sequences analyses indicate that about 64.4% and 35.6% of the SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Tri-nucleotide repeats (57.3%) were preferentially present in the CDS. The rate of successful amplification and polymorphism were investigated using selected primers among 18 black gram accessions. Genic-SSR markers developed from the Illumina paired end sequencing of black gram immature seed transcriptome will provide a valuable resource for genetic diversity, evolution, linkage mapping, comparative genomics and marker-assisted selection in black gram.

Development and Characterization of Transcription Factor Gene-Derived Microsatellite (TFGM) Markers in Medicago truncatula and Their Transferability in Leguminous and Non-Leguminous Species

Transcription factors (TFs) are critical adaptor molecules that regulate many plant processes by controlling gene expression. The recent increase in the availability of TF data has made TFs a valuable resource for genic functional microsatellite marker development. In the present study, we developed TF gene-derived microsatellite (TFGM) markers for Medicago truncatula and assessed their cross-species transferability. A total of 203 SSRs were identified from 1467 M. truncatula TF coding sequences, 87.68% of which were trinucleotide repeats, followed by mono- (4.93%) and hexanucleotide repeats (1.48%). Further, 142 TFGM markers showed a high level of transferability to the leguminous (55.63%-85.21%) and non-leguminous (28.17%-50.00%) species. Polymorphisms of 27 TFGM markers were evaluated in 44 alfalfa accessions. The allele number per marker ranged from two to eight with an average of 4.41, and the PIC values ranged from 0.08 to 0.84 with an average of 0.60. Considering the high polymorphism, these TFGM markers developed in our study will be valuable for genetic relationship assessments, marker-assisted selection and comparative genomic studies in leguminous and non-leguminous species.

Characterization of EST-SSR markers in durum wheat EST library and functional analysis of SSR-containing EST fragments

The goal of this study is to identify characterization of expressed sequence tag (EST)-simple sequence repeats (SSR) markers from EST library of durum wheat and functional analysis of SSR-containing EST sequences for application in comparative genomics and breeding. 19,141 sequences were analyzed among which 18,937 ESTs were selected. Consistent with MISA results, 313 EST-SSRs were yielded. The final EST-SSRs were compared to the GenBank non-redundant database using BLASTX and classified based on these functions. Results indicated that the perfect EST-SSRs are the most frequent. The TTG/CTG imperfect EST-SSR had gamma-gliadin putative function that can be appropriate for durum wheat. Also, the mononucleotides and trinucleotides were the most frequent. Findings suggested that the identified EST-SSRs could be categorized into 83 types. Motifs TTG in trinucleotides and TC in dinucleotides had the highest frequency. TTG is the new motif in durum wheat identified in this study. We identified new EST-SSRs with more than trinucleotide and detected motifs that have potential to code amino acids. Arginine was the most frequent amino acid. Enzymes had the highest frequency among predicted functions. EST-SSRs have been identified in this study can be used for developing ESS-SSR-based detection tool for durum wheat in future studies and will be a useful resource for molecular breeding, genetics, genomics, and environmental stress studies. Motifs coding amino acids could be used as a new source of functional markers and biological study. In addition to, designed new PCR primer pairs are new resources for to identify useful alleles in transcription factors, storage proteins, and enzymes which incorporated them again into the cultivated material.

Development and characterization of simple sequence repeat (SSR) markers based on RNA-sequencing of Medicago sativa and in silico mapping onto the M. truncatula genome

Sufficient codominant genetic markers are needed for various genetic investigations in alfalfa since the species is an outcrossing autotetraploid. With the newly developed next generation sequencing technology, a large amount of transcribed sequences of alfalfa have been generated and are available for identifying SSR markers by data mining. A total of 54,278 alfalfa non-redundant unigenes were assembled through the Illumina HiSeqTM 2000 sequencing technology. Based on 3,903 unigene sequences, 4,493 SSRs were identified. Tri-nucleotide repeats (56.71%) were the most abundant motif class while AG/CT (21.7%), AGG/CCT (19.8%), AAC/GTT (10.3%), ATC/ATG (8.8%), and ACC/GGT (6.3%) were the subsequent top five nucleotide repeat motifs. Eight hundred and thirty- seven EST-SSR primer pairs were successfully designed. Of these, 527 (63%) primer pairs yielded clear and scored PCR products and 372 (70.6%) exhibited polymorphisms. High transferability was observed for ssp falcata at 99.2% (523) and 71.7% (378) in M. truncatula. In addition, 313 of 527 SSR marker sequences were in silico mapped onto the eight M. truncatula chromosomes. Thirty-six polymorphic SSR primer pairs were used in the genetic relatedness analysis of 30 Chinese alfalfa cultivated accessions generating a total of 199 scored alleles. The mean observed heterozygosity and polymorphic information content were 0.767 and 0.635, respectively. The codominant markers not only enriched the current resources of molecular markers in alfalfa, but also would facilitate targeted investigations in marker-trait association, QTL mapping, and genetic diversity analysis in alfalfa.

Development of novel microsatellite markers for the BBCC Oryza genome (Poaceae) using high-throughput sequencing technology

Wild species of Oryza are extremely valuable sources of genetic material that can be used to broaden the genetic background of cultivated rice, and to increase its resistance to abiotic and biotic stresses. Until recently, there was no sequence information for the BBCC Oryza genome; therefore, no special markers had been developed for this genome type. The lack of suitable markers made it difficult to search for valuable genes in the BBCC genome. The aim of this study was to develop microsatellite markers for the BBCC genome. We obtained 13,991 SSR-containing sequences and designed 14,508 primer pairs. The most abundant was hexanuclelotide (31.39%), followed by trinucleotide (27.67%) and dinucleotide (19.04%). 600 markers were selected for validation in 23 accessions of Oryza species with the BBCC genome. A set of 495 markers produced clear amplified fragments of the expected sizes. The average number of alleles per locus (Na) was 2.5, ranging from 1 to 9. The genetic diversity per locus (He) ranged from 0 to 0.844 with a mean of 0.333. The mean polymorphism information content (PIC) was 0.290, and ranged from 0 to 0.825. Of the 495 markers, 12 were only found in the BB genome, 173 were unique to the CC genome, and 198 were also present in the AA genome. These microsatellite markers could be used to evaluate the phylogenetic relationships among different Oryza genomes, and to construct a genetic linkage map for locating and identifying valuable genes in the BBCC genome, and would also for marker-assisted breeding programs that included accessions with the AA genome, especially Oryza sativa.

Transcriptome sequencing of two parental lines of cabbage (Brassica oleracea L. var. capitata L.) and construction of an EST-based genetic map

Background

Expressed sequence tag (EST)-based markers are preferred because they reflect transcribed portions of the genome. We report the development of simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers derived from transcriptome sequences in cabbage, and their utility for map construction.

Results

Transcriptome sequences were obtained from two cabbage parental lines, C1184 and C1234, which are susceptible and resistant to black rot disease, respectively, using the 454 platform. A total of 92,255 and 127,522 reads were generated and clustered into 34,688 and 40,947 unigenes, respectively. We identified 2,405 SSR motifs from the unigenes of the black rot-resistant parent C1234. Trinucleotide motifs were the most abundant (66.15%) among the repeat motifs. In addition, 1,167 SNPs were detected between the two parental lines. A total of 937 EST-based SSR and 97 SNP-based dCAPS markers were designed and used for detection of polymorphism between parents. Using an F₂ population, we built a genetic map comprising 265 loci, and consisting of 98 EST-based SSRs, 21 SNP-based dCAPS, 55 IBP markers derived from B. rapa genome sequence and 91 public SSRs, distributed on nine linkage groups spanning a total of 1,331.88 cM with an average distance of 5.03 cM between adjacent loci. The parental lines used in this study are elite breeding lines with little genetic diversity; therefore, the markers that mapped in our genetic map will have broad spectrum utility.

Conclusions

This genetic map provides additional genetic information to the existing B. oleracea map. Moreover, the new set of EST-based SSR and dCAPS markers developed herein is a valuable resource for genetic studies and will facilitate cabbage breeding. Additionally, this study demonstrates the usefulness of NGS transcriptomes for the development of genetic maps even with little genetic diversity in the mapping population.

Development of microsatellite markers by transcriptome sequencing in two species of Amorphophallus (Araceae)

Background

Amorphophallus is a genus of perennial plants widely distributed in the tropics or subtropics of West Africa and South Asia. Its corms contain a high level of water-soluble glucomannan; therefore, it has long been used as a medicinal herb and food source. Genetic studies of Amorphophallus have been hindered by a lack of genetic markers. A large number of molecular markers are required for genetic diversity study and improving disease resistance in Amorphophallus. Here, we report large scale of transcriptome sequencing of two species: Amorphophallus konjac and Amorphophallus bulbifer using deep sequencing technology, and microsatellite (SSR) markers were identified based on these transcriptome sequences.

Results

cDNAs of A. konjac and A. bulbifer were sequenced using Illumina HiSeq™ 2000 sequencing technology. A total of 135,822 non-redundant unigenes were assembled from about 9.66 gigabases, and 19,596 SSRs were identified in 16,027 non-redundant unigenes. Di-nucleotide SSRs were the most abundant motif (61.6%), followed by tri- (30.3%), tetra- (5.6%), penta- (1.5%), and hexa-nucleotides (1%) repeats. The top di- and tri-nucleotide repeat motifs included AG/CT (45.2%) and AGG/CCT (7.1%), respectively. A total of 10,754 primer pairs were designed for marker development. Of these, 320 primers were synthesized and used for validation of amplification and assessment of polymorphisms in 25 individual plants. The total of 275 primer pairs yielded PCR amplification products, of which 205 were polymorphic. The number of alleles ranged from 2 to 14 and the polymorphism information content valued ranged from 0.10 to 0.90. Genetic diversity analysis was done using 177 highly polymorphic SSR markers. A phenogram based on Jaccard’s similarity coefficients was constructed, which showed a distinct cluster of 25 Amorphophallus individuals.

Conclusion

A total of 10,754 SSR markers have been identified in Amorphophallus using transcriptome sequencing. One hundred and seventy-seven polymorphic markers were successfully validated in 25 individuals. The large number of genetic markers developed in the present study should contribute greatly to research into genetic diversity and germplasm characterization in Amorphophallus.

Development of microsatellite markers by transcriptome sequencing in two species of Amorphophallus (Araceae)

BACKGROUND

Amorphophallus is a genus of perennial plants widely distributed in the tropics or subtropics of West Africa and South Asia. Its corms contain a high level of water-soluble glucomannan; therefore, it has long been used as a medicinal herb and food source. Genetic studies of Amorphophallus have been hindered by a lack of genetic markers. A large number of molecular markers are required for genetic diversity study and improving disease resistance in Amorphophallus. Here, we report large scale of transcriptome sequencing of two species: Amorphophallus konjac and Amorphophallus bulbifer using deep sequencing technology, and microsatellite (SSR) markers were identified based on these transcriptome sequences.

RESULTS

cDNAs of A. konjac and A. bulbifer were sequenced using Illumina HiSeq™ 2000 sequencing technology. A total of 135,822 non-redundant unigenes were assembled from about 9.66 gigabases, and 19,596 SSRs were identified in 16,027 non-redundant unigenes. Di-nucleotide SSRs were the most abundant motif (61.6%), followed by tri- (30.3%), tetra- (5.6%), penta- (1.5%), and hexa-nucleotides (1%) repeats. The top di- and tri-nucleotide repeat motifs included AG/CT (45.2%) and AGG/CCT (7.1%), respectively. A total of 10,754 primer pairs were designed for marker development. Of these, 320 primers were synthesized and used for validation of amplification and assessment of polymorphisms in 25 individual plants. The total of 275 primer pairs yielded PCR amplification products, of which 205 were polymorphic. The number of alleles ranged from 2 to 14 and the polymorphism information content valued ranged from 0.10 to 0.90. Genetic diversity analysis was done using 177 highly polymorphic SSR markers. A phenogram based on Jaccard's similarity coefficients was constructed, which showed a distinct cluster of 25 Amorphophallus individuals.

CONCLUSION

A total of 10,754 SSR markers have been identified in Amorphophallus using transcriptome sequencing. One hundred and seventy-seven polymorphic markers were successfully validated in 25 individuals. The large number of genetic markers developed in the present study should contribute greatly to research into genetic diversity and germplasm characterization in Amorphophallus.

Development of simple sequence repeat markers and diversity analysis in alfalfa (Medicago sativa L.)

Efficient and robust molecular markers are essential for molecular breeding in plant. Compared to dominant and bi-allelic markers, multiple alleles of simple sequence repeat (SSR) markers are particularly informative and superior in genetic linkage map and QTL mapping in autotetraploid species like alfalfa. The objective of this study was to enrich SSR markers directly from alfalfa expressed sequence tags (ESTs). A total of 12,371 alfalfa ESTs were retrieved from the National Center for Biotechnology Information. Total 774 SSR-containing ESTs were identified from 716 ESTs. On average, one SSR was found per 7.7 kb of EST sequences. Tri-nucleotide repeats (48.8 %) was the most abundant motif type, followed by di-(26.1 %), tetra-(11.5 %), penta-(9.7 %), and hexanucleotide (3.9 %). One hundred EST-SSR primer pairs were successfully designed and 29 exhibited polymorphism among 28 alfalfa accessions. The allele number per marker ranged from two to 21 with an average of 6.8. The PIC values ranged from 0.195 to 0.896 with an average of 0.608, indicating a high level of polymorphism of the EST-SSR markers. Based on the 29 EST-SSR markers, assessment of genetic diversity was conducted and found that Medicago sativa ssp. sativa was clearly different from the other subspecies. The high transferability of those EST-SSR markers was also found for relative species.

Development and cross-species/genera transferability of microsatellite markers discovered using 454 genome sequencing in chokecherry (Prunus virginiana L.)

Chokecherry (Prunus virginiana L.) (2n = 4x = 32) is a unique Prunus species for both genetics and disease-resistance research due to its tetraploid nature and X-disease resistance. However, no genetic and genomic information on chokecherry is available. A partial chokecherry genome was sequenced using Roche 454 sequencing technology. A total of 145,094 reads covering 4.8 Mbp of the chokecherry genome were generated and 15,113 contigs were assembled, of which 11,675 contigs were larger than 100 bp in size. A total of 481 SSR loci were identified from 234 (out of 11,675) contigs and 246 polymerase chain reaction (PCR) primer pairs were designed. Of 246 primers, 212 (86.2 %) effectively produced amplification from the genomic DNA of chokecherry. All 212 amplifiable chokecherry primers were used to amplify genomic DNA from 11 other rosaceous species (sour cherry, sweet cherry, black cherry, peach, apricot, plum, apple, crabapple, pear, juneberry, and raspberry). Thus, chokecherry SSR primers can be transferable across Prunus species and other rosaceous species. An average of 63.2 and 58.7 % of amplifiable chokecherry primers amplified DNA from cherry and other Prunus species, respectively, while 47.2 % of amplifiable chokecherry primers amplified DNA from other rosaceous species. Using random genome sequence data generated from next-generation sequencing technology to identify microsatellite loci appears to be rapid and cost-efficient, particularly for species with no sequence information available. Sequence information and confirmed transferability of the identified chokecherry SSRs among species will be valuable for genetic research in Prunus and other rosaceous species. Key message A total of 246 SSR primers were identified from chokecherry genome sequences. Of which, 212 were confirmed amplifiable both in chokecherry and other 11 other rosaceous species.