A lesson from polybutylene succinate plastisphere to the discovery of novel plastic degrading enzyme genes in marine vibrios

Polybutylene succinate (PBS) is an eco-friendly green plastic. However, PBS was shown as being non-biodegradable in marine environments, and up until now, only a limited number of PBS-degrading marine microbes have been discovered. We first set up in vitro PBS- and PBSA (polybutylene succinate adipate)-plastispheres to characterize novel PBS-degrading marine microbes. Microbial growth and oxygen consumption were observed in both PBS- and PBSA-plastispheres enriched with natural seawater collected from Usujiri, Hokkaido, Japan, and Vibrionaceae and Pseudoalteromonadaceae were significantly enriched on these films. Further gene identification indicated that vibrios belonging to the Gazogenes clade possess genes related to a PBS degrading enzyme (PBSase). The PBS degradation assay for six Gazogenes clade vibrios identified Vibrio ruber, Vibrio rhizosphaerae, and Vibrio spartinae as being capable of degrading PBS. We further identified the gene responsible for PBSase from the type strain of V. ruber, and the purified recombinant vibrio PBSase was found to have low-temperature adaptation and was active under high NaCl concentrations. We also provided docking models between the vibrio PBSase and PBS and PBSA units to show how vibrio PBSase interacts with each substrate compared to the Acidovorax PBSase. These results could contribute to a more sustainable society through further utilization of PBS in marine environments and plastic recycling.     

On the Origins of a Vibrio Species

Thirty-two genome sequences of various Vibrionaceae members are compared, with emphasis on what makes V. cholerae unique. As few as 1,000 gene families are conserved across all the Vibrionaceae genomes analysed; this fraction roughly doubles for gene families conserved within the species V. cholerae. Of these, approximately 200 gene families that cluster on various locations of the genome are not found in other sequenced Vibrionaceae; these are possibly unique to the V. cholerae species. By comparing gene family content of the analysed genomes, the relatedness to a particular species is identified for two unspeciated genomes. Conversely, two genomes presumably belonging to the same species have suspiciously dissimilar gene family content. We are able to identify a number of genes that are conserved in, and unique to, V. cholerae. Some of these genes may be crucial to the niche adaptation of this species.     

Genome-Wide SNP-Genotyping Array to Study the Evolution of the Human Pathogen Vibrio vulnificus Biotype 3

Vibrio vulnificus is an aquatic bacterium and an important human pathogen. Strains of V. vulnificus are classified into three different biotypes. The newly emerged biotype 3 has been found to be clonal and restricted to Israel. In the family Vibrionaceae, horizontal gene transfer is the main mechanism responsible for the emergence of new pathogen groups. To better understand the evolution of the bacterium, and in particular to trace the evolution of biotype 3, we performed genome-wide SNP genotyping of 254 clinical and environmental V. vulnificus isolates with worldwide distribution recovered over a 30-year period, representing all phylogeny groups. A custom single-nucleotide polymorphism (SNP) array implemented on the Illumina GoldenGate platform was developed based on 570 SNPs randomly distributed throughout the genome. In general, the genotyping results divided the V. vulnificus species into three main phylogenetic lineages and an additional subgroup, clade B, consisting of environmental and clinical isolates from Israel. Data analysis suggested that 69% of biotype 3 SNPs are similar to SNPs from clade B, indicating that biotype 3 and clade B have a common ancestor. The rest of the biotype 3 SNPs were scattered along the biotype 3 genome, probably representing multiple chromosomal segments that may have been horizontally inserted into the clade B recipient core genome from other phylogroups or bacterial species sharing the same ecological niche. Results emphasize the continuous evolution of V. vulnificus and support the emergence of new pathogenic groups within this species as a recurrent phenomenon. Our findings contribute to a broader understanding of the evolution of this human pathogen.     

Comprehensive analysis of the complete mitochondrial genome of Melanoplus differentialis (Acrididae: Melanoplinae) captured in Korea

In this study, we determined the complete mitochondrial genome of the invasive insect species Melanoplus differentialis captured in Korea. The complete mitochondrial genome of M. differentialis is 15,625 bp long and comprises 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNAs, with a GC ratio of 25.2%. In total, 353 SNPs and 11 INDEL regions (total length 67 bp) were found against the previously sequenced M. differentialis mitochondrial genome recorded as public genome data. The number of interspecific variations was greater than the number of intraspecific variations in this insect. Phylogenetic tree analysis showed that the mitochondrial genome clustered the Melanoplus clade with two previously reported Melanoplus sequences. However, the sequences were not divided at the species-level clade possibly as a consequence of misidentification caused by an error in the public database. Our results extend the molecular database status of Melanoplus by providing a novel complete mitochondrial genome sequence for M. differentialis that could serve as reference for further molecular studies.     

Comparative sequence analyses indicate that Coffea (Asterids) and Vitis (Rosids) derive from the same paleo-hexaploid ancestral genome

The complete sequence of Vitis vinifera revealed that the rosid clade derives from a hexaploid ancestor. At present, no analysis of complete genome sequence is available for an asterid, the other large eudicot clade, which includes the economically important species potato, tomato and coffee. To elucidate the genomic history of asterids, we compared the sequence of an 800 kb region of diploid Coffea genome to the orthologous regions of V. vinifera, Populus trichocarpa and Arabidopsis thaliana. We found a very high level of collinearity between around 80 genes of the three rosid species and Coffea. Collinearity comparisons between orthologous and paralogous regions indicates that (1) the Coffea (and consequently all asterids) and rosids share the same hexaploid ancestor; (2) the diploidization process (loss of duplicated and redundant copies from the whole genome duplication) was very advanced in the most recent common ancestor of rosids and asterids. Finally, no additional polyploidization events were detected in the Coffea lineage. Differences in gene loss rates were detected among the three rosid species and linked to the divergence in protein sequences.     

Competence and natural transformation in vibrios

Natural transformation is a major mechanism of horizontal gene transfer in bacteria. By incorporating exogenous DNA elements into chromosomes, bacteria are able to acquire new traits that can enhance their fitness in different environments. Within the past decade, numerous studies have revealed that natural transformation is prevalent among members of the Vibrionaceae, including the pathogen Vibrio cholerae. Four environmental factors: (i) nutrient limitation, (ii) availability of extracellular nucleosides, (iii) high cell density and (iv) the presence of chitin, promote genetic competence and natural transformation in Vibrio cholerae by co‐ordinating expression of the regulators CRP, CytR, HapR and TfoX respectively. Studies of other Vibrionaceae members highlight the general importance of natural transformation within this bacterial family.     

Beyond an AFLP genome scan towards the identification of immune genes involved in plague resistance in Rattus rattus from Madagascar

Genome scans using amplified fragment length polymorphism (AFLP) markers became popular in nonmodel species within the last 10 years, but few studies have tried to characterize the anonymous outliers identified. This study follows on from an AFLP genome scan in the black rat (Rattus rattus), the reservoir of plague (Yersinia pestis infection) in Madagascar. We successfully sequenced 17 of the 22 markers previously shown to be potentially affected by plague‐mediated selection and associated with a plague resistance phenotype. Searching these sequences in the genome of the closely related species Rattus norvegicus assigned them to 14 genomic regions, revealing a random distribution of outliers in the genome (no clustering). We compared these results with those of an in silico AFLP study of the R. norvegicus genome, which showed that outlier sequences could not have been inferred by this method in R. rattus (only four of the 15 sequences were predicted). However, in silico analysis allowed the prediction of AFLP markers distribution and the estimation of homoplasy rates, confirming its potential utility for designing AFLP studies in nonmodel species. The 14 genomic regions surrounding AFLP outliers (less than 300 kb from the marker) contained 75 genes encoding proteins of known function, including nine involved in immune function and pathogen defence. We identified the two interleukin 1 genes (Il1a and Il1b) that share homology with an antigen of Y. pestis, as the best candidates for genes subject to plague‐mediated natural selection. At least six other genes known to be involved in proinflammatory pathways may also be affected by plague‐mediated selection.     

Advanced Microbial Taxonomy Combined with Genome-Based-Approaches Reveals that Vibrio astriarenae sp. nov., an Agarolytic Marine Bacterium,Forms a New Clade in Vibrionaceae

Advances in genomic microbial taxonomy have opened the way to create a more universal and transparent concept of species but is still in a transitional stage towards becoming a defining robust criteria for describing new microbial species with minimum features obtained using both genome and classical polyphasic taxonomies. Here we performed advanced microbial taxonomies combined with both genome-based and classical approaches for new agarolytic vibrio isolates to describe not only a novel Vibrio species but also a member of a new Vibrio clade. Two novel vibrio strains (Vibrio astriarenae sp. nov. C7^T and C20) showing agarolytic, halophilic and fermentative metabolic activity were isolated from a seawater sample collected in a coral reef in Okinawa. Intraspecific similarities of the isolates were identical in both sequences on the 16S rRNA and pyrH genes, but the closest relatives on the molecular phylogenetic trees on the basis of 16S rRNA and pyrH gene sequences were V. hangzhouensis JCM 15146^T (97.8% similarity) and V. agarivorans CECT 5085^T (97.3% similarity), respectively. Further multilocus sequence analysis (MLSA) on the basis of 8 protein coding genes (ftsZ, gapA, gyrB, mreB, pyrH, recA, rpoA, and topA) obtained by the genome sequences clearly showed the V. astriarenae strain C7^T and C20 formed a distinct new clade protruded next to V. agarivorans CECT 5085^T. The singleton V. agarivorans has never been included in previous MLSA of Vibrionaceae due to the lack of some gene sequences. Now the gene sequences are completed and analysis of 100 taxa in total provided a clear picture describing the association of V. agarivorans into pre-existing concatenated network tree and concluded its relationship to our vibrio strains. Experimental DNA-DNA hybridization (DDH) data showed that the strains C7^T and C20 were conspecific but were separated from all of the other Vibrio species related on the basis of both 16S rRNA and pyrH gene phylogenies (e.g., V. agarivorans CECT 5085^T, V. hangzhouensis JCM 15146^T V. maritimus LMG 25439^T, and V. variabilis LMG 25438^T). In silico DDH data also supported the genomic relationship. The strains C7^T also had less than 95% average amino acid identity (AAI) and average nucleotide identity (ANI) towards V. maritimus C210, V. variabilis C206, and V. mediterranei AK1^T, V. brasiliensis LMG 20546^T, V. orientalis ATCC 33934^T, and V. sinaloensis DSM 21326. The name Vibrio astriarenae sp. nov. is proposed with C7 as the type strains. Both V. agarivorans CECT 5058^T and V. astriarenae C7^T are members of the newest clade of Vibrionaceae named Agarivorans.     

Multiple nuclear genes stabilize the phylogenetic backbone of the genus Quercus

Phylogenetic relationships among 108 oak species (genus Quercus L.) were inferred using DNA sequences of six nuclear genes selected from the existing genomic resources of the genus. Previous phylogenetic reconstructions based on traditional molecular markers are inconclusive at the deeper nodes. Overall, weak phylogenetic signals were obtained for each individual gene analysis, but stronger signals were obtained when gene sequences were concatenated. Our data support the recognition of six major intrageneric groups Cyclobalanopsis, Cerris, Ilex, Quercus, Lobatae and Protobalanus. Our analyses provide resolution at deeper nodes but with moderate support and a more robust infrageneric classification within the two major clades, the ‘Old World Oaks’ (Cyclobalanopsis, Cerris, Ilex) and ‘New World Oaks’ (Quercus, Lobatae, Protobalanus). However, depending on outgroup choice, our analysis yielded two alternative placements of the Cyclobalanopsis clade within the genus Quercus. When Castanea Mill. was chosen as outgroup, our data suggested that the genus Quercus comprised two clades corresponding to two subgenera as traditionally recognized by Camus: subgenus Euquercus Hickel and Camus and subgenus Cyclobalanopsis Øersted (Schneider). However, when Notholithocarpus Manos, Cannon and S. Oh was chosen as an outgroup subgenus Cyclobalanopsis clustered with Cerris and Ilex groups to form the Old World clade. To assess the placement of the root, we complemented our dataset with published data of ITS and CRC sequences. Based on the concatenated eight gene sequences, the most likely root position is at the split between the ‘Old World Oaks’ and the ‘New World Oaks’, which is one of the alternative positions suggested by our six gene analysis. Using a dating approach, we inferred an Eocene age for the primary divergences in Quercus and a root age of about 50–55 Ma, which agrees with palaeobotanical evidence. Finally, irrespective of the outgroup choice, our data boost the topology within the New World clade, where (Protobalanus + Quercus) is a sister clade of Lobatae. Inferred divergence ages within this clade and the Cerris–Ilex clade are generally younger than could be expected from the fossil record, indicating that morphological differentiation pre-dates genetic isolation in this clade.     

Sequencing of Australian wild rice genomes reveals ancestral relationships with domesticated rice

The related A genome species of the Oryza genus are the effective gene pool for rice. Here, we report draft genomes for two Australian wild A genome taxa: O. rufipogon‐like population, referred to as Taxon A, and O. meridionalis‐like population, referred to as Taxon B. These two taxa were sequenced and assembled by integration of short‐ and long‐read next‐generation sequencing (NGS) data to create a genomic platform for a wider rice gene pool. Here, we report that, despite the distinct chloroplast genome, the nuclear genome of the Australian Taxon A has a sequence that is much closer to that of domesticated rice (O. sativa) than to the other Australian wild populations. Analysis of 4643 genes in the A genome clade showed that the Australian annual, O. meridionalis, and related perennial taxa have the most divergent (around 3 million years) genome sequences relative to domesticated rice. A test for admixture showed possible introgression into the Australian Taxon A (diverged around 1.6 million years ago) especially from the wild indica/O. nivara clade in Asia. These results demonstrate that northern Australia may be the centre of diversity of the A genome Oryza and suggest the possibility that this might also be the centre of origin of this group and represent an important resource for rice improvement.     

The complete chloroplast genome sequence of the endangered species Syringa pinnatifolia (Oleaceae)

Syringa pinnatifolia is an endangered endemic species in China with important ornamental and medicinal value, and it needs urgent protection. Here, we report the complete chloroplast (cp) genome structure of S. pinnatifolia and its evolution is inferred through comparative studies with related species. The S. pinnatifolia cp genome was 155 326 bp and contained a large single copy region (LSC) of 86 167 bp and a small single copy region (SSC) of 17 775 bp, as well as a pair of inverted repeat regions (IRs) of 25 692 bp. A total of 113 unique genes were annotated, including 79 protein‐coding genes, 30 tRNA genes and four rRNA genes. The GC content of the S. pinnatifolia cp genome was 37.9%, and the corresponding values in the LSC, SSC and IR regions were 36.0, 32.1, 43.2% respectively. Repetitive sequences analysis revealed that the S. pinnatifolia cp genome contained 38 repeats. Microsatellite marker detection analysis identified 253 simple sequence repeats (SSRs), which provides opportunities for future studies of the population genetics and phylogenetic relationships of Syringa. Phylogenetic analysis of 29 selected cp genomes revealed that S. pinnatifolia is closely related to Syringa vulgaris and all 27 Lamiales species formed a clade separate from the two outgroup species. This newly characterized S. pinnatifolia chloroplast genome will provide a useful genomic resource of phylogenetic inference and the development of more genetic markers for species discrimination and population studies in the genus Syringa.     

Evolution and microsynteny of the apyrase gene family in three legume genomes

Apyrases have been suggested to play important roles in plant nutrition, photomorphogenesis, and nodulation. To help trace the evolution of these genes in the legumes—and possibly, the acquisition of new functions for nodulation—apyrase-containing BACs were sequenced from three legume genomes. Genomic sequences from Medicago truncatula, Glycine max and Lotus japonicus were compared to one another and to corresponding regions in Arabidopsis thaliana. A phylogenetic analysis of apyrase homologs from these regions and sequences from other legume species, as well as other plant families, identified a potentially legume-specific clade that contains a well-characterized soybean ( G. soja) apyrase, Gs52, as well as homologs from Dolichos, Lotus , Medicago and Pisum. Sister clades contain homologs from members of Brassicaceae, Solanaceae, Poaceae and Fabaceae. Comparisons of rates of change at synonymous and nonsynonymous sites in the Gs52 and sister clades show rapid evolution in the potentially legume-specific Gs52 clade. The genomic organization of the apyrase-containing BACs shows evidence of gene duplication, genomic rearrangement, and gene conversion among Gs52 homologs. Taken together, these results suggest a scenario of local apyrase gene duplication in an ancestor of the legumes, followed by functional diversification and increased rates of change in the new genes, and further duplications in the Galegae (which include the genera Medicago and Pisum). The study also provides a detailed comparison of genomic regions between two model genomes which are now being sequenced ( M. truncatula and L. japonicus), and a genome from an economically important legume species ( G. max).Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by A. Kondorosi     

Genomic Adaptation of the Lactobacillus casei Group

