A genome-wide analysis of single nucleotide polymorphism diversity in the world's major cereal crops

Development and validation of an SNP genotyping array and construction of a high-density linkage map in castor

Castor is a commercially important oilseed crop that provides raw materials for several industries. Currently, the availability of genomic resources for castor is very limited. In this study, genome-wide SNPs were discovered in castor via whole-genome sequencing of 14 diverse lines to an average of 34X coverage. A total of 2,179,759 putative SNPs were detected, and a genotyping array was designed with 6,000 high-quality SNPs representing 2,492 scaffolds of the draft castor genome (87.5% genome coverage). The array was validated by genotyping a panel of 314 inbred castor lines, which resulted in 5,025 scorable SNPs with a high call rate (98%) and reproducibility (100%). Using this array, a consensus linkage map consisting of 1,978 SNP loci was constructed with an average inter-marker distance of 0.55 cM. The genome-wide SNP data, the genotyping array and the dense linkage map are valuable genomic tools for promoting high-throughput genomic research and molecular breeding in castor.

High Transferability of Homoeolog-Specific Markers between Bread Wheat and Newly Synthesized Hexaploid Wheat Lines

Bread wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) has a complex allohexaploid genome, which makes it difficult to differentiate between the homoeologous sequences and assign them to the chromosome A, B, or D subgenomes. The chromosome-based draft genome sequence of the ‘Chinese Spring’ common wheat cultivar enables the large-scale development of polymerase chain reaction (PCR)-based markers specific for homoeologs. Based on high-confidence ‘Chinese Spring’ genes with known functions, we developed 183 putative homoeolog-specific markers for chromosomes 4B and 7B. These markers were used in PCR assays for the 4B and 7B nullisomes and their euploid synthetic hexaploid wheat (SHW) line that was newly generated from a hybridization between Triticum turgidum (AABB) and the wild diploid species Aegilops tauschii (DD). Up to 64% of the markers for chromosomes 4B or 7B in the SHW background were confirmed to be homoeolog-specific. Thus, these markers were highly transferable between the ‘Chinese Spring’ bread wheat and SHW lines. Homoeolog-specific markers designed using genes with known functions may be useful for genetic investigations involving homoeologous chromosome tracking and homoeolog expression and interaction analyses.

Comparative transcriptomics and comprehensive marker resource development in mulberry

BACKGROUND

High potential of Morus laevigata and Morus serrata has been proposed in the breeding programs for Morus sp. However, due to the lack of dense molecular markers this goal is still in its nascent stage and not yet realized. We thus, sequenced the transcriptomes of these two wild Morus species and utilized the data for marker development.

RESULTS

We generated 87.0 and 80.3 Mb of transcriptome data from M. laevigata and M. serrata, respectively. The transcriptomes from M. laevigata and M. serrata, were assembled into 95,181 and 85,269 transcripts, respectively, and annotated. We identified around 24,049 Simple Sequence Repeats (SSRs), 1,201,326 Single Nucleotide Polymorphisms (SNPs) and 67,875 Insertion-Deletions (InDels). The variants having a higher impact were also identified and their effect was further investigated.

CONCLUSIONS

The transcriptome resource from the wildly growing mulberry species developed in this study can find wide applicability in gene identification and/or characterization. It can also contribute immensely in the existing mulberry improvement programs.

Comparison of SNP and CAPS markers application in genetic research in wheat and barley

BACKGROUND

Barley and bread wheat show large differences in frequencies of Single Nucleotide Polymorphism (SNP) as determined from genome-wide studies. These frequencies have been estimated as 2.4-3 times higher in the entire barley genome than within each diploid genomes of wheat (A, B or D). However, barley SNPs within individual genes occur significantly more frequently than quoted. Differences between wheat and barley are based on the origin and evolutionary history of the species. Bread wheat contains rarer SNPs due to the double genetic 'bottle-neck' created by natural hybridisation and spontaneous polyploidisation. Furthermore, wheat has the lowest level of useful SNP-derived markers while barley is estimated to have the highest level of polymorphism.

RESULTS

Different strategies are required for the development of suitable molecular markers in these cereal species. For example, SNP markers based on high-throughput technology (Infinium or KASP) are very effective and useful in both barley and bread wheat. In contrast, Cleaved Amplified Polymorphic Sequences (CAPS) are more widely and successfully employed in small-scale experiments with highly polymorphic genetic regions containing multiple SNPs in barley, but not in wheat. However, preliminary 'in silico' search databases for assessing the potential value of SNPs have yet to be developed.

CONCLUSIONS

This mini-review summarises results supporting the development of different strategies for the application of effective SNP and CAPS markers in wheat and barley.

De novo transcriptome sequence and identification of major bast-related genes involved in cellulose biosynthesis in jute (Corchorus capsularis L.)

Background

