Next Article in Journal
Meeting the Challenge of Eliminating Chronic Hepatitis B Infection
Next Article in Special Issue
Development and Validation of a Gene-Targeted dCAPS Marker for Marker-Assisted Selection of Low-Alkaloid Content in Seeds of Narrow-Leafed Lupin (Lupinus angustifolius L.)
Previous Article in Journal
On the Importance to Acknowledge Transposable Elements in Epigenomic Analyses
Previous Article in Special Issue
Convergent Evolution of the Seed Shattering Trait
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Chromosomal Rearrangements on the Interpretation of Lupin Karyotype Evolution

1
Department of Genomics, Institute of Plant Genetics, Polish Academy of Sciences, 60-479 Poznan, Poland
2
Department of Biometry and Bioinformatics, Institute of Plant Genetics, Polish Academy of Sciences, 60-479 Poznan, Poland
3
Department of Plant Anatomy and Cytology, University of Silesia in Katowice, 40-032 Katowice, Poland
4
Center for Applied Genetic Technologies, University of Georgia, Athens, GA 30602, USA
*
Author to whom correspondence should be addressed.
Genes 2019, 10(4), 259; https://doi.org/10.3390/genes10040259
Submission received: 21 March 2019 / Revised: 27 March 2019 / Accepted: 27 March 2019 / Published: 1 April 2019
(This article belongs to the Special Issue Genomics of Plant Domestication and Crop Evolution)

Abstract

:
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.

1. Introduction

Knowing how genomes are organized and distinct at the chromosomal level is fundamental for an understanding of the dynamics of chromosomal structure and karyotypes [1]. Chromosome number variations have been studied as an intermediary of karyotype changes. Traditionally, chromosome numbers, genome size, and the number of nucleolar organizers have been used to establish the basis for further studies of chromosome organization, especially when data on whole genome sequences or differential chromosome markers are not available. Progress in molecular cytogenetics and the availability of bacterial artificial chromosome (BAC) libraries have advanced plant genome studies [2], to provide more precise identification of karyotypes and to subsequently integrate these with linkage groups [3,4]. Whole-genome sequencing was a breakthrough that has improved the quality and accuracy of genome maps, and has made it possible to assign chromosomes to pseudomolecules in various plant species [5,6,7]. Importantly, reference karyotypes have become a key tool for comparative chromosome mapping, e.g., among the legume genera Arachis [8], Glycine [9], and Phaseolus [10,11], as well as the genera of other families, including Solanum [12] and Daucus [13].
Plants of the Lupinus genus (lupins; Genisteae tribe) belong to the important grain legumes, along with those of the genera Glycine, Phaseolus, and Arachis. Lupins are a source of valuable protein, oils, and bioactive compounds, such as alkaloids, isoflavones, and oligosaccharides. Additionally, lupins contribute to the productivity and quality of agricultural soils, and thus they have substantial roles in crop rotation. The Lupinus genus consists of approximately 270 species, which can be divided into two groups, known as Old World lupins (OWLs) and New World lupins (NWLs). OWLs comprise 12 autogamous species, which include three crops: Lupinus angustifolius (narrow-leafed lupin), Lupinus albus (white lupin), and Lupinus luteus (yellow lupin) [14]. The OWLs are geographically and phenotypically diverse, and they are mainly represented in Mediterranean regions and North Africa, which highlights their adaptation to different environments. Based on their seed coat structure, the OWLs are traditionally divided into two groups. The first group includes the smooth-seeded lupins, which is divided into four sections: Angustifolius (2n = 40), Luteus (2n = 52), Albus (2n = 50), and Micranthus (2n = 52). These are distributed around the Mediterranean basin, and there are major genetic barriers between most species. Phylogenetic systematics using internal transcribed spacer data have shown that smooth-seeded species include at least two lineages, i.e., Angustifolius–Luteus and Micranthus–Albus [15]. The second group of the OWLs comprises the rough-seeded lupins, which have been divided into three sections: Pilosus (2n = 42), Atlanticus (2n = 32, 36, 38) and Princei (2n = 38). These are eco-geographically isolated and are mainly in North Africa.
At the chromosome level, OWLs are a complex group within legumes. They have chromosome numbers that range from 2n = 32 to 2n = 52, and a varied basic chromosome number x = 5 − 9, 13 [16], which indicates an obscure evolutionary pattern, as has been exemplified across L. angustifolius and Lupinus cryptanthus, Lupinus micranthus, Lupinus pilosus, and Lupinus cosentinii [17]. Recently, differences in chromatin modification and DNA methylation have also highlighted variations in the organization of the lupin genomes [18]. Progress in genetic and molecular studies in lupins has mostly arisen from crops, and has been linked, for example, with the development of genetic maps [19], whole-genome sequencing [20,21], and transcriptomic analyses [22]. Of note, L. angustifolius and its genomic resources is considered as the reference species in comparative studies. The assignment of chromosomes to linkage groups of L. angustifolius [23,24] established the starting point for the tracking of cytogenetic differences among crops and wild lupins.
Here, we present the first complete genome-based cytogenetic map of L. angustifolius. For this purpose, we assigned the reference karyotype [24] to its whole genome assembly [20]. We explored the possibilities of its application in lupin genome studies by providing the first chromosome markers for related species. Finally, we report on a first global view of the chromosome structure in lupin karyotypes and assume their complex patterns of evolutionary changes. A considerable part of this study was addressed to the fundamental question of whether a common model of lupin karyotype evolution exists, according to the most likely explanations of the observations. Our data provide a foundation for lupin chromosome and karyotype evolution studies.

2. Materials and Methods

2.1. Plant Material

Twelve species of both crop and wild lupins were used in this study: The crop L. angustifolius and its wild botanical form L. cryptanthus; the crop L. luteus and Lupinus hispanicus; the crop L. albus and its wild form Lupinus graecus; and six other wild lupins. Their names, basic taxonomy, and cytogenetic characteristics are provided in Table 1. Their origins and geographic distributions are shown in Supplementary Figure S1.
The seeds were germinated in Petri dishes at 25 °C, to obtain root tips that were suitable for mitotic chromosome preparations.

2.2. BAC DNA Isolation, Sequencing, and Labeling

BAC clones were originated from the genomic BAC library constructed by Kasprzak et al. [25]. The BACs used for assignment of linkage groups to chromosome of L. angustifolius [23,24] and other available in NCBI database are listed in Table 2 (Results section). BAC DNA was isolated using standard miniprep kits (QIAprep Spin; Qiagen, Hilden, Germany). BAC DNA sequencing was carried out using the long-read sequencing (PacBio), as described by Susek et al. [17], or by the next-generation sequencing (NGS) approach (Illumina) [24]. Whole BAC sequences (WBS) were imported into Geneious 8.1.6. (Biomatters, Ltd., Auckland, New Zealand) and aligned to the L. angustifolius pseudomolecules and/or scaffolds of Hane et al. [20]. The BAC-end sequences (BES) were also used when the WBS were not available. For generation of the fluorescence in situ hybridization (FISH) probe, BAC DNA was labeled by nick translation (Sigma-Aldrich, St. Louis, MI, USA), either with digoxygenin-11-dUTP (Sigma-Aldrich), or with tetramethylrhodamine-5-dUTP (Sigma-Aldrich).

2.3. Fluorescence In Situ Hybridization

Mitotic chromosome preparations and FISH using the BAC-based probes (BAC-FISH) were carried out according to Lesniewska et al. [23] for L. angustifolius, and as described by Susek et al. [17] for other lupins, with minor modifications. Root tips were treated in an enzyme solution of 40% (v/v) pectinase (Sigma-Aldrich), 3% (w/v) cellulase (Sigma-Aldrich), and 1.5% (w/v) cellulase ‘Onozuka R-10’ (Serva, Heidelberg, Germany), at 37 °C for 150 min for the primary roots of the crops, and for 60 min to 90 min for both the primary and the lateral roots of the wild lupins. The same set of BAC markers was used for all of the species in a given experiment, to maintain the same conditions of the BAC-FISH reaction. To localize the BAC-FISH signals, the chromosomes were counterstained with 4′,6-diamidino-2-phenylindole (DAPI) in Vectashield (Vector Laboratories, Burlingame, CA, USA).

2.4. Microscopy, Image Acquisition and Processing, Data Presentation

The FISH data were analyzed using an epifluorescence microscope (BX-60; Olympus, Tokyo, Japan). Images were acquired with a high-sensitivity monochromatic camera (Ikegami 47E), separately for each of three fluorochromes, using the appropriate excitation and emission filters. Grey images were pseudocolored (Wasabi; Hamamatsu Photonics v. 1.1.) and superimposed (Micrografx; Corel Picture Publisher 10.1 software). At least five images were analyzed to confirm the intensities and chromosomal distributions of the individual BAC-FISH signals. All of the lupin karyotypes were visualised (ChromDraw package; included in Bioconductor R packages; www.bioconductor.org/packages), as described by Janečka and Lysak [26].

3. Results

3.1. Assignment of Chromosomes to Pseudomolecules of the L. angustifolius Genome

