Development of polymorphic markers from expressed sequence tags of Manihot esculenta Crantz

Six Express Sequence Tag-Simple Sequence Repeat Primers Reveal Genetic Diversity in the Cultivars of Three Zanthoxylum Species

Zanthoxylum (Sichuan pepper), with its rich cultivars, has long been widely cultivated in China for its unique seasoning and medicinal uses, but most of its cultivars have similar morphological characteristics. Therefore, we hypothesized that the genetic diversity of Zanthoxylum cultivars is low because of their apomixis and long cultivation history. In this study, we aimed to investigate the genetic diversity of three Zanthoxylum species on the cultivar level based on express sequence tag-simple sequence repeat (EST-SSR) primers. In total, 121 samples of three Zanthoxylum species (Z. bungeanum, Z. armatum and Z. piperitum) were collected from different areas in China for genetic diversity analysis. A total of six specificity and polymorphism EST-SSR primers, which we selected from among 120 primers based on two transcriptomes (Z. bungeanum, Z. armatum) in our earlier study, were used to evaluate genetic diversity based on capillary electrophoresis technology. The results of our analysis using the unweighted pair group method with arithmetic mean (UPGMA) indicated that most of the samples are clustered in one clade in the UPGMA dendrogram, and the average genetic distance was 0.6409. Principal component analysis (PCA) showed that Z. piperitum may have a closer genetic relationship with Z. bungeanum than with Z. armatum. An analysis of molecular variation (AMOVA) showed that the genetic variation mainly stemmed from individuals within populations; the genetic differentiation coefficient (PhiPT) was 0.429, the gene flow (Nm) between populations was 0.333, and the differences among populations were not significant (p > 0.001). For the intraspecific populations of ZB, the percentage of genetic variation was 53% among populations and 47% within populations, with non-significant differences between populations (p > 0.001). The genetic differentiation coefficient (PhiT) was 0.529, and the gene flow (Nm) was 0.223. For the intraspecific populations of ZA, the results indicated that the percentage of genetic variation was 29% among populations and 71% within populations, with non-significant differences between populations (p > 0.001); the genetic differentiation coefficient (PhiPT) was 0.293, and the gene flow (Nm) was 0.223. Through genetic structure analysis (GSA), we predicted that these 121 samples belonged to two optimal subgroups, which means that all the samples probably originated from two gene pools. Above all, this indicated that the genetic diversity of the 121 Zanthoxylum samples was relatively low at both the species and cultivar levels, a finding which was consistent with our initial assumptions. This study provides a reference, with molecular-level data, for the further identification of Zanthoxylum species.

Conventional breeding, marker-assisted selection, genomic selection and inbreeding in clonally propagated crops: a case study for cassava

Consolidates relevant molecular and phenotypic information on cassava to demonstrate relevance of heterosis, and alternatives to exploit it by integrating different tools. Ideas are useful to other asexually reproduced crops. Asexually propagated crops offer the advantage that all genetic effects can be exploited in farmers' production fields. However, non-additive effects complicate selection because, while influencing the performance of the materials under evaluation, they cannot be transmitted efficiently to the following cycle of selection. Cassava can be used as a model crop for asexually propagated crops because of its diploid nature and the absence of (known) incompatibility effects. New technologies such as genomic selection (GS), use of inbred progenitors based on doubled haploids and induction of flowering can be employed for accelerating genetic gains in cassava. Available information suggests that heterosis, non-additive genetic effects and within-family variation are relatively large for complex traits such as fresh root yield, moderate for dry matter or starch content in the roots, and low for defensive traits (pest and disease resistance) and plant architecture. The present article considers the potential impact of different technologies for maximizing gains for key traits in cassava, and highlights the advantages of integrating them. Exploiting heterosis would be optimized through the implementation of reciprocal recurrent selection. The advantages of using inbred progenitors would allow shifting the current cassava phenotypic recurrent selection method into line improvement, which in turn would allow designing outstanding hybrids rather than finding them by trial and error.

Identification and characterization of microsatellites in expressed sequence tags and their cross transferability in different plants

