Insect speciation rules: unifying concepts in speciation research

Fundamental Patterns of Structural Evolution Revealed by Chromosome-Length Genomes of Cactophilic Drosophila

A thorough understanding of adaptation and speciation requires model organisms with both a history of ecological and phenotypic study as well as a complete set of genomic resources. In particular, high-quality genome assemblies of ecological model organisms are needed to assess the evolution of genome structure and its role in adaptation and speciation. Here, we generate new genomes of cactophilic Drosophila, a crucial model clade for understanding speciation and ecological adaptation in xeric environments. We generated chromosome-level genome assemblies and complete annotations for seven populations across Drosophila mojavensis, Drosophila arizonae, and Drosophila navojoa. We use these data first to establish the most robust phylogeny for this clade to date, and to assess patterns of molecular evolution across the phylogeny, showing concordance with a priori hypotheses regarding adaptive genes in this system. We then show that structural evolution occurs at constant rate across the phylogeny, varies by chromosome, and is correlated with molecular evolution. These results advance the understanding of the D. mojavensis clade by demonstrating core evolutionary genetic patterns and integrating those patterns to generate new gene-level hypotheses regarding adaptation. Our data are presented in a new public database (cactusflybase.arizona.edu), providing one of the most in-depth resources for the analysis of inter- and intraspecific evolutionary genomic data. Furthermore, we anticipate that the patterns of structural evolution identified here will serve as a baseline for future comparative studies to identify the factors that influence the evolution of genome structure across taxa.

Genomic regions underlying the species-specific mating songs of green lacewings

Rapid species radiations provide insight into the process of speciation and diversification. The radiation of Chrysoperla carnea-group lacewings seems to be driven, at least in part, by their species-specific pre-mating vibrational duets. We associated genetic markers from across the genome with courtship song period in the offspring of a laboratory cross between Chrysoperla plorabunda and Chrysoperla adamsi, two species primarily differentiated by their mating songs. Two genomic regions were strongly associated with the song period phenotype. Large regions of chromosomes one and two were associated with song phenotype, as fewer recombination events occurred on these chromosomes relative to the other autosomes. Candidate genes were identified by functional annotation of proteins from the C. carnea reference genome. The majority of genes that are associated with vibrational courtship signals in other insects were found within QTL for lacewing song phenotype. Together these findings suggest that decreased recombination may be acting to keep together loci important to reproductive isolation between these species. Using wild-caught individuals from both species, we identified signals of genomic divergence across the genome. We identified several candidate genes both in song-associated regions and near divergence outliers including nonA, fruitless, paralytic, period, and doublesex. Together these findings bring us one step closer to identifying the genomic basis of a mating song trait critical to the maintenance of species boundaries in green lacewings.

Evolutionary history and diversity in the ball roller beetle Canthon cyanellus

To understand the evolutionary history of species, it is necessary to know the mechanisms for reproductive isolation, divergence-time between populations, and the relative action of the evolutionary forces (e.g., mutation, genetic drift, gene flow) within and between populations of the same, or closely related species. Although Canthon is one of the more diverse genera of neotropical beetles, insufficient research has been done to comprehend the divergent patterns that explain its speciation process. The absence of diagnostic morphological characters and the wide geographic variation of qualitative traits in Scarabaeinae obscures species delimitation, genealogical limits between populations, and its taxonomy. Canthon cyanellus is one of the best-known species in ecological and evolutionary aspects. It is a widely distributed species in the tropical forests of America. Also, the current deforestation has facilitated its incursion into open areas. Individuals from different populations have similar morphological characters but show wide variation in body color throughout their distribution, which makes it difficult to delimit the subspecies that comprise it. Recently, studies have been carried out to elucidate the pre-and postzygotic isolation mechanisms between populations and the historical biogeographical processes favoring cladogenesis events during the Pleistocene. Morphological variation of the male genitalia does not correspond to the phylogeographic structure. However, the morphological differences in one of the pieces of the endophallic sclerites have allowed a preliminary delimitation of some genetically differentiated clades. Finally, we consider that the joint analysis of traditional morphological taxonomy and phylogeography is important to understand the speciation process in the C. cyanellus complex.

The Genetic Basis of Gene Expression Divergence in Antennae of Two Closely Related Moth Species, Helicoverpa armigera and Helicoverpa assulta