Jute fiber, extracted from stem bast, is called golden fiber. It is essential for fiber improvement to discover the genes associated with jute development at the vegetative growth stage. However, only 858 EST sequences of jute were deposited in the GenBank database. Obviously, the public available data is far from sufficient to understand the molecular mechanism of the fiber biosynthesis. It is imperative to conduct transcriptomic sequence for jute, which can be used for the discovery of a number of new genes, especially genes involved in cellulose biosynthesis.

Results

A total of 79,754,600 clean reads (7.98 Gb) were generated using Illumina paired-end sequencing. De novo assembly yielded 48,914 unigenes with an average length of 903 bp. By sequence similarity searching for known proteins, 27,962 (57.16 %) unigenes were annotated for their function. Out of these annotated unigenes, 21,856 and 11,190 unigenes were assigned to gene ontology (GO) and euKaryotic Ortholog Groups (KOG), respectively. Searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG) indicated that 14,216 unigenes were mapped to 268 KEGG pathways. Moreover, 5 Susy, 3 UGPase, 9 CesA, 18 CSL, 2 Kor (Korrigan), and 12 Cobra unigenes involving in cellulose biosynthesis were identified. Among these unigenes, the unigenes of comp11264_c0 (SuSy), comp24568_c0 (UGPase), comp11363_c0 (CesA), comp11363_c1 (CesA), comp24217_c0 (CesA), and comp23531_c0 (CesA), displayed relatively high expression level in stem bast using FPKM and RT-qPCR, indicating that they may have potential value of dissecting mechanism on cellulose biosynthesis in jute. In addition, a total of 12,518 putative gene-associate SNPs were called from these assembled uingenes.

Conclusion

We characterized the transcriptome of jute, discovered a broad survey of unigenes associated with vegetative growth and development, developed large-scale SNPs, and analyzed the expression patterns of genes involved in cellulose biosynthesis for bast fiber. All these provides a valuable genomics resource, which will accelerate the understanding of the mechanism of fiber development in jute.

Electronic supplementary material

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

Asymmetric somatic hybridization induces point mutations and indels in wheat

Background

Allopolyploid genome needs wide structural variation to deal with genomic shock. The introgression line, generated via asymmetric somatic hybridization, is introgressed with a minimum of exogenous chromatin, which also leads to genomic shock to induce genetic variation. However, the extent of its genomic variation and its difference from allopolyploidies remains unknown.

Methods

Here, we explored this issue using the bread wheat cultivar SR3, a derivative of an asymmetric somatic hybrid between the cultivar JN177 and an accession of tall wheatgrass (Thinopyrum elongatum). The ESTs (expressed sequence taqs) were large-scale sequenced using the cDNA library constructed in each of SR3 and JN177. Point mutations and indels (insertions and deletions) of SR3 were calculated, and their difference from the genetic variation of bread wheat and its ancestors were compared, with aim to analyze the extent and pattern of sequence variation induced by somatic hybridization.

Results

Both point mutations and indels (insertions and deletions) were frequently induced by somatic hybridization in the coding sequences. While the genomic shock caused by allopolyploidization tends to favor deletion over insertion, there was no evidence for such a preference following asymmetric somatic hybridization. The GC content of sequence adjacent to indel sites was also distinct from what has been observed in allopolyploids.

Conclusions

This study demonstrates that asymmetric somatic hybridization induces high frequency of genetic variation in a manner partially different from allopolipoidization. Asymmetric somatic hybridization provides appropriate material to comprehensively explore the nature of the genetic variation induced by genomic shock.

Electronic supplementary material

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

High density linkage mapping of genomic and transcriptomic SNPs for synteny analysis and anchoring the genome sequence of chickpea

This study presents genome-wide discovery of SNPs through next generation sequencing of the genome of Cicer reticulatum. Mapping of the C. reticulatum sequenced reads onto the draft genome assembly of C. arietinum (desi chickpea) resulted in identification of 842,104 genomic SNPs which were utilized along with an additional 36,446 genic SNPs identified from transcriptome sequences of the aforementioned varieties. Two new chickpea Oligo Pool All (OPAs) each having 3,072 SNPs were designed and utilized for SNP genotyping of 129 Recombinant Inbred Lines (RILs). Using Illumina GoldenGate Technology genotyping data of 5,041 SNPs were generated and combined with the 1,673 marker data from previously published studies, to generate a high resolution linkage map. The map comprised of 6698 markers distributed on eight linkage groups spanning 1083.93 cM with an average inter-marker distance of 0.16 cM. Utility of the present map was demonstrated for improving the anchoring of the earlier reported draft genome sequence of desi chickpea by ~30% and that of kabuli chickpea by 18%. The genetic map reported in this study represents the most dense linkage map of chickpea , with the potential to facilitate efficient anchoring of the draft genome sequences of desi as well as kabuli chickpea varieties.

A novel noninvasive procedure for high-throughput screening of major seed traits

