Whole Genomes Reveal Evolutionary Relationships and Mechanisms Underlying Gene-Tree Discordance in Neodiprion Sawflies.

Rapidly evolving taxa are excellent models for understanding the mechanisms that give rise to biodiversity. However, developing an accurate historical framework for comparative analysis of such lineages remains a challenge due to ubiquitous incomplete lineage sorting and introgression. Here, we use a whole-genome alignment, multiple locus-sampling strategies, and summary-tree and SNP-based species-tree methods to infer a species tree for eastern North American Neodiprion species, a clade of pine-feeding sawflies (Order: Hymenopteran; Family: Diprionidae). We recovered a well-supported species tree that-except for three uncertain relationships-was robust to different strategies for analyzing whole-genome data. Nevertheless, underlying gene-tree discordance was high. To understand this genealogical variation, we used multiple linear regression to model site concordance factors estimated in 50-kb windows as a function of several genomic predictor variables. We found that site concordance factors tended to be higher in regions of the genome with more parsimony-informative sites, fewer singletons, less missing data, lower GC content, more genes, lower recombination rates, and lower D-statistics (less introgression). Together, these results suggest that incomplete lineage sorting, introgression, and genotyping error all shape the genomic landscape of gene-tree discordance in Neodiprion. More generally, our findings demonstrate how combining phylogenomic analysis with knowledge of local genomic features can reveal mechanisms that produce topological heterogeneity across genomes.

Phylogenomics and biogeography of sawflies and woodwasps (hymenoptera, symphyta).

Phylogenomic approaches have recently helped elucidate various insect relationships, but large-scale comprehensive analyses on relationships within sawflies and woodwasps are still lacking. Here, we infer the relationships and long-term biogeographic history of these hymenopteran groups using a large dataset of 354 UCE loci collected from 385 species that represent all major lineages. Early Hymenoptera started diversifying during the Early Triassic ∼249 Ma and spread all over the ancient supercontinent Pangaea. We recovered Xyeloidea as a monophyletic sister group to other Hymenoptera and Pamphilioidea as sister to Unicalcarida. Within the diverse family Tenthredinidae, our taxonomically and geographically expanded taxon sampling highlights the non-monophyly of several traditionally defined subfamilies. In addition, the recent removal of Athalia and related genera from the Tenthredinidae into the separate family Athaliidae is supported. The deep historical biogeography of the group is characterised by independent dispersals and re-colonisations between the northern (Laurasia) and southern (Gondwana) palaeocontinents. The breakup of these landmasses led to ancient vicariance in several Gondwanan lineages, while interchange across the Northern Hemisphere has continued until the Recent. The little-studied African sawfly fauna is likewise a diverse mixture of groups with varying routes of colonization. Our results reveal interesting parallels in the evolution and biogeography of early hymenopterans and other ancient insect groups.

The Ecology of Hybrid Incompatibilities.

Ecologically mediated selection against hybrids, caused by hybrid phenotypes fitting poorly into available niches, is typically viewed as distinct from selection caused by epistatic Dobzhansky-Muller hybrid incompatibilities. Here, we show how selection against transgressive phenotypes in hybrids manifests as incompatibility. After outlining our logic, we summarize current approaches for studying ecology-based selection on hybrids. We then quantitatively review QTL-mapping studies and find traits differing between parent taxa are typically polygenic. Next, we describe how verbal models of selection on hybrids translate to phenotypic and genetic fitness landscapes, highlighting emerging approaches for detecting polygenic incompatibilities. Finally, in a synthesis of published data, we report that trait transgression-and thus possibly extrinsic hybrid incompatibility in hybrids-escalates with the phenotypic divergence between parents. We discuss conceptual implications and conclude that studying the ecological basis of hybrid incompatibility will facilitate new discoveries about mechanisms of speciation.

Evolution of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae).

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

Multiple mechanisms contribute to isolation by environment in the redheaded pine sawfly, Neodiprion lecontei.

