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

The Diversity of Melia azedarach L. from China Based on Transcriptome-Developed SSR Marker

Melia azedarach L. is a native tree species that can be used in a comprehensive way and is widely distributed in all provinces south of the Yellow River in China. Genetic diversity analysis of different M. azedarach germplasm sources is an important basic work for the selection, evaluation, and genetic improvement of M. azedarach germplasm resources. In this study, 100 pairs of SSR primers were designed and synthesized based on M. azedarach transcriptome data, and 16 pairs of reliable SSR primers were finally selected. The developed primers were used to analyze the genetic diversity of M. azedarach from 15 sources in 10 provinces in East, Central, and South China. The results showed that the frequency of the M. azedarach transcriptome SSR loci was high, and the distribution density was high. There were 15 sources of M. azedarach genetic diversity at a moderate level, and genetic variation was mainly present within the sources. The present study further enriches the existing SSR marker database of the M. azedarach family and can provide a reference for genetic diversity analysis and molecularly assisted breeding of M. azedarach plants at the genomic level.

A successful defense of the narrow-leafed lupin against anthracnose involves quick and orchestrated reprogramming of oxidation-reduction, photosynthesis and pathogenesis-related genes

Narrow-leafed lupin (NLL, Lupinus angustifolius L.) is a legume plant cultivated for grain production and soil improvement. Worldwide expansion of NLL as a crop attracted various pathogenic fungi, including Colletotrichum lupini causing a devastating disease, anthracnose. Two alleles conferring improved resistance, Lanr1 and AnMan, were exploited in NLL breeding, however, underlying molecular mechanisms remained unknown. In this study, European NLL germplasm was screened with Lanr1 and AnMan markers. Inoculation tests in controlled environment confirmed effectiveness of both resistance donors. Representative resistant and susceptible lines were subjected to differential gene expression profiling. Resistance to anthracnose was associated with overrepresentation of "GO:0006952 defense response", "GO:0055114 oxidation-reduction process" and "GO:0015979 photosynthesis" gene ontology terms. Moreover, the Lanr1 (83A:476) line revealed massive transcriptomic reprogramming quickly after inoculation, whereas other lines showed such a response delayed by about 42 h. Defense response was associated with upregulation of TIR-NBS, CC-NBS-LRR and NBS-LRR genes, pathogenesis-related 10 proteins, lipid transfer proteins, glucan endo-1,3-beta-glucosidases, glycine-rich cell wall proteins and genes from reactive oxygen species pathway. Early response of 83A:476, including orchestrated downregulation of photosynthesis-related genes, coincided with the successful defense during fungus biotrophic growth phase, indicating effector-triggered immunity. Mandelup response was delayed and resembled general horizontal resistance.

Evaluation of the Major Seed Storage Proteins, the Conglutins, Across Genetically Diverse Narrow-Leafed Lupin Varieties

Lupin seeds have an excellent nutritional profile, including a high proportion of protein and dietary fiber. These qualities make lupin seeds an ideal candidate to help meet the growing global demand for complementary sources of protein. Of consequence to this application, there are nutritional and antinutritional properties assigned to the major lupin seed storage proteins—referred to as α-, β-, δ- and γ-conglutins The variation in the abundance of these protein families can impact the nutritional and bioactive properties of different lupin varieties. Hence, exploring the conglutin protein profiles across a diverse range of lupin varieties will yield knowledge that can facilitate the selection of superior genotypes for food applications or lupin crop improvement. To support this knowledge generation, discovery proteomics was applied for the identification of the 16 known conglutin subfamilies from 46 domestic and wild narrow-leafed lupin (NLL) genotypes. Consequently, the diversity of abundance of these proteins was evaluated using liquid chromatography–multiple reaction monitoring-mass spectrometry (LC–MRM-MS). This comparative study revealed a larger variability for the β- and δ-conglutin content across the lines under study. The absence/lower abundance of the β2- to β6-conglutin subfamilies in a subset of the domesticated cultivars led to substantially lower overall levels of the allergenic β-conglutin content in these NLLs, for which the elevation of the other conglutin families were observed. The diversity of the conglutin profiles revealed through this study—and the identification of potential hypoallergenic genotypes—will have great significance for lupin allergic consumers, food manufactures as well as grain breeders through the future development of lupin varieties with higher levels of desirable bioactive proteins and lower allergen content.

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.

Development and application of a virus-induced gene silencing protocol for the study of gene function in narrow-leafed lupin

BACKGROUND

Lupins are promising protein crops with an increasing amount of genomic and transcriptomic resources. The new resources facilitate the in silico identification of candidate genes controlling important agronomic traits. However, a major bottleneck for lupin research and crop improvement is the in planta characterization of gene function. Here, we present an efficient protocol for virus-induced gene silencing (VIGS) to down-regulate endogenous genes in narrow-leafed lupin (NLL) using the apple latent spherical virus (ALSV).

RESULTS

We identified ALSV as an appropriate VIGS vector able to infect NLL without causing a discernible phenotype. We created improved ALSV vectors to allow for efficient cloning of gene fragments into the viral genome and for easier viral propagation via agroinfiltration of Nicotiana benthamiana. Using this system, we silenced the visual marker gene phytoene desaturase (PDS), which resulted in systemic, homogenous silencing as indicated by bleaching of newly produced tissues. Furthermore, by silencing lysine decarboxylase (LaLDC)-a gene likely to be involved in toxic alkaloid biosynthesis-we demonstrate the applicability of our VIGS method to silence a target gene alone or alongside PDS in a 'PDS co-silencing' approach. The co-silencing approach allows the visual identification of tissues where silencing is actively occurring, which eases tissue harvesting and downstream analysis, and is useful where the trait under study is not affected by PDS silencing. Silencing LaLDC resulted in a ~ 61% or ~ 67% decrease in transcript level, depending on whether LaLDC was silenced alone or alongside PDS. Overall, the silencing of LaLDC resulted in reduced alkaloid levels, providing direct evidence of its involvement in alkaloid biosynthesis in NLL.

CONCLUSIONS

We provide a rapid and efficient VIGS method for validating gene function in NLL. This will accelerate the research and improvement of this underutilized crop.

Methyl jasmonate-induced resistance to Delia platura (Diptera: Anthomyiidae) in Lupinus mutabilis

BACKGROUND

Andean lupin (Lupinus mutabilis Sweet) is an important leguminous crop from South America with a high protein content. In Ecuador, lupin yields are severely affected by the infestation of Delia platura larvae on germinating seeds. The application of elicitor molecules with activity against herbivorous insects to control D. platura infestation constitutes an unexplored and promising alternative for chemical insecticides. In this study, methyl jasmonate (MeJA), hexanoic acid, menadione sodium bisulfite, and DL-β-aminobutyric acid were evaluated for their ability to induce resistance against D. platura in three commercial lupin cultivars.

RESULTS

Only seeds pretreated with MeJA significantly impaired insect performance during choice and no-choice assays. Additionally, fitness indicators such as seed germination and growth were not affected by MeJA treatment. To investigate the molecular mechanisms behind the MeJA-mediated resistance, RT-qPCR assays were performed. First, RT-qPCR reference genes were validated, showing that LmUBC was the most stable reference gene. Next, expression analysis over time revealed that MeJA application up-regulated the activity of the jasmonic acid biosynthetic genes LmLOX2 and LmAOS, together with other jasmonate-related defense genes, such as LmTPS1, LmTPS4, LmPI2, LmMBL, LmL/ODC, LmCSD1, and LmPOD.

