Resolving interspecific relationships within evolutionarily young lineages using RNA-seq data: An example from Pedicularis section Cyathophora (Orobanchaceae).

Phylogenomics has shown great potential in resolving evolutionary relationships at different taxonomical levels. However, it remains controversial whether all orthologous genes under different selective pressures can be concatenated for phylogenomic reconstruction. Here we used sect. Cyathophora of Pedicularis, one of the most species-rich genera of angiosperms in the alpine and arctic regions of the Northern Hemisphere, as a model to investigate the efficiency of RNA-seq in resolving relationships of closely related congeneric species. Flower transcriptomes were sequenced for all species of sect. Cyathophora and two outgroup species. Forty-one highly conserved single-/low-copy nuclear genes and 1553 orthologous groups (OGs) were identified and concatenated into five datasets based on gene copy numbers and Ka/Ks values to reconstruct the phylogeny of section Cyathophora. We also tested how many genes minimally can resolve the interspecific relationships, and reconstructed the evolutionary history of some floral characters in sect. Cyathophora. The results showed that the five different datasets consistently resolved the interspecific relationships of sect. Cyathophora, and the interspecific relationships can be robustly reconstructed with maximal support when ⩾20 single-/low-copy nuclear genes or 25 OGs are used. Our study suggests that all OGs under different selective pressures can be concatenated for phylogenomic reconstruction, and provides a successful and efficient use of RNA-seq in reconstructing interspecific relationships of a non-model plant group with recent radiations.