The closely related species Helicoverpa armigera (H. armigera) and Helicoverpa assulta (H. assulta) have different host plant ranges and share two principal components of sex pheromones but with reversed ratios. The antennae are the main olfactory organ of insects and play a crucial role in host plant selection and mate seeking. However, the genetic basis for gene expression divergence in the antennae of the two species is unclear. We performed an allele-specific expression (ASE) analysis in the antennal transcriptomes of the two species and their F₁ hybrids, examining the connection between gene expression divergence and phenotypic differences. The results show that the proportion of genes classified as all cis was higher than that of all trans in males and reversed in females. The contribution of regulatory patterns to gene expression divergence in males was less than that in females, which explained the functional differentiation of male and female antennae. Among the five groups of F₁ hybrids, the fertile males from the cross of H. armigera female and H. assulta male had the lowest proportion of misexpressed genes, and the inferred regulatory patterns were more accurate. By using this group of F₁ hybrids, we discovered that cis-related regulations play a crucial role in gene expression divergence of sex pheromone perception-related proteins. These results are helpful for understanding how specific changes in the gene expression of olfactory-related genes can contribute to rapid evolutionary changes in important olfactory traits in closely related moths.

Reproductive isolation is a heuristic, not a measure: a commentary on Westram et al., 2022

Reproductive isolation is the heuristic basis of the biological species concept, but what is it? Westram et al. (this issue) propose that it is a measurable quantity, "barrier strength," that prevents gene flow among populations. However, their attempt to make the concept of reproductive isolation more scientific is unlikely to satisfy the diverse opinions of all evolutionary biologists. There are many different opinions about the nature of species, even under the biological species concept. Complete reproductive isolation, where gene flow is effectively zero, is regarded by some biologists as an important end point of speciation. Others, including Westram et al., argue for a more nuanced approach, and they also suggest that reproductive isolation may differ in different parts of the genome due to variation in genetic linkage to divergently selected loci. In contrast to both these approaches, we favour as a key criterion of speciation the stable coexistence of divergent populations in sympatry. Obviously, such populations must be reproductively isolated in some sense, but neither the fraction of the genome that is exchanged, nor measures of overall barrier strength acting on neutral variation will yield very precise predictions as to species status. Although an overall measure of reproductive isolation is virtually unattainable for these reasons, its early generation components, such as assortative mating, divergent selection, or hybrid inviability and sterility are readily measurable and remain informative. For example, we can make the prediction that to remain divergent in sympatry, almost all sexual species will require strong assortative mating, as well as some sort of ecological or intrinsic selection against hybrids and introgressed variants.

Using a pupal exuvia to designate the undamaged neotype of a species belonging to a complex of sibling species - the case of Aedes coluzzii (Diptera, Culicidae)

The mosquito species Aedes (Ochlerotatus) coluzzii Rioux, Guilvard & Pasteur, 1998 was distinguished from its sibling species Aedes detritus (Haliday, 1833) using an isoenzymatic method that required the destruction of the entire specimen, therefore no holotype was designated by the species authors. We aimed to designate a neotype for Ae. coluzzii from specimens collected from the type-locality and individually reared up to adult stage. Genomic DNA was extracted from pupal exuvia and ITS2 was sequenced, enabling verification of the identity of each specimen as Ae. coluzzii or Ae. detritus. Among the series of Ae. coluzzii, a male was designated as neotype and deposited in a collection. To our knowledge, this is the first time the type of a mosquito species is deposited thanks to its molecular identification from its pupal exuvia. The set of identified specimens allowed additional phylogenetic and morphologic studies.

Genome and cuticular hydrocarbon-based species delimitation shed light on potential drivers of speciation in a Neotropical ant species complex

Geographic separation that leads to the evolution of reproductive isolation between populations generally is considered the most common form of speciation. However, speciation may also occur in the absence of geographic barriers due to phenotypic and genotypic factors such as chemical cue divergence, mating signal divergence, and mitonuclear conflict. Here, we performed an integrative study based on two genome-wide techniques (3RAD and ultraconserved elements) coupled with cuticular hydrocarbon (CHC) and mitochondrial (mt) DNA sequence data, to assess the species limits within the Ectatomma ruidum species complex, a widespread and conspicuous group of Neotropical ants for which heteroplasmy (i.e., presence of multiple mtDNA variants in an individual) has been recently discovered in some populations from southeast Mexico. Our analyses indicate the existence of at least five distinct species in this complex: two widely distributed across the Neotropics, and three that are restricted to southeast Mexico and that apparently have high levels of heteroplasmy. We found that species boundaries in the complex did not coincide with geographic barriers. We therefore consider possible roles of alternative drivers that may have promoted the observed patterns of speciation, including mitonuclear incompatibility, CHC differentiation, and colony structure. Our study highlights the importance of simultaneously assessing different sources of evidence to disentangle the species limits of taxa with complicated evolutionary histories.

Limited dispersal and local adaptation promote allopatric speciation in a biodiversity hotspot