CONCLUSION

This study indicates that MeJA can be used as an environmentally friendly elicitor molecule to protect Andean lupin from D. platura attack without fitness cost. MeJA application induces plant defense responses to insects in Andean lupin that may be modulated by the onset of terpenoid biosynthesis, proteinase inhibitors, lectins, polyamines, and antioxidative enzymes. © 2021 Society of Chemical Industry.

A Trimethylguanosine Synthase1-like (TGS1) homologue is implicated in vernalisation and flowering time control

KEY MESSAGE

A plant-specific Trimethylguanosine Synthase1-like homologue was identified as a candidate gene for the efl mutation in narrow-leafed lupin, which alters phenology by reducing vernalisation requirement. The vernalisation pathway is a key component of flowering time control in plants from temperate regions but is not well understood in the legume family. Here we examined vernalisation control in the temperate grain legume species, narrow-leafed lupin (Lupinus angustifolius L.), and discovered a candidate gene for an ethylene imine mutation (efl). The efl mutation changes phenology from late to mid-season flowering and additionally causes transformation from obligate to facultative vernalisation requirement. The efl locus was mapped to pseudochromosome NLL-10 in a recombinant inbred line (RIL) mapping population developed by accelerated single seed descent. Candidate genes were identified in the reference genome, and a diverse panel of narrow-leafed lupins was screened to validate mutations specific to accessions with efl. A non-synonymous SNP mutation within an S-adenosyl-L-methionine-dependent methyltransferase protein domain of a Trimethylguanosine Synthase1-like (TGS1) orthologue was identified as the candidate mutation giving rise to efl. This mutation caused substitution of an amino acid within an established motif at a position that is otherwise highly conserved in several plant families and was perfectly correlated with the efl phenotype in F2 and F6 genetic population and a panel of diverse accessions, including the original efl mutant. Expression of the TGS1 homologue did not differ between wild-type and efl genotypes, supporting altered functional activity of the gene product. This is the first time a TGS1 orthologue has been associated with vernalisation response and flowering time control in any plant species.

Evaluation of protein extraction methods for in-depth proteome analysis of narrow-leafed lupin (Lupinus angustifolius) seeds

Lupin is slated as a potential contributor towards future food security. Lupin possesses several nutritional and nutraceutical attributes, many linked to seed proteins. For in-depth characterisation of the lupin proteome, liquid chromatography-tandem mass spectrometry was used to evaluate four protein extraction procedures. The proteomes of three narrow-leafed lupin were qualitatively evaluated using protein/peptide identifications and further quantitatively assessed by data-independent proteome measurement. Each extraction buffer led to unique protein identifications; altogether yielding 2,760 protein identifications from lupin varieties. The analysis of protein abundance data highlighted distinct differences between Tris-HCl and urea extracted proteomes, while also revealing variation amongst the cultivar proteomes with the wild accession (P27255) distinctly different from the domesticated cultivars (Tanjil, Unicrop). The extraction buffer used influenced the proteome coverage, downstream functional annotation results and consequently the biological interpretation demonstrating the need to optimise and understand the impact of protein extraction conditions.

Alkaloids of narrow-leaved lupine as a factor determining alternative ways of the crop's utilization and breeding

Narrow-leaved lupine (Lupinus angustifolius L.), a valuable leguminous crop adapted to a wide range of climatic conditions, has a very short history of domestication. For many centuries it was used mainly as a green manure, since the success and prospects of the multi-purpose use of the species depend on its breeding improvement, in particular, on a particular concentration of alkaloids in seeds and green mass. The first varieties of scientific breeding were created only in the 1930s after the appearance of low-alkaloid mutants. Despite wide prospects for use in various areas of the national economy, unstable productivity and susceptibility to diseases hinder the production of this crop. Obviously, breeders deal only with a small part of the gene pool of the species and limited genetic resources, using mainly low-alkaloid (sweet) genotypes to create new varieties. The genetic potential of the species can be used more efficiently. At the same time, it is rational to create highly alkaloid (bitter) varieties for green manure, while food and feed varieties should not lose their adaptive potential, in particular, resistance to pathogens, due to the elimination of alkaloids. In this regard, it seems to be a productive idea to create ‘bitter/sweet’ varieties combining a high content of alkaloids in the vegetative organs and low in seeds, which can be achieved by regulating the synthesis/transport of alkaloids in the plant. The paper discusses the current state of use of the species as a green manure, fodder, food plant. Information is given on the quantity and qualitative composition of narrow-leaved lupine alkaloids, their applied value, in particular, fungicidal, antibacterial, insecticidal, the use of lupine alkaloids as active principles of drugs. Along with promising breeding considerations, the possibility of using technologies for processing raw high-alkaloid materials with the accompanying extraction of valuable ingredients for pharmaceuticals is discussed. Information is briefly presented about the genomic resources of the species and the prospects for their use in marker-assistant selection and genome editing.

Whole-genome assembly and resequencing reveal genomic imprint and key genes of rapid domestication in narrow-leafed lupin

Shifting from a livestock-based protein diet to a plant-based protein diet has been proposed as an essential requirement to maintain global food sustainability, which requires the increased production of protein-rich crops for direct human consumption. Meanwhile, the lack of sufficient genetic diversity in crop varieties is an increasing concern for sustainable food supplies. Countering this concern requires a clear understanding of the domestication process and dynamics. Narrow-leafed lupin (Lupinus angustifolius L.) has experienced rapid domestication and has become a new legume crop over the past century, with the potential to provide protein-rich seeds. Here, using long-read whole-genome sequencing, we assembled the third-generation reference genome for the narrow-leafed lupin cultivar Tanjil, comprising 20 chromosomes with a total genome size of 615.8 Mb and contig N50 = 5.65 Mb. We characterized the original mutation and putative biological pathway resulting in low seed alkaloid level that initiated the recent domestication of narrow-leafed lupin. We identified a 1133-bp insertion in the cis-regulatory region of a putative gene that may be associated with reduced pod shattering (lentus). A comparative analysis of genomic diversity in cultivars and wild types identified an apparent domestication bottleneck, as precisely predicted by the original model of the bottleneck effect on genetic variability in populations. Our results identify the key domestication genetic loci and provide direct genomic evidence for a domestication bottleneck, and open up the possibility of knowledge-driven de novo domestication of wild plants as an avenue to broaden crop plant diversity to enhance food security and sustainable low-carbon emission agriculture.

A functional genomics approach to dissect spotted alfalfa aphid resistance in Medicago truncatula

