The USDA barley core collection: genetic diversity, population structure, and potential for genome-wide association studies

Root system architecture variation among barley (Hordeum vulgare L.) accessions at seedling stage under soil acidity condition

MAIN CONCLUSION

Soil acidity in Ethiopian highlands impacts barley production, affecting root system architecture. Study on 300 accessions showed significant trait variability, with potential for breeding enhancement. Soil acidity poses a significant challenge to crop production in the highland regions of Ethiopia, particularly impacting barley, a crucial staple crop. This acidity serves as a key stressor affecting the root system architecture (RSA) of this crop. Hence, the objective of this study was to assess the RSA traits variability under acidic soil conditions using 300 barley accessions in a greenhouse experiment. The analysis of variance indicated substantial variations among the accessions across all traits studied. The phenotypic coefficient of variation ranged from 24.4% for shoot dry weight to 11.1% for root length, while the genotypic coefficient variation varied between 18.83 and 9.2% for shoot dry weight and root length, respectively. The broad-sense heritability ranged from 36.7% for leaf area to 69.9% for root length, highlighting considerable heritability among multiple traits. The genetic advances as a percent of the mean ranged from 13.63 to 29.9%, suggesting potential for enhancement of these traits through breeding efforts. Principal component analysis and cluster analysis grouped the genotypes into two major clusters, each containing varying numbers of genotypes with contrasting traits. This diverse group presents an opportunity to access a wide range of potential parent candidates to enhance genetic variablity in breeding programs. The Pearson correlation analysis revealed significant negative associations between root angle (RA) and other RSA traits. This helps indirect selection of accessions for further improvement in soil acidity. In conclusion, this study offers valuable insights into the RSA characteristics of barley in acidic soil conditions, aiding in the development of breeding strategies to enhance crop productivity in acidic soil environments.

Genome-Wide Association Study of Host Resistance to Hessian Fly in Barley

The Hessian fly (HF), Mayetiola destructor (Diptera: Cecidomyiidae), is one of the most devastating insect pests of cereals including wheat, barley, and rye. Although wheat is the preferred host for HF, this continuously evolving pest has been emerging as a threat to barley production. However, characterization and identification of genetic resistance to HF has not been conducted in barley. In the present study, we used a genome-wide association study (GWAS) to identify barley resistance loci to HF using a geographically diverse set of 234 barley accessions. The results showed that around 90% of barley lines were highly susceptible, indicating a significant vulnerability to HF in barley, and a total of 29 accessions were resistant, serving as potential resistance resources. GWAS with a mixed linear model revealed two marker-trait associations, both on chromosome 4H. The resistance loci and associated markers will facilitate barley improvement and development for breeders. In addition, our results are fundamental for genetic studies to understand the HF resistance mechanism in barley.

Unlocking the genetic diversity and population structure of the newly introduced two-row spring European HerItage Barley collecTion (ExHIBiT)

In the last century, breeding programs have traditionally favoured yield-related traits, grown under high-input conditions, resulting in a loss of genetic diversity and an increased susceptibility to stresses in crops. Thus, exploiting understudied genetic resources, that potentially harbour tolerance genes, is vital for sustainable agriculture. Northern European barley germplasm has been relatively understudied despite its key role within the malting industry. The European Heritage Barley collection (ExHIBiT) was assembled to explore the genetic diversity in European barley focusing on Northern European accessions and further address environmental pressures. ExHIBiT consists of 363 spring-barley accessions, focusing on two-row type. The collection consists of landraces (~14%), old cultivars (~18%), elite cultivars (~67%) and accessions with unknown breeding history (~1%), with 70% of the collection from Northern Europe. The population structure of the ExHIBiT collection was subdivided into three main clusters primarily based on the accession's year of release using 26,585 informative SNPs based on 50k iSelect single nucleotide polymorphism (SNP) array data. Power analysis established a representative core collection of 230 genotypically and phenotypically diverse accessions. The effectiveness of this core collection for conducting statistical and association analysis was explored by undertaking genome-wide association studies (GWAS) using 24,876 SNPs for nine phenotypic traits, four of which were associated with SNPs. Genomic regions overlapping with previously characterised flowering genes (HvZTLb) were identified, demonstrating the utility of the ExHIBiT core collection for locating genetic regions that determine important traits. Overall, the ExHIBiT core collection represents the high level of untapped diversity within Northern European barley, providing a powerful resource for researchers and breeders to address future climate scenarios.

Deciphering the genetic diversity and population structure of wild barley germplasm against corn leaf aphid, Rhopalosiphum maidis (Fitch)

Corn-leaf aphid (CLA-Rhopalosiphum maidis) is a major insect pest of barley (Hordeum vulgare) causing yield loss upto 30% under severe infestation. Keeping in view of the availability of very few sources of CLA resistance in barley, the present investigation was framed to assess the genetic diversity and population structure of 43 wild barley (H. vulgare subsp. spontaneum) genotypes using eight microsatellite markers against R. maidis. Three statistical methods viz. multivariate-hierarchical clustering, Bayesian clustering and PCoA, unanimously grouped genotypes into three subpopulations (K = 3) with 25.58% (SubPop1-Red), 39.53% (SubPop2-Green) and 34.88% (SubPop3-Blue) genotypes including admixtures. Based on Q ≥ 66.66%, 37.20% genotypes formed a superficial "Mixed/Admixture" subpopulation. All polymorphic SSR markers generated 36 alleles, averaging to 4.5 alleles/locus (2-7 range). The PIC and H were highest in MS31 and lowest in MS28, with averages of 0.66 and 0.71. MAF and mean genetic diversity were 0.16 and 89.28%, respectively. All these parameters indicated the presence of predominant genetic diversity and population structure amongst the studied genotypes. Based on AII, only 6 genotypes were found to be R. maidis resistant. SubPop3 had 91.66% (11) of the resistant or moderately resistant genotypes. SubPop3 also had the most pure genotypes (11), the least aphid infestation (8.78), and the highest GS (0.88), indicating its suitability for future R. maidis resistance breeding initiatives.

Amand P, Diggle AJ, Rengel Z. Genetic characterization of root architectural traits in barley (Hordeum vulgare L.) using SNP markers

Increasing attention is paid to providing new tools to breeders for targeted breeding for specific root traits that are beneficial in low-fertility, drying soils; however, such information is not available for barley (Hordeum vulgare L.). A panel of 191 barley accessions (originating from Australia, Europe, and Africa) was phenotyped for 26 root and shoot traits using the semi-hydroponic system and genotyped using 21 062 high-quality single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing (GBS). The population structure analysis of the barley panel identified six distinct groups. We detected 1199 significant (P<0.001) marker-trait associations (MTAs) with r2 values up to 0.41. The strongest MTAs were found for root diameter in the top 20 cm and the longest root length. Based on the physical locations of these MTAs in the barley reference genome, we identified 37 putative QTLs for the root traits, and three QTLs for shoot traits, with nine QTLs located in the same physical regions. The genomic region 640-653 Mb on chromosome 7H was significant for five root length-related traits, where 440 annotated genes were located. The putative QTLs for various root traits identified in this study may be useful for genetic improvement regarding the adaptation of new barley cultivars to suboptimal environments and abiotic stresses.

Genome-wide association and targeted transcriptomic analyses reveal loci and candidate genes regulating preharvest sprouting in barley

KEY MESSAGE

Genome-wide association study of diverse barley genotypes identified loci, single nucleotide polymorphisms and candidate genes that control seed dormancy and therefore enhance resistance to preharvest sprouting. Preharvest sprouting (PHS) causes significant yield and quality loss in barley and it is strongly associated with the level of seed dormancy. This study performed genome-wide association study using a collection of 255 diverse barley genotypes grown over four environments to identify loci controlling dormancy/PHS. Our phenotypic analysis revealed substantial variation in germination index/dormancy levels among the barley genotypes. Marker-trait association and linkage disequilibrium (LD) decay analyses identified 16 single nucleotide polymorphisms (SNPs) and two QTLs associated with dormancy/PHS, respectively, on chromosome 3H and 5H explaining 6.9% to 11.1% of the phenotypic variation. QTL.5H consist of 14 SNPs of which 12 SNPs satisfy the FDR threshold of α = 0.05, and it may represent the SD2 locus. The QTL on 3H consists of one SNP that doesn't satisfy FDR (α = 0.05). Genes harbouring the significant SNPs were analyzed for their expression pattern in the seeds of selected dormant and non-dormant genotypes. Of these genes, HvRCD1, HvPSRP1 and HvF3H exhibited differential expression between the dormant and non-dormant seed samples, suggesting their role in controlling seed dormancy/PHS. Three SNPs located within the differentially expressed genes residing in QTL.5H explained considerable phenotypic variation (≥ 8.6%), suggesting their importance in regulating PHS resistance. Analysis of the SNP marker data in QTL.5H identified a haplotype for PHS resistance. Overall, the study identified loci, SNPs and candidate genes that control dormancy and therefore play important roles in enhancing PHS resistance in barley through marker-assisted breeding.

