Genomic selection for target traits in the Australian lentil breeding program

Genomic selection (GS) uses associations between markers and phenotypes to predict the breeding values of individuals. It can be applied early in the breeding cycle to reduce the cross-to-cross generation interval and thereby increase genetic gain per unit of time. The development of cost-effective, high-throughput genotyping platforms has revolutionized plant breeding programs by enabling the implementation of GS at the scale required to achieve impact. As a result, GS is becoming routine in plant breeding, even in minor crops such as pulses. Here we examined 2,081 breeding lines from Agriculture Victoria’s national lentil breeding program for a range of target traits including grain yield, ascochyta blight resistance, botrytis grey mould resistance, salinity and boron stress tolerance, 100-grain weight, seed size index and protein content. A broad range of narrow-sense heritabilities was observed across these traits (0.24-0.66). Genomic prediction models were developed based on 64,781 genome-wide SNPs using Bayesian methodology and genomic estimated breeding values (GEBVs) were calculated. Forward cross-validation was applied to examine the prediction accuracy of GS for these targeted traits. The accuracy of GEBVs was consistently higher (0.34-0.83) than BLUP estimated breeding values (EBVs) (0.22-0.54), indicating a higher expected rate of genetic gain with GS. GS-led parental selection using early generation breeding materials also resulted in higher genetic gain compared to BLUP-based selection performed using later generation breeding lines. Our results show that implementing GS in lentil breeding will fast track the development of high-yielding cultivars with increased resistance to biotic and abiotic stresses, as well as improved seed quality traits

Transcriptome sequencing of field pea and faba bean for discovery and validation of SSR genetic markers

BACKGROUND: Field pea (Pisum sativum L.) and faba bean (Vicia faba L.) are cool-season grain legume species that provide rich sources of food for humans and fodder for livestock. To date, both species have been relative ‘genomic orphans’ due to limited availability of genetic and genomic information. A significant enrichment of genomic resources is consequently required in order to understand the genetic architecture of important agronomic traits, and to support germplasm enhancement, genetic diversity, population structure and demographic studies. RESULTS: cDNA samples obtained from various tissue types of specific field pea and faba bean genotypes were sequenced using 454 Roche GS FLX Titanium technology. A total of 720,324 and 304,680 reads for field pea and faba bean, respectively, were de novo assembled to generate sets of 70,682 and 60,440 unigenes. Consensus sequences were compared against the genome of the model legume species Medicago truncatula Gaertn., as well as that of the more distantly related, but better-characterised genome of Arabidopsis thaliana L.. In comparison to M. truncatula coding sequences, 11,737 and 10,179 unique hits were obtained from field pea and faba bean. Totals of 22,057 field pea and 18,052 faba bean unigenes were subsequently annotated from GenBank. Comparison to the genome of soybean (Glycine max L.) resulted in 19,451 unique hits for field pea and 16,497 unique hits for faba bean, corresponding to c. 35% and 30% of the known gene space, respectively. Simple sequence repeat (SSR)- containing expressed sequence tags (ESTs) were identified from consensus sequences, and totals of 2,397 and 802 primer pairs were designed for field pea and faba bean. Subsets of 96 EST-SSR markers were screened for validation across modest panels of field pea and faba bean cultivars, as well as related non-domesticated species. For field pea, 86 primer pairs successfully obtained amplification products from one or more template genotypes, of which 59% revealed polymorphism between 6 genotypes. In the case of faba bean, 81 primer pairs displayed successful amplification, of which 48% detected polymorphism. CONCLUSIONS: The generation of EST datasets for field pea and faba bean has permitted effective unigene identification and functional sequence annotation. EST-SSR loci were detected at incidences of 14-17%, permitting design of comprehensive sets of primer pairs. The subsets from these primer pairs proved highly useful for polymorphism detection within Pisum and Vicia germplasm.Sukhjiwan Kaur, Luke W. Pembleton, Noel O.I. Cogan, Keith W. Savin, Tony Leonforte, Jeffrey Paull, Michael Materne and John W. Forste

Evidence and Consequence of a Highly Adapted Clonal Haplotype within the Australian Ascochyta rabiei Population

The Australian Ascochyta rabiei (Pass.) Labr. (syn. Phoma rabiei) population has low genotypic diversity with only one mating type detected to date, potentially precluding substantial evolution through recombination. However, a large diversity in aggressiveness exists. In an effort to better understand the risk from selective adaptation to currently used resistance sources and chemical control strategies, the population was examined in detail. For this, a total of 598 isolates were quasi-hierarchically sampled between 2013 and 2015 across all major Australian chickpea growing regions and commonly grown host genotypes. Although a large number of haplotypes were identified (66) through short sequence repeat (SSR) genotyping, overall low gene diversity (Hexp = 0.066) and genotypic diversity (D = 0.57) was detected. Almost 70% of the isolates assessed were of a single dominant haplotype (ARH01). Disease screening on a differential host set, including three commonly deployed resistance sources, revealed distinct aggressiveness among the isolates, with 17% of all isolates identified as highly aggressive. Almost 75% of these were of the ARH01 haplotype. A similar pattern was observed at the host level, with 46% of all isolates collected from the commonly grown host genotype Genesis090 (classified as “resistant” during the term of collection) identified as highly aggressive. Of these, 63% belonged to the ARH01 haplotype. In conclusion, the ARH01 haplotype represents a significant risk to the Australian chickpea industry, being not only widely adapted to the diverse agro-geographical environments of the Australian chickpea growing regions, but also containing a disproportionately large number of aggressive isolates, indicating fitness to survive and replicate on the best resistance sources in the Australian germplasm.The research was funded by the Grains Research and Development Cooperation within project UM00052