Aphids are virus-spreading insect pests affecting crops worldwide and their fast population build-up and insecticide resistance make them problematic to control. Here, we aim to understand the molecular basis of spotted alfalfa aphid (SAA) or Therioaphis trifolii f. maculata resistance in Medicago truncatula, a model organism for legume species. We compared susceptible and resistant near isogenic Medicago lines upon SAA feeding via transcriptome sequencing. Expression of genes involved in defense and stress responses, protein kinase activity and DNA binding were enriched in the resistant line. Potentially underlying some of these changes in gene expression was the finding that members of the MYB, NAC, AP2 domain and ERF transcription factor gene families were differentially expressed in the resistant versus susceptible lines. A TILLING population created in the resistant cultivar was screened using exome capture sequencing and served as a reverse genetics tool to functionally characterise genes involved in the aphid resistance response. This screening revealed three transcription factors (a NAC, AP2 domain and ERF) as important regulators in the defence response, as a premature stop-codon in the resistant background led to a delay in aphid mortality and enhanced plant susceptibility. This combined functional genomics approach will facilitate the future development of pest resistant crops by uncovering candidate target genes that can convey enhanced aphid resistance.

The Puzzling Fate of a Lupin Chromosome Revealed by Reciprocal Oligo-FISH and BAC-FISH Mapping

Old World lupins constitute an interesting model for evolutionary research due to diversity in genome size and chromosome number, indicating evolutionary genome reorganization. It has been hypothesized that the polyploidization event which occurred in the common ancestor of the Fabaceae family was followed by a lineage-specific whole genome triplication (WGT) in the lupin clade, driving chromosome rearrangements. In this study, chromosome-specific markers were used as probes for heterologous fluorescence in situ hybridization (FISH) to identify and characterize structural chromosome changes among the smooth-seeded (Lupinus angustifolius L., Lupinus cryptanthus Shuttlew., Lupinus micranthus Guss.) and rough-seeded (Lupinus cosentinii Guss. and Lupinus pilosus Murr.) lupin species. Comparative cytogenetic mapping was done using FISH with oligonucleotide probes and previously published chromosome-specific bacterial artificial chromosome (BAC) clones. Oligonucleotide probes were designed to cover both arms of chromosome Lang06 of the L. angustifolius reference genome separately. The chromosome was chosen for the in-depth study due to observed structural variability among wild lupin species revealed by BAC-FISH and supplemented by in silico mapping of recently released lupin genome assemblies. The results highlighted changes in synteny within the Lang06 region between the lupin species, including putative translocations, inversions, and/or non-allelic homologous recombination, which would have accompanied the evolution and speciation.

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.

Photoperiod and Vernalization Control of Flowering-Related Genes: A Case Study of the Narrow-Leafed Lupin (Lupinus angustifolius L.)

Narrow-leafed lupin (Lupinus angustifolius L.) is a moderate-yielding legume crop known for its high grain protein content and contribution to soil improvement. It is cultivated under photoperiods ranging from 9 to 17 h, as a spring-sown (in colder locations) or as an autumn-sown crop (in warmer regions). Wild populations require a prolonged cold period, called vernalization, to induce flowering. The key achievement of L. angustifolius domestication was the discovery of two natural mutations (named Ku and Jul) conferring vernalization independence. These mutations are overlapping deletion variants in the promoter of LanFTc1, a homolog of the Arabidopsis thaliana FLOWERING LOCUS T (FT) gene. The third deletion, named here as Pal, was recently found in primitive germplasm. In this study, we genotyped L. angustifolius germplasm that differs in domestication status and geographical origin for LanFTc1 alleles, which we then phenotyped to establish flowering time and vernalization responsiveness. The Ku and Jul lines were vernalization-independent and early flowering, wild (ku) lines were vernalization-dependent and late flowering, whereas the Pal line conferred intermediate phenotype. Three lines representing ku, Pal, and Ku alleles were subjected to gene expression surveys under 8- and 16-h photoperiods. FT homologs (LanFTa1, LanFTa2, LanFTc1, and LanFTc2) and some genes selected by recent expression quantitative trait loci mapping were analyzed. Expression profiles of LanFTc1 and LanAGL8 (AGAMOUS-like 8) matched observed differences in flowering time between genotypes, highlighted by high induction after vernalization in the ku line. Moreover, these genes revealed altered circadian clock control in Pal line under short days. LanFD (FD) and LanCRLK1 (CALCIUM/CALMODULIN-REGULATED RECEPTOR-LIKE KINASE 1) were negatively responsive to vernalization in Ku and Pal lines but positively responsive or variable in ku, whereas LanUGT85A2 (UDP-GLUCOSYL TRANSFERASE 85A2) was significantly suppressed by vernalization in all lines. Such a pattern suggests the opposite regulation of these gene pairs in the vernalization pathway. LanCRLK1 and LanUGT85A2 are homologs of A. thaliana genes involved in the FLOWERING LOCUS C (FLC) vernalization pathway. Lupins, like many other legumes, do not have any FLC homologs. Therefore, candidate genes surveyed in this study, namely LanFTc1, LanAGL8, LanCRLK1, and LanUGT85A2, may constitute anchors for further elucidation of molecular components contributing to vernalization response in legumes.

Proteomic Characterisation of Lupin (Lupinus angustifolius) Milk as Influenced by Extraction Techniques, Seed Coat and Cultivars

Lupin seeds are rich in proteins and other essential ingredients that can help to improve human health. The protein contents in both whole and split seeds of two lupin cultivars (Mandleup and PBA Jurien) were used to produce the lupin milk using the cheesecloth and centrifuge method. Proteins were extracted from the lupin milk using thiourea/urea solubilization. The proteins were separated by a two-dimensional polyacrylamide gel electrophoresis and then identified with mass spectrometry. A total of 230 protein spots were identified, 60 of which showed differential abundances. The cheesecloth separation showed protein extractability much better than that of the centrifuge method for both the cultivars. The results from this study could offer guidance for future comparative analysis and identification of lupin milk protein and provide effective separation technique to determine specific proteins in the cheese-making process.

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.

Photoperiod and Vernalization Control of Flowering-Related Genes: A Case Study of the Narrow-Leafed Lupin (Lupinus angustifolius L.)

Narrow-leafed lupin (Lupinus angustifolius L.) is a moderate-yielding legume crop known for its high grain protein content and contribution to soil improvement. It is cultivated under photoperiods ranging from 9 to 17 h, as a spring-sown (in colder locations) or as an autumn-sown crop (in warmer regions). Wild populations require a prolonged cold period, called vernalization, to induce flowering. The key achievement of L. angustifolius domestication was the discovery of two natural mutations (named Ku and Jul) conferring vernalization independence. These mutations are overlapping deletion variants in the promoter of LanFTc1, a homolog of the Arabidopsis thaliana FLOWERING LOCUS T (FT) gene. The third deletion, named here as Pal, was recently found in primitive germplasm. In this study, we genotyped L. angustifolius germplasm that differs in domestication status and geographical origin for LanFTc1 alleles, which we then phenotyped to establish flowering time and vernalization responsiveness. The Ku and Jul lines were vernalization-independent and early flowering, wild (ku) lines were vernalization-dependent and late flowering, whereas the Pal line conferred intermediate phenotype. Three lines representing ku, Pal, and Ku alleles were subjected to gene expression surveys under 8- and 16-h photoperiods. FT homologs (LanFTa1, LanFTa2, LanFTc1, and LanFTc2) and some genes selected by recent expression quantitative trait loci mapping were analyzed. Expression profiles of LanFTc1 and LanAGL8 (AGAMOUS-like 8) matched observed differences in flowering time between genotypes, highlighted by high induction after vernalization in the ku line. Moreover, these genes revealed altered circadian clock control in Pal line under short days. LanFD (FD) and LanCRLK1 (CALCIUM/CALMODULIN-REGULATED RECEPTOR-LIKE KINASE 1) were negatively responsive to vernalization in Ku and Pal lines but positively responsive or variable in ku, whereas LanUGT85A2 (UDP-GLUCOSYL TRANSFERASE 85A2) was significantly suppressed by vernalization in all lines. Such a pattern suggests the opposite regulation of these gene pairs in the vernalization pathway. LanCRLK1 and LanUGT85A2 are homologs of A. thaliana genes involved in the FLOWERING LOCUS C (FLC) vernalization pathway. Lupins, like many other legumes, do not have any FLC homologs. Therefore, candidate genes surveyed in this study, namely LanFTc1, LanAGL8, LanCRLK1, and LanUGT85A2, may constitute anchors for further elucidation of molecular components contributing to vernalization response in legumes.

