High-density SNP genotyping array for hexaploid wheat and its secondary and tertiary gene pool

Detection of ploidy and chromosomal aberrations in commercial oil palm using high-throughput SNP markers

Karyotyping using high-density genome-wide SNP markers identified various chromosomal aberrations in oil palm (Elaeis guineensis Jacq.) with supporting evidence from the 2C DNA content measurements (determined using FCM) and chromosome counts. Oil palm produces a quarter of the world's total vegetable oil. In line with its global importance, an initiative to sequence the oil palm genome was carried out successfully, producing huge amounts of sequence information, allowing SNP discovery. High-capacity SNP genotyping platforms have been widely used for marker-trait association studies in oil palm. Besides genotyping, a SNP array is also an attractive tool for understanding aberrations in chromosome inheritance. Exploiting this, the present study utilized chromosome-wide SNP allelic distributions to determine the ploidy composition of over 1,000 oil palms from a commercial F₁ family, including 197 derived from twin-embryo seeds. Our method consisted of an inspection of the allelic intensity ratio using SNP markers. For palms with a shifted or abnormal distribution ratio, the SNP allelic frequencies were plotted along the pseudo-chromosomes. This method proved to be efficient in identifying whole genome duplication (triploids) and aneuploidy. We also detected several loss of heterozygosity regions which may indicate small chromosomal deletions and/or inheritance of identical by descent regions from both parents. The SNP analysis was validated by flow cytometry and chromosome counts. The triploids were all derived from twin-embryo seeds. This is the first report on the efficiency and reliability of SNP array data for karyotyping oil palm chromosomes, as an alternative to the conventional cytogenetic technique. Information on the ploidy composition and chromosomal structural variation can help to better understand the genetic makeup of samples and lead to a more robust interpretation of the genomic data in marker-trait association analyses.

The role of gene flow and chromosomal instability in shaping the bread wheat genome

Bread wheat (Triticum aestivum) is one of the world's most important crops; however, a low level of genetic diversity within commercial breeding accessions can significantly limit breeding potential. In contrast, wheat relatives exhibit considerable genetic variation and so potentially provide a valuable source of novel alleles for use in breeding new cultivars. Historically, gene flow between wheat and its relatives may have contributed novel alleles to the bread wheat pangenome. To assess the contribution made by wheat relatives to genetic diversity in bread wheat, we used markers based on single nucleotide polymorphisms to compare bread wheat accessions, created in the past 150 years, with 45 related species. We show that many bread wheat accessions share near-identical haplotype blocks with close relatives of wheat's diploid and tetraploid progenitors, while some show evidence of introgressions from more distant species and structural variation between accessions. Hence, introgressions and chromosomal rearrangements appear to have made a major contribution to genetic diversity in cultivar collections. As gene flow from relatives to bread wheat is an ongoing process, we assess the impact that introgressions might have on future breeding strategies.

New resources for genetic studies in maize (Zea mays L.): a genome-wide Maize6H-60K single nucleotide polymorphism array and its application

Despite the availability of numerous molecular markers in maize, effective evaluation of all types of germplasm resources, accurate identification of varieties and analysis of a large number of materials in a timely, low-cost manner is challenging. Here, we present Maize6H-60K, a genome-wide single nucleotide polymorphism (SNP) array to facilitate maize genotyping. We first identified 160 million variants by sequencing data of 388 representative inbreds and then tiled 200 000 high-quality variants on a screening array. These variants were further narrowed down to 61 282 using stringent filtering criteria. Among the 60 000 markers, 21 460 SNPs (35%) were within genic regions and 12 835 (21%) were located in coding regions. To assess their effectiveness, 329 inbreds, 221 hybrids, 34 parent-offspring sets and six breeding samples were genotyped. Overall, 48 972 SNPs (80%) were categorized into the highest quality class, that of 'poly high resolution'. A total of 54 658 (89.29%) and 53 091 (86.73%) SNPs had minor allele frequency values ≥ 0.20 in inbreds and hybrids respectively. A linkage disequilibrium (LD) analysis revealed that LD decline was in equilibrium when r² was between 0.10 and 0.15, which corresponds to a physical distance of 400-600 kb. UPGMA clustering analysis divided the 329 inbred lines into nine groups that were consistent with known pedigrees. A background analysis of breeding materials indicated that the 60 000 markers were suitable for evaluation of breeding populations constructed by materials between or within heterotic groups. The developed Maize6H-60K array should be an important tool in maize genetic studies, variety identification and molecular breeding.

Genomic resources in plant breeding for sustainable agriculture

Climate change during the last 40 years has had a serious impact on agriculture and threatens global food and nutritional security. From over half a million plant species, cereals and legumes are the most important for food and nutritional security. Although systematic plant breeding has a relatively short history, conventional breeding coupled with advances in technology and crop management strategies has increased crop yields by 56 % globally between 1965-85, referred to as the Green Revolution. Nevertheless, increased demand for food, feed, fiber, and fuel necessitates the need to break existing yield barriers in many crop plants. In the first decade of the 21st century we witnessed rapid discovery, transformative technological development and declining costs of genomics technologies. In the second decade, the field turned towards making sense of the vast amount of genomic information and subsequently moved towards accurately predicting gene-to-phenotype associations and tailoring plants for climate resilience and global food security. In this review we focus on genomic resources, genome and germplasm sequencing, sequencing-based trait mapping, and genomics-assisted breeding approaches aimed at developing biotic stress resistant, abiotic stress tolerant and high nutrition varieties in six major cereals (rice, maize, wheat, barley, sorghum and pearl millet), and six major legumes (soybean, groundnut, cowpea, common bean, chickpea and pigeonpea). We further provide a perspective and way forward to use genomic breeding approaches including marker-assisted selection, marker-assisted backcrossing, haplotype based breeding and genomic prediction approaches coupled with machine learning and artificial intelligence, to speed breeding approaches. The overall goal is to accelerate genetic gains and deliver climate resilient and high nutrition crop varieties for sustainable agriculture.

Evaluation of genetic structure in European wheat cultivars and advanced breeding lines using high-density genotyping-by-sequencing approach

BACKGROUND

The genetic diversity and gene pool characteristics must be clarified for efficient genome-wide association studies, genomic selection, and hybrid breeding. The aim of this study was to evaluate the genetic structure of 509 wheat accessions representing registered varieties and advanced breeding lines via the high-density genotyping-by-sequencing approach.

RESULTS

More than 30% of 13,499 SNP markers representing 2162 clusters were mapped to genes, whereas 22.50% of 26,369 silicoDArT markers overlapped with coding sequences and were linked in 3527 blocks. Regarding hexaploidy, perfect sequence matches following BLAST searches were not sufficient for the unequivocal mapping to unique loci. Moreover, allelic variations in homeologous loci interfered with heterozygosity calculations for some markers. Analyses of the major genetic changes over the last 27 years revealed the selection pressure on orthologs of the gibberellin biosynthesis-related GA2 gene and the senescence-associated SAG12 gene. A core collection representing the wheat population was generated for preserving germplasm and optimizing breeding programs.

CONCLUSIONS

Our results confirmed considerable differences among wheat subgenomes A, B and D, with D characterized by the lowest diversity but the highest LD. They revealed genomic regions that have been targeted by breeding.

Linkage mapping identifies a non-synonymous mutation in FLOWERING LOCUS T (FT-B1) increasing spikelet number per spike

Total spikelet number per spike (TSN) is a major component of spike architecture in wheat (Triticumaestivum L.). A major and consistent quantitative trait locus (QTL) was discovered for TSN in a doubled haploid spring wheat population grown in the field over 4 years. The QTL on chromosome 7B explained up to 20.5% of phenotypic variance. In its physical interval (7B: 6.37–21.67 Mb), the gene FLOWERINGLOCUST (FT-B1) emerged as candidate for the observed effect. In one of the parental lines, FT-B1 carried a non-synonymous substitution on position 19 of the coding sequence. This mutation modifying an aspartic acid (D) into a histidine (H) occurred in a highly conserved position. The mutation was observed with a frequency of ca. 68% in a set of 135 hexaploid wheat varieties and landraces, while it was not found in other plant species. FT-B1 only showed a minor effect on heading and flowering time (FT) which were dominated by a major QTL on chromosome 5A caused by segregation of the vernalization gene VRN-A1. Individuals carrying the FT-B1 allele with amino acid histidine had, on average, a higher number of spikelets (15.1) than individuals with the aspartic acid allele (14.3) independent of their VRN-A1 allele. We show that the effect of TSN is not mainly related to flowering time; however, the duration of pre-anthesis phases may play a major role.

Genome-wide association study of yield and related traits in common wheat under salt-stress conditions

BACKGROUND

Soil salinization is a major threat to wheat production. It is essential to understand the genetic basis of salt tolerance for breeding and selecting new salt-tolerant cultivars that have the potential to increase wheat yield.

RESULT