Physical mapping of chromosomes is a robust and efficient way for comparative studies of karyotype structure and evolution based on the synteny and collinearity of related genomes. As a framework, we used the cytogenetic map of L. angustifolius that was developed by Wyrwa et al. [24] to assign the L. angustifolius chromosomes to pseudomolecules [20]. Different abbreviations are used for chromosome names in the literature, such as LANG [23] and Lang [24], while both linkage groups and pseudomolecules are named as NLL. Here, for clarity, we use the following abbreviations: Lang for L. angustifolius chromosomes, and NLL for pseudomolecules. As a consequence of this assignment, we propose to name the chromosomes and pseudomolecules as Lang01 to Lang20.
We selected 53 BAC clone-based cytogenetic markers of L. angustifolius that were previously generated by Lesniewska et al. [23] and Wyrwa et al. [24] to assign the chromosomes to linkage groups of L. angustifolius (see detailed data in Supplementary Table S1). Among these, there were some BAC clones (i.e., 1M23, 2B03, 8C03) that remained unassigned. We sequenced 22 BACs, as 18 by long-read sequencing (PacBio) and four with an NGS platform (Illumina). These resources were used to accurately integrate the chromosomes and pseudomolecules in the Lang vs. NLL map of L. angustifolius. This set of new WBS, along with other lupin BAC sequences available at the NCBI database, established efficient tools for reciprocal integration of chromosomes and pseudomolecules (Table 2). We assigned Lang vs. NLL using varied numbers of WBS for particular chromosomes: Five WBS for Lang06; four for Lang 08 and 13; three for Lang 03, 11, 17, and 20; and two for Lang 05, 10, 14, and 19. For Lang 01, 02, 04, 07, 09, 12, 15, 16, and 18, we used one WBS. In the cases where WBS and/or BES provided misassembled information between chromosomes and pseudomolecules, we used BAC-FISH to integrate the L. angustifolius chromosomes and pseudomolecules. As a consequence, seven chromosomes were re-assigned (Table 2, Lang 01, 07, 08, 13, 14, 19, 20).
The chromosome Lang01 was verified by BAC 43C18, although WBS 43C18 was identified in both pseudomolecules NLL01 and NLL13. BAC-FISH mapping precluded assignment to NLL13 (Figure 1A, Table 2). Interestingly, Lang07 was previously marked by the clone 74I10, and has now been anchored by BAC 2B03 on a different chromosome to BAC 74I10 (Figure 1B). Additionally, the alignment of 74I10 WBS to the whole genome sequence, along with BAC-FISH mapping, defined this clone as a marker for Lang19 (Figure 1C). Lang19 was described by two other BACs, 67C07 and 83F23, with the WBS of 67C07 positioned in NLL19 (Table 2). Unexpectedly, four BACs (84D22, 111B08, 142C04, 142D13) were assigned to Lang08, even though they all corresponded to NLL17 (Table 2). Furthermore, we verified that BACs 11G20, 8C03, and 51F15 belong to the same chromosome Lang13, although 11G20 WBS was aligned to many NLLs (Table 2). BAC-FISH using both clones 11G20 and 51F15 (Figure 1D) showed that these BACs belong to the same chromosome Lang13. Moreover, Lang13 was enriched with the marker 8C03, assigned to chromosome Lang13 (Figure 1E). Similar to Lang13, Lang14 with clones 138N02, 115C21, and 9K06 illustrated some miss-mapping, such as the alignment of 5′BES and 3′BES of 115C21 to two different pseudomolecules, Lang14 and Lang12, respectively (Table 2). Nevertheless, BAC-FISH with 9K06 and 115C21 confirmed that Lang14 corresponded to NLL14 (not shown). Within Lang16, five BACs were identified (115G22, 112E01, 72O21, 8A03, 87N22) in various NLLs based on BESs (Table 2). However, the whole sequence of BAC 87N22 was in NLL16. Considering that all of these BACs were mapped to the same chromosome by Wyrwa et al. [24], all five BACs were linked to the same chromosome, Lang16. Our integrated studies also led to the assignment of clone 1M23 from cluster-2 to pseudomolecule NLL20 (Table 2), and to its colocalization with BAC 17B07, which is a marker of chromosome Lang20 (Figure 1F). Finally, 13 chromosomes (i.e., Lang02–06, Lang09–12, Lang15, Lang16–18) were assigned to pseudomolecules of L. angustifolius, in agreement with former chromosome assignments to linkage groups [23,24]. Finally, combined analyses of multiple BAC-FISH reactions, BAC sequencing, and whole genome assignment led us to construct the first comprehensive ideogram of the L. angustifolius chromosomes (Figure 2).

3.2. Chromosome Variation in Lupins

The ideogram obtained for the L. angustifolius chromosomes (Figure 2) provides the foundation for comparative cytogenomic mapping among lupins. Such chromosome analyses have been conducted using heterologous BAC-FISH mapping among 11 species of the genus: L. cryptanthus (2n = 40, Lcry), L. luteus (2n = 52, Llut), L. hispanicus (2n = 52, Lhis), L. albus (2n = 50, Lalb), L. graecus (2n = 50, Lgra), L. micranthus (2n = 52, Lmic), L. atlanticus (2n = 38, Latl), Lupinus digitatus (2n = 36, Ldig), L. cosentinii (2n = 32, Lcos), L. pilosus (2n = 42, Lpil), and Lupinus palaestinus (2n = 42, Lpal). We comparatively mapped 52 BAC-clone-based chromosome markers of L. angustifolius karyotype to provide cytogenetic landmarks for related lupins. These BACs were either mapped to one unique locus (so-called ‘single’ BAC; indicated as S here) or to many loci along the chromosomes of the analyzed lupins (referred to as ‘repetitive’ BAC; indicated as R here). Particular clones were mapped either as S or R markers; e.g., clone 84A06 was mapped to one locus in L. cryptanthus, L. albus, L. graecus, L. micranthus, L. cosentinii, L. pilosus, and L. palaestinus, but had many loci in the remaining lupins. Their FISH patterns are summarized in Table 3. Among the 52 clones, we have provided a set of 41 markers for Lcry, 25 for Llut, 19 for Lhis, 25 for Lalb, 24 for Lgra, 26 for Lmic, 25 for Latl, 19 for Ldig, 22 for Lcos, 27 for Lpil, and 24 for Lpal (Figure 3). Eight BACs (i.e., 80B11, 59J08, 115G22, 112E01, 72O21, 8A03, 136C16, 17B07) were shown to be unique markers for all lupins.
The set of markers for each species led us to identify chromosomes in the 11 OWLs analyzed (Supplementary Figures S2–S12). These had varied numbers of Lang-like chromosomes—19 in Lpil (Supplementary Figure S11); 18 in Lmic (Supplementary Figure S7); 17 in Lcry (Supplementary Figure S2); 16 in Llut (Supplementary Figure S3) and Lcos (Supplementary Figure S10); 14 in Lalb (Supplementary Figure S5), Lgra (Supplementary Figure S6), and Lpal (Supplementary Figure S12); 13 in Latl (Supplementary Figure S8); 12 in Ldig (Supplementary Figure S9); and 11 in Lhis (Supplementary Figure S4). Due to no reference information about other lupin genomes, the chromosome sizes given here were based on L. angustifolius. Also, the chromosomes of the lupins studied were numbered according to the chromosomes of the reference species; e.g., Lcry02 corresponds to Lang02, and so on.
The ‘repetitive’ BACs were excluded from the chromosomal marker analyses, although they broaden our view about chromosome variations in lupins. Combining the S and R clones, we illustrated Lang-like chromosomes in related species, to identify different levels of synteny between the lupins (Supplementary Figure S13). We identified 11 R clones for Lcry, 27 for Llut, 33 for Lhis, 26 for Lalb, 27 for Lgra, 25 for Lmic, 27 for Latl, 30 for Ldig, 29 for Lcos, 24 for Lpil, and 27 for Lpal (Figure 3). Three clones were not detected in Ldig, and another one was not identified in Lalb, Lgra, Lmic, Lcos, Lpil, and Lpal (Table 3). Additionally, we estimated that eight clones (43C18, 28O01, 15P08, 11G20, 115C21, 5L11, 87N22, 68H10) hybridized to many loci in the chromosomes of all of the lupins analyzed. On the other hand, all of the clones from Lang 05, 11, and 13 were ‘repetitive’ only for Lhis, Ldig, and Lpal, respectively.
Based on the FISH patterns for both ‘single’ and ‘repetitive’ BACs, we distinguished four types of Lang-like chromosomes in related species. The first reflects the most varied BAC-FISH pattern, as Lang 03, 04, 06, 08, 10, 11, 13, 17, 19, and 20 illustrated that BACs from particular chromosomes were localized differently in most species. The second type was recognized as the most conserved pattern, e.g., two clones of Lang05 were localized as S in the chromosomes across all species, with the exception of L. hispanicus (for an example of comparative BAC-FISH mapping, see Figure 4). Moreover, BACs 15P08 and 59J08 from Lang09 were mapped as R and S in all of the lupins. Also, clones of Lang16 were mapped with the same pattern in all of the lupins. The third type comprised five Lang-like chromosomes, e.g., Lang02, where BAC 120E23 was mapped as S in most lupin species, except for L. albus, L. graecus, and L. digitatus, and Lang07, where the clone 2B03 was mapped as R only in L. micranthus and L. digitatus. Additionally, Lang14 and Lang15 showed two kinds of BAC-FISH patterns. The Lang12-like chromosomes were also included with this type, as the clone 94P05 was mapped as R in seven species, with the exceptions of L. luteus, L. albus, L. graecus, and L. digitatus. In Lang14, we recognized three clones as R in five species. In the other five lupins, two clones (138N02, 115C21) were mapped as R, while BAC 9K06 was mapped as S. Exceptionally, for L. cryptanthus, 138N02 and 9K06 were mapped as S, but BAC 115C21 hybridized as R. The same pattern was also noted for Lang15, where both clones of this chromosome were R in six lupins, but clone 5L11 and BAC 134F01 were seen as R and S, respectively, in five species. The last type of chromosomes contains the clones represented as ‘repetitive’ in Lang01-like and Lang18-like chromosomes in all species.
Surprisingly, Lang17 was covered by BAC markers that all were ‘single’ in the lupins analyzed, with two exceptions as L. cryptanthus and L. palaestinus, where BACs 111G03 and 3B118, respectively, were mapped as R. We also found that ‘repetitive’ clones enabled us to distinguish Lang10-like chromosome between crop lupin and its wild form, e.g., marker 77C13 mapped as S in Lalb and R in Lgra. This was similar to the Lang19-like chromosomes of Llut and Lhis, where clone 83F23 was mapped as R and S, respectively. Interestingly, two other clones from Lang19-like chromosomes were mapped with the same pattern (Table 3).

3.3. Chromosome Rearrangements in Lupins

