Analysis of molecular diversity, population structure and linkage disequilibrium in a worldwide survey of cultivated barley germplasm (Hordeum vulgare L.)

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.

Genetic diversity and population structure assessment of Western Canadian barley cooperative trials

Studying the population structure and genetic diversity of historical datasets is a proposed use for association analysis. This is particularly important when the dataset contains traits that are time-consuming or costly to measure. A set of 96 elite barley genotypes, developed from eight breeding programs of the Western Canadian Cooperative Trials were used in the current study. Genetic diversity, allelic variation, and linkage disequilibrium (LD) were investigated using 5063 high-quality SNP markers via the Illumina 9K Barley Infinium iSelect SNP assay. The distribution of SNPs markers across the barley genome ranged from 449 markers on chromosome 1H to 1111 markers on chromosome 5H. The average polymorphism information content (PIC) per locus was 0.275 and ranged from 0.094 to 0.375. Bayesian clustering in STRUCTURE and principal coordinate analysis revealed that the populations are differentiated primarily due to the different breeding program origins and ear-row type into five subpopulations. Analysis of molecular variance based on PhiPT values suggested that high values of genetic diversity were observed within populations and accounted for 90% of the total variance. Subpopulation 5 exhibited the most diversity with the highest values of the diversity indices, which represent the breeding program gene pool of AFC, AAFRD, AU, and BARI. With increasing genetic distance, the LD values, expressed as r², declined to below the critical r^{2 =} 0.18 after 3.91 cM, and the same pattern was observed on each chromosome. Our results identified an important pattern of genetic diversity among the Canadian barley panel that was proposed to be representative of target breeding programs and may have important implications for association mapping in the future. This highlight, that efforts to identify novel variability underlying this diversity may present practical breeding opportunities to develop new barley genotypes.

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.

The Population Structure of a Globe Artichoke Worldwide Collection, as Revealed by Molecular and Phenotypic Analyzes

The knowledge of the organization of the domesticated gene pool of crop species is an essential requirement to understand crop evolution, to rationalize conservation programs, and to support practical decisions in plant breeding. Here, we integrate simple sequence repeat (SSR) analysis and phenotypic characterization to investigate a globe artichoke collection that comprises most of the varieties cultivated worldwide. We show that the cultivated gene pool of globe artichoke includes five distinct genetic groups associated with the major phenotypic typologies: Catanesi (which based on our analysis corresponds to Violetti di Provenza), Spinosi, Violetti di Toscana, Romaneschi, and Macau. We observed that 17 and 11% of the molecular and phenotypic variance, respectively, is between these groups, while within groups, strong linkage disequilibrium and heterozygote excess are evident. The divergence between groups for quantitative traits correlates with the average broad-sense heritability within the groups. The phenotypic divergence between groups for both qualitative and quantitative traits is strongly and positively correlated with SSR divergence (F_ST) between groups. All this implies a low population size and strong bottleneck effects, and indicates a long history of clonal propagation and selection during the evolution of the domesticated gene pool of globe artichoke. Moreover, the comparison between molecular and phenotypic population structures suggests that harvest time, plant architecture (i.e., plant height, stem length), leaf spininess, head morphology (i.e., head shape, bract shape, spininess) together with the number of heads per plant were the main targets of selection during the evolution of the cultivated germplasm. We emphasize our findings in light of the potential exploitation of this collection for association mapping studies.

Modeling allelic diversity of multiparent mapping populations affects detection of quantitative trait loci

The search for quantitative trait loci that explain complex traits such as yield and drought tolerance has been ongoing in all crops. Methods such as biparental quantitative trait loci mapping and genome-wide association studies each have their own advantages and limitations. Multiparent advanced generation intercross populations contain more recombination events and genetic diversity than biparental mapping populations and are better able to estimate effect sizes of rare alleles than association mapping populations. Here, we discuss the results of using a multiparent advanced generation intercross population of doubled haploid maize lines created from 16 diverse founders to perform quantitative trait loci mapping. We compare 3 models that assume bi-allelic, founder, and ancestral haplotype allelic states for quantitative trait loci. The 3 methods have differing power to detect quantitative trait loci for a variety of agronomic traits. Although the founder approach finds the most quantitative trait loci, all methods are able to find unique quantitative trait loci, suggesting that each model has advantages for traits with different genetic architectures. A closer look at a well-characterized flowering time quantitative trait loci, qDTA8, which contains vgt1, highlights the strengths and weaknesses of each method and suggests a potential epistatic interaction. Overall, our results reinforce the importance of considering different approaches to analyzing genotypic datasets, and shows the limitations of binary SNP data for identifying multiallelic quantitative trait loci.

A Genome-Wide Association Study Pinpoints Quantitative Trait Genes for Plant Height, Heading Date, Grain Quality, and Yield in Rye (Secale cereale L.)

Rye is the only cross-pollinating Triticeae crop species. Knowledge of rye genes controlling complex-inherited traits is scarce, which, currently, largely disables the genomics assisted introgression of untapped genetic variation from self-incompatible germplasm collections in elite inbred lines for hybrid breeding. We report on the first genome-wide association study (GWAS) in rye based on the phenotypic evaluation of 526 experimental hybrids for plant height, heading date, grain quality, and yield in 2 years and up to 19 environments. We established a cross-validated NIRS calibration model as a fast, effective, and robust analytical method to determine grain quality parameters. We observed phenotypic plasticity in plant height and tiller number as a resource use strategy of rye under drought and identified increased grain arabinoxylan content as a striking phenotype in osmotically stressed rye. We used DArTseq™ as a genotyping-by-sequencing technology to reduce the complexity of the rye genome. We established a novel high-density genetic linkage map that describes the position of almost 19k markers and that allowed us to estimate a low genome-wide LD based on the assessed genetic diversity in elite germplasm. We analyzed the relationship between plant height, heading date, agronomic, as well as grain quality traits, and genotype based on 20k novel single-nucleotide polymorphism markers. In addition, we integrated the DArTseq™ markers in the recently established 'Lo7' reference genome assembly. We identified cross-validated SNPs in 'Lo7' protein-coding genes associated with all traits studied. These include associations of the WUSCHEL-related homeobox transcription factor DWT1 and grain yield, the DELLA protein gene SLR1 and heading date, the Ethylene overproducer 1-like protein gene ETOL1 and thousand-grain weight, protein and starch content, as well as the Lectin receptor kinase SIT2 and plant height. A Leucine-rich repeat receptor protein kinase and a Xyloglucan alpha-1,6-xylosyltransferase count among the cross-validated genes associated with water-extractable arabinoxylan content. This study demonstrates the power of GWAS, hybrid breeding, and the reference genome sequence in rye genetics research to dissect and identify the function of genes shaping genetic diversity in agronomic and grain quality traits of rye. The described links between genetic causes and phenotypic variation will accelerate genomics-enabled rye improvement.

Characterization of the complete chloroplast genome of Hordeum vulgare L. var. trifurcatum with phylogenetic analysis

In the present study, the complete chloroplast genome of Hordeum vulgare L. var. trifurcatum was sequenced, assembled and compared with closely related species. The chloroplast genome of Hordeum vulgare L. var. trifurcatum was composed of 84 protein-coding genes (PCG), 8 ribosomal RNA (rRNA) genes, and 38 transfer RNA (tRNA) genes. The Hordeum vulgare L. var. trifurcatum chloroplast genome is 136,485 bp in size, with the GC content of 38.32%. Phylogenetic analysis based on the combined chloroplast gene dataset indicated that the Hordeum vulgare L. var. trifurcatum exhibited a close relationship with Hordeum vulgare subsp. spontaneum and Hordeum vulgare subsp. vulgare.

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.

Genomics Analyses Reveal Unique Classification, Population Structure and Novel Allele of Neo-Tetraploid Rice

BACKGROUND

Neo-tetraploid rice (NTR) is a useful new germplasm that developed from the descendants of the autotetraploid rice (ATR) hybrids. NTR showed improved fertility and yield potential, and produced high yield heterosis when crossed with indica ATR for commercial utilization. However, their classification, population structure and genomic feature remain elusive.

RESULTS

Here, high-depth genome resequencing data of 15 NTRs and 18 ATRs, together with 38 publicly available data of diploid rice accessions, were analyzed to conduct classification, population structure and haplotype analyses. Five subpopulations were detected and NTRs were clustered into one independent group that was adjacent to japonica subspecies, which maybe the reason for high heterosis when NTRs crossed with indica ATRs. Haplotype patterns of 717 key genes that associated with yield and other agronomic traits were revealed in these NTRs. Moreover, a novel specific SNP variation was detected in the first exon of HSP101, a known heat-inducible gene, which was conserved in all NTRs but absent in ATRs, 3KRG and RiceVarMap2 databases. The novel allele was named as HSP101-1, which was confirmed to be a heat response factor by qRT-PCR, and knockout of HSP101-1 significantly decreased the thermotolerance capacity of NTR. Interestingly, HSP101-1 was also specifically expressed in the anthers of NTR at pre-meiotic and meiosis stages under optimal environment without heat stress, and its loss-of-function mutant showed significant decrease in fertility of NTR.

CONCLUSION

The construction of first genomic variation repository and the revelation of population structure provide invaluable information for optimizing the designs of tetraploid rice breeding. The detection of specific genomic variations offered useful genomic markers and new directions to resolve high fertility mechanism of NTR.

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.

Hybridisation-based target enrichment of phenology genes to dissect the genetic basis of yield and adaptation in barley

Barley (Hordeum vulgare L.) is a major cereal grain widely used for livestock feed, brewing malts and human food. Grain yield is the most important breeding target for genetic improvement and largely depends on optimal timing of flowering. Little is known about the allelic diversity of genes that underlie flowering time in domesticated barley, the genetic changes that have occurred during breeding, and their impact on yield and adaptation. Here, we report a comprehensive genomic assessment of a worldwide collection of 895 barley accessions based on the targeted resequencing of phenology genes. A versatile target-capture method was used to detect genome-wide polymorphisms in a panel of 174 flowering time-related genes, chosen based on prior knowledge from barley, rice and Arabidopsis thaliana. Association studies identified novel polymorphisms that accounted for observed phenotypic variation in phenology and grain yield, and explained improvements in adaptation as a result of historical breeding of Australian barley cultivars. We found that 50% of genetic variants associated with grain yield, and 67% of the plant height variation was also associated with phenology. The precise identification of favourable alleles provides a genomic basis to improve barley yield traits and to enhance adaptation for specific production areas.

Genetic variation of HvXYN1 associated with endoxylanase activity and TAX content in barley (Hordeum vulgare L.)

BACKGROUND

Endo-β-1,4-xylanase1 (EA), the key endoxylanase in plants, is involved in the degradation of arabinoxylan during grain germination. In barley (Hordeum vulgare L.), one gene (HvXYN-1) that encode a endo-beta-1,4-xylanase, has been cloned. However, the single nucleotide polymorphisms (SNPs) that affect the endoxylanase activity and total arabinoxylan (TAX) content have yet to be characterized. The investigation of genetic variation in HvXYN1 may facilitate a better understanding of the relationship between TAX content and EA activity in barley.

RESULTS

In the current study, 56 polymorphisms were detected in HvXYN1 among 210 barley accessions collected from 34 countries, with 10 distinct haplotypes identified. The SNPs at positions 110, 305, 1045, 1417, 1504, 1597, 1880 bp in the genomic region of HvXYN1 were significantly associated with EA activity (P < 0.0001), and the sites 110, 305, and 1045 were highly significantly associated with TAX content. The amount of phenotypic variation in a given trait explained by each associated polymorphism ranged from 6.96 to 9.85%. Most notably, we found two variants at positions 1504 bp and 1880 bp in the second exon that significantly (P < 0.0001) affected EA activity; this result could be used in breeding programs to improve beer quality. In addition, African accessions had the highest EA activity and TAX content, and the richest germplasm resources were from Asia, indicating the high potential value of Asian barley.

CONCLUSION

This study provided insight into understanding the relationship, EA activity, TAX content with the SNPs of HvXYN1 in barley. These SNPs can be applied as DNA markers in breeding programs to improve the quality of barley for beer brewing after further validation.

Identification of Ideal Allele Combinations for the Adaptation of Spring Barley to Northern Latitudes

The northwards expansion of barley production requires adaptation to longer days, lower temperatures and stronger winds during the growing season. We have screened 169 lines of the current barley breeding gene pool in the Nordic region with regards to heading, maturity, height, and lodging under different environmental conditions in nineteen field trials over 3 years at eight locations in northern and central Europe. Through a genome-wide association scan we have linked phenotypic differences observed in multi-environment field trials (MET) to single nucleotide polymorphisms (SNP). We have identified an allele combination, only occurring among a few Icelandic lines, that affects heat sum to maturity and requires 214 growing degree days (GDD) less heat sum to maturity than the most common allele combination in the Nordic spring barley gene pool. This allele combination is beneficial in a cold environment, where autumn frost can destroy a late maturing harvest. Despite decades of intense breeding efforts relying heavily on the same germplasm, our results show that there still exists considerable variation within the current breeding gene pool and we identify ideal allele combinations for regional adaptation, which can facilitate the expansion of cereal cultivation even further northwards.

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.

Genetic Variation and Alleviation of Salinity Stress in Barley (Hordeum vulgare L.)

