Next Article in Journal
Hydrogen Peroxide-Induced Root Ca2+ and K+ Fluxes Correlate with Salt Tolerance in Cereals: Towards the Cell-Based Phenotyping
Next Article in Special Issue
The Complete Plastome Sequence of an Antarctic Bryophyte Sanionia uncinata (Hedw.) Loeske
Previous Article in Journal
Anti-Hypertensive Effects of Acacia Polyphenol in Spontaneously Hypertensive Rats
Previous Article in Special Issue
Exploring the History of Chloroplast Capture in Arabis Using Whole Chloroplast Genome Sequencing
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 3
10.3390/ijms19030701
ijms-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 Sequences of Four Meliaceae Species and Comparative Analyses

by
Malte Mader
1,
Birte Pakull
1,
Céline Blanc-Jolivet
1,
Maike Paulini-Drewes
1,
Zoéwindé Henri-Noël Bouda
1,
Bernd Degen
1,
Ian Small
2 and
Birgit Kersten
1,*
1
Thünen Institute of Forest Genetics, Sieker Landstrasse 2, D-22927 Grosshansdorf, Germany
2
Australian Research Centre of Excellence in Plant Energy Biology, School of Molecular Sciences, The University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(3), 701; https://doi.org/10.3390/ijms19030701
Submission received: 31 January 2018 / Revised: 26 February 2018 / Accepted: 26 February 2018 / Published: 1 March 2018
(This article belongs to the Special Issue Chloroplast)
Browse Figures
Versions Notes

Abstract

:
The Meliaceae family mainly consists of trees and shrubs with a pantropical distribution. In this study, the complete chloroplast genomes of four Meliaceae species were sequenced and compared with each other and with the previously published Azadirachta indica plastome. The five plastomes are circular and exhibit a quadripartite structure with high conservation of gene content and order. They include 130 genes encoding 85 proteins, 37 tRNAs and 8 rRNAs. Inverted repeat expansion resulted in a duplication of rps19 in the five Meliaceae species, which is consistent with that in many other Sapindales, but different from many other rosids. Compared to Azadirachta indica, the four newly sequenced Meliaceae individuals share several large deletions, which mainly contribute to the decreased genome sizes. A whole-plastome phylogeny supports previous findings that the four species form a monophyletic sister clade to Azadirachta indica within the Meliaceae. SNPs and indels identified in all complete Meliaceae plastomes might be suitable targets for the future development of genetic markers at different taxonomic levels. The extended analysis of SNPs in the matK gene led to the identification of four potential Meliaceae-specific SNPs as a basis for future validation and marker development.

Graphical Abstract

1. Introduction

The Meliaceae or mahogany family is a flowering plant family of mainly trees and shrubs (and a few mangroves and herbaceous plants) in the order Sapindales. The species of the Meliaceae family, which are contained in The Plant List [1] belong to 52 plant genera, all showing a pantropical distribution. The Plant List includes 3198 scientific plant names of species rank for the family Meliaceae. Of these, 669 are accepted species names.
Cedrela odorata L., a fast-growing deciduous tree species, is the most commercially important and widely distributed species in the genus Cedrela, and one of the world’s most important timber species. It is found from Mexico southwards throughout central America to northern Argentina, as well as in the Caribbean. The aromatic wood is in high demand in the American tropics because it is naturally termite- and rot-resistant. It contains an aromatic and insect-repelling resin that is the source of one of its popular names, Spanish-cedar (it resembles the aroma of true cedars, Cedrus spp.). Other common names include Cuban cedar or cedro in Spanish. It is used for a wide range of purposes, including the production of furniture and craft items as well as different medicinal uses of the bark. It is an excellent choice for use in reforestation, because it is considered a pioneer species.
Entandrophragma cylindricum (Sprague) Sprague, commonly known as sapele or sapelli, is a large tree native to tropical Africa. The commercially important wood is reminiscent of mahogany (Swietenia macrophylla), a member of the same family. Demand for sapelli as a mahogany substitute, often traded as “African mahogany”, has increased sharply in recent years. It is sold both in lumber and veneer form. Among other applications is its use in musical instruments.
Khaya senegalensis (Desv.) A. Juss. represents a Khaya species of the African riparian woodlands and savanna zone. Common names include African mahogany, dry zone mahogany, Gambia mahogany, or Senegal mahogany, among others. The wood is used for a variety of purposes, such as carpentry, interior trim, and construction. The bitter tasting bark is utilized for a variety of medical purposes.
Carapa guianensis Aubl. is a tree from the tropical regions of Southern Central America, the Amazon region, and the Caribbean. The wood resembles mahogany and is used in quality furniture. The seed oil is used in traditional medicine and as an insect repellent.
All four Meliaceae species analyzed in this study (Cedrela odorata, Entandrophragma cylindricum, Khaya senegalensis, Carapa guianensis) are listed as vulnerable on the IUCN red list of threatened species [2]. Furthermore, Cedrela odorata is one of the six Meliaceae species that are included in the list of CITES-protected species to ensure that international trade will not threaten the survival of this species [3]. Detecting violations of CITES regulations in the tropical timber trade requires accurate genus and species identification of wood, which is often difficult or even impossible using anatomical methods, especially if the wood is processed. Thus, there is a high demand for the development of genetic markers for these purposes.
For genus and species identification, chloroplast DNA (cpDNA) markers are often useful. In contrast to nuclear genomes, chloroplast genomes are inherited uniparentally (maternally in most seed plants) [4], show a dense gene content and a slower evolutionary rate of change (e.g., [5]). The double-stranded cpDNA is present in many copies per cell, leading to the convenient situation that cpDNA can be retrieved relatively easily from low-quantity and/or degraded DNA samples, including wood samples. These advantages of cpDNA are reflected in the recommendation of the Barcode of Life Consortia to apply molecular markers, mainly based on organelle DNA, to genetically differentiate all eukaryotic species [6]. For the genetic differentiation of vascular plants, molecular markers which are mainly based on DNA variations in two chloroplast regions (rbcL and matK) are used as a two-locus barcode [7].
Chloroplast-derived DNA sequences have been widely used for taxonomic purposes and phylogenetic studies (e.g., [8,9]). Complete cp genome sequences provide valuable data sets to resolve complex evolutionary relationships [10] and have been shown to improve resolution at lower taxonomic levels (e.g., [11]). The application of Next-Generation Sequencing (NGS) technologies [12], especially using genome-skimming strategies [13,14] has made it relatively easy to obtain complete cpDNA sequences for low cost. The high abundance of cpDNA compared to nuclear DNA allows the use of total DNA for genome skimming without prior purification of chloroplasts or cpDNA [15,16].
The comparative analysis of complete plastomes allows the extension of two-locus barcoding to next-generation barcoding (whole-plastome barcoding), gene-based phylogenetics to genome-based phylogenomics, and the development of molecular markers for taxonomic and phylogeographic purposes [15,16,17,18].
In Meliaceae, the development of cpDNA markers has so far been restricted to specific loci [19,20]. Although the Meliaceae form a large family, only the plastome sequence of one species, Azadirachta indica, has been previously published in GenBank (KF986530.1). In this study, we sequenced the plastomes of four additional species and compared them to that of Azadirachta indica as well as those of other species in the order Sapindales to analyze interspecific variation within the Meliaceae and uncover Meliaceae-specific genome features.

