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.

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.

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.].

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.

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.

Inter-genomic sexual conflict drives antagonistic coevolution in harvester ants.

The reproductive interests of males and females are not always aligned, leading to sexual conflict over parental investment, rate of reproduction and mate choice. Traits that increase the genetic interests of one sex often occur at the expense of the other, selecting for counter-adaptations leading to antagonistic coevolution. Reproductive conflict is not limited to intraspecific interactions; interspecific hybridization can produce pronounced sexual conflict between males and females of different species, but it is unclear whether such conflict can drive sexually antagonistic coevolution between reproductively isolated genomes. We tested for hybridization-driven sexually antagonistic adaptations in queens and males of the socially hybridogenetic 'J' lineages of Pogonomyrmex harvester ants, whose mating system promotes hybridization in queens but selects against it in males. We conducted no-choice mating assays to compare patterns of mating behaviour and sperm transfer between inter- and intra-lineage pairings. There was no evidence for mate discrimination on the basis of pair type, and the total quantity of sperm transferred did not differ between intra- and inter-lineage pairs; however, further dissection of the sperm transfer process into distinct mechanistic components revealed significant, and opposing, cryptic manipulation of copulatory investment by both sexes. Males of both lineages increased their rate of sperm transfer to high-fitness intra-lineage mates, with a stronger response in the rarer lineage for whom mating mistakes are the most likely. By contrast, the total duration of copulation for intra-lineage mating pairs was significantly shorter than for inter-lineage crosses, suggesting that queens respond to prevent excessive sperm loading by prematurely terminating copulation. These findings demonstrate that sexual conflict can lead to antagonistic coevolution in both intra-genomic and inter-genomic contexts. Indeed, the resolution of sexual conflict may be a key determinant of the long-term evolutionary potential of host-dependent reproductive strategies, counteracting the inherent instabilities arising from such systems.

A novel family of terminal-repeat retrotransposon in miniature (TRIM) in the genome of the red harvester ant, Pogonomyrmex barbatus.

We report the first described non-plant family of TRIMs (terminal-repeat retrotransposons in miniature), which are small nonautonomous LTR retrotransposons, from the whole-genome sequence of the red harvester ant, Pogonomyrmex barbatus (Hymenoptera: Myrmicinae). Members of this retrotransposon family, named PbTRIM, have typical features of plant TRIMs in length and structure, although they share no overall sequence similarity. PbTRIM elements and their solo-LTRs are abundant in the host genome and exhibit an uneven distribution pattern. Elements are preferentially inserted into TA-rich regions with ATAT as the most common pattern of target site duplication (TSD). PbTRIM is most likely mobile as indicated by the young age of many complete elements, the high degree of sequence similarity among elements at different genomic locations, the abundance of elements in the host genome, and the presence of 4-bp target site duplications (TSDs) flanking the elements and solo-LTRs. Many PbTRIM elements and their solo-LTRs are located within or near genes, suggesting their potential roles in restructuring the host genes and genome. Database search, PCR and sequencing analysis revealed the presence of homologous PbTRIM elements in other ant species. The high sequence similarity between elements from distantly related ant species, the incongruence between the phylogenies of PbTRIM and its hosts, and the patchy distribution of the retroelement within the Myrmicinae subfamily indicate possible horizontal transfer events of the retroelement.

Intergenerational effect of juvenile hormone on offspring in Pogonomyrmex harvester ants.

Parents can influence the phenotypes of their offspring via a number of mechanisms. In harvester ants, whether female progeny develop into workers or daughter queens is strongly influenced by the age and temperature conditions experienced by their mother, which is associated with variation in maternal ecdysteroid deposition in fertilized eggs. In many insects, juvenile hormone (JH) is antagonistic to ecdysteroid release, suggesting that seasonal and age-based variation in maternal JH titers may explain maternal effects on offspring size and reproductive caste. To test this hypothesis, we artificially increased maternal JH titers with methoprene, a JH analog, in laboratory colonies of two Pogonomyrmex populations exhibiting genetic caste determination. Increasing maternal JH resulted in a 50% increase in worker body size, as well as a sharp reduction in total number of progeny reared, but did not alter the genotype of progeny reared to adulthood. The intergenerational effect of JH manipulation was not mediated by a reduction in ecdysteroid deposition into eggs; instead, changes in egg size, trophic egg availability or brood/worker ratio may have altered the nutritional environment of developing larvae. Egg ecdysteroid content was significantly negatively correlated with natural variation in worker body size, however, suggesting that there are multiple independent routes by which queens can modify offspring phenotypes.

