Reproductive barriers between two sympatric beetle species specialized on different host plants

Coevolutionary Interactions between Sexual and Habitat Isolation during Reinforcement

Speciation often involves the evolution of multiple genetic-based barriers to gene flow (i.e., "coupling"). However, barriers may exhibit a diversity of evolutionary interactions during speciation. These dynamics are important in reinforcement, where selection may favor different prezygotic isolating barriers to avoid maladaptive hybridization. Here we study the interaction between evolution of sexual and habitat isolation. We first review the empirical literature where both barriers were explicitly considered, and then develop a population genetic model of reinforcement. Most studies of both sexual and habitat isolation were found in phytophagous insect systems. In 76% of these studies, both barriers coevolved; the remaining cases either showed only habitat isolation (21%) or only sexual isolation (3%). Our two-allele genetic mechanism model of each barrier also found that these often coevolved, but habitat isolation was generally more effective during reinforcement. Depending on the fitness of hybrids (e.g., Dobzhansky-Muller incompatibilities) and initial migration rate, these barriers could either facilitate, curtail, or have no effect on each other. This indicates that basic parameters will alter the underlying evolutionary dynamics, and thus the nature of "speciation coupling" will be highly variable in natural systems. Finally, we studied initially asymmetrical migration rates and found that populations with higher initial emigration evolved stronger habitat isolation, while populations that initially received more immigrants exhibited stronger sexual isolation. These results are in line with observations in some empirical studies, but more data is needed to test their generality.

Understanding the joint evolution of dispersal and host specialisation using phytophagous arthropods as a model group

Theory generally predicts that host specialisation and dispersal should evolve jointly. Indeed, many models predict that specialists should be poor dispersers to avoid landing on unsuitable hosts while generalists will have high dispersal abilities. Phytophagous arthropods are an excellent group to test this prediction, given extensive variation in their host range and dispersal abilities. Here, we explore the degree to which the empirical literature on this group is in accordance with theoretical predictions. We first briefly outline the theoretical reasons to expect such a correlation. We then report empirical studies that measured both dispersal and the degree of specialisation in phytophagous arthropods. We find a correlation between dispersal and levels of specialisation in some studies, but with wide variation in this result. We then review theoretical attributes of species and environment that may blur this correlation, namely environmental grain, temporal heterogeneity, habitat selection, genetic architecture, and coevolution between plants and herbivores. We argue that theoretical models fail to account for important aspects, such as phenotypic plasticity and the impact of selective forces stemming from other biotic interactions, on both dispersal and specialisation. Next, we review empirical caveats in the study of this interplay. We find that studies use different measures of both dispersal and specialisation, hampering comparisons. Moreover, several studies do not provide independent measures of these two traits. Finally, variation in these traits may occur at scales that are not being considered. We conclude that this correlation is likely not to be expected from large-scale comparative analyses as it is highly context dependent and should not be considered in isolation from the factors that modulate it, such as environmental scale and heterogeneity, intrinsic traits or biotic interactions. A stronger crosstalk between theoretical and empirical studies is needed to understand better the prevalence and basis of the correlation between dispersal and specialisation.

Integrated analysis of miRNA profiles and gut bacterial changes in Altica viridicyanea following antibiotic treatment

The gut bacteria involves in insect homeostasis by playing essential roles in host physiology, metabolism, innate immunity, and so forth. microRNAs (miRNAs) are endogenous small noncoding RNAs that posttranscriptionally regulate gene expression to affect immune or metabolic processes in insects. For several non-model insects, the available knowledge on the relationship between changes in the gut bacteria and miRNA profiles is limited. In this study, we investigated the gut bacterial diversity, composition, and function from Altica viridicyanea feeding on normal- and antibiotic-treated host plants using 16S rRNA amplicon sequencing; antibiotics have been shown to affect the body weight and development time in A. viridicyanea, suggesting that the gut bacteria of the normal sample were more diverse and abundant than those of the antibiotic-fed group, and most of them were involved in various physical functions by enrichment analysis. Furthermore, we executed small RNA transcriptome sequencing using the two experimental groups to obtain numerous sRNAs, such as piRNAs, siRNAs, and known and novel miRNAs, by data mapping and quality control, and furthermore, a total of 224 miRNAs were identified as significantly differentially expressed miRNAs, of which some DEMs and their target genes participated in immune- and metabolism-related pathways based on GO and KEGG annotation. Besides, regarding the regulatory roles of miRNA and target genes, a interaction network of DEM-target gene pairs from eight immune- or metabolism-related signaling pathways were constructed. Finally, we discovered that DEMs from above pathways were significantly positively or negatively correlated with gut bacterial alterations following antibiotic treatment. Collectively, the observations of this study expand our understanding of how the disturbance of gut bacteria affects miRNA profiles in A. viridicyanea and provide new valuable resources from extreme ranges for future studies on the adaptive evolution in insects.