Genome-wide association study of resistance to anthracnose in pepper (Capsicum chinense) germplasm

BACKGROUND

Anthracnose is a fungal disease caused by Colletotrichum spp. that has a significant impact on worldwide pepper production. Colletotrichum scovillei is the most common pathogenic anthracnose-causing species in the Republic of Korea.

RESULTS

The resistances of 197 pepper (Capsicum chinense) accessions deposited in Korea's National Agrobiodiversity Center were evaluated for their response against the virulent pathogens Colletotrichum acutatum isolate 'KSCa-1' and C. scovillei isolate 'Hana') in the field and in vitro methods for three consecutive years (2018 to 2020). The severity of the disease was recorded and compared between inoculation methods. Six phenotypically resistant pepper accessions were selected based on three years of disease data. All of the selected resistant pepper accessions outperformed the control resistant pepper in terms of resistance (PI 594,137). A genome-wide association study (GWAS) was carried out to identify single nucleotide polymorphisms (SNPs) associated with anthracnose resistance. An association analysis was performed using 53,518 SNPs and the disease score of the 2020 field and in vitro experiment results. Both field and in vitro experiments revealed 25 and 32 significantly associated SNPs, respectively. These SNPs were found on all chromosomes except Ch06 and Ch07 in the field experiment, whereas in the in vitro experiment they were found on all chromosomes except Ch04 and Ch11.

CONCLUSION

In this study, six resistant C. chinense accessions were selected. Additionally, in this study, significantly associated SNPs were found in a gene that codes for a protein kinase receptor, such as serine/threonine-protein kinase, and other genes that are known to be involved in disease resistance. This may strengthen the role of these genes in the development of anthracnose resistance in Capsicum spp. As a result, the SNPs discovered to be strongly linked in this study can be used to identify a potential marker for selecting pepper material resistant to anthracnose, which will assist in the development of resistant varieties.

Population structure analysis and genome-wide association study of a hexaploid oat landrace and cultivar collection

INTRODUCTION

Oat (Avena sativa L.) is an important cereal crop grown worldwide for grain and forage, owing to its high adaptability to diverse environments. However, the genetic and genomics research of oat is lagging behind that of other staple cereal crops.

METHODS

In this study, a collection of 288 oat lines originating worldwide was evaluated using 2,213 single nucleotide polymorphism (SNP) markers obtained from an oat iSelect 6K-beadchip array to study its genetic diversity, population structure, and linkage disequilibrium (LD) as well as the genotype-phenotype association for hullessness and lemma color.

RESULTS

The average gene diversity and polymorphic information content (PIC) were 0.324 and 0.262, respectively. The first three principal components (PCs) accounted for 30.33% of the genetic variation, indicating that the population structure of this panel of oat lines was stronger than that reported in most previous studies. In addition, accessions could be classified into two subpopulations using a Bayesian clustering approach, and the clustering pattern of accessions was closely associated with their region of origin. Additionally, evaluation of LD decay using 2,143 mapped markers revealed that the intrachromosomal whole-genome LD decayed rapidly to a critical r² value of 0.156 for marker pairs separated by a genetic distance of 1.41 cM. Genome-wide association study (GWAS) detected six significant associations with the hullessness trait. Four of these six markers were located on the Mrg21 linkage group between 194.0 and 205.7 cM, while the other two significant markers mapped to Mrg05 and Mrg09. Three significant SNPs, showing strong association with lemma color, were located on linkage groups Mrg17, Mrg18, and Mrg20.

DISCUSSION

Our results discerned relevant patterns of genetic diversity, population structure, and LD among members of a worldwide collection of oat landraces and cultivars proposed to be 'typical' of the Qinghai-Tibetan Plateau. These results have important implications for further studies on association mapping and practical breeding in high-altitude oat.

Single nucleotide polymorphism (SNP) markers for genetic diversity and population structure study in Ethiopian barley (Hordeum vulgare L.) germplasm

BACKGROUND

High-density single nucleotide polymorphisms (SNPs) are the most abundant and robust form of genetic variants and hence make highly favorable markers to determine the genetic diversity and relationship, enhancing the selection of breeding materials and the discovery of novel genes associated with economically important traits. In this study, a total of 105 barley genotypes were sampled from various agro-ecologies of Ethiopia and genotyped using 10 K single nucleotide polymorphism (SNP) markers. The refined dataset was used to assess genetic diversity and population structure.

RESULTS

The average gene diversity was 0.253, polymorphism information content (PIC) of 0.216, and minor allelic frequency (MAF) of 0.118 this revealed a high genetic variation in barley genotypes. The genetic differentiation also showed the existence of variations, ranging from 0.019 to 0.117, indicating moderate genetic differentiation between barley populations. Analysis of molecular variance (AMOVA) revealed that 46.43% and 52.85% of the total genetic variation occurred within the accessions and populations, respectively. The heat map, principal components and population structure analysis further confirm the presence of four distinct clusters.

CONCLUSIONS

This study confirmed that there is substantial genetic variation among the different barley genotypes. This information is useful in genomics, genetics and barley breeding.

Wild barley (Hordeum spontaneum) and landraces (Hordeum vulgare) from Turkey contain an abundance of novel Rhynchosporium commune resistance loci

Rhynchosporium commune is a globally devastating pathogen of barley. Wild and landrace barley are underutilized, however, contain an abundance of loci that can be used as potential sources of resistance. Rhynchosporium commune, the causal agent of the disease scald or leaf blotch of barley, is a hemibiotrophic fungal pathogen of global importance, responsible for yield losses ranging from 30 to 40% on susceptible varieties. To date, over 150 resistance loci have been characterized in barley. However, due to the suspected location of the R. commune host jump in Europe, European germplasm has been the primary source used to screen for R. commune resistance leaving wild (Hordeum spontaneum) and landrace (H. vulgare) barley populations from the center of origin largely underutilized. A diverse population consisting of 94 wild and 188 barley landraces from Turkey were genotyped using PCR-GBS amplicon sequencing and screened with six Turkish R. commune isolates. The isolates were collected from distinct geographic regions of Turkey with two from the Aegean region, two from central Turkey and two from the Fertile Crescent region. The data set was utilized for association mapping analysis with a total of 21 loci identified, of which 12 were novel, indicating that these diverse primary barley gene pools contain an abundance of novel R. commune resistances that could be utilized for resistance breeding.

Genetic diversity and population structure of barley landraces from Southern Ethiopia's Gumer district: Utilization for breeding and conservation

Barley is the fifth most important food crop in Ethiopia. The genetic relationship and population structure studies of barley are limited to gene bank collections. Therefore, this study fills a gap by investigating the selection, consumption, economic value, genetic diversity, and population structure of farm-collected barley from the Gumer district of the Gurage Zone, which has received little attention. The information on the use of barley in the study area was collected using semi-structured interviews and questionnaires. 124 households of 11 kebeles, the smallest community unit, were interviewed. Barley landraces collected were compared with those collected from Japan, the United States (USA), and other Ethiopian locations. Illumina iSelect (50K genotyping platform) was used to identify single nucleotide polymorphisms (SNP) (20,367). Thirty landraces were found in Gumer. Burdaenadenber had the highest on-farm Shannon index estimate (2.0), followed by Aselecha (1.97) and Enjefo (1.95). Aselecha and Fetazer had the highest (44%) and the lowest (29%) richness values, respectively. High and low Simpson index values were found in Aselecha (84%) and Wulbaragenateretero (79%), respectively. The neighbor-joining tree revealed that Gumer landraces formed a separate subcluster with a common ancestral node; a sister subcluster contained barley landraces from Japan. According to the population structure analysis, barley landraces from Gumer differed from Japan and the United States. The principal component analysis revealed that US barley was the most distant group from Gumer barley. The markers' allele frequencies ranged from 0.10 to 0.50, with an average value of 0.28. The mean values of Nei's gene diversity (0.38) and the polymorphic information content (0.30) indicated the presence of high genetic diversity in the samples. The clustering of accessions was not based on geographic origin. Significant genetic diversity calls for additional research and analysis of local barley diversity because the selection and use of barley in Ethiopia would have been affected by the preference of ethnic groups.

CGIAR Barley Breeding Toolbox: A diversity panel to facilitate breeding and genomic research in the developing world