Isolation by environment (IBE) is a population genomic pattern that arises when ecological barriers reduce gene flow between populations. Although current evidence suggests IBE is common in nature, few studies have evaluated the underlying mechanisms that generate IBE patterns. In this study, we evaluate five proposed mechanisms of IBE (natural selection against immigrants, sexual selection against immigrants, selection against hybrids, biased dispersal, and environment-based phenological differences) that may give rise to host-associated differentiation within a sympatric population of the redheaded pine sawfly, Neodiprion lecontei, a species for which IBE has previously been detected. We first characterize the three pine species used by N. lecontei at the site, finding morphological and chemical differences among the hosts that could generate divergent selection on sawfly host-use traits. Next, using morphometrics and ddRAD sequencing, we detect modest phenotypic and genetic differentiation among sawflies originating from different pines that is consistent with recent, in situ divergence. Finally, via a series of laboratory assays-including assessments of larval performance on different hosts, adult mate and host preferences, hybrid fitness, and adult eclosion timing-we find evidence that multiple mechanisms contribute to IBE in N. lecontei. Overall, our results suggest IBE can emerge quickly, possibly due to multiple mechanisms acting in concert to reduce migration between different environments.

A Test of Haldane's Rule in Neodiprion Sawflies and Implications for the Evolution of Postzygotic Isolation in Haplodiploids.

AbstractHaldane's rule-a pattern in which hybrid sterility or inviability is observed in the heterogametic sex of an interspecific cross-is one of the most widely obeyed rules in nature. Because inheritance patterns are similar for sex chromosomes and haplodiploid genomes, Haldane's rule may apply to haplodiploid taxa, predicting that haploid male hybrids will evolve sterility or inviability before diploid female hybrids. However, there are several genetic and evolutionary mechanisms that may reduce the tendency of haplodiploids to obey Haldane's rule. Currently, there are insufficient data from haplodiploids to determine how frequently they adhere to Haldane's rule. To help fill this gap, we crossed a pair of haplodiploid hymenopteran species (Neodiprion lecontei and Neodiprion pinetum) and evaluated the viability and fertility of female and male hybrids. Despite considerable divergence, we found no evidence of reduced fertility in hybrids of either sex, consistent with the hypothesis that hybrid sterility evolves slowly in haplodiploids. For viability, we found a pattern opposite to that of Haldane's rule: hybrid females, but not males, had reduced viability. This reduction was most pronounced in one direction of the cross, possibly due to a cytoplasmic-nuclear incompatibility. We also found evidence of extrinsic postzygotic isolation in hybrids of both sexes, raising the possibility that this form or reproductive isolation tends to emerge early in speciation in host-specialized insects. Our work emphasizes the need for more studies on reproductive isolation in haplodiploids, which are abundant in nature but underrepresented in the speciation literature.

Whole-genome resequencing data support a single introduction of the invasive white pine sawfly, Diprion similis.

Biological introductions are unintended "natural experiments" that provide unique insights into evolutionary processes. Invasive phytophagous insects are of particular interest to evolutionary biologists studying adaptation, as introductions often require rapid adaptation to novel host plants. However, adaptive potential of invasive populations may be limited by reduced genetic diversity-a problem known as the "genetic paradox of invasions." One potential solution to this paradox is if there are multiple invasive waves that bolster genetic variation in invasive populations. Evaluating this hypothesis requires characterizing genetic variation and population structure in the invaded range. To this end, we assemble a reference genome and describe patterns of genetic variation in the introduced white pine sawfly, Diprion similis. This species was introduced to North America in 1914, where it has rapidly colonized the thin-needled eastern white pine (Pinus strobus), making it an ideal invasion system for studying adaptation to novel environments. To evaluate evidence of multiple introductions, we generated whole-genome resequencing data for 64 D. similis females sampled across the North American range. Both model-based and model-free clustering analyses supported a single population for North American D. similis. Within this population, we found evidence of isolation-by-distance and a pattern of declining heterozygosity with distance from the hypothesized introduction site. Together, these results support a single-introduction event. We consider implications of these findings for the genetic paradox of invasion and discuss priorities for future research in D. similis, a promising model system for invasion biology.

Body size as a magic trait in two plant-feeding insect species.

