Molecular cytogenetic analysis of genome-specific repetitive elements in Citrus clementina Hort. Ex Tan. and its taxonomic implications

The genome of Citrus australasica reveals disease resistance and other species specific genes

BACKGROUND

The finger lime (Citrus australasica), one of six Australian endemic citrus species shows a high natural phenotypic diversity and novel characteristics. The wide variation and unique horticultural features have made this lime an attractive candidate for domestication. Currently no haplotype resolved genome is available for this species. Here we present a high quality, haplotype-resolved reference genome for this species using PacBio HiFi and Hi-C sequencing.

RESULTS

Hifiasm assembly and SALSA scaffolding resulted in a collapsed genome size of 344.2 Mb and 321.1 Mb and 323.2 Mb size for the two haplotypes. The nine pseudochromosomes of the collapsed genome had an N50 of 35.2 Mb, 99.1% genome assembly completeness and 98.9% gene annotation completeness (BUSCO). A total of 41,304 genes were predicted in the nuclear genome. Comparison with C. australis revealed that 13,661 genes in pseudochromosomes were unique in C. australasica. These were mainly involved in plant-pathogen interactions, stress response, cellular metabolic and developmental processes, and signal transduction. The two genomes showed a syntenic arrangement at the chromosome level with large structural rearrangements in some chromosomes. Genetic variation among five C. australasica cultivars was analysed. Genes related to defense, synthesis of volatile compounds and red/yellow coloration were identified in the genome. A major expansion of genes encoding thylakoid curvature proteins was found in the C. australasica genome.

CONCLUSIONS

The genome of C. australasica present in this study is of high quality and contiguity. This genome helps deepen our understanding of citrus evolution and reveals disease resistance and quality related genes with potential to accelerate the genetic improvement of citrus.

Five Species of Taxus Karyotype Based on Oligo-FISH for 5S rDNA and (AG3T3)3

As a relict plant, Taxus is used in a variety of medicinal ingredients, for instance to treat a variety of cancers. Taxus plants are difficult to distinguish from one another due to their similar morphology; indeed, some species of Taxus cytogenetic data still are unclear. Oligo-FISH can rapidly and efficiently provide insight into the genetic composition and karyotype. This is important for understanding the organization and evolution of chromosomes in Taxus species. We analysed five Taxus species using two oligonucleotide probes. (AG₃T₃)₃ signals were distributed at the chromosome ends and the centromere of five species of Taxus. The 5S rDNA signal was displayed on two chromosomes of five species of Taxus. In addition to Taxus wallichiana var. mairei, 5S rDNA signals were found proximal in the remaining four species, which signals a difference in its location. The karyotype formula of Taxus wallichiana was 2n = 2x = 24m, its karyotype asymmetry index was 55.56%, and its arm ratio was 3.0087. Taxus × media's karyotype formula was 2n = 2x = 24m, its karyotype asymmetry index was 55.09%, and its arm ratio was 3.4198. The karyotype formula of Taxus yunnanensis was 2n = 2x = 24m, its karyotype asymmetry index was 55.56%, and its arm ratio was 2.6402. The karyotype formula of Taxus cuspidate was 2n = 2x = 24m, its karyotype asymmetry index was 54.67%, its arm ratio was 3.0135, and two chromosomes exhibited the 5S rDNA signal. The karyotype formula of T. wallichiana var. mairei was 2n= 2x = 22m + 2sm, its karyotype asymmetry index was 54.33%, and its arm ratio was 2.8716. Our results provide the karyotype analysis and physical genetic map of five species of Taxus, which contributes to providing molecular cytogenetics data for Taxus.

Accurate Chromosome Identification in the Prunus Subgenus Cerasus (Prunus pseudocerasus) and its Relatives by Oligo-FISH