2. Results

2.1. The Structure of Chloroplast Genomes from Four Meliaceae Species

The complete cp genomes of Cedrela odorata, Entandrophragma cylindricum, Khaya senegalensis and Carapa guianensis were sequenced using Illumina sequencing technology in a genome-skimming approach (Table 1) and compared to the previously published plastome of Azadirachta indica (GenBank KF986530.1).
The plastomes of these four Meliaceae species (and of Azadirachta indica) are small circular DNA molecules of sizes in the range of 158,558 bp to 160,737 bp, with the typical quadripartite structure of land plant cp genomes consisting of two inverted repeats (IRa and IRb) separated by large (LSC) and small (SSC) single copy regions, respectively (Table 2).
In each of the cp genomes of the four newly sequenced Meliaceae species (Table 1), 112 different genes were annotated (78 protein-coding, 30 tRNA and 4 rRNA genes), 18 of which were duplicated in the IR regions, giving a total of 130 genes encoding 85 proteins, 37 tRNAs and 8 rRNAs (Table 2 and Table 3). Among the 112 unique genes, 18 included one or two introns. The intron sizes are ranging from 532 bp for trnL-UAA to 2535 bp for trnK-UUU when considering the plastome of Cedrela odorata. One gene, rps12, is presumed to require trans-splicing (Table 3; [21]).
The gene maps of the newly sequenced Meliaceae species are provided in Figures 1 (Cedrela odorata), and in Figure S1–S3 (Entandrophragma cylindricum, Khaya senegalensis and Carapa guianensis). The gene content and gene order in the cp genomes of the four species are nearly identical to the previously published plastome of Azadirachta indica (GenBank KF986530.1) with the exception that the ycf1 gene at the IRa/SSC border was not annotated as a pseudogene in Azadirachta indica resulting in a total number of 131 genes in this species compared to 130 genes in the other four Meliaceae species (Table 2). Another exception is that the gene names of trnG-UCC and trnC-GCC are swapped due to misannotation in Azadirachta indica.
As mentioned above, 18 unique genes were annotated to include introns in the four newly sequenced Meliaceae species (Table 3), whereas introns are missing in the annotations of two of these genes in Azadirachta indica (GenBank KF986530.1), namely petD and rps12 (intron in the 3′ part of the gene). The annotations of these genes are probably not correct in Azadirachta indica; both exons are relatively short (exon 1 of petD: 8 bp; exon 3 of rps12: 26 bp), thus making their annotation difficult.
The gene rps19 is included in the inverted repeats (close to the IRa/LSC or IRb/LSC border, respectively; Figure 1, Figures S1–S3, Table 3) in the four Meliaceae species sequenced in this study as well as in Azadirachta indica (KF986530.1) thus resulting in a duplication of rps19.

2.2. Diversity of the Meliaceae Plastome Sequences

The newly sequenced individuals of the four different Meliaceae species share several large regions, located in intergenic regions and in the rpl16 intron (LSC), which show low similarity to Azadirachta indica (Figure 2). These regions are related to large deletions which become obvious in a multiple alignment of related whole plastomes (Figure S4). The largest deletion which is in the psbE-petL linker is of a size of about 199 bp (Figure S4; at position 69531 bp of the Azadirachta indica sequence). These deletions mainly contribute to the smaller genome sizes of the four individuals compared to Azadirachta indica (Table 2). The Cedrela odorata individual, which is the individual with the smallest cp genome size (158,558 bp; Table 2), shows, compared to the four other individuals, additional large deletions in different intergenic linkers, e.g., the ycf3-rps4 linker (Figure 2). A large exclusive deletion was also detected for the Carapa guianensis individual in the psbM-psbD linker region in the LSC (Figure 2, Figure S4).
Based on the multiple whole-plastome alignment of the five Meliaceae individuals (Figure S4), 7635 positions that showed DNA sequence variations (SNPs or indels) in at least one of the plastomes (compared to the consensus sequence), were called by the SNiPlay tool (SNiPlay genotyping table in Table S5). The following regions of the consensus sequence showed the highest frequencies of variations (considering intervals of 10,000 bp): 1–10,000 bp (1–9681 bp in Azadirachta indica) with 923 variable positions (top1); 120,000–130,000 bp (117,660–127,364 bp in A. indica) with 771 positions (top2); and 130,000–140,000 bp (127,365–137,199 bp in A. indica) with 735 positions (top3). The top1-region is located at the 5-prime part of the LSC including psbA, matK, rps16, psbK and psbI (Figure 2). The top2- and top3-region are connected and represent the SSC downstream of ndhF, the SSC/IRb border and parts of the rRNA cluster of the IRb (Figure 2).
A phylogenetic tree based on a multiple alignment of the five complete Meliaceae cpDNA sequences (Table 2) and Acer buergerianum (family Sapindaceae in the Sapindales; NC_034744.1) as an outgroup, shows that the analyzed Meliaceae form a monophyletic sister clade to Azadirachta indica within the Meliaceae. Within this clade, Cedrela odorata and Entandrophragma cylindricum group together, as do Khaya senegalensis and Carapa guianensis (Figure 3).

2.3. Comparative Analyses for the Identification of Potential Meliaceae-Specific Plastome DNA Variations in One Barcoding Region

