AFLP Linkage Map of Hybridizing Swallowtail Butterflies, Papilio glaucus and Papilio canadensis

Insights into the Structure of the Spruce Budworm (Choristoneura fumiferana) Genome, as Revealed by Molecular Cytogenetic Analyses and a High-Density Linkage Map

Genome structure characterization can contribute to a better understanding of processes such as adaptation, speciation, and karyotype evolution, and can provide useful information for refining genome assemblies. We studied the genome of an important North American boreal forest pest, the spruce budworm, Choristoneura fumiferana, through a combination of molecular cytogenetic analyses and construction of a high-density linkage map based on single nucleotide polymorphism (SNP) markers obtained through a genotyping-by-sequencing (GBS) approach. Cytogenetic analyses using fluorescence in situ hybridization methods confirmed the haploid chromosome number of n = 30 in both sexes of C. fumiferana and showed, for the first time, that this species has a WZ/ZZ sex chromosome system. Synteny analysis based on a comparison of the Bombyx mori genome and the C. fumiferana linkage map revealed the presence of a neo-Z chromosome in the latter species, as previously reported for other tortricid moths. In this neo-Z chromosome, we detected an ABC transporter C2 (ABCC2) gene that has been associated with insecticide resistance. Sex-linkage of the ABCC2 gene provides a genomic context favorable to selection and rapid spread of resistance against Bacillus thuringiensis serotype kurstaki (Btk), the main insecticide used in Canada to control spruce budworm populations. Ultimately, the linkage map we developed, which comprises 3586 SNP markers distributed over 30 linkage groups for a total length of 1720.41 cM, will be a valuable tool for refining our draft assembly of the spruce budworm genome.

The genetic architecture of ecological adaptation: intraspecific variation in host plant use by the lepidopteran crop pest Chloridea virescens

Intraspecific variation in ecologically important traits is a cornerstone of Darwin's theory of evolution by natural selection. The evolution and maintenance of this variation depends on genetic architecture, which in turn determines responses to natural selection. Some models suggest that traits with complex architectures are less likely to respond to selection than those with simple architectures, yet rapid divergence has been observed in such traits. The simultaneous evolutionary lability and genetic complexity of host plant use in the Lepidopteran subfamily Heliothinae suggest that architecture may not constrain ecological adaptation in this group. Here we investigate the response of Chloridea virescens, a generalist that feeds on diverse plant species, to selection for performance on a novel host, Physalis angulata (Solanaceae). P. angulata is the preferred host of Chloridea subflexa, a narrow specialist on the genus Physalis. In previous experiments, we found that the performance of C. subflexa on P. angulata depends on many loci of small effect distributed throughout the genome, but whether the same architecture would be involved in the generalist's adoption of P. angulata was unknown. Here we report a rapid response to selection in C. virescens for performance on P. angulata, and establish that the genetic architecture of intraspecific variation is quite similar to that of the interspecific differences in terms of the number, distribution, and effect sizes of the QTL involved. We discuss the impact of genetic architecture on the ability of Heliothine moths to respond to varying ecological selection pressures.

Molecular evidence for hybridization in Colias (Lepidoptera: Pieridae): are Colias hybrids really hybrids?

Gene flow and hybridization among species dramatically affect our understanding of the species as a biological unit, species relationships, and species adaptations. In North American Colias eurytheme and Colias eriphyle, there has been historical debate over the extent of hybridization occurring and the identity of phenotypically intermediate individuals as genetic hybrids. This study assesses the population structure of these two species to measure the extent of hybridization and the genetic identity of phenotypic intermediates as hybrids. Amplified fragment length polymorphism (AFLP) marker analysis was performed on 378 specimens collected from northern California and Nevada. Population structure was inferred using a Bayesian/Markov chain Monte Carlo method, which probabilistically assigns individuals to genetic clusters. Three genetic clusters provided the best fit for the data. C. eurytheme individuals were primarily assigned to two closely related clusters, and C. eriphyle individuals were mostly assigned to a third, more distantly related cluster. There appeared to be significant hybridization between the two species. Individuals of intermediate phenotype (putative hybrids) were found to be genetically indistinguishable from C. eriphyle, indicating that previous work based on the assumption that these intermediate forms are hybrids may warrant reconsideration.

A new Heraclides swallowtail (Lepidoptera, Papilionidae) from North America is recognized by the pattern on its neck