Barley (Hordeum vulgare L.) represents one of the most important cereals cultivated worldwide. Investigating genetic variability and structure of barley is important for enhancing the crop productivity. This study aimed to investigate the diversity and structure of 40 barley genotypes originated from three European countries (France, the Netherlands, Poland) using amplified fragment length polymorphisms (AFLPs). It also aimed to study 5-aminolevulinic acid (ALA) effect on salinity tolerance of six barley genotypes. The expected heterozygosity (H_e) diverged from 0.126 to 0.501, with a mean of 0.348. Polymorphic information content (PIC) diverged from 0.103 to 0.482 across barley genotypes, with a mean of 0.316, indicating that barley genotypes are rich in a considerable level of genetic diversity. The 40 barley genotypes were further studied based on their geographical origin (Western Europe and Eastern Europe). The Eastern European region (Poland) has a higher barley variability than the Western European region (France and the Netherlands). Nei’s distance-based cluster tree divided the 40 barley accessions into two major clusters; one cluster comprised all the varieties originated from the Eastern European region, while the other major cluster included all accessions originated from the Western European region. Structure analysis results were in a complete concordance with our cluster analysis results. Slaski 2, Damseaux and Urbanowicki genotypes have the highest diversity level, whereas Carmen, Bigo and Cambrinus genotypes have the lowest level. The response of these six varieties to NaCl stress was also investigated. Salt stress (100 mM NaCl) slightly decreased levels of chlorophyll, carotenoid and osmolytes (proteins, soluble sugars, phenolics and flavonoids) in the leaves of Slaski 2, Damseaux and Urbanowicki genotypes at non-significant level, as compared to control samples. However, pigment contents and osmolytes in leaves of Carmen, Bigo and Cambrinus genotypes were significantly decreased by salt stress. Antioxidant enzyme activities were significantly increased in Slaski 2 genotype, but non-significantly increased in Carmen by salt stress. Priming Slaski 2 and Carmen cultivars with ALA under salt stress significantly induced pigment contents, antioxidants enzymes activity and stress-responsive genes expression, relative to NaCl-stressed plants. In conclusion, this study suggested a correlation between variability percentage and degree of salinity resistance. ALA improved salt tolerance in barley.

GWAS analysis in spring barley (Hordeum vulgare L.) for morphological traits exposed to drought

Association analysis based on linkage disequilibrium has become a common and powerful approach for detection of QTLs underlying complex agronomic traits including drought tolerance. To determine marker/trait association, 148 modern European spring barley cultivars were evaluated under drought stress. Associations of morphological traits with AFLP/SSR markers were investigated based on the mixed linear model using the TASSEL_3.0. Population structure was estimated using various methods including Bayesian clustering model by STRUCTURE software, PCoA analysis, NJ dendrogram and Hierarchical Clustering. Linkage disequilibrium patterns were explored among the whole genome and each chromosome separately. All the analysis for population structure divided the population into two sub-groups. Linkage disequilibrium analysis showed that by increasing genetic distance, LD decreases. Totally, 167 significant marker trait associations were found which delineated into 65 QTLs in both treatments. Two stable QTLs on 5H at 86.880 cM were detected for Internode Length and on 3H at 126.421 cM for flag leaf length in drought stress treatment. Fourteen QTLs were co-localized with previously reported QTLs and others were novel. The results indicate that these putative genomic regions contain genes that have pleiotropic effects on morphological traits in drought condition.

Unlocking historical phenotypic data from an ex situ collection to enhance the informed utilization of genetic resources of barley (Hordeum sp.)

Key message Historical data generated during seed regeneration are valuable to populate a bio-digital resource center for barley (Hordeum sp.). Precise estimates of trait performance of genetic resources are considered as an intellectually challenging, complex, costly and time-consuming step needed to exploit the phenotypic and genetic diversity maintained in genebanks for breeding and research. Using barley (Hordeum sp.) as a model, we examine strategies to tap into historical data available from regeneration trials. This is a first step toward extending the Federal ex situ Genebank into a bio-digital resource center facilitating an informed choice of barley accessions for research and breeding. Our study is based on historical data of seven decades collected for flowering time, plant height, and thousand grain weight during the regeneration of 12,872 spring and winter barley accessions. Linear mixed models were implemented in conjunction with routines for assessment of data quality. A resampling study highlights the potential risk of biased estimates in second-order statistics when grouping accessions for regeneration according to the year of collection or geographic origin. Based on rigorous quality assessment, we obtained high heritability estimates for the traits under consideration exceeding 0.8. Thus, the best linear unbiased estimations for the three traits are a valuable source to populate a bio-digital resource center for the IPK barley collection. The proposed strategy to leverage historical data from regeneration trials is not crop specific and can be used as a blueprint for other ex situ collections.

Association of functional markers with flowering time in lentil

In the present study, a diverse panel of 96 accessions of lentil germplasm was used to study flowering time over environments and to identify simple sequence repeat markers associated with flowering time through association mapping. The study showed high broad sense heritability estimate (h ² bs=0.93) for flowering time in lentil. Screening of 534 SSR markers resulted in an identification of 75 SSR polymorphic markers (13.9%) across studied genotypes. These markers amplified 266 loci and generated 697 alleles ranging from two to 16 alleles per locus. Model-based cluster analysis used for the determination of population structure resulted in the identification of two distinct subpopulations. Distribution of flowering time was ranged from 40 to 70 days in subpopulation I and from 54 to 69 days in subpopulation II and did not skew either late or early flowering time within a subpopulation. No admixture was observed within the subpopulations. Use of the most accepted maximum likelihood model (P3D mixed linear model with optimum compression) of MTA analysis showed significant association of 26 SSR markers with flowering time at <0.05 probability. The percent of phenotypic explained by each associated marker with flowering time ranged from 2.1 to 21.8% and identified QTLs for flowering time explaining high phenotypic variation across the environments (10.7-21.8%) or in a particular environment (10.2-21.4%). In the present study, 13 EST-SSR showed significant association with flowering time and explained large phenotypic variation (2.3-21.8%) compared to genomic SSR markers (2.1-10.2%). Hence, these markers can be used as functional markers in the lentil breeding program to develop short duration cultivars.

Genetic diversity, linkage disequilibrium, and association mapping analyses of Gossypium barbadense L. germplasm

Limited polymorphism and narrow genetic base, due to genetic bottleneck through historic domestication, highlight a need for comprehensive characterization and utilization of existing genetic diversity in cotton germplasm collections. In this study, 288 worldwide Gossypium barbadense L. cotton germplasm accessions were evaluated in two diverse environments (Uzbekistan and USA). These accessions were assessed for genetic diversity, population structure, linkage disequilibrium (LD), and LD-based association mapping (AM) of fiber quality traits using 108 genome-wide simple sequence repeat (SSR) markers. Analyses revealed structured population characteristics and a high level of intra-variability (67.2%) and moderate interpopulation differentiation (32.8%). Eight percent and 4.3% of markers revealed LD in the genome of the G. barbadense at critical values of r2 ≥ 0.1 and r2 ≥ 0.2, respectively. The LD decay was on average 24.8 cM at the threshold of r2 ≥ 0.05. LD retained on average distance of 3.36 cM at the threshold of r2 ≥ 0.1. Based on the phenotypic evaluations in the two diverse environments, 100 marker loci revealed a strong association with major fiber quality traits using mixed linear model (MLM) based association mapping approach. Fourteen marker loci were found to be consistent with previously identified quantitative trait loci (QTLs), and 86 were found to be new unreported marker loci. Our results provide insights into the breeding history and genetic relationship of G. barbadense germplasm and should be helpful for the improvement of cotton cultivars using molecular breeding and omics-based technologies.

