Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity¹. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value². Faba bean (Vicia fabaL.) has a high yield potential and is well suited for cultivation in temperate regions, but genomic resources are scarce. Here, we report a high-quality chromosome-scale assembly of the faba bean genome and show that it has expanded to a massive 13 Gb in size through an imbalance between the rates of amplification and elimination of retrotransposons and satellite repeats. Genes and recombination events are evenly dispersed across chromosomes and the gene space is remarkably compact considering the genome size, although with substantial copy number variation driven by tandem duplication. Demonstrating practical application of the genome sequence, we develop a targeted genotyping assay and use high-resolution genome-wide association analysis to dissect the genetic basis of seed size and hilum colour. The resources presented constitute a genomics-based breeding platform for faba bean, enabling breeders and geneticists to accelerate the improvement of sustainable protein production across the Mediterranean, subtropical and northern temperate agroecological zones.

Lentil (Lens culinaris Medik.) is a self-pollinating, diploid, annual, cool-season, food legume crop that is cultivated throughout the world. Ascochyta blight (AB), caused by Ascochyta lentis Vassilievsky, is an economically important and widespread disease of lentil. Development of cultivars with high levels of durable resistance provides an environmentally acceptable and economically feasible method for AB control. A detailed understanding of the genetic basis of AB resistance is hence highly desirable, in order to obtain insight into the number and influence of resistance genes. Genetic linkage maps based on single nucleotide polymorphisms (SNP) and simple sequence repeat (SSR) markers have been developed from three recombinant inbred line (RIL) populations. The IH x NF map contained 460 loci across 1461.6 cM, while the IH x DIG map contained 329 loci across 1302.5 cM and the third map, NF x DIG contained 330 loci across 1914.1 cM. Data from these maps were combined with a map from a previously published study through use of bridging markers to generate a consensus linkage map containing 689 loci distributed across 7 linkage groups (LGs), with a cumulative length of 2429.61 cM at an average density of one marker per 3.5 cM. Trait dissection of AB resistance was performed for the RIL populations, identifying totals of two and three quantitative trait loci (QTLs) explaining 52% and 69% of phenotypic variation for resistance to infection in the IH x DIG and IH x NF populations, respectively. Presence of common markers in the vicinity of the AB_IH1- and AB_IH2.1/AB_IH2.2-containing regions on both maps supports the inference that a common genomic region is responsible for conferring resistance and is associated with the resistant parent, Indianhead. The third QTL was derived from Northfield. Evaluation of markers associated with AB resistance across a diverse lentil germplasm panel revealed that the identity of alleles associated with AB_IH1 predicted the phenotypic responses with high levels of accuracy (~86%), and therefore have the potential to be widely adopted in lentil breeding programs. The availability of RIL-based maps, a consensus map, and validated markers linked to AB resistance provide important resources for lentil improvement.

RNA-Seq using second-generation sequencing technologies permits generation of a reference unigene set for a given species, in the absence of a well-annotated genome sequence, supporting functional genomics studies, gene characterisation and detailed expression analysis for specific morphophysiological or environmental stress response traits. A reference unigene set for lentil has been developed, consisting of 58,986 contigs and scaffolds with an N50 length of 1719 bp. Comparison to gene complements from related species, reference protein databases, previously published lentil transcriptomes and a draft genome sequence validated the current dataset in terms of degree of completeness and utility. A large proportion (98%) of unigenes were expressed in more than one tissue, at varying levels. Candidate genes associated with mechanisms of tolerance to both boron toxicity and time of flowering were identified, which can eventually be used for the development of gene-based markers. This study has provided a comprehensive, assembled and annotated reference gene set for lentil that can be used for multiple applications, permitting identification of genes for pathway-specific expression analysis, genetic modification approaches, development of resources for genotypic analysis, and assistance in the annotation of a future lentil genome sequence.