Breeding programs in developing countries still cannot afford the new genotyping technologies, hindering their research. We aimed to assemble an Association Mapping panel to serve as CGIAR Barley Breeding Toolbox (CBBT), especially for the Developing World. The germplasm had to be representative of the one grown in the Developing World; with high genetic variability and be of public domain. For it, we genotyped with the Infinium iSelect 50K chip, a Global Barley Panel (GBP) of 530 genotypes representing a wide range of row-types, end-uses, growth habits, geographical origins and environments. 40,342 markers were polymorphic with an average polymorphism information content of 0.35 and 66% of them exceeding 0.25. The analysis of the population structure identified 8 subpopulations mostly linked to geographical origin, four of them with significant ICARDA origin. The 16 allele combinations at 4 major flowering genes (HvVRN-H3, HvPPD-H1, HvVRN-H1 and HvCEN) explained 11.07% genetic variation and were linked to the geographic origins of the lines. ICARDA material showed the widest diversity as revealed by the highest number of polymorphic loci (99.76% of all polymorphic SNPs in GBP), number of private alleles and the fact that ICARDA lines were present in all 8 subpopulations and carried all 16 allelic combinations. Due to their genetic diversity and their representativity of the germplasm adapted to the Developing World, ICARDA-derived lines and cultivated landraces were pre-selected to form the CBBT. Using the Mean of Transformed Kinships method, we assembled a panel capturing most of the allelic diversity in the GBP. The CBBT (N=250) preserves good balance between row-types and good representation of both phenology allelic combinations and subpopulations of the GBP. The CBBT and its genotypic data is available to researchers worldwide as a collaborative tool to underpin the genetic mechanisms of traits of interest for barley cultivation.

Novel Single-Nucleotide Variants for Morpho-Physiological Traits Involved in Enhancing Drought Stress Tolerance in Barley

Barley (Hordeum vulgare L.) thrives in the arid and semi-arid regions of the world; nevertheless, it suffers large grain yield losses due to drought stress. A panel of 426 lines of barley was evaluated in Egypt under deficit (DI) and full irrigation (FI) during the 2019 and 2020 growing seasons. Observations were recorded on the number of days to flowering (NDF), total chlorophyll content (CH), canopy temperature (CAN), grain filling duration (GFD), plant height (PH), and grain yield (Yield) under DI and FI. The lines were genotyped using the 9K Infinium iSelect single nucleotide polymorphisms (SNP) genotyping platform, which resulted in 6913 high-quality SNPs. In conjunction with the SNP markers, the phenotypic data were subjected to a genome-wide association scan (GWAS) using Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK). The GWAS results indicated that 36 SNPs were significantly associated with the studied traits under DI and FI. Furthermore, eight markers were significant and common across DI and FI water regimes, while 14 markers were uniquely associated with the studied traits under DI. Under DI and FI, three (11_10326, 11_20042, and 11_20170) and five (11_20099, 11_10326, 11_20840, 12_30298, and 11_20605) markers, respectively, had pleiotropic effect on at least two traits. Among the significant markers, 24 were annotated to known barley genes. Most of these genes were involved in plant responses to environmental stimuli such as drought. Overall, nine of the significant markers were previously reported, and 27 markers might be considered novel. Several markers identified in this study could enable the prediction of barley accessions with optimal agronomic performance under DI and FI.

Host and pathogen genetics reveal an inverse gene-for-gene association in the P. teres f. maculata-barley pathosystem

Pathogen and host genetics were used to uncover an inverse gene-for-gene interaction where virulence genes from the pathogen Pyrenophora teres f. maculata target barley susceptibility genes, resulting in disease. Although models have been proposed to broadly explain how plants and pathogens interact and coevolve, each interaction evolves independently, resulting in various scenarios of host manipulation and plant defense. Spot form net blotch is a foliar disease of barley caused by Pyrenophora teres f. maculata. We developed a barley population (Hockett × PI 67381) segregating for resistance to a diverse set of P. teres f. maculata isolates. Quantitative trait locus analysis identified major loci on barley chromosomes (Chr) 2H and 7H associated with resistance/susceptibility. Subsequently, we used avirulent and virulent P. teres f. maculata isolates to develop a pathogen population, identifying two major virulence loci located on Chr1 and Chr2. To further characterize this host-pathogen interaction, progeny from the pathogen population harboring virulence alleles at either the Chr1 or Chr2 locus was phenotyped on the Hockett × PI 67381 population. Progeny harboring only the Chr1 virulence allele lost the barley Chr7H association but maintained the 2H association. Conversely, isolates harboring only the Chr2 virulence allele lost the barley Chr2H association but maintained the 7H association. Hockett × PI 67381 F₂ individuals showed susceptible/resistant ratios not significantly different than 15:1 and results from F₂ inoculations using the single virulence genotypes were not significantly different from a 3:1 (S:R) ratio, indicating two dominant susceptibility genes. Collectively, this work shows that P. teres f. maculata virulence alleles at the Chr1 and Chr2 loci are targeting the barley 2H and 7H susceptibility alleles in an inverse gene-for-gene manner to facilitate colonization.

Genomic-environmental associations in wild cranberry (Vaccinium macrocarpon Ait.)

Understanding the genetic basis of local adaptation in natural plant populations, particularly crop wild relatives, may be highly useful for plant breeding. By characterizing genetic variation for adaptation to potentially stressful environmental conditions, breeders can make targeted use of crop wild relatives to develop cultivars for novel or changing environments. This is especially appealing for improving long-lived woody perennial crops such as the American cranberry (Vaccinium macrocarpon Ait.), the cultivation of which is challenged by biotic and abiotic stresses. In this study, we used environmental association analyses in a collection of 111 wild cranberry accessions to identify potentially adaptive genomic regions for a range of bioclimatic and soil conditions. We detected 126 significant associations between SNP marker loci and environmental variables describing temperature, precipitation, and soil attributes. Many of these markers tagged genes with functional annotations strongly suggesting a role in adaptation to biotic or abiotic conditions. Despite relatively low genetic variation in cranberry, our results suggest that local adaptation to divergent environments is indeed present, and the identification of potentially adaptive genetic variation may enable a selective use of this germplasm for breeding more stress-tolerant cultivars.

A reference-guided TILLING by amplicon-sequencing platform supports forward and reverse genetics in barley

Barley is a diploid species with a genome smaller than those of other members of the Triticeae tribe, making it an attractive model for genetic studies in Triticeae crops. The recent development of barley genomics has created a need for a high-throughput platform to identify genetically uniform mutants for gene function investigations. In this study, we report an ethyl methanesulfonate (EMS)-mutagenized population consisting of 8525 M₃ lines in the barley landrace "Hatiexi" (HTX), which we complement with a high-quality de novo assembly of a reference genome for this genotype. The mutation rate within the population ranged from 1.51 to 4.09 mutations per megabase, depending on the treatment dosage of EMS and the mutation discrimination platform used for genotype analysis. We implemented a three-dimensional DNA pooling strategy combined with multiplexed amplicon sequencing to create a highly efficient and cost-effective TILLING (targeting induced locus lesion in genomes) platform in barley. Mutations were successfully identified from 72 mixed amplicons within a DNA pool containing 64 individual mutants and from 56 mixed amplicons within a pool containing 144 individuals. We discovered abundant allelic mutants for dozens of genes, including the barley Green Revolution contributor gene Brassinosteroid insensitive 1 (BRI1). As a proof of concept, we rapidly determined the causal gene responsible for a chlorotic mutant by following the MutMap strategy, demonstrating the value of this resource to support forward and reverse genetic studies in barley.

A pan-genome data structure induced by pooled sequencing facilitates variant mining in heterogeneous germplasm

Valuable genetic variation lies unused in gene banks due to the difficulty of exploiting heterogeneous germplasm accessions. Advances in molecular breeding, including transgenics and genome editing, present the opportunity to exploit hidden sequence variation directly. Here we describe the pan-genome data structure induced by whole-genome sequencing of pooled individuals from wild populations of Patellifolia spp., a source of disease resistance genes for the related crop species sugar beet (Beta vulgaris). We represent the pan-genome as a map of reads from pooled sequencing of a heterogeneous population sample to a reference genome, plus a BLAST data base of the mapped reads. We show that this basic data structure can be queried by reference genome position or homology to identify sequence variants present in the wild relative, at genes of agronomic interest in the crop, a process known as allele or variant mining. Further we demonstrate the possibility of cataloging variants in all Patellifolia genomic regions that have corresponding single copy orthologous regions in sugar beet. The data structure, termed a "pooled read archive," can be produced, altered, and queried using standard tools to facilitate discovery of agronomically-important sequence variation.

Supplementary information

The online version contains supplementary material available at 10.1007/s11032-022-01308-6.

The genetic diversity of Ethiopian barley genotypes in relation to their geographical origin

Ethiopia is recognized as a center of diversity for barley, and its landraces are known for the distinct genetic features compared to other barley collections. The genetic diversity of Ethiopian barley likely results from the highly diverse topography, altitude, climate conditions, soil types, and farming systems. To get detailed information on the genetic diversity a panel of 260 accessions, comprising 239 landraces and 21 barley breeding lines, obtained from the Ethiopian biodiversity institute (EBI) and the national barley improvement program, respectively were studied for their genetic diversity using the 50k iSelect single nucleotide polymorphism (SNP) array. A total of 983 highly informative SNP markers were used for structure and diversity analysis. Three genetically distinct clusters were obtained from the structure analysis comprising 80, 71, and 109 accessions, respectively. Analysis of molecular variance (AMOVA) revealed the presence of higher genetic variation (89%) within the clusters than between the clusters (11%), with moderate genetic differentiation (PhiPT = 0.11) and five accessions were detected as first-generation migrants using Monte Carlo resampling methods. The Mantel test revealed that the genetic distance between accessions is poorly associated with their geographical distance. Despite the observed weak correlation between geographic distance and genetic differentiation, for some regions like Gonder, Jimma, Gamo-Gofa, Shewa, and Welo, more than 50% of the landraces derived from these regions are assigned to one of the three clusters.