Transcriptome sequencing of lentil based on second-generation technology permits large-scale unigene assembly and SSR marker discovery

<p>Abstract</p> <p>Background</p> <p>Lentil (<it>Lens culinaris </it>Medik.) is a cool-season grain legume which provides a rich source of protein for human consumption. In terms of genomic resources, lentil is relatively underdeveloped, in comparison to other Fabaceae species, with limited available data. There is hence a significant need to enhance such resources in order to identify novel genes and alleles for molecular breeding to increase crop productivity and quality.</p> <p>Results</p> <p>Tissue-specific cDNA samples from six distinct lentil genotypes were sequenced using Roche 454 GS-FLX Titanium technology, generating c. 1.38 × 10⁶expressed sequence tags (ESTs). <it>De novo </it>assembly generated a total of 15,354 contigs and 68,715 singletons. The complete unigene set was sequence-analysed against genome drafts of the model legume species <it>Medicago truncatula </it>and <it>Arabidopsis thaliana </it>to identify 12,639, and 7,476 unique matches, respectively. When compared to the genome of <it>Glycine max</it>, a total of 20,419 unique hits were observed corresponding to c. 31% of the known gene space. A total of 25,592 lentil unigenes were subsequently annoated from GenBank. Simple sequence repeat (SSR)-containing ESTs were identified from consensus sequences and a total of 2,393 primer pairs were designed. A subset of 192 EST-SSR markers was screened for validation across a panel 12 cultivated lentil genotypes and one wild relative species. A total of 166 primer pairs obtained successful amplification, of which 47.5% detected genetic polymorphism.</p> <p>Conclusions</p> <p>A substantial collection of ESTs has been developed from sequence analysis of lentil genotypes using second-generation technology, permitting unigene definition across a broad range of functional categories. As well as providing resources for functional genomics studies, the unigene set has permitted significant enhancement of the number of publicly-available molecular genetic markers as tools for improvement of this species.</p

Genetic analysis of global faba bean diversity, agronomic traits and selection signatures

Faba bean (Vicia faba L.) is a high-protein grain legume crop with great potential for sustainable protein production. However, little is known about the genetics underlying trait diversity. In this study, we used 21,345 high-quality SNP markers to genetically characterize 2678 faba bean genotypes. We performed genome-wide association studies of key agronomic traits using a seven-parent-MAGIC population and detected 238 significant marker-trait associations linked to 12 traits of agronomic importance. Sixty-five of these were stable across multiple environments. Using a non-redundant diversity panel of 685 accessions from 52 countries, we identified three subpopulations differentiated by geographical origin and 33 genomic regions subjected to strong diversifying selection between subpopulations. We found that SNP markers associated with the differentiation of northern and southern accessions explained a significant proportion of agronomic trait variance in the seven-parent-MAGIC population, suggesting that some of these traits were targets of selection during breeding. Our findings point to genomic regions associated with important agronomic traits and selection, facilitating faba bean genomics-based breeding

The giant diploid faba genome unlocks variation in a global protein crop

Publisher Copyright: © 2023, The Author(s).Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity1. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value2. Faba bean (Vicia faba L.) 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.Peer reviewe

Identification of boron tolerance in Brassica rapa