Recently diverged or diverging populations can offer unobstructed insights into early barriers to gene flow during the initial stages of speciation. The current study utilised a novel insect system (order Mantophasmatodea) to shed light on the early drivers of speciation. The members of this group have limited dispersal abilities, small allopatric distributions and strong habitat associations in the Cape Floristic Region biodiversity hotspot in South Africa. Sister taxa from the diverse family Austrophasmatidae were chosen as focal species (Karoophasma biedouwense, K. botterkloofense). Population genetics and Generalized Dissimilarity Modelling (GDM) were used to characterise spatial patterns of genetic variation and evaluate the contribution of environmental factors to population divergence and speciation. Extensive sampling confirmed the suspected allopatry of these taxa. However, hybrids were identified in a narrow region occurring between the species' distributions. Strong population structure was found over short geographic distances; particularly in K. biedouwense in which geographic distance accounted for 32% of genetic variation over a scale of 50 km (r = .56, p < .001). GDM explained 42%-78% of the deviance in observed genetic dissimilarities. Geographic distance was consistently indicated to be important for between species and within population differentiation, suggesting that limited dispersal ability may be an important neutral driver of divergence. Temperature, altitude, precipitation and vegetation were also indicated as important factors, suggesting the possible role of adaptation to local environmental conditions for species divergence. The discovery of the hybrid-zone, and the multiple allopatric species pairs in Austrophasmatidae support the idea that this could be a promising group to further our understanding of speciation modes.

Behavioral and Physiological Plasticity Provides Insights into Molecular Based Adaptation Mechanism to Strain Shift in Spodoptera frugiperda

How herbivorous insects adapt to host plants is a key question in ecological and evolutionary biology. The fall armyworm, (FAW) Spodoptera frugiperda (J.E. Smith), although polyphagous and a major pest on various crops, has been reported to have a rice and corn (maize) feeding strain in its native range in the Americas. The species is highly invasive and has recently established in China. We compared behavioral changes in larvae and adults of a corn population (Corn) when selected on rice (Rice) and the molecular basis of these adaptational changes in midgut and antennae based on a comparative transcriptome analysis. Larvae of S. frugiperda reared on rice plants continuously for 20 generations exhibited strong feeding preference for with higher larval performance and pupal weight on rice than on maize plants. Similarly, females from the rice selected population laid significantly more eggs on rice as compared to females from maize population. The most highly expressed DEGs were shown in the midgut of Rice vs. Corn. A total of 6430 DEGs were identified between the populations mostly in genes related to digestion and detoxification. These results suggest that potential adaptations for feeding on rice crops, may contribute to the current rapid spread of fall armyworm on rice crops in China and potentially elsewhere. Consistently, highly expressed DEGs were also shown in antennae; a total of 5125 differentially expressed genes (DEGs) s were identified related to the expansions of major chemosensory genes family in Rice compared to the Corn feeding population. These results not only provide valuable insight into the molecular mechanisms in host plants adaptation of S. frugiperda but may provide new gene targets for the management of this pest.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Reproductive isolation caused by azoospermia in sterile male hybrids of Drosophila

Recently diverged populations in the early stages of speciation offer an opportunity to understand mechanisms of isolation and their relative contributions. Drosophila willistoni is a tropical species with broad distribution from Argentina to the southern United States, including the Caribbean islands. A postzygotic barrier between northern populations (North America, Central America, and the northern Caribbean islands) and southern populations (South American and the southern Caribbean islands) has been recently documented and used to propose the existence of two different subspecies. Here, we identify premating isolation between populations regardless of their subspecies status. We find no evidence of postmating prezygotic isolation and proceeded to characterize hybrid male sterility between the subspecies. Sterile male hybrids transfer an ejaculate that is devoid of sperm but causes elongation and expansion of the female uterus. In sterile male hybrids, bulging of the seminal vesicle appears to impede the movement of the sperm toward the sperm pump, where sperm normally mixes with accessory gland products. Our results highlight a unique form of hybrid male sterility in Drosophila that is driven by a mechanical impediment to transfer sperm rather than by an abnormality of the sperm itself. Interestingly, this form of sterility is reminiscent of a form of infertility (azoospermia) that is caused by lack of sperm in the semen due to blockages that impede the sperm from reaching the ejaculate.

Transcriptome Profiling Revealed Potentially Critical Roles for Digestion and Defense-Related Genes in Insects' Use of Resistant Host Plants: A Case Study with Sitobion Avenae

