Transcriptomic Analysis Pipeline (TAP) for quality control and functional assessment of transcriptomes.

Background: RNA-sequencing (RNA-seq) has revolutionized the exploration of biological mechanisms, shedding light on the roles of non-coding RNAs, including long non-coding RNAs (lncRNAs), across various biological processes, including stress responses. Despite these advancements, there remains a gap in our understanding of the implications of different RNA-seq library protocols on comprehensive lncRNA expression analysis, particularly in non-mammalian organisms. Results: In this study, we sought to bridge this knowledge gap by investigating lncRNA expression patterns in Drosophila melanogaster under thermal stress conditions. To achieve this, we conducted a comparative analysis of two RNA-seq library protocols: polyA + RNA capture and rRNA-depletion. Our approach involved the development and application of a Transcriptome Analysis Pipeline (TAP) designed to systematically assess both the technical and functional dimensions of RNA-seq, facilitating a robust comparison of these library protocols. Our findings underscore the efficacy of the polyA + protocol in capturing the majority of expressed lncRNAs within the Drosophila melanogaster transcriptome. In contrast, rRNA-depletion exhibited limited advantages in the context of D. melanogaster studies. Notably, the polyA + protocol demonstrated superior performance in terms of usable read yield and the accurate detection of splice junctions. Conclusions: Our study introduces a versatile transcriptomic analysis pipeline, TAP, designed to uniformly process RNA-seq data from any organism with a reference genome. It also highlights the significance of selecting an appropriate RNA-seq library protocol tailored to the specific research context. Background: Advances in next generation sequencing (NGS) technologies enable the comprehensive analysis of genetic sequences of organisms in a relatively cost-effective manner [1, 2]. Among these technologies, RNA-sequencing (RNA-seq) has emerged as a preeminent method to study fundamental biological mechanisms at the level of cells, tissues, and whole organisms. RNA-seq enables the detection and quantification of various RNA populations, including messenger RNA (mRNA) and various species of non-coding RNA, such as long non-coding RNA (lncRNA), as well as an assessment of features including splice junctions in RNA.

Population genomics supports multiple hybrid zone origins of socially hybridogenetic lineages of Pogonomyrmex harvester ants.

Reproductive division of labor in the social insects is typically determined by environmental cues; however, genetic effects on caste have been discovered in a growing set of ant taxa. An extreme form of genetic caste determination is "social hybridogenesis," in which co-occurring genetic lineages obligately interbreed to produce workers, whereas daughter queens are of pure-lineage ancestry. In this study, we tested the hypothesis that social hybridogenesis in the genus Pogonomyrmex resulted from one or more interspecific hybridization events, and if so, whether individual lineages were of hybrid ancestry. We reconstructed evolutionary relationships of four lineage pairs to populations of two closely related non-hybridogenetic species, Pogonomyrmex barbatus and Pogonomyrmex rugosus, using nuclear SNP loci and mitochondrial sequencing. The nuclear phylogeny supported a hybridization hypothesis, with one member of each pair nested within P. rugosus, whereas the other was nested within P. barbatus. The source populations corresponded with two distinct geographic areas at the eastern and western edges of a zone of contact. Relatively little gene flow was detected between interbreeding lineages, either historically or currently. This suggests that shifts in reproductive caste determination may reinforce reproductive incompatibility, in a manner similar to the evolution of hybridogenesis in nonsocial systems.

Insights from a comprehensive study of Trypanosoma cruzi: A new mitochondrial clade restricted to North and Central America and genetic structure of TcI in the region.

