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

Wheat improvement through advances in single nucleotide polymorphism (SNP) detection and genotyping with a special emphasis on rust resistance

KEY MESSAGE

Single nucleotide polymorphism (SNP) markers in wheat and their prospects in breeding with special reference to rust resistance. Single nucleotide polymorphism (SNP)-based markers are increasingly gaining momentum for screening and utilizing vital agronomic traits in wheat. To date, more than 260 million SNPs have been detected in modern cultivars and landraces of wheat. This rapid SNP discovery was made possible through the release of near-complete reference and pan-genome assemblies of wheat and its wild relatives, coupled with whole genome sequencing (WGS) of thousands of wheat accessions. Further, genotyping customized SNP sites were facilitated by a series of arrays (9 to 820Ks), a cost effective substitute WGS. Lately, germplasm-specific SNP arrays have been introduced to characterize novel traits and detect closely linked SNPs for marker-assisted breeding. Subsequently, the kompetitive allele-specific PCR (KASP) assay was introduced for rapid and large-scale screening of specific SNP markers. Moreover, with the advances and reduction in sequencing costs, ample opportunities arise for generating SNPs artificially through mutations and in combination with next-generation sequencing and comparative genomic analyses. In this review, we provide historical developments and prospects of SNP markers in wheat breeding with special reference to rust resistance where over 50 genetic loci have been characterized through SNP markers. Rust resistance is one of the most essential traits for wheat breeding as new strains of the Puccinia fungus, responsible for rust diseases, evolve frequently and globally.

Wheat Water-Soluble Carbohydrate Remobilisation under Water Deficit by 1-FEH w3

Fructan 1-exohydrolase (1-FEH) is one of the major enzymes in water-soluble carbohydrate (WSC) remobilisation for grains in wheat. We investigated the functional role of 1-FEH w1, w2, and w3 isoforms in WSC remobilisation under post-anthesis water deficit using mutation lines derived from the Australian wheat variety Chara. F1 seeds, developed by backcrossing the 1-FEH w1, w2, and w3 mutation lines with Chara, were genotyped using the Infinium 90K SNP iSelect platform to characterise the mutated region. Putative deletions were identified in FEH mutation lines encompassing the FEH genomic regions. Mapping analysis demonstrated that mutations affected significantly longer regions than the target FEH gene regions. Functional roles of the non-target genes were carried out utilising bioinformatics and confirmed that the non-target genes were unlikely to confound the effects considered to be due to the influence of 1-FEH gene functions. Glasshouse experiments revealed that the 1-FEH w3 mutation line had a slower degradation and remobilisation of fructans than the 1-FEH w2 and w1 mutation lines and Chara, which reduced grain filling and grain yield. Thus, 1-FEH w3 plays a vital role in reducing yield loss under drought. This insight into the distinct role of the 1-FEH isoforms provides new gene targets for water-deficit-tolerant wheat breeding.

Genetic approaches to exploit landraces for improvement of Triticum turgidum ssp. durum in the age of climate change

Addressing the challenges of climate change and durum wheat production is becoming an important driver for food and nutrition security in the Mediterranean area, where are located the major producing countries (Italy, Spain, France, Greece, Morocco, Algeria, Tunisia, Turkey, and Syria). One of the emergent strategies, to cope with durum wheat adaptation, is the exploration and exploitation of the existing genetic variability in landrace populations. In this context, this review aims to highlight the important role of durum wheat landraces as a useful genetic resource to improve the sustainability of Mediterranean agroecosystems, with a focus on adaptation to environmental stresses. We described the most recent molecular techniques and statistical approaches suitable for the identification of beneficial genes/alleles related to the most important traits in landraces and the development of molecular markers for marker-assisted selection. Finally, we outline the state of the art about landraces genetic diversity and signature of selection, already identified from these accessions, for adaptability to the environment.

Significance of the Diversification of Wheat Species for the Assembly and Functioning of the Root-Associated Microbiome

Wheat, one of the major crops in the world, has had a complex history that includes genomic hybridizations between Triticum and Aegilops species and several domestication events, which resulted in various wild and domesticated species (especially Triticum aestivum and Triticum durum), many of them still existing today. The large body of information available on wheat-microbe interactions, however, was mostly obtained without considering the importance of wheat evolutionary history and its consequences for wheat microbial ecology. This review addresses our current understanding of the microbiome of wheat root and rhizosphere in light of the information available on pre- and post-domestication wheat history, including differences between wild and domesticated wheats, ancient and modern types of cultivars as well as individual cultivars within a given wheat species. This analysis highlighted two major trends. First, most data deal with the taxonomic diversity rather than the microbial functioning of root-associated wheat microbiota, with so far a bias toward bacteria and mycorrhizal fungi that will progressively attenuate thanks to the inclusion of markers encompassing other micro-eukaryotes and archaea. Second, the comparison of wheat genotypes has mostly focused on the comparison of T. aestivum cultivars, sometimes with little consideration for their particular genetic and physiological traits. It is expected that the development of current sequencing technologies will enable to revisit the diversity of the wheat microbiome. This will provide a renewed opportunity to better understand the significance of wheat evolutionary history, and also to obtain the baseline information needed to develop microbiome-based breeding strategies for sustainable wheat farming.

Development of GBTS and KASP Panels for Genetic Diversity, Population Structure, and Fingerprinting of a Large Collection of Broccoli (Brassica oleracea L. var. italica) in China

Broccoli (Brassica oleracea var. italica) is one of the most important and nutritious vegetables widely cultivated in China. In the recent four decades, several improved varieties were bred and developed by Chinese breeders. However, the efforts for improvement of broccoli are hindered by limited information of genetic diversity and genetic relatedness contained within the available germplasms. This study evaluated the genetic diversity, genetic relationship, population structure, and fingerprinting of 372 accessions of broccoli representing most of the variability of broccoli in China. Millions of SNPs were identified by whole-genome sequencing of 23 representative broccoli genotypes. Through several stringent selection criteria, a total of 1,167 SNPs were selected to characterize genetic diversity and population structure. Of these markers, 1,067 SNPs were genotyped by target sequencing (GBTS), and 100 SNPs were genotyped by kompetitive allele specific PCR (KASP) assay. The average polymorphism information content (PIC) and expected heterozygosity (gene diversity) values were 0.33 and 0.42, respectively. Diversity analysis revealed the prevalence of low to moderate genetic diversity in the broccoli accessions indicating a narrow genetic base. Phylogenetic and principal component analyses revealed that the 372 accessions could be clustered into two main groups but with weak groupings. STRUCTURE analysis also suggested the presence of two subpopulations with weak genetic structure. Analysis of molecular variance (AMOVA) identified 13% variance among populations and 87% within populations revealing very low population differentiation, which could be attributed to massive gene flow and the reproductive biology of the crop. Based on high resolving power, a set of 28 KASP markers was chosen for DNA fingerprinting of the broccoli accessions for seed authentication and varietal identification. To the best of our knowledge, this is the first comprehensive study to measure diversity and population structure of a large collection of broccoli in China and also the first application of GBTS and KASP techniques in genetic characterization of broccoli. This work broadens the understanding of diversity, phylogeny, and population structure of a large collection of broccoli, which may enhance future breeding efforts to achieve higher productivity.

High-throughput sequencing techniques to flax genetics and breeding

Flax (Linum usitatissimum L.) is an important oil and fiber crop. Using modern methods for flax breeding allows accelerating the introduction of some desired genes into the genotypes of future varieties. Today, an important condition for their creation is the development of research, that is based on next-generation sequencing (NGS). This review summarizes the results obtained using NGS in flax research. To date, a linkage map with a high marker density has been obtained for L. usitatissimum, which is already being used for a more efficient search for quantitative traits loci. Comparative studies of transcriptomes and miRNomes of flax under stress and in control conditions elucidated molecular-genetic mechanisms of abiotic and biotic stress responses. The very accurate model for genomic selection of flax resistant to pasmo was constructed. Based on NGS-sequencing also some details of the genus Linum evolution were clarified. The knowledge systematized in the review can be useful for researchers working in flax breeding and whereas fundamental interest for understanding the phylogenetic relationships within the genus Linum, the ontogenesis, and the mechanisms of the response of flax plants to various stress factors.

