A view of the pan-genome of domesticated Cowpea (Vigna unguiculata [L.] Walp.)

Cowpea, Vigna unguiculata L. Walp., is a diploid warm-season legume of critical importance as both food and fodder in sub-Saharan Africa. This species is also grown in Northern Africa, Europe, Latin America, North America, and East to Southeast Asia. To capture the genomic diversity of domesticates of this important legume, de novo genome assemblies were produced for representatives of six subpopulations of cultivated cowpea identified previously from genotyping of several hundred diverse accessions. In the most complete assembly (IT97K-499-35), 26,026 core and 4963 noncore genes were identified, with 35,436 pan genes when considering all seven accessions. GO terms associated with response to stress and defense response were highly enriched among the noncore genes, while core genes were enriched in terms related to transcription factor activity, and transport and metabolic processes. Over 5 million single nucleotide polymorphisms (SNPs) relative to each assembly and over 40 structural variants >1 Mb in size were identified by comparing genomes. Vu10 was the chromosome with the highest frequency of SNPs, and Vu04 had the most structural variants. Noncore genes harbor a larger proportion of potentially disruptive variants than core genes, including missense, stop gain, and frameshift mutations; this suggests that noncore genes substantially contribute to diversity within domesticated cowpea.

Urine and Dried Blood Spots From Children and Pregnant Women Reveal Phytochemicals, Amino Acids, and Carnitine Metabolites as Cowpea Consumption Biomarkers

SCOPE

Legumes consumption has been proven to promote health across the lifespan; cowpeas have demonstrated efficacy in combating childhood malnutrition and growth faltering, with an estimated malnutrition prevalence of 35.6% of children in Ghana. This cowpea feeding study aimed to identify a suite of metabolic consumption biomarkers in children and adults.

METHODS AND RESULTS

Urine and dried blood spots (DBS) from 24 children (9-21 months) and 21 pregnant women (>18 years) in Northern Ghana are collected before and after dose-escalated consumption of four cowpea varieties for 15 days. Untargeted metabolomics identified significant increases in amino acids, phytochemicals, and lipids. The carnitine metabolism pathway is represented by 137 urine and 43 DBS metabolites, with significant changes to tiglylcarnitine and acetylcarnitine. Additional noteworthy candidate biomarkers are mansouramycin C, N-acetylalliin, proline betaine, N2, N5-diacetylornithine, S-methylcysteine, S-methylcysteine sulfoxide, and cis-urocanate. S-methylcysteine and S-methylcysteine sulfoxide are targeted and quantified in urine.

CONCLUSION

This feeding study for cowpea biomarkers supports the utility of a suite of key metabolites classified as amino acids, lipids, and phytochemicals for dietary legume and cowpea-specific food exposures of global health importance.

Corrigendum: Seed coat pattern QTL and development in cowpea (Vigna unguiculata [L.] Walp.)

[This corrects the article DOI: 10.3389/fpls.2019.01346.].

Non-targeted metabolomics of cooked cowpea (Vigna unguiculata) and pigeon pea (Cajanus cajan) from Ghana using two distinct and complementary analytical platforms

•

Cowpea varieties represent diverse staple foods in Sub-Saharan Africa.

•

Cowpea metabolite composition is different from pigeon pea and common bean.

•

Cowpea metabolites included tonkinelin, pheophytin A, and linoleoyl ethanolamide.

•

Pheophytin A was only detected in the cowpea variety Sangyi.

•

Pipecolic acid identification was confirmed for all three legumes.

Legumes are global staple foods with multiple human health properties that merit detailed composition analysis in cooked forms. This study analyzed cowpea [Vigna unguiculata] (three varieties: Dagbantuya, Sangyi, and Tukara), pigeon pea [Cajanus cajan], and common bean [Phaseolus vulgaris] using two distinct ultra-performance liquid chromatography mass spectrometry (UPLC-MS) platforms and analytical workflows. Comparisons between cowpea and pigeon pea consumed in Ghana, and common bean (navy bean) from USA, revealed 75 metabolites that differentiated cowpeas. Metabolite fold-change comparisons resulted in 142 metabolites with significantly higher abundance in cowpea, and 154 higher in abundance from pigeon pea. 3-(all-trans-nonaprenyl)benzene-1,2-diol, N-tetracosanoylphytosphingosine, and sitoindoside II are novel identifications in cowpea, with notably higher abundance than other legumes tested. Cowpea variety specific markers were tonkinelin (Dagbantuya), pheophytin A (Sangyi), and linoleoyl ethanolamide (Tukara). This study identified novel cowpea and pigeon pea food metabolites that warrant continued investigation as bioactive food components following consumption in people.

