Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 27 May 2020
Sec. Plant Breeding

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

\r\nHonghong Deng,Honghong Deng1,2Guohao TangGuohao Tang1Nuo XuNuo Xu1Zhijian GaoZhijian Gao1Lijin LinLijin Lin1Dong LiangDong Liang1Hui XiaHui Xia1Qunxian DengQunxian Deng1Jin WangJin Wang1Zexi CaiZexi Cai3Guolu Liang*Guolu Liang2*Xiulan Lv*Xiulan Lv1*
  • 1Institute of Pomology and Olericulture, Sichuan Agricultural University, Chengdu, China
  • 2College of Horticulture and Landscape Architecture, Southwest University, Chongqing, China
  • 3National Maize Improvement Center, College of Agronomy and Biotechnology, China Agricultural University, Beijing, China

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.

Introduction

Citrus is one of the most horticulturally and economically important fruit crops globally (Wu et al., 2018), extensively cultivated in more than 140 countries and regions in the tropics and subtropics (Cuenca et al., 2018). Citrus cultivars are grafted on the rootstock for commercial production worldwide (Oustric et al., 2017, 2019; Bowman and Joubert, 2020). Citrus scion cultivars may respond differently to vegetative growth, canopy size, fruit quality and production, resistance to pests and diseases, as well as tolerance to various abiotic and biotic stresses when grown on diverse rootstocks (Warschefsky et al., 2016; Wang et al., 2017; Oustric et al., 2019; Bowman and Joubert, 2020). Thus, the use of proper rootstock is utmost importance to the success in growing attractive and productive citrus fruits (Albrecht et al., 2019; Bowman and Joubert, 2020).

Carrizo citrange [Citrus sinensis (L.) Osbeck × Poncirus trifoliata (L.) Raf., CC] is a generic hybrid (Deng, 2008) of a Washington navel orange and a trifoliate orange that was obtained almost a century ago (Huerta et al., 2009). Currently, CC remains one of the most commonly used rootstocks in citriculture worldwide, especially in China (Deng, 2008), Spain (Ruiz et al., 2016b), the United States (Albrecht et al., 2012), and several other major citrus-producing countries. The main advantages of CC are excellent tolerance to disease and flooding, intermediate resistance to Phytophthora root rot, great contribution to fruit quality and yield, and good compatibility with citrus scion cultivars (Bowman and Joubert, 2020). For instance, one of the most important strategies for controlling citrus tristeza virus (CTV) employs the utilization of tolerant CC rootstock (Castle et al., 2009). Additionally, CC rootstock has been described as having good cold hardiness for citrus scion cultivars (Oustric et al., 2017). The multibillion-dollar citrus industry is facing a century-old conundrum from Huanglongbing (HLB) (Clark et al., 2018). Albrecht et al. (2012) evaluated the influence of 15 different rootstock varieties on HLB disease development in Florida field-grown Valencia and Early Gold sweet orange [C. sinensis (L.) Osb.], and they found that the highest fruit production, under high HLB pressure, were obtained from the combinations of citrus scion cultivars on US-802 and CC rootstock.

In a citrus rootstock breeding program of our team, an unexpected tetraploid CC genotype was found via seedling screening in our previous study. Tetraploid CC plants exhibit many characteristics that are superior to those of diploid plants with respect to diverse performance, and this is probably due to chromosome doubling (Oustric et al., 2019). Tetraploid CC rootstock improved the chilling stress tolerance of common Clementine mandarin (C. clementina Hort. ex Tan.) scion when compared to the graft of Clementine onto diploid CC rootstock (Oustric et al., 2017). Tetraploid CC trees also present more tolerance toward salinization, hence conferring greater salt tolerance to citrus scion cultivars (Ruiz et al., 2016b). Tetraploidy is associated with an increased boron-excess tolerance in CC rootstock (Ruiz et al., 2016a). In addition, tetraploid CC rootstock is more tolerant to water deficit than the corresponding diploid one (Oliveira et al., 2017). Given the importance of diploid CC plant in citrus production system and various tetraploid superiority, a comprehensive investigation into them is critically needed and of particular interest.

Each species has a characteristic chromosome complement, the karyotype (Jiang, 2019), which represents the highest level of structural and functional organization of genome (Weiss-Schneeweiss and Schneeweiss, 2013; Liang and Chen, 2015). Karyotypic features have greatly facilitated taxonomic and systematic studies, and they have also provided important insight into genome size estimations as well as genome structure and organization at the chromosomal level in numerous plant species (Chèvre et al., 2018), including radish (Raphanus sativus L.) (He et al., 2015), Coix lacryma-jobi L. and Coix aquatica Roxb. (Cai et al., 2014), peanut (Arachis hypogaea) (Zhang et al., 2016), sweet potato (Sun et al., 2019), Oleaceae plants (Luo and Liu, 2019), Clementine mandarin (Deng et al., 2019b), and blood orange [C. sinensis (L.) Osb.] (Deng et al., 2019a). Elucidating chromosome constitution and the karyotype including the number, absolute and relative size, symmetry, and centromere position of the chromosome complement in nucleus cell of individual eukaryotic species, is a fundamental question in plant biology (Liang and Chen, 2015; Jiang, 2019). Nevertheless, up to now, the chromosomal characterization of CC rootstock remains unexplored, partly because of its inherent small size and morphological similarity of mitotic chromosomes as well as the paucity of accessible chromosome landmarks in the Aurantioideae subfamily (Krug, 1943; Silva et al., 2015).

Fluorescence in situ hybridization (FISH) is one of the most powerful and widely used modern molecular cytogenetic techniques with which to perform chromosomal identification and karyotype analysis (Jiang, 2019). The utility of FISH depends largely on the availability of probes that allow the identification of chromosome-specific regions (Jiang, 2019; Sun et al., 2019). Repetitive DNA sequences, including tandem repeats (also known as satellite repeats), 45S ribosomal DNA (rDNA), 5S rDNA, and some oligonucleotides (Garrido-Ramos, 2017), are valuable cytogenetic markers that have been extensively developed as representative FISH probes for chromosome identification in diverse plant species (Cai et al., 2014; He et al., 2015; Jiang, 2019; Sun et al., 2019).