Four Consistent Loci Confer Adult Plant Resistance to Leaf Rust in the Durum Wheat Lines Heller#1 and Dunkler

The durum wheat lines Heller#1 and Dunkler from the International Maize and Wheat Improvement Center Global Wheat Program showed moderate and stable adult plant resistance to leaf rust under high disease pressure over field environments in northwestern Mexico. Leaf rust phenotyping was performed on two recombinant inbred line (RIL) populations derived from crosses of Heller#1 and Dunkler with the susceptible parent Atred#2, conducted under artificially induced Puccinia triticina epidemics in 2013, 2014, 2015, and 2016. The Atred#2 × Heller#1 and Atred#2 × Dunkler populations were genotyped by single nucleotide polymorphism (SNP) platforms and diversity arrays technology markers, respectively. Four leaf rust resistance quantitative trait loci were detected simultaneously in the two RIL populations: Lr46, QLr.cim-2BC, QLr.cim-5BL, and QLr.cim-6BL based on phenotypic data across all four crop seasons. They explained 11.7 to 46.8%, 7.2 to 26.1%, 8.4 to 24.1%, and 12.4 to 28.5%, respectively, of the phenotypic variation for leaf rust resistance in Atred#2 × Heller#1 and 16.3 to 56.6%, 6.7 to 15.7%, 4.1 to 10.1%, and 5.1 to 20.2% of the variation in the Atred#2 × Dunkler population. Only the resistance allele of QLr.cim-2BC was from the susceptible parent Atred#2, and resistance alleles at other loci came from the resistant parents Heller#1 and Dunkler. The SNP markers closely linked to Lr46 and QLr.cim-2BC were converted to kompetitive allele specific PCR markers for use in marker-assisted selection to improve leaf rust resistance through crosses with Heller#1 and Dunkler sources.

Genetic diversity and linkage disequilibrium using SNP (KASP) and AFLP markers in a worldwide durum wheat (Triticum turgidum L. var durum) collection

The aim of this work was to analyze the genetic diversity and linkage disequilibrium in a collection of 168 durum wheat accessions (Triticum turgidum L. var. durum) of different origins. Our collection was mainly composed of released and unreleased Argentinian germplasm, with additional genotypes from Italy, Chile, France, CIMMYT, Cyprus, USA and WANA region. To this end, the entire collection was characterized with 85 Single Nucleotide Polymorphism (SNP) markers obtained by Kompetitive Allele Specific PCR (KASP), giving a heterozygosity (He) mean value of 0.183 and a coefficient of genetic differentiation (Gst) value of 0.139. A subset of 119 accessions was characterized with six Amplified Fragment Length Polymorphism (AFLP) primer combinations. A total of 181 polymorphic markers (125 AFLP and 56 SNP) amplified across this subset revealed He measures of 0.352 and 0.182, respectively. Of these, 134 were selected to estimate the genome-wide linkage disequilibrium obtaining low significant values (r2 = 0.11) in the subset, indicating its suitability for future genome-wide association studies (GWAS). The structure analysis conducted in the entire collection with SNP detected two subpopulations. However, the structure analysis conducted with AFLP markers in the subset of 119 accessions proved to have greater degree of resolution and detect six subpopulations. The information provided by both marker types was complementary and showed a strong association between old Argentinian and Italian germplasm and a contribution of CIMMYT germplasm to modern Argentinian, Chilean and Cypriot accessions. The influence of Mediterranean germplasm, mainly from Italy, on part of the modern Argentinian cultivars or breeding lines was also clearly evidenced. Although our analysis yields conclusive results and useful information for association mapping studies, further analyses are needed to refine the number of subpopulations present in the germplasm collection analyzed.

High-Throughput Genotype, Morphology, and Quality Traits Evaluation for the Assessment of Genetic Diversity of Wheat Landraces from Sicily

During the XX Century, the widespread use of modern wheat cultivars drastically reduced the cultivation of ancient landraces, which nowadays are confined to niche cultivation areas. Several durum wheat landraces adapted to the extreme environments of the Mediterranean region, are still being cultivated in Sicily, Italy. Detailed knowledge of the genetic diversity of this germplasm could lay the basis for their efficient management in breeding programs, for a wide-range range of traits. The aim of the present study was to characterize a collection of durum wheat landraces from Sicily, using single nucleotide polymorphisms (SNP) markers, together with agro-morphological, phenological and quality-related traits. Two modern cv. Simeto, Claudio, and the hexaploid landrace, Cuccitta, were used as outgroups. Cluster analysis and Principal Coordinates Analysis (PCoA) allowed us to identify four main clusters across the analyzed germplasm, among which a cluster included only historical and modern varieties. Likewise, structure analysis was able to distinguish the ancient varieties from the others, grouping the entries in seven cryptic genetic clusters. Furthermore, a Principal Component Analysis (PCA) was able to separate the modern testers from the ancient germplasm. This approach was useful to classify and evaluate Sicilian ancient wheat germplasm, supporting their safeguard and providing a genetic fingerprint that is necessary for avoiding commercial frauds to sustaining the economic profits of farmers resorting to landraces cultivation.

High density genetic mapping of Fusarium head blight resistance QTL in tetraploid wheat

Breeding for Fusarium head blight (FHB) resistance in durum wheat is complicated by the quantitative trait expression and narrow genetic diversity of available resources. High-density mapping of the FHB resistance quantitative trait loci (QTL), evaluation of their co-localization with plant height and maturity QTL and the interaction among the identified QTL are the objectives of this study. Two doubled haploid (DH) populations, one developed from crosses between Triticum turgidum ssp. durum lines DT707 and DT696 and the other between T. turgidum ssp. durum cv. Strongfield and T. turgidum ssp. carthlicum cv. Blackbird were genotyped using the 90K Infinium iSelect chip and evaluated phenotypically at multiple field FHB nurseries over years. A moderate broad-sense heritability indicated a genotype-by-environment interaction for the expression of FHB resistance in both populations. Resistance QTL were identified for the DT707 × DT696 population on chromosomes 1B, 2B, 5A (two loci) and 7A and for the Strongfield × Blackbird population on chromosomes 1A, 2A, 2B, 3A, 6A, 6B and 7B with the QTL on chromosome 1A and those on chromosome 5A being more consistently expressed over environments. FHB resistance co-located with plant height and maturity QTL on chromosome 5A and with a maturity QTL on chromosome 7A for the DT707 × DT696 population. Resistance also co-located with plant height QTL on chromosomes 2A and 3A and with maturity QTL on chromosomes 1A and 7B for the Strongfield × Blackbird population. Additive × additive interactions were identified, for example between the two FHB resistance QTL on chromosome 5A for the DT707 × DT696 population and the FHB resistance QTL on chromosomes 1A and 7B for the Strongfield × Blackbird population. Application of the Single Nucleotide Polymorphic (SNP) markers associated with FHB resistance QTL identified in this study will accelerate combining genes from the two populations.

High-density genetic mapping of a major QTL for resistance to multiple races of loose smut in a tetraploid wheat cross

Loose smut, caused by Ustilago tritici (Pers.) Rostr., is a systemic disease of tetraploid durum wheat (Triticum turgidum L.). Loose smut can be economically controlled by growing resistant varieties, making it important to find and deploy new sources of resistance. Blackbird, a variety of T. turgidum L. subsp. carthlicum (Nevski) A. Love & D. Love, carries a high level of resistance to loose smut. Blackbird was crossed with the loose smut susceptible durum cultivar Strongfield to produce a doubled haploid (DH) mapping population. The parents and progenies were inoculated with U. tritici races T26, T32 and T33 individually and as a mixture at Swift Current, Canada in 2011 and 2012 and loose smut incidence (LSI) was assessed. Genotyping of the DH population and parents using an Infinium iSelect 90K single nucleotide polymorphism (SNP) array identified 12,952 polymorphic SNPs. The SNPs and 426 SSRs (previously genotyped in the same population) were mapped to 16 linkage groups spanning 3008.4 cM at an average inter-marker space of 0.2 cM in a high-density genetic map. Composite interval mapping analysis revealed three significant quantitative trait loci (QTL) for loose smut resistance on chromosomes 3A, 6B and 7A. The loose smut resistance QTL on 6B (QUt.spa-6B.2) and 7A (QUt.spa-7A.2) were derived from Blackbird. Strongfield contributed the minor QTL on 3A (QUt.spa-3A.2). The resistance on 6B was a stable major QTL effective against all individual races and the mixture of the three races; it explained up to 74% of the phenotypic variation. This study is the first attempt in durum wheat to identify and map loose smut resistance QTL using a high-density genetic map. The QTL QUt.spa-6B.2 would be an effective source for breeding resistance to multiple races of the loose smut pathogen because it provides near-complete broad resistance to the predominant virulence on the Canadian prairies.

