International Journal of
Genetics and Molecular Biology

  • Abbreviation: Int. J. Genet. Mol. Biol.
  • Language: English
  • ISSN: 2006-9863
  • DOI: 10.5897/IJGMB
  • Start Year: 2009
  • Published Articles: 139

Full Length Research Paper

The characterization of Citrus sp. from Parang Island Karimunjawa based on morphological, DNA barcoding and nutritional analysis

Hermin Pancasakti Kusumaningrum1*,
  • Hermin Pancasakti Kusumaningrum1*,
  • Department of Biology, Faculty of Sciences and Mathematic, Diponegoro University, Jl. Prof Sudharto SH. Tembalang Semarang 50275, Phone/Fax +62-02476480923, Indonesia.
  • Google Scholar
Anto Budiharjo
  • Anto Budiharjo
  • Department of Biology, Faculty of Sciences and Mathematic, Diponegoro University, Jl. Prof Sudharto SH. Tembalang Semarang 50275, Phone/Fax +62-02476480923, Indonesia.
  • Google Scholar
Agung Suprihadi
  • Agung Suprihadi
  • Department of Biology, Faculty of Sciences and Mathematic, Diponegoro University, Jl. Prof Sudharto SH. Tembalang Semarang 50275, Phone/Fax +62-02476480923, Indonesia.
  • Google Scholar
Yuriza Eshananda
  • Yuriza Eshananda
  • Department of Biology, Faculty of Sciences and Mathematic, Diponegoro University, Jl. Prof Sudharto SH. Tembalang Semarang 50275, Phone/Fax +62-02476480923, Indonesia.
  • Google Scholar
Annisa Fadillah
  • Annisa Fadillah
  • Department of Biology, Faculty of Sciences and Mathematic, Diponegoro University, Jl. Prof Sudharto SH. Tembalang Semarang 50275, Phone/Fax +62-02476480923, Indonesia.
  • Google Scholar
Dina Rahayuning Pangestuti
  • Dina Rahayuning Pangestuti
  • Faculty of Public Health, Diponegoro University, Jl. Prof Sudharto SH. Tembalang Semarang 50275, Indonesia.
  • Google Scholar


  •  Received: 21 June 2018
  •  Accepted: 04 September 2018
  •  Published: 30 September 2018

 ABSTRACT

The Citrus sp. from Parang Island Karimunjawa is a wild type of Citrus plant that grow on salt area and exhibits a high level of vitamin C.  Due to its ability, morphological, molecular and nutritional characterization needs to be carried out to improve its potential. The study was carried out at the Biotechnology Laboratory and the Integrated Laboratory of Diponegoro University from 2016 to 2017. The characterization was based on the morphplogical appearence of the tree, its fuits, and leaves. Its DNA barcoding consists of 18S ribosomal RNA and the ITS region on the plants is dispersed along the coastal ends and the centre of the island. Furthermore, the nutritional characterization consists of an edible fruit part, a high vitamin C level and protein content. The research showed that the Citrus sp. on the coastal ends and at the center of Parang Island Karimunjawa exhibits a close relation with the C. hystrix and members of Papeda clade. However, the morphology of the Citrus is quite different from the common C. hystrix in Indonesia. The nutritional content also shows that the Vitamin C content is higher than that of the C. hystrix by almost four times. The features of the Citrus sp. in Parang Island Karimunjawa have potential and showed possibility to improve its superiority in industrial applications and breeding programs.

 

Key words: Karimunjawa, Parang, Citrus, 18S rRNA, ITS, nutrition.


 INTRODUCTION

Citrus is one of the most popular world fruits. It also contains important nutritional elements for health.  Citrus is a good source of  vitamin  C  (ascorbic  acid),  phenolic compound, flavonoid, folic acid, potassium, pectin and antioxidant properties (Chiba et al., 2003; Abirami et al., 2014; Gosslau et al., 2014; Rafiq et al., 2016; Samraj and Rajamurgugan, 2017). As humans lack the ability to synthesize and store ascorbic acid, their daily requirement  depend on fresh fruits.and vegetables. Citrus fruit production in Indonesia reaches 2 million tons each year while its consumption in 2018 totaled almost 2.76 to 2.45 million tons. Commonly, it was planted in the highland and lowland (Direktorat Jenderal Holtikultura, 2015). Citrus plants in Indonesia that can grow in coastal areas or areas with high salt content are very rare, except Swing Orange [Limnocitrus littoralis (Mig) Swing] in the area of ​​Jepara  which is considered as endangered species and has the strength of disease resistance (World Conservation Monitoring Centre, 1998). We found the indigenous citrus plant that grows on salinity area on Parang Island Karimunjawa which is an archipelago located about 80 km to the north of the island of Java with a total area of about 78 km2 and the largest island covering an area of 2700 ha (Tomascir et al., 1997). Sea salinity in the western season is 32.6 ppm and 32.2 ppm in the eastern season. Its territory consists of several islands including Parang Island with an area of 690 ha. Their geographical position is at 5°42'-6°00'S, 110°07'-110°37'E, with air temperatures reaching 23 to 32°C and altitude of about 0 to 605 mabove sea level (Forestry Department, 2004). The wild Citrus in Indonesia is generally recognized by their morphology and nutrition content as conducted on Swing orange (Nuryandani, 2012; Adelina and Adelina, 2017). Characterization of some Citrus in Indonesia had been performed by molecular methods using isozymes, RAPD and ISSR related to their high genetic diversity and environmental adaptation (Fang et al., 1998; Karsinah et al., 2002; Agisimanto et al., 2007; Novelli et al., 2006; Bayer et al., 2009;  Morton,  2009;  Uzun  et  al.,  2009;  Penjor  et  al., 2014; Shimizu et al., 2016; Uchoi et al., 2017). Identification of  Citrus sp. from Parang Island has never been found. Molecular markers provide abundant information compared to morphological data; they are more efficient and are insensitive to environmental factors. Therefore, characterization based on a combination of  them has become important in the identification of this indigenous species. Until now, phylogenetic data for oranges that grow in salt area is very limited, especially oranges from Indonesia. Citrus sp. from Parang Island  is able to live in coastal area and it can be used as a commercial orange rootstock that is susceptible to salinity. Since the sea area of Indonesia is around 63%, salinity is a major agricultural problem in lands that decreases growth and productivity, it is also possible to develop varieties with enhanced salt tolerance that  can expand citrus cultivation into the salt affected marginal lands (Singh et al., 2003; Vijayan, 2009; Syvertsen et al., 2010; Hanin et al., 2013).


 MATERIALS AND METHODS