A precise, rapid and straightforward approach to chromosome identification is fundamental for cytogenetics studies. However, the identification of individual chromosomes was not previously possible for Chinese cherry or other Prunus species due to the small size and similar morphology of their chromosomes. To address this issue, we designed a pool of oligonucleotides distributed across specific pseudochromosome regions of Chinese cherry. This oligonucleotide pool was amplified through multiplex PCR with specific internal primers to produce probes that could recognize specific chromosomes. External primers modified with red and green fluorescence tags could produce unique signal barcoding patterns to identify each chromosome concomitantly. The same oligonucleotide pool could also discriminate all chromosomes in other Prunus species. Additionally, the 5S/45S rDNA probes and the oligo pool were applied in two sequential rounds of fluorescence in situ hybridization (FISH) localized to chromosomes and showed different distribution patterns among Prunus species. At the same time, comparative karyotype analysis revealed high conservation among P. pseudocerasus, P. avium, and P. persica. Together, these findings establish this oligonucleotide pool as the most effective tool for chromosome identification and the analysis of genome organization and evolution in the genus Prunus.

Establishment and Optimization of Molecular Cytogenetic Techniques (45S rDNA-FISH, GISH, and Fiber-FISH) in Kiwifruit (Actinidia Lindl.)

The kiwifruit (Actinidia chinensis) has long been regarded as "the king of fruits" for its nutritional importance. However, the molecular cytogenetics of kiwifruit has long been hampered because of the large number of basic chromosome (x = 29), the inherent small size and highly similar morphology of metaphase chromosomes. Fluorescence in situ hybridization (FISH) is an indispensable molecular cytogenetic technique widely used in many plant species. Herein, the effects of post-hybridization washing temperature on FISH, blocking DNA concentration on genomic in situ hybridization (GISH), extraction method on nuclei isolation and the incubation time on the DNA fiber quality in kiwifruit were evaluated. The post-hybridization washing in 2 × saline sodium citrate (SSC) solution for 3 × 5 min at 37^°C ensured high stringency and distinct specific FISH signals in kiwifruit somatic chromosomes. The use of 50 × blocking DNA provided an efficient and reliable means of discriminating between chromosomes derived from in the hybrids of A. chinensis var. chinensis (2n = 2x = 58) × A. eriantha (2n = 2x = 58), and inferring the participation of parental genitors. The chopping method established in the present study were found to be very suitable for preparation of leaf nuclei in kiwifruit. A high-quality linear DNA fiber was achieved by an incubation of 20 min. The physical size of 45S rDNA signals was approximately 0.35-0.40 μm revealed by the highly reproducible fiber-FISH procedures established and optimized in this study. The molecular cytogenetic techniques (45S rDNA-FISH, GISH, and high-resolution fiber-FISH) for kiwifruit was for the first time established and optimized in the present study, which is the foundation for the future genomic and evolutionary studies and provides chromosomal characterization for kiwifruit breeding programs.

Oligo-FISH Can Identify Chromosomes and Distinguish Hippophaë rhamnoides L. Taxa

Oligo-fluorescence in situ hybridization (FISH) facilitates precise chromosome identification and comparative cytogenetic analysis. Detection of autosomal chromosomes of Hippophaë rhamnoides has not been achieved using oligonucleotide sequences. Here, the chromosomes of five H. rhamnoides taxa in the mitotic metaphase and mitotic metaphase to anaphase were detected using the oligo-FISH probes (AG3T3)3, 5S rDNA, and (TTG)6. In total, 24 small chromosomes were clearly observed in the mitotic metaphase (0.89–3.03 μm), whereas 24–48 small chromosomes were observed in the mitotic metaphase to anaphase (0.94–3.10 μm). The signal number and intensity of (AG3T3)3, 5S rDNA, and (TTG)6 in the mitotic metaphase to anaphase chromosomes were nearly consistent with those in the mitotic metaphase chromosomes when the two split chromosomes were integrated as one unit. Of note, 14 chromosomes (there is a high chance that sex chromosomes are included) were exclusively identified by (AG3T3)3, 5S rDNA, and (TTG)6. The other 10 also showed a terminal signal with (AG3T3)3. Moreover, these oligo-probes were able to distinguish one wild H. rhamnoides taxon from four H. rhamnoides taxa. These chromosome identification and taxa differentiation data will help in elucidating visual and elaborate physical mapping and guide breeders’ utilization of wild resources of H. rhamnoides.