Toward Identification of Black Lemma and Pericarp Gene Blp1 in Barley Combining Bulked Segregant Analysis and Specific-Locus Amplified Fragment Sequencing

Black barley is caused by phytomelanin synthesized in lemma and/or pericarp and the trait is controlled by one dominant gene Blp1. The gene is mapped on chromosome 1H by molecular markers, but it is yet to be isolated. Specific-locus amplified fragment sequencing (SLAF-seq) is an effective method for large-scale de novo single nucleotide polymorphism (SNP) discovery and genotyping. In the present study, SLAF-seq with bulked segregant analysis (BSA) was employed to obtain sufficient markers to fine mapping Blp1 gene in an F2 population derived from Hatiexi No.1 × Zhe5819. Based on SNP screening criteria, a total of 77,542 polymorphic SNPs met the requirements for association analysis. Combining two association analysis methods, the overlapped region with a size of 32.41 Mb on chromosome 1H was obtained as the candidate region of Blp1 gene. According to SLAF-seq data, markers were developed in the target region and were used for mapping the Blp1 gene. Linkage analysis showed that Blp1 co-segregated with HZSNP34 and HZSNP36, and was delimited by two markers (HZSNP35 and HZSNP39) spanning 8.1 cM in 172 homozygous yellow grain F2 plants of Hatiexi No.1 × Zhe5819. More polymorphic markers were screened in the reduced target region and were used to genotype the population. As a result, Blp1 was delimited within a 1.66 Mb on chromosome 1H by the upstream marker HZSNP63 and the downstream marker HZSNP59. Our results demonstrated the utility of SLAF-seq-BSA approach to identify the candidate region and discover polymorphic markers at the specific targeted genomic region.

High-throughput SNP genotyping of modern and wild emmer wheat for yield and root morphology using a combined association and linkage analysis

Durum wheat (Triticum turgidum var. durum Desf.) is a major world crop that is grown primarily in areas of the world that experience periodic drought, and therefore, breeding climate-resilient durum wheat is a priority. High-throughput single nucleotide polymorphism (SNP) genotyping techniques have greatly increased the power of linkage and association mapping analyses for bread wheat, but as yet there is no durum wheat-specific platform available. In this study, we evaluate the new 384HT Wheat Breeders Array for its usefulness in tetraploid wheat breeding by genotyping a breeding population of F₆ hybrids, derived from multiple crosses between T. durum cultivars and wild and cultivated emmer wheat accessions. Using a combined linkage and association mapping approach, we generated a genetic map including 1345 SNP markers, and identified markers linked to 6 QTLs for coleoptile length (2), heading date (1), anthocyanin accumulation (1) and osmotic stress tolerance (2). We also developed a straightforward approach for combining genetic data from multiple families of limited size that will be useful for evaluating and mapping pre-existing breeding material.

Genetic analysis and mapping of adult plant resistance loci to leaf rust in durum wheat cultivar Bairds

New leaf rust adult plant resistance (APR) QTL QLr.cim - 6BL was mapped and confirmed the known pleotropic APR gene Lr46 effect on leaf rust in durum wheat line Bairds. CIMMYT-derived durum wheat line Bairds displays an adequate level of adult plant resistance (APR) to leaf rust in Mexican field environments. A recombinant inbred line (RIL) population developed from a cross of Bairds with susceptible parent Atred#1 was phenotyped for leaf rust response at Ciudad Obregon, Mexico, during 2013, 2014, 2015 and 2016 under artificially created epidemics of Puccinia triticina (Pt) race BBG/BP. The RIL population and its parents were genotyped with the 50 K diversity arrays technology (DArT) sequence system and simple sequence repeat (SSR) markers. A genetic map comprising 1150 markers was used to map the resistance loci. Four significant quantitative trait loci (QTLs) were detected on chromosomes 1BL, 2BC (centromere region), 5BL and 6BL. These QTLs, named Lr46, QLr.cim-2BC, QLr.cim-5BL and QLr.cim-6BL, respectively, explained 13.5-60.8%, 9.0-14.3%, 2.8-13.9%, and 11.6-29.4%, respectively, of leaf rust severity variation by the inclusive composite interval mapping method. All of these resistance loci were contributed by the resistant parent Bairds, except for QLr.cim-2BC, which came from susceptible parent Atred#1. Among these, the QTL on chromosome 1BL was the known pleiotropic APR gene Lr46, whereas QLr.cim-6BL, a consistently detected locus, should be a new leaf rust resistance locus in durum wheat. The mean leaf rust severity of RILs carrying all four QTLs ranged from 8.0 to 17.5%, whereas it ranged from 10.9 to 38.5% for three QTLs (Lr46 + 5BL + 6BL) derived from the resistant parent Bairds. Two RILs with four QTLs combinations can be used as sources of complex APR in durum wheat breeding.

SNP-based association analysis for seedling traits in durum wheat (Triticum turgidum L. durum (Desf.))

In the present study, 150 accessions of worldwide originated durum wheat germplasm (Triticum turgidum spp. durum) were observed for major seedling traits and their growth. The accessions were evaluated for major seedling traits under controlled conditions of hydroponics at the 13th, 20th, 27th and 34th day-after germination. Biomass traits were measured at the 34th day-after germination. Correlation analysis was conducted among the seedling traits and three field traits at maturity, plant height, grain weight and 1000-grain weight observed in four consecutive years. Associations of the measured seedling traits and SNP markers were analyzed based on the mixed linear model (MLM). The results indicated that highly significant genetic variation and robust heritability were found for the seedling and field mature traits. In total, 259 significant associations were detected for all the traits and four growth stages. The phenotypic variation explained (R2) by a single SNP marker is higher than 10% for most (84%) of the significant SNP markers. Forty-six SNP markers associated with multiple traits, indicating non-neglectable pleiotropy in seedling stage. The associated SNP markers could be helpful for genetic analysis of seedling traits, and marker-assisted breeding of new wheat varieties with strong seedling vigor.

Combined Genomic and Genetic Data Integration of Major Agronomical Traits in Bread Wheat (Triticum aestivum L.)

The high resolution integration of bread wheat genetic and genomic resources accumulated during the last decades offers the opportunity to unveil candidate genes driving major agronomical traits to an unprecedented scale. We combined 27 public quantitative genetic studies and four genetic maps to deliver an exhaustive consensus map consisting of 140,315 molecular markers hosting 221, 73, and 82 Quantitative Trait Loci (QTL) for respectively yield, baking quality, and grain protein content (GPC) related traits. Projection of the consensus genetic map and associated QTLs onto the wheat syntenome made of 99,386 genes ordered on the 21 chromosomes delivered a complete and non-redundant repertoire of 18, 8, 6 metaQTLs for respectively yield, baking quality and GPC, altogether associated to 15,772 genes (delivering 28,630 SNP-based makers) including 37 major candidates. Overall, this study illustrates a translational research approach in transferring information gained from grass relatives to dissect the genomic regions hosting major loci governing key agronomical traits in bread wheat, their flanking markers and associated candidate genes to be now considered as a key resource for breeding programs.

Development of a Genetic Map for Onion (Allium cepa L.) Using Reference-Free Genotyping-by-Sequencing and SNP Assays

