Highly Variable Chloroplast Markers for Evaluating Plant Phylogeny at Low Taxonomic Levels and for DNA Barcoding

BACKGROUND: At present, plant molecular systematics and DNA barcoding techniques rely heavily on the use of chloroplast gene sequences. Because of the relatively low evolutionary rates of chloroplast genes, there are very few choices suitable for molecular studies on angiosperms at low taxonomic levels, and for DNA barcoding of species. METHODOLOGY/PRINCIPAL FINDINGS: We scanned the entire chloroplast genomes of 12 genera to search for highly variable regions. The sequence data of 9 genera were from GenBank and 3 genera were of our own. We identified nearly 5% of the most variable loci from all variable loci in the chloroplast genomes of each genus, and then selected 23 loci that were present in at least three genera. The 23 loci included 4 coding regions, 2 introns, and 17 intergenic spacers. Of the 23 loci, the most variable (in order from highest variability to lowest) were intergenic regions ycf1-a, trnK, rpl32-trnL, and trnH-psbA, followed by trnS(UGA)-trnG(UCC), petA-psbJ, rps16-trnQ, ndhC-trnV, ycf1-b, ndhF, rpoB-trnC, psbE-petL, and rbcL-accD. Three loci, trnS(UGA)-trnG(UCC), trnT-psbD, and trnW-psaJ, showed very high nucleotide diversity per site (π values) across three genera. Other loci may have strong potential for resolving phylogenetic and species identification problems at the species level. The loci accD-psaI, rbcL-accD, rpl32-trnL, rps16-trnQ, and ycf1 are absent from some genera. To amplify and sequence the highly variable loci identified in this study, we designed primers from their conserved flanking regions. We tested the applicability of the primers to amplify target sequences in eight species representing basal angiosperms, monocots, eudicots, rosids, and asterids, and confirmed that the primers amplified the desired sequences of these species. SIGNIFICANCE/CONCLUSIONS: Chloroplast genome sequences contain regions that are highly variable. Such regions are the first consideration when screening the suitable loci to resolve closely related species or genera in phylogenetic analyses, and for DNA barcoding

Molecular evolution and phylogenomic analysis of complete chloroplast genomes of Cotinus (Anacardiaceae)

Abstract Cotinus is an oligo‐specific ornamentally valuable genus with a disjunct distribution in the Northern Hemisphere. Traditionally, the taxonomy of Cotinus was mainly based on leaf morphological characteristics. However, the limited availability of genomic information greatly hindered the study of molecular evolution and phylogeny of this genus. This study sequenced the chloroplast (cp) genomes of all currently recognized taxa of Cotinus, including three species and four varieties. A comparative analysis was performed to investigate their cp genome characteristics and evolution. Furthermore, we inferred the phylogenetic relationships of Cotinus based on whole cp genomes, protein‐coding genes, and nuclear ITS data. All cp genomes exhibited a typical quadripartite structure with genome sizes ranging from 158,865 to 160,155 bp. A total of 113–114 genes were identified in the genomes. Seven non‐coding and four coding regions were identified as the most divergent hotspots for potential molecular barcodes and phylogenetic markers. Selection pressure analysis showed that there had been positive selection on genes matK and rps8 in the Cotinus cp genomes. Phylogenetic results confirmed that Cotinus is a monophyletic group but the widely distributed species Cotinus coggygria is not monophyletic. The divergence‐time analysis suggested that Cotinus underwent an evolutionary divergence from the middle Eocene and rapid adaptive radiation from the middle Miocene. This study revealed new insights into the cp genome evolution and phylogeny of Cotinus and related taxa

Complete chloroplast genome of Sedum sarmentosum and chloroplast genome evolution in Saxifragales