Plant material
 
Fresh leaves and fruits of Citrus sp. were collected from Parang Island Karimunjawa. Location, state and name of species are provided in Figure 1.
 
 
Determination of nutrition and ascorbic acid/vitamin C using iodometric titration method
 
The nutritional analysis and iodometric titration method was conducted according to another study with some modifications (Ciancaglini  et  al.,  2001;  AOAC,  2005;  Ywassaki and Canniattii- Brazaca, 2011; Spínola et al., 2013 ). Citrus samples was ground in a  mortar and pestle followed by adding distilled water several times while grinding the sample, each time decanting off the liquid extract into a 100 mL volumetric flask.  The solution was centrifuged to obtain  the filtrate. About 5 ml of filtrate was placed  in  the erlenmeyer and a solution of starch was added as much as 1% or about 2 ml. About 20 ml of distilled water was added into the place and the solution was titrated with 0.01 N iodine solution. The endpoint of the titration was identified as the first permanent trace of a dark blue-black colour due to the starch-iodine complex.
 
DNA extraction
 
Total genomic DNA was extracted from 0.5 to 1.0 g of fresh  leaf material. Genomic DNA of Citrus species was extracted through cetyl trimethyl ammonium bromide (CTAB) method (Doyle  and Doyle, 1987) with some modification. Quantification of DNA isolation product was conducted spectrophotometrically using Nanodrop.  The quality of DNA was determined by electrophoresis on 0.8% agarose gel.
 
DNA barcoding analysis
 
The DNA barcoding analysis was conducted using the 18S ribosomal RNA  (Region I) and sequences of  ITS1,5.8S ribosomal RNA gene, ITS2, sequence of 28S rRNA gene (Region II). Both regions are DNA barcoding commonly used for phylogenetic analysis by researchers. The 18S ribosomal RNA fragment was amplified using the primer pair of Primer F (5′-GTA GTC ATA TGC TTG TCT-3′) and Primer R (5′-GCT GGC ACC ASA CTT GCC CT-3′) (Kusumaningrum, 2008). The final PCR cocktail of 25 μL ml contained 50 ng of genomic DNA, 2.5 μL PCR buffer (KAPA), 10 mM concentration of deoxynucleoside triphosphate mix (KAPA), 2.5 pmol of each primer and 0.625 U of Taq Extra Hotstart DNA Polymerase and ddH2O. The amplifying reactions were run for 25 cycles for 3 min of pre-denaturation at 94°C, 25 s of denaturation at 94°C, 30 s of primer annealing at 55°C, 50s of elongation at 72°C, and 1 min of final elongation at 72°C. The PCR primers ITS 4 and ITS 5 (White et al., 2016) were used to amplify the ITS region (ITS 1,5.8S, and ITS 2) utilizing same primers for sequencing. The amplification program consisted of one cycle of initial denaturation at 94°C for 4 min followed by 25 cycles of 94°C for 1 min, 55°C for 3 min and 72°C for 1 min. This was followed by a 7 min extension at 72°C to allow completion of unfinished DNA strands, which in turn links to a soak file at 4°C. The PCR mixture of ITS primer contains 50 ng of genomic DNA, 2.5 μL PCR buffer, 10 mM concentration of deoxynucleoside triphosphate mix, 2.5 pmol of forward primer, 2.5 pmol of reverse primer and 0.625 U of Taq Extra Hotstart DNA Polymerase and Nuclease Free Water until reaching volume of 50 μL. DNA amplification was performed in a thermal cycler system. Amplified PCR products were purified using QIAquick gelextraction kit.
 
Phylogenetic  analysis
 
The sequence characteristics of the ribosomal RNA and ITS region were calculated using MEGA version 6 (Tamura et al., 2013). For data analysis, published sequences of some members from the genus Citrus were downloaded from GenBank. Sequence data of the Citrus 18S rRNA partial sequences and ITS are listed in Table 1. Juke-Cantor method (Jukes and Cantor, 1969)  was used to analyze the aligned sequence data. The phylogenetic tree was constructed using Phylip (Felsenstein, 2004). Bootstrap analysis was carried out with 999 random seed and 1000 replicates to examine the relative level of  support  for  individual  clades  on  the cladograms of each search. Genetic relationship was analyzed by phylogenetic tree construction.
 


 RESULTS

Morphological appearance of the Citrus sp. on Parang Island Karimunjawa
 
The Citrus sp. on Parang Island Karimunjawa is grown in the center and coastal areas of the land and based on our observation, the morphology of the plant exhibits some differences in some parts of the plant. As depicted in Figure 2, the fruit of Citrus plant grown on land area were big, globose, ovoid, very irregularly bumpy, glabrous with scattered  glandular  dots. The young green fruits become ripe and turns yellow with an average size of 7 to 10 cm in diameter. The peel is thick with its exterior layer of ±0.3 cm thickness and yellowish green. The  inner part was white and the pulp was yellowish green. The taste is very sour and slightly bitter with a faint fragrance. The fruiting pedicel was about 0.3 to 0.5 cm long. The seeds of the Citrus fruit are numerous, ovoid-oblong, ridged, 1.5 to 1.7cm long, 1 to 1.1 cm wide and 0.4 cm thick. 
 
The morphological characteristics of the Citrus fruit that grow on coastal areas were small, a little bumpy, globose, light green when ripe, feebly shining, with 5-7 cm diameter; thick peel, an exterior layer of ±0.2 cm thickness, light green, the inner part white; pulp yellowish green, very sour and slightly bitter, and releases a faint fragrance. The fruiting pedicel was about 0.3-0.5 cm long. The seeds of the Citrus fruit are numerous, ovoid-oblong, ridged, 1.2-1.5 cm long, 0.5-0.8 cm wide and 0.3 cm thick.
 
The characteristics of the  Citrus  leaf  which grows   on the land area as showed in Figure 3 were thin. The leaf surfaces was glabrous on both surfaces, light green adaxial, dull, light green or yellowish-green abaxial, densely pellucid dotted, fragrant when bruised, alternate, stalked unifoliate, broadly orbicular-ovate or ovate-oblong, lanceolate; base cuneate, obtuse or rounded, rarely subcordate, apex obtuse, rounded or slightly acuminate, often notched; patently serrate-crenate, coriaceous,  9.5-10 cm long and 3.5-4 cm wide.
 
 
The Citrus leaves that grows on the coastal area were also thin. The leaf surfaces were glabrous on both surfaces, light green adaxial, dull, light green or yellowish-green abaxial, dull, densely pellucid dotted, fragrant when bruised, alternate, stalked unifoliate, broadly orbicular-ovate or ovate-oblong, lanceolate; base cuneate, obtuse or rounded, rarely subcordate, apex obtuse, rounded or slightly acuminate, often notched; patently serrate-crenate, coriaceous,  9.5-10 cm long, and 3.5-4 cm wide. The petiole of Citrus leaves that grows on the coastal  areas were  almost  the  same  with
 
the plants found in the center of the island. It is long and expands into prominent wings, the winged part are obovate or obcordate-oblong, with an acute, cuneote, obtuse or rounded base and an obtuse, truncate, rounded or slightly emarginate apex, patently crenate-exsculptate, coriaceous, 0.3-0.5 cm above the base upwards with large, foliaceous wings, light green adaxial, shining, yellowish green abaxial,  the wings 1-8 cm long and 1-4.5 cm wide.
 
The Citrus trees that grows on the inland area are tall with a height of about 2-10 m (Figure 4). The trunk crooked with glabrous spiny branches, asymmetric or angular, thick, branched near the base; irregular crown, densely branched; branchlets rather thin. The tree branches were compressed-acutangular when young and become terete as they grow older. The branches are dark green,     glabrous     with     scattered     glandular    dots, accomplished with axillary spines; spines long, stiff, subulate, green with hard brown or orange-coloured tips, obliquely erect, solitary, glabrous, 1.5-2 cm long.
 
 
The height of the Citrus tree which grows on the coastal area is almost 15 m high. The trunk is also crooked with glabrous spiny branches, asymmetric or angular, thick, branched near the base; irregular crown, densely branched; branchlets are rather thin, when young compressed-acutangular, when older terete, dark green, glabrous with scattered glandular dots. The tree branch is coupled with axillary spines; very long spines, stiff, subulate, green with hard brown or orange-coloured tips, obliquely erect, solitary, glabrous, and 3.5-4 cm long. Based on the characteristics, the clear morphological differences between the coastal and land Citrus was the decrease of leaf thickness and greenness, folllowed by longer  and  bigger  stem  spines  on  the  coastal  plants.
 
Nutritional analysis of Citrus sp. of Parang Karimunjawa Island
 
The Citrus sp. of Parang Island Karimunjawa nutritional analysis displayed in Table 2 shows a remarkable amount of Vitamin C on the Citrus fruit which grows on coastal region about 66.37 mg/100 g BDD. This amount was almost two times higher than Citrus which was grown on the center of the island. The quantity of Citrus nutrition showed various content in the fruit. High level of ascorbic acid  in Citrus sp. of  Parang Island will increase resistance to a lot of biotic and abiotic stresses like high salt concentration (Davey et al., 2006).
 
The juice contains the total amount of soluble sugar (15.43 mg/100 ml) and free amino acids (15.8 mg/100 ml) as other species of Citrus, whereas C. limon contains the least amount of the same compounds that is 4.37 and 3.6 mg/100 ml, respectively. Another researcher reported the lower (Kumar et al., 2013) and almost similar findings with the Citrus sp. inland of  Parang  Island  Karimunjawa
(Mahmud et al., 2009).
 
DNA barcoding of Karimunjawa Island Citrus using fragment  of 18S rRNA gene, and ITS 1,5.8S ribosomal RNA gene, ITS2, fragment of 28S rRNA gene
 
The DNA barcoding to characterize the Citrus sp. on Parang Island from Karimunjawa archipelago consists of two regions, which are the 18S rRNA (region I) and  fragment of 18S rRNA gene, ITS 1,5.8S ribosomal RNA gene, ITS2 and fragment of 28S rRNA gene (region II).  Figure 5 shows the products of the DNA barcoding of the Citrus sp. in Parang Island.
 