The large numbers of samples processed in breeding and biodiversity programmes require the development of efficient methods for the nondestructive evaluation of basic seed properties. Near-infrared spectroscopy is the state-of-the-art solution for this analytical demand, but it also has some limitations. Here, we present a novel, rapid, accurate procedure based on time domain-nuclear magnetic resonance (TD-NMR), designed to simultaneously quantify a number of basic seed traits without any seed destruction. Using a low-field, benchtop (1) H-NMR instrument, the procedure gives a high-accuracy measurement of oil content (R(2) = 0.98), carbohydrate content (R(2) = 0.99), water content (R(2) = 0.98) and both fresh and dry weight of seeds/grains (R(2) = 0.99). The method requires a minimum of ~20 mg biomass per sample and thus enables to screen individual, intact seeds. When combined with an automated sample delivery system, a throughput of ~1400 samples per day is achievable. The procedure has been trialled as a proof of concept on cereal grains (collection of ~3000 accessions of Avena spp. curated at the IPK genebank). A mathematical multitrait selection approach has been designed to simplify the selection of outlying (most contrasting) accessions. To provide deeper insights into storage oil topology, some oat accessions were further analysed by three-dimensional seed modelling and lipid imaging. We conclude that the novel TD-NMR-based screening tool opens perspectives for breeding and plant biology in general.

Molecular marker applications in plants

Individuals within a population of a sexually reproducing species will have some degree of heritable genomic variation caused by mutations, insertion/deletions (INDELS), inversions, duplications, and translocations. Such variation can be detected and screened using molecular, or genetic, markers. By definition, molecular markers are genetic loci that can be easily tracked and quantified in a population and may be associated with a particular gene or trait of interest. This chapter will review the current major applications of molecular markers in plants.

Molecular marker databases

The detection and analysis of genetic variation plays an important role in plant breeding and this role is increasing with the continued development of genome sequencing technologies. Molecular genetic markers are important tools to characterize genetic variation and assist with genomic breeding. Processing and storing the growing abundance of molecular marker data being produced requires the development of specific bioinformatics tools and advanced databases. Molecular marker databases range from species specific through to organism wide and often host a variety of additional related genetic, genomic, or phenotypic information. In this chapter, we will present some of the features of plant molecular genetic marker databases, highlight the various types of marker resources, and predict the potential future direction of crop marker databases.

Genome-Wide Linkage Mapping of QTL for Yield Components, Plant Height and Yield-Related Physiological Traits in the Chinese Wheat Cross Zhou 8425B/Chinese Spring

Identification of genes for yield components, plant height (PH), and yield-related physiological traits and tightly linked molecular markers is of great importance in marker-assisted selection (MAS) in wheat breeding. In the present study, 246 F8 RILs derived from the cross of Zhou 8425B/Chinese Spring were genotyped using the high-density Illumina iSelect 90K single nucleotide polymorphism (SNP) assay. Field trials were conducted at Zhengzhou and Zhoukou of Henan Province, during the 2012-2013 and 2013-2014 cropping season under irrigated conditions, providing data for four environments. Analysis of variance (ANOVA) of agronomic and physiological traits revealed significant differences (P < 0.01) among RILs, environments, and RILs × environments interactions. Broad-sense heritabilities of all traits including thousand kernel weight (TKW), PH, spike length (SL), kernel number per spike (KNS), spike number/m(2) (SN), normalized difference in vegetation index at anthesis (NDVI-A) and at 10 days post-anthesis (NDVI-10), SPAD value of chlorophyll content at anthesis (Chl-A) and at 10 days post-anthesis (Chl-10) ranged between 0.65 and 0.94. A linkage map spanning 3609.4 cM was constructed using 5636 polymorphic SNP markers, with an average chromosome length of 171.9 cM and marker density of 0.64 cM/marker. A total of 866 SNP markers were newly mapped to the hexaploid wheat linkage map. Eighty-six QTL for yield components, PH, and yield-related physiological traits were detected on 18 chromosomes except 1D, 5D, and 6D, explaining 2.3-33.2% of the phenotypic variance. Ten stable QTL were identified across four environments, viz. QTKW.caas-6A.1, QTKW.caas-7AL, QKNS.caas-4AL, QSN.caas-1AL.1, QPH.caas-4BS.2, QPH.caas-4DS.1, QSL.caas-4AS, QSL.caas-4AL.1, QChl-A.caas-5AL, and QChl-10.caas-5BL. Meanwhile, 10 QTL-rich regions were found on chromosome 1BS, 2AL (2), 3AL, 4AL (2), 4BS, 4DS, 5BL, and 7AL exhibiting pleiotropic effects. These QTL or QTL clusters are tightly linked to SNP markers, with genetic distances to the closest SNPs ranging from 0 to 1.5 cM, and could serve as target regions for fine mapping, candidate gene discovery, and MAS in wheat breeding.

Using genic sequence capture in combination with a syntenic pseudo genome to map a deletion mutant in a wheat species

