Modulation of the heat shock response is associated with acclimation to novel temperatures but not adaptation to climatic variation in the ants Aphaenogaster picea and A. rudis.

Ecological diversification into thermally divergent habitats can push species toward their physiological limits, requiring them to accommodate temperature extremes through plastic or evolutionary changes that increase persistence under the local thermal regime. One way to withstand thermal stress is to increase production of heat shock proteins, either by maintaining higher baseline abundance within cells or by increasing the magnitude of induction in response to heat stress. We evaluated whether environmental variation was associated with expression of three heat shock protein genes in two closely-related species of woodland ant, Aphaenogaster picea and A. rudis. We compared adult workers from colonies collected from 25 sites across their geographic ranges. Colonies were maintained at two different laboratory temperatures, and tested for the independent effects of environment, phylogeny, and acclimation temperature on baseline and heat-induced gene expression. The annual maximum temperature at each collection site (Tmax) was not a significant predictor of either baseline expression or magnitude of induction of any of the heat shock protein genes tested. A phylogenetic effect was detected only for basal expression of Hsp40, which was lower in the most southern populations of A. rudis and higher in a mid-range population of possible hybrid ancestry. In contrast, a higher acclimation temperature significantly increased baseline expression of Hsc70-4, and increased induction of Hsp40 and Hsp83. Thus, physiological acclimation to temperature variation appears to involve modulation of the heat shock response, whereas other mechanisms are likely to be responsible for evolutionary shifts in thermal performance associated with large-scale climate gradients.

How Do Genomes Create Novel Phenotypes? Insights from the Loss of the Worker Caste in Ant Social Parasites.

A central goal of biology is to uncover the genetic basis for the origin of new phenotypes. A particularly effective approach is to examine the genomic architecture of species that have secondarily lost a phenotype with respect to their close relatives. In the eusocial Hymenoptera, queens and workers have divergent phenotypes that may be produced via either expression of alternative sets of caste-specific genes and pathways or differences in expression patterns of a shared set of multifunctional genes. To distinguish between these two hypotheses, we investigated how secondary loss of the worker phenotype in workerless ant social parasites impacted genome evolution across two independent origins of social parasitism in the ant genera Pogonomyrmex and Vollenhovia. We sequenced the genomes of three social parasites and their most-closely related eusocial host species and compared gene losses in social parasites with gene expression differences between host queens and workers. Virtually all annotated genes were expressed to some degree in both castes of the host, with most shifting in queen-worker bias across developmental stages. As a result, despite >1 My of divergence from the last common ancestor that had workers, the social parasites showed strikingly little evidence of gene loss, damaging mutations, or shifts in selection regime resulting from loss of the worker caste. This suggests that regulatory changes within a multifunctional genome, rather than sequence differences, have played a predominant role in the evolution of social parasitism, and perhaps also in the many gains and losses of phenotypes in the social insects.

Global invasion history of the tropical fire ant: a stowaway on the first global trade routes.

Biological invasions are largely thought to be contemporary, having recently increased sharply in the wake of globalization. However, human commerce had already become global by the mid-16th century when the Spanish connected the New World with Europe and Asia via their Manila galleon and West Indies trade routes. We use genetic data to trace the global invasion of one of the world's most widespread and invasive pest ants, the tropical fire ant, Solenopsis geminata. Our results reveal a pattern of introduction of Old World populations that is highly consistent with historical trading routes suggesting that Spanish trade introduced the tropical fire ant to Asia in the 16th century. We identify southwestern Mexico as the most likely source for the invasive populations, which is consistent with the use of Acapulco as the major Spanish port on the Pacific Ocean. From there, the Spanish galleons brought silver to Manila, which served as a hub for trade with China. The genetic data document a corresponding spread of S. geminata from Mexico via Manila to Taiwan and from there, throughout the Old World. Our descriptions of the worldwide spread of S. geminata represent a rare documented case of a biological invasion of a highly invasive and globally distributed pest species due to the earliest stages of global commerce.

Catch me if you can: Under-detection of Trypanosoma cruzi (Kinetoplastea: Trypanosomatida) infections in Triatoma dimidiata s.l. (Hemiptera: Reduviidae) from Central America.

Assays for parasite detection in insect vectors provide important information for disease control. American Trypanosomiasis (Chagas disease) is the most devastating vector-borne illness and the fourth most common in Central America behind HIV/AIDS and acute respiratory and diarrheal infections (Peterson et al., 2019). Under-detection of parasites is a general problem which may be influenced by parasite genetic variation; however, little is known about the genetic variation of the Chagas parasite, especially in this region. In this study we compared six assays for detecting the Chagas parasite, Trypanosoma cruzi: genomic reduced representation sequencing (here referred to as genotype-by-sequencing or GBS), two with conventional PCR (i.e., agarose gel detection), two with qPCR, and microscopy. Our results show that, compared to GBS genomic analysis, microscopy and PCR under-detected T. cruzi in vectors from Central America. Of 94 samples, 44% (50/94) were positive based on genomic analysis. The lowest detection, 9% (3/32) was in a subset assayed with microscopy. Four PCR assays, two with conventional PCR and two with qPCR showed intermediate levels of detection. Both qPCR tests and one conventional PCR test targeted the 195 bp repeat of satellite DNA while the fourth test targeted the 18S gene. Statistical analyses of the genomic and PCR results indicate that the PCR assays significantly under detect infections of Central American T. cruzi genotypes.

