Anchoring and ordering NGS contig assemblies by population sequencing (POPSEQ)

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

The genetics of Ug99 stem rust resistance in spring wheat variety 'Linkert'

Wheat stem rust caused by Puccinia graminis f. sp. tritici (Pgt) threatens wheat production worldwide. The objective of this study was to characterize wheat stem rust resistance in 'Linkert', a variety with adult plant resistance effective to emerging wheat stem rust pathogen strain Ug99. Two doubled haploid (DH) populations and one recombinant inbred line (RIL) population were developed with 'Linkert' as a stem rust resistant parent. Hard red spring wheat variety 'Forefront' and genetic stock 'LMPG' were used as stem rust susceptible parents of the DH populations. Breeding line 'MN07098-6' was used as a susceptible parent of the RIL population. Both DH and RIL populations with their parents were evaluated both at the seedling stage and in the field against Pgt races. Genotyping data of the DH populations were generated using the wheat iSelect 90k SNP assay. The RIL population was genotyped by genotyping-by-sequencing. We found QTL consistently associated with wheat stem rust resistance on chromosome 2BS for the Linkert/Forefront DH population and the Linkert/MN07098-6 RIL population both in Ethiopia and Kenya. Additional reliable QTL were detected on chromosomes 5BL (125.91 cM) and 4AL (Sr7a) for the Linkert/LMPG population in Ethiopia and Kenya. Different QTL identified in the populations reflect the importance of examining the genetics of resistance in populations derived from adapted germplasm (Forefront and MN07098-6) in addition to a genetic stock (LMPG). The associated markers in this study could be used to track and select for the identified QTL in wheat breeding programs.

High-resolution mapping of Ryd4Hb, a major resistance gene to Barley yellow dwarf virus from Hordeum bulbosum

KEY MESSAGE

We mapped Ryd4Hb in a 66.5 kbp interval in barley and dissociated it from a sublethality factor. These results will enable a targeted selection of the resistance in barley breeding. Virus diseases are causing high yield losses in crops worldwide. The Barley yellow dwarf virus (BYDV) complex is responsible for one of the most widespread and economically important viral diseases of cereals. While no gene conferring complete resistance (immunity) has been uncovered in the primary gene pool of barley, sources of resistance were searched and identified in the wild relative Hordeum bulbosum, representing the secondary gene pool of barley. One such locus, Ryd4Hb, has been previously introgressed into barley, and was allocated to chromosome 3H, but is tightly linked to a sublethality factor that prevents the incorporation and utilization of Ryd4Hb in barley varieties. To solve this problem, we fine-mapped Ryd4Hb and separated it from this negative factor. We narrowed the Ryd4Hb locus to a corresponding 66.5 kbp physical interval in the barley 'Morex' reference genome. The region comprises a gene from the nucleotide-binding and leucine-rich repeat immune receptor family, typical of dominant virus resistance genes. The closest homolog to this Ryd4Hb candidate gene is the wheat Sr35 stem rust resistance gene. In addition to the fine mapping, we reduced the interval bearing the sublethality factor to 600 kbp in barley. Aphid feeding experiments demonstrated that Ryd4Hb provides a resistance to BYDV rather than to its vector. The presented results, including the high-throughput molecular markers, will permit a more targeted selection of the resistance in breeding, enabling the use of Ryd4Hb in barley varieties.

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

KEY MESSAGE

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

Identifying developmental QTL alleles with favorable effect on grain yield components under late-terminal drought in spring barley MAGIC population

Barley is the fourth most cultivated cereal worldwide, and drought is a major cause of its yield loss by negatively affecting its development. Hence, better understanding developmental mechanisms that control complex polygenic yield-related traits under drought is essential to uncover favorable yield regulators. This study evaluated seven above-ground yield-related traits under well-watered (WW) and late-terminal drought (TD) treatment using 534 spring barley multiparent advanced generation intercross double haploid (DH) lines. The analysis of quantitative trait loci (QTL) for WW, TD, marker by treatment interaction, and drought stress tolerance identified 69, 64, 25, and 25 loci, respectively, for seven traits from which 15 loci were common for at least three traits and 17 were shared by TD and drought stress tolerance. Evaluation of allelic effects for a QTL revealed varying effect of parental alleles. Results showed prominent QTL located on major flowering time gene Ppd-H1 with favorable effects for grain weight under TD when flowering time was not significantly affected, suggesting that this gene might be linked with increasing grain weight by ways other than timing of flowering under late-terminal drought stress. Furthermore, a desirable novel QTL allele was identified on chromosome 5H for grain number under TD nearby sucrose transporter gene HvSUT2. The findings indicated that spring barley multiparent advanced generation intercross population can provide insights to improve yield under complex condition of drought.

Genome-wide association analysis provides insights into the genetic basis of photosynthetic responses to low-temperature stress in spring barley

Low-temperature stress (LTS) is among the major abiotic stresses affecting the geographical distribution and productivity of the most important crops. Understanding the genetic basis of photosynthetic variation under cold stress is necessary for developing more climate-resilient barley cultivars. To that end, we investigated the ability of chlorophyll fluorescence parameters (F_VF_M, and F_VF₀) to respond to changes in the maximum quantum yield of Photosystem II photochemistry as an indicator of photosynthetic energy. A panel of 96 barley spring cultivars from different breeding zones of Canada was evaluated for chlorophyll fluorescence-related traits under cold acclimation and freeze shock stresses at different times. Genome-wide association studies (GWAS) were performed using a mixed linear model (MLM). We identified three major and putative genomic regions harboring 52 significant quantitative trait nucleotides (QTNs) on chromosomes 1H, 3H, and 6H for low-temperature tolerance. Functional annotation indicated several QTNs were either within the known or close to genes that play important roles in the photosynthetic metabolites such as abscisic acid (ABA) signaling, hydrolase activity, protein kinase, and transduction of environmental signal transduction at the posttranslational modification levels. These outcomes revealed that barley plants modified their gene expression profile in response to decreasing temperatures resulting in physiological and biochemical modifications. Cold tolerance could influence a long-term adaption of barley in many parts of the world. Since the degree and frequency of LTS vary considerably among production sites. Hence, these results could shed light on potential approaches for improving barley productivity under low-temperature stress.

The TCP transcription factor HvTB2 heterodimerizes with VRS5 and controls spike architecture in barley