Candidate Domestication-Related Genes Revealed by Expression Quantitative Trait Loci Mapping of Narrow-Leafed Lupin (Lupinus angustifolius L.)

The last century has witnessed rapid domestication of the narrow-leafed lupin (Lupinus angustifolius L.) as a grain legume crop, exploiting discovered alleles conferring low-alkaloid content (iucundus), vernalization independence (Ku and Julius), and reduced pod shattering (lentus and tardus). In this study, a L. angustifolius mapping population was subjected to massive analysis of cDNA ends (MACE). The MACE yielded 4185 single nucleotide polymorphism (SNP) markers for linkage map improvement and 30,595 transcriptomic profiles for expression quantitative trait loci (eQTL) mapping. The eQTL highlighted a high number of cis- and trans-regulated alkaloid biosynthesis genes with gene expression orchestrated by a regulatory agent localized at iucundus locus, supporting the concept that ETHYLENE RESPONSIVE TRANSCRIPTION FACTOR RAP2-7 may control low-alkaloid phenotype. The analysis of Ku shed light on the vernalization response via FLOWERING LOCUS T and FD regulon in L. angustifolius, providing transcriptomic evidence for the contribution of several genes acting in C-repeat binding factor (CBF) cold responsiveness and in UDP-glycosyltransferases pathways. Research on lentus selected a DUF1218 domain protein as a candidate gene controlling the orientation of the sclerified endocarp and a homolog of DETOXIFICATION14 for purplish hue of young pods. An ABCG transporter was identified as a hypothetical contributor to sclerenchyma fortification underlying tardus phenotype.

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%).

The first genetic map for yellow lupin enables genetic dissection of adaptation traits in an orphan grain legume crop

BACKGROUND

Yellow lupin (Lupinus luteus L.) is a promising grain legume for productive and sustainable crop rotations. It has the advantages of high tolerance to soil acidity and excellent seed quality, but its current yield potential is poor, especially in low rainfall environments. Key adaptation traits such as phenology and enhanced stress tolerance are often complex and controlled by several genes. Genomic-enabled technologies may help to improve our basic understanding of these traits and to provide selective markers in breeding. However, in yellow lupin there are very limited genomic resources to support research and no published information is available on the genetic control of adaptation traits.

RESULTS

We aimed to address these deficiencies by developing the first linkage map for yellow lupin and conducting quantitative trait locus (QTL) analysis of yield under well-watered (WW) and water-deficit (WT) conditions. Two next-generation sequencing marker approaches - genotyping-by-sequencing (GBS) and Diversity Array Technology (DArT) sequencing - were employed to genotype a recombinant inbred line (RIL) population developed from a bi-parental cross between wild and domesticated parents. A total of 2,458 filtered single nucleotide polymorphism (SNP) and presence / absence variation (PAV) markers were used to develop a genetic map comprising 40 linkage groups, the first reported for this species. A number of significant QTLs controlling total biomass and 100-seed weight under two water (WW and WD) regimes were found on linkage groups YL-03, YL-09 and YL-26 that together explained 9 and 28% of total phenotypic variability. QTLs associated with length of the reproductive phase and time to flower were found on YL-01, YL-21, YL-35 and YL-40 that together explained a total of 12 and 44% of total phenotypic variation.

CONCLUSION

These genomic resources and the QTL information offer significant potential for use in marker-assisted selection in yellow lupin.

Development of gene-based molecular markers tagging low alkaloid pauper locus in white lupin (Lupinus albus L.)

White lupin (Lupinus albus L.) is a legume grain crop cultivated since ancient Greece and Egypt. Modern white lupin cultivars are appreciated as a source of protein with positive nutraceutical impact. However, white lupins produce anti-nutritional compounds, quinolizidine alkaloids, which provide bitter taste and have a negative influence on human health. During domestication of this species, several recessive alleles at unlinked loci controlling low alkaloid content were selected. One of these loci, pauper, was exploited worldwide providing numerous low-alkaloid cultivars. However, molecular tracking of pauper has been hampered due to the lack of diagnostic markers. In the present study, the synteny-based approach was harnessed to target pauper locus. Single-nucleotide polymorphisms flanking pauper locus on white lupin linkage map as well as candidate gene sequences elucidated from the narrow-leafed lupin (L. angustifolius L.) chromosome segment syntenic to the pauper linkage group region were transformed to PCR-based molecular markers. These markers were analyzed both in the mapping population and world germplasm collection. From fourteen markers screened, eleven were localized at a distance below 1.5 cM from this locus, including five co-segregating with pauper. The linkage of these markers was confirmed by high LOD values (up to 58.4). Validation performed in the set of 127 bitter and 23 sweet accessions evidenced high applicability of one marker, LAGI01_35805_F1_R1, for pauper locus selection, highlighted by the low ratio of false-positive scores (2.5%). LAGI01_35805 represents a homolog of L. angustifolius acyltransferase-like (LaAT) gene which might hypothetically participate in the alkaloid biosynthesis process in lupins.

The role of jasmonate signalling in quinolizidine alkaloid biosynthesis, wounding and aphid predation response in narrow-leafed lupin

Quinolizidine alkaloids (QAs) are toxic secondary metabolites produced in lupin species that protect the plant against insects. They form in vegetative tissues and accumulate to a different extent in the grains: high levels in 'bitter' narrow-leafed lupin (NLL) and low levels in 'sweet' NLL. Grain QA levels vary considerably, and sometimes exceed the industry limit for food and feed purposes. We hypothesised that jasmonates regulate QA biosynthesis in response to environmental stresses such as wounding and aphid predation, which may explain non-genetic variability in grain QA levels. Methyl jasmonate (MeJA)-inducible genes were identified and verified in NLL. Exogenous MeJA application-induced expression of QA biosynthetic genes and QA levels for bitter, but not sweet NLL. Although MeJA-inducible genes responded to wounding, the expression of QA biosynthetic genes was not induced for bitter and sweet NLL. We assessed the effect of aphid predation on QA production for two cultivars - one moderately resistant and one susceptible to aphid predation. Although MeJA-inducible genes responded to aphid predation, no change in QA levels was found for either cultivar. These findings offer insights into the regulation of QA biosynthesis in bitter and sweet NLL and concludes that aphids are not a concern for increasing grain QAs in NLL cultivars.

