Rapid development of molecular markers by next-generation sequencing linked to a gene conferring phomopsis stem blight disease resistance for marker-assisted selection in lupin (Lupinus angustifolius L.) breeding

Development and validation of PCR marker array for molecular selection towards spring, vernalization-independent and winter, vernalization-responsive ecotypes of white lupin (Lupinus albus L.)

White lupin (Lupinus albus L.) is an ancient grain legume that is still undergoing improvement of domestication traits, including vernalization-responsiveness, providing frost tolerance and preventing winter flowering in autumn-sowing agriculture, and vernalization-independence, conferring drought escape by rapid flowering in spring-sowing. A recent genome-wide association study highlighted several loci significantly associated with the most contrasting phenotypes, including deletions in the promoter of the FLOWERING LOCUS T homolog, LalbFTc1, and some DArT-seq/silicoDArT loci. The present study aimed to develop and validate a versatile PCR marker array enabling molecular selection of spring- and winter-type white lupin ecotypes. Candidate DArT-seq and silicoDArT loci were transformed into cleaved amplified polymorphic sequence (CAPS) or derived CAPS markers. Developed markers, together with those previously published for LalbFTc1 INDELs and quantitative trait loci from linkage maps, were implemented for screening of white lupin germplasm panel subjected to 2-year phenotyping of phenology traits. Three DArT-seq, two silicoDArT and seven LalbFTc1 INDEL markers were positively validated, constituting a convenient PCR-based marker assay for rapid and accurate reselection of white lupin germplasm towards early flowering and thermoneutrality or late flowering and vernalization-responsiveness, as well as for tracking high genetic and phenotypic diversity within white lupin landraces, revealed in the present study.

Development of SNP markers linked to purple blotch resistance for marker-assisted selection in onion (Allium cepa L.) breeding

Purple blotch (PB), caused by Alternaria porri (Ellis) Cifferi, is one of the most destructive diseases of onion worldwide. Rapid development and deployment of resistant onion varieties is the most effective approach to control this disease. A single dominant gene, ApR1 was previously linked to PB resistance in onion cultivar 'Arka Kalyan'. In this study, an advanced RIL population derived from a cross between the resistant (Arka Kalyan) and susceptible (Agrifound Rose) cultivar of onion was used to fine map the resistant locus with SNP markers. Twenty plants from the RIL population, ten each with disease resistance and susceptibility trait, were subjected to restriction site-associated DNA sequencing (RAD-Seq) and generated 7388 single nucleotide polymorphisms (SNPs). Correlation analysis between marker genotypes and PB disease phenotype on the 20 plants identified 27 SNPs as candidate markers linked to ApR1 gene for PB resistance. Six candidate SNPs were converted to Kompetitive Allele-Specific PCR (KASP) markers designated as ApRsnip5, ApRsnip8, ApRsnip14, ApRsnip21, ApRsnip23 and ApRsnip25. Marker-trait association based on disease phenotyping and KASP genotyping data on 153 RILs confirmed that all six KASP markers were tightly associated with ApR1 gene within the genetic distance of 1.3 CentiMorgan (cM). ApRsnip14 co-segregated with the ApR1 locus. Further, the six KASP markers were tested on 27 onion lines with different genetic backgrounds. ApRsnip14, ApRsnip21, ApRsnip5 and ApRsnip23 not only showed the correct resistance allele in 3 resistance genotypes, but also clustered together in the remaining 24 susceptible lines. Alternatively, ApRsnip8 and ApRsnip25 exhibited false positives in two onion lines which do not have the R-gene. Overall, our results suggest that ApRsnip14 and ApRsnip23 with their close linkage to ApR1 locus and greater applicability on breeding germplasm are recommended in marker-assisted selection for PB resistance in onion breeding program.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03562-7.

Biodiscoveries within the Australian plant genus Eremophila based on international and interdisciplinary collaboration: results and perspectives on outstanding ethical dilemmas

In a cross-continental research initiative, including researchers working in Australia and Denmark, and based on joint external funding by a 3-year grant from the Novo Nordisk Foundation, we have used DNA sequencing, extensive chemical profiling and molecular networking analyses across the entire Eremophila genus to provide new knowledge on the presence of natural products and their bioactivities using polypharmocological screens. Sesquiterpenoids, diterpenoids and dimers of branched-chain fatty acids with previously unknown chemical structures were identified. The collection of plant material from the Eremophila genus was carried out according to a 'bioprospecting agreement' with the Government of Western Australia. We recognize that several Eremophila species hold immense cultural significance to Australia's First Peoples. In spite of our best intentions to ensure that new knowledge gained about the genus Eremophila and any potential future benefits are shared in an equitable manner, in accordance with the Nagoya Protocol, we encounter serious dilemmas and potential conflicts in making benefit sharing with Australia's First Peoples a reality.

Progress of Genomics-Driven Approaches for Sustaining Underutilized Legume Crops in the Post-Genomic Era

Legume crops, belonging to the Fabaceae family, are of immense importance for sustaining global food security. Many legumes are profitable crops for smallholder farmers due to their unique ability to fix atmospheric nitrogen and their intrinsic ability to thrive on marginal land with minimum inputs and low cultivation costs. Recent progress in genomics shows promise for future genetic gains in major grain legumes. Still it remains limited in minor legumes/underutilized legumes, including adzuki bean, cluster bean, horse gram, lathyrus, red clover, urd bean, and winged bean. In the last decade, unprecedented progress in completing genome assemblies of various legume crops and resequencing efforts of large germplasm collections has helped to identify the underlying gene(s) for various traits of breeding importance for enhancing genetic gain and contributing to developing climate-resilient cultivars. This review discusses the progress of genomic resource development, including genome-wide molecular markers, key breakthroughs in genome sequencing, genetic linkage maps, and trait mapping for facilitating yield improvement in underutilized legumes. We focus on 1) the progress in genomic-assisted breeding, 2) the role of whole-genome resequencing, pangenomes for underpinning the novel genomic variants underlying trait gene(s), 3) how adaptive traits of wild underutilized legumes could be harnessed to develop climate-resilient cultivars, 4) the progress and status of functional genomics resources, deciphering the underlying trait candidate genes with putative function in underutilized legumes 5) and prospects of novel breeding technologies, such as speed breeding, genomic selection, and genome editing. We conclude the review by discussing the scope for genomic resources developed in underutilized legumes to enhance their production and play a critical role in achieving the "zero hunger" sustainable development goal by 2030 set by the United Nations.

