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

Diversity of Gibberellin 2-oxidase genes in the barley genome offers opportunities for genetic improvement

INTRODUCTION

Gibberellin (GA) is a vital phytohormone in regulating plant growth and development. During the "Green Revolution", modification of GA-related genes created semi-dwarfing phenotype in cereal crops but adversely affected grain weight. Gibberellin 2-oxidases (GA2oxs) in barley act as key catabolic enzymes in deactivating GA, but their functions are still less known.

OBJECTIVES

This study investigates the physiological function of two HvGA2ox genes in barley and identifies novel semi-dwarf alleles with minimum impacts on other agronomic traits.

METHODS

Virus-induced gene silencing and CRISPR/Cas9 technology were used to manipulate gene expression of HvGA2ox9 and HvGA2ox8a in barley and RNA-seq was conducted to compare the transcriptome between wild type and mutants. Also, field trials in multiple environments were performed to detect the functional haplotypes.

RESULTS

There were ten GA2oxs that distinctly expressed in shoot, tiller, inflorescence, grain, embryo and root. Knockdown of HvGA2ox9 did not affect plant height, while ga2ox8a mutants generated by CRISPR/Cas9 increased plant height and significantly altered seed width and weight due to the increased bioactive GA4 level. RNA-seq analysis revealed that genes involved in starch and sucrose metabolism were significantly decreased in the inflorescence of ga2ox8a mutants. Furthermore, haplotype analysis revealed one naturally occurring HvGA2ox8a haplotype was associated with decreased plant height, early flowering and wider and heavier seed.

CONCLUSION

Our results demonstrate the potential of manipulating GA2ox genes to fine tune GA signalling and biofunctions in desired plant tissues and open a promising avenue for minimising the trade-off effects of Green Revolution semi-dwarfing genes on grain size and weight. The knowledge will promote the development of next generation barley cultivars with better adaptation to a changing climate.

Sense in sensitivity: difference in the meaning of photoperiod-insensitivity between wheat and barley

The description of long photoperiod sensitivity in wheat and barley is a cause of confusion for researchers working in these crops, usually accustomed to free exchange of physiological and genetic knowledge of such similar crops. Indeed, wheat and barley scientists customarily quote studies of either crop species when researching one of them. Among their numerous similarities the main gene controlling that response is the same in both crops (PPD1; PPD-H1 in barley and PPD-D1 in hexaploid wheat). However, the photoperiod responses are different: (i) the main dominant allele inducing shorter time to anthesis is the insensitive allele in wheat (Ppd-D1a) but the sensitive allele in barley (Ppd-H1) (i.e. sensitivity to photoperiod produces opposite effects on time to heading in wheat and barley), (ii) the main "insensitive" allele in wheat, Ppd-D1a, does confer insensitivity, whilst that of barley reduces the sensitivity but still responds to photoperiod. The different behaviour of PPD1 genes in wheat and barley is put in a common framework based on the similarities and differences of the molecular bases of their mutations, which include polymorphism at gene expression levels, copy number variation, and sequence of coding regions. This common perspective sheds light on a source on confusion for cereal researchers, and prompts us to recommend accounting for the photoperiod sensitivity status of the plant materials when doing research on genetic control of phenology. Finally, we provide advice to facilitate the management of natural PPD1 diversity in breeding programs and suggest targets for further modification through gene editing, based on mutual knowledge on the two crops.

New semi-dwarfing alleles with increased coleoptile length by gene editing of gibberellin 3-oxidase 1 using CRISPR-Cas9 in barley (Hordeum vulgare L.)

The green revolution was based on genetic modification of the gibberellin (GA) hormone system with "dwarfing" gene mutations that reduces GA signals, conferring shorter stature, thus enabling plant adaptation to modern farming conditions. Strong GA-related mutants with shorter stature often have reduced coleoptile length, discounting yield gain due to their unsatisfactory seedling emergence under drought conditions. Here we present gibberellin (GA) 3-oxidase1 (GA3ox1) as an alternative semi-dwarfing gene in barley that combines an optimal reduction in plant height without restricting coleoptile and seedling growth. Using large-scale field trials with an extensive collection of barley accessions, we showed that a natural GA3ox1 haplotype moderately reduced plant height by 5-10 cm. We used CRISPR/Cas9 technology, generated several novel GA3ox1 mutants and validated the function of GA3ox1. We showed that altered GA3ox1 activities changed the level of active GA isoforms and consequently increased coleoptile length by an average of 8.2 mm, which could provide essential adaptation to maintain yield under climate change. We revealed that CRISPR/Cas9-induced GA3ox1 mutations increased seed dormancy to an ideal level that could benefit the malting industry. We conclude that selecting HvGA3ox1 alleles offers a new opportunity for developing barley varieties with optimal stature, longer coleoptile and additional agronomic traits.