Several excellent reports have emphasized the important role of rDNAs, and some specific genomic fragments, such as bacterial artificial chromosomes (BACs) and rDNAs, have proved particularly useful in identifying chromosomes in some species of Citrus and closely related genera (Brasileiro-Vidal et al., 2007; Moraes et al., 2008; Marques et al., 2011; Mendes et al., 2011; Silva et al., 2013, 2015). As sequencing technology and associated bioinformatics tools have advanced, citrus genomics have followed suit. Recently, we analyzed the repeat sequence composition in the Clementine mandarin genome and successfully developed a set of FISH probes, encompassing a centromeric-like repeat (CL17), four major satellite repeats (CL1, CL2, CL3, and CL4), two rDNAs (45S and 5S rDNA), and a telomeric repeat (TTTAGGG)n (Deng et al., 2019b). In that study, the physical distribution of these repetitive DNAs on the same somatic metaphase chromosomes of Clementine mandarin and other several citrus cultivars were also analyzed (Deng et al., 2019b). Subsequently, a detailed karyotype of diploid blood orange was established by multicolor FISH with these excellent probes, and we did the same for its spontaneously occurring tetraploid plant (Deng et al., 2019a).

The present study is a continuation of our recent two papers (Deng et al., 2019a, b). This study aimed to test the feasibility of repetitive DNA-based FISH on chromosome discrimination of citrus rootstock cultivars to establish the integrated somatic metaphase chromosome karyotypes of diploid and tetraploid CC rootstocks and to analyze the chromosome complement of the tetraploid CC plant. Herein, this follow-up study applied sequential multicolor FISH with the aforementioned repetitive DNAs and proposed the FISH-based karyotypes of diploid CC rootstock and its spontaneously occurring tetraploid. The results presented extend the limited knowledge of molecular cytological characterization of this extensively used citrus rootstock. Moreover, the present study also provides the fundamental basis for future exploration of genome organization of CC rootstock and intricate polyploid formation.

Materials and Methods

Plant Materials

The tetraploid CC plant used in this study was originally obtained from a seedling screening in the annual laboratory activity of our team. The diploid and tetraploid plants were transplanted into the field simultaneously and maintained under identical natural field conditions at the experimental base of the College of Horticulture and Landscape Architecture, Southwest University, Chongqing, China (latitude 23.39°N, longitude 34.95°E). The ploidy status of the plant materials was verified by flow cytometry in combination with somatic chromosome counting before the seedlings were transplanted.

Somatic Chromosome Preparation

Somatic chromosome preparation was made following our previously published protocol with minor optimizations (Deng et al., 2019a). Briefly, at least 20 fresh and fine roots (approximately 3–4 cm proximal to the root tip) were carefully dissected from plants with a tweezer on a sunny morning. These excised root tips were rinsed in tap water and pretreated in ice-cold water for 24 h in dark, and this was followed by fixation in a freshly prepared 3:1 Carnoy’s fixative solution of absolute ethanol (EtOH) to glacial acetic acid for at least 24 h at room temperature (RT). The treated samples were then stored in 70% EtOH at −20°C until further use. The fixed root tips were thoroughly washed with distilled water and citrate buffer (0.01 M citric acid and 0.01 M sodium citrate buffer, pH 4.8). Thereafter, actively growing root apexes (approximately 1.5–2.5 mm long) with meristematic cells were transferred into an enzyme solution consisting of 3% (w/v) cellulose Onozuka R-10 (Yakult Pharmaceutical Industry Co., Ltd., Tokyo, Japan) plus 0.3% (w/v) pectolyase Y-23 (Kikkoman Corp., Tokyo, Japan) and digested for 1.5–2.0 h at 37°C. After enzyme treatment, the softened material was carefully rinsed in distilled water to remove excess enzymes and transferred into a drop of freshly prepared 3:1 Carnoy’s fixative solution. Each root was carefully transferred to a microscope slide and immediately squashed using a tweezer. Finally, a small drop of fixative was added to the slide and flame dried. The chromosomes of mitotic metaphase cells were examined under phase contrast microscope, and good-quality chromosomes at metaphase stage were selected for further sequential multicolor FISH experiments.

DNA Probe Labeling and Sequential Multicolor FISH

The following seven probes, encompassing 45S rDNA, 5S rDNA, four satellite repeats (CL1, CL2, CL3, and CL4), and a centromere-like repeat (CL17), were used for sequential multicolor FISH. These seven probes, representing repetitive chromosomal markers, were previously amplified from genomic DNA of Clementine mandarin by PCR using primer pairs as detailed by Deng et al. (2019b). A digoxigenin-labeled oligonucleotide of (TTTAGGG)3 (Nelson et al., 2014) was directly synthesized to map telomeric sequences. Details of DNA probe labeling via nick translation and sequential multicolor FISH targeting all repetitive elements followed our reported protocol without modification (Deng et al., 2019a). Three rounds of sequential FISH experiments were performed on the same well-spread mitotic metaphase chromosomes of diploid and tetraploid CC plants. Probes (2 μL) were equally mixed before hybridization and added to the hybridization mixture, including 10 μL deionized formamide, 4 μL 50% dextran sulfate, 2 μL 20 × SSC, and 2 μL sheared salmon sperm blocking DNA (10 mg/mL; Sigma-Aldrich, St. Louis, MO, United States). A total volume of 20 μL hybridization mixture was added to each slide, and the slides were covered with coverslips of appropriate sizes. Probes and chromosomes were denatured together on a hot plate at 85°C for 1.5 min and incubated in a 37°C moist chamber at least 6 h or overnight. The procedure of FISH signal detection and posthybridization washes described by Deng et al. (2019a) without modifications were used. The slides were mounted with Vectashield§ antifade mounting medium containing 1.5 μg/mL 4′,6-diamidino-2-phenylindole (DAPI) (Vector Laboratories, Inc., Burlingame, CA, United States).