Using host plant resistance (HPR) in management of insect pests is often environmentally friendly and suitable for sustainable development of agricultural industries. However, this strategy can be limited by rapid evolution of insect populations that overcome HPR, for which the underlying molecular factors and mechanisms are not well understood. To address this issue, we analyzed transcriptomes of two distinct biotypes of the grain aphid, Sitobion avenae (Fabricius), on wheat and barley. This analysis revealed a large number of differentially expressed genes (DEGs) between biotypes 1 and 3 on wheat and barley. The majority of them were common DEGs occurring on both wheat and barley. GO and KEGG enrichment analyses for these common DEGs demonstrated significant expression divergence between both biotypes in genes associated with digestion and defense. Top defense-related common DEGs with the most significant expression changes included three peroxidases, two UGTs (UDP-glycosyltransferase), two cuticle proteins, one glutathione S-transferases (GST), one superoxide dismutase, and one esterase, suggesting their potentially critical roles in the divergence of S. avenae biotypes. A relatively high number of specific DEGs on wheat were identified for peroxidases (9) and P450s (8), indicating that phenolic compounds and hydroxamic acids may play key roles in resistance of wheat against S. avenae. Enrichment of specific DEGs on barley for P450s and ABC transporters suggested their key roles in this aphid's detoxification against secondary metabolites (e.g., alkaloids) in barley. Our results can provide insights into the molecular factors and functions that explain biotype adaptation in insects and their use of resistant plants. This study also has significant implications for developing new resistant cultivars, developing strategies that limit rapid development of insect biotypes, and extending resistant crop cultivars' durability and sustainability in integrated management programs.

Review of Molecular Identification Techniques for Forensically Important Diptera

The medico-legal section of forensic entomology focuses on the analysis of insects associated with a corpse. Such insects are identified, and their life history characteristics are evaluated to provide information related to the corpse, such as postmortem interval and time of colonization. Forensically important insects are commonly identified using dichotomous keys, which rely on morphological characteristics. Morphological identifications can pose a challenge as local keys are not always available and can be difficult to use, especially when identifying juvenile stages. If a specimen is damaged, certain keys cannot be used for identification. In contrast, molecular identification can be a better instrument to identify forensically important insects, regardless of life stage or specimen completeness. Despite more than 20 yr since the first use of molecular data for the identification of forensic insects, there is little overlap in gene selection or phylogenetic methodology among studies, and this inconsistency reduces efficiency. Several methods such as genetic distance, reciprocal monophyly, or character-based methods have been implemented in forensic identification studies. It can be difficult to compare the results of studies that employ these different methods. Here we present a comprehensive review of the published results for the molecular identification of Diptera of forensic interest, with an emphasis on evaluating variation among studies in gene selection and phylogenetic methodology.

Fast Evolution and Lineage-Specific Gene Family Expansions of Aphid Salivary Effectors Driven by Interactions with Host-Plants

Effector proteins play crucial roles in plant-parasite interactions by suppressing plant defenses and hijacking plant physiological responses to facilitate parasite invasion and propagation. Although effector proteins have been characterized in many microbial plant pathogens, their nature and role in adaptation to host plants are largely unknown in insect herbivores. Aphids rely on salivary effector proteins injected into the host plants to promote phloem sap uptake. Therefore, gaining insight into the repertoire and evolution of aphid effectors is key to unveiling the mechanisms responsible for aphid virulence and host plant specialization. With this aim in mind, we assembled catalogues of putative effectors in the legume specialist aphid, Acyrthosiphon pisum, using transcriptomics and proteomics approaches. We identified 3,603 candidate effector genes predicted to be expressed in A. pisum salivary glands (SGs), and 740 of which displayed up-regulated expression in SGs in comparison to the alimentary tract. A search for orthologs in 17 arthropod genomes revealed that SG-up-regulated effector candidates of A. pisum are enriched in aphid-specific genes and tend to evolve faster compared with the whole gene set. We also found that a large fraction of proteins detected in the A. pisum saliva belonged to three gene families, of which certain members show evidence consistent with positive selection. Overall, this comprehensive analysis suggests that the large repertoire of effector candidates in A. pisum constitutes a source of novelties promoting plant adaptation to legumes.

Identifying genomic hotspots of differentiation and candidate genes involved in the adaptive divergence of pea aphid host races