Preadapted to adapt: underpinnings of adaptive plasticity revealed by the downy brome genome

Bromus tectorum L. is arguably the most successful invasive weed in the world. It has fundamentally altered arid ecosystems of the western United States, where it now found on an excess of 20 million hectares. Invasion success is related to avoidance of abiotic stress and human management. Early flowering is a heritable trait utilized by B. tectorum, enabling the species to temporally monopolize limited resources and outcompete the native plant community. Thus, understanding the genetic underpinning of flowering time is critical for the design of integrated management strategies. To study flowering time traits in B. tectorum, we assembled a chromosome scale reference genome for B. tectorum. To assess the utility of the assembled genome, 121 diverse B. tectorum accessions are phenotyped and subjected to a genome wide association study (GWAS). Candidate genes, representing homologs of genes that have been previously associated with plant height or flowering phenology traits in related species are located near QTLs we identified. This study uses a high-resolution GWAS to identify reproductive phenology genes in a weedy species and represents a considerable step forward in understanding the mechanisms underlying genetic plasticity in one of the most successful invasive weed species.

Genome architecture and diverged selection shaping pattern of genomic differentiation in wild barley

Divergent selection of populations in contrasting environments leads to functional genomic divergence. However, the genomic architecture underlying heterogeneous genomic differentiation remains poorly understood. Here, we de novo assembled two high-quality wild barley (Hordeum spontaneum K. Koch) genomes and examined genomic differentiation and gene expression patterns under abiotic stress in two populations. These two populations had a shared ancestry and originated in close geographic proximity but experienced different selective pressures due to their contrasting micro-environments. We identified structural variants that may have played significant roles in affecting genes potentially associated with well-differentiated phenotypes such as flowering time and drought response between two wild barley genomes. Among them, a 29-bp insertion into the promoter region formed a cis-regulatory element in the HvWRKY45 gene, which may contribute to enhanced tolerance to drought. A single SNP mutation in the promoter region may influence HvCO5 expression and be putatively linked to local flowering time adaptation. We also revealed significant genomic differentiation between the two populations with ongoing gene flow. Our results indicate that SNPs and small SVs link to genetic differentiation at the gene level through local adaptation and are maintained through divergent selection. In contrast, large chromosome inversions may have shaped the heterogeneous pattern of genomic differentiation along the chromosomes by suppressing chromosome recombination and gene flow. Our research offers novel insights into the genomic basis underlying local adaptation and provides valuable resources for the genetic improvement of cultivated barley.

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

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

Comparative transcriptomic and epigenetic analyses reveal conserved and divergent regulatory pathways in barley response to temperature stresses

DNA methylation and histone modification enable plants to rapidly adapt to adverse temperature stresses, including low temperature (LT) and high temperature (HT) stress. In this study, we conducted physiological, epigenetic, and transcriptomic analyses of barley seedlings grown under control (22°C), mild low temperature (MLT, 14°C) and HT (38°C) conditions to elucidate the underlying molecular mechanisms. Compared to MLT, HT implies greater deleterious effects on barley seedlings' growth. The methylation-sensitive amplification polymorphism analysis showed that MLT induced more DNA methylation and HT more DNA demethylation compared to control. Besides, the higher levels of H3K9ac and H3K4me3 under HT compared to MLT stresses might lead to the loosening of chromatin and, subsequently, the activation of gene expression. Consistently, the transcriptome analysis revealed that there were more differentially expressed genes (DEGs) in plants subjected to HT stress than MLT stress compared to control. The common and unique pathways of these DEGs between MLT and HT were also analyzed. Transcription factors, such as ERF, bHLH, NAC, HSF, and MYB, were most involved in MLT and HT stress. The underlying gene regulation networks of epigenetic modulation-related genes were further explored by weight gene co-expression network analysis. Our study provides new insights into the understanding of epigenetic regulation responses to temperature stress in barley, which will lead to improved strategies for the development of cold- and heat-tolerant barley varieties for sustainable barley production in a climate-changing world.