Expressed sequence tags (EST) are potential source for the development of genic microsatellite markers, gene discovery, comparative genomics, and other genomic studies. In the present study, 7630 ESTs were examined from NCBI for SSR identification and characterization. A total of 263 SSRs were identified with an average density of one SSR/4.2 kb (3.4% frequency). Analysis revealed that trinucleotide repeats (47.52%) were most abundant followed by tetranucleotide (19.77%), dinucleotide (19.01%), pentanucleotide (9.12%), and hexanucleotide repeats (4.56%). Functional annotation was done through homology search and gene ontology, and 35 EST-SSRs were selected. Primer pairs were designed for evaluation of cross transferability and polymorphism among 11 plants belonging to five different families. Total 402 alleles were generated at 155 loci with an average of 2.6 alleles/locus and the polymorphic information content (PIC) ranged from 0.15 to 0.92 with an average of 0.75. The cross transferability ranged from 34.84% to 98.06% in different plants, with an average of 67.86%. Thus, the validation study of annotated 35 EST-SSR markers which correspond to particular metabolic activity revealed polymorphism and evolutionary nature in different families of Angiospermic plants.

In silico genome comparison and distribution analysis of simple sequences repeats in cassava

We conducted a SSRs density analysis in different cassava genomic regions. The information obtained was useful to establish comparisons between cassava's SSRs genomic distribution and those of poplar, flax, and Jatropha. In general, cassava has a low SSR density (~50 SSRs/Mbp) and has a high proportion of pentanucleotides, (24,2 SSRs/Mbp). It was found that coding sequences have 15,5 SSRs/Mbp, introns have 82,3 SSRs/Mbp, 5' UTRs have 196,1 SSRs/Mbp, and 3' UTRs have 50,5 SSRs/Mbp. Through motif analysis of cassava's genome SSRs, the most abundant motif was AT/AT while in intron sequences and UTRs regions it was AG/CT. In addition, in coding sequences the motif AAG/CTT was also found to occur most frequently; in fact, it is the third most used codon in cassava. Sequences containing SSRs were classified according to their functional annotation of Gene Ontology categories. The identified SSRs here may be a valuable addition for genetic mapping and future studies in phylogenetic analyses and genomic evolution.

Development of simple sequence repeat markers from expressed sequence tags of the black tiger shrimp (Penaeus monodon)

In this study, we describe the development of expressed sequence tag-simple sequence repeat markers from expressed sequence tags of the black tiger shrimp (Penaeus monodon) deposited in public sequence databases. A total of 46 primer pairs were designed and screened on 26 individuals of P. monodon from a natural population. Of these, 16 primer pairs showed polymorphic profiles with between two and five alleles per locus. The average unbiased and direct count heterozygosities were 0.4662 and 0.3516, respectively. Cross-amplification was tested with five individuals of Penaeus vannamei and polymorphic products were detected at five loci.

The Cassava Genome: Current Progress, Future Directions

The starchy swollen roots of cassava provide an essential food source for nearly a billion people, as well as possibilities for bioenergy, yet improvements to nutritional content and resistance to threatening diseases are currently impeded. A 454-based whole genome shotgun sequence has been assembled, which covers 69% of the predicted genome size and 96% of protein-coding gene space, with genome finishing underway. The predicted 30,666 genes and 3,485 alternate splice forms are supported by 1.4 M expressed sequence tags (ESTs). Maps based on simple sequence repeat (SSR)-, and EST-derived single nucleotide polymorphisms (SNPs) already exist. Thanks to the genome sequence, a high-density linkage map is currently being developed from a cross between two diverse cassava cultivars: one susceptible to cassava brown streak disease; the other resistant. An efficient genotyping-by-sequencing (GBS) approach is being developed to catalog SNPs both within the mapping population and among diverse African farmer-preferred varieties of cassava. These resources will accelerate marker-assisted breeding programs, allowing improvements in disease-resistance and nutrition, and will help us understand the genetic basis for disease resistance.

Identification, validation and high-throughput genotyping of transcribed gene SNPs in cassava

The availability of genomic resources can facilitate progress in plant breeding through the application of advanced molecular technologies for crop improvement. This is particularly important in the case of less researched crops such as cassava, a staple and food security crop for more than 800 million people. Here, expressed sequence tags (ESTs) were generated from five drought stressed and well-watered cassava varieties. Two cDNA libraries were developed: one from root tissue (CASR), the other from leaf, stem and stem meristem tissue (CASL). Sequencing generated 706 contigs and 3,430 singletons. These sequences were combined with those from two other EST sequencing initiatives and filtered based on the sequence quality. Quality sequences were aligned using CAP3 and embedded in a Windows browser called HarvEST:Cassava which is made available. HarvEST:Cassava consists of a Unigene set of 22,903 quality sequences. A total of 2,954 putative SNPs were identified. Of these 1,536 SNPs from 1,170 contigs and 53 cassava genotypes were selected for SNP validation using Illumina's GoldenGate assay. As a result 1,190 SNPs were validated technically and biologically. The location of validated SNPs on scaffolds of the cassava genome sequence (v.4.1) is provided. A diversity assessment of 53 cassava varieties reveals some sub-structure based on the geographical origin, greater diversity in the Americas as opposed to Africa, and similar levels of diversity in West Africa and southern, eastern and central Africa. The resources presented allow for improved genetic dissection of economically important traits and the application of modern genomics-based approaches to cassava breeding and conservation.

Transcriptome sequencing of Hevea brasiliensis for development of microsatellite markers and construction of a genetic linkage map

To obtain more information on the Hevea brasiliensis genome, we sequenced the transcriptome from the vegetative shoot apex yielding 2 311 497 reads. Clustering and assembly of the reads produced a total of 113 313 unique sequences, comprising 28 387 isotigs and 84 926 singletons. Also, 17 819 expressed sequence tag (EST)-simple sequence repeats (SSRs) were identified from the data set. To demonstrate the use of this EST resource for marker development, primers were designed for 430 of the EST-SSRs. Three hundred and twenty-three primer pairs were amplifiable in H. brasiliensis clones. Polymorphic information content values of selected 47 SSRs among 20 H. brasiliensis clones ranged from 0.13 to 0.71, with an average of 0.51. A dendrogram of genetic similarities between the 20 H. brasiliensis clones using these 47 EST-SSRs suggested two distinct groups that correlated well with clone pedigree. These novel EST-SSRs together with the published SSRs were used for the construction of an integrated parental linkage map of H. brasiliensis based on 81 lines of an F1 mapping population. The map consisted of 97 loci, consisting of 37 novel EST-SSRs and 60 published SSRs, distributed on 23 linkage groups and covered 842.9 cM with a mean interval of 11.9 cM and ∼4 loci per linkage group. Although the numbers of linkage groups exceed the haploid number (18), but with several common markers between homologous linkage groups with the previous map indicated that the F1 map in this study is appropriate for further study in marker-assisted selection.

A genome scan for quantitative trait loci affecting cyanogenic potential of cassava root in an outbred population

Background

Cassava (Manihot esculenta Crantz) can produce cyanide, a toxic compound, without self-injury. That ability was called the cyanogenic potential (CN). This project aimed to identify quantitative trait loci (QTL) associated with the CN in an outbred population derived from 'Hanatee' × 'Huay Bong 60', two contrasting cultivars. CN was evaluated in 2008 and in 2009 at Rayong province, and in 2009 at Lop Buri province, Thailand. CN was measured using a picrate paper kit. QTL analysis affecting CN was performed with 303 SSR markers.

Results

The phenotypic values showed continuous variation with transgressive segregation events with more (115 ppm) and less CN (15 ppm) than either parent ('Hanatee' had 33 ppm and 'Huay Bong 60' had 95 ppm). The linkage map consisted of 303 SSR markers, on 27 linkage groups with a map that encompassed 1,328 cM. The average marker interval was 5.8 cM. Five QTL underlying CN were detected. CN08R1from 2008 at Rayong, CN09R1and CN09R2 from 2009 at Rayong, and CN09L1 and CN09L2 from 2009 at Lop Buri were mapped on linkage group 2, 5, 10 and 11, respectively. Among all the identified QTL, CN09R1 was the most significantly associated with the CN trait with LOD score 5.75 and explained the greatest percentage of phenotypic variation (%Expl.) of 26%.

Conclusions

Five new QTL affecting CN were successfully identified from 4 linkage groups. Discovery of these QTL can provide useful markers to assist in cassava breeding and studying genes affecting the trait.

SSR and EST-SSR-based genetic linkage map of cassava (Manihot esculenta Crantz)

Simple sequence repeat (SSR) markers provide a powerful tool for genetic linkage map construction that can be applied for identification of quantitative trait loci (QTL). In this study, a total of 640 new SSR markers were developed from an enriched genomic DNA library of the cassava variety 'Huay Bong 60' and 1,500 novel expressed sequence tag-simple sequence repeat (EST-SSR) loci were developed from the Genbank database. To construct a genetic linkage map of cassava, a 100 F(1) line mapping population was developed from the cross Huay Bong 60 by 'Hanatee'. Polymorphism screening between the parental lines revealed that 199 SSRs and 168 EST-SSRs were identified as novel polymorphic markers. Combining with previously developed SSRs, we report a linkage map consisted of 510 markers encompassing 1,420.3 cM, distributed on 23 linkage groups with a mean distance between markers of 4.54 cM. Comparison analysis of the SSR order on the cassava linkage map and the cassava genome sequences allowed us to locate 284 scaffolds on the genetic map. Although the number of linkage groups reported here revealed that this F(1) genetic linkage map is not yet a saturated map, it encompassed around 88% of the cassava genome indicating that the map was almost complete. Therefore, sufficient markers now exist to encompass most of the genomes and efficiently map traits in cassava.

Developing new SSR markers from ESTs of pea (Pisum sativum L.)

The development of expressed sequence tags (ESTs) from pea has provided a useful source for mining novel simple sequence repeat (SSR) markers. In the present research, in order to find EST-derived SSR markers, 18 552 pea ESTs from the National Center for Biotechnology Information (NCBI) database were downloaded and assembled into 10 086 unigenes. A total of 586 microsatellites in 530 unigenes were identified, indicating that merely 5.25% of sequences contained SSRs. The most abundant SSRs within pea were tri-nucleotide repeat motifs, and among all the tri-nucleotide repeats, the motif GAA was the most abundant type. In total, 49 SSRs were used for primer design. EST-SSR loci were subsequently screened on 10 widely adapted varieties in China. Of these, nine loci showed polymorphic profiles that revealed two to three alleles per locus. The polymorphism information content value ranged from 0.18 to 0.58 with an average of 0.41. Furthermore, transferable analysis revealed that some of these loci showed transferability to faba bean. Because of their polymorphism and transferability, these nine novel EST-SSRs will be valuable tools for marker-assisted breeding and comparative mapping of pea in the future.

Gene-based microsatellites for cassava (Manihot esculenta Crantz): prevalence, polymorphisms, and cross-taxa utility

BACKGROUND

Cassava (Manihot esculenta Crantz), a starchy root crop grown in tropical and subtropical climates, is the sixth most important crop in the world after wheat, rice, maize, potato and barley. The repertoire of simple sequence repeat (SSR) markers for cassava is limited and warrants a need for a larger number of polymorphic SSRs for germplasm characterization and breeding applications.

RESULTS

A total of 846 putative microsatellites were identified in silico from an 8,577 cassava unigene set with an average density of one SSR every 7 kb. One hundred and ninety-two candidate SSRs were screened for polymorphism among a panel of cassava cultivars from Africa, Latin America and Asia, four wild Manihot species as well as two other important taxa in the Euphorbiaceae, leafy spurge (Euphorbia esula) and castor bean (Ricinus communis). Of 168 markers with clean amplification products, 124 (73.8%) displayed polymorphism based on high resolution agarose gels. Of 85 EST-SSR markers screened, 80 (94.1%) amplified alleles from one or more wild species (M epruinosa, M glaziovii, M brachyandra, M tripartita) whereas 13 (15.3%) amplified alleles from castor bean and 9 (10.6%) amplified alleles from leafy spurge; hence nearly all markers were transferable to wild relatives of M esculenta while only a fraction was transferable to the more distantly related taxa. In a subset of 20 EST-SSRs assessed by fluorescence-based genotyping the number of alleles per locus ranged from 2 to 10 with an average of 4.55 per locus. These markers had a polymorphism information content (PIC) from 0.19 to 0.75 with an average value of 0.55 and showed genetic relationships consistent with existing information on these genotypes.

CONCLUSION

A set of 124 new, unique polymorphic EST-SSRs was developed and characterized which extends the repertoire of SSR markers for cultivated cassava and its wild relatives. The markers show high PIC values and therefore will be useful for cultivar identification, taxonomic studies, and genetic mapping. The study further shows that mining ESTs is a highly efficient strategy for polymorphism detection within the cultivated cassava gene pool.