Lactobacillus casei, L. paracasei, and L. rhamnosus form a closely related taxonomic group (Lactobacillus casei group) within the facultatively heterofermentative lactobacilli. Here, we report the complete genome sequences of L. paracasei JCM 8130 and L. casei ATCC 393, and the draft genome sequence of L. paracasei COM0101, all of which were isolated from daily products. Furthermore, we re-annotated the genome of L. rhamnosus ATCC 53103 (also known as L. rhamnosus GG), which we have previously reported. We confirmed that ATCC 393 is distinct from other strains previously described as L. paracasei. The core genome of 10 completely sequenced strains of the L. casei group comprised 1,682 protein-coding genes. Although extensive genome-wide synteny was found among the L. casei group, the genomes of ATCC 53103, JCM 8130, and ATCC 393 contained genomic islands compared with L. paracasei ATCC 334. Several genomic islands, including carbohydrate utilization gene clusters, were found at the same loci in the chromosomes of the L. casei group. The spaCBA pilus gene cluster, which was first identified in GG, was also found in other strains of the L. casei group, but several L. paracasei strains including COM0101 contained truncated spaC gene. ATCC 53103 encoded a higher number of proteins involved in carbohydrate utilization compared with intestinal lactobacilli, and extracellular adhesion proteins, several of which are absent in other strains of the L. casei group. In addition to previously fully sequenced L. rhamnosus and L. paracasei strains, the complete genome sequences of L. casei will provide valuable insights into the evolution of the L. casei group.     

Continuation and replacement of Vibrio cholerae non-O1 clonal genomic groups isolated from Plecoglossus altivelis fish in freshwaters

The dissemination and abundances of Vibrio species in aquatic environments are of interest, as some species cause emerging diseases in humans and in aquatic organisms like fish. It is suggested that Vibrio cholerae non-O1 infections of Plecoglossus altivelis (‘ayu’) were spread to various parts of Japan through the annual transplantation of juvenile fish. To investigate this, we used genome-aided tracing of 17 V. cholerae strains isolated from ayu between the 1970s and 1990s in different Japanese freshwater systems. The strains formed a genomic clade distinct from all known clades, which we designate as the Ayu clade. Two clonal genomic groups identified within the clade, Ayu-1 and Ayu-2, persisted for a few years (between 1977 to 1979 and 1987 to 1990, respectively), and clonal replacement of Ayu-1 by Ayu-2 took place over an 8-year period. Despite the high similarity between Ayu-1 and Ayu-2 (> 99.9% identity and > 97% fraction of genomes shared), differences in their gene repertoires were found, raising the possibility that they are phenotypically distinct. These results highlight the importance of genome-based studies for understanding the long-term dynamics of populations over the timescale of years.     

Genomic organization and gene expression of the multiple globins in Atlantic cod: conservation of globin-flanking genes in chordates infers the origin of the vertebrate globin clusters

Background  

The vertebrate globin genes encoding the α- and β-subunits of the tetrameric hemoglobins are clustered at two unlinked loci. The highly conserved linear order of the genes flanking the hemoglobins provides a strong anchor for inferring common ancestry of the globin clusters. In fish, the number of α-β-linked globin genes varies considerably between different sublineages and seems to be related to prevailing physico-chemical conditions. Draft sequences of the Atlantic cod genome enabled us to determine the genomic organization of the globin repertoire in this marine species that copes with fluctuating environments of the temperate and Arctic regions.     

An unusual overrepresentation of genetic factors related to iron homeostasis in the genome of the fluorescent Pseudomonas sp. ABC1

Members of the genus Pseudomonas inhabit diverse environments, such as soil, water, plants and humans. The variability of habitats is reflected in the diversity of the structure and composition of their genomes. This cosmopolitan bacterial genus includes species of biotechnological, medical and environmental importance. In this study, we report on the most relevant genomic characteristics of Pseudomonas sp. strain ABC1, a siderophore-producing fluorescent strain recently isolated from soil. Phylogenomic analyses revealed that this strain corresponds to a novel species forming a sister clade of the recently proposed Pseudomonas kirkiae. The genomic information reveals an overrepresented repertoire of mechanisms to hoard iron when compared to related strains, including a high representation of fecI-fecR family genes related to iron regulation and acquisition. The genome of the Pseudomonas sp. ABC1 contains the genes for non-ribosomal peptide synthetases (NRPSs) of a novel putative Azotobacter-related pyoverdine-type siderophore, a yersiniabactin-type siderophore and an antimicrobial betalactone; the last two are found only in a limited number of Pseudomonas genomes. Strain ABC1 can produce siderophores in a low-cost medium, and the supernatants from cultures of this strain promote plant growth, highlighting their biotechnological potential as a sustainable industrial microorganism.     