Bringing more players into play: Leveraging stress in genome wide association studies

In order to meet the demand of the burgeoning human population as well as to adapt crops to the enhanced abiotic and biotic stress caused by the global climatic change, breeders focus on identifying valuable genes to improve both crop stress tolerance and crop quality. Recently, with the development of next-generation sequencing methods, millions of high quality single-nucleotide polymorphisms (SNPs) have been made available and genome-wide association studies (GWAS) are widely used in crop improvement studies to identify the associations between genetic variants of genomes and relevant crop agronomic traits. Here, we review classic cases of use of GWAS to identify genetic variants associated with valuable traits such as geographic adaptation, crop quality and metabolites. We discuss the power of stress GWAS to identify further associations including those with genes that are not, or only lowly, expressed during optimal growth conditions. Finally, we emphasize recent demonstrations of the efficiency and accuracy of time-resolved dynamic stress GWAS and GWAS based on genomic gene expression and structural variations, which can be applied to resolve more comprehensively the genetic regulation mechanisms of complex traits.

Responses of Barley to High Ambient Temperature Are Modulated by Vernalization

Ambient temperatures are increasing due to climate change. Cereal crops development and production will be affected consequently. Flowering time is a key factor for adaptation of small grain cereals and, therefore, exploring developmental responses of barley to rising temperatures is required. In this work, we studied phasic growth, and inflorescence traits related to yield, in eight near isogenic lines of barley (Hordeum vulgare L.) differing at the VRN-H1, VRN-H2 and PPD-H1 genes, representing different growth habits. The lines were grown in contrasting vernalization treatments, under two temperature regimes (18 and 25°C), in long days. Lines with recessive ppd-H1 presented delayed development compared to lines with the sensitive PPD-H1 allele, across the two growth phases considered. High temperature delayed flowering in all unvernalized plants, and in vernalized spring barleys carrying the insensitive ppd-H1 allele, whilst it accelerated flowering in spring barleys with the sensitive PPD-H1 allele. This finding evidenced an interaction between PPD-H1, temperature and vernalization. At the high temperature, PPD-H1 lines in spring backgrounds (VRN-H1-7) yielded more, whereas lines with ppd-H1 were best in vrn-H1 background. Our study revealed new information that will support breeding high-yielding cultivars with specific combinations of major adaptation genes tailored to future climatic conditions.

Novel Design of Imputation-Enabled SNP Arrays for Breeding and Research Applications Supporting Multi-Species Hybridization

