Next Article in Journal
Domesticating Commercially Important Native Tree Species in the Philippines: Early Growth Performance Level
Next Article in Special Issue
Profiling of Widely Targeted Metabolomics for the Identification of Secondary Metabolites in Heartwood and Sapwood of the Red-Heart Chinese Fir (Cunninghamia Lanceolata)
Previous Article in Journal
Diplodia fraxini and Diplodia subglobosa: The Main Species Associated with Cankers and Dieback of Fraxinus excelsior in North-Eastern Italy
Previous Article in Special Issue
Physiological Characterization and Transcriptome Analysis of Camellia oleifera Abel. during Leaf Senescence
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 11
Issue 8
10.3390/f11080884
forests-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Complete Chloroplast Genome of Japanese Larch (Larix kaempferi): Insights into Intraspecific Variation with an Isolated Northern Limit Population

by
Shufen Chen
1,
Wataru Ishizuka
2,
Toshihiko Hara
3 and
Susumu Goto
1,*
1
Education and Research Center, The University of Tokyo Forests, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-8657, Japan
2
Forestry Research Institute, Hokkaido Research Organization, Koushunai, Bibai, Hokkaido 079-0166, Japan
3
Institute of Low Temperature Science, Hokkaido University, Sapporo-city, Hokkaido 060-0819, Japan
*
Author to whom correspondence should be addressed.
Forests 2020, 11(8), 884; https://doi.org/10.3390/f11080884
Submission received: 25 July 2020 / Revised: 10 August 2020 / Accepted: 11 August 2020 / Published: 14 August 2020
(This article belongs to the Special Issue Genetics and Improvement of Forest Trees)
Browse Figures
Review Reports Versions Notes

Abstract

:
Research Highlights: The complete chloroplast genome for eight individuals of Japanese larch, including from the isolated population at the northern limit of the range (Manokami larch), revealed that Japanese larch forms a monophyletic group, within which Manokami larch can be phylogenetically placed in Japanese larch. We detected intraspecific variation for possible candidate cpDNA markers in Japanese larch. Background and Objectives: The natural distribution of Japanese larch is limited to the mountainous range in the central part of Honshu Island, Japan, with an isolated northern limit population (Manokami larch). In this study, we determined the phylogenetic position of Manokami larch within Japanese larch, characterized the chloroplast genome of Japanese larch, detected intraspecific variation, and determined candidate cpDNA markers. Materials and Methods: The complete genome sequence was determined for eight individuals, including Manokami larch, in this study. The genetic position of the northern limit population was evaluated using phylogenetic analysis. The chloroplast genome of Japanese larch was characterized by comparison with eight individuals. Furthermore, intraspecific variations were extracted to find candidate cpDNA markers. Results: The phylogenetic tree showed that Japanese larch forms a monophyletic group, within which Manokami larch can be phylogenetically placed, based on the complete chloroplast genome, with a bootstrap value of 100%. The value of nucleotide diversity (π) was calculated at 0.00004, based on SNP sites for Japanese larch, suggesting that sequences had low variation. However, we found three hyper-polymorphic regions within the cpDNA. Finally, we detected 31 intraspecific variations, including 19 single nucleotide polymorphisms, 8 simple sequence repeats, and 4 insertions or deletions. Conclusions: Using a distant genotype in a northern limit population (Manokami larch), we detected sufficient intraspecific variation for the possible candidates of cpDNA markers in Japanese larch.

1. Introduction

