Diversity arrays technology (DArT) for high-throughput profiling of the hexaploid wheat genome

Novel SNP markers for flowering and seed quality traits in faba bean (Vicia faba L.): characterization and GWAS of a diversity panel

Faba bean (Vicia faba L.) is a legume crop grown in diverse climates worldwide. It has a high potential for increased cultivation to meet the need for more plant-based proteins in human diets, a prerequisite for a more sustainable food production system. Characterization of diversity panels of crops can identify variation in and genetic markers for target traits of interest for plant breeding. In this work, we collected a diversity panel of 220 accessions of faba bean from around the world consisting of gene bank material and commercially available cultivars. The aims of this study were to quantify the phenotypic diversity in target traits to analyze the impact of breeding on these traits, and to identify genetic markers associated with traits through a genome-wide association study (GWAS). Characterization under field conditions at Nordic latitude across two years revealed a large genotypic variation and high broad-sense heritability for eleven agronomic and seed quality traits. Pairwise correlations showed that seed yield was positively correlated to plant height, number of seeds per plant, and days to maturity. Further, susceptibility to bean weevil damage was significantly higher for early flowering accessions and accessions with larger seeds. In this study, no yield penalty was found for higher seed protein content, but protein content was negatively correlated to starch content. Our results showed that while breeding advances in faba bean germplasm have resulted in increased yields and number of seeds per plant, they have also led to a selection pressure towards delayed onset of flowering and maturity. DArTseq genotyping identified 6,606 single nucleotide polymorphisms (SNPs) by alignment to the faba bean reference genome. These SNPs were used in a GWAS, revealing 51 novel SNP markers significantly associated with ten of the assessed traits. Three markers for days to flowering were found in predicted genes encoding proteins for which homologs in other plant species regulate flowering. Altogether, this work enriches the growing pool of phenotypic and genotypic data on faba bean as a valuable resource for developing efficient breeding strategies to expand crop cultivation.

Favorable Loci Identified for Stripe Rust Resistance in Chinese Winter Wheat Accessions via Genome-Wide Association Study

Stripe rust (or yellow rust), caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most devastating diseases of wheat worldwide. Currently, the utilization of resistant cultivars is the most viable way to reduce yield losses. In this study, a panel of 188 wheat accessions from China was evaluated for stripe rust resistance, and genome-wide association studies were performed using high-quality Diversity Arrays Technology markers. According to the phenotype and genotype data, a total of 26 significant marker-trait associations were identified, representing 18 quantitative trait loci (QTLs) on chromosomes 1B, 2A, 2B, 3A, 3B, 5A, 5B, 6B, 7B, and 7D. Of the 18 QTLs, almost all were associated with adult plant resistance (APR) except QYr.nwsuaf-6B.2, which was associated with all-stage resistance (also known as seedling resistance). Three of the 18 QTLs were mapped far from previously identified Pst resistance genes and QTLs and were considered potentially new loci. The other 15 QTLs were mapped close to known resistance genes and QTLs. Subsequent haplotype analysis for QYr.nwsuaf-2A and QYr.nwsuaf-7B.3 revealed the degrees of resistance of the panel in the APR stage. In summary, the favorable alleles identified in this study may be useful in breeding for disease resistance to stripe rust.

Genome-Wide Association Study of Preharvest Sprouting in Wheat

Genome-wide association study (GWAS) is widely used to characterize genes or quantitative trait loci (QTLs) associated with preharvest sprouting and seed dormancy. GWAS can identify both previously discovered and novel QTLs across diverse genetic panels. The high-throughput SNP arrays or next-generation sequencing technologies have facilitated the identification of numerous genetic markers, thereby significantly enhancing the resolution of GWAS. Although various methods have been developed, the fundamental principles underlying these techniques remain constant. Here, we provide a basic technological flow to perform seed dormancy assay, followed by GWAS using population structure control, and compared it with previous identified QTLs and genes.

Genetic architecture of adult-plant resistance to stripe rust in bread wheat (Triticum aestivum L.) association panel

Stripe rust, caused by Puccinia striiformis f. sp. tritici, is a severe disease in wheat worldwide, including Ethiopia, causing up to 100% wheat yield loss in the worst season. The use of resistant cultivars is considered to be the most effective and durable management technique for controlling the disease. Therefore, the present study targeted the genetic architecture of adult plant resistance to yellow rust in 178 wheat association panels. The panel was phenotyped for yellow rust adult-plant resistance at three locations. Phonological, yield, yield-related, and agro-morphological traits were recorded. The association panel was fingerprinted using the genotyping-by-sequencing (GBS) platform, and a total of 6,788 polymorphic single nucleotide polymorphisms (SNPs) were used for genome-wide association analysis to identify effective yellow rust resistance genes. The marker-trait association analysis was conducted using the Genome Association and Prediction Integrated Tool (GAPIT). The broad-sense heritability for the considered traits ranged from 74.52% to 88.64%, implying the presence of promising yellow rust resistance alleles in the association panel that could be deployed to improve wheat resistance to the disease. The overall linkage disequilibrium (LD) declined within an average physical distance of 31.44 Mbp at r2 = 0.2. Marker-trait association (MTA) analysis identified 148 loci significantly (p = 0.001) associated with yellow rust adult-plant resistance. Most of the detected resistance quantitative trait loci (QTLs) were located on the same chromosomes as previously reported QTLs for yellow rust resistance and mapped on chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 3B, 3D, 4A, 4B, 4D, 5A, 5B, 6A, 6B, 7A, and 7D. However, 12 of the discovered MTAs were not previously documented in the wheat literature, suggesting that they could represent novel loci for stripe rust resistance. Zooming into the QTL regions in IWGSC RefSeq Annotation v1 identified crucial disease resistance-associated genes that are key in plants' defense mechanisms against pathogen infections. The detected QTLs will be helpful for marker-assisted breeding of wheat to increase resistance to stripe rust. Generally, the present study identified putative QTLs for field resistance to yellow rust and some important agronomic traits. Most of the discovered QTLs have been reported previously, indicating the potential to improve wheat resistance to yellow rust by deploying the QTLs discovered by marker-assisted selection.

Using plants in forensics: State-of-the-art and prospects

The increasing use of plant evidence in forensic investigations gave rise to a powerful new discipline - Forensic Botany - that analyses micro- or macroscopic plant materials, such as the totality or fragments of an organ (i.e., leaves, stems, seeds, fruits, roots) and tissue (i.e., pollen grains, spores, fibers, cork) or its chemical composition (i. e., secondary metabolites, isotopes, DNA, starch grains). Forensic botanists frequently use microscopy, chemical analysis, and botanical expertise to identify and interpret evidence crucial to solving civil and criminal issues, collaborating in enforcing laws or regulations, and ensuring public health safeguards. The present work comprehensively examines the current state and future potential of Forensic Botany. The first section conveys the critical steps of plant evidence collection, documentation, and preservation, emphasizing the importance of these initial steps in maintaining the integrity of the items. It explores the different molecular analyses, covering the identification of plant species and varieties or cultivars, and discusses the limitations and challenges of these techniques in forensics. The subsequent section covers the diversity of Forensic Botany approaches, examining how plant evidence exposes food and pharmaceutical frauds, uncovers insufficient or erroneous labeling, traces illegal drug trafficking routes, and combats the illegal collection or trade of protected species and derivatives. National and global security issues, including the implications of biological warfare, bioterrorism, and biocrime are addressed, and a review of the contributions of plant evidence in crime scene investigations is provided, synthesizing a comprehensive overview of the diverse facets of Forensic Botany.