Understanding the molecular network, including protein-protein interactions, of VRS5 provide new routes towards the identification of other key regulators of plant architecture in barley. The TCP transcriptional regulator TEOSINTE BRANCHED 1 (TB1) is a key regulator of plant architecture. In barley, an important cereal crop, HvTB1 (also referred to as VULGARE SIX-ROWED spike (VRS) 5), inhibits the outgrowth of side shoots, or tillers, and grains. Despite its key role in barley development, there is limited knowledge on the molecular network that is utilized by VRS5. In this work, we performed protein-protein interaction studies of VRS5. Our analysis shows that VRS5 potentially interacts with a diverse set of proteins, including other class II TCP's, NF-Y TF, but also chromatin remodelers. Zooming in on the interaction capacity of VRS5 with other TCP TFs shows that VRS5 preferably interacts with other class II TCP TFs in the TB1 clade. Induced mutagenesis through CRISPR-Cas of one of the putative VRS5 interactors, HvTB2 (also referred to as COMPOSITUM 1 and BRANCHED AND INDETERMINATE SPIKELET 1), resulted in plants that have lost their characteristic unbranched spike architecture. More specifically, hvtb2 mutants exhibited branches arising at the main spike, suggesting that HvTB2 acts as inhibitor of branching. Our protein-protein interaction studies of VRS5 resulted in the identification of HvTB2 as putative interactor of VRS5, another key regulator of spike architecture in barley. The study presented here provides a first step to underpin the protein-protein interactome of VRS5 and to identify other, yet unknown, key regulators of barley plant architecture.

Measuring the frequency and distribution of meiotic crossovers in homozygous barley inbred lines

We report a novel approach for establishing the number and position of CO events in individual homozygous inbred plants by combining low level EMS mutagenesis, speed breeding, whole genome shotgun sequencing and sliding window analysis of the induced molecular variant data. We demonstrate the approach by exploring CO frequency and distribution in self-fertilised progeny of the inbred barley cultivar Bowman and compare these observations to similar data obtained from a Bowman nearly isogenic line (BW230 Hvmlh3) containing a mutation in the DNA mismatch repair gene HvMLH3. We have previously shown that Hvmlh3 decreases both plant fertility and recombination by ~50%. We compare our results to those from previously published traditional genetic analysis of F3 families derived from multiple F2 lines containing WT or mutant alleles of HvMLH3, revealing a high level of correspondence between analyses. We discuss possible applications of the approach in streamlining the assessment of recombination in plant meiosis research.

Genotyping by Sequencing Advancements in Barley

Barley is considered an ideal crop to study cereal genetics due to its close relationship with wheat and diploid ancestral genome. It plays a crucial role in reducing risks to global food security posed by climate change. Genetic variations in the traits of interest in crops are vital for their improvement. DNA markers have been widely used to estimate these variations in populations. With the advancements in next-generation sequencing, breeders could access different types of genetic variations within different lines, with single-nucleotide polymorphisms (SNPs) being the most common type. However, genotyping barley with whole genome sequencing (WGS) is challenged by the higher cost and computational demand caused by the large genome size (5.5GB) and a high proportion of repetitive sequences (80%). Genotyping-by-sequencing (GBS) protocols based on restriction enzymes and target enrichment allow a cost-effective SNP discovery by reducing the genome complexity. In general, GBS has opened up new horizons for plant breeding and genetics. Though considered a reliable alternative to WGS, GBS also presents various computational difficulties, but GBS-specific pipelines are designed to overcome these challenges. Moreover, a robust design for GBS can facilitate the imputation to the WGS level of crops with high linkage disequilibrium. The complete exploitation of GBS advancements will pave the way to a better understanding of crop genetics and offer opportunities for the successful improvement of barley and its close relatives.

The barley mutant multiflorus2.b reveals quantitative genetic variation for new spikelet architecture

Spikelet indeterminacy and supernumerary spikelet phenotypes in barley multiflorus2.b mutant show polygenic inheritance. Genetic analysis of multiflorus2.b revealed major QTLs for spikelet determinacy and supernumerary spikelet phenotypes on 2H and 6H chromosomes. Understanding the genetic basis of yield forming factors in small grain cereals is of extreme importance, especially in the wake of stagnation of further yield gains in these crops. One such yield forming factor in these cereals is the number of grain-bearing florets produced per spikelet. Wild-type barley (Hordeum vulgare L.) spikelets are determinate structures, and the spikelet axis (rachilla) degenerates after producing single floret. In contrast, the rachilla of wheat (Triticum ssp.) spikelets, which are indeterminate, elongates to produce up to 12 florets. In our study, we characterized the barley spikelet determinacy mutant multiflorus2.b (mul2.b) that produced up to three fertile florets on elongated rachillae of lateral spikelets. Apart from the lateral spikelet indeterminacy (LS-IN), we also characterized the supernumerary spikelet phenotype in the central spikelets (CS-SS) of mul2.b. Through our phenotypic and genetic analyses, we identified two major QTLs on chromosomes 2H and 6H, and two minor QTLs on 3H for the LS-IN phenotype. For, the CS-SS phenotype, we identified one major QTL on 6H, and a minor QTL on 5H chromosomes. Notably, the 6H QTLs for CS-SS and LS-IN phenotypes co-located with each other, potentially indicating that a single genetic factor might regulate both phenotypes. Thus, our in-depth phenotyping combined with genetic analyses revealed the quantitative nature of the LS-IN and CS-SS phenotypes in mul2.b, paving the way for cloning the genes underlying these QTLs in the future.

Scripting Analyses of Genomes in Ensembl Plants