When gene flow accompanies speciation, recombination can decouple divergently selected loci and loci conferring reproductive isolation. This barrier to sympatric divergence disappears when assortative mating and disruptive selection involve the same "magic" trait. Although magic traits could be widespread, the relative importance of different types of magic traits to speciation remains unclear. Because body size frequently contributes to host adaptation and assortative mating in plant-feeding insects, we evaluated several magic trait predictions for this trait in a pair of sympatric Neodiprion sawfly species adapted to different pine hosts. A large morphological dataset revealed that sawfly adults from populations and species that use thicker-needled pines are consistently larger than those that use thinner-needled pines. Fitness data from recombinant backcross females revealed that egg size is under divergent selection between the preferred pines. Lastly, mating assays revealed strong size-assortative mating within and between species in three different crosses, with the strongest prezygotic isolation between populations that have the greatest interspecific size differences. Together, our data support body size as a magic trait in pine sawflies and possibly many other plant-feeding insects. Our work also demonstrates how intraspecific variation in morphology and ecology can cause geographic variation in the strength of prezygotic isolation.

Faster-haplodiploid evolution under divergence-with-gene-flow: Simulations and empirical data from pine-feeding hymenopterans.

Although haplodiploidy is widespread in nature, the evolutionary consequences of this mode of reproduction are not well characterized. Here, we examine how genome-wide hemizygosity and a lack of recombination in haploid males affects genomic differentiation in populations that diverge via natural selection while experiencing gene flow. First, we simulated diploid and haplodiploid "genomes" (500-kb loci) evolving under an isolation-with-migration model with mutation, drift, selection, migration and recombination; and examined differentiation at neutral sites both tightly and loosely linked to a divergently selected site. As long as there is divergent selection and migration, sex-limited hemizygosity and recombination cause elevated differentiation (i.e., produce a "faster-haplodiploid effect") in haplodiploid populations relative to otherwise equivalent diploid populations, for both recessive and codominant mutations. Second, we used genome-wide single nucleotide polymorphism data to model divergence history and describe patterns of genomic differentiation between sympatric populations of Neodiprion lecontei and N. pinetum, a pair of pine sawfly species (order: Hymenoptera; family: Diprionidae) that are specialized on different pine hosts. These analyses support a history of continuous gene exchange throughout divergence and reveal a pattern of heterogeneous genomic differentiation that is consistent with divergent selection on many unlinked loci. Third, using simulations of haplodiploid and diploid populations evolving according to the estimated divergence history of N. lecontei and N. pinetum, we found that divergent selection would lead to higher differentiation in haplodiploids. Based on these results, we hypothesize that haplodiploids undergo divergence-with-gene-flow and sympatric speciation more readily than diploids.

The impact of life stage and pigment source on the evolution of novel warning signal traits.

Our understanding of how novel warning color traits evolve in natural populations is largely based on studies of reproductive stages and organisms with endogenously produced pigmentation. In these systems, genetic drift is often required for novel alleles to overcome strong purifying selection stemming from frequency-dependent predation and positive assortative mating. Here, we integrate data from field surveys, predation experiments, population genomics, and phenotypic correlations to explain the origin and maintenance of geographic variation in a diet-based larval pigmentation trait in the redheaded pine sawfly (Neodiprion lecontei), a pine-feeding hymenopteran. Although our experiments confirm that N. lecontei larvae are indeed aposematic-and therefore likely to experience frequency-dependent predation-our genomic data do not support a historical demographic scenario that would have facilitated the spread of an initially deleterious allele via drift. Additionally, significantly elevated differentiation at a known color locus suggests that geographic variation in larval color is currently maintained by selection. Together, these data suggest that the novel white morph likely spread via selection. However, white body color does not enhance aposematic displays, nor is it correlated with enhanced chemical defense or immune function. Instead, the derived white-bodied morph is disproportionately abundant on a pine species with a reduced carotenoid content relative to other pine hosts, suggesting that bottom-up selection via host plants may have driven divergence among populations. Overall, our results suggest that life stage and pigment source can have a substantial impact on the evolution of novel warning signals, highlighting the need to investigate diverse aposematic taxa to develop a comprehensive understanding of color variation in nature.