Mechanisms of reproductive isolation between an ant species of hybrid origin and one of its parents.

The establishment of new species by hybridization is difficult because it requires the development of reproductive isolation (RI) in sympatry to escape the homogenizing effects of gene flow from the parental species. Here we investigated the role of two pre- and two postzygotic mechanisms of RI in a system comprising two interdependent Pogonomyrmex harvester ant lineages (the H1 and H2 lineages) of hybrid origin and one of their parental species (P. rugosus). Similar to most other ants, P. rugosus is characterized by an environmental system of caste determination with female brood developing either into queens or workers depending on nongenetic factors. By contrast, there is a strong genetic component to caste determination in the H1 and H2 lineages because the developmental fate of female brood depends on the genetic origin of the parents, with interlineage eggs developing into workers and intralineage eggs developing into queens. The study of a mixed mating aggregation revealed strong differences in mating flight timing between P. rugosus and the two lineages as a first mechanism of RI. A second important prezygotic mechanism was assortative mating. Laboratory experiments also provided support for one of the two investigated mechanisms of postzygotic isolation. The majority of offspring produced from the few matings between P. rugosus and the lineages aborted at the egg stage. This hybrid inviability was under maternal influence, with hybrids produced by P. rugosus queens being always inviable whereas a small proportion of H2 lineage queens produced large numbers of adult hybrid offspring. Finally, we found no evidence that genetic caste determination acted as a second postzygotic mechanism reducing gene flow between P. rugosus and the H lineages. The few viable P. rugosus-H hybrids were not preferentially shunted into functionally sterile workers but developed into both workers and queens. Overall, these results reveal that the nearly complete (99.5%) RI between P. rugosus and the two hybrid lineages stems from the combination of two typical prezygotic mechanisms (mating time divergence and assortative mating) and one postzygotic mechanism (hybrid inviability).

Maternal effect on female caste determination in a social insect.

Caste differentiation and division of labor are the hallmarks of social insect colonies [1, 2]. The current dogma for female caste differentiation is that female eggs are totipotent, with morphological and physiological differences between queens and workers stemming from a developmental switch during the larval stage controlled by nutritional and other environmental factors (e.g., [3-8]). In this study, we tested whether maternal effects influence caste differentiation in Pogonomyrmex harvester ants. By conducting crossfostering experiments we identified two key factors in the process of caste determination. New queens were produced only from eggs laid by queens exposed to cold. Moreover, there was a strong age effect, with development into queens occurring only in eggs laid by queens that were at least two years old. Biochemical analyses further revealed that the level of ecdysteroids was significantly lower in eggs developing into queens than workers. By contrast, we found no significant effect of colony size or worker exposure to cold, suggesting that the trigger for caste differentiation may be independent of the quantity and quality of resources provided to larvae. Altogether these data demonstrate that the developmental fate of female brood is strongly influenced by maternal effects in ants of the genus Pogonomyrmex.

Two alternate mechanisms contribute to the persistence of interdependent lineages in Pogonomyrmex harvester ants.