Transcriptome-derived investigation of biosynthesis of quinolizidine alkaloids in narrow-leafed lupin (Lupinus angustifolius L.) highlights candidate genes linked to iucundus locus

Unravelling the biosynthetic pathway of quinolizidine alkaloids (QAs), regarded as antinutritional compounds of narrow-leafed lupin (NLL) seeds, is fundamental to best exploit NLL as food or feed. We investigated 12 candidate genes connected to QA biosynthesis, selecting them by transcriptomic and genomic approaches, from the landscape of genes differentially expressed in leaves of the high- and low-alkaloid NLL accessions. Linkage analysis enabled the assessment of the location of the candidate genes in relation to iucundus, a major locus of unknown identity, that confers reduced QA content in seeds. The key finding was the identification of APETALA2/ethylene response transcription factor, RAP2-7, cosegregating with the iucundus locus and located within a region with highly significant QTLs that affect QA composition. We additionally identified a 4-hydroxy-tetrahydrodipicolinate synthase (DHDPS) gene involved in L-lysine biosynthesis as being closely linked to iucundus. The distributed location of other remaining candidates (including previously known QA genes) across different linkage groups, also indirectly supports the transcription factor as a possible regulator of lupin alkaloid biosynthesis. Our findings provide crucial insight into QA biosynthesis in NLL. Additionally, we evaluated and selected appropriate reference genes for qRT-PCRs to analyse the expression levels of QA genes in NLL.

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.

Identification and profiling of narrow-leafed lupin (Lupinus angustifolius) microRNAs during seed development

BACKGROUND

Whilst information regarding small RNAs within agricultural crops is increasing, the miRNA composition of the nutritionally valuable pulse narrow-leafed lupin (Lupinus angustifolius) remains unknown.

RESULTS

By conducting a genome- and transcriptome-wide survey we identified 7 Dicer-like and 16 Argonaute narrow-leafed lupin genes, which were highly homologous to their legume counterparts. We identified 43 conserved miRNAs belonging to 16 families, and 13 novel narrow-leafed lupin-specific miRNAs using high-throughput sequencing of small RNAs from foliar and root and five seed development stages. We observed up-regulation of members of the miRNA families miR167, miR399, miR156, miR319 and miR164 in narrow-leafed lupin seeds, and confirmed expression of miR156, miR166, miR164, miR1507 and miR396 using quantitative RT-PCR during five narrow-leafed lupin seed development stages. We identified potential targets for the conserved and novel miRNAs and were able to validate targets of miR399 and miR159 using 5' RLM-RACE. The conserved miRNAs are predicted to predominately target transcription factors and 93% of the conserved miRNAs originate from intergenic regions. In contrast, only 43% of the novel miRNAs originate from intergenic regions and their predicted targets were more functionally diverse.

CONCLUSION

This study provides important insights into the miRNA gene regulatory networks during narrow-leafed lupin seed development.

INDEL variation in the regulatory region of the major flowering time gene LanFTc1 is associated with vernalization response and flowering time in narrow-leafed lupin (Lupinus angustifolius L.)

Narrow-leafed lupin (Lupinus angustifolius L.) cultivation was transformed by 2 dominant vernalization-insensitive, early flowering time loci known as Ku and Julius (Jul), which allowed expansion into shorter season environments. However, reliance on these loci has limited genetic and phenotypic diversity for environmental adaptation in cultivated lupin. We recently predicted that a 1,423-bp deletion in the cis-regulatory region of LanFTc1, a FLOWERING LOCUS T (FT) homologue, derepressed expression of LanFTc1 and was the underlying cause of the Ku phenotype. Here, we surveyed diverse germplasm for LanFTc1 cis-regulatory variation and identified 2 further deletions of 1,208 and 5,162 bp in the 5' regulatory region, which overlap the 1,423-bp deletion. Additionally, we confirmed that no other polymorphisms were perfectly associated with vernalization responsiveness. Phenotyping and gene expression analyses revealed that Jul accessions possessed the 5,162-bp deletion and that the Jul and Ku deletions were equally capable of removing vernalization requirement and up-regulating gene expression. The 1,208-bp deletion was associated with intermediate phenology, vernalization responsiveness, and gene expression and therefore may be useful for expanding agronomic adaptation of lupin. This insertion/deletion series may also help resolve how the vernalization response is mediated at the molecular level in legumes.

Adapting legume crops to climate change using genomic approaches

Our agricultural system and hence food security is threatened by combination of events, such as increasing population, the impacts of climate change, and the need to a more sustainable development. Evolutionary adaptation may help some species to overcome environmental changes through new selection pressures driven by climate change. However, success of evolutionary adaptation is dependent on various factors, one of which is the extent of genetic variation available within species. Genomic approaches provide an exceptional opportunity to identify genetic variation that can be employed in crop improvement programs. In this review, we illustrate some of the routinely used genomics-based methods as well as recent breakthroughs, which facilitate assessment of genetic variation and discovery of adaptive genes in legumes. Although additional information is needed, the current utility of selection tools indicate a robust ability to utilize existing variation among legumes to address the challenges of climate uncertainty.

Legume Cytosolic and Plastid Acetyl-Coenzyme-A Carboxylase Genes Differ by Evolutionary Patterns and Selection Pressure Schemes Acting before and after Whole-Genome Duplications

Acetyl-coenzyme A carboxylase (ACCase, E.C.6.4.1.2) catalyzes acetyl-coenzyme A carboxylation to malonyl coenzyme A. Plants possess two distinct ACCases differing by cellular compartment and function. Plastid ACCase contributes to de novo fatty acid synthesis, whereas cytosolic enzyme to the synthesis of very long chain fatty acids, phytoalexins, flavonoids, and anthocyanins. The narrow leafed lupin (Lupinus angustifolius L.) represents legumes, a plant family which evolved by whole-genome duplications (WGDs). The study aimed on the contribution of these WGDs to the multiplication of ACCase genes and their further evolutionary patterns. The molecular approach involved bacterial artificial chromosome (BAC) library screening, fluorescent in situ hybridization, linkage mapping, and BAC sequencing. In silico analysis encompassed sequence annotation, comparative mapping, selection pressure calculation, phylogenetic inference, and gene expression profiling. Among sequenced legumes, the highest number of ACCase genes was identified in lupin and soybean. The most abundant plastid ACCase subunit genes were accB. ACCase genes in legumes evolved by WGDs, evidenced by shared synteny and Bayesian phylogenetic inference. Transcriptional activity of almost all copies was confirmed. Gene duplicates were conserved by strong purifying selection, however, positive selection occurred in Arachis (accB2) and Lupinus (accC) lineages, putatively predating the WGD event(s). Early duplicated accA and accB genes underwent transcriptional sub-functionalization.

Characterization of the genetic factors affecting quinolizidine alkaloid biosynthesis and its response to abiotic stress in narrow-leafed lupin (Lupinus angustifolius L.)

