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

Variation in Ribosomal DNA in the Genus Trifolium (Fabaceae)

The genus Trifolium L. is characterized by basic chromosome numbers 8, 7, 6, and 5. We conducted a genus-wide study of ribosomal DNA (rDNA) structure variability in diploids and polyploids to gain insight into evolutionary history. We used fluorescent in situ hybridization to newly investigate rDNA variation by number and position in 30 Trifolium species. Evolutionary history among species was examined using 85 available sequences of internal transcribed spacer 1 (ITS1) of 35S rDNA. In diploid species with ancestral basic chromosome number (x = 8), one pair of 5S and 26S rDNA in separate or adjacent positions on a pair of chromosomes was prevalent. Genomes of species with reduced basic chromosome numbers were characterized by increased number of signals determined on one pair of chromosomes or all chromosomes. Increased number of signals was observed also in diploids Trifolium alpestre and Trifolium microcephalum and in polyploids. Sequence alignment revealed ITS1 sequences with mostly single nucleotide polymorphisms, and ITS1 diversity was greater in diploids with reduced basic chromosome numbers compared to diploids with ancestral basic chromosome number (x = 8) and polyploids. Our results suggest the presence of one 5S rDNA site and one 26S rDNA site as an ancestral state.

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.

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.

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

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.

Impact of Chromosomal Rearrangements on the Interpretation of Lupin Karyotype Evolution

Plant genome evolution can be very complex and challenging to describe, even within a genus. Mechanisms that underlie genome variation are complex and can include whole-genome duplications, gene duplication and/or loss, and, importantly, multiple chromosomal rearrangements. Lupins (Lupinus) diverged from other legumes approximately 60 mya. In contrast to New World lupins, Old World lupins show high variability not only for chromosome numbers (2n = 32–52), but also for the basic chromosome number (x = 5–9, 13) and genome size. The evolutionary basis that underlies the karyotype evolution in lupins remains unknown, as it has so far been impossible to identify individual chromosomes. To shed light on chromosome changes and evolution, we used comparative chromosome mapping among 11 Old World lupins, with Lupinus angustifolius as the reference species. We applied set of L. angustifolius-derived bacterial artificial chromosome clones for fluorescence in situ hybridization. We demonstrate that chromosome variations in the species analyzed might have arisen from multiple changes in chromosome structure and number. We hypothesize about lupin karyotype evolution through polyploidy and subsequent aneuploidy. Additionally, we have established a cytogenomic map of L. angustifolius along with chromosome markers that can be used for related species to further improve comparative studies of crops and wild lupins.

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.

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.

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.

Epigenomic diversification within the genus Lupinus

Deciphering the various chemical modifications of both DNA and the histone compound of chromatin not only leads to a better understanding of the genome-wide organisation of epigenetic landmarks and their impact on gene expression but may also provide some insights into the evolutionary processes. Although both histone modifications and DNA methylation have been widely investigated in various plant genomes, here we present the first study for the genus Lupinus. Lupins, which are members of grain legumes (pulses), are beneficial for food security, nutrition, health and the environment. In order to gain a better understanding of the epigenetic organisation of genomes in lupins we applied the immunostaining of methylated histone H3 and DNA methylation as well as whole-genome bisulfite sequencing. We revealed variations in the patterns of chromatin modifications at the chromosomal level among three crop lupins, i.e. L. angustifolius (2n = 40), L. albus (2n = 50) and L. luteus (2n = 52), and the legume model plant Medicago truncatula (2n = 16). Different chromosomal patterns were found depending on the specific modification, e.g. H3K4me2 was localised in the terminal parts of L. angustifolius and M. truncatula chromosomes, which is in agreement with the results that have been obtained for other species. Interestingly, in L. albus and L. luteus this modification was limited to one arm in the case of all of the chromosomes in the complement. Additionally, H3K9me2 was detected in all of the analysed species except L. luteus. DNA methylation sequencing (CG, CHG and CHH contexts) of aforementioned crop but also wild lupins such as L. cosentinii (2n = 32), L. digitatus (2n = 36), L. micranthus (2n = 52) and L. pilosus (2n = 42) supported the range of interspecific diversity. The examples of epigenetic modifications illustrate the diversity of lupin genomes and could be helpful for elucidating further epigenetic changes in the evolution of the lupin genome.

Potential Uses of Wild Germplasms of Grain Legumes for Crop Improvement

Challenged by population increase, climatic change, and soil deterioration, crop improvement is always a priority in securing food supplies. Although the production of grain legumes is in general lower than that of cereals, the nutritional value of grain legumes make them important components of food security. Nevertheless, limited by severe genetic bottlenecks during domestication and human selection, grain legumes, like other crops, have suffered from a loss of genetic diversity which is essential for providing genetic materials for crop improvement programs. Illustrated by whole-genome-sequencing, wild relatives of crops adapted to various environments were shown to maintain high genetic diversity. In this review, we focused on nine important grain legumes (soybean, peanut, pea, chickpea, common bean, lentil, cowpea, lupin, and pigeonpea) to discuss the potential uses of their wild relatives as genetic resources for crop breeding and improvement, and summarized the various genetic/genomic approaches adopted for these purposes.

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.