Mapping-by-sequencing analyses have largely required a complete reference sequence and employed whole genome re-sequencing. In species such as wheat, no finished genome reference sequence is available. Additionally, because of its large genome size (17 Gb), re-sequencing at sufficient depth of coverage is not practical. Here, we extend the utility of mapping by sequencing, developing a bespoke pipeline and algorithm to map an early-flowering locus in einkorn wheat (Triticum monococcum L.) that is closely related to the bread wheat genome A progenitor. We have developed a genomic enrichment approach using the gene-rich regions of hexaploid bread wheat to design a 110-Mbp NimbleGen SeqCap EZ in solution capture probe set, representing the majority of genes in wheat. Here, we use the capture probe set to enrich and sequence an F2 mapping population of the mutant. The mutant locus was identified in T. monococcum, which lacks a complete genome reference sequence, by mapping the enriched data set onto pseudo-chromosomes derived from the capture probe target sequence, with a long-range order of genes based on synteny of wheat with Brachypodium distachyon. Using this approach we are able to map the region and identify a set of deleted genes within the interval.

Application of next-generation sequencing technology to study genetic diversity and identify unique SNP markers in bread wheat from Kazakhstan

BACKGROUND

New SNP marker platforms offer the opportunity to investigate the relationships between wheat cultivars from different regions and assess the mechanism and processes that have led to adaptation to particular production environments. Wheat breeding has a long history in Kazakhstan and the aim of this study was to explore the relationship between key varieties from Kazakhstan and germplasm from breeding programs for other regions.

RESULTS

The study revealed 5,898 polymorphic markers amongst ten cultivars, of which 2,730 were mapped in the consensus genetic map. Mapped SNP markers were distributed almost equally across the A and B genomes, with between 279 and 484 markers assigned to each chromosome. Marker coverage was approximately 10-fold lower in the D genome. There were 863 SNP markers identified as unique to specific cultivars, and clusters of these markers (regions containing more than three closely mapped unique SNPs) showed specific patterns on the consensus genetic map for each cultivar. Significant intra-varietal genetic polymorphism was identified in three cultivars (Tzelinnaya 3C, Kazakhstanskaya rannespelaya and Kazakhstanskaya 15). Phylogenetic analysis based on inter-varietal polymorphism showed that the very old cultivar Erythrospermum 841 was the most genetically distinct from the other nine cultivars from Kazakhstan, falling in a clade together with the American cultivar Sonora and genotypes from Central and South Asia. The modern cultivar Kazakhstanskaya 19 also fell into a separate clade, together with the American cultivar Thatcher. The remaining eight cultivars shared a single sub-clade but were categorised into four clusters.

CONCLUSION

The accumulated data for SNP marker polymorphisms amongst bread wheat genotypes from Kazakhstan may be used for studying genetic diversity in bread wheat, with potential application for marker-assisted selection and the preparation of a set of genotype-specific markers.

Functional screening of willow alleles in Arabidopsis combined with QTL mapping in willow (Salix) identifies SxMAX4 as a coppicing response gene

Willows (Salix spp.) are important biomass crops due to their ability to grow rapidly with low fertilizer inputs and ease of cultivation in short-rotation coppice cycles. They are relatively undomesticated and highly diverse, but functional testing to identify useful allelic variation is time-consuming in trees and transformation is not yet possible in willow. Arabidopsis is heralded as a model plant from which knowledge can be transferred to advance the improvement of less tractable species. Here, knowledge and methodologies from Arabidopsis were successfully used to identify a gene influencing stem number in coppiced willows, a complex trait of key biological and industrial relevance. The strigolactone-related More AXillary growth (MAX) genes were considered candidates due to their role in shoot branching. We previously demonstrated that willow and Arabidopsis show similar response to strigolactone and that transformation rescue of Arabidopsis max mutants with willow genes could be used to detect allelic differences. Here, this approach was used to screen 45 SxMAX1, SxMAX2, SxMAX3 and SxMAX4 alleles cloned from 15 parents of 11 mapping populations varying in shoot-branching traits. Single-nucleotide polymorphism (SNP) frequencies were locus dependent, ranging from 29.2 to 74.3 polymorphic sites per kb. SxMAX alleles were 98%-99% conserved at the amino acid level, but different protein products varying in their ability to rescue Arabidopsis max mutants were identified. One poor rescuing allele, SxMAX4D, segregated in a willow mapping population where its presence was associated with increased shoot resprouting after coppicing and colocated with a QTL for this trait.

SNP discovery and high-density genetic mapping in faba bean (Vicia faba L.) permits identification of QTLs for ascochyta blight resistance