More than 100 years since the first description of Chagas Disease and with over 29,000 new cases annually due to vector transmission (in 2010), American Trypanosomiasis remains a Neglected Tropical Disease (NTD). This study presents the most comprehensive Trypanosoma cruzi sampling in terms of geographic locations and triatomine species analyzed to date and includes both nuclear and mitochondrial genomes. This addresses the gap of information from North and Central America. We incorporate new and previously published DNA sequence data from two mitochondrial genes, Cytochrome oxidase II (COII) and NADH dehydrogenase subunit 1 (ND1). These T. cruzi samples were collected over a broad geographic range including 111 parasite DNA samples extracted from triatomines newly collected across North and Central America, all of which were infected with T. cruzi in their natural environment. In addition, we present parasite reduced representation (Restriction site Associated DNA markers, RAD-tag) genomic nuclear data combined with the mitochondrial gene sequences for a subset of the triatomines (27 specimens) collected from Guatemala and El Salvador. Our mitochondrial phylogenetic reconstruction revealed two of the major mitochondrial lineages circulating across North and Central America, as well as the first ever mitochondrial data for TcBat from a triatomine collected in Central America. Our data also show that within mtTcIII, North and Central America represent an independent, distinct clade from South America, named here as mtTcIIINA-CA, geographically restricted to North and Central America. Lastly, the most frequent lineage detected across North and Central America, mtTcI, was also an independent, distinct clade from South America, noted as mtTcINA-CA. Furthermore, nuclear genome data based on Single Nucleotide Polymorphism (SNP) showed genetic structure of lineage TcI from specimens collected in Guatemala and El Salvador supporting the hypothesis that genetic diversity at a local scale has a geographical component. Our multiscale analysis contributes to the understanding of the independent and distinct evolution of T. cruzi lineages in North and Central America regions.

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.

Integrating GWAS and Transcriptomics to Identify the Molecular Underpinnings of Thermal Stress Responses in Drosophila melanogaster.

Thermal tolerance of an organism depends on both the ability to dynamically adjust to a thermal stress and preparatory developmental processes that enhance thermal resistance. However, the extent to which standing genetic variation in thermal tolerance alleles influence dynamic stress responses vs. preparatory processes is unknown. Here, using the model species Drosophila melanogaster, we used a combination of Genome Wide Association mapping (GWAS) and transcriptomic profiling to characterize whether genes associated with thermal tolerance are primarily involved in dynamic stress responses or preparatory processes that influence physiological condition at the time of thermal stress. To test our hypotheses, we measured the critical thermal minimum (CTmin) and critical thermal maximum (CTmax) of 100 lines of the Drosophila Genetic Reference Panel (DGRP) and used GWAS to identify loci that explain variation in thermal limits. We observed greater variation in lower thermal limits, with CTmin ranging from 1.81 to 8.60°C, while CTmax ranged from 38.74 to 40.64°C. We identified 151 and 99 distinct genes associated with CTmin and CTmax, respectively, and there was strong support that these genes are involved in both dynamic responses to thermal stress and preparatory processes that increase thermal resistance. Many of the genes identified by GWAS were involved in the direct transcriptional response to thermal stress (72/151 for cold; 59/99 for heat), and overall GWAS candidates were more likely to be differentially expressed than other genes. Further, several GWAS candidates were regulatory genes that may participate in the regulation of stress responses, and gene ontologies related to development and morphogenesis were enriched, suggesting many of these genes influence thermal tolerance through effects on development and physiological status. Overall, our results suggest that thermal tolerance alleles can influence both dynamic plastic responses to thermal stress and preparatory processes that improve thermal resistance. These results also have utility for directly comparing GWAS and transcriptomic approaches for identifying candidate genes associated with thermal tolerance.

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.

Trade-Offs in Cold Resistance at the Northern Range Edge of the Common Woodland Ant Aphaenogaster picea (Formicidae).

Geographic variation in low temperatures at poleward range margins of terrestrial species often mirrors population variation in cold resistance, suggesting that range boundaries may be set by evolutionary constraints on cold physiology. The northeastern woodland ant Aphaenogaster picea occurs up to approximately 45°N in central Maine. We combined presence/absence surveys with classification tree analysis to characterize its northern range limit and assayed two measures of cold resistance operating on different timescales to determine whether and how marginal populations adapt to environmental extremes. The range boundary of A. picea was predicted primarily by temperature, but low winter temperatures did not emerge as the primary correlate of species occurrence. Low summer temperatures and high seasonal variability predicted absence above the boundary, whereas high mean annual temperature (MAT) predicted presence in southern Maine. In contrast, assays of cold resistance across multiple sites were consistent with the hypothesis of local cold adaptation at the range edge: among populations, there was a 4-min reduction in chill coma recovery time across a 2° reduction in MAT. Baseline resistance and capacity for additional plastic cold hardening shifted in opposite directions, with hardening capacity approaching zero at the coldest sites. This trade-off between baseline resistance and cold-hardening capacity suggests that populations at range edges may adapt to colder temperatures through genetic assimilation of plastic responses, potentially constraining further adaptation and range expansion.