In this study, a panel of 191 wheat accessions was subjected to genome wide association study (GWAS) to identify SNP markers linked with adult-stage characters. The population was genotyped by Wheat660K SNP array and eight phenotype traits were investigated under low and high salinity environments for three consecutive years. A total of 389 SNPs representing 11 QTLs were significantly associated with plant height, spike number, spike length, grain number, thousand kernels weight, yield and biological mass under different salt treatments, with the phenotypic explanation rate (R²) ranging from 9.14 to 50.45%. Of these, repetitive and pleiotropic loci on chromosomes 4A, 5A, 5B and 7A were significantly linked to yield and yield related traits such as thousand kernels weight, spike number, spike length, grain number and so on under low salinity conditions. Spike length-related loci were mainly located on chromosomes 1B, 3B, 5B and 7A under different salt treatments. Two loci on chromosome 4D and 5A were related with plant height in low and high salinity environment, respectively. Three salt-tolerant related loci were confirmed to be important in two bi-parental populations. Distribution of favorable haplotypes indicated that superior haplotypes of pleiotropic loci on group-5 chromosomes were strongly selected and had potential for increasing wheat salt tolerance. A total of 14 KASP markers were developed for nine loci associating with yield and related traits to improve the selection efficiency of wheat salt-tolerance breeding.

CONCLUSION

Utilizing a Wheat660K SNPs chip, QTLs for yield and its related traits were detected under salt treatment in a natural wheat population. Important salt-tolerant related loci were validated in RIL and DH populations. This study provided reliable molecular markers that could be crucial for marker-assisted selection in wheat salt tolerance breeding programs.

Generation of Doubled Haploid Wheat-Triticum urartu Introgression Lines and Their Characterisation Using Chromosome-Specific KASP Markers

Wheat is one of the most important food and protein sources in the world and although, in recent years wheat breeders have achieved yield gains, they are not sufficient to meet the demands of an ever-growing population. Development of high yielding wheat varieties, resilient to abiotic and biotic stress resulting from climate change, has been limited by wheat's narrow genetic base. In contrast to wheat, the wild relatives of wheat provide a vast reservoir of genetic variation for most, if not all, agronomic traits. Previous studies by the authors have shown the transfer of genetic variation from T. urartu into bread wheat. However, before the introgression lines can be exploited for trait analysis, they are required to have stable transmission of the introgressions to the next generation. In this work, we describe the generation of 86 doubled haploid (DH) wheat-T. urartu introgression lines that carry homozygous introgressions which are stably inherited. The DH lines were characterised using the Axiom^® Wheat Relative Genotyping Array and 151 KASP markers to identify 65 unique T. urartu introgressions in a bread wheat background. DH production has helped accelerate the breeding process and facilitated the early release of homozygous wheat-T. urartu introgression lines. Together with the KASP markers, this valuable resource could greatly advance identification of beneficial alleles that can be used in wheat improvement.

Mapping QTL for agronomic traits under two levels of salt stress in a new constructed RIL wheat population

QTL for 15 agronomic traits under two levels of salt stress in dry salinity field were mapped in a new constructed RIL population utilizing a Wheat55K SNP array. Furthermore, eight QTL were validated in a collected natural population. Soil salinity is one of the major abiotic stresses causing serious impact on crop growth, development and yield. As one of the three most important crops in the world, bread wheat (Triticum aestivum L.) is severely affected by salinity, too. In this study, an F₇ recombinant inbred line (RIL) population derived from a cross between high-yield wheat cultivar Zhongmai 175 and salt-tolerant cultivar Xiaoyan 60 was constructed. The adult stage performances of the RIL population and their parent lines under low and high levels of salt stress were evaluated for three consecutive growing seasons. Utilizing a Wheat55K SNP array, a high-density genetic linkage map spinning 3250.71 cM was constructed. QTL mapping showed that 90 stable QTL for 15 traits were detected, and they were distributed on all wheat chromosomes except 4D, 6B and 7D. These QTL individually explained 2.34-32.43% of the phenotypic variation with LOD values ranging from 2.68 to 47.15. It was found that four QTL clusters were located on chromosomes 2D, 3D, 4B and 6A, respectively. Notably, eight QTL from the QTL clusters were validated in a collected natural population. Among them, QPh-4B was deduced to be an allele of Rht-B1. In addition, three kompetitive allele-specific PCR (KASP) markers derived from SNPs were successfully designed for three QTL clusters. This study provides an important base for salt-tolerant QTL (gene) cloning in wheat, and the markers, especially the KASP markers, will be useful for marker-assisted selection in salt-tolerant wheat breeding.

Genomic prediction of agronomic traits in wheat using different models and cross-validation designs

Genomic predictions across environments and within populations resulted in moderate to high accuracies but across-population genomic prediction should not be considered in wheat for small population size. Genomic selection (GS) is a marker-based selection suggested to improve the genetic gain of quantitative traits in plant breeding programs. We evaluated the effects of training population (TP) composition, cross-validation design, and genetic relationship between the training and breeding populations on the accuracy of GS in spring wheat (Triticum aestivum L.). Two populations of 231 and 304 spring hexaploid wheat lines that were phenotyped for six agronomic traits and genotyped with the wheat 90 K array were used to assess the accuracy of seven GS models (RR-BLUP, G-BLUP, BayesB, BL, RKHS, GS + de novo GWAS, and reaction norm) using different cross-validation designs. BayesB outperformed the other models for within-population genomic predictions in the presence of few quantitative trait loci (QTL) with large effects. However, including fixed-effect marker covariates gave better performance for an across-population prediction when the same QTL underlie traits in both populations. The accuracy of prediction was highly variable based on the cross-validation design, which suggests the importance to use a design that resembles the variation within a breeding program. Moderate to high accuracies were obtained when predictions were made within populations. In contrast, across-population genomic prediction accuracies were very low, suggesting that the evaluated models are not suitable for prediction across independent populations. On the other hand, across-environment prediction and forward prediction designs using the reaction norm model resulted in moderate to high accuracies, suggesting that GS can be applied in wheat to predict the performance of newly developed lines and lines in incomplete field trials.

Identification and cross-validation of genetic loci conferring resistance to Septoria nodorum blotch using a German multi-founder winter wheat population

We identified allelic variation at two major loci, QSnb.nmbu-2A.1 and QSnb.nmbu-5A.1, showing consistent and additive effects on SNB field resistance. Validation of QSnb.nmbu-2A.1 across genetic backgrounds further highlights its usefulness for marker-assisted selection. Septoria nodorum blotch (SNB) is a disease of wheat (Triticum aestivum and T. durum) caused by the necrotrophic fungal pathogen Parastagonospora nodorum. SNB resistance is a typical quantitative trait, controlled by multiple quantitative trait loci (QTL) of minor effect. To achieve increased plant resistance, selection for resistance alleles and/or selection against susceptibility alleles must be undertaken. Here, we performed genetic analysis of SNB resistance using an eight-founder German Multiparent Advanced Generation Inter-Cross (MAGIC) population, termed BMWpop. Field trials and greenhouse testing were conducted over three seasons in Norway, with genetic analysis identifying ten SNB resistance QTL. Of these, two QTL were identified over two seasons: QSnb.nmbu-2A.1 on chromosome 2A and QSnb.nmbu-5A.1 on chromosome 5A. The chromosome 2A BMWpop QTL co-located with a robust SNB resistance QTL recently identified in an independent eight-founder MAGIC population constructed using varieties released in the United Kingdom (UK). The validation of this SNB resistance QTL in two independent multi-founder mapping populations, regardless of the differences in genetic background and agricultural environment, highlights the value of this locus in SNB resistance breeding. The second robust QTL identified in the BMWpop, QSnb.nmbu-5A.1, was not identified in the UK MAGIC population. Combining resistance alleles at both loci resulted in additive effects on SNB resistance. Therefore, using marker assisted selection to combine resistance alleles is a promising strategy for improving SNB resistance in wheat breeding. Indeed, the multi-locus haplotypes determined in this study provide markers for efficient tracking of these beneficial alleles in future wheat genetics and breeding activities.

Development of a minimal KASP marker panel for distinguishing genotypes in apple collections

Accurate identification of named accessions in germplasm collections is extremely important, especially for vegetatively propagated crops which are expensive to maintain. Thus, an inexpensive, reliable, and rapid genotyping method is essential because it avoids the need for laborious and time-consuming morphological comparisons. Single Nucleotide Polymorphism (SNP) marker panels containing large numbers of SNPs have been developed for many crop species, but such panels are much too large for basic cultivar identification. Here, we have identified a minimum set of SNP markers sufficient to distinguish apple cultivars held in the English and Welsh national collections providing a cheaper and automatable alternative to the markers currently used by the community. We show that SNP genotyping with a small set of well selected markers is equally efficient as microsatellites for the identification of apple cultivars and has the added advantage of automation and reduced cost when screening large numbers of samples.

Identification and validation of quantitative trait loci for kernel traits in common wheat (Triticum aestivum L.)

BACKGROUND

Kernel weight and morphology are important traits affecting cereal yields and quality. Dissecting the genetic basis of thousand kernel weight (TKW) and its related traits is an effective method to improve wheat yield.

RESULTS