Genetic, anatomical, and environmental patterns related to pod shattering resistance in domesticated cowpea [Vigna unguiculata (L.) Walp]

Pod shattering, which causes the explosive release of seeds from the pod, is one of the main sources of yield losses in cowpea in arid and semi-arid areas. Reduction of shattering has therefore been a primary target for selection during domestication and improvement of cowpea, among other species. Using a mini-core diversity panel of 368 cowpea accessions, four regions with a statistically significant association with pod shattering were identified. Two genes (Vigun03g321100 and Vigun11g100600), involved in cell wall biosynthesis, were identified as strong candidates for pod shattering. Microscopic analysis was conducted on a subset of accessions representing the full spectrum of shattering phenotypes. This analysis indicated that the extent of wall fiber deposition was highly correlated with shattering. The results from this study also demonstrate that pod shattering in cowpea is exacerbated by arid environmental conditions. Finally, using a subset of West African landraces, patterns of historical selection for shattering resistance related to precipitation in the environment of origin were identified. Together, these results shed light on sources of resistance to pod shattering, which will, in turn, improve climate resilience of a major global nutritional staple.

Genome-Wide Association Study Reveals Candidate Genes for Flowering Time in Cowpea (Vigna unguiculata [L.] Walp.)

Cowpea (Vigna unguiculata [L.] Walp., diploid, 2n = 22) is a major crop used as a protein source for human consumption as well as a quality feed for livestock. It is drought and heat tolerant and has been bred to develop varieties that are resilient to changing climates. Plant adaptation to new climates and their yield are strongly affected by flowering time. Therefore, understanding the genetic basis of flowering time is critical to advance cowpea breeding. The aim of this study was to perform genome-wide association studies (GWAS) to identify marker trait associations for flowering time in cowpea using single nucleotide polymorphism (SNP) markers. A total of 368 accessions from a cowpea mini-core collection were evaluated in Ft. Collins, CO in 2019 and 2020, and 292 accessions were evaluated in Citra, FL in 2018. These accessions were genotyped using the Cowpea iSelect Consortium Array that contained 51,128 SNPs. GWAS revealed seven reliable SNPs for flowering time that explained 8-12% of the phenotypic variance. Candidate genes including FT, GI, CRY2, LSH3, UGT87A2, LIF2, and HTA9 that are associated with flowering time were identified for the significant SNP markers. Further efforts to validate these loci will help to understand their role in flowering time in cowpea, and it could facilitate the transfer of some of this knowledge to other closely related legume species.

BAC- and oligo-FISH mapping reveals chromosome evolution among Vigna angularis, V. unguiculata, and Phaseolus vulgaris

Cytogenomic resources have accelerated synteny and chromosome evolution studies in plant species, including legumes. Here, we established the first cytogenetic map of V. angularis (Va, subgenus Ceratotropis) and compared this new map with those of V. unguiculata (Vu, subgenus Vigna) and P. vulgaris (Pv) by BAC-FISH and oligopainting approaches. We mapped 19 Vu BACs and 35S rDNA probes to the 11 chromosome pairs of Va, Vu, and Pv. Vigna angularis shared a high degree of macrosynteny with Vu and Pv, with five conserved syntenic chromosomes. Additionally, we developed two oligo probes (Pv2 and Pv3) used to paint Vigna orthologous chromosomes. We confirmed two reciprocal translocations (chromosomes 2 and 3 and 1 and 8) that have occurred after the Vigna and Phaseolus divergence (~9.7 Mya). Besides, two inversions (2 and 4) and one translocation (1 and 5) have occurred after Vigna and Ceratotropis subgenera separation (~3.6 Mya). We also observed distinct oligopainting patterns for chromosomes 2 and 3 of Vigna species. Both Vigna species shared similar major rearrangements compared to Pv: one translocation (2 and 3) and one inversion (chromosome 3). The sequence synteny identified additional inversions and/or intrachromosomal translocations involving pericentromeric regions of both orthologous chromosomes. We propose chromosomes 2 and 3 as hotspots for chromosomal rearrangements and de novo centromere formation within and between Vigna and Phaseolus. Our BAC- and oligo-FISH mapping contributed to physically trace the chromosome evolution of Vigna and Phaseolus and its application in further studies of both genera.

Population Structure and Genetic Diversity in Korean Cowpea Germplasm Based on SNP Markers