Identifying the genomic bases of adaptation to novel environments is a long-term objective in evolutionary biology. Because genetic differentiation is expected to increase between locally adapted populations at the genes targeted by selection, scanning the genome for elevated levels of differentiation is a first step towards deciphering the genomic architecture underlying adaptive divergence. The pea aphid Acyrthosiphon pisum is a model of choice to address this question, as it forms a large complex of plant-specialized races and cryptic species, resulting from recent adaptive radiation. Here, we characterized genomewide polymorphisms in three pea aphid races specialized on alfalfa, clover and pea crops, respectively, which we sequenced in pools (poolseq). Using a model-based approach that explicitly accounts for selection, we identified 392 genomic hotspots of differentiation spanning 47.3 Mb and 2,484 genes (respectively, 9.12% of the genome size and 8.10% of its genes). Most of these highly differentiated regions were located on the autosomes, and overall differentiation was weaker on the X chromosome. Within these hotspots, high levels of absolute divergence between races suggest that these regions experienced less gene flow than the rest of the genome, most likely by contributing to reproductive isolation. Moreover, population-specific analyses showed evidence of selection in every host race, depending on the hotspot considered. These hotspots were significantly enriched for candidate gene categories that control host-plant selection and use. These genes encode 48 salivary proteins, 14 gustatory receptors, 10 odorant receptors, five P450 cytochromes and one chemosensory protein, which represent promising candidates for the genetic basis of host-plant specialization and ecological isolation in the pea aphid complex. Altogether, our findings open new research directions towards functional studies, for validating the role of these genes on adaptive phenotypes.

The Genetics of a Behavioral Speciation Phenotype in an Island System

Mating behavior divergence can make significant contributions to reproductive isolation and speciation in various biogeographic contexts. However, whether the genetic architecture underlying mating behavior divergence is related to the biogeographic history and the tempo and mode of speciation remains poorly understood. Here, we use quantitative trait locus (QTL) mapping to infer the number, distribution, and effect size of mating song rhythm variations in the crickets Laupala eukolea and Laupala cerasina, which occur on different islands (Maui and Hawaii). We then compare these results with a similar study of an independently evolving species pair that diverged within the same island. Finally, we annotate the L. cerasina transcriptome and test whether the QTL fall in functionally enriched genomic regions. We document a polygenic architecture behind the song rhythm divergence in the inter-island species pair that is remarkably similar to that previously found for an intra-island species pair in the same genus. Importantly, the QTL regions were significantly enriched for potential homologs of the genes involved in pathways that may be modulating the cricket song rhythm. These clusters of loci could constrain the spatial genomic distribution of the genetic variation underlying the cricket song variation and harbor several candidate genes that merit further study.

Testing the adaptive hypothesis of Batesian mimicry among hybridizing North American admiral butterflies

Batesian mimicry is characterized by phenotypic convergence between an unpalatable model and a palatable mimic. However, because convergent evolution may arise via alternative evolutionary mechanisms, putative examples of Batesian mimicry must be rigorously tested. Here, we used artificial butterfly facsimiles (N = 4000) to test the prediction that (1) palatable Limenitis lorquini butterflies should experience reduced predation when in sympatry with their putative model, Adelpha californica, (2) protection from predation on L. lorquini should erode outside of the geographical range of the model, and (3) mimetic color pattern traits are more variable in allopatry, consistent with relaxed selection for mimicry. We find support for these predictions, implying that this convergence is the result of selection for Batesian mimicry. Additionally, we conducted mark-recapture studies to examine the effect of mimicry and found that mimics survive significantly longer at sites where the model is abundant. Finally, in contrast to theoretical predictions, we found evidence that the Batesian model (A. californica) is protected from predation outside of its geographic range. We discuss these results considering the ongoing hybridization between L. lorquini and its sister species, L. weidemeyerii, and growing evidence that selection for mimicry predictably leads to a reduction in gene flow between nascent species.

Widespread plant specialization in the polyphagous planthopper Hyalesthes obsoletus (Cixiidae), a major vector of stolbur phytoplasma: Evidence of cryptic speciation

The stolbur phytoplasma vector Hyalesthes obsoletus is generally considered as a polyphagous species associated with numerous wild and cultivated plants. However, recent research in southeastern Europe, the distribution centre of H. obsoletus and the area of most stolbur-inflicted crop diseases, points toward specific host-plant associations of the vector, indicating specific vector-based transmission routes. Here, we study the specificity of populations associated with four host-plants using mitochondrial and nuclear genetic markers, and we evaluate the evolution of host-shifts in H. obsoletus. Host-plant use was confirmed for Convolvulus arvensis, Urtica dioica, Vitex agnus-castus and Crepis foetida. Mitochondrial genetic analysis showed sympatric occurrence of three phylogenetic lineages that were ecologically delineated by host-plant preference, but were morphologically inseparable. Nuclear data supported the existence of three genetic groups (Evanno’s ΔK(3) = 803.72) with average genetic membership probabilities > 90%. While populations associated with C. arvensis and U. dioica form a homogenous group, populations affiliated with V. agnus-castus and C. foetida constitute two independent plant-associated lineages. The geographical signal permeating the surveyed populations indicated complex diversification processes associated with host-plant selection and likely derived from post-glacial refugia in the eastern Mediterranean. This study provides evidence for cryptic species diversification within H. obsoletus sensu lato: i) consistent mitochondrial differentiation (1.1–1.5%) among host-associated populations in syntopy and in geographically distant areas, ii) nuclear genetic variance supporting mitochondrial data, and iii) average mitochondrial genetic distances among host-associated meta-populations are comparable to the most closely related, morphologically distinguishable species, i.e., Hyalesthes thracicus (2.1–3.3%).