Interstitial Telomeric-like Repeats (ITR) in Seed Plants as Assessed by Molecular Cytogenetic Techniques: A Review

The discovery of telomeric repeats in interstitial regions of plant chromosomes (ITRs) through molecular cytogenetic techniques was achieved several decades ago. However, the information is scattered and has not been critically evaluated from an evolutionary perspective. Based on the analysis of currently available data, it is shown that ITRs are widespread in major evolutionary lineages sampled. However, their presence has been detected in only 45.6% of the analysed families, 26.7% of the sampled genera, and in 23.8% of the studied species. The number of ITR sites greatly varies among congeneric species and higher taxonomic units, and range from one to 72 signals. ITR signals mostly occurs as homozygous loci in most species, however, odd numbers of ITR sites reflecting a hemizygous state have been reported in both gymnosperm and angiosperm groups. Overall, the presence of ITRs appears to be poor predictors of phylogenetic and taxonomic relatedness at most hierarchical levels. The presence of ITRs and the number of sites are not significantly associated to the number of chromosomes. The longitudinal distribution of ITR sites along the chromosome arms indicates that more than half of the ITR presences are between proximal and terminal locations (49.5%), followed by proximal (29.0%) and centromeric (21.5%) arm regions. Intraspecific variation concerning ITR site number, chromosomal locations, and the differential presence on homologous chromosome pairs has been reported in unrelated groups, even at the population level. This hypervariability and dynamism may have likely been overlooked in many lineages due to the very low sample sizes often used in cytogenetic studies.

Distribution of FISH oligo-5S rDNA and oligo-(AGGGTTT)3 in Hibiscus mutabilis L

Hibiscus exhibits high variation in chromosome number both within and among species. The Hibiscus mutabilis L. karyotype was analyzed in detail using fluorescence in situ hybridization (FISH) with oligonucleotide probes for (AG₃T₃)₃ and 5S rDNA, which were tested here for the first time. In total, 90 chromosomes were counted in prometaphase and metaphase, and all exhibited similarly intense (AG₃T₃)₃ signals at both ends. (AG₃T₃)₃ showed little variation and thus did not allow discrimination among H. mutabilis chromosomes, but its location at both ends confirmed the integrity of each chromosome, thus contributing to accurate counting of the numerous, small chromosomes. Oligo-5S rDNA marked the proximal/distal regions of six chromosomes: weak signals on chromosomes 7 and 8, slightly stronger signals on chromosomes 15 and 16, and very strong signals on chromosomes 17 and 18. Therefore, 5S rDNA could assist in chromosome identification in H. mutabilis. Metaphase chromosome lengths ranged from 3.00 to 1.18 μm, indicating small chromosomes. The ratios of longest to shortest chromosome length in prometaphase and metaphase were 2.58 and 2.54, respectively, indicating karyotype asymmetry in H. mutabilis. These results provide an exact chromosome number and a physical map, which will be useful for genome assembly and contribute to molecular cytogenetics in the genus Hibiscus.

Localization and characterization of Citrus centromeres by combining half-tetrad analysis and CenH3-associated sequence profiling

The physical locations of citrus centromere are revealed by combining genetic and immunological assays for the first time and nine citrus centromere-specific markers for cytogenetics are mined. Centromere localization is challenging, because highly redundant repetitive sequences in centromeric regions make sequence assembly difficult. Although several citrus genomes have been released, the centromeric regions and their characteristics remain to be elucidated. Here, we mapped citrus centromeres through half-tetrad analysis (HTA) that included the genotyping of 54 tetraploid hybrids derived from 2n megagametophytes of Nadorcott tangor with 212 single nucleotide polymorphism (SNP) markers. The sizes of centromeric regions, which estimated based on the heterozygosity restitution rate pattern along the chromosomes, ranged from 1.12 to 18.19 Mb. We also profiled the binding sequences with the centromere-specific histone variant CenH3 by chromatin immunoprecipitation sequencing (ChIP-seq). Based on the positions of the top ten CenH3-enriched contigs, the sizes of centromeric regions were estimated to range from 0.01 to 7.60 Mb and were either adjacent to or included in the centromeric regions identified by HTA. We used DNA probes from two repeats selected from the centromeric regions and seven CenH3-binding centromeric repeats to verify centromeric locations by fluorescence in situ hybridization (FISH). Centromere localization in citrus will contribute to the mining of centromeric/pericentromeric markers, thus to facilitate the rapid identification of mechanisms underlying 2n gamete formation and serve the polyploidy breeding.