Support for the adaptive decoupling hypothesis from whole-transcriptome profiles of a hypermetamorphic and sexually dimorphic insect, Neodiprion lecontei.

Though seemingly bizarre, the dramatic morphological and ecological transformation that occurs when immature life stages metamorphose into reproductive adults is one of the most successful developmental strategies on the planet. The adaptive decoupling hypothesis (ADH) proposes that metamorphosis is an adaptation for breaking developmental links between traits expressed in different life stages, thereby facilitating their independent evolution when exposed to opposing selection pressures. Here, we draw inspiration from the ADH to develop a conceptual framework for understanding changes in gene expression across ontogeny. We hypothesized that patterns of stage-biased and sex-biased gene expression are the product of both decoupling mechanisms and selection history. To test this hypothesis, we characterized transcriptome-wide patterns of gene-expression traits for three ecologically distinct larval stages (all male) and adult males and females of a hypermetamorphic insect (Neodiprion lecontei). We found that stage-biased gene expression was most pronounced between larval and adult males, which is consistent with the ADH. However, even in the absence of a metamorphic transition, considerable stage-biased expression was observed among morphologically and behaviourally distinct larval stages. Stage-biased expression was also observed across ecologically relevant Gene Ontology categories and genes, highlighting the role of ecology in shaping patterns of gene expression. We also found that the magnitude and prevalence of stage-biased expression far exceeded adult sex-biased expression. Overall, our results highlight how the ADH can shed light on transcriptome-wide patterns of gene expression in organisms with complex life cycles. For maximal insight, detailed knowledge of organismal ecology is also essential.

The population genetics of crypsis in vertebrates: recent insights from mice, hares, and lizards.

By combining well-established population genetic theory with high-throughput sequencing data from natural populations, major strides have recently been made in understanding how, why, and when vertebrate populations evolve crypsis. Here, we focus on background matching, a particular facet of crypsis that involves the ability of an organism to conceal itself through matching its color to the surrounding environment. While interesting in and of itself, the study of this phenotype has also provided fruitful population genetic insights into the interplay of strong positive selection with other evolutionary processes. Specifically, and predicated upon the findings of previous candidate gene association studies, a primary focus of this recent literature involves the realization that the inference of selection from DNA sequence data first requires a robust model of population demography in order to identify genomic regions which do not conform to neutral expectations. Moreover, these demographic estimates provide crucial information about the origin and timing of the onset of selective pressures associated with, for example, the colonization of a novel environment. Furthermore, such inference has revealed crypsis to be a particularly useful phenotype for investigating the interplay of migration and selection-with examples of gene flow constraining rates of adaptation, or alternatively providing the genetic variants that may ultimately sweep through the population. Here, we evaluate the underlying evidence, review the strengths and weaknesses of the many population genetic methodologies used in these studies, and discuss how these insights have aided our general understanding of the evolutionary process.

Predicting evolutionary predictability.

The observation that phenotypic convergence and genetic convergence are widespread in nature implies that evolution is at least somewhat predictable. But to what extent and under what circumstances? In other words, how predictable is evolutionary predictability? Answering this question requires going beyond documenting examples of repeated evolution to actually quantifying predictability at different hierarchical levels. At present, few such studies exist. In this issue of Molecular Ecology, Chaturvedi et al. () quantify the predictability of genomewide changes that accompany shifts to an introduced host plant (alfalfa) in populations of the Melissa blue butterfly (Lycaeides melissa). They evaluate predictability in two contexts: (i) overlap in host-associated loci among populations that have independently colonized alfalfa, and (ii) overlap between host-associated loci in nature and loci associated with host performance in laboratory experiments. Overall, they find that the genomic changes that accompany host shifts in this system are indeed somewhat predictable. However, the degree of predictability depends on the type of comparison (among natural populations vs. between natural and experimental populations), type of convergence (specific genomic locations vs. direction of allele frequency change), geographic scale (rangewide vs. specific population pairs) and location in the genome (autosomes vs. sex chromosomes). Together with a handful of comparable data sets, Chaturvedi et al.'s () work suggests that the relative contribution of stochastic and deterministic processes to genomewide responses to novel selection pressures may be highly variable, but possibly predictably so.