The matK gene, one of the cpDNA barcoding regions (see Introduction), was selected for comparative analyses aiming at the identification of potential Meliaceae-specific cpDNA variations. Based on a multiple alignment (Figure S5) of the extracted matK gene sequences of the five Meliaceae species listed in Table 2 together with five non-Meliaceae species of the order Sapindales (Boswellia sacra Flueck (Burseraceae), NC_029420.1; Anacardium occidentale L. (Anacardiaceae), NC_035235.1; Leitneria floridana Chapm. (Simaroubaceae), NC_030482.1; Citrus sinensis (L.) Osbeck (Rutaceae), NC_008334.1 [23]; Acer buergerianum Miq. (Aceraceae), NC_034744.1), 16 SNP positions were identified where the five Meliaceae individuals showed the same nucleotide, which, however, differed from the nucleotide(s) of the non-Meliaceae species at the same position (SNP positions summarized in Table 4).
These 16 SNP positions were further analyzed in a multiple alignment of matK gene sequences from 100 Meliaceae individuals downloaded from GenBank (Figure S6). Only SNPs where all the downloaded Meliaceae sequences showed the same nucleotide at the given position were further validated with more matK sequences from members of other taxonomic groups (other Sapindales, other Rosids, and other land plants; Table 4).
The SNPs at the following four positions of the matK gene were selected as potential Meliaceae-specific: 346 bp/328 bp (consensus sequence/Cedrela odorata sequence), 1318 bp/1270 bp, 1478 bp/1430 bp and 1494 bp/1446 bp (Table 4). At position 1318 bp/1270 bp, e.g., all Meliaceae individuals considered include a G (Table 4, Figure S6) in contrast to members of other Sapindales families, other Rosids, as well as other land plants, where the considered individuals show an A (Table 4). This SNP is at the first position of the codon encoding for MatK amino acid 424 in Cedrela odorata and will result in an amino acid exchange (Meliaceae analyzed: Gly (Figure S6); other Sapindales, Rosids and land plants analyzed: Ser or Asp).
The other barcoding region, the rbcL gene (see Introduction), showed less nucleotide variation compared to matK when analyzed in a similar way (Figure S7). In a multiple alignment of the rbcL sequences of five Meliaceae species and five non-Meliaceae species of the order Sapindales (Figure S7), only 2 SNP positions were identified where the five Meliaceae individuals showed another nucleotide than the non-Meliaceae species at the same position (SNP positions 225 and 582 related to the consensus sequence). However, at these two positions, some other Meliaceae individuals showed differing nucleotides at GenBank; thus, these positions were not further analyzed.

3. Discussion

In this study, the complete chloroplast genomes of four Meliaceae species (order Sapindales) were sequenced by NGS (genome-skimming), annotated (Figure 1, Figures S1–S3, Table 3) and compared with each other and with the only previously published Meliaceae plastome sequence, the sequence of Azadirachta indica (KF986530.1; Table 2). Across eight other families within the Sapindales, there are currently 38 complete cpDNA sequences available at the Organellar Genome Resource at NCBI [24].
Gene and intron content are highly conserved among land plant plastomes, although losses have been identified in several angiosperm lineages (reviewed, e.g., in [25,26]). The following genes known to be lost in some other Angiosperm species [25] are also missing in the four newly sequenced Meliaceae plastomes (Table 3) as well as in Azadirachta indica (KF986530.1): chlB, chlL, chlN, infA, trnP-GGG. The losses of chlB, chlL, chlN, and trnP-GGG represent synapomorphies for flowering plants [25]. The most common gene loss involves infA [25,27]. There is evidence that infA has been transferred to the nucleus in some species [27], which has not yet been investigated in Meliaceae. The low depth of coverage (related to the nuclear genome) of the genome skimming data used in this study did not allow to perform such an investigation.
Intron content is also highly conserved across angiosperms with most genomes containing 18 genes with introns [25]. In the four Meliaceae species sequenced in this study, 18 unique genes with one or two introns were annotated (Table 3). The apparent lack of the petD intron and the intron in the 3′-rps12 locus in Azadirachta indica (KF986530.1) are annotation errors in our opinion.
Hotspots for structural rearrangements within plastomes include the IRs, which are frequently subjected to expansion, contraction, or even complete loss (summarized in [21]). Inverted repeat expansion resulted in a duplication of rps19 in the five Meliaceae plastomes analyzed so far (Figure 1, Figures S1–S3, Table 3, KF986530.1). This is consistent with research considering the plastomes of Nitraria sibirica (Nitrariaceae [28]) and 38 other Sapindales species [24] not belonging to the Meliaceae. However, some of these plastomes have only incomplete second rps19 copies (Anacardium occidentale, NC_035235.1; Acer miaotaiense, NC_030343.1; Acer buergerianum, NC_034744.1; Mangifera indica, NC_035239.1) or the second copy is even completely missing, such as in Aesculus wangii (NC_035955.1) or Pistacia vera (NC_034998.1), where the first rps19 copy is in addition incomplete. The rps19 gene is completely absent in Rhus chinensis (NC_033535.1). In contrast to the Sapindales, many other rosid species include only one rps19 gene in the LSC of the cpDNA sequence [28]. To answer the question if the duplication of rps19 is a general plastome feature in Meliaceae, the plastomes of other member species of the 52 plant genera within the Meliaceae must be sequenced and annotated.
Complete plastome sequences are valuable for deciphering phylogenetic relationships especially between closely related taxa, or where recent divergence, rapid speciation or slow genome evolution has resulted in limited sequence variation [18,29,30]. Considering that most species-level diversity of Meliaceae in rainforests is recent [31], a whole-plastome phylogeny is highly desired for this family. The whole-plastome phylogeny that was generated in this study (Figure 3) based on the whole-plastome alignment for five Meliaceae species (Figure S4), is an initial step in this direction. The result that Azadirachta indica (member of the subfamily Melioidaeae) belongs to another subclade than the four other Meliaceae species (members of the subfamily Swietenioideae) as well as the further sub-grouping of the four Meliaceae species, agrees with previous studies (subfamilies according to [31,32,33,34]). In the future, more complete cpDNA sequences of member individuals of other Meliaceae species are needed to exploit the power of the whole-plastome phylogenies, especially for differentiation at lower taxonomic levels. Additional integration of complete cpDNA sequences with small amplicon datasets could further improve the phylogenetic resolution, as recently shown in Acacia, where the greatest support has been achieved when using a whole-plastome phylogeny as a constraint on the amplicon-derived phylogeny [30].
Whole-plastome alignments are also very useful to develop cpDNA markers for the genetic identification of species or other taxonomic categories [17]. Especially, large indels can be easily identified in whole-plastome alignments and used, e.g., for the development of robust PCR-based markers, as recently shown for a Populus tremula-specific marker that has been developed based on a 96 bp-indel [17,35]. In the present study, we identified and compactly visualized large indels between the analyzed Meliaceae individuals based on pairwise alignments using VISTA (Figure 2). Exclusive large deletions were identified for the Cedrela odorata and the Carapa guianensis individuals, respectively (Figure 2). Before effective markers can be developed, it must be further validated with more Meliaceae species and individuals whether these deletions are genus-, species- or individual-specific. PCR-marker development based on indels in small intergenic regions might be particularly robust, because the primers could be placed into conserved genic regions adjacent to the indel. The SNPs and indels identified between the five Meliaceae plastomes (Table S5), based on a multiple alignment, may also serve—after further validation—as targets for future marker development at different taxonomic levels.
Aiming at the identification of potential Meliaceae-specific SNPs, the matK gene, one of the barcoding regions [7]—was selected and the SNPs identified in this gene were analyzed in extended sets of previously published matK sequences. In general, nucleotide variation per site in matK is three times higher than in rbcL (the large subunit of RUBISCO; [36]), and the amino acid substitution rate is six times that of rbcL [37]. Here, four potential Meliaceae-specific SNPs were identified (Table 4), three of which are in the 3′-terminus of the matK gene encoding for the C-terminus of MatK in a region that is homolog to domain X of mitochondrial group II intron maturases [38] and has a significantly higher amount of basic amino acids compared to the N-terminal region [39]. The potential Meliaceae-specific SNPs—once further validated—should be attractive for the taxonomic differentiation of samples from wood or wood products.