Ascochyta blight, caused by the fungus Ascochyta fabae Speg., is a common and destructive disease of faba bean (Vicia faba L.) on a global basis. Yield losses vary from typical values of 35-40% to 90% under specific environmental conditions. Several sources of resistance have been identified and used in breeding programs. However, introgression of the resistance gene determinants into commercial cultivars as a gene pyramiding approach is reliant on selection of closely linked genetic markers. A total of 14,552 base variants were identified from a faba bean expressed sequence tag (EST) database, and were further quality assessed to obtain a set of 822 high-quality single nucleotide polymorphisms (SNPs). Sub-sets of 336 EST-derived simple sequence repeats (SSRs) and 768 SNPs were further used for high-density genetic mapping of a biparental faba bean mapping population (Icarus×Ascot) that segregates for resistance to ascochyta blight. The linkage map spanned a total length of 1216.8 cM with 12 linkage groups (LGs) and an average marker interval distance of 2.3 cM. Comparison of map structure to the genomes of closely related legume species revealed a high degree of conserved macrosynteny, as well as some rearrangements. Based on glasshouse evaluation of ascochyta blight resistance performed over two years, four genomic regions controlling resistance were identified on Chr-II, Chr-VI and two regions on Chr-I.A. Of these, one (QTL-3) may be identical with quantitative trait loci (QTLs) identified in prior studies, while the others (QTL-1, QTL-2 and QTL-4) may be novel. Markers in close linkage to ascochyta blight resistance genes identified in this study can be further validated and effectively implemented in faba bean breeding programs.

SNP marker discovery, linkage map construction and identification of QTLs for enhanced salinity tolerance in field pea (Pisum sativum L.)

BACKGROUND

Field pea (Pisum sativum L.) is a self-pollinating, diploid, cool-season food legume. Crop production is constrained by multiple biotic and abiotic stress factors, including salinity, that cause reduced growth and yield. Recent advances in genomics have permitted the development of low-cost high-throughput genotyping systems, allowing the construction of saturated genetic linkage maps for identification of quantitative trait loci (QTLs) associated with traits of interest. Genetic markers in close linkage with the relevant genomic regions may then be implemented in varietal improvement programs.

RESULTS

In this study, single nucleotide polymorphism (SNP) markers associated with expressed sequence tags (ESTs) were developed and used to generate comprehensive linkage maps for field pea. From a set of 36,188 variant nucleotide positions detected through in silico analysis, 768 were selected for genotyping of a recombinant inbred line (RIL) population. A total of 705 SNPs (91.7%) successfully detected segregating polymorphisms. In addition to SNPs, genomic and EST-derived simple sequence repeats (SSRs) were assigned to the genetic map in order to obtain an evenly distributed genome-wide coverage. Sequences associated with the mapped molecular markers were used for comparative genomic analysis with other legume species. Higher levels of conserved synteny were observed with the genomes of Medicago truncatula Gaertn. and chickpea (Cicer arietinum L.) than with soybean (Glycine max [L.] Merr.), Lotus japonicus L. and pigeon pea (Cajanus cajan [L.] Millsp.). Parents and RIL progeny were screened at the seedling growth stage for responses to salinity stress, imposed by addition of NaCl in the watering solution at a concentration of 18 dS m-1. Salinity-induced symptoms showed normal distribution, and the severity of the symptoms increased over time. QTLs for salinity tolerance were identified on linkage groups Ps III and VII, with flanking SNP markers suitable for selection of resistant cultivars. Comparison of sequences underpinning these SNP markers to the M. truncatula genome defined genomic regions containing candidate genes associated with saline stress tolerance.

CONCLUSION

The SNP assays and associated genetic linkage maps developed in this study permitted identification of salinity tolerance QTLs and candidate genes. This constitutes an important set of tools for marker-assisted selection (MAS) programs aimed at performance enhancement of field pea cultivars.

The genetic study utility of a hexaploid wheat DH population with non-recombinant A- and B-genomes

To study the D-genome of the wild wheat relative Aegilops tauschii Cosson at the hexaploid level, we developed a synthetic doubled-haploid (DH) hexaploid wheat population, SynDH3. This population was derived from the spontaneous chromosome doubling of triploid F1 hybrid plants obtained from a cross between Triticum turgidum ssp. dicoccon PI377655 and A. tauschii ssp. strangulata AS66 × ssp. tauschii AS87. SynDH3 is a diploidization-hexaploid DH population containing recombinant D chromosomes from two different A. tauschii genotypes, with A and B chromosomes from T. turgidum being homogenous across the entire population. Using this population, we constructed a genetic map. Of the 440 markers used to construct the map, 421 (96%) were assigned to 12 linkage groups; these included 346 Diversity Arrays Technology (DArT) and 75 simple sequence repeat (SSR) markers. The total map length of the seven D chromosomes spanned 916.27 cM, with an average length of 130.90 cM per chromosome and an average distance between markers of 3.47 cM. Seven segregation distortion regions were detected on seven linkage groups. Out of 50 markers shared with those on a common wheat map, 37 showed a consistent order. The utility of the diploidization-hexaploid DH population for mapping qualitative trait genes was confirmed using the dominant glaucousness-inhibiting gene W2 (I) as an example.

Genotype-specific SNP map based on whole chromosome 3B sequence information from wheat cultivars Arina and Forno