Genotyping Analysis by RAD-Seq Reads Is Useful to Assess the Genetic Identity and Relationships of Breeding Lines in Lavender Species Aimed at Managing Plant Variety Protection

Lavender species are widely distributed in their wild forms around the Mediterranean Basin and they are also cultivated worldwide as improved and registered clonal varieties. The economic interest of the species belonging to the Lavandula genus is determined by their use as ornamental plants and important source of essential oils that are destinated to the production of cosmetics, pharmaceuticals and foodstuffs. Because of the increasing number of cases of illegal commercialization of selected varieties, the protection of plant breeders’ rights has become of main relevance for the recognition of breeding companies’ royalties. With this aim, genomic tools based on molecular markers have been demonstrated to be very reliable and transferable among laboratories, and also much more informative than morphological descriptors. With the rising of the next-generation sequencing (NGS) technologies, several genotyping-by-sequencing approaches are now available. This study deals with a deep characterization of 15 varietal clones, belonging to two distinct Lavandula species, by means of restriction-site associated DNA sequencing (RAD-Seq). We demonstrated that this technology screens single nucleotide variants that enable to assess the genetic identity of individual accessions, to reconstruct genetic relationships among related breeding lines, to group them into genetically distinguishable main subclusters, and to assign their molecular lineages to distinct ancestors. Moreover, a number of polymorphic sites were identified within genes putatively involved in biosynthetic pathways related to both tissue pigmentation and terpene production, useful for breeding and/or protecting newly registered varieties. Overall, the results highlighted the presence of pure ancestries and interspecific hybrids for the analyzed Lavandula species, and demonstrated that RAD-Seq analysis is very informative and highly reliable for characterizing Lavandula clones and managing plant variety protection.

Identification of Gray Leaf Spot Disease Candidate Gene in Narrow-Leafed Lupin (Lupinus angustifolius L.)

Selection for resistance against gray leaf spot (GLS) is a major objective in the lupin breeding programs. A segregation ratio of 1:1 (resistant:susceptible) in F8 recombinant inbred lines (RIL8) derived from a cross between a breeding line 83A:476 (resistant to GLS) and a wild accession P27255 (susceptible to GLS) indicated that GLS was controlled by a single major gene. To develop molecular markers linked to GLS, in the beginning, only 11 resistant lines and six susceptible lines from the 83A:476 and P27255 population were genotyped with MFLP markers, and three MFLP markers were identified to be co-segregated with GLS. This method was very efficient, but the markers were located outside of the gene, and could not be used in other germplasms. Then QTL analysis and fine mapping were conducted to identify the gene. Finally, the gene was narrowed down to a 241-kb region containing two disease resistance genes. To further identify the candidate gene, DNA variants between accessions Tanjil (resistant to GLS) and Unicrop (susceptible to GLS) were analyzed. The results indicated that only one SNP was detected in the 241 kb region. This SNP was located in the TMV resistance protein N-like gene region and also identified between 83A:476 and P27255. Genotyping the Tanjil/Unicrop RIL8 population showed that this SNP co-segregated with GLS resistance. The phylogenetic tree analysis of this gene among 18 lupin accessions indicates that Australian resistant breeding line/varieties were clustered into one group and carry two resistant alleles, while susceptible accessions were clustered into different groups.

The Resistance of Narrow-Leafed Lupin to Diaporthe toxica Is Based on the Rapid Activation of Defense Response Genes

Narrow-leafed lupin (Lupinus angustifolius L.) is a grain legume crop that is advantageous in animal nutrition due to its high protein content; however, livestock grazing on stubble may develop a lupinosis disease that is related to toxins produced by a pathogenic fungus, Diaporthe toxica. Two major unlinked alleles, Phr1 and PhtjR, confer L. angustifolius resistance to this fungus. Besides the introduction of these alleles into modern cultivars, the molecular mechanisms underlying resistance remained unsolved. In this study, resistant and susceptible lines were subjected to differential gene expression profiling in response to D. toxica inoculation, spanning the progress of the infection from the early to latent phases. High-throughput sequencing of stem transcriptome and PCR quantification of selected genes were performed. Gene Ontology term analysis revealed that an early (24 h) response in the resistant germplasm encompassed activation of genes controlling reactive oxygen species and oxylipin biosynthesis, whereas in the susceptible germplasm, it comprised induction of xyloglucan endotransglucosylases/hydrolases. During the first five days of the infection, the number of genes with significantly altered expressions was about 2.6 times higher in resistant lines than in the susceptible line. Global transcriptome reprogramming involving the activation of defense response genes occurred in lines conferring Phr1 and PhtjR resistance alleles about 4–8 days earlier than in the susceptible germplasm.

Innovative transcriptome-based genotyping highlights environmentally responsive genes for phenology, growth and yield in a non-model grain legume

The narrow-leafed lupin, Lupinus angustifolius L., is a grain legume crop, cultivated both as a green manure and as a source of protein for animal feed and human food production. During its domestication process, numerous agronomic traits were improved, however, only two trait-related genes were identified hitherto, both by linkage mapping. Genome-wide association studies (GWAS), exploiting genomic sequencing, did not select any novel candidate gene. In the present study, an innovative method of 3'-end reduced representation transcriptomic profiling, a massive analysis of cDNA ends, has been used for genotyping of 126 L. angustifolius lines surveyed by field phenotyping. Significant genotype × environment interactions were identified for all phenology and yield traits analysed. Principal component analysis of population structure evidenced European domestication bottlenecks, visualized by clustering of breeding materials and cultivars. GWAS provided contribution towards deciphering vernalization pathway in legumes, and, apart from highlighting known domestication loci (Ku/Julius and mol), designated novel candidate genes for L. angustifolius traits. Early phenology was associated with genes from vernalization, cold-responsiveness and phosphatidylinositol signalling pathways whereas high yield with genes controlling photosynthesis performance and abiotic stress (drought or heat) tolerance. PCR-based toolbox was developed and validated to enable tracking desired alleles in marker-assisted selection. Narrow-leafed lupin was genotyped with an innovative method of transcriptome profiling and phenotyped for phenology, growth and yield traits in field. Early phenology was found associated with genes from cold-response, vernalization and phosphatidylinositol signalling pathways, whereas high yield with genes running photosystem II and drought or heat stress response. Key loci were supplied with PCR-based toolbox for marker-assisted selection.