4. Materials and Methods

4.1. Sampling, DNA Extraction and Sequencing

The individual trees analyzed are described in detail in Table 1. Genomic DNA was isolated from leaves or cambium according to Dumolin et al. 1995. Genomic library generation and sequencing on the Illumina MiSeq v3 (2 × 300 bp paired-end reads) was done by Eurofins Genomics (Ebersberg, Germany).

4.2. Assemblies of Chloroplast Genome Sequences and Annotation

If not otherwise stated, CLC Genomics Workbench (CLC-GWB; v8.5.1 and v9.5.3; CLC bio, A QIAGEN company; Aarhus, Denmark) was used for data processing. Initial quality control of the NGS reads was done with FastQC [40].

4.2.1. Assembly of cpDNA Sequences of Khaya Senegalensis and Entandrophragma Cylindricum

All reads were trimmed with CLC-GWB including adapter trimming, quality trimming (quality limit of 0.01), trimming of 10 nucleotides at the 5′-end and removing reads of less than 50 bp in length. All other options were set to default. Potential chloroplast reads were extracted by mapping all trimmed reads against the chloroplast sequence of Azadirachta indica (KF986530.1) using the tool MITObim [41]. Multiple de novo assemblies based on different word sizes were performed on extracted chloroplast reads. Overlapping contigs of one or more assemblies were used to combine the contigs to complete cpDNA sequences.

4.2.2. Assembly of cpDNA Sequences of Carapa Guianensis and Cedrela Odorata

Adapter sequences were removed by the trimming software Trimmomatic [42] using simple and palindromic trim. Further trimming was done by CLC-GWB including quality trimming (quality limit of 0.01), trimming of 10 nucleotides at the 3′- and 15 nucleotides 5′-end and removing reads shorter than 50 bp. All other options were set to default. A first set of contigs was generated by de novo assembly of all trimmed reads, using a length fraction of 0.9 and a similarity fraction of 0.95. All resulting contigs were blasted with the command line tool blastn [43] against the nucleotide BLAST database downloaded from GenBank [44] to identify and select chloroplast contigs. All resulting Blast hits were filtered for the keyword “chloroplast” and finally validated with the web Blast tool at NCBI [45]. Trimmed reads were mapped to the chloroplast contigs, mapped reads were extracted and stored in a separate read set as chloroplast reads. Multiple de novo assemblies based on different word sizes were performed on these chloroplast reads. Overlapping contigs of one or more assemblies were used to assemble the complete cpDNA sequence.

4.3. Annotation of the cpDNA Sequences and Preparation of GenBank Submission Files

