Comparative Analysis of Chloroplast Genomes of 19 Saxifraga Species, Mostly from the European Alps

Complete chloroplast genome sequences are widely used in the analyses of phylogenetic relationships among angiosperms. As a species-rich genus, species diversity centers of Saxifraga L. include mountainous regions of Eurasia, such as the Alps and the Qinghai-Tibetan Plateau (QTP) sensu lato. However, to date, datasets of chloroplast genomes of Saxifraga have been concentrated on the QTP species; those from European Alps are largely unavailable, which hinders comprehensively comparative and evolutionary analyses of chloroplast genomes in this genus. Here, complete chloroplast genomes of 19 Saxifraga species were de novo sequenced, assembled and annotated, and of these 15 species from Alps were reported for the first time. Subsequent comparative analysis and phylogenetic reconstruction were also conducted. Chloroplast genome length of the 19 Saxifraga species range from 149,217 bp to 152,282 bp with a typical quadripartite structure. All individual chloroplast genome included in this study contains 113 unique genes, including 79 protein-coding genes, four rRNAs and 30 tRNAs. The IR boundaries keep relatively conserved with minor expansion in S. consanguinea. mVISTA analysis and identification of polymorphic loci for molecular markers shows that six intergenic regions (ndhC-trnV, psbE-petL, rpl32-trnL, rps16-trnQ, trnF-ndhJ, trnS-trnG) can be selected as the potential DNA barcodes. A total of 1204 SSRs, 433 tandem repeats and 534 Large sequence repeats were identified in the 19 Saxifraga chloroplast genomes. The codon usage analysis revealed that Saxifraga chloroplast genome codon prefers to end in A/T. Phylogenetic reconstruction of 33 species (31 Saxifraga species included) based on 75 common protein coding genes received high bootstrap support values for nearly all identified nodes, and revealed a tree topology similar to previous studies.