Quinolizidine alkaloids (QAs) are toxic secondary metabolites that complicate the end use of narrow-leafed lupin (NLL; Lupinus angustifolius L.) grain, as levels sometimes exceed the industry limit for its use as a food and feed source. The genotypic and environmental influences on QA production in NLL are poorly understood. Here, the expression of QA biosynthetic genes was analysed in vegetative and reproductive tissues of bitter (high QA) and sweet (low QA) accessions. It was demonstrated that sweet accessions are characterized by lower QA biosynthetic gene expression exclusively in leaf and stem tissues than bitter NLL, consistent with the hypothesis that QAs are predominantly produced in aerial tissues and transported to seeds, rather than synthesized within the seed itself. This analysis informed our identification of additional candidate genes involved in QA biosynthesis. Drought and temperature stress are two major abiotic stresses that often occur during NLL pod set. Hence, we assessed the effect of drought, increased temperature, and their combination, on QA production in three sweet NLL cultivars. A cultivar-specific response to drought and temperature in grain QA levels was observed, including the identification of a cultivar where alkaloid levels did not change with these stress treatments.

Construction of an ultra-high density consensus genetic map, and enhancement of the physical map from genome sequencing in Lupinus angustifolius

An ultra-high density genetic map containing 34,574 sequence-defined markers was developed in Lupinus angustifolius. Markers closely linked to nine genes of agronomic traits were identified. A physical map was improved to cover 560.5 Mb genome sequence. Lupin (Lupinus angustifolius L.) is a recently domesticated legume grain crop. In this study, we applied the restriction-site associated DNA sequencing (RADseq) method to genotype an F₉ recombinant inbred line population derived from a wild type × domesticated cultivar (W × D) cross. A high density linkage map was developed based on the W × D population. By integrating sequence-defined DNA markers reported in previous mapping studies, we established an ultra-high density consensus genetic map, which contains 34,574 markers consisting of 3508 loci covering 2399 cM on 20 linkage groups. The largest gap in the entire consensus map was 4.73 cM. The high density W × D map and the consensus map were used to develop an improved physical map, which covered 560.5 Mb of genome sequence data. The ultra-high density consensus linkage map, the improved physical map and the markers linked to genes of breeding interest reported in this study provide a common tool for genome sequence assembly, structural genomics, comparative genomics, functional genomics, QTL mapping, and molecular plant breeding in lupin.

Transcript profiling of a bitter variety of narrow-leafed lupin to discover alkaloid biosynthetic genes

Lupins (Lupinus spp.) are nitrogen-fixing legumes that accumulate toxic alkaloids in their protein-rich beans. These anti-nutritional compounds belong to the family of quinolizidine alkaloids (QAs), which are of interest to the pharmaceutical and chemical industries. To unleash the potential of lupins as protein crops and as sources of QAs, a thorough understanding of the QA pathway is needed. However, only the first enzyme in the pathway, lysine decarboxylase (LDC), is known. Here, we report the transcriptome of a high-QA variety of narrow-leafed lupin (L. angustifolius), obtained using eight different tissues and two different sequencing technologies. In addition, we present a list of 33 genes that are closely co-expressed with LDC and that represent strong candidates for involvement in lupin alkaloid biosynthesis. One of these genes encodes a copper amine oxidase able to convert the product of LDC, cadaverine, into 1-piperideine, as shown by heterologous expression and enzyme assays. Kinetic analysis revealed a low KM value for cadaverine, supporting a role as the second enzyme in the QA pathway. Our transcriptomic data set represents a crucial step towards the discovery of enzymes, transporters, and regulators involved in lupin alkaloid biosynthesis.

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.

Legume isoflavone synthase genes have evolved by whole-genome and local duplications yielding transcriptionally active paralogs

Isoflavone synthase (IFS) is the key enzyme of isoflavonoid biosynthesis. IFS genes were identified in numerous species, although their evolutionary patterns have not yet been reconstructed. To address this issue, we performed structural and functional genomic analysis. Narrow leafed lupin, Lupinus angustifolius L., was used as a reference species for the genus, because it has the most developed molecular tools available. Nuclear genome BAC library clones carrying IFS homologs were localized by linkage mapping and fluorescence in situ hybridization in three chromosome pairs. Annotation of BAC, scaffold and transcriptome sequences confirmed the presence of three full-length IFS genes in the genome. Microsynteny analysis and Bayesian inference provided clear evidence that IFS genes in legumes have evolved by lineage-specific whole-genome and tandem duplications. Gene expression profiling and RNA-seq data mining showed that the vast majority of legume IFS copies have maintained their transcriptional activity. L. angustifolius IFS homologs exhibited organ-specific expression patterns similar to those observed in other Papilionoideae. Duplicated lupin IFS homologs retained non-negligible levels of substitutions in conserved motifs, putatively due to positive selection acting during early evolution of the genus, before the whole-genome duplication. Strong purifying selection preserved newly arisen IFS duplicates from further nonsynonymous changes.

Toward Genomics-Based Breeding in C3 Cool-Season Perennial Grasses

Most important food and feed crops in the world belong to the C3 grass family. The future of food security is highly reliant on achieving genetic gains of those grasses. Conventional breeding methods have already reached a plateau for improving major crops. Genomics tools and resources have opened an avenue to explore genome-wide variability and make use of the variation for enhancing genetic gains in breeding programs. Major C3 annual cereal breeding programs are well equipped with genomic tools; however, genomic research of C3 cool-season perennial grasses is lagging behind. In this review, we discuss the currently available genomics tools and approaches useful for C3 cool-season perennial grass breeding. Along with a general review, we emphasize the discussion focusing on forage grasses that were considered orphan and have little or no genetic information available. Transcriptome sequencing and genotype-by-sequencing technology for genome-wide marker detection using next-generation sequencing (NGS) are very promising as genomics tools. Most C3 cool-season perennial grass members have no prior genetic information; thus NGS technology will enhance collinear study with other C3 model grasses like Brachypodium and rice. Transcriptomics data can be used for identification of functional genes and molecular markers, i.e., polymorphism markers and simple sequence repeats (SSRs). Genome-wide association study with NGS-based markers will facilitate marker identification for marker-assisted selection. With limited genetic information, genomic selection holds great promise to breeders for attaining maximum genetic gain of the cool-season C3 perennial grasses. Application of all these tools can ensure better genetic gains, reduce length of selection cycles, and facilitate cultivar development to meet the future demand for food and fodder.

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.

Quinolizidine Alkaloid Biosynthesis in Lupins and Prospects for Grain Quality Improvement

Quinolizidine alkaloids (QAs) are toxic secondary metabolites found within the genus Lupinus, some species of which are commercially important grain legume crops including Lupinus angustifolius (narrow-leafed lupin, NLL), L. luteus (yellow lupin), L. albus (white lupin), and L. mutabilis (pearl lupin), with NLL grain being the most largely produced of the four species in Australia and worldwide. While QAs offer the plants protection against insect pests, the accumulation of QAs in lupin grain complicates its use for food purposes as QA levels must remain below the industry threshold (0.02%), which is often exceeded. It is not well understood what factors cause grain QA levels to exceed this threshold. Much of the early work on QA biosynthesis began in the 1970-1980s, with many QA chemical structures well-characterized and lupin cell cultures and enzyme assays employed to identify some biosynthetic enzymes and pathway intermediates. More recently, two genes associated with these enzymes have been characterized, however, the QA biosynthetic pathway remains only partially elucidated. Here, we review the research accomplished thus far concerning QAs in lupin and consider some possibilities for further elucidation and manipulation of the QA pathway in lupin crops, drawing on examples from model alkaloid species. One breeding strategy for lupin is to produce plants with high QAs in vegetative tissues while low in the grain in order to confer insect resistance to plants while keeping grain QA levels within industry regulations. With the knowledge achieved on alkaloid biosynthesis in other plant species in recent years, and the recent development of genomic and transcriptomic resources for NLL, there is considerable scope to facilitate advances in our knowledge of QAs, leading to the production of improved lupin crops.

