Comparative analysis of chloroplast genomes of 19 Saxifraga species, mostly from Europe Alps

Complete chloroplast genome sequences are wildly 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, dataset of chloroplast genomes of Saxifraga has concentrated in the QTP species, those from European Alps are largely unavailable, which hinders comprehensive 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, of which 15 species from Alps are first reported. 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.