Organellar genome assembly methods and comparative analysis of horticultural plants

Abstract

Although organellar genomes (including chloroplast and mitochondrial genomes) are smaller than nuclear genomes in size and gene number, organellar genomes are very important for the investigation of plant evolution and molecular ecology mechanisms. Few studies have focused on the organellar genomes of horticultural plants. Approximately 1193 chloroplast genomes and 199 mitochondrial genomes of land plants are available in the National Center for Biotechnology Information (NCBI), of which only 39 are from horticultural plants. In this paper, we report an innovative and efficient method for high-quality horticultural organellar genome assembly from next-generation sequencing (NGS) data. Sequencing reads were first assembled by Newbler, Amos, and Minimus software with default parameters. The remaining gaps were then filled through BLASTN search and PCR. The complete DNA sequence was corrected based on Illumina sequencing data using BWA (Burrows–Wheeler Alignment tool) software. The advantage of this approach is that there is no need to isolate organellar DNA from total DNA during sample preparation. Using this procedure, the complete mitochondrial and chloroplast genomes of an ornamental plant, Salix suchowensis, and a fruit tree, Ziziphus jujuba, were identified. This study shows that horticultural plants have similar mitochondrial and chloroplast sequence organization to other seed plants. Most horticultural plants demonstrate a slight bias toward A+T rich features in the mitochondrial genome. In addition, a phylogenetic analysis of 39 horticultural plants based on 15 protein-coding genes showed that some mitochondrial genes are horizontally transferred from chloroplast DNA. Our study will provide an important reference for organellar genome assembly in other horticultural plants. Furthermore, phylogenetic analysis of the organellar genomes of horticultural plants could accurately clarify the unanticipated relationships among these plants.

References

[1] Richman, V., Bennett, J., Jackson, R.S. et al. Horticulture-Plant needs, Horticultural plants. Science Encyclopedia.Web. 20 Dec 2017. http://science.jrank.org/pages/3392/Horticulture.html.

[2] Pyke, K. A. Plastid division and development. Plant Cell 11, 549-556 (1999).

[3] Mccoy, S. R., Kuehl, J. V., Boore, J. L. & Raubeson, L. A. The complete plastid genome sequence of Welwitschia mirabilis: an unusually compact plastome with accelerated divergence rates. BMC Evol. Biol. 8, 130 (2008).

[4] Wu, C. S., Wang, Y. N., Liu, S. M., & Chaw, S. M. Chloroplast Genome (cpDNA) of Cycas taitungensis and 56 cp protein-coding genes of Gnetum parvifolium: insights into cpDNA evolution and phylogeny of extant seed plants. Mol. Biol. Evol. 24, 1366-1379 (2007).

[5] Yang, M. et al. The complete chloroplast genome sequence of date palm (Phoenix dactylifera L.). PLoS ONE 5, e12762 (2012).

[6] Kubo, T. & Newton, K. J. Angiosperm mitochondrial genomes and mutations. Mitochondrion 8, 5-14 (2008).

[7] Tanaka, Y., Tsuda, M., Yasumoto, K., Yamagishi, H. & Terachi, T. A complete mitochondrial genome sequence of Ogura-type male-sterile cytoplasm and its comparative analysis with that of normal cytoplasm in radish (Raphanus sativus L.). BMC Genom. 13, 1-12 (2012).

[8] Alverson, A. J. et al. Insights into the evolution of mitochondrial genome size from complete sequences of Citrullus lanatus and Cucurbita pepo (Cucurbitaceae). Mol. Biol. Evol. 27, 1436-1448 (2010).

[9] Alverson, A. J., Zhuo, S., Rice, D. W., Sloan, D. B. & Palmer, J. D. The mitochondrial genome of the legume Vigna radiata and the analysis of recombination across short mitochondrial repeats. PLoS ONE 6, e16404 (2011).

[10] Jeong, Y. M. et al. The complete mitochondrial genome of cultivated radish WK10039 (Raphanus sativus L.). Mitochondrial DNA A DNA Mapp. Seq. Anal. 27, 1-2 (2014).