Residual survival and local dispersal drive reinfestation by Triatoma dimidiata following insecticide application in Guatemala.

Chagas disease is caused by the protozoan parasite Trypanosoma cruzi and transmitted by triatomine insect vectors. In Guatemala, insecticide spraying is an integral part of management of the main vector, Triatoma dimidiata. Spraying typically has low efficacy, which may be due to incomplete elimination from infested houses, within-village dispersal, or influx from other villages or sylvan environments. To evaluate how these mechanisms contribute to reinfestation, we conducted a time-course analysis of T. dimidiata infestation, abundance and household genetic structure in two nearby villages in Jutiapa, Guatemala; houses in the first village were surveyed, treated with insecticide if infested and then re-surveyed at eight and 22 months following spraying, while the second village served as an untreated control to quantify changes associated with seasonal dispersal. Insects were genotyped at 2-3000 SNP loci for kinship and population genetic analyses. Insecticide application reduced overall infestation and abundance, while the untreated village was stable over time. Nevertheless, within two years 35.5% of treated houses were reinfested and genetic diversity had largely recovered. Insects collected from reinfested houses post-spraying were most closely related to pre-spray collections from the same house, suggesting that infestations had not been fully eliminated. Immigration by unrelated insects was also detected within a year of spraying; when it occurred, dispersal was primarily local from neighboring houses. Similar dispersal patterns were observed following the annual dispersal season in the untreated village, with high-infestation houses serving as sources for neighboring homes. Our findings suggest that the efficacy of pyrethroid application is rapidly diminished by both within-house breeding by survivors and annual cycles of among-house movement. Given these patterns, we conclude that house structural improvements, an integral part of the Ecohealth approach that makes houses refractory to vector colonization and persistence, are critical for long-term reduction of T. dimidiata infestation.

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.

Correction: A single Danio rerio hars gene encodes both cytoplasmic and mitochondrial histidyl-tRNA synthetases.

[This corrects the article DOI: 10.1371/journal.pone.0185317.].

Vectors of diversity: Genome wide diversity across the geographic range of the Chagas disease vector Triatoma dimidiata sensu lato (Hemiptera: Reduviidae).

To date, the phylogeny of Triatoma dimidiata sensu lato (s. l.) (Hemiptera: Reduviidae: Triatominae), the epidemiologically most important Chagas disease vector in Central America and a secondary vector in Mexico and northern South America, has only been investigated by one multi-copy nuclear gene (Internal Transcribed Spacer - 2) and a few mitochondrial genes. We examined 450 specimens sampled across most of its native range from Mexico to Ecuador using reduced representation next-generation sequencing encompassing over 16,000 single nucleotide polymorphisms (SNPs). Using a combined phylogenetic and species delimitation approach we uncovered two distinct species, as well as a well-defined third group that may contain multiple species. The findings are discussed with respect to possible drivers of diversification and the epidemiological importance of the distinct species and groups.

A single Danio rerio hars gene encodes both cytoplasmic and mitochondrial histidyl-tRNA synthetases.