Core Collections: Is There Any Value for Cotton Breeding?

Global plant breeding activities are reliant on the available genetic variation held in extant varieties and germplasm collections. Throughout the mid- to late 1900s, germplasm collecting efforts were prioritized for breeding programs to archive precious material before it disappeared and led to the development of the numerous large germplasm resources now available in different countries. In recent decades, however, the maintenance and particularly the expansion of these germplasm resources have come under threat, and there has been a significant decline in investment in further collecting expeditions, an increase in global biosecurity restrictions, and restrictions placed on the open exchange of some commercial germplasm between breeders. The large size of most genebank collections, as well as constraints surrounding the availability and reliability of accurate germplasm passport data and physical or genetic characterization of the accessions in collections, limits germplasm utilization by plant breeders. To overcome these constraints, core collections, defined as a representative subset of the total germplasm collection, have gained popularity. Core collections aim to increase germplasm utilization by containing highly characterized germplasm that attempts to capture the majority of the variation in a whole collection. With the recent availability of many new genetic tools, the potential to unlock the value of these resources can now be realized. The Commonwealth Scientific and Industrial Research Organisation (CSIRO) cotton breeding program supplies 100% of the cotton cultivars grown in Australia. The program is reliant on the use of plant genetic resources for the development of improved cotton varieties to address emerging challenges in pest and disease resistance as well as the global changes occurring in the climate. Currently, the CSIRO germplasm collection is actively maintained but underutilized by plant breeders. This review presents an overview of the Australian cotton germplasm resources and discusses the appropriateness of a core collection for cotton breeding programs.

A high-resolution consensus linkage map for barley based on GBS-derived genotypes

As genotyping-by-sequencing (GBS) is widely used in barley genetic studies, the translation of the physical position of GBS-derived SNPs into accurate genetic positions has become relevant. The main aim of this study was to develop a high-resolution consensus linkage map based on GBS-derived SNPs. The construction of this integrated map involved 11 bi-parental populations composed of 3743 segregating progenies. We adopted a uniform set of SNP-calling and filtering conditions to identify 50 875 distinct SNPs segregating in at least one population. These SNPs were grouped into 18 580 non-redundant SNPs (bins). The resulting consensus linkage map spanned 1050.1 cM, providing an average density of 17.7 bins and 48.4 SNPs per cM. The consensus map is characterized by the absence of large intervals devoid of marker coverage (significant gaps), the largest interval between bins was only 3.7 cM and the mean distance between adjacent bins was 0.06 cM. This high-resolution linkage map will contribute to several applications in genomic research, such as providing useful information on the recombination landscape for QTLs/genes identified via GWAS or ensuring a uniform distribution of SNPs when developing low-cost genotyping tools offering a limited number of markers.

Identification and Validation of a Core Single-Nucleotide Polymorphism Marker Set for Genetic Diversity Assessment, Fingerprinting Identification, and Core Collection Development in Bottle Gourd

Germplasm collections are indispensable resources for the mining of important genes and variety improvement. To preserve and utilize germplasm collections in bottle gourd, we identified and validated a highly informative core single-nucleotide polymorphism (SNP) marker set from 1,100 SNPs. This marker set consisted of 22 uniformly distributed core SNPs with abundant polymorphisms, which were established to have strong representativeness and discriminatory power based on analyses of 206 bottle gourd germplasm collections and a multiparent advanced generation inter-cross (MAGIC) population. The core SNP markers were used to assess genetic diversity and population structure, and to fingerprint important accessions, which could provide an optimized procedure for seed authentication. Furthermore, using the core SNP marker set, we developed an accessible core population of 150 accessions that represents 100% of the genetic variation in bottle gourds. This core population will make an important contribution to the preservation and utilization of bottle gourd germplasm collections, cultivar identification, and marker-assisted breeding.

Adult resistance genes to barley powdery mildew confer basal penetration resistance associated with broad-spectrum resistance

Powdery mildew isa major disease of barley (Hordeum vulgare L.) for which breeders have traditionally relied on dominant, pathogen race-specific resistance genes for genetic control. Directional selection pressures in extensive monocultures invariably result in such genes being overcome as the pathogen mutates to evade recognition. This has led to a widespread reliance on fungicides and a single broad-spectrum recessive resistance provided by the mlo gene. The range of resistance genes and alleles found in wild crop relatives and landraces has been reduced in agricultural cultivars through an erosion of genetic diversity during domestication and selective breeding. Three novel major-effect adult plant resistance (APR) genes from landraces, designated Resistance to Blumeria graminis f. sp. hordei (Rbgh1 to Rbgh3), were identified in the terminal regions of barley chromosomes 5HL, 7HS, and 1HS, respectively. The phenotype of the new APR genes showed neither pronounced penetration resistance, nor the spontaneous necrosis and mesophyll cell death typical of mlo resistance, nor a whole epidermal cell hypersensitive response, typical of race-specific resistance. Instead, resistance was localized to the site of attempted penetration in an epidermal cell and was associated with cell wall appositions and cytosolic vesicle-like bodies, and lacked strong induction of reactive oxygen species. The APR genes exhibited differences in vesicle-like body sizes, their distribution, and the extent of localized 3,3-diaminobenzidine staining in individual doubled haploid lines. The results revealed a set of unique basal penetration resistance genes that offer opportunities for combining different resistance mechanisms in breeding programs for robust mildew resistance.

Genome-wide association mapping of Pyrenophora teres f. maculata and Pyrenophora teres f. teres resistance loci utilizing natural Turkish wild and landrace barley populations

Unimproved landraces and wild relatives of crops are sources of genetic diversity that were lost post domestication in modern breeding programs. To tap into this rich resource, genome-wide association studies in large plant genomes have enabled the rapid genetic characterization of desired traits from natural landrace and wild populations. Wild barley (Hordeum spontaneum), the progenitor of domesticated barley (Hordeum vulgare), is dispersed across Asia and North Africa, and has co-evolved with the ascomycetous fungal pathogens Pyrenophora teres f. teres and P. teres f. maculata, the causal agents of the diseases net form of net blotch and spot form of net blotch, respectively. Thus, these wild and local adapted barley landraces from the region of origin of both the host and pathogen represent a diverse gene pool to identify new sources of resistance, due to millions of years of co-evolution. The barley—P. teres pathosystem is governed by complex genetic interactions with dominant, recessive, and incomplete resistances and susceptibilities, with many isolate-specific interactions. Here, we provide the first genome-wide association study of wild and landrace barley from the Fertile Crescent for resistance to both forms of P. teres. A total of 14 loci, four against P. teres f. maculata and 10 against P. teres f. teres, were identified in both wild and landrace populations, showing that both are genetic reservoirs for novel sources of resistance. We also highlight the importance of using multiple algorithms to both identify and validate additional loci.

Population Structure and Genetic Diversity of Two-Rowed Barley Accessions from Kazakhstan Based on SNP Genotyping Data

The genetic relationship and population structure of two-rowed barley accessions from Kazakhstan were assessed using single-nucleotide polymorphism (SNP) markers. Two different approaches were employed in the analysis: (1) the accessions from Kazakhstan were compared with barley samples from six different regions around the world using 1955 polymorphic SNPs, and (2) 94 accessions collected from six breeding programs from Kazakhstan were studied using 5636 polymorphic SNPs using a 9K Illumina Infinium assay. In the first approach, the neighbor-joining tree showed that the majority of the accessions from Kazakhstan were grouped in a separate subcluster with a common ancestral node; there was a sister subcluster that comprised mainly barley samples that originated in Europe. The Pearson’s correlation analysis suggested that Kazakh accessions were genetically close to samples from Africa and Europe. In the second approach, the application of the STRUCTURE package using 5636 polymorphic SNPs suggested that Kazakh barley samples consisted of five subclusters in three major clusters. The principal coordinate analysis plot showed that, among six breeding origins in Kazakhstan, the Krasnovodopad (KV) and Karaganda (KA) samples were the most distant groups. The assessment of the pedigrees in the KV and KA samples showed that the hybridization schemes in these breeding stations heavily used accessions from Ethiopia and Ukraine, respectively. The comparative analysis of the KV and KA samples allowed us to identify 214 SNPs with opposite allele frequencies that were tightly linked to 60 genes/gene blocks associated with plant adaptation traits, such as the heading date and plant height. The identified SNP markers can be efficiently used in studies of barley adaptation and deployed in breeding projects to develop new competitive cultivars.

Major flowering time genes of barley: allelic diversity, effects, and comparison with wheat