Complex hybrid origin of genetic caste determination in harvester ants.

Caste differentiation and division of labour are the hallmarks of insect societies and at the root of their ecological success. Kin selection predicts that caste determination should result from environmentally induced differences in gene expression, a prediction largely supported by empirical data. However, two exceptional cases of genetically determined caste differentiation have recently been found in harvester ants. Here we show that genetic caste determination evolved in these populations after complex hybridization events. We identified four distinct genetic lineages, each consisting of unique blends of the genomes of the parental species, presumably Pogonomyrmex barbatus and P. rugosus. Crosses between lineages H1 and H2 and between J1 and J2 give rise to workers, whereas queens develop from within-lineage matings. Although historical gene flow is evident, genetic exchange among lineages and between lineages and the parental species no longer occurs. This unusual system of caste determination seems to be evolutionarily stable.

Novel Evolutionary Algorithm Identifies Interactions Driving Infestation of Triatoma dimidiata, a Chagas Disease Vector.

Chagas disease is a lethal, neglected tropical disease. Unfortunately, aggressive insecticide-spraying campaigns have not been able to eliminate domestic infestation of Triatoma dimidiata, the native vector in Guatemala. To target interventions toward houses most at risk of infestation, comprehensive socioeconomic and entomologic surveys were conducted in two towns in Jutiapa, Guatemala. Given the exhaustively large search space associated with combinations of risk factors, traditional statistics are limited in their ability to discover risk factor interactions. Two recently developed statistical evolutionary algorithms, specifically designed to accommodate risk factor interactions and heterogeneity, were applied to this large combinatorial search space and used in tandem to identify sets of risk factor combinations associated with infestation. The optimal model includes 10 risk factors in what is known as a third-order disjunctive normal form (i.e., infested households have chicken coops AND deteriorated bedroom walls OR an accumulation of objects AND dirt floors AND total number of occupants ≥ 5 AND years of electricity ≥ 5 OR poor hygienic condition ratings AND adobe walls AND deteriorated walls AND dogs). Houses with dirt floors and deteriorated walls have been reported previously as risk factors and align well with factors currently targeted by Ecohealth interventions to minimize infestation. However, the tandem evolutionary algorithms also identified two new socioeconomic risk factors (i.e., households having many occupants and years of electricity ≥ 5). Identifying key risk factors may help with the development of new Ecohealth interventions and/or reduce the survey time needed to identify houses most at risk.

Draft Aphaenogaster genomes expand our view of ant genome size variation across climate gradients.

Given the abundance, broad distribution, and diversity of roles that ants play in many ecosystems, they are an ideal group to serve as ecosystem indicators of climatic change. At present, only a few whole-genome sequences of ants are available (19 of >16,000 species), mostly from tropical and sub-tropical species. To address this limited sampling, we sequenced genomes of temperate-latitude species from the genus Aphaenogaster, a genus with important seed dispersers. In total, we sampled seven colonies of six species: Aphaenogaster ashmeadi, Aphaenogaster floridana, Aphaenogaster fulva, Aphaenogaster miamiana, Aphaenogaster picea, and Aphaenogaster rudis. The geographic ranges of these species collectively span eastern North America from southern Florida to southern Canada, which encompasses a latitudinal gradient in which many climatic variables are changing rapidly. For the six genomes, we assembled an average of 271,039 contigs into 47,337 scaffolds. The Aphaenogaster genomes displayed high levels of completeness with 96.1% to 97.6% of Hymenoptera BUSCOs completely represented, relative to currently sequenced ant genomes which ranged from 88.2% to 98.5%. Additionally, the mean genome size was 370.5 Mb, ranging from 310.3 to 429.7, which is comparable to that of other sequenced ant genomes (212.8-396.0 Mb) and flow cytometry estimates (210.7-690.4 Mb). In an analysis of currently sequenced ant genomes and the new Aphaenogaster sequences, we found that after controlling for both spatial autocorrelation and phylogenetics ant genome size was marginally correlated with sample site climate similarity. Of all examined climate variables, minimum temperature, and annual precipitation had the strongest correlations with genome size, with ants from locations with colder minimum temperatures and higher levels of precipitation having larger genomes. These results suggest that climate extremes could be a selective force acting on ant genomes and point to the need for more extensive sequencing of ant genomes.