Development of PCR-based markers and whole-genome selection model for anthracnose resistance in white lupin (Lupinus albus L.)

White lupin (Lupinus albus L.) is a high-protein grain legume crop, grown since ancient Greece and Rome. Despite long domestication history, its cultivation remains limited, partly because of susceptibility to anthracnose. Only some late-flowering, bitter, low-yielding landraces from Ethiopian mountains displayed resistance to this devastating disease. The resistance is controlled by various genes, thereby complicating the breeding efforts. The objective of this study was developing tools for molecular tracking of Ethiopian resistance genes based on genotyping-by-sequencing (GBS) data, envisaging linkage mapping and genomic selection approaches. Twenty GBS markers from two major quantitative trait loci (QTLs), antr04_1/antr05_1 and antr04_2/antr05_2, were converted to PCR-based markers using assigned transcriptome sequences. Newly developed markers improved mapping resolution around both anthracnose resistance loci, providing more precise QTL estimation. PCR-based screening of diversified domesticated and primitive germplasm revealed the high specificity of two markers for the antr04_1/antr05_1 locus (TP222136 and TP47110) and one for the antr04_2/antr05_2 locus (TP338761), highlighted by simple matching coefficients of 0.96 and 0.89, respectively. Moreover, a genomic selection approach based on GBS data of a recombinant inbred line mapping population was assessed, providing an average predictive ability of 0.56. These tools can be used for preselection of candidate white lupin germplasm for anthracnose resistance assays.

A Tale of Two Families: Whole Genome and Segmental Duplications Underlie Glutamine Synthetase and Phosphoenolpyruvate Carboxylase Diversity in Narrow-Leafed Lupin (Lupinus angustifolius L.)

Narrow-leafed lupin (Lupinus angustifolius L.) has recently been supplied with advanced genomic resources and, as such, has become a well-known model for molecular evolutionary studies within the legume family—a group of plants able to fix nitrogen from the atmosphere. The phylogenetic position of lupins in Papilionoideae and their evolutionary distance to other higher plants facilitates the use of this model species to improve our knowledge on genes involved in nitrogen assimilation and primary metabolism, providing novel contributions to our understanding of the evolutionary history of legumes. In this study, we present a complex characterization of two narrow-leafed lupin gene families—glutamine synthetase (GS) and phosphoenolpyruvate carboxylase (PEPC). We combine a comparative analysis of gene structures and a synteny-based approach with phylogenetic reconstruction and reconciliation of the gene family and species history in order to examine events underlying the extant diversity of both families. Employing the available evidence, we show the impact of duplications on the initial complement of the analyzed gene families within the genistoid clade and posit that the function of duplicates has been largely retained. In terms of a broader perspective, our results concerning GS and PEPC gene families corroborate earlier findings pointing to key whole genome duplication/triplication event(s) affecting the genistoid lineage.

Validation of Diaporthe toxica resistance markers in European Lupinus angustifolius germplasm and identification of novel resistance donors for marker-assisted selection

The fungus, Diaporthe toxica, anamorph Phomopsis sp., previously classified as P. leptostromiformis, is a plant endophyte and occasional pathogen, causing Phomopsis stem blight. This disease is damaging not only to lupins but also to the animals grazing on infected plants, due to the toxic secondary metabolites called phomopsins. The aim of this work was to validate markers for resistance to Phomopsis stem blight in narrow-leafed lupins and identify novel germplasm with increased levels of resistance to the disease. Plant inoculations were performed using ten isolates of D. toxica, originating from Australia and Poland. The European core collection of L. angustifolius was evaluated both in a controlled environment and with field experiments to classify the accessions based on their resistance to the disease. Simultaneously, the accessions were assayed with disease resistance markers to identify donors of hypothetical resistance alleles. We have found that the European lupin germplasm collection preserves wild and domesticated donors of at least two resistance genes to Phomopsis stem blight, including Phr1 and PhtjR. Molecular markers PhtjM7, InDel2, and InDel10, tagging PhtjR gene, were applicable for marker-assisted selection targeting the European gene pool with an expected accuracy of 95%. None of diagnostic markers for the Phr1 locus was found useful for European breeding programs; two existing markers Ph258M1 and Ph258M2 were unreliable, due to a high percentage of false-positive results (up to 58%) and a high recombination rate between markers (~ 30%).

Genome-wide RAD sequencing to identify a sex-specific marker in Chinese giant salamander Andrias davidianus

Background

Chinese giant salamander Andrias davidianus is an endangered species. The success of artificial breeding provides a useful way to protect this species. However, the method to identify the sex and mechanism of sex determination were unclear which hinder the improvement of the artificial breeding. Detection of a sex specific marker provides an effective approach to identify genetic sex and investigate the sex determination mechanism.

Results

We used restriction-site-associated DNA (RAD) sequencing to isolate a sex-specific genetic marker in A. davidianus to expand knowledge of the sex determination mechanism. Four male and four female specimens were subjected to RAD sequencing, which generated 934,072,989 reads containing approximately 134.4 Gb of sequences. The first round of comparison of the assembled sequence against the opposite sex raw reads revealed 19,097 female and 17,994 male unmatched sequences. Subsequently, 19,097 female sequences were subjected to a BLAST search against male genomic data, which revealed 308 sequences unmapped to the male genome. One hundred of these were randomly selected and validated by PCR in five male and five female specimens, and four putative sex-specific sequences were produced. Further validation was performed by PCR in another 24 females and 24 males, and all female individuals exhibited the expected specific bands, while the males did not. To apply the sex-specific marker, three specimens reversed from genetic female to physiological male were found in a group exposed to elevated temperature, and 13 individuals reversed from genetic male to physiological female were obtained in a 17β-estradiol exposed group.

Conclusion