The loss of vernalization requirement in narrow-leafed lupin is associated with a deletion in the promoter and de-repressed expression of a Flowering Locus T (FT) homologue

Adaptation of Lupinus angustifolius (narrow-leafed lupin) to cropping in southern Australian and northern Europe was transformed by a dominant mutation (Ku) that removed vernalization requirement for flowering. The Ku mutation is now widely used in lupin breeding to confer early flowering and maturity. We report here the identity of the Ku mutation. We used a range of genetic, genomic and gene expression approaches to determine whether Flowering Locus T (FT) homologues are associated with the Ku locus. One of four FT homologues present in the narrow-leafed lupin genome, LanFTc1, perfectly co-segregated with the Ku locus in a reference mapping population. Expression of LanFTc1 in the ku (late-flowering) parent was strongly induced by vernalization, in contrast to the Ku (early-flowering) parent, which showed constitutively high LanFTc1 expression. Co-segregation of this expression phenotype with the LanFTc1 genotype indicated that the Ku mutation impairs cis-regulation of LanFTc1. Sequencing of LanFTc1 revealed a 1.4-kb deletion in the promoter region, which was perfectly predictive of vernalization response in 216 wild and domesticated accessions. Linkage disequilibrium rapidly decayed around LanFTc1, suggesting that this deletion caused the loss of vernalization response. This is the first time a legume FTc subclade gene has been implicated in the vernalization response.

Narrow-Leafed Lupin (Lupinus angustifolius) β1- and β6-Conglutin Proteins Exhibit Antifungal Activity, Protecting Plants against Necrotrophic Pathogen Induced Damage from Sclerotinia sclerotiorum and Phytophthora nicotianae

Vicilins (7S globulins) are seed storage proteins and constitute the main protein family in legume seeds, particularly in narrow-leafed lupin (Lupinus angustifolius L.; NLL), where seven vicilin genes, called β1- to β7-conglutin have been identified. Vicilins are involved in germination processes supplying amino acids for seedling growth and plant development, as well as in some cases roles in plant defense and protection against pathogens. The roles of NLL β-conglutins in plant defense are unknown. Here the potential role of five NLL β-conglutin family members in protection against necrotrophic fungal pathogens was investigated and it was demonstrated that recombinant purified 6xHis-tagged β1- and β6-conglutin proteins exhibited the strongest in vitro growth inhibitory activity against a range of necrotrophic fungal pathogens compared to β2, β3, and β4 conglutins. To examine activity in vivo, two representative necrotrophic pathogens, the fungus Sclerotinia sclerotiorum and oomycete Phytophthora nicotianae were used. Transient expression of β1- and β6-conglutin proteins in Nicotiana benthamiana leaves demonstrated in vivo growth suppression of both of these pathogens, resulting in low percentages of hyphal growth and elongation in comparison to control treated leaves. Cellular studies using β1- and β6-GFP fusion proteins showed these conglutins localized to the cell surface including plasmodesmata. Analysis of cellular death following S. sclerotiorum or P. nicotianae revealed both β1- and β6-conglutins suppressed pathogen induced cell death in planta and prevented pathogen induced suppression of the plant oxidative burst as determined by protein oxidation in infected compared to mock-inoculated leaves.

Expansion of the phosphatidylethanolamine binding protein family in legumes: a case study of Lupinus angustifolius L. FLOWERING LOCUS T homologs, LanFTc1 and LanFTc2

BACKGROUND

The Arabidopsis FLOWERING LOCUS T (FT) gene, a member of the phosphatidylethanolamine binding protein (PEBP) family, is a major controller of flowering in response to photoperiod, vernalization and light quality. In legumes, FT evolved into three, functionally diversified clades, FTa, FTb and FTc. A milestone achievement in narrow-leafed lupin (Lupinus angustifolius L.) domestication was the loss of vernalization responsiveness at the Ku locus. Recently, one of two existing L. angustifolius homologs of FTc, LanFTc1, was revealed to be the gene underlying Ku. It is the first recorded involvement of an FTc homologue in vernalization. The evolutionary basis of this phenomenon in lupin has not yet been deciphered.

RESULTS

Bacterial artificial chromosome (BAC) clones carrying LanFTc1 and LanFTc2 genes were localized in different mitotic chromosomes and constituted sequence-specific landmarks for linkage groups NLL-10 and NLL-17. BAC-derived superscaffolds containing LanFTc genes revealed clear microsyntenic patterns to genome sequences of nine legume species. Superscaffold-1 carrying LanFTc1 aligned to regions encoding one or more FT-like genes whereas superscaffold-2 mapped to a region lacking such a homolog. Comparative mapping of the L. angustifolius genome assembly anchored to linkage map localized superscaffold-1 in the middle of a 15 cM conserved, collinear region. In contrast, superscaffold-2 was found at the edge of a 20 cM syntenic block containing highly disrupted collinearity at the LanFTc2 locus. 118 PEBP-family full-length homologs were identified in 10 legume genomes. Bayesian phylogenetic inference provided novel evidence supporting the hypothesis that whole-genome and tandem duplications contributed to expansion of PEBP-family genes in legumes. Duplicated genes were subjected to strong purifying selection. Promoter analysis of FT genes revealed no statistically significant sequence similarity between duplicated copies; only RE-alpha and CCAAT-box motifs were found at conserved positions and orientations.

CONCLUSIONS

Numerous lineage-specific duplications occurred during the evolution of legume PEBP-family genes. Whole-genome duplications resulted in the origin of subclades FTa, FTb and FTc and in the multiplication of FTa and FTb copy number. LanFTc1 is located in the region conserved among all main lineages of Papilionoideae. LanFTc1 is a direct descendant of ancestral FTc, whereas LanFTc2 appeared by subsequent duplication.

Root trait diversity, molecular marker diversity, and trait-marker associations in a core collection of Lupinus angustifolius