Genetic Basis of Body Color and Spotting Pattern in Redheaded Pine Sawfly Larvae (Neodiprion lecontei).

Pigmentation has emerged as a premier model for understanding the genetic basis of phenotypic evolution, and a growing catalog of color loci is starting to reveal biases in the mutations, genes, and genetic architectures underlying color variation in the wild. However, existing studies have sampled a limited subset of taxa, color traits, and developmental stages. To expand the existing sample of color loci, we performed QTL mapping analyses on two types of larval pigmentation traits that vary among populations of the redheaded pine sawfly (Neodiprion lecontei): carotenoid-based yellow body color and melanin-based spotting pattern. For both traits, our QTL models explained a substantial proportion of phenotypic variation and suggested a genetic architecture that is neither monogenic nor highly polygenic. Additionally, we used our linkage map to anchor the current N. lecontei genome assembly. With these data, we identified promising candidate genes underlying (1) a loss of yellow pigmentation in populations in the mid-Atlantic/northeastern United States [C locus-associated membrane protein homologous to a mammalian HDL receptor-2 gene (Cameo2) and lipid transfer particle apolipoproteins II and I gene (apoLTP-II/I)], and (2) a pronounced reduction in black spotting in Great Lakes populations [members of the yellow gene family, tyrosine hydroxylase gene (pale), and dopamine N-acetyltransferase gene (Dat)]. Several of these genes also contribute to color variation in other wild and domesticated taxa. Overall, our findings are consistent with the hypothesis that predictable genes of large effect contribute to color evolution in nature.

The Evolutionary History of Nebraska Deer Mice: Local Adaptation in the Face of Strong Gene Flow.

The interplay of gene flow, genetic drift, and local selective pressure is a dynamic process that has been well studied from a theoretical perspective over the last century. Wright and Haldane laid the foundation for expectations under an island-continent model, demonstrating that an island-specific beneficial allele may be maintained locally if the selection coefficient is larger than the rate of migration of the ancestral allele from the continent. Subsequent extensions of this model have provided considerably more insight. Yet, connecting theoretical results with empirical data has proven challenging, owing to a lack of information on the relationship between genotype, phenotype, and fitness. Here, we examine the demographic and selective history of deer mice in and around the Nebraska Sand Hills, a system in which variation at the Agouti locus affects cryptic coloration that in turn affects the survival of mice in their local habitat. We first genotyped 250 individuals from 11 sites along a transect spanning the Sand Hills at 660,000 single nucleotide polymorphisms across the genome. Using these genomic data, we found that deer mice first colonized the Sand Hills following the last glacial period. Subsequent high rates of gene flow have served to homogenize the majority of the genome between populations on and off the Sand Hills, with the exception of the Agouti pigmentation locus. Furthermore, mutations at this locus are strongly associated with the pigment traits that are strongly correlated with local soil coloration and thus responsible for cryptic coloration.

Variation in DNA Methylation Is Not Consistently Reflected by Sociality in Hymenoptera.

Changes in gene regulation that underlie phenotypic evolution can be encoded directly in the DNA sequence or mediated by chromatin modifications such as DNA methylation. It has been hypothesized that the evolution of eusocial division of labor is associated with enhanced gene regulatory potential, which may include expansions in DNA methylation in the genomes of Hymenoptera (bees, ants, wasps, and sawflies). Recently, this hypothesis garnered support from analyses of a commonly used metric to estimate DNA methylation in silico, CpG content. Here, we test this hypothesis using direct, nucleotide-level measures of DNA methylation across nine species of Hymenoptera. In doing so, we generated new DNA methylomes for three species of interest, including one solitary and one facultatively eusocial halictid bee and a sawfly. We demonstrate that the strength of correlation between CpG content and DNA methylation varies widely among hymenopteran taxa, highlighting shortcomings in the utility of CpG content as a proxy for DNA methylation in comparative studies of taxa with sparse DNA methylomes. We observed strikingly high levels of DNA methylation in the sawfly relative to other investigated hymenopterans. Analyses of molecular evolution suggest the relatively distinct sawfly DNA methylome may be associated with positive selection on functional DNMT3 domains. Sawflies are an outgroup to all ants, bees, and wasps, and no sawfly species are eusocial. We find no evidence that either global expansions or variation within individual ortholog groups in DNA methylation are consistently associated with the evolution of social behavior.