Single nucleotide polymorphisms (SNPs) play important roles as molecular markers in plant genomics and breeding studies. Although onion (Allium cepa L.) is an important crop globally, relatively few molecular marker resources have been reported due to its large genome and high heterozygosity. Genotyping-by-sequencing (GBS) offers a greater degree of complexity reduction followed by concurrent SNP discovery and genotyping for species with complex genomes. In this study, GBS was employed for SNP mining in onion, which currently lacks a reference genome. A segregating F₂ population, derived from a cross between 'NW-001' and 'NW-002,' as well as multiple parental lines were used for GBS analysis. A total of 56.15 Gbp of raw sequence data were generated and 1,851,428 SNPs were identified from the de novo assembled contigs. Stringent filtering resulted in 10,091 high-fidelity SNP markers. Robust SNPs that satisfied the segregation ratio criteria and with even distribution in the mapping population were used to construct an onion genetic map. The final map contained eight linkage groups and spanned a genetic length of 1,383 centiMorgans (cM), with an average marker interval of 8.08 cM. These robust SNPs were further analyzed using the high-throughput Fluidigm platform for marker validation. This is the first study in onion to develop genome-wide SNPs using GBS. The resulting SNP markers and developed linkage map will be valuable tools for genetic mapping of important agronomic traits and marker-assisted selection in onion breeding programs.

Genotyping by Sequencing Using Specific Allelic Capture to Build a High-Density Genetic Map of Durum Wheat

Targeted sequence capture is a promising technology which helps reduce costs for sequencing and genotyping numerous genomic regions in large sets of individuals. Bait sequences are designed to capture specific alleles previously discovered in parents or reference populations. We studied a set of 135 RILs originating from a cross between an emmer cultivar (Dic2) and a recent durum elite cultivar (Silur). Six thousand sequence baits were designed to target Dic2 vs. Silur polymorphisms discovered in a previous RNAseq study. These baits were exposed to genomic DNA of the RIL population. Eighty percent of the targeted SNPs were recovered, 65% of which were of high quality and coverage. The final high density genetic map consisted of more than 3,000 markers, whose genetic and physical mapping were consistent with those obtained with large arrays.

Genome-wide linkage mapping of flour color-related traits and polyphenol oxidase activity in common wheat

KEY MESSAGE

Fifty-six QTL for flour color-related traits and polyphenol oxidase activity were identified using a genome-wide linkage mapping of data from a RIL population derived from a Gaocheng 8901/Zhoumai 16 cross.

ABSTRACT

Flour color-related traits, including L*, a*, b*, yellow pigment content (YPC), and polyphenol oxidase (PPO) activity are important parameters influencing the quality of wheat end-use products. Mapping quantitative trait loci (QTL) for these traits and characterization of candidate genes are important for improving wheat quality. The aims of this study were to identify QTL for flour color-related traits and PPO activity and to characterize candidate genes using a high-density genetic linkage map in a common wheat recombinant inbred line (RIL) population derived from a cross between Gaocheng 8901 and Zhoumai 16. A linkage map was constructed by genotyping the RILs with the wheat 90 K iSelect array. Fifty-six QTL were mapped on 35 chromosome regions on homoeologous groups 1, 2, 5 and 7 chromosomes, and chromosomes 3B, 4A, 4B and 6B. Four QTL were for PPO activity, and the others were for flour color-related traits. Compared with previous studies, five QTL for a*, two for b*, one for L*, one for YPC and one for PPO activity were new. The new QTL on chromosome 2DL was involved in both a* and YPC, and another on chromosome 7DS affected both a* and L*. The scan for SNP sequences tightly linked to QTL for flour color-related traits against the wheat and/or related cereals genomes identified six candidate genes significantly related to these traits, and five of them were associated with the terpenoid backbone biosynthesis pathway. The high-density genetic linkage map of Gaocheng 8901/Zhoumai 16 represents a useful tool to identify QTL for important quality traits and candidate genes.

TaXA21-A1 on chromosome 5AL is associated with resistance to multiple pests in wheat

KEY MESSAGE

The wheat ortholog of the rice gene OsXA21 against bacterial leaf blight showed resistance to multiple pests in bread wheat but different interacting proteins.

ABSTRACT

A quantitative trait locus QYr.osu-5A on the long arm of chromosome 5A in bread wheat (Triticum aestivum L., 2n = 6x = 42; AABBDD) was previously reported to confer consistent resistance in adult plants to predominant stripe rust races, but the gene causing the quantitative trait locus (QTL) is not known. Single-nucleotide polymorphism (SNP) markers were used to saturate the QTL region. Comparative and syntenic regions between wheat and rice (Oryza sativa) were applied to identify candidate genes for QYr.osu-5A. TaXA21-A1, which is referred to as a wheat ortholog of OsXA21-like gene on chromosome 9 in rice, was mapped under the peak of the QYr.osu-5A. TaXA21-A1 not only explained the phenotypic variation in reaction to different stripe rust races but also showed significant effects on resistance to powdery mildew and Hessian fly biotype BP. The natural allelic variation resulted in the alternations of four amino acids in deduced TaXA21-A1 proteins. The interacting proteins of TaXA21-A1 were different from those identified by OsXA21 on rice chromosome 11 against bacterial leaf blight. TaXA21-A1 confers unique resistance against multiple pests in wheat but might not have common protein interactors or thus overlapping functions with OsXA21 in rice. XA21 function has diverged during evolution of cereal crops. The molecular marker developed for TaXA21-A1 would accelerate its application of the candidate gene at the QYr.osu-5A locus in wheat breeding 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.

Genetic Map of Triticale Integrating Microsatellite, DArT and SNP Markers

Triticale (×Triticosecale Wittm) is an economically important crop for fodder and biomass production. To facilitate the identification of markers for agronomically important traits and for genetic and genomic characteristics of this species, a new high-density genetic linkage map of triticale was constructed using doubled haploid (DH) population derived from a cross between cultivars 'Hewo' and 'Magnat'. The map consists of 1615 bin markers, that represent 50 simple sequence repeat (SSR), 842 diversity array technology (DArT), and 16888 DArTseq markers mapped onto 20 linkage groups assigned to the A, B, and R genomes of triticale. No markers specific to chromosome 7R were found, instead mosaic linkage group composed of 1880 highly distorted markers (116 bins) from 10 wheat chromosomes was identified. The genetic map covers 4907 cM with a mean distance between two bins of 3.0 cM. Comparative analysis in respect to published maps of wheat, rye and triticale revealed possible deletions in chromosomes 4B, 5A, and 6A, as well as inversion in chromosome 7B. The number of bin markers in each chromosome varied from 24 in chromosome 3R to 147 in chromosome 6R. The length of individual chromosomes ranged between 50.7 cM for chromosome 2R and 386.2 cM for chromosome 7B. A total of 512 (31.7%) bin markers showed significant (P < 0.05) segregation distortion across all chromosomes. The number of 8 the segregation distorted regions (SDRs) were identified on 1A, 7A, 1B, 2B, 7B (2 SDRs), 5R and 6R chromosomes. The high-density genetic map of triticale will facilitate fine mapping of quantitative trait loci, the identification of candidate genes and map-based cloning.

Association of Agronomic Traits with SNP Markers in Durum Wheat (Triticum turgidum L. durum (Desf.))

Association mapping is a powerful approach to detect associations between traits of interest and genetic markers based on linkage disequilibrium (LD) in molecular plant breeding. In this study, 150 accessions of worldwide originated durum wheat germplasm (Triticum turgidum spp. durum) were genotyped using 1,366 SNP markers. The extent of LD on each chromosome was evaluated. Association of single nucleotide polymorphisms (SNP) markers with ten agronomic traits measured in four consecutive years was analyzed under a mix linear model (MLM). Two hundred and one significant association pairs were detected in the four years. Several markers were associated with one trait, and also some markers were associated with multiple traits. Some of the associated markers were in agreement with previous quantitative trait loci (QTL) analyses. The function and homology analyses of the corresponding ESTs of some SNP markers could explain many of the associations for plant height, length of main spike, number of spikelets on main spike, grain number per plant, and 1000-grain weight, etc. The SNP associations for the observed traits are generally clustered in specific chromosome regions of the wheat genome, mainly in 2A, 5A, 6A, 7A, 1B, and 6B chromosomes. This study demonstrates that association mapping can complement and enhance previous QTL analyses and provide additional information for marker-assisted selection.

A high-density, SNP-based consensus map of tetraploid wheat as a bridge to integrate durum and bread wheat genomics and breeding