Agronomically important traits are frequently controlled by rare, genotype-specific alleles. Such genes can only be mapped in a population derived from the donor genotype. This requires the development of a specific genetic map, which is difficult in wheat because of the low level of polymorphism among elite cultivars. The absence of sufficient polymorphism, the complexity of the hexaploid wheat genome as well as the lack of complete sequence information make the construction of genetic maps with a high density of reproducible and polymorphic markers challenging. We developed a genotype-specific genetic map of chromosome 3B from winter wheat cultivars Arina and Forno. Chromosome 3B was isolated from the two cultivars and then sequenced to 10-fold coverage. This resulted in a single-nucleotide polymorphisms (SNP) database of the complete chromosome. Based on proposed synteny with the Brachypodium model genome and gene annotation, sequences close to coding regions were used for the development of 70 SNP-based markers. They were mapped on a Arina × Forno Recombinant Inbred Lines population and found to be spread over the complete chromosome 3B. While overall synteny was well maintained, numerous exceptions and inversions of syntenic gene order were identified. Additionally, we found that the majority of recombination events occurred in distal parts of chromosome 3B, particularly in hot-spot regions. Compared with the earlier map based on SSR and RFLP markers, the number of markers increased fourfold. The approach presented here allows fast development of genotype-specific polymorphic markers that can be used for mapping and marker-assisted selection.

Comparative analysis of kabuli chickpea transcriptome with desi and wild chickpea provides a rich resource for development of functional markers

Chickpea (Cicer arietinum L.) is an important crop legume plant with high nutritional value. The transcriptomes of desi and wild chickpea have already been sequenced. In this study, we sequenced the transcriptome of kabuli chickpea, C. arietinum (genotype ICCV2), having higher commercial value, using GS-FLX Roche 454 and Illumina technologies. The assemblies of both Roche 454 and Illumina datasets were optimized using various assembly programs and parameters. The final optimized hybrid assembly generated 43,389 transcripts with an average length of 1065 bp and N50 length of 1653 bp representing 46.2 Mb of kabuli chickpea transcriptome. We identified a total of 5409 simple sequence repeats (SSRs) in these transcript sequences. Among these, at least 130 and 493 SSRs were polymorphic with desi (ICC4958) and wild (PI489777) chickpea, respectively. In addition, a total of 1986 and 37,954 single nucleotide polymorphisms (SNPs) were predicted in kabuli/desi and kabuli/wild genotypes, respectively. The SNP frequency was 0.043 SNP per kb for kabuli/desi and 0.821 SNP per kb for kabuli/wild, reflecting very low genetic diversity in chickpea. Further, SSRs and SNPs present in tissue-specific and transcription factor encoding transcripts have been identified. The experimental validation of a selected set of polymorphic SSRs and SNPs exhibited high intra-specific polymorphism potential between desi and kabuli chickpea, suggesting their utility in large-scale genotyping applications. The kabuli chickpea gene index assembled, and SSRs and SNPs identified in this study will serve as useful genomic resource for genetic improvement of chickpea.

Genome wide SNP discovery in flax through next generation sequencing of reduced representation libraries

BACKGROUND

Flax (Linum usitatissimum L.) is a significant fibre and oilseed crop. Current flax molecular markers, including isozymes, RAPDs, AFLPs and SSRs are of limited use in the construction of high density linkage maps and for association mapping applications due to factors such as low reproducibility, intense labour requirements and/or limited numbers. We report here on the use of a reduced representation library strategy combined with next generation Illumina sequencing for rapid and large scale discovery of SNPs in eight flax genotypes. SNP discovery was performed through in silico analysis of the sequencing data against the whole genome shotgun sequence assembly of flax genotype CDC Bethune. Genotyping-by-sequencing of an F6-derived recombinant inbred line population provided validation of the SNPs.

RESULTS

Reduced representation libraries of eight flax genotypes were sequenced on the Illumina sequencing platform resulting in sequence coverage ranging from 4.33 to 15.64X (genome equivalents). Depending on the relatedness of the genotypes and the number and length of the reads, between 78% and 93% of the reads mapped onto the CDC Bethune whole genome shotgun sequence assembly. A total of 55,465 SNPs were discovered with the largest number of SNPs belonging to the genotypes with the highest mapping coverage percentage. Approximately 84% of the SNPs discovered were identified in a single genotype, 13% were shared between any two genotypes and the remaining 3% in three or more. Nearly a quarter of the SNPs were found in genic regions. A total of 4,706 out of 4,863 SNPs discovered in Macbeth were validated using genotyping-by-sequencing of 96 F6 individuals from a recombinant inbred line population derived from a cross between CDC Bethune and Macbeth, corresponding to a validation rate of 96.8%.

CONCLUSIONS

Next generation sequencing of reduced representation libraries was successfully implemented for genome-wide SNP discovery from flax. The genotyping-by-sequencing approach proved to be efficient for validation. The SNP resources generated in this work will assist in generating high density maps of flax and facilitate QTL discovery, marker-assisted selection, phylogenetic analyses, association mapping and anchoring of the whole genome shotgun sequence.

Single nucleotide polymorphism discovery from wheat next-generation sequence data