Gregariousness does not vary with geography, developmental stage, or group relatedness in feeding redheaded pine sawfly larvae.

Aggregations are widespread across the animal kingdom, yet the underlying proximate and ultimate causes are still largely unknown. An ideal system to investigate this simple, social behavior is the pine sawfly genus Neodiprion, which is experimentally tractable and exhibits interspecific variation in larval gregariousness. To assess intraspecific variation in this trait, we characterized aggregative tendency within a single widespread species, the redheaded pine sawfly (N. lecontei). To do so, we developed a quantitative assay in which we measured interindividual distances over a 90-min video. This assay revealed minimal behavioral differences: (1) between early-feeding and late-feeding larval instars, (2) among larvae derived from different latitudes, and (3) between groups composed of kin and those composed of nonkin. Together, these results suggest that, during the larval feeding period, the benefits individuals derive from aggregating outweigh the costs and that this cost-to-benefit ratio does not vary dramatically across space (geography) or ontogeny (developmental stage). In contrast to the feeding larvae, our assay revealed a striking reduction in gregariousness following the final larval molt in N. lecontei. We also found some intriguing interspecific variation: While N. lecontei and N. maurus feeding larvae exhibit significant aggregative tendencies, feeding N. compar larvae do not aggregate at all. These results set the stage for future work investigating the proximate and ultimate mechanisms underlying developmental and interspecific variation in larval gregariousness across Neodiprion.

Evolutionary genetics of host shifts in herbivorous insects: insights from the age of genomics.

Adaptation to different host taxa is a key driver of insect diversification. Herbivorous insects are classic models for ecological and evolutionary research, but it is recent advances in sequencing, statistics, and molecular technologies that have cleared the way for investigations into the proximate genetic mechanisms underlying host shifts. In this review, we discuss how genome-scale data are revealing-at resolutions previously unimaginable-the genetic architecture of host-use traits, the causal loci underlying host shifts, and the predictability of host-use evolution. Collectively, these studies are providing novel insights into longstanding questions about host-use evolution. On the basis of this synthesis, we suggest that different host-use traits are likely to differ in their genetic architecture (number of causal loci and the nature of their genetic correlations) and genetic predictability (extent of gene or mutation reuse), indicating that any conclusions about the causes and consequences of host-use evolution will depend heavily on which host-use traits are investigated. To draw robust conclusions and identify general patterns in host-use evolution, we argue that investigation of diverse host-use traits and identification of causal genes and mutations should be the top priorities for future studies on the evolutionary genetics of host shifts.

Oviposition traits generate extrinsic postzygotic isolation between two pine sawfly species.

BACKGROUND: Although empirical data indicate that ecological speciation is prevalent in nature, the relative importance of different forms of reproductive isolation and the traits generating reproductive isolation remain unclear. To address these questions, we examined a pair of ecologically divergent pine-sawfly species: while Neodiprion pinetum specializes on a thin-needled pine (Pinus strobus), N. lecontei utilizes thicker-needled pines. We hypothesized that extrinsic postzygotic isolation is generated by oviposition traits. To test this hypothesis, we assayed ovipositor morphology, oviposition behavior, and host-dependent oviposition success in both species and in F1 and backcross females. RESULTS: Compared to N. lecontei, N. pinetum females preferred P. strobus more strongly, had smaller ovipositors, and laid fewer eggs per needle. Additionally, we observed host- and trait-dependent reductions in oviposition success in F1 and backcross females. Hybrid females that had pinetum-like host preference (P. strobus) and lecontei-like oviposition traits (morphology and egg pattern) fared especially poorly. CONCLUSIONS: Together, these data indicate that maladaptive combinations of oviposition traits in hybrids contribute to extrinsic postzygotic isolation between N. lecontei and N. pinetum, suggesting that oviposition traits may be an important driver of divergence in phytophagous insects.