The ITS of the 18S-26S nuclear ribosomal RNA separates the three gene region coding for the 18S, 5.8S, and 26S ribosomal subunits respectively. The ITS1 spacer is located between the 18S and 5.8S regions, the ITS2 spacer is between 5.8S and 26S. The annealing temperature for region I showed positive bands at 51 to 52°C and exhibited the best single DNA band. These data confirmed the primer we chose in this study was suitable for amplifying the conserved region of 18S ribosomal RNA in this Citrus. The size of the PCR product of 18S rRNA fragment was about 1000 bp. The best amplifying result of region II was achieved at the annealing temperature of 52°C. The size of the PCR product of ITS1, 5.8S rRNA gene, ITS2 region were about 750 bp.
 
The identity of our sequencing results using 18S rRNA is high and the value ranged from 92 to 98% compared to the existing sequence sources of existent Citrus species in GenBank database. This result suggested that the ITS universal primers has been successfully applied for the genus Citrus plants; the nrDNA ITS region could be successfully amplified using ITS universal primer sets. This result also in accordance with others (Sun et al., 2015). The sequences of region I showed homologies about 98 to 99% with 18S rRNA partial sequences of 30 Citrus species in the GenBank. In this study, as exhibited in Figure 6 and phylogenetic tree on Figure 7, the two Citrus sp. of Parang  Island  displayed  high similarity with C. trifoliata (99%) with accession number KJ740219.1. C. tachibana. Phylogenetic tree exhibited most of the features that is a species of Citrus that originated from China. A similarity analysis shows that there are 2 deletion bases and 2 subtitution bases between 18S rRNA partial sequences of the Citrus sp. on Parang Island with C. trifoliata marked with a red box.
 
 
 
 
The deletion and subtitution bases potentially indicated the nature itself for the Citrus sp. of Parang Island since 18S rRNA was the region with conserved sequences. The differences between the coastal and inland Citrus sp. of Parang Island was displayed by base substitution which exhibited a replacement of  cytosine (C) on the Citrus sp. Parang Island land with thymine (T) on the same plant which lives on the coast. This subtitution was a turnover among the same base pairs (transitions) of pyrimidine. We also found a guanine (G) base insertion on the inland Citrus sp. of Parang Island comparing that of the coastal region.
 
Different results were reflected among the 18S rRNA partial sequences of the Citrus sp. on Parang  Island Karimunjawa and their ITS region. These regions gained high similarity (97%) with C. hystrix as exhibited in Figure 8 and the phylogenetic tree in Figure 9. This result also is supported with the similarity on little bumpy morphology in their fruit skin as shown in Figure 1. The results study about high homology between Citrus sp. Parang Island with C. hystrix was not correlated with their vitamin C content. The  higher vitamin C content  of Citrus sp.  Parang Island than C. hystrix obtained from this study exhibited the unique and specific character of this native plant. Furthermore, the sequences in the ITS1-ITS2 region also exhibited higher differences among the coastal and inland Citrus sp. of Parang Island Karimunjawa. In contrast with 18S rRNA partial sequences, the sequences of ITS1-ITS2 region showed lower homologies of about 89-95% with other Citrus from Vietnam, India, Japan, and China. This result was in agreement with the similarity analysis in showing variable sequences at the end of 18S rRNA sequences. The ITS1 sequences showed lowest differences compared to the 5.8S rRNA and  ITS2. The  several subtitution in the ITS2 was considered to be the specific characteristic of the Citrus sp. on the coastal area of Parang Island Karimunjawa because the land area Citrus sp. of Parang Island and C. hystrix did not show this subtitution in this region. The reasons for the difference in the homology obtained from phylogenetic analysis and ITS region sequencing for the both Citrus sp. from Parang Island which live inland and coastal is having the same anchestor. Nei and Kumar (2000) stated that when two DNA sequences come from the same ancestral sequence, the sequence of descendants will gradually be differentiated by nucleotide substitution. The variation in sequences and differences of ITS2 among the coastal and inland plants has been shown to be valuable in identifying both of them. This study is also in good agreement with the previous report (Alvarez and Wendel, 2003; Sun et al., 2015). Amongst the two Citrus sp. on Parang Island Karimunjawa investigated in this study, there were some found variations in nucleotide substitution, deletion, or insertion. Number of base substitution found between them were 8 transitions and 14 transversions.
 
The probability level of transversion is greater than transition which indicates that the changes in the ITS region can potentially alter the genes to ensure its position in classification at Citrus group. The phylogenetic and genetic distance analysis using 18S rRNA showed that Citrus sp. of Parang Island was in the same branch with C. tachibana. The ITS region showed different result in gained close genetic relationship with C. montana, C. hystrix and C. macroptera.
 
 
 

 


 DISCUSSION

In this study, Citrus sp. ws characterized which grow on  Parang Island Karimunjawa based on morphogical, nutritional analysis, and DNA barcoding. In some studies, morphological characterization fails to obtain the representative identity of the Citrus. The study using morphologic  appearence   proved   these  characteristics were independent and related to diverse pressures and evolutionary factors and require other methods to complete the assessment. Therefore, it still raises a lot of questions in attempt to further improve breeding programs  (Koehler-Santos  et  al.,  2003;  Yahata  et  al.,
2017). Morphological characterization of the Citrus sp. on Parang Island Karimunjawa also showed some differences in some parts of the plant compared with other Citrus plants in Indonesia. These characteristics seems quite specific for the Citrus on Parang Island Karimunjawa. The fruit size and weight on the coastal region tends to be smaller but higher in Vitamin C content. It have greener leaf, narrower, bigger and stronger stem pine. Presumably, these character was related to salinity susceptibility and is essential for plant survival and maintained growth rate. This characteristics were also verified from other studies (Moya et al., 1999; Murkute et al., 2005; Syvertsen et al., 2010; Hussain et al., 2012; Acosta-Motos et al., 2017). Some study explain that the effects of high salt environment had induce plant adaptation and rise the level of ascorbic acid in several plants (Hernandez et al., 2000; Davey et al., 2006; Acosta-Motos et al., 2017). These studies supported the nutritional analysis of the Citrus sp. from Parang Island which showed that its Vitamin C content is higher than C. hystrix almost by four times, although they have close genetic relationship and similarity in morphology (Murkute et al., 2005).  Vitamin C is one of the known forms of ascorbic acid which is a water-soluble chemical in fruits (Kumar et al., 2013). This substance amount has been shown to be valuable since vitamin C is one of important components of citrus. Another study reported that it contributes to the antioxidant activity about 56 to 77% from Citrus extract (Kumar et al., 2013; Abirami et al., 2014). In higher plants, ascorbic acid  biosynthesis is from glucose using L-galactose pathway. Overexpression of several structural ascorbic acid-related genes from various ascorbic acid metabolic pathways is not succesfull in most species in order to enrich ascorbic acid levels (Mellidou and Kannelis, 2017). The role of two genes of Guanosin diphosphate (GDP) in the L-galactose pathway, GDP-D-mannose pyrophosphorylase (VTC1 or GMP) and GDP-D-mannose-3,5-epimerase (GME), has been correlated with ascorbic acid concentrations in fruit of  some species such as kiwifruit (Bulley et al., 2009), apple (Li et al., 2010) and blueberry (Liu et al., 2015). It is also related with ascorbic acid accumulation under salt stress (Zhang et al., 2011). Several studies on genetic factors indicated that ascorbic acid accumulation showed influence of heritability (Davey et al., 2006). The high concentration of absisic acid in the Citrus sp. from Costal area of Parang Island offers a potency to be exploited to enhance it accumulation in fruit due to its importance as bioactive nutrients and dietary antioxidants.
 
The DNA barcoding approach involving the 18S rRNA region evolves relatively slowly compared to the ITS region due to their high conservation. Despite their weakness, the use of ribosomal RNA as an alignment tool and similarity analysis in plant phylogenetic studies still offers advantages because it is not influenced by environmental changes compared with other methods such as SRAP, matK, ISSR (Rogers and  Bendich,  1987; Fang et al., 1998; Johnson et al., 1999; Uzun et al., 2009; Penjor et al., 2014). The rapidly evolving ITS spacer sequences have been used extensively in phylogenetic studies due to its application not only in lower levels, but also help resolve intra-family relationships (Rogers and Bendich, 1987, Johnson et al., 1999; Kyndt et al., 2010). This was exhibited by the inconsistent results of similarity analysis which showed that the Citrus sp. on Parang Island Karimunjawa was the closest to C. aurantium but the phylogenetic tree of 18S rRNA shows its relationship with C. tachibana. Interestingly, Hirai et  al. (1990) and Yamaji et al. (2013) found that C. tachibana was known as wild species of Citrus with intraspecific nuclear rDNA (nrDNA) variation in ITS that was thought to be fixed into a single ribotype. The study with ITS region found that based on the number of base substitution, deletion and insertion between the coastal and inland Citrus sp. on Parang island, it was clustered together with C. montana, C. hystrix and C. macroptera. This study of ITS was supported with other in showing that C. hystrix, C. macroptera, and C. montana had belong to one clade Papeda respectively (Li et al., 2010, Hynniewta et al., 2014, Yamaji et al., 2013). Papeda is a common name for a group of native tropical Asian citrus. Some species of Papeda have been used as genomic sources for breeding disease-resistance (Wang et al., 2017). Based on molecular studies, papeda is one of the ancestors  of many types of commercial limes (Xu et al., 2013; Wang et al., 2017; Wu et al., 2018). The morphological, nutritional and DNA barcoding analysis showed the possibility that Citrus sp. on Parang Island Karimunjawa was considered to have the specific characters and belongs to the Papeda clade. The ITS analysis showed  that the longer branch of Citrus sp. which grows in coastal area suggests a possibility of a more advanced evolution. This study showed that Citrus sp. from Parang Island with spesific emphasis on the ones that grow on the coast has genetic potentials and a remarkable vitamin C content and can live on salinity of 32.2-32.6 ppm. This potential will increase the genetic resources of Citrus and increase the possibility to elevate their superiority for breeding programs and industrial applications. Previous studies have shown different levels of ascorbic acid in fruits during their growth in saline soil (Bulley et al., 2009; Li et al., 2010; Liu et al., 2015). It raised open question whether the character of Citrus sp. Parang Island is inherited from the parent or is the result of adaptation to the environment. Further experiment will deal with salt treatment to the Citrus on the laboratory and the examinination on the effect of the EC or TDS of water and the soil on Citrus sp. Parang Island vitamin C content.

 


 CONCLUSION

Characterization of Citrus sp. Pulau Parang Karimunjawa based on morphology, DNA barcoding and nutritional analysis shows the closest relationship with C. hystrix and members of the Citrus species in Papepa clade. However, some variations in the ITS region sequence and vitamin C content indicate the specific character of Citrus sp. of Parang Island Karimunjawa, especially in the coastal area.


 ACKNOWLEDGEMENT