Single nucleotide polymorphisms (SNPs) are the most abundant type of molecular genetic marker and can be used for producing high-resolution genetic maps, marker-trait association studies and marker-assisted breeding. Large polyploid genomes such as wheat present a challenge for SNP discovery because of the potential presence of multiple homoeologs for each gene. AutoSNPdb has been successfully applied to identify SNPs from Sanger sequence data for several species, including barley, rice and Brassica, but the volume of data required to accurately call SNPs in the complex genome of wheat has prevented its application to this important crop. DNA sequencing technology has been revolutionized by the introduction of next-generation sequencing, and it is now possible to generate several million sequence reads in a timely and cost-effective manner. We have produced wheat transcriptome sequence data using 454 sequencing technology and applied this for SNP discovery using a modified autoSNPdb method, which integrates SNP and gene annotation information with a graphical viewer. A total of 4,694,141 sequence reads from three bread wheat varieties were assembled to identify a total of 38 928 candidate SNPs. Each SNP is within an assembly complete with annotation, enabling the selection of polymorphism within genes of interest.

Targeted re-sequencing of the allohexaploid wheat exome

Bread wheat, Triticum aestivum, is an allohexaploid composed of the three distinct ancestral genomes, A, B and D. The polyploid nature of the wheat genome together with its large size has limited our ability to generate the significant amount of sequence data required for whole genome studies. Even with the advent of next-generation sequencing technology, it is still relatively expensive to generate whole genome sequences for more than a few wheat genomes at any one time. To overcome this problem, we have developed a targeted-capture re-sequencing protocol based upon NimbleGen array technology to capture and characterize 56.5 Mb of genomic DNA with sequence similarity to over 100 000 transcripts from eight different UK allohexaploid wheat varieties. Using this procedure in conjunction with a carefully designed bioinformatic procedure, we have identified more than 500 000 putative single-nucleotide polymorphisms (SNPs). While 80% of these were variants between the homoeologous genomes, A, B and D, a significant number (20%) were putative varietal SNPs between the eight varieties studied. A small number of these latter polymorphisms were experimentally validated using KASPar technology and 94% proved to be genuine. The procedures described here to sequence a large proportion of the wheat genome, and the various SNPs identified should be of considerable use to the wider wheat community.

The organelle genomes of Hassawi rice (Oryza sativa L.) and its hybrid in saudi arabia: genome variation, rearrangement, and origins

Hassawi rice (Oryza sativa L.) is a landrace adapted to the climate of Saudi Arabia, characterized by its strong resistance to soil salinity and drought. Using high quality sequencing reads extracted from raw data of a whole genome sequencing project, we assembled both chloroplast (cp) and mitochondrial (mt) genomes of the wild-type Hassawi rice (Hassawi-1) and its dwarf hybrid (Hassawi-2). We discovered 16 InDels (insertions and deletions) but no SNP (single nucleotide polymorphism) is present between the two Hassawi cp genomes. We identified 48 InDels and 26 SNPs in the two Hassawi mt genomes and a new type of sequence variation, termed reverse complementary variation (RCV) in the rice cp genomes. There are two and four RCVs identified in Hassawi-1 when compared to 93–11 (indica) and Nipponbare (japonica), respectively. Microsatellite sequence analysis showed there are more SSRs in the genic regions of both cp and mt genomes in the Hassawi rice than in the other rice varieties. There are also large repeats in the Hassawi mt genomes, with the longest length of 96,168 bp and 96,165 bp in Hassawi-1 and Hassawi-2, respectively. We believe that frequent DNA rearrangement in the Hassawi mt and cp genomes indicate ongoing dynamic processes to reach genetic stability under strong environmental pressures. Based on sequence variation analysis and the breeding history, we suggest that both Hassawi-1 and Hassawi-2 originated from the Indonesian variety Peta since genetic diversity between the two Hassawi cultivars is very low albeit an unknown historic origin of the wild-type Hassawi rice.

Reference genome-independent assessment of mutation density using restriction enzyme-phased sequencing

BACKGROUND

The availability of low cost sequencing has spurred its application to discovery and typing of variation, including variation induced by mutagenesis. Mutation discovery is challenging as it requires a substantial amount of sequencing and analysis to detect very rare changes and distinguish them from noise. Also challenging are the cases when the organism of interest has not been sequenced or is highly divergent from the reference.

RESULTS

We describe the development of a simple method for reduced representation sequencing. Input DNA was digested with a single restriction enzyme and ligated to Y adapters modified to contain a sequence barcode and to provide a compatible overhang for ligation. We demonstrated the efficiency of this method at SNP discovery using rice and arabidopsis. To test its suitability for the discovery of very rare SNP, one control and three mutagenized rice individuals (1, 5 and 10 mM sodium azide) were used to prepare genomic libraries for Illumina sequencers by ligating barcoded adapters to NlaIII restriction sites. For genome-dependent discovery 15-30 million of 80 base reads per individual were aligned to the reference sequence achieving individual sequencing coverage from 7 to 15×. We identified high-confidence base changes by comparing sequences across individuals and identified instances consistent with mutations, i.e. changes that were found in a single treated individual and were solely GC to AT transitions. For genome-independent discovery 70-mers were extracted from the sequence of the control individual and single-copy sequence was identified by comparing the 70-mers across samples to evaluate copy number and variation. This de novo "genome" was used to align the reads and identify mutations as above. Covering approximately 1/5 of the 380 Mb genome of rice we detected mutation densities ranging from 0.6 to 4 per Mb of diploid DNA depending on the mutagenic treatment.