Adaptive Divergence in a Defense Symbiosis Driven from the Top Down

Most studies of adaptive radiation in animals focus on resource competition as the primary driver of trait divergence. The roles of other ecological interactions in shaping divergent phenotypes during such radiations have received less attention. We evaluate natural enemies as primary agents of diversifying selection on the phenotypes of an actively diverging lineage of gall midges on tall goldenrod. In this system, the gall of the midge consists of a biotrophic fungal symbiont that develops on host-plant leaves and forms distinctly variable protective carapaces over midge larvae. Through field studies, we show that fungal gall morphology, which is induced by midges (i.e., it is an extended phenotype), is under directional and diversifying selection by parasitoid enemies. Overall, natural enemies disruptively select for either small or large galls, mainly along the axis of gall thickness. These results imply that predators are driving the evolution of phenotypic diversity in symbiotic defense traits in this system and that divergence in defensive morphology may provide ecological opportunities that help to fuel the adaptive radiation of this genus of midges on goldenrods. This enemy-driven phenotypic divergence in a diversifying lineage illustrates the potential importance of consumer-resource and symbiotic species interactions in adaptive radiation.

Genetic Variability of Two Leaffooted Bugs, Leptoglossus clypealis and Leptoglossus zonatus (Hemiptera: Coreidae) in the Central Valley of California

Leaffooted plant bugs (LFPBs) (Leptoglossus spp., Guérin-Méneville) (Hemiptera: Coreidae) are large seed-feeding bugs native to the Western Hemisphere. In California, several Leptoglossus spp. feed on almonds, pistachios, and pomegranate and are occasional pests. The objective of this study was to survey the different species of Leptoglossus present in almond, pistachio, and pomegranate orchards in the Central Valley of California. We used two molecular markers, amplified fragment length polymorphisms (AFLPs) and mitochondrial DNA COI, to determine the number of species or strains of each species, and to infer whether individuals of each species move and possibly interbreed with populations from the other host plants. Two species of leaffooted bugs were abundant, Leptoglossus clypealis Heidemann, and Leptoglossus zonatus (Dallas). L. clypealis was collected in almond and pistachio, while L. zonatus was found on all three host plants, but was the dominant species in pomegranate. The AFLP results indicated that L. clypealis consisted of one species, which suggests it moves between almonds and pistachios during the growing season. Mitochondrial DNA COI for L. clypealis found 1-2% divergence between sequences, and a high haplotype diversity of 0.979 with 17 haplotypes. The AFLP results for L. zonatus found two genetically divergent populations which were morphologically similar. The mtDNA COI sequences for L. zonatus were used for haplotype analysis; three haplotypes were found in California, with one haplotype shared with collections from Brazil. The importance of genetic variability and cryptic species for pest management are discussed.

Multiple strong postmating and intrinsic postzygotic reproductive barriers isolate florally diverse species of Jaltomata (Solanaceae)

Divergence in phenotypic traits often contributes to premating isolation between lineages, but could also promote isolation at postmating stages. Phenotypic differences could directly result in mechanical isolation or hybrids with maladapted traits; alternatively, when alleles controlling these trait differences pleiotropically affect other components of development, differentiation could indirectly produce genetic incompatibilities in hybrids. Here, we determined the strength of nine postmating and intrinsic postzygotic reproductive barriers among 10 species of Jaltomata (Solanaceae), including species with highly divergent floral traits. To evaluate the relative importance of floral trait diversification for the strength of these postmating barriers, we assessed their relationship to floral divergence, genetic distance, geographical context, and ecological differences, using conventional tests and a new linear-mixed modeling approach. Despite close evolutionary relationships, all species pairs showed moderate to strong isolation. Nonetheless, floral trait divergence was not a consistent predictor of the strength of isolation; instead this was best explained by genetic distance, although we found evidence for mechanical isolation in one species, and a positive relationship between floral trait divergence and fruit set isolation across species pairs. Overall, our data indicate that intrinsic postzygotic isolation is more strongly associated with genome-wide genetic differentiation, rather than floral divergence.

Difference in diel mating time contributes to assortative mating between host plant-associated populations of Chilo suppressalis