Cowpea is one of the most essential legume crops providing inexpensive dietary protein and nutrients. The aim of this study was to understand the genetic diversity and population structure of global and Korean cowpea germplasms. A total of 384 cowpea accessions from 21 countries were genotyped with the Cowpea iSelect Consortium Array containing 51,128 single-nucleotide polymorphisms (SNPs). After SNP filtering, a genetic diversity study was carried out using 35,116 SNPs within 376 cowpea accessions, including 229 Korean accessions. Based on structure and principal component analysis, a total of 376 global accessions were divided into four major populations. Accessions in group 1 were from Asia and Europe, those in groups 2 and 4 were from Korea, and those in group 3 were from West Africa. In addition, 229 Korean accessions were divided into three major populations (Q1, Jeonra province; Q2, Gangwon province; Q3, a mixture of provinces). Additionally, the neighbor-joining tree indicated similar results. Further genetic diversity analysis within the global and Korean population groups indicated low heterozygosity, a low polymorphism information content, and a high inbreeding coefficient in the Korean cowpea accessions. The population structure analysis will provide useful knowledge to support the genetic potential of the cowpea breeding program, especially in Korea.

Identification of QTL for perenniality and floral scent in cowpea (Vigna unguiculata [L.] Walp.)

Perennial habit and floral scent are major traits that distinguish domesticated cowpeas from their wild relatives. However, the genetic basis of these two important traits remains largely unknown in cowpea. Plant longevity, a perenniality-related trait, and floral scent, an outcrossing trait, were investigated using a RIL population derived from a cross between a domesticated and a wild cowpea. QTL analysis revealed three significant loci, one on chromosome 8 associated with plant longevity and two, on chromosomes 1 and 11, for floral scent. Genes within the QTL regions were identified. Genes encoding an F-box protein (Vigun08g215300) and two kinases (Vigun08g217000, Vigun08g217800), and involved in physiological processes including regulation of flowering time and plant longevity, were identified within the perenniality QTL region. A cluster of O-methyltransferase genes (Vigun11g096800, Vigun11g096900, Vigun11g097000, Vigun11g097600, and Vigun11g097800) was identified within the floral scent QTL region. These O-methyltransferase cowpea genes are orthologs of the Arabidopsis N-acetylserotonin O-methyltransferase (ASMT) gene, which is involved in the biosynthesis of melatonin. Melatonin is an indole derivative, which is an essential molecule for plant interactions with pollinators. These findings lay the foundation for further exploration of the genetic mechanisms of perenniality and floral scent in cowpea. Knowledge from this study can help in the development of new extended-growth cycle lines with increased yield or lines with increased outcrossing for population breeding.

Perspectives on Low Temperature Tolerance and Vernalization Sensitivity in Barley: Prospects for Facultative Growth Habit

One option to achieving greater resiliency for barley production in the face of climate change is to explore the potential of winter and facultative growth habits: for both types, low temperature tolerance (LTT) and vernalization sensitivity are key traits. Sensitivity to short-day photoperiod is a desirable attribute for facultative types. In order to broaden our understanding of the genetics of these phenotypes, we mapped quantitative trait loci (QTLs) and identified candidate genes using a genome-wide association studies (GWAS) panel composed of 882 barley accessions that was genotyped with the Illumina 9K single-nucleotide polymorphism (SNP) chip. Fifteen loci including 5 known and 10 novel QTL/genes were identified for LTT-assessed as winter survival in 10 field tests and mapped using a GWAS meta-analysis. FR-H1, FR-H2, and FR-H3 were major drivers of LTT, and candidate genes were identified for FR-H3. The principal determinants of vernalization sensitivity were VRN-H1, VRN-H2, and PPD-H1. VRN-H2 deletions conferred insensitive or intermediate sensitivity to vernalization. A subset of accessions with maximum LTT were identified as a resource for allele mining and further characterization. Facultative types comprised a small portion of the GWAS panel but may be useful for developing germplasm with this growth habit.

A genome-wide association and meta-analysis reveal regions associated with seed size in cowpea [Vigna unguiculata (L.) Walp]