Heraclidesrumiko Shiraiwa & Grishin, sp. n. is described from southwestern United States, Mexico, and Central America (type locality: USA, Texas, Duval County). It is closely allied to Heraclidescresphontes (Cramer, 1777) and the two species are sympatric in central Texas. The new species is diagnosed by male genitalia and exhibits a nearly 3% difference from Heraclidescresphontes in the COI DNA barcode sequence of mitochondrial DNA. The two Heraclides species can usually be told apart by the shape and size of yellow spots on the neck, by the wing shape, and the details of wing patterns. "Western Giant Swallowtail" is proposed as the English name for Heraclidesrumiko. To stabilize nomenclature, neotype for Papiliocresphontes Cramer, 1777, an eastern United States species, is designated from Brooklyn, New York, USA; and lectotype for Papiliothoas Linnaeus, 1771 is designated from Suriname. We sequenced DNA barcodes and ID tags of nearly 400 Papilionini specimens completing coverage of all Heraclides species. Comparative analyses of DNA barcodes, genitalia, and facies suggest that Heraclidesoviedo (Gundlach, 1866), reinstated status, is a species-level taxon rather than a subspecies of Heraclidesthoas (Linnaeus, 1771); and Heraclidespallas (G. Gray, [1853]), reinstated status, with its subspecies HeraclidesPapiliobajaensis (J. Brown & Faulkner, 1992), comb. n., and Heraclidesanchicayaensis Constantino, Le Crom & Salazar, 2002, stat. n., are not conspecific with Heraclidesastyalus (Godart, 1819).

Genome-wide characterization of adaptation and speciation in tiger swallowtail butterflies using de novo transcriptome assemblies

Hybrid speciation appears to be rare in animals, yet characterization of possible examples offers to shed light on the genomic consequences of this unique phenomenon, as well as more general processes such as the role of adaptation in speciation. Here, we first generate transcriptome assemblies for a putative hybrid butterfly species, Papilio appalachiensis, its parental species, P. glaucus and P. canadensis, and an outgroup, P. polytes. Then, we use these data to infer genome-wide patterns of introgression and genomic mosaicism using both phylogenetic and population genetic approaches. Our results reveal that there is little genetic divergence among all three of the focal species, but the subset of gene trees that strongly support a specific tree topology suggest widespread sharing of genetic variation between P. appalachiensis and both parental species, likely as a result of hybrid speciation. We also find evidence for substantial shared genetic variation between P. glaucus and P. canadensis, which may be due to gene flow or ancestral variation. Consistent with previous work, we show that P. applachiensis is more similar to P. canadensis at Z-linked genes and more similar to P. glaucus at mitochondrial genes. We also identify a variety of targets of adaptive evolution, which appear to be enriched for traits that are likely to be important in the evolution of this butterfly system, such as pigmentation, hormone sensitivity, developmental processes, and cuticle formation. Overall, our results provide a genome-wide portrait of divergence and introgression associated with adaptation and speciation in an iconic butterfly radiation.

Genetic mapping of two components of reproductive isolation between two sibling species of moths, Ostrinia nubilalis and O. scapulalis

We report the quantitative trait loci (QTL) mapping of reproductive isolation traits between Ostrinia nubilalis (the European corn borer) and its sibling species O. scapulalis (the Adzuki bean borer), focusing on two traits: mating isolation (mi) and pheromone production (Pher). Four genetic maps were generated from two backcross families, with two maps (one chromosomal map and one linkage map) per backcross. We located 165-323 AFLP markers on these four maps, resulting in the identification of 27-31 linkage groups, depending on the map considered. No-choice mating experiments with the offspring of each backcross led to the detection of at least two QTLs for mi in different linkage groups. QTLs underlying Pher were located in a third linkage group. The Z heterochromosome was identified by a specific marker (Tpi) and did not carry any of these QTLs. Finally, we considered the global divergence between the two sibling species, distortions of segregation throughout the genome, and the location and effect of mi and Pher QTLs in light of the known candidate genes for reproductive isolation within the genus Ostrinia and, more broadly, in phytophagous insects.

Linkage map of the peppered moth, Biston betularia (Lepidoptera, Geometridae): a model of industrial melanism

We have constructed a linkage map for the peppered moth (Biston betularia), the classical ecological genetics model of industrial melanism, aimed both at localizing the network of loci controlling melanism and making inferences about chromosome dynamics. The linkage map, which is based primarily on amplified fragment length polymorphisms (AFLPs) and genes, consists of 31 linkage groups (LGs; consistent with the karyotype). Comparison with the evolutionarily distant Bombyx mori suggests that the gene content of chromosomes is highly conserved. Gene order is conserved on the autosomes, but noticeably less so on the Z chromosome, as confirmed by physical mapping using bacterial artificial chromosome fluorescence in situ hybridization (BAC-FISH). Synteny mapping identified three pairs of B. betularia LGs (11/29, 23/30 and 24/31) as being orthologous to three B. mori chromosomes (11, 23 and 24, respectively). A similar finding in an outgroup moth (Plutella xylostella) indicates that the B. mori karyotype (n=28) is a phylogenetically derived state resulting from three chromosome fusions. As with other Lepidoptera, the B. betularia W chromosome consists largely of repetitive sequence, but exceptionally we found a W homolog of a Z-linked gene (laminin A), possibly resulting from ectopic recombination between the sex chromosomes. The B. betularia linkage map, featuring the network of known melanization genes, serves as a resource for melanism research in Lepidoptera. Moreover, its close resemblance to the ancestral lepidopteran karyotype (n=31) makes it a useful reference point for reconstructing chromosome dynamic events and ancestral genome architectures. Our study highlights the unusual evolutionary stability of lepidopteran autosomes; in contrast, higher rates of intrachromosomal rearrangements support a special role of the Z chromosome in adaptive evolution and speciation.