Extraordinarily conserved chromosomal synteny of Citrus species revealed by chromosome-specific painting

Reliable identification of individual chromosomes in eukaryotic species is the foundation for comparative chromosome synteny and evolutionary studies. Unfortunately, chromosome identification has been a major challenge for plants with small chromosomes, such as the Citrus species. We developed oligonucleotide-based chromosome painting probes for all nine chromosomes in Citrus maxima (Pummelo). We were able to identify all C. maxima chromosomes in the same metaphase cells using multiple rounds of sequential fluorescence in situ hybridization with the painting probes. We conducted comparative chromosome painting analysis in six different Citrus and related species. We found that each painting probe hybridized to only a single chromosome in all other five species, suggesting that the six species have maintained a complete chromosomal synteny after more than 9 million years of divergence. No interchromosomal rearrangement was identified in any species. These results support the hypothesis that karyotypes of woody species are more stable than herbaceous plants because woody plants need a longer period to fix chromosome structural variants in natural populations.

Integration of mandarin (Citrus reticulata) cytogenetic map with its genome sequence

Citrus is an extremely important genus in terms of world fruit production. Despite its economic importance and the small genome sizes of its species (2n = 18, 1C = 430 ± 68 Mbp), entire genomic assemblies have only recently become available for some of its representatives. Together with the previous CMA/DAPI banding and fluorescence in situ hybridization (FISH) in the group, these data are important for understanding the complex relationships between its species and for assisting breeding programs. To anchor genomic data with the cytogenetic map of mandarin (Citrus reticulata), the parental species of several economically important hybrids such as sweet orange and clementine, 18 BAC (bacterial artificial chromosome) clones were used. Eleven clementine BACs were positioned by BAC-FISH, doubling the number of chromosome markers so far available for BAC-FISH in citrus. Additionally, six previously mapped BACs were end-sequenced, allowing, together with one BAC previously sequenced, their assignment to scaffolds and the subsequent integration of chromosomes and the genome assembly. This study therefore established correlations between mandarin scaffolds and chromosomes, allowing further structural genomic and comparative study with the sweet orange genome, as well as insights into the chromosomal evolution of the group.

Distinguishing Sichuan Walnut Cultivars and Examining Their Relationships with Juglans regia and J. sigillata by FISH, Early-Fruiting Gene Analysis, and SSR Analysis

Walnuts are economically important tree species in Sichuan Province (China) that provide heathy nuts. Fluorescence in situ hybridization (FISH) and analyses of an early-fruiting gene fragment and simple sequence repeats (SSRs) were used to distinguish Sichuan walnut cultivars and examine their relationships with Juglans regia L. and Juglans sigillata Dode. Thirty-four small chromosomes were counted in four Sichuan walnut cultivars. In the four cultivars, 5S rDNA was located in the proximal regions of two chromosomes (5 and 6), while (AG₃T₃)₃ was located at both ends of each chromosome. The existence of the signal at both chromosome ends ensured accurate chromosome counts. 5S rDNA and (AG₃T₃)₃ were not effective in identifying Sichuan walnut cultivars. Evolutionary analysis involving 32 early-fruiting nucleotide sequences from Sichuan walnut materials were performed with the maximum likelihood method. There were a total of 602 positions. All positions with gaps and missing data were eliminated, resulting in a final dataset of 562 positions. The ML tree with the highest log likelihood (-1607.82) revealed two obvious groups: one including materials of J. regia, which fruits 1 year after grafting, and another including materials of J. sigillata, which fruits >3 years after grafting. The early-fruiting gene fragment divided 22 walnut materials (10 walnut cultivars and 12 walnut accessions) into two groups, indicating that it was somewhat effective for distinguishing Sichuan walnut cultivars. Furthermore, 22 SSR loci were revealed to identify nine walnut cultivars. Eight cultivars were exclusively discerned by one SSR locus each: Chuanzao 1 [CUJRB307 (116) or CUJRA206a (182)], Chuanzao 2 [JSI-73 (154)], Shuangzao [CUJRB103a (123), CUJRB218 (144), JSI-71 (146), or CUJRA206a (176)], Shimianju [ZMZ11 (138)], Meigupao [CUJRB218 (149), CUJRB103a (151), or CUJRA206a (190)], Muzhilinhe [CUJRB220 (136), ZMZ11 (147), CUJRC310 (156), or JSI-73 (166)], Maerkang [CUJRA124 (154), CUJRB218 (159), or CUJRA123 (182)], Yanyuanzao [CUJRA124 (150) or CUJRA206a (192)]. The Shuling cultivar was identified by the combination of ZMZ11 (148) and other SSR loci, which distinguished and excluded the Chuanzao 1 and Yanyuanzao cultivars. Our results will guide the identification and breeding of Sichuan walnut cultivars.