The chloroplast genome is highly conserved and has a much lower mutation rate than the nuclear genome [1]. Chloroplast DNA (cpDNA) has been widely used to clarify interspecific relationships, and to evaluate the magnitude of intraspecific variation [2,3]. The cpDNA of gymnosperms, particularly of the conifers, is characterized by high levels of intraspecific variation [4,5] and paternal inheritance [6]. A high-resolution chloroplast-specific polymorphic assay would facilitate the analysis of population differentiation and gene flow in gymnosperms [7].
The next-generation sequencing (NGS) technique enables the sequencing of whole chloroplast genomes. The chloroplast genome has a circular molecular structure, with a length ranging from 120 to 160 kbp in most plants. The cpDNA contains a pair of inverted repeats (IRs), a large single-copy region (LSC), and a small single-copy region (SSC) [8]. IRs are a crucial feature of the chloroplast genome in most plants, likely contributing to the maintenance of a conserved arrangement of cpDNA sequences. Previous studies have reported that the length of observed short IRs is roughly consistent among gymnosperm species [9]. The whole chloroplast genome is of significant use for phylogenetic studies [2,3]; Parks et al. [10] presented complete chloroplast genomes for 37 pine species and documented a notable degree of variation at several loci (particularly at ycf1 and ycf2). Intraspecific variation in whole chloroplast genomes derived from multiple individuals can clarify the phylogenic lineage of target individuals [11]. In particular, single nucleotide polymorphisms (SNPs) have been efficiently used in the fields of phylogeography and conservation biology [12].
Japanese larch (Larix kaempferi (Lamb.) Carr.) is a deciduous coniferous tree species endemic to Japan, and integral to the country’s forestry efforts. The natural distribution of Japanese larch is limited to the mountainous range in the central part of Honshu Island, Japan [13]. An isolated population with ten mature trees was discovered at Manokami (hereafter Manokami larch), in the Zao Mountains in 1932 [14], extending the known northern limit of the species. Manokami larch was initially believed to be Larix gmelinii var. japonica, based on the morphological traits. An analysis of partial cpDNA sequences, and random amplified polymorphic DNA (RAPD) analysis, indicated that the Manokami larch population was actually Japanese larch [15]. However, the phylogenic position of Manokami larch, with relation to Japanese larch, has not yet been sufficiently defined [16,17].
The complete chloroplast genome of the genus Larix has been reported for several species [11,18,19], and the complete chloroplast genome of Japanese larch introduced in Korea, was reported by Kim et al. [20]. However, the intraspecific variations of Japanese larch have not yet been examined based on the complete chloroplast genome.
In this study, we identified the complete chloroplast genome for eight individuals of Japanese larch, including from the isolated population at the northern limit of the range (Manokami larch), to (1) determine the phylogenetic position of Manokami larch within Japanese larch, (2) characterize the chloroplast genome of Japanese larch, with included chloroplast data from Manokami larch, and (3) detect intraspecific variation and determine candidate cpDNA markers.

2. Materials and Methods

