Bacterial Artificial Chromosome Clones Research Articles

Improved protocols for BAC insert DNA isolation, BAC end sequencing and FISH for construction of BAC based physical map of genes on the chromosomes.

Bacterial artificial chromosome (BAC) library is an important genomic resource useful in targeted marker development, positional cloning, physical mapping and a substrate for genome sequencing for better understanding the genome organization of a species. The present manuscript elucidates the improvement in protocols for economical and efficient BAC insert DNA isolation, BAC end sequencing and FISH for physical localization on the metaphase chromosome complements. BAC clones of Clarias magur, maintained in 384-well plate format in our laboratory, were used in this study. The protocols gave consistent and efficient results. We use routinely these protocols for BAC insert DNA extraction, generating end sequence data of the clone and constructing DNA probes to hybridize on the metaphase spreads of C. magur using FISH for physical their localization.

BAC selection of the large white butterfly, Pieris brassicae (Lepidoptera: Pieridae) containing orthologs of the silkworm, Bombyx mori (Lepidoptera: Bombycidae)

We constructed a bacterial artificial chromosome (BAC) library from the large white butterfly, Pieris brassicae (L.) (Lepidoptera: Pieridae), which composed of 15,360 clones with an average insert size of 103 kb. We designed 39 degenerate oligonucleotide primer pairs and cloned the PCR products for sequencing. From these plus pierisin data, we designed single tagged sequence primers and selected 40 BAC clones using a 3-step PCR protocol. Additional 40 BACs were sequenced to improve identification of P. brassicae BACs carrying orthologs of the silkworm, Bombyx mori (L.) (Lepidoptera: Bombycidae). Among 34 BAC clones containing multiple B. mori orthologs, sequences in 27 BACs correspond to a block (less than 1 Mb) of every B. mori chromosome. While orthologs in 5 BACs showed collinearity to each B. mori chromosome with a distance (1.4–6.8 Mb), orthologs in two BACs corresponded to genes in more than two different B. mori chromosomes.

Fluorescent In Situ Hybridization Using Oligonucleotide-Based Probes.

Efficient and consistent chromosome identification is the foundation for successful cytogenetic studies. Fluorescent in situ hybridization (FISH) has been the most popular technique for chromosome identification in plants. Large insert genomic DNA clones, such as bacterial artificial chromosome (BAC) clones, and repetitive DNA sequences have been the most commonly used DNA probes for FISH. However, most of such traditional probes can only be used to identify a single chromosome or are too polymorphic to consistently identify the same chromosome in the target species. In contrast, FISH using oligonucleotide (oligo)-based probes is highly versatile. In this procedure, a large number of oligos specific to a chromosomal region, to an entire chromosome, or to multiple chromosomes are computationally identified, synthesized in parallel, and labeled as probes. In addition, each oligo probe can be used for thousands of FISH experiments and represents an infinite resource. In this chapter we describe a detailed protocol for amplification and labeling of oligo-based probes, relevant chromosome preparation, and FISH procedures.

Fine mapping of a major locus representing the lack of prickles in eggplant revealed the availability of a 0.5-kb insertion/deletion for marker-assisted selection.

As prickles cause labour inefficiency during cultivation and scratches on the skin of fruits during transportation, they are considered undesirable traits of eggplant (Solanum melongena L.). Because the molecular basis of prickle emergence has not been entirely revealed in plants, we mapped an eggplant semi-dominant Prickle (Pl) gene locus, which causes the absence of prickles, on chromosome 6 of a linkage map of the F2 population derived from crossing the no-prickly cultivar ‘Togenashi-senryo-nigo’ and the prickly line LS1934. By performing synteny mapping with tomato, the genomic region corresponding to the eggplant Pl locus was identified. Through bacterial artificial chromosome (BAC) screening, positive BAC clones and the contig sequence that harbour the Pl locus in the prickly eggplant genome were revealed. The BAC contig length was 133 kb, and it contained 16 predicted genes. Among them, a characteristic 0.5-kb insertion/deletion was detected. As the 0.5-kb insertion was commonly identified with the prickly phenotype worldwide, a primer pair that amplifies the insertion/deletion could be used for marker-assisted selection of the no-prickly phenotype. Such findings contribute to map-based-cloning of the Pl gene and the understanding of gene function, ultimately providing new insights into the regulatory molecular mechanisms underlying prickle emergence in plants.