Draft annotations were generated with the web-based software CPGAVAS [46,47] to check gene content and order. These draft annotations were corrected where necessary, guided by alignments to other well-characterized eudicot plastomes including those of Arabidopsis thaliana (NC_000932), Nicotiana tabacum (NC_001879) and Spinacia oleracea (NC_002202).
The file resulting from the fine annotation of the plastome of Khaya senegalensis (GB-format) was transferred to a draft SQN-file using the CHLOROBOX-GenBank2Sequin-tool [48] and edited using the Sequin tool (v13.05; [49]). The error-corrected SQN-file was submitted to GenBank. Because gene content and order were the same in all species analyzed (according to the fine annotation), the GenBank submission files (SQN-format) for the other species were created by updating the sqn-file of Khaya senegalensis with the sequences of the other three species (using the “update sequence” function in Sequin) and subsequent manual editing of shifted annotations in Sequin.
The circular gene maps of the four Meliaceae plastomes (Figure 1; Figures S1–S3) were obtained using the OrganellarGenomeDRAW software (OGDRAW v1.2; [22,50].

4.4. Alignments and Construction of a Phylogenetic Tree

Pairwise alignments of complete cpDNA sequences were performed using the VISTA tool mVISTA (“AVID” as alignment program) at VISTA [51,52]. Identity plots of each pairwise alignment were downloaded from the related VISTA-point results page.
Multiple alignments of complete cpDNA sequences were run using CLC-GWB (v. 8.5.1.; parameters: gap open cost = 10.0; gap extension cost = 1.0; end gap cost = as any other; alignment mode = very accurate; redo alignments = no; use fixpoints = no).
The phylogenetic tree was constructed based on the alignment of five Meliaceae species together with one outgroup (Acer buergerianum, NC_034744.1; family Sapindaceae in the order Sapindales) using the “Maximum likelihood phylogeny” tool of CLC-GWB including bootstrap analysis with 100 replicates (other parameters: construction method for the start tree = UPGMA; existing start tree = not set; nucleotide substitution model = Jukes Cantor; protein substitution model = WAG; transition/transversion ratio = 2.0; include rate variation = No; number of substitution rate categories = 4; gamma distribution parameter = 1.0; estimate substitution rate parameter(s) = Yes; estimate topology = Yes; estimate gamma distribution parameter = no).

4.5. SNP and Indel Detection in Multiple Alignments Using SNiPlay

To identify SNPs between the complete cpDNA sequences of the 5 Meliaceae species (Table 2), the alignment-FASTA of the related multiple alignment was exported from CLC-GWB and further edited (replacement of alignment spaces “-“—if present—by “?”at the 5′- or 3′- end of the sequences). The edited alignment-FASTA was used as an input for the web tool SNiPlay (pipeline v2; [53]) to run SNP/indel discovery (default parameters) [54].

4.6. NCBI-Blast Analyses of matK and Download of matK Gene Sequences of Different Taxonomic Groups for Multiple Alignments

The matK gene sequence of Cedrela odorata (MG724915) was used as a query in different BlastN searches (parameters: optimized for highly similar sequences, maximal target sequences: 100) at GenBank (NCBI; [45]): (i) restrict to Meliaceae (taxid:43707); (ii) restrict to Sapindales (taxid:41937) and exclude Meliaceae; (iii) restrict to Rosids (taxid:71275) and exclude Sapindales; and (iv) no restriction, but exclude Rosids. After each Blast analysis, all 100 hits were selected and a FASTA was downloaded with the aligned sequences. Each FASTA file was used as input for a multiple alignment together with the matK sequence of Cedrela odorata (used as reference). In the case of the Meliaceae alignment, sequences of wrong orientation in the alignment were reverse complemented and the alignment was repeated. In the case of other alignments only sequences in the right orientation (CDS orientation) were considered in the further analyses.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/19/3/701/s1. Figure S1. Gene map of the complete cpDNA sequence of Entandrophragma cylindricum (KY923074). Figure S2. Gene map of the complete cpDNA sequence of Khaya senegalensis (KX364458). Figure S3. Gene map of the complete cpDNA sequence of Carapa guianensis (MF401522). Figure S4. Whole-plastome alignment of five Meliaceae species. Figure S5. Alignment of the matK gene sequences of 5 Meliaceae species (Table 2) and 5 non-Meliaceae members of the order Sapindales (Boswellia sacra, NC_029420.1; Anacardium occidentale, NC_035235.1; Leitneria floridana, NC_030482.1; Citrus sinensis, NC_008334.1; Acer buergerianum, NC_034744.1). Figure S6. Alignment of the matK gene sequences of 100 Meliaceae individuals (GenBank) and Cedrela odorata (MG724915; reference). Figure S7. Alignment of the rbcL gene sequences of 5 Meliaceae species (Table 2) and 5 non-Meliaceae members of the order Sapindales (Boswellia sacra, NC_029420.1; Anacardium occidentale, NC_035235.1; Leitneria floridana, NC_030482.1; Citrus sinensis, NC_008334.1; Acer buergerianum, NC_034744.1). Table S1. SNPs and Indels identified by SNIplay in the whole-plastome alignment of 5 Meliaceae species.

Acknowledgments

This work was funded by the Federal Ministry of Food and Agriculture (BMEL), Germany in the scope of the “Large scale project on genetic timber verification” and by the International Tropical Timber Organization (project PD 620/11 Rev.1). We would like to thank Susanne Bein for technical assistance in the laboratory, Stephen Cavers and Niklas Tysklind for providing samples from South America. Open access costs were covered by the Thuenen Institute, Braunschweig, Germany.

Author Contributions

B.D. initiated the project; B.P., M.M., C.B.-J., B.D. and B.K. conceived and designed the experiments; B.P. and M.P.-D. performed the experiments; M.M., B.K. and I.S. analyzed the data; Z.H.-N.B. collected and contributed sample material; B.K., M.M. and I.S. wrote the manuscript; and all authors edited and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Abbreviations

cpchloroplast
cpDNAchloroplast DNA
CITESConvention on International Trade in Endangered Species of Wild Fauna and Flora
IRaInverted Repeat a
IRbInverted Repeat b
IUCNInternational Union for Conservation of Nature
LSCLarge Single-Copy Region
NGSNext-Generation Sequencing
SSCSmall Single-Copy Region

References

  1. The Plant List. Available online: http://www.theplantlist.org/ (accessed on 15 January 2018).
  2. The IUCN Red List of Threatened Species. Available online: http://www.iucnredlist.org/ (accessed on 11 January 2018).
  3. CITES—Convention on International Trade in Endangered Species of Wild Fauna and Flora/the Cites Species. Available online: https://cites.org/eng/disc/species.php (accessed on 11 January 2018).
  4. Birky, C.W. Uniparental inheritance of mitochondrial and chloroplast genes—Mechanisms and evolution. Proc. Natl. Acad. Sci. USA 1995, 92, 11331–11338. [Google Scholar] [CrossRef] [PubMed]
  5. Drouin, G.; Daoud, H.; Xia, J. Relative rates of synonymous substitutions in the mitochondrial, chloroplast and nuclear genomes of seed plants. Mol. Phylogenet. Evol. 2008, 49, 827–831. [Google Scholar] [CrossRef] [PubMed]
  6. Barcode of Life. Available online: http://www.barcodeoflife.org/ (accessed on 1 January 2018).
  7. Hollingsworth, P.M.; Forrest, L.L.; Spouge, J.L.; Hajibabaei, M.; Ratnasingham, S.; van der Bank, M.; Chase, M.W.; Cowan, R.S.; Erickson, D.L.; Fazekas, A.J.; et al. A DNA barcode for land plants. Proc. Natl. Acad. Sci. USA 2009, 106, 12794–12797. [Google Scholar] [CrossRef]
  8. Shaw, J.; Shafer, H.L.; Leonard, O.R.; Kovach, M.J.; Schorr, M.; Morris, A.B. Chloroplast DNA sequence utility for the lowest phylogenetic and phylogeographic inferences in angiosperms: The tortoise and the hare IV. Am. J. Bot. 2014, 101, 1987–2004. [Google Scholar] [CrossRef] [PubMed]
  9. Wang, W.; Li, H.L.; Chen, Z.D. Analysis of plastid and nuclear DNA data in plant phylogenetics-evaluation and improvement. Sci. China Life Sci. 2014, 57, 280–286. [Google Scholar] [CrossRef] [PubMed]
  10. Moore, M.J.; Soltis, P.S.; Bell, C.D.; Burleigh, J.G.; Soltis, D.E. Phylogenetic analysis of 83 plastid genes further resolves the early diversification of eudicots. Proc. Natl. Acad. Sci. USA 2010, 107, 4623–4628. [Google Scholar] [CrossRef] [PubMed]
  11. Carbonell-Caballero, J.; Alonso, R.; Ibanez, V.; Terol, J.; Talon, M.; Dopazo, J. A phylogenetic analysis of 34 chloroplast genomes elucidates the relationships between wild and domestic species within the genus Citrus. Mol. Biol. Evol. 2015, 32, 2015–2035. [Google Scholar] [CrossRef] [PubMed]
  12. Goodwin, S.; McPherson, J.D.; McCombie, W.R. Coming of age: Ten years of next-generation sequencing technologies. Nat. Rev. Genet. 2016, 17, 333–351. [Google Scholar] [CrossRef] [PubMed]
  13. Dodsworth, S. Genome skimming for next-generation biodiversity analysis. Trends Plant Sci. 2015, 20, 525–527. [Google Scholar] [CrossRef] [PubMed]
  14. Straub, S.C.K.; Parks, M.; Weitemier, K.; Fishbein, M.; Cronn, R.C.; Liston, A. Navigating the tip of the genomic iceberg: Next-generation sequencing for plant systematics. Am. J. Bot. 2012, 99, 349–364. [Google Scholar] [CrossRef] [PubMed]
  15. Pakull, B.; Mader, M.; Kersten, B.; Ekue, M.R.M.; Dipelet, U.G.B.; Paulini, M.; Bouda, Z.H.N.; Degen, B. Development of nuclear, chloroplast and mitochondrial SNP markers for Khaya sp. Conserv. Genet. Resour. 2016, 8, 283–297. [Google Scholar] [CrossRef]
  16. Schroeder, H.; Cronn, R.; Yanbaev, Y.; Jennings, T.; Mader, M.; Degen, B.; Kersten, B. Development of molecular markers for determining continental origin of wood from white oaks (Quercus l. Sect. Quercus). PLoS ONE 2016, 11, e0158221. [Google Scholar] [CrossRef] [PubMed]
  17. Kersten, B.; Rampant, P.F.; Mader, M.; Le Paslier, M.C.; Bounon, R.; Berard, A.; Vettori, C.; Schroeder, H.; Leple, J.C.; Fladung, M. Genome sequences of Populus tremula chloroplast and mitochondrion: Implications for holistic poplar breeding. PLoS ONE 2016, 11, e0147209. [Google Scholar] [CrossRef] [PubMed]
  18. Tonti-Filippini, J.; Nevill, P.G.; Dixon, K.; Small, I. What can we do with 1000 plastid genomes? Plant J. 2017, 90, 808–818. [Google Scholar] [CrossRef] [PubMed]
  19. Duminil, J.; Kenfack, D.; Viscosi, V.; Grumiau, L.; Hardy, O.J. Testing species delimitation in sympatric species complexes: The case of an african tropical tree, Carapa spp. (Meliaceae). Mol. Phylogenet. Evol. 2012, 62, 275–285. [Google Scholar] [CrossRef] [PubMed]
  20. Holtken, A.M.; Schroder, H.; Wischnewski, N.; Degen, B.; Magel, E.; Fladung, M. Development of DNA-based methods to identify cites-protected timber species: A case study in the Meliaceae family. Holzforschung 2012, 66, 97–104. [Google Scholar] [CrossRef]
  21. Wicke, S.; Schneeweiss, G.M.; dePamphilis, C.W.; Muller, K.F.; Quandt, D. The evolution of the plastid chromosome in land plants: Gene content, gene order, gene function. Plant Mol. Biol. 2011, 76, 273–297. [Google Scholar] [CrossRef] [PubMed]
  22. 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]
  23. Bausher, M.G.; Singh, N.D.; Lee, S.B.; Jansen, R.K.; Daniell, H. The complete chloroplast genome sequence of Citrus sinensis (L.) Osbeck var ‘ridge pineapple’: Organization and phylogenetic relationships to other angiosperms. BMC Plant Biol. 2006, 6, 21. [Google Scholar] [CrossRef] [PubMed]
  24. Organelle Resources at NCBI. Available online: http://www.ncbi.nlm.nih.gov/genome/organelle/ (accessed on 1 February 2018).
  25. Jansen, R.K.; Cai, Z.; Raubeson, L.A.; Daniell, H.; Depamphilis, C.W.; Leebens-Mack, J.; Muller, K.F.; Guisinger-Bellian, M.; Haberle, R.C.; Hansen, A.K.; et al. Analysis of 81 genes from 64 plastid genomes resolves relationships in angiosperms and identifies genome-scale evolutionary patterns. Proc. Natl. Acad. Sci. USA 2007, 104, 19369–19374. [Google Scholar] [CrossRef] [PubMed]
  26. Green, B.R. Chloroplast genomes of photosynthetic eukaryotes. Plant J. 2011, 66, 34–44. [Google Scholar] [CrossRef] [PubMed]
  27. Millen, R.S.; Olmstead, R.G.; Adams, K.L.; Palmer, J.D.; Lao, N.T.; Heggie, L.; Kavanagh, T.A.; Hibberd, J.M.; Giray, J.C.; Morden, C.W.; et al. Many parallel losses of infa from chloroplast DNA during angiosperm evolution with multiple independent transfers to the nucleus. Plant Cell 2001, 13, 645–658. [Google Scholar] [CrossRef] [PubMed]
  28. Lu, L.; Li, X.; Hao, Z.; Yang, L.; Zhang, J.; Peng, Y.; Xu, H.; Lu, Y.; Zhang, J.; Shi, J.; et al. Phylogenetic studies and comparative chloroplast genome analyses elucidate the basal position of halophyte Nitraria sibirica (Nitrariaceae) in the sapindales. Mitochondrial DNA Part A DNA Mapp. Seq. Anal. 2017, 1–11. [Google Scholar] [CrossRef] [PubMed]
  29. Daniell, H.; Lin, C.S.; Yu, M.; Chang, W.J. Chloroplast genomes: Diversity, evolution, and applications in genetic engineering. Genome Biol. 2016, 17, 134. [Google Scholar] [CrossRef] [PubMed]
  30. Williams, A.V.; Miller, J.T.; Small, I.; Nevill, P.G.; Boykin, L.M. Integration of complete chloroplast genome sequences with small amplicon datasets improves phylogenetic resolution in Acacia. Mol. Phylogenet. Evol. 2016, 96, 1–8. [Google Scholar] [CrossRef] [PubMed]
  31. Koenen, E.J.M.; Clarkson, J.J.; Pennington, T.D.; Chatrou, L.W. Recently evolved diversity and convergent radiations of rainforest mahoganies (Meliaceae) shed new light on the origins of rainforest hyperdiversity. New Phytol. 2015, 207, 327–339. [Google Scholar] [CrossRef] [PubMed]
  32. Pennington, T.D.; Styles, B.T. A generic monograph of the meliaceae. Blumea 1975, 22, 419–540. [Google Scholar]
  33. Muellner, A.N.; Samuel, R.; Johnson, S.A.; Cheek, M.; Pennington, T.D.; Chase, M.W. Molecular phylogenetics of Meliaceae (Sapindales) based on nuclear and plastid DNA sequences. Am. J. Bot. 2003, 90, 471–480. [Google Scholar] [CrossRef] [PubMed]
  34. Muellner, A.N.; Pennington, T.D.; Chase, M.W. Molecular phylogenetics of neotropical Cedreleae (mahogany family, Meliaceae) based on nuclear and plastid DNA sequences reveal multiple origins of “cedrela odorata”. Mol. Phylogenet. Evol. 2009, 52, 461–469. [Google Scholar] [CrossRef] [PubMed]
  35. Schroeder, H.; Kersten, B.; Fladung, M. Development of multiplexed marker sets to identify the most relevant poplar species for breeding. Forests 2017, 8, 492. [Google Scholar] [CrossRef]
  36. Soltis, D.E.; Soltis, P.S. Choosing an approach and an appropriate gene for phylogenetic analysis. In Molecular Systematics of Plants II; Soltis, D.E., Soltis, P.S., Doyle, J.J., Eds.; Kluwer Academic Publishers: Boston, MA, USA, 1998; pp. 1–42. [Google Scholar]
  37. Olmstead, R.G.; Palmer, J.D. Chloroplast DNA systematics—A review of methods and data-analysis. Am. J. Bot. 1994, 81, 1205–1224. [Google Scholar] [CrossRef]
  38. Neuhaus, H.; Link, G. The chloroplast tRNALys(UUU) gene from mustard (Sinapis alba) contains a class-II intron potentially coding for a maturase-related polypeptide. Curr. Genet. 1987, 11, 251–257. [Google Scholar] [CrossRef] [PubMed]
  39. Barthet, M.M.; Hilu, K.W. Evaluating evolutionary constraint on the rapidly evolving gene matK using protein composition. J. Mol. Evol. 2008, 66, 85–97. [Google Scholar] [CrossRef] [PubMed]
  40. Braham Bioinformatics/Fastqc. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 1 December 2017).
  41. Hahn, C.; Bachmann, L.; Chevreux, B. Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads-a baiting and iterative mapping approach. Nucleic Acids Res. 2013, 41. [Google Scholar] [CrossRef] [PubMed]
  42. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  43. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. Blast+: Architecture and applications. BMC Bioinf. 2009, 10, 421. [Google Scholar] [CrossRef] [PubMed]
  44. Genbank. Available online: http://www.ncbi.nlm.nih.gov/genbank/ (accessed on 1 January 2018).
  45. Ncbi/Blast. Available online: https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastn&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome (accessed on 19 January 2018).
  46. Liu, C.; Shi, L.C.; Zhu, Y.J.; Chen, H.M.; Zhang, J.H.; Lin, X.H.; Guan, X.J. Cpgavas, an integrated web server for the annotation, visualization, analysis, and genbank submission of completely sequenced chloroplast genome sequences. BMC Genom. 2012, 13, 715. [Google Scholar] [CrossRef] [PubMed]
  47. CpGAVAS: Chloroplast Genome Annotation, Visualization, Analysis and Genbank. Available online: http://www.herbalgenomics.org/0506/cpgavas (accessed on 1 February 2018).
  48. Chlorobox/genbank2sequin. Available online: https://chlorobox.mpimp-golm.mpg.de/GenBank2Sequin.html (accessed on 1 February 2018).
  49. NCBI/sequin. Available online: https://www.ncbi.nlm.nih.gov/Sequin/ (accessed on 1 February 2018).
  50. OrganellarGenomeDRAW. Available online: http://ogdraw.mpimp-golm.mpg.de/ (accessed on 1 February 2018).
  51. VISTA. Available online: http://genome.lbl.gov/vista/index.shtml (accessed on 1 February 2018).
  52. Mayor, C.; Brudno, M.; Schwartz, J.R.; Poliakov, A.; Rubin, E.M.; Frazer, K.A.; Pachter, L.S.; Dubchak, I. VISTA: Visualizing global DNA sequence alignments of arbitrary length. Bioinformatics 2000, 16, 1046–1047. [Google Scholar] [CrossRef] [PubMed]
  53. SNiPlay Pipeline v2. Available online: http://sniplay.cirad.fr/cgi-bin/analysis.cgi (accessed on 1 February 2018).
  54. Dereeper, A.; Homa, F.; Andres, G.; Sempere, G.; Sarah, G.; Hueber, Y.; Dufayard, J.F.; Ruiz, M. SNiPlay3: A web-based application for exploration and large scale analyses of genomic variations. Nucleic Acids Res. 2015, 43, W295–W300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 19 00701 g001 550