Histidyl tRNA Synthetase (HARS) is a member of the aminoacyl tRNA synthetase (ARS) family of enzymes. This family of 20 enzymes is responsible for attaching specific amino acids to their cognate tRNA molecules, a critical step in protein synthesis. However, recent work highlighting a growing number of associations between ARS genes and diverse human diseases raises the possibility of new and unexpected functions in this ancient enzyme family. For example, mutations in HARS have been linked to two different neurological disorders, Usher Syndrome Type IIIB and Charcot Marie Tooth peripheral neuropathy. These connections raise the possibility of previously undiscovered roles for HARS in metazoan development, with alterations in these functions leading to complex diseases. In an attempt to establish Danio rerio as a model for studying HARS functions in human disease, we characterized the Danio rerio hars gene and compared it to that of human HARS. Using a combination of bioinformatics, molecular biology, and cellular approaches, we found that while the human genome encodes separate genes for cytoplasmic and mitochondrial HARS protein, the Danio rerio genome encodes a single hars gene which undergoes alternative splicing to produce the respective cytoplasmic and mitochondrial versions of Hars. Nevertheless, while the HARS genes of humans and Danio differ significantly at the genomic level, we found that they are still highly conserved at the amino acid level, underscoring the potential utility of Danio rerio as a model organism for investigating HARS function and its link to human diseases in vivo.

Heat tolerance predicts the importance of species interaction effects as the climate changes.

SYNOPSIS: Few studies have quantified the relative importance of direct effects of climate change on communities versus indirect effects that are mediated thorough species interactions, and the limited evidence is conflicting. Trait-based approaches have been popular in studies of climate change, but can they be used to estimate direct versus indirect effects? At the species level, thermal tolerance is a trait that is often used to predict winners and losers under scenarios of climate change. But thermal tolerance might also inform when species interactions are likely to be important because only subsets of species will be able to exploit the available warmer climatic niche space, and competition may intensify in the remaining, compressed cooler climatic niche space. Here, we explore the relative roles of the direct effects of temperature change and indirect effects of species interactions on forest ant communities that were heated as part of a large-scale climate manipulation at high- and low-latitude sites in eastern North America. Overall, we found mixed support for the importance of negative species interactions (competition), but found that the magnitude of these interaction effects was predictable based on the heat tolerance of the focal species. Forager abundance and nest site occupancy of heat-intolerant species were more often influenced by negative interactions with other species than by direct effects of temperature. Our findings suggest that measures of species-specific heat tolerance may roughly predict when species interactions will influence responses to global climate change.

Effects of desiccation and starvation on thermal tolerance and the heat-shock response in forest ants.

Temperature increases associated with global climate change are likely to be accompanied by additional environmental stressors such as desiccation and food limitation, which may alter how temperature impacts organismal performance. To investigate how interactions between stressors influence thermal tolerance in the common forest ant, Aphaenogaster picea, we compared the thermal resistance of workers to heat shock with and without pre-exposure to desiccation or starvation stress. Knockdown (KD) time at 40.5 °C of desiccated ants was reduced 6% compared to controls, although longer exposure to desiccation did not further reduce thermal tolerance. Starvation, in contrast, had an increasingly severe effect on thermal tolerance: at 21 days, average KD time of starved ants was reduced by 65% compared to controls. To test whether reduction in thermal tolerance results from impairment of the heat-shock response, we measured basal gene expression and transcriptional induction of two heat-shock proteins (hsp70 and hsp40) in treated and control ants. We found no evidence that either stressor impaired the Hsp response: both desiccation and starvation slightly increased basal Hsp expression under severe stress conditions and did not affect the magnitude of induction under heat shock. These results suggest that the co-occurrence of multiple environmental stressors predicted by climate change models may make populations more vulnerable to future warming than is suggested by the results of single-factor heating experiments.

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.

Climatic warming destabilizes forest ant communities.

How will ecological communities change in response to climate warming? Direct effects of temperature and indirect cascading effects of species interactions are already altering the structure of local communities, but the dynamics of community change are still poorly understood. We explore the cumulative effects of warming on the dynamics and turnover of forest ant communities that were warmed as part of a 5-year climate manipulation experiment at two sites in eastern North America. At the community level, warming consistently increased occupancy of nests and decreased extinction and nest abandonment. This consistency was largely driven by strong responses of a subset of thermophilic species at each site. As colonies of thermophilic species persisted in nests for longer periods of time under warmer temperatures, turnover was diminished, and species interactions were likely altered. We found that dynamical (Lyapunov) community stability decreased with warming both within and between sites. These results refute null expectations of simple temperature-driven increases in the activity and movement of thermophilic ectotherms. The reduction in stability under warming contrasts with the findings of previous studies that suggest resilience of species interactions to experimental and natural warming. In the face of warmer, no-analog climates, communities of the future may become increasingly fragile and unstable.