This paper combined GWAS, meta-analysis and sequence homology comparison with common bean to identify regions associated with seed size variation in domesticated cowpea. Seed size is an important trait for yield and commercial value in dry-grain cowpea. Seed size varies widely among different cowpea accessions, and the genetic basis of such variation is not yet well understood. To better decipher the genetic basis of seed size, a genome-wide association study (GWAS) and meta-analysis were conducted on a panel of 368 cowpea diverse accessions from 51 countries. Four traits, including seed weight, length, width and density were evaluated across three locations. Using 51,128 single nucleotide polymorphisms covering the cowpea genome, 17 loci were identified for these traits. One locus was common to weight, width and length, suggesting pleiotropy. By integrating synteny-based analysis with common bean, six candidate genes (Vigun05g036000, Vigun05g039600, Vigun05g204200, Vigun08g217000, Vigun11g187000, and Vigun11g191300) which are implicated in multiple functional categories related to seed size such as endosperm development, embryo development, and cell elongation were identified. These results suggest that a combination of GWAS meta-analysis with synteny comparison in a related plant is an efficient approach to identify candidate gene (s) for complex traits in cowpea. The identified loci and candidate genes provide useful information for improving cowpea varieties and for molecular investigation of seed size.

Seed Coat Pattern QTL and Development in Cowpea (Vigna unguiculata [L.] Walp.)

The appearance of the seed is an important aspect of consumer preference for cowpea (Vigna unguiculata [L.] Walp.). Seed coat pattern in cowpea has been a subject of study for over a century. This study makes use of newly available resources, including mapping populations, a reference genome and additional genome assemblies, and a high-density single nucleotide polymorphism genotyping platform, to map various seed coat pattern traits to three loci, concurrent with the Color Factor (C), Watson (W), and Holstein (H) factors identified previously. Several gene models encoding proteins involved in regulating the later stages of the flavonoid biosynthesis pathway have been identified as candidate genes, including a basic helix-loop-helix gene (Vigun07g110700) for the C locus, a WD-repeat gene (Vigun09g139900) for the W locus and an E3 ubiquitin ligase gene (Vigun10g163900) for the H locus. A model of seed coat development, consisting of six distinct stages, is described to explain some of the observed pattern phenotypes.

Identification of QTL controlling domestication-related traits in cowpea (Vigna unguiculata L. Walp)

Cowpea (Vigna unguiculata L. Walp) is a warm-season legume with a genetically diverse gene-pool composed of wild and cultivated forms. Cowpea domestication involved considerable phenotypic changes from the wild progenitor, including reduction of pod shattering, increased organ size, and changes in flowering time. Little is known about the genetic basis underlying these changes. In this study, 215 recombinant inbred lines derived from a cross between a cultivated and a wild cowpea accession were used to evaluate nine domestication-related traits (pod shattering, peduncle length, flower color, days to flowering, 100-seed weight, pod length, leaf length, leaf width and seed number per pod). A high-density genetic map containing 17,739 single nucleotide polymorphisms was constructed and used to identify 16 quantitative trait loci (QTL) for these nine traits. Based on annotations of the cowpea reference genome, genes within these regions are reported. Four regions with clusters of QTL were identified, including one on chromosome 8 related to increased organ size. This study provides new knowledge of the genomic regions controlling domestication-related traits in cowpea as well as candidate genes underlying those QTL. This information can help to exploit wild relatives in cowpea breeding programs.

A multi-parent advanced generation inter-cross (MAGIC) population for genetic analysis and improvement of cowpea (Vigna unguiculata L. Walp.)

Multi-parent advanced generation inter-cross (MAGIC) populations are an emerging type of resource for dissecting the genetic structure of traits and improving breeding populations. We developed a MAGIC population for cowpea (Vigna unguiculata L. Walp.) from eight founder parents. These founders were genetically diverse and carried many abiotic and biotic stress resistance, seed quality and agronomic traits relevant to cowpea improvement in the United States and sub-Saharan Africa, where cowpea is vitally important in the human diet and local economies. The eight parents were inter-crossed using structured matings to ensure that the population would have balanced representation from each parent, followed by single-seed descent, resulting in 305 F₈ recombinant inbred lines each carrying a mosaic of genome blocks contributed by all founders. This was confirmed by single nucleotide polymorphism genotyping with the Illumina Cowpea Consortium Array. These lines were on average 99.74% homozygous but also diverse in agronomic traits across environments. Quantitative trait loci (QTLs) were identified for several parental traits. Loci with major effects on photoperiod sensitivity and seed size were also verified by biparental genetic mapping. The recombination events were concentrated in telomeric regions. Due to its broad genetic base, this cowpea MAGIC population promises breakthroughs in genetic gain, QTL and gene discovery, enhancement of breeding populations and, for some lines, direct releases as new varieties.

Genetic diversity and structure of Iberian Peninsula cowpeas compared to world-wide cowpea accessions using high density SNP markers

