Development and mapping of EST-derived simple sequence repeat markers for hexaploid wheat

Prebreeding studies of leaf rust resistant Triticum aestivum/T. timopheevii line L624

Triticum timopheevii Zhuk. attracts the attention of bread wheat breeders with its high immunity to the leaf rust pathogen. However, introgressions from this species in Triticum aestivum L. are little used in practical breeding. In the presented study, the agronomic value of T. aestivum/T. timopheevii line L624 was studied in comparison with the parent cultivars Saratovskaya 68, Dobrynya and the standard cultivar Favorit during 2017-2022. Introgressions from T. timopheevii in L624 were detected by the FISH method with probes pSc119.2, pAs1 and Spelt1, as well as microsatellite markers Xgwm312, Xgpw4480 and Xksum73. Translocations of 2AS.2AL-2AtL and on 2DL were detected as well. Line L624 is highly resistant to Puccinia triticina both under the background of natural epiphytotics and under laboratory conditions. PCR analysis with the DNA marker of the LrTt1 gene (Xgwm312) revealed that it is not identical to the Lr gene(s) in L624. According to a five-year study, the grain yield of L624 was, on average, higher than that of Favorit and Dobrynya, but lower than that of Saratovskaya 68. Line L624 had a lower weight of 1000 grains than the recipients, and was at the same level with the standard cultivar Favorit. Introgressions from T. timopheevii in L624 increased the grain protein content by comparison with Saratovskaya 68 and Favorit, but it was at the same level as in Dobrynya. As for parameters of flour and bread, L624 was not inferior to the recipient cultivars, but by volume and porosity of bread, it surpassed Saratovskaya 68. Moreover, L624 surpassed Favorit by the elasticity of the dough, the ratio of the elasticity of the dough to the extensibility and the strength of the flour. Thus, the results obtained suggest that introgressions in chromosomes 2A and 2D in L624 do not impair baking properties.

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

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

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

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

SUPPLEMENTARY INFORMATION

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

Fine mapping of wheat powdery mildew resistance gene Pm6 using 2B/2G homoeologous recombinants induced by the ph1b mutant

Using the ph1b mutant, the recombination frequency between the homoeologous region of 2B and 2G was significantly increased. By this, we narrowed Pm6 to a 0.9 Mb physical region. The powdery mildew (Pm) resistance gene Pm6 from Triticum timopheevii (2n = 48, AAGG) was mapped to the long arm of chromosome 2G and introduced into common wheat in the form of 2B-2G introgressions. The introgression line IGV1-465 has the shortest 2G segment, which is estimated 37 Mb in size when referring to 2BL genome reference of Chinese Spring (CS). The further fine mapping of Pm6 was impeded by the inhibition of allogeneic chromosome recombination between 2B and 2G in the Pm6 region. In the present study, to overcome 2B/2G recombination suppression, a ph1b-based strategy was employed to produce introgressions with reduced 2G fragments for the fine mapping of Pm6. IGV1-465 was crossed and backcrossed to the CSph1b mutant to produce plants with increased 2B/2G chromosome pairing frequency at the Pm6 region. A total of 182 allogeneic recombinants were obtained through two-round screening, i.e., first round of screening of 820 BC₁F_2:3 progenies using the flanking markers CIT02g-14/CIT02g-19 and second round of screening of 642 BC₁F_2:4 progenies using the flanking markers CIT02g-13/CIT02g-18, respectively. Through marker analysis using 30 chromosome 2G-specific markers located in the Pm6 region, the identified recombinants were divided into 14 haplotypes. Pm resistance evaluation of these haplotypes enabled us to narrow Pm6 to a 0.9 Mb physical region of 2BL, flanked by markers CIT02g-20 and CIT02g-18. Six wheat varieties containing Pm6 were identified from a natural population, and they showed increased Pm resistance. This implied Pm6 is still effective, especially when used in combination with other Pm resistance genes.

Development and validation of whole genome-wide and genic microsatellite markers in oil palm (Elaeis guineensis Jacq.): First microsatellite database (OpSatdb)

The availability of large expressed sequence tag (EST) and whole genome databases of oil palm enabled the development of a data base of microsatellite markers. For this purpose, an EST database consisting of 40,979 EST sequences spanning 27 Mb and a chromosome-wise whole genome databases were downloaded. A total of 3,950 primer pairs were identified and developed from EST sequences. The tri and tetra nucleotide repeat motifs were most prevalent (each 24.75%) followed by di-nucleotide repeat motifs. Whole genome-wide analysis found a total of 245,654 SSR repeats across the 16 chromosomes of oil palm, of which 38,717 were compound microsatellite repeats. A web application, OpSatdb, the first microsatellite database of oil palm, was developed using the PHP and MySQL database ( https://ssr.icar.gov.in/index.php ). It is a simple and systematic web-based search engine for searching SSRs based on repeat motif type, repeat type, and primer details. High synteny was observed between oil palm and rice genomes. The mapping of ESTs having SSRs by Blast2GO resulted in the identification of 19.2% sequences with gene ontology (GO) annotations. Randomly, a set of ten genic SSRs and five genomic SSRs were used for validation and genetic diversity on 100 genotypes belonging to the world oil palm genetic resources. The grouping pattern was observed to be broadly in accordance with the geographical origin of the genotypes. The identified genic and genome-wide SSRs can be effectively useful for various genomic applications of oil palm, such as genetic diversity, linkage map construction, mapping of QTLs, marker-assisted selection, and comparative population studies.

Molecular and Cytogenetic Characterization of Six Wheat-Aegilops markgrafii Disomic Addition Lines and Their Resistance to Rusts and Powdery Mildew

Aegilops markgrafii (Greuter) Hammer is an important source of genes for resistance to abiotic stresses and diseases in wheat (Triticum aestivum L.). A series of six wheat 'Alcedo'-Ae. markgrafii chromosome disomic addition lines, designated as AI(B), AII(C), AIII(D), AV(E), AIV(F), and AVIII(G) carrying the Ae. markgrafii chromosomes B, C, D, E, F, and G, respectively, were tested with SSR markers to establish homoeologous relationships to wheat and identify markers useful in chromosome engineering. The addition lines were evaluated for resistance to rust and powdery mildew diseases. The parents Alcedo and Ae. markgrafii accession 'S740-69' were tested with 1500 SSR primer pairs and 935 polymorphic markers were identified. After selecting for robust markers and confirming the polymorphisms on the addition lines, 132 markers were considered useful for engineering and establishing homoeologous relationships. Based on the marker analysis, we concluded that the chromosomes B, C, D, E, F, and G belong to wheat homoeologous groups 2, 5, 6, 7, 3, and 4, respectively. Also, we observed chromosomal rearrangements in several addition lines. When tested with 20 isolates of powdery mildew pathogen (Blumeria graminis f. sp. tritici) from five geographic regions of the United States, four addition lines [AIII(D), AV(E), AIV(F), and AVIII(G)] showed resistance to some isolates, with addition line AV(E) being resistant to 19 of 20 isolates. The addition lines were tested with two races (TDBJ and TNBJ) of the leaf rust pathogen (Puccinia triticina), and only addition line AI(B) exhibited resistance at a level comparable to the Ae. markgrafii parent. Addition lines AII(C) and AIII(D) had been previously identified as resistant to the Ug99 race group of the stem rust pathogen (Puccinia graminis f. sp. tritici). The addition lines were also tested for resistance to six United States races (PSTv-4, PSTv-14, PSTv-37, PSTv-40, PSTv-51, and PSTv-198) of the stripe rust pathogen (Puccinia striiformis f. sp. tritici); we found no resistance either in Alcedo or any of the addition lines. The homoeologous relationships of the chromosomes in the addition lines, molecular markers located on each chromosome, and disease resistance associated with each chromosome will allow for chromosome engineering of the resistance genes.

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.

Molecular Resources from Transcriptomes in the Brassicaceae Family

The rapidly falling costs and the increasing availability of large DNA sequence data sets facilitate the fast and affordable mining of large molecular markers data sets for comprehensive evolutionary studies. The Brassicaceae (mustards) are an important species-rich family in the plant kingdom with taxa distributed worldwide and a complex evolutionary history. We performed Simple Sequence Repeats (SSRs) mining using de novo assembled transcriptomes from 19 species across the Brassicaceae in order to study SSR evolution and provide comprehensive sets of molecular markers for genetic studies within the family. Moreover, we selected the genus Cochlearia to test the transferability and polymorphism of these markers among species. Additionally, we annotated Cochlearia pyrenaica transcriptome in order to identify the position of each of the mined SSRs. While we introduce a new set of tools that will further enable evolutionary studies across the Brassicaceae, we also discuss some broader aspects of SSR evolution. Overall, we developed 2012 ready-to-use SSR markers with their respective primers in 19 Brassicaceae species and a high quality annotated transcriptome for C. pyrenaica. As indicated by our transferability test with the genus Cochlearia these SSRs are transferable to species within the genus increasing exponentially the number of targeted species. Also, our polymorphism results showed substantial levels of variability for these markers. Finally, despite its complex evolutionary history, SSR evolution across the Brassicaceae family is highly conserved and we found no deviation from patterns reported in other Angiosperms.

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

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

ESAP plus: a web-based server for EST-SSR marker development

Background

Simple sequence repeats (SSRs) have become widely used as molecular markers in plant genetic studies due to their abundance, high allelic variation at each locus and simplicity to analyze using conventional PCR amplification. To study plants with unknown genome sequence, SSR markers from Expressed Sequence Tags (ESTs), which can be obtained from the plant mRNA (converted to cDNA), must be utilized. With the advent of high-throughput sequencing technology, huge EST sequence data have been generated and are now accessible from many public databases. However, SSR marker identification from a large in-house or public EST collection requires a computational pipeline that makes use of several standard bioinformatic tools to design high quality EST-SSR primers. Some of these computational tools are not users friendly and must be tightly integrated with reference genomic databases.

Results

A web-based bioinformatic pipeline, called EST Analysis Pipeline Plus (ESAP Plus), was constructed for assisting researchers to develop SSR markers from a large EST collection. ESAP Plus incorporates several bioinformatic scripts and some useful standard software tools necessary for the four main procedures of EST-SSR marker development, namely 1) pre-processing, 2) clustering and assembly, 3) SSR mining and 4) SSR primer design. The proposed pipeline also provides two alternative steps for reducing EST redundancy and identifying SSR loci. Using public sugarcane ESTs, ESAP Plus automatically executed the aforementioned computational pipeline via a simple web user interface, which was implemented using standard PHP, HTML, CSS and Java scripts. With ESAP Plus, users can upload raw EST data and choose various filtering options and parameters to analyze each of the four main procedures through this web interface. All input EST data and their predicted SSR results will be stored in the ESAP Plus MySQL database. Users will be notified via e-mail when the automatic process is completed and they can download all the results through the web interface.