Array-based single nucleotide polymorphism (SNP) genotyping platforms have low genotype error and missing data rates compared to genotyping-by-sequencing technologies. However, design decisions used to create array-based SNP genotyping assays for both research and breeding applications are critical to their success. We describe a novel approach applicable to any animal or plant species for the design of cost-effective imputation-enabled SNP genotyping arrays with broad utility and demonstrate its application through the development of the Illumina Infinium Wheat Barley 40K SNP array Version 1.0. We show that the approach delivers high quality and high resolution data for wheat and barley, including when samples are jointly hybridised. The new array aims to maximally capture haplotypic diversity in globally diverse wheat and barley germplasm while minimizing ascertainment bias. Comprising mostly biallelic markers that were designed to be species-specific and single-copy, the array permits highly accurate imputation in diverse germplasm to improve the statistical power of genome-wide association studies (GWAS) and genomic selection. The SNP content captures tetraploid wheat (A- and B-genome) and Aegilops tauschii Coss. (D-genome) diversity and delineates synthetic and tetraploid wheat from other wheat, as well as tetraploid species and subgroups. The content includes SNP tagging key trait loci in wheat and barley, as well as direct connections to other genotyping platforms and legacy datasets. The utility of the array is enhanced through the web-based tool, Pretzel (https://plantinformatics.io/) which enables the content of the array to be visualized and interrogated interactively in the context of numerous genetic and genomic resources to be connected more seamlessly to research and breeding. The array is available for use by the international wheat and barley community.

A global alfalfa diversity panel reveals genomic selection signatures in Chinese varieties and genomic associations with root development

Alfalfa (Medicago sativa L.) is an important forage crop worldwide. However, little is known about the effects of breeding status and different geographical populations on alfalfa improvement. Here, we sequenced 220 alfalfa core germplasms and determined that Chinese alfalfa cultivars form an independent group, as evidenced by comparisons of F_ST values between different subgroups, suggesting that geographical origin plays an important role in group differentiation. By tracing the influence of geographical regions on the genetic diversity of alfalfa varieties in China, we identified 350 common candidate genetic regions and 548 genes under selection. We also defined 165 loci associated with 24 important traits from genome-wide association studies. Of those, 17 genomic regions closely associated with a given phenotype were under selection, with the underlying haplotypes showing significant differences between subgroups of distinct geographical origins. Based on results from expression analysis and association mapping, we propose that 6-phosphogluconolactonase (MsPGL) and a gene encoding a protein with NHL domains (MsNHL) are critical candidate genes for root growth. In conclusion, our results provide valuable information for alfalfa improvement via molecular breeding.

Genomic structural equation modelling provides a whole-system approach for the future crop breeding

Using genomic structural equation modelling, this research demonstrates an efficient way to identify genetically correlating traits and provides an effective proxy for multi-trait selection to consider the joint genetic architecture of multiple interacting traits in crop breeding. Breeding crop cultivars with optimal value across multiple traits has been a challenge, as traits may negatively correlate due to pleiotropy or genetic linkage. For example, grain yield and grain protein content correlate negatively with each other in cereal crops. Future crop breeding needs to be based on practical yet accurate evaluation and effective selection of beneficial trait to retain genes with the best agronomic score for multiple traits. Here, we test the framework of whole-system-based approach using structural equation modelling (SEM) to investigate how one trait affects others to guide the optimal selection of a combination of agronomically important traits. Using ten traits and genome-wide SNP profiles from a worldwide barley panel and SEM analysis, we revealed a network of interacting traits, in which tiller number contributes positively to both grain yield and protein content; we further identified common genetic factors affecting multiple traits in the network of interaction. Our method demonstrates an efficient way to identify genetically correlating traits and underlying pleiotropic genetic factors and provides an effective proxy for multi-trait selection within a whole-system framework that considers the joint genetic architecture of multiple interacting traits in crop breeding. Our findings suggest the promise of a whole-system approach to overcome challenges such as the negative correlation of grain yield and protein content to facilitating quantitative and objective breeding decisions in future crop breeding.

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.

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.

Highly efficient and genotype-independent barley gene editing based on anther culture

Recalcitrance to tissue culture and genetic transformation is the major bottleneck for gene manipulation in crops. In barley, immature embryos of Golden Promise have typically been used as explants for transformation. However, the genotype dependence of this approach limits the genetic modification of commercial varieties. Here, we developed an anther culture-based system that permits the effective creation of transgenic and gene-edited plants from commercial barley varieties. The protocol was tested in Golden Promise and four Australian varieties, which differed in phenology, callus induction, and green plant regeneration responses. Agrobacterium-mediated transformation was performed on microspore-derived callus to target the HvPDS gene, and T0 albinos with targeted mutations were successfully obtained from commercial varieties. Further editing of three targets was achieved with an average mutation rate of 53% in the five varieties. In 51 analyzed T0 individuals, Cas9 induced a large proportion (69%) of single-base indels and two-base deletions in the target sites, with variable mutation rates among targets and varieties. Both on-target and off-target activities were detected in T1 progenies. Compared with immature embryo protocols, this genotype-independent platform can deliver a high editing efficiency and more regenerant plants within a similar time frame. It shows promise for functional genomics and the application of CRISPR technologies for the precise improvement of commercial varieties.

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

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

Genome-Wide Association Study and Identification of Candidate Genes for Nitrogen Use Efficiency in Barley (Hordeum vulgare L.)

Nitrogen (N) fertilizer is largely responsible for barley grain yield potential and quality, yet excessive application leads to environmental pollution and high production costs. Therefore, efficient use of N is fundamental for sustainable agriculture. In the present study, we investigated the performance of 282 barley accessions through hydroponic screening using optimal and low NH4NO3 treatments. Low-N treatment led to an average shoot dry weight reduction of 50%, but there were significant genotypic differences among the accessions. Approximately 20% of the genotypes showed high (>75%) relative shoot dry weight under low-N treatment and were classified as low-N tolerant, whereas 20% were low-N sensitive (≤55%). Low-N tolerant accessions exhibited well-developed root systems with an average increase of 60% in relative root dry weight to facilitate more N absorption. A genome-wide association study (GWAS) identified 66 significant marker trait associations (MTAs) conferring high nitrogen use efficiency, four of which were stable across experiments. These four MTAs were located on chromosomes 1H(1), 3H(1), and 7H(2) and were associated with relative shoot length, relative shoot and root dry weight. Genes corresponding to the significant MTAs were retrieved as candidate genes, including members of the asparagine synthetase gene family, several transcription factor families, protein kinases, and nitrate transporters. Most importantly, the high-affinity nitrate transporter 2.7 (HvNRT2.7) was identified as a promising candidate on 7H for root and shoot dry weight. The identified candidate genes provide new insights into our understanding of the molecular mechanisms driving nitrogen use efficiency in barley and represent potential targets for genetic improvement.

Genome-wide association mapping reveals novel genes associated with coleoptile length in a worldwide collection of barley

BACKGROUND

Drought is projected to become more frequent and severe in a changing climate, which requires deep sowing of crop seeds to reach soil moisture. Coleoptile length is a key agronomic trait in cereal crops such as barley, as long coleoptiles are linked to drought tolerance and improved seedling establishment under early water-limited growing conditions.

RESULTS

In this study, we detected large genetic variation in a panel of 328 diverse barley (Hordeum vulgare L.) accessions. To understand the overall genetic basis of barley coleoptile length, all accessions were germinated in the dark and phenotyped for coleoptile length after 2 weeks. The investigated barleys had significant variation for coleoptile length. We then conducted genome-wide association studies (GWASs) with more than 30,000 molecular markers and identified 8 genes and 12 intergenic loci significantly associated with coleoptile length in our barley panel. The Squamosa promoter-binding-like protein 3 gene (SPL3) on chromosome 6H was identified as a major candidate gene. The missense variant on the second exon changed serine to alanine in the conserved SBP domain, which likely impacted its DNA-binding activity.

CONCLUSION

This study provides genetic loci for seedling coleoptile length along with candidate genes for future potential incorporation in breeding programmes to enhance early vigour and yield potential in water-limited environments.

Characterization of Genetic Diversity and Genome-Wide Association Mapping of Three Agronomic Traits in Qingke Barley (Hordeum Vulgare L.) in the Qinghai-Tibet Plateau

Barley (Hordeum vulgare L.) is one of the most important cereal crops worldwide. In the Qinghai-Tibet Plateau, six-rowed hulless (or naked) barley, called “qingke” in Chinese or “nas” in Tibetan, is produced mainly in Tibet. The complexity of the environment in the Qinghai-Tibet Plateau has provided unique opportunities for research on the breeding and adaptability of qingke barley. However, the genetic architecture of many important agronomic traits for qingke barley remains elusive. Heading date (HD), plant height (PH), and spike length (SL) are three prominent agronomic traits in barley. Here, we used genome-wide association (GWAS) mapping and GWAS with eigenvector decomposition (EigenGWAS) to detect quantitative trait loci (QTL) and selective signatures for HD, PH, and SL in a collection of 308 qingke barley accessions. The accessions were genotyped using a newly-developed, proprietary genotyping-by-sequencing (tGBS) technology, that yielded 14,970 high quality single nucleotide polymorphisms (SNPs). We found that the number of SNPs was higher in the varieties than in the landraces, which suggested that Tibetan varieties and varieties in the Tibetan area may have originated from different landraces in different areas. We have identified 62 QTLs associated with three important traits, and the observed phenotypic variation is well-explained by the identified QTLs. We mapped 114 known genes that include, but are not limited to, vernalization, and photoperiod genes. We found that 83.87% of the identified QTLs are located in the non-coding regulatory regions of annotated barley genes. Forty-eight of the QTLs are first reported here, 28 QTLs have pleotropic effects, and three QTL are located in the regions of the well-characterized genes HvVRN1, HvVRN3, and PpD-H2. EigenGWAS analysis revealed that multiple heading-date-related loci bear signatures of selection. Our results confirm that the barley panel used in this study is highly diverse, and showed a great promise for identifying the genetic basis of adaptive traits. This study should increase our understanding of complex traits in qingke barley, and should facilitate genome-assisted breeding for qingke barley improvement.

Genome-Wide Association Study of Salinity Tolerance During Germination in Barley (Hordeum vulgare L.)

Barley seeds need to be able to germinate and establish seedlings in saline soils in Mediterranean-type climates. Despite being a major cereal crop, barley has few reported quantitative trait loci (QTL) and candidate genes underlying salt tolerance at the germination stage. Breeding programs targeting salinity tolerance at germination require an understanding of genetic loci and alleles in the current germplasm. In this study, we investigated seed-germination-related traits under control and salt stress conditions in 350 diverse barley accessions. A genome-wide association study, using ~24,000 genetic markers, was undertaken to detect marker-trait associations (MTA) and the underlying candidate genes for salinity tolerance during germination. We detected 19 loci containing 52 significant salt-tolerance-associated markers across all chromosomes, and 4 genes belonging to 4 family functions underlying the predicted MTAs. Our results provide new genetic resources and information to improve salt tolerance at germination in future barley varieties via genomic and marker-assisted selection and to open up avenues for further functional characterization of the identified candidate genes.

Gene-set association and epistatic analyses reveal complex gene interaction networks affecting flowering time in a worldwide barley collection

Single-marker genome-wide association studies (GWAS) have successfully detected associations between single nucleotide polymorphisms (SNPs) and agronomic traits such as flowering time and grain yield in barley. However, the analysis of individual SNPs can only account for a small proportion of genetic variation, and can only provide limited knowledge on gene network interactions. Gene-based GWAS approaches provide enormous opportunity both to combine genetic information and to examine interactions among genetic variants. Here, we revisited a previously published phenotypic and genotypic data set of 895 barley varieties grown in two years at four different field locations in Australia. We employed statistical models to examine gene-phenotype associations, as well as two-way epistasis analyses to increase the capability to find novel genes that have significant roles in controlling flowering time in barley. Genetic associations were tested between flowering time and corresponding genotypes of 174 putative flowering time-related genes. Gene-phenotype association analysis detected 113 genes associated with flowering time in barley, demonstrating the unprecedented power of gene-based analysis. Subsequent two-way epistasis analysis revealed 19 pairs of gene×gene interactions involved in controlling flowering time. Our study demonstrates that gene-based association approaches can provide higher capacity for future crop improvement to increase crop performance and adaptation to different environments.

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.

Targeted enrichment by solution-based hybrid capture to identify genetic sequence variants in barley

In barley and other cereal crops, phenological diversity drives adaptation to different cultivation areas. Improvement of barley yield and quality traits requires adaptation to specific production areas with introgression of favorable alleles dependent upon precise identification of the underlying genes. Combining targeted sequence capture systems with next-generation sequencing provides an efficient approach to explore target genetic regions at high resolution, and allows rapid discovery of thousands of genetic polymorphisms. Here, we apply a versatile target-capture method to detect genome-wide polymorphisms in 174 flowering time-related genes, chosen based on prior knowledge from barley, rice, and Arabidopsis thaliana. Sequences were generated across a phenologically diverse panel of 895 barley varieties, resulting a high mean depth coverage of ~25x allowing reliable discovery and calling of insertion-deletion (InDel) and single nucleotide polymorphisms (SNPs). Sequences of InDel and SNPs from the targeted enrichment were utilized to develop 67 Kompetitive Allele Specific PCR (KASP) markers for validation. This work provides researchers and breeders a comprehensive molecular toolkit for the selection of phenology-related traits in barley.

Design Type(s)	genetic polymorphism analysis objective • strain comparison design
Measurement Type(s)	genetic variation
Technology Type(s)	DNA sequencing
Factor Type(s)	cultivar
Sample Characteristic(s)	Hordeum vulgare subsp. vulgare • leaf

Machine-accessible metadata file describing the reported data (ISA-Tab format)