Figure 1. Gene map of the complete chloroplast genome of Cedrela odorata (GenBank MG724915). The grey arrows indicate the direction of transcription of the two DNA strands. A GC-content graph is depicted within the inner circle. The circle inside the GC content graph marks the 50% threshold. The maps were created using OrganellarGenomeDRAW [22].
Figure 1. Gene map of the complete chloroplast genome of Cedrela odorata (GenBank MG724915). The grey arrows indicate the direction of transcription of the two DNA strands. A GC-content graph is depicted within the inner circle. The circle inside the GC content graph marks the 50% threshold. The maps were created using OrganellarGenomeDRAW [22].
Ijms 19 00701 g001
Ijms 19 00701 g002 550
Figure 2. Visualization of pairwise alignments of complete cpDNA sequences of four Meliaceae species each with the cpDNA sequence of Azadirachta indica (reference). VISTA-based similarity plots portraying the sequence identity of each of the four Meliacea species with the reference Azadirachta indica are shown. The annotation (protein-encoding genes) is provided for Azadirachta indica on top (based on the related GenBank file; KF986530.1). Plastome regions with the highest diversity between the 5 Meliaceae individuals are marked by blue arrows (top1–3). Further details are provided in the text below.
Figure 2. Visualization of pairwise alignments of complete cpDNA sequences of four Meliaceae species each with the cpDNA sequence of Azadirachta indica (reference). VISTA-based similarity plots portraying the sequence identity of each of the four Meliacea species with the reference Azadirachta indica are shown. The annotation (protein-encoding genes) is provided for Azadirachta indica on top (based on the related GenBank file; KF986530.1). Plastome regions with the highest diversity between the 5 Meliaceae individuals are marked by blue arrows (top1–3). Further details are provided in the text below.
Ijms 19 00701 g002
Ijms 19 00701 g003 550
Figure 3. Phylogenetic tree (maximum likelihood) based on whole-plastome sequences of five Meliaceae species and Acer buergerianum (outgroup). Bootstrap values (%) are shown above branches. The bootstrap value on the branch separating Azadirachta indica from the other Meliaceae is below 70% and was not shown for this reason. GenBank accession numbers of the plastomes are given in Table 2.
Figure 3. Phylogenetic tree (maximum likelihood) based on whole-plastome sequences of five Meliaceae species and Acer buergerianum (outgroup). Bootstrap values (%) are shown above branches. The bootstrap value on the branch separating Azadirachta indica from the other Meliaceae is below 70% and was not shown for this reason. GenBank accession numbers of the plastomes are given in Table 2.
Ijms 19 00701 g003
Table
Table 1. Information on samples and NGS data of four Meliaceae individuals sequenced in this study.
Table 1. Information on samples and NGS data of four Meliaceae individuals sequenced in this study.
SpeciesIndividual (Thuenen-ID)Origin/LocationLongitudeLatitudeTrimmed ReadsCoverage *
Cedrela odorataCEODO_205_2Cuba, population Guisa−76.6820.16254214101×
Entandrophragma cylindricumc-5-ENTC-46Cameroon, FBR, Parc National de Lobeke15.64422.262862206300165×
Khaya senegalensisKSUnknown/Green house, Thünen Institute Hamburg-Lohbrügge 2783117346×
Carapa guianensisCAGUI_332_1French Guiana, region Rorota−52.2623924.87761422759135×
* The coverage is based on a final mapping of the trimmed reads to the assembled cpDNA sequence.
Table
Table 2. Summary of Meliaceae chloroplast genome features.
Table 2. Summary of Meliaceae chloroplast genome features.
Cedrela Odorata (MG724915)Entandrophragma Cylindricum (KY923074.1)Khaya Senegalensis (KX364458.1)Carapa Guianensis (MF401522.1)Azadirachta Indica * (KF986530.1)
Genome size (bp)158,558159,609159,787159,483160,737
LSC length (bp)86,39087,11787,40487,05488,137
SSC length (bp)18,38018,53218,31118,27718,636
IR length (bp)26,89426,98027,03627,07626,982
Number of genes130130130130131
* Not sequenced in this study. Identifiers (in parenthesis) under the species name refer to GenBank accession numbers.
Table
Table 3. List of genes annotated in the cp genomes of four Meliaceae species sequenced in this study (Table 2).
Table 3. List of genes annotated in the cp genomes of four Meliaceae species sequenced in this study (Table 2).
FunctionGenes
RNAs, ribosomalrrn23, rrn16, rrn5, rrn4.5
RNAs, transfertrnA-UGC *, trnC-GCA, trnD-GUC, trnE-UUC, trnF-GAA, trnG-GCC, trnG-UCC *, trnH-GUG, trnI-CAU, trnI-GAU *, trnK-UUU *, trnL-CAA, trnL-UAA *, trnL-UAG, trnM-CAU, trnfM-CAU, trnN-GUU, trnP-UGG, trnQ-UUG, trnR-ACG, trnR-UCU, trnS-GCU, trnS-GGA, trnS-UGA, trnT-GGU, trnT-UGU, trnV-GAC, trnV-UAC *, trnW-CCA, trnY-GUA
Transcription and splicingrpoA, rpoB, rpoC1 *, rpoC2, matK
Translation, ribosomal proteins
Small subunitrps2, rps3, rps4, rps7, rps8, rps11, rps12 **,T, rps14, rps15, rps16 *, rps18, rps19
Large subunitrpl2 *, rpl14, rpl16 *, rpl20, rpl22, rpl23, rpl32, rpl33, rpl36
Photosynthesis
ATP synthaseatpA, atpB, atpE, atpF *, atpH, atpI
Photosystem IpsaA, psaB, psaC, psaI, psaJ, ycf3 **, ycf4
Photosystem IIpsbA, psbB, psbC, psbD, psbE, psbF, psbH, psbI, psbJ, psbK, psbL, psbM, psbN, psbT, psbZ
Calvin cyclerbcL
Cytochrome complexpetA, petB *, petD *, petG, petL, petN
NADH dehydrogenaseNdhA *, ndhB *, ndhC, ndhD, ndhE, ndhF, ndhG, ndhH, ndhI, ndhJ, ndhK
OthersclpP1 **, accD, cemA, ccsA, ycf1, ycf2
* Genes containing one intron; ** genes containing two introns; T trans-splicing of the related gene. Genes in bold are located in inverted repeats (two gene copies in the genome).
Table
Table 4. Identification of potential Meliaceae-specific SNPs (highlighted in grey) in the matK gene.
Table 4. Identification of potential Meliaceae-specific SNPs (highlighted in grey) in the matK gene.
Position (Consensus)Position (Cedrela odorata)MeliaceaeSapindales without MeliaceaeRosids without SapindalesAll without Rosids
208208TA or C or TC or AC or A
280262T or C 1
346328GC or TC or AC or A
402378AC or T or A or GA or CA or C
574550C or T 3
618588C or A 1
639609GT or G or A or CT or G or A or CG or A or T
861819TCC or T or GC
995953C or T 4
11941146TC or TC or TC
12371189G or A 2 or C 1
12391191G or A 3
13181270GAAA
13891341C or T 1
14781430CG or AG or A or TA or G
14941446TA or GG or AG or C
The SNP positions listed were validated in different multiple alignments of matK sequences from member individuals of different taxonomic groups downloaded from GenBank (100 sequences each). The “position (consensus)” refers to the position in the consensus sequence in the alignment in Figure S5. 1,2,3,4 Only 1/2/3/4 individual(s) show the indicated nucleotide.

