The Lomatium foeniculaceum complex of the Perennial North American Clade (PENA) in Apiaceae is present across much of the Great Plains and Intermountain Region, and as far north as northeast British Columbia. It is comprised of five subspecies and is the most widely distributed complex in the PENA clade. As delineated by previous taxonomic treatments, these infrataxa, which were once recognized at the specific rank, are separated by the presence or absence of variable characters such as petal ciliation, petal color, and involucel geometry. Recent phylogenetic research has found that morphology alone in the PENA clade does not always reflect evolutionary history and is often misleading when used for species delimitation. Fifty populations of L. foeniculaceum were selected based on general distribution, and to include all named subspecies and areas where sympatry or peripatry of subspecies occurs. The Angiosperm353 bait kit is being employed in conjunction with next generation sequencing for this project. Phylogenetic analyses will be used to test congruence with associated historic climate data and soil sample data from representative sites. Results from this project will enhance the current understanding of species relations within the PENA clade and will help verify the taxonomic validity and subspecific boundaries within the complex.