The draft genome of the specialist flea beetle Altica viridicyanea (Coleoptera: Chrysomelidae)

BACKGROUND

Altica (Coleoptera: Chrysomelidae) is a highly diverse and taxonomically challenging flea beetle genus that has been used to address questions related to host plant specialization, reproductive isolation, and ecological speciation. To further evolutionary studies in this interesting group, here we present a draft genome of a representative specialist, Altica viridicyanea, the first Alticinae genome reported thus far.

RESULTS

The genome is 864.8 Mb and consists of 4490 scaffolds with a N50 size of 557 kb, which covered 98.6% complete and 0.4% partial insect Benchmarking Universal Single-Copy Orthologs. Repetitive sequences accounted for 62.9% of the assembly, and a total of 17,730 protein-coding gene models and 2462 non-coding RNA models were predicted. To provide insight into host plant specialization of this monophagous species, we examined the key gene families involved in chemosensation, detoxification of plant secondary chemistry, and plant cell wall-degradation.

CONCLUSIONS

The genome assembled in this work provides an important resource for further studies on host plant adaptation and functionally affiliated genes. Moreover, this work also opens the way for comparative genomics studies among closely related Altica species, which may provide insight into the molecular evolutionary processes that occur during ecological speciation.

Analysis of longevity in Chordata identifies species with exceptional longevity among taxa and points to the evolution of longer lifespans

Animals have a considerable variation in their longevity. This fundamental life-history trait is shaped by both intrinsic and extrinsic mortality pressures, influenced by multiple parameters including ecological variables and mode-of-life traits. Here, we examined the distribution of maximum age at multiple taxonomic ranks (class, order and family) in Chordata, and identified species with exceptional longevity within various taxa. We used a curated dataset of maximum longevity of animals from AnAge database, containing a total of 2542 chordates following our filtering criteria. We determined shapes of maximum age distributions at class, order and family taxonomic ranks, and calculated skewness values for each distribution, in R programming environment. We identified species with exceptional longevity compared to other species belonging to the same taxa, based on our definition of outliers. We collected data on ecological variables and mode-of-life traits which might possibly contribute, at least in part, to the exceptional lifespans of certain chordates. We found that 23, 12 and 4 species have exceptional longevity when we grouped chordates by their class, order and family, respectively. Almost all distributions of maximum age among taxa were positively skewed (towards increased longevity), possibly showing the emergence of longer lifespans in contrast to shorter lifespans, through the course of evolution. However, potential biases in the collection of data should be taken into account. Most of the identified species in the current study have not been previously studied in the context of animal longevity. Our analyses point that certain chordates may have evolved to have longer lifespans compared to other species belonging to the same taxa, and that among taxa, outliers in terms of maximum age have always longer lifespans, not shorter. Future research is required to understand how and why increased longevity have arose in certain species.

Patterns, Predictors, and Consequences of Dominance in Hybrids

AbstractCompared to those of their parents, are the traits of first-generation (F₁) hybrids typically intermediate, biased toward one parent, or mismatched for alternative parental phenotypes? To address this empirical gap, we compiled data from 233 crosses in which traits were measured in a common environment for two parent taxa and their F₁ hybrids. We find that individual traits in F₁s are halfway between the parental midpoint and one parental value. Considering pairs of traits together, a hybrid's bivariate phenotype tends to resemble one parent (parent bias) about 50% more than the other, while also exhibiting a similar magnitude of mismatch due to different traits having dominance in conflicting directions. Using data from an experimental field planting of recombinant hybrid sunflowers, we illustrate that parent bias improves fitness, whereas mismatch reduces fitness. Our study has three major conclusions. First, hybrids are not phenotypically intermediate but rather exhibit substantial mismatch. Second, dominance is likely determined by the idiosyncratic evolutionary trajectories of individual traits and populations. Finally, selection against hybrids likely results from selection against both intermediate and mismatched phenotypes.