Mapping of a Stripe Rust Resistance Gene Yr72 in the Common Wheat Landraces AUS27506 and AUS27894 from the Watkins Collection

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is among the major threats to global wheat production. The common wheat landraces AUS27506 and AUS27894 displayed stripe rust resistance against several commercially prevailing Pst pathotypes. These genotypes were crossed with a stripe-rust-susceptible landrace AUS27229 to understand the inheritance of resistance and to determine the genomic location(s) of underlying gene(s). F3 generations of crosses AUS27506/AUS27229 and AUS27894/AUS27229 showed monogenic segregation for stripe rust resistance under greenhouse conditions. The absence of segregation for stripe rust response among the AUS27506/AUS27894-derived F3 population suggested that both genotypes carry the same gene. The stripe rust resistance gene carried by AUS27506 and AUS27894 was tentatively named YrAW4. A bulked segregant analysis placed YrAW4 in the long arm of chromosome 2B. The AUS27506/AUS27229 F3 population was enhanced to develop an F6 recombinant inbred line (RIL) population for detailed mapping of chromosome 2BL. DArT-based SSR, STS and SNP markers were employed to enrich the 2BL map. DArT-based STS markers sun481 and SNP marker IWB12294 flanked YrAW4 proximally (1.8 cM) and distally (1.2 cM), respectively. Deletion mapping placed sun481 in the deletion bin 2BL-5. All stripe rust resistance genes, previously located on chromosome 2BL, neither express an infection type like YrAW4, nor are they mapped in the deletion bin 2BL-5. Hence, YrAW4 represented a new locus and was formally named Yr72. The usefulness of the markers IWB12294 and sun481 in marker-assisted selection was demonstrated by the amplification of alleles that are different to that linked with Yr72 in 19 common wheat and two durum wheat cultivars.

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.

Transfer of the high-temperature adult-plant stripe rust resistance gene Yr62 in four Chinese wheat cultivars

Wheat stripe rust is one of the diseases that seriously affect wheat production worldwide. Breeding resistant cultivars is an effective way to control this disease. The wheat stripe rust resistance gene Yr62 has high-temperature adult-plant resistance (HTAP). In this study, PI 660,060, a single Yr62 gene line, was crossed with four Chinese wheat cultivars, LunXuan987 (LX987), Bainongaikang58 (AK58), ZhengMai9023 (ZM9023), and HanMai6172 (H6172). F1 seeds of four cross combinations were planted and self-crossed to develop the advance generations in the field. The seeds of each cross were mixed harvested and about 2400 to 3000 seeds were sown in each generation for F1 to F4 to maintain the maximum possible genotypes. Forty-five lines were selected and evaluated for resistance to stripe rust and agronomic traits, including plant height, number of grains per spike, and tiller number, in F5 and F6. Then, 33 lines with good agronomic traits and high disease resistance were developed to F9 generation. SSR markers Xgwm251 and Xgwm192 flank linked with the Yr62 were used to detect the presence of Yr62 in these 33 F9 lines. Of these, 22 lines were confirmed with the resistance gene Yr62. Finally, nine lines with good agronomic traits and disease resistance were successfully selected. The selected wheat lines in this study provide material support for the future breeding of wheat for stripe rust resistance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01393-1.

Mapping and characterization of the recessive leaf rust resistance gene Lr83 on wheat chromosome arm 1DS

KEY MESSAGE

The leaf rust resistance gene in Thatcher wheat derivative 78-1 was mapped to chromosome 1DS with SNP markers and designated as Lr83. 'Thatcher' wheat near isogenic line RL6149, a putative derivative of Triticum dicoccoides, was previously determined to carry leaf rust resistance gene Lr64 on chromosome arm 6AL and a second gene temporarily named LrX on chromosome arm 1DS. The objective of this study was to map and characterize LrX in a population of recombinant inbred lines (RILs) that segregated for a single gene. Thatcher line 78-1 with LrX was crossed with Thatcher and individual F₂ seedlings and F₆ RILs were evaluated for leaf rust response. The 208 F₂ plants segregated for a single recessive gene and 148 F₆ lines for a single gene. The RILs and parents were characterized by genotyping by sequencing (GBS). Six GBS markers and five Kompetitive Allele-Specific PCR (KASP) markers were used to map LrX on the distal region of chromosome arm 1DS. LrX was 1 centiMorgan (cM) proximal to marker K-IWB38437 and 0.4 cM distal to GBS marker 1D_9037138. Line 78-1 was crossed with Thatcher wheat lines with Lr21, Lr42, and Lr60 for allelism tests. LrX mapped 19.49 cM from Lr21 and 11.93 cM from Lr42. In the cross of line 78-1 with the Thatcher line with Lr60, one recombinant in 1,003 F₂ plants was found. LrX and Lr60 are at tightly linked loci on the distal region of chromosome arm 1DS. The gene in line 78-1 was designated as Lr83. Cytological examination of RL6149 provided no evidence of transfer of a chromosome segment of an A- or B-genome chromosome to chromosome 1D.

Comparison of linear and semi-parametric models incorporating genomic, pedigree, and associated loci information for the prediction of resistance to stripe rust in an Austrian winter wheat breeding program

We used a historical dataset on stripe rust resistance across 11 years in an Austrian winter wheat breeding program to evaluate genomic and pedigree-based linear and semi-parametric prediction methods. Stripe rust (yellow rust) is an economically important foliar disease of wheat (Triticum aestivum L.) caused by the fungus Puccinia striiformis f. sp. tritici. Resistance to stripe rust is controlled by both qualitative (R-genes) and quantitative (small- to medium-effect quantitative trait loci, QTL) mechanisms. Genomic and pedigree-based prediction methods can accelerate selection for quantitative traits such as stripe rust resistance. Here we tested linear and semi-parametric models incorporating genomic, pedigree, and QTL information for cross-validated, forward, and pairwise prediction of adult plant resistance to stripe rust across 11 years (2008-2018) in an Austrian winter wheat breeding program. Semi-parametric genomic modeling had the greatest predictive ability and genetic variance overall, but differences between models were small. Including QTL as covariates improved predictive ability in some years where highly significant QTL had been detected via genome-wide association analysis. Predictive ability was moderate within years (cross-validated) but poor in cross-year frameworks.

Salt tolerance QTLs of an endemic rice landrace, Horkuch at seedling and reproductive stages

Salinity has a significant negative impact on production of rice. To cope with the increased soil salinity due to climate change, we need to develop salt tolerant rice varieties that can maintain their high yield. Rice landraces indigenous to coastal Bangladesh can be a great resource to study the genetic basis of salt adaptation. In this study, we implemented a QTL analysis framework with a reciprocal mapping population developed from a salt tolerant landrace Horkuch and a high yielding rice variety IR29. Our aim was to detect genetic loci that contributes to the salt adaptive responses of the two different developmental stages of rice which are very sensitive to salinity stress. We identified 14 QTLs for 9 traits and found that most are unique to specific developmental stages. In addition, we detected a significant effect of the cytoplasmic genome on the QTL model for some traits such as leaf total potassium and filled grain weight. This underscores the importance of considering cytoplasm-nuclear interaction for breeding programs. Finally, we identified QTLs co-localization for multiple traits that highlights the possible constraint of multiple QTL selection for breeding programs due to different contributions of a donor allele for different traits.

Genetic incorporation of genes for the optimal plant architecture in common wheat

Wheat grain yield is affected by plant height, which is the total length of spike, the uppermost internode, and other elongated internodes. In this study, a population of recombinant inbred lines generated from a cross between two advanced winter wheat breeding lines were phenotyped over four locations/years and genotyped by using markers of genotyping-by-sequencing (GBS) and Diversity Array Technology (DArT) for mapping of genes for three traits, including spike length, the uppermost internode length, and plant height. Five genomic regions or quantitative trait loci (QTLs) were associated with candidate genes for these traits. A major QTL was associated with Q5A, and two novel haplotypes of Q5A were identified, one for a single nucleotide polymorphism (SNP) at position -2,149 in promoter region and the other for copy number variation. Compared with one copy Q5A on chromosome 5A in Chinese Spring, the novel haplotype of Q5A with two copies Q5A was found to generate spikes that are extremely compacted. A major QTL was associated with allelic variation in the recessive vrn-A1 alleles involving in protein sequences, and this QTL was associated with increased uppermost internode length but not with plant height. A major QTL for plant height was associated with Rht-B1b on chromosome 4B, but its effects could be compromised by two new minor QTLs on chromosome 7. Collectively, the favorable alleles from the four loci can be used to establish the optimal plant height in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01336-2.

Multi-locus genome-wide association studies reveal novel alleles for flowering time under vernalisation and extended photoperiod in a barley MAGIC population