Integrated Karyotypes of Diploid and Tetraploid Carrizo Citrange (Citrus sinensis L. Osbeck × Poncirus trifoliata L. Raf.) as Determined by Sequential Multicolor Fluorescence in situ Hybridization With Tandemly Repeated DNA Sequences

Carrizo citrange [Citrus sinensis (L.) Osbeck × Poncirus trifoliata (L.) Raf., CC] is one of the most widely used rootstocks in citriculture worldwide, but its cytogenetic study has been hampered by its inherent small size, morphological similarity to mitotic chromosomes, and lack of accessible cytological landmarks. In our previous study, a spontaneously occurring tetraploid CC seedling was discovered. The main goals of this study were to elucidate the chromosome constitution and construct the karyotypes of diploid CC rootstock and its corresponding spontaneously occurring tetraploid. To accomplish these, the chromosomal characteristics were investigated by sequential multicolor fluorescence in situ hybridization (FISH) with eight properly labeled repetitive DNA sequences, including a centromere-like repeat, four satellite repeats, two rDNAs, and an oligonucleotide of telomeric (TTTAGGG) _n repeat. The results nicely demonstrated that these repetitive DNAs are reliable cytogenetic markers that collectively facilitate simultaneous and unequivocal identification of homologous chromosome pairs. Based on chromosome size and morphology together with FISH patterns of repetitive DNAs, an integrated karyotype of CC rootstock was constructed, consisting of 2n = 2x = 12m (1sat) + 6sm with karyotype asymmetry degree being divided into 2B category. Cytogenetically speaking, the variable and asymmetric distribution patterns of these repetitive DNAs were fully confirmed the hybrid nature of CC rootstock. In addition, comparative distribution patterns and chromosomal localizations of these repetitive DNAs convincingly showed that this tetraploid CC material arose from somatic chromosome doubling of diploid CC rootstock. This study revealed, for the first time, the integrated karyotype and chromosomal characteristics of this important citrus rootstock as well as its spontaneously occurring tetraploid plant. Furthermore, this study is a good prospective model for study species with morphologically indistinguishable small chromosomes.

Physical Map of FISH 5S rDNA and (AG3T3)3 Signals Displays Chimonanthus campanulatus R.H. Chang & C.S. Ding Chromosomes, Reproduces its Metaphase Dynamics and Distinguishes Its Chromosomes