Comparative chloroplast genome analyses are mostly carried out at lower taxonomic levels, such as the family and genus levels. At higher taxonomic levels, chloroplast genomes are generally used to reconstruct phylogenies. However, little attention has been paid to chloroplast genome evolution within orders. Here, we present the chloroplast genome of Sedum sarmentosum and take advantage of several available (or elucidated) chloroplast genomes to examine the evolution of chloroplast genomes in Saxifragales. The chloroplast genome of S. sarmentosum is 150,448 bp long and includes 82,212 bp of a large single-copy (LSC) region, 16.670 bp of a small single-copy (SSC) region, and a pair of 25,783 bp sequences of inverted repeats (IRs).The genome contains 131 unique genes, 18 of which are duplicated within the IRs. Based on a comparative analysis of chloroplast genomes from four representative Saxifragales families, we observed two gene losses and two pseudogenes in Paeonia obovata, and the loss of an intron was detected in the rps16 gene of Penthorum chinense. Comparisons among the 72 common protein-coding genes confirmed that the chloroplast genomes of S. sarmentosum and Paeonia obovata exhibit accelerated sequence evolution. Furthermore, a strong correlation was observed between the rates of genome evolution and genome size. The detected genome size variations are predominantly caused by the length of intergenic spacers, rather than losses of genes and introns, gene pseudogenization or IR expansion or contraction. The genome sizes of these species are negatively correlated with nucleotide substitution rates. Species with shorter duration of the life cycle tend to exhibit shorter chloroplast genomes than those with longer life cycles

Plastome evolution and phylogeny of the tribe Ruteae (Rutaceae)

Abstract Rutaceae is a large family, and the genus‐level classification in the subfamilies or tribes of this family is not unified based on different taxonomic treatments. Until now, phylogenetic relationships of some genera in traditional tribe Ruteae have not been clearly resolved. In this study, seven new complete plastomes of this tribe were sequenced, and a comparative analysis was performed to investigate their plastome characteristics and evolution. In addition, we inferred the phylogenetic relationships of Ruteae based on complete plastome and nuclear ITS data. All plastomes exhibited a typical quadripartite structure and were relatively conserved in their structure and gene arrangement. Their genome sizes ranged from 154,656 bp to 160,677 bp, and the size variation was found to be associated with differences in IR expansion and gene loss. A total of 112 to 114 genes were identified in the genomes, including 78 to 79 protein‐coding genes, 30 tRNA genes, 4 rRNA genes, and 2 pseudogenes. Sequence divergence analysis indicated that non‐coding regions exhibited a higher percentage of variable characters, and nine non‐coding and six coding regions were identified as divergent hotspots. Phylogenetic results based on different datasets showed that this tribe was divided into three reciprocally exclusive groups. The phylogenetic analyses between plastome and nuclear ITS data were partly incongruent with each other. This study provides new insights into plastome evolution of Ruteae as well as Rutaceae. The availability of these plastomes provides useful genomic resources for molecular DNA barcodes and phylogenetically informative markers and deepens our understanding of the phylogeny in Ruteae

Phylogenetic Resolution in Juglans Based on Complete Chloroplast Genomes and Nuclear DNA Sequences

Walnuts (Juglans of the Juglandaceae) are well-known economically important resource plants for the edible nuts, high-quality wood, and medicinal use, with a distribution from tropical to temperate zones and from Asia to Europe and Americas. There are about 21 species in Juglans. Classification of Juglans at section level is problematic, because the phylogenetic position of Juglans cinerea is disputable. Lacking morphological and DNA markers severely inhibited the development of related researches. In this study, the complete chloroplast genomes and two nuclear DNA regions (the internal transcribed spacer and ubiquitin ligase gene) of 10 representative taxa of Juglans were used for comparative genomic analyses in order to deepen the understanding on the application value of genetic information for inferring the phylogenetic relationship of the genus. The Juglans chloroplast genomes possessed the typical quadripartite structure of angiosperms, consisting of a pair of inverted repeat regions separated by a large single-copy region and a small single-copy region. All the 10 chloroplast genomes possessed 112 unique genes arranged in the same order, including 78 protein-coding, 30 tRNA, and 4 rRNA genes. A combined sequence data set from two nuclear DNA regions revealed that Juglans plants could be classified into three branches: (1) section Juglans, (2) section Cardiocaryon including J. cinerea which is closer to J. mandshurica, and (3) section Rhysocaryon. However, three branches with a different phylogenetic topology were recognized in Juglans using the complete chloroplast genome sequences: (1) section Juglans, (2) section Cardiocaryon, and (3) section Rhysocaryon plus J. cinerea. The molecular taxonomy of Juglans is almost compatible to the morphological taxonomy except J. cinerea (section Trachycaryon). Based on the complete chloroplast genome sequence data, the divergence time between section Juglans and section Cardiocaryon was 44.77 Mya, while section Rhysocaryon diverged from other sections in the genus Juglans was 47.61 Mya. Eleven of the 12 small inversions in the chloroplast genomes provided valuable phylogenetic information for classification of walnut plants at section and species levels. Our results are valuable for future studies on Juglans genetic diversity and will enhance the understanding on the phylogenetic evolution of Juglandaceae

