777
Views
3
CrossRef citations to date
0
Altmetric
Mitogenome Announcement

The complete chloroplast genome sequence of Poncirus polyandra (Rutaceae), an endangered species endemic to Yunnan Province, China

, , &
Pages 766-768 | Received 10 Dec 2018, Accepted 23 Dec 2018, Published online: 08 Feb 2019

Abstract

Poncirus polyandra, a plant species with extremely small populations endemic to Fumin county of Yunnan province, has been classified as a national second-class protected wild plant. In this study, we assembled its complete chloroplast genome. The total genome size of P. polyandra was 160,211 bp in length, containing a large single-copy region of 87,406 bp, a small single-copy region of 18,771 bp, and a pair of inverted repeat regions of 27,017 bp. The all GC content of P. polyandra chloroplast genome was 38.4%. It contains 137 function genes, including 92 protein-coding genes, 37 tRNA genes, and eight rRNA ribosomal genes, 14 genes contain a single intron, and three genes have two introns. Phylogenetic analysis results strongly supported that P. polyandra was closely related to the genera of Citrus.

Poncirus Raf. is a genus in the Rutaceae of only two species: P. trifoliate Raf. and P. polyandra S. Q. Ding (Fang Citation1993; Fang et al. Citation1999). Of them, P. polyandra, also known as wild orange and used as a dwarf stock of Citrus reticulate, with extremely small populations endemic to Fumin county of Yunnan province, had been classified as a national second-class protected wild plant in the Information System of Chinese Rare and Endangered Plants (ISCRPE) (http://rep.iplant.cn/protlist). However, it has been extinct in the wild due to the drastic changes in the natural environment (Li et al. Citation2010) and the excessive excavation of human beings (Gao et al. Citation2017). So it is necessary to protect P. polyandra germplasm resources.

Chloroplast genomes are widely used for phylogeny (Xue et al. Citation2012), DNA barcoding (Dong et al. Citation2012, Citation2014), genome evolution and species conservation (Dong et al. Citation2013). So far, the chloroplast genome such as Clauseneae tribe (Shivakumar et al. Citation2017), genus Citrus (Bausher et al. Citation2006; Yang et al. Citation2016; Liu and Shi Citation2017), Phellodendron amurense (Chen Citation2018) and Zanthoxylum bungeanum (Liu and Wei Citation2017) within the family of Rutaceae have been reported, while the plastome of P. polyandra has not been reported. Here, we reported the complete chloroplast genome sequence of P. polyandra based on the next-generation sequencing, and the annotated genomic sequence was submitted to GenBank under accession number MK250977.

The fresh leaves of P. polyandra were collected from three-year-old seedling planted in the greenhouse of Southwest Forestry University, Kunming, China. Total genome DNA was extracted with the Ezup plant genomic DNA prep kit (Sangon Biotech, Shanghai, China), and DNA samples were properly stored at the Key Laboratory of State Forestry Administration on Biodiversity Conservation in Southwest China, Southwest Forestry University, Kunming, China.

Total DNA was used to generate libraries with an average insert size of 350 bp, which were sequenced using the Illumina HiSeq X platform. Approximately 14.0 GB of raw data were generated with 150 bp paired-end read lengths. Then, the raw data were used to assemble the complete cp genome using GetOrganelle software (Jin et al. Citation2018) with Citrus aurantifolia as the reference. Genome annotation was performed with the program Geneious R8 (Biomatters Ltd, Auckland, New Zealand) by comparing the sequences with the cp genome of Citrus aurantifolia. The tRNA genes were further confirmed through online tRNAscan-SE web servers (Schattner et al. Citation2005). A gene map of the annotated P. polyandra cp genome was drawn by OGdraw online (Lohse et al. Citation2013).

The complete chloroplast genome of P. polyandra is a circular molecule of 160,211 bp in length comprising a large single copy (LSC) region of 87,406 bp and a small single copy (SSC) region of 18,771 bp separated by a pair of inverted repeats (IRs), each 27,017 bp in length. The all GC content of P. polyandra chloroplast genome was 38.4%, while the GC contents of the LSC, SSC, and IR regions were 36.8%, 33.2%, 42.9%, respectively.

The chloroplast genome of P. polyandra contains 137 function genes, including 92 protein-coding genes, 37 tRNA genes, and eight rRNA ribosomal genes. Ten protein-coding genes, 7 tRNA genes, and all rRNA genes were duplicated in the IR regions. Among 116 unique genes, 14 genes (atpF, ndhA, ndhB, petB, petD, rpl2, rpoC1, rps16, trnA-UGC, trnG-GCC, trnI-GAU, trnK-UUU, trnL-UAA, and trnV-UAC) contain a single intron, and three genes (clpP, rps12, and ycf3) have two introns.

To confirm the phylogenetic location of P. polyandra within the family of Rutaceae, a total of 16 complete cp genomes of Rutaceae were obtained from GenBank, and two species, Dimocarpus longan and Acer buergerianum in the family of Sapindaceae were used as out-groups. The 19 complete chloroplast sequences were aligned by the MAFFT version 7 software (Katoh and Standley Citation2013). Phylogenetic analysis was conducted based on maximum likelihood (ML) analyses were implemented in IQ-TREE 1.5.5 (Nguyen et al. Citation2015) under the GTR + F+R3 nucleotide substitution model, which was selected by ModelFinder (Kalyaanamoorthy et al. Citation2017). Support for the inferred ML tree was inferred by bootstrapping with 1000 replicates. Phylogenetic analysis results strongly supported that P. polyandra was closely related to the genera of Citrus (). The chloroplast genome of P. polyandra will provide useful genetic information for further study on genetic diversity and conservation of Poncirus species.

Figure 1. Phylogenetic relationships among 17 complete chloroplast genomes of Rutaceae. Bootstrap support values are given at the nodes. Chloroplast genome accession number used in this phylogeny analysis: Acer buergerianum: KY419137; Citrus aurantiifolia: KJ865401; Citrus depressa: LC147381; Citrus platymamma: KR259987; Citrus sinensis: DQ864733; Clausena excavata: KU949003; Dimocarpus longan: MG214255; Glycosmis mauritiana: KU949004; Glycosmis pentaphylla: KU949005; Merrillia caloxylon: KU949006; Micromelum minutum: KU949007; Murraya koenigii: KU949002; Phellodendron amurense: KY707335; Poncirus polyandra: MK250977; Zanthoxylum bungeanum: KX497031; Zanthoxylum piperitum: KT153018; Zanthoxylum schinifolium: KT321318; Zanthoxylum simulans: MF716524; Zanthoxylum sp. NH-2018: MF716521.

Figure 1. Phylogenetic relationships among 17 complete chloroplast genomes of Rutaceae. Bootstrap support values are given at the nodes. Chloroplast genome accession number used in this phylogeny analysis: Acer buergerianum: KY419137; Citrus aurantiifolia: KJ865401; Citrus depressa: LC147381; Citrus platymamma: KR259987; Citrus sinensis: DQ864733; Clausena excavata: KU949003; Dimocarpus longan: MG214255; Glycosmis mauritiana: KU949004; Glycosmis pentaphylla: KU949005; Merrillia caloxylon: KU949006; Micromelum minutum: KU949007; Murraya koenigii: KU949002; Phellodendron amurense: KY707335; Poncirus polyandra: MK250977; Zanthoxylum bungeanum: KX497031; Zanthoxylum piperitum: KT153018; Zanthoxylum schinifolium: KT321318; Zanthoxylum simulans: MF716524; Zanthoxylum sp. NH-2018: MF716521.

Disclosure statement

No potential conflict of interest was reported by the authors.

Additional information

Funding

The study was supported by the National Natural Science Foundation of China under Grant [Number 31460205, 31860219, 31360184].

References

  • Bausher MG, Singh ND, Lee SB, Jansen RK, Daniell H. 2006. The complete chloroplast genome sequence of Citrus sinensis (L.) Osbeck var 'Ridge Pineapple': organization and phylogenetic relationships to other angiosperms. BMC Plant Biol. 6:21.
  • Chen KK. 2018. Characterization of the complete chloroplast genome of the Tertiary relict tree Phellodendron amurense (Sapindales: Rutaceae) using Illumina sequencing technology. Conservation Genet Resour. 10:43–46.
  • Dong W, Liu H, Xu C, Zuo Y, Chen Z, Zhou S. 2014. A chloroplast genomic strategy for designing taxon specific DNA mini-barcodes: a case study on ginsengs. BMC Genet. 15:138.
  • Dong W, Liu J, Yu J, Wang L, Zhou S. 2012. Highly variable chloroplast markers for evaluating plant phylogeny at low taxonomic levels and for DNA barcoding. PLoS One. 7:e35071.
  • Dong W, Xu C, Cheng T, Zhou S. 2013. Complete chloroplast genome of Sedum sarmentosum and chloroplast genome evolution in Saxifragales. PLoS One. 8:e77965.
  • Fang DQ. 1993. Intra and intergeneric relationships of Poncirus polyandra: investigation by leaf isozymes. J Wuhan Botanic Res. 11:34–40.
  • Fang ZF, Zhao SD, Skvortsov AK. 1999. Flora of China (English version). Vol 43. Beijing: Science Press.
  • Gao XS, Zhang YP, Tang JW, Gu WH. 2017. Discussion on survival status, protection and development of Poncirus polyandra. Prot Forest Sci Technol. 12:83–85.
  • Jin JJ, Yu WB, Yang JB, Song Y, Yi TS, Li DZ. 2018. GetOrganelle: a simple and fast pipeline for de novo assemble of a complete circular chloroplast genome using genome skimming data. BioRxiv [in Preprint]. http://doi.org/10.1101/256479.
  • Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. 2017. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 14:587–589.
  • Katoh K, Standley DM. 2013. MAFFT Multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 30:772–780.
  • Li G, Ouyang ZQ, Qi Y, Lu LJ. 2010. Protection measures on rare and endangered plants in Kunming region. Environ Sci Survey. 29:27–29.
  • Liu J, Shi C. 2017. The complete chloroplast genome of wild shaddock, Citrus maxima (Burm.) Merr. Conservation Genet Resour. 9:599–601.
  • Liu YL, Wei AZ. 2017. The complete chloroplast genome sequence of an economically important plant, Zanthoxylum bungeanum (Rutaceae). Conservation Genet Resour. 9:25–27.
  • Lohse M, Drechsel O, Kahlau S, Bock R. 2013. Organellar GenomeDRAW-a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucleic Acids Res. 41:575–581.
  • Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ. 2015. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 32:268–274.
  • Schattner P, Brooks AN, Lowe TM. 2005. The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res. 33:686–689.
  • Shivakumar VS, Appelhans MS, Johnson G, Carlsen M, Zimmer EA. 2017. Analysis of whole chloroplast genomes from the genera of the Clauseneae, the curry tribe (Rutaceae, Citrus family). Mol Phylogenet Evol. 117:135–140.
  • Xue JH, Dong WP, Cheng T, Zhou SL. 2012. Nelumbonaceae: systematic position and species diversification revealed by the complete chloroplast genome. J Syst Evol. 50:477–487.
  • Yang K, Lee J, Lee SC, Yang TJ, Kim H, Chung MY, Nou IS. 2016. The complete chloroplast genome of Korean popular Citrus hybrid Hallabong mandarin [(Citrus unshiu × C. sinensis) × C. reticulate] (Rutaceae). Mitochondrial DNA Part B. 1:29–30.