A Novel QTL for Powdery Mildew Resistance in Nordic Spring Barley (Hordeum vulgare L. ssp. vulgare) Revealed by Genome-Wide Association Study

The powdery mildew fungus, Blumeria graminis f. sp. hordei is a worldwide threat to barley (Hordeum vulgare L. ssp. vulgare) production. One way to control the disease is by the development and deployment of resistant cultivars. A genome-wide association study was performed in a Nordic spring barley panel consisting of 169 genotypes, to identify marker-trait associations significant for powdery mildew. Powdery mildew was scored during three years (2012-2014) in four different locations within the Nordic region. There were strong correlations between data from all locations and years. In total four QTLs were identified, one located on chromosome 4H in the same region as the previously identified mlo locus and three on chromosome 6H. Out of these three QTLs identified on chromosome 6H, two are in the same region as previously reported QTLs for powdery mildew resistance, whereas one QTL appears to be novel. The top NCBI BLASTn hit of the SNP markers within the novel QTL predicted the responsible gene to be the 26S proteasome regulatory subunit, RPN1, which is required for innate immunity and powdery mildew-induced cell death in Arabidopsis. The results from this study have revealed SNP marker candidates that can be exploited for use in marker-assisted selection and stacking of genes for powdery mildew resistance in barley.

Genetic variations of HvP5CS1 and their association with drought tolerance related traits in barley (Hordeum vulgare L.)

Delta-1-pyrroline-5-carboxylate synthase gene1 (P5CS1) is the key gene involved in the biosynthesis of proline and is significantly induced by drought stress. The exploration of genetic variation in HvP5CS1 may facilitate a better understanding of the mechanism of drought adaptation in barley. In the current study, 41 polymorphisms including 16 single nucleotide polymorphisms (SNPs) and 25 insertions/deletions (indels) were detected in HvP5CS1 among 287 barley (Hordeum vulgare L.) accessions collected worldwide, with 13 distinct haplotypes identified in the barley collection. Five polymorphisms in HvP5CS1 were significantly (P < 0.001) associated with drought tolerance related traits in barley. The phenotypic variation of a given trait explained by each associated polymorphism ranged from 4.43% to 9.81%. Two sequence variations that were significantly (P < 0.0001) associated with grain yield had marginally significant positive Tajima’s D values in the sliding window, so they might have been selected for environmental adaptation. Meanwhile, two haplotypes HvP5CS1_H1 and HvP5CS1_H4, which contained desired alleles of the two variations mentioned above, were significantly (P < 0.001) associated with drought tolerance related traits, and explained 5.00~11.89% of the phenotypic variations. These variations associated with drought tolerance related traits can be used as potential markers for improving drought tolerance in barley.

Association mapping analysis of fiber yield and quality traits in Upland cotton (Gossypium hirsutum L.)

Fiber yield and quality are the most important traits for Upland cotton (Gossypium hirsutum L.). Identifying high yield and good fiber quality genes are the prime concern of researchers in cotton breeding. Association mapping offers an alternative and powerful method for detecting those complex agronomic traits. In this study, 198 simple sequence repeats (SSRs) were used to screen markers associated with fiber yield and quality traits with 302 elite Upland cotton accessions that were evaluated in 12 locations representing the Yellow River and Yangtze River cotton growing regions of China. Three subpopulations were found after the estimation of population structure. The pair-wise kinship values varied from 0 to 0.867. Only 1.59% of the total marker locus pairs showed significant linkage disequilibrium (LD, p < 0.001). The genome-wide LD decayed within the genetic distance of ~30 to 32 cM at r ² = 0.1, and decreased to ~1 to 2 cM at r ² = 0.2, indicating the potential for association mapping. Analysis based on a mixed linear model detected 57 significant (p < 0.01) marker-trait associations, including seven associations for fiber length, ten for fiber micronaire, nine for fiber strength, eight for fiber elongation, five for fiber uniformity index, five for fiber uniformity ratio, six for boll weight and seven for lint percent, for a total of 35 SSR markers, of which 11 markers were associated with more than one trait. Among marker-trait associations, 24 associations coincided with the previously reported quantitative trait loci (QTLs), the remainder were newly identified QTLs/genes. The QTLs identified in this study will potentially facilitate improvement of fiber yield and quality in the future cotton molecular breeding programs.

Comparative assessment of genetic diversity in Sesamum indicum L. using RAPD and SSR markers

Sesame (Sesamum indicum L.) is an ancient oilseed crop known for its nutty seeds and high-quality edible oil. It is an unexplored crop with a great economic potential. The present study deals with assessment of genetic diversity in the crop. Twenty two RAPD and 18 SSR primers were used for analysis of the 47 different sesame accessions grown in different agroclimatic zones of India. A total of 256 bands were obtained with RAPD primers, of which 191 were polymorphic. SSR primers gave 64 DNA bands, of which all of were polymorphic. The Jaccardʼs similarity coefficient of RAPD, SSR, and pooled RAPD and SSR data ranged from 0.510 to 0.885, 0.167 to 0.867, and 0.505 to 0.853, respectively. Maximum polymorphic information content was reported with SSRs (0.194) compared to RAPDs (0.186). Higher marker index was observed with RAPDs (1.426) than with SSRs (0.621). Similarly, maximum resolving power was found with RAPD (4.012) primers than with SSRs (0.884). The RAPD primer RPI-B11 and SSR primer S16 were the most informative in terms of describing genetic variability among the varieties under study. At a molecular level, the seed coat colour was distinguishable by the presence and absence of a group of marker amplicon/s. White and brown seeded varieties clustered close to each other, while black seeded varieties remained distanced from the cluster. In the present study, we found higher variability in Sesamum indicum L. using RAPD and SSR markers and these could assist in DNA finger printing, conservation of germplasm, and crop improvement.

Patterns of Evolutionary Trajectories and Domestication History within the Genus Hordeum Assessed by REMAP Markers

