Chromatin characterization by banding techniques, in situ hybridization, and nuclear DNA content in Cicer L. (Leguminosae)

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.

Brazilian Kayabi Indian accessions of peanut, Arachis hypogaea (Fabales, Fabaceae): origin, diversity and evolution

Peanut is a crop of the Kayabi tribe, inhabiting the Xingu Indigenous Park, Brazil. Morphological analysis of Xingu accessions showed variation exceeding that described for cultivated peanuts. This raised questions as to the origin of the Xingu accessions: are they derived from different species, or is their diversity a result of different evolutionary and selection processes? To answer these questions, cytogenetic and genotyping analyses were conducted. The karyotypes of Xingu accessions analyzed are very similar to each other, to an A. hypogaea subsp. fastigiata accession and to the wild allotetraploid A. monticola. The accessions share the number and general morphology of the chromosomes; DAPI+ bands; 5S and 45S rDNA loci distribution and a high genomic affinity with A. duranensis and A. ipaënsis genomic probes. However, the number of CMA3+ bands differs from those determined for A. hypogaea and A. monticola, which are also different from each other. SNP genotyping grouped all Arachis allotetraploids into four taxonomic groups: Xingu accessions were closer to A. monticola and A. hypogaea subsp. hypogaea. Our data suggests that the morphological diversity within these accessions is not associated with a different origin and can be attributed to morphological plasticity and different selection by the Indian tribes.

Chromosome diversity in species of the genus Arachis, revealed by FISH and CMA/DAPI banding, and inferences about their karyotype differentiation

The species of the genus Arachis (Leguminosae) are ordered into nine sections. The assignment of genome types in this genus has been based on cross-compatibility analysis and molecular cytogenetic studies. The latter has also allowed karyotypically establishing well-defined genomes and reassigning the genome of several species. However, most of these studies have been focused mainly on the sections Arachis and Rhizomatosae. To increase the knowledge about the chromosome diversity of the whole genus, here we performed a detailed karyotype characterization of representative species of most of the sections and genomes of Arachis. This characterization included chromosome morphology, CMA/DAPI chromosome banding, and chromosome marker localization (rDNAloci and one satDNA sequence) by fluorescent in situ hybridization (FISH). Based on the data obtained and other previously published data, we established the karyotype similarities by cluster analysis and defined eleven karyotype groups. The grouping was partly coincident with the traditional genome assignment, except for some groups and some individual species. Karyotype similarities among some genomes were also found. The main characteristics of each karyotype group of Arachis were summarized. Together, our results provide information that may be beneficial for future cytogenetic and evolutionary studies, and also contribute to the identification of interspecific hybrids.

Karyotype diversity and 2C DNA content in species of the Caesalpinia group

BACKGROUND

The Leguminosae family is the third-largest family of angiosperms, and Caesalpinioideae is its second-largest subfamily. A great number of species (approximately 205) are found in the Caesalpinia group within this subfamily; together with these species' phenotypic plasticity and the similarities in their morphological descriptors, make this a complex group for taxonomic and phylogenetic studies. The objective of the present work was to evaluate the karyotypic diversity and the 2C DNA content variation in 10 species of the Caesalpinia group, representing six genera: Paubrasilia, Caesalpinia, Cenostigma, Poincianella, Erythrostemon and Libidibia. The GC-rich heterochromatin and 45S rDNA sites (which are used as chromosome markers) were located to evaluate the karyotype diversity in the clade. The variation in the 2C DNA content was determined through flow cytometry.

RESULTS

The fluorochrome banding indicated that the chromomycin A₃⁺/4',6-diamidino-2-phenylindole^- blocks were exclusively in the terminal regions of the chromosomes, coinciding with 45S rDNA sites in all analyzed species. Physical mapping of the species (through fluorescence in situ hybridization) revealed variation in the size of the hybridization signals and in the number and distribution of the 45S rDNA sites. All hybridization sites were in the terminal regions of the chromosomes. In addition, all species had a hybridization site in the fourth chromosome pair. The 2C DNA content ranged from 1.54 pg in Erythrostemon calycina to 2.82 pg in the Paubrasilia echinata large-leaf variant. The Pa. echinata small-leaf variant was isolated from the other leaf variants through Scoot-Knott clustering.