Gut bacterial communities and their contribution to performance of specialist Altica flea beetles

Host plant shifts are a common mode of speciation in herbivorous insects. Although insects can evolve adaptations to successfully incorporate a new host plant, it is becoming increasingly recognized that the gut bacterial community may play a significant role in allowing insects to detoxify novel plant chemical defenses. Here, we examined differences in gut bacterial communities between Altica flea beetle species that feed on phylogenetically unrelated host plants in sympatry. We surveyed the gut bacterial communities of three closely related flea beetles from multiple locations using 16S rRNA amplicon sequencing. The results showed that the beetle species shared a high proportion (80.7%) of operational taxonomic units. Alpha-diversity indicators suggested that gut bacterial diversity did not differ among host species, whereas geography had a significant effect on bacterial diversity. In contrast, analyses of beta-diversity showed significant differences in gut bacterial composition among beetle species when we used species composition and relative abundance metrics, but there was no difference in composition when species presence/absence and phylogenetic distance indices were used. Within host beetle species, gut bacterial composition varied significantly among sites. A metagenomic functionality analysis predicted that the gut microbes had functions involved in xenobiotic biodegradation and metabolism as well as metabolism of terpenoids and polyketides. These predictions, however, did not differ among beetle host species. Antibiotic curing experiments showed that development time was significantly prolonged, and there was a significant decline in body weight of newly emerged adults in beetles lacking gut bacteria, suggesting the beetles may receive a potential benefit from the gut microbe-insect interaction. On the whole, our results suggest that although the gut bacterial community did not show clear host-specific patterns among Altica species, spatiotemporal variability is an important determinant of gut bacterial communities. Furthermore, the similarity of communities among these beetle species suggests that microbial facilitation may not be a determinant of host plant shifts in Altica.

Genetic architecture underlying host choice differentiation in the sympatric host races of Lochmaea capreae leaf beetles

Speciation in herbivorous insects has received considerable attention during the last few decades. Much of this group's diversity originates from adaptive population divergence onto different host plants, which often involves the evolution of specialized patterns of host choice behaviour. Differences in host choice often translates directly into divergence in mating sites, and therefore positive assortative mating will be created which will act as a strong barrier to gene flow. In this study, we first explored whether host choice is a genetically determined trait in the sympatric willow and birch host races of the leaf feeding beetle Lochmaea capreae, or whether larval experience influences adult host choice. Once we had established that host choice is a genetically based trait we determined its genetic architecture. To achieve this, we employed a reciprocal transplant design in which offspring from pure willow and birch cross-types, F1, F2 and backcrosses were raised on each host plant and their preference was determined upon reaching adulthood. We then applied joint-scaling analysis to uncover the genetic architecture of host preference. Our results suggest that rearing host does not have a pronounced effect on adult's host choice; rather the segregation pattern implies the existence of genetic loci affecting host choice in these host races. The joint-scaling analysis revealed that population differences in host choice are mainly influenced by the contribution of additive genetic effects and also maternally inherited cytoplasmic effects. We explore the implications of our findings for evolutionary dynamics of sympatric host race formation and speciation.

Evidence of a Genetic Basis for Differences in Parasitization Success between Strains of Avetianella longoi (Siscaro)

When the cerambycid, Phoracantha recurva, invaded California in the mid 1990’s a parasitoid wasp was imported from its native range in Australia as part of a biological control program. The wasp was later identified to be Avetianella longoi, which had already been released years earlier to control the congener longhorned beetle, Phoracantha semipunctata. Despite being recognized as the same species, the two wasps exhibited differential success on P. recurva eggs, indicating the presence of two separate strains. Here we determine if the differentiating factor between the two strains of A. longoi is a heritable genetic trait. All four pairings between the two strains were conducted, resulting in two homogenous and two heterogeneous crosses. All crosses except one produced viable F1 female offspring. F1 females were allowed to oviposit on P. recurva eggs and the survival of their offspring was compared to determine if survival can be affected by paternal contributions. The result was that the offspring of females with fathers from the second introduced strain showed significantly increased survival compared to F1 females with parents from the first introduced strain. This increased survival demonstrated that there is a heritable dominant trait that is associated with increased survival on P. recurva host eggs.

