- 1Department of Biology, University of Indonesia, Depok, Indonesia
- 2Research Center for Applied Botany, National Research and Innovation Agency, Bogor, Indonesia
- 3Research Center for Biota Systems, National Research and Innovation Agency, Bogor, Indonesia
- 4Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
- 5Research Center for Biosystematics and Evolution, National Research and Innovation Agency, Bogor, Indonesia
Introduction
Dehaasia pugerensis Koord. & Valeton is an endemic Indonesian species restricted to Jember Regency, East Java, particularly in the Gunung Watangan area (Bijmoer et al., 2020). This species belongs to the Lauraceae family and has been classified as Critically Endangered (CR) on the IUCN Red List based on criteria B1ab (i, ii, iii, iv, v) and C2a (i) (Helmanto et al., 2022). Its distribution is highly limited, occurring at elevations between 61 and 391 meters above sea level, and its population is threatened by habitat degradation and overexploitation, particularly for charcoal production (Helmanto et al., 2022). These anthropogenic pressures have led to a significant population decline, highlighting the urgent need for science-based conservation efforts, including genomic studies as a foundation for sustainable management of rare species.
One of the key approaches in plant genomics is the characterization of the chloroplast genome. The chloroplast genome in land plants typically exhibits a conserved circular structure, consisting of four major regions: the large single-copy (LSC), the small single-copy (SSC), and two inverted repeat (IR) regions. The total chloroplast genome size in land plants ranges from approximately 120 to 200 kb, with IR regions typically spanning 20–26 kb (Xiao-Ming et al., 2017). In recent decades, advances in sequencing technologies have enabled faster and more accurate assembly of chloroplast genomes. These developments have opened new opportunities to explore the structure, variation, and evolution of chloroplast genomes across plant species, contributing to phylogenetic studies and the conservation of rare genetic resources (An et al., 2022).
To date, the complete chloroplast genome of D. pugerensis has not been reported. The absence of such genomic information limits molecular-based conservation initiatives for this species. Within the Lauraceae family, chloroplast genome sizes exhibit considerable variation, as recorded in Alseodaphne semecarpifolia (153,051 bp), Eusideroxylon zwageri (157,535–157,577 bp), and Neocinnamomum spp. (150,753–150,956 bp) (Song et al., 2017; Chen et al., 2017; Li et al., 2017; Cao et al., 2023; Zhu et al., 2023). These data underscore the importance of generating plastome sequences from rare species to deepen our understanding of evolutionary patterns and genetic diversity within Lauraceae. Therefore, this study aims to characterize and reconstruct the complete chloroplast genome of D. pugerensis as a foundational resource for genomics-based conservation.
Methods
Plant materials
Fresh young leaves were collected from a healthy D. pugerensis seedling approximately 1 meter tall, originally sourced from its natural habitat in the Perhutani Forest, Puger, East Java (8°24’10.58” S, 113°30’33.74” E; accession number P.3.1.41). The research and sampling activities were conducted under formal permission from the Perhutani Forestry Institute (PeFI), as stated in Letter No. 0392/001.6/PeFI/2024. The seedling has been conserved at the Cibinong Botanic Gardens (CBG) for ex situ conservation purposes.
DNA extraction, library preparation, and next-generation sequencing
Genomic DNA was extracted from plant leaves using the cetyltrimethylammonium bromide (CTAB) method as described by Doyle and Doyle (1987). The initial DNA concentration and purity were assessed using a Nanodrop 2000 spectrophotometer (Thermo Scientific, MA, USA). DNA integrity was visualized by agarose gel electrophoresis, and quantification was performed with the Qubit dsDNA HS Assay Kit (Thermo Scientific, MA, USA). Further evaluation of DNA integrity was conducted using the 4150 TapeStation system (Agilent Technologies, CA, USA). High-quality genomic DNA was then utilized for library preparation. Sequencing was carried out on the Illumina NextSeq 2000 platform (Genetika Science Lab, Tangerang, Indonesia) employing a paired-end 150 bp strategy, targeting a total data output of 10 Gb.
Chloroplast genome assembly, annotation and analysis
The reads were assessed for quality using FASTQC software version 0.11.8 (Andrews, 2010). Filtering and trimming were performed using Trimmomatic version 0.39 to remove low-quality bases (less than 30), adapters, nucleotide position bias at the 3’ and 5’ ends, and sequence contamination. The parameters applied included TruSeq3-PE.fa:2:30:10, SLIDINGWINDOW:4:28, LEADING:28, TRAILING:28, and MINLEN:20 (Bolger et al., 2014). The trimmed reads results were subsequently assembled utilizing GetOrganelle version 1.7.7.1 (Jin et al., 2020). The annotation of complete chloroplast genome of D. pugerensis was conducted utilizing CPGAVAS2 (http://47.96.249.172:16019/analyzer/annotate) (Shi et al., 2019; Lestari et al., 2024), with the cp genome of D. hainanensis (accession number: OP374101.1) serving as the reference. The annotation process faced challenges such as resolving ambiguous gene regions and validating gene boundaries, which were addressed via manual verification using Unipro UGENE v. 45.1 (Okonechnikov et al., 2012) and NCBI Genomic Workbench v. 3.8.2 (Kuznetsov and Bollin, 2021). To ensure the cp genome sequence contained no N bases and had 21 amino acids, Unipro UGENEv. 45.1 was employed. Genes without a start codon were manually edited using the edit menu in NCBI Genomic Workbench v. 3.8.2. The circular genome visualization was performed using Organellar Genome DRAW (OGDRAW) accessed through MPI-MP Chlorobox (Greiner et al., 2019).
Characterization of simple sequence repeats
Simple sequence repeats (SSRs) of D. pugerensis chloroplas genome were identified using the MIcroSAtellite (MISA) web tool (Beier et al., 2017). Search parameters were configured to detect perfect mono-, di-, tri-, tetra-, penta-, and hexa-nucleotide motifs, with minimum repeat thresholds of 10, 5, 3, 3, 3, and 3, respectively. Compound SSRs were allowed when two adjacent repeat motifs were separated by no more than 100 bp.
Codon usage analysis
Codon usage patterns and Relative Synonymous Codon Usage (RSCU) values were analyzed using MEGA X software (Kumar et al., 2018). Visualization of codon frequency distributions was subsequently performed using the “ggpubr” package in R version 4.2.3.
Results
The complete chloroplast genome of D. pugerensis spans 153,111 bp, exhibiting the typical quadripartite structure of angiosperms (Figure 1A), comprising a large single-copy (LSC) region of 93,852 bp, a small single-copy (SSC) region of 18,699 bp, and two inverted repeats (IRs) of 20,280 bp each. The GC content is 39.07%, consistent with other Lauraceae species. Genome annotation identified 128 functional genes, including 85 protein-coding genes, 8 rRNA genes, and 36 tRNA genes (Table 1). Seven genes are duplicated in the IRs, and 16 genes contain introns, reflecting a level of structural and regulatory complexity typical of Lauraceae plastomes.
Figure 1. Chloroplast genome features of Dehaasia pugerensis. (A) Circular chloroplast genome map generated with OGDRAW. Genes on the inner circle are transcribed counterclockwise, while those on the outer circle are transcribed clockwise. Functional gene groups are color-coded, and GC/AT content is shown in grey and light grey. (B) Diagram of cis-splicing genes, showing genes that contain introns within a single, continuous locus. (C) Diagram of trans-splicing genes, showing genes whose exons are separated across different genomic regions and joined through post-transcriptional splicing. (D) Distribution of SSR motif types identified in the chloroplast genome. (E) Codon usage patterns for amino acids based on all protein-coding genes in the chloroplast genome.
The conserved structure of the plastome provides a robust molecular framework for phylogenetic resolution within Lauraceae, where morphological convergence often obscures species boundaries (Liu et al., 2017; Tian et al., 2021). Duplicated genes and intron-rich regions offer potential molecular markers for evaluating genetic diversity, population structure, and gene flow (Han et al., 2022). This is critical for D. pugerensis, a narrowly endemic and critically endangered species, as chloroplast genomic data inform historical biogeography, demographic shifts, and conservation prioritization (Crawford and Stuessy, 2016).
The chloroplast genome of D. pugerensis harbors 11 genes with cis-splicing introns, where exons and introns reside on the same transcript (Figure 1B). These include rps16, atpF, rpoC1, ycf3, clpP, petB, petD, rpl2, ndhA, and ndhB (the latter two located in IRs and thus duplicated). Gene structures, including exons (black) and introns (white), are annotated with genomic coordinates in the corresponding figure. Most genes, such as atpF and rpoC1, contain one intron, while ycf3 and clpP have two, indicating complex splicing regulation. Genes located in the IR regions, such as rpl2 and ndhB, appear twice in the genome.
The presence of introns in genes like clpP and ycf3 is conserved among angiosperms and may reflect regulatory or evolutionary functions (Rogalski et al., 2015). The cis-splicing profile in D. pugerensis mirrors that of other Lauraceae species, indicating a conserved regulatory mechanism across the family (Song et al., 2017).
The chloroplast genome of D. pugerensis contains a single trans-splicing gene, rps12 (Figure 1C). Unlike cis-splicing, trans-splicing joins exons located in separate genomic regions. In D. pugerensis, rps12 is split into three exons: exon 1 resides in the LSC region, while exons 2 and 3 are duplicated in the IRs. These exons are spliced post-transcriptionally to form a functional mRNA. This complex splicing pattern is highly conserved in land plants and is essential for proper chloroplast gene expression. rps12 encodes a component of the small ribosomal subunit, critical for translation of chloroplast-encoded proteins. Accurate trans-splicing is thus fundamental to chloroplast function and plant development, underscoring the evolutionary significance of maintaining structural and functional integrity in organelle genomes (Oldenburg and Bendich, 2015). The stability of this arrangement also makes rps12 a reliable phylogenetic marker, offering taxonomic utility in resolving species relationships within Lauraceae (Jacobs et al., 2010; Horiuchi and Aigaki, 2006).
The chloroplast genome of D. pugerensis contains 149 simple sequence repeats (SSRs), comprising 53 mononucleotide, 11 dinucleotide, 72 trinucleotide, 10 tetranucleotide, one pentanucleotide, and two hexanucleotide repeats. Trinucleotide repeats—particularly TTA and TAA—are most abundant, followed by A/T-rich mononucleotide repeats (Figure 1D). In contrast, the plastome of D. hainanensis (NC_068504) is dominated by mononucleotide SSRs with relatively few trinucleotide repeats (n=10), reflecting interspecific variation in SSR profiles (Gao et al., 2018). Such patterns may indicate lineage-specific mutation rates or demographic processes such as historical bottlenecks or prolonged population isolation (Dobrogojski et al., 2020). The elevated proportion of trinucleotide SSRs in D. pugerensis may therefore signal unique evolutionary pressures acting on its small, fragmented populations, making these markers valuable for future population-genetic and conservation studies.
Synonymous codon usage analysis in the D. pugerensis chloroplast genome provides insights into translational dynamics and selective pressures shaping plastome evolution. Relative Synonymous Codon Usage (RSCU) analysis revealed a distinct codon bias among protein-coding genes (Figure 1E). AGA (arginine) showed the highest RSCU value (1.79), while CGG (arginine) had the lowest (0.35), indicating a strong preference among synonymous codons. Leucine was the most abundant amino acid, and tryptophan the least (Figure 1E).
A marked bias toward codons ending in A or U was observed, consistent with the AT-rich nature of plastid genomes. Most A/U-ending codons had RSCU > 1, enhancing translational efficiency, while codons such as CGA, GGU, and AGC deviated from this trend. AUG (methionine) and UGG (tryptophan) had neutral RSCU values (1.0), reflecting their lack of synonymous alternatives.
Genomic analyses demonstrate that D. pugerensis shares strong plastome conservation with D. hainanensis, supporting the need for coordinated conservation strategies across their overlapping native range in East Java. Despite this similarity, distinct SSR profiles and codon usage signatures provide powerful genomic markers for species identification and for evaluating adaptive potential. Combined with ongoing conservation research—including population surveys, DNA barcoding using rbcL, matK, trnH–psbA, and ITS (Widjaya et al., 2025), SSR and ISSR-based genetic diversity assessments, and vegetative propagation—the plastome data strengthen our understanding of evolutionary resilience and inform integrated in situ and ex situ conservation planning.
Recent conservation efforts have resulted in the collection of 23 seeds and 15 seedlings of D. pugerensis, now cultivated at the Bogor Botanic Gardens. These ex situ collections serve as essential material for research, propagation trials, and public education (Williams et al., 2015; Westwood et al., 2021). They also support in situ recovery efforts by supplying seeds or planting stock for population reinforcement (Heywood, 2017; Abeli et al., 2020). Field assessments by Helmanto et al. (2022) confirm the species’ extremely restricted distribution and high vulnerability, underscoring the need to incorporate genomic evidence into conservation management to enhance long-term survival.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found at the National Center for Biotechnology Information (NCBI) using accession number PQ560536.1. We also wish to inform you that the chloroplast genome sequence of Dehaasia pugerensis has been published in NCBI under accession number PQ560536.1 GI: 2844831064. (inc ase 18485643).
Author contributions
AHW: Formal analysis, Methodology, Project administration, Writing – original draft. AS: Conceptualization, Funding acquisition, Supervision, Writing – review & editing. TM: Investigation, Supervision, Writing – review & editing. IR: Resources, Supervision, Writing – review & editing. WS: Data curation, Writing – review & editing. MM: Software, Writing – review & editing. IM: Formal analysis, Writing – review & editing. MRH: Visualization, Writing – review & editing. AN: Validation, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This research was funded by the UI PUTI grant in 2024 (No. NKB-376/UN2.RST/HKP.05.00/2024) and the Research Program of the Research Organization for Life Sciences and Environment 2024, BRIN (No. B-1539/III.5/PR.03.06/5/2024).
Acknowledgments
Gratitude is extended to the Perhutani Forestry Institute (PeFI) for providing the materials utilized in this research activity.
Conflict of interest
The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that Generative AI was used in the creation of this manuscript. Regarding AI usage, we confirm that AI tools were used only for grammar refinement and sentence structure improvement, not for data analysis, interpretation, or content generation.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
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.
References
Abeli, T., Dalrymple, S., Godefroid, S., Mondoni, A., Müller, J. V., Rossi, G., et al. (2020). Ex situ collections and their potential for the restoration of extinct plants. Conserv. Biol. 34, 303–313. doi: 10.1111/cobi.13391
An, Y., Wang, Y., Wang, X., and Xiao, J. (2022). Development of chloroplast transformation and gene expression regulation technology in land plants. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.1037038
Andrews, S. (2010). FastQC: a Quality Control Tool for High Throughput Sequence Data. Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (Accessed September 05, 2024).
Beier, S., Thiel, T., Münch, T., Scholz, U., and Mascher, M. (2017). MISA-web: a web server for microsatellite prediction. Bioinformatics 33 (16), 2583–2585. doi: 10.1093/bioinformatics/btx198
Bijmoer, R., Scherrenberg, M., and Creuwels, J. (2020). Naturalis Biodiversity Center (Nl)-Botany. Naturalis Biodiversity Center. Available online at: https://www.gbif.org/occurrence/2514084514 (Accessed September 05, 2024). Occurrence Dataset 10.15468/Ib5ypt.
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Cao, Z., Yang, L., Xin, Y., Xu, W., Li, Q., Zhang, H., et al. (2023). Comparative and phylogenetic analysis of complete chloroplast genomes from seven Neocinnamomum taxa (Lauraceae). Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1205051
Chen, C., Zheng, Y., Liu, S., Zhong, Y., Wu, Y., Li, J., et al. (2017). The complete chloroplast genome of Cinnamomum camphora and its comparison with related Lauraceae species. PeerJ 5, e3820. doi: 10.7717/peerj.3820
Crawford, D. J. and Stuessy, T. F. (2016). Cryptic variation, molecular data, and the challenge of conserving plant diversity in oceanic archipelagos: the critical role of plant systematics. Korean. J. Pl. Taxon. 46, 129–148. doi: 10.11110/kjpt.2016.46.2.129
Dobrogojski, J., Adamiec, M., and Luciński, R. (2020). The chloroplast genome: a review. Acta Physiol. Plant. 42, 98. doi: 10.1007/s11738-020-03089-x
Doyle, J. J. and Doyle, J. L. (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem. Bull. 19 (1), 11–15.
Gao, X., Zhang, X., Meng, H., Li, J., Zhang, D., and Liu, C. (2018). Comparative chloroplast genomes of Paris Sect. Marmorata: insights into repeat regions and evolutionary implications. BMC Genomics 19, 133–144. doi: 10.1186/s12864-018-5281-x
Greiner, S., Lehwark, P., and Bock, R. (2019). OrganellarGenomeDRAW (OGDRAW) version 1.3.1: expanded toolkit for the graphical visualization of organellar genomes. Nucleic Acids Res. 47, W59eW64. doi: 10.1093/nar/gkz238
Han, S., Ding, H., Bi, D., Zhang, S., Yi, R., Gao, J., et al. (2022). Structural diversities and phylogenetic signals in plastomes of the early-divergent angiosperms: a case study in Saxifragales. Plants 11, 3544. doi: 10.3390/plants11243544
Helmanto, H., Primananda, E., Rinandio, D. S., and Robiansyah, I. (2022). Conservation of dehaasia pugerensis (Lauraceae), A tree endemic to east java (Indonesia) and last collected in 1940. J. Nat. Conserv. 65, 1–6. doi: 10.1016/J.Jnc.2021.126096
Heywood, V. H. (2017). The future of plant conservation and the role of botanic gardens. Plant Diversity 39, 309–313. doi: 10.1016/j.pld.2017.12.002
Horiuchi, T. and Aigaki, T. (2006). Alternative trans-splicing: a novel mode of pre-mRNA processing. Biol. Cell 98, 135–140. doi: 10.1042/BC20050002
Jacobs, J., Glanz, S., Bunse-Grassmann, A., Kruse, O., and Kück, U. (2010). RNA trans-splicing: identification of components of a putative chloroplast spliceosome. Eur. J. Cell Biol. 89, 932–939. doi: 10.1016/j.ejcb.2010.06.015
Jin, J. J., Yu, W. B., Yang, J. B., Song, Y., DePamphilis, C. W., Yi, T. S., et al. (2020). GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 21, 1–31. doi: 10.1186/s13059-020-02154-5
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096
Kuznetsov, A. and Bollin, C. J. (2021). NCBI genome workbench: desktop software for comparative genomics, visualization, and GenBank data submission, in Multiple Sequence Alignment: Methods and Protocols. Ed. Katoh, K. (New York: Humana Press), 261–295. doi: 10.1007/978-1-0716-1036-7_16
Lestari, R., Magandhi, M., Hariri, M. R., Noviady, I., Nugroho, A., and Indriani, F. (2024). Characterization of the complete chloroplast genome of the endangered and endemic bornean fruit Artocarpus tamaran Becc. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1513364
Li, Y., Xu, W., Zou, W., Jiang, D., and Liu, X. (2017). Complete chloroplast genome sequences of two endangered Phoebe (Lauraceae) species. Bot. Stud. 58, 1–10. doi: 10.1186/s40529-017-0192-8
Liu, Z. F., Ci, X. Q., Li, L., Li, H. W., Conran, J. G., and Li, J. (2017). DNA barcoding evaluation and implications for phylogenetic relationships in Lauraceae from China. PloS One 12, e0175788. doi: 10.1371/journal.pone.0175788
Okonechnikov, K., Golosova, O., Fursov, M., and Ugene Team (2012). Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics 28, 1166–1167. doi: 10.1093/bioinformatics/bts091
Oldenburg, D. J. and Bendich, A. J. (2015). DNA maintenance in plastids and mitochondria of plants. Frontiers in Plant Science 6, 883. doi: 10.3389/fpls.2015.00883
Rogalski, M., do Nascimento Vieira, L., Fraga, H. P., and Guerra, M. P. (2015). Plastid genomics in horticultural species: importance and applications for plant population genetics, evolution, and biotechnology. Frontiers in plant science 6, 586. doi: 10.3389/fpls.2015.00586
Shi, L., Chen, H., Jiang, M., Wang, L., Wu, X., Huang, L., et al. (2019). CPGAVAS2, an integrated plastome sequence annotator and analyzer. Nucleic Acids Res. 47, W65–W73. doi: 10.1093/nar/gkz345
Song, Y., Yao, X., Liu, B., Tan, Y., and Corlett, R. T. (2017). Complete plastid genome sequence of three tropical Alseodaphne trees in the family Lauraceae. Holzforschung 72 (4), 337–345. doi: 10.1515/hf-2017-0065
Tian, Y., Zhou, J., Zhang, Y., Wang, S., Wang, Y., Liu, H., et al. (2021). Research progress in plant molecular systematics of Lauraceae. Biology 10, 391. doi: 10.3390/biology10050391
Westwood, M., Cavender, N., Meyer, A., and Smith, P. (2021). Botanic garden solutions to the plant extinction crisis. Plant. People. Planet. 3, 22–32. doi: 10.1002/ppp3.10134
Widjaya, A. H., Salamah, A., Maryenti, T., Robiansyah, I., Martiansyah, I., Hariri, M. R., et al. (2025). Phylogenetic and DNA barcoding study on Dehaasia pugerensis Koord. & Valeton, an endemic and critically endangered species from East Java, Indonesia. Beni-Suef. Univ. J. Basic. Appl. Sci. 14, 58. doi: 10.1186/s43088-025-00641-x
Williams, S. J., Jones, J. P., Gibbons, J. M., and Clubbe, C. (2015). Botanic gardens can positively influence visitors’ environmental attitudes. Biodivers. Conserv. 24, 1609–1620. doi: 10.1007/s10531-015-0879-7
Xiao-Ming, Z., Junrui, W., Li, F., Sha, L., Hongbo, P., Lan, Q., et al. (2017). Inferring the evolutionary mechanism of the chloroplast genome size by comparing whole-chloroplast genome sequences in seed plants. Sci. Rep. 7, 1555. doi: 10.1038/s41598-017-01518-5
Keywords: chloroplast genome, conservation, East Java, endemic species, Indonesia
Citation: Widjaya AH, Salamah A, Maryenti T, Robiansyah I, Sun W, Magandhi M, Martiansyah I, Hariri MR and Nugroho A (2026) Plastome of Dehaasia pugerensis Koord. & Valeton: a critically endangered Lauraceae species. Front. Plant Sci. 16:1632459. doi: 10.3389/fpls.2025.1632459
Received: 21 May 2025; Accepted: 11 December 2025; Revised: 25 November 2025;
Published: 09 January 2026.
Edited by:
Yan Cheng, Fujian Agriculture and Forestry University, ChinaReviewed by:
Lihong Xiao, Zhejiang Agriculture and Forestry University, ChinaHoang Dang Khoa Do, Nguyen Tat Thanh University, Vietnam
Copyright © 2026 Widjaya, Salamah, Maryenti, Robiansyah, Sun, Magandhi, Martiansyah, Hariri and Nugroho. 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: Andi Salamah, c2FsYW1haEBzY2kudWkuYWMuaWQ=
Tety Maryenti1