Thermal reactionomes reveal divergent responses to thermal extremes in warm and cool-climate ant species.

BACKGROUND: The distributions of species and their responses to climate change are in part determined by their thermal tolerances. However, little is known about how thermal tolerance evolves. To test whether evolutionary extension of thermal limits is accomplished through enhanced cellular stress response (enhanced response), constitutively elevated expression of protective genes (genetic assimilation) or a shift from damage resistance to passive mechanisms of thermal stability (tolerance), we conducted an analysis of the reactionome: the reaction norm for all genes in an organism's transcriptome measured across an experimental gradient. We characterized thermal reactionomes of two common ant species in the eastern U.S, the northern cool-climate Aphaenogaster picea and the southern warm-climate Aphaenogaster carolinensis, across 12 temperatures that spanned their entire thermal breadth. RESULTS: We found that at least 2 % of all genes changed expression with temperature. The majority of upregulation was specific to exposure to low temperatures. The cool-adapted A. picea induced expression of more genes in response to extreme temperatures than did A. carolinensis, consistent with the enhanced response hypothesis. In contrast, under high temperatures the warm-adapted A. carolinensis downregulated many of the genes upregulated in A. picea, and required more extreme temperatures to induce down-regulation in gene expression, consistent with the tolerance hypothesis. We found no evidence for a trade-off between constitutive and inducible gene expression as predicted by the genetic assimilation hypothesis. CONCLUSIONS: These results suggest that increases in upper thermal limits may require an evolutionary shift in response mechanism away from damage repair toward tolerance and prevention.

The evolution of heat shock protein sequences, cis-regulatory elements, and expression profiles in the eusocial Hymenoptera.

BACKGROUND: The eusocial Hymenoptera have radiated across a wide range of thermal environments, exposing them to significant physiological stressors. We reconstructed the evolutionary history of three families of Heat Shock Proteins (Hsp90, Hsp70, Hsp40), the primary molecular chaperones protecting against thermal damage, across 12 Hymenopteran species and four other insect orders. We also predicted and tested for thermal inducibility of eight Hsps from the presence of cis-regulatory heat shock elements (HSEs). We tested whether Hsp induction patterns in ants were associated with different thermal environments. RESULTS: We found evidence for duplications, losses, and cis-regulatory changes in two of the three gene families. One member of the Hsp90 gene family, hsp83, duplicated basally in the Hymenoptera, with shifts in HSE motifs in the novel copy. Both copies were retained in bees, but ants retained only the novel HSE copy. For Hsp70, Hymenoptera lack the primary heat-inducible orthologue from Drosophila melanogaster and instead induce the cognate form, hsc70-4, which also underwent an early duplication. Episodic diversifying selection was detected along the branch predating the duplication of hsc70-4 and continued along one of the paralogue branches after duplication. Four out of eight Hsp genes were heat-inducible and matched the predictions based on presence of conserved HSEs. For the inducible homologues, the more thermally tolerant species, Pogonomyrmex barbatus, had greater Hsp basal expression and induction in response to heat stress than did the less thermally tolerant species, Aphaenogaster picea. Furthermore, there was no trade-off between basal expression and induction. CONCLUSIONS: Our results highlight the unique evolutionary history of Hsps in eusocial Hymenoptera, which has been shaped by gains, losses, and changes in cis-regulation. Ants, and most likely other Hymenoptera, utilize lineage-specific heat inducible Hsps, whose expression patterns are associated with adaptive variation in thermal tolerance between two ant species. Collectively, our analyses suggest that Hsp sequence and expression patterns may reflect the forces of selection acting on thermal tolerance in ants and other social Hymenoptera.

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.