There has been increasing interest in developing canola quality B. juncea for low rainfall areas across Australia over the past two decades. However, B. juncea genotypes are susceptible to high levels of boron in Western Victorian soils. An understanding of the genetics and the molecular basis of boron tolerance may enable fast and accurate tolerance selection and lead to improved boron tolerance. Being an allotetraploid species, B. juncea is difficult to understand at the genetic level because of chromosomal duplication and the potential presence of multiple copies of the loci of interest. Therefore, once the tolerance genes or chromosomal loci governing tolerance are identified in the diploid progenitor genomes, B. rapa and B. nigra, boron tolerant B. juncea lines may be resynthesized. Thus, as an initial step in this process, this thesis aimed to understand the physiological, genomic and molecular mechanisms involved in boron tolerance in B. rapa. Initially, B. rapa genotypes were screened for tolerance to boron toxicity using hydroponic and soil assays. On the basis of primary root length, severity of leaf toxicity symptoms, dry matter accumulation and shoot boron uptake, the B. rapa genotypes WWY Sarson and Local were identified as the most tolerant and the B. rapa genotypes Shillong and Kaga the most susceptible to toxic boron concentrations (1000 ?M B in hydroponic assay; 54 mg B kg-1 soil in soil assay). The main mechanism of tolerance to boron toxicity in B. rapa involved reduced net boron uptake by roots, with some boron accumulation in the tap roots and partial exclusion of boron from shoots. Furthermore, boron uptake was much lower in the WWY Sarson and Local genotypes than in the Shillong genotype, despite higher rates of transpiration. This implied that an active boron efflux mechanism may be operating in the tolerant genotypes. The inheritance pattern of tolerance to boron toxicity in B. rapa genotype, WWY Sarson best fitted a Mendelian model of two major dominant and epistatic genes. A B. rapa linkage map was constructed from an intraspecific F2 population (WWY Sarson X Shillong) with ISSR, RAPD, SRAP and SSR marker loci. The linkage map spanned a total length of 874.1 cM and contained 12 linkage groups. Chisquare analysis (P < 0.05) revealed 25 dominant markers that showed segregation distortion in the F2 progeny. QTL analysis using composite interval analysis identified three significant peaks on LG2 and LG8 that were associated with primary root length and which accounted for 17% of the trait variation. Differential transcript analysis of SRAP markers following exposure to a toxic boron concentration identified up-regulation of me4+em2570bp, me2+em2650bp, me2+em1 1600bp, me2+em1800bp and me4+em2500bp genes in Shillong and Kaga and down-regulation of me2+em2650bp, me2+em1 1600bp, me2+em1800bp and me1+em21200bp genes in WWY Sarson and Local. Of these, a UDP-glycosyltransferase gene (sharing 80% similarity to the Arabidopsis thaliana homolog) was highly transcribed only in the sensitive genotype, Shillong, and may be involved in excessive boron cross-linking to the glycosyl groups present in the cell walls and/or membranes eventually causing the observed reductions in shoot and root growth

Characterisation of Faba Bean (Vicia faba L.) Transcriptome Using RNA-Seq: Sequencing, De Novo Assembly, Annotation, and Expression Analysis

RNA sequencing (RNA-Seq) is a deep sequencing method used for transcriptome profiling. RNA-Seq assemblies have successfully been used for a broad variety of applications, such as gene characterisation, functional genomic studies, and gene expression analysis, particularly useful in the absence of a well-studied genome reference sequence. This study reports on the development of reference unigene sets from faba bean using RNA-Seq. Two Australian faba bean cultivars (Doza and Farah) that differ in terms of disease resistance, breeding habit, and adaptation characteristics, and have been extensively used in breeding programs, were utilised in this study. The de novo assembly resulted in a total of 58,962 and 53,275 transcripts with approximately 67 Mbp (1588 bp N50) and 61 Mbp (1629 bp N50) for Doza and Farah, respectively. The generated transcripts have been compared to the protein and nucleotide databases of NCBI, as well as to the gene complements of several related legume species such as Medicago truncatula, soybean, and chickpea. Both assemblies were compared to previously-published faba bean transcriptome reference sets for the degree of completeness and utility. Annotation of unigenes has been performed, and patterns of tissue-specific expression identified. The gene complement derived from this comprehensive transcriptome analysis shows that faba bean, despite its complex 13 Gbp genome, compares well to other legumes in expressed gene content. This study in faba bean represents the most comprehensive reference transcriptomes from two different Australian cultivars available to date and it provides a valuable resource for future genomics-assisted breeding activities in this species

Characterisation of Faba Bean (Vicia faba L.) Transcriptome Using RNA-Seq: Sequencing, De Novo Assembly, Annotation, and Expression Analysis

RNA sequencing (RNA-Seq) is a deep sequencing method used for transcriptome profiling. RNA-Seq assemblies have successfully been used for a broad variety of applications, such as gene characterisation, functional genomic studies, and gene expression analysis, particularly useful in the absence of a well-studied genome reference sequence. This study reports on the development of reference unigene sets from faba bean using RNA-Seq. Two Australian faba bean cultivars (Doza and Farah) that differ in terms of disease resistance, breeding habit, and adaptation characteristics, and have been extensively used in breeding programs, were utilised in this study. The de novo assembly resulted in a total of 58,962 and 53,275 transcripts with approximately 67 Mbp (1588 bp N50) and 61 Mbp (1629 bp N50) for Doza and Farah, respectively. The generated transcripts have been compared to the protein and nucleotide databases of NCBI, as well as to the gene complements of several related legume species such as Medicago truncatula, soybean, and chickpea. Both assemblies were compared to previously-published faba bean transcriptome reference sets for the degree of completeness and utility. Annotation of unigenes has been performed, and patterns of tissue-specific expression identified. The gene complement derived from this comprehensive transcriptome analysis shows that faba bean, despite its complex 13 Gbp genome, compares well to other legumes in expressed gene content. This study in faba bean represents the most comprehensive reference transcriptomes from two different Australian cultivars available to date and it provides a valuable resource for future genomics-assisted breeding activities in this species