Some populations of Pogonomyrmex harvester ants comprise pairs of highly differentiated lineages with queens mating at random with several males of their own and of the alternate lineage. These queens produce two types of diploid offspring, those fertilized by males of the queens' lineage which develop into new queens and those fertilized by males of the other lineage which mostly develop into functionally sterile workers. This unusual mode of genetic caste determination has been found in 26 populations and a total of four lineage pairs (F(1)-F(2), G(1)-G(2), H(1)-H(2) and J(1)-J(2)) have been described in these populations. Despite the fact that a few interlineage queens are produced, previous studies revealed that there is a complete lack of genetic introgression between lineages. Here we quantify the proportion of interlineage queens produced in each of the four lineage pairs and determine the fate of these queens. In the F(1)-F(2), G(1)-G(2) and H(1)-H(2) lineage pairs, interlineage queens were produced by a minority of colonies. These colonies exclusively produced interlineage queens and workers, suggesting that interlineage eggs can develop into queens in these three pairs of lineages in the absence of competition with pure-lineage brood. An analysis of three key stages of the colony life cycle revealed that colonies headed by interlineage queens failed to grow sufficiently to produce reproductive individuals. In laboratory comparisons, interlineage queens produced fewer viable eggs, with the effect that they raised fewer workers and lost more weight per worker produced than pure-lineage queens. In the J(1)-J(2) lineage pair, we did not find a single interlineage queen, raising the possibility that interlineage eggs have completely lost the ability to develop into queens in this lineage pair. Hence, two distinct mechanisms seem to account for the complete lack of between-lineage gene flow in the F(1)-F(2), G(1)-G(2), H(1)-H(2) and J(1)-J(2) lineage pairs.

Characterization and distribution of Pogonomyrmex harvester ant lineages with genetic caste determination.

Genetic caste determination has been described in two populations of Pogonomyrmex harvester ants, each comprising a pair of interbreeding lineages. Queens mate with males of their own and of the alternate lineage and produce two types of diploid offspring, those fertilized by males of the queens' lineage which develop into queens and those fertilized by males of the other lineage which develop into workers. Each of the lineages has been shown to be itself of hybrid origin between the species Pogonomyrmex barbatus and Pogonomyrmex rugosus, which both have typical, environmentally determined caste differentiation. In a large scale genetic survey across 35 sites in Arizona, New Mexico and Texas, we found that genetic caste determination associated with pairs of interbreeding lineages occurred frequently (in 26 out of the 35 sites). Overall, we identified eight lineages with genetic caste determination that always co-occurred in the same complementary lineage pairs. Three of the four lineage pairs appear to have a common origin while their relationship with the fourth remains unclear. The level of genetic differentiation among these eight lineages was significantly higher than the differentiation between P. rugosus and P. barbatus, which questions the appropriate taxonomic status of these genetic lineages. In addition to being genetically isolated from one another, all lineages with genetic caste determination were genetically distinct from P. rugosus and P. barbatus, even when colonies of interbreeding lineages co-occurred with colonies of either putative parent at the same site. Such nearly complete reproductive isolation between the lineages and the species with environmental caste determination might prevent the genetic caste determination system to be swept away by gene flow.

Reproductive isolation between Pogonomyrmex rugosus and two lineages with genetic caste determination.

Hybrid speciation occurs when combination of two interspecific genomes results in individuals that are of high fitness but reproductively incompatible with the parental species. Although hybrid speciation is a relatively common source of new species in plants, it appears to be a much rarer occurrence in animal taxa. Here we report on reproductive isolation and range overlap between the rough harvester ant Pogonomyrmex rugosus and two lineages with hybrid genotypes (H 1 and H2). Both lineages obligately interbreed and produce genetically distinct queen and worker offspring, a phenomenon referred to as genetic caste determination (GCD). Diploid offspring produced by gametes of the same lineage develop only into queens, whereas diploid offspring derived from gametes of distinct lineages develop into workers. We investigated small-scale patterns of gene flow between the parent and the two H lineages by sampling along an 80-km transect between a pure P. rugosus population and a two-lineage population. Microsatellite and mitochondrial markers both indicated virtually no gene flow between the parent species and either lineage even at sites where parental and H-lineage colonies co-occurred. The geographic ranges of the parental species and the two-lineage population were essentially parapatric, with a surprisingly narrow band of overlap and evidence of spatial structuring even at microgeographic scales within the transition zone. This suggests that ecological competition with the parent species plays a significant role in determining the evolutionary persistence and current distribution of the hybrid lineages and the genetic caste system.