Share and Cite

MDPI and ACS Style

Mader, M.; Pakull, B.; Blanc-Jolivet, C.; Paulini-Drewes, M.; Bouda, Z.H.-N.; Degen, B.; Small, I.; Kersten, B. Complete Chloroplast Genome Sequences of Four Meliaceae Species and Comparative Analyses. Int. J. Mol. Sci. 2018, 19, 701. https://doi.org/10.3390/ijms19030701

AMA Style

Mader M, Pakull B, Blanc-Jolivet C, Paulini-Drewes M, Bouda ZH-N, Degen B, Small I, Kersten B. Complete Chloroplast Genome Sequences of Four Meliaceae Species and Comparative Analyses. International Journal of Molecular Sciences. 2018; 19(3):701. https://doi.org/10.3390/ijms19030701

Chicago/Turabian Style

Mader, Malte, Birte Pakull, Céline Blanc-Jolivet, Maike Paulini-Drewes, Zoéwindé Henri-Noël Bouda, Bernd Degen, Ian Small, and Birgit Kersten. 2018. "Complete Chloroplast Genome Sequences of Four Meliaceae Species and Comparative Analyses" International Journal of Molecular Sciences 19, no. 3: 701. https://doi.org/10.3390/ijms19030701

APA Style

Mader, M., Pakull, B., Blanc-Jolivet, C., Paulini-Drewes, M., Bouda, Z. H. -N., Degen, B., Small, I., & Kersten, B. (2018). Complete Chloroplast Genome Sequences of Four Meliaceae Species and Comparative Analyses. International Journal of Molecular Sciences, 19(3), 701. https://doi.org/10.3390/ijms19030701

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