This is the first report of a sex-specific marker in A. davidianus and may have potential for elucidation of its sex determination mechanism and, hence, its conservation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5771-5) contains supplementary material, which is available to authorized users.

The Use of Lupin as a Source of Protein in Animal Feeding: Genomic Tools and Breeding Approaches

Livestock production in the European Union EU is highly dependent on imported soybean, exposing the livestock farming system to risks related to the global trade of soybean. Lupin species could be a realistic sustainable alternative source of protein for animal feeding. Lupinus is a very diverse genus with many species. However, only four of them—namely, L. albus, L. angustifolius, L. luteus and L. mutabilis—are cultivated. Their use in livestock farming systems has many advantages in relation to economic and environmental impact. Generally, lupin grains are characterized by high protein content, while their oil content is relatively low but of high quality. On the other hand, the presence of quinolizidine alkaloids and their specific carbohydrate composition are the main antinutritional factors that prevent their use in animal feeding. This research is mainly related to L. albus and to L. angustifolius, and to a lesser extent, to L. lauteus and L. mutabilis. The breeding efforts are mostly focused on yield stabilization, resistance to biotic and abiotic stresses, biochemical structure associated with seed quality and late maturing. Progress is made in improving lupin with respect to the seed quality, as well as the tolerance to biotic and abiotic stress. It has to be noted that modern cultivars, mostly of L. albus and L. angustifolius, contain low levels of alkaloids. However, for future breeding efforts, the implementation of marker-assisted selection and the available genomic tools is of great importance.

A high-density consensus linkage map of white lupin highlights synteny with narrow-leafed lupin and provides markers tagging key agronomic traits

White lupin (Lupinus albus L.) is a valuable source of seed protein, carbohydrates and oil, but requires genetic improvement to attain its agronomic potential. This study aimed to (i) develop a new high-density consensus linkage map based on new, transcriptome-anchored markers; (ii) map four important agronomic traits, namely, vernalization requirement, seed alkaloid content, and resistance to anthracnose and Phomopsis stem blight; and, (iii) define regions of synteny between the L. albus and narrow-leafed lupin (L. angustifolius L.) genomes. Mapping of white lupin quantitative trait loci (QTLs) revealed polygenic control of vernalization responsiveness and anthracnose resistance, as well as a single locus regulating seed alkaloid content. We found high sequence collinearity between white and narrow-leafed lupin genomes. Interestingly, the white lupin QTLs did not correspond to previously mapped narrow-leafed lupin loci conferring vernalization independence, anthracnose resistance, low alkaloids and Phomopsis stem blight resistance, highlighting different genetic control of these traits. Our suite of allele-sequenced and PCR validated markers tagging these QTLs is immediately applicable for marker-assisted selection in white lupin breeding. The consensus map constitutes a platform for synteny-based gene cloning approaches and can support the forthcoming white lupin genome sequencing efforts.

A comprehensive draft genome sequence for lupin (Lupinus angustifolius), an emerging health food: insights into plant-microbe interactions and legume evolution

Lupins are important grain legume crops that form a critical part of sustainable farming systems, reducing fertilizer use and providing disease breaks. It has a basal phylogenetic position relative to other crop and model legumes and a high speciation rate. Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is gaining popularity as a health food, which is high in protein and dietary fibre but low in starch and gluten-free. We report the draft genome assembly (609 Mb) of NLL cultivar Tanjil, which has captured >98% of the gene content, sequences of additional lines and a dense genetic map. Lupins are unique among legumes and differ from most other land plants in that they do not form mycorrhizal associations. Remarkably, we find that NLL has lost all mycorrhiza-specific genes, but has retained genes commonly required for mycorrhization and nodulation. In addition, the genome also provided candidate genes for key disease resistance and domestication traits. We also find evidence of a whole-genome triplication at around 25 million years ago in the genistoid lineage leading to Lupinus. Our results will support detailed studies of legume evolution and accelerate lupin breeding programmes.

Genotyping-by-sequencing approaches to characterize crop genomes: choosing the right tool for the right application

In the last decade, the revolution in sequencing technologies has deeply impacted crop genotyping practice. New methods allowing rapid, high-throughput genotyping of entire crop populations have proliferated and opened the door to wider use of molecular tools in plant breeding. These new genotyping-by-sequencing (GBS) methods include over a dozen reduced-representation sequencing (RRS) approaches and at least four whole-genome resequencing (WGR) approaches. The diversity of methods available, each often producing different types of data at different cost, can make selection of the best-suited method seem a daunting task. We review the most common genotyping methods used today and compare their suitability for linkage mapping, genomewide association studies (GWAS), marker-assisted and genomic selection and genome assembly and improvement in crops with various genome sizes and complexity. Furthermore, we give an outline of bioinformatics tools for analysis of genotyping data. WGR is well suited to genotyping biparental cross populations with complex, small- to moderate-sized genomes and provides the lowest cost per marker data point. RRS approaches differ in their suitability for various tasks, but demonstrate similar costs per marker data point. These approaches are generally better suited for de novo applications and more cost-effective when genotyping populations with large genomes or high heterozygosity. We expect that although RRS approaches will remain the most cost-effective for some time, WGR will become more widespread for crop genotyping as sequencing costs continue to decrease.

Identification of Genome-Wide Variants and Discovery of Variants Associated with Brassica rapa Clubroot Resistance Gene Rcr1 through Bulked Segregant RNA Sequencing