This review summarizes the allelic series, effects, interactions between genes and with the environment, for the major flowering time genes that drive phenological adaptation of barley. The optimization of phenology is a major goal of plant breeding addressing the production of high-yielding varieties adapted to changing climatic conditions. Flowering time in cereals is regulated by genetic networks that respond predominately to day length and temperature. Allelic diversity at these genes is at the basis of barley wide adaptation. Detailed knowledge of their effects, and genetic and environmental interactions will facilitate plant breeders manipulating flowering time in cereal germplasm enhancement, by exploiting appropriate gene combinations. This review describes a catalogue of alleles found in QTL studies by barley geneticists, corresponding to the genetic diversity at major flowering time genes, the main drivers of barley phenological adaptation: VRN-H1 (HvBM5A), VRN-H2 (HvZCCTa-c), VRN-H3 (HvFT1), PPD-H1 (HvPRR37), PPD-H2 (HvFT3), and eam6/eps2 (HvCEN). For each gene, allelic series, size and direction of QTL effects, interactions between genes and with the environment are presented. Pleiotropic effects on agronomically important traits such as grain yield are also discussed. The review includes brief comments on additional genes with large effects on phenology that became relevant in modern barley breeding. The parallelisms between flowering time allelic variation between the two most cultivated Triticeae species (barley and wheat) are also outlined. This work is mostly based on previously published data, although we added some new data and hypothesis supported by a number of studies. This review shows the wide variety of allelic effects that provide enormous plasticity in barley flowering behavior, which opens new avenues to breeders for fine-tuning phenology of the barley crop.

Combined effects of a glycine-rich RNA-binding protein and a NAC transcription factor extend grain fill duration and improve malt barley agronomic performance

Two key barley genes independently control anthesis and senescence timing, enabling the manipulation of grain fill duration, grain size/plumpness, and grain protein concentration. Plant developmental processes such as flowering and senescence have direct effects on cereal yield and quality. Previous work highlighted the importance of two tightly linked genes encoding a glycine-rich RNA-binding protein (HvGR-RBP1) and a NAC transcription factor (HvNAM1), controlling barley anthesis timing, senescence, and percent grain protein. Varieties that differ in HvGR-RBP1 expression, 'Karl'(low) and 'Lewis'(high), also differ in sequence 1 KB upstream of translation start site, including an ~ 400 bp G rich insertion in the 5'-flanking region of the 'Karl' allele, which could disrupt gene expression. To improve malt quality, the (low-grain protein, delayed-senescence) 'Karl' HvNAM1 allele was introgressed into Montana germplasm. After several seasons of selection, the resulting germplasm was screened for the allelic combinations of HvGR-RBP1 and HvNAM1, finding lines combining 'Karl' alleles for both genes (-/-), lines combining 'Lewis' (functional, expressed) HvGR-RBP1 with 'Karl' HvNAM1 alleles ( ±), and lines combining 'Lewis' alleles for both genes (+ / +). Field experiments indicate that the functional ('Lewis,'  +) HvGR-RBP1 allele is associated with earlier anthesis and with slightly shorter plants, while the 'Karl' (-) HvNAM1 allele delays maturation. Genotypes carrying the ± allele combination therefore had a significantly (3 days) extended grain fill duration, leading to a higher percentage of plump kernels, slightly enhanced test weight, and lower grain protein concentration when compared to the other allele combinations. Overall, our data suggest an important function for HvGR-RBP1 in the control of barley reproductive development and set the stage for a more detailed functional analysis of this gene.

The influence of breeding history, origin and growth type on population structure of barley as revealed by SSR markers

Natural and mass selection during domestication and cultivation favored particular traits of interest in barley. In the present study, population structure, and genetic relationships among 144 accessions of barley landraces and breeding materials from various countries were studied using a set of 77 and 72 EST-SSR and gSSR markers, respectively distributed on seven chromosomes of barley. In total, 262 and 429 alleles were amplified in 77 EST-SSRs and 72 gSSR loci, respectively. Out of which, 185 private/group-specific alleles were identified in the landraces compared with 14 in "cultivar and advanced breeding lines", indicating the possibility to introgress favorite alleles from landraces into breeding materials. Comparative analysis of genetic variation among breeding materials, Iranian landraces, and exotic landraces revealed higher genetic diversity in Iranian landraces compared with others. A total of 37, 15, and 14 private/group-specific alleles were identified in Iranian landraces, exotic landraces, and breeding materials, respectively. The most likely groups for 144 barley genotypes were three as inferred using model- and distance-based clustering as well as principal coordinate analysis which assigned the landraces and breeding materials into separate groups. The distribution of alleles was found to be correlated with population structure, domestication history and eco-geographical factors. The high allelic richness in the studied set of barley genotype provides insights into the available diversity and allows the construction of core groups based on maximizing allelic diversity for use in barley breeding programs.

Introgression Breeding in Barley: Perspectives and Case Studies

Changing production scenarios resulting from unstable climatic conditions are challenging crop improvement efforts. A deeper and more practical understanding of plant genetic resources is necessary if these assets are to be used effectively in developing improved varieties. In general, current varieties and potential varieties have a narrow genetic base, making them prone to suffer the consequences of new and different abiotic and biotic stresses that can reduce crop yield and quality. The deployment of genomic technologies and sophisticated statistical analysis procedures has generated a dramatic change in the way we characterize and access genetic diversity in crop plants, including barley. Various mapping strategies can be used to identify the genetic variants that lead to target phenotypes and these variants can be assigned coordinates in reference genomes. In this way, new genes and/or new alleles at known loci present in wild ancestors, germplasm accessions, land races, and un-adapted introductions can be located and targeted for introgression. In principle, the introgression process can now be streamlined and linkage drag reduced. In this review, we present an overview of (1) past and current efforts to identify diversity that can be tapped to improve barley yield and quality, and (2) case studies of our efforts to introgress resistance to stripe and stem rust from un-adapted germplasm. We conclude with a description of a modified Nested Association Mapping (NAM) population strategy that we are implementing for the development of multi-use naked barley for organic systems and share perspectives on the use of genome editing in introgression breeding.

A Genome-Wide Association Study of Resistance to Puccinia striiformis f. sp. hordei and P. graminis f. sp. tritici in Barley and Development of Resistant Germplasm

Stripe rust (incited by Puccinia striiformis f. sp. hordei) and stem rust (incited by P. graminis f. sp. tritici) are two of the most important diseases affecting barley. Building on prior work involving the introgression of the resistance genes rpg4/Rpg5 into diverse genetic backgrounds and the discovery of additional quantitative trait locus (QTLs) for stem rust resistance, we generated an array of germplasm in which we mapped resistance to stripe rust and stem rust. Stem rust races TTKSK and QCCJB were used for resistance mapping at the seedling and adult plant stages, respectively. Resistance to stripe rust, at the adult plant stage, was determined by QTLs on chromosomes 1H, 4H, and 5H that were previously reported in the literature. The rpg4/Rpg5 complex was validated as a source of resistance to stem rust at the seedling stage. Some parental germplasm, selected as potentially resistant to stem rust or susceptible but having other positive attributes, showed resistance at the seedling stage, which appears to be allelic to rpg4/Rpg5. The rpg4/Rpg5 complex, and this new allele, were not sufficient for adult plant resistance to stem rust in one environment. A QTL on 5H, distinct from Rpg5 and a previously reported resistance QTL, was required for resistance at the adult plant stage in all environments. This QTL is coincident with the QTL for stripe rust resistance. Germplasm with mapped genes/QTLs conferring resistance to stripe and stem rust was identified and is available as a resource to the research and breeding communities.

Natural Genetic Variation Underlying Tiller Development in Barley (Hordeum vulgare L)

In barley (Hordeum vulgare L.), lateral branches called tillers contribute to grain yield and define shoot architecture, but genetic control of tiller number and developmental rate are not well characterized. The primary objectives of this work were to examine relationships between tiller number and other agronomic and morphological traits and identify natural genetic variation associated with tiller number and rate, and related traits. We grew 768 lines from the USDA National Small Grain Collection in the field and collected data over two years for tiller number and rate, and agronomic and morphological traits. Our results confirmed that spike row-type and days to heading are correlated with tiller number, and as much as 28% of tiller number variance was associated with these traits. In addition, negative correlations between tiller number and leaf width and stem diameter were observed, indicating trade-offs between tiller development and other vegetative growth. Thirty-three quantitative trait loci (QTL) were associated with tiller number or rate. Of these, 40% overlapped QTL associated with days to heading and 22% overlapped QTL associated with spike row-type, further supporting that tiller development is associated with these traits. Some QTL associated with tiller number or rate, including the major QTL on chromosome 3H, were not associated with other traits, suggesting that some QTL may be directly related to rate of tiller development or axillary bud number. These results enhance our knowledge of the genetic control of tiller development in barley, which is important for optimizing tiller number and rate for yield improvement.

Research advances in the Pyrenophora teres-barley interaction