Microscopy and Image Processing

The fluorescence signals were visualized using an Olympus epifluorescence microscope (BX61; Tokyo, Japan) at the same magnification times, and the photomicrographs of FISH chromosomes were captured using a Sensys CCD camera (Qimaging RetigaTM SRV Fast 1394, Vancouver, BC, Canada) under the same exposure intensity and time. Images were preudocolored (blue for DAPI, red for anti-digoxigenin-rhodamin, green for anti-avidin antibody conjugated with FITC, yellow for diethylaminocoumarin-5-dUTP), inverted and merged using ImageJ software (National Institutes of Health, Wayne Rasband, MD, United States), and further cropped and processed with Adobe Photoshop CS6.0 software (Adobe Systems Incorporated, San Jose, CA, United States) using only functions affecting the entire image simultaneously and equally.

Statistical Analysis of Karyomorphological Indices

The chromosome number was estimated based on at least 30 split phases with well-separated and good-quality metaphase spreads. All the chromosome indices and karyotype parameters were calculated based on the mean values of at least five better split phase through KaryoType software (Altinordu et al., 2016). A total of 13 independent karyotype parameters, including long arm length, short arm length, arm ratio, relative length of chromosome length, and centromeric index, were calculated for both diploid and tetraploid CC plants. Chromosome classifications were determined based on definitions provided by Levan et al. (1964). A karyotype idiogram was produced based on the decreasing relative chromosome lengths and centromere positions, or a combination of these two aspects of chromosome morphology, as well as the position and distribution of other chromosomal landmarks, in reference to previous reports of Cai et al. (2014), He et al. (2015), and (Deng et al., 2019a, b).

Results

Physical Localization of Highly Repetitive DNA Sequences on Somatic Chromosomes of Carrizo Citrange