Key genes controlling flowering and interactions of different photoperiod alleles with various environments were identified in a barley MAGIC population. A new candidate gene for vernalisation requirements was also detected. Optimal flowering time has a major impact on grain yield in crop species, including the globally important temperate cereal crop barley (Hordeum vulgare L.). Understanding the genetics of flowering is a key avenue to enhancing yield potential. Although bi-parental populations were used intensively to map genes controlling flowering, their lack of genetic diversity requires additional work to obtain desired gene combinations in the selected lines, especially when the two parental cultivars did not carry the genes. Multi-parent mapping populations, which use a combination of four or eight parental cultivars, have higher genetic and phenotypic diversity and can provide novel genetic combinations that cannot be achieved using bi-parental populations. This study uses a Multi-parent advanced generation intercross (MAGIC) population from four commercial barley cultivars to identify genes controlling flowering time in different environmental conditions. Genome-wide association studies (GWAS) were performed using 5,112 high-quality markers from Diversity Arrays Technology sequencing (DArT-seq), and Kompetitive allele-specific polymerase chain reaction (KASP) genetic markers were developed. Phenotypic data were collected from fifteen different field trials for three consecutive years. Planting was conducted at various sowing times, and plants were grown with/without additional vernalisation and extended photoperiod treatments. This study detected fourteen stable regions associated with flowering time across multiple environments. GWAS combined with pangenome data highlighted the role of CEN gene in flowering and enabled the prediction of different CEN alleles from parental lines. As the founder lines of the multi-parental population are elite germplasm, the favourable alleles identified in this study are directly relevant to breeding, increasing the efficiency of subsequent breeding strategies and offering better grain yield and adaptation to growing conditions.

SilicoDArT and SNP markers for genetic diversity and population structure analysis of Trema orientalis; a fodder species

Establishing the genetic diversity and population structure of a species can guide the selection of appropriate conservation and sustainable utilization strategies. Next-generation sequencing (NGS) approaches are increasingly being used to generate multi-locus data for genetic structure determination. This study presents the genetic structure of a fodder species -Trema orientalis based on two genome-wide high-throughput diversity array technology (DArT) markers; silicoDArT and single nucleotide polymorphisms (SNPs). Genotyping of 119 individuals generated 40,650 silicoDArT and 4767 SNP markers. Both marker types had a high average scoring reproducibility (>99%). Genetic relationships explored by principal coordinates analysis (PCoA) showed that the first principal coordinate axis explained most of the variation in both the SilicoDArT (34.2%) and SNP (89.6%) marker data. The average polymorphic information content did not highly differ between silicoDArT (0.22) and SNPs (0.17) suggesting minimal differences in informativeness in the two groups of markers. The, mean observed (H_o) and expected (H_e) heterozygosity were low and differed between the silicoDArT and SNPs respectively, estimated at H_o = 0.08 and H_e = 0.05 for silicoDArT and H_o = 0.23 and H_e = 0.19 for SNPs. The population of T. orientalis was moderately differentiated (F_ST  =  0.20–0.53) and formed 2 distinct clusters based on maximum likelihood and principal coordinates analysis. Analysis of molecular variance revealed that clusters contributed more to the variation (46.3–60.8%) than individuals (32.9–31.2%). Overall, the results suggest a high relatedness of the individuals sampled and a threatened genetic potential of T. orientalis in the wild. Therefore, genetic management activities such as ex-situ germplasm management are required for the sustainability of the species. Ex-situ conservation efforts should involve core collection of individuals from different populations to capture efficient diversity. This study demonstrates the importance of silicoDArT and SNP makers in population structure and genetic diversity analysis of Trema orientalis, useful for future genome wide studies in the species.

Transpirational Leaf Cooling Effect Did Not Contribute Equally to Biomass Retention in Wheat Genotypes under High Temperature

High temperature and water deficit are the most critical yield-limiting environmental factors for wheat in rainfed environments. It is important to understand the heat avoidance mechanisms and their associations with leaf morpho-physiological traits that allow crops to stay cool and retain high biomass under warm and dry conditions. We examined 20 morpho-physiologically diverse wheat genotypes under ambient and elevated temperatures (T_air) to investigate whether increased water use leads to high biomass retention due to increased leaf cooling. An experiment was conducted under well-watered conditions in two partially controlled glasshouses. We measured plant transpiration (T_r), leaf temperature (T_leaf), vapor pressure deficit (VPD), and associated leaf morpho-physiological characteristics. High water use and leaf cooling increased biomass retention under high temperatures, but increased use did not always increase biomass retention. Some genotypes maintained biomass, irrespective of water use, possibly through mechanisms other than leaf cooling, indicating their adaptation under water shortage. Genotypic differences in leaf cooling capacity did not always correlate with T_r (VPD) response. In summary, the contribution of high water use or the leaf cooling effect on biomass retention under high temperature is genotype-dependent and possibly due to variations in leaf morpho-physiological traits. These findings are useful for breeding programs to develop climate resilient wheat cultivars.

Genetic and genomic diversity in the sorghum gene bank collection of Uganda

BACKGROUND

The Plant Genetic Resources Centre at the Uganda National Gene Bank houses has over 3000 genetically diverse landraces and wild relatives of Sorghum bicolor accessions. This genetic diversity resource is untapped, under-utilized, and has not been systematically incorporated into sorghum breeding programs. In this study, we characterized the germplasm collection using whole-genome SNP markers (DArTseq). Discriminant analysis of principal components (DAPC) was implemented to study the racial ancestry of the accessions in comparison to a global sorghum diversity set and characterize the sub-groups present in the Ugandan (UG) germplasm.

RESULTS

Population structure and phylogenetic analysis revealed the presence of five subgroups among the Ugandan accessions. The samples from the highlands of the southwestern region were genetically distinct as compared to the rest of the population. This subset was predominated by the caudatum race and unique in comparison to the other sub-populations. In this study, we detected QTL for juvenile cold tolerance by genome-wide association studies (GWAS) resulting in the identification of 4 markers associated (-log10p > 3) to survival under cold stress under both field and climate chamber conditions, located on 3 chromosomes (02, 06, 09). To our best knowledge, the QTL on Sb09 with the strongest association was discovered for the first time.

CONCLUSION

This study demonstrates how genebank genomics can potentially facilitate effective and efficient usage of valuable, untapped germplasm collections for agronomic trait evaluation and subsequent allele mining. In face of adverse climate change, identification of genomic regions potentially involved in the adaptation of Ugandan sorghum accessions to cooler climatic conditions would be of interest for the expansion of sorghum production into temperate latitudes.

Exploration of the Genetic Diversity of Solina Wheat and Its Implication for Grain Quality

Different Solina wheat accessions (n = 24) collected in the Abruzzo region (Italy) were studied using 45,000 SNP markers generated from the DarTseq platform. The structure of genetic data was analyzed by Principal Component Analysis and Hierarchical Cluster analysis that revealed the existence of two main clusters (Clu1 and Clu2) characterized by samples with different geographical origin. The Solina genetic dataset was further merged and analyzed with a public genetic one provided by CIMMYT containing 25,963 genotypes from all over the world. The Solina accessions occupied a vast space, thus confirming a high heterogeneity of this landrace that, nevertheless, is considerably unique and placed quite far from other clusters. Clu1 and Clu2 divergence were clearly visible. Solina clusters were genetically closer to landraces from Turkey and the central fertile crescent than to the Italian genotypes present in the dataset. Selected commercial quality traits of accessions of the two Solina clusters were analyzed (yield, thousand kernel weight, test weight, and protein content), and significant differences were found between clusters. The results of this investigation did not highlight any relationships of Solina with Italian genotypes, and confirmed its wide genetic diversity by permitting to identify two genetic groups with distinct origin and quality traits.

Elucidating the Drought Responsiveness in Wheat Genotypes

Drought is one of the major abiotic constraints on wheat yields and also for sustainability of production levels around the world. In the near future, the occurrence likelihood of droughts is predicted to become more common, due to changing climatic conditions, thereby posing a serious threat to the food security system. Heterogeneity, in its time of occurrence and severity levels, is likely to further augment the complexity of drought conditions. Although wheat crop growth has progressively risen to good levels, as evident by notable increases in both area and production, the expected wheat demand for the ever-growing population is quite high. Besides crop yield volatility in the era of climate change and dwindling resources, “trait-based” breeding programs are required, so as to develop high yielding, climate resilient and stable genotypes, at a faster pace. For this to happen, a broad genetic base and wider adaptability to suit varied agro-ecologies would provide enough scope for their quicker spread. The current review places emphasis on making distinct categories of the wheat cultivars/advanced breeding lines, as tolerant, moderately tolerant or susceptible to drought stresses, duly supported by an extensive up-to-date literature base and will be useful for wheat researchers, in order to choose the best potential donors as parents, coupled with the associated traits for the development of drought-tolerant wheat varieties, and also to facilitate molecular studies.