Conclusions

ESAP Plus is a comprehensive and convenient web-based bioinformatic tool for SSR marker development. ESAP Plus offers all necessary EST-SSR development processes with various adjustable options that users can easily use to identify SSR markers from a large EST collection. With familiar web interface, users can upload the raw EST using the data submission page and visualize/download the corresponding EST-SSR information from within ESAP Plus. ESAP Plus can handle considerably large EST datasets. This EST-SSR discovery tool can be accessed directly from: http://gbp.kku.ac.th/esap_plus/.

Electronic supplementary material

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

Development of transcriptome-derived SSR markers for Hoya ledongensis (Apocynaceae) and cross-amplification in a congener

PREMISE OF THE STUDY

To examine patterns of genetic diversity and test possible hybridization events, microsatellite markers were identified and characterized in Hoya ledongensis (Apocynaceae), and cross-amplification was tested in a congener, H. jianfenglingensis.

METHODS AND RESULTS

Based on the transcriptome data of H. ledongensis, 46 microsatellite primer pairs were randomly selected for initial validation. From these, 28 primer pairs were successfully amplified, 12 of which were polymorphic in 36 individuals across three populations of H. ledongensis. The number of alleles per microsatellite locus ranged from two to 11. The observed and expected heterozygosities for the 12 loci ranged from 0.133 to 0.867 and 0.128 to 0.894, respectively. Cross-species amplification was successful for these 12 loci in the congeneric species H. jianfenglingensis.

CONCLUSIONS

These polymorphic transcriptome-derived simple sequence repeat markers have the potential to be used as multilocus molecular markers to study the population genetics and natural hybridization in species of Hoya.

Genetic analysis of a novel broad-spectrum powdery mildew resistance gene from the wheat-Agropyron cristatum introgression line Pubing 74

A novel broad-spectrum powdery mildew resistance gene PmPB74 was identified in wheat- Agropyron cristatum introgression line Pubing 74. Development of wheat cultivars with broad-spectrum, durable resistance to powdery mildew has been restricted by lack of superior genetic resources. In this study, a wheat-A. cristatum introgression line Pubing 74, originally selected from a wide cross between the common wheat cultivar Fukuhokomugi (Fukuho) and Agropyron cristatum (L.) Gaertn (2n = 4x = 28; genome PPPP), displayed resistance to powdery mildew at both the seedling and adult stages. The putative alien chromosomal fragment in Pubing 74 was below the detection limit of genomic in situ hybridization (GISH), but evidence for other non-GISH-detectable introgressions was provided by the presence of three STS markers specific to A. cristatum. Genetic analysis indicated that Pubing 74 carried a single dominant gene for powdery mildew resistance, temporarily designated PmPB74. Molecular mapping showed that PmPB74 was located on wheat chromosome arm 5DS, and flanked by markers Xcfd81 and HRM02 at genetic distances of 2.5 and 1.7 cM, respectively. Compared with other lines with powdery mildew resistance gene(s) on wheat chromosome arm 5DS, Pubing 74 was resistant to all 28 Blumeria graminis f. sp tritici (Bgt) isolates from different wheat-producing regions of northern China. Allelism tests indicated that PmPB74 was not allelic to PmPB3558 or Pm2. Our work showed that PmPB74 is a novel gene with broad resistance to powdery mildew, and hence will be helpful in broadening the genetic basis of powdery mildew resistance in wheat.

In silico search, characterization and validation of new EST-SSR markers in the genus Prunus

BACKGROUND

Simple sequence repeats (SSRs) are defined as sequence repeat units between 1 and 6 bp that occur in both coding and non-coding regions abundant in eukaryotic genomes, which may affect the expression of genes. In this study, expressed sequence tags (ESTs) of eight Prunus species were analyzed for in silico mining of EST-SSRs, protein annotation, and open reading frames (ORFs), and the identification of codon repetitions.

RESULTS

A total of 316 SSRs were identified using MISA software. Dinucleotide SSR motifs (26.31 %) were found to be the most abundant type of repeats, followed by tri- (14.58 %), tetra- (0.53 %), and penta- (0.27 %) nucleotide motifs. An attempt was made to design primer pairs for 316 identified SSRs but these were successful for only 175 SSR sequences. The positions of SSRs with respect to ORFs were detected, and annotation of sequences containing SSRs was performed to assign function to each sequence. SSRs were also characterized (in terms of position in the reference genome and associated gene) using the two available Prunus reference genomes (mei and peach). Finally, 38 SSR markers were validated across peach, almond, plum, and apricot genotypes. This validation showed a higher transferability level of EST-SSR developed in P. mume (mei) in comparison with the rest of species analyzed.