Genomic and Phenotypic Diversity of Coastal Vibrio cholerae Strains Is Linked to Environmental Factors

Studies of Vibrio cholerae diversity have focused primarily on pathogenic isolates of the O1 and O139 serotypes. However, autochthonous environmental isolates of this species routinely display more extensive genetic diversity than the primarily clonal pathogenic strains. In this study, genomic and metabolic profiles of 41 non-O1/O139 environmental isolates from central California coastal waters and four clinical strains are used to characterize the core genome and metabolome of V. cholerae. Comparative genome hybridization using microarrays constructed from the fully sequenced V. cholerae O1 El Tor N16961 genome identified 2,787 core genes that approximated the projected species core genome within 1.6%. Core genes are almost universally present in strains with widely different niches, suggesting that these genes are essential for persistence in diverse aquatic environments. In contrast, the dispensable genes and phenotypic traits identified in this study should provide increased fitness for certain niche environments. Environmental parameters, measured in situ during sample collection, are correlated to the presence of specific dispensable genes and metabolic capabilities, including utilization of mannose, sialic acid, citrate, and chitosan oligosaccharides. These results identify gene content and metabolic pathways that are likely selected for in certain coastal environments and may influence V. cholerae population structure in aquatic environments.     

Detection of quorum sensing signal molecules in the family Vibrionaceae

Aims: The aim of this study was to detect the production of three kinds of quorum sensing (QS) signal molecules, i.e. the N‐acyl‐homoserine lactone (AHL), the autoinducer‐2 (AI‐2) and the cholerae autoinducer‐1‐like (CAI‐1‐like) molecules in 25 Vibrionaceae strains. Methods and Results: The QS signal molecules in 25 Vibrionaceae strains were detected with different biosensors. Except Salinivibrio costicola VIB288 and Vibrio natriegens VIB299, all the other 23 Vibrionaceae strains could produce one or more kinds of detectable QS signal molecules. Twenty‐one of the 25 strains were found to produce AHL signal molecules by using Vibrio harveyi JMH612 and Agrobacterium tumefaciens KYC55 (pJZ372; pJZ384; pJZ410) as biosensors. The AHL fingerprints of eight strains were detected by thin‐layer chromatography with Ag. tumefaciens KYC55, and two of them, i.e. V. mediterranei VIB296 and Aliivibrio logei VIB414 had a high diversity of AHLs. Twenty of the 25 strains were found to have the AI‐2 activity, and the luxS gene sequences in 18 strains were proved to be conserved by PCR amplification and sequencing. Only six (five Vibrio strains and A. logei VIB414) of the 25 strains possessed the CAI‐1‐like activity. A. logei VIB414, V. campbellii VIB285, V. furnissii VIB293, V. pomeroyi LMG20537 and two V. harveyi strains VIB571 and VIB645 were found to produce all the three kinds of QS signal molecules. Conclusions: The results indicated that the QS signal molecules, especially AHL and AI‐2 molecules, were widespread in the family Vibrionaceae. Significance and Impact of the Study: In response to a variety of environmental conditions and selection forces, the family Vibrionaceae produced QS signal molecules with great diversity and complexity. The knowledge we obtained from this study will be useful for further research on the roles of different QS signal molecules in this family.     

Single nucleotide polymorphism-based species phylogeny of greater fritillary butterflies (Lepidoptera: Nymphalidae: Speyeria) demonstrates widespread mitonuclear discordance

The systematics of Speyeria butterflies has historically been complicated by intraspecific variability that has challenged efforts to delimit species and reconstruct phylogenies. Our study presents a phylogenetic comparison of genomic single nucleotide polymorphisms (SNPs) and mitochondrial COI gene sequences, with comprehensive taxon sampling that includes 15 species and 46 subspecies. Increased sampling of genetic markers and taxa improved the match between genetic clusters, obtained with both phylogenetic and cluster-based analyses, and species previously detected using morphology, as well as showing two species delimitations that may need revision. We also recovered extensive mitonuclear discordance between genomic SNPs and the COI gene, confirming that mitochondrial DNA does not reliably identify several species at broad geographic scales. Resolution of the relationships of Speyeria species demonstrates the importance of sampling variation across the whole genome, and provides an essential foundation for understanding the evolution of this charismatic clade of North American butterflies.