Clubroot, caused by Plasmodiophora brassicae, is an important disease on Brassica species worldwide. A clubroot resistance gene, Rcr1, with efficacy against pathotype 3 of P. brassicae, was previously mapped to chromosome A03 of B. rapa in pak choy cultivar "Flower Nabana". In the current study, resistance to pathotypes 2, 5 and 6 was shown to be associated with Rcr1 region on chromosome A03. Bulked segregant RNA sequencing was performed and short read sequences were assembled into 10 chromosomes of the B. rapa reference genome v1.5. For the resistant (R) bulks, a total of 351.8 million (M) sequences, 30,836.5 million bases (Mb) in length, produced 120-fold coverage of the reference genome. For the susceptible (S) bulks, 322.9 M sequences, 28,216.6 Mb in length, produced 109-fold coverage. In total, 776.2 K single nucleotide polymorphisms (SNPs) and 122.2 K insertion / deletion (InDels) in R bulks and 762.8 K SNPs and 118.7 K InDels in S bulks were identified; each chromosome had about 87% SNPs and 13% InDels, with 78% monomorphic and 22% polymorphic variants between the R and S bulks. Polymorphic variants on each chromosome were usually below 23%, but made up 34% of the variants on chromosome A03. There were 35 genes annotated in the Rcr1 target region and variants were identified in 21 genes. The numbers of poly variants differed significantly among the genes. Four out of them encode Toll-Interleukin-1 receptor / nucleotide-binding site / leucine-rich-repeat proteins; Bra019409 and Bra019410 harbored the higher numbers of polymorphic variants, which indicates that they are more likely candidates of Rcr1. Fourteen SNP markers in the target region were genotyped using the Kompetitive Allele Specific PCR method and were confirmed to associate with Rcr1. Selected SNP markers were analyzed with 26 recombinants obtained from a segregating population consisting of 1587 plants, indicating that they were completely linked to Rcr1. Nine SNP markers were used for marker-assisted introgression of Rcr1 into B. napus canola from B. rapa, with 100% accuracy in this study.

Identification of Genome-Wide Variants and Discovery of Variants Associated with Brassica rapa Clubroot Resistance Gene Rcr1 through Bulked Segregant RNA Sequencing

Clubroot, caused by Plasmodiophora brassicae, is an important disease on Brassica species worldwide. A clubroot resistance gene, Rcr1, with efficacy against pathotype 3 of P. brassicae, was previously mapped to chromosome A03 of B. rapa in pak choy cultivar “Flower Nabana”. In the current study, resistance to pathotypes 2, 5 and 6 was shown to be associated with Rcr1 region on chromosome A03. Bulked segregant RNA sequencing was performed and short read sequences were assembled into 10 chromosomes of the B. rapa reference genome v1.5. For the resistant (R) bulks, a total of 351.8 million (M) sequences, 30,836.5 million bases (Mb) in length, produced 120-fold coverage of the reference genome. For the susceptible (S) bulks, 322.9 M sequences, 28,216.6 Mb in length, produced 109-fold coverage. In total, 776.2 K single nucleotide polymorphisms (SNPs) and 122.2 K insertion / deletion (InDels) in R bulks and 762.8 K SNPs and 118.7 K InDels in S bulks were identified; each chromosome had about 87% SNPs and 13% InDels, with 78% monomorphic and 22% polymorphic variants between the R and S bulks. Polymorphic variants on each chromosome were usually below 23%, but made up 34% of the variants on chromosome A03. There were 35 genes annotated in the Rcr1 target region and variants were identified in 21 genes. The numbers of poly variants differed significantly among the genes. Four out of them encode Toll-Interleukin-1 receptor / nucleotide-binding site / leucine-rich-repeat proteins; Bra019409 and Bra019410 harbored the higher numbers of polymorphic variants, which indicates that they are more likely candidates of Rcr1. Fourteen SNP markers in the target region were genotyped using the Kompetitive Allele Specific PCR method and were confirmed to associate with Rcr1. Selected SNP markers were analyzed with 26 recombinants obtained from a segregating population consisting of 1587 plants, indicating that they were completely linked to Rcr1. Nine SNP markers were used for marker-assisted introgression of Rcr1 into B. napus canola from B. rapa, with 100% accuracy in this study.

Application of whole genome re-sequencing data in the development of diagnostic DNA markers tightly linked to a disease-resistance locus for marker-assisted selection in lupin (Lupinus angustifolius)

BACKGROUND

Molecular marker-assisted breeding provides an efficient tool to develop improved crop varieties. A major challenge for the broad application of markers in marker-assisted selection is that the marker phenotypes must match plant phenotypes in a wide range of breeding germplasm. In this study, we used the legume crop species Lupinus angustifolius (lupin) to demonstrate the utility of whole genome sequencing and re-sequencing on the development of diagnostic markers for molecular plant breeding.

RESULTS

Nine lupin cultivars released in Australia from 1973 to 2007 were subjected to whole genome re-sequencing. The re-sequencing data together with the reference genome sequence data were used in marker development, which revealed 180,596 to 795,735 SNP markers from pairwise comparisons among the cultivars. A total of 207,887 markers were anchored on the lupin genetic linkage map. Marker mining obtained an average of 387 SNP markers and 87 InDel markers for each of the 24 genome sequence assembly scaffolds bearing markers linked to 11 genes of agronomic interest. Using the R gene PhtjR conferring resistance to phomopsis stem blight disease as a test case, we discovered 17 candidate diagnostic markers by genotyping and selecting markers on a genetic linkage map. A further 243 candidate diagnostic markers were discovered by marker mining on a scaffold bearing non-diagnostic markers linked to the PhtjR gene. Nine out from the ten tested candidate diagnostic markers were confirmed as truly diagnostic on a broad range of commercial cultivars. Markers developed using these strategies meet the requirements for broad application in molecular plant breeding.

CONCLUSIONS

We demonstrated that low-cost genome sequencing and re-sequencing data were sufficient and very effective in the development of diagnostic markers for marker-assisted selection. The strategies used in this study may be applied to any trait or plant species. Whole genome sequencing and re-sequencing provides a powerful tool to overcome current limitations in molecular plant breeding, which will enable plant breeders to precisely pyramid favourable genes to develop super crop varieties to meet future food demands.

No evidence for Fabaceae Gametophytic self-incompatibility being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes

BACKGROUND

Fabaceae species are important in agronomy and livestock nourishment. They have a long breeding history, and most cultivars have lost self-incompatibility (SI), a genetic barrier to self-fertilization. Nevertheless, to improve legume crop breeding, crosses with wild SI relatives of the cultivated varieties are often performed. Therefore, it is fundamental to characterize Fabaceae SI system(s). We address the hypothesis of Fabaceae gametophytic (G)SI being RNase based, by recruiting the same S-RNase lineage gene of Rosaceae, Solanaceae or Plantaginaceae SI species.

RESULTS