CONCLUSIONS

Findings will aid analysis of functionally important molecular markers and facilitate the analysis of genetic diversity.

Expressed Sequence Tags Analysis and Design of Simple Sequence Repeats Markers from a Full-Length cDNA Library in Perilla frutescens (L.)

Perilla frutescens is valuable as a medicinal plant as well as a natural medicine and functional food. However, comparative genomics analyses of P. frutescens are limited due to a lack of gene annotations and characterization. A full-length cDNA library from P. frutescens leaves was constructed to identify functional gene clusters and probable EST-SSR markers via analysis of 1,056 expressed sequence tags. Unigene assembly was performed using basic local alignment search tool (BLAST) homology searches and annotated Gene Ontology (GO). A total of 18 simple sequence repeats (SSRs) were designed as primer pairs. This study is the first to report comparative genomics and EST-SSR markers from P. frutescens will help gene discovery and provide an important source for functional genomics and molecular genetic research in this interesting medicinal plant.

Mining microsatellite markers from public expressed sequence tags databases for the study of threatened plants

BACKGROUND

Simple Sequence Repeats (SSRs) are widely used in population genetic studies but their classical development is costly and time-consuming. The ever-increasing available DNA datasets generated by high-throughput techniques offer an inexpensive alternative for SSRs discovery. Expressed Sequence Tags (ESTs) have been widely used as SSR source for plants of economic relevance but their application to non-model species is still modest.

METHODS

Here, we explored the use of publicly available ESTs (GenBank at the National Center for Biotechnology Information-NCBI) for SSRs development in non-model plants, focusing on genera listed by the International Union for the Conservation of Nature (IUCN). We also search two model genera with fully annotated genomes for EST-SSRs, Arabidopsis and Oryza, and used them as controls for genome distribution analyses. Overall, we downloaded 16 031 555 sequences for 258 plant genera which were mined for SSRsand their primers with the help of QDD1. Genome distribution analyses in Oryza and Arabidopsis were done by blasting the sequences with SSR against the Oryza sativa and Arabidopsis thaliana reference genomes implemented in the Basal Local Alignment Tool (BLAST) of the NCBI website. Finally, we performed an empirical test to determine the performance of our EST-SSRs in a few individuals from four species of two eudicot genera, Trifolium and Centaurea.

RESULTS

We explored a total of 14 498 726 EST sequences from the dbEST database (NCBI) in 257 plant genera from the IUCN Red List. We identify a very large number (17 102) of ready-to-test EST-SSRs in most plant genera (193) at no cost. Overall, dinucleotide and trinucleotide repeats were the prevalent types but the abundance of the various types of repeat differed between taxonomic groups. Control genomes revealed that trinucleotide repeats were mostly located in coding regions while dinucleotide repeats were largely associated with untranslated regions. Our results from the empirical test revealed considerable amplification success and transferability between congenerics.

CONCLUSIONS

The present work represents the first large-scale study developing SSRs by utilizing publicly accessible EST databases in threatened plants. Here we provide a very large number of ready-to-test EST-SSR (17 102) for 193 genera. The cross-species transferability suggests that the number of possible target species would be large. Since trinucleotide repeats are abundant and mainly linked to exons they might be useful in evolutionary and conservation studies. Altogether, our study highly supports the use of EST databases as an extremely affordable and fast alternative for SSR developing in threatened plants.

Identification and analysis of novel salt responsive candidate gene based SSRs (cgSSRs) from rice (Oryza sativa L.)

BACKGROUND

Majority of the Asian people depend on rice for nutritional energy. Rice cultivation and yield are severely affected by soil salinity stress worldwide. Marker assisted breeding is a rapid and efficient way to develop improved variety for salinity stress tolerance. Genomic microsatellite markers are an elite group of markers, but there is possible uncertainty of linkage with the important genes. In contrast, there are better possibilities of linkage detection with important genes if SSRs are developed from candidate genes. To the best of our knowledge, there is no such report on SSR markers development from candidate gene sequences in rice. So the present study was aimed to identify and analyse SSRs from salt responsive candidate genes of rice.

RESULTS

In the present study, based on the comprehensive literature survey, we selected 220 different salt responsive genes of rice. Out of them, 106 genes were found to contain 180 microsatellite loci with, tri-nucleotide motifs (56%) being most abundant, followed by di-(41%) and tetra nucleotide (2.8%) motifs. Maximum loci were found in the coding sequences (37.2%), followed by in 5'UTR (26%), intron (21.6%) and 3'UTR (15%). For validation, 19 primer sets were evaluated to detect polymorphism in diversity analysis among the two panels consisting of 17 salt tolerant and 17 susceptible rice genotypes. Except one, all primer sets exhibited polymorphic nature with an average of 21.8 alleles/primer and with a mean PIC value of 0.28. Calculated genetic similarity among genotypes was ranged from 19%-89%. The generated dendrogram showed 3 clusters of which one contained entire 17 susceptible genotypes and another two clusters contained all tolerant genotypes.

CONCLUSION

The present study represents the potential of salt responsive candidate gene based SSR (cgSSR) markers to be utilized as novel and remarkable candidate for diversity analysis among rice genotypes differing in salinity response.

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.

Identification, validation and cross-species transferability of novel Lavandula EST-SSRs

We identified and characterized EST-SSRs with strong discrimination power against Lavandula angustifolia and Lavandula x intermedia . The markers also showed considerable cross-species transferability rate into six related Lavandula species. Lavenders (Lavandula) are important economical crops grown around the globe for essential oil production. In an attempt to develop genetic markers for these plants, we analyzed over 13,000 unigenes developed from L. angustifolia and L. x intermedia EST databases, and identified 3,459 simple sequence repeats (SSR), which were dominated by trinucleotides (41.2 %) and dinucleotides (31.45 %). Approximately, 19 % of the unigenes contained at least one SSR marker, over 60 % of which were localized in the UTRs. Only 252 EST-SSRs were 18 bp or longer from which 31 loci were validated, and 24 amplified discrete fragments with 85 % polymorphism in L. x intermedia and L. angustifolia. The average number of alleles in L. x intermedia and L. angustifolia were 3.42 and 3.71 per marker with average PIC values of 0.47 and 0.52, respectively. These values suggest a moderate to strong level of informativeness for the markers, with some loci producing unique fingerprints. The cross-species transferability rate of the markers ranges 50-100 % across eight species. The utility of these markers was assessed in eight Lavandula species and 15 L. angustifolia and L. x intermedia cultivars, and the dendrogram deduced from their similarity indexes successfully delineated the species into their respective sections and the cultivars into their respective species. These markers have potential for application in fingerprinting, diversity studies and marker-assisted breeding of Lavandula.

QTL mapping of adult-plant resistance to leaf rust in a RIL population derived from a cross of wheat cultivars Shanghai 3/Catbird and Naxos

Six QTL for adult plant resistance to leaf rust, including two QTL effective against additional diseases, were identified in a RIL population derived from a cross between Shanghai 3/Catbird and Naxos. Leaf rust is an important wheat disease and utilization of adult-plant resistance (APR) may be the best approach to achieve long-term protection from the disease. The CIMMYT spring wheat line Shanghai 3/Catbird (SHA3/CBRD) showed a high level of APR to Chinese Puccinia triticina pathotypes in the field. To identify APR genes in this line, a mapping population of 164 recombinant inbred lines (RILs) was developed from a cross of this line and Naxos, a moderately susceptible German cultivar. The RILs were evaluated for final disease severity (FDS) at Baoding, Hebei province, and Zhoukou, Henan province, in the 2010-2011 and 2011-2012 cropping seasons. QTL analysis detected one major QTL derived from SHA3/CBRD on chromosome 2BS explaining from 15 to 37 % of the phenotypic variance across environments. In addition one minor resistance QTL on chromosome 1AL from SHA3/CBRD and four minor QTL from Naxos on chromosomes 2DL, 5B, 7BS, and 7DS were also detected. SHA3/CBRD also possessed seedling resistance gene Lr26, and Naxos contained Lr1 based on gene postulation following tests with an array of P. triticina pathotypes and molecular marker assays. These seedling resistance and APR genes and their closely linked molecular markers are potentially useful for improving leaf rust resistance in wheat breeding programs.

