Bryan M. T. Brunet

Genome-wide SNPs resolve phylogenetic relationships in the North American spruce budworm (Choristoneura fumiferana) species complex

High throughput sequencing technologies have revolutionized the potential to reconcile incongruence between gene and species trees, and numerous approaches have been developed to take advantage of these advances. Genotyping-by-sequencing is becoming a regular tool for gathering phylogenetic data, yet comprehensive evaluations of phylogenetic methods using these data are sparse. Here we use multiple phylogenetic and population genetic methods for genotyping-by-sequencing data to assess species relationships in a group of forest insect pests, the spruce budworm (Choristoneura fumiferana) species complex. With few exceptions, all methods agree on the same relationships, most notably placing C. pinus as basal to the remainder of the group, rather than C. fumiferana as previously suggested. We found strong support for the monophyly of C. pinus, C. fumiferana, and C. retiniana, but more ambiguous relationships and signatures of introgression in a clade of western lineages, including C. carnana, C. lambertiana, C. occidentalis occidentalis, C. occidentalis biennis, and C. orae. This represents the most taxonomically comprehensive genomic treatment of the spruce budworm species group, which is further supported by the broad agreement among multiple methodologies.

Two's company, three's a crowd: new insights on spruce budworm species boundaries using genotyping-by-sequencing in an integrative species assessment (Lepidoptera: Tortricidae)

Species delimitation requires an assessment of varied traits that can contribute to reproductive isolation, as well as of the permanence of evolutionary differentiation among closely related lineages. Integrative taxonomy, including the combination of genome‐wide molecular data with ecological data, offers an effective approach to this issue. We use genotyping‐by‐sequencing together with a review of ecological divergence to assess the traditionally recognized species status of three closely related members of the spruce budworm species complex, Choristoneura fumiferana (Clemens), C. occidentalis Freeman (=C. freemani Razowski) and C. biennis Freeman, each of which is a major defoliator of conifer forests. We sampled a broad region of overlap between these three taxa in Alberta and British Columbia (Canada) where potential for gene flow provides a strong test of the durability of divergence among lineages. A total of 2218 single nucleotide polymorphisms (SNPs) were assayed, and patterns of differentiation were evaluated under the biological, ecological, genotypic cluster and phylogenetic species concepts. Choristoneura fumiferana was genetically distinct with substantial barriers to genetic exchange with C. occidentalis and C. biennis. Conversely, divergence between C. occidentalis and C. biennis was limited to a small subset of outlier loci and was within the range observed within any one of the taxa. Considering both population genetic and ecological patterns of divergence, C. fumiferana should continue to be recognized as a distinct species, and C. biennis (syn.n.) should be treated as a subspecies (C. occidentalis biennis Freeman, 1967) of C. occidentalis, thereby automatically establishing the nominate name C. occidentalis occidentalis Freeman, 1967 for univoltine populations of this species.

Life‐stage differences in spatial genetic structure in an irruptive forest insect: implications for dispersal and spatial synchrony

Dispersal determines the flux of individuals, energy and information and is therefore a key determinant of ecological and evolutionary dynamics. Yet, it remains difficult to quantify its importance relative to other factors. This is particularly true in cyclic populations in which demography, drift and dispersal contribute to spatio‐temporal variability in genetic structure. Improved understanding of how dispersal influences spatial genetic structure is needed to disentangle the multiple processes that give rise to spatial synchrony in irruptive species. In this study, we examined spatial genetic structure in an economically important irruptive forest insect, the spruce budworm (Choristoneura fumiferana) to better characterize how dispersal, demography and ecological context interact to influence spatial synchrony in a localized outbreak. We characterized spatial variation in microsatellite allele frequencies using 231 individuals and seven geographic locations. We show that (i) gene flow among populations is likely very high (Fst ≈ 0); (ii) despite an overall low level of genetic structure, important differences exist between adult (moth) and juvenile (larvae) life stages; and (iii) the localized outbreak is the likely source of moths captured elsewhere in our study area. This study demonstrates the potential of using molecular methods to distinguish residents from migrants and for understanding how dispersal contributes to spatial synchronization. In irruptive populations, the strength of genetic structure depends on the timing of data collection (e.g. trough vs. peak), location and dispersal. Taking into account this ecological context allows us to make more general characterizations of how dispersal can affect spatial synchrony in irruptive populations.

Characterization of EST-based SSR loci in the spruce budworm, Choristoneura fumiferana (Lepidoptera: Tortricidae)

After identifying 114 microsatellite loci from Choristoneura fumiferana expressed sequence tags, 87 loci were assayed in a panel of 11 wild-caught individuals, giving 29 polymorphic loci. Further analysis of 20 of these loci on 31 individuals collected from a single population in northern Minnesota identified 14 in Hardy–Weinberg equilibrium.

Convergent herbivory on conifers by Choristoneura moths after boreal forest formation

