Resolving robust phylogenetic relationships of core Brassicaceae using genome skimming data

The Brassicaceae is an economically and scientifically important family distributed globally, including oilseed rape and the model plant, Arabidopsis thaliana. Although growing molecular data have been used in phylogenetic studies, the relationships among major clades and tribes of Brassicaceae are still controversial. Here, we investigated the core Brassicaceae phylogenetics using 222 plastomes and 235 nrDNA cistrons, including 106 plastomes and 112 nrDNA cistrons assembled from newly sequenced genome skimming data of 112 taxa. The sampling covered 73 genera from 61.5% tribes and four unassigned genera and species. Three well supported lineages LI, LII, and LIII were revealed in our plastomic analyses, with LI sister to LII + LIII. In addition, the monophyly of the newly delimitated LII was strongly supported by three different partition strategies, concatenated methods under Bayesian and Maximum Likelihood analyses. LII comprised 13 tribes, including four tribes previously unassigned to any lineage, that is Biscutelleae as the earliest diverging clade and Cochlearieae as the sister to Megacarpaeeae + Anastaticeae. Within LII, the intertribal relationships were also well resolved, except that a conflicting position of Orychophragmus was detected among different datasets. In LIII, Shehbazia was resolved as a member of Chorisproreae, but Chorisproreae, Dontostemoneae, and Euclidieae were all resolved as paraphyletic, which was also confirmed by nrDNA analyses. Moreover, the loss of the rps16 gene was detected as likely to be a synapomorphy of the tribes Arabideae and Alysseae. Overall, using genome skimming data, we resolved robust phylogenetic relationships of core Brassicaceae and shed new light on the complex evolutionary history of this family.