Identification, cross-taxon transferability and application of full-length cDNA SSR markers in Phyllostachys pubescens

Current databases of Phyllostachys pubescens full-length cDNAs (FL-cDNAs) provide a rich source of sequences for the development of potential FL-cDNA simple sequence repeat (SSR) markers. We screened 10,608 P. pubescens cDNAs, discovering 1614 SSRs in 1382 SSR-containing FL-cDNAs. The SSRs were more abundant within transposable elements (TEs) than expressed sequence tags (ESTs) and genome survey sequences (GSSs), and specific dinucleotide repeats tended to associate with particular TE families: (TA)n with En/Spm and (CT)n with Mutator. A selected panel of 100 FL-cDNAs containing type I SSRs yielded 68 functional SSR markers with an average polymorphism information content (PIC) value of 0.12, among which 22 loci contained polymorphisms. These markers became less transferrable (83.1% → 69.9% → 49.3%) but more polymorphic (79.4% → 92.3% → 92.8%) with increasing phylogenetic distance (intra-genus → intra-subtribe → intra-family). Transferability and polymorphism also depended on the location of the marker, with those located in the coding region being more transferrable (69.1%) and less polymorphic (89.4%) than those in the 5'-UTR (63.4% transferable, 90.7% polymorphic) and the 3'-UTR (61.8% transferable, 91.4% polymorphic). As proof of principle, we were able to use our FL-cDNA SSR markers to identify the parental stocks in interspecific hybrids of bamboo within and beyond P. pubescens, and estimate the outcrossing rate for P. pubescens. Our research should facilitate molecular breeding in bamboo species where original genetic markers are scarce.

Comparative genomics and association mapping approaches for blast resistant genes in finger millet using SSRs

The major limiting factor for production and productivity of finger millet crop is blast disease caused by Magnaporthe grisea. Since, the genome sequence information available in finger millet crop is scarce, comparative genomics plays a very important role in identification of genes/QTLs linked to the blast resistance genes using SSR markers. In the present study, a total of 58 genic SSRs were developed for use in genetic analysis of a global collection of 190 finger millet genotypes. The 58 SSRs yielded ninety five scorable alleles and the polymorphism information content varied from 0.186 to 0.677 at an average of 0.385. The gene diversity was in the range of 0.208 to 0.726 with an average of 0.487. Association mapping for blast resistance was done using 104 SSR markers which identified four QTLs for finger blast and one QTL for neck blast resistance. The genomic marker RM262 and genic marker FMBLEST32 were linked to finger blast disease at a P value of 0.007 and explained phenotypic variance (R²) of 10% and 8% respectively. The genomic marker UGEP81 was associated to finger blast at a P value of 0.009 and explained 7.5% of R². The QTLs for neck blast was associated with the genomic SSR marker UGEP18 at a P value of 0.01, which explained 11% of R². Three QTLs for blast resistance were found common by using both GLM and MLM approaches. The resistant alleles were found to be present mostly in the exotic genotypes. Among the genotypes of NW Himalayan region of India, VHC3997, VHC3996 and VHC3930 were found highly resistant, which may be effectively used as parents for developing blast resistant cultivars in the NW Himalayan region of India. The markers linked to the QTLs for blast resistance in the present study can be further used for cloning of the full length gene, fine mapping and their further use in the marker assisted breeding programmes for introgression of blast resistant alleles into locally adapted cultivars.

A novel genetic map of wheat: utility for mapping QTL for yield under different nitrogen treatments

Background

Common wheat (Triticum aestivum L.) is one of the most important food crops worldwide. Wheat varieties that maintain yield (YD) under moderate or even intense nitrogen (N) deficiency can adapt to low input management systems. A detailed genetic map is necessary for both wheat molecular breeding and genomics research. In this study, an F_6:7 recombinant inbred line population comprising 188 lines was used to construct a novel genetic map and subsequently to detect quantitative trait loci (QTL) for YD and response to N stress.

Results

A genetic map consisting of 591 loci distributed across 21 wheat chromosomes was constructed. The map spanned 3930.7 cM, with one marker per 6.7 cM on average. Genomic simple sequence repeat (g-SSR), expressed sequence tag-derived microsatellite (e-SSR), diversity arrays technology (DArT), sequence-tagged sites (STS), sequence-related amplified polymorphism (SRAP), and inter-simple sequence repeat (ISSR) molecular markers were included in the map. The linear relationships between loci found in the present map and in previously compiled physical maps were presented, which were generally in accordance. Information on the genetic and physical positions and allele sizes (when possible) of 17 DArT, 50 e-SSR, 44 SRAP, five ISSR, and two morphological markers is reported here for the first time. Seven segregation distortion regions (SDR) were identified on chromosomes 1B, 3BL, 4AL, 6AS, 6AL, 6BL, and 7B. A total of 22 and 12 QTLs for YD and yield difference between the value (YDDV) under HN and the value under LN were identified, respectively. Of these, QYd-4B-2 and QYddv-4B, two major stable QTL, shared support interval with alleles from KN9204 increasing YD in LN and decreasing YDDV. We probe into the use of these QTLs in wheat breeding programs. Moreover, factors affecting the SDR and total map length are discussed in depth.

Conclusions

This novel map may facilitate the use of novel markers in wheat molecular breeding programs and genomics research. Moreover, QTLs for YD and YDDV provide useful markers for wheat molecular breeding programs designed to increase yield potential under N stress.

Molecular and comparative mapping of genes governing spike compactness from wild emmer wheat

The development and morphology of the wheat spike is important because the spike is where reproduction occurs and it holds the grains until harvest. Therefore, genes that influence spike morphology are of interest from both theoretical and practical stand points. When substituted for the native chromosome 2A in the tetraploid Langdon (LDN) durum wheat background, the Triticum turgidum ssp. dicoccoides chromosome 2A from accession IsraelA confers a short, compact spike with fewer spikelets per spike compared to LDN. Molecular mapping and quantitative trait loci (QTL) analysis of these traits in a homozygous recombinant population derived from LDN × the chromosome 2A substitution line (LDNIsA-2A) indicated that the number of spikelets per spike and spike length were controlled by linked, but different, loci on the long arm of 2A. A QTL explaining most of the variation for spike compactness coincided with the QTL for spike length. Comparative mapping indicated that the QTL for number of spikelets per spike overlapped with a previously mapped QTL for Fusarium head blight susceptibility. The genes governing spike length and compactness were not orthologous to either sog or C, genes known to confer compact spikes in diploid and hexaploid wheat, respectively. Mapping and sequence analysis indicated that the gene governing spike length and compactness derived from wild emmer could be an ortholog of the barley Cly1/Zeo gene, which research indicates is an AP2-like gene pleiotropically affecting cleistogamy, flowering time, and rachis internode length. This work provides researchers with knowledge of new genetic loci and associated markers that may be useful for manipulating spike morphology in durum wheat.

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.

Genetic rearrangements of six wheat-agropyron cristatum 6P addition lines revealed by molecular markers