Based on comparative analysis of chromosome variations, we distinguished two main types of chromosome rearrangements. The first type consists of structural changes that lead to either increases or decreases in chromosome numbers and involve chromosomes Lang06-like and Lang17-like. We tracked Lang06 based on five clones (44J16, 76K16, 80B11, 51D03, 127N17) that marked the two chromosome arms (here named as arm A, carrying clones 44J16, 76K16, and 80B11, and arm B, with BACs 51D03 and 127N17). We observed different mapping patterns (Figure 5), e.g., clones 44J16 and 80B11 hybridized to one chromosome of L. luteus (Llut06), and BACs 76K16 and 80B11 to one chromosome of L. albus (Lalb06) and L. graecus (Lgra06). However, clone 51D03 was identified on another chromosome (additional chromosomes were marked by ′ or ′′) of these three species, referred to as Llut06′ (e.g., Figure 5A, Supplementary Figure S3, karyotype of Llut), Lalb06′ (Supplementary Figure S5, karyotype of Llalb), and Lgra06′ (Supplementary Figure S6, karyotype of Lgra). In L. pilosus, clones 76K16 and 80B11 were mapped in two different chromosomes (Figure 5B), as Lpil06 and Lpil06′, respectively. In addition, clone 51D03 together with BAC 127N17 were mapped in different chromosomes to Lpil06 (for an example of multiple BAC-FISH reactions, see Figure 5C) and Lpil06’, named as Lpil06′′ (Supplementary Figure S11, karyotype of Lpil). In L. palaestinus, we observed two chromosomes that carried clones only from one arm of Lang06, one called Lpal06 that hybridized with both BACs 76K16 and 80B11, and the other called Lpal06′ that hybridized with the clone 76K16 (Figure 5D, Supplementary Figure S12, karyotype of Lpal). Of note, one chromosome arm of Lang06 was also detected in Lhis (Supplementary Figure S4, karyotype of Lhis), Ldig (Supplementary Figure S9, karyotype of Ldig), and Lcos (Supplementary Figure S10, karyotype of Lcos), with the difference that Ldig06 carried both clones 76K16 and 80B11 (Figure 5E), while Lhis06 (e.g., Figure 5F) and Lcos06 carried only BAC 80B11. On the other hand, in L. micranthus, we visualized two chromosomes (Figure 6G,H) that corresponded to one chromosome Lang06. In this species, clones 76K16 and 80B11 and clones 51D03 and 127N17 were mapped in two chromosomes, Lmic06 and Lmic06′, respectively (Supplementary Figure S7, karyotype of Lmic). The entire set of markers of Lang06 was identified in the Lcry06 chromosome of L. cryptanthus and the Latl06 chromosome of L. atlanticus. However, in L. atlanticus, only BAC 127N17 was mapped, which is one of two clones that originated from the same arm of Lang06 (Supplementary Figure S2, karyotype of Lcry; Supplementary Figure S8, karyotype of Latl).
We also found that clones 111G03, 136C16, and 3B18 of Lang17 mapped differently in the lupins analysed. In L. digitatus, three chromosomes carried these BACs, of which chromosome Ldig17 carried only BAC 3B18. However, clone 136C16 was mapped in another two chromosomes, in Ldig17’ and Lang17”, together with clone 111G03. In L. pilosus, two chromosomes were identified by clones that originated from Lang17. Clone 3B18 was the only one that mapped in chromosome Lpil17, while both 111G03 and 136C16 hybridized to chromosome Lpil17’. In L. palaestinus, one chromosome Lpal17 was marked by clones 111G03 and 136C16. Finally, in another seven lupin species (L. luteus, L. hispanicus, L. albus, L. graecus, L. micranthus, L. atlanticus, L. cosentinii), the BAC-FISH pattern of Lang17-like chromosome revealed its similarity to Lang17, with the exception of L. cryptanthus, where Lcry17 carried only two clones (3B18, 136C16).
The second group of chromosomal rearrangements involved insertions and/or deletions, as reflected by the differences in the BAC presence/absence in a given chromosome (for example, see Supplementary Figure S14). We found that eight chromosomes (Lang03–04, Lang08, Lang10–11, Lang13, Lang19–20) showed the greatest changes in the BAC-FISH patterns. All of these three clones (84A06, 57J20, 137N08) from Lang03 were mapped in L. cryptanthus, L. albus, and L. palaestinus, in agreement with L. angustifolius (Supplementary Figure S13). However, clones 84A06 and 57J20 were localized in the same chromosome in L. graceus (Lgra03), with 84A06 and 137N08 in L. micranthus (Lmic03), and L. cosentinii (Lcos03) and L. pilosus (Lpil03) (Supplementary Figure S13). In the case of Lang08, we determined that clone 111B08 was mapped together with either BAC 142D13 in L. luteus (Llut08) or 84D22 in L. cosentinii (Lcos08) and L. pilosus (Lpil08). The clones 84D22 and 142C04 were colocalized in one chromosome in L. micranthus (Lmic08) and L. cosentinii (Lcos08) (Supplementary Figures S7, S10). Additionally, we found another chromosome (Lmic08’) that carried BAC 111B08 in L. micranthus (Supplementary Figure S7). Tracking the localization of clones from Lang11 in related OWL species, we noted that one (BAC 36L23) of four clones from chromosome Lang11 was detected in Lang11-like chromosomes of L. luteus, L. micranthus, L. cosentinii, and L. pilosus (Supplementary Figures S3, S7, S10, S11). The second, BAC 97D16, was mapped in L. atlanticus (Latl11) and L. palaestinus (Lpal11) (Supplementary Figures S3, S7, S10, and S11). Additionally, clone 97D16 was identified along with clone 36L23 in L. hispanicus and L. graecus, and with 20A06 in L. albus (Supplementary Figures S4–S6).
The Lang13 chromosome was marked by three clones (11G20, 8C03, 51F15) and was visualized in other lupins, either by the clone 8C03 in L. micranthus (Lmic13) or 51F15 in seven other species: L. hispanicus, L. albus, L. graecus, L. atlanticus, L. digitatus, L. cosentinii, and L. pilosus. However, in L. luteus, we identified two chromosomes (Llut13, Llut13′) that carried 51F15 (Supplementary Figure S3). Chromosome markers of Lang19 have been used to show that BACs 83F23 and 67C07 were mapped in one Lang19-like chromosome in seven species, as L. cryptanthus, L. hispanicus, L. micranthus, L. digitatus, L. cosentinii, L. pilosus, and L. palaestinus, and one BAC 83F23 was found in L. luteus. Surprisingly, all of the clones in L. atlanticus were mapped with patterns that were typical for L. angustifolius. Also, in the case of the Lang20-specific clones, we observed that both 17B07 and 83F23 were localized in one Lang20-like chromosome of species such as L. albus, L. micranthus, L. atlanticus, L. digitatus, L. cosentinii, and L. pilosus, but only BAC 17B07 hybridized with Lhis20 and Lgra20 in L. hispanicus and L. graecus, respectively. All of the three Lang20-specific clones (17B07, 83F23, 1M23) were represented in L. cryptanthus and L. luteus (Supplementary Figures S2 and S3). In the Lang04-like chromosome, four lupins had a chromosome that carried both clones, but the other four only had one of them, BAC 47P22. Furthermore, within Lang10-like chromosomes, clone 77C13 was mapped in one chromosome of six lupin species, while its colocalization with clone 57K22 was detected only in two species: L. cryptanthus and L. micranthus (Supplementary Figures S2 and S7).

4. Discussion

4.1. Lupin Chromosomes Reveal Their Puzzling Evolution