The Genome Resequencing of TCR Loci in Gallus gallus Revealed Their Distinct Evolutionary Features in Avians.

The TCR is consisted of four chains: α (TCRα), β (TCRβ), γ (TCRγ), and δ (TCRδ) that are present in all jawed vertebrates. Birds are very important in terms of evolutionary aspects of the adaptive immune system, in which it bridges the evolutionary gap between mammals and other vertebrates. To gain better understanding into the genomic organization and complexity of birds' TCR loci, we applied cross-reference error-correction sequencing approach by using Illumina and single-molecule real-time sequencing technology to resequence genomic regions of chicken TCR loci based on 10 mapped bacterial artificial chromosome clones. We did de novo classification of V and J genes for all four chains of the TCR loci according to our sequencing results using the Immunogenetics nomenclature. In sum, we identified 85, 8, and 37 TCR V gene segments in the chicken TCRα/TCRδ, TCRβ, and TCRγ loci, respectively. The phylogenetic analysis showed the Vα 7 and Vα family 4 gene sequences shared greater sequence similarity with mammalian species, whereas the other Vα segment sequences are evolutionary closer with sequences from bony fishes. The organization of chicken TCRβ locus is more similar to fish TCRβ locus over mammalian species, as chicken TCRβ locus has a single translocon of its V-D-J-C and exhibits significantly fewer Vβ gene segments. In this study, we present a highly precise genomic map for chicken TCR loci and phylogenetic relationships of TCR variable gene segments against other animal species and verified the relative stability of the receptor structure during evolutional process.

Resistance to Sharka in Apricot: Comparison of Phase-Reconstructed Resistant and Susceptible Haplotypes of 'Lito' Chromosome 1 and Analysis of Candidate Genes.

Sharka, a common disease among most stone fruit crops, is caused by the Plum Pox Virus (PPV). Resistant genotypes have been found in apricot (Prunus armeniaca L.), one of which—the cultivar ‘Lito’ heterozygous for the resistance—has been used to map a major quantitative trait locus (QTL) on linkage group 1, following a pseudo-test-cross mating design with 231 individuals. In addition, 19 SNP markers were selected from among the hundreds previously developed, which allowed the region to be limited to 236 kb on chromosome 1. A ‘Lito’ bacterial artificial chromosome (BAC) library was produced, screened with markers of the region, and positive BAC clones were sequenced. Resistant (R) and susceptible (S) haplotypes were assembled independently. To refine the assembly, the whole genome of ‘Lito’ was sequenced to high coverage (98×) using PacBio technology, enabling the development of a detailed assembly of the region that was able to predict and annotate the genes in the QTL region. The selected cultivar ‘Lito’ allowed not only to discriminate structural variants between the two haplotypic regions but also to distinguish specific allele expression, contributing towards mining the PPVres locus. In light of these findings, genes previously indicated (i.e., MATHd genes) to have a possible role in PPV resistance were further analyzed, and new candidates were discussed. Although the results are not conclusive, the accurate and independent assembly of R and S haplotypes of ‘Lito’ is a valuable resource to predict and test alternative transcription and regulation mechanisms underpinning PPV resistance.

Comparatively Barcoded Chromosomes of Brachypodium Perennials Tell the Story of Their Karyotype Structure and Evolution.

The Brachypodium genus is an informative model system for studying grass karyotype organization. Previous studies of a limited number of species and reference chromosomes have not provided a comprehensive picture of the enigmatic phylogenetic relationships in the genus. Comparative chromosome barcoding, which enables the reconstruction of the evolutionary history of individual chromosomes and their segments, allowed us to infer the relationships between putative ancestral karyotypes of extinct species and extant karyotypes of current species. We used over 80 chromosome-specific BAC (bacterial artificial chromosome) clones derived from five reference chromosomes of B. distachyon as probes against the karyotypes of twelve accessions representing five diploid and polyploid Brachypodium perennials. The results showed that descending dysploidy is common in Brachypodium and occurs primarily via nested chromosome fusions. Brachypodium distachyon was rejected as a putative ancestor for allotetraploid perennials and B. stacei for B. mexicanum. We propose two alternative models of perennial polyploid evolution involving either the incorporation of a putative x = 5 ancestral karyotype with different descending dysploidy patterns compared to B. distachyon chromosomes or hybridization of two x = 9 ancestors followed by genome doubling and descending dysploidy. Details of the karyotype structure and evolution in several Brachypodium perennials are revealed for the first time.

Small-Scale Isolation of Bacterial Artificial Chromosome (BAC) DNA and Verification by Polymerase Chain Reaction (PCR).

In this protocol, small amounts of selected bacterial artificial chromosome (BAC) DNAs are prepared from 5-mL cultures of the Escherichia coli host transformed with the BAC clone. DNA is isolated by an adaptation of the alkaline lysis method. The yield is 0.1-0.4 µg, and the BAC DNA is suitable for analysis by polymerase chain reaction (PCR). Considerations for clone selection are also addressed.

ZW Sex Chromosomes in Australian Dragon Lizards (Agamidae) Originated from a Combination of Duplication and Translocation in the Nucleolar Organising Region.

Sex chromosomes in some reptiles share synteny with distantly related amniotes in regions orthologous to squamate chromosome 2. The latter finding suggests that chromosome 2 was formerly part of a larger ancestral (amniote) super-sex chromosome and raises questions about how sex chromosomes are formed and modified in reptiles. Australian dragon lizards (Agamidae) are emerging as an excellent model for studying these processes. In particular, they exhibit both genotypic (GSD) and temperature-dependent (TSD) sex determination, show evidence of transitions between the two modes and have evolved non-homologous ZW sex microchromosomes even within the same evolutionary lineage. They therefore represent an excellent group to probe further the idea of a shared ancestral super-sex chromosome and to investigate mechanisms for transition between different sex chromosome forms. Here, we compare sex chromosome homology among eight dragon lizard species from five genera to identify key cytological differences and the mechanisms that may be driving sex chromosome evolution in this group. We performed fluorescence in situ hybridisation to physically map bacterial artificial chromosome (BAC) clones from the bearded dragon, Pogona vitticeps’ ZW sex chromosomes and a nucleolar organising region (NOR) probe in males and females of eight Agamid species exhibiting either GSD or TSD. We show that the sex chromosome derived BAC clone hybridises near the telomere of chromosome 2q in all eight species examined. This clone also hybridises to the sex microchromosomes of three species (P vitticeps, P. barbata and Diporiphora nobbi) and a pair of microchromosomes in three others (Ctenophorus pictus, Amphibolurus norrisi and Amphibolurus muricatus). No other chromosomes are marked by the probe in two species from the closely related genus Physignathus. A probe bearing nucleolar organising region (NOR) sequences maps close to the telomere of chromosome 2q in all eight species, and to the ZW pair in P. vitticeps and P. barbata, the W microchromosome in D. nobbi, and several microchromosomes in P. cocincinus. Our findings provide evidence of sequence homology between chromosome 2 and the sex chromosomes of multiple agamids. These data support the hypothesis that there was an ancestral sex chromosome in amniotes that gave rise to squamate chromosome 2 and raises the prospect that some particular property of this chromosome has favoured its role as a sex chromosome in amniotes. It is likely that the amplification of repetitive sequences associated with this region has driven the high level of heterochromatinisation of the sex-specific chromosomes in three species of agamid. Our data suggest a possible mechanism for chromosome rearrangement, including inversion and duplication near the telomeric regions of the ancestral chromosome 2 and subsequent translocation to the ZW sex microchromosomes in three agamid species. It is plausible that these chromosome rearrangements involving sex chromosomes also drove speciation in this group.

Lipopolysaccharide-induced expansion of histidine decarboxylase-expressing Ly6G+ myeloid cells identified by exploiting histidine decarboxylase BAC-GFP transgenic mice

Histamine is a biogenic amine that is chiefly produced in mast cells and basophils and elicits an allergic response upon stimulation. Histidine decarboxylase (HDC) is a unique enzyme that catalyzes the synthesis of histamine. Therefore, the spatiotemporally specific Hdc gene expression profile could represent the localization of histamine-producing cells under various pathophysiological conditions. Although the bioactivity of histamine is well defined, the regulatory mechanism of Hdc gene expression and the distribution of histamine-producing cell populations in various disease contexts remains unexplored. To address these issues, we generated a histidine decarboxylase BAC (bacterial artificial chromosome) DNA-directed GFP reporter transgenic mouse employing a 293-kb BAC clone containing the entire Hdc gene locus and extended flanking sequences (Hdc-GFP). We found that the GFP expression pattern in the Hdc-GFP mice faithfully recapitulated that of conventional histamine-producing cells and that the GFP expression level mirrored the increased Hdc expression in lipopolysaccharide (LPS)-induced septic lungs. Notably, a CD11b+Ly6G+Ly6Clow myeloid cell population accumulated in the lung during sepsis, and most of these cells expressed high levels of GFP and indeed contain histamine. This study reveals the accumulation of a histamine-producing myeloid cell population during sepsis, which likely participates in the immune process of sepsis.

An Improved Recombineering Toolset for Plants.

Gene functional studies often rely on the expression of a gene of interest as transcriptional and translational fusions with specialized tags. Ideally, this is done in the native chromosomal contexts to avoid potential misexpression artifacts. Although recent improvements in genome editing have made it possible to directly modify the target genes in their native chromosomal locations, classical transgenesis is still the preferred experimental approach chosen in most gene tagging studies because of its time efficiency and accessibility. We have developed a recombineering-based tagging system that brings together the convenience of the classical transgenic approaches and the high degree of confidence in the results obtained by direct chromosomal tagging using genome-editing strategies. These simple, scalable, customizable recombineering toolsets and protocols allow a variety of genetic modifications to be generated. In addition, we developed a highly efficient recombinase-mediated cassette exchange system to facilitate the transfer of the desired sequences from a bacterial artificial chromosome clone to a transformation-compatible binary vector, expanding the use of the recombineering approaches beyond Arabidopsis (Arabidopsis thaliana). We demonstrated the utility of this system by generating more than 250 whole-gene translational fusions and 123 Arabidopsis transgenic lines corresponding to 62 auxin-related genes and characterizing the translational reporter expression patterns for 14 auxin biosynthesis genes.

Characterization and fine mapping of Arabidopsis RLL3, a locus required for DNA demethylation pathway.

DNA methylation pattern is found to be established by the combined actions of DNA methylation and demethylation. Compared to the DNA methylation pathway, DNA demethylation pathway, however, remains largely unknown. To better understand the DNA demethylation pathway, we performed genetic screening for Arabidopsis mutants with increased genomic DNA methylation levels through a 2 × 35S:LUC (LUC, luciferase) reporter system. A mutant with reduced LUC luminescence was identified by such a system, therefore named rll3-1 (for reduced LUC luminescence 3-1). The rll3-1 mutant exhibited pleiotropic developmental defects, such as delayed bolting as well as flowering, more branches, etc. By map-based cloning approach, rll3 locus that contains a single nuclear recessive mutation as revealed by the genetic analysis was mapped to a region between molecular markers CL102_B1 M1 and CL102_B3M1, which are located in bacterial artificial chromosome (BAC) clones F9P14 and F12K11, respectively, on chromosome 1. Chop-PCR analysis indicated that a total of seven tested loci displayed elevated DNA methylation levels. Whole-genome bisulfite sequencing further revealed 1536 loci exhibiting increased DNA methylation levels relative to Col-LUC control, among which there are 507 such loci overlapping between the rll3-1 and ros1-7 mutants, suggestive of a functional association between RLL3 and REPRESSOR OF SILENCING 1 (ROS1). Further investigations demonstrated that the expression levels of a few genes (like ROS1, IDM1, etc.), which are involved in DNA demethylation pathway, remained unchanged in the rll3-1 mutant, indicating that the increased DNA methylation levels in rll3-1 mutant are not attributable to downregulation of such genes. Taken together, our studies provide a demonstration of the involvement of RLL3 in the DNA demethylation pathway.

Optimized Expression of Duck Tembusu Virus E Gene Delivered by a Vectored Duck Enteritis Virus In Vitro.

In our previous study, a recombinant duck enteritis virus (DEV) delivering codon-optimized E gene (named as E-ch) of duck Tembusu virus (DTMUV) optimized referring to chicken's codon bias has been obtained based on the infectious bacterial artificial chromosome (BAC) clone of duck enteritis virus vaccine strain pDEV-EF1, but the expression level of E-ch in recombinant virus rDEV-E-ch-infected cells was very low. To optimize DTMUV E gene expression delivered by the vectored DEV, different forms of E gene (collectively called EG) including origin E gene (E-ori), truncated E451-ori gene, codon-optimized E-dk gene optimized referring to duck's codon bias, as well as the truncated E451-ch and E451-dk, Etpa-ori and Etpa-451-ori, which contain prefixing chick TPA signal peptide genes, were cloned into transfer vector pEP-BGH-end, and several recombinant plasmids pEP-BGH-EG were constructed. Then the expression cassettes pCMV-EG-polyABGH amplified from pEP-BGH-EG by PCR were inserted into US7/US8 gene intergenicregion of pDEV-EF1 by two-step Red/ET recombination, 7 strain recombinant mutated BAC clones pDEV-EG carrying different E genes were constructed. Next, the recombinant viruses rDEV-EG were reconstituted from chicken embryo fibroblasts (CEFs) by calcium phosphate precipitation. Western blot analysis showed that E or E451 protein is expressed in rDEV-E-ori, rDEV-E-ch, rDEV-Etpa-ori, rDEV-E451-ori, rDEV-E451-dk, and rDEV-E451-ch-infected CEFs, and protein expression level in rDEV-E451-dk-infected CEFs is the highest. These studies have laid a foundation for developing bivalent vaccine controlling DEV and DTMUV infection.

Identification of passion fruit (Passiflora edulis) chromosomes using BAC-FISH.

Passiflora edulis, the yellow passion fruit, is the main crop from the Passiflora genus, which comprises 525 species with its diversity center in South America. Genetic maps and a BAC (bacterial artificial chromosome) genomic library are available, but the nine chromosome pairs of similar size and morphology (2n = 18) hamper chromosome identification, leading to different proposed karyotypes. Thus, the aim of this study was to establish chromosome-specific markers for the yellow passion fruit using single-copy and repetitive sequences as probes in fluorescent in situ hybridizations (FISH) to allow chromosome identification and future integration with whole genome data. Thirty-six BAC clones harboring genes and three retrotransposons (Ty1-copy, Ty3-gypsy, and LINE) were selected. Twelve BACs exhibited a dispersed pattern similar to that revealed by retroelements, and one exhibited subtelomeric distribution. Twelve clones showed unique signals in terminal or subterminal regions of the chromosomes, allowing their genes to be anchored to six chromosome pairs that can be identified with single-copy markers. The markers developed herein will provide an important tool for genomic and evolutionary studies in the Passiflora genus.

Fine mapping and identification of the rice blast-resistance locus Pi-kf2(t) as a new member of the Pi2/Pi9 multigene family

The identification and utilization of broad-spectrum resistance genes is the most effective and economical strategy of controlling rice blast disease. Kangfeng B (KFB), an elite Chinese rice cultivar, has been shown to exhibit broad-spectrum resistance to 53 isolates of Magnaporthe oryzae, the rice blast causative agent, from different regions of China. In this study, we identified a dominant blast-resistance gene at the Pi2/Pi9 locus from cultivar KFB, through genetic analysis and physical mapping. Allele-specific marker-based assessment revealed that Pi2, Pi9, and Piz-t are not the blast-resistance genes in KFB. By combining bulked segregation analysis (BSA) and recessive class analysis (RCA), the blast-resistance gene was fine-mapped to an approximately 249-kb interval between markers InDel-22 and Rm7213 on chromosome 6. Three bacterial artificial chromosome (BAC) clones spanning the region were identified. This region contains 19 predicted genes, including 7 nucleotide binding site–leucine-rich repeat (NBS-LRR) genes at the Pi2/Pi9 locus in japonica cv. Nipponbare genome. Further sequence comparison of the four functional NBS-LRR genes revealed that NBS-LRR2 and NBS-LRR4, as evidenced by their allelic/orthologous to Pi9 or Pi2, had significant differences of 9 to 43 and 14 to 48 amino acids in KFB, respectively, unlike the other known Pi2/Pi9 alleles, suggesting that KFB carries a hitherto undocumented member of the Pi2/Pi9 multigene family. It was tentatively designated as Pi-kf2(t). Our results provide essential information for the isolation of the Pi-kf2(t) gene and will facilitate both map-based cloning and marker-assisted selection of this gene for rice blast-resistance breeding.

Genome and Phylogenetic Analysis of Genes Involved in the Immune System of Solea senegalensis - Potential Applications in Aquaculture.

Global aquaculture production continues to increase rapidly. One of the most important species of marine fish currently cultivated in Southern Europe is Solea senegalensis, reaching more than 300 Tn in 2017. In the present work, 14 Bacterial Artificial Chromosome (BAC) clones containing candidate genes involved in the immune system (b2m, il10, tlr3, tap1, tnfα, tlr8, trim25, lysg, irf5, hmgb2, calr, trim16, and mx), were examined and compared with other species using multicolor Fluorescence in situ Hybridization (mFISH), massive sequencing and bioinformatic analysis to determine the genomic surroundings and syntenic chromosomal conservation of the genomic region contained in each BAC clone. The mFISH showed that the groups of genes hmgb2-trim25-irf5-b2m; tlr3-lysg; tnfα-tap1, and il10-mx-trim16 were co-localized on the same chromosomes. Synteny results suggested that the studied BACs are placed in a smaller number of chromosomes in S. senegalensis that in other species. Phylogenetic analyses suggested that the evolutionary rate of immune system genes studied is similar among the taxa studied, given that the clustering obtained was in accordance with the accepted phylogenetic relationships among these species. This study contributes to a better understanding of the structure and function of the immune system of the Senegalese sole, which is essential for the development of new technologies and products to improve fish health and productivity.

A high-density BAC physical map covering the entire MHC region of addax antelope genome

BackgroundThe mammalian major histocompatibility complex (MHC) harbours clusters of genes associated with the immunological defence of animals against infectious pathogens. At present, no complete MHC physical map is available for any of the wild ruminant species in the world.ResultsThe high-density physical map is composed of two contigs of 47 overlapping bacterial artificial chromosome (BAC) clones, with an average of 115 Kb for each BAC, covering the entire addax MHC genome. The first contig has 40 overlapping BAC clones covering an approximately 2.9 Mb region of MHC class I, class III, and class IIa, and the second contig has 7 BAC clones covering an approximately 500 Kb genomic region that harbours MHC class IIb. The relative position of each BAC corresponding to the MHC sequence was determined by comparative mapping using PCR screening of the BAC library of 192,000 clones, and the order of BACs was determined by DNA fingerprinting. The overlaps of neighboring BACs were cross-verified by both BAC-end sequencing and co-amplification of identical PCR fragments within the overlapped region, with their identities further confirmed by DNA sequencing.ConclusionsWe report here the successful construction of a high-quality physical map for the addax MHC region using BACs and comparative mapping. The addax MHC physical map we constructed showed one gap of approximately 18 Mb formed by an ancient autosomal inversion that divided the MHC class II into IIa and IIb. The autosomal inversion provides compelling evidence that the MHC organizations in all of the ruminant species are relatively conserved.

Accessing a Russian Wheat Aphid Resistance Gene in Bread Wheat by Long-Read Technologies.

Russian wheat aphid (RWA) ( Kurdjumov) is a serious invasive pest of small-grain cereals and many grass species. An efficient strategy to defy aphid attacks is to identify sources of natural resistance and transfer resistance genes into susceptible crop cultivars. Revealing the genes helps understand plant defense mechanisms and engineer plants with durable resistance to the pest. To date, more than 15 RWA resistance genes have been identified in wheat ( L.) but none of them has been cloned. Previously, we genetically mapped the RWA resistance gene into an interval of 0.83 cM on the short arm of chromosome 7D and spanned it with five bacterial artificial chromosome (BAC) clones. Here, we used a targeted strategy combining traditional approaches toward gene cloning (genetic mapping and sequencing of BAC clones) with novel technologies, including optical mapping and long-read nanopore sequencing. The latter, with reads spanning the entire length of a BAC insert, enabled us to assemble the whole region, a task that was not achievable with short reads. Long-read optical mapping validated the DNA sequence in the interval and revealed a difference in the locus organization between resistant and susceptible genotypes. The complete and accurate sequence of the region facilitated the identification of new markers and precise annotation of the interval, revealing six high-confidence genes. Identification of as the most likely candidate opens an avenue for its validation through functional genomics approaches.

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.

Development of a Sequence-Based Reference Physical Map of Pea (Pisum sativum L.).

Whole genome profiling (WGP) is a sequence-based physical mapping technology and uses sequence tags generated by next generation sequencing for construction of bacterial artificial chromosome (BAC) contigs of complex genomes. The physical map provides a framework for assembly of genome sequence and information for localization of genes that are difficult to find through positional cloning. To address the challenges of accurate assembly of the pea genome (∼4.2 GB of which approximately 85% is repetitive sequences), we have adopted the WGP technology for assembly of a pea BAC library. Multi-dimensional pooling of 295,680 BAC clones and sequencing the ends of restriction fragments of pooled DNA generated 1,814 million high quality reads, of which 825 million were deconvolutable to 1.11 million unique WGP sequence tags. These WGP tags were used to assemble 220,013 BACs into contigs. Assembly of the BAC clones using the modified Fingerprinted Contigs (FPC) program has resulted in 13,040 contigs, consisting of 213,719 BACs, and 6,294 singleton BACs. The average contig size is 0.33 Mbp and the N50 contig size is 0.62 Mbp. WGPTM technology has proved to provide a robust physical map of the pea genome, which would have been difficult to assemble using traditional restriction digestion based methods. This sequence-based physical map will be useful to assemble the genome sequence of pea. Additionally, the 1.1 million WGP tags will support efficient assignment of sequence scaffolds to the BAC clones, and thus an efficient sequencing of BAC pools with targeted genome regions of interest.