Behavioral isolation in animals can be mediated by inherent mating preferences and assortative traits, such as divergence in the diel timing of mating activity. Although divergence in the diel mating time could, in principle, promote the reproductive isolation of sympatric, conspecific populations, there is currently no unequivocal evidence of this. We conducted different mate-choice experiments to investigate the contribution of differences in diel mating activity to the reproductive isolation of the rice and water-oat populations of Chilo suppressalis. The results show that inter-population difference in diel mating activity contributes to assortative mating in these populations. In the rice population, most mating activity occurred during the first half of the scotophase, whereas in the water-oat population virtually all mating activity was confined to the second half of the scotophase. However, when the photoperiod of individuals from the water-oat population was altered to more closely align their mating activity with that of the rice population, mate choice was random. We conclude that inter-population differences in diel mating time contribute to assortative mating, and thereby the partial reproductive isolation, of these host-associated populations of C. suppressalis.

A Natural Population Derived from Species Hybridizationin the Drosophila ananassae Species Complexon Penang Island, Malaysia

We surveyed natural population of the Drosophila ananassae species complex on Penang Island, Malaysia. Analyses of phenotypic traits, chromosome arrangements, molecular markers, and reproductive isolation suggest the existence of two species: D. ananassae and D. cf. parapallidosa. Molecular marker analysis indicates that D. cf. parapallidosa carries chromosome Y and 4 introgressions from D. ananassae. Thus, D. cf. parapallidosa seems to be a hybrid descendant that recently originated from a natural D. parapallidosa♀× D. ananassae♂ cross. Furthermore, D. cf. parapallidosa behaves differently from authentic D. parapallidosa with respect to its reproductive isolation from D. ananassae. Premating isolation is usually seen in only the D. ananassae♀× D. parapallidosa♂ cross, but we observed it in crosses of both directions between D. ananassae and D. cf. parapallidosa. In addition, hybrid males from the D. ananassae♀× D. parapallidosa♂ cross are usually sterile, but they were fertile when D. ananassae♀ were mated with D. cf. parapallidosa ♂. We attempted an artificial reconstruction of the hybrid species to simulate the evolutionary process(es) that produced D. cf. parapallidosa. This is a rare case of natural hybrid population in Drosophila and may be a useful system for elucidating speciation with gene flow.

Preadult life history variation determines adult transcriptome expression

Preadult determinants of adult fitness and behaviour have been documented in a variety of organisms with complex life cycles, but little is known about expression patterns of genes underlying these adult traits. We explored the effects of differences in egg-to-adult development time on adult transcriptome and cuticular hydrocarbon variation in order to understand the nature of the genetic correlation between preadult development time and premating isolation between populations of Drosophila mojavensis reared in different host cactus environments. Transcriptome variation was analysed separately in flies reared on each host and revealed that hundreds of genes in adults were differentially expressed (FDR P < 0.05) due to development time differences. For flies reared on pitaya agria cactus, longer preadult development times caused increased expression of genes in adults enriched for ribosome production, protein metabolism, chromatin remodelling and regulation of alternate splicing and transcription. Baja California flies reared on organ pipe cactus showed fewer differentially expressed genes in adults due to longer preadult development time, but these were enriched for ATP synthesis and the TCA cycle. Mainland flies reared on organ pipe cactus with shorter development times showed increased transcription of genes enriched for mitochondria and energy production, protein synthesis and glucose metabolism: adults with longer development times had increased expression of genes enriched for adult life span, cuticle proteins and ion binding, although most differentially expressed genes were unannotated. Differences due to population, sex, mating status and their interactions were also assessed. Adult cuticular hydrocarbon profiles also showed shifts due to egg-to-adult development time and were influenced by population and mating status. These results help to explain why preadult life history variation determines subsequent expression of the adult transcriptome along with traits involved with reproductive isolation and revealed previously undocumented connections between genetic and environmental influences over the entire life cycle in this desert insect.

Reproductive isolation among allopatric Drosophila montana populations

An outstanding goal in speciation research is to trace the mode and tempo of the evolution of barriers to gene flow. Such research benefits from studying incipient speciation, in which speciation between populations has not yet occurred, but where multiple potential mechanisms of reproductive isolation (RI: i.e., premating, postmating-prezygotic (PMPZ), and postzygotic barriers) may act. We used such a system to investigate these barriers among allopatric populations of Drosophila montana. In all heteropopulation crosses we found premating (sexual) isolation, which was either symmetric or asymmetric depending on the population pair compared. Postmating isolation was particularly strong in crosses involving males from one of the study populations, and while sperm were successfully transferred, stored, and motile, we experimentally demonstrated that the majority of eggs produced were unfertilized. Thus, we identified the nature of a PMPZ incompatibility. There was no evidence of intrinsic postzygotic effects. Measures of absolute and relative strengths of pre- and postmating barriers showed that populations differed in the mode and magnitude of RI barriers. Our results indicate that incipient RI among populations can be driven by different contributions of both premating and PMPZ barriers occurring between different population pairs and without the evolution of postzygotic barriers.

