Population genomic evidence of selection on structural variants in a natural hybrid zone.

Structural variants (SVs) can promote speciation by directly causing reproductive isolation or by suppressing recombination across large genomic regions. Whereas examples of each mechanism have been documented, systematic tests of the role of SVs in speciation are lacking. Here, we take advantage of long-read (Oxford nanopore) whole-genome sequencing and a hybrid zone between two Lycaeides butterfly taxa (L. melissa and Jackson Hole Lycaeides) to comprehensively evaluate genome-wide patterns of introgression for SVs and relate these patterns to hypotheses about speciation. We found >100,000 SVs segregating within or between the two hybridizing species. SVs and SNPs exhibited similar levels of genetic differentiation between species, with the exception of inversions, which were more differentiated. We detected credible variation in patterns of introgression among SV loci in the hybrid zone, with 562 of 1419 ancestry-informative SVs exhibiting genomic clines that deviated from null expectations based on genome-average ancestry. Overall, hybrids exhibited a directional shift towards Jackson Hole Lycaeides ancestry at SV loci, consistent with the hypothesis that these loci experienced more selection on average than SNP loci. Surprisingly, we found that deletions, rather than inversions, showed the highest skew towards excess ancestry from Jackson Hole Lycaeides. Excess Jackson Hole Lycaeides ancestry in hybrids was also especially pronounced for Z-linked SVs and inversions containing many genes. In conclusion, our results show that SVs are ubiquitous and suggest that SVs in general, but especially deletions, might disproportionately affect hybrid fitness and thus contribute to reproductive isolation.

A continental-scale survey of Wolbachia infections in blue butterflies reveals evidence of interspecific transfer and invasion dynamics.

Infections by maternally inherited bacterial endosymbionts, especially Wolbachia, are common in insects and other invertebrates but infection dynamics across species ranges are largely under studied. Specifically, we lack a broad understanding of the origin of Wolbachia infections in novel hosts, and the historical and geographical dynamics of infections that are critical for identifying the factors governing their spread. We used Genotype-by-Sequencing data from previous population genomics studies for range-wide surveys of Wolbachia presence and genetic diversity in North American butterflies of the genus Lycaeides. As few as one sequence read identified by assembly to a Wolbachia reference genome provided high accuracy in detecting infections in host butterflies as determined by confirmatory PCR tests, and maximum accuracy was achieved with a threshold of only 5 sequence reads per host individual. Using this threshold, we detected Wolbachia in all but 2 of the 107 sampling localities spanning the continent, with infection frequencies within populations ranging from 0% to 100% of individuals, but with most localities having high infection frequencies (mean = 91% infection rate). Three major lineages of Wolbachia were identified as separate strains that appear to represent 3 separate invasions of Lycaeides butterflies by Wolbachia. Overall, we found extensive evidence for acquisition of Wolbachia through interspecific transfer between host lineages. Strain wLycC was confined to a single butterfly taxon, hybrid lineages derived from it, and closely adjacent populations in other taxa. While the other 2 strains were detected throughout the rest of the continent, strain wLycB almost always co-occurred with wLycA. Our demographic modeling suggests wLycB is a recent invasion. Within strain wLycA, the 2 most frequent haplotypes are confined almost exclusively to separate butterfly taxa with haplotype A1 observed largely in Lycaeides melissa and haplotype A2 observed most often in Lycaeides idas localities, consistent with either cladogenic mode of infection acquisition from a common ancestor or by hybridization and accompanying mutation. More than 1 major Wolbachia strain was observed in 15 localities. These results demonstrate the utility of using resequencing data from hosts to quantify Wolbachia genetic variation and infection frequency and provide evidence of multiple colonizations of novel hosts through hybridization between butterfly lineages and complex dynamics between Wolbachia strains.

Additive genetic effects in interacting species jointly determine the outcome of caterpillar herbivory.

Plant-insect interactions are common and important in basic and applied biology. Trait and genetic variation can affect the outcome and evolution of these interactions, but the relative contributions of plant and insect genetic variation and how these interact remain unclear and are rarely subject to assessment in the same experimental context. Here, we address this knowledge gap using a recent host-range expansion onto alfalfa by the Melissa blue butterfly. Common garden rearing experiments and genomic data show that caterpillar performance depends on plant and insect genetic variation, with insect genetics contributing to performance earlier in development and plant genetics later. Our models of performance based on caterpillar genetics retained predictive power when applied to a second common garden. Much of the plant genetic effect could be explained by heritable variation in plant phytochemicals, especially saponins, peptides, and phosphatidyl cholines, providing a possible mechanistic understanding of variation in the species interaction. We find evidence of polygenic, mostly additive effects within and between species, with consistent effects of plant genotype on growth and development across multiple butterfly species. Our results inform theories of plant-insect coevolution and the evolution of diet breadth in herbivorous insects and other host-specific parasites.

Genomic time-series data show that gene flow maintains high genetic diversity despite substantial genetic drift in a butterfly species.

Effective population size affects the efficacy of selection, rate of evolution by drift and neutral diversity levels. When species are subdivided into multiple populations connected by gene flow, evolutionary processes can depend on global or local effective population sizes. Theory predicts that high levels of diversity might be maintained by gene flow, even very low levels of gene flow, consistent with species long-term effective population size, but tests of this idea are mostly lacking. Here, we show that Lycaeides butterfly populations maintain low contemporary (variance) effective population sizes (e.g. ~200 individuals) and thus evolve rapidly by genetic drift. However, populations harboured high levels of genetic diversity consistent with an effective population size several orders of magnitude larger. We hypothesized that the differences in the magnitude and variability of contemporary versus long-term effective population sizes were caused by gene flow of sufficient magnitude to maintain diversity but only subtly affect evolution on generational timescales. Consistent with this hypothesis, we detected low but nontrivial gene flow among populations. Furthermore, using short-term population-genomic time-series data, we documented patterns consistent with predictions from this hypothesis, including a weak but detectable excess of evolutionary change in the direction of the mean (migrant gene pool) allele frequencies across populations and consistency in the direction of allele frequency change over time. The documented decoupling of diversity levels and short-term change by drift in Lycaeides has implications for our understanding of contemporary evolution and the maintenance of genetic variation in the wild.

Recent hybrids recapitulate ancient hybrid outcomes.

Genomic outcomes of hybridization depend on selection and recombination in hybrids. Whether these processes have similar effects on hybrid genome composition in contemporary hybrid zones versus ancient hybrid lineages is unknown. Here we show that patterns of introgression in a contemporary hybrid zone in Lycaeides butterflies predict patterns of ancestry in geographically adjacent, older hybrid populations. We find a particularly striking lack of ancestry from one of the hybridizing taxa, Lycaeides melissa, on the Z chromosome in both the old and contemporary hybrids. The same pattern of reduced L. melissa ancestry on the Z chromosome is seen in two other ancient hybrid lineages. More generally, we find that patterns of ancestry in old or ancient hybrids are remarkably predictable from contemporary hybrids, which suggests selection and recombination affect hybrid genomes in a similar way across disparate time scales and during distinct stages of speciation and species breakdown.

Genomic evidence of genetic variation with pleiotropic effects on caterpillar fitness and plant traits in a model legume.

Plant-insect interactions are ubiquitous, and have been studied intensely because of their relevance to damage and pollination in agricultural plants, and to the ecology and evolution of biodiversity. Variation within species can affect the outcome of these interactions. Specific genes and chemicals that mediate these interactions have been identified, but genome- or metabolome-scale studies might be necessary to better understand the ecological and evolutionary consequences of intraspecific variation for plant-insect interactions. Here, we present such a study. Specifically, we assess the consequences of genome-wide genetic variation in the model plant Medicago truncatula for Lycaeides melissa caterpillar growth and survival (larval performance). Using a rearing experiment and a whole-genome SNP data set (>5 million SNPs), we found that polygenic variation in M. truncatula explains 9%-41% of the observed variation in caterpillar growth and survival. Genetic correlations among caterpillar performance and other plant traits, including structural defences and some anonymous chemical features, suggest that multiple M. truncatula alleles have pleiotropic effects on plant traits and caterpillar performance (or that substantial linkage disequilibrium exists among distinct loci affecting subsets of these traits). A moderate proportion of the genetic effect of M. truncatula alleles on L. melissa performance can be explained by the effect of these alleles on the plant traits we measured, especially leaf toughness. Taken together, our results show that intraspecific genetic variation in M. truncatula has a substantial effect on the successful development of L. melissa caterpillars (i.e., on a plant-insect interaction), and further point toward traits potentially mediating this genetic effect.

The predictability of genomic changes underlying a recent host shift in Melissa blue butterflies.

Despite accumulating evidence that evolution can be predictable, studies quantifying the predictability of evolution remain rare. Here, we measured the predictability of genome-wide evolutionary changes associated with a recent host shift in the Melissa blue butterfly (Lycaeides melissa). We asked whether and to what extent genome-wide patterns of evolutionary change in nature could be predicted (i) by comparisons among instances of repeated evolution and (ii) from SNP × performance associations in a laboratory experiment. We delineated the genetic loci (SNPs) most strongly associated with host use in two L. melissa lineages that colonized alfalfa. Whereas most SNPs were strongly associated with host use in none or one of these lineages, we detected a an approximately twofold excess of SNPs associated with host use in both lineages. Similarly, we found that host-associated SNPs in nature could also be partially predicted from SNP × performance (survival and weight) associations in a laboratory rearing experiment. But the extent of overlap, and thus degree of predictability, was somewhat reduced. Although we were able to predict (to a modest extent) the SNPs most strongly associated with host use in nature (in terms of parallelism and from the experiment), we had little to no ability to predict the direction of evolutionary change during the colonization of alfalfa. Our results show that different aspects of evolution associated with recent adaptation can be more or less predictable and highlight how stochastic and deterministic processes interact to drive patterns of genome-wide evolutionary change.

Genetic constraints on wing pattern variation in Lycaeides butterflies: A case study on mapping complex, multifaceted traits in structured populations.

Patterns of phenotypic variation within and among species can be shaped and constrained by trait genetic architecture. This is particularly true for complex traits, such as butterfly wing patterns, that consist of multiple elements. Understanding the genetics of complex trait variation across species boundaries is difficult, as it necessitates mapping in structured populations and can involve many loci with small or variable phenotypic effects. Here, we investigate the genetic architecture of complex wing pattern variation in Lycaeides butterflies as a case study of mapping multivariate traits in wild populations that include multiple nominal species or groups. We identify conserved modules of integrated wing pattern elements within populations and species. We show that trait covariances within modules have a genetic basis and thus represent genetic constraints that can channel evolution. Consistent with this, we find evidence that evolutionary changes in wing patterns among populations and species occur in the directions of genetic covariances within these groups. Thus, we show that genetic constraints affect patterns of biological diversity (wing pattern) in Lycaeides, and we provide an analytical template for similar work in other systems.

Sources of Variation in the Gut Microbial Community of Lycaeides melissa Caterpillars.

Microbes can mediate insect-plant interactions and have been implicated in major evolutionary transitions to herbivory. Whether microbes also play a role in more modest host shifts or expansions in herbivorous insects is less clear. Here we evaluate the potential for gut microbial communities to constrain or facilitate host plant use in the Melissa blue butterfly (Lycaeides melissa). We conducted a larval rearing experiment where caterpillars from two populations were fed plant tissue from two hosts. We used 16S rRNA sequencing to quantify the relative effects of sample type (frass versus whole caterpillar), diet (plant species), butterfly population and development (caterpillar age) on the composition and diversity of the caterpillar gut microbial communities, and secondly, to test for a relationship between microbial community and larval performance. Gut microbial communities varied over time (that is, with caterpillar age) and differed between frass and whole caterpillar samples. Diet (host plant) and butterfly population had much more limited effects on microbial communities. We found no evidence that gut microbe community composition was associated with caterpillar weight, and thus, our results provide no support for the hypothesis that variation in microbial community affects performance in L. melissa.

Genetic architecture, biochemical underpinnings and ecological impact of floral UV patterning.

Floral attraction traits can significantly affect pollinator visitation patterns, but adaptive evolution of these traits may be constrained by correlations with other traits. In some cases, molecular pathways contributing to floral attraction are well characterized, offering the opportunity to explore loci potentially underlying variation among individuals. Here, we quantify the range of variation in floral UV patterning (i.e. UV 'bulls-eye nectar guides) among crop and wild accessions of Brassica rapa. We then use experimental crosses to examine the genetic architecture, candidate loci and biochemical underpinnings of this patterning as well as phenotypic manipulations to test the ecological impact. We find qualitative variation in UV patterning between wild (commonly lacking UV patterns) and crop (commonly exhibiting UV patterns) accessions. Similar to the majority of crops, recombinant inbred lines (RILs) derived from an oilseed crop × WI fast-plant® cross exhibit UV patterns, the size of which varies extensively among genotypes. In RILs, we further observe strong statistical-genetic and QTL correlations within petal morphological traits and within measurements of petal UV patterning; however, correlations between morphology and UV patterning are weak or nonsignificant, suggesting that UV patterning is regulated and may evolve independently of overall petal size. HPLC analyses reveal a high concentration of sinapoyl glucose in UV-absorbing petal regions, which, in concert with physical locations of UV-trait QTLs, suggest a regulatory and structural gene as candidates underlying observed quantitative variation. Finally, insects prefer flowers with UV bulls-eye patterns over those that lack patterns, validating the importance of UV patterning in pollen-limited populations of B. rapa.

Morphological outcomes of gynandromorphism in Lycaeides butterflies (Lepidoptera: Lycaenidae).

The genitalia of male insects have been widely used in taxonomic identification and systematics and are potentially involved in maintaining reproductive isolation between species. Although sexual selection has been invoked to explain patterns of morphological variation in genitalia among populations and species, developmental plasticity in genitalia likely contributes to observed variation but has been rarely examined, particularly in wild populations. Bilateral gynandromorphs are individuals that are genetically male on one side of the midline and genetically female on the other, while mosaic gynandromorphs have only a portion of their body developing as the opposite sex. Gynandromorphs might offer unique insights into developmental plasticity because individuals experience abnormal cellular interactions at the genitalic midline. In this study, we compare the genitalia and wing patterns of gynandromorphic Anna and Melissa blue butterflies, Lycaeides anna (Edwards) (formerly L. idas anna) and L. melissa (Edwards) (Lepidoptera: Lycaenidae), to the morphology of normal individuals from the same populations. Gynandromorph wing markings all fell within the range of variation of normal butterflies; however, a number of genitalic measurements were outliers when compared with normal individuals. From these results, we conclude that the gynandromorphs' genitalia, but not wing patterns, can be abnormal when compared with normal individuals and that the gynandromorphic genitalia do not deviate developmentally in a consistent pattern across individuals. Finally, genetic mechanisms are considered for the development of gynandromorphism in Lycaeides butterflies.

The evolution of novel host use is unlikely to be constrained by trade-offs or a lack of genetic variation.

The genetic and ecological factors that shape the evolution of animal diets remain poorly understood. For herbivorous insects, the expectation has been that trade-offs exist, such that adaptation to one host plant reduces performance on other potential hosts. We investigated the genetic architecture of alternative host use by rearing individual Lycaeides melissa butterflies from two wild populations in a crossed design on two hosts (one native and one introduced) and analysing the genetic basis of differences in performance using genomic approaches. Survival during the experiment was highest when butterfly larvae were reared on their natal host plant, consistent with local adaptation. However, cross-host correlations in performance among families (within populations) were not different from zero. We found that L. melissa populations possess genetic variation for larval performance and variation in performance had a polygenic basis. We documented very few genetic variants with trade-offs that would inherently constrain diet breadth by preventing the optimization of performance across hosts. Instead, most genetic variants that affected performance on one host had little to no effect on the other host. In total, these results suggest that genetic trade-offs are not the primary cause of dietary specialization in L. melissa butterflies.

Pregnancy and lactation alter biomarkers of biotin metabolism in women consuming a controlled diet.

BACKGROUND: Biotin functions as a cofactor for several carboxylase enzymes with key roles in metabolism. At present, the dietary requirement for biotin is unknown and intake recommendations are provided as Adequate Intakes (AIs). The biotin AI for adults and pregnant women is 30 µg/d, whereas 35 µg/d is recommended for lactating women. However, pregnant and lactating women may require more biotin to meet the demands of these reproductive states. OBJECTIVE: The current study sought to quantify the impact of reproductive state on biotin status response to a known dietary intake of biotin. METHODS: To achieve this aim, we measured a panel of biotin biomarkers among pregnant (gestational week 27 at study entry; n = 26), lactating (postnatal week 5 at study entry; n = 28), and control (n = 21) women who participated in a 10- to 12-wk feeding study providing 57 µg of dietary biotin/d as part of a mixed diet. RESULTS: Over the course of the study, pregnant women excreted 69% more (vs. control; P < 0.001) 3-hydroxyisovaleric acid (3-HIA), a metabolite that accumulates during the catabolism of leucine when the activity of biotin-dependent methylcrotonyl-coenzyme A carboxylase is impaired. Interestingly, urinary excretion of 3-hydroxyisovaleryl-carnitine (3-HIA-carnitine), a downstream metabolite of 3-HIA, was 27% lower (P = 0.05) among pregnant (vs. control) women, a finding that may arise from carnitine inadequacy during gestation. No differences (P > 0.05) were detected in plasma biotin, urinary biotin, or urinary bisnorbiotin between pregnant and control women. Lactating women excreted 76% more (vs. control; P = 0.001) of the biotin catabolite bisnorbiotin, indicating that lactation accelerates biotin turnover and loss. Notably, with respect to control women, lactating women excreted 23% less (P = 0.04) urinary 3-HIA and 26% less (P = 0.05) urinary 3-HIA-carnitine, suggesting that lactation reduces leucine catabolism and that these metabolites may not be useful indicators of biotin status during lactation. CONCLUSIONS: Overall, these data demonstrate significant alterations in markers of biotin metabolism during pregnancy and lactation and suggest that biotin intakes exceeding current recommendations are needed to meet the demands of these reproductive states. This trial was registered at clinicaltrials.gov as NCT01127022.

Admixture and the organization of genetic diversity in a butterfly species complex revealed through common and rare genetic variants.

Detailed information about the geographic distribution of genetic and genomic variation is necessary to better understand the organization and structure of biological diversity. In particular, spatial isolation within species and hybridization between them can blur species boundaries and create evolutionary relationships that are inconsistent with a strictly bifurcating tree model. Here, we analyse genome-wide DNA sequence and genetic ancestry variation in Lycaeides butterflies to quantify the effects of admixture and spatial isolation on how biological diversity is organized in this group. We document geographically widespread and pervasive historical admixture, with more restricted recent hybridization. This includes evidence supporting previously known and unknown instances of admixture. The genome composition of admixed individuals varies much more among than within populations, and tree- and genetic ancestry-based analyses indicate that multiple distinct admixed lineages or populations exist. We find that most genetic variants in Lycaeides are rare (minor allele frequency <0.5%). Because the spatial and taxonomic distributions of alleles reflect demographic and selective processes since mutation, rare alleles, which are presumably younger than common alleles, were spatially and taxonomically restricted compared with common variants. Thus, we show patterns of genetic variation in this group are multifaceted, and we argue that this complexity challenges simplistic notions concerning the organization of biological diversity into discrete, easily delineated and hierarchically structured entities.

Genome divergence and the genetic architecture of barriers to gene flow between Lycaeides idas and L. Melissa.

Genome divergence during speciation is a dynamic process that is affected by various factors, including the genetic architecture of barriers to gene flow. Herein we quantitatively describe aspects of the genetic architecture of two sets of traits, male genitalic morphology and oviposition preference, that putatively function as barriers to gene flow between the butterfly species Lycaeides idas and L. melissa. Our analyses are based on unmapped DNA sequence data and a recently developed Bayesian regression approach that includes variable selection and explicit parameters for the genetic architecture of traits. A modest number of nucleotide polymorphisms explained a small to large proportion of the variation in each trait, and average genetic variant effects were nonnegligible. Several genetic regions were associated with variation in multiple traits or with trait variation within- and among-populations. In some instances, genetic regions associated with trait variation also exhibited exceptional genetic differentiation between species or exceptional introgression in hybrids. These results are consistent with the hypothesis that divergent selection on male genitalia has contributed to heterogeneous genetic differentiation, and that both sets of traits affect fitness in hybrids. Although these results are encouraging, we highlight several difficulties related to understanding the genetics of speciation.

Hybrid speciation and independent evolution in lineages of alpine butterflies.

The power of hybridization between species to generate variation and fuel adaptation is poorly understood despite long-standing interest. There is, however, increasing evidence that hybridization often generates biodiversity, including via hybrid speciation. We tested the hypothesis of hybrid speciation in butterflies occupying extreme, high-altitude habitats in four mountain ranges in western North America with an explicit, probabilistic model, and genome-wide DNA sequence data. Using this approach, in concert with ecological experiments and observations and morphological data, we document three lineages of hybrid origin. These lineages have different genome admixture proportions and distinctive trait combinations that suggest unique and independent evolutionary histories.

Geographically multifarious phenotypic divergence during speciation.

Speciation is an important evolutionary process that occurs when barriers to gene flow evolve between previously panmictic populations. Although individual barriers to gene flow have been studied extensively, we know relatively little regarding the number of barriers that isolate species or whether these barriers are polymorphic within species. Herein, we use a series of field and lab experiments to quantify phenotypic divergence and identify possible barriers to gene flow between the butterfly species Lycaeides idas and Lycaeides melissa. We found evidence that L. idas and L. melissa have diverged along multiple phenotypic axes. Specifically, we identified major phenotypic differences in female oviposition preference and diapause initiation, and more moderate divergence in mate preference. Multiple phenotypic differences might operate as barriers to gene flow, as shown by correlations between genetic distance and phenotypic divergence and patterns of phenotypic variation in admixed Lycaeides populations. Although some of these traits differed primarily between species (e.g., diapause initiation), several traits also varied among conspecific populations (e.g., male mate preference and oviposition preference).

Genomic regions with a history of divergent selection affect fitness of hybrids between two butterfly species.

Speciation is the process by which reproductively isolated lineages arise, and is one of the fundamental means by which the diversity of life increases. Whereas numerous studies have documented an association between ecological divergence and reproductive isolation, relatively little is known about the role of natural selection in genome divergence during the process of speciation. Here, we use genome-wide DNA sequences and Bayesian models to test the hypothesis that loci under divergent selection between two butterfly species (Lycaeides idas and L. melissa) also affect fitness in an admixed population. Locus-specific measures of genetic differentiation between L. idas and L. melissa and genomic introgression in hybrids varied across the genome. The most differentiated genetic regions were characterized by elevated L. idas ancestry in the admixed population, which occurs in L. idas-like habitat, consistent with the hypothesis that local adaptation contributes to speciation. Moreover, locus-specific measures of genetic differentiation (a metric of divergent selection) were positively associated with extreme genomic introgression (a metric of hybrid fitness). Interestingly, concordance of differentiation and introgression was only partial. We discuss multiple, complementary explanations for this partial concordance.

Secondary contact between Lycaeides idas and L. melissa in the Rocky Mountains: extensive admixture and a patchy hybrid zone.

Studies of hybridization have increased our understanding of the nature of species boundaries, the process of speciation, and the effects of hybridization on the evolution of populations and species. In the present study we use genetic and morphological data to determine the outcome and consequences of secondary contact and hybridization between the butterfly species Lycaeides idas and L. melissa in the Rocky Mountains. Admixture proportions estimated from structure and geographical cline analysis indicate L. idas and L. melissa have hybridized extensively in the Rocky Mountains and that reproductive isolation was insufficient to prevent introgression for much of the genome. Geographical patterns of admixture suggest that hybridization between L. idas and L. melissa has led to the formation of a hybrid zone. The hybrid zone is relatively wide, given estimates of dispersal for Lycaeides butterflies, and does not show strong evidence of cline concordance among characters. We believe the structure of the Lycaeides hybrid zone might be best explained by the patchy distribution of Lycaeides, local extinction and colonization of habitat patches, environmental variation and weak overall selection against hybrids. We found no evidence that hybridization in the Rocky Mountains has resulted in the formation of independent hybrid species, in contrast to the outcome of hybridization between L. idas and L. melissa in the Sierra Nevada. Finally, our results suggest that differences in male morphology between L. idas and L. melissa might contribute to isolation, or perhaps even that selection has favoured the spread of L. melissa male genitalia alleles.