The complete mitochondrial genome of <i>Arabidopsis thaliana</i> (Brassicaceae) isolated in Korea

Jongsun PARK; Hong XI; Yongsung KIM

doi:10.11110/kjpt.2021.51.2.176

Korean J. Pl. Taxon > Volume 51(2); 2021 > Article

/home/virtual/kjpt/journal//../xmls/kjpt-51-2-176.xml PARK, XI, and KIM: The complete mitochondrial genome of Arabidopsis thaliana (Brassicaceae) isolated in Korea

Short Communication

Korean Journal of Plant Taxonomy 2021; 51(2): 176-180.

Published online: June 30, 2021

DOI: https://doi.org/10.11110/kjpt.2021.51.2.176

The complete mitochondrial genome of Arabidopsis thaliana (Brassicaceae) isolated in Korea

Jongsun PARK^*

, Hong XI, Yongsung KIM¹

InfoBoss Inc. and InfoBoss Research Center, Seoul 06278, Korea

¹Honam National Institute of Biological Resources, Mokpo 58762, Korea

^*Author for correspondence: starflr@infoboss.co.kr

Received May 07, 2021 Revised May 21, 2021 Accepted June 05, 2021

(open-access):

This is an open-access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Arabidopsis thaliana (L.) Heynh. is a small plant species that serves as a model organism of plant biology and genetics. Here, we present the first complete mitochondrial genome of Korean A. thaliana natural isolate (named as 180404IB4), which is 368,875 bp long and contains 58 genes (33 protein-coding genes, 22 tRNAs, and three rRNAs), with a GC ratio of 44.8%. Sixty-four single-nucleotide polymorphisms and 11 insertion and deletion regions (1,089 bp in length) are identified against the Col-0 ecotype, showing one large insertion of 1,069 bp without structural variation. Phylogenetic trees constructed from 30 conserved genes indicate that the 180404IB4 mitochondrial genome is clustered with Col-0 and three East Asian ecotypes.

Keywords: Arabidopsis thaliana, Brassicaceae, Korean individual, intraspecific variations, mitochondrial genome

Arabidopsis thaliana (L.) Heynh. is a small plant species distributed in Eurasia and Africa (Hoffmann, 2002). It is a model organism of plant biology and genetics (Rensink and Buell, 2004) because of its rapid life cycle (ca. 10 weeks) (Mandel and Yanofsky, 1995) and relatively small genome size (ca. 119 Mb) (The Arabidopsis Genome Initiative, 2000). As next-generation sequencing technologies have been established, above one thousand whole genomes have been sequenced (Gan et al., 2011; Schmitz et al., 2013): e.g., 1,135 natural isolates from Eurasian and North African (The 1001 Genomes Consortium, 2016) and 118 strains from Yangtze River in China (Zou et al., 2017). In contrast, only six complete chloroplast (Park et al., 2020a) and three mitochondrial genomes (mitogenomes) (Unseld et al., 1997; Davila et al., 2011; Sloan et al., 2018) are available even though they are valuable to investigate phylogenetic relationship (Park et al., 2020a, 2020c, 2021). Here, we presented the complete mitogenome of A. thaliana isolated in Korea together with three East Asian isolates to understand their phylogenetic relationship.

Materials and Methods

Plant material

We collected an individual accession of A. thaliana in the population located in Yeonggwang-gun, Jeollanam-do province in Korea (180404IB4). A voucher and isolated DNA were deposited in the InfoBoss Cyber Herbarium (IN, the voucher number IB-00925).

DNA extraction and mitochondrial genome determination