CONCLUSIONS

The chromosome diversity and the variation in the 2C DNA content reinforce that the actual taxonomy and clustering of the analyzed taxa requires more genera that were previously proposed. This fact indicates that taxonomy, phylogeny and cytoevolutionary inference related to the complex Caesalpinia group have to be done through integrative evaluation.

Heterochromatin distribution and comparative karyo-morphological studies in Vignaumbellata Thunberg, 1969 and V.aconitifolia Jacquin, 1969 (Fabaceae) accessions

Chromosome studies along with heterochromatin distribution pattern analysis have been carried out in two domesticated species of Vigna Savi, 1824 which grow in contrasting geo-climatic conditions of India: Vignaumbellata Thunberg, 1969, a legume well acclimatized to subtropical hilly regions of North-east India and Vignaaconitifolia Jacquin, 1969, a species of arid and semi-arid regions in desert plains of Western India. Karyo-morphological studies in both species reveal 2n = 22 chromosomes without any evidence of numerical variation and the overall karyotype symmetry in chromosome morphology suggest that the diversification at intraspecific level in genus Vigna has occurred through structural alteration of chromosomes, rather than numerical changes. Heterochromatin distribution as revealed by fluorochrome binding pattern using CMA3 and DAPI, confirms the occurrence of relatively more GC content in Vignaaconitifolia as compared to Vignaumbellata. However, AT content was found to be comparatively higher in Vignaumbellata which perhaps play a role in species interrelationships.

Heterochromatin characterization through differential fluorophore binding pattern in some species of Vigna Savi

Heterochromatin regions are the most intensively studied and best known chromosome markers in plants. In Vigna species, blocks of constitutive heterochromatin were found either in the terminal or interstitial region of the chromosomes. The number and distribution of CMA(+) and DAPI(+) binding sites exhibit high chromosomal variability with characteristic unique banding patterns in all the eight taxa. A predominant feature was observed, i.e., most of the CMA(+) binding sites were in the terminal region of the short arm of some chromosomes while DAPI(+) binding sites were found mostly in the intercalary region of the chromosomes. The higher divergence in the heterochromatin blocks, as revealed by chromomycin A3 (CMA) binding pattern, in a few taxa, viz. Vigna glabrescens, Vigna khandalensis, and Vigna mungo, suggests that the processes of divergent evolution of repetitive sequences in genomic DNA involve a guanine-cytosine (GC)-rich region. On the contrary, Vigna dalzelliana had shown a prominent adenine-thymine (AT)-rich repetitive DNA sequence in terminal regions in the short arm of chromosomes while Vigna umbellata had shown in interstitial regions. The presence of prominent heterochromatic-rich regions, either GC- or AT-rich regions, does facilitate the rate of chromosomal rearrangements leading to restructuring of the karyotypes and thereby helping the species to attempt structural alterations as means of speciation.

Comparative analysis of kabuli chickpea transcriptome with desi and wild chickpea provides a rich resource for development of functional markers

Chickpea (Cicer arietinum L.) is an important crop legume plant with high nutritional value. The transcriptomes of desi and wild chickpea have already been sequenced. In this study, we sequenced the transcriptome of kabuli chickpea, C. arietinum (genotype ICCV2), having higher commercial value, using GS-FLX Roche 454 and Illumina technologies. The assemblies of both Roche 454 and Illumina datasets were optimized using various assembly programs and parameters. The final optimized hybrid assembly generated 43,389 transcripts with an average length of 1065 bp and N50 length of 1653 bp representing 46.2 Mb of kabuli chickpea transcriptome. We identified a total of 5409 simple sequence repeats (SSRs) in these transcript sequences. Among these, at least 130 and 493 SSRs were polymorphic with desi (ICC4958) and wild (PI489777) chickpea, respectively. In addition, a total of 1986 and 37,954 single nucleotide polymorphisms (SNPs) were predicted in kabuli/desi and kabuli/wild genotypes, respectively. The SNP frequency was 0.043 SNP per kb for kabuli/desi and 0.821 SNP per kb for kabuli/wild, reflecting very low genetic diversity in chickpea. Further, SSRs and SNPs present in tissue-specific and transcription factor encoding transcripts have been identified. The experimental validation of a selected set of polymorphic SSRs and SNPs exhibited high intra-specific polymorphism potential between desi and kabuli chickpea, suggesting their utility in large-scale genotyping applications. The kabuli chickpea gene index assembled, and SSRs and SNPs identified in this study will serve as useful genomic resource for genetic improvement of chickpea.

Transcriptome sequencing of wild chickpea as a rich resource for marker development

The transcriptome of cultivated chickpea (Cicer arietinum L.), an important crop legume, has recently been sequenced. Here, we report sequencing of the transcriptome of wild chickpea, C. reticulatum (PI489777), the progenitor of cultivated chickpea, by GS-FLX 454 technology. The optimized assembly of C. reticulatum transcriptome generated 37 265 transcripts in total with an average length of 946 bp. A total of 4072 simple sequence repeats (SSRs) could be identified in these transcript sequences, of which at least 561 SSRs were polymorphic between C. arietinum and C. reticulatum. In addition, a total of 36 446 single-nucleotide polymorphisms (SNPs) were identified after optimization of probability score, quality score, read depth and consensus base ratio. Several of these SSRs and SNPs could be associated with tissue-specific and transcription factor encoding transcripts. A high proportion (92-94%) of polymorphic SSRs and SNPs identified between the two chickpea species were validated successfully. Further, the estimation of synonymous substitution rates of orthologous transcript pairs suggested that the speciation event for divergence of C. arietinum and C. reticulatum may have happened approximately 0.53 million years ago. The results of our study provide a rich resource for exploiting genetic variations in chickpea for breeding programmes.

EST-derived genic molecular markers: development and utilization for generating an advanced transcript map of chickpea

Well-saturated linkage maps especially those based on expressed sequence tag (EST)-derived genic molecular markers (GMMs) are a pre-requisite for molecular breeding. This is especially true in important legumes such as chickpea where few simple sequence repeats (SSR) and even fewer GMM-based maps have been developed. Therefore, in this study, 2,496 ESTs were generated from chickpea seeds and utilized for the development of 487 novel EST-derived functional markers which included 125 EST-SSRs, 151 intron targeted primers (ITPs), 109 expressed sequence tag polymorphisms (ESTPs), and 102 single nucleotide polymorphisms (SNPs). Whereas ESTSSRs, ITPs, and ESTPs were developed by in silico analysis of the developed EST sequences, SNPs were identified by allele resequencing and their genotyping was performedusing the Illumina GoldenGate Assay. Parental polymorphism was analyzed between C. arietinum ICC4958 and C. reticulatum PI489777, parents of the reference chickpea mapping population, using a total of 872 markers: 487 new gene-based markers developed in this study along with 385 previously published markers, of which 318 (36.5%) were found to be polymorphic and were used for genotyping. The genotypic data were integrated with the previously published data of 108 markers and an advanced linkage map was generated that contained 406 loci distributed on eight linkage groups that spanned 1,497.7 cM. The average marker density was 3.68 cM and the average number of markers per LG was 50.8. Among the mapped markers, 303 new genomic locations were defined that included 177 gene-based and 126 gSSRs (genomic SSRs) thereby producing the most advanced gene-rich map of chickpea solely based on co-dominant markers.

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.

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.

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

Procedures for flow cytometric analysis and sorting of mitotic chromosomes (flow cytogenetics) have been developed for chickpea (Cicer arietinum). Suspensions of intact chromosomes were prepared from root tips treated to achieve a high degree of metaphase synchrony. The optimal protocol consisted of a treatment of roots with 2 mmol/L hydroxyurea for 18 h, a 4.5-h recovery in hydroxyurea-free medium, 2 h incubation with 10 micromol/L oryzalin, and ice-water treatment overnight. This procedure resulted in an average metaphase index of 47%. Synchronized root tips were fixed in 2% formaldehyde for 20 min, and chromosome suspensions prepared by mechanical homogenization of fixed root tips. More than 4 x 10(5) morphologically intact chromosomes could be isolated from 15 root tips. Flow cytometric analysis of DAPI-stained chromosomes resulted in histograms of relative fluorescence intensity (flow karyotypes) containing eight peaks, representing individual chromosomes and/or groups of chromosomes with a similar relative DNA content. Five peaks could be assigned to individual chromosomes (A, B, C, G, H). The parity of sorted chromosome fractions was high, and chromosomes B and H could be sorted with 100% purity. PCR on flow-sorted chromosome fractions with primers for sequence-tagged microsatellite site (STMS) markers permitted assignment of the genetic linkage group LG8 to the smallest chickpea chromosome H. This study extends the number of legume species for which flow cytogenetics is available, and demonstrates the potential of flow cytogenetics for genome mapping in chickpea.

An intraspecific linkage map of the chickpea ( Cicer arietinum L.) genome based on sequence tagged microsatellite site and resistance gene analog markers

An intraspecific linkage map of the chickpea genome based on STMS as anchor markers, was established using an F(2) population of chickpea cultivars with contrasting disease reactions to Ascochyta rabiei (Pass.) Lab. At a LOD-score of 2.0 and a maximum recombination distance of 20 cM, 51 out of 54 chickpea-STMS markers (94.4%), three ISSR markers (100%) and 12 RGA markers (57.1%) were mapped into eight linkage groups. The chickpea-derived STMS markers were distributed throughout the genome, while the RGA markers clustered with the ISSR markers on linkage groups LG I, II and III. The intraspecific linkage map spanned 534.5 cM with an average interval of 8.1 cM between markers. Sixteen markers (19.5%) were unlinked, while l1 chickpea-STMS markers (20.4%) deviated significantly ( P < 0.05) from the expected Mendelian segregation ratio and segregated in favor of the maternal alleles. However, ten of the distorted chickpea-STMS markers were mapped and clustered mostly on LG VII, suggesting the association of these loci in the preferential transmission of the maternal germ line. Preliminary comparative mapping revealed that chickpea may have evolved from Cicer reticulatum, possibly via inversion of DNA sequences and minor chromosomal translocation. At least three linkage groups that spanned a total of approximately 79.2 cM were conserved in the speciation process.

MicroMeasure: A new computer program for the collection and analysis of cytogenetic data

The ability to identify individual chromosomes in cytological preparations is an essential component of many investigations. While several computer software applications have been used to facilitate such quantitative karyotype analysis, most of these programs are limited by design for specific types of analyses, or can be used only with specific hardware configurations. MicroMeasure is a new image analysis application that may be used to collect data for a wide variety of chromosomal parameters from electronically captured or scanned images. Unlike similar applications, MicroMeasure may be individually configured by the end user to suit a wide variety of research needs. This program can be used with most common personal computers, and requires no unusual or specific hardware. MicroMeasure is made available to the research community without cost by the Department of Biology at Colorado State University via the World Wide Web at http://www.biology.colostate.edu/MicroMeasure.Key words: MicroMeasure, computer program, chromosome measurement, cytogenetics.

Chromosomal localization and distribution of simple sequence repeats and the Arabidopsis-type telomere sequence in the genome of Cicer arietinum L

We used fluorescence in situ hybridization to probe the physical organization of five simple sequence repeat motifs and the Arabidopsis-type telomeric repeat in metaphase chromosomes and interphase nuclei of chickpea (Cicer arietinum L.). Hybridization signals were observed with the whole set of probes and on all chromosomes, but the distribution and intensity of signals varied depending on the motif. On root-tip metaphase chromosomes, CA and GATA repeats were mainly restricted to centromeric areas, with additional GATA signals along some chromosomes. TA, A and AAC repeats were organized in a more dispersed manner, with centromeric regions being largely excluded. In interphase nuclei of the inner integument, CA and GATA signals predominantly occurred in the heterochromatic endochromocentres, whereas the other motifs were found both in eu- and heterochromatin. The distribution of the Arabidopsis-type telomeric repeat (TTTAGGG)n on metaphase chromosomes was found to be quite exceptional. One major cluster of repeats was spread along the short arm of chromosome B, whereas a second, weaker signal occurred interstitially on chromosome A. Only faint and inconsistent hybridization signals were visualized with the same probe at the chromosomal termini.