Chimonanthus campanulatus R.H. Chang & C.S. Ding is a good horticultural tree because of its beautiful yellow flowers and evergreen leaves. In this study, fluorescence in situ hybridization (FISH) was used to analyse mitotic metaphase chromosomes of Ch. campanulatus with 5S rDNA and (AG₃T₃)₃ oligonucleotides. Twenty-two small chromosomes were observed. Weak 5S rDNA signals were observed only in proximal regions of two chromosomes, which were adjacent to the (AG₃T₃)₃ proximal signals. Weak (AG₃T₃)₃ signals were observed on both chromosome ends, which enabled accurate chromosome counts. A pair of satellite bodies was observed. (AG₃T₃)₃ signals displayed quite high diversity, changing in intensity from weak to very strong as follows: far away from the chromosome ends (satellites), ends, subtelomeric regions, and proximal regions. Ten high-quality spreads revealed metaphase dynamics from the beginning to the end and the transition to anaphase. Chromosomes gradually grew larger and thicker into linked chromatids, which grew more significantly in width than in length. Based on the combination of 5S rDNA and (AG₃T₃)₃ signal patterns, ten chromosomes were exclusively distinguished, and the remaining twelve chromosomes were divided into two distinct groups. Our physical map, which can reproduce dynamic metaphase progression and distinguish chromosomes, will powerfully guide cytogenetic research on Chimonanthus and other trees.

Fluorescence In Situ Hybridization (FISH) Analysis of the Locations of the Oligonucleotides 5S rDNA, (AGGGTTT)3, and (TTG)6 in Three Genera of Oleaceae and Their Phylogenetic Framework

We report the cytogenetic map for a collection of species in the Oleaceae, and test similarities among the karyotypes relative to their known species phylogeny. The oligonucleotides 5S ribosomal DNA (rDNA), (AGGGTTT)₃, and (TTG)₆ were used as fluorescence in situ hybridization (FISH) probes to locate the corresponding chromosomes in three Oleaceae genera: Fraxinus pennsylvanica, Syringa oblata, Ligustrum lucidum, and Ligustrum × vicaryi. Forty-six small chromosomes were identified in four species. (AGGGTTT)₃ signals were observed on almost all chromosome ends of four species, but (AGGGTTT)₃ played no role in distinguishing the chromosomes but displayed intact chromosomes and could thus be used as a guide for finding chromosome counts. (TTG)₆ and 5S rDNA signals discerned several chromosomes located at subterminal or central regions. Based on the similarity of the signal pattern (mainly in number and location and less in intensity) of the four species, the variations in the 5S rDNA and (TTG)₆ distribution can be ordered as L. lucidum < L. × vicaryi < F. pennsylvanica < S. oblata. Variations have observed in the three genera. The molecular cytogenetic data presented here might serve as a starting point for further larger-scale elucidation of the structure of the Oleaceae genome, and comparison with the known phylogeny of Oleaceae family.

Karyotype Analysis of Diploid and Spontaneously Occurring Tetraploid Blood Orange [Citrus sinensis (L.) Osbeck] Using Multicolor FISH With Repetitive DNA Sequences as Probes

Blood orange [Citrus sinensis (L.) Osbeck] has been increasingly appreciated by consumers worldwide owing to its brilliant red color, abundant anthocyanin and other health-promoting compounds. However, there is still relatively little known about its cytogenetic characteristics, probably because of the small size and similar morphology of metaphase chromosomes and the paucity of chromosomal landmarks. In our previous study, a naturally occurring tetraploid blood orange plant was obtained via seedling screening. Before this tetraploid germplasm can be manipulated into a citrus triploid seedless breeding program, it is of great importance to determine its chromosome characterization and composition. In the present study, an integrated karyotype of blood orange was constructed using sequential multicolor fluorescence in situ hybridization (FISH) with four satellite repeats, two ribosomal DNAs (rDNAs), a centromere-like repeat and an oligonucleotide of telomere repeat (TTTAGGG)3 as probes. Satellite repeats were preferentially located at the terminal regions of the chromosomes of blood orange. Individual somatic chromosome pairs of blood orange were unambiguously identified by repetitive DNA-based multicolor FISH. These probes proved to be effective chromosomal landmarks. The karyotype was formulated as 2n = 2x = 18 = 16m+2sm (1sat) with the karyotype asymmetry degree belonging to 2B. The chromosomal distribution pattern of these repetitive DNAs in this spontaneously occurring tetraploid was identical to that of the diploid, but the tetraploid carried twice the number of hybridization sites as the diploid, indicating a possible pathway involving the spontaneous duplication of chromosome sets in nucellar cells. Our work may facilitate the molecular cytogenetic study of blood orange and provide chromosomal characterization for the future utilization of this tetraploid germplasm in the service of seedless breeding programs.