Agropyron cristatum (L.) Gaertn. (2n = 4x = 28, PPPP) not only is cultivated as pasture fodder but also could provide many desirable genes for wheat improvement. It is critical to obtain common wheat–A. cristatum alien disomic addition lines to locate the desired genes on the P genome chromosomes. Comparative analysis of the homoeologous relationships between the P genome chromosome and wheat genome chromosomes is a key step in transferring different desirable genes into common wheat and producing the desired alien translocation line while compensating for the loss of wheat chromatin. In this study, six common wheat–A. cristatum disomic addition lines were produced and analyzed by phenotypic examination, genomic in situ hybridization (GISH), SSR markers from the ABD genomes and STS markers from the P genome. Comparative maps, six in total, were generated and demonstrated that all six addition lines belonged to homoeologous group 6. However, chromosome 6P had undergone obvious rearrangements in different addition lines compared with the wheat chromosome, indicating that to obtain a genetic compensating alien translocation line, one should recombine alien chromosomal regions with homoeologous wheat chromosomes. Indeed, these addition lines were classified into four types based on the comparative mapping: 6P_I, 6P_II, 6P_III, and 6P_IV. The different types of chromosome 6P possessed different desirable genes. For example, the 6P_I type, containing three addition lines, carried genes conferring high numbers of kernels per spike and resistance to powdery mildew, important traits for wheat improvement. These results may prove valuable for promoting the development of conventional chromosome engineering techniques toward molecular chromosome engineering.

Construction of an integrative linkage map and QTL mapping of grain yield-related traits using three related wheat RIL populations

A novel high-density consensus wheat genetic map was obtained based on three related RIL populations, and the important chromosomal regions affecting yield and related traits were specified. A prerequisite for mapping quantitative trait locus (QTL) is to build a genetic linkage map. In this study, three recombinant inbred line populations (represented by WL, WY, and WJ) sharing one common parental line were used for map construction and subsequently for QTL detection of yield-related traits. PCR-based and diversity arrays technology markers were screened in the three populations. The integrated genetic map contains 1,127 marker loci, which span 2,976.75 cM for the whole genome, 985.93 cM for the A genome, 922.16 cM for the B genome, and 1,068.65 cM for the D genome. Phenotypic values were evaluated in four environments for populations WY and WJ, but three environments for population WL. Individual and combined phenotypic values across environments were used for QTL detection. A total of 165 putative additive QTL were identified, 22 of which showed significant additive-by-environment interaction effects. A total of 65 QTL (51.5%) were stable across environments, and 23 of these (35.4%) were common stable QTL that were identified in at least two populations. Notably, QTkw-5B.1, QTkw-6A.2, and QTkw-7B.1 were common major stable QTL in at least two populations, exhibiting 11.28-16.06, 5.64-18.69, and 6.76-21.16% of the phenotypic variance, respectively. Genetic relationships between kernel dimensions and kernel weight and between yield components and yield were evaluated. Moreover, QTL or regions that commonly interact across genetic backgrounds were discussed by comparing the results of the present study with those of previous similar studies. The present study provides useful information for marker-assisted selection in breeding wheat varieties with high yield.

Genetic analysis and molecular mapping of crown rust resistance in common wheat

This is the first report on genetic analysis and genome mapping of major dominant genes for near non-host resistance to barley crown rust ( Puccinia coronata var. hordei ) in common wheat. Barley crown rust, caused by Puccinia coronata var. hordei, primarily occurs on barley (Hordeum vulgare L.) in the Great Plain regions of the United States. However, a few genotypes of common wheat (Triticum aestivum L.) were susceptible to this pathogen among 750 wheat accessions evaluated. To investigate the genetics of crown rust resistance in wheat, a susceptible winter wheat accession PI 350005 was used in crosses with two resistant wheat varieties, Chinese Spring and Chris. Analysis of F1 plants and F2 populations from these two crosses indicated that crown rust resistance is controlled by one and two dominant genes in Chris and Chinese Spring, respectively. To determine the chromosome location of the resistance gene Cr1 in Chris, a set of 21 monosomic lines derived from Chris was used as female parents to cross with a susceptible spring type selection (SSTS35) derived from the PI 350005/Chris cross. Monosomic analysis indicated that Cr1 is located on chromosome 5D in Chris and one of the crown rust resistance genes is located on chromosome 2D in Chinese Spring. The other gene in Chinese Spring is not on 5D and thus is different from Cr1. Molecular linkage analysis and QTL mapping using a population of 136 doubled haploid lines derived from Chris/PI 350005 further positioned Cr1 between SSR markers Xwmc41-2 and Xgdm63 located on the long arm of chromosome 5D. Our study suggests that near non-host resistance to crown rust in these different common wheat genotypes is simply inherited.

Development of chloroplast simple sequence repeats (cpSSRs) for the intraspecific study of Gracilaria tenuistipitata (Gracilariales, Rhodophyta) from different populations

Background

Gracilaria tenuistipitata is an agarophyte with substantial economic potential because of its high growth rate and tolerance to a wide range of environment factors. This red seaweed is intensively cultured in China for the production of agar and fodder for abalone. Microsatellite markers were developed from the chloroplast genome of G. tenuistipitata var. liui to differentiate G. tenuistipitata obtained from six different localities: four from Peninsular Malaysia, one from Thailand and one from Vietnam. Eighty G. tenuistipitata specimens were analyzed using eight simple sequence repeat (SSR) primer-pairs that we developed for polymerase chain reaction (PCR) amplification.

Findings

Five mononucleotide primer-pairs and one trinucleotide primer-pair exhibited monomorphic alleles, whereas the other two primer-pairs separated the G. tenuistipitata specimens into two main clades. G. tenuistipitata from Thailand and Vietnam were grouped into one clade, and the populations from Batu Laut, Middle Banks and Kuah (Malaysia) were grouped into another clade. The combined dataset of these two primer-pairs separated G. tenuistipitata obtained from Kelantan, Malaysia from that obtained from other localities.

Conclusions

Based on the variations in repeated nucleotides of microsatellite markers, our results suggested that the populations of G. tenuistipitata were distributed into two main geographical regions: (i) populations in the west coast of Peninsular Malaysia and (ii) populations facing the South China Sea. The correct identification of G. tenuistipitata strains with traits of high economic potential will be advantageous for the mass cultivation of seaweeds.

In-silico mining, type and frequency analysis of genic microsatellites of finger millet (Eleusine coracana (L.) Gaertn.): a comparative genomic analysis of NBS-LRR regions of finger millet with rice

In recent years, the increased availability of the DNA sequences has given the possibility to develop and explore the expressed sequence tags (ESTs) derived SSR markers. In the present study, a total of 1956 ESTs of finger millet were used to find the microsatellite type, distribution, frequency and developed a total of 545 primer pairs from the ESTs of finger millet. Thirty-two EST sequences had more than two microsatellites and 1357 sequences did not have any SSR repeats. The most frequent type of repeats was trimeric motif, however the second place was occupied by dimeric motif followed by tetra-, hexa- and penta repeat motifs. The most common dimer repeat motif was GA and in case of trimeric SSRs, it was CGG. The EST sequences of NBS-LRR region of finger millet and rice showed higher synteny and were found on nearly same positions on the rice chromosome map. A total of eight, out of 15 EST based SSR primers were polymorphic among the selected resistant and susceptible finger millet genotypes. The primer FMBLEST5 could able to differentiate them into resistant and susceptible genotypes. The alleles specific to the resistant and susceptible genotypes were sequenced using the ABI 3130XL genetic analyzer and found similarity to NBS-LRR regions of rice and finger millet and contained the characteristic kinase-2 and kinase 3a motifs of plant R-genes belonged to NBS-LRR region. The In-silico and comparative analysis showed that the genes responsible for blast resistance can be identified, mapped and further introgressed through molecular breeding approaches for enhancing the blast resistance in finger millet.