Background

Cowpea (Vigna unguiculata L. Walp) is an important legume crop due to its high protein content, adaptation to heat and drought and capacity to fix nitrogen. Europe has a deficit of cowpea production. Knowledge of genetic diversity among cowpea landraces is important for the preservation of local varieties and is the basis to obtain improved varieties. The aims of this study were to explore diversity and the genetic structure of a set of Iberian Peninsula cowpea accessions in comparison to a worldwide collection and to infer possible dispersion routes of cultivated cowpea.

Results

The Illumina Cowpea iSelect Consortium Array containing 51,128 SNPs was used to genotype 96 cowpea accessions including 43 landraces and cultivars from the Iberian Peninsula, and 53 landraces collected worldwide. Four subpopulations were identified. Most Iberian Peninsula accessions clustered together with those from other southern European and northern African countries. Only one accession belonged to another subpopulation, while two accessions were ‘admixed’. A lower genetic diversity level was found in the Iberian Peninsula accessions compared to worldwide cowpeas.

Conclusions

The genetic analyses performed in this study brought some insights into worldwide genetic diversity and structure and possible dispersion routes of cultivated cowpea. Also, it provided an in-depth analysis of genetic diversity in Iberian Peninsula cowpeas that will help guide crossing strategies in breeding programs.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4295-0) contains supplementary material, which is available to authorized users.

Construction of a map-based reference genome sequence for barley, Hordeum vulgare L

Barley (Hordeum vulgare L.) is a cereal grass mainly used as animal fodder and raw material for the malting industry. The map-based reference genome sequence of barley cv. ‘Morex’ was constructed by the International Barley Genome Sequencing Consortium (IBSC) using hierarchical shotgun sequencing. Here, we report the experimental and computational procedures to (i) sequence and assemble more than 80,000 bacterial artificial chromosome (BAC) clones along the minimum tiling path of a genome-wide physical map, (ii) find and validate overlaps between adjacent BACs, (iii) construct 4,265 non-redundant sequence scaffolds representing clusters of overlapping BACs, and (iv) order and orient these BAC clusters along the seven barley chromosomes using positional information provided by dense genetic maps, an optical map and chromosome conformation capture sequencing (Hi-C). Integrative access to these sequence and mapping resources is provided by the barley genome explorer (BARLEX).

Genome resources for climate-resilient cowpea, an essential crop for food security

Cowpea (Vigna unguiculata L. Walp.) is a legume crop that is resilient to hot and drought-prone climates, and a primary source of protein in sub-Saharan Africa and other parts of the developing world. However, genome resources for cowpea have lagged behind most other major crops. Here we describe foundational genome resources and their application to the analysis of germplasm currently in use in West African breeding programs. Resources developed from the African cultivar IT97K-499-35 include a whole-genome shotgun (WGS) assembly, a bacterial artificial chromosome (BAC) physical map, and assembled sequences from 4355 BACs. These resources and WGS sequences of an additional 36 diverse cowpea accessions supported the development of a genotyping assay for 51 128 SNPs, which was then applied to five bi-parental RIL populations to produce a consensus genetic map containing 37 372 SNPs. This genetic map enabled the anchoring of 100 Mb of WGS and 420 Mb of BAC sequences, an exploration of genetic diversity along each linkage group, and clarification of macrosynteny between cowpea and common bean. The SNP assay enabled a diversity analysis of materials from West African breeding programs. Two major subpopulations exist within those materials, one of which has significant parentage from South and East Africa and more diversity. There are genomic regions of high differentiation between subpopulations, one of which coincides with a cluster of nodulin genes. The new resources and knowledge help to define goals and accelerate the breeding of improved varieties to address food security issues related to limited-input small-holder farming and climate stress.

Fine mapping of barley locus Rps6 conferring resistance to wheat stripe rust

KEY MESSAGE

Barley resistance to wheat stripe rust has remained effective for a long time and, therefore, the genes underlying this resistance can be a valuable tool to engineer durable resistance in wheat. Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a major disease of wheat that is causing large economic losses in many wheat-growing regions of the world. Deployment of Pst resistance genes has been an effective strategy for controlling this pathogen, but many of these genes have been defeated by new Pst races. In contrast, genes providing resistance to this wheat pathogen in other grass species (nonhost resistance) have been more durable. Barley varieties (Hordeum vulgare ssp. vulgare) are predominately immune to wheat Pst, but we identified three accessions of wild barley (Hordeum vulgare ssp. spontaneum) that are susceptible to Pst. Using these accessions, we mapped a barley locus conferring resistance to Pst on the distal region of chromosome arm 7HL and designated it as Rps6. The detection of the same locus in the cultivated barley 'Tamalpais' and in the Chinese barley 'Y12' by an allelism test suggests that Rps6 may be a frequent component of barley intermediate host resistance to Pst. Using a high-density mapping population (>10,000 gametes) we precisely mapped Rps6 within a 0.14 cM region (~500 kb contig) that is colinear to regions in Brachypodium (<94 kb) and rice (<9 kb). Since no strong candidate gene was identified in these colinear regions, a dedicated positional cloning effort in barley will be required to identify Rps6. The identification of this and other barley genes conferring resistance to Pst can contribute to our understanding of the mechanisms for durable resistance against this devastating wheat pathogen.

The Relationships between Development and Low Temperature Tolerance in Barley Near Isogenic Lines Differing for Flowering Behavior

Flowering time, vernalization requirement, photoperiod sensitivity and low temperature tolerance are key traits in the Triticeae. We characterized a set of isogenic genetic stocks-representing single and pairwise substitutions of spring alleles at the VRN-H1, VRN-H2 and VRN-H3 loci in a winter barley background-at the structural, functional and phenotypic levels. High density mapping with reference to the barley genome sequence confirmed that in all cases target VRN alleles were present in the near isogenic lines (NILs) and allowed estimates of introgression size (at the genetic and physical levels) and gene content. Expression data corroborated the structural and phenotypic results. The latter confirmed that substitution of a spring allele at any of the VRN loci is sufficient to eliminate vernalization requirement. There was no significant change in low temperature tolerance with substitution of a spring allele at VRN-H2, but there were significant losses in cold tolerance with substitutions at VRN-H1 and VRN-H3. Reductions in cold tolerance are ascribed to an accelerated transition from the vegetative to reproductive state. The set of NILs will be a rich resource for understanding the genetics of vernalization, low temperature tolerance and other traits encoded/regulated by genes within the introgressed intervals.

Sequencing of 15 622 gene-bearing BACs clarifies the gene-dense regions of the barley genome

Barley (Hordeum vulgare L.) possesses a large and highly repetitive genome of 5.1 Gb that has hindered the development of a complete sequence. In 2012, the International Barley Sequencing Consortium released a resource integrating whole-genome shotgun sequences with a physical and genetic framework. However, because only 6278 bacterial artificial chromosome (BACs) in the physical map were sequenced, fine structure was limited. To gain access to the gene-containing portion of the barley genome at high resolution, we identified and sequenced 15 622 BACs representing the minimal tiling path of 72 052 physical-mapped gene-bearing BACs. This generated ~1.7 Gb of genomic sequence containing an estimated 2/3 of all Morex barley genes. Exploration of these sequenced BACs revealed that although distal ends of chromosomes contain most of the gene-enriched BACs and are characterized by high recombination rates, there are also gene-dense regions with suppressed recombination. We made use of published map-anchored sequence data from Aegilops tauschii to develop a synteny viewer between barley and the ancestor of the wheat D-genome. Except for some notable inversions, there is a high level of collinearity between the two species. The software HarvEST:Barley provides facile access to BAC sequences and their annotations, along with the barley-Ae. tauschii synteny viewer. These BAC sequences constitute a resource to improve the efficiency of marker development, map-based cloning, and comparative genomics in barley and related crops. Additional knowledge about regions of the barley genome that are gene-dense but low recombination is particularly relevant.

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

New sources of genetic diversity must be incorporated into plant breeding programs if they are to continue increasing grain yield and quality, and tolerance to abiotic and biotic stresses. Germplasm collections provide a source of genetic and phenotypic diversity, but characterization of these resources is required to increase their utility for breeding programs. We used a barley SNP iSelect platform with 7,842 SNPs to genotype 2,417 barley accessions sampled from the USDA National Small Grains Collection of 33,176 accessions. Most of the accessions in this core collection are categorized as landraces or cultivars/breeding lines and were obtained from more than 100 countries. Both STRUCTURE and principal component analysis identified five major subpopulations within the core collection, mainly differentiated by geographical origin and spike row number (an inflorescence architecture trait). Different patterns of linkage disequilibrium (LD) were found across the barley genome and many regions of high LD contained traits involved in domestication and breeding selection. The genotype data were used to define 'mini-core' sets of accessions capturing the majority of the allelic diversity present in the core collection. These 'mini-core' sets can be used for evaluating traits that are difficult or expensive to score. Genome-wide association studies (GWAS) of 'hull cover', 'spike row number', and 'heading date' demonstrate the utility of the core collection for locating genetic factors determining important phenotypes. The GWAS results were referenced to a new barley consensus map containing 5,665 SNPs. Our results demonstrate that GWAS and high-density SNP genotyping are effective tools for plant breeders interested in accessing genetic diversity in large germplasm collections.