Assortative mating between two sympatric closely-related specialists: inferred from molecular phylogenetic analysis and behavioral data

Host plant shifting of phytophagous insects can lead to the formation of host associated differentiation and ultimately speciation. In some cases, host plant specificity alone acts as a nearly complete pre-mating isolating barrier among insect populations. We here test whether effective pre-mating isolation and host-independent behavioral isolation have evolved under the condition of extreme host specilization using two sympatric flea beetles with incomplete post-mating isolation under laboratory conditions. Phylogenetic analysis and coalescent simulation results showed that there is a limited interspecific gene flow, indicating effctive isolation between these species. Three types of mating tests in the absence of host plant cues showed that strong host-independent behavioral isolation has evolved between them. We conclude that almost perfect assortative mating between these two extreme host specialists results from a combination of reduced encounter rates due to differential host preference and strong sexual isolation.

Mode of inheritance of increased host acceptance in a seed beetle

Colonization of a novel plant by herbivorous insects is frequently accompanied by genetic changes that progressively improve larval or adult performance on the new host. This study examined the genetic basis of adaptation to a marginal host (lentil) by the seed beetle Callosobruchus maculatus (F.). Quasi-natural selection in the laboratory rapidly increased the tendency to oviposit on lentil. The mode of inheritance of this increase in host acceptance was determined from crosses between three lentil-adapted lines and a line maintained on the ancestral host, mung bean. In each set of crosses, females from the lentil lines laid two to three times more eggs on lentil than did females from the mung-bean line. Hybrid females consistently displayed an intermediate level of host acceptance, which did not differ between reciprocal crosses. Alleles promoting greater oviposition on lentil thus were inherited additively, with no evidence of sex-linkage or cytoplasmic effects. In a time-course study, hybrid females initially resembled the parent from the mung-bean line, as few eggs were laid on lentil during the first 24 h. However, oviposition rates on lentil after 72 h were closer to the rate observed in the lentil-line parent. Inferences about additivity vs. dominance in genes affecting oviposition may, therefore, depend on experimental protocol. Comparison with earlier work suggests that inheritance patterns observed in crosses between recently derived selection lines (as in this study) may differ from those obtained in crosses between long-divergent geographic populations.

Recent speciation in three closely related sympatric specialists: inferences using multi-locus sequence, post-mating isolation and endosymbiont data

Shifting between unrelated host plants is relatively rare for phytophagous insects, and distinct host specificity may play crucial roles in reproductive isolation. However, the isolation status and the relationship between parental divergence and post-mating isolation among closely related sympatric specialists are still poorly understood. Here, multi-locus sequence were used to estimate the relationship among three host plant–specific closely related flea beetles, Altica cirsicola, A. fragariae and A. viridicyanea (abbreviated as AC, AF and AV respectively). The tree topologies were inconsistent using different gene or different combinations of gene fragments. The relationship of AF+(AC+AV) was supported, however, by both gene tree and species tree based on concatenated data. Post-mating reproductive data on the results of crossing these three species are best interpreted in the light of a well established phylogeny. Nuclear-induced but not Wolbachia-induced unidirectional cytoplasmic incompatibility, which was detected in AC-AF and AF-AV but not in AC-AV, may also suggest more close genetic affinity between AC and AV. Prevalence of Wolbachia in these three beetles, and the endosymbiont in most individuals of AV and AC sharing a same wsp haplotype may give another evidence of AF+(AC+AV). Our study also suggested that these three flea beetles diverged in a relative short time (0.94 My), which may be the result of shifting between unrelated host plants and distinct host specificity. Incomplete post-mating isolation while almost complete lineage sorting indicated that effective pre-mating isolation among these three species should have evolved.