Evolutionary directions of single nucleotide substitutions and structural mutations in the chloroplast genomes of the family Calycanthaceae

Background Chloroplast genome sequence data is very useful in studying/addressing the phylogeny of plants at various taxonomic ranks. However, there are no empirical observations on the patterns, directions, and mutation rates, which are the key topics in chloroplast genome evolution. In this study, we used Calycanthaceae as a model to investigate the evolutionary patterns, directions and rates of both nucleotide substitutions and structural mutations at different taxonomic ranks. Results There were 2861 polymorphic nucleotide sites on the five chloroplast genomes, and 98% of polymorphic sites were biallelic. There was a single-nucleotide substitution bias in chloroplast genomes. A -> T or T -> A (2.84%) and G -> C or C -> G (3.65%) were found to occur significantly less frequently than the other four transversion mutation types. Synonymous mutations kept balanced pace with nonsynonymous mutations, whereas biased directions appeared between transition and transversion mutations and among transversion mutations. Of the structural mutations, indels and repeats had obvious directions, but microsatellites and inversions were non-directional. Structural mutations increased the single nucleotide mutations rates. The mutation rates per site per year were estimated to be 0.14-0.34 x 10(- 9)for nucleotide substitution at different taxonomic ranks, 0.64 x 10(- 11)for indels and 1.0 x 10(- 11)for repeats. Conclusions Our direct counts of chloroplast genome evolution events provide raw data for correctly modeling the evolution of sequence data for phylogenetic inferences

Determination of a criminal suspect using environmental plant DNA metabarcoding technology

There are criminal cases that no frequently used evidence, for example, human DNAs from the criminal, is available. Such cases usually are unresolvable. With the advent of DNA metabarcoding, evidences are mined from environmental DNA and such cases become resolvable. This study reports how a criminal suspect was determined by environmental plant DNA metabarcoding technology. A girl was killed in a rural wet area in China without a witness or video record. Pants with dried mud was found from one of her classmate's house. The mud was removed from the pants and 11 more mud or soil samples surrounding murder scene were collected. DNA was extracted from the soil. Chloroplast rbcL gene were amplified and sequenced on a next generation sequencing platform. After bioinformatics analysis, ZOTU composition of 12 samples de-monstrated that the mud on the suspect's pants was from the criminal scene. The suspect finally made a clean breast of his crime. This case implies that plant DNA in the environment soil is a new source of evidence in determination of suspects using DNA metabarcoding technology and has high potentials of extensive applications in criminal cases. (C) 2021 Elsevier B.V. All rights reserved

Chloroplast Genomes Evolution and Phylogenetic Relationships of <i>Caragana</i> species

Caragana sensu lato (s.l.) includes approximately 100 species that are mainly distributed in arid and semi-arid regions. Caragana species are ecologically valuable for their roles in windbreaking and sand fixation. However, the taxonomy and phylogenetic relationships of the genus Caragana are still unclear. In this study, we sequenced and assembled the chloroplast genomes of representative species of Caragana and reconstructed robust phylogenetic relationships at the section level. The Caragana chloroplast genome has lost the inverted repeat region and wascategorized in the inverted repeat loss clade (IRLC). The chloroplast genomes of the eight species ranged from 128,458 bp to 135,401 bp and contained 110 unique genes. All the Caragana chloroplast genomes have a highly conserved structure and gene order. The number of long repeats and simple sequence repeats (SSRs) showed significant variation among the eight species, indicating heterogeneous evolution in Caragana. Selective pressure analysis of the genes revealed that most of the protein-coding genes evolved under purifying selection. The phylogenetic analyses indicated that each section forms a clade, except the section Spinosae, which was divided into two clades. This study elucidated the evolution of the chloroplast genome within the widely distributed genus Caragana. The detailed information obtained from this study can serve as a valuable resource for understanding the molecular dynamics and phylogenetic relationships within Caragana

Editorial : Rise to the challenges in plastome phylogenomics

Non peer reviewe