Consensus linkage maps are important tools in crop genomics. We have assembled a high-density tetraploid wheat consensus map by integrating 13 data sets from independent biparental populations involving durum wheat cultivars (Triticum turgidum ssp. durum), cultivated emmer (T. turgidum ssp. dicoccum) and their ancestor (wild emmer, T. turgidum ssp. dicoccoides). The consensus map harboured 30 144 markers (including 26 626 SNPs and 791 SSRs) half of which were present in at least two component maps. The final map spanned 2631 cM of all 14 durum wheat chromosomes and, differently from the individual component maps, all markers fell within the 14 linkage groups. Marker density per genetic distance unit peaked at centromeric regions, likely due to a combination of low recombination rate in the centromeric regions and even gene distribution along the chromosomes. Comparisons with bread wheat indicated fewer regions with recombination suppression, making this consensus map valuable for mapping in the A and B genomes of both durum and bread wheat. Sequence similarity analysis allowed us to relate mapped gene-derived SNPs to chromosome-specific transcripts. Dense patterns of homeologous relationships have been established between the A- and B-genome maps and between nonsyntenic homeologous chromosome regions as well, the latter tracing to ancient translocation events. The gene-based homeologous relationships are valuable to infer the map location of homeologs of target loci/QTLs. Because most SNP and SSR markers were previously mapped in bread wheat, this consensus map will facilitate a more effective integration and exploitation of genes and QTL for wheat breeding purposes.

Genetic diversity in tef [Eragrostis tef (Zucc.) Trotter]

Tef [Eragrostis tef (Zucc.) Trotter] is a cereal crop resilient to adverse climatic and soil conditions, and possessing desirable storage properties. Although tef provides high quality food and grows under marginal conditions unsuitable for other cereals, it is considered to be an orphan crop because it has benefited little from genetic improvement. Hence, unlike other cereals such as maize and wheat, the productivity of tef is extremely low. In spite of the low productivity, tef is widely cultivated by over six million small-scale farmers in Ethiopia where it is annually grown on more than three million hectares of land, accounting for over 30% of the total cereal acreage. Tef, a tetraploid with 40 chromosomes (2n = 4x = 40), belongs to the family Poaceae and, together with finger millet (Eleusine coracana Gaerth.), to the subfamily Chloridoideae. It was originated and domesticated in Ethiopia. There are about 350 Eragrostis species of which E. tef is the only species cultivated for human consumption. At the present time, the gene bank in Ethiopia holds over five thousand tef accessions collected from geographical regions diverse in terms of climate and elevation. These germplasm accessions appear to have huge variability with regard to key agronomic and nutritional traits. In order to properly utilize the variability in developing new tef cultivars, various techniques have been implemented to catalog the extent and unravel the patterns of genetic diversity. In this review, we show some recent initiatives investigating the diversity of tef using genomics, transcriptomics and proteomics and discuss the prospect of these efforts in providing molecular resources that can aid modern tef breeding.

Development and use of molecular markers: past and present

Molecular markers, due to their stability, cost-effectiveness and ease of use provide an immensely popular tool for a variety of applications including genome mapping, gene tagging, genetic diversity diversity, phylogenetic analysis and forensic investigations. In the last three decades, a number of molecular marker techniques have been developed and exploited worldwide in different systems. However, only a handful of these techniques, namely RFLPs, RAPDs, AFLPs, ISSRs, SSRs and SNPs have received global acceptance. A recent revolution in DNA sequencing techniques has taken the discovery and application of molecular markers to high-throughput and ultrahigh-throughput levels. Although, the choice of marker will obviously depend on the targeted use, microsatellites, SNPs and genotyping by sequencing (GBS) largely fulfill most of the user requirements. Further, modern transcriptomic and functional markers will lead the ventures onto high-density genetic map construction, identification of QTLs, breeding and conservation strategies in times to come in combination with other high throughput techniques. This review presents an overview of different marker technologies and their variants with a comparative account of their characteristic features and applications.

Analysis of agronomic and domestication traits in a durum × cultivated emmer wheat population using a high-density single nucleotide polymorphism-based linkage map

Development of a high-density SNP map and evaluation of QTL shed light on domestication events in tetraploid wheat and the potential utility of cultivated emmer wheat for durum wheat improvement. Cultivated emmer wheat (Triticum turgidum ssp. dicoccum) is tetraploid and considered as one of the eight founder crops that spawned the Agricultural Revolution about 10,000 years ago. Cultivated emmer has non-free-threshing seed and a somewhat fragile rachis, but mutations in genes governing these and other agronomic traits occurred that led to the formation of today's fully domesticated durum wheat (T. turgidum ssp. durum). Here, we evaluated a population of recombinant inbred lines (RILs) derived from a cross between a cultivated emmer accession and a durum wheat variety. A high-density single nucleotide polymorphism (SNP)-based genetic linkage map consisting of 2,593 markers was developed for the identification of quantitative trait loci. The major domestication gene Q had profound effects on spike length and compactness, rachis fragility, and threshability as expected. The cultivated emmer parent contributed increased spikelets per spike, and the durum parent contributed higher kernel weight, which led to the identification of some RILs that had significantly higher grain weight per spike than either parent. Threshability was governed not only by the Q locus, but other loci as well including Tg-B1 on chromosome 2B and a putative Tg-A1 locus on chromosome 2A indicating that mutations in the Tg loci occurred during the transition of cultivated emmer to the fully domesticated tetraploid. These results not only shed light on the events that shaped wheat domestication, but also demonstrate that cultivated emmer is a useful source of genetic variation for the enhancement of durum varieties.

A consensus framework map of durum wheat (Triticum durum Desf.) suitable for linkage disequilibrium analysis and genome-wide association mapping

BACKGROUND

Durum wheat (Triticum durum Desf.) is a tetraploid cereal grown in the medium to low-precipitation areas of the Mediterranean Basin, North America and South-West Asia. Genomics applications in durum wheat have the potential to boost exploitation of genetic resources and to advance understanding of the genetics of important complex traits (e.g. resilience to environmental and biotic stresses). A dense and accurate consensus map specific for T. durum will greatly facilitate genetic mapping, functional genomics and marker-assisted improvement.

RESULTS

High quality genotypic data from six core recombinant inbred line populations were used to obtain a consensus framework map of 598 simple sequence repeats (SSR) and Diversity Array Technology® (DArT) anchor markers (common across populations). Interpolation of unique markers from 14 maps allowed us to position a total of 2,575 markers in a consensus map of 2,463 cM. The T. durum A and B genomes were covered in their near totality based on the reference SSR hexaploid wheat map. The consensus locus order compared to those of the single component maps showed good correspondence, (average Spearman's rank correlation rho ρ value of 0.96). Differences in marker order and local recombination rate were observed between the durum and hexaploid wheat consensus maps. The consensus map was used to carry out a whole-genome search for genetic differentiation signatures and association to heading date in a panel of 183 accessions adapted to the Mediterranean areas. Linkage disequilibrium was found to decay below the r2 threshold=0.3 within 2.20 cM, on average. Strong molecular differentiations among sub-populations were mapped to 87 chromosome regions. A genome-wide association scan for heading date from 27 field trials in the Mediterranean Basin and in Mexico yielded 50 chromosome regions with evidences of association in multiple environments.

CONCLUSIONS

The consensus map presented here was used as a reference for genetic diversity and mapping analyses in T. durum, providing nearly complete genome coverage and even marker density. Markers previously mapped in hexaploid wheat constitute a strong link between the two species. The consensus map provides the basis for high-density single nucleotide polymorphic (SNP) marker implementation in durum wheat.

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.

Transcriptome and complexity-reduced, DNA-based identification of intraspecies single-nucleotide polymorphisms in the polyploid Gossypium hirsutum L

Varietal single nucleotide polymorphisms (SNPs) are the differences within one of the two subgenomes between different tetraploid cotton varieties and have not been practically used in cotton genetics and breeding because they are difficult to identify due to low genetic diversity and very high sequence identity between homeologous genes in cotton. We have used transcriptome and restriction site-associated DNA sequencing to identify varietal SNPs among 18 G. hirsutum varieties based on the rationale that varietal SNPs can be more confidently called when flanked by subgenome-specific SNPs. Using transcriptome data, we successfully identified 37,413 varietal SNPs and, of these, 22,121 did not have an additional varietal SNP within their 20-bp flanking regions so can be used in most SNP genotyping assays. From restriction site-associated DNA sequencing data, we identified an additional 3090 varietal SNPs between two of the varieties. Of the 1583 successful SNP assays achieved using different genotyping platforms, 1363 were verified. Many of the SNPs behaved as dominant markers because of coamplification from homeologous loci, but the number of SNPs acting as codominant markers increased when one or more subgenome-specific SNP(s) were incorporated in their assay primers, giving them greater utility for breeding applications. A G. hirsutum genetic map with 1244 SNP markers was constructed covering 5557.42 centiMorgan and used to map qualitative and quantitative traits. This collection of G. hirsutum varietal SNPs complements existing intra-specific SNPs and provides the cotton community with a valuable marker resource applicable to genetic analyses and breeding programs.

Genetic and molecular bases of yield-associated traits: a translational biology approach between rice and wheat

Transferring the knowledge bases between related species may assist in enlarging the yield potential of crop plants. Being cereals, rice and wheat share a high level of gene conservation; however, they differ at metabolic levels as a part of the environmental adaptation resulting in different yield capacities. This review focuses on the current understanding of genetic and molecular regulation of yield-associated traits in both crop species, highlights the similarities and differences and presents the putative knowledge gaps. We focus on the traits associated with phenology, photosynthesis, and assimilate partitioning and lodging resistance; the most important drivers of yield potential. Currently, there are large knowledge gaps in the genetic and molecular control of such major biological processes that can be filled in a translational biology approach in transferring genomics and genetics informations between rice and wheat.

A high-density genetic map with array-based markers facilitates structural and quantitative trait locus analyses of the common wheat genome

The large genome and allohexaploidy of common wheat have complicated construction of a high-density genetic map. Although improvements in the throughput of next-generation sequencing (NGS) technologies have made it possible to obtain a large amount of genotyping data for an entire mapping population by direct sequencing, including hexaploid wheat, a significant number of missing data points are often apparent due to the low coverage of sequencing. In the present study, a microarray-based polymorphism detection system was developed using NGS data obtained from complexity-reduced genomic DNA of two common wheat cultivars, Chinese Spring (CS) and Mironovskaya 808. After design and selection of polymorphic probes, 13,056 new markers were added to the linkage map of a recombinant inbred mapping population between CS and Mironovskaya 808. On average, 2.49 missing data points per marker were observed in the 201 recombinant inbred lines, with a maximum of 42. Around 40% of the new markers were derived from genic regions and 11% from repetitive regions. The low number of retroelements indicated that the new polymorphic markers were mainly derived from the less repetitive region of the wheat genome. Around 25% of the mapped sequences were useful for alignment with the physical map of barley. Quantitative trait locus (QTL) analyses of 14 agronomically important traits related to flowering, spikes, and seeds demonstrated that the new high-density map showed improved QTL detection, resolution, and accuracy over the original simple sequence repeat map.

Wheat in the Mediterranean revisited--tetraploid wheat landraces assessed with elite bread wheat Single Nucleotide Polymorphism markers

Background

Single Nucleotide Polymorphism (SNP) panels recently developed for the assessment of genetic diversity in wheat are primarily based on elite varieties, mostly those of bread wheat. The usefulness of such SNP panels for studying wheat evolution and domestication has not yet been fully explored and ascertainment bias issues can potentially affect their applicability when studying landraces and tetraploid ancestors of bread wheat. We here evaluate whether population structure and evolutionary history can be assessed in tetraploid landrace wheats using SNP markers previously developed for the analysis of elite cultivars of hexaploid wheat.

Results

We genotyped more than 100 tetraploid wheat landraces and wild emmer wheat accessions, some of which had previously been screened with SSR markers, for an existing SNP panel and obtained publically available genotypes for the same SNPs for hexaploid wheat varieties and landraces. Results showed that quantification of genetic diversity can be affected by ascertainment bias but that the effects of ascertainment bias can at least partly be alleviated by merging SNPs to haplotypes. Analyses of population structure and genetic differentiation show strong subdivision between the tetraploid wheat subspecies, except for durum and rivet that are not separable. A more detailed population structure of durum landraces could be obtained than with SSR markers. The results also suggest an emmer, rather than durum, ancestry of bread wheat and with gene flow from wild emmer.

Conclusions

SNP markers developed for elite cultivars show great potential for inferring population structure and can address evolutionary questions in landrace wheat. Issues of marker genome specificity and mapping need, however, to be addressed. Ascertainment bias does not seem to interfere with the ability of a SNP marker system developed for elite bread wheat accessions to detect population structure in other types of wheat.

SNP Discovery for mapping alien introgressions in wheat

Background

Monitoring alien introgressions in crop plants is difficult due to the lack of genetic and molecular mapping information on the wild crop relatives. The tertiary gene pool of wheat is a very important source of genetic variability for wheat improvement against biotic and abiotic stresses. By exploring the 5M^g short arm (5M^gS) of Aegilops geniculata, we can apply chromosome genomics for the discovery of SNP markers and their use for monitoring alien introgressions in wheat (Triticum aestivum L).

Results

The short arm of chromosome 5M^g of Ae. geniculata Roth (syn. Ae. ovata L.; 2n = 4x = 28, U^gU^gM^gM^g) was flow-sorted from a wheat line in which it is maintained as a telocentric chromosome. DNA of the sorted arm was amplified and sequenced using an Illumina Hiseq 2000 with ~45x coverage. The sequence data was used for SNP discovery against wheat homoeologous group-5 assemblies. A total of 2,178 unique, 5M^gS-specific SNPs were discovered. Randomly selected samples of 59 5M^gS-specific SNPs were tested (44 by KASPar assay and 15 by Sanger sequencing) and 84% were validated. Of the selected SNPs, 97% mapped to a chromosome 5M^g addition to wheat (the source of t5M^gS), and 94% to 5M^g introgressed from a different accession of Ae. geniculata substituting for chromosome 5D of wheat. The validated SNPs also identified chromosome segments of 5M^gS origin in a set of T5D-5M^g translocation lines; eight SNPs (25%) mapped to TA5601 [T5DL · 5DS-5M^gS(0.75)] and three (8%) to TA5602 [T5DL · 5DS-5M^gS (0.95)]. SNPs (gsnp_5ms83 and gsnp_5ms94), tagging chromosome T5DL · 5DS-5M^gS(0.95) with the smallest introgression carrying resistance to leaf rust (Lr57) and stripe rust (Yr40), were validated in two released germplasm lines with Lr57 and Yr40 genes.

Conclusion

This approach should be widely applicable for the identification of species/genome-specific SNPs. The development of a large number of SNP markers will facilitate the precise introgression and monitoring of alien segments in crop breeding programs and further enable mapping and cloning novel genes from the wild relatives of crop plants.

Transcriptome sequencing for high throughput SNP development and genetic mapping in Pea

BACKGROUND

Pea has a complex genome of 4.3 Gb for which only limited genomic resources are available to date. Although SNP markers are now highly valuable for research and modern breeding, only a few are described and used in pea for genetic diversity and linkage analysis.

RESULTS

We developed a large resource by cDNA sequencing of 8 genotypes representative of modern breeding material using the Roche 454 technology, combining both long reads (400 bp) and high coverage (3.8 million reads, reaching a total of 1,369 megabases). Sequencing data were assembled and generated a 68 K unigene set, from which 41 K were annotated from their best blast hit against the model species Medicago truncatula. Annotated contigs showed an even distribution along M. truncatula pseudochromosomes, suggesting a good representation of the pea genome. 10 K pea contigs were found to be polymorphic among the genetic material surveyed, corresponding to 35 K SNPs.We validated a subset of 1538 SNPs through the GoldenGate assay, proving their ability to structure a diversity panel of breeding germplasm. Among them, 1340 were genetically mapped and used to build a new consensus map comprising a total of 2070 markers. Based on blast analysis, we could establish 1252 bridges between our pea consensus map and the pseudochromosomes of M. truncatula, which provides new insight on synteny between the two species.

CONCLUSIONS

Our approach created significant new resources in pea, i.e. the most comprehensive genetic map to date tightly linked to the model species M. truncatula and a large SNP resource for both academic research and breeding.

Identification and characterization of nucleotide variations in the genome of Ziziphus jujuba (Rhamnaceae) by next generation sequencing

In this study, single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) in the genome of Ziziphus jujuba were identified using sequences generated by the Roche 454 GS-FLX sequencer. A total of, 573,141 reads were produced with an average read length of 360 bp. After quality control, 258,754 of the filtered reads were assembled into 23,864 contigs, and 293,458 remained as singletons. Using the contig assemblies as a reference, 17,160 SNPs and 478 InDels were identified. Among the SNPs, transitions occurred three times more frequently than transversions. In transitions, the number of C/T and G/A transitions was similar. Among the transversions, A/T was the most abundant, and C/G was much rarer than any of the other types of transversions, accounting for only about half the numbers of A/C, A/T and G/T transversions. For the InDels, mononucleotide changes amounted to 64.4% of the total number of InDels. In general, the frequency of detected InDels decreased as the length of the InDels increased. This study provides valuable marker resources for future genetic studies of Ziziphus spp.