Field pea (Pisum sativum L.) is a cool-season grain legume that is cultivated world-wide for both human consumption and stock-feed purposes. Enhancement of genetic and genomic resources for field pea will permit improved understanding of the control of traits relevant to crop productivity and quality. Advances in second-generation sequencing and associated bioinformatics analysis now provide unprecedented opportunities for the development of such resources. The objective of this study was to perform transcriptome sequencing and characterisation from two genotypes of field pea that differ in terms of seed and plant morphological characteristics. Transcriptome sequencing was performed with RNA templates from multiple tissues of the field pea genotypes Kaspa and Parafield. Tissue samples were collected at various growth stages, and a total of 23 cDNA libraries were sequenced using Illumina high-throughput sequencing platforms. A total of 407 and 352 million paired-end reads from the Kaspa and Parafield transcriptomes, respectively were assembled into 129,282 and 149,272 contigs, which were filtered on the basis of known gene annotations, presence of open reading frames (ORFs), reciprocal matches and degree of coverage. Totals of 126,335 contigs from Kaspa and 145,730 from Parafield were subsequently selected as the reference set. Reciprocal sequence analysis revealed that c. 87 % of contigs were expressed in both cultivars, while a small proportion were unique to each genotype. Reads from different libraries were aligned to the genotype-specific assemblies in order to identify and characterise expression of contigs on a tissue-specific basis, of which 87 % were expressed in more than one tissue, while others showed distinct expression patterns in specific tissues, providing unique transcriptome signatures. This study provided a comprehensive assembled and annotated transcriptome set for field pea that can be used for development of genetic markers, in order to assess genetic diversity, construct linkage maps, perform trait-dissection and implement whole-genome selection strategies in varietal improvement programs, as well to identify target genes for genetic modification approaches on the basis of annotation and expression analysis. In addition, the reference field pea transcriptome will prove highly valuable for comparative genomics studies and construction of a finalised genome sequence.

Within the Brassicaceae, six species from the genus Brassica are widely cultivated throughout the world as oilseed, condiment, fodder or vegetable crops. The genetic relationships among the six Brassica species are described by U’s triangle model. Extensive shared traits and diverse morphotypes among Brassica species make identification and classification based on phenotypic data alone challenging and unreliable, especially when dealing with large germplasm collections. Consequently, a major issue for genebank collections is ensuring the correct identification of species. Molecular genotyping based on simple sequence repeat (SSR) marker sequencing or the Illumina Infinium Brassica napus 60K single nucleotide polymorphism (SNP) array has been used to identify species and assess genetic diversity of Brassica collections. However, these methods are technically challenging, expensive and time-consuming, making them unsuitable for routine or rapid screening of Brassica accessions for germplasm management. A cheaper, faster and simpler method for Brassica species identification is described here. A multiplex polymerase chain reaction (MPCR) consisting of new and existing primers specific to the Brassica A, B and C genomes was able to reliably distinguish all six Brassica species in the triangle of U with 16 control samples of known species identity. Further validation against 120 Brassica accessions previously genotyped showed that the MPCR is highly accurate and comparable to more advanced techniques such as SSR marker sequencing or the Illumina Infinium B. napus 60K SNP array. In addition, the MPCR was sensitive enough to detect seed contaminations in pooled seed samples of Brassica accessions. A cheap and fast multiplex PCR assay for identification of Brassica species in the triangle of U was developed and validated in this study. The MPCR assay can be readily implemented in any basic molecular laboratory and should prove useful for the management of Brassica germplasm collections in genebanks.

Lentil (Lens culinaris Medik.) is a diploid (2n = 2x = 14), self-pollinating, cool-season, grain legume that is cultivated worldwide and is highly valuable due to its high protein content. However, lentil production is constrained by many factors including biotic stresses, majority of which are fungal diseases such as ascochyta blight (AB), fusarium wilt, rust, stemphylium blight, anthracnose and botrytis grey mould. Among various diseases, AB is a major -problem in many lentil-producing countries and can significantly reduce crop production. Breeding for AB resistance has been a priority for breeding programs across the globe and consequently, a number of resistance sources have been identified and extensively exploited. In order to increase the efficiency of combining genes from different genetic backgrounds, molecular genetic tools can be integrated with conventional breeding methods. A range of genetic linkage maps have been generated based on DNA-based markers, and quantitative trait loci (QTLs) for AB resistance have been identified. Molecular markers linked to these QTLs may potentially be used for efficient pyramiding of the AB disease resistance genes. Significant genomic resources have been established to identify and characterise resistance genes, including an integrated genetic map, expressed sequence tag libraries, gene based markers and draft genome sequences. These resources are already being utilised for lentil crop improvement, to more effectively select for disease resistance, as a case study of the Australian breeding program will show. The combination of genomic resources, effective molecular genetic tools and high resolution phenotyping tools will improve the efficiency of selection for ascochyta blight resistance and accelerate varietal development of global lentil breeding programs.