Developing repetitive oligo probes from whole genome
Genome sequences of Tifrunner and A. ipaensis were analyzed using the Tandem Repeats Finder (TRF) [40], resulting in 4,595 and 894 repetitive sequences, respectively. The length of these sequences varied between 4 bp and 723 bp, and the copy number between 50 and 29,162. After CD-HIT [41] elimination, a total of 80 and 35 TRs of Tifrunner and A. ipaensis, respectively, were selected for further development of oligos.
A total of 249 oligos were designed using the Oligo 7 [42], and in silico mapped to the reference sequences of peanut with B2DSC [32]. These oligos were first labeled via a random primer labeling method [28] and located on the chromosomes of peanut. 114 oligos produced clear signals in different positions of the chromosomes of Tifrunner (Table S1). Based on their unique patterns and positions after dual-color and sequential FISH, the 114 oligo probes were categorized into 28 types. Oligos with the same position and overlapping signals were classified as the same type. For each of the 28 types, a single oligo was selected and further modified with 6-carboxyfluorescein (FAM) or 6-carboxytetramethylrhodamine (TAMRA) (Table S2). Figure 1 showed the results of oligo Ipa-1463 after in silico mapping and FISH. 1,260 copies were observed in chromosome plots mapped in the region of 63–77 Mbp of chromosome B9 (Fig. 1a). The FISH analysis confirmed that the signals of this oligo were only present in one pair of chromosomes at a similar region. Thus, this chromosome actually corresponded with chromosome B9 in the reference sequence (Fig. 1b).
Development of a genome map-based karyotype of Tifrunner
Based on the unique patterns and sequence composition of 28 oligo probes, a new multiplex #3 oligo probe cocktail was developed with eight oligos, including FAM-modified TIF-439, TIF-185-1, TIF-134-3, and TIF-165-3, and TAMRA-modified Ipa-1162, Ipa-1137, DP-1, and DP-5. Both DP-1 and DP-5 were derived from a previous study by Du et al. [17]. Following the sequential FISH/GISH with multiplex #3, total genomic DNA of A. duranensis and A. ipaensis, and 45S and 5S rDNA assays, a robust karyotype of Tifrunner were established (Fig. 3a–e).
Comparisons of the distributions of the eight oligos in the karyotype and chromosome plots, we found that most signal sites and intensities in the actual chromosomes corresponded well with their positions and copy numbers in the reference sequences (Fig. S1). Finally, a genome map-based karyotype of Tifrunner was established. Each of the actual chromosomes in the karyotype were renumbered as A1~A10 and B1~B10, according to their pseudomolecule number in the genome map of Tifrunner [36]. However, in this karyotype, nine significant non-correspondent signals were observed in seven chromosomes (Fig. 3e–f). For example, the oligos Ipa-1162 and Ipa-1137 evidently had distribution sites on chromosomes A1 and B1 in the chromosome plots, but no signal was observed in the actual chromosomes. In contrast, oligos TIF-439, TIF-185-1, TIF-134-3, and TIF-165-3 produced strong signals in the centromeric regions of chromosomes B2 and B10, but were not in silico mapped in the chromosome plots (Fig. 3d–f, Fig. S1).
To validate the genome map-based karyotype, two chromosome-specific single-copy sequence oligo libraries, L1A-1 and L3A-1, from the upper arm of chromosomes A1 and A3, were used for sequential FISH/GISH analysis, combined with the multiplex #3 and total genomic DNA of A. duranensis and A. ipaensis as probes (Fig. S2). The specific signals of the two chromosome-specific oligo libraries were clearly shown in the two expected chromosomes, indicating considerable correspondence between the actual chromosomes of the karyotype and the genome maps.
Among the A subgenome of this karyotype, chromosomes A1 and A8 (the smallest chromosome) both contained intense green signals in the centromeric regions, while A1 additionally had strong red signals in the terminal region of its short arm. Chromosomes A6, A7, and A10 had 45S or 5S rDNA sites. Chromosomes A2, A3, and A4 had green signals in the terminal regions of the short arms, while A3 and A4 had green signals in the centromeric regions. Chromosome A5 had red signals in the short arm, while A9 had red signals at the subtelomeric region of the long arm. Among the B genome, chromosomes B6, B7, and B8 showed signals with the probes for 45S or 5S rDNA. Chromosome B9 had red signals, and B3 had green signals in the centromeric regions. The other chromosomes had green signals either in the centromeric or telomeric regions with varying intensity (Fig. 3e). Based on the unique patterns observed, all chromosomes could be clearly differentiated in the karyotype.
Chromosome allocation of the oligo probe
To map TRs in the genome map-based karyotype of Tifrunner, 28 representative oligos were analyzed by both FISH and in silico mapping (Fig. 4). The 28 oligos produced more signals in chromosomes of the B genome than in those of the A genome (Fig. 4). Among them, six oligos (TIF-165-3, TIF-439, TIF-556, TIF-198-1, TIF-384-3, and TIF-185-1) produced signals in the interstitial or terminal regions of the chromosomes. Four oligos (TIF-198-2, TIF-416-3, TIF-497, and TIF-342-2) had signals exclusively on the B genome, indicating that these oligos are specific for the B genome. Two oligos (Ipa-1137 and Ipa-1162) had signals only at the secondary constrictions, which fully overlapped with the signals of 45S rDNA, following sequential FISH. Oligo Ipa-1463 exclusively showed signals on one pair of chromosomes, which indicated that it is chromosome-specific. The other 15 oligos had signals in the pericentric regions (Fig. S3).
Among the 28 oligos physically mapped via FISH, the distributions of 22 oligos were same or similar to those in the in silico mapping results. However, six oligos (TIF-89-3, TIF-155-5, TIF-198-1, TIF-359-3, TIF-76-1, and Ipa-1757) showed significant differences between the two maps. For example, TIF-89-3 and TIF-155-5 were in silico mapped in just two pairs of chromosomes with a high number of copies. However, eight pairs of chromosomes evidently showed FISH signals. Similarly, TIF-198-1 was in silico mapped onto one pair of chromosomes alone, but produced signals on 16 pairs of chromosomes. In contrast, obvious sites of TIF-76-1 were mapped to 13 pairs of chromosomes, but produced signals on only five pairs of chromosomes following FISH (Fig. 4, Fig. S4). This may indicated that not all TRs were unambiguously assembled in all chromosomes of Tifrunner.
Identification of chromosomal variations of Silihong (SLH) induced by radiation
To check chromosomal variations in peanut, sequential FISH using multiplex #3 assays, and GISH were conducted on Chinese variety SLH and 70 radiation-induced M1 plants of SLH (Fig. 5). Fourteen M1 plants showed chromosomal variations. For example, in plant 161-1a, one reciprocal translocation was evident, based on unique patterns, which probably occurred between chromosomes 1A and 3B. The segment of one chromosome translocation from the A genome had extensive green signals covering the arm and red signals at the terminal. This pattern was similar to that of the signals on the centromere and the lower arm (L) of 1A, and the signal patterns of another segment from the B genome was like the upper arm (S) of 3B with the centromere, which indicated that they were T 3BL-1AS·1AL and T 1AS·3BL-3BS. Further FISH was performed using the single-copy oligo library probes L1A-1 and L3A-1 which are exclusively hybridize to the upper arms of 1A and 3A. FISH signals confirmed the identification of the two translocated chromosomes (Fig. 5f–h).
A total of 13 other plants were identified with 17 translocations, one deletion, and eight monosomic chromosomes (Fig. S5). Among these chromosomal variations, eight translocations were observed between homoeologous chromosomes, and nine translocations were observed between non-homologous chromosomes. Chromosomes 1, 3, and 5 showed a greater number of translocations (Fig. 6).