We first identify SSK1 like genes (described only in species having RNase based GSI), in the Trifolium pratense, Medicago truncatula, Cicer arietinum, Glycine max, and Lupinus angustifolius genomes. Then, we characterize the S-lineage T2-RNase genes in these genomes. In T. pratense, M. truncatula, and C. arietinum we identify S-RNase lineage genes that in phylogenetic analyses cluster with Pyrinae S-RNases. In M. truncatula and C. arietinum genomes, where large scaffolds are available, these sequences are surrounded by F-box genes that in phylogenetic analyses also cluster with S-pollen genes. In T. pratense the S-RNase lineage genes show, however, expression in tissues not involved in GSI. Moreover, levels of diversity are lower than those observed for other S-RNase genes. The M. truncatula and C. arietinum S-RNase and S-pollen like genes phylogenetically related to Pyrinae S-genes, are also expressed in tissues other than those involved in GSI. To address if other T2-RNases could be determining Fabaceae GSI, here we obtained a style with stigma transcriptome of Cytisus striatus, a species that shows significant difference on the percentage of pollen growth in self and cross-pollinations. Expression and polymorphism analyses of the C. striatus S-RNase like genes revealed that none of these genes, is the S-pistil gene.

CONCLUSION

We find no evidence for Fabaceae GSI being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes. There is no evidence that T2-RNase lineage genes could be determining GSI in C. striatus. Therefore, to characterize the Fabaceae S-pistil gene(s), expression analyses, levels of diversity, and segregation analyses in controlled crosses are needed for those genes showing high expression levels in the tissues where GSI occurs.

Sequencing consolidates molecular markers with plant breeding practice

Plenty of molecular markers have been developed by contemporary sequencing technologies, whereas few of them are successfully applied in breeding, thus we present a review on how sequencing can facilitate marker-assisted selection in plant breeding. The growing global population and shrinking arable land area require efficient plant breeding. Novel strategies assisted by certain markers have proven effective for genetic gains. Fortunately, cutting-edge sequencing technologies bring us a deluge of genomes and genetic variations, enlightening the potential of marker development. However, a large gap still exists between the potential of molecular markers and actual plant breeding practices. In this review, we discuss marker-assisted breeding from a historical perspective, describe the road from crop sequencing to breeding, and highlight how sequencing facilitates the application of markers in breeding practice.

Analysis of conglutin seed storage proteins across lupin species using transcriptomic, protein and comparative genomic approaches

BACKGROUND

The major proteins in lupin seeds are conglutins that have primary roles in supplying carbon, sulphur and nitrogen and energy for the germinating seedling. They fall into four families; α, β, γ and δ. Interest in these conglutins is growing as family members have been shown to have beneficial nutritional and pharmaceutical properties.

RESULTS

An in-depth transcriptome and draft genome from the narrow-leafed lupin (NLL; Lupinus angustifolius) variety, Tanjil, were examined and 16 conglutin genes were identified. Using RNAseq data sets, the structure and expression of these 16 conglutin genes were analysed across eight lupin varieties from five lupin species. Phylogenic analysis suggest that the α and γ conglutins diverged prior to lupin speciation while β and δ members diverged both prior and after speciation. A comparison of the expression of the 16 conglutin genes was performed, and in general the conglutin genes showed similar levels of RNA expression among varieties within species, but quite distinct expression patterns between lupin species. Antibodies were generated against the specific conglutin families and immunoblot analyses were used to compare the levels of conglutin proteins in various tissues and during different stages of seed development in NLL, Tanjil, confirming the expression in the seed. This analysis showed that the conglutins were expressed highly at the mature seed stage, in all lupin species, and a range of polypeptide sizes were observed for each conglutin family.

CONCLUSIONS

This study has provided substantial information on the complexity of the four conglutin families in a range of lupin species in terms of their gene structure, phylogenetic relationships as well as their relative RNA and protein abundance during seed development. The results demonstrate that the majority of the heterogeneity of conglutin polypeptides is likely to arise from post-translational modification from a limited number of precursor polypeptides rather than a large number of different genes. Overall, the results demonstrate a high degree of plasticity for conglutin expression during seed development in different lupin species.

Whole genome sequences in pulse crops: a global community resource to expedite translational genomics and knowledge-based crop improvement

Unprecedented developments in legume genomics over the last decade have resulted in the acquisition of a wide range of modern genomic resources to underpin genetic improvement of grain legumes. The genome enabled insights direct investigators in various ways that primarily include unearthing novel structural variations, retrieving the lost genetic diversity, introducing novel/exotic alleles from wider gene pools, finely resolving the complex quantitative traits and so forth. To this end, ready availability of cost-efficient and high-density genotyping assays allows genome wide prediction to be increasingly recognized as the key selection criterion in crop breeding. Further, the high-dimensional measurements of agronomically significant phenotypes obtained by using new-generation screening techniques will empower reference based resequencing as well as allele mining and trait mapping methods to comprehensively associate genome diversity with the phenome scale variation. Besides stimulating the forward genetic systems, accessibility to precisely delineated genomic segments reveals novel candidates for reverse genetic techniques like targeted genome editing. The shifting paradigm in plant genomics in turn necessitates optimization of crop breeding strategies to enable the most efficient integration of advanced omics knowledge and tools. We anticipate that the crop improvement schemes will be bolstered remarkably with rational deployment of these genome-guided approaches, ultimately resulting in expanded plant breeding capacities and improved crop performance.

Identification of sex-linked SNP markers using RAD sequencing suggests ZW/ZZ sex determination in Pistacia vera L

Background

Pistachio (Pistacia vera L.) is a dioecious species that has a long juvenility period. Therefore, development of marker-assisted selection (MAS) techniques would greatly facilitate pistachio cultivar-breeding programs. The sex determination mechanism is presently unknown in pistachio. The generation of sex-linked markers is likely to reduce time, labor, and costs associated with breeding programs, and will help to clarify the sex determination system in pistachio.

Results

