Genotyping-by-sequencing as a tool to infer phylogeny and ancestral hybridization: a case study in Carex (Cyperaceae)

Determining phylogenetic relationships among very closely related species has remained a challenge for evolutionary biologists due to interlocus phylogenetic discordance and the difficulty of obtaining variable markers. Here, we used a Genotyping-by-Sequencing (GBS) approach to sample a reduced representation genomic data set and infer the phylogeny of seven closely related species in the genus Carex (Cyperaceae). Past attempts to reconstruct phylogenetic relationships among these species produced conflicting and poorly-supported results. We inferred a robust phylogeny based on >3000 GBS loci and >1300 SNPs (with a minimum sequence depth within individuals of 10) using maximum likelihood and Bayesian inference. We also tested for historical introgression using the D-statistic test. We compared these analyses with partitioned RAD analysis, which is designed to identify suboptimal trees reflecting secondary phylogenetic signal that may be obscured by the dominant signal in the data. Phylogenetic analyses yielded fully resolved trees with high support. We found two main clades, one grouping Carex scoparia populations and C. waponahkikensis, and a second clade grouping C. longii, C. vexans, C. suberecta and C. albolutescens. We detected marginally significant signals of introgression between C. scoparia and C. suberecta or C. albolutescens, and we rejected a hybrid origin hypothesis for C. waponahkikensis. Our results demonstrate the power of NGS data sets for resolving some of the most difficult phylogenetic challenges where traditional phylogenetic markers have failed.