This research was funding by DIPA Fakultas Sains dan Matematika Universitas Diponegoro Year 2017 according to the Letter of Agreement   Surat Perjanjian Tugas Pelaksanaan Penelitian Para Dosen Departemen Biologi Fakultas Sains dan Matematika Universitas Diponegoro Number: 1643G/UN7.5.8/PP/2017 date 3 April 2017 which is gratefully acknowledged.


 CONFLICT OF INTERESTS

The authors declare that they have no conflict of interest.

 



 REFERENCES

Abirami A, Nagarani G, Siddhuraju P (2014). In vitro antioxidant, anti-diabetic, cholinesterase and tyrosinase inhibitory potential of fresh juice from Citrus hystrix and C. maxima fruits. Food Science and Human Wellness 3:6-25.
Crossref

 

Acosta-Motos J, Ortu-o MF, Bernal-Vicente A, Diaz-Vivancos P, Sanchez-Blanco M, Hernandez JA (2017). Plant responses to salt stress: adaptive mechanisms. Agronomy Journal 7(18):1-38.

 
 

Adelina SO, Adelina E (2017). Identifikasi morfologi dan anatomi jeruk lokal (Citrus sp) di desa Doda dan desa Lempe Kecamatan Lore. Jurnal Agrotekbis (In Indonesian) 5(1):58-65.

 
 

Agisimanto D, Martasari C, Supriyanto A (2007). Perbedaan primer RAPD dan ISSR dalam identifikasi hubungan kekerabatan genetik jeruk siam (Citrus suhuniensis L. Tan) Indonesia. Jurnal Hortikultura (In Indonesian) 17(1):101-110.

 
 

Álvarez I, Wendel JF (2003). Ribosomal ITS sequences and plant phylogenetic inference. Molecular Phylogenetics and Evolution 29:417-434.
Crossref

 
 

AOAC (Association of Official Analitycal Chemists) (2005). Official method of analysis (18th Edition). Washington.DC: Association of Officiating Analytical Chemists.

 
 

Bayer RJ, Mabberley DJ, Morton C, Miller CH, Sharma IK, Pfeil BE, Rich S, Hitchcock R, Sykes S (2009). A molecular phylogeny of the orange subfamily (Rutaceae: Aurantioideae) using nine cpDNA sequences. American Journal of Botany 96(3):668-685.
Crossref

 
 

Bulley SM, Rassam M, Hoser D, Otto W, Schünemann N, Wright M, MacRae E, Gleave A, Laing W. (2009). Gene expression studies in kiwifruit and gene over-expression in Arabidopsis indicates that GDP-L-galactose guanyltransferase is a major control point of vitamin C biosynthesis. Journal of Experimental Botany 60(3):765-778.
Crossref

 
 

Chiba H, Uehara M, Wu J, Wang X, Masuyama R, Suzuki K, Kanazawa K, Ishimi, Y (2003). Hesperidin, a citrus flavonoid, inhibits bone loss and decreases serum and hepatic lipids in ovariectomized mice. Journal of Nutrition 133(6):1892-1897.
Crossref

 
 

Ciancaglini P, Santos HL, Daghastanli KRP, Thedei G (2001). Using a classical method of vitamin C quantification as a tool for discussion of its role in the body. Biochemistry and Molecular Biology Education 29(3):110-114.

 
 

Davey MW, Kenis K, Keulemans J (2006) Genetic control of fruit vitamin C contents. Plant Physiology 142(1):343-351.
Crossref

 
 

Direktorat Jenderal Holtikultura (2015). Sub Sektor Holtikultura,  (In Indonesian).

 
 