Our analysis shows that multiple chromosomal rearrangements were responsible for the chromosome variations in OWLs. We propose a model that involves complex changes that shaped the variations in chromosome number and/or structure. We assume that the mechanisms driving these differences involved DNA gain or loss through insertion/duplication or deletion. In addition, chromosome changes can also be derived from rearrangements within or between chromosomes, through insertion or translocation [27]. In particular, we hypothesized that at least four types (I–IV) of rearrangements have impacted on lupin chromosome changes. Possible changes of Lang06-like and Lang17-like chromosomes (Figure 6) might be explained by translocations of both arms of Lang chromosomes, to lead to two Lang06-like chromosomes in L. luteus, L. albus, L. graecus and L. micranthus, and two Lang17-like chromosomes in L. pilosus. However, these chromosomes might also have arisen from chromosome breakage, and have then undergone chromosome fusion-fission. We have called these possible rearrangements type I. Chromosomal changes within type II might encompass translocation of arm A and arm B to three chromosomes, which would result in two chromosomes carrying arm A, and one with arm B. However, translocation of these arms might have involved two chromosomes (one with arm A, the other with arm B), with the subsequent duplication of chromosome-carrying clones from arm A. Additionally, chromosome breakage might occur as well, together with the translocation of arm A or the duplication of the entire chromosome-carrying clones from arm A. This type of change was exclusively seen for the Lang06-like chromosome in L. pilosus and the Lang17-like chromosomes in L. digitatus, where three chromosomes (Lpil06, Lpil06′, Lpil06′′ Ldig17, Ldig17′, Ldig′′; respectively) were identified. Type III rearrangements were seen for the Lang06-like chromosome in L. hispanicus, L. digitatus and L. cosentinii, and for the Lang17-like chromosome in L. palaestinus. They show changes that underlie either translocation of one of the chromosome arms (arm A) or chromosome breakage, with the loss of one of the chromosome arm (arm B). The fourth type of rearrangement (type IV) was detected uniquely in the Lang06-like chromosomes of L. palaestinus. It might have evolved via multiple changes, such as translocation of arm A to two chromosomes, translocation of arm A to one chromosome and then duplication of the entire chromosome, or chromosome breakage and loss of one chromosome carrying clones from B and then duplication of the entire chromosome with translocation of arm A. In L. cryptanthus (2n = 40) and L. atlanticus (2n = 38), Lang06-like chromosomes were identified with patterns typical to L. angustifolius, beyond ‘repetitive’ clone 51D03 in L. atlanticus. Changes within the Lang17-like chromosome were noted only in three species: L. digitatus, L. pilosus and L. palaestinus. In other lupins, the BAC-FISH patterns were visualised as for L. angustifolius, with the exception of the ‘repetitive’ BAC 111G03 in L. cryptanthus.
Chromosome translocation changes between lupins can be explained in a similar way to those identified between Phaseolus vulgaris and Vigna unguiculata, where two BACs from two chromosomes of P. vulgaris were mapped to one chromosome of V. unguiculata [28]. However, the possibility of chromosome breakage and then chromosome fusion-fission in lupins, which could well occur in centromere regions and lead to chromosome reduction, needs to be further demonstrated using additional centromeric markers. Despite the tendency of chromosome reduction as a general evolutionary event, the mechanisms of these events are unclear in lupins [17]. Dysploidy via chromosome fusion-fission has been widely recognized in grasses, such as the Brachypodium genus [29,30], as well as in Brassicaceae [31], but in legumes it has only been seen in Phaseolus leptostachyus, with its chromosome number 2n = 20 [32]. On the other hand, ascending or descending aneuploidy via duplication or losses of individual chromosomes might also have a pivotal role in plants [1]. It has been shown that descending aneuploidy occurred mostly in early diverging papilionoid lineages (to which lupins belong), especially in species with low basic chromosome numbers [33]. Aneuploidy is exemplified in the Arachis genus (2n = 2x = 18) [34], Lathyrus [35], and Medicago sections [36]. Both dysploidy and aneuploidy in lupins cannot be ruled out, considering that polyploidy is associated with chromosome number reduction in core papilionoids [33,37]. However, aneuploidy complex series [38] and duplication and triplication [39] have been already suggested as crucial players in lupin genome variations. Importantly, Drummond [38] noted that an additional evolutionary mechanism like allopolyploidy might be unique to OWLs. They also assumed that the diversification within NWLs with stable chromosome numbers of 2n = 36 or 2n = 48 resulted from ecological speciation chromosome changes related to ploidy, with the exception of eastern North American lupins, with chromosome 2n = 52.
Thus, we address the question whether the course of lupin karyotypes followed from karyotype 2n = 52 to 2n = 36. In theory, species with 2n = 52 can evolve from ‘ancestral’ species with 2n = 54. There is no data about lupin with 2n = 54 chromosomes; however, they might be extinct. Furthermore, there have not been species with 2n = 44, 46 and 48 described within OWLs, while in NWLs, species with 2n = 48 are widely seen [40,41]. It can be suggested that species with 2n = 54 were the ‘ancestral’ ones, whereby the species with lower chromosome numbers might have undergone three rounds of whole genome duplication (WGD), to result in 2n = 48, 2n = 42 and 2n = 36, with the assistance of aneuploidy giving chromosome numbers 2n = 52, 50, 40, 38, 32. In this sense, lupin karyotypes might also have arisen due to aneuploidy, as their chromosome number is often different from the simple multiple of the basic chromosome number, like in L. cosentinii (2n = 32), L. atlanticus (2n = 38), L. angustifolius (2n = 40), L. albus (2n = 50), and L. luteus (2n = 52). We now present here our hypothetic model of lupin chromosome evolution, which is based on the polyploidy events supplemented by the likely aneuploidy events (Figure 7). We assume that the basic chromosome number in OWLs is x = 6. The idea of the basic chromosome number x = 6 is consistent the same base number considered for NWLs [41] and the ancestral rosid karyotype model [42], but contrary to comparative phylogenetic analyses, where the ancestral basic chromosome number for lupins as a member of genistoids was postulated to be x = 9 [43].
Of importance, multiple and various rearrangements in chromosome structures support our hypothesis regarding variations in the lupin chromosomes. Lupin chromosome variation can be explained by different genome organization and evolutionary relationships between the reference L. angustifolius and the related lupins. Of note, the identification of chromosome changes was not directly linked to the number of clones used for comparative mapping. Four clones were enough to detect structural changes for Lang03 and Lang11 in all of the lupins. However, the localization of four clones in the Lang08-like chromosomes illustrated structural changes in the species analyzed, and chromosome number rearrangements only in L. micranthus. Furthermore, it can also be considered that an individual chromosome might have undergone its specific and independent evolutionary changes.
We realize that a cytogenetically based perspective can underestimate the evolutionary events that have shaped plant genomes. Lupins belong to the genistoid clade, which is known as one of the most complex groups in the facet of polyploidy, and one of the weakest supported groups in core papilionoids [33]. The advent of lupin genome sequencing data and the availability of extensive genomic resources should be accomplished in the future. Additionally, the power of informative cytogenetic mapping to track genome duplication or extensive synteny among species can be increased while suitable model species are adopted, like model soybean and common bean [44], or Arabidopsis thaliana and other species within Brassicaceae [45,46].
Thus, for further comparative chromosome analyses, we designated L. digitatus, L. pilosus, and L. palaestinus as the reference species to track chromosome evolution within lupins. These are species showing chromosome changes of both Lang06- and Lang17-like chromosomes (Figure 6). Interestingly, these species have chromosome numbers that directly correspond to the basic chromosome number x = 6. Phylogeny analyses indicated that L. digitatus and L. palaestinus are the youngest ones, with diversification times estimated at about 0.5 mya. The diversification time of L. pilosus was estimated as about 3 mya [47]. Although L. pilosus and L. palaestinus share the same chromosome number, they differ in the divergence times and different types of chromosome changes among the Lang06-like and Lang17-like chromosomes. This sheds light on different evolutionary patterns of these two species with 2n = 42, and also L. digitatus with 2n = 36. It can be hypothesized that L. digitatus and L. palaestinus have undergone at least one common round of WGD, while L. pilosus might have been driven by different WGD.
On the other hand, there is a lupin group with chromosome number 2n = 52 or 2n = 50. In contrast to L. digitatus, L. micranthus (2n = 52), and L. albus (2n = 50) are two of the oldest species within the OWLs, with divergence about 7 mya and 7.5 mya, respectively [47]. It is intriguing the way species carrying chromosome 2n = 50-52 might have evolved. We wonder whether species with chromosome number 2n = 52, like L. micranthus, L. luteus, and L. hispanicus, might have evolved through different evolutionary events, as L. luteus and L. hispanicus diverged about 2.0 mya [47,48]. Consider also that L. micranthus and L. luteus have the smallest (0.98 pg/2C DNA) and largest (2.44 pg/2C DNA) genomes, respectively, in OWLs. On the other hand, L. hispanicus showed Lang06-like chromosome changes that were the same as L. digitatus. To highlight the complexity of this lupin genome variation, we note that L. digitatus and L. hispanicus share the highest number of clones that were mapped as ‘repetitive’. They have different chromosome numbers and diversification rates, and they are geographically separated. L. digitatus is distributed mainly in the northern part of Africa (Sahara and Egypt), and L. hispanicus is in north-western Spain, Portugal, Greece, and Turkey [14].

4.2. Cytogenetic Resources for Lupin Genome Analyses

Here we have demonstrated that the integration of whole genome sequencing and chromosome mapping through BAC-FISH can be successfully used to determine the genome assembly of the genome of L. angustifolius. This integrative approach has allowed us to construct a comprehensive chromosomal map that might be erroneously interpreted, if only the sequence assembly was to be applied, as it was applied to precisely assign the chromosomes and pseudomolecules of Brachypodium distachyon [49] and cotton [50]. As a consequence, we can finally integrate the unassigned linkage groups, as was shown for cluster-2, which has QTLs for flowering time [51] and flowering time gene homologue markers [39]. This ideogram developed for L. angustifolius defines the complete genomic map through the high number of chromosome markers that were assigned using whole BAC sequences.
For example, in our previous study [17], Lang08 carrying clones BACs 84D22, 111B88, 142C04, and 142D13 were mapped to scaffold AOCW01048841of the L. angutifolius genome published by Yang, et al. [52], while in the present study it was localized in the pseudomolecule NLL17 of the L. angutifolius genome [20]. Even so, implementation of FISH with these BACs as probes, excludes them as markers of Lang17 [23]. In this case, we established them as markers for Lang08 by the supportive ideogram provided by Wyrwa et al. [24]. Conversely, some BACs were localized in pseudomolecules, while they were not mapped in linkage groups. Similarly, we were able to saturate the pseudomolecules NLL13 by clones 8C03, which remained unmapped in the genetic map [23]. In this way, the chromosome order was also reshuffled, as was done for Lang19 and Lang07.
As the number of sequenced genomes is rising rapidly, there is the need to produce the appropriate tools to establish high quality reference genomes among plants studied that can be used for analyses of genomic variation. Furthermore, and especially in the case of some complex genomes, BAC sequencing can still contribute to improved genome sequencing carried out by NGS methods [53], and can remain as an additional ‘check point’, even if advanced techniques like BioNano optical mapping and Hi-C sequencing are implemented for the comparative assembly of genomes [54].
In Lupinus species, only the whole genome sequence of L. angustifolius is available so far. Thus, the sequence-comparison-based evolutionary research in this group of plants is impossible, to date. In addition, the other genetic resources, such as genetic maps and reliable chromosome-specific markers, are also rather limited. Hence, the unambiguous BAC markers developed can help guide comparative research in lupins and can also be useful to distinguish crop and wild genotypes (e.g., L. albus from L. graecus). In addition, we have provided markers for individual chromosomes of L. cosentinii (2n = 32). Further reciprocal BAC-FISH within this species should provide the completed ideogram of its karyotype. Additionally, in L. albus (2n = 50), these markers can be used as anchors to integrate the pseudomolecules and chromosomes of the already sequenced white lupin genome (Peret laboratory; personal communication).
A few BAC clones per chromosome, similar to the number we used for particular lupin chromosomes, were successful used to compare karyotypes of Phaseolus [10] and to discover dysploidy in this genus, as mentioned above [32]. Considering the high chromosome numbers in lupins, there is the suggestion to track particular chromosome of complement in more detailed pathways, either by providing higher numbers of BAC probes for FISH, as were generated for chromosome 6 in tomato [12] and chromosome 9 in maize [55]. The most advanced comparative studies require the use of high numbers of clones assigned to particular chromosomes, to track conserved chromosomal syntenies [56], reveal hypothesized ancestral chromosome numbers, and reconstruct the course of chromosome rearrangements during evolution [57,58]. An alternative approach is to implement the promising oligo-painting method [59] to investigate changes at different chromosome or chromosome-arm scales.

5. Conclusions

We have shown that some Lang-like chromosomes retain the same BAC-FISH patterns, while others expand with multiple rearrangements. We cannot determine the direction of chromosomal changes that shaped these lupin genomes, as we have assumed that any individual chromosome might undergone independent evolution. However, a tendency of the chromosome number to be reduced through translocation and/or breakage in the lupin genus was seen (mainly in the Lang06-like and Lang17-like chromosomes, and a few other chromosomes). We propose L. digitatus (2n = 36), L. pilosus (2n = 42) and L. palaestinus (2n = 42) as reference species for comparative evolutionary analyses here. Given possibly three independent WGD in lupins, together with the predominant chromosomal data of chromosome numbers and phylogenetic analyses, we hypothesize that the lupin basic chromosomes with x = 6 arose from a reduction from 2n = 54, with assistance of aneuploidy. Lupinus are considered as polyploids, like for many agricultural crops. However, the evolutionary relationships between lupins in terms of their complex polyploidy origin are neither evident nor simple.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/10/4/259/s1, Table S1: Characterization of the BAC clone position in the pseudomolecules of the L. angustifolius genome (Hane et al. 2017), Figure S1: Origin and distribution of Old World lupins. The geographical range is based on Gladstones (1998), with minor editorial modifications, Figure S2–S12: Ideograms of individual chromosomes in lupins. (S2) L. cryptanthus. (S3) L. luteus. (S4) L. hispanicus (S5) L. albus. (S6) L. graecus. (S7) L. micranthus. (S8) L. atlanticus. (S9) L. digitatus. (S10) L. cosentinii. (S11) L. pilosus. (S12) L. palaestinus. BAC markers are indicated by rectangles with unique colors, Figure S13: Schematic visualization of the results of comparative BAC-FISH among lupins in this study. BAC clones mapped either as ‘single’ or ‘repetitive’ were indicated by single rectangle or triple rectangle with its unique color, respectively. ‘Repetitive’ clones were described by the R at the end of their names, for example the clone named as 43C18_R is ‘repetitive’, while the clone 43C18 is a ‘single’ one. The colors of rectangles were used according to the RGB code. All chromosomes were drawn to scale to L. angustifolius chromosomes and BAC positions in L. angustifolius genome (Hane et al. 2017), Figure S14: Diagrams with examples of the changes in chromosome structure, deduced according to the presence/absence of ‘single’ BACs in particular Lang-like chromosomes.