Unraveling the Genetic Basis of Key Agronomic Traits of Wrinkled Vining Pea (Pisum sativum L.) for Sustainable Production

Estimating the allelic variation and exploring the genetic basis of quantitatively inherited complex traits are the two foremost breeding scenarios for sustainable crop production. The current study utilized 188 wrinkled vining pea genotypes comprising historical varieties and breeding lines to evaluate the existing genetic diversity and to detect molecular markers associated with traits relevant to vining pea production, such as wrinkled vining pea yield (YTM100), plant height (PH), earliness (ERL), adult plant resistance to downy mildew (DM), pod length (PDL), numbers of pods per plant (PDP), number of peas per pod (PPD), and percent of small wrinkled vining peas (PSP). Marker-trait associations (MTAs) were conducted using 6902 quality single nucleotide polymorphism (SNP) markers generated from the diversity arrays technology sequencing (DArTseq) and Genotyping-by-sequencing (GBS) sequencing methods. The best linear unbiased prediction (BLUP) values were estimated from the two-decades-long (1999-2020) unbalanced phenotypic data sets recorded from two private breeding programs, the Findus and the Birds eye, now owned by Nomad Foods. Analysis of variance revealed a highly significant variation between genotypes and genotype-by-environment interactions for the ten traits. The genetic diversity and population structure analyses estimated an intermediate level of genetic variation with two optimal sub-groups within the current panel. A total of 48 significant (P < 0.0001) MTAs were identified for eight different traits, including five for wrinkled vining pea yield on chr2LG1, chr4LG4, chr7LG7, and scaffolds (two), and six for adult plant resistance to downy mildew on chr1LG6, chr3LG5 (two), chr6LG2, and chr7LG7 (two). We reported several novel MTAs for different crucial traits with agronomic importance in wrinkled vining pea production for the first time, and these candidate markers could be easily validated and integrated into the active breeding programs for marker-assisted selection.

Genetic Diversity and Population Structure of Myanmar Rice (Oryza sativa L.) Varieties Using DArTseq-Based SNP and SilicoDArT Markers

Myanmar is well known as a primary center of plant genetic resources for rice. A considerable number of genetic diversity studies have been conducted in Myanmar using various DNA markers. However, this is the first report using DArTseq technology for exploring the genetic diversity of Myanmar rice. In our study, two ultra-high-throughput diversity array technology markers were employed to investigate the genetic diversity and population structure of local rice varieties in the Ayeyarwady delta, the major region of rice cultivation. The study was performed using 117 rice genotypes with 7643 SNP and 4064 silicoDArT markers derived from the DArT platform. Genetic variance among the genotypes ranged from 0 to 0.753 in SNPs, and from 0.001 to 0.954 in silicoDArT. Two distinct population groups were identified from SNP data analysis. Cluster analysis with both markers clearly separated traditional Pawsan varieties and modern high-yielding varieties. A significant divergence was found between populations according to the Fst values (0.737) obtained from the analysis of molecular variance, which revealed 74% genetic variation at the population level. These findings support rice researchers in identifying useful DNA polymorphisms in genes and pinpointing specific genes conferring desirable phenotypic traits for further genome-wide association studies and parental selection for recombination breeding to enhance rice varietal development and release.

Genome-Wide Association Study Reveals Novel Genetic Loci for Quantitative Resistance to Septoria Tritici Blotch in Wheat (Triticum aestivum L.)

Septoria tritici blotch, caused by the fungus Zymoseptoria titici, poses serious and persistent challenges to wheat cultivation in Ethiopia and worldwide. Deploying resistant cultivars is a major component of controlling septoria tritici blotch (STB). Thus, the objective of this study was to elucidate the genomic architecture of STB resistance in an association panel of 178 bread wheat genotypes. The association panel was phenotyped for STB resistance, phenology, yield, and yield-related traits in three locations for 2 years. The panel was also genotyped for single nucleotide polymorphism (SNP) markers using the genotyping-by-sequencing (GBS) method, and a total of 7,776 polymorphic SNPs were used in the subsequent analyses. Marker-trait associations were also computed using a genome association and prediction integrated tool (GAPIT). The study then found that the broad-sense heritability for STB resistance ranged from 0.58 to 0.97 and 0.72 to 0.81 at the individual and across-environment levels, respectively, indicating the presence of STB resistance alleles in the association panel. Population structure and principal component analyses detected two sub-groups with greater degrees of admixture. A linkage disequilibrium (LD) analysis in 338,125 marker pairs also detected the existence of significant (p ≤ 0.01) linkage in 27.6% of the marker pairs. Specifically, in all chromosomes, the LD between SNPs declined within 2.26-105.62 Mbp, with an overall mean of 31.44 Mbp. Furthermore, the association analysis identified 53 loci that were significantly (false discovery rate, FDR, <0.05) associated with STB resistance, further pointing to 33 putative quantitative trait loci (QTLs). Most of these shared similar chromosomes with already published Septoria resistance genes, which were distributed across chromosomes 1B, 1D, 2A, 2B, 2D, 3A,3 B, 3D, 4A, 5A, 5B, 6A, 7A, 7B, and 7D. However, five of the putative QTLs identified on chromosomes 1A, 5D, and 6B appeared to be novel. Dissecting the detected loci on IWGSC RefSeq Annotation v2.1 revealed the existence of disease resistance-associated genes in the identified QTL regions that are involved in plant defense responses. These putative QTLs explained 2.7-13.2% of the total phenotypic variation. Seven of the QTLs (R ² = 2.7-10.8%) for STB resistance also co-localized with marker-trait associations (MTAs) for agronomic traits. Overall, this analysis reported on putative QTLs for adult plant resistance to STB and some important agronomic traits. The reported and novel QTLs have been identified previously, indicating the potential to improve STB resistance by pyramiding QTLs by marker-assisted selection.

Genomic signatures of selection for resistance to stripe rust in Austrian winter wheat

We combined quantitative and population genetic methods to identify loci under selection for adult plant resistance to stripe rust in an Austrian winter wheat breeding population from 2008 to 2018. Resistance to stripe rust, a foliar disease caused by the fungus P. striiformis f. sp. tritici, in wheat (Triticum aestivum L.) is both qualitatively and quantitatively controlled. Resistance genes confer complete, race-specific resistance but are easily overcome by evolving pathogen populations, while quantitative resistance is controlled by many small- to medium-effect loci that provide incomplete yet more durable protection. Data on resistance loci can be applied in marker-assisted selection and genomic prediction frameworks. We employed genome-wide association to detect loci associated with stripe rust and selection testing to identify regions of the genome that underwent selection for stripe rust resistance in an Austrian winter wheat breeding program from 2008 to 2018. Genome-wide association mapping identified 150 resistance loci, 62 of which showed significant evidence of selection over time. The breeding population also demonstrated selection for resistance at the genome-wide level.

High-density DArT-based SilicoDArT and SNP markers for genetic diversity and population structure studies in cassava (Manihot esculenta Crantz)

Cassava (Manihot esculenta Crantz) is an important industrial and staple crop due to its high starch content, low input requirement, and resilience which makes it an ideal crop for sustainable agricultural systems and marginal lands in the tropics. However, the lack of genomic information on local genetic resources has impeded efficient conservation and improvement of the crop and the exploration of its full agronomic and breeding potential. This work was carried out to obtain information on population structure and extent of genetic variability among some local landraces conserved at the Plant Genetic Resources Research Institute, Ghana and exotic cassava accessions with Diversity Array Technology based SilicoDArT and SNP markers to infer how the relatedness in the genetic materials can be used to enhance germplasm curation and future breeding efforts. A total of 10521 SilicoDArT and 10808 SNP markers were used with varying polymorphic information content (PIC) values. The average PIC was 0.36 and 0.28 for the SilicoDArT and SNPs respectively. Population structure and average linkage hierarchical clustering based on SNPs revealed two distinct subpopulations and a large number of admixtures. Both DArT platforms identified 22 landraces as potential duplicates based on Gower's genetic dissimilarity. The expected heterozygosity which defines the genetic variation within each subpopulation was 0.008 for subpop1 which were mainly landraces and 0.391 for subpop2 indicating the homogeneous and admixture nature of the two subpopulations. Further analysis upon removal of the duplicates increased the expected heterozygosity of subpop1 from 0.008 to 0.357. A mantel test indicated strong interdependence (r = 0.970; P < 0.001) between SilicoDArT and DArTSeq SNP genotypic data suggesting both marker platforms as a robust system for genomic studies in cassava. These findings provide important information for efficient ex-situ conservation of cassava, future heterosis breeding, and marker-assisted selection (MAS) to enhance cassava improvement.

Using chlorate as an analogue to nitrate to identify candidate genes for nitrogen use efficiency in barley

Nitrogen (N) is one of the most important macronutrients for crop growth and development. Large amounts of N fertilizers are applied exogenously to improve grain yield and quality, which has led to environmental pollution and high cost of production. Therefore, improvement of N use efficiency (NUE) is a very important aspect for sustainable agriculture. Here, a pilot experiment was firstly conducted with a set of barley genotypes with confirmed NUE to validate the fast NUE screening, using chlorate as an analogue to nitrate. High NUE genotypes were susceptible to chlorate-induced toxicity whereas the low NUE genotypes were tolerant. A total of 180 barley RILs derived from four parents (Compass, GrangeR, Lockyer and La Trobe) were further screened for NUE. Leaf chlorosis induced by chlorate toxicity was the key parameter observed which was later related to low-N tolerance of the RILs. There was a distinct variation in chlorate susceptibility of the RILs with leaf chlorosis in the oldest leaf ranging from 10 to 80%. A genome-wide association study (GWAS) identified 9 significant marker-trait associations (MTAs) conferring high chlorate sensitivity on chromosomes 2H (2), 3H (1), 4H (4), 5H (1) and Un (1). Genes flanking with these markers were retrieved as potential targets for genetic improvement of NUE. Genes encoding Ferredoxin 3, leucine-rich receptor-like protein kinase family protein and receptor kinase are responsive to N stress. MTA4H5468 which exhibits concordance with high NUE phenotype can further be explored under different genetic backgrounds and successfully applied in marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01239-8.

QTL Mapping and Phenotypic Variation for Seedling Vigour Traits in Barley (Hordeum vulgare L.)

Seed vigour is considered a critical stage for barley production, and cultivars with early seedling vigour (ESV) facilitate rapid canopy formation. In this study, QTLs for 12 ESV-related traits were mapped using 185 RILs derived from a Xena x H94061120 evaluated across six independent environments. DArT markers were used to develop a genetic map (1075.1 cM; centimorgans) with an average adjacent-marker distance of 3.28 cM. In total, 46 significant QTLs for ESV-related traits were detected. Fourteen QTLs for biomass yield were found on all chromosomes, two of them co-localized with QTLs on 1H for grain yield. The related traits: length of the first and second leaves and dry weight of the second leaf, biomass yield and grain yield, had high heritability (>30%). Meanwhile, a significant correlation was observed between grain yield and biomass yield, which provided a clear image of these traits in the selection process. Our results demonstrate that a pleiotropic QTL related to the specific leaf area of the second leaf, biomass yield, and grain yield was linked to the DArT markers bPb-9280 and bPb-9108 on 1H, which could be used to significantly improve seed vigour by marker-assisted selection and facilitate future map-based cloning efforts.

DArTseq-Based High-Throughput SilicoDArT and SNP Markers Applied for Association Mapping of Genes Related to Maize Morphology

Today, agricultural productivity is essential to meet the needs of a growing population, and is also a key tool in coping with climate change. Innovative plant breeding technologies such as molecular markers, phenotyping, genotyping, the CRISPR/Cas method and next-generation sequencing can help agriculture meet the challenges of the 21st century more effectively. Therefore, the aim of the research was to identify single-nucleotide polymorphisms (SNPs) and SilicoDArT markers related to select morphological features determining the yield in maize. The plant material consisted of ninety-four inbred lines of maize of various origins. These lines were phenotyped under field conditions. A total of 14 morphological features was analyzed. The DArTseq method was chosen for genotyping because this technique reduces the complexity of the genome by restriction enzyme digestion. Subsequently, short fragment sequencing was used. The choice of a combination of restrictases allowed the isolation of highly informative low copy fragments of the genome. Thanks to this method, 90% of the obtained DArTseq markers are complementary to the unique sequences of the genome. All the observed features were normally distributed. Analysis of variance indicated that the main effect of lines was statistically significant (p < 0.001) for all 14 traits of study. Thanks to the DArTseq analysis with the use of next-generation sequencing (NGS) in the studied plant material, it was possible to identify 49,911 polymorphisms, of which 33,452 are SilicoDArT markers and the remaining 16,459 are SNP markers. Among those mentioned, two markers associated with four analyzed traits deserved special attention: SNP (4578734) and SilicoDArT (4778900). SNP marker 4578734 was associated with the following features: anthocyanin coloration of cob glumes, number of days from sowing to anthesis, number of days from sowing to silk emergence and anthocyanin coloration of internodes. SilicoDArT marker 4778900 was associated with the following features: number of days from sowing to anthesis, number of days from sowing to silk emergence, tassel: angle between the axis and lateral branches and plant height. Sequences with a length of 71 bp were used for physical mapping. The BLAST and EnsemblPlants databases were searched against the maize genome to identify the positions of both markers. Marker 4578734 was localized on chromosome 7, the closest gene was Zm00001d022467, approximately 55 Kb apart, encoding anthocyanidin 3-O-glucosyltransferase. Marker 4778900 was located on chromosome 7, at a distance of 45 Kb from the gene Zm00001d045261 encoding starch synthase I. The latter observation indicated that these flanking SilicoDArT and SNP markers were not in a state of linkage disequilibrium.

A global barley panel revealing genomic signatures of breeding in modern Australian cultivars

The future of plant cultivar improvement lies in the evaluation of genetic resources from currently available germplasm. Today's gene pool of crop genetic diversity has been shaped during domestication and more recently by breeding. Recent efforts in plant breeding have been aimed at developing new and improved varieties from poorly adapted crops to suit local environments. However, the impact of these breeding efforts is poorly understood. Here, we assess the contributions of both historical and recent breeding efforts to local adaptation and crop improvement in a global barley panel by analysing the distribution of genetic variants with respect to geographic region or historical breeding category. By tracing the impact that breeding had on the genetic diversity of Hordeum vulgare (barley) released in Australia, where the history of barley production is relatively young, we identify 69 candidate regions within 922 genes that were under selection pressure. We also show that modern Australian barley varieties exhibit 12% higher genetic diversity than historical cultivars. Finally, field-trialling and phenotyping for agriculturally relevant traits across a diverse range of Australian environments suggests that genomic regions under strong breeding selection and their candidate genes are closely associated with key agronomic traits. In conclusion, our combined data set and germplasm collection provide a rich source of genetic diversity that can be applied to understanding and improving environmental adaptation and enhanced yields.

Detection of grain yield QTLs in the durum population Lahn/Cham1 tested in contrasting environments

Durum wheat (Triticum turgidum L. var durum) is tetraploid wheat (AABB); it is the main source of semolina and other pasta products. Grain yield in wheat is quantitatively inherited and influenced by the environment. The genetic map construction constitutes the essential step in identifying quantitative trait loci (QTLs) linked to complex traits, such as grain yield. The study aimed to construct a genetic linkage map of two parents that are widely grown durum cultivars (Lahn and Cham1) in the Mediterranean basin, which is characterized by varying climate changes. The genetic linkage map of Lahn/Cham1 population consisted of 112 recombinant inbred lines (RILs) and was used to determine QTLs linked to the grain yield in 11 contrasting environments (favorable, cold, dry, and hot). Simple sequence repeat (SSR) molecular markers were used to construct an anchor map, which was later enriched with single nucleotide polymorphisms (SNPs). The map was constructed with 247 SSRs and enriched with 1425 SNPs. The map covered 6122.22 cM. One hundred and twenty-six QTLs were detected on different chromosomes. Chromosomes 2A and 4B harbored the most significant grain yield QTLs. Furthermore, by comparison with several wheat mapping populations, all the A and B chromosomes of Lahn/Cham1 QTLs contributed to grain yield. The results showed that the detected QTLs can be used as a potential candidate for marker-assisted selection in durum breeding programs.