Doyle JJ, Doyle JL (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical Bulletin 19:11-15.

 
 

Fang D, Krueger RR, Roose ML (1998). Phylogenetic relationships among selected Citrus germplasm accessions revealed by inter-simple sequence repeat (ISSR) markers. Journal of the American Society for Horticultural Science 123(4):612-617.

 
 

Felsenstein J (2004). PHYLIP (Phylogeny Inference Package) Version 3.6. (Distributed by the author) Department of Genome Sciences, University of Washington, Seattle.

 
 

Forestry Department (2004). Zonatıon of Karımunjawa natıonal park. Jepara Regency. Central Java Provınce.

 
 

Gosslau A, Chen KY, Ho CT, Li S (2014). Anti-inflammatory effects of characterized orange peel extracts enriched with bioactive polymethoxyflavones. Food Science and Human Wellness 3:26-35.
Crossref

 
 

Hanin M, Ebel C, Ngom M, Laplaze L, Masmoudi K (2016). New insights on plant salt tolerance mechanisms and their potential use for breeding. Frontiers in Plant Science pp. 1-17.
Crossref

 
 

Hernandez JA, Jimenez A, Mullineaux P, Sevilla F (2000). Tolerance of Pea (Pisum sativum L.) to long-term salt stree is associated with antioxidant defences. Plant, Cell and Environment 23(8):853-862.
Crossref

 
 

Hirai M, Mitsue S, Kita K, Kajiura I (1990). A survey and isozyme analysis of wild mandarin, Tachibana (Citrus tachibana (Mak.) Tanaka) growing in Japan. Japanese Society for Horticultural Science 59:1-7.
Crossref

 
 

Hussain S, Luro F, Costantino G, Ollitrault P, Morillon R (2012). Physiological analysis of salt stress behaviour of Citrus species and genera: Low chloride accumulation as an indicator of salt tolerance. South African Journal of Botany 81:103-112.
Crossref

 
 

Hynniewta M, Malik SK, Rao SR (2014). Genetic diversity and phylogenetic analysis of Citrus (L) from north-east India as revealed by meiosis, and molecular analysis of internal transcribed spacer region of rDNA. Meta Gene 2:237-251.
Crossref

 
 

Johnson LA, Soltis DE, Soltis PS (1999). Phylogenetic relationships of Polemoniaceae inferred from 18S ribosomal DNA sequences. Plant Systematics and Evolution 214:65-89.
Crossref

 
 

Jukes TH, Cantor CR (1969). Evolution of protein molecules. In Munro HN, editor, Mammalian Protein Metabolism, Academic Press, New York 3(21):132.

 
 

Karsinah, Sudarsono, Setyobudi L, Aswidinnoor H (2002). Keragaman genetik plasma nutfah jeruk berdasarkan analisis penanda RAPD. Jurnal Bioteknologi Pertanian (In Indonesian) 7(1):8-16.

 
 

Koehler-Santos P, Dornelles ALC, De Freitas LB (2003). Characterization of mandarin citrus germplasm from Southern Brazil by morphological and molecular analyses. Pesquisa Agropecuária Brasileira 38(7):797-806.
Crossref

 
 

Kumar R, Vijay S, Khan N (2013). Comparative nutritional analysis and antioxidant activity of fruit juices of some Citrus spp. Octa Journal of Biosciences 1(1):44-53.

 
 

Kusumanıngrum HP (2008). Karakterisasi alga hijau Dunaliella sp. dan isolat sianobakteria serta deteksi gen DXS penyandi enzim kunci biosintesis karotenoid [Characterization of green microalgae Dunaliella sp. and cyanobacteria isolate and detection of DXS gene encoding key enzyme of carotenoid biosynthetic] [PhD Dissertation]. Yogyakarta: Gadjah Mada Universitas (In Indonesian).

 
 

Kyndt T, Dung TN, Goetghebeur P, Toan HT, Gheysen G (2010). Analysis of ITS of the rDNA to infer phylogenetic relationships among Vietnamese Citrus accessions. Genetic Resources and Crop Evolution 57(2):183-192.
Crossref

 
 

Li X, Xie R, Lu Z, Zhou Z (2010). The origin of cultivated Citrus as Inferred from internal transcribed spacer and chloroplast DNA sequence and amplified fragment length polymorphism fingerprints. Journal of The American Society for Horticultural Science 135:341-350.

 
 

Liu F, Wang L, Gu L, Zhao W, Su H, and Cheng X. (2015). Higher transcription levels in ascorbic acid biosynthetic and recycling genes were associated with higher ascorbic acid accumulation in blueberry. Food Chemistry 188:399-405.
Crossref

 
 

Mahmud MK, Hermana N, Zulfianto A, Rozanna R, Apriyantono NI, Hartati B, Bernardus T (2009). Table of Indonesian food composition. Jakarta. Indonesian Nutritionist Association (In Indonesian).

 
 

Mellidou I, Kanellis AK (2017). Genetic control of ascorbic acid biosynthesis and recycling in horticultural crops. Frontiers in Chemistry Mini Review 5(40):1-8.
Crossref

 
 

Morton CM (2009). Phylogenetic relationships of the Aurantioideae (Rutaceae) based on the nuclear ribosomal DNA ITS region and three noncoding chloroplast DNA regions, atpB-rbcL spacer, rps16, and trnL-trnF. Organisms Diversity and Evolution 9(1):52-68.
Crossref

 
 

Moya JL, Primo-Millo E, Talon M (1999). Morphological factors determining salt tolerance in Citrus seedlings: The shoot to root ratio modulates passive root uptake of chloride ions and their accumulation in leaves. Plant Cell and Environment 22:1425-1433.
Crossref

 
 

Murkute A, Sharma S, Singh SK (2005). Citrus in terms of soil and water salinity: A review. Journal of Scientific and Industrial Research 64(06):393-402.

 
 

Nei M, Kumar S (2000) Molecular evolution and Phylogenetics. Oxford University Press pp. 33-46

 
 

Novelli VM, Cristofani M, Souza A, Machado M (2006). Development and characterization of polymorphic microsatellite markers for the sweet orange (Citrus sinensis L. Osbeck). Genetics and Molecular Biology 29(1):90-96.
Crossref

 
 

Nuryandani, E (2012). Persebaran dan karakterisasi induk jeruk keprok Tawangmangu asli (Citrus reticulata Blanco ssp Tawangmangu) (In Indonesian) 13(1):33-42.

 
 

Penjor T, Mimura T, Matsumoto R, Yamamoto M, Nagano Y (2014). Characterization of limes (Citrus aurantifolia) grown in Bhutan and Indonesia using high-throughput sequencing. Scientific Reports 4:1-9.

 
 

Rafiq S, Kaul R, Sofi SA, Bashir N, Nazir F, Nayik GA (2016). Citrus peel as a source of functional ingredient: A review. Journal of The Saudi Society of Agricultural Sciences 7:1-8.
Crossref

 
 

Rogers SO, Bendich AJ (1987). Ribosomal RNA genes in plants: variability in copy number and in the intergenic spacer. Plant Molecular Biology 9(5):509-520.
Crossref

 
 

Samraj S, Rajamurgugan S (2017). Qualitative and quantitative estimation of bioactive compounds and antioxidant activity in Citrus hystrix. International Journal of Engineering Science and Computing 7(6):13154-13163.

 
 

Shimizu T, Kitajima A, Nonaka K, Yoshioka T, Ohta S, Goto S, Toyoda A, Fujiyama A, Mochizuki T, Nagasaki H, Kaminuma E, Nakamura Y (2016). Hybrid origins of citrus varieties inferred from DNA marker analysis of nuclear and organelle genomes. PLoS ONE 11(11): 1-58.
Crossref

 
 

Singh A, Saini ML, Behl RK (2003). Screening of citrus rootstocks for salt tolerance in semi-arid climates-A review. Tropics 13:53-66.
Crossref

 
 

Spínola V, Mendes B, Câmara JS, Castilho PC (2013). Effect of time and temperature on vitamin C stability in horticultural extracts. UHPLC-PDA vs iodometric titration as analytical methods. LWT - Food Science and Technology 50:489-495.
Crossref

 
 

Sun YL, Kang HM, Han SH, Park YC, Hong SK (2015). Taxonomy and phylogeny of the genus Citrus based on the nuclear ribosomal and its region sequence. Pakistan Journal of Botany 47(1):95-101.

 
 

Syvertsen JP, Melgar JC, García-Sánchez F (2010). Salinity tolerance and leaf water use efficiency in Citrus. Journal of the American Society for Horticultural Science 135(1):33-39.

 
 

Tamura K, Glen S, Daniel P, Alan F, Sudhir K (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Molecular Biology and Evolution 30(12):2725-2729.
Crossref

 
 

Tomascir T, Mah AJ, Nontji A, Moosa MK (1970). The ecology of the Indonesian Seas, Part Two. Hong Kong: Eric Oey, Periplus Editions Ltd pp. 685-686.

 
 

Uchoi A, Malik SK, Choudhary R, Kumar S, Pal D, Rohini MR, Chaudhury R (2017). Molecular markers in assessing genetic variation of Indian citron (Citrus medica L.) cultivars collected from different parts of India. Indian Journal of Biotechnology 16:346-356.

 
 

Uzun A, Yesiloglu T, Aka-Kacar Y, Tuzcu O, Gulsen O (2009). Genetic diversity and relationships within Citrus and related genera based on sequence related amplified polymorphism markers (SRAPs). Scientia Horticulturae 121(3):306-312.
Crossref

 
 

Vijayan K (2009). Approaches for enhancing salt tolerance in mulberry (Morus L) -A review. Review Article. Plant Omics Journal 2(1):41-59.

 
 

Wang X, Xu Y, Zhang S, Cao L, Huang Y, Cheng J, Wu G, Tian S, Chen C, Liu Y, Yu H, Yang X, Lan H, Wang N, Wang L, Xu J, Jiang X, Xie Z, Tan M, Larkin RM, Chen L-L, Ma B-G, Ruan Y, Deng X, Xu X (2017). Genomic analyses of primitive, wild and cultivated citrus provide insights into asexual reproduction. Nature 49(5):765-775.

 
 

White TJ, Burns T, Lee S, Taylor J (1990). Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. M.A. Innis. D.H. Gelfand, J.J. Sninsky, T.J. White (Eds.), PCR Protocols: A Guide to Methods and Amplifications, Academic Press, San Diego & London pp. 315-322.

 
 

World Conservation Monitoring Centre (1998). Limnocitrus littoralis. The IUCN Red List of Threatened Species 1998: e.T37412A10051534. 
Crossref

 
 

Wu GA, Terol J, Ibanez V, Lopez-Garcia A, Perez-Roman E Borredá C, Domingo CR, Tadeo FR, Carbonell-Caballero J, Alonso R, Curk F, Du D, Ollitrault P, Roose ML, Dopazo J, Gmitter Jr FG, Rokhsar DS, Talon M (2018). Genomics of the origin and evolution of Citrus. Nature 554:311-316.
Crossref

 
 

Xu Q, Chen L, Ruan X, Chen D, Zhu A, Chen C, Bertrand D, Jiao W, Hao B, Lyon MP, Chen J, Gao S, Xing F, Lan H, Chang J, Ge X, Lei Y, Hu Q, Miao Y, Wang L, Xiao S, Biswas MK, Zeng W, Guo F, Cao H, Yang X, Xu X, Cheng Y, Xu J, Liu J, Luo OJ, Tang Z, Guo W, Kuang H, Zhang H, Roose ML, Nagarajan N, Deng X, Ruan Y (2013). The draft genome of sweet orange (Citrus sinensis). Nature Genetics 45(1):59-66.
Crossref

 
 

Yahata M, Kunitake H, Komatsu H (2017). Morphological characterization and evaluation of reproductive function in a haploid pummelo [Citrus maxima (Burm.) Merr.]. Japan Agricultural Research Quarterly 51(4):293-298.
Crossref

 
 

Yamaji H, Kondo K, Kuniga T, Nesumi H, Yoshida T, Hashimoto K, Takheda O (2013). Origin of cultivated Citrus (Rutaceae) documented by the contents of internal transcribed spacer sequences (ITS) in nuclear ribosomal DNA. Journal of Japanese Botany 88:222-238.

 
 

Ywassaki LA, Canniatti-Brazaca SG (2011). Ascorbic acid and pectin in different sizes and parts of citric fruits. Ciência e Tecnologia de Alimentos 31(2):319-326.
Crossref

 
 

Zhang C, Liu J, Zhang Y, Cai X, Gong P, Zhang J, Wang T, Li H, Ye Z. (2011). Overexpression of SlGMEs leads to ascorbate accumulation with enhanced oxidative stress, cold and salt tolerance in tomato. Plant Cell Report 30(3):389-398.
Crossref

 

 




          */?>