Author Contributions

Conceptualization, K.S.; methodology, K.S.; investigation, K.S., W.B., and K.B.C.; data acquisition, K.S. and W.B.; analysis and interpretation of data, K.S.; data curation, K.S.; writing—original draft preparation, K.S.; writing/editing, K.S., R.H., and S.A.J.; visualization of data, K.S.; supervision of the research, K.S.; project administration, K.S.; funding acquisition, K.S.; final reading, K.S., W.B., K.B.C., R.H, S.A.J., B.W., and B.N.

Funding

This study was supported by the National Science Centre, Poland (grant no. OPUS 2 2011/03/B/NZ2/01420 and HARMONIA 7 2015/18/M/NZ2/00422 to KS).

Acknowledgments

We thank Dorota Narożna (Poznań University of Life Sciences, Poznań, Poland) for providing BAC 94P05 sequences and Jan Pomorski (Museum and Institute of Zoology, of the Polish Academy of Sciences, Warsaw, Poland) for BAC sequencing.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Schubert, I.; Lysak, M.A. Interpretation of karyotype evolution should consider chromosome structural constraints. Trends Genet. 2011, 27, 207–216. [Google Scholar] [CrossRef] [PubMed]
  2. Doležel, J.; Vrána, J.; Cápal, P.; Kubaláková, M.; Burešová, V.; Šimková, H. Advances in plant chromosome genomics. Biotechnol. Adv. 2014, 32, 122–136. [Google Scholar] [CrossRef] [PubMed]
  3. Iwata-Otsubo, A.; Lin, J.Y.; Gill, N.; Jackson, S.A. Highly distinct chromosomal structures in cowpea (Vigna unguiculata), as revealed by molecular cytogenetic analysis. Chromosome Res. 2016, 24, 197–216. [Google Scholar] [CrossRef] [PubMed]
  4. Kim, H.; Choi, S.R.; Bae, J.; Hong, C.P.; Lee, S.Y.; Hossain, M.J.; Van Nguyen, D.; Jin, M.; Park, B.-S.; Bang, J.-W.; et al. Sequenced BAC anchored reference genetic map that reconciles the ten individual chromosomes of Brassica rapa. BMC Genom. 2009, 10, 432. [Google Scholar] [CrossRef] [PubMed]
  5. Sun, J.; Zhang, Z.; Zong, X.; Huang, S.; Li, Z.; Han, Y. A high-resolution cucumber cytogenetic map integrated with the genome assembly. BMC Genom. 2013, 14, 461. [Google Scholar] [CrossRef]
  6. Chao, Y.-T.; Chen, W.-C.; Chen, C.-Y.; Ho, H.-Y.; Yeh, C.-H.; Kuo, Y.-T.; Su, C.-L.; Yen, S.-H.; Hsueh, H.-Y.; Yeh, J.-H.; et al. Chromosome-level assembly, genetic and physical mapping of Phalaenopsis aphrodite genome provides new insights into species adaptation and resources for orchid breeding. Plant Biotechnol. J. 2018, 16, 2027–2041. [Google Scholar] [CrossRef]
  7. Shearer, L.A.; Anderson, L.K.; de Jong, H.; Smit, S.; Goicoechea, J.L.; Roe, B.A.; Hua, A.; Giovannoni, J.J.; Stack, S.M. Fluorescence in situ hybridization and optical mapping to correct scaffold arrangement in the tomato genome. G3 (Bethesda) 2014, 4, 1395–1405. [Google Scholar] [CrossRef]
  8. Zhang, L.; Yang, X.; Tian, L.; Chen, L.; Yu, W. Identification of peanut (Arachis hypogaea) chromosomes using a fluorescence in situ hybridization system reveals multiple hybridization events during tetraploid peanut formation. New Phytol. 2016, 211, 1424–1439. [Google Scholar] [CrossRef]
  9. Findley, S.D.; Cannon, S.; Varala, K.; Du, J.; Ma, J.; Hudson, M.E.; Birchler, J.A.; Stacey, G. A fluorescence in situ hybridization system for karyotyping soybean. Genetics 2010, 185, 727–744. [Google Scholar] [CrossRef] [PubMed]
  10. Fonsêca, A.; Ferreira, J.; dos Santos, T.R.B.; Mosiolek, M.; Bellucci, E.; Kami, J.; Gepts, P.; Geffroy, V.; Schweizer, D.; dos Santos, K.G.B.; et al. Cytogenetic map of common bean (Phaseolus vulgaris L.). Chromosome Res. 2010, 18, 487–502. [Google Scholar] [CrossRef]
  11. Fonsêca, A.; Pedrosa-Harand, A. Karyotype stability in the genus Phaseolus evidenced by the comparative mapping of the wild species Phaseolus microcarpus. Genome 2013, 56, 335–343. [Google Scholar] [CrossRef]
  12. Szinay, D.; Wijnker, E.; van den Berg, R.; Visser, R.G.F.; Jong, H.D.; Bai, Y. Chromosome evolution in Solanum traced by cross-species BAC-FISH. New Phytol. 2012, 195, 688–698. [Google Scholar] [CrossRef]
  13. Iovene, M.; Cavagnaro, P.F.; Senalik, D.; Buell, C.R.; Jiang, J.; Simon, P.W. Comparative FISH mapping of Daucus species (Apiaceae family). Chromosome Res. 2011, 19, 493–506. [Google Scholar] [CrossRef] [PubMed]
  14. Gladstones, J.S. Distribution, origin, taxonomy, history and importance. In Lupins as Crop Plants: Biology, Production, and Utilization; Gladstones, J.S., Atkins, C.A., Hamblin, J., Eds.; CAB International: Wallingford, UK, 1998; pp. 1–36. [Google Scholar]
  15. Aïnouche, A.; Bayer, R.J.; Misset, M.T. Molecular phylogeny, diversification and character evolution in Lupinus (Fabaceae) with special attention to Mediterranean and African lupines. Plant Syst. Evol. 2004, 246, 211–222. [Google Scholar] [CrossRef]
  16. Naganowska, B.; Wolko, B.; Sliwinska, E.; Kaczmarek, Z. Nuclear DNA content variation and species relationships in the genus Lupinus (Fabaceae). Ann. Bot. 2003, 92, 349–355. [Google Scholar] [CrossRef]
  17. Susek, K.; Bielski, W.K.; Hasterok, R.; Naganowska, B.; Wolko, B. A first glimpse of wild lupin karyotype variation as revealed by comparative cytogenetic mapping. Front. Plant Sci. 2016, 7, 1152. [Google Scholar] [CrossRef] [PubMed]
  18. Susek, K.; Braszewska-Zalewska, A.; Bewick, A.J.; Hasterok, R.; Schmitz, R.J.; Naganowska, B. Epigenomic diversification within the genus Lupinus. PLoS ONE 2017, 12, e0179821. [Google Scholar] [CrossRef] [PubMed]
  19. Nelson, M.N.; Moolhuijzen, P.M.; Boersma, J.G.; Chudy, M.; Lesniewska, K.; Bellgard, M.; Oliver, R.P.; Swiecicki, W.; Wolko, B.; Cowling, W.A.; et al. Aligning a new reference genetic map of Lupinus angustifolius with the genome sequence of the model legume, Lotus japonicus. DNA Res. 2010, 17, 73–83. [Google Scholar] [CrossRef] [PubMed]
  20. Hane, J.K.; Ming, Y.; Kamphuis, L.G.; Nelson, M.N.; Garg, G.; Atkins, C.A.; Bayer, P.E.; Bravo, A.; Bringans, S.; Cannon, S.; et al. A comprehensive draft genome sequence for lupin (Lupinus angustifolius), an emerging health food: Insights into plant–microbe interactions and legume evolution. Plant Biotechnol. J. 2017, 15, 318–330. [Google Scholar] [CrossRef] [PubMed]
  21. Książkiewicz, M.; Nazzicari, N.; Yang, H.A.; Nelson, M.N.; Renshaw, D.; Rychel, S.; Ferrari, B.; Carelli, M.; Tomaszewska, M.; Stawiński, S.; et al. A high-density consensus linkage map of white lupin highlights synteny with narrow-leafed lupin and provides markers tagging key agronomic traits. Sci. Rep. 2017, 7, 15335. [Google Scholar] [CrossRef]
  22. Parra-González, L.B.; Aravena-Abarzúa, G.A.; Navarro-Navarro, C.S.; Udall, J.; Maughan, J.; Peterson, L.M.; Salvo-Garrido, H.E.; Maureira-Butler, I.J. Yellow lupin (Lupinus luteus L.) transcriptome sequencing: Molecular marker development and comparative studies. BMC Genom. 2012, 13, 1–15. [Google Scholar] [CrossRef]
  23. Lesniewska, K.; Ksiazkiewicz, M.; Nelson, M.N.; Mahe, F.; Ainouche, A.; Wolko, B.; Naganowska, B. Assignment of 3 genetic linkage groups to 3 chromosomes of narrow-leafed lupin. J. Hered. 2011, 102, 228–236. [Google Scholar] [CrossRef]
  24. Wyrwa, K.; Ksiazkiewicz, M.; Szczepaniak, A.; Susek, K.; Podkowinski, J.; Naganowska, B. Integration of Lupinus angustifolius L. (narrow-leafed lupin) genome maps and comparative mapping within legumes. Chromosome Res. 2016, 24, 355–378. [Google Scholar] [CrossRef]
  25. Kasprzak, A.; Safar, J.; Janda, J.; Dolezel, J.; Wolko, B.; Naganowska, B. The bacterial artificial chromosome (BAC) library of the narrow-leafed lupin (Lupinus angustifolius L.). Cell. Mol. Biol. Lett. 2006, 11, 396–407. [Google Scholar] [CrossRef]
  26. Janečka, J.; Lysak, M.A. chromDraw: An R package for visualization of linear and circular karyotypes. Chromosome Res. 2016, 24, 217–223. [Google Scholar] [CrossRef]
  27. Schubert, I. Chromosome evolution. Curr. Opin. Plant Biol. 2007, 10, 109–115. [Google Scholar] [CrossRef]
  28. Vasconcelos, E.V.; de Andrade Fonsêca, A.F.; Pedrosa-Harand, A.; de Andrade Bortoleti, K.C.; Benko-Iseppon, A.M.; da Costa, A.F.; Brasileiro-Vidal, A.C. Intra- and interchromosomal rearrangements between cowpea [Vigna unguiculata (L.) Walp.] and common bean (Phaseolus vulgaris L.) revealed by BAC-FISH. Chromosome Res. 2015, 23, 253–266. [Google Scholar] [CrossRef]
  29. Wolny, E.; Lesniewska, K.; Hasterok, R.; Langdon, T. Compact genomes and complex evolution in the genus Brachypodium. Chromosoma 2011, 120, 199–212. [Google Scholar] [CrossRef]
  30. Lusinska, J.; Majka, J.; Betekhtin, A.; Susek, K.; Wolny, E.; Hasterok, R. Chromosome identification and reconstruction of evolutionary rearrangements in Brachypodium distachyon, B. stacei and B. hybridum. Ann. Bot. 2018, 122, 445–459. [Google Scholar] [CrossRef]
  31. Lysak, M.A.; Mandáková, T.; Schranz, M.E. Comparative paleogenomics of crucifers: Ancestral genomic blocks revisited. Curr. Opin. Plant Biol. 2016, 30, 108–115. [Google Scholar] [CrossRef]
  32. Fonsêca, A.; Ferraz, M.E.; Pedrosa-Harand, A. Speeding up chromosome evolution in Phaseolus: Multiple rearrangements associated with a one-step descending dysploidy. Chromosoma 2016, 125, 413–421. [Google Scholar] [CrossRef]
  33. Doyle, J.J. Polyploidy in Legumes. In Polyploidy and Genome Evolution; Soltis, P.S., Soltis, D.E., Eds.; Springer: Berlin/Heidelberg, Germany, 2012; pp. 147–180. [Google Scholar] [CrossRef]
  34. Leal-Bertioli, S.C.M.; Santos, S.P.; Dantas, K.M.; Inglis, P.W.; Nielen, S.; Araujo, A.C.G.; Silva, J.P.; Cavalcante, U.; Guimarães, P.M.; Brasileiro, A.C.M.; et al. Arachis batizocoi: A study of its relationship to cultivated peanut (A. hypogaea) and its potential for introgression of wild genes into the peanut crop using induced allotetraploids. Ann. Bot. 2015, 115, 237–249. [Google Scholar] [CrossRef]
  35. Talukdar, D. Meiotic consequences of selfing in grass pea (Lathyrus sativus L.) autotetraploids in the advanced generations: Cytogenetics of chromosomal rearrangement and detection of aneuploids. Nucleus 2012, 55, 73–82. [Google Scholar] [CrossRef]
  36. Steele, K.P.; Ickert-Bond, S.M.; Zarre, S.; Wojciechowski, M.F. Phylogeny and character evolution in Medicago (Leguminosae): Evidence from analyses of plastid trnK/matK and nuclear GA3ox1 sequences. Am. J. Bot. 2010, 97, 1142–1155. [Google Scholar] [CrossRef]
  37. Murat, F.; Zhang, R.; Guizard, S.; Gavranovic, H.; Flores, R.; Steinbach, D.; Quesneville, H.; Tannier, E.; Salse, J. Karyotype and gene order evolution from reconstructed extinct ancestors highlight contrasts in genome plasticity of modern rosid crops. Genome Biol. Evol. 2015, 7, 735–749. [Google Scholar] [CrossRef]
  38. Drummond, C.S. Diversification of Lupinus (Leguminosae) in the western New World: Derived evolution of perennial life history and colonization of montane habitats. Mol. Phylogenet. Evol. 2008, 48, 408–421. [Google Scholar] [CrossRef]
  39. Kroc, M.; Koczyk, G.; Swiecicki, W.; Kilian, A.; Nelson, M.N. New evidence of ancestral polyploidy in the Genistoid legume Lupinus angustifolius L. (narrow-leafed lupin). Theor. Appl. Genet. 2014, 127, 1237–1249. [Google Scholar] [CrossRef]
  40. Conterato, I.F.; Schifino-Wittmann, M.T. New chromosome numbers, meiotic behaviour and pollen fertility in American taxa of Lupinus (Leguminosae): Contributions to taxonomic and evolutionary studies. Bot. J. Linn. Soc. 2006, 150, 229–240. [Google Scholar] [CrossRef]
  41. Naganowska, B.; Wolko, B.; Śliwińska, E.; Kaczmarek, Z.; Schifino-Wittmann, M. 2C DNA variation and relationships among New World species of the genus Lupinus (Fabaceae). Plant Syst. Evol. 2006, 256, 147–157. [Google Scholar] [CrossRef]
  42. Murat, F.; Zhang, R.; Guizard, S.; Flores, R.; Armero, A.; Pont, C.; Steinbach, D.; Quesneville, H.; Cooke, R.; Salse, J. Shared subgenome dominance following polyploidization explains grass genome evolutionary plasticity from a seven protochromosome ancestor with 16K protogenes. Genome Biol. Evol. 2014, 6, 12–33. [Google Scholar] [CrossRef]
  43. Cannon, S.B.; McKain, M.R.; Harkess, A.; Nelson, M.N.; Dash, S.; Deyholos, M.K.; Peng, Y.; Joyce, B.; Stewart, C.N., Jr.; Rolf, M.; et al. Multiple polyploidy events in the early radiation of nodulating and nonnodulating legumes. Mol. Biol. Evol. 2015, 32, 193–210. [Google Scholar] [CrossRef]
  44. Lin, J.Y.; Stupar, R.M.; Hans, C.; Hyten, D.L.; Jackson, S.A. Structural and functional divergence of a 1-Mb duplicated region in the soybean (Glycine max) genome and comparison to an orthologous region from Phaseolus vulgaris. Plant Cell 2010, 22, 2545–2561. [Google Scholar] [CrossRef] [PubMed]
  45. Mandáková, T.; Guo, X.; Özüdoğru, B.; Mummenhoff, K.; Lysak, M.A. Hybridization-facilitated genome merger and repeated chromosome fusion after 8 million years. Plant J. 2018, 96, 748–760. [Google Scholar] [CrossRef]
  46. Mandáková, T.; Kovařík, A.; Zozomová-Lihová, J.; Shimizu-Inatsugi, R.; Shimizu, K.K.; Mummenhoff, K.; Marhold, K.; Lysak, M.A. The more the merrier: Recent hybridization and polyploidy in Cardamine. Plant Cell 2013, 25, 3280–3295. [Google Scholar] [CrossRef] [PubMed]
  47. Drummond, C.S.; Eastwood, R.J.; Miotto, S.T.S.; Hughes, C.E. Multiple continental radiations and correlates of diversification in Lupinus (Leguminosae): Testing for key innovation with incomplete taxon sampling. Syst. Biol. 2012, 61, 443–460. [Google Scholar] [CrossRef] [PubMed]
  48. Mahe, F.; Markova, D.; Pasquet, R.; Misset, M.T.; Ainouche, A. Isolation, phylogeny and evolution of the SymRK gene in the legume genus Lupinus L. Mol. Phylogenet. Evol. 2011, 60, 49–61. [Google Scholar] [CrossRef]
  49. International Brachypodium Initiative. Genome sequencing and analysis of the model grass Brachypodium distachyon. Nature 2010, 463, 763–768. [Google Scholar] [CrossRef] [PubMed]
  50. Saski, C.A.; Scheffler, B.E.; Hulse-Kemp, A.M.; Liu, B.; Song, Q.; Ando, A.; Stelly, D.M.; Scheffler, J.A.; Grimwood, J.; Jones, D.C.; et al. Sub genome anchored physical frameworks of the allotetraploid Upland cotton (Gossypium hirsutum L.) genome, and an approach toward reference-grade assemblies of polyploids. Sci. Rep. 2017, 7, 15274. [Google Scholar] [CrossRef] [PubMed]
  51. Boersma, J.G.; Li, C.; Leśniewska, K.; Sivasithamparam, K.; Yang, H. Identification of quantitative trait loci (QTLs) influencing early vigour, height, flowering date, and seed size and their implications for breeding of narrow-leafed lupin (Lupinus angustifolius L.). Aust. J. Agr. Res. 2008, 59, 527–535. [Google Scholar] [CrossRef]
  52. Yang, H.; Tao, Y.; Zheng, Z.; Zhang, Q.; Zhou, G.; Sweetingham, M.W.; Howieson, J.G.; Li, C. Draft genome sequence, and a sequence-defined genetic linkage map of the legume crop species Lupinus angustifolius L. PLoS ONE 2013, 8, e64799. [Google Scholar] [CrossRef]
  53. Visendi, P.; Berkman, P.J.; Hayashi, S.; Golicz, A.A.; Bayer, P.E.; Ruperao, P.; Hurgobin, B.; Montenegro, J.; Chan, C.-K.K.; Staňková, H.; et al. An efficient approach to BAC based assembly of complex genomes. Plant Methods 2016, 12, 2. [Google Scholar] [CrossRef]
  54. Kolmogorov, M.; Armstrong, J.; Raney, B.J.; Streeter, I.; Dunn, M.; Yang, F.; Odom, D.; Flicek, P.; Keane, T.M.; Thybert, D.; et al. Chromosome assembly of large and complex genomes using multiple references. Genome Res. 2018, 28, 1720–1732. [Google Scholar] [CrossRef]
  55. Danilova, T.V.; Birchler, J.A. Integrated cytogenetic map of mitotic metaphase chromosome 9 of maize: Resolution, sensitivity, and banding paint development. Chromosoma 2008, 117, 345–356. [Google Scholar] [CrossRef]
  56. Mandáková, T.; Lysak, M.A. Healthy roots and leaves: Comparative genome structure of horseradish and watercress. Plant Physiol. 2019, 179, 66–73. [Google Scholar] [CrossRef] [PubMed]
  57. Mandakova, T.; Pouch, M.; Harmanova, K.; Zhan, S.H.; Mayrose, I.; Lysak, M.A. Multispeed genome diploidization and diversification after an ancient allopolyploidization. Mol. Ecol. 2017, 26, 6445–6462. [Google Scholar] [CrossRef] [PubMed]
  58. Lysak, M.A.; Berr, A.; Pecinka, A.; Schmidt, R.; McBreen, K.; Schubert, I. Mechanisms of chromosome number reduction in Arabidopsis thaliana and related Brassicaceae species. Proc. Natl. Acad. Sci. USA 2006, 103, 5224. [Google Scholar] [CrossRef] [PubMed]
  59. Han, Y.; Zhang, T.; Thammapichai, P.; Weng, Y.; Jiang, J. Chromosome-specific painting in Cucumis species using bulked oligonucleotides. Genetics 2015, 200, 771–779. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Fluorescence in situ hybridization verification of BAC localization in chromosomes of L. angustifolius. (A) BAC 43C18 (red signals)—Lang 01 and BAC 51F15 (green signals)—Lang 13. (B) BAC 2B03 (red)—Lang07 and 74I10 (green)—Lang19. (C) 67C07 (red)—Lang19 and 74I10 (green)—Lang19. (D) BAC 11G20 (red)—Lang13 and BAC 51F15 (green)—Lang13. (E) BAC 51F15 (red)—Lang13 and 8C03 (green)—Lang13. (F) 17B07 (red)—Lang20 and 1M23 (green)—Lang20. Chromosomes were stained with DAPI (blue). Scale bars = 5 µm.