Its total DNA was extracted from fresh leaves by using a DNeasy Plant Mini Kit (QIAGEN, Hilden, Germany). Genome sequencing was performed using HiSeqX at Macrogen Inc., Korea, and de novo assembly was done by Velvet v1.2.10 (Zerbino and Birney, 2008) and GapCloser v1.12 (Zhao et al., 2011). Assembled sequences were modified and confirmed by BWA v0.7.17 (Li, 2013) and SAMtools v1.9 (Li et al., 2009). All bioinformatic analyses were conducted in the Genome Information System (http://geis.infoboss.co.kr/) like the previous studies (Kwon et al., 2019a, 2019b; Park et al., 2019a, 2019d, 2020b; Min et al., 2020; Choi et al., 2021a, 2021b; Kim et al., 2021). The same method was applied to assemble mitogenomes of three East Asian A. thaliana isolates based on public raw reads (Table 1).

Genome annotation was conducted based on Col-0 mitogenome (NC_037304) (Sloan et al., 2018) with Geneious Prime 2020.2.4 (Biomatters Ltd., Auckland, New Zealand). A circular map of A. thaliana mitogenome was drawn using OGDRAW v1.31 (Greiner et al., 2019).

Phylogenetic analysis

Maximum-Likelihood (ML) and Bayesian inference (BI) phylogenetic trees were constructed based on the concatenated alignment of the 30 aligned conserved genes from the seven A. thaliana including Korean individual and three ecotypes (Table 1) and Arabis alpina (Xu and Bi, 2018) mitogenomes (Table 1) by MAFFT v7.450 (Katoh and Standley, 2013). The ML tree was reconstructed in MEGA X with 1,000 bootstrap repeats (Kumar et al., 2018). In the ML analysis, a heuristic search was used with nearest-neighbor interchange branch swapping, TVM + F + R4 model determined by jModelTest v2.1.5 (Darriba et al., 2012), and uniform rates among sites. All other options used the default settings. The posterior probability of each node was estimated by BI using MrBayes v3.2.6 (Huelsenbeck and Ronquist, 2001) plug-in implemented in Geneious Prime 2020.2.4 (Biomatters Ltd.). The HKY85 model with gamma rates was used as a molecular model. A Markov chain Monte Carlo algorithm was employed for 1,100,000 generations, sampling trees every 200 generations, with four chains running simultaneously. Trees from the first 100,000 generations were discarded as burn-in.

Data availability

Mitochondrial genome sequence of Korean A. thaliana can be accessed via accession number MK358445 in NCBI GenBank. The associated BioProject, SRA, and Bio-Sample numbers are PRJNA727818, SAMN19040823, and SRR14458670, respectively.

Results and Discussion

The complete mitogenome of A. thaliana 180404IB4 (MK358445) is 368,875 bp long (Fig. 1A), longer than three previously assembled and one newly assembled A. thaliana mitogenomes (Table 1). Overall GC content is 44.8% and it contains 58 genes (33 protein-coding genes, 22 tRNAs, and three rRNAs) (Fig. 1A), which are the same as the Col-0 mitogenome.

Sixty-four single nucleotide polymorphisms (SNPs; 0.017%) and 11 insertion and deletion (INDEL) regions (1,089 bp in length; 0.30%) were identified against Col-0 mitogenome. Interestingly, one 1,069-bp insertion was found, similar to the insertions identified in chloroplast genomes of Coffea arabica L. (Park et al., 2019c) and Duchesnea chrysantha (Zoll. & Moritzi) Miq. (Park et al., 2019b) and mitogenomes of Liriodendron tulipifera L. (Park et al., 2019a) and Populus alba x Populus glandulosa (Park et al., 2019d). Proportion of intraspecific variations identified between 180404IB4 and Col-0 is similar to that of Rosa rugosa Thunb. (124 SNPs, 0.041% and 769-bp INDELs, 0.25%) (Park et al., 2020b), is less than those of Liriodendron tulipifera L. (365 SNPs, 0.07% and 2,117-bp INDELs, 0.38%) (Park et al., 2019a) and Glycyrrhiza uralensis Fisch. ex DC. (1,099 SNPs, 0.24% and 1,736-bp INDELs, 0.38%) (Baek et al., under review), but is greater in number than that of Malus x domestica Borkh. (140 SNPs, 0.035% and 6-bp INDELs, 0.0015%) (Goremykin et al., 2012; Ge et al., 2020).

No structural variation among 180404IB4, 15-11, and Col-0 mitogenomes was found, while the remaining A. thaliana genomes (Table 1) presented structural variations (Davila et al., 2011). It is also congruent to the recent analysis which presents a dynamic structure of A. thaliana mitogenome along with different ecotypes (Masutani et al., 2021).

Phylogenetic trees constructed from the conserved genes of the seven A. thaliana and Arabis alpina mitogenomes (Table 1) show that four East Asian Arabidopsis including 180404IB4 were clustered with the Col-0 (USA) with high supportive values and Ler-0 (Germany) and C-24 (Tanzania) were clustered in an independent clade (Fig. 1B). It is different from the phylogenetic tree constructed from the complete chloroplast genomes, displaying the intraspecific topology of the six A. thaliana ecotypes except for C-24 with high supportive values (Park et al., 2020a). It might be caused by the low coverage of the 30 conserved genes (6.98%) against the complete mitogenome in comparison to that of the chloroplast genomes. To improve the resolution of phylogenetic analysis using mitogenomes, more conserved areas except genic regions should be rescued and addition mitogenomes of East Asian A. thaliana are needed to be assembled to investigate their phylogenetic relationship. Throughout this study, this mitogenome presents that Korean A. thaliana is phylogenetically related to the other East Asian isolates and Col-0.

ACKNOWLEDGMENTS

This work was supported by InfoBoss Research Grant (IBG-0023).

NOTES

Conflict of Interest

The authors declare that there are no conflicts of interest.

Fig. 1.

A. Circular map of mitochondrial genome of Arabidopsis thaliana isolated in Korea. Gene shown outside are transcribed clockwise, and inside the circle are transcribed counter clockwise. Genes are colour-coded to distinguish different functional groups. The dark grey and the light grey plot in the inner circle correspond to the GC content and AT content, respectively. B. Maximum-Likelihood (ML) and Bayesian inference (BI) phylogenetic trees of 30 conserved genes from seven Arabidopsis and one Arabis complete mitochondrial genomes. Phylogenetic tree was drawn based on ML tree. The numbers above branches indicate bootstrap support value of ML tree and posterior probability of BI tree, respectively.

Table 1.

List of the mitochondrial genomes used in this study

Strain species	Origin	Complete genome	Accession number	Length (bp)	GC ratio (%)	Reference / Source
180404IB4	Korea	Y	MK358445	368,875	44.8	This study
Col-0	USA	Y	NC_037304	367,808	44.8	Sloan et al. 2018
Ler-0	Germany	Y	JF729202	345,418	44.8	Davila et al. 2011
C24	Tanzania	Y	JF729200	366,922	44.8	Davila et al. 2011
Kyoto	Japan	Y	N/A	350,265	44.8	SRR492307 (Schmitz et al. 2013)
Tsu-0	Japan	N	N/A	389,442	44.7	ERR031555 (Gan et al. 2011)
15-11	China	Y	N/A	368,875	44.8	SRR2204166 (Zou et al. 2017)
Arabis alpina		Y	NC_037070	323,159	41.3	Xu and Bi 2018

Literature Cited

Choi, SS. Bakalin, VA. Kwon, W and Park, J. 2021a. The complete mitochondrial genome of Douinia plicata (Lindb.) Konstant. et. Vilnet (Scapaniaceae, Jungermanniales). Mitochondrial DNA Part B 6: 789-791.

Choi, SS. Min, J. Kwon, W and Park, J. 2021b. The complete mitochondrial genome of Scapania ampliata Steph., 1897 (Scapaniaceae, Jungermanniales). Mitochondrial DNA Part B Resources 6: 686-688.

Darriba, D. Taboada, GL. Doallo, R and Posada, D. 2012. jModelTest 2: more models, new heuristics and parallel computing. Nature Methods 9: 772.

Davila, JI. Arrieta-Montiel, MP. Wamboldt, Y. Cao, J. Hagmann, J. Shedge, V. Xu, Y.-Z. Weigel, D and Mackenzie, SA. 2011. Double-strand break repair processes drive evolution of the mitochondrial genome in Arabidopsis . BMC Biology 9: 64.

Gan, X. Stegle, O. Behr, J. Steffen, JG. Drewe, P. Hildebrand, KL. Lyngsoe, R. Schultheiss, SJ. Osborne, EJ. Sreedharan, VT. Kahles, A. Bohnert, R. Jean, G. Derwent, P. Kersey, P. Belfield, EJ. Harberd, NP. Kemen, E. Toomajian, C. Kover, PX. Clark, RM. Rätsch, G and Mott, R. 2011. Multiple reference genomes and transcriptomes for Arabidopsis thaliana . Nature 477: 419-423.

Ge, D. Dong, J. Guo, L. Yan, M. Zhao, X and Yuan, Z. 2020. The complete mitochondrial genome sequence of cultivated apple (Malus domestica cv.‘Yantai Fuji 8’). Mitochondrial DNA Part B Resources 5: 1317-1318.

Goremykin, VV. Lockhart, PJ. Viola, R and Velasco, R. 2012. The mitochondrial genome of Malus domestica and the import-driven hypothesis of mitochondrial genome expansion in seed plants. The Plant Journal 71: 615-626.

Greiner, S. Lehwark, P and Bock, R. 2019. OrganellarGenome-DRAW (OGDRAW) version 1.3. 1: expanded toolkit for the graphical visualization of organellar genomes. Nucleic Acids Research 47: W59-W64.

Hoffmann, MH. 2002. Biogeography of Arabidopsis thaliana (L.) Heynh. (Brassicaceae). Journal of Biogeography 29: 125-134.

Huelsenbeck, JP and Ronquist, F. 2001. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17: 754-755.

Katoh, K and Standley, DM. 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Molecular Biology and Evolution 30: 772-780.

Kim, M. Xi, H and Park, J. 2021. Genome-wide comparative analyses of GATA transcription factors among 19 Arabidopsis ecotype genomes: Intraspecific characteristics of GATA transcription factors. PLoS One 16: e0252181.

Kumar, S. Stecher, G. Li, M. Knyaz, C and Tamura, K. 2018. MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Molecular Biology and Evolution 35: 1547-1549.

Kwon, W. Kim, Y and Park, J. 2019a. The complete mitochondrial genome of Dumortiera hirsuta (Sw.) Nees (Dumortieraceae, Marchantiophyta). Mitochondrial DNA Part B Resources 4: 1586-1587.

Kwon, W. Kim, Y and Park, J. 2019b. The complete mitochondrial genome of Korean Marchantia polymorpha subsp. rude-ralis Bischl. & Boisselier: inverted repeats on mitochondrial genome between Korean and Japanese isolates. Mitochondrial DNA Part B Resources 4: 769-770.

Li, H. 2013. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. Preprint at: https://arxiv.org/abs/1303.3997 . 2013.

Li, H. Handsaker, B. Wysoker, A. Fennell, T. Ruan, J. Homer, N. Marth, G. Abecasis, G and Durbin, R. 1000. Genome Project Processing Subgroup. 2009. The sequence alignment/map format and SAMtools. Bioinformatics 25: 2078-2079.

Mandel, MA and Yanofsky, MF. 1995. A gene triggering flower formation in Arabidopsis . Nature 377: 522-524.

Masutani, B. Arimura, S.-I and Morishita, S. 2021. Investigating the mitochondrial genomic landscape of Arabidopsis thaliana by long-read sequencing. PLoS Computational Biology 17: e1008597.

Min, J. Kwon, W. Xi, H and Park, J. 2020. The complete mitochondrial genome of Riccia fluitans L. (Ricciaceae, Marchantiophyta): Investigation of intraspecific variations on mitochondrial genomes of R. fluitans . Mitochondrial DNA Part B Resources 5: 1220-1222.

Park, J. Kim, Y and Kwon, M. 2019a. The complete mitochondrial genome of tulip tree, Liriodendron tulifipera L. (Magnoliaceae): Intra-species variations on mitochondrial genome. Mitochondrial DNA Part B Resources 4: 1308-1309.

Park, J. Kim, Y and Lee, K. 2019b. The complete chloroplast genome of Korean mock strawberry, Duchesnea chrysantha (Zoll. & Moritzi) Miq. (Rosoideae). Mitochondrial DNA Part B Resources 4: 864-865.

Park, J. Kim, Y. Xi, H and Heo, K.-I. 2019c. The complete chloroplast genome of ornamental coffee tree, Coffea arabica L. (Rubiaceae). Mitochondrial DNA Part B Resources 4: 1059-1060.

Park, J. Kim, Y. Xi, H. Kwon, W and Kwon, M. 2019d. The complete chloroplast and mitochondrial genomes of Hyunsasi tree, Populus alba x Populus glandulosa (Salicaceae). Mitochondrial DNA Part B Resources 4: 2521-2522.

Park, J. Min, J. Kim, Y and Chung, Y. 2021. The comparative analyses of six complete chloroplast genomes of morphologically diverse Chenopodium album L. (Amaranthaceae) collected in Korea. International Journal of Genomics 2021: 6643444.

Park, J. Xi, H and Kim, Y. 2020a. The complete chloroplast genome of Arabidopsis thaliana isolated in Korea (Brassicaceae): an investigation of intraspecific variations of the chloroplast genome of Korean A. thaliana . International Journal of Genomics 2020: 3236461.

Park, J. Xi, H. Kim, Y. Nam, S and Heo, K.-I. 2020b. The complete mitochondrial genome of new species candidate of Rosa rugosa (Rosaceae). Mitochondrial DNA Part B Resources 5: 3435-3437.

Park, J. Xi, H and Oh, S.-H. 2020c. Comparative chloroplast genomics and phylogenetic analysis of the Viburnum dilatatum complex (Adoxaceae) in Korea. Korean Journal of Plant Taxonomy 50: 8-16.

Rensink, WA and Buell, CR. 2004. Arabidopsis to rice: applying knowledge from a weed to enhance our understanding of a crop species. Plant Physiology 135: 622-629.

Schmitz, RJ. Schultz, MD. Urich, MA. Nery, JR. Pelizzola, M. Libiger, O. Alix, A. McCosh, RB. Chen, H. Schork, NJ and Ecker, JR. 2013. Patterns of population epigenomic diversity. Nature 495: 193-198.

Sloan, DB. Wu, Z and Sharbrough, J. 2018. Correction of persistent errors in Arabidopsis reference mitochondrial genomes. The Plant Cell 30: 525-527.

The 1001 Genomes Consortium. 2016. 1,135 genomes reveal the global pattern of polymorphism in Arabidopsis thaliana . Cell 166: 481-491.

The Arabidopsis Genome Initiative. 2000. Analysis of the genome sequence of the flowering plant Arabidopsis thaliana . Nature 408: 796-815.

Unseld, M. Marienfeld, P. Brandt, P and Brennicke, A. 1997. The mitochondrial genome of Arabidopsis thaliana contains 57 genes in 366,924 nucleotides. Nature Genetics 15: 57-61.

Xu, Y and Bi, C. 2018. The complete mitochondrial genome sequence of an alpine plant Arabis alpina . Mitochondrial DNA Part B Resources 3: 725-727.

Zerbino, DR and Birney, E. 2008. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Research 18: 821-829.

Zhao, Q.-Y. Wang, Y. Kong, Y.-M. Luo, D. Li, X and Hao, P. 2011. Optimizing de novo transcriptome assembly from shortread RNA-Seq data: a comparative study. BMC Bioinformatics 12 Suppl 14: S2.

Zou, Y.-P. Hou, X.-H. Wu, Q. Chen, J.-F. Li, Z.-W. Han, T.-S. Niu, X.-M. Yang, L. Xu, Y.-C. Zhang, Jx. Zhang, J. Tan, D. Tian, Z. Gu, H and Guo, Y.-L. 2017. Adaptation of Arabidopsis thaliana to the Yangtze River basin. Genome Biology 18: 239.