Narrow-leafed lupin (Lupinus angustifolius L.) is the predominant grain legume crop in southern Australia, contributing half of the total grain legume production of Australia. Its yield in Australia is hampered by a range of subsoil constraints. The adaptation of lupin genotypes to subsoil constraints may be improved by selecting for optimal root traits from new and exotic germplasm sources. We assessed root trait diversity and genetic diversity of a core collection of narrow-leafed lupin (111 accessions) using 191 Diversity Arrays Technology (DArT) markers. The genetic relationship among accessions was determined using the admixture model in STRUCTURE. Thirty-eight root-associated traits were characterized, with 21 having coefficient of variation values >0.5. Principal coordinate analysis and cluster analysis of the DArT markers revealed broad diversity among the accessions. An ad hoc statistics calculation resulted in 10 distinct populations with significant differences among and within them (P < 0.001). The mixed linear model test in TASSEL showed a significant association between all root traits and some DArT markers, with the numbers of markers associated with an individual trait ranging from 2 to 13. The percentage of phenotypic variation explained by any one marker ranged from 6.4 to 21.8%, with 15 associations explaining >10% of phenotypic variation. The genetic variation values ranged from 0 to 7994, with 23 associations having values >240. Root traits such as deeper roots and lateral root proliferation at depth would be useful for this species for improved adaptation to drier soil conditions. This study offers opportunities for discovering useful root traits that can be used to increase the yield of Australian cultivars across variable environmental conditions.

Integration of Lupinus angustifolius L. (narrow-leafed lupin) genome maps and comparative mapping within legumes

Narrow-leafed lupin (Lupinus angustifolius L.) has recently been considered a reference genome for the Lupinus genus. In the present work, genetic and cytogenetic maps of L. angustifolius were supplemented with 30 new molecular markers representing lupin genome regions, harboring genes involved in nitrogen fixation during the symbiotic interaction of legumes and soil bacteria (Rhizobiaceae). Our studies resulted in the precise localization of bacterial artificial chromosomes (BACs) carrying sequence variants for early nodulin 40, nodulin 26, nodulin 45, aspartate aminotransferase P2, asparagine synthetase, cytosolic glutamine synthetase, and phosphoenolpyruvate carboxylase. Together with previously mapped chromosomes, the integrated L. angustifolius map encompasses 73 chromosome markers, including 5S ribosomal DNA (rDNA) and 45S rDNA, and anchors 20 L. angustifolius linkage groups to corresponding chromosomes. Chromosomal identification using BAC fluorescence in situ hybridization identified two BAC clones as narrow-leafed lupin centromere-specific markers, which served as templates for preliminary studies of centromere composition within the genus. Bioinformatic analysis of these two BACs revealed that centromeric/pericentromeric regions of narrow-leafed lupin chromosomes consisted of simple sequence repeats ordered into tandem repeats containing the trinucleotide and pentanucleotide simple sequence repeats AGG and GATAC, structured into long arrays. Moreover, cross-genus microsynteny analysis revealed syntenic patterns of 31 single-locus BAC clones among several legume species. The gene and chromosome level findings provide evidence of ancient duplication events that must have occurred very early in the divergence of papilionoid lineages. This work provides a strong foundation for future comparative mapping among legumes and may facilitate understanding of mechanisms involved in shaping legume chromosomes.

Proteomics for exploiting diversity of lupin seed storage proteins and their use as nutraceuticals for health and welfare

Lupins have a variety of both traditional and modern uses. In the last decade, reports assessing the benefits of lupin seed proteins have proliferated and, nowadays, the pharmaceutical industry is interested in lupin proteins for human health. Modern genomics and proteomics have hugely contributed to describing the diversity of lupin storage genes and, above all, proteins. Most of these studies have been centered on few edible lupin species. However, Lupinus genus comprises hundreds of species spread throughout the Old and New Worlds, and these resources have been scarcely explored and exploited. We present here a detailed review of the literature on the potential of lupin seed proteins as nutraceuticals, and the use of -omic tools to analyze seed storage polypeptides in main edible lupins and their diversity at the Lupinus inter- and intra-species level. In this sense, proteomics, more than any other, has been a key approach. Proteomics has shown that lupin seed protein diversity, where post-translational modifications yield a large number of peptide variants with a potential concern in bioactivity, goes far beyond gene diversity. The future extended use of second and third generation proteomics should definitely help to go deeper into coverage and characterization of lupin seed proteome.

BIOLOGICAL SIGNIFICANCE

Some important topics concerning storage proteins from lupin seeds are presented and analyzed in an integrated way in this review. Proteomic approaches have been essential in characterizing lupin seed protein diversity, which goes far beyond gene diversity since the protein level adds to the latter differential proteolytic cleavage of conglutin pro-proteins and a diverse array of glycosylation forms and sites. Proteomics has also proved helpful for screening and studying Lupinus germplasm with the future aim of exploiting and improving food production, quality, and nutritional values.

Identifying Stable Reference Genes for qRT-PCR Normalisation in Gene Expression Studies of Narrow-Leafed Lupin (Lupinus angustifolius L.)

Quantitative Reverse Transcription PCR (qRT-PCR) is currently one of the most popular, high-throughput and sensitive technologies available for quantifying gene expression. Its accurate application depends heavily upon normalisation of gene-of-interest data with reference genes that are uniformly expressed under experimental conditions. The aim of this study was to provide the first validation of reference genes for Lupinus angustifolius (narrow-leafed lupin, a significant grain legume crop) using a selection of seven genes previously trialed as reference genes for the model legume, Medicago truncatula. In a preliminary evaluation, the seven candidate reference genes were assessed on the basis of primer specificity for their respective targeted region, PCR amplification efficiency, and ability to discriminate between cDNA and gDNA. Following this assessment, expression of the three most promising candidates [Ubiquitin C (UBC), Helicase (HEL), and Polypyrimidine tract-binding protein (PTB)] was evaluated using the NormFinder and RefFinder statistical algorithms in two narrow-leafed lupin lines, both with and without vernalisation treatment, and across seven organ types (cotyledons, stem, leaves, shoot apical meristem, flowers, pods and roots) encompassing three developmental stages. UBC was consistently identified as the most stable candidate and has sufficiently uniform expression that it may be used as a sole reference gene under the experimental conditions tested here. However, as organ type and developmental stage were associated with greater variability in relative expression, it is recommended using UBC and HEL as a pair to achieve optimal normalisation. These results highlight the importance of rigorously assessing candidate reference genes for each species across a diverse range of organs and developmental stages. With emerging technologies, such as RNAseq, and the completion of valuable transcriptome data sets, it is possible that other potentially more suitable reference genes will be identified for this species in future.

Characterization and mapping of LanrBo: a locus conferring anthracnose resistance in narrow-leafed lupin (Lupinus angustifolius L.)

A novel and highly effective source of anthracnose resistance in narrow-leafed lupin was identified. Resistance was shown to be governed by a single dominant locus. Molecular markers have been developed, which can be used for selecting resistant genotypes in lupin breeding. A screening for anthracnose resistance of a set of plant genetic resources of narrow-leafed lupin (Lupinus angustifolius L.) identified the breeding line Bo7212 as being highly resistant to anthracnose (Colletotrichum lupini). Segregation analysis indicated that the resistance of Bo7212 is inherited by a single dominant locus. The corresponding resistance gene was given the designation LanrBo. Previously published molecular anchor markers allowed us to locate LanrBo on linkage group NLL-11 of narrow-leafed lupin. Using information from RNAseq data obtained with inoculated resistant vs. susceptible lupin entries as well as EST-sequence information from the model genome Lotus japonicus, additional SNP and EST markers linked to LanrBo were derived. A bracket of two LanrBo-flanking markers allows for precise marker-assisted selection of the novel resistance gene in narrow-leafed lupin breeding programs.

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.