Figure 1. Fluorescence in situ hybridization verification of BAC localization in chromosomes of L. angustifolius. (A) BAC 43C18 (red signals)—Lang 01 and BAC 51F15 (green signals)—Lang 13. (B) BAC 2B03 (red)—Lang07 and 74I10 (green)—Lang19. (C) 67C07 (red)—Lang19 and 74I10 (green)—Lang19. (D) BAC 11G20 (red)—Lang13 and BAC 51F15 (green)—Lang13. (E) BAC 51F15 (red)—Lang13 and 8C03 (green)—Lang13. (F) 17B07 (red)—Lang20 and 1M23 (green)—Lang20. Chromosomes were stained with DAPI (blue). Scale bars = 5 µm.
Genes 10 00259 g001
Figure 2. Ideogram of L. angustifolius chromosomes (Lang01 to Lang20). The karyotype was created based on reciprocal assignment to the whole genome of L. angustifolius. The BAC clones are indicated by rectangles with a unique color, according to their position in the genome. The overlapping rectangles are shown overlapping BAC positions in L. angustifolius. All of the chromosomes are drawn to scale, whereby the Mb units refer to the L. angustifolius genome sequence.
Figure 2. Ideogram of L. angustifolius chromosomes (Lang01 to Lang20). The karyotype was created based on reciprocal assignment to the whole genome of L. angustifolius. The BAC clones are indicated by rectangles with a unique color, according to their position in the genome. The overlapping rectangles are shown overlapping BAC positions in L. angustifolius. All of the chromosomes are drawn to scale, whereby the Mb units refer to the L. angustifolius genome sequence.
Genes 10 00259 g002
Figure 3. Representation of ‘single’ (S) and ‘repetitive’ (R) BAC clones in 11 lupin genomes. Lang, L. angustifolius (reference species). Lcry, L. cryptanthus. Llut, L. luteus. Lhis, L. hispanicus. Lalb, L. albus. Lgra, L. graecus. Lmic, L. micranthus. Latl, L. atlanticus. Ldig, L. digitatus. Lcos, L. cosentini. Lpil, L. pilosus. Lpal, L. palaestinus. N/D, not detected.
Figure 3. Representation of ‘single’ (S) and ‘repetitive’ (R) BAC clones in 11 lupin genomes. Lang, L. angustifolius (reference species). Lcry, L. cryptanthus. Llut, L. luteus. Lhis, L. hispanicus. Lalb, L. albus. Lgra, L. graecus. Lmic, L. micranthus. Latl, L. atlanticus. Ldig, L. digitatus. Lcos, L. cosentini. Lpil, L. pilosus. Lpal, L. palaestinus. N/D, not detected.
Genes 10 00259 g003
Figure 4. Identification of the Lang05-like chromosome in related lupins using fluorescence in situ hybridization (FISH) with BAC 131K22 (red signals) and BAC 4G15 (green signals). (A) L. albus. (B) L. pilosus. (C) L. palaestinus. (D) L. atlanticus. (E) L. digitatus. (F) L. cosentinii. Chromosomes were stained with DAPI (blue). Scale bars = 5 µm.
Figure 4. Identification of the Lang05-like chromosome in related lupins using fluorescence in situ hybridization (FISH) with BAC 131K22 (red signals) and BAC 4G15 (green signals). (A) L. albus. (B) L. pilosus. (C) L. palaestinus. (D) L. atlanticus. (E) L. digitatus. (F) L. cosentinii. Chromosomes were stained with DAPI (blue). Scale bars = 5 µm.
Genes 10 00259 g004
Figure 5. Localization of BAC clones from chromosome Lang06 in chromosomes of related lupins. (A) L. luteus—BAC 44J16 (red signals) and BAC 51D03 (green signals). (B) L. pilosus—BAC 80B11 (red) and BAC 76K16 (green). (C) L. pilosus—BAC 51D03 (red) and BAC 76K16 (green). (D) L. palaestinus—BAC 80B11 (red) and BAC 76K16 (green). (E) L. digitatus – BAC 80B11 (red) and BAC 76K16 (green). (F) L. hispanicus—BAC 80B11 (red). (G) L. micranthus – BAC 80B11 (red) and BAC 127N17 (green). (H) L. micranthus—BAC 51D03 (red) and BAC 127N17 (green). Chromosomes were stained with DAPI (blue). Arrows in (H) show the localization of BAC 80B11 (red). Scale bars, 5 µm.
Figure 5. Localization of BAC clones from chromosome Lang06 in chromosomes of related lupins. (A) L. luteus—BAC 44J16 (red signals) and BAC 51D03 (green signals). (B) L. pilosus—BAC 80B11 (red) and BAC 76K16 (green). (C) L. pilosus—BAC 51D03 (red) and BAC 76K16 (green). (D) L. palaestinus—BAC 80B11 (red) and BAC 76K16 (green). (E) L. digitatus – BAC 80B11 (red) and BAC 76K16 (green). (F) L. hispanicus—BAC 80B11 (red). (G) L. micranthus – BAC 80B11 (red) and BAC 127N17 (green). (H) L. micranthus—BAC 51D03 (red) and BAC 127N17 (green). Chromosomes were stained with DAPI (blue). Arrows in (H) show the localization of BAC 80B11 (red). Scale bars, 5 µm.
Genes 10 00259 g005
Figure 6. Scheme of the structural rearrangements within Lang06 and Lang17 of L. angustifolius and the related lupins. The pink and green rectangles correspond to arm A and arm B of both of the Lang06 and Lang17 chromosomes in related lupins.
Figure 6. Scheme of the structural rearrangements within Lang06 and Lang17 of L. angustifolius and the related lupins. The pink and green rectangles correspond to arm A and arm B of both of the Lang06 and Lang17 chromosomes in related lupins.
Genes 10 00259 g006
Figure 7. Hypothetical pattern of direction of lupin chromosome evolution through polyploidy and aneuploidy events. † describes extinct species; considering x = 6 the 1 × 6 corresponds to one round of whole genome duplication (WGD), while 3 × 6 three rounds of WGD; number 2 (4) indicates the aneuploidy events.
Figure 7. Hypothetical pattern of direction of lupin chromosome evolution through polyploidy and aneuploidy events. † describes extinct species; considering x = 6 the 1 × 6 corresponds to one round of whole genome duplication (WGD), while 3 × 6 three rounds of WGD; number 2 (4) indicates the aneuploidy events.
Genes 10 00259 g007
Table 1. General characteristics of the Lupinus species used in this study.
Table 1. General characteristics of the Lupinus species used in this study.
GroupSectionSpeciesAccessionChromosome Number (2n)Genome Size (pg/2C DNA)
Smooth-seededAngustifoliusL. angustifoliuscv. ‘Sonet’ *401.89
L. cryptanthus96361 *401.86
LuteusL. luteuscv. ‘Talar’ *522.44
L. hispanicus96385 *522.15
AlbusL. albuscv. ‘Boros’ *501.16
L. graecus95601 *501.13
MicranthusL. micranthus98552 *520.98
Rough-seededAtlanticusL. atlanticus98401 *381.61
L. digitatusPI 660697 **361.37
L. cosentinii98452 *321.42
PilosusL. pilosus98653 *421.36
L. palaestinus98605 *421.39
* Polish Lupinus Gene Bank, Breeding Station Wiatrowo, Poznan Plant Breeders Ltd., Poland ** US Department of Agriculture, USA.
Table 2. Bacterial artificial chromosome (BAC) clones used for the chromosome and pseudomolecule integration in L. angustifolius.
Table 2. Bacterial artificial chromosome (BAC) clones used for the chromosome and pseudomolecule integration in L. angustifolius.
AssignmentChromosome (Lang)Pseudomolecule (NLL)Pseudomolecule (bp)BAC IDBAC DataGenBank Accession
Lang0101011..36 457 58143C18WBSMK650088
13
Lang0202021..24 697 652120E23WBSMK650090
Lang0303031..30 153 01984A06WBSMK650076
0303 57J20WBSHE804810.1
0303 28O013′ STSGF112056.1
0303 137N08WBSMK650084
Lang0404041..27 333 97547P22WBSMK650072
0404 131K225′ BESKU678257.1
0404 3′ BESKU678256.1
Lang0505051..26376911131K15WBSMK650081
05054G15WBSHE804808.1
Lang0606061..33 111 45044J16WBSKX298066
0606 76K16WBSMK650073
0606 80B11WBSHE804812.1
0606 51D03WBSMK650071
0606 127N17WBSKU678223.1
Lang07Not assigned071..19 782 1702B03WBSKX298069
Lang0808171...25 521 64684D22WBSKX298065
0817 111B08WBSKX298071
0817 142C04WBSKX298073
0817 142D13WBSKX298074
Lang0909091..21 753 71215P085′ BESAB809174.1
09 3′ BESAB809173.1
0909 59J08WBSMK650082
Lang1010101...16 341 95577C13WBSMK650077
1010 57K22WBSMK650069
Lang1111111..35 963 95897D16WBSMK650085
1111 36L23WBSMK650075
1111 20A06WBSMK650080
1111 60F02WBSMK045265
Lang1212121..19 065 70194P05WBSMK650086
Lang13Not analyzedMany1..17 820 68011G20WBSMK650079
Not assigned13 8C03WBSKX298063
Not analyzed13 51F15WBSMK045272
Lang1414031...16 251 777138N02WBSMK650089
1414 115C215′ BESHR864196.1
12 3′ BESAB809318.1
Not analyzed14 9K06WBSMK045273
Lang1515151...20 964 7035L113′ STSGF112057.1
15 134F01WBSMK650070
Lang1616Many1..20 786 881115G225′ BESAB811328.1
Many 3′ BESAB811323.1
16Many 112E015′ BESAB809308.1
Many 3′ BESAB809307.1
16Many 72O215′ BESAB809272.1
1620 8A035′ BESAB809167.1
Many 3′ BESAB809166.1
1616 87N22WBSMK650087
Lang1717171..21 299 880111G03WBSKX298064
1717 136C16WBSKX298072
1717 3B18WBSKX298070
Lang1818181..16 588 00768H10WBSKU678221
Lang1919191..18 159 81267C07WBSMK650078
1919 83F235′ BESKU678304.1
19 3′ BESKU678303.1
0719 74I10WBSMK650074
Lang2020201..21 988 59017B07WBSHF937076.1
2020 83C06WBSMK650083
Cluster-220 1M23WBSKX298067
BAC: bacterial artificial chromosome; BES: BAC-end sequences; WBS: whole BAC sequences; in bold when provided in this study. Lang, NLL, results of chromosome, and pseudomolecule integration, respectively, with reassigned chromosome numbers in italics.
Table 3. Characteristics and localizations of the BAC clones used for the comparative BAC-FISH.
Table 3. Characteristics and localizations of the BAC clones used for the comparative BAC-FISH.
No.BAC IDLcryLlutLhisLalbLgraLmicLatlLdigLcosLpilLpal
Lang011-R1-R1-R1-R1-R1-R1-R1-R1-R1-R1-R
143C18RRRRRRRRRRR
Lang02SSS2-R2-RSS2-RSSS
2120E23SSSRRSSRSSS
Lang033=4=6; 5-R3-4-5-6-R3-4-5-6-R3=4=6; 5-R3=4; 5-6-R3=6; 4-5-R3-4-5-6-R3-5-6-R3=6; 4-5-R3=6; 4-5-R3=4=6; 5-R
384A06SRRSSSRRSSS
457J20SRRSSRRN/DRRS
528O01RRRRRRRRRRR
6137N08SRRSRSRRSSS
Lang047=87-8-R7=87; 8-R7=87; 8-R7-8-R7-8-R7; 8-R7; 8-R7=8
747P22SRSSSSRRSSS
8131K22SRSRSRRRRRS
Lang059=109=109-10-R9=109=109=109=109=109=109=109=10
9131K15SSRSSSSSSSS
104G15SSRSSSSSSSS
Lang0611=12=13=15=1611=13 #14; 12-15-R13; 11-12-14-15-R12=13 # 14; 11-15-R12=13 # 14; 11-15-R12=13 # 14=1511=12=13=15; 14-R12=13; 11-14-15-R13; 12-14-15-R12 # 13 # 14=1512=13#12′; 11-14-15-R
1144J16SSRRRN/DSRN/DN/DR
1276K16SRRSSSSSRSS, S′
1380B11SSSSSSSSSSS
1451D03SSRSSSRRRSR
15127N17SRRRN/DSSRRSR
Lang07SSSSS16-RS16-RSSS
162B03SSSSSRSRSSS
Lang0817=18=19=2018=20; 17-19-R17-18-19-20-R17-18-19-20-R17-18-19-20-R17=19 # 18; 20-R17-18-19-20-R17-18-19-20-R17=18; 19-20-R17=18; 19-20-R17-18-19-20-R
1784D22SRRRRSRRSSR
18111B08SSRRRSRRSSR
19142C04SRRRRSRRRRR
20142D13SSRRRRRRRRR
Lang0922; 21-R22; 21-R22; 21-R22; 21-R22; 21-R22; 21-R22; 21-R22; 21-R22; 21-R22; 21-R22; 21-R
2115P08RRRRRRRRRRR
2259J08SSSSSSSSSSS
Lang1023=2423; 24-R23-24-R23; 24-R23-24-R23=2423; 24-R23; 24-R23-24-R23-24-R23-24-R
2377C13SSRSRSSSRRR
2457K22SRRRRSRRRRR
Lang1125=26=27=2826; 25-27-28-R25=26; 27-28-R25=27; 26-28-R25=26; 27-28-R26; 25-27-28-R25; 26-27-28-R26-27-28-R26; 25-27-28-R26; 25-27-28-R25; 27-28-R
2597D16SRSSSRSN/DRRS
2636L23SSSRSSRRSSN/D
2720A06SRRSRRRRRRR
2860F02SRRRRRRRRRR
Lang1229-RS29-RSS29-R29-R29-R29-R29-RS
2994P05RSRSSRRRRRS
Lang1331=32; 30-R32 # 32′; 30-31-R32; 30-31-R32; 30-31-R32; 30-31-R31; 30-32-R32; 30-31-R32; 30-31-R32; 30-31-R32; 30-31-R30-31-32-R
3011G20RRRRRRRRRRR
318C03SRRRRSRN/DRRR
3251F15SS, S′SSSRSSSSR
Lang1433=35; 34-R33-34-35-R33-34-35-R33-34-35-R35; 33-34-R35; 33-34-R35; 33-34-R33-34-35-R35; 33-34-R35; 33-34-R33-34-35-R
33138N02SRRRRRRRRRR
34115C21RRRRRRRRRRR
359K06SRRRSSSRSSR
Lang1537; 36-R36-37-R36-37-R36-37-R36-37-R37; 36-R36-37-R37; 36-R37; 36-R37; 36-R36-37-R
365L11RRRRRRRRRRR
37134F01SRRRRSRSSSR
Lang1638=39=40=41; 42-R38=39=40=41; 42-R38=39=40=41; 42-R38=39=40=41; 42-R38=39=40=41; 42-R38=39=40=41; 42-R38=39=40=41; 42-R38=39=40=41; 42-R38=39=40=41; 42-R38=39=40=41; 42-R38=39=40=41; 42-R
38115G22SSSSSSSSSSS
39112E01SSSSSSSSSSS
4072O21SSSSSSSSSSS
418A03SSSSSSSSSSS
4287N22RRRRRRRRRRR
Lang1743=45; 44-R43=44=4543=44=4543=44=4543=44=4543=44=4543=44=4543 # 44=45 # 45′43=44=4543 # 44=4544=45; 43-R
433B18SSSSSSSSSSR
44111G03RSSSSSSSSSS
45136C16SSSSSSSS, S′SSS
Lang1846-R46-R46-R46-R46-R46-R46-R46-R46-R46-R46-R
4668H10RRRRRRRRRRR
Lang1947=48; 49-R47; 48-49-R47=48; 49-R47-49-R47-48-49-R47=48; 49-R47=48=4947=48; 49-R47=48; 49-R47=48; 49-R47=48; 49-R
4767C07SSSRRSSSSSS
4883F23SRSN/DRSSSSSS
4974I10 RRRRRRSRRRR
Lang2050=51=5250=51=5250; 51-52-R50=51; 52-R50; 51-52-R50=51; 52-R50=51; 52-R50=51; 52-R50=51; 52-R50=51; 52-R50=51; 52-R
5017B07SSSSSSSSSSS
5183C06SSRSRSSSSSS
521M23SSRRRRRRRRR
BAC clones are individually numbered from 1 to 52 and grouped with correspondence to L. angustifolius chromosomes (Lang). There are four possible BAC-FISH mapping results: S, ‘single’ BAC; S, S′, BAC mapped into two loci of two chromosomes; R, ‘repetitive’ BAC; N/D, not detected. For example: (a) BAC clone no. 2—120E23 corresponding to Lang02 is ‘single’ (S) in Lcry, but ‘repetitive’ (R) in Lalb; (b) 3=4; 5-6-R means that BACs 84A06 and 57J20 were mapped as ‘single’ in the same chromosome, but both clones 28O01 and 137N08 were ‘repetitive’; (c) 11=13 #14; 12-15-R means that BACs 11 and 13 were mapped in the same chromosome, while BAC 14 was mapped in different chromosome, both BACs 12 and 15 were ‘repetitive’.