Pyrenophora teres f. teres and P. teres f. maculata are significant pathogens that cause net blotch of barley. An increased number of loci involved in P. teres resistance or susceptibility responses of barley as well as interacting P. teres virulence effector loci have recently been identified through biparental and association mapping studies of both the pathogen and host. Characterization of the resistance/susceptibility loci in the host and the interacting effector loci in the pathogen will provide a path for targeted gene validation for better-informed release of resistant barley cultivars. This review assembles concise consensus maps for all loci published for both the host and pathogen, providing a useful resource for the community to be used in pathogen characterization and barley breeding for resistance to both forms of P. teres.

Evaluation of genetic diversity, agronomic traits, and anthracnose resistance in the NPGS Sudan Sorghum Core collection

Background

The United States Department of Agriculture (USDA) National Plant Germplasm System (NPGS) sorghum core collection contains 3011 accessions randomly selected from 77 countries. Genomic and phenotypic characterization of this core collection is necessary to encourage and facilitate its utilization in breeding programs and to improve conservation efforts. In this study, we examined the genome sequences of 318 accessions belonging to the NPGS Sudan sorghum core set, and characterized their agronomic traits and anthracnose resistance response.

Results

We identified 183,144 single nucleotide polymorphisms (SNPs) located within or in proximity of 25,124 annotated genes using the genotyping-by-sequencing (GBS) approach. The core collection was genetically highly diverse, with an average pairwise genetic distance of 0.76 among accessions. Population structure and cluster analysis revealed five ancestral populations within the Sudan core set, with moderate to high level of genetic differentiation. In total, 171 accessions (54%) were assigned to one of these populations, which covered 96% of the total genomic variation. Genome scan based on Tajima’s D values revealed two populations under balancing selection. Phenotypic analysis showed differences in agronomic traits among the populations, suggesting that these populations belong to different ecogeographical regions. A total of 55 accessions were resistant to anthracnose; these accessions could represent multiple resistance sources. Genome-wide association study based on fixed and random model Circulating Probability (farmCPU) identified genomic regions associated with plant height, flowering time, panicle length and diameter, and anthracnose resistance response. Integrated analysis of the Sudan core set and sorghum association panel indicated that a large portion of the genetic variation in the Sudan core set might be present in breeding programs but remains unexploited within some clusters of accessions.

Conclusions

The NPGS Sudan core collection comprises genetically and phenotypically diverse germplasm with multiple anthracnose resistance sources. Population genomic analysis could be used to improve screening efforts and identify the most valuable germplasm for breeding programs. The new GBS data set generated in this study represents a novel genomic resource for plant breeders interested in mining the genetic diversity of the NPGS sorghum collection.

Non-Invasive Phenotyping Reveals Genomic Regions Involved in Pre-Anthesis Drought Tolerance and Recovery in Spring Barley

With ongoing climate change, drought events are becoming more frequent and will affect biomass formation when occurring during pre-flowering stages. We explored growth over time under such a drought scenario, via non-invasive imaging and revealed the underlying key genetic factors in spring barley. By comparing with well-watered conditions investigated in an earlier study and including information on timing, QTL could be classified as constitutive, drought or recovery-adaptive. Drought-adaptive QTL were found in the vicinity of genes involved in dehydration tolerance such as dehydrins (Dhn4, Dhn7, Dhn8, and Dhn9) and aquaporins (e.g. HvPIP1;5, HvPIP2;7, and HvTIP2;1). The influence of phenology on biomass formation increased under drought. Accordingly, the main QTL during recovery was the region of HvPPD-H1. The most important constitutive QTL for late biomass was located in the vicinity of HvDIM, while the main locus for seedling biomass was the HvWAXY region. The disappearance of QTL marked the genetic architecture of tiller number. The most important constitutive QTL was located on 6HS in the region of 1-FEH. Stage and tolerance specific QTL might provide opportunities for genetic manipulation to stabilize biomass and tiller number under drought conditions and thereby also grain yield.

Environmental Association Identifies Candidates for Tolerance to Low Temperature and Drought

Barley (Hordeum vulgare ssp. vulgare) is cultivated from the equator to the Arctic Circle. The wild progenitor species, Hordeum vulgare ssp. spontaneum, occupies a relatively narrow latitudinal range (∼30 - 40° N) primarily at low elevation (< 1,500 m). Adaptation to the range of cultivation has occurred over ∼8,000 years. The genetic basis of adaptation is amenable to study through environmental association. An advantage of environmental association in a well-characterized crop is that many loci that contribute to climatic adaptation and abiotic stress tolerance have already been identified. This provides the opportunity to determine if environmental association approaches effectively identify these loci of large effect. Using published genotyping from 7,864 SNPs in 803 barley landraces, we examined allele frequency differentiation across multiple partitions of the data and mixed model associations relative to bioclimatic variables. Using newly generated resequencing data from a subset of these landraces, we tested for linkage disequilibrium (LD) between SNPs queried in genotyping and SNPs in neighboring loci. Six loci previously reported to contribute to adaptive differences in flowering time and abiotic stress in barley and six loci previously identified in other plant species were identified in our analyses. In many cases, patterns of LD are consistent with the causative variant occurring in the immediate vicinity of the queried SNP. The identification of barley orthologs to well-characterized genes may provide a new understanding of the nature of adaptive variation and could permit a more targeted use of potentially adaptive variants in barley breeding and germplasm improvement.

Development of a Multiparent Population for Genetic Mapping and Allele Discovery in Six-Row Barley

Germplasm collections hold valuable allelic diversity for crop improvement and genetic mapping of complex traits. To gain access to the genetic diversity within the USDA National Small Grain Collection (NSGC), we developed the Barley Recombinant Inbred Diverse Germplasm Population (BRIDG6), a six-row spring barley multiparent population (MPP) with 88 cultivated accessions crossed to a common parent (Rasmusson). The parents were randomly selected from a core subset of the NSGC that represents the genetic diversity of landrace and breeding accessions. In total, we generated 6160 F5 recombinant inbred lines (RILs), with an average of 69 and a range of 37-168 RILs per family, that were genotyped with 7773 SNPs, with an average of 3889 SNPs segregating per family. We detected 23 quantitative trait loci (QTL) associated with flowering time with five QTL found coincident with previously described flowering time genes. A major QTL was detected near the flowering time gene, HvPpd-H1 which affects photoperiod. Haplotype-based analysis of HvPpd-H1 identified private alleles to families of Asian origin conferring both positive and negative effects, providing the first observation of flowering time-related alleles private to Asian accessions. We evaluated several subsampling strategies to determine the effect of sample size on the power of QTL detection, and found that, for flowering time in barley, a sample size >50 families or 3000 individuals results in the highest power for QTL detection. This MPP will be useful for uncovering large and small effect QTL for traits of interest, and identifying and utilizing valuable alleles from the NSGC for barley improvement.

Exome sequences and multi-environment field trials elucidate the genetic basis of adaptation in barley

Broadening the genetic base of crops is crucial for developing varieties to respond to global agricultural challenges such as climate change. Here, we analysed a diverse panel of 371 domesticated lines of the model crop barley to explore the genetics of crop adaptation. We first collected exome sequence data and phenotypes of key life history traits from contrasting multi-environment common garden trials. Then we applied refined statistical methods, including some based on exomic haplotype states, for genotype-by-environment (G×E) modelling. Sub-populations defined from exomic profiles were coincident with barley's biology, geography and history, and explained a high proportion of trial phenotypic variance. Clear G×E interactions indicated adaptation profiles that varied for landraces and cultivars. Exploration of circadian clock-related genes, associated with the environmentally adaptive days to heading trait (crucial for the crop's spread from the Fertile Crescent), illustrated complexities in G×E effect directions, and the importance of latitudinally based genic context in the expression of large-effect alleles. Our analysis supports a gene-level scientific understanding of crop adaption and leads to practical opportunities for crop improvement, allowing the prioritisation of genomic regions and particular sets of lines for breeding efforts seeking to cope with climate change and other stresses.

Association mapping for agronomic traits in six-rowed spring barley from the USA harvested in Kazakhstan