Ensembl Plants ( http://plants.ensembl.org ) offers genome-scale information for plants, with four releases per year. As of release 47 (April 2020) it features 79 species and includes genome sequence, gene models, and functional annotation. Comparative analyses help reconstruct the evolutionary history of gene families, genomes, and components of polyploid genomes. Some species have gene expression baseline reports or variation across genotypes. While the data can be accessed through the Ensembl genome browser, here we review specifically how our plant genomes can be interrogated programmatically and the data downloaded in bulk. These access routes are generally consistent across Ensembl for other non-plant species, including plant pathogens, pests, and pollinators.

Mutations in a barley cytochrome P450 gene enhances pathogen induced programmed cell death and cutin layer instability

Disease lesion mimic mutants (DLMMs) are characterized by the spontaneous development of necrotic spots with various phenotypes designated as necrotic (nec) mutants in barley. The nec mutants were traditionally considered to have aberrant regulation of programmed cell death (PCD) pathways, which have roles in plant immunity and development. Most barley nec3 mutants express cream to orange necrotic lesions contrasting them from typical spontaneous DLMMs that develop dark pigmented lesions indicative of serotonin/phenolics deposition. Barley nec3 mutants grown under sterile conditions did not exhibit necrotic phenotypes until inoculated with adapted pathogens, suggesting that they are not typical DLMMs. The F₂ progeny of a cross between nec3-γ1 and variety Quest segregated as a single recessive susceptibility gene post-inoculation with Bipolaris sorokiniana, the causal agent of the disease spot blotch. Nec3 was genetically delimited to 0.14 cM representing 16.5 megabases of physical sequence containing 149 annotated high confidence genes. RNAseq and comparative analysis of the wild type and five independent nec3 mutants identified a single candidate cytochrome P450 gene (HORVU.MOREX.r2.6HG0460850) that was validated as nec3 by independent mutations that result in predicted nonfunctional proteins. Histology studies determined that nec3 mutants had an unstable cutin layer that disrupted normal Bipolaris sorokiniana germ tube development.

At the site of pathogen infection, plant defense mechanisms rely on controlled programmed cell death (PCD) to sequester biotrophic pathogens that require living cells to extract nutrients from the host. However, these defense mechanisms are hijacked by necrotrophic plant pathogens that purposefully induce PCD to feed on the dead cells, thus facilitating further disease development. Thus, understanding PCD responses is important for resistance to both classes of pathogens. We characterized five independent disease lesion mimic mutants of barley designated necrotic 3 (nec3) that show aberrant regulation of PCD responses upon pathogen challenge. A cytochrome P450 gene was identified as Nec3 encoding a Tryptamine 5-Hydroxylase that functions as a terminal serotonin biosynthetic enzyme in the Tryptophan pathway of plants. We posit that nec3 mutants have disrupted serotonin biosynthesis resulting in expanded PCD, necrotrophic pathogen susceptibility and cutin layer instability. The nec3 mutants show expanded PCD and disease susceptibility of pathogen-induced necrotic lesions, suggesting a role of serotonin to sequester PCD and suppress pathogen colonization. The identification of Nec3 will facilitate functional analysis to elucidate the role that serotonin plays in the elicitation or suppression of PCD immunity responses to diverse pathogens and the effects it has on cutin layer biosynthesis.

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

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

The barley immune receptor Mla recognizes multiple pathogens and contributes to host range dynamics

Crop losses caused by plant pathogens are a primary threat to stable food production. Stripe rust (Puccinia striiformis) is a fungal pathogen of cereal crops that causes significant, persistent yield loss. Stripe rust exhibits host species specificity, with lineages that have adapted to infect wheat and barley. While wheat stripe rust and barley stripe rust are commonly restricted to their corresponding hosts, the genes underlying this host specificity remain unknown. Here, we show that three resistance genes, Rps6, Rps7, and Rps8, contribute to immunity in barley to wheat stripe rust. Rps7 cosegregates with barley powdery mildew resistance at the Mla locus. Using transgenic complementation of different Mla alleles, we confirm allele-specific recognition of wheat stripe rust by Mla. Our results show that major resistance genes contribute to the host species specificity of wheat stripe rust on barley and that a shared genetic architecture underlies resistance to the adapted pathogen barley powdery mildew and non-adapted pathogen wheat stripe rust.

Untargeted metabotyping to study phenylpropanoid diversity in crop plants

Plant genebanks constitute a key resource for breeding to ensure crop yield under changing environmental conditions. Because of their roles in a range of stress responses, phenylpropanoids are promising targets. Phenylpropanoids comprise a wide array of metabolites; however, studies regarding their diversity and the underlying genes are still limited for cereals. The assessment of barley diversity via genotyping-by-sequencing is in rapid progress. Exploring these resources by integrating genetic association studies to in-depth metabolomic profiling provides a valuable opportunity to study barley phenylpropanoid metabolism; but poses a challenge by demanding large-scale approaches. Here, we report an LC-PDA-MS workflow for barley high-throughput metabotyping. Without prior construction of a species-specific library, this method produced phenylpropanoid-enriched metabotypes with which the abundance of putative metabolic features was assessed across hundreds of samples in a single-processed data matrix. The robustness of the analytical performance was tested using a standard mix and extracts from two selected cultivars: Scarlett and Barke. The large-scale analysis of barley extracts showed (1) that barley flag leaf profiles were dominated by glycosylation derivatives of isovitexin, isoorientin, and isoscoparin; (2) proved the workflow's capability to discriminate within genotypes; (3) highlighted the role of glycosylation in barley phenylpropanoid diversity. Using the barley S42IL mapping population, the workflow proved useful for metabolic quantitative trait loci purposes. The protocol can be readily applied not only to explore the barley phenylpropanoid diversity represented in genebanks but also to study species whose profiles differ from those of cereals: the crop Helianthus annuus (sunflower) and the model plant Arabidopsis thaliana.

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

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

Identification of genomic loci conferring broad-spectrum resistance to multiple nematode species in exotic soybean accession PI 567305

KEY MESSAGE

Genetic analysis identified a unique combination of major QTL for resistance to important soybean nematodes concurrently present in a single soybean accession, which has not been reported earlier. An exotic soybean [Glycine max (L.) Merr.] accession, PI 567305, was reported to be highly resistant to three important nematode species, soybean cyst (SCN), root-knot (RKN), and reniform (RN) nematodes. However, genetic basis controlling broad-spectrum resistance in this germplasm has not been investigated. We report results of genetic analysis to identify genomic loci conferring resistance to these nematode species. A bi-parental population consisting of 242 F8-derived recombinant inbred lines (RILs) was developed from a cross of a nematode susceptible cultivar, Magellan, and resistant accession, PI 567305. The RILs were phenotyped for nematode resistance to three SCN HG types. They were genotyped using the Infinium SoySNP6K BeadChips and genotype-by-sequencing (GBS) methods in an attempt to evaluate the cost-effectiveness and efficiency of these two genotyping platforms. Genetic analysis confirmed the major QTL on chromosomes (Chrs) 10 and 18 with broad-spectrum resistance to the three nematodes present in this germplasm. Haplotype and copy number variation analyses of SCN resistance QTL indicated that PI 567305 has a different haplotype, which is associated with likely a unique SCN resistance mechanism different from Peking- or PI 88788-type resistance. The evaluations of both Infinium Beadchip- and GBS-based genotyping technologies provided comprehensive insights for researchers to choose a cost-effective and efficient platform for QTL mapping and for other genomic studies in soybeans.

Genome-Wide Association Study of Waterlogging Tolerance in Barley (Hordeum vulgare L.) Under Controlled Field Conditions

Waterlogging is one of the main abiotic stresses severely reducing barley grain yield. Barley breeding programs focusing on waterlogging tolerance require an understanding of genetic loci and alleles in the current germplasm. In this study, 247 worldwide spring barley genotypes grown under controlled field conditions were genotyped with 35,926 SNPs with minor allele frequency (MAF) > 0.05. Significant phenotypic variation in each trait, including biomass, spikes per plant, grains per plant, kernel weight per plant, plant height and chlorophyll content, was observed. A genome-wide association study (GWAS) based on linkage disequilibrium (LD) for waterlogging tolerance was conducted. Population structure analysis divided the population into three subgroups. A mixed linkage model using both population structure and kinship matrix (Q+K) was performed. We identified 17 genomic regions containing 51 significant waterlogging-tolerance-associated markers for waterlogging tolerance response, accounting for 5.8-11.5% of the phenotypic variation, with a majority of them localized on chromosomes 1H, 2H, 4H, and 5H. Six novel QTL were identified and eight potential candidate genes mediating responses to abiotic stresses were located at QTL associated with waterlogging tolerance. To our awareness, this is the first GWAS for waterlogging tolerance in a worldwide barley collection under controlled field conditions. The marker-trait associations could be used in the marker-assisted selection of waterlogging tolerance and will facilitate barley breeding.

Long-read sequence assembly: a technical evaluation in barley

Sequence assembly of large and repeat-rich plant genomes has been challenging, requiring substantial computational resources and often several complementary sequence assembly and genome mapping approaches. The recent development of fast and accurate long-read sequencing by circular consensus sequencing (CCS) on the PacBio platform may greatly increase the scope of plant pan-genome projects. Here, we compare current long-read sequencing platforms regarding their ability to rapidly generate contiguous sequence assemblies in pan-genome studies of barley (Hordeum vulgare). Most long-read assemblies are clearly superior to the current barley reference sequence based on short-reads. Assemblies derived from accurate long reads excel in most metrics, but the CCS approach was the most cost-effective strategy for assembling tens of barley genomes. A downsampling analysis indicated that 20-fold CCS coverage can yield very good sequence assemblies, while even five-fold CCS data may capture the complete sequence of most genes. We present an updated reference genome assembly for barley with near-complete representation of the repeat-rich intergenic space. Long-read assembly can underpin the construction of accurate and complete sequences of multiple genomes of a species to build pan-genome infrastructures in Triticeae crops and their wild relatives.

Building pan-genome infrastructures for crop plants and their use in association genetics

Pan-genomic studies aim at representing the entire sequence diversity within a species to provide useful resources for evolutionary studies, functional genomics and breeding of cultivated plants. Cost reductions in high-throughput sequencing and advances in sequence assembly algorithms have made it possible to create multiple reference genomes along with a catalogue of all forms of genetic variations in plant species with large and complex or polyploid genomes. In this review, we summarize the current approaches to building pan-genomes as an in silico representation of plant sequence diversity and outline relevant methods for their effective utilization in linking structural with phenotypic variation. We propose as future research avenues (i) transcriptomic and epigenomic studies across multiple reference genomes and (ii) the development of user-friendly and feature-rich pan-genome browsers.

Genomic region associated with pod color variation in pea (Pisum sativum)

Pea (Pisum sativum) was chosen as the research material by Gregor Mendel to discover the laws of inheritance. Out of seven traits studied by Mendel, genes controlling three traits including pod shape, pod color, and flower position have not been identified to date. With the aim of identifying the genomic region controlling pod color, we determined the genome sequence of a pea line with yellow pods. Genome sequence reads obtained using a Nanopore sequencing technology were assembled into 117,981 contigs (3.3 Gb), with an N50 value of 51.2 kb. A total of 531,242 potential protein-coding genes were predicted, of which 519,349 (2.8 Gb) were located within repetitive sequences (2.8 Gb). The assembled sequences were ordered using a reference as a guide to build pseudomolecules. Subsequent genetic and association analyses led to the identification of a genomic region that controls pea pod color. DNA sequences at this genomic location and transcriptome profiles of green and yellow pod lines were analyzed, and genes encoding 3' exoribonucleases were selected as potential candidates controlling pod color. The results presented in this study are expected to accelerate pan-genome studies in pea and facilitate the identification of the gene controlling one of the traits studied by Mendel.

Identification of QTL for barley grain size

Background

Barley grain size is one of the key factors determining storage capacity during grain filling. Large, well-filled grains also have a high malt extract potential. Grain size is a complex quantitative trait and can be easily affected by environmental factors thus the identification of genes controlling the trait and the use of molecular markers linked to the genes in breeding program is the most effective way of improving grain size.

Methods

Grain sizes of 188 doubled-haploid (DH) lines derived from the cross of a Japanese malting barley variety (Naso Nijo) and a Chinese feed barley variety (TX9425) were obtained from three different sites in two consecutive years. The average data were used for identifying QTL for grain size.

Results

A total of four significant QTL were identified for grain length (GL) and three for grain width (GW). The two major GL QTL are located at similar positions to the QTL for malt extract on 2H and uzu gene on 3H, respectively. However, the GL QTL on 2H is more likely a different one from the malt extract QTL as most of the candidate genes are located outside the fine mapped QTL region for malt extract. The GL QTL on 3H is closely linked with uzu gene but not due to a pleiotropic effect of uzu. The three QTL for grain width on 1H, 2H and 5H, respectively, were located at same position to those for GL.

Genome-wide association study of multiple traits linked to heat tolerance in emmer-derived hexaploid wheat genotypes

Heat stress tolerance in plants is a complex trait controlled by multiple genes of minor effect which are influenced by the environment and this makes breeding and selection complicated. Emmer wheat (Triticum dicoccon Schrank) carries valuable diversity that can be used to improve the heat tolerance of modern bread wheat. A diverse set of emmer-based genotypes was developed by crossing emmer wheat with hexaploid wheat. These materials, along with their hexaploid recurrent parents and commercial cultivars, were evaluated at optimum (E1) and heat stressed (E2) sowing times in the field for three consecutive years (2014-2016). The material was genotyped using the Infinium iSelect SNP 90K SNP Assay. The phenotypic data were combined across years within each sowing time and best linear unbiased estimators calculated for each genotype in each environment. These estimates were used for GWAS analysis. Significant phenotypic and genotypic variation was observed for all traits. A total of 125 and 142 marker-trait associations (MTAs) were identified in E1 and E2, respectively. The highest number of MTAs were observed on the A genome (106), followed by the B (105) and D (56) genomes. MTAs with pleiotropic effects within and across the environments were observed. Many of the MTAs found were reported previously for various traits, and a few significant MTAs under heat stress were new and linked to emmer genome. Genomic regions identified on chromosomes 2B and 3A had a significant positive impact on grain yield under stress with a 7% allelic effect. Genomic regions on chromosomes 1A and 4B contributed 11% and 9% of the variation for thousand kernel weight (TKW) under heat stress respectively. Following fine mapping, these regions could be used for marker-assisted selection to improve heat tolerance in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01222-3.

Genome-wide association mapping identifies HvNIP2;2/HvLsi6 accounting for efficient boron transport in barley

Boron (B) is an essential mineral element for plant growth, and the seed B pool of crops can be crucial when seedlings need to establish on low-B soils. To date, it is poorly understood how B accumulation in grain crops is genetically controlled. Here, we assessed the genotypic variation of the B concentration in grains of a spring barley (Hordeum vulgare L.) association panel that represents broad genetic diversity. We found a large genetic variation of the grain B concentration and detected in total 23 quantitative trait loci (QTLs) using genome-wide association mapping. HvNIP2;2/HvLsi6, encoding a potential B-transporting membrane protein, mapped closely to a major-effect QTL accounting for the largest proportion of grain B variation. Based on transport studies using heterologous expression systems and gene expression analysis, we demonstrate that HvNIP2;2/HvLsi6 represents a functional B channel and that expression variation in its transcript level associates with root and shoot B concentrations as well as with root dry mass formation under B-deficient conditions.

[NGS sequencing in barley breeding and genetic studies]

Barley (Hordeum vulgare L.) is the one of the most important cereal species used as food and feed crops, as well as for malting and alcohol production. At the end of the last century, traditional breeding techniques were complemented by the use of DNA markers. Molecular markers have also been used extensively for molecular genetic mapping and QTL analysis. In 2012, the barley genome sequencing was completed, which provided a broad range of new opportunities - from a more efficient search for candidate genes controlling economically important traits to genomic selection. The review summarizes the results of the studies performed after barley genome sequencing, which discovered new areas of barley genetics and breeding with high throughput screening and genotyping methods. During this period, intensive studies aimed at identification of barley genomic loci associated with economically important traits have been carried out; online databases and tools for working with barley genomic data and their deposition have appeared and are being replenished. In recent years, GWAS analysis has been used for large-scale phenotypegenotype association studies, which has been widely used in barley since 2010 due to the developed SNP-arrays, as well as genotyping methods based on direct NGS sequencing of selected fractions of the genome. To date, more than 80 papers have been published that describe the results of the GWAS analysis in barley. SNP identification associated with economically important traits and their transformation into CAPS or KASP markers convenient for screening selection material significantly expands the possibilities of marker-assisted selection of barley. In addition, the currently available information on potential target genes and the quality of the whole barley genome sequence provides a good base for applying genome editing technologies to create material for the creation of varieties with desired properties.

High-resolution mapping of Rym14Hb, a wild relative resistance gene to barley yellow mosaic disease

We mapped the Rym14^Hb resistance locus to barley yellow mosaic disease in a 2Mbp interval. The co-segregating markers will be instrumental for marker-assisted selection in barley breeding. Barley yellow mosaic disease is caused by Barley yellow mosaic virus and Barley mild mosaic virus and leads to severe yield losses in barley (Hordeum vulgare) in Central Europe and East-Asia. Several resistance loci are used in barley breeding. However, cases of resistance-breaking viral strains are known, raising concerns about the durability of those genes. Rym14^Hb is a dominant major resistance gene on chromosome 6HS, originating from barley's secondary genepool wild relative Hordeum bulbosum. As such, the resistance mechanism may represent a case of non-host resistance, which could enhance its durability. A susceptible barley variety and a resistant H. bulbosum introgression line were crossed to produce a large F₂ mapping population (n = 7500), to compensate for a ten-fold reduction in recombination rate compared to intraspecific barley crosses. After high-throughput genotyping, the Rym14^Hb locus was assigned to a 2Mbp telomeric interval on chromosome 6HS. The co-segregating markers developed in this study can be used for marker-assisted introgression of this locus into barley elite germplasm with a minimum of linkage drag.

The Barley HvWRKY6 Transcription Factor Is Required for Resistance Against Pyrenophora teres f. teres

Barley is an important cereal crop worldwide because of its use in the brewing and distilling industry. However, adequate supplies of quality malting barley are threatened by global climate change due to drought in some regions and excess precipitation in others, which facilitates epidemics caused by fungal pathogens. The disease net form net blotch caused by the necrotrophic fungal pathogen Pyrenophora teres f. teres (Ptt) has emerged as a global threat to barley production and diverse populations of Ptt have shown a capacity to overcome deployed genetic resistances. The barley line CI5791 exhibits remarkably effective resistance to diverse Ptt isolates from around the world that maps to two major QTL on chromosomes 3H and 6H. To identify genes involved in this effective resistance, CI5791 seed were γ-irradiated and two mutants, designated CI5791-γ3 and CI5791-γ8, with compromised Ptt resistance were identified from an M₂ population. Phenotyping of CI5791-γ3 and -γ8 × Heartland F₂ populations showed three resistant to one susceptible segregation ratios and CI5791-γ3 × -γ8 F₁ individuals were susceptible, thus these independent mutants are in a single allelic gene. Thirty-four homozygous mutant (susceptible) CI5791-γ3 × Heartland F₂ individuals, representing 68 recombinant gametes, were genotyped via PCR genotype by sequencing. The data were used for single marker regression mapping placing the mutation on chromosome 3H within an approximate 75 cM interval encompassing the 3H CI5791 resistance QTL. Sequencing of the mutants and wild-type (WT) CI5791 genomic DNA following exome capture identified independent mutations of the HvWRKY6 transcription factor located on chromosome 3H at ∼50.7 cM, within the genetically delimited region. Post transcriptional gene silencing of HvWRKY6 in barley line CI5791 resulted in Ptt susceptibility, confirming that it functions in NFNB resistance, validating it as the gene underlying the mutant phenotypes. Allele analysis and transcript regulation of HvWRKY6 from resistant and susceptible lines revealed sequence identity and upregulation upon pathogen challenge in all genotypes analyzed, suggesting a conserved transcription factor is involved in the defense against the necrotrophic pathogen. We hypothesize that HvWRKY6 functions as a conserved signaling component of defense mechanisms that restricts Ptt growth in barley.

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

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

Quantitative Trait Loci Mapping for Earliness, Fruit, and Seed Related Traits Using High Density Genotyping-by-Sequencing-Based Genetic Map in Bitter Gourd (Momordica charantia L.)

Bitter gourd (Momordica charantia L.) is an important vegetable crop having numerous medicinal properties. Earliness and yield related traits are main aims of bitter gourd breeding program. High resolution quantitative trait loci (QTLs) mapping can help in understanding the molecular basis of phenotypic variation of these traits and thus facilitate marker-assisted breeding. The aim of present study was to identify genetic loci controlling earliness, fruit, and seed related traits. To achieve this, genotyping-by-sequencing (GBS) approach was used to genotype 101 individuals of F₄ population derived from a cross between an elite cultivar Punjab-14 and PAUBG-6. This population was phenotyped under net-house conditions for three years 2018, 2019, and 2021. The linkage map consisting of 15 linkage groups comprising 3,144 single nucleotide polymorphism (SNP) markers was used to detect the QTLs for nine traits. A total of 50 QTLs for these traits were detected which were distributed on 11 chromosomes. The QTLs explained 5.09-29.82% of the phenotypic variance. The highest logarithm of the odds (LOD) score for a single QTL was 8.68 and the lowest was 2.50. For the earliness related traits, a total of 22 QTLs were detected. For the fruit related traits, a total of 16 QTLs and for seed related traits, a total of 12 QTLs were detected. Out of 50 QTLs, 20 QTLs were considered as frequent QTLs (FQ-QTLs). The information generated in this study is very useful in the future for fine-mapping and marker-assisted selection for these traits in bitter gourd improvement program.

LDscaff: LD-based scaffolding of de novo genome assemblies

Background

Genome assembly is fundamental for de novo genome analysis. Hybrid assembly, utilizing various sequencing technologies increases both contiguity and accuracy. While such approaches require extra costly sequencing efforts, the information provided millions of existed whole-genome sequencing data have not been fully utilized to resolve the task of scaffolding. Genetic recombination patterns in population data indicate non-random association among alleles at different loci, can provide physical distance signals to guide scaffolding.

Results

In this paper, we propose LDscaff for draft genome assembly incorporating linkage disequilibrium information in population data. We evaluated the performance of our method with both simulated data and real data. We simulated scaffolds by splitting the pig reference genome and reassembled them. Gaps between scaffolds were introduced ranging from 0 to 100 KB. The genome misassembly rate is 2.43% when there is no gap. Then we implemented our method to refine the Giant Panda genome and the donkey genome, which are purely assembled by NGS data. After LDscaff treatment, the resulting Panda assembly has scaffold N50 of 3.6 MB, 2.5 times larger than the original N50 (1.3 MB). The re-assembled donkey assembly has an improved N50 length of 32.1 MB from 23.8 MB.

Conclusions

Our method effectively improves the assemblies with existed re-sequencing data, and is an potential alternative to the existing assemblers required for the collection of new data.

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.

Assembly of whole-chromosome pseudomolecules for polyploid plant genomes using outbred mapping populations

Despite advances in sequencing technologies, assembly of complex plant genomes remains elusive due to polyploidy and high repeat content. Here we report PolyGembler for grouping and ordering contigs into pseudomolecules by genetic linkage analysis. Our approach also provides an accurate method with which to detect and fix assembly errors. Using simulated data, we demonstrate that our approach is of high accuracy and outperforms three existing state-of-the-art genetic mapping tools. Particularly, our approach is more robust to the presence of missing genotype data and genotyping errors. We used our method to construct pseudomolecules for allotetraploid lawn grass utilizing PacBio long reads in combination with restriction site-associated DNA sequencing, and for diploid Ipomoea trifida and autotetraploid potato utilizing contigs assembled from Illumina reads in combination with genotype data generated by single-nucleotide polymorphism arrays and genotyping by sequencing, respectively. We resolved 13 assembly errors for a published I. trifida genome assembly and anchored eight unplaced scaffolds in the published potato genome.

Estimating the contribution of Greenland Halibut (Reinhardtius hippoglossoides) stocks to nurseries by means of genotyping-by-sequencing: Sex and time matter

Identification of stocks and quantification of their relative contribution to recruitment are major objectives toward improving the management and conservation of marine exploited species. Next-generation sequencing allows for thousands of genomic markers to be analyzed, which provides the resolution needed to address these questions in marine species with weakly differentiated populations. Greenland Halibut (Reinhardtius hippoglossoides) is one of the most important exploited demersal species throughout the North Atlantic, in particular in the Gulf of St. Lawrence, Canada. There, two nurseries are known, the St. Lawrence Estuary and the northern Anticosti Island, but their contribution to the renewal of stocks remains unknown. The goals of this study were (a) to document the genetic structure and (b) to estimate the contribution of the different identified breeding stocks to nurseries. We sampled 100 juveniles per nursery and 50 adults from seven sites ranging from Saguenay Fjord to offshore Newfoundland, with some sites sampled over two consecutive years in order to evaluate the temporal stability of the contribution. Our results show that after removing sex-linked markers, the Estuary/Gulf of St. Lawrence represents a single stock which is genetically distinct from the Atlantic around Newfoundland (F ST = 0.00146, p-value = .001). Population assignment showed that recruitment in both nurseries is predominantly associated with the St. Lawrence stock. However, we found that the relative contribution of both stocks to the nurseries is temporally variable with 1% contribution of the Newfoundland stock one year but up to 33% for the second year, which may be caused by year-to-year variation in larval transport into the Gulf of St. Lawrence. This study serves as a model for the identification of stocks for fisheries resources in a context where few barriers to dispersal occur, in addition to demonstrating the importance of considering sex-linked markers and temporal replicates in studies of population genomics.

A high-density SNP-based linkage map using genotyping-by-sequencing and its utilization for improved genome assembly of chickpea (Cicer arietinum L.)

Genotyping-by-sequencing (GBS) allows rapid identification of markers for use in development of linkage maps, which expedite efficient breeding programs. In the present study, we have utilized GBS approach to identify and genotype single-nucleotide polymorphism (SNP) markers in an inter-specific RIL population of Cicer arietinum L. X C. reticulatum. A total of 141,639 raw SNPs were identified using the TASSEL-GBS pipeline. After stringent filtering, 8208 candidate SNPs were identified of which ~ 37% were localized in the intragenic regions followed by genic regions (~ 30%) and intergenic regions (~ 27%). We then utilized 6920 stringent selected SNPs from present study and 6714 SNPs and microsatellite markers available from previous studies for construction of linkage map. The resulting high-density linkage map comprising of eight linkage groups contained 13,590 markers which spanned 1299.14 cM of map length with an average marker density of 0.095 cM. Further, the derived linkage map was used to improve the available assembly of desi chickpea genome by anchoring 443 previously unplaced scaffolds onto eight linkage groups. The present efforts have refined anchoring of the desi chickpea genome assembly to 55.57% of the ~ 520 Mb of assembled desi genome. To the best of our knowledge, the linkage map generated in the present study represents one of the most dense linkage map developed for the crop till date. It will serve as a valuable resource for fine mapping and positional cloning of important quantitative trait loci (QTLs) associated with agronomical traits and also for anchoring and ordering of future genome sequence assemblies.

History and future perspectives of barley genomics

Barley (Hordeum vulgare), one of the most widely cultivated cereal crops, possesses a large genome of 5.1 Gbp. Through various international collaborations, the genome has recently been sequenced and assembled at the chromosome-scale by exploiting available genetic and genomic resources. Many wild and cultivated barley accessions have been collected and preserved around the world. These accessions are crucial to obtain diverse natural and induced barley variants. The barley bioresource project aims to investigate the diversity of this crop based on purified seed and DNA samples of a large number of collected accessions. The long-term goal of this project is to analyse the genome sequences of major barley accessions worldwide. In view of technical limitations, a strategy has been employed to establish the exome structure of a selected number of accessions and to perform high-quality chromosome-scale assembly of the genomes of several major representative accessions. For the future project, an efficient annotation pipeline is essential for establishing the function of genomes and genes as well as for using this information for sequence-based digital barley breeding. In this article, the author reviews the existing barley resources along with their applications and discuss possible future directions of research in barley genomics.

European maize genomes highlight intraspecies variation in repeat and gene content

The diversity of maize (Zea mays) is the backbone of modern heterotic patterns and hybrid breeding. Historically, US farmers exploited this variability to establish today’s highly productive Corn Belt inbred lines from blends of dent and flint germplasm pools. Here, we report de novo genome sequences of four European flint lines assembled to pseudomolecules with scaffold N50 ranging from 6.1 to 10.4 Mb. Comparative analyses with two US Corn Belt lines explains the pronounced differences between both germplasms. While overall syntenic order and consolidated gene annotations reveal only moderate pangenomic differences, whole-genome alignments delineating the core and dispensable genome, and the analysis of heterochromatic knobs and orthologous long terminal repeat retrotransposons unveil the dynamics of the maize genome. The high-quality genome sequences of the flint pool complement the maize pangenome and provide an important tool to study maize improvement at a genome scale and to enhance modern hybrid breeding.

De novo genome assemblies of four European flint maize lines and comparison with two US Corn Belt genomes provide insights into the dynamics of intraspecies variation in repeat and gene content in maize genomes.

The Battle to Sequence the Bread Wheat Genome: A Tale of the Three Kingdoms

In the year 2018, the world witnessed the finale of the race to sequence the genome of the world's most widely grown crop, the common wheat. Wheat has been known to bear a notoriously large and complicated genome of a polyploidy nature. A decade competition to sequence the wheat genome initiated with a single consortium of multiple countries, taking a conventional strategy similar to that for sequencing Arabidopsis and rice, became ferocious over time as both sequencing technologies and genome assembling methodologies advanced. At different stages, multiple versions of genome sequences of the same variety (e.g., Chinese Spring) were produced by several groups with their special strategies. Finally, 16 years after the rice genome was finished and 9 years after that of maize, the wheat research community now possesses its own reference genome. Armed with these genomics tools, wheat will reestablish itself as a model for polyploid plants in studying the mechanisms of polyploidy evolution, domestication, genetic and epigenetic regulation of homoeolog expression, as well as defining its genetic diversity and breeding on the genome level. The enhanced resolution of the wheat genome should also help accelerate development of wheat cultivars that are more tolerant to biotic and/or abiotic stresses with better quality and higher yield.

An improved high-quality genome assembly and annotation of Tibetan hulless barley

Hulless barley (Hordeum vulgare L. var. nudum) is a barley variety that has loose husk cover of the caryopses. Because of the ease in processing and edibility, hulless barley has been locally cultivated and used as human food. For example, in Tibetan Plateau, hulless barley is the staple food for human and essential livestock feed. Although the draft genome of hulless barley has been sequenced, the assembly remains fragmented. Here, we reported an improved high-quality assembly and annotation of the Tibetan hulless barley genome using more than 67X PacBio long-reads. The N50 contig length of the new assembly is at least more than 19 times larger than other available barley assemblies. The new genome assembly also showed high gene completeness and high collinearity of genome synteny with the previously reported barley genome. The new genome assembly and annotation will not only remove major hurdles in genetic analysis and breeding of hulless barley, but will also serve as a key resource for studying barley genomics and genetics.

Genome-Wide Mapping of Adult Plant Resistance to Leaf Rust and Stripe Rust in CIMMYT Wheat Line Arableu#1

Leaf (brown) rust (LR) and stripe (yellow) rust (YR), caused by Puccinia triticina and P. striiformis f. sp. tritici, respectively, significantly reduce wheat production worldwide. Disease-resistant wheat varieties offer farmers one of the most effective ways to manage these diseases. The common wheat (Triticum aestivum L.) Arableu#1, developed by the International Maize and Wheat Improvement Center and released as Deka in Ethiopia, shows susceptibility to both LR and YR at the seedling stage but a high level of adult plant resistance (APR) to the diseases in the field. We used 142 F₅ recombinant inbred lines (RILs) derived from Apav#1 × Arableu#1 to identify quantitative trait loci (QTLs) for APR to LR and YR. A total of 4,298 genotyping-by-sequencing markers were used to construct a genetic linkage map. The study identified four LR resistance QTLs and six YR resistance QTLs in the population. Among these, QLr.cim-1BL.1/QYr.cim-1BL.1 was located in the same location as Lr46/Yr29, a known pleiotropic resistance gene. QLr.cim-1BL.2 and QYr.cim-1BL.2 were also located on wheat chromosome 1BL at 37 cM from Lr46/Yr29 and may represent a new segment for pleiotropic resistance to both rusts. QLr.cim-7BL is likely Lr68 given its association with the tightly linked molecular marker cs7BLNLRR. In addition, QLr.cim-3DS, QYr.cim-2AL, QYr.cim-4BL, QYr.cim-5AL, and QYr.cim-7DS are probably new resistance loci based on comparisons with published QTLs for resistance to LR and YR. Our results showed the diversity of minor resistance QTLs in Arableu#1 and their role in conferring near-immune levels of APR to both LR and YR, when combined with the pleiotropic APR gene Lr46/Yr29.

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

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

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

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

High Resolution Mapping of a Hordeum bulbosum-Derived Powdery Mildew Resistance Locus in Barley Using Distinct Homologous Introgression Lines

Powdery mildew caused by Blumeria graminis f. sp. hordei (Bgh) is one of the main foliar diseases in barley (Hordeum vulgare L.; Hv). Naturally occurring resistance genes used in barley breeding are a cost effective and environmentally sustainable strategy to minimize the impact of pathogens, however, the primary gene pool of H. vulgare contains limited diversity owing to recent domestication bottlenecks. To ensure durable resistance against this pathogen, more genes are required that could be unraveled by investigation of secondary barley gene-pool. A large set of Hordeum bulbosum (Hb) introgression lines (ILs) harboring a diverse set of desirable resistance traits have been developed and are being routinely used as source of novel diversity in gene mapping studies. Nevertheless, this strategy is often compromised by a lack of recombination between the introgressed fragment and the orthologous chromosome of the barley genome. In this study, we fine-mapped a Hb gene conferring resistance to barley powdery mildew. The initial genotyping of two Hb ILs mapping populations with differently sized 2HS introgressions revealed severely reduced interspecific recombination in the region of the introgressed segment, preventing precise localization of the gene. To overcome this difficulty, we developed an alternative strategy, exploiting intraspecific recombination by crossing two Hv/Hb ILs with collinear Hb introgressions, one of which carries a powdery mildew resistance gene, while the other doesn't. The intraspecific recombination rate in the Hb-introgressed fragment of 2HS was approximately 20 times higher than it was in the initial simple ILs mapping populations. Using high-throughput genotyping-by-sequencing (GBS), we allocated the resistance gene to a 1.4 Mb interval, based on an estimate using the Hv genome as reference, in populations of only 103 and 146 individuals, respectively, similar to what is expected at this locus in barley. The most likely candidate resistance gene within this interval is part of the coiled-coil nucleotide-binding-site leucine-rich-repeat (CC-NBS-LLR) gene family, which is over-represented among genes conferring strong dominant resistance to pathogens. The reported strategy can be applied as a general strategic approach for identifying genes underlying traits of interest in crop wild relatives.

Identification of quantitative trait loci for net form net blotch resistance in contemporary barley breeding germplasm from the USA using genome-wide association mapping

Association mapping study conducted in a population of 3490 elite barley breeding lines from ten barley breeding programs of the USA identified 12 QTLs for resistance/susceptibility to net form of net blotch. Breeding resistant varieties is the best management strategy for net form of net blotch (NFNB) in barley (Hordeum vulgare L.) caused by Pyrenophora teres f. teres (Ptt). Several resistance QTL have been previously identified in barley via linkage mapping and genome-wide association studies (GWAS). A GWAS conducted in a collection of advanced breeding lines (n = 3490) representing elite germplasm from ten barley breeding programs of the USA identified 42 unique marker-trait associations (MTA) for NFNB resistance. The lines were genotyped with 3072 SNP markers and phenotyped with four Ptt isolates in controlled environment. The lines were used to construct 13 different GWAS panels. Efficient mixed model association method with principal components and kinship was used for GWAS. Significance threshold for MTA was set at a false discovery rate of 0.05. Two, eight, six, one and 25 MTA were identified in chromosomes 1H, 3H, 4H, 5H and 6H, respectively. Based on genetic positions and linkage disequilibrium, these MTA's correspond to two, three, two, one and four QTLs in chromosome 1H, 3H, 4H, 5H and 6H, respectively. A comparison with previous linkage and GWAS studies revealed several previously identified and novel QTLs. Moreover, different genomic regions were found to be responsible for NFNB resistance in two-row versus six-row germplasm. The germplasm-specific SNP markers with additive effects and allelic distribution is reported to facilitate breeders in selection of markers for MAS to introgress novel net blotch resistance.

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

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

GWAS: Fast-forwarding gene identification and characterization in temperate Cereals: lessons from Barley - A review

Understanding the genetic complexity of traits is an important objective of small grain temperate cereals yield and adaptation improvements. Bi-parental quantitative trait loci (QTL) linkage mapping is a powerful method to identify genetic regions that co-segregate in the trait of interest within the research population. However, recently, association or linkage disequilibrium (LD) mapping using a genome-wide association study (GWAS) became an approach for unraveling the molecular genetic basis underlying the natural phenotypic variation. Many causative allele(s)/loci have been identified using the power of this approach which had not been detected in QTL mapping populations. In barley (Hordeum vulgare L.), GWAS has been successfully applied to define the causative allele(s)/loci which can be used in the breeding crop for adaptation and yield improvement. This promising approach represents a tremendous step forward in genetic analysis and undoubtedly proved it is a valuable tool in the identification of candidate genes. In this review, we describe the recently used approach for genetic analyses (linkage mapping or association mapping), and then provide the basic genetic and statistical concepts of GWAS, and subsequently highlight the genetic discoveries using GWAS. The review explained how the candidate gene(s) can be detected using state-of-art bioinformatic tools.

Effect of epistasis and environment on flowering time in barley reveals a novel flowering-delaying QTL allele

Flowering time is a complex trait and has a key role in crop yield and adaptation to environmental stressors such as heat and drought. This study aimed to better understand the interconnected dynamics of epistasis and environment and look for novel regulators. We investigated 534 spring barley MAGIC DH lines for flowering time at various environments. Analysis of quantitative trait loci (QTLs), epistatic interactions, QTL × environment (Q×E) interactions, and epistasis × environment (E×E) interactions were performed with single SNP and haplotype approaches. In total, 18 QTLs and 2420 epistatic interactions were detected, including intervals harboring major genes such as Ppd-H1, Vrn-H1, Vrn-H3, and denso/sdw1. Epistatic interactions found in field and semi-controlled conditions were distinctive. Q×E and E×E interactions revealed that temperature influenced flowering time by triggering different interactions between known and newly detected regulators. A novel flowering-delaying QTL allele was identified on chromosome 1H (named 'HvHeading') and was shown to be engaged in epistatic and environment interactions. Results suggest that investigating epistasis, environment, and their interactions, rather than only single QTLs, is an effective approach for detecting novel regulators. We assume that barley can adapt flowering time to the environment via alternative routes within the pathway.

Analysis of early-flowering genes at barley chromosome 2H expands the repertoire of mutant alleles at the Mat-c locus

Analyses of barley mat-c loss of function mutants reveal deletions, splice-site mutations and nonsynonymous substitutions in a key gene regulating early flowering. Optimal timing of flowering is critical for reproductive success and crop yield improvement. Several major quantitative trait loci for flowering time variation have been identified in barley. In the present study, we analyzed two near-isogenic lines, BW507 and BW508, which were reported to carry two independent early-flowering mutant loci, mat-b.7 and mat-c.19, respectively. Both introgression segments are co-localized in the pericentromeric region of chromosome 2H. We mapped the mutation in BW507 to a 31 Mbp interval on chromosome 2HL and concluded that BW507 has a deletion of Mat-c, which is an ortholog of Antirrhinum majus CENTRORADIALIS (AmCEN) and Arabidopsis thaliana TERMINAL FLOWER1 (AtTFL1). Contrary to previous reports, our data showed that both BW507 and BW508 are Mat-c deficient and none of them are mat-b.7 derived. This work complements previous studies by identifying the uncharacterized mat-c.19 mutant and seven additional mat-c mutants. Moreover, we explored the X-ray structure of AtTFL1 for prediction of the functional effects of nonsynonymous substitutions caused by mutations in Mat-c.

TRITEX: chromosome-scale sequence assembly of Triticeae genomes with open-source tools

Chromosome-scale genome sequence assemblies underpin pan-genomic studies. Recent genome assembly efforts in the large-genome Triticeae crops wheat and barley have relied on the commercial closed-source assembly algorithm DeNovoMagic. We present TRITEX, an open-source computational workflow that combines paired-end, mate-pair, 10X Genomics linked-read with chromosome conformation capture sequencing data to construct sequence scaffolds with megabase-scale contiguity ordered into chromosomal pseudomolecules. We evaluate the performance of TRITEX on publicly available sequence data of tetraploid wild emmer and hexaploid bread wheat, and construct an improved annotated reference genome sequence assembly of the barley cultivar Morex as a community resource.