Mitogenomes are useful markers for phylogenetic studies across a range of taxonomic levels. Here, we focus on mitogenome variation across the tortricid moth genus Choristoneura and particularly the spruce budworm (Choristoneura fumiferana) species complex, a notorious pest group of North American conifer forests. Phylogenetic relationships of Tortricidae, representing two subfamilies, four tribes and nine genera, were analyzed using 21 mitogenomes. These included six newly-sequenced mitogenomes for species in the spruce budworm complex plus three additional Choristoneura species and 12 previously published mitogenomes from other tortricids and one from the Cossidae. We evaluated the phylogenetic informativeness of the mitogenomes and reconstructed a time-calibrated tree with fossil and secondary calibrations. We found that tortricid mitogenomes had conserved protein and ribosomal regions, and analysis of all protein-coding plus ribosomal genes together provided an efficient marker at any taxonomic rank. The time-calibrated phylogeny showed evolutionary convergence of conifer feeding within Choristoneura, with two independent lineages, the Nearctic spruce budworm complex and the Palearctic species Choristoneura murinana, both shifting onto conifers about 11 million years ago from angiosperms. These two host-plant shifts both occurred after the formation of boreal forest in the late Miocene. Haplotype diversification within the spruce budworm complex occurred in the last 4 million years, and is probably linked to the initial cooling cycles of the Northern Hemisphere in the Pliocene.

Genus delimitation, biogeography and diversification of Choristoneura Lederer (Lepidoptera: Tortricidae) based on molecular evidence: Genus Choristoneura delimitation

Widely known for pest species that include major modulators of temperate forests, the genus Choristoneura is part of the species‐rich tribe Archipini of leafroller moths (Tortricidae). Delimitation of the genus has remained unresolved because no phylogeny has included species endemic to Africa and studies have often omitted the type species of the genus. Further taxonomic confusion has been generated by the transfer of Archips occidentalis (Walsingham) to Choristoneura, creating a homonym with Choristoneura occidentalis Freeman, an important defoliator of North American forests. To define the limits of the genus, we reconstructed a phylogeny using DNA sequences for mitochondrial cytochrome oxidase subunit I and nuclear ribosomal 28S genes. Our ingroup included 23 Choristoneura species‐level taxa, complemented by a large sample of outgroups comprising 82 species of Archipini and other Tortricidae. We generated a time‐calibrated tree using fossil and secondary calibrations and we inferred biogeographic and diversification processes in Choristoneura. Our analysis recovered the genus as polyphyletic, with Archips occidentalis, Choristoneura simonyi and Choristoneura evanidana excluded from the main clade. Based on the recovered phylogenies and a redefinition, we restrict Choristoneura primarily to species with a northern hemisphere distribution. Our analysis supports A. occidentalis as the sister group of Cacoecimorpha pronubana, C. simonyi as the sister of ‘Xenotemna’ pallorana, and C. evanidana as the sister of Archips purpurana. A new combination is proposed: Archips evanidana comb.n.; the availability of ‘Xenotemna’ as a valid name is discussed and A. occidentalis is considered as an orphaned name within the Archipini. We found support for a Holarctic origin of Choristoneura about 23 Ma, followed by early divergence in the Palearctic region. The main divergence occurred at 16 Ma, with one clade in the Nearctic and another in the Palearctic. Subsequent cladogenetic events were synchronous and related to herbivorous specialization, with each clade divided into coniferophagous and polyphagous lineages. Their specialization as conifer feeders temporally matched the expansion of boreal forest during the Miocene.

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.

Genomic data indicate ubiquitous evolutionary distinctiveness among populations of California metalmark butterflies

Conservation geneticists have argued that evolutionarily significant units (ESUs) must be both genetically distinct and adaptively significant to be recognized for conservation protection. High-throughput DNA approaches can greatly increase the power to identify genetic distinctiveness, even if inferring adaptive significance remains a challenge. Here we present the first genomic evaluation of Lange’s metalmark, Apodemia mormo langei (Lepidoptera: Riodinidae), a U.S. federally endangered subspecies restricted to sand dune habitats in a single National Wildlife Refuge in California. Previous work based on very few genetic markers detected little genetic distinction for Lange’s metalmark. We use several thousand genome-wide single nucleotide polymorphisms to characterize the population structure of the A. mormo complex across California and determine if Lange’s metalmark qualifies as an ESU. We found that Lange’s metalmark is genetically identifiable, but is no more distinct than many other isolated populations across the study area. It remains unclear whether this genetic variation is adaptive, and so conservation efforts would benefit from more ecological characterization to determine conservation priorities.

Distinct sources of gene flow produce contrasting population genetic dynamics at different range boundaries of a Choristoneura budworm

Populations are often exposed to multiple sources of gene flow, but accounts are lacking of the population genetic dynamics that result from these interactions or their effects on local evolution. Using a genomic clines framework applied to 1,195 single nucleotide polymorphisms, we documented genomewide, locus‐specific patterns of introgression between Choristoneura occidentalis biennis spruce budworms and two ecologically divergent relatives, C. o. occidentalis and Choristoneura fumiferana, that it interacts with at alternate boundaries of its range. We observe contrasting hybrid indexes between the two hybrid zones, no overlap in “gene‐flow outliers” (clines showing relatively extreme extents or rates of locus‐specific introgression) and variable linkage disequilibrium among those outliers. At the same time, correlated genomewide rates of introgression between zones suggest the presence of processes common to both boundaries. These findings highlight the contrasting population genetic dynamics that can occur at separate frontiers of a single population, while also suggesting that shared patterns may frequently accompany cases of divergence‐with‐gene‐flow that involve a lineage in common. Our results point to potentially complex evolutionary outcomes for populations experiencing multiple sources of gene flow.