In order to precisely identify each pair of homologous chromosomes and obtain detailed karyotypes of diploid and tetraploid CC plants, the physical distribution of highly repetitive DNAs on somatic chromosomes was determined by sequential multicolor FISH on the mitotic metaphase chromosomes coming from the same metaphase spreads (Figures 1, 2). Some other additional metaphase plates produced consistent results with respect to chromosome counts, and both signal numbers and locations are used for karyotype analysis, as shown in Supplementary Table S1. Somatic chromosome images at mitotic metaphase stages of the analyzed diploid and tetraploid CC plants are shown in Figures 1, 2, respectively. In general, the different probes distributed to various positions and showed great variation in the numbers and fluorescence intensities of hybridization sites (Figures 1, 2). FISH signals of highly repetitive DNAs in the tetraploid CC plants (Figure 2) corresponded quantitatively to the double signals of its metaphase counterpart in the diploid CC rootstock (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1. Physical localization of different highly repetitive DNA sequences on root tip metaphase chromosomes of diploid Carrizo Citrange rootstock. The merged images of (a), (e), and (i) came from the hybridization sites in (b–d), (f–h), and (j,k), respectively. (b) Satellite repeat CL1 (green), (c) a centromere-like repeat CL17 (red), (d) satellite repeat CL3 (yellow), (f) satellite repeat CL2 (green), (g) 5S rDNA (red), (h) satellite repeat CL4 (yellow), (j) 45S rDNA (green), and (k) telomere repeat (red) were used as FISH probes. Signal color corresponds to probe name in (b–d,f–h,j,k). (l) DAPI-stained chromosomes were converted into a blue image. Bar = 2 μm.

FIGURE 2
www.frontiersin.org

Figure 2. Physical localization of different highly repetitive DNA sequences on root tip metaphase chromosomes of tetraploid Carrizo Citrange rootstock. The merged images of (a), (e), and (i) came from the hybridization sites in (b–d), (f–h), and (j,k), respectively. (b) Satellite repeat CL1 (green), (c) a centromere-like repeat CL17 (red), (d) satellite repeat CL3 (yellow), (f) satellite repeat CL2 (green), (g) 5S rDNA (red), (h) satellite repeat CL4 (yellow), (j) 45S rDNA (green), and (k) telomere repeat (red) were used as FISH probes. Signal color corresponds to probe name in (b–d,f–h,j,k). (l) DAPI-stained chromosomes were converted into a blue image. Bar = 2 μm.

CL1 (Figure 1B), CL2 (Figure 1F), CL3 (Figure 1D), and CL4 (Figure 1H) signals were located predominantly in subtelomeric regions of several homologous chromosome pairs, with 11, 12, 11, and 4 hybridization sites being observed in diploid CC rootstock, respectively (Figure 1). Specifically, FISH with a CL1 probe revealed an asymmetrical pattern that the CL1 loci were localized in the subtelomeric regions of the long arms of the chromosome pairs 2, 3, 4, 5, and 6, and only one short arm of the chromosome pair 3 (Figures 1B, 2B). Positive signals of CL2 were visible at the subtelomeric regions of the long arms of the chromosome pairs 2, 3, 4, 5, and 6, and the two short arms of chromosome pair 3 (Figures 1F, 2F). CL3 FISH signals (Figures 1D, 2D) were co-localized with CL2 locus. CL4 FISH signals were detected in the subtelomeric regions of two long arms of chromosome pair 4, only one short arm of the chromosome pair 3 co-localized with other three satellite repeats, and only one long arm of the chromosome pair 7 (Figures 1H, 2H).

The CL17 probe hybridized to sites on chromosome centromeric regions (Figures 1C, 2C). FISH with the oligonucleotides of (TTTAGGG)3 showed noticeable twin fluorescent signals at the terminal regions of all chromosomes in both diploid (Figure 1K) and tetraploid CC rootstock (Figure 2K). A few weak signals of the oligonucleotides of (TTTAGGG)3 were also detected in interstitial regions of several chromosome pairs, but there were much less intense than those in the terminal regions (Figures 1K, 2K). FISH with 45S rDNA showed four hybridization sites in diploid CC rootstock, with two sites on the short arms of chromosome pair 2 corresponding to one pair of nucleolus organizer regions (NORs), one site on one short arm of chromosome pair 3 co-localized with other four satellite repeats, and one site occurred on the chromosomal satellite position next to the short arm of chromosome pair 8 (Figure 1J). Two 5S rDNA sites were identified in the pericentromeric regions of the short arms of chromosome pair 2 co-carrying 45S rDNA locus, and another two 5S rDNA FISH signals were found in the pericentromeric regions of the short arms of chromosome pair 8 (Figure 1G). The tetraploid CC rootstock had the same rDNA distribution pattern but twice as many loci as the diploid CC rootstock (Figures 2J, 2G). Chromosomes stained by DAPI are shown in Figures 1L, 2L, from which we can see it is hard to recognize the chromosome pair for the small size and high similarity in morphology without the aids of chromosome landmarks.

Characterizing Karyotype of Diploid and Tetraploid Carrizo Citrange Chromosomes

Chromosome numbers were confirmed by chromosome counts in at least 30 mitotic metaphases. All of the examined mitotic metaphase plates of CC rootstock showed a fundamental number x = 9 (Figures 1, 2). Karyotypes of CC individuals possessed the same standard diploid constitution 2n = 18 (Figure 1) and tetraploid constitution 4n = 36 (Figure 2). The karyotype parameters, chromosome classification, and morphological features were detailed in Supplementary Table S1. The karyotype formula of CC rootstock was 2n = 2x = 12m (1sat) + 6sm. The chromosome complement consisted of three pairs of submetacentric chromosomes and six pairs of metacentric chromosomes, with one chromosome being satellite chromosome. Concerning Stebbin’s classification of karyotypes, CC rootstock falls into the 2B category (Supplementary Table S1). The homologous chromosomes were paired based on chromosome size and morphology, centromere position, and the presence and distribution of repetitive DNA loci, as shown in Figures 1, 2, collectively facilitating the integration of karyotype of CC rootstock (Figure 3 and Supplementary Table S1).

FIGURE 3
www.frontiersin.org

Figure 3. Idiogram of diploid Carrizo Citrange rootstock showing the chromosome complement and relative positions of highly repetitive DNA sequences. Numbers on the x- and y-axis indicate the numbers of homologous chromosome pairs and the relative arm length of chromosome pairs, respectively. The ideogram is construed based on the sequential multicolor FISH-based karyotype. Chromosomes are arranged in matched, homologous pairs and ordered by the entire length from the longest to shortest and centromere positions, and a combination of these two aspects of chromosome morphology. The color scheme at right side shows the color of each probe as represented in this idiogram.

From Figure 3, the karyotype results convincingly showed that diploid CC rootstock is a hybrid. In the karyotype, chromosomes were ordered by a gradual decrease excluding the satellite length in computing arm ratios. Thus, chromosome pairs 1 and 9 can be recognized by the largest and smallest size among all the homologous chromosome pairs, respectively. It is evident that chromosome pair 2 carried a pair of strong 45S rDNA loci co-located with 5S rDNA at the NOR of short arms and three pairs of satellite DNA repeat (CL1, CL2, and CL3) loci at the terminal regions of long arms. Chromosome pair 3 was characterized by a pair of robust CL2 and CL3 loci on the terminal positions of both short and long arms; one 45S rDNA locus, one CL1 locus, and one CL4 locus were terminally positioned in only one short arm of chromosome pair 3 near the CL2 and CL3 loci. Chromosome pair 4 could be easily distinguished by four pairs of satellite repeat loci on the long arms. Morphologically, chromosome pairs 5 and 6 resembled each other and presented symmetrical CL1, CL2, and CL3 signals at the subtelomeric regions of long arms. Only one long arm of chromosome pair 7 contained a CL4 locus. Chromosome 8 showed strong 5S rDNA signals in the pericentromeric regions of long arms, with one satellite chromosome being labeling by intense 45S rDNA signal next to the long arm. The corresponding tetraploid CC chromosome pairs could also be identified for these chromosomal landmarks, and, therefore, only the idiogram of diploid CC is presented.

Discussion

The small size (between 2 and 4 μm) and similar morphology of Citrus chromosomes (Krug, 1943) has been long challenging unequivocal identification of homologous chromosome pairs (Silva et al., 2015). Krug (1943) investigated the first chromosome characterization of Citrus. Since then, chromosome banding techniques, such as C-banding (Liang, 1988) and CMA+/DAPE fluorescence banding (Carvalho et al., 2005; Brasileiro-Vidal et al., 2007; Guerra and Souza, 2009; Silva et al., 2010; Marques et al., 2011; Mendes et al., 2011), and the combination of CMA banding and rDNA- and BAC-based FISH (Carvalho et al., 2005; Brasileiro-Vidal et al., 2007; de Moraes et al., 2007; Guerra and Souza, 2009; Silva et al., 2010, 2013; Marques et al., 2011; Mendes et al., 2011) have gradually been developed for citrus cytogenetic study over time.

Taking advantage of the availability of citrus reference genome sequences and using the newly developed bioinformatic tools, we have investigated the content and the types as well as the proportional contribution of repetitive DNAs in C. clementina mandarin genome (Deng et al., 2019b). The chromosome distribution patterns of these highly repetitive DNAs on mitotic chromosomes of C. clementina mandarin and several citrus scion cultivars within Citrus, Poncirus, and Fortunella were analyzed through FISH (Deng et al., 2019a, b). In the present study, we extended the FISH-based karyotype analysis to one of the most widely used citrus rootstocks and its spontaneously occurring tetraploid plant. While the previous karyotype studies focused on citrus scion variety, the karyotype of citrus rootstock and its corresponding tetraploid plants were studied for the first time in the present study.

FISH technique has been successfully applied to identify specific chromosome pair or chromosome regions in a number of plant species (Cai et al., 2014; He et al., 2015; Zhang et al., 2016; Deng et al., 2019a, b; Jiang, 2019; Luo and Liu, 2019; Sun et al., 2019). The results in the present study confirmed the importance and advantage of using repetitive DNA-based FISH for chromosome discrimination and cytogenetic characterization in CC rootstock. Accordingly, the chromosome complement of CC rootstock was clearly unraveled (Figures 13). The chromosome numbers (2n = 18, Supplementary Table S1) were in accordance with previously published citrus diploid data (Krug, 1943; Silva et al., 2010; Marques et al., 2011).

The repetitive DNA probes developed in our previous work (Deng et al., 2019b) were of great assistance in investigating the cytological characterization of diverse citrus scion cultivars (Deng et al., 2019a). Here, when this set of eight such probes was applied to a sequential multicolor FISH-based karyotype of CC rootstock, and the results collaborated that eight such repetitive DNAs are good chromosomal landmarks for citrus rootstock species (Figures 13). In this cytogenetic methodology, it is able to simultaneously reveal the diverse features of the small or difficult-to-distinguish chromosomes of citrus rootstock. Repetitive DNAs are responsible for the main components of the structural and functional organization of higher plant genomes, resulting in the genome size variation (Guerra et al., 2000; Garrido-Ramos, 2017). The tandemly repetitive DNA sequences exhibit characteristic chromosomal locations, usually at pericentromeric, centromeric, and subtelomeric regions of eukaryotic chromosomes (Cai et al., 2014; He et al., 2015). CL17 loci are valuable chromosome centromere markers for CC rootstock, and the physical mapping of CL17 loci using in situ hybridization allows for the determination of centromere position of chromosome pairs in CC rootstock (Figures 1C, 2C).

The repetitive DNAs could be specific for a family or genus taxonomically or may be specific for a plant species, genome or even a chromosome pair (Jiang, 2019). Thus, the distribution of repetitive DNAs may vary qualitatively and quantitatively between citrus taxa (Guerra et al., 2000). The genome changes may be intriguingly associated with elimination and amplification of some types of repetitive DNAs, such as tandemly arranged DNA sequences (rDNA and satellite repeats) or dispersed (DNA transposons and retrotransposons) (Garrido-Ramos, 2017). In this study, satellite repeats (CL1, CL2, CL3, and CL4) showed subtelomeric preferences in the root tip metaphase chromosome material of both diploid (Figures 1B,D,F,H) and tetraploid (Figures 2B,D,F,H) CC rootstock, which has been considered as a common and typical feature of satellite repeat in Citrus L (Silva et al., 2010; Marques et al., 2011; Deng et al., 2019a, b).

Subtelomere, the regions proximal to telomeres, of most eukaryotes harbor fast-evolving gene family usually in relation to adaptive processes (Brown et al., 2010). Specifically speaking, the subtelomere of common bean (Phaseolus vulgaris) favors the rapid evolution of genes associated with resistance (Chen et al., 2018). Satellite repeat represents a fast-evolving portion of various plant species (Cai et al., 2014; He et al., 2015). Therefore, it must be extremely interesting to understand what the possible biological functions behind this preferential subtelomere localizations of satellite repeats in citrus in the near future and how the satellites have undergone rapid evolution within Citrus L. genus.

Another major tandem repetitive DNA in CC genome is rDNA (Figures 1G,J, 2G,J). Results of the cytogenetic studies have been consistent in previous views that rDNA is a very useful marker to identify chromosomes and species (Silva et al., 2010; Marques et al., 2011; Deng et al., 2019a, b), as rDNAs are presented in numerous repeated units, generally organized in clusters and are thus easily visualized in plant chromosomes (Garrido-Ramos, 2017). In our study, the rDNA hybridization signals presented in only one of the homologous chromosomes (Figures 1G,J, 2G,J) explicitly supported the generic hybrid nature of CC rootstock (Deng, 2008) from cytological aspect. The 45S rDNA, generally confined to NOR (Marques et al., 2011; Lan et al., 2016), was consistently located proximally on the short arms of chromosome pair 2 in our study (Figures 1J, 2J), which may have potential to illustrate several issues of chromosomal structure and composition. The plant-type telomere probe (TTTAGGG)3 (Nelson et al., 2014) efficiently hybridized to the ends of all chromosomes (Figures 1K, 2K), indicating the presence of TTTAGGG repeat in CC telomere.

Chromosome doubling may generate during some period of seed formation in Citrus L., and tetraploid seedling may arise spontaneously in the nucellar cells. These clones are regarded to be genetically identical because they have the same genome expression profile (Cameron and Soost, 1970; Deng et al., 2019a). The distribution pattern of repetitive DNAs in the tetraploid metaphase chromosomes (Figure 2) is concordant with those in the diploid CC plant (Figure 1). The number of repetitive DNAs sites strictly correlated with the ploidy, and the tetraploid (Figure 2) possessed twice the number of sites as the diploid (Figure 1). Previous molecular study using SSR markers has been applied to hundreds of spontaneous tetraploid citrus lines, and showed that in the absence of tetraploid pollinator spontaneous tetraploid lines resulted from somatic chromosome doubling and that the rates were affected by both genotype and the environment (Aleza et al., 2011). These results presented herein confirmed that the naturally occurring tetraploid CC material arose from somatic chromosome doubling from the cytogenetic aspect. In apomictic citrus germplasm, tetraploidization is a common event, and the autotetraploids arising from chromosome doubling are extensively used as male parents in citrus improvement (Aleza et al., 2011). Further molecular cytogenetic studies, including the parental species of CC rootstock and comparative studies regarding the genome sequencing, are highly needed to elucidate a more accurate genome structure of this genetic hybrid. This results in the present study have provided valuable cytogenetic information for CC rootstock, which could be a good starting point for future chromosome engineering-based citrus rootstock improvement.

Conclusion

In conclusion, sequential multicolor FISH with three sets of efficient probes containing a centromere-like repeat, four satellite repeats, two rDNAs, and an oligonucleotide of telomeric (TTTAGGG)3 repeat were conducted on the same root tip metaphase chromosomes, by which each chromosome pair was unequivocally identified. The asymmetric distribution patterns of the repetitive DNA sequences demonstrated the hybrid origin of CC rootstock. This study revealed, for the first time, the integrated karyotype and chromosomal characteristics of one most widely used citrus rootstock, Carrizo citrange (Citrus sinensis L. Osbeck × Poncirus trifoliata L. Raf.), as well as its spontaneously occurring tetraploid plant.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Author Contributions

HD, GL, and XL conceived the research, reviewed and edited the manuscript, and acquired funding. HD, LL, DL, HX, QD, JW, and ZC performed the methodology optimization and data curation. HD, GT, NX, and ZG conducted the experiment and managed the plant material collection. HD wrote the manuscript. All authors read and approved the final manuscript.

Funding

This work was supported by the Program for Chongqing Graduate Student Innovation (CYB14050) and the China Scholarship Council (No. 201606990048) to HD and Pujiang (Sichuan, China) Modern Agriculture Industrialization Exemplary Project to QD and XL.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

HD expresses her sincere heartfelt thanks to her Ph.D. co-supervisor, Dr. Weiwei Jin at the China Agricultural University, for his excellent supervision, generous assistance, great friendship, and continuous support on various aspects of her study, work, and life. Thanks for technical assistance kindly provided by Dr. Huijun Liu at China Agricultural University. We would like to thank the reviewers and the editor for their critical review of the manuscript and beneficial comments.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.00569/full#supplementary-material

References

Albrecht, U., McCollum, G., and Bowman, K. D. (2012). Influence of rootstock variety on Huanglongbing disease development in field-grown sweet orange (Citrus sinensis [L.] Osbeck) trees. Sci. Hortic. 138, 210–220. doi: 10.1016/j.scienta.2012.02.027

CrossRef Full Text | Google Scholar

Albrecht, U., Tripathi, I., Kim, H., and Bowman, K. D. (2019). Rootstock effects on metabolite composition in leaves and roots of young navel orange (Citrus sinensis L. Osbeck) and pummelo (C. grandis L. Osbeck) trees. Trees 33, 243–256. doi: 10.1007/s00468-018-1773-1

CrossRef Full Text | Google Scholar

Aleza, P., Froelicher, Y., Schwarz, S., Agusti, M., Hernandez, M., Juarez, J., et al. (2011). Tetraploidization events by chromosome doubling of nucellar cells are frequent in apomictic citrus and are dependant of genotype and environment. Ann. Bot. 108, 37–50. doi: 10.1093/aob/mcr099

CrossRef Full Text | Google Scholar

Altinordu, F., Peruzzi, L., Yu, Y., and He, X. (2016). A tool for the analysis of chromosomes: karyoType. Taxon 65, 586–592. doi: 10.12705/653.9

CrossRef Full Text | Google Scholar

Bowman, K. D., and Joubert, J. (2020). “Citrus rootstocks,” in The Genus Citurs, eds M. Talon, M. Caruso, and F. G. Gmitter (Amsterdam: Elsevier Inc), 105–127.

Google Scholar

Brasileiro-Vidal, A. C., Dos Santos-Serejo, J. A., Soares Filho, W. D. S., and Guerra, M. (2007). A simple chromosomal marker can reliably distinguishes Poncirus from Citrus species. Genetica 129, 273–279. doi: 10.1007/s10709-006-0007-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, C. A., Murray, A. W., and Verstrepen, K. J. (2010). Rapid expansion and functional divergence of subtelomeric gene families in yeasts. Curr. Biol. 20, 895–903. doi: 10.1016/J.CUB.2010.04.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, Z., Liu, H., He, Q., Pu, M., Chen, J., Lai, J., et al. (2014). Differential genome evolution and speciation of Coix lacryma-jobi L. and Coix aquatica Roxb. hybrid guangxi revealed by repetitive sequence analysis and fine karyotyping. BMC Genomics 15:1025. doi: 10.1186/1471-2164-15-1025

CrossRef Full Text | Google Scholar

Cameron, J. W., and Soost, R. K. (1970). “Characters of new populations of Citrus polyploids, and the relation between tetraploidy in the pollen parent and hybrid tetraploid progeny,” in Proceedings of the 1st International Citrus Symposium, Riverside, CA, 199–205.

Google Scholar

Carvalho, R., Soares Filho, W. S., Brasileiro-Vidal, A. C., and Guerra, M. (2005). The relationships among lemons, limes and citron: a chromosomal comparison. Cytogenet. Genome Res. 109, 276–282. doi: 10.1159/000082410

PubMed Abstract | CrossRef Full Text | Google Scholar

Castle, W. S., Nunnallee, J., and Manthey, J. A. (2009). Screening citrus rootstocks and related selections in soil and solution culture for tolerance to low-iron stress. HortScience 44, 638–645. doi: 10.21273/HORTSCI.44.3.638

CrossRef Full Text | Google Scholar

Chen, N. W. G., Thareau, V., Ribeiro, T., Magdelenat, G., Ashfield, T., Innes, R. W., et al. (2018). Common bean subtelomeres are hot spots of recombination and favor resistance gene evolution. Front. Plant Sci. 9:1185. doi: 10.3389/fpls.2018.01185

PubMed Abstract | CrossRef Full Text | Google Scholar

Chèvre, A.-M., Mason, A. S., Coriton, O., Grandont, L., Jenczewski, E., and Lysak, M. A. (2018). “Cytogenetics, a science linking genomics and breeding: the Brassica model,” in The Brassica napus Genome, eds S. Liu, R. Snowdon, and B. Chalhoub (Cham: Springer), 21–39.

Google Scholar

Clark, K., Franco, J. Y., Schwizer, S., Pang, Z., Hawara, E., Liebrand, T. W. H., et al. (2018). An effector from the huanglongbing-associated pathogen targets citrus proteases. Nat. Commun. 9:1718. doi: 10.1038/s41467-018-04140-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuenca, J., Garcia-Lor, A., Navarro, L., and Aleza, P. (2018). “Citrus genetics and breeding,” in Advances in Plant Breeding Strategies: Fruits, eds J. M. Al-Khayri, S. M. Jain, and D. V. Johnson (Cham: Springer International Publishing), 403–436.

Google Scholar

de Moraes, A. P., Filho, W. D. S. S., and Guerra, M. (2007). Karyotype diversity and the origin of grapefruit. Chromosom. Res. 15, 115–121. doi: 10.1007/s10577-006-1101-2

CrossRef Full Text | Google Scholar

Deng, H., Cai, Z., Xiang, S., Guo, Q., Huang, W., and Liang, G. (2019a). Karyotype analysis of diploid and spontaneously occurring tetraploid blood orange (Citrus sinensis (L.) Osbeck) using multicolor FISH with repetitive DNA sequences as probes. Front. Plant Sci. 10:331. doi: 10.3389/FPLS.2019.00331

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, H., Xiang, S., Guo, Q., Jin, W., Cai, Z., and Liang, G. (2019b). Molecular cytogenetic analysis of genome-specific repetitive elements in Citrus clementina Hort. Ex Tan. and its taxonomic implications. BMC Plant Biol. 19:77. doi: 10.1186/s12870-019-1676-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, X. (2008). “Trifoliate orange and ity hybrids with Citrus genus,” in Citrus varieties in China (in Chinese), eds X. Deng, C. Peng, Z. Chen, Z. Deng, J. Xu, and J. Li (Beijing: China Agriculture Press), 186–189.

Google Scholar

Garrido-Ramos, M. (2017). Satellite DNA: an evolving topic. Genes 8, 230. doi: 10.3390/genes8090230

CrossRef Full Text | Google Scholar

Guerra, M., Galvão Bezerra Dos Santos, K., Barros, E., Silva, A. E., and Ehrendorfer, F. (2000). Heterochromatin banding patterns in rutaceae-aurantioideae - A case of parallel chromosomal evolution. Am. J. Bot. 87, 735–747. doi: 10.2307/2656860

PubMed Abstract | CrossRef Full Text | Google Scholar

Guerra, M., and Souza, M. (2009). “Chromosomal variability and the origin of Citrus species,” in Genetic Diversity, ed. C. L. Mohoney and D. A. Springer (New York, NY: Nova Science Publisher, Inc), 51–68.

Google Scholar

He, Q., Cai, Z., Hu, T., Liu, H., Bao, C., Mao, W., et al. (2015). Repetitive sequence analysis and karyotyping reveals centromere-associated DNA sequences in radish (Raphanus sativus L.). BMC Plant Biol. 15:105. doi: 10.1186/s12870-015-0480-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Huerta, L., Garcia-Lor, A., and Garcia-Martinez, J. L. (2009). Characterization of gibberellin 20-oxidases in the citrus hybrid Carrizo citrange. Tree Physiol. 29, 569–577. doi: 10.1093/treephys/tpn049

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, J. (2019). Fluorescence in situ hybridization in plants: recent developments and future applications. Chromosom. Res. 27, 153–165. doi: 10.1007/s10577-019-09607-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Krug, C. A. (1943). Chromosome numbers in the subfamily aurantioideae with special reference to the genus Citrus. Botanical Gazette 104, 602–611. doi: 10.1086/335173

CrossRef Full Text | Google Scholar

Lan, H., Chen, C.-L., Miao, Y., Yu, C.-X., Guo, W.-W., Xu, Q., et al. (2016). Fragile sites of ‘Valencia’ sweet orange (Citrus sinensis) chromosomes are related with active 45s rDNA. PLoS One 11:e0151512. doi: 10.1371/journal.pone.0151512

CrossRef Full Text | Google Scholar

Levan, A., Fredga, K., and Sandberg, A. A. (1964). Nomenclature for centromeric position on chromosomes. Hereditas 52, 201–220. doi: 10.1111/j.1601-5223.1964.tb01953.x

CrossRef Full Text | Google Scholar

Liang, G. (1988). Studies on the giemsa C-banding patterns of some Citrus and its related genera. Acta Genet. Sin. 15, 409–415.

Google Scholar

Liang, G., and Chen, H. (2015). Scaling chromosomes for an evolutionary karyotype: a chromosomal tradeoff between size and number across woody species. PLoS One 10:e0144669. doi: 10.1371/journal.pone.0144669

CrossRef Full Text | Google Scholar

Luo, X., and Liu, J. (2019). 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. Genes 10:375. doi: 10.3390/genes10050375

CrossRef Full Text | Google Scholar

Marques, A., Fuchs, J., Ma, L., Heckmann, S., Guerra, M., and Houben, A. (2011). Characterization of eu-and heterochromatin of Citrus with a focus on the condensation behavior of 45S rDNA chromatin. Cytogenet. Genome Res. 134, 72–82. doi: 10.1159/000323971

CrossRef Full Text | Google Scholar

Mendes, S., Moraes, A. P., Mirkov, T. E., and Pedrosa-Harand, A. (2011). Chromosome homeologies and high variation in heterochromatin distribution between Citrus L. and Poncirus Raf. as evidenced by comparative cytogenetic mapping. Chromosom. Res. 19, 521–530. doi: 10.1007/s10577-011-9203-x

CrossRef Full Text | Google Scholar

Moraes, A. P., Mirkov, T. E., and Guerra, M. (2008). Mapping the chromosomes of Poncirus trifoliata Raf. by BAC-FISH. Cytogenet. Genome Res. 281, 277–281. doi: 10.1159/000138897

CrossRef Full Text | Google Scholar

Nelson, A. D. L., Beilstein, M. A., and Shippen, D. E. (2014). “Plant telomeres and telomerase,” in Molecular Biology, The Plant Sciences 2, ed. S. H. Howell (New York, NY: Springer), 25–49.

Google Scholar

Oliveira, T. M., Yahmed, J. B., Dutra, J., Maserti, B. E., Talon, M., and Navarro, L. (2017). Better tolerance to water deficit in doubled diploid ‘Carrizo citrange’ compared to diploid seedlings is associated with more limited water consumption. Acta Physiol. Plant 39:204. doi: 10.1007/s11738-017-2497-3

CrossRef Full Text | Google Scholar

Oustric, J., Morillon, R., Luro, F., Herbette, S., Lourkisti, R., Giannettini, J., et al. (2017). Tetraploid Carrizo citrange rootstock (Citrus sinensis Osb. × Poncirus trifoliata L. Raf.) enhances natural chilling stress tolerance of common Clementine (Citrus clementina Hort. ex Tan). J. Plant Physiol. 214, 108–115. doi: 10.1016/j.jplph.2017.04.014

CrossRef Full Text | Google Scholar

Oustric, J., Morillon, R., Luro, F., Herbette, S., Martin, P., Giannettini, J., et al. (2019). Nutrient deficiency tolerance in citrus is dependent on genotype or ploidy level. Front. Plant Sci. 10:127. doi: 10.3389/fpls.2019.00127

CrossRef Full Text | Google Scholar

Ruiz, M., Quiñones, A., Martínez-Alcántara, B., Aleza, P., Morillon, R., Navarro, L., et al. (2016a). Tetraploidy enhances boron-excess tolerance in Carrizo Citrange (Citrus sinensis L. Osb. × Poncirus trifoliata L. Raf.). Front. Plant Sci. 7:701. doi: 10.3389/fpls.2016.00701

CrossRef Full Text | Google Scholar

Ruiz, M., Quiñones, A., Martínez-Cuenca, M. R., Aleza, P., Morillon, R., Navarro, L., et al. (2016b). Tetraploidy enhances the ability to exclude chloride from leaves in carrizo citrange seedlings. J. Plant Physiol. 205, 1–10. doi: 10.1016/j.jplph.2016.08.002

CrossRef Full Text | Google Scholar

Silva, A. E., Barros, S. E., Filho, W. S. S., and Guerra, M. (2013). Linked 5S and 45S rDNA sites are highly conserved through the subfamily Aurantioideae (Rutaceae). Cytogenet. Genome Res. 140, 62–69. doi: 10.1159/000350695

CrossRef Full Text | Google Scholar

Silva, A. E. B., Marques, A., Santos, K. G. B., and Guerra, M. (2010). The evolution of CMA bands in Citrus and related genera. Chromosom. Res. 18, 503–514. doi: 10.1007/s10577-010-9130-2

CrossRef Full Text | Google Scholar

Silva, S. C., Mendes, S., Filho, W. S. S., and Pedrosa-Harand, A. (2015). Chromosome homologies between Citrus and Poncirus—the comparative cytogenetic map of mandarin (Citrus reticulata). Tree Genet. Genomes 11:811. doi: 10.1007/s11295-014-0811-4

CrossRef Full Text | Google Scholar

Sun, J., Yu, L., Cai, Z., Zhang, A., Jin, W., Han, Y., et al. (2019). Comparative karyotype analysis among six Ipomoea species based on two newly identified repetitive sequences. Genome. 64, 243–253. doi: 10.1139/gen-2018-0169

CrossRef Full Text | Google Scholar

Wang, J., Jiang, L., and Wu, R. (2017). Plant grafting: how genetic exchange promotes vascular reconnection. New Phytol. 214, 56–65. doi: 10.1111/nph.14383

CrossRef Full Text | Google Scholar

Warschefsky, E. J., Klein, L. L., Frank, M. H., Chitwood, D. H., Londo, J. P., von Wettberg, E. J. B., et al. (2016). Rootstocks: diversity, domestication, and impacts on shoot phenotypes. Trends Plant Sci. 21, 418–437. doi: 10.1016/J.TPLANTS.2015.11.008

CrossRef Full Text | Google Scholar

Weiss-Schneeweiss, H., and Schneeweiss, G. M. (2013). “Karyotype diversity and evolutionary trends in angiosperms,” in Plant Genome Diversity (Volume 2), Plant Structure, Beheavior and Evolution of Plant Genmones, eds I. J. Leitch, J. Greilhuber, J. Dolezel, and J. F. Wendel (Berlin: Springer), 209–230.

Google Scholar

Wu, G. A., Terol, J., Ibanez, V., López-García, A., Pérez-Román, E., Borredá, C., et al. (2018). Genomics of the origin and evolution of Citrus. Nature 554, 311–316. doi: 10.1038/nature25447

CrossRef Full Text | Google Scholar

Zhang, L., Yang, X., Tian, L., Chen, L., and Yu, W. (2016). Identification of peanut (Arachis hypogaea) chromosomes using a fluorescence in situ hybridization system reveals multiple hybridization events during tetraploid peanut formation. New Phytol. 211, 1424–1439. doi: 10.1111/nph.13999

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: citrus rootstock, satellite repeat, ribosomal DNA, oligonucleotide probe, cytogenetic markers, chromosome discrimination

Citation: Deng H, Tang G, Xu N, Gao Z, Lin L, Liang D, Xia H, Deng Q, Wang J, Cai Z, Liang G and Lv X (2020) 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. Front. Plant Sci. 11:569. doi: 10.3389/fpls.2020.00569

Received: 29 November 2019; Accepted: 16 April 2020;
Published: 27 May 2020.

Edited by:

Sergio J. Ochatt, INRA UMR 1347 Agroécologie, France

Reviewed by:

Wen-Wu Guo, Huazhong Agricultural University, China
Ki-Byung Lim, Kyungpook National University, South Korea

Copyright © 2020 Deng, Tang, Xu, Gao, Lin, Liang, Xia, Deng, Wang, Cai, Liang and Lv. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Guolu Liang, lianggl@swu.edu.cn; Xiulan Lv, xllvjj@163.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.