Barley genetic variation: implications for crop improvement

Genetic variation is crucial for successful barley improvement. Genomic technologies are improving dramatically and are providing access to the genetic diversity within this important crop species. Diverse collections of barley germplasm are being assembled and mined via genome-wide association studies and the identified variation can be linked to the barley sequence assembly. Introgression of favorable alleles via marker-assisted selection is now faster and more efficient due to the availability of single nucleotide polymorphism platforms. High-throughput genotyping is also making genomic selection an essential tool in modern barley breeding. Contemporary plant breeders now benefit from publicly available user-friendly databases providing genotypic and phenotypic information on large numbers of barley accessions. These resources facilitate access to allelic variation. In this review we explore how the most recent genomics and molecular breeding advances are changing breeding practices. The Coordinated Agricultural Projects (CAPs), Barley CAP and Triticeae CAP coupled with international collaborations, are discussed in detail as examples of a collaborative approach to exploit diverse germplasm resources for barley improvement.

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

Next-generation whole-genome shotgun assemblies of complex genomes are highly useful, but fail to link nearby sequence contigs with each other or provide a linear order of contigs along individual chromosomes. Here, we introduce a strategy based on sequencing progeny of a segregating population that allows de novo production of a genetically anchored linear assembly of the gene space of an organism. We demonstrate the power of the approach by reconstructing the chromosomal organization of the gene space of barley, a large, complex and highly repetitive 5.1 Gb genome. We evaluate the robustness of the new assembly by comparison to a recently released physical and genetic framework of the barley genome, and to various genetically ordered sequence-based genotypic datasets. The method is independent of the need for any prior sequence resources, and will enable rapid and cost-efficient establishment of powerful genomic information for many species.

The detection of QTLs in barley associated with endosperm hardness, grain density, grain size and malting quality using rapid phenotyping tools

Using a barley mapping population, 'Vlamingh' × 'Buloke' (V × B), whole grain analyses were undertaken for physical seed traits and malting quality. Grain density and size were predicted by digital image analysis (DIA), while malt extract and protein content were predicted using near infrared (NIR) analysis. Validation of DIA and NIR algorithms confirmed that data for QTL analysis was highly correlated (R (2) > 0.82), with high RPD values (the ratio of the standard error of prediction to the standard deviation, 2.31-9.06). Endosperm hardness was measured on this mapping population using the single kernel characterisation system. Grain density and endosperm hardness were significantly inter-correlated in all three environments (r > 0.22, P < 0.001); however, other grain components were found to interact with the traits. QTL for these traits were also found on different genomic regions, for example, grain density QTLs were found on chromosomes 2H and 6H, whereas endosperm hardness QTLs were found on 1H, 5H, and 7H. In this study, the majority of the genomic regions associated with grain texture were also coincident with QTLs for grain size, yield, flowering date and/or plant development genes. This study highlights the complexity of genomic regions associated with the variation of endosperm hardness and grain density, and their relationships with grain size traits, agronomic-related traits, and plant development loci.

Distribution, functional impact, and origin mechanisms of copy number variation in the barley genome

Background

There is growing evidence for the prevalence of copy number variation (CNV) and its role in phenotypic variation in many eukaryotic species. Here we use array comparative genomic hybridization to explore the extent of this type of structural variation in domesticated barley cultivars and wild barleys.

Results

A collection of 14 barley genotypes including eight cultivars and six wild barleys were used for comparative genomic hybridization. CNV affects 14.9% of all the sequences that were assessed. Higher levels of CNV diversity are present in the wild accessions relative to cultivated barley. CNVs are enriched near the ends of all chromosomes except 4H, which exhibits the lowest frequency of CNVs. CNV affects 9.5% of the coding sequences represented on the array and the genes affected by CNV are enriched for sequences annotated as disease-resistance proteins and protein kinases. Sequence-based comparisons of CNV between cultivars Barke and Morex provided evidence that DNA repair mechanisms of double-strand breaks via single-stranded annealing and synthesis-dependent strand annealing play an important role in the origin of CNV in barley.