In this study, we performed quantitative trait loci (QTL) analysis using recombinant inbred lines derived from the cross 'PuBing3228 × Gao8901' (PG-RIL) to dissect the genetic basis of kernel traits. A total of 17 stable QTLs related to kernel traits were identified, notably, two stable QTLs QTkw.cas-1A.2 and QTkw.cas-4A explained the largest portion of the phenotypic variance for TKW and kernel length (KL), and the other two stable QTLs QTkw.cas-6A.1 and QTkw.cas-7D.2 contributed more effects on kernel width (KW). Conditional QTL analysis revealed that the stable QTLs for TKW were mainly affected by KW. The QTLs QTkw.cas-7D.2 and QKw.cas-7D.1 associated with TKW and KW were delimited to the physical interval of approximately 3.82 Mb harboring 47 candidate genes. Among them, the candidate gene TaFT-D1 had a 1 bp insertions/deletion (InDel) within the third exon, which might be the reason for diversity in TKW and KW between the two parents. A Kompetitive Allele-Specific PCR (KASP) marker of TaFT-D1 allele was developed and verified by PG-RIL and a natural population consisted of 141 cultivar/lines. It was found that the favorable TaFT-D1 (G)-allele has been positively selected during Chinese wheat breeding. Thus, these results can be used for further positional cloning and marker-assisted selection in wheat breeding programs.

CONCLUSIONS

Seventeen stable QTLs related to kernel traits were identified. The stable QTLs for thousand kernel weight were mainly affected by kernel width. TaFT-D1 could be the candidate gene for QTLs QTkw.cas-7D.2 and QKw.cas-7D.1.

A novel adult plant leaf rust resistance gene Lr2K38 mapped on wheat chromosome 1AL

Soft red winter wheat (SRWW) cultivar AGS 2038 has a high level of seedling and adult plant leaf rust (LR) resistance. To map and characterize LR resistance in AGS 2038, a recombinant inbred line (RIL) population consisting of 225 lines was developed from a cross between AGS 2038 and moderately resistant line UGA 111729. The parents and RIL population were phenotyped for LR response in three field environments at Plains and Griffin, GA, in the 2017-2018 and 2018-2019 growing seasons, one greenhouse environment at the adult-plant stage, and at seedling stage. The RIL population was genotyped with the Illumina iSelect 90K SNP marker array, and a total of 7667 polymorphic markers representing 1513 unique loci were used to construct a linkage map. Quantitative trait loci (QTL) analysis detected six QTL, QLr.ags-1AL, QLr.ags-2AS, QLr.ags-2BS1, QLr.ags-2BS2, QLr.ags-2BS3, and QLr.ags-2DS, for seedling and adult plant LR resistance. Of these, the major adult plant leaf rust resistance QTL, QLr.ags-1AL, was detected on all field and greenhouse adult plant tests and explained up to 34.45% of the phenotypic variation. QLr.ags-1AL, tightly flanked by IWB20487 and IWA4022 markers, was contributed by AGS 2038. Molecular marker analysis using a diagnostic marker linked to Lr59 showed that QLr.ags-1AL was different from Lr59, the only known LR resistance gene on 1AL. Therefore, the QTL was temporarily designated as Lr2K38. Lr2K38-linked marker IWB20487 was highly polymorphic among 30 SRWW lines and should be useful for selecting the Lr2K38 in wheat breeding programs.

Genome-wide identification and characterization of novel non-coding RNA-derived SSRs in wheat

Expression of eukaryotic genes is largely regulated by non-coding RNAs (ncRNA). Sequence variations in the regulatory RNAs may have critical biological consequences including transcriptional and post-transcriptional gene regulation. ncRNA-derived markers thus can be proved useful in molecular breeding, QTL mapping and association studies for trait dissection. In present study, we identified a total of 661 SSRs dwelling in pre-miRNA (15), small nuclear RNA (25) and lncRNA (621). Of these, 46 were validated and 100% amplification success was observed in selected wheat genotypes. A set of 36 ncRNA-SSRs markers was utilized for genetic variability assessment in forty-eight Indian wheat genotypes (which includes bread wheat, durum wheat and relatives). Number of alleles ranged from 1 to 4 with an average of two alleles per SSR locus. Mean PIC, observed heterozygosity and Shannon information index were found to be 0.258, 0.37 and 0.476 which suggests ncRNA-SSRs show higher polymorphism compared to genic SSRs but lower polymorphism compared to genomic SSRs. Thirty-six ncRNA-SSRs showed transferability ranging from 42.1% to 100%. Average genetic dissimilarity among wheat genotypes was found to be 0.29 based on Jaccard's dissimilarity. This is the first report of ncRNA-SSRs in wheat which will be useful for molecular breeding and genetic improvement of wheat.

Multi-locus genome-wide association studies reveal novel genomic regions associated with vegetative stage salt tolerance in bread wheat (Triticum aestivum L.)

Soil salinity is one of the typical abiotic stresses affecting sustainability of wheat production worldwide. In the present study, we performed a 35 K SNP genotyping assay on association panel of 135 diverse wheat genotypes evaluated for vegetative stage tolerance in hydroponics. Association analyses using five multi-locus GWAS models revealed 42 reliable QTNs for 10 salt tolerance associated traits. Among these 42 reliable QTNs, 9, 17 and 16 QTNs were associated with physiological, biomass and shoot ionic traits respectively. Novel major QTNs were identified for chlorophyll content, shoot fresh weight, seedling total biomass, Na⁺ and K⁺ concentration and Na⁺/K⁺ ratio in shoots. Further, 10 major QTNs showed significant effect on the corresponding salt tolerance traits. Gene ontology analysis of the associated genomic regions identified 58 candidate genes. The information generated in this study will be of potential value for improvement of salt tolerance of wheat cultivars using marker assisted selection.

Appraising the Genetic Architecture of Kernel Traits in Hexaploid Wheat Using GWAS

Kernel morphology is one of the major yield traits of wheat, the genetic architecture of which is always important in crop breeding. In this study, we performed a genome-wide association study (GWAS) to appraise the genetic architecture of the kernel traits of 319 wheat accessions using 22,905 single nucleotide polymorphism (SNP) markers from a wheat 90K SNP array. As a result, 111 and 104 significant SNPs for Kernel traits were detected using four multi-locus GWAS models (mrMLM, FASTmrMLM, FASTmrEMMA, and pLARmEB) and three single-locus models (FarmCPU, MLM, and MLMM), respectively. Among the 111 SNPs detected by the multi-locus models, 24 SNPs were simultaneously detected across multiple models, including seven for kernel length, six for kernel width, six for kernels per spike, and five for thousand kernel weight. Interestingly, the five most stable SNPs (RAC875_29540_391, Kukri_07961_503, tplb0034e07_1581, BS00074341_51, and BobWhite_049_3064) were simultaneously detected by at least three multi-locus models. Integrating these newly developed multi-locus GWAS models to unravel the genetic architecture of kernel traits, the mrMLM approach detected the maximum number of SNPs. Furthermore, a total of 41 putative candidate genes were predicted to likely be involved in the genetic architecture underlining kernel traits. These findings can facilitate a better understanding of the complex genetic mechanisms of kernel traits and may lead to the genetic improvement of grain yield in wheat.

Diverse Wheat-Alien Introgression Lines as a Basis for Durable Resistance and Quality Characteristics in Bread Wheat

Wheat productivity has been significantly improved worldwide through the incorporation of novel genes from various gene pools, not least from wild relatives of wheat, into the commonly cultivated bread and durum wheat. Here, we present and summarize results obtained from a diverse set of wheat-alien introgression lines with mainly introgressions of rye, but also of Leymus spp. and Thinopyrum junceiforme into bread-wheat (Triticum aestivum L.). From this material, lines carrying 2RL were found with good agronomic performance and multiple resistance not least towards several races of powdery mildew. A novel resistance gene, one of few showing resistance towards all today identified stem rust races, designated Sr59, was also found originating from 2RL. Lines with multiple introgressions from 4R, 5R, and 6R were found resistant towards the majority of the stripe rust races known today. Due to lack of agricultural adaptation in these lines, transfer of useful genes into more adapted wheat material is a necessity, work which is also in progress through crosses with the CSph1b mutant, to be able to only transfer small chromosome segments that carry the target gene. Furthermore, resistance towards Russian wheat aphid was found in lines having a substitution of 1R (1D) and translocations of 3DL.3RS and 5AL.5RS. The rye chromosomes 1R, 2R, and 6R were found responsible for resistance towards the Syrian Hessian fly. High levels of especially zinc was found in several lines obtained from crosses with Leymus racemosus and Leymus mollis, while also some lines with 1R, 2R, or 5R showed increased levels of minerals and in particular of iron and zinc. Moreover, lines with 1R, 2R, 3R, and Leymus spp. introgressions were also found to have a combination of high iron and zinc and low cadmium concentrations. High variation was found both in grain protein concentration and gluten strength, measured as %UPP, within the lines, indicating large variation in bread-making quality. Thus, our study emphasizes the impact that wheat-alien introgression lines can contribute to current wheat lines and shows large opportunities both to improve production, resistance, and quality. To obtain such improvements, novel plant breeding tools, as discussed in this paper, opens unique opportunities, to transfer suitable genes into the modern and adapted wheat cultivars.

Tool for genomic selection and breeding to evolutionary adaptation: Development of a 100K single nucleotide polymorphism array for the honey bee

High-throughput high-density genotyping arrays continue to be a fast, accurate, and cost-effective method for genotyping thousands of polymorphisms in high numbers of individuals. Here, we have developed a new high-density SNP genotyping array (103,270 SNPs) for honey bees, one of the most ecologically and economically important pollinators worldwide. SNPs were detected by conducting whole-genome resequencing of 61 honey bee drones (haploid males) from throughout Europe. Selection of SNPs for the chip was done in multiple steps using several criteria. The majority of SNPs were selected based on their location within known candidate regions or genes underlying a range of honey bee traits, including hygienic behavior against pathogens, foraging, and subspecies. Additionally, markers from a GWAS of hygienic behavior against the major honey bee parasite Varroa destructor were brought over. The chip also includes SNPs associated with each of three major breeding objectives-honey yield, gentleness, and Varroa resistance. We validated the chip and make recommendations for its use by determining error rates in repeat genotypings, examining the genotyping performance of different tissues, and by testing how well different sample types represent the queen's genotype. The latter is a key test because it is highly beneficial to be able to determine the queen's genotype by nonlethal means. The array is now publicly available and we suggest it will be a useful tool in genomic selection and honey bee breeding, as well as for GWAS of different traits, and for population genomic, adaptation, and conservation questions.

Exploiting the genome of Thinopyrum elongatum to expand the gene pool of hexaploid wheat

One hundred and thirty four introgressions from Thinopyrum elongatum have been transferred into a wheat background and were characterised using 263 SNP markers. Species within the genus Thinopyrum have been shown to carry genetic variation for a very wide range of traits including biotic and abiotic stresses and quality. Research has shown that one of the species within this genus, Th. elongatum, has a close relationship with the genomes of wheat making it a highly suitable candidate to expand the gene pool of wheat. Homoeologous recombination, in the absence of the Ph1 gene, has been exploited to transfer an estimated 134 introgressions from Th. elongatum into a hexaploid wheat background. The introgressions were detected and characterised using 263 single nucleotide polymorphism markers from a 35 K Axiom^® Wheat-Relative Genotyping Array, spread across seven linkage groups and validated using genomic in situ hybridisation. The genetic map had a total length of 187.8 cM and the average chromosome length was 26.8 cM. Comparative analyses of the genetic map of Th. elongatum and the physical map of hexaploid wheat confirmed previous work that indicated good synteny at the macro-level, although Th. elongatum does not contain the 4A/5A/7B translocation found in wheat.

Role of NGS and SNP genotyping methods in sugarcane improvement programs

Sugarcane (Saccharum spp.) is one of the most economically significant crops because of its high sucrose content and it is a promising biomass feedstock for biofuel production. Sugarcane genome sequencing and analysis is a difficult task due to its heterozygosity and polyploidy. Long sequence read technologies, PacBio Single-Molecule Real-Time (SMRT) sequencing, the Illumina TruSeq, and the Oxford Nanopore sequencing could solve the problem of genome assembly. On the applications side, next generation sequencing (NGS) technologies played a major role in the discovery of single nucleotide polymorphism (SNP) and the development of low to high throughput genotyping platforms. The two mainstream high throughput genotyping platforms are the SNP microarray and genotyping by sequencing (GBS). This paper reviews the NGS in sugarcane genomics, genotyping methodologies, and the choice of these methods. Array-based SNP genotyping is robust, provides consistent SNPs, and relatively easier downstream data analysis. The GBS method identifies large scale SNPs across the germplasm. A combination of targeted GBS and array-based genotyping methods should be used to increase the accuracy of genomic selection and marker-assisted breeding.

Recommendations for Choosing the Genotyping Method and Best Practices for Quality Control in Crop Genome-Wide Association Studies

High-throughput genotyping boosts genome-wide association studies (GWAS) in crop species, leading to the identification of single-nucleotide polymorphisms (SNPs) associated with economically important traits. Choosing a cost-effective genotyping method for crop GWAS requires careful examination of several aspects, namely, the purpose and the scale of the study, crop-specific genomic features, and technical and economic matters associated with each genotyping option. Once genotypic data have been obtained, quality control (QC) procedures must be applied to avoid bias and false signals in genotype–phenotype association tests. QC for human GWAS has been extensively reviewed; however, QC for crop GWAS may require different actions, depending on the GWAS population type. Here, we review most popular genotyping methods based on next-generation sequencing (NGS) and array hybridization, and report observations that should guide the investigator in the choice of the genotyping method for crop GWAS. We provide recommendations to perform QC in crop species, and deliver an overview of bioinformatics tools that can be used to accomplish all needed tasks. Overall, this work aims to provide guidelines to harmonize those procedures leading to SNP datasets ready for crop GWAS.

Genome-Wide Association Studies in Diverse Spring Wheat Panel for Stripe, Stem, and Leaf Rust Resistance

Among several important wheat foliar diseases, Stripe rust (YR), Leaf rust (LR), and Stem rust (SR) have always been an issue of concern to the farmers and wheat breeders. Evolution of virulent pathotypes of these rusts has posed frequent threats to an epidemic. Pyramiding rust-resistant genes are the most economical and environment-friendly approach in postponing this inevitable threat. To achieve durable long term resistance against the three rusts, an attempt in this study was made searching for novel sources of resistant alleles in a panel of 483 spring wheat genotypes. This is a unique and comprehensive study where evaluation of a diverse panel comprising wheat germplasm from various categories and adapted to different wheat agro-climatic zones was challenged with 18 pathotypes of the three rusts with simultaneous screening in field conditions. The panel was genotyped using 35K SNP array and evaluated for each rust at two locations for two consecutive crop seasons. High heritability estimates of disease response were observed between environments for each rust type. A significant effect of population structure in the panel was visible in the disease response. Using a compressed mixed linear model approach, 25 genomic regions were found associated with resistance for at least two rusts. Out of these, seven were associated with all the three rusts on chromosome groups 1 and 6 along with 2B. For resistance against YR, LR, and SR, there were 16, 18, and 27 QTL (quantitative trait loci) identified respectively, associated at least in two out of four environments. Several of these regions got annotated with resistance associated genes viz. NB-LRR, E3-ubiquitin protein ligase, ABC transporter protein, etc. Alien introgressed (on 1B and 3D) and pleiotropic (on 7D) resistance genes were captured in seedling and adult plant disease responses, respectively. The present study demonstrates the use of genome-wide association for identification of a large number of favorable alleles for leaf, stripe, and stem rust resistance for broadening the genetic base. Quick conversion of these QTL into user-friendly markers will accelerate the deployment of these resistance loci in wheat breeding programs.

Association Analysis Identifies New Loci for Resistance to Chinese Yr26-Virulent Races of the Stripe Rust Pathogen in a Diverse Panel of Wheat Germplasm

Stripe rust caused by Puccinia striiformis f. sp. tritici (Pst) is one of the most destructive fungal diseases of wheat worldwide. The expanding Yr26-virulent Pst race (V26) group overcomes almost all currently deployed resistance genes in China and has continued to accumulate new virulence. Investigating the genetic architecture of stripe rust resistance in common wheat is an important basis for a successful utilization of resistance in breeding programs. A panel of 410 exotic wheat germplasms was used for characterizing new stripe rust resistance loci. This panel was genotyped using high-density wheat 660K single-nucleotide polymorphism (SNP) array, and phenotypic evaluation of seedlings for stripe rust resistance was performed using multiple Pst races. Thirty-five loci conferring resistance were identified through genome-wide association mapping, and explained phenotypic variances ranged from 53 to 75%. Of these, 14 were colocated in the proximity of the known loci, including cataloged Yr genes Yr9, Yr10, Yr26, Yr33, Yr47, Yr56, Yr57, Yr64, Yr67, Yr72, and Yr81 and three temporarily designated as YrCen, YrNP63, and YrRC detected in our quantitative trait locus (QTL) mapping studies. Seven of them (Yr9, Yr10, Yr24/26, Yr81, YrCEN, YrNP63, and YrRC) were confirmed by molecular detection or genetic analysis. New loci that were identified to be different from reported Yr genes need further confirmation. Nine QTL with significantly large phenotypic effect on resistance to all tested races were considered as major loci for effective resistance. The identified loci enrich our stripe rust resistance gene pool, and the linked SNPs should be useful for marker-assisted selection in breeding programs.

The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat

The rapid development and application of molecular marker assays have facilitated genomic selection and genome-wide linkage and association studies in wheat breeding. Although PCR-based markers (e.g. simple sequence repeats and functional markers) and genotyping by sequencing have contributed greatly to gene discovery and marker-assisted selection, the release of a more accurate and complete bread wheat reference genome has resulted in the design of single-nucleotide polymorphism (SNP) arrays based on different densities or application targets. Here, we evaluated seven types of wheat SNP arrays in terms of their SNP number, distribution, density, associated genes, heterozygosity and application. The results suggested that the Wheat 660K SNP array contained the highest percentage (99.05%) of genome-specific SNPs with reliable physical positions. SNP density analysis indicated that the SNPs were almost evenly distributed across the whole genome. In addition, 229 266 SNPs in the Wheat 660K SNP array were located in 66 834 annotated gene or promoter intervals. The annotated genes revealed by the Wheat 660K SNP array almost covered all genes revealed by the Wheat 35K (97.44%), 55K (99.73%), 90K (86.9%) and 820K (85.3%) SNP arrays. Therefore, the Wheat 660K SNP array could act as a substitute for other 6 arrays and shows promise for a wide range of possible applications. In summary, the Wheat 660K SNP array is reliable and cost-effective and may be the best choice for targeted genotyping and marker-assisted selection in wheat genetic improvement.

Mapping Quantitative Trait Loci for 1000-Grain Weight in a Double Haploid Population of Common Wheat

Thousand-grain weight (TGW) is a very important yield trait of crops. In the present study, we performed quantitative trait locus (QTL) analysis of TGW in a doubled haploid population obtained from a cross between the bread wheat cultivar "Superb" and the breeding line "M321" using the wheat 55-k single-nucleotide polymorphism (SNP) genotyping assay. A genetic map containing 15,001 SNP markers spanning 2209.64 cM was constructed, and 9 QTLs were mapped to chromosomes 1A, 2D, 4B, 4D, 5A, 5D, 6A, and 6D based on analyses conducted in six experimental environments during 2015-2017. The effects of the QTLs qTgw.nwipb-4DS and qTgw.nwipb-6AL were shown to be strong and stable in different environments, explaining 15.31-32.43% and 21.34-29.46% of the observed phenotypic variance, and they were mapped within genetic distances of 2.609 cM and 5.256 cM, respectively. These novel QTLs may be used in marker-assisted selection in wheat high-yield breeding.

High-Resolution Mapping of the Novel Early Leaf Senescence Gene Els2 in Common Wheat

Early leaf senescence negatively impacts the grain yield in wheat (Triticum aestivum L.). Induced mutants provide an important resource for mapping and cloning of genes for early leaf senescence. In our previous study, Els2, a single incomplete dominance gene, that caused early leaf senescence phenotype in the wheat mutant LF2099, had been mapped on the long arm of chromosome 2B. The objective of this study was to develop molecular markers tightly linked to the Els2 gene and construct a high-resolution map surrounding the Els2 gene. Three tightly linked single-nucleotide polymorphism (SNP) markers were obtained from the Illumina Wheat 90K iSelect SNP genotyping array and converted to Kompetitive allele-specific polymerase chain reaction (KASP) markers. To saturate the Els2 region, the Axiom® Wheat 660K SNP array was used to screen bulked extreme phenotype DNA pools, and 9 KASP markers were developed. For fine mapping of the Els2 gene, these KASP markers and previously identified polymorphic markers were analyzed in a large F₂ population of the LF2099 × Chinese Spring cross. The Els2 gene was located in a 0.24-cM genetic region flanked by the KASP markers AX-111643885 and AX-111128667, which corresponded to a physical interval of 1.61 Mb in the Chinese Spring chromosome 2BL containing 27 predicted genes with high confidence. The study laid a foundation for a map-based clone of the Els2 gene controlling the mutation phenotype and revealing the molecular regulatory mechanism of wheat leaf senescence.

Combining QTL Analysis and Genomic Predictions for Four Durum Wheat Populations Under Drought Conditions

Durum wheat is an important crop for the human diet and its consumption is gaining popularity. In order to ensure that durum wheat production maintains the pace with the increase in demand, it is necessary to raise productivity by approximately 1.5% per year. To deliver this level of annual genetic gain the incorporation of molecular strategies has been proposed as a key solution. Here, four RILs populations were used to conduct QTL discovery for grain yield (GY) and 1,000 kernel weight (TKW). A total of 576 individuals were sown at three locations in Morocco and one in Lebanon. These individuals were genotyped by sequencing with 3,202 high-confidence polymorphic markers, to derive a consensus genetic map of 2,705.7 cM, which was used to impute any missing data. Six QTLs were found to be associated with GY and independent from flowering time on chromosomes 2B, 4A, 5B, 7A and 7B, explaining a phenotypic variation (PV) ranging from 4.3 to 13.4%. The same populations were used to train genomic prediction models incorporating the relationship matrix, the genotype by environment interaction, and marker by environment interaction, to reveal significant advantages for models incorporating the marker effect. Using training populations (TP) in full sibs relationships with the validation population (VP) was shown to be the only effective strategy, with accuracies reaching 0.35–0.47 for GY. Reducing the number of markers to 10% of the whole set, and the TP size to 20% resulted in non-significant changes in accuracies. The QTLs identified were also incorporated in the models as fixed effects, showing significant accuracy gain for all four populations. Our results confirm that the prediction accuracy depends considerably on the relatedness between TP and VP, but not on the number of markers and size of TP used. Furthermore, feeding the model with information on markers associated with QTLs increased the overall accuracy.

Identification of candidate chromosome region of Sbwm1 for Soil-borne wheat mosaic virus resistance in wheat

Soil-borne wheat mosaic virus (SBWMV) causes a serious viral disease that can significantly reduce grain yield in winter wheat worldwide. Using resistant cultivars is the only feasible strategy to reduce the losses caused by SBWMV. To fine map the resistance gene Sbwm1, 205 wheat accessions was genotyped using wheat Infinium iSelect Beadchips with 90 K SNPs. Association analysis identified 35 SNPs in 12 wheat genes and one intergenic SNP in the Sbwm1 region that showed a significant association with SBWMV resistance. Those SNPs were converted into Kompetitive Allele-Specific Polymerase assays (KASP) and analyzed in two F₆-derived recombinant inbred line (RIL) populations derived from the crosses between two resistant cultivars 'Wesley' and 'Deliver' and a susceptible line 'OK03825-5403-6'. Linkage analysis mapped this gene on chromosome 5D at intervals of 5.1 cM and 3.4 cM in the two populations, respectively. The two flanking markers in both populations delimited the gene to a 620 kb region where 19 genes were annotated. Comparative analysis identified a syntenic region of 660 kb in Ae. tauschii with 18 annotated genes and a syntenic region in chromosome 1 of B. distachyon. The candidate region includes several disease resistance related genes and we identified a PTI1-like tyrosine-protein kinase 1 gene as a putative candidate gene for Sbwm1. The two flanking SNPs for Sbwm1 can effectively separate the resistant and susceptible lines in a new diversity panel of 159 wheat germplasm. The results from this study lay a solid foundation for the cloning, functional characterization and marker-assisted selection of Sbwm1.

Bymovirus-induced yellow mosaic diseases in barley and wheat: viruses, genetic resistances and functional aspects

Bymovirus-induced yellow mosaic diseases seriously threaten global production of autumn-sown barley and wheat, which are two of the presently most important crops around the world. Under natural field conditions, the diseases are caused by infection of soil-borne plasmodiophorid Polymyxa graminis-transmitted bymoviruses of the genus Bymovirus of the family Potyviridae. Focusing on barley and wheat, this article summarizes the achievements on taxonomy, geography and host specificity of these disease-conferring viruses, as well as the genetics of resistance in barley, wheat and wild relatives. Moreover, based on recent progress of barley and wheat genomics, germplasm resources and large-scale sequencing, the exploration and isolation of corresponding resistant genes from wheat and barley as well as relatives, no matter what a large and complicated genome is present, are becoming feasible and are discussed. Furthermore, the foreseen advances on cloning of the resistance or susceptibility-encoding genes, which will provide the possibility to explore the functional interaction between host plants and soil-borne viral pathogens, are discussed as well as the benefits for marker-assisted resistance breeding in barley and wheat.

A roadmap for gene functional characterisation in crops with large genomes: Lessons from polyploid wheat

Understanding the function of genes within staple crops will accelerate crop improvement by allowing targeted breeding approaches. Despite their importance, a lack of genomic information and resources has hindered the functional characterisation of genes in major crops. The recent release of high-quality reference sequences for these crops underpins a suite of genetic and genomic resources that support basic research and breeding. For wheat, these include gene model annotations, expression atlases and gene networks that provide information about putative function. Sequenced mutant populations, improved transformation protocols and structured natural populations provide rapid methods to study gene function directly. We highlight a case study exemplifying how to integrate these resources. This review provides a helpful guide for plant scientists, especially those expanding into crop research, to capitalise on the discoveries made in Arabidopsis and other plants. This will accelerate the improvement of crops of vital importance for food and nutrition security.

A 62K genic-SNP chip array for genetic studies and breeding applications in pigeonpea (Cajanus cajan L. Millsp.)

Pigeonpea is the second most important pulse legume crop for food and nutritional security of South Asia that requires accelerated breeding using high throughput genomic tools. Single nucleotide polymorphisms (SNPs) are highly suitable markers for this purpose because of their bi-allelic nature, reproducibility and high abundance in the genome. Here we report on development and use of a pigeonpea 62 K SNP chip array ‘CcSNPnks’ for Affymetrix GeneTitan^® platform. The array was designed after filtering 645,662 genic-SNPs identified by re-sequencing of 45 diverse genotypes and has 62,053 SNPs from 9629 genes belonging to five different categories, including 4314 single-copy genes unique to pigeonpea, 4328 single-copy genes conserved between soybean and pigeonpea, 156 homologs of agronomically important cloned genes, 746 disease resistance and defense response genes and 85 multi-copy genes of pigeonpea. This fully genic chip has 28.94% exonic, 33.04% intronic, 27.56% 5′UTR and 10.46% 3′UTR SNPs and incorporates multiple SNPs per gene allowing gene haplotype network analysis. It was used successfully for the analysis of genetic diversity and population structure of 95 pigeonpea varieties and high resolution mapping of 11 yield related QTLs for number of branches, pod bearing length and number of seeds per pod in a biparental RIL population. As an accurate high-density genotyping tool, ‘CcSNPnks’ chip array will be useful for high resolution fingerprinting, QTL mapping and genome wide as well as gene-based association studies in pigeonpea.

The Resurgence of Introgression Breeding, as Exemplified in Wheat Improvement

Breeding progress in most crops has relied heavily on the exploitation of variation within the species' primary gene pool, a process which is destined to fail once the supply of novel variants has been exhausted. Accessing a crop's secondary gene pool, as represented by its wild relatives, has the potential to greatly expand the supply of usable genetic variation. The crop in which this approach has been most strongly championed is bread wheat (Triticum aestivum), a species which is particularly tolerant of the introduction of chromosomal segments of exotic origin thanks to the genetic buffering afforded by its polyploid status. While the process of introgression can be in itself cumbersome, a larger problem is that linkage drag and/or imperfect complementation frequently impose a yield and/or quality penalty, which explains the reluctance of breeders to introduce such materials into their breeding populations. Thanks to the development of novel strategies to induce introgression and of genomic tools to facilitate the selection of desirable genotypes, introgression breeding is returning as a mainstream activity, at least in wheat. Accessing variation present in progenitor species has even been able to drive genetic advance in grain yield. The current resurgence of interest in introgression breeding can be expected to result in an increased deployment of exotic genes in commercial wheat cultivars.

Scale development and utilization of universal PCR-based and high-throughput KASP markers specific for chromosome arms of rye (Secale cereale L.)

Background

Rye (Secale cereale L., 2n = 2x = 14, RR), a relative of common wheat, is a large gene resource pool for wheat improvement. Accurate and convenient identification of the rye chromatin in wheat background will facilitate the transfer and utilization of elite genes derived from rye in wheat breeding.

Results

In the present study, five rye cultivars including Imperial, German White, Jingzhouheimai, Baili and Guyuan were sequenced by specific-locus amplified fragment sequencing (SLAF-seq) to develop large-scale rye-specific markers. Based on SLAF-seq and bioinformatics analyses, a total of 404 universal PCR-based and a whole set of Kompetitive allele-specific PCR (KASP) markers specific for the 14 individual rye chromosome arms were developed and validated. Additionally, two KASP markers specific for 1RS and 2RL were successfully applied in the detection of 1RS translocations in a natural population and 2RL chromosome arms in wheat-rye derived progenies that conferred adult resistance to powdery mildew.

Conclusion

The 404 PCR-based markers and 14 KASP markers specific for the 14 individual rye chromosome arms developed in this study can enrich the marker densities for gene mapping and accelerate the utilization of rye-derived genes in wheat improvement. Especially, the KASP markers achieved high-throughput and accurate detection of rye chromatin in wheat background, thus can be efficiently used in marker-assisted selection (MAS). Besides, the strategy of rye-specific PCR-based markers converting into KASP markers was high-efficient and low-cost, which will facilitate the tracing of alien genes, and can also be referred for other wheat relatives.

Incorporating Genome-Wide Association Mapping Results Into Genomic Prediction Models for Grain Yield and Yield Stability in CIMMYT Spring Bread Wheat

Untangling the genetic architecture of grain yield (GY) and yield stability is an important determining factor to optimize genomics-assisted selection strategies in wheat. We conducted in-depth investigation on the above using a large set of advanced bread wheat lines (4,302), which were genotyped with genotyping-by-sequencing markers and phenotyped under contrasting (irrigated and stress) environments. Haplotypes-based genome-wide-association study (GWAS) identified 58 associations with GY and 15 with superiority index Pi (measure of stability). Sixteen associations with GY were "environment-specific" with two on chromosomes 3B and 6B with the large effects and 8 associations were consistent across environments and trials. For Pi, 8 associations were from chromosomes 4B and 7B, indicating 'hot spot' regions for stability. Epistatic interactions contributed to an additional 5-9% variation on average. We further explored whether integrating consistent and robust associations identified in GWAS as fixed effects in prediction models improves prediction accuracy. For GY, the model accounting for the haplotype-based GWAS loci as fixed effects led to up to 9-10% increase in prediction accuracy, whereas for Pi this approach did not provide any advantage. This is the first report of integrating genetic architecture of GY and yield stability into prediction models in wheat.

Dmc1 is a candidate for temperature tolerance during wheat meiosis

KEY MESSAGE

The meiotic recombination gene Dmc1 on wheat chromosome 5D has been identified as a candidate for the maintenance of normal chromosome synapsis and crossover at low and possibly high temperatures. We initially assessed the effects of low temperature on meiotic chromosome synapsis and crossover formation in the hexaploid wheat (Triticum aestivum L.) variety 'Chinese Spring'. At low temperatures, asynapsis and chromosome univalence have been observed before in Chinese Spring lines lacking the long arm of chromosome 5D (5DL), which led to the proposal that 5DL carries a gene (Ltp1) that stabilises wheat chromosome pairing at low temperatures. In the current study, Chinese Spring wild type and 5DL interstitial deletion mutant plants were exposed to low temperature in a controlled environment room during a period from premeiotic interphase to early meiosis I. A 5DL deletion mutant was identified whose meiotic chromosomes exhibit extremely high levels of asynapsis and chromosome univalence at metaphase I after 7 days at 13 °C, suggesting that Ltp1 is deleted in this mutant. Immunolocalisation of the meiotic proteins ASY1 and ZYP1 on ltp1 mutants showed that low temperature results in a failure to complete synapsis at pachytene. KASP genotyping revealed that the ltp1 mutant has a 4-Mb deletion in 5DL. Of 41 genes within this deletion region, the strongest candidate for the stabilisation of chromosome pairing at low temperatures is the meiotic recombination gene Dmc1. The ltp1 mutants were subsequently treated at 30 °C for 24 h during meiosis and exhibited a reduced number of crossovers and increased univalence, though to a lesser extent than at 13 °C. We therefore renamed our ltp1 mutant 'ttmei1' (temperature-tolerant meiosis 1) to reflect this additional loss of high temperature tolerance.

Aegilops umbellulata introgression carrying leaf rust and stripe rust resistance genes Lr76 and Yr70 located to 9.47-Mb region on 5DS telomeric end through a combination of chromosome sorting and sequencing

Lr76 and Yr70 have been fine mapped using the sequence of flow-sorted recombinant 5D chromosome from wheat-Ae. umbellulata introgression line. The alien introgression has been delineated to 9.47-Mb region on short arm of wheat chromosome 5D. Leaf rust and stripe rust are among the most damaging diseases of wheat worldwide. Wheat cultivation based on limited number of rust resistance genes deployed over vast areas expedites the emergence of new pathotypes warranting a continuous deployment of new resistance genes. In this paper, fine mapping of Aegilops umbellulata-derived leaf rust and stripe rust resistance genes Lr76 and Yr70 is being reported. We flow sorted and paired-end sequenced 5U chromosome of Ae. umbellulata, recombinant chromosome 5D (5D^IL) from wheat-Ae. umbellulata introgression line pau16057 and 5D^RP of recurrent parent WL711. Chromosome 5U reads were mapped against the reference Chinese Spring chromosome 5D sequence, and alien-specific SNPs were identified. Chromosome 5D^IL and 5D^RP sequences were de novo assembled, and alien introgression-specific markers were designed by selecting 5U- and 5D-specific SNPs. Overall, 27 KASP markers were mapped in high-resolution population consisting of 1404 F₅ RILs. The mapping population segregated for single gene each for leaf rust and stripe rust resistance. The physical order of the SNPs in pau16057 was defined by projecting the 27 SNPs against the IWGSC RefSeq v1.0 sequence. Based on this physical map, the size of Ae. umbellulata introgression was determined to be 9.47 Mb on the distal most end of the short arm of chromosome 5D. This non-recombining alien segment carries six NB-LRR encoding genes based on NLR annotation of assembled chromosome 5D^IL sequence and IWGSC RefSeq v1.1 gene models. The presence of SNPs and other sequence variations in these genes between pau16057 and WL711 suggested that they are candidates for Lr76 and Yr70.

Yield reduction historically associated with the Aegilops ventricosa 7DV introgression is genetically and physically distinct from the eyespot resistance gene Pch1

Yield penalty and increased grain protein content traits associated with Aegilops ventricosa 7D introgression have been mapped for the first time, and they are physically distinct from the eyespot resistance locus Pch1. Wheat wild relatives represent an important source of genetic variation, but introgression of agronomically relevant genes, such as for disease resistance, may lead to the simultaneous introduction of genetically linked deleterious traits. Pch1 is a dominant gene, conferring resistance to eyespot and was introgressed to wheat from Aegilops ventricosa as part of a large segment of the 7D^V chromosome. This introgression has been associated with a significant yield reduction and a concomitant increase in grain protein content. In this study, we evaluated both traits and their relationship to the location of the Pch1 gene. We found that both QTLs were clearly distinct from the Pch1 gene, being located on a different linkage group to Pch1. In addition, we found that the QTL for increased grain protein content was strong and consistent across field trials, whereas the yield penalty QTL was unstable and environmentally dependent. The yield and grain protein content QTLs were genetically linked and located in the same linkage group. This finding is due in part to the small size of the population, and to the restricted recombination between wheat 7D and Ae. ventricosa 7D^v chromosomes. Although recombination in this interval is rare, it does occur. A recombinant line containing Pch1 and 7D_KASP6, the marker associated with increase in grain protein content, but not Xwmc221, the marker associated with the yield penalty effect, was identified.

Rapid identification of homozygosity and site of wild relative introgressions in wheat through chromosome-specific KASP genotyping assays

For future food security, it is important that wheat, one of the most widely consumed crops in the world, can survive the threat of abiotic and biotic stresses. New genetic variation is currently being introduced into wheat through introgressions from its wild relatives. For trait discovery, it is necessary that each introgression is homozygous and hence stable. Breeding programmes rely on efficient genotyping platforms for marker-assisted selection (MAS). Recently, single nucleotide polymorphism (SNP)-based markers have been made available on high-throughput Axiom^® SNP genotyping arrays. However, these arrays are inflexible in their design and sample numbers, making their use unsuitable for long-term MAS. SNPs can potentially be converted into Kompetitive allele-specific PCR (KASP™) assays that are comparatively cost-effective and efficient for low-density genotyping of introgression lines. However, due to the polyploid nature of wheat, KASP assays for homoeologous SNPs can have difficulty in distinguishing between heterozygous and homozygous hybrid lines in a backcross population. To identify co-dominant SNPs, that can differentiate between heterozygotes and homozygotes, we PCR-amplified and sequenced genomic DNA from potential single-copy regions of the wheat genome and compared them to orthologous copies from different wild relatives. A panel of 620 chromosome-specific KASP assays have been developed that allow rapid detection of wild relative segments and provide information on their homozygosity and site of introgression in the wheat genome. A set of 90 chromosome-nonspecific assays was also produced that can be used for genotyping introgression lines. These multipurpose KASP assays represent a powerful tool for wheat breeders worldwide.

Development, identification, and characterization of blue-grained wheat-Triticum boeoticum substitution lines

Diploid wild einkorn wheat, Triticum boeoticum Boiss (A^bA^b, 2n = 2x = 14), is a wheat-related species with a blue aleurone layer. In this study, six blue-grained wheat lines were developed from F₈ progeny of crosses between common wheat and T. boeoticum. The chromosome constitutions of these lines were characterized by fluorescence in situ hybridization (FISH) using the oligonucleotide probes Oligo-pTa535-1, Oligo-pSc119.2-1, Oligo-pTa71-2, and (AAC)₇. Multicolor FISH using Oligo-pTa535-1, Oligo-pSc119.2-1, and Oligo-pTa71-2 identified all 42 common wheat chromosomes, while Oligo-pTa535-1 and (AAC)₇ discriminated the 14 chromosomes of T. boeoticum. FISH revealed that all six blue-grained lines were wheat-T. boeoticum 4A^b (4B) disomic substitution lines. The substitution lines were validated by genotyping using the wheat 55 K single nucleotide polymorphism (SNP) array containing 53,063 markers. These 4A^b (4B) substitution lines represent novel germplasm for blue-grained wheat breeding. The FISH probes and SNP markers used here can be applied in the development of blue-grained wheat-Triticum boeoticum translocation lines.

Suitability of Single-Nucleotide Polymorphism Arrays Versus Genotyping-By-Sequencing for Genebank Genomics in Wheat

Genebank genomics promises to unlock valuable diversity for plant breeding but first, one key question is which marker system is most suitable to fingerprint entire genebank collections. Using wheat as model species, we tested for the presence of an ascertainment bias and investigated its impact on estimates of genetic diversity and prediction ability obtained using three marker platforms: simple sequence repeat (SSR), genotyping-by-sequencing (GBS), and array-based SNP markers. We used a panel of 378 winter wheat genotypes including 190 elite lines and 188 plant genetic resources (PGR), which were phenotyped in multi-environmental trials for grain yield and plant height. We observed an ascertainment bias for the array-based SNP markers, which led to an underestimation of the molecular diversity within the population of PGR. In contrast, the marker system played only a minor role for the overall picture of the population structure and precision of genome-wide predictions. Interestingly, we found that rare markers contributed substantially to the prediction ability. This combined with the expectation that valuable novel diversity is most likely rare suggests that markers with minor allele frequency deserve careful consideration in the design of a pre-breeding program.

Genetic diversity and population structure analysis based on the high density SNP markers in Ethiopian durum wheat (Triticum turgidum ssp. durum)

Background

Ethiopia has been considered as a center of diversity and the second possible center of domestication of durum wheat. Genetic diversity and population structure analysis in the existing Ethiopian durum wheat germplasm have enormous importance in enhancing breeding effort and for sustainable conservation. Hence, 192 Ethiopian durum wheat accessions comprising 167 landraces collected from major wheat-growing areas of the country and 25 improved varieties released from Debre Zeit and Sinana Agricultural Research Centers, Ethiopia in different years (1994–2010) were assembled for the current study.

Results

The panel was genotyped with a High-density 90 K wheat SNP array by Illumina and generated 15,338 polymorphic SNPs that were used to analyze the genetic diversity and to estimate the population structure. Varied values of genetic diversity indices were scored across chromosomes and genomes. Genome-wide mean values of Nei’s gene diversity (0.246) and polymorphism information content (0.203) were recorded signifying the presence of high genetic diversity within this collection. Minor allele frequency of the genome varied with a range of 0.005 to 0.5 scoring a mean value of 0.175. Improved varieties clustered separately to landraces in population structure analysis resulted from STRUCTURE, PCA and neighbor joining tree. Landraces clustering was irrespective of their geographical origin signifying the presence of higher admixture that could arise due to the existence of historical exchanges of seeds through informal seed system involving regional and countrywide farming communities in Ethiopia.

Conclusions

Sustainable utilization and conservation of this rich Ethiopian durum wheat genetic resource is an irreplaceable means to cope up from the recurrent climate changes and biotic stresses happening widely and thereby able to keep meeting the demand of durum productivity for the ever-growing human population.

Molecular Mapping of Loci Conferring Susceptibility to Spot Blotch and Resistance to Powdery Mildew in Barley Using the Sequencing-Based Genotyping Approach

Spot blotch (SB) caused by Bipolaris sorokiniana and powdery mildew (PM) caused by Blumeria graminis f. sp. hordei are two important diseases of barley. To map genetic loci controlling susceptibility and resistance to these diseases, a mapping population consisting of 138 recombinant inbred lines (RILs) was developed from the cross between Bowman and ND5883. A genetic map was constructed for the population with 852 unique single nucleotide polymorphism markers generated by sequencing-based genotyping. Bowman and ND5883 showed distinct infection responses at the seedling stage to two isolates (ND90Pr and ND85F) of Bipolaris sorokiniana and one isolate (Race I) of Blumeria graminis f. sp. hordei. Genetic analysis of the RILs revealed that one major gene (Scs6) controls susceptibility to Bipolaris sorokiniana isolate ND90Pr, and another major gene (Mla8) confers resistance to Blumeria graminis f. sp. hordei isolate Race I, respectively. Scs6 was mapped on chromosome 1H of Bowman, as previously reported. Mla8 was also mapped to the short arm of 1H, which was tightly linked but not allelic to the Rcs6/Scs6 locus. Quantitative trait locus (QTL) analysis identified two QTLs, QSbs-1H-P1 and QSbs-7H-P1, responsible for susceptibility to spot blotch caused by Bipolaris sorokiniana isolate ND85F in ND5883, which are located on chromosome 1H and 7H, respectively. QSbs-7H-P1 was mapped to the same region as Rcs5, whereas QSbs-1H-P1 may represent a novel allele conferring seedling stage susceptibility to isolate ND85F. Identification and molecular mapping of the loci for SB susceptibility and PM resistance will facilitate development of barley cultivars with resistance to the diseases.

An Axiom SNP genotyping array for Douglas-fir

BACKGROUND

In forest trees, genetic markers have been used to understand the genetic architecture of natural populations, identify quantitative trait loci, infer gene function, and enhance tree breeding. Recently, new, efficient technologies for genotyping thousands to millions of single nucleotide polymorphisms (SNPs) have finally made large-scale use of genetic markers widely available. These methods will be exceedingly valuable for improving tree breeding and understanding the ecological genetics of Douglas-fir, one of the most economically and ecologically important trees in the world.

RESULTS

We designed SNP assays for 55,766 potential SNPs that were discovered from previous transcriptome sequencing projects. We tested the array on ~ 2300 related and unrelated coastal Douglas-fir trees (Pseudotsuga menziesii var. menziesii) from Oregon and Washington, and 13 trees of interior Douglas-fir (P. menziesii var. glauca). As many as ~ 28 K SNPs were reliably genotyped and polymorphic, depending on the selected SNP call rate. To increase the number of SNPs and improve genome coverage, we developed protocols to 'rescue' SNPs that did not pass the default Affymetrix quality control criteria (e.g., 97% SNP call rate). Lowering the SNP call rate threshold from 97 to 60% increased the number of successful SNPs from 20,669 to 28,094. We used a subset of 395 unrelated trees to calculate SNP population genetic statistics for coastal Douglas-fir. Over a range of call rate thresholds (97 to 60%), the median call rate for SNPs in Hardy-Weinberg equilibrium ranged from 99.2 to 99.7%, and the median minor allele frequency ranged from 0.198 to 0.233. The successful SNPs also worked well on interior Douglas-fir.

CONCLUSIONS

Based on the original transcriptome assemblies and comparisons to version 1.0 of the Douglas-fir reference genome, we conclude that these SNPs can be used to genotype about 10 K to 15 K loci. The Axiom genotyping array will serve as an excellent foundation for studying the population genomics of Douglas-fir and for implementing genomic selection. We are currently using the array to construct a linkage map and test genomic selection in a three-generation breeding program for coastal Douglas-fir.

Examining the Effects of Temperature on Recombination in Wheat

Meiotic recombination plays a crucial role in the generation of new varieties. The effectiveness of recombination is limited by the distribution of crossover events, which in wheat and many other crops is skewed toward the distal regions of the chromosomes. Whole-genome sequencing of wheat has revealed that there are numerous important genes in the pericentromeric regions, which are inaccessible to manipulation due to the lack of crossover events. Studies in barley have shown that the distribution of recombination events can be shifted toward the centromeres by increasing temperature during meiosis. Here we present an analysis of the effects of temperature on the distribution and frequency of recombination events in wheat. Our data show that although increased temperature during meiosis does cause an inward shift in recombination distribution for some chromosomes, its overall utility is limited, with many genes remaining highly linked.

CerealsDB-new tools for the analysis of the wheat genome: update 2020

CerealsDB (www.cerealsdb.uk.net) is an online repository of mainly hexaploid wheat (Triticum aestivum) single nucleotide polymorphisms (SNPs) and genotyping data. The CerealsDB website has been designed to enable wheat breeders and scientists to select the appropriate markers for research breeding tasks, such as marker-assisted selection. We report a large update of genotyping information for over 6000 wheat accessions and describe new webtools for exploring and visualizing the data. We also describe a new database of quantitative trait loci that links phenotypic traits to CerealsDB SNP markers and allelic scores for each of those markers. CerealsDB is an open-access website that hosts information on wheat SNPs considered useful for both plant breeders and research scientists. The latest CerealsDB database is available at https://www.cerealsdb.uk.net/cerealgenomics/CerealsDB/indexNEW.php.

Amand P, Le HQ, Su Z, Bai G. Multiplex restriction amplicon sequencing: a novel next-generation sequencing-based marker platform for high-throughput genotyping

To enable rapid selection of traits in marker-assisted breeding, markers must be technically simple, low-cost, high-throughput and randomly distributed in a genome. We developed such a technology, designated as Multiplex Restriction Amplicon Sequencing (MRASeq), which reduces genome complexity by polymerase chain reaction (PCR) amplification of amplicons flanked by restriction sites. The first PCR primers contain restriction site sequences at 3'-ends, preceded by 6-10 bases of specific or degenerate nucleotide sequences and then by a unique M13-tail sequence which serves as a binding site for a second PCR that adds sequencing primers and barcodes to allow sample multiplexing for sequencing. The sequences of restriction sites and adjacent nucleotides can be altered to suit different species. Physical mapping of MRASeq SNPs from a biparental population of allohexaploid wheat (Triticum aestivum L.) showed a random distribution of SNPs across the genome. MRASeq generated thousands of SNPs from a wheat biparental population and natural populations of wheat and barley (Hordeum vulgare L.). This novel, next-generation sequencing-based genotyping platform can be used for linkage mapping to screen quantitative trait loci (QTL), background selection in breeding and many other genetics and breeding applications of various species.

Identification of New Leaf Rust Resistance Loci in Wheat and Wild Relatives by Array-Based SNP Genotyping and Association Genetics

Leaf rust caused by Puccinia triticina is the most widespread rust disease of wheat. As pathogen populations are constantly evolving, identification of novel sources of resistance is necessary to maintain disease resistance and stay ahead of this plant-pathogen evolutionary arms race. The wild genepool of wheat is a rich source of genetic diversity, accounting for 44% of the Lr genes identified. Here we performed a genome-wide association study (GWAS) on a diverse germplasm of 385 accessions, including 27 different Triticum and Aegilops species. Genetic characterization using the wheat 90 K array and subsequent filtering identified a set of 20,501 single nucleotide polymorphic (SNP) markers. Of those, 9,570 were validated using exome capture and mapped onto the Chinese Spring reference sequence v1.0. Phylogenetic analyses illustrated four major clades, clearly separating the wild species from the T. aestivum and T. turgidum species. GWAS was conducted using eight statistical models for infection types against six leaf rust isolates and leaf rust severity rated in field trials for 3-4 years at 2-3 locations in Canada. Functional annotation of genes containing significant quantitative trait nucleotides (QTNs) identified 96 disease-related loci associated with leaf rust resistance. A total of 21 QTNs were in haplotype blocks or within flanking markers of at least 16 known Lr genes. The remaining significant QTNs were considered loci that putatively harbor new Lr resistance genes. Isolation of these candidate genes will contribute to the elucidation of their role in leaf rust resistance and promote their usefulness in marker-assisted selection and introgression.

Grain protein content and thousand kernel weight QTLs identified in a durum × wild emmer wheat mapping population tested in five environments

Genetic dissection of GPC and TKW in tetraploid durum × WEW RIL population, based on high-density SNP genetic map, revealed 12 GPC QTLs and 11 TKW QTLs, with favorable alleles for 11 and 5 QTLs, respectively, derived from WEW. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, WEW) was shown to exhibit high grain protein content (GPC) and therefore possess a great potential for improvement of cultivated wheat nutritional value. Genetic dissection of thousand kernel weight (TKW) and grain protein content (GPC) was performed using a high-density genetic map constructed based on a recombinant inbred line (RIL) population derived from a cross between T. durum var. Svevo and WEW acc. Y12-3. Genotyping of 208 F₆ RILs with a 15 K wheat single nucleotide polymorphism (SNP) array yielded 4166 polymorphic SNP markers, of which 1510 were designated as skeleton markers. A total map length of 2169 cM was obtained with an average distance of 1.5 cM between SNPs. A total of 12 GPC QTLs and 11 TKW QTLs were found under five different environments. No significant correlations were found between GPC and TKW across all environments. Four major GPC QTLs with favorable alleles from WEW were found on chromosomes 4BS, 5AS, 6BS and 7BL. The 6BS GPC QTL coincided with the physical position of the NAC transcription factor TtNAM-B1, underlying the cloned QTL, Gpc-B1. Comparisons of the physical intervals of the GPC QTLs described here with the results previously reported in other durum × WEW RIL population led to the discovery of seven novel GPC QTLs. Therefore, our research emphasizes the importance of GPC QTL dissection in diverse WEW accessions as a source of novel alleles for improvement of GPC in cultivated wheat.

Dissection of genetic factors underlying grain size and fine mapping of QTgw.cau-7D in common wheat (Triticum aestivum L.)

Thirty environmentally stable QTL controlling grain size and/or plant height were identified, among which QTgw.cau-7D was delimited into the physical interval of approximately 4.4 Mb. Grain size and plant height (PHT) are important agronomic traits in wheat breeding. To dissect the genetic basis of these traits, we conducted a quantitative trait locus (QTL) analysis using recombinant inbred lines (RILs). In total, 30 environmentally stable QTL for thousand grain weight (TGW), grain length (GL), grain width (GW) and PHT were detected. Notably, one major pleiotropic QTL on chromosome arm 3DS explained the highest phenotypic variance for TGW, GL and PHT, and two stable QTL (QGw.cau-4B and QGw.cau-7D) on chromosome arms 4BS and 7DS contributed greater effects for GW. Furthermore, the stable QTL controlling grain size (QTgw.cau-7D and QGw.cau-7D) were delimited into the physical interval of approximately 4.4 Mb harboring 56 annotated genes. The elite NILs of QTgw.cau-7D increased TGW by 12.79-21.75% and GW by 4.10-8.47% across all three environments. Collectively, these results provide further insight into the genetic basis of TGW, GL, GW and PHT, and the fine-mapped QTgw.cau-7D will be an attractive target for positional cloning and marker-assisted selection in wheat breeding programs.