Uncovering vector, parasite, blood meal and microbiome patterns from mixed-DNA specimens of the Chagas disease vector Triatoma dimidiata.

Chagas disease, considered a neglected disease by the World Health Organization, is caused by the protozoan parasite Trypanosoma cruzi, and transmitted by >140 triatomine species across the Americas. In Central America, the main vector is Triatoma dimidiata, an opportunistic blood meal feeder inhabiting both domestic and sylvatic ecotopes. Given the diversity of interacting biological agents involved in the epidemiology of Chagas disease, having simultaneous information on the dynamics of the parasite, vector, the gut microbiome of the vector, and the blood meal source would facilitate identifying key biotic factors associated with the risk of T. cruzi transmission. In this study, we developed a RADseq-based analysis pipeline to study mixed-species DNA extracted from T. dimidiata abdomens. To evaluate the efficacy of the method across spatial scales, we used a nested spatial sampling design that spanned from individual villages within Guatemala to major biogeographic regions of Central America. Information from each biotic source was distinguished with bioinformatics tools and used to evaluate the prevalence of T. cruzi infection and predominant Discrete Typing Units (DTUs) in the region, the population genetic structure of T. dimidiata, gut microbial diversity, and the blood meal history. An average of 3.25 million reads per specimen were obtained, with approximately 1% assigned to the parasite, 20% to the vector, 11% to bacteria, and 4% to putative blood meals. Using a total of 6,405 T. cruzi SNPs, we detected nine infected vectors harboring two distinct DTUs: TcI and a second unidentified strain, possibly TcIV. Vector specimens were sufficiently variable for population genomic analyses, with a total of 25,710 T. dimidiata SNPs across all samples that were sufficient to detect geographic genetic structure at both local and regional scales. We observed a diverse microbiotic community, with significantly higher bacterial species richness in infected T. dimidiata abdomens than those that were not infected. Unifrac analysis suggests a common assemblage of bacteria associated with infection, which co-occurs with the typical gut microbial community derived from the local environment. We identified vertebrate blood meals from five T. dimidiata abdomens, including chicken, dog, duck and human; however, additional detection methods would be necessary to confidently identify blood meal sources from most specimens. Overall, our study shows this method is effective for simultaneously generating genetic data on vectors and their associated parasites, along with ecological information on feeding patterns and microbial interactions that may be followed up with complementary approaches such as PCR-based parasite detection, 18S eukaryotic and 16S bacterial barcoding.

Reproductive division of labor between hybrid and nonhybrid offspring in a fire ant hybrid zone.

Interspecific hybridization can often impose a substantial fitness cost due to reduced hybrid viability or fecundity. In social insects, however, such costs disproportionately impact reproductive offspring, whereas hybrids who become sterile workers can be functional, and even beneficial, colony members. Genomic imprinting of the paternal genome in reproductive, but not worker female offspring has been proposed as a mechanism to avoid genomic incompatibilities in hybrid queens in a hybrid zone between two fire ant species, Solenopsis geminata and S. xyloni. A study of allozyme variation demonstrated differences between the worker caste displaying a hybrid phenotype, and the winged queen caste displaying only the mother's phenotype. In this study, we investigate whether these differences are caused by genomic imprinting or genetic differences between castes by comparing variability of proteins to that of microsatellite markers. Workers and winged queens differed genetically at both classes of marker, indicating that allozyme differences were caused by underlying genetic differences between castes rather than differences in gene expression due to imprinting. Workers were F1 S. geminata x S. xyloni hybrids, whereas nearly all winged queens were of pure S. xyloni ancestry. Thus, S. xyloni within the hybrid zone appears to have evolved social hybridogenesis, in which the loss of worker potential in pure-species offspring necessitates hybridization for worker production, but prevents hybrids from being represented in the reproductive caste.

Extreme genetic differences between queens and workers in hybridizing Pogonomyrmex harvester ants.

The process of reproductive caste determination in eusocial insect colonies is generally understood to be mediated by environmental, rather than genetic factors. We present data demonstrating unexpected genetic differences between reproductive castes in a variant of the rough harvester ant, Pogonomyrmex rugosus var. fuscatus. Across multiple loci, queens were consistently more homozygous than expected, while workers were more heterozygous. Adult colony queens were divided into two highly divergent genetic groups, indicating the presence of two cryptic species, rather than a single population. The observed genetic differences between castes reflect differential representation of heterospecific and conspecific patrilines in these offspring groups. All workers were hybrids; by contrast, winged queens were nearly all pure-species. The complete lack of pure-species workers indicates a loss of worker potential in pure-species female offspring. Hybrids appear to be bipotential, but do not normally develop into reproductives because they are displaced by pure-species females in the reproductive pool. Genetic differences between reproductive castes are expected to be rare in non-hybridizing populations, but within hybrid zones they may be evolutionarily stable and thus much more likely to occur.