The patterns of genetic diversity related to the taxonomy and domestication history of 85 accessions representing the main four species of the genus Hordeum were examined by retrotransposon-microsatellite amplified polymorphism (REMAP) markers based on the retrotransposon BARE-1. A substantial level of genetic polymorphisms at among- and within-species level was observed showing that this retrotransposon family and its adjacent genomic regions has been a target for genome dynamics during the evolution and domestication of barley. The obtained data are consistent with the current taxonomic status within the genus Hordeum. Similar level of genetic diversity was observed between the wild and the domesticated barley accessions suggesting that transposable elements` activity and accumulation may counteract the decrease of genome-wide diversity following domestication. In addition, eco-geographical sub-genome pools of the cultivated barley were identified in support to the theory of multiple origins of domestication within the genus Hordeum. We also provide conclusions about the relationship between accessions of different species and the putative routes of barley domestication. In conclusion, the retrotransposon BARE-1 stands as a reliable and perspective DNA marker for the assessment of the phylogenetic and domestication history in the genus Hordeum and other crop species.

Genome-Wide Association Study of Spot Form of Net Blotch Resistance in the Upper Midwest Barley Breeding Programs

Pyrenophora teres f. maculata, the causal agent of spot form of net blotch (SFNB), is an emerging pathogen of barley in the United States and Australia. Compared with net form of net blotch (NFNB), less is known in the U.S. Upper Midwest barley breeding programs about host resistance and quantitative trait loci (QTL) associated with SFNB in breeding lines. The main objective of this study was to identify QTL associated with SFNB resistance in the Upper Midwest two-rowed and six-rowed barley breeding programs using a genome-wide association study approach. A total of 376 breeding lines of barley were evaluated for SFNB resistance at the seedling stage in the greenhouse in Fargo in 2009. The lines were genotyped with 3,072 single nucleotide polymorphism (SNP) markers. Phenotypic evaluation showed a wide range of variability among populations from the four breeding programs and the two barley-row types. The two-rowed barley lines were more susceptible to SFNB than the six-rowed lines. Continuous distributions of SFNB severity indicate the quantitative nature of SFNB resistance. The mixed linear model (MLM) analysis, which included both population structure and kinship matrices, was used to identify significant SNP-SFNB associations. Principal component analysis was used to control false marker-trait association. The linkage disequilibrium (LD) estimates varied among chromosomes (10 to 20 cM). The MLM analysis identified 10 potential QTL in barley: SFNB-2H-8-10, SFNB-2H-38.03, SFNB-3H-58.64, SFNB-3H-78.53, SFNB-3H-91.88, SFNB-3H-117.1, SFNB-5H-155.3, SFNB-6H-5.4, SFNB-6H-33.74, and SFNB-7H-34.82. Among them, four QTL (SFNB-2H-8-10, SFNB-2H-38.03 SFNB-3H-78.53, and SFNB-3H-117.1) have not previously been published. Identification of SFNB resistant lines and QTL associated with SFNB resistance in this study will be useful in the development of barley genotypes with better SFNB resistance.

Genomic Prediction of Seed Quality Traits Using Advanced Barley Breeding Lines

Genomic selection was recently introduced in plant breeding. The objective of this study was to develop genomic prediction for important seed quality parameters in spring barley. The aim was to predict breeding values without expensive phenotyping of large sets of lines. A total number of 309 advanced spring barley lines tested at two locations each with three replicates were phenotyped and each line was genotyped by Illumina iSelect 9Kbarley chip. The population originated from two different breeding sets, which were phenotyped in two different years. Phenotypic measurements considered were: seed size, protein content, protein yield, test weight and ergosterol content. A leave-one-out cross-validation strategy revealed high prediction accuracies ranging between 0.40 and 0.83. Prediction across breeding sets resulted in reduced accuracies compared to the leave-one-out strategy. Furthermore, predicting across full and half-sib-families resulted in reduced prediction accuracies. Additionally, predictions were performed using reduced marker sets and reduced training population sets. In conclusion, using less than 200 lines in the training set can result in low prediction accuracy, and the accuracy will then be highly dependent on the family structure of the selected training set. However, the results also indicate that relatively small training sets (200 lines) are sufficient for genomic prediction in commercial barley breeding. In addition, our results indicate a minimum marker set of 1,000 to decrease the risk of low prediction accuracy for some traits or some families.

Assessment of genetic diversity in Brazilian barley using SSR markers

Barley is a major cereal grown widely and used in several food products, beverage production and animal fodder. Genetic diversity is a key component in breeding programs. We have analyzed the genetic diversity of barley accessions using microsatellite markers. The accessions were composed of wild and domesticated barley representing genotypes from six countries and three breeding programs in Brazil. A total of 280 alleles were detected, 36 unique to Brazilian barley. The marker Bmag120 showed the greatest polymorphism information content (PIC), with the highest mean value found on chromosome three, and the lowest on chromosomes four and six. The wild accessions presented the highest diversity followed by the foreign genotypes. Genetic analysis was performed using Principal Coordinates Analysis, UPGMA clustering, and Bayesian clustering analysis implemented in Structure. All results obtained by the different methods were similar. Loss of genetic diversity has occurred in Brazilian genotypes. The number of alleles detected in genotypes released in 1980s was higher, whereas most of the cultivars released thereafter showed lower PIC and clustered in separate subgroups from the older cultivars. The use of a more diverse panel of genotypes should be considered in order to exploit novel alleles in Brazilian barley breeding programs.

Genetic Diversity and Association of EST-SSR and SCoT Markers with Rust Traits in Orchardgrass (Dactylis glomerata L.)

Orchardgrass (Dactylis glomerata L.), is a well-known perennial forage species; however, rust diseases have caused a noticeable reduction in the quality and production of orchardgrass. In this study, genetic diversity was assessed and the marker-trait associations for rust were examined using 18 EST-SSR and 21 SCoT markers in 75 orchardgrass accessions. A high level of genetic diversity was detected in orchardgrass with an average genetic diversity index of 0.369. For the EST-SSR and SCoT markers, 164 and 289 total bands were obtained, of which 148 (90.24%) and 272 (94.12%) were polymorphic, respectively. Results from an AMOVA analysis showed that more genetic variance existed within populations (87.57%) than among populations (12.43%). Using a parameter marker index, the efficiencies of the EST-SSR and SCoT markers were compared to show that SCoTs have higher marker efficiency (8.07) than EST-SSRs (4.82). The results of a UPGMA cluster analysis and a STRUCTURE analysis were both correlated with the geographic distribution of the orchardgrass accessions. Linkage disequilibrium analysis revealed an average r² of 0.1627 across all band pairs, indicating a high extent of linkage disequilibrium in the material. An association analysis between the rust trait and 410 bands from the EST-SSR and SCoT markers using TASSEL software revealed 20 band panels were associated with the rust trait in both 2011 and 2012. The 20 bands obtained from association analysis could be used in breeding programs for lineage selection to prevent great losses of orchardgrass caused by rust, and provide valuable information for further association mapping using this collection of orchardgrass.

Identification of Genetic Loci Associated with Quality Traits in Almond via Association Mapping

To design an appropriate association study, we need to understand population structure and the structure of linkage disequilibrium within and among populations as well as in different regions of the genome in an organism. In this study, we have used a total of 98 almond accessions, from five continents located and maintained at the Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA; Spain), and 40 microsatellite markers. Population structure analysis performed in 'Structure' grouped the accessions into two principal groups; the Mediterranean (Western-Europe) and the non-Mediterranean, with K = 3, being the best fit for our data. There was a strong subpopulation structure with linkage disequilibrium decaying with increasing genetic distance resulting in lower levels of linkage disequilibrium between more distant markers. A significant impact of population structure on linkage disequilibrium in the almond cultivar groups was observed. The mean r2 value for all intra-chromosomal loci pairs was 0.040, whereas, the r2 for the inter-chromosomal loci pairs was 0.036. For analysis of association between the markers and phenotypic traits, five models comprising both general linear models and mixed linear models were selected to test the marker trait associations. The mixed linear model (MLM) approach using co-ancestry values from population structure and kinship estimates (K model) as covariates identified a maximum of 16 significant associations for chemical traits and 12 for physical traits. This study reports for the first time the use of association mapping for determining marker-locus trait associations in a world-wide almond germplasm collection. It is likely that association mapping will have the most immediate and largest impact on the tier of crops such as almond with the greatest economic value.

A genome-wide association study of malting quality across eight U.S. barley breeding programs

We report malt quality QTLs relevant to breeding with greater precision than previous mapping studies. The distribution of favorable alleles suggests strategies for marker-assisted breeding and germplasm exchange. This study leverages the breeding data of 1,862 barley breeding lines evaluated in 97 field trials for genome-wide association study of malting quality traits in barley. The mapping panel consisted of six-row and two-row advanced breeding lines from eight breeding populations established at six public breeding programs across the United States. A total of 4,976 grain samples were subjected to micro-malting analysis and mapping of nine quality traits was conducted with 3,072 SNP markers distributed throughout the genome. Association mapping was performed for individual breeding populations and for combined six-row and two-row populations. Only 16% of the QTL we report here had been detected in prior bi-parental mapping studies. Comparison of the analyses of the combined two-row and six-row panels identified only two QTL regions that were common to both. In total, 108 and 107 significant marker-trait associations were identified in all six-row and all two-row breeding programs, respectively. A total of 102 and 65 marker-trait associations were specific to individual six-row and two-row breeding programs, respectively indicating that most marker-trait associations were breeding population specific. Combining datasets from different breeding program resulted in both the loss of some QTL that were apparent in the analyses of individual programs and the discovery of new QTL not identified in individual programs. This suggests that simply increasing sample size by pooling samples with different breeding history does not necessarily increase the power to detect associations. The genetic architecture of malting quality and the distribution of favorable alleles suggest strategies for marker-assisted selection and germplasm exchange.

Mapping of yellow mosaic virus (YMV) resistance in soybean (Glycine max L. Merr.) through association mapping approach

Yellow Mosaic Virus (YMV) is a serious disease of soybean. Resistance to YMV was mapped in 180 soybean genotypes through association mapping approach using 121 simple sequence repeats (SSR) and four resistance gene analogue (RGA)-based markers. The association mapping population (AMP) (96 genotypes) and confirmation population (CP) (84 genotypes) was tested for resistance to YMV at hot-spot consecutively for 3 years (2007-2009). The genotypes exhibited significant variability for YMV resistance (P < 0.01). Molecular genotyping and population structure analysis with 'admixture' co-ancestry model detected seven optimal sub-populations in the AMP. Linkage disequilibrium (LD) between the markers extended up to 35 and 10 cM with r2 > 0.15, and >0.25, respectively. The 4 RGA-based markers showed no association with YMV resistance. Two SSR markers, Satt301 and GMHSP179 on chromosome 17 were found to be in significant LD with YMV resistance. Contingency Chi-square test confirmed the association (P < 0.01) and the utility of the markers was validated in the CP. It would pave the way for marker assisted selection for YMV resistance in soybean. This is the first report of its kind in soybean.

Association mapping for kernel phytosterol content in almond

Almond kernels are a rich source of phytosterols, which are important compounds for human nutrition. The genetic control of phytosterol content has not yet been documented in almond. Association mapping (AM), also known as linkage disequilibrium (LD), was applied to an almond germplasm collection in order to provide new insight into the genetic control of total and individual sterol contents in kernels. Population structure analysis grouped the accessions into two principal groups, the Mediterranean and the non-Mediterranean. There was a strong subpopulation structure with LD decaying with increasing genetic distance, resulting in lower levels of LD between more distant markers. A significant impact of population structure on LD in the almond cultivar groups was observed. The mean r(2) -value for all intra-chromosomal loci pairs was 0.040, whereas, the r(2) for the inter-chromosomal loci pairs was 0.036. For analysis of association between the markers and phenotypic traits five models were tested. The mixed linear model (MLM) approach using co-ancestry values from population structure and kinship estimates (K model) as covariates identified a maximum of 13 significant associations. Most of the associations found appeared to map within the interval where many candidate genes involved in the sterol biosynthesis pathway are predicted in the peach genome. These findings provide a valuable foundation for quality gene identification and molecular marker assisted breeding in almond.

Genetic diversity and population structure in a legacy collection of spring barley landraces adapted to a wide range of climates

Global environmental change and increasing human population emphasize the urgent need for higher yielding and better adapted crop plants. One strategy to achieve this aim is to exploit the wealth of so called landraces of crop species, representing diverse traditional domesticated populations of locally adapted genotypes. In this study, we investigated a comprehensive set of 1485 spring barley landraces (Lrc1485) adapted to a wide range of climates, which were selected from one of the largest genebanks worldwide. The landraces originated from 5° to 62.5° N and 16° to 71° E. The whole collection was genotyped using 42 SSR markers to assess the genetic diversity and population structure. With an average allelic richness of 5.74 and 372 alleles, Lrc1485 harbours considerably more genetic diversity than the most polymorphic current GWAS panel for barley. Ten major clusters defined most of the population structure based on geographical origin, row type of the ear and caryopsis type – and were assigned to specific climate zones. The legacy core reference set Lrc648 established in this study will provide a long-lasting resource and a very valuable tool for the scientific community. Lrc648 is best suited for multi-environmental field testing to identify candidate genes underlying quantitative traits but also for allele mining approaches.

Population structure and association mapping studies for important agronomic traits in soybean

The present study was carried out with a set of 96 diverse soybean genotypes with the objectives of analysing the population structure and to identify molecular markers associated with important agronomic traits. Large phenotypic variability was observed for the agronomic traits under study indicating suitability of the genotypes for association studies. The maximum values for plant height, pods per plant, seeds per pod, 100-seed weight and seed yield per plant were approximately two and half to three times more than the minimum values for the genotypes. Seed yield per plant was found to be significantly correlated with pods per plant (r = 0.77), 100-seed weight (r = 0.35) and days to maturity (r = 0.23). The population structure studies depicted the presence of seven subpopulations which nearly corresponded with the source of geographical origin of the genotypes. Linkage disequilibrium (LD) between the linked markers decreased with the increased distance, and a substantial drop in LD decay values was observed between 30 and 35 cM. Genomewide marker-traits association analysis carried out using general linear (GLM) and mixed linear models (MLM) identified six genomic regions (two of them were common in both) on chromosomes 6, 7, 8, 13, 15 and 17, which were found to be significantly associated with various important traits viz., plant height, pods per plant, 100-seed weight, plant growth habit, average number of seeds per pod, days to 50% flowering and days to maturity. The phenotypic variation explained by these loci ranged from 6.09 to 13.18% and 4.25 to 9.01% in the GLM and MLM studies, respectively. In conclusion, association mapping (AM) in soybean could be a viable alternative to conventional QTL mapping approach.

Detection of QTLs for salt tolerance in Asian barley (Hordeum vulgare L.) by association analysis with SNP markers

Two hundred ninety-six Asian barley (Hordeum vulgare L.) accessions were assessed to detect QTLs underlying salt tolerance by association analysis using a 384 single nucleotide polymorphism (SNP) marker system. The experiment was laid out at the seedling stage in a hydroponic solution under control and 250 mM NaCl solution with three replications of four plants each. Salt tolerance was assessed by leaf injury score (LIS) and salt tolerance indices (STIs) of the number of leaves (NL), shoot length (SL), root length (RL), shoot dry weight (SDW) and root dry weight (RDW). LIS was scored from 1 to 5 according to the severity of necrosis and chlorosis observed on leaves. There was a wide variation in salt tolerance among Asian barley accessions. LIS and STI (SDW) were the most suitable traits for screening salt tolerance. Association was estimated between markers and traits to detect QTLs for LIS and STI (SDW). Seven significant QTLs were located on chromosomes 1H (2 QTLs), 2H (2 QTLs), 3H (1 QTL), 4H (1 QTL) and 5H (1 QTL). Five QTLs were associated with LIS and 2 QTLs with STI (SDW). Two QTLs associated with LIS were newly identified on chromosomes 3H and 4H.

Promoting utilization of Saccharum spp. genetic resources through genetic diversity analysis and core collection construction

Sugarcane (Saccharum spp.) and other members of Saccharum spp. are attractive biofuel feedstocks. One of the two World Collections of Sugarcane and Related Grasses (WCSRG) is in Miami, FL. This WCSRG has 1002 accessions, presumably with valuable alleles for biomass, other important agronomic traits, and stress resistance. However, the WCSRG has not been fully exploited by breeders due to its lack of characterization and unmanageable population. In order to optimize the use of this genetic resource, we aim to 1) genotypically evaluate all the 1002 accessions to understand its genetic diversity and population structure and 2) form a core collection, which captures most of the genetic diversity in the WCSRG. We screened 36 microsatellite markers on 1002 genotypes and recorded 209 alleles. Genetic diversity of the WCSRG ranged from 0 to 0.5 with an average of 0.304. The population structure analysis and principal coordinate analysis revealed three clusters with all S. spontaneum in one cluster, S. officinarum and S. hybrids in the second cluster and mostly non-Saccharum spp. in the third cluster. A core collection of 300 accessions was identified which captured the maximum genetic diversity of the entire WCSRG which can be further exploited for sugarcane and energy cane breeding. Sugarcane and energy cane breeders can effectively utilize this core collection for cultivar improvement. Further, the core collection can provide resources for forming an association panel to evaluate the traits of agronomic and commercial importance.

Farmers without borders-genetic structuring in century old barley (Hordeum vulgare)

The geographic distribution of genetic diversity can reveal the evolutionary history of a species. For crop plants, phylogeographic patterns also indicate how seed has been exchanged and spread in agrarian communities. Such patterns are, however, easily blurred by the intense seed trade, plant improvement and even genebank conservation during the twentieth century, and discerning fine-scale phylogeographic patterns is thus particularly challenging. Using historical crop specimens, these problems are circumvented and we show here how high-throughput genotyping of historical nineteenth century crop specimens can reveal detailed geographic population structure. Thirty-one historical and nine extant accessions of North European landrace barley (Hordeum vulgare L.), in total 231 individuals, were genotyped on a 384 single nucleotide polymorphism assay. The historical material shows constant high levels of within-accession diversity, whereas the extant accessions show more varying levels of diversity and a higher degree of total genotype sharing. Structure, discriminant analysis of principal components and principal component analysis cluster the accessions in latitudinal groups across country borders in Finland, Norway and Sweden. FST statistics indicate strong differentiation between accessions from southern Fennoscandia and accessions from central or northern Fennoscandia, and less differentiation between central and northern accessions. These findings are discussed in the context of contrasting historical records on intense within-country south to north seed movement. Our results suggest that although seeds were traded long distances, long-term cultivation has instead been of locally available, possibly better adapted, genotypes.

Genomewide association studies for 50 agronomic traits in peanut using the 'reference set' comprising 300 genotypes from 48 countries of the semi-arid tropics of the world

Peanut is an important and nutritious agricultural commodity and a livelihood of many small-holder farmers in the semi-arid tropics (SAT) of world which are facing serious production threats. Integration of genomics tools with on-going genetic improvement approaches is expected to facilitate accelerated development of improved cultivars. Therefore, high-resolution genotyping and multiple season phenotyping data for 50 important agronomic, disease and quality traits were generated on the ‘reference set’ of peanut. This study reports comprehensive analyses of allelic diversity, population structure, linkage disequilibrium (LD) decay and marker-trait association (MTA) in peanut. Distinctness of all the genotypes can be established by using either an unique allele detected by a single SSR or a combination of unique alleles by two or more than two SSR markers. As expected, DArT features (2.0 alleles/locus, 0.125 PIC) showed lower allele frequency and polymorphic information content (PIC) than SSRs (22.21 alleles /locus, 0.715 PIC). Both marker types clearly differentiated the genotypes of diploids from tetraploids. Multi-allelic SSRs identified three sub-groups (K = 3) while the LD simulation trend line based on squared-allele frequency correlations (r²) predicted LD decay of 15–20 cM in peanut genome. Detailed analysis identified a total of 524 highly significant MTAs (pvalue >2.1×10–6) with wide phenotypic variance (PV) range (5.81–90.09%) for 36 traits. These MTAs after validation may be deployed in improving biotic resistance, oil/ seed/ nutritional quality, drought tolerance related traits, and yield/ yield components.

Genome-wide association studies of agronomic and quality traits in a set of German winter barley (Hordeum vulgare L.) cultivars using Diversity Arrays Technology (DArT)

A set of about 100 winter barley (Hordeum vulgare L.) cultivars, comprising diverse and economically important German barley elite germplasm released during the last six decades, was previously genotypically characterized by single nucleotide polymorphism (SNP) markers using the Illumina GoldenGate BeadArray Technology to detect associations with phenotypic data estimated in three-year field trials at 12 locations. In order to identify further associations and to obtain information on whether the marker type influences the outcome of association genetics studies, the set of winter barley cultivars was re-analyzed using Diversity Arrays Technology (DArT) markers. As with the analysis of the SNPs, only polymorphic markers present at an allele frequency >5% were included to detect associations in a mixed linear model (MLM) approach using the TASSEL software (P ≤ 0.001). The population structure and kinship matrix were estimated on 72 simple sequence repeats (SSRs) covering the whole barley genome. The respective average linkage disequilibrium (LD) analyzed with DArT markers was estimated at 5.73 cM. A total of 52 markers gave significant associations with at least one of the traits estimated which, therefore, may be suitable for marker-assisted breeding. In addition, by comparing the results to those generated using the Illumina GoldenGate BeadArray Technology, it turned out that a different number of associations for respective traits is detected, depending on the marker system. However, as only a few of the respective DArT and Illumina markers are present in a common map, no comprehensive comparison of the detected associations was feasible, but some were probably detected in the same chromosomal regions. Because of the identification of additional marker-trait associations, it may be recommended to use both marker techniques in genome-wide association studies.

Association analysis of fiber quality traits and exploration of elite alleles in Upland cotton cultivars/accessions (Gossypium hirsutum L.)

Exploring the elite alleles and germplasm accessions related to fiber quality traits will accelerate the breeding of cotton for fiber quality improvement. In this study, 99 Gossypium hirsutum L. accessions with diverse origins were used to perform association analysis of fiber quality traits using 97 polymorphic microsatellite marker primer pairs. A total of 107 significant marker-trait associations were detected for three fiber quality traits under three different environments, with 70 detected in two or three environments and 37 detected in only one environment. Among the 70 significant marker-trait associations, 52.86% were reported previously, implying that these are stable loci for target traits. Furthermore, we detected a large number of elite alleles associated simultaneously with two or three traits. These elite alleles were mainly from accessions collected in China, introduced to China from the United States, or rare alleles with a frequency of less than 5%. No one cultivar contained more than half of the elite alleles, but 10 accessions were collected from China and the two introduced from the United States did contain more than half of these alleles. Therefore, there is great potential for mining elite alleles from germplasm accessions for use in fiber quality improvement in modern cotton breeding.

Development of SNP markers for genes of the phenylpropanoid pathway and their association to kernel and malting traits in barley

BACKGROUND

Flavonoids are an important class of secondary compounds in angiosperms. Next to certain biological functions in plants, they play a role in the brewing process and have an effect on taste, color and aroma of beer. The aim of this study was to reveal the haplotype diversity of candidate genes involved in the phenylpropanoid biosynthesis pathway in cultivated barley varieties (Hordeum vulgare L.) and to determine associations to kernel and malting quality parameters.

RESULTS

Five genes encoding phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), chalcone synthase (CHS), flavanone 3-hydroxylase (F3H) and dihydroflavonol reductase (DFR) of the phenylpropanoid biosynthesis pathway were partially resequenced in 16 diverse barley reference genotypes. Their localization in the barley genome, their genetic structure, and their genetic variation e.g. single nucleotide polymorphism (SNP) and Insertion/Deletion (InDel) patterns were revealed. In total, 130 SNPs and seven InDels were detected. Of these, 21 polymorphisms were converted into high-throughput pyrosequencing markers. The resulting SNP and haplotype patterns were used to calculate associations with kernel and malting quality parameters.

CONCLUSIONS

SNP patterns were found to be highly variable for the investigated genes. The developed high-throughput markers are applicable for assessing the genetic variability and for the determination of haplotype patterns in a set of barley accessions. The candidate genes PAL, C4H and F3H were shown to be associated to several malting properties like glassiness (PAL), viscosity (C4H) or to final attenuation (F3H).

Association mapping of salt tolerance in barley (Hordeum vulgare L.)

A spring barley collection of 192 genotypes from a wide geographical range was used to identify quantitative trait loci (QTLs) for salt tolerance traits by means of an association mapping approach using a thousand SNP marker set. Linkage disequilibrium (LD) decay was found with marker distances spanning 2-8 cM depending on the methods used to account for population structure and genetic relatedness between genotypes. The association panel showed large variation for traits that were highly heritable under salt stress, including biomass production, chlorophyll content, plant height, tiller number, leaf senescence and shoot Na(+), shoot Cl(-) and shoot, root Na(+)/K(+) contents. The significant correlations between these traits and salt tolerance (defined as the biomass produced under salt stress relative to the biomass produced under control conditions) indicate that these traits contribute to (components of) salt tolerance. Association mapping was performed using several methods to account for population structure and minimize false-positive associations. This resulted in the identification of a number of genomic regions that strongly influenced salt tolerance and ion homeostasis, with a major QTL controlling salt tolerance on chromosome 6H, and a strong QTL for ion contents on chromosome 4H.

Association mapping combined with linkage analysis for aluminum tolerance among soybean cultivars released in Yellow and Changjiang River Valleys in China

Association mapping (AM) combined with linkage mapping (LM) was executed to identify molecular markers and QTL regions associated with aluminum (Al) tolerance using relative root elongation (RRE) in hydroponics as an indicator. A set of 188 soybean cultivars released in Yellow and Changjiang River Valleys and 184 recombinant inbred lines (RIL) derived from a cross KF No. 1 (tolerant) × NN1138-2 (susceptible) was used in the study. Inheritance analysis of the RIL population suggested four major genes and polygenes controlled Al-tolerance. Further, LM indicated four additive and four epistatic QTL pairs plus a collective unmapped minor QTL were responsible for Al-tolerance and explained 29.39, 18.75 and 43.07 % of the phenotypic variation (PV), respectively. In the set of released cultivars, AM identified 11 markers significant at P < 0.03 that explained 85.2 % of PV with six of which at P < 0.01 accounted for 57.9 % of PV. Ten of these eleven AM marker-QTL were mapped within range of ~2.0 cM to ~43.0 cM outside confidence interval of respective Al-tolerance QTL in previous studies. Five markers, Satt209, Sat_364, Sat_240, Sct_190 and Satt284, were located near Al-tolerance QTL regions in this and previous LM studies. Thus, the two methods confirmed these markers as being the most likely candidate regions for Al-tolerance. Allele effects relative to the population mean for the 11 QTL were estimated, and the allele A210 of Satt209 showed greatest phenotypic effect on Al-tolerance. The two most favorable alleles from each of the 11 marker loci and their carriers were identified, and accordingly the genetic constitution of Al-tolerance for the 188 cultivars was dissected as a QTL-allele matrix. Therefore, marker-assisted pairing of crosses and marker-assisted selection of progenies can be carried out to pyramid favorable alleles of all the 11 loci. This marker-assisted breeding procedure was designated as breeding by design using a QTL-allele matrix.

Polymorphism of TaSAP1-A1 and its association with agronomic traits in wheat

TaSAP1, a member of the stress association protein (SAP) gene family from wheat, is involved in response to several abiotic stresses, including drought, salt and cold. In this study, TaSAP1-A1, a TaSAP1 member on chromosome 7A, and its flanking sequences were isolated. Polymorphism analysis indicated that the average nucleotide diversity (π) of the whole region was 0.00296. The highest nucleotide diversity occurred in the promoter region (π = 0.00631) and no polymorphism was identified in the coding region. Three markers T7AM5, T7AM2606 and T7AM39 located in the promoter region, were developed from sequence variations (InDel5-1810, SNP-2606 and InDel39-1637). Six haplotypes were identified among 300 accessions based on the three markers. TaSAP1-A1 was located on chromosome 7A using marker T7AM39 and was flanked by markers Xwmc530 and Xbarc174. QTL for yield-related traits, including 1,000-grain weight, number of grains per spike and grain yield, were located in the same region. In marker- and haplotype-trait association analyses, TaSAP1-A1 was significantly associated with 1,000-grain weight, number of grains per spike, spike length, peduncle length and total number of spikelets per spike in multiple environments. These results provide useful information for marker-assisted selection for yield-related traits under well-watered and drought-stressed conditions.