In barley, six-rowed barley is advantageous over two-rowed barley for feed due to the larger number of seeds per spike and the higher seed protein content. The growth of six-rowed barley is potentially important for breeding in agriculturally oriented countries, such as Kazakhstan. Nevertheless, until recently, very little attention was given to six-rowed barley in breeding projects in Kazakhstan, one of the largest countries in the world. In this study, phenotyping and single nucleotide polymorphism (SNP) genotyping data were generated from 275 accessions originating from six different breeding organizations in the USA as well as 9 accessions from Kazakhstan in field trials at six breeding institutions. The USA six-rowed barley was tested in comparison to local accessions over three years (2009–2011) based on analyses of key agronomic traits. It was determined that the average yield in the USA accessions in comparison to local lines showed heavier yield in all six tested sites. Principal Coordinate Analysis based on 1618 polymorphic SNP markers separated Kazakh lines from six USA barley origin groups based on PC1 (77.9%), and Montana lines from the remaining five USA groups based on PC2 (15.1%). A genome-wide association study based on eighteen field trials allowed the identification of 47 stable marker-trait associations (MTA) for ten agronomic traits, including key yield related characters such as yield per square meter, thousand grain weight, number of kernels per spike, and productive tillers. The comparison of chromosomal positions of identified MTA with positions of known genes and quantitative trait loci suggests that 25 out of those 47 MTAs are presumably novel. The analysis of 42 SNPs associated with 47 MTAs in the Ensemble genome annotation system (http://ensemblgenomes.org) suggested that 40 SNPs were in genic positions of the genome, as their sequences successfully aligned with corresponding Gen ID.

Worldwide Research Trends on Wheat and Barley: A Bibliometric Comparative Analysis

Grain cereals such as wheat, barley, rice, and maize are the nutritional basis of humans and animals worldwide. Thus, these crop plants are essential in terms of global food security. We conducted a bibliometric assessment of scientific documents and patents related to wheat and barley through the Scopus database. The number of documents published per year, their affiliation and corresponding scientific areas, the publishing journals, document types and languages were metricized. The main keywords included in research publications concerning these crops were also analysed globally and clustered in thematic groups. In the case of keywords related to agronomy or genetics and molecular biology, we considered documents dated up to 1999, and from 2000 to 2018, separately. Comparison of the results obtained for wheat and barley revealed some remarkable different trends, for which the underlying reasons are further discussed.

Deciphering the genetic basis for vitamin E accumulation in leaves and grains of different barley accessions

Tocopherols and tocotrienols, commonly referred to as vitamin E, are essential compounds in food and feed. Due to their lipophilic nature they protect biomembranes by preventing the propagation of lipid-peroxidation especially during oxidative stress. Since their synthesis is restricted to photosynthetic organisms, plant-derived products are the major source of natural vitamin E. In the present study the genetic basis for high vitamin E accumulation in leaves and grains of different barley (Hordeum vulgare L.) accessions was uncovered. A genome wide association study (GWAS) allowed the identification of two genes located on chromosome 7H, homogentisate phytyltransferase (HPT-7H) and homogentisate geranylgeranyltransferase (HGGT) that code for key enzymes controlling the accumulation of tocopherols in leaves and tocotrienols in grains, respectively. Transcript profiling showed a correlation between HPT-7H expression and vitamin E content in leaves. Allele sequencing allowed to decipher the allelic variation of HPT-7H and HGGT genes corresponding to high and low vitamin E contents in the respective tissues. Using the obtained sequence information molecular markers have been developed which can be used to assist smart breeding of high vitamin E barley varieties. This will facilitate the selection of genotypes more tolerant to oxidative stress and producing high-quality grains.

Characterization of genes required for both Rpg1 and rpg4-mediated wheat stem rust resistance in barley

Background

Puccinia graminis f. sp. tritici (Pgt) race TTKSK and its lineage pose a threat to barley production world-wide justifying the extensive efforts to identify, clone, and characterize the rpg4-mediated resistance locus (RMRL), the only effective resistance to virulent Pgt races in the TTKSK lineage. The RMRL contains two nucleotide-binding domain and leucine-rich repeat (NLR) resistance genes, Rpg5 and HvRga1, which are required for resistance. The two NLRs have head-to-head genome architecture with one NLR, Rpg5, containing an integrated C-terminal protein kinase domain, characteristic of an “integrated sensory domain” resistance mechanism. Fast neutron mutagenesis of line Q21861 was utilized in a forward genetics approach to identify genetic components that function in the RMRL or Rpg1 resistance mechanisms, as Q21861 contains both genes. A mutant was identified that compromises both RMRL and Rpg1-mediated resistances and had stunted seedling roots, designated required for P. graminis resistance 9 (rpr9).

Results

The rpr9 mutant generated in the Q21861 background was crossed with the Swiss landrace Hv584, which carries RMRL but contains polymorphism across the genome compared to Q21861. To map Rpr9, a Hv584 x rpr9 F_6:7 recombinant inbred line (RIL) population was developed. The RIL population was phenotyped with Pgt race QCCJB. The Hv584 x rpr9 RIL population was genotyped with the 9 k Illumina Infinium iSelect marker panel, producing 2701 polymorphic markers. A robust genetic map consisting of 563 noncosegregating markers was generated and used to map Rpr9 to an ~ 3.4 cM region on barley chromosome 3H. The NimbleGen barley exome capture array was utilized to capture rpr9 and wild type Q21861 exons, followed by Illumina sequencing. Comparative analysis, resulting in the identification of a 1.05 Mbp deletion at the chromosome 3H rpr9 locus. The identified deletion contains ten high confidence annotated genes with the best rpr9 candidates encoding a SKP1-like 9 protein and a F-box family protein.

Conclusion

Genetic mapping and exome capture rapidly identified candidate gene/s that function in RMRL and Rpg1 mediated resistance pathway/s. One or more of the identified candidate rpr9 genes are essential in the only two known effective stem rust resistance mechanisms, present in domesticated barley.

Electronic supplementary material

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

Role of genomics in promoting the utilization of plant genetic resources in genebanks

Global efforts have seen the world's plant genetic resources (PGRs) conserved in about 1625 germ plasm repositories. Utility of these resources is important in increasing the resilience and productivity of agricultural production systems. However, despite their importance, utility of these resources has been poor. This article reviews the real and potential application of the current advances in genomic technologies in improving the utilization of these resources. The actual and potential application of these genomic approaches in plant identification, phylogenetic analysis, analysing the genetic value of germ plasm, facilitating germ plasm selection in genebanks as well as instilling confidence in international germ plasm exchange system is discussed. We note that if genebanks are to benefit from this genomic revolution, there is need for fundamental changes in the way genebanks are managed, perceived, organized and funded. Increased collaboration between genebank managers and the user community is also recommended.

Introgression of rpg4/Rpg5 Into Barley Germplasm Provides Insights Into the Genetics of Resistance to Puccinia graminis f. sp. tritici Race TTKSK and Resources for Developing Resistant Cultivars

Stem rust (incited by Puccinia graminis f. sp. tritici) is a devastating disease of wheat and barley in many production areas. The widely virulent African P. graminis f. sp. tritici race TTKSK is of particular concern, because most cultivars are susceptible. To prepare for the possible arrival of race TTKSK in North America, we crossed a range of barley germplasm-representing different growth habits and end uses-with donors of stem rust resistance genes Rpg1 and rpg4/Rpg5. The former confers resistance to prevalent races of P. graminis f. sp. tritici in North America, and the latter confers resistance to TTKSK and other closely related races from Africa. We produced doubled haploids from these crosses and determined their allele type at the Rpg loci and haplotype at 7,864 single-nucleotide polymorphism loci. The doubled haploids were phenotyped for TTKSK resistance at the seedling stage. Integration of genotype and phenotype data revealed that (i) Rpg1 was not associated with TTKSK resistance, (ii) rpg4/Rpg5 was necessary but was not sufficient for resistance, and (iii) specific haplotypes at two quantitative trait loci were required for rpg4/Rpg5 to confer resistance to TTKSK. To confirm whether lines found resistant to TTKSK at the seedling resistance were also resistant at the adult plant stage, a subset of doubled haploids was evaluated in Kenya. Additionally, adult plant resistance to leaf rust and stripe rust (incited by Puccinia hordei and Puccinia striiformis f. sp. hordei, respectively) was also assessed on the doubled haploids in field trials at three locations in the United States over a 2-year period. Doubled haploids were identified with adult plant resistance to all three rusts, and this germplasm is available to the research and breeding communities.

Harnessing Novel Diversity From Landraces to Improve an Elite Barley Variety

The Spanish Barley Core Collection (SBCC) is a source of genetic variability of potential interest for breeding, particularly for adaptation to Mediterranean environments. Two backcross populations (BC2F5) were developed using the elite cultivar Cierzo as the recurrent parent. The donor parents, namely SBCC042 and SBCC073, were selected from the SBCC lines due to their outstanding yield in drought environments. Flowering time, yield and drought-related traits were evaluated in two field trials in Zaragoza (Spain) during the 2014-15 and 2015-16 seasons and validated in the 2017-18 season. Two hundred sixty-four lines of each population were genotyped with the Barley Illumina iSelect 50k SNP chip. Genetic maps for each population were generated. The map for SBCC042 × Cierzo contains 12,893 SNPs distributed in 9 linkage groups. The map for SBCC073 × Cierzo includes 12,026 SNPs in 7 linkage groups. Both populations shared two QTL hotspots. There are QTLs for flowering time, thousand-kernel weight (TKW), and hectoliter weight on a segment of 23 Mb at ~515 Mb on chromosome 1H, which encompasses the HvFT3 gene. In both populations, flowering was accelerated by the landrace allele, which also increased the TKW. In the same region, better soil coverage was contributed by SBCC042 but coincident with a lower hectoliter weight. The second large hotspot was on chromosome 6H and contained QTLs with wide intervals for grain yield, plant height and TKW. Landrace alleles contributed to increased plant height and TKW and reduced grain yield. Only SBCC042 contributed favorable alleles for "green area," with three significant QTLs that increased ground coverage after winter, which might be exploited as an adaptive trait of this landrace. Some genes of interest found in or very close to the peaks of the QTLs are highlighted. Strategies to deploy the QTLs found for breeding and pre-breeding are proposed.

Genetic association with high-resolution climate data reveals selection footprints in the genomes of barley landraces across the Iberian Peninsula

Landraces are local populations of crop plants adapted to a particular environment. Extant landraces are surviving genetic archives, keeping signatures of the selection processes experienced by them until settling in their current niches. This study intends to establish relationships between genetic diversity of barley (Hordeum vulgare L.) landraces collected in Spain and the climate of their collection sites. A high-resolution climatic data set (5 × 5 km spatial, 1-day temporal grid) was computed from over 2,000 temperature and 7,000 precipitation stations across peninsular Spain. This data set, spanning the period 1981-2010, was used to derive agroclimatic variables meaningful for cereal production at the collection sites of 135 barley landraces. Variables summarize temperature, precipitation, evapotranspiration, potential vernalization and frost probability at different times of the year and time scales (season and month). SNP genotyping of the landraces was carried out combining Illumina Infinium assays and genotyping-by-sequencing, yielding 9,920 biallelic markers (7,479 with position on the barley reference genome). The association of these SNPs with agroclimatic variables was analysed at two levels of genetic diversity, with and without taking into account population structure. The whole data sets and analysis pipelines are documented and available at https://eead-csic-compbio.github.io/barley-agroclimatic-association. We found differential adaptation of the germplasm groups identified to be dominated by reactions to cold temperature and late-season frost occurrence, as well as to water availability. Several significant associations pointing at specific adaptations to agroclimatic features related to temperature and water availability were observed, and candidate genes underlying some of the main regions are proposed.

Development of Decreased-Gluten Wheat Enabled by Determination of the Genetic Basis of lys3a Barley

Celiac disease is the most common food-induced enteropathy in humans, with a prevalence of approximately 1% worldwide. It is induced by digestion-resistant, proline- and glutamine-rich seed storage proteins, collectively referred to as gluten, found in wheat (Triticum aestivum). Related prolamins are present in barley (Hordeum vulgare) and rye (Secale cereale). The incidence of both celiac disease and a related condition called nonceliac gluten sensitivity is increasing. This has prompted efforts to identify methods of lowering gluten in wheat, one of the most important cereal crops. Here, we used bulked segregant RNA sequencing and map-based cloning to identify the genetic lesion underlying a recessive, low-prolamin mutation (lys3a) in diploid barley. We confirmed the mutant identity by complementing the lys3a mutant with a transgenic copy of the wild-type barley gene and then used targeting-induced local lesions in genomes to identify induced single-nucleotide polymorphisms in the three homeologs of the corresponding wheat gene. Combining inactivating mutations in the three subgenomes of hexaploid bread wheat in a single wheat line lowered gliadin and low-molecular-weight glutenin accumulation by 50% to 60% and increased free and protein-bound lysine by 33%.

A tiered approach to genome-wide association analysis for the adherence of hulls to the caryopsis of barley seeds reveals footprints of selection

BACKGROUND

Seeds of domesticated barley are grouped into two distinct types, which differ in morphology. Caryopses covered by adaxial (palea) and abaxial (lemma) hulls that tightly adhere to the pericarp at maturity give rise to hulled seeds whereas caryopses without adhering hulls give rise to naked seeds. The naked caryopsis character is an essential trait regarding the end use of barley.

RESULTS

To uncover the genetic basis of the trait, a genome-wide association study (GWAS) has been performed in a panel comprising 222 2-rowed and 303 6-rowed spring barley landrace accessions. In addition to the well-described Nud locus on chromosome 7H, three novel loci showed strong associations with the trait: the first locus on 2H was specifically detected in 6-rowed accessions, the second locus on 3H was found in 2-rowed accessions from Eurasia and the third locus on 6H was revealed in 6-rowed accessions from Ethiopia. PCR analysis of naked accessions also confirmed the absence of a 17 kb region harboring the Nud gene on chromosome 7H for all but one naked accession. The latter was characterized by a slightly variant phenotype of the caryopsis.

CONCLUSION

Our findings provide evidence of the pervasiveness of the 17 kb deletion in spring barley from different geographic regions and at the same time reveal genomic footprints of selection in naked barley, which follow both geographic and morphological patterns.

Prospects of pan-genomics in barley

The concept of a pan-genome refers to intraspecific diversity in genome content and structure, encompassing both genes and intergenic space. Pan-genomic studies employ a combination of de novo sequence assembly and reference-based alignment to discover and genotype structural variants. The large size and complex structure of Triticeae genomes were for a long time an obstacle for genomic research in barley and its relatives. Now that a reference genome is available, computational pipelines for high-quality sequence assembly are in place, and sequence costs continue to drop, investigations into the structural diversity of the barley genome seem within reach. Here, we review the recent progress on pan-genomics in the model grass Brachypodium distachyon, and the cereal crops rice and maize, and devise a multi-tiered strategy for a pan-genome project in barley. Our design involves: (1) the construction of high-quality de novo sequence assemblies for a small core set of representative genotypes, (2) short-read sequencing of a large diversity panel of genebank accessions to medium coverage and (3) the use of complementary methods such as chromosome-conformation capture sequencing and k-mer-based association genetics. The in silico representation of the barley pan-genome may inform about the mechanisms of structural genome evolution in the Triticeae and supplement quantitative genetics models of crop performance for better accuracy and predictive ability.

Pyramiding rpg4- and Rpg1-Mediated Stem Rust Resistance in Barley Requires the Rrr1 Gene for Both to Function

Stem rust, caused by Puccinia graminis f. sp. tritici (Pgt) is an economically important disease of wheat and barley. Rpg1 is the only resistance gene deployed in Midwestern US barley varieties and provides remarkable resistance to most North American races, except Pgt race QCCJB. Rpg1 is also ineffective against Pgt race TTKSK and its lineage that originated in Africa. The barley rpg4-mediated resistance locus (RMRL) conferring resistance to Pgt races QCCJB and TTKSK was isolated from line Q21861, which is resistant to all known Pgt races due to Rpg1 and RMRL. To develop elite barley varieties RMRL was pyramided into the varieties, Pinnacle and Conlon (both contain Rpg1), producing the near isogenic lines (NILs), Pinnacle RMRL-NIL (PRN) and Conlon RMRL-NIL (CRN). The CRN was resistant to Pgt races QCCJB (RMRL specific) and HKHJC (Rpg1 specific) at the seedling stage and Pgt race TTKSK (RMRL specific) at the adult stage. In contrast, PRN was susceptible to QCCJB and HKHJC at the seedling stage and TTKSK at the adult stage. Interestingly, PRN's susceptibility to QCCJB and HKHJC showed that RMRL was non-functional in the Pinnacle background but its presence also suppressed Rpg1-mediated resistance. Thus, in the absence of a gene/s found in the Q21861 background, Rpg1 becomes non-functional if RMRL is present, suggesting that another polymorphic gene, that we designated Rrr1 (required for rpg4-mediated resistance 1), is required for RMRL resistance and Rpg1-mediated resistance in the presence of RMRL. Utilizing a PRN/Q21861 derived recombinant inbred line (RIL) population, Rrr1 was delimited to a ∼0.5 MB physical region, slightly proximal (∼1.8 MB) of RMRL on barley chromosome 5H. A second gene, designated required for Rpg1-mediated resistance 2 (Rrr2), with duplicate gene action to Rrr1 in Rpg1-mediated resistance function, was genetically delimited to a physical region of ∼0.7 MB, slightly distal (∼3.1 MB) to Rpg1 on the short arm of barley chromosome 7H. Thus, Rrr1 is required for RMRL resistance and Rrr1 or Rrr2 is required for functional Rpg1-mediated resistance in the presence of the RMRL introgression. Candidate Rrr1 and Rrr2 genes were identified that need to be considered when pyramiding Rpg1 and RMRL in barley.

Unbalanced historical phenotypic data from seed regeneration of a barley ex situ collection

The scarce knowledge on phenotypic characterization restricts the usage of genetic diversity of plant genetic resources in research and breeding. We describe original and ready-to-use processed data for approximately 60% of ~22,000 barley accessions hosted at the Federal ex situ Genebank for Agricultural and Horticultural Plant Species. The dataset gathers records for three traits with agronomic relevance: flowering time, plant height and thousand grain weight. This information was collected for seven decades for winter and spring barley during the seed regeneration routine. The curated data represent a source for research on genetics and genomics of adaptive and yield related traits in cereals due to the importance of barley as model organism. This data could be used to predict the performance of non-phenotyped individuals in other collections through genomic prediction. Moreover, the dataset empowers the utilization of phenotypic diversity of genetic resources for crop improvement.