Restriction site-associated DNA (RAD) markers were used to identify sex-linked markers and to elucidate the sex determination system in pistachio. Eight male and eight female F₁ progenies from a Pistacia vera L. Siirt × Bağyolu cross, along with the parents, were subjected to RAD sequencing in two lanes of a Hi-Seq 2000 sequencing platform. This generated 449 million reads, comprising approximately 37.7 Gb of sequences. There were 33,757 polymorphic single nucleotide polymorphism (SNP) loci between the parents. Thirty-eight of these, from 28 RAD reads, were detected as putative sex-associated loci in pistachio. Validation was performed by SNaPshot analysis in 42 mature F₁ progenies and in 124 cultivars and genotypes in a germplasm collection. Eight loci could distinguish sex with 100% accuracy in pistachio. To ascertain cost-effective application of markers in a breeding program, high-resolution melting (HRM) analysis was performed; four markers were found to perfectly separate sexes in pistachio. Because of the female heterogamety in all candidate SNP loci, we report for the first time that pistachio has a ZZ/ZW sex determination system. As the reported female-to-male segregation ratio is 1:1 in all known segregating populations and there is no previous report of super-female genotypes or female heteromorphic chromosomes in pistachio, it appears that the WW genotype is not viable.

Conclusion

Sex-linked SNP markers were identified and validated in a large germplasm and proved their suitability for MAS in pistachio. HRM analysis successfully validated the sex-linked markers for MAS. For the first time in dioecious pistachio, a female heterogamety ZW/ZZ sex determination system is suggested.

Transcriptome sequencing of different narrow-leafed lupin tissue types provides a comprehensive uni-gene assembly and extensive gene-based molecular markers

Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is an important grain legume crop that is valuable for sustainable farming and is becoming recognized as a human health food. NLL breeding is directed at improving grain production, disease resistance, drought tolerance and health benefits. However, genetic and genomic studies have been hindered by a lack of extensive genomic resources for the species. Here, the generation, de novo assembly and annotation of transcriptome datasets derived from five different NLL tissue types of the reference accession cv. Tanjil are described. The Tanjil transcriptome was compared to transcriptomes of an early domesticated cv. Unicrop, a wild accession P27255, as well as accession 83A:476, together being the founding parents of two recombinant inbred line (RIL) populations. In silico predictions for transcriptome-derived gene-based length and SNP polymorphic markers were conducted and corroborated using a survey assembly sequence for NLL cv. Tanjil. This yielded extensive indel and SNP polymorphic markers for the two RIL populations. A total of 335 transcriptome-derived markers and 66 BAC-end sequence-derived markers were evaluated, and 275 polymorphic markers were selected to genotype the reference NLL 83A:476 × P27255 RIL population. This significantly improved the completeness, marker density and quality of the reference NLL genetic map.

Remnants of the Legume Ancestral Genome Preserved in Gene-Rich Regions: Insights from Lupinus angustifolius Physical, Genetic, and Comparative Mapping

The narrow-leafed lupin (Lupinus angustifolius) was recently considered as a legume reference species. Genetic resources have been developed, including a draft genome sequence, linkage maps, nuclear DNA libraries, and cytogenetic chromosome-specific landmarks. Here, we used a complex approach, involving DNA fingerprinting, sequencing, genetic mapping, and molecular cytogenetics, to localize and analyze L. angustifolius gene-rich regions (GRRs). A L. angustifolius genomic bacterial artificial chromosome (BAC) library was screened with short sequence repeat (SSR)-based probes. Selected BACs were fingerprinted and assembled into contigs. BAC-end sequence (BES) annotation allowed us to choose clones for sequencing, targeting GRRs. Additionally, BESs were aligned to the scaffolds of the genome sequence. The genetic map was supplemented with 35 BES-derived markers, distributed in 14 linkage groups and tagging 37 scaffolds. The identified GRRs had an average gene density of 19.6 genes/100 kb and physical-to-genetic distance ratios of 11 to 109 kb/cM. Physical and genetic mapping was supported by multi-BAC-fluorescence in situ hybridization (FISH), and five new linkage groups were assigned to the chromosomes. Syntenic links to the genome sequences of five legume species (Medicago truncatula, Glycine max, Lotus japonicus, Phaseolus vulgaris, and Cajanus cajan) were identified. The comparative mapping of the two largest lupin GRRs provides novel evidence for ancient duplications in all of the studied species. These regions are conserved among representatives of the main clades of Papilionoideae. Furthermore, despite the complex evolution of legumes, some segments of the nuclear genome were not substantially modified and retained their quasi-ancestral structures. Cytogenetic markers anchored in these regions constitute a platform for heterologous mapping of legume genomes.

Application of in silico bulked segregant analysis for rapid development of markers linked to Bean common mosaic virus resistance in common bean

BACKGROUND

Common bean was one of the first crops that benefited from the development and utilization of molecular marker-assisted selection (MAS) for major disease resistance genes. Efficiency of MAS for breeding common bean is still hampered, however, due to the dominance, linkage phase, and loose linkage of previously developed markers. Here we applied in silico bulked segregant analysis (BSA) to the BeanCAP diversity panel, composed of over 500 lines and genotyped with the BARCBEAN_3 6K SNP BeadChip, to develop codominant and tightly linked markers to the I gene controlling resistance to Bean common mosaic virus (BCMV).

RESULTS

We physically mapped the genomic region underlying the I gene. This locus, in the distal arm of chromosome Pv02, contains seven putative NBS-LRR-type disease resistance genes. Two contrasting bulks, containing BCMV host differentials and ten BeanCAP lines with known disease reaction to BCMV, were subjected to in silico BSA for targeting the I gene and flanking sequences. Two distinct haplotypes, containing a cluster of six single nucleotide polymorphisms (SNP), were associated with resistance or susceptibility to BCMV. One-hundred and twenty-two lines, including 115 of the BeanCAP panel, were screened for BCMV resistance in the greenhouse, and all of the resistant or susceptible plants displayed distinct SNP haplotypes as those found in the two bulks. The resistant/susceptible haplotypes were validated in 98 recombinant inbred lines segregating for BCMV resistance. The closest SNP (~25-32 kb) to the distal NBS-LRR gene model for the I gene locus was targeted for conversion to codominant KASP (Kompetitive Allele Specific PCR) and CAPS (Cleaved Amplified Polymorphic Sequence) markers. Both marker systems accurately predicted the disease reaction to BCMV conferred by the I gene in all screened lines of this study.

CONCLUSIONS