Adult Plant Leaf Rust Resistance in AC Taber Wheat Maps to Chromosomes 2BS and 3BS

AC Taber is a hard red spring wheat cultivar that has had long-lasting resistance to the leaf rust fungus Puccinia triticina. The objective of this study was to determine the chromosome location of the leaf rust resistance genes in AC Taber. The leaf rust-susceptible cultivar Thatcher was crossed with AC Taber to develop an F₆ recombinant inbred line (RIL) population. The RILs and parents were evaluated for segregation of leaf rust resistance in five field plot tests and in two seedling tests to race BBBDB of P. triticina. A genetic map of the RIL population was developed using 90,000 single nucleotide polymorphism markers with the Illumina Infinium iSelect 90K wheat bead array. Quantitative trait loci (QTLs) with significant effects for lower leaf rust severity in the field plot tests were found on chromosomes 2BS and 3BS. The same QTLs also had significant effects for lower infection type in seedlings to leaf rust race BBBDB. The gene on 2BS was the adult plant resistance gene Lr13, and the gene on 3BS mapped to the same region as the adult plant resistance gene Lr74 and other QTLs for leaf rust resistance. Kompetitive allele-specific PCR assay markers linked to the 2BS and 3BS regions were developed and should be useful for marker-based selection of these genes.

Genome-wide association analysis of Mexican bread wheat landraces for resistance to yellow and stem rust

Deploying under-utilized landraces in wheat breeding has been advocated to accelerate genetic gains in current era of genomics assisted breeding. Mexican bread wheat landraces (Creole wheats) represent an important resource for the discovery of novel alleles including disease resistance. A core set of 1,098 Mexican landraces was subjected to multi-location testing for rust diseases in India, Mexico and Kenya. The landrace core set showed a continuous variation for yellow (YR) and stem rust (SR) disease severity. Principal component analysis differentiated Mexican landraces into three groups based on their respective collection sites. Linkage disequilibrium (LD) decay varied from 10 to 32 Mb across chromosomes with an averge of 23Mb across whole genome. Genome-wide association analysis revealed marker-trait associations for YR resistance in India and Mexico as well as for SR resistance in Kenya. In addition, significant additive-additive interaction effects were observed for both YR and SR resistance including genomic regions on chromosomes 1BL and 3BS, which co-locate with pleiotropic genes Yr29/Lr46/Sr58/Pm39/Ltn2 and Sr2/Yr30/Lr27, respectively. Study reports novel genomic associations for YR (chromosomes 1AL, 2BS, and 3BL) and SR (chromosomes 2AL, 4DS, and 5DS). The novel findings in Creole wheat landraces can be efficiently utilized for the wheat genetic improvement.

The barley pan-genome reveals the hidden legacy of mutation breeding

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

Morphological and genome-wide evidence for natural hybridisation within the genus Stipa (Poaceae)

Hybridisation in the wild between closely related species is a common mechanism of speciation in the plant kingdom and, in particular, in the grass family. Here we explore the potential for natural hybridisation in Stipa (one of the largest genera in Poaceae) between genetically distant species at their distribution edges in Mountains of Central Asia using integrative taxonomy. Our research highlights the applicability of classical morphological and genome reduction approaches in studies on wild plant species. The obtained results revealed a new nothospecies, Stipa × lazkovii, which exhibits intermediate characters to S. krylovii and S. bungeana. A high-density DArTseq assay disclosed that S. × lazkovii is an F1 hybrid, and established that the plastid and mitochondrial DNA was inherited from S. bungeana. In addition, molecular markers detected a hybridisation event between morphologically and genetically distant species S. bungeana and probably S. glareosa. Moreover, our findings demonstrated an uncertainty on the taxonomic status of S. bungeana that currently belongs to the section Leiostipa, but it is genetically closer to S. breviflora from the section Barbatae. Finally, we noticed a discrepancy between the current molecular data with the previous findings on S. capillata and S. sareptana.

Genome-wide association studies of Ca and Mn in the seeds of the common bean (Phaseolus vulgaris L.)

SNP markers linked to genes controlling Ca and Mn uptake were identified in the common bean seeds using DArT-based association mapping (AM). The Ca concentration in the seeds varied between 475 and 3,100 mg kg^-1 with an average of 1,280.9 mg kg^-1 and the Mn concentration ranged from 4.87 to 27.54 mg kg^-1 with a mean of 11.76 mg kg^-1. A total of 19,204 SNP markers were distributed across 11 chromosomes that correspond to the haploid genome number of the common bean. The highest value of ΔK was determined as K = 2, and 173 common bean genotypes were split into two main subclusters as POP1 (Mesoamerican) and POP2 (Andean). The results of the UPGMA dendrogram and PCA confirmed those of STRUCTURE analysis. MLM based on the Q + K model identified a large number of markers-trait associations. Of the 19,204 SNPs, five (on Pv2, 3, 8, 10 and 11) and four (on Pv2, 3, 8 and 11) SNPs were detected to be significantly related to the Ca content of the beans grown in Bornova and Menemen, respectively in 2015. In 2016, six SNPs (on Pv1-4, 8 and 10) were identified to be significantly associated with the Ca content of the seeds obtained from Bornova and six SNPs (on Pv1-4, 8 and 10) from Menemen. Eight (on Pv3, 5 and 11) and four (on Pv2, 5 and 11) SNPs had a significant association with Mn content in Bornova in 2015 and 2016, respectively. In Menemen, eight (on Pv3, 5, 8 and 11) and 11 (on Pv1, 2, 5, 10 and 11) SNPs had a significant correlation with Mn content in 2015 and 2016, respectively.

The development of genome-wide single nucleotide polymorphisms in blue wildebeest using the DArTseq platform

Blue wildebeest (Connochaetes taurinus taurinus) are economically important antelope that are widely utilised in the South African wildlife industry. However, very few genomic resources are available for blue wildebeest that can assist in breeding management and facilitate research. This study aimed to develop a set of genome-wide single nucleotide polymorphism (SNP) markers for blue wildebeest. The DArTseq genotyping platform, commonly used in polyploid plant species, was selected for SNP discovery. A limited number of published articles have described the use of the DArTseq platform in animals and, therefore, this study also provided a unique opportunity to assess the performance of the DArTseq platform in an animal species. A total of 20,563 SNPs, each located within a 69 bp sequence, were generated. The developed SNP markers had a high average scoring reproducibility (>99%) and a low percentage missing data (~9.21%) compared to other reduced representation sequencing approaches that have been used in animal studies. Furthermore, the number of candidate SNPs per nucleotide position decreased towards the 3' end of sequence reads, and the ratio of transitions (Ts) to transversions (Tv) remained similar for each read position. These observations indicate that there was no read position bias, such as the identification of false SNPs due to low sequencing quality, towards the tail-end of sequencing reads. The DArTseq platform was also successful in identifying a large number of informative SNPs with desirable polymorphism parameters such as a high minor allele frequency (MAF). The Bos taurus genome was used for the in silico mapping of the marker sequences and a total of 6020 (29.28%) sequences were successfully mapped against the bovine genome. The marker sequences mapped to all of the bovine chromosomes establishing the genome-wide distribution of the SNPs. Moreover, the high observed Ts:Tv ratio (2.84:1) indicate that the DArTseq platform targeted gene-rich regions of the blue wildebeest genome. Finally, functional annotation of the marker sequences revealed a wide range of different putative functions indicating that these SNP markers can be useful in functional gene studies. The DArTseq platform, therefore, represents a high-throughput, robust and cost-effective genotyping platform, which may find adoption in several other African antelope and animal species.

Genetic and transcriptional dissection of resistance to Claviceps purpurea in the durum wheat cultivar Greenshank