Mixed signals? Morphological and molecular evidence suggest a color polymorphism in some neotropical polythore damselflies

The study of color polymorphisms (CP) has provided profound insights into the maintenance of genetic variation in natural populations. We here offer the first evidence for an elaborate wing polymorphism in the Neotropical damselfly genus Polythore, which consists of 21 described species, distributed along the eastern slopes of the Andes in South America. These damselflies display highly complex wing colors and patterning, incorporating black, white, yellow, and orange in multiple wing bands. Wing colors, along with some components of the male genitalia, have been the primary characters used in species description; few other morphological traits vary within the group, and so there are few useful diagnostic characters. Previous research has indicated the possibility of a cryptic species existing in P. procera in Colombia, despite there being no significant differences in wing color and pattern between the populations of the two putative species. Here we analyze the complexity and diversity of wing color patterns of individuals from five described Polythore species in the Central Amazon Basin of Peru using a novel suite of morphological analyses to quantify wing color and pattern: geometric morphometrics, chromaticity analysis, and Gabor wavelet transformation. We then test whether these color patterns are good predictors of species by recovering the phylogenetic relationships among the 5 species using the barcode gene (COI). Our results suggest that, while highly distinct and discrete wing patterns exist in Polythore, these “wingforms” do not represent monophyletic clades in the recovered topology. The wingforms identified as P. victoria and P. ornata are both involved in a polymorphism with P. neopicta; also, cryptic speciation may have taking place among individuals with the P. victoria wingform. Only P. aurora and P. spateri represent monophyletic species with a single wingform in our molecular phylogeny. We discuss the implications of this polymorphism, and the potential evolutionary mechanisms that could maintain it.

Host association drives genetic divergence in the bed bug, Cimex lectularius

Genetic differentiation may exist among sympatric populations of a species due to long-term associations with alternative hosts (i.e. host-associated differentiation). While host-associated differentiation has been documented in several phytophagus insects, there are far fewer cases known in animal parasites. The bed bug, Cimex lectularius, a wingless insect, represents a potential model organism for elucidating the processes involved in host-associated differentiation in animal parasites with relatively limited mobility. In conjunction with the expansion of modern humans from Africa into Eurasia, it has been speculated that bed bugs extended their host range from bats to humans in their shared cave domiciles throughout Eurasia. C. lectularius that associate with humans have a cosmopolitan distribution, whereas those associated with bats occur across Europe, often in human-built structures. We assessed genetic structure and gene flow within and among populations collected in association with each host using mtDNA, microsatellite loci and knock-down resistance gene variants. Both nuclear and mitochondrial data support a lack of significant contemporary gene flow between host-specific populations. Within locations human-associated bed bug populations exhibit limited genetic diversity and elevated levels of inbreeding, likely due to human-mediated movement, infrequent additional introduction events per infestation, and pest control. In contrast, populations within bat roosts exhibit higher genetic diversity and lower levels of relatedness, suggesting populations are stable with temporal fluctuations due to host dispersal and bug mortality. In concert with previously published evidence of morphological and behavioural differentiation, the genetic data presented here suggest C. lectularius is currently undergoing lineage divergence through host association.

Widespread host-dependent hybrid unfitness in the pea aphid species complex

Linking adaptive divergence to hybrid unfitness is necessary to understand the ecological factors contributing to reproductive isolation and speciation. To date, this link has been demonstrated in few model systems, most of which encompass ecotypes that occupy relatively early stages in the speciation process. Here we extend these studies by assessing how host-plant adaptation conditions hybrid fitness in the pea aphid, Acyrthosiphon pisum. We made crosses between and within five pea aphid biotypes adapted to different host plants and representing various stages of divergence within the complex. Performance of F1 hybrids and nonhybrids was assessed on a "universal" host that is favorable to all pea aphid biotypes in laboratory conditions. Although hybrids performed equally well as nonhybrids on the universal host, their performance was much lower than nonhybrids on the natural hosts of their parental populations. Hence, hybrids, rather than being intrinsically deficient, are maladapted to their parents' hosts. Interestingly, the impact of this maladaptation was stronger in certain hybrids from crosses involving the most divergent biotype, suggesting that host-dependent postzygotic isolation has continued to evolve late in divergence. Even though host-independent deficiencies are not excluded, hybrid maladaptation to parental hosts supports the hypothesis of ecological speciation in this complex.