Assessment of genetic diversity in the sorghum reference set using EST-SSR markers

Selection and use of genetically diverse genotypes are key factors in any crop breeding program to develop cultivars with a broad genetic base. Molecular markers play a major role in selecting diverse genotypes. In the present study, a reference set representing a wide range of sorghum genetic diversity was screened with 40 EST-SSR markers to validate both the use of these markers for genetic structure analyses and the population structure of this set. Grouping of accessions is identical in distance-based and model-based clustering methods. Genotypes were grouped primarily based on race within the geographic origins. Accessions derived from the African continent contributed 88.6 % of alleles confirming the African origin of sorghum. In total, 360 alleles were detected in the reference set with an average of 9 alleles per marker. The average PIC value was 0.5230 with a range of 0.1379-0.9483. Sub-race, guinea margaritiferum (Gma) from West Africa formed a separate cluster in close proximity to wild accessions suggesting that the Gma group represents an independent domestication event. Guineas from India and Western Africa formed two distinct clusters. Accessions belongs to the kafir race formed the most homogeneous group as observed in earlier studies. This analysis suggests that the EST-SSR markers used in the present study have greater discriminating power than the genomic SSRs. Genetic variance within the subpopulations was very high (71.7 %) suggesting that the germplasm lines included in the set are more diverse. Thus, this reference set representing the global germplasm is an ideal material for the breeding community, serving as a community resource for trait-specific allele mining as well as genome-wide association mapping.

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 Reproducible EST-derived SSR Markers and Assessment of Genetic Diversity in Panax ginseng Cultivars and Related Species

Little is known about the genetics or genomics of Panax ginseng. In this study, we developed 70 expressed sequence tag-derived polymorphic simple sequence repeat markers by trials of 140 primer pairs. All of the 70 markers showed reproducible polymorphism among four Panax speciesand 19 of them were polymorphic in six P. ginseng cultivars. These markers segregated 1:2:1 manner of Mendelian inheritance in an F2 population of a cross between two P. ginseng cultivars, 'Yunpoong' and 'Chunpoong', indicating that these are reproducible and inheritable mappable markers. A phylogenetic analysis using the genotype data showed three distinctive groups: a P. ginseng-P. japonicus clade, P. notoginseng and P. quinquefolius, with similarity coefficients of 0.70. P. japonicus was intermingled with P. ginseng cultivars, indicating that both species have similar genetic backgrounds. P. ginseng cultivars were subdivided into three minor groups: an independent cultivar 'Chunpoong', a subgroup with three accessions including two cultivars, 'Gumpoong' and 'Yunpoong' and one landrace 'Hwangsook' and another subgroup with two accessions including one cultivar, 'Gopoong' and one landrace 'Jakyung'. Each primer pair produced 1 to 4 bands, indicating that the ginseng genome has a highly replicated paleopolyploid genome structure.

Development, cross-species/genera transferability of novel EST-SSR markers and their utility in revealing population structure and genetic diversity in sugarcane

Sugarcane (Saccharum spp. hybrid) with complex polyploid genome requires a large number of informative DNA markers for various applications in genetics and breeding. Despite the great advances in genomic technology, it is observed in several crop species, especially in sugarcane, the availability of molecular tools such as microsatellite markers are limited. Now-a-days EST-SSR markers are preferred to genomic SSR (gSSR) as they represent only the functional part of the genome, which can be easily associated with desired trait. The present study was taken up with a new set of 351 EST-SSRs developed from the 4085 non redundant EST sequences of two Indian sugarcane cultivars. Among these EST-SSRs, TNR containing motifs were predominant with a frequency of 51.6%. Thirty percent EST-SSRs showed homology with annotated protein. A high frequency of SSRs was found in the 5'UTR and in the ORF (about 27%) and a low frequency was observed in the 3'UTR (about 8%). Two hundred twenty-seven EST-SSRs were evaluated, in sugarcane, allied genera of sugarcane and cereals, and 134 of these have revealed polymorphism with a range of PIC value 0.12 to 0.99. The cross transferability rate ranged from 87.0% to 93.4% in Saccharum complex, 80.0% to 87.0% in allied genera, and 76.0% to 80.0% in cereals. Cloning and sequencing of EST-SSR size variant amplicons revealed that the variation in the number of repeat-units was the main source of EST-SSR fragment polymorphism. When 124 sugarcane accessions were analyzed for population structure using model-based approach, seven genetically distinct groups or admixtures thereof were observed in sugarcane. Results of principal coordinate analysis or UPGMA to evaluate genetic relationships delineated also the 124 accessions into seven groups. Thus, a high level of polymorphism adequate genetic diversity and population structure assayed with the EST-SSR markers not only suggested their utility in various applications in genetics and genomics in sugarcane but also enriched the microsatellite marker resources in sugarcane.

SNP discovery and chromosome anchoring provide the first physically-anchored hexaploid oat map and reveal synteny with model species

A physically anchored consensus map is foundational to modern genomics research; however, construction of such a map in oat (Avena sativa L., 2n = 6x = 42) has been hindered by the size and complexity of the genome, the scarcity of robust molecular markers, and the lack of aneuploid stocks. Resources developed in this study include a modified SNP discovery method for complex genomes, a diverse set of oat SNP markers, and a novel chromosome-deficient SNP anchoring strategy. These resources were applied to build the first complete, physically-anchored consensus map of hexaploid oat. Approximately 11,000 high-confidence in silico SNPs were discovered based on nine million inter-varietal sequence reads of genomic and cDNA origin. GoldenGate genotyping of 3,072 SNP assays yielded 1,311 robust markers, of which 985 were mapped in 390 recombinant-inbred lines from six bi-parental mapping populations ranging in size from 49 to 97 progeny. The consensus map included 985 SNPs and 68 previously-published markers, resolving 21 linkage groups with a total map distance of 1,838.8 cM. Consensus linkage groups were assigned to 21 chromosomes using SNP deletion analysis of chromosome-deficient monosomic hybrid stocks. Alignments with sequenced genomes of rice and Brachypodium provide evidence for extensive conservation of genomic regions, and renewed encouragement for orthology-based genomic discovery in this important hexaploid species. These results also provide a framework for high-resolution genetic analysis in oat, and a model for marker development and map construction in other species with complex genomes and limited resources.

Generation and analysis of expressed sequence tags in the extreme large genomes Lilium and Tulipa

BACKGROUND

Bulbous flowers such as lily and tulip (Liliaceae family) are monocot perennial herbs that are economically very important ornamental plants worldwide. However, there are hardly any genetic studies performed and genomic resources are lacking. To build genomic resources and develop tools to speed up the breeding in both crops, next generation sequencing was implemented. We sequenced and assembled transcriptomes of four lily and five tulip genotypes using 454 pyro-sequencing technology.

RESULTS

Successfully, we developed the first set of 81,791 contigs with an average length of 514 bp for tulip, and enriched the very limited number of 3,329 available ESTs (Expressed Sequence Tags) for lily with 52,172 contigs with an average length of 555 bp. The contigs together with singletons covered on average 37% of lily and 39% of tulip estimated transcriptome. Mining lily and tulip sequence data for SSRs (Simple Sequence Repeats) showed that di-nucleotide repeats were twice more abundant in UTRs (UnTranslated Regions) compared to coding regions, while tri-nucleotide repeats were equally spread over coding and UTR regions. Two sets of single nucleotide polymorphism (SNP) markers suitable for high throughput genotyping were developed. In the first set, no SNPs flanking the target SNP (50 bp on either side) were allowed. In the second set, one SNP in the flanking regions was allowed, which resulted in a 2 to 3 fold increase in SNP marker numbers compared with the first set. Orthologous groups between the two flower bulbs: lily and tulip (12,017 groups) and among the three monocot species: lily, tulip, and rice (6,900 groups) were determined using OrthoMCL. Orthologous groups were screened for common SNP markers and EST-SSRs to study synteny between lily and tulip, which resulted in 113 common SNP markers and 292 common EST-SSR. Lily and tulip contigs generated were annotated and described according to Gene Ontology terminology.