Four QTL for ergot resistance (causal pathogen Claviceps purpurea) have been identified in the durum wheat cultivar Greenshank. Claviceps purpurea is a pathogen of grasses that infects flowers, replacing the seed with an ergot sclerotium. Ergot presents a significant problem to rye, barley and wheat, in particular hybrid seed production systems. In addition, there is evidence that the highly toxic alkaloids that accumulate within sclerotia can cross-contaminate otherwise healthy grain. Host resistance to C. purpurea is rare, few resistance loci having been identified. In this study, four ergot resistance loci are located on chromosomes 1B, 2A, 5A and 5B in the durum wheat cv. Greenshank. Ergot resistance was assessed through analysis of phenotypes associated with C. purpurea infection, namely the number of inoculated flowers that produced sclerotia, or resulted in ovary death but no sclerotia, the levels of honeydew produced, total sclerotia weight and average sclerotia weight and size per spike. Ergot testing was undertaken in Canada and the UK. A major effect QTL, QCp.aafc.DH-2A, was detected in both the Canadian and UK experiments and had a significant effect on honeydew production levels. QCp.aafc.DH-5B had the biggest influence on total sclerotia weight per spike. QCp.aafc.DH-1B was only detected in the Canadian experiments and QCp.aafc.DH-5A in the UK experiment. An RNASeq analysis, undertaken to identify wheat differentially expressed genes associated with different combinations of the four ergot resistance QTL, revealed a disproportionate number of DEGs locating to the QCp.aafc.DH-1B, QCp.aafc.DH-2A and QCp.aafc.DH-5B QTL intervals.

In-depth genome diversity, population structure and linkage disequilibrium analysis of worldwide diverse safflower (Carthamus tinctorius L.) accessions using NGS data generated by DArTseq technology

Safflower (Carthamus tinctorius L.) is one of the most important oilseed crops for its seed oil rich in unsaturated fatty acids. Precise utilization of diverse genetic resources is fundamental in breeding programs to improve high yield genotypes with desirable traits. In this study, for the first time we report successful application of DArTseq technology; an efficient genotyping-by-sequencing (NGS); to analysis genetic diversity and population structure of 89 safflower accessions from worldwide origins. Totally, 19,639 DArTseq markers (10,130 SilicoDArTs and 9509 SNPs) generated through DArTseq genotyping. After filtering the data, 3431 polymorphic DArTseq markers (1136 SilicoDArTs and 2295 SNPs) used for genetic diversity, population structure and linkage disequilibrium analysis in safflower genotypes. All the SilicoDArT and SNP markers showed high reproducibility and call rate. Polymorphism information content (PIC) values ranged from 0.1 to 0.5, while ≥ 0.50% of SilicoDArTs and ≥ 0.64% SNPs showed PIC values more than median. Genotypes grouping using DArTseq markers resulted in three distinct clusters. Results showed weak correlation between safflower diversity pattern and origins. Analysis of molecular variance revealed that the majority of genetic variation was attributed to the differences among varieties within cluster populations and there was no significant molecular variance between origins. However, safflower of accessions belonged to Iran, Turkey, Pakistan and India indeed appear to be genetically similar and grouped close in referred cluster, while the accessions from Near East (Afghanistan, China) being distinct. Our results were in agreement with hypothesis that safflower domesticated in somewhere west of Fertile Crescent and then expanded through Africa and Europe. Present study using a panel of globally diverse safflower accessions and large number of DArTseq markers set the stage for future analysis of safflower domestication using large germplasm from proposed domestication centers. Also, studied germplasm in this study can be used as a valuable source for future genomic studies in safflower for mapping desirable traits through genome-wide association mapping studies.

Whole-genome diversity, population structure and linkage disequilibrium analysis of globally diverse wheat genotypes using genotyping-by-sequencing DArTseq platform

In this study, 129 wheat genotypes from globally diverse origins were genotyped using DArTseq (SilicoDArT and SNP) markers. After filtering markers for quality-filtering, 14,270 SilicoDArTs and 6484 SNPs were retained and used for genetic diversity, population structure and linkage disequilibrium analyses. The highest number of SilicoDArT and SNP markers mapped on genome A and B compared to genome D. In both marker types, polymorphism information content (PIC) values ranged from 0.1 to 0.5, while > 0.80% of SilicoDArTs and > 0.44% SNPs showed PIC value more than median (0.25%). Un-weighted Neighbor Joining cluster analysis and Bayesian-based model population structure grouped wheat genotypes into three and four clusters, respectively. Principal component analysis and discriminant analysis of principal component results showed highly match with cluster and population structure analysis. Linkage disequilibrium (LD) was more extensive in both marker types, while graphical display of LD decay for both marker types showed that LD declined in the region close to 15 kbp, where r 2-values corresponded to r 2 = 0.16. Overall, our genetic diversity analysis showed high level of variation in studied wheat genotypes, even though there was no relationship between wheat grouping and origins. This might be attributed to admixture level that occurred during long-term natural selection of wheat genotypes in different parts of the world. Highly diverse wheat genotypes used in this study may possess unique genes and are useful sources in breeding programs to improve grain yield and quality.

New Male Specific Markers for Hop and Application in Breeding Program

Male specific DNA sequences were selected from a Diversity Arrays Technology (DArT) mapping study to evaluate their suitability for determination of the sex phenotype among young seedlings in a hop (Humulus lupulus L.) breeding program. Ten male specific DArT markers showed complete linkage with male sex phenotype in three crossing families. Following optimization, four were successfully converted into PCR markers and a multiplex PCR approach for their use was developed. Among 197 plants (97 from the world collection; 100 from three segregating families), 94-100% positive correlation with sex phenotypic data was achieved for the single PCR amplification, whereas the multiplex approach showed 100% correlation. To develop a fast and low-cost method, crude sample multiplex PCR was evaluated in 253 progenies from 14 segregating populations without losing accuracy. The study describes, for the first time, the routine application of molecular markers linked to male sex in an intensive Slovenian hop breeding program. The methods described could be employed for screening of sex at the seedling stage in other hop programs worldwide, thereby saving resources for desirable female plants.

Thatcher wheat line RL6149 carries Lr64 and a second leaf rust resistance gene on chromosome 1DS

The leaf rust resistance gene Lr64 in the Thatcher wheat RL6149 was mapped to chromosome 6AL with SNP and KASP markers and a second leaf rust resistance gene was mapped to chromosome 1DS. RL6149, a near-isogenic line of Thatcher wheat, carries leaf rust resistance gene Lr64 on chromosome arm 6AL. The objective of this study was to develop molecular markers that can be easily used to select wheat lines with Lr64. RL6149 was crossed with Thatcher and F₂ plants derived from a single F₁ plant were advanced to F₆ lines by single seed descent. The 100 F₇ recombinant inbred lines (RIL) were inoculated with two races of P.triticina that differed widely for virulence in order to identify resistant and susceptible RIL. Thirty RIL that differed for resistance and the parental lines were genotyped with the 90 K Infinium iSelect single nucleotide polymorphism (SNP) array to find closely linked markers with Lr64. Seven linked SNPs on chromosome arm 6AL were converted into Kompetitive Allele Specific PCR (KASP) markers that were genotyped on the 100 RIL. A genetic linkage map for the seven KASP markers spanned 19.1 cM on chromosome arm 6AL. KASP marker K-IWB59855 was tightly linked to Lr64. A second unexpected gene for leaf rust resistance also segregated in the F₇ lines. Four KASP markers that spanned 18.6 cM located the gene on chromosome 1DS. The KASP marker K-IWB38437 was tightly linked to the second leaf rust resistance gene.

Interspecific Genetic Differences and Historical Demography in South American Arowanas (Osteoglossiformes, Osteoglossidae, Osteoglossum)

The South American arowanas (Osteoglossiformes, Osteoglossidae, Osteoglossum) are emblematic species widely distributed in the Amazon and surrounding basins. Arowana species are under strong anthropogenic pressure as they are extensively exploited for ornamental and food purposes. Until now, limited genetic and cytogenetic information has been available, with only a few studies reporting to their genetic diversity and population structure. In the present study, cytogenetic and DArTseq-derived single nucleotide polymorphism (SNP) data were used to investigate the genetic diversity of the two Osteoglossum species, the silver arowana O. bicirrhosum, and the black arowana O. ferreirai. Both species differ in their 2n (with 2n = 54 and 56 for O. ferreirai and O. bicirrhosum, respectively) and in the composition and distribution of their repetitive DNA content, consistent with their taxonomic status as different species. Our genetic dataset was coupled with contemporary and paleogeographic niche modeling, to develop concurrent demographic models that were tested against each other with a deep learning approach in O. bicirrhosum. Our genetic results reveal that O. bicirrhosum colonized the Tocantins-Araguaia basin from the Amazon basin about one million years ago. In addition, we highlighted a higher genetic diversity of O. bicirrhosum in the Amazon populations in comparison to those from the Tocantins-Araguaia basin.