We demonstrated the utility of the in silico BSA approach using genetically diverse germplasm, genotyped with a high-density SNP chip array, to discover SNP variation at a specific targeted genomic region. In common bean, many disease resistance genes are mapped and their physical genomic position can now be determined, thus the application of this approach will facilitate further development of codominant and tightly linked markers for use in MAS.

Application of in silico bulked segregant analysis for rapid development of markers linked to Bean common mosaic virus resistance in common bean

Background

Common bean was one of the first crops that benefited from the development and utilization of molecular marker-assisted selection (MAS) for major disease resistance genes. Efficiency of MAS for breeding common bean is still hampered, however, due to the dominance, linkage phase, and loose linkage of previously developed markers. Here we applied in silico bulked segregant analysis (BSA) to the BeanCAP diversity panel, composed of over 500 lines and genotyped with the BARCBEAN_3 6K SNP BeadChip, to develop codominant and tightly linked markers to the I gene controlling resistance to Bean common mosaic virus (BCMV).

Results

We physically mapped the genomic region underlying the I gene. This locus, in the distal arm of chromosome Pv02, contains seven putative NBS-LRR-type disease resistance genes. Two contrasting bulks, containing BCMV host differentials and ten BeanCAP lines with known disease reaction to BCMV, were subjected to in silico BSA for targeting the I gene and flanking sequences. Two distinct haplotypes, containing a cluster of six single nucleotide polymorphisms (SNP), were associated with resistance or susceptibility to BCMV. One-hundred and twenty-two lines, including 115 of the BeanCAP panel, were screened for BCMV resistance in the greenhouse, and all of the resistant or susceptible plants displayed distinct SNP haplotypes as those found in the two bulks. The resistant/susceptible haplotypes were validated in 98 recombinant inbred lines segregating for BCMV resistance. The closest SNP (~25-32 kb) to the distal NBS-LRR gene model for the I gene locus was targeted for conversion to codominant KASP (Kompetitive Allele Specific PCR) and CAPS (Cleaved Amplified Polymorphic Sequence) markers. Both marker systems accurately predicted the disease reaction to BCMV conferred by the I gene in all screened lines of this study.

Conclusions

We demonstrated the utility of the in silico BSA approach using genetically diverse germplasm, genotyped with a high-density SNP chip array, to discover SNP variation at a specific targeted genomic region. In common bean, many disease resistance genes are mapped and their physical genomic position can now be determined, thus the application of this approach will facilitate further development of codominant and tightly linked markers for use in MAS.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-903) contains supplementary material, which is available to authorized users.

Genomics and molecular breeding in lesser explored pulse crops: current trends and future opportunities

Pulses are multipurpose crops for providing income, employment and food security in the underprivileged regions, notably the FAO-defined low-income food-deficit countries. Owing to their intrinsic ability to endure environmental adversities and the least input/management requirements, these crops remain central to subsistence farming. Given their pivotal role in rain-fed agriculture, substantial research has been invested to boost the productivity of these pulse crops. To this end, genomic tools and technologies have appeared as the compelling supplement to the conventional breeding. However, the progress in minor pulse crops including dry beans (Vigna spp.), lupins, lablab, lathyrus and vetches has remained unsatisfactory, hence these crops are often labeled as low profile or lesser researched. Nevertheless, recent scientific and technological breakthroughs particularly the next generation sequencing (NGS) are radically transforming the scenario of genomics and molecular breeding in these minor crops. NGS techniques have allowed de novo assembly of whole genomes in these orphan crops. Moreover, the availability of a reference genome sequence would promote re-sequencing of diverse genotypes to unlock allelic diversity at a genome-wide scale. In parallel, NGS has offered high-resolution genetic maps or more precisely, a robust genetic framework to implement whole-genome strategies for crop improvement. As has already been demonstrated in lupin, sequencing-based genotyping of the representative sample provided access to a number of functionally-relevant markers that could be deployed straight away in crop breeding programs. This article attempts to outline the recent progress made in genomics of these lesser explored pulse crops, and examines the prospects of genomics assisted integrated breeding to enhance and stabilize crop yields.

Draft genome sequence, and a sequence-defined genetic linkage map of the legume crop species Lupinus angustifolius L

Lupin (Lupinus angustifolius L.) is the most recently domesticated crop in major agricultural cultivation. Its seeds are high in protein and dietary fibre, but low in oil and starch. Medical and dietetic studies have shown that consuming lupin-enriched food has significant health benefits. We report the draft assembly from a whole genome shotgun sequencing dataset for this legume species with 26.9x coverage of the genome, which is predicted to contain 57,807 genes. Analysis of the annotated genes with metabolic pathways provided a partial understanding of some key features of lupin, such as the amino acid profile of storage proteins in seeds. Furthermore, we applied the NGS-based RAD-sequencing technology to obtain 8,244 sequence-defined markers for anchoring the genomic sequences. A total of 4,214 scaffolds from the genome sequence assembly were aligned into the genetic map. The combination of the draft assembly and a sequence-defined genetic map made it possible to locate and study functional genes of agronomic interest. The identification of co-segregating SNP markers, scaffold sequences and gene annotation facilitated the identification of a candidate R gene associated with resistance to the major lupin disease anthracnose. We demonstrated that the combination of medium-depth genome sequencing and a high-density genetic linkage map by application of NGS technology is a cost-effective approach to generating genome sequence data and a large number of molecular markers to study the genomics, genetics and functional genes of lupin, and to apply them to molecular plant breeding. This strategy does not require prior genome knowledge, which potentiates its application to a wide range of non-model species.

Stacks: an analysis tool set for population genomics

Massively parallel short-read sequencing technologies, coupled with powerful software platforms, are enabling investigators to analyse tens of thousands of genetic markers. This wealth of data is rapidly expanding and allowing biological questions to be addressed with unprecedented scope and precision. The sizes of the data sets are now posing significant data processing and analysis challenges. Here we describe an extension of the Stacks software package to efficiently use genotype-by-sequencing data for studies of populations of organisms. Stacks now produces core population genomic summary statistics and SNP-by-SNP statistical tests. These statistics can be analysed across a reference genome using a smoothed sliding window. Stacks also now provides several output formats for several commonly used downstream analysis packages. The expanded population genomics functions in Stacks will make it a useful tool to harness the newest generation of massively parallel genotyping data for ecological and evolutionary genetics.