[11] Alverson, A. J. et al. Insights into the evolution of mitochondrial genome size from complete sequences of Citrullus lanatus and Cucurbita pepo (Cucurbitaceae). Mol. Biol. Evol. 27, 1436 (2010).

[12] Atherton, R. A. et al. Whole genome sequencing of enriched chloroplast DNA using the Illumina GAII platform. Plant Methods 6, 1-6 (2010).

[13] Unseld, M., Marienfeld, J. R., Brandt, P. & Brennicke, A. The mitochondrial genome of Arabidopsis thaliana contains 57 genes in 366,924 nucleotides. Nat. Genet. 15, 57-61 (1997).

[14] Handa, H. The complete nucleotide sequence and RNA editing content of the mitochondrial genome of rapeseed (Brassica napus L.): comparative analysis of the mitochondrial genomes of rapeseed and Arabidopsis thaliana. Nucleic Acids Res. 31, 5907 (2003).

[15] Ogihara, Y. et al. Structural dynamics of cereal mitochondrial genomes as revealed by complete nucleotide sequencing of the wheat mitochondrial genome. Nucleic Acids Res. 33, 6235-6250 (2005).

[16] Jansen, R. K. et al. Methods for obtaining and analyzing whole chloroplast genome sequences. Methods Enzymol. 395, 348-384 (2010).

[17] Cronn, R. et al. Multiplex sequencing of plant chloroplast genomes using Solexa sequencing-by-synthesis technology. Nucleic Acids Res. 36, e122- e122 (2008).

[18] Moore, M. J. et al. Rapid and accurate pyrosequencing of angiosperm plastid genomes. BMC Plant. Biol. 6, 1-13 (2006).

[19] Tangphatsornruang, S. et al. The chloroplast genome sequence of mungbean (Vigna radiata) determined by high-throughput pyrosequencing: structural organization and phylogenetic relationships. DNA Res. 17, 11-22 (2010).

[20] Kumar, S., Stecher, G. & Tamura, K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870 (2016).

[21] Simon, P. W. et al. De novo assembly and characterization of the carrot mitochondrial genome using next generation sequencing data from whole genomic DNA provides first evidence of DNA transfer into an angiosperm plastid genome. BMC Plant Biol. 12, 1-17 (2012).

[22] Ma, Q. et al. Complete chloroplast genome sequence of a major economic species, Ziziphus jujuba (Rhamnaceae). Curr. Genet. 63, 1-13 (2017).

[23] Ma, Q. et al. Identification and characterization of nucleotide variations in the genome of Ziziphus jujuba (Rhamnaceae) by next generation sequencing. Mol. Biol. Rep. 41, 3219-3223 (2014).

[24] Camacho C. et al. BLAST plus: architecture and applications. BMC Bioinformatics 10, 421 (2009).

[25] Treangen, T. J., Sommer, D. D., Angly, F. E., Sergey, K. & Mihai, P. Next generation sequence assembly with AMOS. Curr. Protoc. Bioinformatics Chapter 11, 11.18.11-11.18.18 (2011).

[26] Sommer, D. D., Delcher, A. L., Salzberg, S. L. & Pop, M. Minimus: a fast, lightweight genome assembler. BMC Bioinformatics 8, 1-64 (2007).

[27] Ewing, B. & Green, P. Base-calling of automated sequencer traces using Phred. II error probabilities. Genome Res. 8, 186-194 (1998).

[28] Peters, D., Qiu, K., Liang, P. Faster short DNA sequence alignment with parallel BWA. AIP Conf. Proc. 1368, 131-134 (2011).

[29] Surhone, L. M., Tennoe, M. T., Henssonow, S. F., Group, T. O., & Done, G. T. OmniGraffle (Betascript Publishing, Beau Bassin, Mauritius, 2010).

[30] Zhao, K. & Chu, X. G-BLASTN: accelerating nucleotide alignment by graphics processors. Bioinformatics 30, 1384-1391 (2014).

[31] Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. Preprint at https://arxiv.org/pdf/1303.3997.pdf (2013).

[32] Tamura, K., Stecher, G., Peterson, D., Filipski, A. & Kumar, S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725-2729 (2013).

[33] Andrews, S. FastQC: a quality control for high throughout sequence data. http://www.bioinformatics.babraham.ac.uk/projects/fastqc (2010).

[34] Sommer, D. D., Delcher, A. L., Salzberg, S. L. & Pop, M. Minimus: a fast, lightweight genome assembler. BMC Bioinformatics 8, 1-64 (2007).

[35] Shao, W. et al. Analysis of 454 sequencing error rate, error sources, and artifact recombination for detection of low-frequency drug resistance mutations in HIV-1 DNA. Retrovirology 10, 1-16 (2013).

[36] Nock, C. J. et al. Chloroplast genome sequences from total DNA for plant identification. Plant Biotechnol. J. 9, 328-333 (2011).

[37] Zhang, T., Zhang, X., Hu, S. & Yu, J. An efficient procedure for plant organellar genome assembly, based on whole genome data from the 454 GS FLX sequencing platform. Plant Methods 7, 1-8 (2011).

[38] Xuelin, W. et al. The whole genome assembly and comparative genomic research of Thellungiella parvula (Extremophile crucifer) mitochondrion. Int. J. Genomics 2016, 5283628 (2016).

[39] Wang, D. et al. Transfer of chloroplast genomic DNA to mitochondrial genome occurred at least 300 MYA. Mol. Biol. Evol. 24, 2040-2048 (2007).

[40] Ye, N. et al. Assembly and comparative analysis of complete mitochondrial genome sequence of an economic plant Salix suchowensis. Peer J. 5, e3148 (2017).

[41] Barbrook, A. C., Howe, C. J., Kurniawan, D. P. & Tarr, S. J. Organization and expression of organellar genomes. Philos. Trans. R. Soc. B Biol. Sci. 365, 785-797 (2010).

[42] Ohyama, K. et al. Gene content, organization and molecular evolution of plant organellar genomes and sex chromosomes: insights from the case of the liverwort Marchantia polymorpha. Proc. Jpn. Acad. 85, 108-124 (2009).

[43] Wei, S. et al. Assembly and analysis of the complete Salix purpurea L. (Salicaceae) mitochondrial genome sequence. Springerplus 5, 1-10 (2016).

[44] Wang, W. & Messing, J. High-throughput sequencing of three Lemnoideae (duckweeds) chloroplast genomes from total DNA. PLoS ONE 6, e24670 (2011).

[45] Knoop V., Volkmar U., Hecht J., & Grewe F. Mitochondrial Genome Evolution in the Plant Lineage 3-29 (Springer, New York, 2011).

[46] Etminan, M., Fitzgerald, J. M., Gleave, M. & Chambers, K . Recombination and the maintenance of plant organelle genome stability. N. Phytol. 186, 299-317 (2010).

[47] Alverson, A. J., Rice, D. W., Dickinson, S., Barry, K. & Palmer, J. D. Origins and recombination of the bacterial-sized multichromosomal mitochondrial genome of cucumber. Plant Cell 23, 2499-2513 (2011).

[48] Chaubey, A. & Rajam, M. V. in Plant Biology and Biotechnology (eds Bahadur B., Venkat Rajam M., Sahijram L., Krishnamurthy K.) 179-204 (Springer, New Delhi, 2015).

[49] Peace, C. P. DNA-informed breeding of rosaceous crops: promises, progress and prospects. Hortic. Res. 4, 17006 (2017).

[50] Nederbragt, A. J. On the middle ground between open source and commercial software-the case of the Newbler program. Genome Biol. 15, 1-2 (2014).

[51] Greene, C. S. & Troyanskaya, O. G. Accurate evaluation and analysis of functional genomics data and methods. Ann. N. Y. Acad. Sci. 1260, 95-100 (2012).