CONCLUSIONS

Two transcriptome sets were built that are valuable resources for marker development, comparative genomic studies and candidate gene approaches. Next generation sequencing of leaf transcriptome is very effective; however, deeper sequencing and using more tissues and stages is advisable for extended comparative studies.

QTL mapping of adult-plant resistance to stripe rust in a population derived from common wheat cultivars Naxos and Shanghai 3/Catbird

Stripe rust, caused by Puccinia striiformis Westend. f. sp. tritici Erikss., is a severe foliar disease of common wheat (Triticum aestivum L.) worldwide. Use of adult-plant resistance (APR) is an efficient approach to provide long-term protection of crops from the disease. The German spring wheat cultivar Naxos showed a high level of APR to stripe rust in the field. To identify the APR genes in this cultivar, a mapping population of 166 recombinant inbred lines (RILs) was developed from a cross between Naxos and Shanghai 3/Catbird (SHA3/CBRD), a moderately susceptible line developed by CIMMYT. The RILs were evaluated for maximum disease severity (MDS) in Sichuan and Gansu in the 2009-2010 and 2010-2011 cropping seasons. Composite interval mapping (CIM) identified four QTL, QYr.caas-1BL.1RS, QYr.caas-1DS, QYr.caas-5BL.3 and QYr.caas-7BL.1, conferring stable resistance to stripe rust across all environments, each explaining 1.9-27.6, 2.1-5.8, 2.5-7.8 and 3.7-9.1 % of the phenotypic variance, respectively. QYr.caas-1DS flanked by molecular markers XUgwm353-Xgdm33b was likely a new QTL for APR to stripe rust. Because the interval between flanking markers for each QTL was less than 6.5 cM, these QTL and their closely linked markers are potentially useful for improving resistance to stripe rust in wheat breeding.

Durable resistance to stripe rust is due to three specific resistance genes in French bread wheat cultivar Apache

Quantitative resistance is postulated to be more durable than qualitative (R-gene mediated) resistance, which is usually quickly overcome by the pathogen population. Despite its wide use for nearly 10 years in France, the French bread wheat cultivar Apache remains resistant to stripe rust. Here, we investigated the genetic architecture of cv. Apache resistance to examine whether its durability could be explained by quantitative characteristics. We identified quantitative trait loci (QTL) by composite interval mapping of disease progress data recorded throughout 4 years of field assays. These assays included inoculation with three different pathotypes on a segregating population originating from a cross between cv. Apache and cv. Taldor, a French cultivar susceptible to stripe rust. Three QTLs derived from Apache, QYr.inra-2AS, QYr.inra-2BL and QYr.inra-4B, were detected. Each of these QTLs contributed between approximately 15 and 69 % of the phenotypic variance and corresponds to a race-specific resistance gene. We showed that QYr.inra-2AS and QYr.inra-2BS map to the positions of Yr17 and Yr7, respectively, whereas QYr.inra-4B corresponds to an adult plant resistance gene. Our results demonstrate that a combination of two or more race-specific resistance genes can confer durable resistance provided that it is properly managed at a continental level. Race-specific resistance genes should not be removed from breeding programs provided that they are properly managed.

Simple sequence repeat marker development from bacterial artificial chromosome end sequences and expressed sequence tags of flax (Linum usitatissimum L.)

Flax is an important oilseed crop in North America and is mostly grown as a fibre crop in Europe. As a self-pollinated diploid with a small estimated genome size of ~370 Mb, flax is well suited for fast progress in genomics. In the last few years, important genetic resources have been developed for this crop. Here, we describe the assessment and comparative analyses of 1,506 putative simple sequence repeats (SSRs) of which, 1,164 were derived from BAC-end sequences (BESs) and 342 from expressed sequence tags (ESTs). The SSRs were assessed on a panel of 16 flax accessions with 673 (58 %) and 145 (42 %) primer pairs being polymorphic in the BESs and ESTs, respectively. With 818 novel polymorphic SSR primer pairs reported in this study, the repertoire of available SSRs in flax has more than doubled from the combined total of 508 of all previous reports. Among nucleotide motifs, trinucleotides were the most abundant irrespective of the class, but dinucleotides were the most polymorphic. SSR length was also positively correlated with polymorphism. Two dinucleotide (AT/TA and AG/GA) and two trinucleotide (AAT/ATA/TAA and GAA/AGA/AAG) motifs and their iterations, different from those reported in many other crops, accounted for more than half of all the SSRs and were also more polymorphic (63.4 %) than the rest of the markers (42.7 %). This improved resource promises to be useful in genetic, quantitative trait loci (QTL) and association mapping as well as for anchoring the physical/genetic map with the whole genome shotgun reference sequence of flax.

Development and validation of genic-SSR markers in sesame by RNA-seq

Background

Sesame (Sesamum indicum L.) is one of the most important oil crops; however, a lack of useful molecular markers hinders current genetic research. We performed transcriptome sequencing of samples from different sesame growth and developmental stages, and mining of genic-SSR markers to identify valuable markers for sesame molecular genetics research.

Results

In this study, 75 bp and 100 bp paired-end RNA-seq was used to sequence 24 cDNA libraries, and 42,566 uni-transcripts were assembled from more than 260 million filtered reads. The total length of uni-transcript sequences was 47.99 Mb, and 7,324 SSRs (SSRs ≥15 bp) and 4,440 SSRs (SSRs ≥18 bp) were identified. On average, there was one genic-SSR per 6.55 kb (SSRs ≥15 bp) or 10.81 kb (SSRs ≥18 bp). Among perfect SSRs (≥18 bp), di-nucleotide motifs (48.01%) were the most abundant, followed by tri- (20.96%), hexa- (25.37%), penta- (2.97%), tetra- (2.12%), and mono-nucleotides (0.57%). The top four motif repeats were (AG/CT)n [1,268 (34.51%)], (CA/TG)n [281 (7.65%)], (AT/AT)n [215 (5.85%)], and (GAA/TTC)n [131 (3.57%)]. A total of 2,164 SSR primer pairs were identified in the 4,440 SSR-containing sequences (≥18 bp), and 300 SSR primer pairs were randomly chosen for validation. These SSR markers were amplified and validated in 25 sesame accessions (24 cultivated accessions, one wild species). 276 (92.0%) primer pairs yielded PCR amplification products in 24 cultivars. Thirty two primer pairs (11.59%) exhibited polymorphisms. Moreover, 203 primer pairs (67.67%) yielded PCR amplicons in the wild accession and 167 (60.51%) were polymorphic between species. A UPGMA dendrogram based on genetic similarity coefficients showed that the correlation between genotype and geographical source was low and that the genetic basis of sesame in China is narrow, as previously reported. The 32 polymorphic primer pairs were validated using an F₂ mapping population; 18 primer pairs exhibited polymorphisms between the parents, and 14 genic-SSRs could be integrated into 9 main linkage groups.

Conclusions

2,164 genic-SSR markers have been developed in sesame using transcriptome sequencing. 276 of 300 validated primer pairs successfully yielded PCR amplicons in 24 cultivated sesame accessions. These markers increase current SSR marker resources and will greatly benefit genetic diversity, qualitative and quantitative trait mapping and marker-assisted selection studies in sesame.

Molecular cytogenetic identification of a wheat (Triticum aestivum)-American dune grass (Leymus mollis) translocation line resistant to stripe rust