Conclusions

We present the first catalog of CNVs in a diploid Triticeae species, which opens the door for future genome diversity research in a tribe that comprises the economically important cereal species wheat, barley, and rye. Our findings constitute a valuable resource for the identification of CNV affecting genes of agronomic importance. We also identify potential mechanisms that can generate variation in copy number in plant genomes.

Transcriptome analysis of a barley breeding program examines gene expression diversity and reveals target genes for malting quality improvement

Background

Advanced cycle breeding utilizes crosses among elite lines and is a successful method to develop new inbreds. However, it results in a reduction in genetic diversity within the breeding population. The development of malting barley varieties requires the adherence to a narrow malting quality profile and thus the use of advanced cycle breeding strategies. Although attention has been focused on diversity in gene expression and its association with genetic diversity, there are no studies performed in a single breeding program examining the implications that consecutive cycles of breeding have on gene expression variation and identifying the variability still available for future improvement.

Results

Fifteen lines representing the historically important six-rowed malting barley breeding program of the University of Minnesota were genotyped with 1,524 SNPs, phenotypically examined for six malting quality traits, and analyzed for transcript accumulation during germination using the Barley1 GeneChip array. Significant correlation was detected between genetic and transcript-level variation. We observed a reduction in both genetic and gene expression diversity through the breeding process, although the expression of many genes have not been fixed. A high number of quality-related genes whose expression was fixed during the breeding process was identified, indicating that much of the diversity reduction was associated with the improvement of the complex phenotype "malting quality", the main goal of the University of Minnesota breeding program. We also identified 49 differentially expressed genes between the most recent lines of the program that were correlated with one or more of the six primary malting quality traits. These genes constitute potential targets for the improvement of malting quality within the breeding program.

Conclusions

The present study shows the repercussion of advanced cycle breeding on gene expression diversity within an important barley breeding program. A reduction in gene expression diversity was detected, although there is diversity still present after forty years of breeding that can exploited for future crop improvement. In addition, the identification of candidate genes for enhancing malting quality may be used to optimize the selection of targets for further improvements in this economically important phenotype.

Structural and functional characterization of a winter malting barley

The development of winter malting barley (Hordeum vulgare L.) varieties is emerging as a worldwide priority due to the numerous advantages of these varieties over spring types. However, the complexity of both malting quality and winter hardiness phenotypes makes simultaneous improvement a challenge. To obtain an understanding of the relationship between loci controlling winter hardiness and malt quality and to assess the potential for breeding winter malting barley varieties, we structurally and functionally characterized the six-row accession "88Ab536", a cold-tolerant line with superior malting quality characteristics that derives from the cross of NE76129/Morex//Morex. We used 4,596 SNPs to construct the haplotype structure of 88Ab536 on which malting quality and winter hardiness loci reported in the literature were aligned. The genomic regions determining malting quality and winter hardiness traits have been defined in this founder germplasm, which will assist breeders in targeting regions for marker-assisted selection. The Barley1 GeneChip array was used to functionally characterize 88Ab536 during malting. Its gene expression profile was similar to that of the archetypical malting variety Morex, which is consistent with their similar malting quality characteristics. The characterization of 88Ab536 has increased our understanding of the genetic relationships of malting quality and winter hardiness, and will provide a genetic foundation for further development of more cold-tolerant varieties that have malt quality characteristics that meet or exceed current benchmarks.

Microspore embryogenesis: assignment of genes to embryo formation and green vs. albino plant production

Plant microspores can be reprogrammed from their normal pollen development to an embryogenic route in a process termed microspore embryogenesis or androgenesis. Stress treatment has a critical role in this process, inducing the dedifferentiation of microspores and conditioning the following androgenic response. In this study, we have used three barley doubled haploid lines with similar genetic background but different androgenic response. The Barley1 GeneChip was used for transcriptome comparison of these lines after mannitol stress treatment, allowing the identification of 213 differentially expressed genes. Most of these genes belong to the functional categories "cell rescue, defense, and virulence"; "metabolism"; "transcription"; and "transport". These genes were grouped into clusters according to their expression profiles among lines. A principal component analysis allowed us to associate specific gene expression clusters to phenotypic variables. Genes associated with the ability of microspores to divide and form embryos were mainly involved in changes in the structure and function of membranes, efficient use of available energy sources, and cell fate. Genes related to stress response, transcription and translation regulation, and degradation of pollen-specific proteins were associated with green plant production, while expression of genes related to plastid development was associated with albino plant regeneration.