CONCLUSIONS

The combination of a simple and cost-effective library construction method, with Illumina sequencing, and the use of a bioinformatic pipeline allows practical SNP discovery regardless of whether a genomic reference is available.

Transcript-specific, single-nucleotide polymorphism discovery and linkage analysis in hexaploid bread wheat (Triticum aestivum L.)

Food security is a global concern and substantial yield increases in cereal crops are required to feed the growing world population. Wheat is one of the three most important crops for human and livestock feed. However, the complexity of the genome coupled with a decline in genetic diversity within modern elite cultivars has hindered the application of marker-assisted selection (MAS) in breeding programmes. A crucial step in the successful application of MAS in breeding programmes is the development of cheap and easy to use molecular markers, such as single-nucleotide polymorphisms. To mine selected elite wheat germplasm for intervarietal single-nucleotide polymorphisms, we have used expressed sequence tags derived from public sequencing programmes and next-generation sequencing of normalized wheat complementary DNA libraries, in combination with a novel sequence alignment and assembly approach. Here, we describe the development and validation of a panel of 1114 single-nucleotide polymorphisms in hexaploid bread wheat using competitive allele-specific polymerase chain reaction genotyping technology. We report the genotyping results of these markers on 23 wheat varieties, selected to represent a broad cross-section of wheat germplasm including a number of elite UK varieties. Finally, we show that, using relatively simple technology, it is possible to rapidly generate a linkage map containing several hundred single-nucleotide polymorphism markers in the doubled haploid mapping population of Avalon × Cadenza.

High-throughput SNP discovery and genotyping in durum wheat (Triticum durum Desf.)

We describe the application of complexity reduction of polymorphic sequences (CRoPS(®)) technology for the discovery of SNP markers in tetraploid durum wheat (Triticum durum Desf.). A next-generation sequencing experiment was carried out on reduced representation libraries obtained from four durum cultivars. SNP validation and minor allele frequency (MAF) estimate were carried out on a panel of 12 cultivars, and the feasibility of genotyping these SNPs in segregating populations was tested using the Illumina Golden Gate (GG) technology. A total of 2,659 SNPs were identified on 1,206 consensus sequences. Among the 768 SNPs that were chosen irrespective of their genomic repetitiveness level and assayed on the Illumina BeadExpress genotyping system, 275 (35.8%) SNPs matched the expected genotypes observed in the SNP discovery phase. MAF data indicated that the overall SNP informativeness was high: a total of 196 (71.3%) SNPs had MAF >0.2, of which 76 (27.6%) showed MAF >0.4. Of these SNPs, 157 were mapped in one of two mapping populations (Meridiano × Claudio and Colosseo × Lloyd) and integrated into a common genetic map. Despite the relatively low genotyping efficiency of the GG assay, the validated CRoPS-derived SNPs showed valuable features for genomics and breeding applications such as a uniform distribution across the wheat genome, a prevailing single-locus codominant nature and a high polymorphism. Here, we report a new set of 275 highly robust genome-wide Triticum SNPs that are readily available for breeding purposes.

Identification of SNP and SSR markers in eggplant using RAD tag sequencing

Background

The eggplant (Solanum melongena L.) genome is relatively unexplored, especially compared to those of the other major Solanaceae crops tomato and potato. In particular, no SNP markers are publicly available; on the other hand, over 1,000 SSR markers were developed and publicly available. We have combined the recently developed Restriction-site Associated DNA (RAD) approach with Illumina DNA sequencing for rapid and mass discovery of both SNP and SSR markers for eggplant.

Results

RAD tags were generated from the genomic DNA of a pair of eggplant mapping parents, and sequenced to produce ~17.5 Mb of sequences arrangeable into ~78,000 contigs. The resulting non-redundant genomic sequence dataset consisted of ~45,000 sequences, of which ~29% were putative coding sequences and ~70% were in common between the mapping parents. The shared sequences allowed the discovery of ~10,000 SNPs and nearly 1,000 indels, equivalent to a SNP frequency of 0.8 per Kb and an indel frequency of 0.07 per Kb. Over 2,000 of the SNPs are likely to be mappable via the Illumina GoldenGate assay. A subset of 384 SNPs was used to successfully fingerprint a panel of eggplant germplasm, producing a set of informative diversity data. The RAD sequences also included nearly 2,000 putative SSRs, and primer pairs were designed to amplify 1,155 loci.

Conclusion

The high throughput sequencing of the RAD tags allowed the discovery of a large number of DNA markers, which will prove useful for extending our current knowledge of the genome organization of eggplant, for assisting in marker-aided selection and for carrying out comparative genomic analyses within the Solanaceae family.