Leymus mollis, a perennial allotetraploid (2n = 4x = 28), known as American dune grass, is a wild relative of wheat that could be useful for cultivar improvement. Shannong0096, developed from interspecific hybridization between common wheat cv. Yannong15 and L. mollis, was analyzed with cytological procedures, genomic in situ hybridization, stripe-rust resistance screening and molecular marker analysis. We found that Shannong0096 has 42 chromosomes in the root-tip cells at mitotic metaphase and 21 bivalents in the pollen mother cells at meiotic metaphase I, demonstrating cytogenetic stability. Genomic in situ hybridization probed with total genomic DNA from L. mollis gave strong hybridization signals in the distal region of two wheat chromosome arms. A single dominant Yr gene, derived from L. mollis and temporarily designated as YrSn0096, was found on the long arm of chromosome 4A of Shannong0096. YrSn0096 should be a novel Yr gene because none of the previously reported Yr genes on chromosome 4A are related to L. mollis. This gene was found to be closely linked to the loci Xbarc236 and Xksum134 with genetic distances of 5.0 and 4.8 cM, respectively. Based on data from 267 F(2) plants of Yannong15/Huixianhong, the linkage map of YrSn0096, using the two molecular markers, was established in the order Xbarc236-YrSn0096-Xksum134. Shannong0096 appeared to be a unique wheat-L. mollis translocation with cryptic alien introgression. Cytogenetic stability, a high level of stripe-rust resistance, the common wheat background, and other positive agronomic traits make it a desirable donor for introducing novel alien resistance genes in wheat breeding programs, with the advantage of molecular markers that can be used to confirm introgression.

Identification of osmotic stress-responsive genes from Leymus mollis, a wild relative of wheat (Triticum aestivum L.)

The tolerance of the dune grass Leymus mollis (Triticeae; Poaceae) to various biotic and abiotic stresses makes it a very useful genetic resource for wheat breeding. Wide hybridization between L. mollis and wheat allows the introduction of Leymus chromosomes into the wheat genetic background and facilitates the integration of useful traits into wheat. However, the genetic basis controlling the physiological tolerance of L. mollis to multiple environmental stresses remains largely unexplored. Using suppression subtractive hybridization, we identified 112 osmotic-stress-responsive genes from L. mollis and confirmed their differential expression under osmotic stress. These genes were categorized into 13 functional categories, including cell defense and stress response, transcriptional regulation, signal transduction, biosynthesis of compatible solutes and cell wall metabolism. Representative genes were validated by northern blot and RT-PCR analyses of expression patterns in response to osmotic stress and abscisic acid treatment. The genes identified here represent a useful source of expressed sequence tags (ESTs) for the analysis and identification of Leymus chromosomes introduced into wheat. Furthermore, being highly conserved, genetically associated with osmotic stress tolerance and transferable to wheat, these ESTs provide significant tools for the development of EST-derived molecular markers for introgression of osmotic stress tolerance genes into wheat.

Whole-genome QTL analysis of Stagonospora nodorum blotch resistance and validation of the SnTox4-Snn4 interaction in hexaploid wheat

Necrotrophic effectors (also known as host-selective toxins) are important determinants of disease in the wheat-Stagonospora nodorum pathosystem. To date, five necrotrophic effector-host gene interactions have been identified in this system. Most of these interactions have additive effects while some are epistatic. The Snn4-SnTox4 interaction was originally identified in a recombinant-inbred population derived from a cross between the Swiss winter wheat cultivars 'Arina' and 'Forno' using the S. nodorum isolate Sn99CH 1A7a. Here, we used a recombinant-inbred population consisting of 121 lines developed from a cross between the hexaploid land race Salamouni and the hexaploid wheat 'Katepwa' (SK population). The SK population was used for the construction of linkage maps and quantitative trait loci (QTL) detection using the Swiss S. nodorum isolate Sn99CH 1A7a. The linkage maps developed in the SK population spanned 3,228 centimorgans (cM) and consisted of 441 simple-sequence repeats, 9 restriction fragment length polymorphisms, 29 expressed sequence tag sequence-tagged site markers, and 5 phenotypic markers. The average marker density was 6.7 cM/marker. Two QTL, designated QSnb.fcu-1A and QSnb.fcu-7A on chromosome arms 1AS and 7AS, respectively, were associated with disease caused by the S. nodorum isolate Sn99CH 1A7a. The effects of QSnb.fcu-1A were determined by the Snn4-SnTox4 interaction and accounted for 23.5% of the phenotypic variation in this population, whereas QSnb.fcu-7A accounted for 16.4% of the phenotypic variation for disease but was not associated with any known effector sensitivity locus. The effects of both QTL were largely additive and collectively accounted for 35.7% of the total phenotypic variation. The results of this research validate the effects of a compatible Snn4-SnTox4 interaction in a different genetic background, and it provides knowledge regarding genomic regions and molecular markers that can be used to improve Stagonospora nodorum blotch resistance in wheat germplasm.

Microsatellite mapping of Ae. speltoides and map-based comparative analysis of the S, G, and B genomes of Triticeae species

The first microsatellite linkage map of Ae. speltoides Tausch (2n = 2x = 14, SS), which is a wild species with a genome closely related to the B and G genomes of polyploid wheats, was developed based on two F(2) mapping populations using microsatellite (SSR) markers from Ae. speltoides, wheat genomic SSRs (g-SSRs) and EST-derived SSRs. A total of 144 different microsatellite loci were mapped in the Ae. speltoides genome. The transferability of the SSRs markers between the related S, B, and G genomes allowed possible integration of new markers into the T. timopheevii G genome chromosomal maps and map-based comparisons. Thirty-one new microsatellite loci assigned to the genetic framework of the T. timopheevii G genome maps were composed of wheat g-SSR (genomic SSR) markers. Most of the used Ae. speltoides SSRs were mapped onto chromosomes of the G genome supporting a close relationship between the G and S genomes. Comparative microsatellite mapping of the S, B, and G genomes demonstrated colinearity between the chromosomes within homoeologous groups, except for intergenomic T6A(t)S.1G, T4AL.5AL.7BS translocations. A translocation between chromosomes 2 and 6 that is present in the T. aestivum B genome was found in neither Ae. speltoides nor in T. timopheevii. Although the marker order was generally conserved among the B, S, and G genomes, the total length of the Ae. speltoides chromosomal maps and the genetic distances between homoeologous loci located in the proximal regions of the S genome chromosomes were reduced compared with the B, and G genome chromosomes.

Genetic map of artichoke × wild cardoon: toward a consensus map for Cynara cardunculus

An integrated consensus linkage map is proposed for globe artichoke. Maternal and paternal genetic maps were constructed on the basis of an F(1) progeny derived from crossing an artichoke genotype (Mola) with its progenitor, the wild cardoon (Tolfa), using EST-derived SSRs, genomic SSRs, AFLPs, ten genes, and two morphological traits. For most genes, mainly belonging to the chlorogenic acid pathway, new markers were developed. Five of these were SNP markers analyzed through high-resolution melt technology. From the maternal (Mola) and paternal (Tolfa) maps, an integrated map was obtained, containing 337 molecular and one morphological markers ordered in 17 linkage groups (LGs), linked between Mola and Tolfa. The integrated map covers 1,488.8 cM, with an average distance of 4.4 cM between markers. The map was aligned with already existing maps for artichoke, and 12 LGs were linked via 31 bridge markers. LG numbering has been proposed. A total of 124 EST-SSRs and two genes were mapped here for the first time, providing a framework for the construction of a functional map in artichoke. The establishment of a consensus map represents a necessary condition to plan a complete sequencing of the globe artichoke genome.