Eight individuals of Japanese larch were used in this study. Five individuals (Lk_Ho1, Lk_Ho2, Lk_Ho3, Lk_Ho4, and Sorachi 3) were collected from test plantations or the Arboretum Garden at the Forestry Research Institute, Hokkaido Research Organization (HRO). Sorachi 3 was selected specifically due to its superior growth in a larch-breeding program in Hokkaido. Two individuals (Lk_Ka1 and Lk_Ka2), were collected from the open-pollinated progeny of artificial plantations of Japanese larch, at the southern edge of Sakhalin Island, Russia. The detailed location of the seed collection of Lk_Ka2 was described in a previous study [11]. These seven individuals were originally derived from the eastern region of Nagano, in the central part of Honshu Island in Japan, near the center of distribution for this species. We added one grafted tree (Manokami 15), as an isolated germplasm of Manokami larch. The ortet of this tree was conserved in situ with the label “No.15” on Mt. Manokami in the Zao Mountains.
Fresh leaves from all eight individuals were collected in 2016 March for Lk_Ka2, June for Lk_Ho2, Lk_Ho3, Lk_Ho4, Sorachi3, 2017 March for Lk_Ka1, June for Lk_Ho1, and 2018 July for Manokami15. Leaf sampling, isolation of purified intact chloroplasts, and extraction of high-concentrate chloroplast DNA were performed as previously described [11]. Briefly, we used a saline Percoll (GE Healthcare, Uppsala, Sweden) gradient for chloroplast isolation and the DNeasy Plant Mini kit (QIAGEN, Hilden, Germany) for DNA extraction.
The cpDNA sequence reads were obtained using the Illumina platform. CLC Genomics Workbench 9.5.3 software (CLC bio, Aarhus, Denmark) was used for genetic analysis. After trimming low-quality sequences from the reads, bulked reads for all eight individuals were used to determine the draft consensus sequence for L. kaempferi. Reference mapping to L. gmelinii var. japonica (LC228570; [11]) was performed with parameter settings of mismatch cost 3, In/Del cost 3, length fraction 0.9, and similarity fraction 0.9. The complete chloroplast genome for each sample was then determined by mapping reads for each sample to our consensus sequence, using the same parameter settings as described above. The initial annotation of the chloroplast genome was performed using DOGMA [21]. Prediction of tRNA genes was performed using tRNAscan (http://lowelab.ucsc.edu/tRNAscan-SE). The annotation was finalized, with reference to that of L. gmelinii var. japonica (LC228570). To estimate the pseudogenes of ndh (subunits of an NADH dehydrogenase), we referred to the Pinus thunbergii chloroplast genome (NC_001631; [9]). REPuter [22] was used to confirm repetitive sequences in the chloroplast genome, (i.e., tandem repeats, duplicated genes, and IR regions). Finally, the gene map of the circular chloroplast genome of Japanese larch was drawn using OrganellarGenomeDRAW [23].
A phylogenetic tree was constructed based on the chloroplast genome sequences of the eight individuals identified in this study, and of MF990369 in the NCBI database (https://www.ncbi.nlm.nih.gov/) for Japanese larch, as well as five reference sequences of related Larix species derived from the database: LC228570 (L. gmelinii var. japonica), MF990370 (L. gmelinii var. olgensis), NC_016058 (L. decidua), KX880508 (L. potaninii), and NC_036811 (L. sibirica). The alignment of these fourteen chloroplast genome sequences was performed in MAFFT [24], and the final alignment was checked using CLC Genomics Workbench 9.5.3. A phylogenetic tree was constructed by MEGA X [25], based on maximum likelihood (ML) methods. A total of 1000 bootstrap replicates were applied to evaluate the branch supports.
The SNP data from the eight Japanese larch individuals was used for subsequent analyses. Haplotype networks have been demonstrated to show alternative genealogical relationships at the intraspecific population level, with low divergence [26]. We estimated haplotype networks for the chloroplast data using the software Network 10 (https://www.fluxus-engineering.com/sharenet.htm). Nucleotide diversity can be used as an inference parameter for evolutionary and demographic forces [12]; here, nucleotide diversity was calculated using DnaSP v6 software [27] and estimated as π. Divergent regions of the chloroplast genomes were identified according to the variation in π, by sliding window analysis, with a 500 bp step size and 10,000 bp window length.
Tandem repeats were identified using the Tandem Repeats Finder website [28]. In addition, simple sequence repeats (SSRs) were identified by MISA-web (https://webblast.ipk-gatersleben.de/misa/) [29], with minimum repetition numbers of 10 for mononucleotides, 6 for dinucleotides, and 5 for trinucleotides, tetranucleotides, pentanucleotides, and hexanucleotides each.

3. Results

3.1. Phylogenetic Analysis

The phylogenetic tree showed that Japanese larch forms a monophyletic group, within which Manokami larch can be phylogenetically placed based on the complete chloroplast genome, with a bootstrap value of 100% (Figure 1). Japanese larch is genetically close to L. decidua and L. gmelinii, but distant from L. sibirica and L. potaninii (Figure 1). The haplotype network among the eight sampled individuals revealed that Manokami larch was genetically distinct from other Japanese larches (Figure 2).

3.2. Characteristics of the Japanese Larch Chloroplast Genome

Japanese larch circular chloroplast genomes were characterized in the range of 122,394–122,409 bp with accession numbers from the DNA Data Bank of Japan (DDBJ) from LC574969 to LC574976. The gene type, number, and order were identical among the Japanese larch chloroplast genomes used in this study. Lk_Ho1 (LC574969) was used as a representative of Japanese larch; Figure 3 illustrates the physical mapping of its chloroplast genome, which contained a pair of IRs (436 bp each) separated by large single copy (LSC), and small single copy (SSC) regions, of 65,398 bp and 56,136 bp, respectively. The trnI-CAU gene was duplicated within inverted repeats, the trnS-GCU and psbl genes were duplicated as another inverted repeat of 457 bp in the LSC region (Figure 3), and two trnT-GGU genes were dispersed in the LSC and SSC regions, respectively. A total of 119 genes were identified (Table S1), including 72 protein genes, 35 transfer RNA genes, 4 ribosomal RNA genes, and 8 pseudogenes. Thirteen genes contained an intron, including trnA-UGC, trnG-UCC, trnI-GAU, trnK-UUU, trnL-UAA, trnV-UAC, rpoC1, rps12, rpl2, rpl16, petB, petD, and atpF. In addition, ycf3 contained two introns. Furthermore, rps12 was a trans-splicing gene with 5′ end and 3′ end exons, located in the LSC region and the SSC region, respectively. The G+C content of the complete chloroplast genome of Japanese larch was 38.7%.

3.3. Nucleotide Diversity Analysis

The value of nucleotide diversity (π) was calculated at 0.00004, based on SNP sites for Japanese larch, suggesting that sequences had low variation. As shown in Figure 4, there were three divergent regions (A, B, and C) in Japanese larch. Two regions (A, C), which were roughly in the range of the rpl16 gene psaB and the rbcL gene psbA, respectively, were classified as moderately variable (π > 0.00004); these regions contained variant sites in psaB, rpl16, ψndhK, atpB, psbK, and matK, and three intergenic spacers (between the rpl23 and psbA-partial genes, between the trnS-GGA and ycf3 genes, and between the trnS-GCU and trnT-GGU genes). The B region (roughly from chlL to rpl32, π > 0.0001) was identified as a hypervariable region, in which mutation occurred twice in the ψndhD and the ycf1.

3.4. Repeat Sequence Analysis of the Japanese Larch Chloroplast Genome

Tandem repeats were detected in approximately 25 sites in the Japanese larch chloroplast genomes. Repeated lengths of tandem repeats varied from 12 to 117 bp, and 64% of all tandem repeats occurred in ycf1, which belongs to a protein-coding region containing 76% of all detected tandem repeats. Nineteen SSR motifs were detected in the Japanese larch chloroplast genome. The majority of the detected SSR motifs were mononucleotide motifs, of which the SSR motif of mononucleotide T was the most frequent, followed by mononucleotide A and mononucleotide G. With the exception of three SSR motifs of dinucleotide AT, no other multiple nucleotide motifs were detected in the Japanese larch chloroplast genome. Furthermore, most (77.8%) of the detected SSR motifs were found in the intergenic region, followed by introns (16.7%), and protein-coding genes (5.5%). Eight cpSSR variants out of nineteen cpSSR motifs were detected in the intergenic region of the chloroplast genome of Japanese larch.

3.5. Genetic Variation among Japanese Larch Chloroplast Genomes

Among the eight individuals sequenced in this study, 31 variants (including 19 SNPs, 8 SSRs, and 4 In/Dels) were detected. For SNP variants, six and thirteen SNPs were identified in the intergenic spacer (IGS) and coding sequence (CDS) regions, respectively. These were detected in the ψndhK (one SNP), ψndhD (two SNPs), and protein-coding ycf1 (two SNPs), all of which belong to the SSC region. Six SSR variants were identified in the IGS, whereas two SSR variants were detected in the CDS region. Four In/Del variants were identified in the ψndhK gene (one In/Del variant) and the ycf1 gene (three In/Del variants), belonging to CDS region. (Table 1)

4. Discussion

The phylogenetic position of Manokami larch, has been discussed by several researchers [14,16,17]. This study clearly indicates that Manokami larch should be phylogenetically categorized into Japanese larch, with a bootstrap value of 100% (Figure 1). Our findings support the assertion by Shiraishi et al. [15] that Manokami larch must be a Japanese larch. Manokami larch is located far from other Japanese larches (Figure 2); genetically divergent genotypes, such as that of Manokami larch, could be used to efficiently detect intraspecific variation in Japanese larch.
The chloroplast genomes of Japanese larch obtained from this study were similar in size and gene order to those of L. gmelinii [11], L. sibirica [19], L. decidua [30], and L. potaninii [18]. The chloroplast structure types were classified in Pinaceae according to their alignment order and the orientation of the F1 (fragment flanked by trnG-UCC and trnE-UUC), F2 (fragment flanked by clpP and trnT-GGU), T1 (type 1 Pinaceae-specific repeat containing trnS-GCU and psbI), and T2 (type 2 Pinaceae-specific repeats in intergenic spacers) fragments in the LSC region, which can produce eight different cpDNA forms, including A, B, C, D, E, F, and G forms [30]. The chloroplast DNA form used in this study was classified into the C form, the same form identified for L. gmelinii, L. decidua, L. griffithiana, and Pinus elliottii [11,30] based on the alignment order and orientation of T1, T2, −F1 (reverse strand), T2, +F2 (forward strand), and T1. Due to this T1 repeat, there were longer inverted repeats (457 bp) in the LSC region than two IRs (436 bp) in Japanese larch. Extremely shortened IRs, with another pair of inverted repeats, is regarded as a common feature in Pinaceae [30,31].
In this study, three hotspots of variation were detected throughout the entire chloroplast genome (Figure 3 and Figure 4). The ycf1 and ψndhD were included in the hypervariable region (region B), and the ψndhK was included in the moderately variable region (region A). Three In/Del variants occurred in the ycf1, and previous research has reported insertions or deletions in the ycf1 of L. gmelinii [11]. Although it was considered a possibility that the ycf1 might be a nonfunctional pseudogene, another study [32] indicated that ycf1 is a functional gene, and encodes a product essential for cell survival. Dong et al. [33] revealed that the divergence of the ycf1 was obvious in gymnosperms. Additionally, Firetti et al. [34] indicated that the ycf1 was more divergent than the non-coding regions in the genus Anemopaegma. Regarding pseudogenes, eleven ndh genes (ndhA—K) have been identified in the cpDNA sequences of photosynthetic land plants [9,35,36]. In our study, five ψndh genes were found only in Japanese larch, of which ψndhD (two SNPs) and ψndhK (one SNP, one SSR, one In/Del variant) belonged to the region of frequent variation; these genes did not, however, exhibit a function consistent with other Pinus and Larix species [9,11].
Repeat sequences may play an important role in chloroplast genome arrangement and sequence divergence. In particular, tandem repeats may induce In/Dels [37,38]. In this study, tandem repeats were primarily identified in the ycf1. Tandem repeats were also located in the ycf1 of other conifers, such as Cryptomeria japonica [39] and L. gmelinii [11].
Among the eight Japanese larch individuals, we detected 31 variants (19 SNPs, 8 SSRs, and 4 In/Dels) located in psaB, ψndhD, ψndhK, psbE, psbK, rpoC1, rpoC2, the intron of rpl16, matK, atpB, ycf1, six intergenic spacers (between the rpl23 and psbA-partial genes, between the trnS-GGA and ycf3 genes, between the trnS-GCU and trnT-GGU genes, between the clpP and trnE-UUC genes twice, between the psbE and petL genes) with SNPs, the intron of atpF, ψndhK, six intergenic spacers (between the trnC-GCA and rpoB genes, between the ycf1 and rps15 genes, between the trnL-CAA and ycf2 genes, between the ψycf2 and trnV-GAC genes, between the trnT-GGU and trnV-UAC genes twice) with SSRs, ψndhK and ycf1 with In/Dels that could prove useful for providing candidate cpDNA markers. Chloroplast simple sequence repeat (cpSSR) markers often contain highly polymorphic variations within a population of conifers (see [7]), although Zhang et al. [40] found only three polymorphic cpSSR markers among 11 candidate markers in Japanese larch. We identified 19 SSR motifs within the chloroplast genome of Japanese larch, preferentially within the intergenic space, and only 8 SSR motifs occurred among 19 SSR motifs in the intergenic region of the chloroplast genome. These results lay a foundation for the development of cpDNA markers for Japanese larch.

5. Conclusions

The complete chloroplast genome of Japanese larch (122,398–122,409 bp) was obtained using next-generation sequencing technology. The comparison of whole chloroplast genomes clearly indicated that the isolated population, forming the northern limit of the species’ range (Manokami larch), should be placed phylogenetically within Japanese larch. The Manokami larch was found to be genetically different from other Japanese larches, indicating that sufficient genetic variation should be detected within the samples used in this study. Based on an analysis of intraspecific variation, 31 variants were detected, including 19 SNPs, 8 SSRs, and 4 In/Dels, all of which can be applied for the development of cpDNA markers. These variations should be useful for paternity analysis and population genetics analysis of Japanese larch in future studies.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4907/11/8/884/s1, Table S1: List of estimated chloroplast genes of Japanese larch. Genes with * have intron(s). Duplicated genes are shown in parenthesis.

Author Contributions

Conceptualization, W.I. and T.H.; Founding acquisition, T.H.; Software, W.I. and S.C.; visualization, W.I; Validation, W.I, S.C., T.H. and S.G.; Formal analysis, S.C. and W.I.; data curation, S.C. and W.I.; Writing—original draft preparation, S.C., S.G.; Writing—review and editing, S.C., W.I., T.H. and S.G. All authors have read and agreed to the published version of the manuscript.

Funding

This study was partly supported by the Japan Society for the Promotion of Science [15K18715] and the Grant for Joint Research Program of the Institute of Low Temperature Science, Hokkaido University.

Acknowledgments

We would like to thank K. Ono for the help of laboratory work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wolfe, K.H.; Li, W.-H.; Sharp, P.M. Rates of nucleotide substitution vary greatly among plant mitochondrial, chloroplast, and nuclear DNAs. Proc. Natl. Acad. Sci. USA 1987, 84, 9054–9058. [Google Scholar] [CrossRef] [Green Version]
  2. Wu, F.H.; Kan, D.P.; Lee, S.B.; Daniell, H.; Lee, Y.W.; Lin, C.C.; Lin, N.S.; Lin, C.S. Complete nucleotide sequence of Dendrocalamus latiflorus and Bambusa oldhamii chloroplast genomes. Tree Physiol. 2009, 29, 847–856. [Google Scholar] [CrossRef] [Green Version]
  3. Nock, C.J.; Waters, D.L.E.; Edwards, M.A.; Bowen, S.G.; Rice, N.; Cordeiro, G.M.; Henry, R.J. Chloroplast genome sequences from total DNA for plant identification. Plant Biotechnol. J. 2011, 9, 328–333. [Google Scholar] [CrossRef]
  4. Wagner, D.B.; Furnier, G.R.; Saghai-Maroof, M.A.; Williams, S.M.; Dancik, B.P.; Allard, R.W. Chloroplast DNA polymorphisms in lodgepole and jack pines and their hybrids. Proc. Natl. Acad. Sci. USA 1987, 84, 2097. [Google Scholar] [CrossRef] [Green Version]
  5. Tsumura, Y.; Suyama, Y.; Taguchi, H.; Ohba, K. Geographical cline of chloroplast DNA variation in Abies mariesii. Theor. Appl. Genet. 1994, 89, 922–926. [Google Scholar] [CrossRef]
  6. Neale, D.B.; Sederoff, R.R. Paternal inheritance of chloroplast DNA and maternal inheritance of mitochondrial DNA in loblolly pine. Theor. Appl. Genet. 1989, 77, 212–216. [Google Scholar] [CrossRef]
  7. Vendramin, G.G.; Lelli, L.; Rossi, P.; Morgante, M. A set of primers for the amplification of 20 chloroplast microsatellites in Pinaceae. Mol. Ecol. 1996, 5, 595–598. [Google Scholar] [CrossRef]
  8. Sugiura, M. The chloroplast chromosomes in land plants. Annu. Rev. Cell Biol. 1989, 5, 51–70. [Google Scholar] [CrossRef]
  9. Wakasugi, T.; Tsudzuki, J.; Ito, S.; Nakashima, K.; Tsudzuki, T.; Sugiura, M. Loss of all ndh genes as determined by sequencing the entire chloroplast genome of the black pine Pinus thunbergii. Proc. Natl. Acad. Sci. USA 1994, 91, 9794–9798. [Google Scholar] [CrossRef] [Green Version]
  10. 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. [Google Scholar] [CrossRef] [Green Version]
  11. Ishizuka, W.; Tabata, A.; Ono, K.; Fukuda, Y.; Hara, T. Draft chloroplast genome of Larix gmelinii var. japonica: Insight into intraspecific divergence. J. For. Res. 2017, 22, 393–398. [Google Scholar] [CrossRef]
  12. Branca, A.; Paape, T.D.; Zhou, P.; Briskine, R.; Farmer, A.D.; Mudge, J.; Bharti, A.K.; Woodward, J.E.; May, G.D.; Gentzbittel, L.; et al. Whole-genome nucleotide diversity, recombination, and linkage disequilibrium in the model legume. Proc. Natl. Acad. Sci. USA 2011, 108, E864. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Ngamitsu, T.; Nagasaka, K.; Yoshimaru, H.; Tsumura, Y. Provenance tests for survival and growth of 50-year-old Japanese larch (Larix kaempferi) trees related to climatic conditions in central Japan. Tree Genet. Genomes 2014, 10, 87–99. [Google Scholar] [CrossRef]
  14. Japan Forest Tree Breeding Association. Characteristics and conservation of isolated and northernmost natural population of larch in Japan-observing the population and trying to reveal what it really is like. In Proceedings of the 1st Forest Genetics Seminar, Miyagi, Japan, July 1994; p. 97. NAID: 10015777891. (In Japanese). [Google Scholar]
  15. Shiraishi, S.; Isoda, K.; Watanabe, A.; Kawasaki, H. DNA systematical study on the Larix relicted at Mt.Manokami, the Zao Mountains. J. Jpn. For. Soc. 1996, 78, 175–182. [Google Scholar]
  16. Oda, H. The northern natural Japanese larch forest at Manokami, the Zao Mountaions. For. Sci. 2013, 38, 52–58. (In Japanese) [Google Scholar]
  17. Tozawa, T. Ecological studies on the natural forest of the Japanese larch in Tohoku districts (Northeast of Honshu Island, Japan), 1: Natural forest of the now as its northern limit in Japan. Univ. For. Rep. Iwate Univ. 1975, 6, 1–18. [Google Scholar]
  18. Han, K.; Li, J.; Zeng, S.; Liu, Z.-L. Complete chloroplast genome sequence of Chinese larch (Larix potaninii var. chinensis), an endangered conifer endemic to China. Conserv. Genet. Resour. 2017, 9, 111–113. [Google Scholar]
  19. Bondar, E.I.; Putintseva, Y.A.; Oreshkova, N.V.; Krutovsky, K.V. Siberian larch (Larix sibirica Ledeb.) chloroplast genome and development of polymorphic chloroplast markers. BMC Bioinform. 2019, 20, 38. [Google Scholar] [CrossRef]
  20. Kim, S.-C.; Lee, J.-W.; Lee, M.-W.; Baek, S.-H.; Hong, K.-N. The complete chloroplast genome sequences of Larix kaempferi and Larix olgensis var. koreana (Pinaceae). Mitochondrial DNA Part B 2018, 3, 36–37. [Google Scholar] [CrossRef] [Green Version]
  21. Wyman, S.K.; Jansen, R.K.; Boore, J.L. Automatic annotation of organellar genomes with DOGMA. Bioinformatics 2004, 20, 3252–3255. [Google Scholar] [CrossRef] [Green Version]
  22. Kurtz, S.; Choudhuri, J.V.; Ohlebusch, E.; Schleiermacher, C.; Stoye, J.; Giegerich, R. REPuter: The manifold applications of repeat analysis on a genomic scale. Nucleic Acids Res. 2001, 29, 4633–4642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. 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, W575–W581. [Google Scholar] [CrossRef] [PubMed]
  24. Kuraku, S.; Zmasek, C.M.; Nishimura, O.; Katoh, K. aLeaves facilitates on-demand exploration of metazoan gene family trees on MAFFT sequence alignment server with enhanced interactivity. Nucleic Acids Res. 2013, 41, W22–W28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  26. Templeton, A.R.; Crandall, K.A.; Sing, C.F. A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram Estim. Genet. 1992, 132, 619. [Google Scholar]
  27. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  28. Benson, G. Tandem repeats finder: A program to analyze DNA sequences. Nucleic Acids Res. 1999, 27, 573–580. [Google Scholar] [CrossRef] [Green Version]
  29. Beier, S.; Thiel, T.; Münch, T.; Scholz, U.; Mascher, M. MISA-web: A web server for microsatellite prediction. Bioinformatics 2017, 33, 2583–2585. [Google Scholar] [CrossRef] [Green Version]
  30. Wu, C.-S.; Lin, C.-P.; Hsu, C.-Y.; Wang, R.-J.; Chaw, S.-M. Comparative chloroplast genomes of Pinaceae: Insights into the mechanism of diversified genomic organizations. Genome Biol. Evol. 2011, 3, 309–319. [Google Scholar] [CrossRef]
  31. Kang, H.-I.; Lee, H.O.; Lee, I.H.; Kim, I.S.; Lee, S.-W.; Yang, T.J.; Shim, D. Complete chloroplast genome of Pinus densiflora Siebold & Zucc. and comparative analysis with five pine trees. Forests 2019, 10, 600. [Google Scholar]
  32. Drescher, A.; Ruf, S.; Calsa, T., Jr.; Carrer, H.; Bock, R. The two largest chloroplast genome-encoded open reading frames of higher plants are essential genes. Plant J. 2000, 22, 97–104. [Google Scholar] [CrossRef]
  33. Dong, W.; Xu, C.; Li, C.; Sun, J.; Zuo, Y.; Shi, S.; Cheng, T.; Guo, J.; Zhou, S. ycf1, the most promising plastid DNA barcode of land plants. Sci. Rep. 2015, 5, 8348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Firetti, F.; Zuntini, A.R.; Gaiarsa, J.W.; Oliveira, R.S.; Lohmann, L.G.; Van Sluys, M.-A. Complete chloroplast genome sequences contribute to plant species delimitation: A case study of the Anemopaegma species complex. Am. J. Bot. 2017, 104, 1493–1509. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Ohyama, K.; Fukuzawa, H.; Kohchi, T.; Shirai, H.; Sano, T.; Sano, S.; Umesono, K.; Shiki, Y.; Takeuchi, M.; Chang, Z. Chloroplast gene organization deduced from complete sequence of liverwort Marchantia polymorpha chloroplast DNA. Nature 1986, 322, 572–574. [Google Scholar] [CrossRef]
  36. Hiratsuka, J.; Shimada, H.; Whittier, R.; Ishibashi, T.; Sakamoto, M.; Mori, M.; Kondo, C.; Honji, Y.; Sun, C.-R.; Meng, B.-Y.; et al. The complete sequence of the rice (Oryza sativa) chloroplast genome: Intermolecular recombination between distinct tRNA genes accounts for a major plastid DNA inversion during the evolution of the cereals. Mol. Gen. Genet. 1989, 217, 185–194. [Google Scholar] [CrossRef] [PubMed]
  37. Guisinger, M.M.; Kuehl, J.V.; Boore, J.L.; Jansen, R.K. Extreme reconfiguration of plastid genomes in the Angiosperm family Geraniaceae: Rearrangements, repeats, and codon usage. Mol. Biol. Evol. 2011, 28, 583–600. [Google Scholar] [CrossRef] [Green Version]
  38. Weng, M.-L.; Blazier, J.C.; Govindu, M.; Jansen, R.K. Reconstruction of the ancestral plastid genome in Geraniaceae reveals a correlation between genome rearrangements, repeats, and nucleotide substitution rates. Mol. Biol. Evol. 1994, 31, 645–659. [Google Scholar] [CrossRef] [Green Version]
  39. Hirao, T.; Watanabe, A.; Kurita, M.; Kondo, T.; Takata, K. A frameshift mutation of the chloroplast matK coding region is associated with chlorophyll deficiency in the Cryptomeria japonica virescent mutant Wogon-Sugi. Curr. Genet. 2009, 55, 311–321. [Google Scholar] [CrossRef] [Green Version]
  40. Zhang, X.Y.; Shiraishi, S.; Huang, M.R. Analysis of genetic structure in population of Larix Kaempferi by chloroplast SSR markers. Hereditas 2004, 26, 486–490. [Google Scholar]
Forests 11 00884 g001 550
Figure 1. The maximum likelihood (ML) phylogenetic tree based on 14 chloroplast genomes of Larix species. The red square represents Japanese larch.
Figure 1. The maximum likelihood (ML) phylogenetic tree based on 14 chloroplast genomes of Larix species. The red square represents Japanese larch.
Forests 11 00884 g001
Forests 11 00884 g002 550
Figure 2. Haplotype network based on SNP sites within eight individuals of Japanese larch. The red dot represents missing haplotypes.
Figure 2. Haplotype network based on SNP sites within eight individuals of Japanese larch. The red dot represents missing haplotypes.
Forests 11 00884 g002
Forests 11 00884 g003 550
Figure 3. A gene map of Japanese larch chloroplast genome (accession number: LC574969). Genes shown outside and inside the circle are transcribed clockwise, and transcribed counterclockwise, respectively. Genes were colour-coded to distinguish different functional groups. The dark and light gray inner circle indicates the GC and AT content of the chloroplast genome, respectively. “†” represents the location of a longer inverted repeat. A, B and C represent hotspot of variation.
Figure 3. A gene map of Japanese larch chloroplast genome (accession number: LC574969). Genes shown outside and inside the circle are transcribed clockwise, and transcribed counterclockwise, respectively. Genes were colour-coded to distinguish different functional groups. The dark and light gray inner circle indicates the GC and AT content of the chloroplast genome, respectively. “†” represents the location of a longer inverted repeat. A, B and C represent hotspot of variation.
Forests 11 00884 g003
Forests 11 00884 g004 550
Figure 4. The nucleotide diversity (π) of Japanese larch chloroplast genomes, based on sliding window analysis.
Figure 4. The nucleotide diversity (π) of Japanese larch chloroplast genomes, based on sliding window analysis.
Forests 11 00884 g004
Table
Table 1. Genetic variation among chloroplast genomes of Japanese larch. CDS: coding sequence; IGS: intergenic spacer; In/Del: insertion or deletion; SNP: single nucleotide polymorphism; SSR: simple sequence repeat.
Table 1. Genetic variation among chloroplast genomes of Japanese larch. CDS: coding sequence; IGS: intergenic spacer; In/Del: insertion or deletion; SNP: single nucleotide polymorphism; SSR: simple sequence repeat.
SequenceNIn/DelSNPSSR
CDSIGSCDSIGSCDSIGS
Lk_Ho1(Reference)———————
Lk_Ho22 2
Lk_Ho331 11
Lk_Ho48 42 2
Lk_Ka14 31
Lk_Ka251 3 1
Sorachi33 3
Manokami15122 2314
Whole sequences314198

Share and Cite

MDPI and ACS Style

Chen, S.; Ishizuka, W.; Hara, T.; Goto, S. Complete Chloroplast Genome of Japanese Larch (Larix kaempferi): Insights into Intraspecific Variation with an Isolated Northern Limit Population. Forests 2020, 11, 884. https://doi.org/10.3390/f11080884

AMA Style

Chen S, Ishizuka W, Hara T, Goto S. Complete Chloroplast Genome of Japanese Larch (Larix kaempferi): Insights into Intraspecific Variation with an Isolated Northern Limit Population. Forests. 2020; 11(8):884. https://doi.org/10.3390/f11080884

Chicago/Turabian Style

Chen, Shufen, Wataru Ishizuka, Toshihiko Hara, and Susumu Goto. 2020. "Complete Chloroplast Genome of Japanese Larch (Larix kaempferi): Insights into Intraspecific Variation with an Isolated Northern Limit Population" Forests 11, no. 8: 884. https://doi.org/10.3390/f11080884

APA Style

Chen, S., Ishizuka, W., Hara, T., & Goto, S. (2020). Complete Chloroplast Genome of Japanese Larch (Larix kaempferi): Insights into Intraspecific Variation with an Isolated Northern Limit Population. Forests, 11(8), 884. https://doi.org/10.3390/f11080884

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop