- Open Access
Siberian larch (Larix sibirica Ledeb.) chloroplast genome and development of polymorphic chloroplast markers
BMC Bioinformaticsvolume 20, Article number: 38 (2019)
The main objectives of this study were sequencing, assembling, and annotation of chloroplast genome of one of the main Siberian boreal forest tree conifer species Siberian larch (Larix sibirica Ledeb.) and detection of polymorphic genetic markers – microsatellite loci or simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs).
We used the data of the whole genome sequencing of three Siberian larch trees from different regions - the Urals, Krasnoyarsk, and Khakassia, respectively. Sequence reads were obtained using the Illumina HiSeq2000 in the Laboratory of Forest Genomics at the Genome Research and Education Center of the Siberian Federal University. The assembling was done using the Bowtie2 mapping program and the SPAdes genomic assembler. The genome annotation was performed using the RAST service. We used the GMATo program for the SSRs search, and the Bowtie2 and UGENE programs for the SNPs detection. Length of the assembled chloroplast genome was 122,561 bp, which is similar to 122,474 bp in the closely related European larch (Larix decidua Mill.). As a result of annotation and comparison of the data with the existing data available only for three larch species - L. decidua, L. potaninii var. chinensis (complete genome 122,492 bp), and L. occidentalis (partial genome of 119,680 bp), we identified 110 genes, 34 of which represented tRNA, 4 rRNA, and 72 protein-coding genes. In total, 13 SNPs were detected; two of them were in the tRNA-Arg and Cell division protein FtsH genes, respectively. In addition, 23 SSR loci were identified.
The complete chloroplast genome sequence was obtained for Siberian larch for the first time. The reference complete chloroplast genomes, such as one described here, would greatly help in the chloroplast resequencing and search for additional genetic markers using population samples. The results of this research will be useful for further phylogenetic and gene flow studies in conifers.
The chloroplast genome in conifers, including larch species , has a unique, strictly paternal inheritance via pollen, unlike angiosperms, where it has a maternal inheritance via seeds . It allows tracing paternal gene flow and lineages separately from maternal (mitochondrial genes) and bi-parental (nuclear genes) ones. Therefore, chloroplast DNA sequences are the most important source of genetic markers to study distribution of paternal genes and paternally based molecular phylogenetic relationships in conifers.
Larch species, as well as many other conifer species are the main boreal forest tree species, which comprise ~ 30% of the world’s forested lands . Boreal forests play a very important ecological role, but are also affected by the global climate change. On one hand, they suffer now from more frequent and drastic droughts, but on the other hand, their area is expanding in the northern regions, and their tree line is moving towards the north creating an ecotone, a highly dynamic transition area . It is important to know how much of paternal associated gene flow by pollen contributes into establishing this zone, compared to the maternal and bi-parental contributions by seeds. Such studies require chloroplast markers. The next generation sequencing (NGS) technique allows whole chloroplast genome sequencing in multiple individuals and makes a search for molecular genetic markers more efficient. For instance, Parks et al.  sequenced chloroplast genomes in 37 pine species using NGS nearly completely. They found a significant amount of variation (especially in two loci ycf1 and ycf2) that provided them with additional data for inferring intrageneric phylogeny of genus Pinus.
Whole chloroplast genome comparison across different species and genera allows also studying organelle evolution and how it is associated with speciation and dispersal. Complete chloroplast genome sequences are available in NCBI Genbank for multiple plant species, including conifers. However, most of them represent the Pinus genus, and only three chloroplast genomes are available for the Larix genus: complete - for European (Larix decidua Mill.; AB501189.1) and Chinese (L. potaninii var. chinensis Beissn.; KX808508) larch and partial - for Western larch (L. occidentalis Nutt.; FJ899578.1).
Variation in the chloroplast genome is effectively used in phylogenetics at different levels. It allowed discriminating different subgenera and genera. For instance, Cronn et al.  compared chloroplast genome sequences of seven pine and one spruce species and found three regions that have deletions corresponding to the subgenera specific deletions in three genes: ycf12 (78 bp at the nucleotide starting position 51,051), psaM (93 bp at position 51,442), and ndhI (371 bp at position 101,988), respectively. These are common deletions in the chloroplast genome in pine species of the subgenus Strobus (i.e., P. gerardiana, P. krempfii, P. lambertiana, P. longaeva, P. monophylla, P. nelsonii, P. koraiensis); the corresponding genes were present in the subgenus Pinus (P. contorta, P. ponderosa, P. thunbergii) and in spruce Picea sitchensis .
Variation in the chloroplast genome can be also effectively used in discriminating different populations of the same species. For instance, Whittall et al.  demonstrated a strong differentiation between the mainland and island populations of Torrey pine (Pinus torreyana) based on 5 SNPs found in the entire chloroplast genome of 120 Kbp.
We used the data of the whole genome sequencing of three Siberian larch trees generated by Illumina HiSeq2000 . DNA samples were isolated from needles and haploid callus of three Siberian larch trees, representing different regions in Russia – the Ural Mountains, Krasnoyarsk Region and Khakassia Republic, respectively. The Larix decidua Mill.  and L. occidentalis Nutt.  chloroplast genomes were used as a reference (NCBI Genbank accession numbers AB501189.1 and FJ899578.1, respectively). We did not use the chloroplast genome of L. potaninii  as a reference, because it was assembled by using the chloroplast genome of L. decidua (NC_016058; ) as a reference, but we used it in the comparative analysis. The paired-end (PE) and mate-pair (MP) libraries with fragment sizes of 400–500 bp (Ural and Krasnoyarsk trees) and 300–400 bp (Khakassia tree), respectively, were used for sequencing via 2 × 100 cycles by Illumina HiSeq2000.
The sequence reads were mapped to the reference chloroplast genomes using the Bowtie2 software , which is good for mapping short sequence reads to medium-sized and large genomes. This software implements an algorithm to derive the FM-index based on the Burrows-Wheeler Transform. The SPAdes genome assembler has been used to assemble the larch genome, which implements the De Bruijn graph approach . The Rapid Annotation service with Subsystem Technology (RAST) has been used for annotation .
The first step in our assembly procedure consisted of mapping short reads to the available chloroplast genome references of L. decidua and L. occidentalis using the Bowtie2 software. Then, the aligned reads were assembled by SPAdes. The obtained contigs were aligned again on the reference of L. decidua using BLAST. In the third step, the selected contigs were verified to get the “trusted” status. Then, the assembly was carried out using SPAdes. The final step of the assembly was scaffolding, which was done using the generated contigs and MP reads using the SSPACE program .
Considering the well-known fact that chloroplast organelle originated from cyanobacteria, and that, therefore, chloroplast genes are still very similar to the bacterial ones, the RAST service, which was designed for annotation of bacterial and archaeal genomes, was used for the larch genome annotation. The annotation obtained by the RAST contained both the confirmed known genes and the predicted genes, potentially coding hypothetical proteins. In order to clarify the roles of these hypothetical coding regions, our annotation was compared with the annotations of two closely related species, L. decidua and L. occidentalis, respectively. In addition, some fragments of the genome have been also selectively aligned with BLAST. The sites of hypothetical proteins confirmed by BLAST were identified and recorded.
The assembled chloroplast genome of L. sibirica has been deposited in the NCBI GenBank with the accession number NC_036811.1 and used as a reference to search for polymorphisms (SNPs and SSRs) among the three above mentioned trees. SNPs were searched using the Bowtie2 and UGENE  software (option Call Variants with SAMtools). First, the reads of the Urals and Khakassian trees were mapped to the finally assembled genome of the Krasnoyarsk tree. Then, the resulting sam-file together with the assembled genome was used by the UGENE program to search for SNPs. The SSR loci were searched using GMATo  with a threshold of minimum 6 repeats for di-, tri-, tetra-, penta-, and hexanucleotide motifs, and 10 minimum repeats for mononucleotide motifs.
The total length of the final Siberian larch chloroplast genome assembly was 122,560 bp, which is very close to 122,474 bp in the closely related European larch (Larix decidua). The annotation through comparison with the available data for L. decidua and L. occidentalis identified 110 genes, from which 34 represented tRNA genes, 4 rRNA, and 72 protein-coding genes. In three trees 13 SNPs were detected. Two of them were found in the coding regions of the tRNA-Arg and Cell division protein ycf2 genes.
We used the available software, such as Bowtie2, BLAST, and SPAdes to assemble the chloroplast genome using the reads generated in the whole genome sequencing of the Siberian larch project. We used SSPACE for scaffolding and the RAST service for annotation of the obtained chloroplast genome. We developed a procedure that allowed us to successfully extract chloroplast genome specific reads and then assemble and annotate the resulting sequences. We identified and verified 110 coding regions representing 38 RNA and 72 protein genes, which is equal to the number of genes in chloroplast sequences of L. decidua and L. potaninii and close to 105 genes in the partial chloroplast genome sequence of L. occidentalis. A gene map of the genome was generated using OGDRAW  and is presented in Fig. 1. The search for SNPs using UGENE revealed a relatively small number of SNPs (Fig. 2; Additional file 1), but it is only preliminary data based on a limited sample size. In addition, 23 SSR loci (16 with mono- and 7 with dinucleotide repeat motifs) were identified in the chloroplast genome (Additional file 2). No SSR loci with tri-, tetra-, penta-, and hexanucleotide repeats were found with the search parameters used.
The chloroplast genome variation in most plants is often limited due to a relatively low frequency of mutations in this organelle. For example, the mutation rate of the chloroplast genome in pines is approximately 0.2–0.4 × 10− 9 synonymous substitutions per nucleotide per year [18, 19]. However, with an average length of 120–160 Kbp and 130 genes, the chloroplast genomes are sufficiently large and complex, and include structural and point mutations that reflect population differentiation and evolutionary divergence .
Unlike angiosperms, the conifer chloroplast DNA (cpDNA) lacks large inverted repeats (IR), but contains dispersed repetitive DNA that is associated with structural rearrangements. In addition to large dispersed repeated sequences, the conifer cpDNA also possesses a number of small repeats. It contains variable numbers of tandem repeats of 124 to 150 bp in size, which are associated with the polymorphic rearranged region near trnK-psbA, where the psbA gene has been duplicated .
Most variation in the chloroplast genome is associated with microsatellite loci [21, 22]. However, these markers have a too high mutation rate that can lead to incorrect phylogenetic inferences [23,24,25]. SNPs could be better markers for phylogenetic inferences, and comparative complete chloroplast genome studies are needed to discover these markers. The reference complete chloroplast genomes, such as the one described here, would greatly help in chloroplast resequencing and search for SNPs using population samples.
The complete chloroplast genome sequence was obtained for Siberian larch for the first time. Annotation and comparison of the obtained data with data available only for two other larch species helped us identify and verify 110 coding regions representing 38 RNA and 72 protein genes. The total of 13 SNPs were detected; two of them were in the coding regions of the genome. The results of this research will be useful for further phylogenetic and gene flow studies in conifers.
Szmidt AE, Aldén T, Hällgren J-E. Paternal inheritance of chloroplast DNA in Larix. Plant Mol Biol. 1987;9:59–64 https://doi.org/10.1007/BF00017987.
Hipkins VD, Krutovskii KV, Strauss SH. Organelle genomes in conifers: structure, evolution, and diversity. Forest Genet. 1995;1:179–89.
Brandt JP, Flannigan MD, Maynard DG, Thompson ID, Volney WJA. An introduction to Canada’s boreal zone: ecosystem processes, health, sustainability, and environmental issues. Environm Rev. 2013;21:207–26 https://doi.org/10.1139/er-2013-0040.
Johnson JS, Gaddis KD, Cairns DM, Konganti K, Krutovsky KV. Landscape genomic insights into the historic migration of mountain hemlock in response to Holocene climate change. Am J Bot. 2017;104(3):439–50 https://doi.org/10.3732/ajb.1600262.
Parks M, Cronn R, Liston A. Increasing phylogenetic resolution at low taxonomic levels using massively parallel sequencing of chloroplast genomes. BMC Biol. 2009;7:84 https://doi.org/10.1186/1741-7007-7-84.
Cronn R, Liston A, Parks M, Gernandt DS, Shen R, Mockler T. Multiplex sequencing of plant chloroplast genomes using Solexa sequencing-by-synthesis technology. Nucleic Acids Res. 2008;36(19):e122 https://doi.org/10.1093/nar/gkn502.
Whittall JB, Syring J, Parks M, Buenrostro J, Dick C, Liston A, Cronn R. Finding a (pine) needle in a haystack: chloroplast genome sequence divergence in rare and widespread pines. Mol Ecol. 2010;19:100–14 https://doi.org/10.1111/j.1365-294X.2009.04474.x.
Kuzmin DA, Feranchuk SI, Sharov VV, Cybin AN, Makolov SV, Putintseva YA, Oreshkova NV, Krutovsky KV. Stepwise large genome assembly approach: A case of Siberian larch (Larix sibirica Ledeb.). BMC Bioinformatics. 2019;20(Suppl. 1). https://doi.org/10.1186/s12859-018-2570-y.
Wu CS, Lin CP, Hsu CY, Wang RJ, Chaw SM. Comparative chloroplast genomes of Pinaceae: insights into the mechanism of diversified genomic organizations. Genome Biol Evol. 2011;3:309–19 https://doi.org/10.1093/gbe/evr026.
Han K, Li J, Zeng S, Liu Z. Complete chloroplast genome sequence of Chinese larch (Larix potaninii var. chinensis), an endangered conifer endemic to China. Conservation Genet Res. 2017;9:111–3 https://doi.org/10.1007/s12686-016-0633-9.
Langmead B, Salzberg S. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9 https://doi.org/10.1038/nmeth.1923.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77 https://doi.org/10.1089/cmb.2012.0021.
Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 2014;42(Database issue):D206–14 https://doi.org/10.1093/nar/gkt1226.
Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011;27:578–9 https://doi.org/10.1093/bioinformatics/btq683.
Okonechnikov K, Golosova O, Fursov M. The UGENE team. Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics. 2012;28:1166–7 https://doi.org/10.1093/bioinformatics/bts091.
Wang X, Lu P, Luo Z. GMATo: a novel tool for the identification and analysis of microsatellites in large genomes. Bioinformation. 2013;9:541–4.
Lohse M, Drechsel O, Kahlau S, Bock R. OrganellarGenomeDRAW – a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucleic Acids Res. 2013;41:575–81 https://doi.org/10.1093/nar/gkt289.
Willyard A, Syring J, Gernandt DS, Liston A, Cronn R. Fossil calibration of molecular divergence infers a moderate mutation rate and recent radiations for Pinus. Mol Biol Evol. 2007;24:90–101 https://doi.org/10.1093/molbev/msl131.
Gernandt DS, Magallon S, Lopez GG, Flores OZ, Willyard A, Liston A. Use of simultaneous analyses to guide fossil-based calibrations of Pinaceae phylogeny. Int J Plant Sci. 2008;169:1086–99 https://doi.org/10.1086/590472.
Lidholm J. Gustafsson P. a three-step model for the rearrangement of the chloroplast trnK-psbA region of the gymnosperm Pinus contorta. Nucleic Acids Res. 1991;19:2881–7.
Provan J, Soranzo N, Wilson NJ, Goldstein DB. Powell W a. a low mutation rate for chloroplast microsatellites. Genetics. 1999;153:943–7.
Ebert D, Peakall R. Chloroplast simple sequence repeats (cpSSRs): technical resources and recommendations for expanding cpSSR discovery and applications to a wide array of plant species. Mol Ecol Res. 2009;9:673–90 https://doi.org/10.1111/j.1755-0998.2008.02319.x.
Afzal-Rafii Z, Dodd A. Chloroplast DNA supports a hypothesis of glacial refugia over postglacial recolonization in disjunct populations of black pine (Pinus nigra) in Western Europe. Mol Ecol. 2007;16:723–36 https://doi.org/10.1111/j.1365-294X.2006.03183.x.
Höhn M, Gugerli F, Abran P, Bisztray G, Buonamici A, Cseke K, et al. Variation in the chloroplast DNA of Swiss stone pine (Pinus cembra L.) reflects contrasting post-glacial history of populations from the Carpathians and the Alps. J Biogeogr. 2009;36:1798–806 https://doi.org/10.1111/j.1365-2699.2009.02122.x.
Moreno-Letelier A, Piñero D. Phylogeographic structure of Pinus strobiformis Engelm. Across the Chihuahuan Desert filter-barrier. J Biogeogr. 2009;36:121–31 https://doi.org/10.1111/j.1365-2699.2008.02001.x.
We thank Dr. M. Sadovsky for discussing results and his valuable comments, Dr. I. N. Tretyakova and M. E. Pak for providing with callus culture, Dr. V. L. Semerikov and K. Deich for help with sampling, Dr. A. Ibe – with laboratory analysis, and Dr. D. Kuzmin and S. Makolov - with computing. This study was supported by a research grant No. 14.Y26.31.0004 from the Government of the Russian Federation.
The presented study was a part of the project “Genomic studies of major boreal coniferous forest tree species and their most dangerous pathogens in the Russian Federation” funded by the Government of the Russian Federation (grant № 14.Y26.31.0004). Publication costs are funded by the BioMed Central Membership of the University of Göttingen.
Availability of data and materials
The annotated chloroplast genome of L. sibirica has been deposited in the NCBI GenBank with the accession number NC_036811.1.
About this supplement
This article has been published as part of BMC Bioinformatics Volume 20 Supplement 1, 2019: Selected articles from BGRS\SB-2018: bioinformatics. The full contents of the supplement are available online at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-20-supplement-1.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.