Share and Cite

MDPI and ACS Style

Susek, K.; Bielski, W.; Czyż, K.B.; Hasterok, R.; Jackson, S.A.; Wolko, B.; Naganowska, B. Impact of Chromosomal Rearrangements on the Interpretation of Lupin Karyotype Evolution. Genes 2019, 10, 259. https://doi.org/10.3390/genes10040259

AMA Style

Susek K, Bielski W, Czyż KB, Hasterok R, Jackson SA, Wolko B, Naganowska B. Impact of Chromosomal Rearrangements on the Interpretation of Lupin Karyotype Evolution. Genes. 2019; 10(4):259. https://doi.org/10.3390/genes10040259

Chicago/Turabian Style

Susek, Karolina, Wojciech Bielski, Katarzyna B. Czyż, Robert Hasterok, Scott A. Jackson, Bogdan Wolko, and Barbara Naganowska. 2019. "Impact of Chromosomal Rearrangements on the Interpretation of Lupin Karyotype Evolution" Genes 10, no. 4: 259. https://doi.org/10.3390/genes10040259

APA Style

Susek, K., Bielski, W., Czyż, K. B., Hasterok, R., Jackson, S. A., Wolko, B., & Naganowska, B. (2019). Impact of Chromosomal Rearrangements on the Interpretation of Lupin Karyotype Evolution. Genes, 10(4), 259. https://doi.org/10.3390/genes10040259

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop