Development of flow cytogenetics and physical genome mapping in chickpea (Cicer arietinum L.)

Chromosome analysis and sorting

Flow cytometric analysis and sorting of plant mitotic chromosomes has been mastered by only a few laboratories worldwide. Yet, it has been contributing significantly to progress in plant genetics, including the production of genome assemblies and the cloning of important genes. The dissection of complex genomes by flow sorting into the individual chromosomes that represent small parts of the genome reduces DNA sample complexity and streamlines projects relying on molecular and genomic techniques. Whereas flow cytometric analysis, that is, chromosome classification according to fluorescence and light scatter properties, is an integral part of any chromosome sorting project, it has rarely been used on its own due to lower resolution and sensitivity as compared to other cytogenetic methods. To perform chromosome analysis and sorting, commercially available electrostatic droplet sorters are suitable. However, in order to resolve and purify chromosomes of interest the instrument must offer high resolution of optical signals as well as stability during long runs. The challenge is thus not the instrumentation, but the adequate sample preparation. The sample must be a suspension of intact mitotic metaphase chromosomes and the protocol, which includes the induction of cell cycle synchrony, accumulation of dividing cells at metaphase, and release of undamaged chromosomes, is time consuming and laborious and needs to be performed very carefully. Moreover, in addition to fluorescent staining chromosomal DNA, the protocol may include specific labelling of DNA repeats to facilitate discrimination of particular chromosomes. This review introduces the applications of chromosome sorting in plants, and discusses in detail sample preparation, chromosome analysis and sorting to achieve the highest purity in flow-sorted fractions, and their suitability for downstream applications.

Karyotype Differentiation in Cultivated Chickpea Revealed by Oligopainting Fluorescence in situ Hybridization

Chickpea (Cicer arietinum L.) is one of the main sources of plant proteins in the Indian subcontinent and West Asia, where two different morphotypes, desi and kabuli, are grown. Despite the progress in genome mapping and sequencing, the knowledge of the chickpea genome at the chromosomal level, including the long-range molecular chromosome organization, is limited. Earlier cytogenetic studies in chickpea suffered from a limited number of cytogenetic landmarks and did not permit to identify individual chromosomes in the metaphase spreads or to anchor pseudomolecules to chromosomes in situ. In this study, we developed a system for fast molecular karyotyping for both morphotypes of cultivated chickpea. We demonstrate that even draft genome sequences are adequate to develop oligo-fluorescence in situ hybridization (FISH) barcodes for the identification of chromosomes and comparative analysis among closely related chickpea genotypes. Our results show the potential of oligo-FISH barcoding for the identification of structural changes in chromosomes, which accompanied genome diversification among chickpea cultivars. Moreover, oligo-FISH barcoding in chickpea pointed out some problematic, most probably wrongly assembled regions of the pseudomolecules of both kabuli and desi reference genomes. Thus, oligo-FISH appears as a powerful tool not only for comparative karyotyping but also for the validation of genome assemblies.

Flow Analysis and Sorting of Plant Chromosomes

Analysis and sorting of plant chromosomes (plant flow cytogenetics) is a special application of flow cytometry in plant genomics and its success depends critically on sample quality. This unit describes the methodology in a stepwise manner, starting with the induction of cell cycle synchrony and accumulation of dividing cells in mitotic metaphase, and continues with the preparation of suspensions of intact mitotic chromosomes, flow analysis and sorting of chromosomes, and finally processing of the sorted chromosomes. Each step of the protocol is described in detail as some procedures have not been used widely. Supporting histograms are presented as well as hints on dealing with plant material; the utility of sorted chromosomes for plant genomics is also discussed. © 2016 by John Wiley & Sons, Inc.

Flow Sorting Plant Chromosomes

Nuclear genomes of many important plant species are tremendously complicated to map and sequence. The ability to isolate single chromosomes, which represent small units of nuclear genome, is priceless in many areas of plant research including cytogenetics, genomics, and proteomics. Flow cytometry is the only technique which can provide large quantities of pure chromosome fractions suitable for downstream applications including physical mapping, preparation of chromosome-specific BAC libraries, sequencing, and optical mapping. Here, we describe step-by-step procedure of preparation of liquid suspensions of intact mitotic metaphase chromosomes and their flow cytometric analysis and sorting.

A chromosomal genomics approach to assess and validate the desi and kabuli draft chickpea genome assemblies

With the expansion of next-generation sequencing technology and advanced bioinformatics, there has been a rapid growth of genome sequencing projects. However, while this technology enables the rapid and cost-effective assembly of draft genomes, the quality of these assemblies usually falls short of gold standard genome assemblies produced using the more traditional BAC by BAC and Sanger sequencing approaches. Assembly validation is often performed by the physical anchoring of genetically mapped markers, but this is prone to errors and the resolution is usually low, especially towards centromeric regions where recombination is limited. New approaches are required to validate reference genome assemblies. The ability to isolate individual chromosomes combined with next-generation sequencing permits the validation of genome assemblies at the chromosome level. We demonstrate this approach by the assessment of the recently published chickpea kabuli and desi genomes. While previous genetic analysis suggests that these genomes should be very similar, a comparison of their chromosome sizes and published assemblies highlights significant differences. Our chromosomal genomics analysis highlights short defined regions that appear to have been misassembled in the kabuli genome and identifies large-scale misassembly in the draft desi genome. The integration of chromosomal genomics tools within genome sequencing projects has the potential to significantly improve the construction and validation of genome assemblies. The approach could be applied both for new genome assemblies as well as published assemblies, and complements currently applied genome assembly strategies.

Accessing complex crop genomes with next-generation sequencing

Many important crop species have genomes originating from ancestral or recent polyploidisation events. Multiple homoeologous gene copies, chromosomal rearrangements and amplification of repetitive DNA within large and complex crop genomes can considerably complicate genome analysis and gene discovery by conventional, forward genetics approaches. On the other hand, ongoing technological advances in molecular genetics and genomics today offer unprecedented opportunities to analyse and access even more recalcitrant genomes. In this review, we describe next-generation sequencing and data analysis techniques that vastly improve our ability to dissect and mine genomes for causal genes underlying key traits and allelic variation of interest to breeders. We focus primarily on wheat and oilseed rape, two leading examples of major polyploid crop genomes whose size or complexity present different, significant challenges. In both cases, the latest DNA sequencing technologies, applied using quite different approaches, have enabled considerable progress towards unravelling the respective genomes. Our ability to discover the extent and distribution of genetic diversity in crop gene pools, and its relationship to yield and quality-related traits, is swiftly gathering momentum as DNA sequencing and the bioinformatic tools to deal with growing quantities of genomic data continue to develop. In the coming decade, genomic and transcriptomic sequencing, discovery and high-throughput screening of single nucleotide polymorphisms, presence-absence variations and other structural chromosomal variants in diverse germplasm collections will give detailed insight into the origins, domestication and available trait-relevant variation of polyploid crops, in the process facilitating novel approaches and possibilities for genomics-assisted breeding.

Chromosomes in the flow to simplify genome analysis

Nuclear genomes of human, animals, and plants are organized into subunits called chromosomes. When isolated into aqueous suspension, mitotic chromosomes can be classified using flow cytometry according to light scatter and fluorescence parameters. Chromosomes of interest can be purified by flow sorting if they can be resolved from other chromosomes in a karyotype. The analysis and sorting are carried out at rates of 10(2)-10(4) chromosomes per second, and for complex genomes such as wheat the flow sorting technology has been ground-breaking in reducing genome complexity for genome sequencing. The high sample rate provides an attractive approach for karyotype analysis (flow karyotyping) and the purification of chromosomes in large numbers. In characterizing the chromosome complement of an organism, the high number that can be studied using flow cytometry allows for a statistically accurate analysis. Chromosome sorting plays a particularly important role in the analysis of nuclear genome structure and the analysis of particular and aberrant chromosomes. Other attractive but not well-explored features include the analysis of chromosomal proteins, chromosome ultrastructure, and high-resolution mapping using FISH. Recent results demonstrate that chromosome flow sorting can be coupled seamlessly with DNA array and next-generation sequencing technologies for high-throughput analyses. The main advantages are targeting the analysis to a genome region of interest and a significant reduction in sample complexity. As flow sorters can also sort single copies of chromosomes, shotgun sequencing DNA amplified from them enables the production of haplotype-resolved genome sequences. This review explains the principles of flow cytometric chromosome analysis and sorting (flow cytogenetics), discusses the major uses of this technology in genome analysis, and outlines future directions.

Integration of genetic and physical maps of the chickpea (Cicer arietinum L.) genome using flow-sorted chromosomes

Cultivated chickpea is the third most important legume after field bean and garden pea worldwide. Despite considerable breeding towards improved yield and resistance to biotic and abiotic stresses, the production of chickpea remained stagnant, but molecular tools are expected to increase the impact of current improvement programs. As a first step towards this goal, various genetic linkage maps have been established and markers linked to resistance genes been identified. However, until now, only one linkage group (LG) has been assigned to a specific chromosome. In the present work, mitotic chromosomes were sorted using flow cytometry and used as template for PCR with primers designed for genomic regions flanking microsatellites. These primers amplify sequence-tagged microsatellite site markers. This approach confirmed the assignment of LG8 to the smallest chromosome H. For the first time, LG5 was linked to the largest chromosome A, LG4 to a medium-sized chromosome E, while LG3 was anchored to the second largest chromosome B. Chromosomes C and D could not be flow-sorted separately and were jointly associated to LG6 and LG7. By the same token, chromosomes F and G were anchored to LG1 and LG2. To establish a set of preferably diagnostic cytogenetic markers, the genomic distribution of various probes was verified using FISH. Moreover, a partial genomic bacterial artificial chromosome (BAC) library was constructed and putative single/low-copy BAC clones were mapped cytogenetically. As a result, two clones were identified localizing specifically to chromosomes E and H, for which no cytogenetic markers were yet available.

Unlocking the barley genome by chromosomal and comparative genomics

We used a novel approach that incorporated chromosome sorting, next-generation sequencing, array hybridization, and systematic exploitation of conserved synteny with model grasses to assign ~86% of the estimated ~32,000 barley (Hordeum vulgare) genes to individual chromosome arms. Using a series of bioinformatically constructed genome zippers that integrate gene indices of rice (Oryza sativa), sorghum (Sorghum bicolor), and Brachypodium distachyon in a conserved synteny model, we were able to assemble 21,766 barley genes in a putative linear order. We show that the barley (H) genome displays a mosaic of structural similarity to hexaploid bread wheat (Triticum aestivum) A, B, and D subgenomes and that orthologous genes in different grasses exhibit signatures of positive selection in different lineages. We present an ordered, information-rich scaffold of the barley genome that provides a valuable and robust framework for the development of novel strategies in cereal breeding.

Assembling a puzzle of dispersed retrotransposable sequences in the genome of chickpea (Cicer arietinum L.)

Several repetitive elements are known to be present in the genome of chickpea (Cicer arietinum L.) including satellite DNA and En/Spm transposons as well as two dispersed, highly repetitive elements, CaRep1 and CaRep2. PCR was used to prove that CaRep1, CaRep2, and previously isolated CaRep3 of C. arietinum represent different segments of a highly repetitive Ty3-gypsy-like retrotransposon (Metaviridae) designated CaRep that makes up large parts of the intercalary heterochromatin. The full sequence of this element including the LTRs and untranslated internal regions was isolated by selective amplification. The restriction pattern of CaRep was different within the annual species of the genus Cicer, suggesting its rearrangement during the evolution of the genus during the last 100 000 years. In addition to CaRep, another LTR and a non-LTR retrotransposon family were isolated, and their restriction patterns and physical localization in the chickpea genome were characterized. The LINE-like element CaLin is only of comparatively low abundance and reveals a considerable heterogeneity. The Ty1-copia-like element (Pseudoviridae) CaTy is located in the distal parts of the intercalary heterochromatin and adjacent euchromatic regions, but it is absent from the centromeric regions. These results together with earlier findings allow to depict the distribution of retroelements on chickpea chromosomes, which extensively resembles the retroelement landscape of the genome of the model legume Medicago truncatula Gaertn.

Isolated chromosomes as a new and efficient source of DArT markers for the saturation of genetic maps

We describe how the diversity arrays technology (DArT) can be coupled with chromosome sorting to increase the density of genetic maps in specific genome regions. Chromosome 3B and the short arm of chromosome 1B (1BS) of wheat were isolated by flow cytometric sorting and used to develop chromosome- and chromosome arm-enriched genotyping arrays containing 2,688 3B clones and 384 1BS clones. Linkage analysis showed that 553 of the 711 polymorphic 3B-derived markers (78%) mapped to chromosome 3B, and 59 of the 68 polymorphic 1BS-derived markers (87%) mapped to chromosome 1BS, confirming the efficiency of the chromosome-sorting approach. To demonstrate the potential for saturation of genetic maps, we constructed a consensus map of chromosome 3B using 19 mapping populations, including some that were genotyped with the 3B-enriched array. The 3B-derived DArT markers doubled the number of genetic loci covered. The resulting consensus map, probably the densest genetic map of 3B available to this date, contains 939 markers (779 DArTs and 160 other markers) that segregate on 304 genetically distinct loci. Importantly, only 2,688 3B-derived clones (probes) had to be screened to obtain almost twice as many polymorphic 3B markers (510) as identified by screening approximately 70,000 whole genome-derived clones (269). Since an enriched DArT array can be developed from less than 5 ng of chromosomal DNA, a quantity which can be obtained within 1 h of sorting, this approach can be readily applied to any crop for which chromosome sorting is available.

Flow cytometry in plant breeding

Since the first report on the flow cytometric study of plant material 35 years ago, analyzing the nuclear DNA content of field bean, an ever increasing number of applications of FCM has been developed and applied in plant science and industry, but a similar length of time elapsed before the appearance of the first complete volume devoted to FCM of plant cells. Most published information on the uses of FCM addresses various aspects of animal (including human) cell biology, thus failing to provide a pertinent substitute. FCM represents an ideal means for the analysis of both cells and subcellular particles, with a potentially large number of parameters analyzed both rapidly, simultaneously, and quantitatively, thereby furnishing statistically exploitable data and allowing for an accurate and facilitated detection of subpopulations. It is, indeed, the summation of these facts that has established FCM as an important, and sometimes essential, tool for the understanding of fundamental mechanisms and processes underlying plant growth, development, and function. In this review, special attention is paid to FCM as applied to plant cells in the context of plant breeding, and some new and less well-known uses of it for plants will be discussed.

Coupling amplified DNA from flow-sorted chromosomes to high-density SNP mapping in barley

Background

Flow cytometry facilitates sorting of single chromosomes and chromosome arms which can be used for targeted genome analysis. However, the recovery of microgram amounts of DNA needed for some assays requires sorting of millions of chromosomes which is laborious and time consuming. Yet, many genomic applications such as development of genetic maps or physical mapping do not require large DNA fragments. In such cases time-consuming de novo sorting can be minimized by utilizing whole-genome amplification.

Results

Here we report a protocol optimized in barley including amplification of DNA from only ten thousand chromosomes, which can be isolated in less than one hour. Flow-sorted chromosomes were treated with proteinase K and amplified using Phi29 multiple displacement amplification (MDA). Overnight amplification in a 20-microlitre reaction produced 3.7 – 5.7 micrograms DNA with a majority of products between 5 and 30 kb. To determine the purity of sorted fractions and potential amplification bias we used quantitative PCR for specific genes on each chromosome. To extend the analysis to a whole genome level we performed an oligonucleotide pool assay (OPA) for interrogation of 1524 loci, of which 1153 loci had known genetic map positions. Analysis of unamplified genomic DNA of barley cv. Akcent using this OPA resulted in 1426 markers with present calls. Comparison with three replicates of amplified genomic DNA revealed >99% concordance. DNA samples from amplified chromosome 1H and a fraction containing chromosomes 2H – 7H were examined. In addition to loci with known map positions, 349 loci with unknown map positions were included. Based on this analysis 40 new loci were mapped to 1H.

Conclusion

The results indicate a significant potential of using this approach for physical mapping. Moreover, the study showed that multiple displacement amplification of flow-sorted chromosomes is highly efficient and representative which considerably expands the potential of chromosome flow sorting in plant genomics.

Coupling amplified DNA from flow-sorted chromosomes to high-density SNP mapping in barley

BACKGROUND

Flow cytometry facilitates sorting of single chromosomes and chromosome arms which can be used for targeted genome analysis. However, the recovery of microgram amounts of DNA needed for some assays requires sorting of millions of chromosomes which is laborious and time consuming. Yet, many genomic applications such as development of genetic maps or physical mapping do not require large DNA fragments. In such cases time-consuming de novo sorting can be minimized by utilizing whole-genome amplification.

RESULTS

Here we report a protocol optimized in barley including amplification of DNA from only ten thousand chromosomes, which can be isolated in less than one hour. Flow-sorted chromosomes were treated with proteinase K and amplified using Phi29 multiple displacement amplification (MDA). Overnight amplification in a 20-microlitre reaction produced 3.7 - 5.7 micrograms DNA with a majority of products between 5 and 30 kb. To determine the purity of sorted fractions and potential amplification bias we used quantitative PCR for specific genes on each chromosome. To extend the analysis to a whole genome level we performed an oligonucleotide pool assay (OPA) for interrogation of 1524 loci, of which 1153 loci had known genetic map positions. Analysis of unamplified genomic DNA of barley cv. Akcent using this OPA resulted in 1426 markers with present calls. Comparison with three replicates of amplified genomic DNA revealed >99% concordance. DNA samples from amplified chromosome 1H and a fraction containing chromosomes 2H - 7H were examined. In addition to loci with known map positions, 349 loci with unknown map positions were included. Based on this analysis 40 new loci were mapped to 1H.

CONCLUSION

The results indicate a significant potential of using this approach for physical mapping. Moreover, the study showed that multiple displacement amplification of flow-sorted chromosomes is highly efficient and representative which considerably expands the potential of chromosome flow sorting in plant genomics.

Functional genomics in chickpea: an emerging frontier for molecular-assisted breeding

Chickpea is a valuable and important agricultural crop, but yield potential is limited by a series of biotic and abiotic stresses, including Ascochyta blight, Fusarium wilt, drought, cold and salinity. To accelerate molecular breeding efforts for the discovery and introgression of stress tolerance genes into cultivated chickpea, functional genomics approaches are rapidly growing. Recently a series of genetic tools for chickpea have become available that have allowed high-powered functional genomics studies to proceed, including a dense genetic map, large insert genome libraries, expressed sequence tag libraries, microarrays, serial analysis of gene expression, transgenics and reverse genetics. This review summarises the development of these genomic tools and the achievements made in initial and emerging functional genomics studies. Much of the initial research focused on Ascochyta blight resistance, and a resistance model has been synthesised based on the results of various studies. Use of the rich comparative genomics resources from the model legumes Medicago truncatula and Lotus japonicus is also discussed. Finally, perspectives on the future directions for chickpea functional genomics, with the goal of developing elite chickpea cultivars, are discussed.

Development of an integrated intraspecific map of chickpea (Cicer arietinum L.) using two recombinant inbred line populations

A composite intraspecific linkage map of chickpea was developed by integrating individual maps developed from two F(8:9) RIL populations with one common parent. Different molecular markers viz. RAPD, ISSR, RGA, SSR and ASAP were analyzed along with three yield related traits: double podding, seeds per pod and seed weight. A total of 273 markers and 186 RILs were used to generate the map with eight linkage groups at a LOD score of >/=3.0 and maximum recombination fraction of 0.4. The map spanned 739.6 cM with 230 markers at an average distance of 3.2 cM between markers. The predominantly used SSR markers facilitated identification of homologous linkage groups from the previously published interspecific linkage map of chickpea and confirmed conservation of the SSR markers across the two maps as well as the variation in terms of marker distance and order. The double podding gene was tagged by the markers NCPGR33 and UBC249z at 2.0 and 1.1 cM, respectively. Whereas, seeds per pod, was tagged by the markers TA2x and UBC465 at 0.1 and 1.8 cM, respectively. Eight QTLs were identified that influence seed weight. The joint map approach allowed mapping a large number of markers with a moderate coverage of the chickpea genome and few linkage gaps.

A linkage map of chickpea (Cicer arietinum L.) based on populations from Kabuli x Desi crosses: location of genes for resistance to fusarium wilt race 0

Two recombinant inbred line (RIL) populations derived from intraspecific crosses with a common parental line (JG62) were employed to develop a chickpea genetic map. Molecular markers, flower colour, double podding, seed coat thickness and resistance to fusarium wilt race 0 (FOC-0) were included in the study. Joint segregation analysis involved a total of 160 markers and 159 RILs. Ten linkage groups (LGs) were obtained that included morphological markers and 134 molecular markers (3 ISSRs, 13 STMSs and 118 RAPDs). Flower colour (B/b) and seed coat thickness (Tt/tt) appeared to be linked to STMS (GAA47). The single-/double-podding locus was located on LG9 jointly with two RAPD markers and STMS TA80. LG3 included a gene for resistance to FOC-0 (Foc0(1)/foc0(1)) flanked by RAPD marker OPJ20(600) and STMS marker TR59. The association of this LG with FOC-0 resistance was confirmed by QTL analysis in the CA2139 x JG62 RIL population where two genes were involved in the resistance reaction. The STMS markers enabled comparison of LGs with preceding maps.

Chromosome sorting in tetraploid wheat and its potential for genome analysis

This study evaluates the potential of flow cytometry for chromosome sorting in durum wheat (Triticum turgidum Desf. var. durum, 2n = 4x = 28). Histograms of fluorescence intensity (flow karyotypes) obtained after the analysis of DAPI-stained chromosomes consisted of three peaks. Of these, one represented chromosome 3B, a small peak corresponded to chromosomes 1A and 6A, and a large peak represented the remaining 11 chromosomes. Chromosomes sorted onto microscope slides were identified after fluorescence in situ hybridization (FISH) with probes for GAA microsatellite, pSc119.2, and Afa repeats. Genomic distribution of these sequences was determined for the first time in durum wheat and a molecular karyotype has been developed for this crop. Flow karyotyping in double-ditelosomic lines of durum wheat revealed that the lines facilitated sorting of any arm of the wheat A- and B-genome chromosomes. Compared to hexaploid wheat, flow karyotype of durum wheat is less complex. This property results in better discrimination of telosomes and high purities in sorted fractions, ranging from 90 to 98%. We have demonstrated that large insert libraries can be created from DNA purified using flow cytometry. This study considerably expands the potential of flow cytogenetics for use in wheat genomics and opens the possibility of sequencing the genome of this important crop one chromosome arm at a time.