Mapping a Resistance Gene to Puccinia graminis f. sp. tritici in the Bread Wheat Cultivar 'Matlabas'

Puccinia graminis f. sp. tritici race TTKSF+ was collected from the South African wheat cultivar 'Matlabas' in 2010. F₂ and F₃ populations derived from a Matlabas × Line 37-07 cross segregated for a single resistance gene to race TTKSF that is avirulent to Matlabas. In screening genomic DNA bulks of susceptible or resistant F₂ plants with simple sequence repeat (SSR) markers, three chromosome arm 2BS markers and one multilocus marker amplified alleles present only in the resistant bulks and Matlabas. Additional 2B-specific SSR markers, incorporating markers spanning regions containing Sr9h, SrWLR, Sr28, and Sr47, were screened in the parental lines and mapped in the F₂ population. Linkage and QTL mapping showed that the gene is located between Xbarc160 in the centromeric region and Xgwm47 on the long arm of chromosome 2B. When 2B-specific SNP markers were mapped, the area of interest was delimited to a 15.3 cM region on chromosome arm 2BL, with XIWA543-HRM and Xgwm47 as flanking loci. Matlabas, Webster, and related Sr9h lines all produced a similar, low infection type to race TTKSF, but were susceptible to race TTKSF+. Phenotypic data and allelic studies suggested that stem rust resistance in Matlabas was derived from an Sr9h source.

Bread wheat: a role model for plant domestication and breeding

BACKGROUND

Bread wheat is one of the most important crops in the world. Its domestication coincides with the beginning of agriculture and since then, it has been constantly under selection by humans. Its breeding has followed millennia of cultivation, sometimes with unintended selection on adaptive traits, and later by applying intentional but empirical selective pressures. For more than one century, wheat breeding has been based on science, and has been constantly evolving due to on farm agronomy and breeding program improvements. The aim of this work is to briefly review wheat breeding, with emphasis on the current advances.

DISCUSSION

Improving yield potential, resistance/tolerance to biotic and abiotic stresses, and baking quality, have been priorities for breeding this cereal, however, new objectives are arising, such as biofortification enhancement. The narrow genetic diversity and complexity of its genome have hampered the breeding progress and the application of biotechnology. Old approaches, such as the introgression from relative species, mutagenesis, and hybrid breeding are strongly reappearing, motivated by an accumulation of knowledge and new technologies. A revolution has taken place regarding the use of molecular markers whereby thousands of plants can be routinely genotyped for thousands of loci. After 13 years, the wheat reference genome sequence and annotation has finally been completed, and is currently available to the scientific community. Transgenics, an unusual approach for wheat improvement, still represents a potential tool, however it is being replaced by gene editing, whose technology along with genomic selection, speed breeding, and high-throughput phenotyping make up the most recent frontiers for future wheat improvement.

FINAL CONSIDERATION

Agriculture and plant breeding are constantly evolving, wheat has played a major role in these processes and will continue through decades to come.

Genome-Wide Genetic Diversity and Population Structure of Tunisian Durum Wheat Landraces Based on DArTseq Technology

Tunisia, being part of the secondary center of diversity for durum wheat, has rich unexploited landraces that are being continuously lost and replaced by high yielding modern cultivars. This study aimed to investigate the genetic diversity and population structure of 196 durum wheat lines issued from landraces collected from Tunisia using Diversity Array Technology sequencing (DArTseq) and to understand possible ways of introduction in comparing them to landraces from surrounding countries. A total of 16,148 polymorphic DArTseq markers covering equally the A and B genomes were effective to assess the genetic diversity and to classify the accessions. Cluster analysis and discriminant analysis of principal components (DAPC) allowed us to distinguish five distinct groups that matched well with the farmer's variety nomenclature. Interestingly, Mahmoudi and Biskri landraces constitute the same gene pool while Jenah Zarzoura constitutes a completely different group. Analysis of molecular variance (AMOVA) showed that the genetic variation was among rather than within the landraces. DAPC analysis of the Tunisian, Mediterranean and West Asian landraces confirmed our previous population structure and showed a genetic similarity between the Tunisian and the North African landraces with the exception of Jenah Zarzoura being the most distant. The genomic characterization of the Tunisian collection will enhance their conservation and sustainable use.

From markers to genome-based breeding in wheat

Recent technological advances in wheat genomics provide new opportunities to uncover genetic variation in traits of breeding interest and enable genome-based breeding to deliver wheat cultivars for the projected food requirements for 2050. There has been tremendous progress in development of whole-genome sequencing resources in wheat and its progenitor species during the last 5 years. High-throughput genotyping is now possible in wheat not only for routine gene introgression but also for high-density genome-wide genotyping. This is a major transition phase to enable genome-based breeding to achieve progressive genetic gains to parallel to projected wheat production demands. These advances have intrigued wheat researchers to practice less pursued analytical approaches which were not practiced due to the short history of genome sequence availability. Such approaches have been successful in gene discovery and breeding applications in other crops and animals for which genome sequences have been available for much longer. These strategies include, (i) environmental genome-wide association studies in wheat genetic resources stored in genbanks to identify genes for local adaptation by using agroclimatic traits as phenotypes, (ii) haplotype-based analyses to improve the statistical power and resolution of genomic selection and gene mapping experiments, (iii) new breeding strategies for genome-based prediction of heterosis patterns in wheat, and (iv) ultimate use of genomics information to develop more efficient and robust genome-wide genotyping platforms to precisely predict higher yield potential and stability with greater precision. Genome-based breeding has potential to achieve the ultimate objective of ensuring sustainable wheat production through developing high yielding, climate-resilient wheat cultivars with high nutritional quality.

Characterization of a dominant mutation for the liguleless trait: Aegilops tauschii liguleless (Lgt)

BACKGROUND

Leaves of Poaceae have a unique morphological feature: they consist of a proximal sheath and a distal blade separated by a ligular region. The sheath provides structural support and protects young developing leaves, whereas the main function of the blade is photosynthesis. The auricles allow the blade to tilt back for optimal photosynthesis and determine the angle of a leaf, whereas the ligule protects the stem from the entry of water, microorganisms, and pests. Liguleless variants have an upright leaf blade that wraps around the culm. Research on liguleless mutants of maize and other cereals has led to identification of genes that are involved in leaf patterning and differentiation.

RESULTS

We characterized an induced liguleless mutant (LM) of Aegilops tauschii Coss., a donor of genome D of bread wheat Triticum aestivum L.. The liguleless phenotype of LM is under dominant monogenic control (Lg^t). To determine precise position of Lg^t on the Ae. tauschii genetic map, highly saturated genetic maps were constructed containing 887 single-nucleotide polymorphism (SNP) markers derived via diversity arrays technology (DArT)seq. The Lg^t gene was mapped to chromosome 5DS. Taking into account coordinates of the SNP markers, flanking Lg^t, on the pseudomolecule 5D, a chromosomal region that contains this gene was determined, and a list of candidate genes was identified. Morphological features of the LM phenotype suggest that Lg^t participates in the control of leaf development, mainly, in leaf proximal-distal patterning, and its dominant mutation causes abnormal ligular region but does not affect reproductive development.

CONCLUSIONS

Here we report characterization of a liguleless Ae. tauschii mutant, whose phenotype is under control of a dominant mutation of Lg^t. The dominant mode of inheritance of the liguleless trait in a Triticeae species is reported for the first time. The position of the Lg^t locus on chromosome 5DS allowed us to identify a list of candidate genes. This list does not contain Ae. tauschii orthologs of any well-characterized cereal genes whose mutations cause liguleless phenotypes. Thus, the characterized Lg^t mutant represents a new model for further investigation of plant leaf patterning and differentiation.