Association mapping for root architectural traits in durum wheat seedlings as related to agronomic performance

Association mapping provides useful insights on the genetic architecture of quantitative traits across a large number of unrelated genotypes, which in turn allows an informed choice of the lines to be crossed for a more accurate characterization of major QTLs in a biparental genetic background. In this study, seedlings of 183 durum wheat elite accessions were evaluated in order to identify QTLs for root system architecture (RSA). The QTLs identified were compared with QTLs detected for grain yield and its component traits, plant height and peduncle length measured in a previous study where the same accessions were evaluated in 15 field trials with a broad range of soil moisture availability and productivity (Maccaferri et al. in J Exp Bot 62:409-438, 2011). The following RSA features were investigated in seedlings at the four-leaf stage: seminal root angle, primary root length, total root length, average root length, root number and shoot length. Highly significant differences among accessions were detected for all traits. The highest repeatability (h² = 0.72) was observed for seminal root angle. Out of the 48 QTLs detected for RSA, 15 overlapped with QTLs for agronomic traits and/or grain yield in two or more environments. The congruency of the effects of RSA traits and agronomic traits was evaluated. Seminal root angle and root number appear the most promising traits for further studies on the adaptive role of RSA plasticity on field performance in environments differing for water availability. Our results provide novel insights on the genetic control of RSA and its implications on field performance of durum wheat.

Genome-wide SNP calling using next generation sequencing data in tomato

The tomato (Solanum lycopersicum L.) is a model plant for genome research in Solanaceae, as well as for studying crop breeding. Genome-wide single nucleotide polymorphisms (SNPs) are a valuable resource in genetic research and breeding. However, to do discovery of genome-wide SNPs, most methods require expensive high-depth sequencing. Here, we describe a method for SNP calling using a modified version of SAMtools that improved its sensitivity. We analyzed 90 Gb of raw sequence data from next-generation sequencing of two resequencing and seven transcriptome data sets from several tomato accessions. Our study identified 4,812,432 non-redundant SNPs. Moreover, the workflow of SNP calling was improved by aligning the reference genome with its own raw data. Using this approach, 131,785 SNPs were discovered from transcriptome data of seven accessions. In addition, 4,680,647 SNPs were identified from the genome of S. pimpinellifolium, which are 60 times more than 71,637 of the PI212816 transcriptome. SNP distribution was compared between the whole genome and transcriptome of S. pimpinellifolium. Moreover, we surveyed the location of SNPs within genic and intergenic regions. Our results indicated that the sufficient genome-wide SNP markers and very sensitive SNP calling method allow for application of marker assisted breeding and genome-wide association studies.

Assessment of the geographic origins of pinewood nematode isolates via single nucleotide polymorphism in effector genes

The pinewood nematode, Bursaphelenchus xylophilus, is native to North America but it only causes damaging pine wilt disease in those regions of the world where it has been introduced. The accurate detection of the species and its dispersal routes are thus essential to define effective control measures. The main goals of this study were to analyse the genetic diversity among B. xylophilus isolates from different geographic locations and identify single nucleotide polymorphism (SNPs) markers for geographic origin, through a comparative transcriptomic approach. The transcriptomes of seven B. xylophilus isolates, from Continental Portugal (4), China (1), Japan (1) and USA (1), were sequenced in the next generation platform Roche 454. Analysis of effector gene transcripts revealed inter-isolate nucleotide diversity that was validated by Sanger sequencing in the genomic DNA of the seven isolates and eight additional isolates from different geographic locations: Madeira Island (2), China (1), USA (1), Japan (2) and South Korea (2). The analysis identified 136 polymorphic positions in 10 effector transcripts. Pairwise comparison of the 136 SNPs through Neighbor-Joining and the Maximum Likelihood methods and 5-mer frequency analysis with the alignment-independent bilinear multivariate modelling approach correlated the SNPs with the isolates geographic origin. Furthermore, the SNP analysis indicated a closer proximity of the Portuguese isolates to the Korean and Chinese isolates than to the Japanese or American isolates. Each geographic cluster carried exclusive alleles that can be used as SNP markers for B. xylophilus isolate identification.

Quantitative and qualitative differences in celiac disease epitopes among durum wheat varieties identified through deep RNA-amplicon sequencing

Background

Wheat gluten is important for the industrial quality of bread wheat (Triticum aestivum L.) and durum wheat (T. turgidum L.). Gluten proteins are also the source of immunogenic peptides that can trigger a T cell reaction in celiac disease (CD) patients, leading to inflammatory responses in the small intestine. Various peptides with three major T cell epitopes involved in CD are derived from alpha-gliadin fraction of gluten. Alpha-gliadins are encoded by a large multigene family and amino acid variation in the CD epitopes is known to influence the immunogenicity of individual gene family members. Current commercial methods of gluten detection are unable to distinguish between immunogenic and non-immunogenic CD epitope variants and thus to accurately quantify the overall CD epitope load of a given wheat variety. Such quantification is indispensable for correct selection of wheat varieties with low potential to cause CD.

Results

A 454 RNA-amplicon sequencing method was developed for alpha-gliadin transcripts encompassing the three major CD epitopes and their variants. The method was used to screen developing grains on plants of 61 different durum wheat cultivars and accessions. A dedicated sequence analysis pipeline returned a total of 304 unique alpha-gliadin transcripts, corresponding to a total of 171 ‘unique deduced protein fragments’ of alpha-gliadins. The numbers of these fragments obtained in each plant were used to calculate quantitative and quantitative differences between the CD epitopes expressed in the endosperm of these wheat plants. A few plants showed a lower fraction of CD epitope-encoding alpha-gliadin transcripts, but none were free of CD epitopes.

Conclusions

The dedicated 454 RNA-amplicon sequencing method enables 1) the grouping of wheat plants according to the genetic variation in alpha-gliadin transcripts, and 2) the screening for plants which are potentially less CD-immunogenic. The resulting alpha-gliadin sequence database will be useful as a reference in proteomics analysis regarding the immunogenic potential of mature wheat grains.

Wheat syntenome unveils new evidences of contrasted evolutionary plasticity between paleo- and neoduplicated subgenomes

Bread wheat derives from a grass ancestor structured in seven protochromosomes followed by a paleotetraploidization to reach a 12 chromosomes intermediate and a neohexaploidization (involving subgenomes A, B and D) event that finally shaped the 21 modern chromosomes. Insights into wheat syntenome in sequencing conserved orthologous set (COS) genes unravelled differences in genomic structure (such as gene conservation and diversity) and genetical landscape (such as recombination pattern) between ancestral as well as recent duplicated blocks. Contrasted evolutionary plasticity is observed where the B subgenome appears more sensitive (i.e. plastic) in contrast to A as dominant (i.e. stable) in response to the neotetraploidization and D subgenome as supra-dominant (i.e. pivotal) in response to the neohexaploidization event. Finally, the wheat syntenome, delivered through a public web interface PlantSyntenyViewer at http://urgi.versailles.inra.fr/synteny-wheat, can be considered as a guide for accelerated dissection of major agronomical traits in wheat.

Next generation characterisation of cereal genomes for marker discovery

Cereal crops form the bulk of the world’s food sources, and thus their importance cannot be understated. Crop breeding programs increasingly rely on high-resolution molecular genetic markers to accelerate the breeding process. The development of these markers is hampered by the complexity of some of the major cereal crop genomes, as well as the time and cost required. In this review, we address current and future methods available for the characterisation of cereal genomes, with an emphasis on faster and more cost effective approaches for genome sequencing and the development of markers for trait association and marker assisted selection (MAS) in crop breeding programs.

Disentangling homeologous contigs in allo-tetraploid assembly: application to durum wheat

BACKGROUND

Using Next Generation Sequencing, SNP discovery is relatively easy on diploid species and still hampered in polyploid species by the confusion due to homeology. We develop HomeoSplitter; a fast and effective solution to split original contigs obtained by RNAseq into two homeologous sequences. It uses the differential expression of the two homeologous genes in the RNA. We verify that the new sequences are closer to the diploid progenitors of the allopolyploid species than the original contig. By remapping original reads on these new sequences, we also verify that the number of valuable detected SNPs has significantly increased.

RESULTS

HomeoSplitter is a fast and effective solution to disentangle homeologous sequences based on a maximum likelihood optimization. On a benchmark set of 2,505 clusters containing homologous sequences of urartu, speltoides and durum, HomeoSplitter was efficient to build sequences closer to the diploid references and increased the number of valuable SNPs from 188 out of 1,360 SNPs detected when mapping the reads on the de novo durum assembly to 762 out of 1,620 SNPs when mapping on HomeoSplitter contigs.

CONCLUSIONS

The HomeoSplitter program is freely available at http://bioweb.supagro.inra.fr/homeoSplitter/. This work provides a practical solution to the complex problem of disentangling homeologous transcripts in allo-tetraploids, which further allows an improved SNP detection.

Sequence-based SNP genotyping in durum wheat

Marker development for marker-assisted selection in plant breeding is increasingly based on next-generation sequencing (NGS). However, marker development in crops with highly repetitive, complex genomes is still challenging. Here we applied sequence-based genotyping (SBG), which couples AFLP®-based complexity reduction to NGS, for de novo single nucleotide polymorphisms (SNP) marker discovery in and genotyping of a biparental durum wheat population. We identified 9983 putative SNPs in 6372 contigs between the two parents and used these SNPs for genotyping 91 recombinant inbred lines (RILs). Excluding redundant information from multiple SNPs per contig, 2606 (41%) markers were used for integration in a pre-existing framework map, resulting in the integration of 2365 markers over 2607 cM. Of the 2606 markers available for mapping, 91% were integrated in the pre-existing map, containing 708 SSRs, DArT markers, and SNPs from CRoPS technology, with a map-size increase of 492 cM (23%). These results demonstrate the high quality of the discovered SNP markers. With this methodology, it was possible to saturate the map at a final marker density of 0.8 cM/marker. Looking at the binned marker distribution (Figure 2), 63 of the 268 10-cM bins contained only SBG markers, showing that these markers are filling in gaps in the framework map. As to the markers that could not be used for mapping, the main reason was the low sequencing coverage used for genotyping. We conclude that SBG is a valuable tool for efficient, high-throughput and high-quality marker discovery and genotyping for complex genomes such as that of durum wheat.

Sequence-based mapping of the polyploid wheat genome

The emergence of new sequencing technologies has provided fast and cost-efficient strategies for high-resolution mapping of complex genomes. Although these approaches hold great promise to accelerate genome analysis, their application in studying genetic variation in wheat has been hindered by the complexity of its polyploid genome. Here, we applied the next-generation sequencing of a wheat doubled-haploid mapping population for high-resolution gene mapping and tested its utility for ordering shotgun sequence contigs of a flow-sorted wheat chromosome. A bioinformatical pipeline was developed for reliable variant analysis of sequence data generated for polyploid wheat mapping populations. The results of variant mapping were consistent with the results obtained using the wheat 9000 SNP iSelect assay. A reference map of the wheat genome integrating 2740 gene-associated single-nucleotide polymorphisms from the wheat iSelect assay, 1351 diversity array technology, 118 simple sequence repeat/sequence-tagged sites, and 416,856 genotyping-by-sequencing markers was developed. By analyzing the sequenced megabase-size regions of the wheat genome we showed that mapped markers are located within 40-100 kb from genes providing a possibility for high-resolution mapping at the level of a single gene. In our population, gene loci controlling a seed color phenotype cosegregated with 2459 markers including one that was located within the red seed color gene. We demonstrate that the high-density reference map presented here is a useful resource for gene mapping and linking physical and genetic maps of the wheat genome.

Searching for novel sources of field resistance to Ug99 and Ethiopian stem rust races in durum wheat via association mapping

Puccinia graminis f. sp. tritici, the causative agent of stem rust in wheat, is a devastating disease of durum wheat. While more than 50 stem rust resistance (Sr) loci have been identified in wheat, only a few of them have remained effective against Ug99 (TTKSK race) and other durum-specific Ethiopian races. An association mapping (AM) approach based on 183 diverse durum wheat accessions was utilized to identify resistance loci for stem rust response in Ethiopia over four field-evaluation seasons and artificial inoculation with Ug99 and a mixture of durum-specific races. The panel was profiled with simple sequence repeat, Diversity Arrays Technology and sequence-tagged site markers (1,253 in total). The resistance turned out to be oligogenic, with twelve QTL-tagging markers that were significant (P < 0.05) across three or four seasons. R (2) values ranged from 1.1 to 11.3 %.Twenty-four additional single-marker/QTL regions were found to be significant over two seasons. The AM results confirmed the role of Sr13, previously described in bi-parental mapping studies, and the role of chromosome regions putatively harbouring Sr9, Sr14, Sr17 and Sr28. Three minor QTLs were coincident with those reported in hexaploid wheat and five overlapped with those recently reported in the Sebatel × Kristal durum mapping population. Thirteen single-marker/QTL regions were located in chromosome regions where no Sr genes/QTLs have been previously reported. The allelic variation identified in this study is readily available and can be exploited for marker-assisted selection, thus providing additional opportunities for a more durable stem rust resistance under field conditions.

Application of large-scale sequencing to marker discovery in plants

Advances in DNA sequencing provide tools for efficient large-scale discovery of markers for use in plants. Discovery options include large-scale amplicon sequencing, transcriptome sequencing, gene-enriched genome sequencing and whole genome sequencing. Examples of each of these approaches and their potential to generate molecular markers for specific applications have been described. Sequencing the whole genome of parents identifies all the polymorphisms available for analysis in their progeny. Sequencing PCR amplicons of sets of candidate genes from DNA bulks can be used to define the available variation in these genes that might be exploited in a population or germplasm collection. Sequencing of the transcriptomes of genotypes varying for the trait of interest may identify genes with patterns of expression that could explain the phenotypic variation. Sequencing genomic DNA enriched for genes by hybridization with probes for all or some of the known genes simplifies sequencing and analysis of differences in gene sequences between large numbers of genotypes and genes especially when working with complex genomes. Examples of application of the above-mentioned techniques have been described.

Identification and mapping of leaf, stem and stripe rust resistance quantitative trait loci and their interactions in durum wheat

Leaf rust (Puccinia triticina Eriks.), stripe rust (Puccinia striiformis f. tritici Eriks.) and stem rust (Puccinia graminis f. sp. tritici) cause major production losses in durum wheat (Triticum turgidum L. var. durum). The objective of this research was to identify and map leaf, stripe and stem rust resistance loci from the French cultivar Sachem and Canadian cultivar Strongfield. A doubled haploid population from Sachem/Strongfield and parents were phenotyped for seedling reaction to leaf rust races BBG/BN and BBG/BP and adult plant response was determined in three field rust nurseries near El Batan, Obregon and Toluca, Mexico. Stripe rust response was recorded in 2009 and 2011 nurseries near Toluca and near Njoro, Kenya in 2010. Response to stem rust was recorded in field nurseries near Njoro, Kenya, in 2010 and 2011. Sachem was resistant to leaf, stripe and stem rust. A major leaf rust quantitative trait locus (QTL) was identified on chromosome 7B at Xgwm146 in Sachem. In the same region on 7B, a stripe rust QTL was identified in Strongfield. Leaf and stripe rust QTL around DArT marker wPt3451 were identified on chromosome 1B. On chromosome 2B, a significant leaf rust QTL was detected conferred by Strongfield, and at the same QTL, a Yr gene derived from Sachem conferred resistance. Significant stem rust resistance QTL were detected on chromosome 4B. Consistent interactions among loci for resistance to each rust type across nurseries were detected, especially for leaf rust QTL on 7B. Sachem and Strongfield offer useful sources of rust resistance genes for durum rust breeding.