Parthenogenetic Stick Insects Exhibit Signatures of Preservation in the Molecular Architecture of Male Reproduction.

After the loss of a trait, theory predicts that the molecular machinery underlying its phenotypic expression should decay. Yet, empirical evidence is contrasting. Here, we test the hypotheses that (i) the molecular ground plan of a lost trait could persist due to pleiotropic effects on other traits and (ii) that gene co-expression network architecture could constrain individual gene expression. Our testing ground has been the Bacillus stick insect species complex, which contains close relatives that are either bisexual or parthenogenetic. After the identification of genes expressed in male reproductive tissues in a bisexual species, we investigated their gene co-expression network structure in two parthenogenetic species. We found that gene co-expression within the male gonads was partially preserved in parthenogens. Furthermore, parthenogens did not show relaxed selection on genes upregulated in male gonads in the bisexual species. As these genes were mostly expressed in female gonads, this preservation could be driven by pleiotropic interactions and an ongoing role in female reproduction. Connectivity within the network also played a key role, with highly connected-and more pleiotropic-genes within male gonad also having a gonad-biased expression in parthenogens. Our findings provide novel insight into the mechanisms which could underlie the production of rare males in parthenogenetic lineages; more generally, they provide an example of the cryptic persistence of a lost trait molecular architecture, driven by gene pleiotropy on other traits and within their co-expression network.

Sodium Channel ß Subunits-An Additional Element in Animal Tetrodotoxin Resistance?

Tetrodotoxin (TTX) is a neurotoxic molecule used by many animals for defense and/or predation, as well as an important biomedical tool. Its ubiquity as a defensive agent has led to repeated independent evolution of tetrodotoxin resistance in animals. TTX binds to voltage-gated sodium channels (VGSC) consisting of α and ß subunits. Virtually all studies investigating the mechanisms behind TTX resistance have focused on the α subunit of voltage-gated sodium channels, where tetrodotoxin binds. However, the possibility of ß subunits also contributing to tetrodotoxin resistance was never explored, though these subunits act in concert. In this study, we present preliminary evidence suggesting a potential role of ß subunits in the evolution of TTX resistance. We gathered mRNA sequences for all ß subunit types found in vertebrates across 12 species (three TTX-resistant and nine TTX-sensitive) and tested for signatures of positive selection with a maximum likelihood approach. Our results revealed several sites experiencing positive selection in TTX-resistant taxa, though none were exclusive to those species in subunit ß1, which forms a complex with the main physiological target of TTX (VGSC Nav1.4). While experimental data validating these findings would be necessary, this work suggests that deeper investigation into ß subunits as potential players in tetrodotoxin resistance may be worthwhile.

Evolutionarily diverse origins of deformed wing viruses in western honey bees.

Novel transmission routes can allow infectious diseases to spread, often with devastating consequences. Ectoparasitic varroa mites vector a diversity of RNA viruses, having switched hosts from the eastern to western honey bees (Apis cerana to Apis mellifera). They provide an opportunity to explore how novel transmission routes shape disease epidemiology. As the principal driver of the spread of deformed wing viruses (mainly DWV-A and DWV-B), varroa infestation has also driven global honey bee health declines. The more virulent DWV-B strain has been replacing the original DWV-A strain in many regions over the past two decades. Yet, how these viruses originated and spread remains poorly understood. Here, we use a phylogeographic analysis based on whole-genome data to reconstruct the origins and demography of DWV spread. We found that, rather than reemerging in western honey bees after varroa switched hosts, as suggested by previous work, DWV-A most likely originated in East Asia and spread in the mid-20th century. It also showed a massive population size expansion following the varroa host switch. By contrast, DWV-B was most likely acquired more recently from a source outside East Asia and appears absent from the original varroa host. These results highlight the dynamic nature of viral adaptation, whereby a vector's host switch can give rise to competing and increasingly virulent disease pandemics. The evolutionary novelty and rapid global spread of these host-virus interactions, together with observed spillover into other species, illustrate how increasing globalization poses urgent threats to biodiversity and food security.

Genome of the Lord Howe Island Stick Insect Reveals a Highly Conserved Phasmid X Chromosome.

We present a chromosome-scale genome assembly for Dryococelus australis, a critically endangered Australian phasmid. The assembly, constructed with Pacific Biosciences continuous long reads and chromatin conformation capture (Omni-C) data, is 3.42 Gb in length with a scaffold N50 of 262.27 Mb and L50 of 5. Over 99% of the assembly is contained in 17 major scaffolds, which corresponds to the species' karyotype. The assembly contains 96.3% of insect Benchmarking Unique Single Copy Ortholog genes in single copy. A custom repeat library identified 63.29% of the genome covered by repetitive elements; most were not identifiable based on similarity to sequences in existing databases. A total of 33,793 putative protein-coding genes were annotated. Despite the high contiguity and single-copy Benchmarking Unique Single Copy Ortholog content of the assembly, over 1 Gb of the flow-cytometry-estimated genome size is not represented, likely due to the large and repetitive nature of the genome. We identified the X chromosome with a coverage-based analysis and searched for homologs of genes known to be X-linked across the genus Timema. We found 59% of these genes on the putative X chromosome, indicating strong conservation of X-chromosomal content across 120 million years of phasmid evolution.

The final frontier: ecological and evolutionary dynamics of a global parasite invasion.

Studying rapid biological changes accompanying the introduction of alien organisms into native ecosystems can provide insights into fundamental ecological and evolutionary theory. While powerful, this quasi-experimental approach is difficult to implement because the timing of invasions and their consequences are hard to predict, meaning that baseline pre-invasion data are often missing. Exceptionally, the eventual arrival of Varroa destructor (hereafter Varroa) in Australia has been predicted for decades. Varroa is a major driver of honeybee declines worldwide, particularly as vectors of diverse RNA viruses. The detection of Varroa in 2022 at over a hundred sites poses a risk of further spread across the continent. At the same time, careful study of Varroa's spread, if it does become established, can provide a wealth of information that can fill knowledge gaps about its effects worldwide. This includes how Varroa affects honeybee populations and pollination. Even more generally, Varroa invasion can serve as a model for evolution, virology and ecological interactions between the parasite, the host and other organisms.

Vector-virus interaction affects viral loads and co-occurrence.

BACKGROUND: Vector-borne viral diseases threaten human and wildlife worldwide. Vectors are often viewed as a passive syringe injecting the virus. However, to survive, replicate and spread, viruses must manipulate vector biology. While most vector-borne viral research focuses on vectors transmitting a single virus, in reality, vectors often carry diverse viruses. Yet how viruses affect the vectors remains poorly understood. Here, we focused on the varroa mite (Varroa destructor), an emergent parasite that can carry over 20 honey bee viruses, and has been responsible for colony collapses worldwide, as well as changes in global viral populations. Co-evolution of the varroa and the viral community makes it possible to investigate whether viruses affect vector gene expression and whether these interactions affect viral epidemiology. RESULTS: Using a large set of available varroa transcriptomes, we identified how abundances of individual viruses affect the vector's transcriptional network. We found no evidence of competition between viruses, but rather that some virus abundances are positively correlated. Furthermore, viruses that are found together interact with the vector's gene co-expression modules in similar ways, suggesting that interactions with the vector affect viral epidemiology. We experimentally validated this observation by silencing candidate genes using RNAi and found that the reduction in varroa gene expression was accompanied by a change in viral load. CONCLUSIONS: Combined, the meta-transcriptomic analysis and experimental results shed light on the mechanism by which viruses interact with each other and with their vector to shape the disease course.

Experimental inheritance of antibiotic acquired dysbiosis affects host phenotypes across generations.

Microbiomes can enhance the health, fitness and even evolutionary potential of their hosts. Many organisms propagate favorable microbiomes fully or partially via vertical transmission. In the long term, such co-propagation can lead to the evolution of specialized microbiomes and functional interdependencies with the host. However, microbiomes are vulnerable to environmental stressors, particularly anthropogenic disturbance such as antibiotics, resulting in dysbiosis. In cases where microbiome transmission occurs, a disrupted microbiome may then become a contagious pathology causing harm to the host across generations. We tested this hypothesis using the specialized socially transmitted gut microbiome of honey bees as a model system. By experimentally passaging tetracycline-treated microbiomes across worker 'generations' we found that an environmentally acquired dysbiotic phenotype is heritable. As expected, the antibiotic treatment disrupted the microbiome, eliminating several common and functionally important taxa and strains. When transmitted, the dysbiotic microbiome harmed the host in subsequent generations. Particularly, naïve bees receiving antibiotic-altered microbiomes died at higher rates when challenged with further antibiotic stress. Bees with inherited dysbiotic microbiomes showed alterations in gene expression linked to metabolism and immunity, among other pathways, suggesting effects on host physiology. These results indicate that there is a possibility that sublethal exposure to chemical stressors, such as antibiotics, may cause long-lasting changes to functional host-microbiome relationships, possibly weakening the host's progeny in the face of future ecological challenges. Future studies under natural conditions would be important to examine the extent to which negative microbiome-mediated phenotypes could indeed be heritable and what role this may play in the ongoing loss of biodiversity.

Gene transcriptional profiles in gonads of Bacillus taxa (Phasmida) with different cytological mechanisms of automictic parthenogenesis.

The evolution of automixis - i.e., meiotic parthenogenesis - requires several features, including ploidy restoration after meiosis and maintenance of fertility. Characterizing the relative contribution of novel versus pre-existing genes and the similarities in their expression and sequence evolution is fundamental to understand the evolution of reproductive novelties. Here we identify gonads-biased genes in two Bacillus automictic stick-insects and compare their expression profile and sequence evolution with a bisexual congeneric species. The two parthenogens restore ploidy through different cytological mechanisms: in Bacillus atticus, nuclei derived from the first meiotic division fuse to restore a diploid egg nucleus, while in Bacillus rossius, diploidization occurs in some cells of the haploid blastula through anaphase restitution. Parthenogens' gonads transcriptional program is found to be largely assembled from genes that were already present before the establishment of automixis. The three species transcriptional profiles largely reflect their phyletic relationships, yet we identify a shared core of genes with gonad-biased patterns of expression in parthenogens which are either male gonads-biased in the sexual species or are not differentially expressed there. At the sequence level, just a handful of gonads-biased genes were inferred to have undergone instances of positive selection exclusively in the parthenogen species. This work is the first to explore the molecular underpinnings of automixis in a comparative framework: it delineates how reproductive novelties can be sustained by genes whose origin precedes the establishment of the novelty itself and shows that different meiotic mechanisms of reproduction can be associated with a shared molecular ground plan.

Genomic signatures of recent convergent transitions to social life in spiders.

The transition from solitary to social life is a major phenotypic innovation, but its genetic underpinnings are largely unknown. To identify genomic changes associated with this transition, we compare the genomes of 22 spider species representing eight recent and independent origins of sociality. Hundreds of genes tend to experience shifts in selection during the repeated transition to social life. These genes are associated with several key functions, such as neurogenesis, behavior, and metabolism, and include genes that previously have been implicated in animal social behavior and human behavioral disorders. In addition, social species have elevated genome-wide rates of molecular evolution associated with relaxed selection caused by reduced effective population size. Altogether, our study provides unprecedented insights into the genomic signatures of social evolution and the specific genetic changes that repeatedly underpin the evolution of sociality. Our study also highlights the heretofore unappreciated potential of transcriptomics using ethanol-preserved specimens for comparative genomics and phylotranscriptomics.

A Genome for Edith's Checkerspot Butterfly: An Insect with Complex Host-Adaptive Suites and Rapid Evolutionary Responses to Environmental Changes.

Insects have been key players in the assessments of biodiversity impacts of anthropogenically driven environmental change, including the evolutionary and ecological impacts of climate change. Populations of Edith's Checkerspot Butterfly (Euphydryas editha) adapt rapidly to diverse environmental conditions, with numerous high-impact studies documenting these dynamics over several decades. However, studies of the underlying genetic bases of these responses have been hampered by missing genomic resources, limiting the ability to connect genomic responses to environmental change. Using a combination of Oxford Nanopore long reads, haplotype merging, HiC scaffolding followed by Illumina polishing, we generated a highly contiguous and complete assembly (contigs n = 142, N50 = 21.2 Mb, total length = 607.8 Mb; BUSCOs n = 5,286, single copy complete = 97.8%, duplicated = 0.9%, fragmented = 0.3%, missing = 1.0%). A total of 98% of the assembled genome was placed into 31 chromosomes, which displayed large-scale synteny with other well-characterized lepidopteran genomes. The E. editha genome, annotation, and functional descriptions now fill a missing gap for one of the leading field-based ecological model systems in North America.

The first steps toward a global pandemic: Reconstructing the demographic history of parasite host switches in its native range.

Host switching allows parasites to expand their niches. However, successful switching may require suites of adaptations and also may decrease performance on the old host. As a result, reductions in gene flow accompany many host switches, driving speciation. Because host switches tend to be rapid, it is difficult to study them in real-time, and their demographic parameters remain poorly understood. As a result, fundamental factors that control subsequent parasite evolution, such as the size of the switching population or the extent of immigration from the original host, remain largely unknown. To shed light on the host switching process, we explored how host switches occur in independent host shifts by two ectoparasitic honey bee mites (Varroa destructor and V. jacobsoni). Both switched to the western honey bee (Apis mellifera) after being brought into contact with their ancestral host (Apis cerana), ~70 and ~12 years ago, respectively. Varroa destructor subsequently caused worldwide collapses of honey bee populations. Using whole-genome sequencing on 63 mites collected in their native ranges from both the ancestral and novel hosts, we were able to reconstruct the known temporal dynamics of the switch. We further found multiple previously undiscovered mitochondrial lineages on the novel host, along with the genetic equivalent of tens of individuals that were involved in the initial host switch. Despite being greatly reduced, some gene flow remains between mites adapted to different hosts. Our findings suggest that while reproductive isolation may facilitate the fixation of traits beneficial for exploiting the new host, ongoing genetic exchange may allow genetic amelioration of inbreeding effects.

Honey Bee Larval and Adult Microbiome Life Stages Are Effectively Decoupled with Vertical Transmission Overcoming Early Life Perturbations.

Microbiomes provide a range of benefits to their hosts which can lead to the coevolution of a joint ecological niche. However, holometabolous insects, some of the most successful organisms on Earth, occupy different niches throughout development, with larvae and adults being physiologically and morphologically highly distinct. Furthermore, transition between the stages usually involves the loss of the gut microbiome since the gut is remodeled during pupation. Most eusocial organisms appear to have evolved a workaround to this problem by sharing their communal microbiome across generations. However, whether this vertical microbiome transmission can overcome perturbations of the larval microbiome remains untested. Honey bees have a relatively simple, conserved, coevolved adult microbiome which is socially transmitted and affects many aspects of their biology. In contrast, larval microbiomes are more variable, with less clear roles. Here, we manipulated the gut microbiome of in vitro-reared larvae, and after pupation of the larvae, we inoculated the emerged bees with adult microbiome to test whether adult and larval microbiome stages may be coupled (e.g., through immune priming). Larval treatments differed in bacterial composition and abundance, depending on diet, which also drove larval gene expression. Nonetheless, adults converged on the typical core taxa and showed limited gene expression variation. This work demonstrates that honey bee adult and larval stages are effectively microbiologically decoupled, and the core adult microbiome is remarkably stable to early developmental perturbations. Combined with the transmission of the microbiome in early adulthood, this allows the formation of long-term host-microbiome associations. IMPORTANCE This work investigated host-microbiome interactions during a crucial developmental stage-the transition from larvae to adults, which is a challenge to both, the insect host and its microbiome. Using the honey bee as a tractable model system, we showed that microbiome transfer after emergence overrides any variation in the larvae, indicating that larval and adult microbiome stages are effectively decoupled. Together with the reliable vertical transfer in the eusocial system, this decoupling ensures that the adults are colonized with a consistent and derived microbiome after eclosion. Taken all together, our data provide additional support that the evolution of sociality, at least in the honey bee system tested here, is linked with host-microbiome relationships.

Co-option of the same ancestral gene family gave rise to mammalian and reptilian toxins.

BACKGROUND: Evolution can occur with surprising predictability when organisms face similar ecological challenges. For most traits, it is difficult to ascertain whether this occurs due to constraints imposed by the number of possible phenotypic solutions or because of parallel responses by shared genetic and regulatory architecture. Exceptionally, oral venoms are a tractable model of trait evolution, being largely composed of proteinaceous toxins that have evolved in many tetrapods, ranging from reptiles to mammals. Given the diversity of venomous lineages, they are believed to have evolved convergently, even though biochemically similar toxins occur in all taxa. RESULTS: Here, we investigate whether ancestral genes harbouring similar biochemical activity may have primed venom evolution, focusing on the origins of kallikrein-like serine proteases that form the core of most vertebrate oral venoms. Using syntenic relationships between genes flanking known toxins, we traced the origin of kallikreins to a single locus containing one or more nearby paralogous kallikrein-like clusters. Additionally, phylogenetic analysis of vertebrate serine proteases revealed that kallikrein-like toxins in mammals and reptiles are genetically distinct from non-toxin ones. CONCLUSIONS: Given the shared regulatory and genetic machinery, these findings suggest that tetrapod venoms evolved by co-option of proteins that were likely already present in saliva. We term such genes 'toxipotent'-in the case of salivary kallikreins they already had potent vasodilatory activity that was weaponized by venomous lineages. Furthermore, the ubiquitous distribution of kallikreins across vertebrates suggests that the evolution of envenomation may be more common than previously recognized, blurring the line between venomous and non-venomous animals.

Functional innovation promotes diversification of form in the evolution of an ultrafast trap-jaw mechanism in ants.

Evolutionary innovations underlie the rise of diversity and complexity-the 2 long-term trends in the history of life. How does natural selection redesign multiple interacting parts to achieve a new emergent function? We investigated the evolution of a biomechanical innovation, the latch-spring mechanism of trap-jaw ants, to address 2 outstanding evolutionary problems: how form and function change in a system during the evolution of new complex traits, and whether such innovations and the diversity they beget are repeatable in time and space. Using a new phylogenetic reconstruction of 470 species, and X-ray microtomography and high-speed videography of representative taxa, we found the trap-jaw mechanism evolved independently 7 to 10 times in a single ant genus (Strumigenys), resulting in the repeated evolution of diverse forms on different continents. The trap mechanism facilitates a 6 to 7 order of magnitude greater mandible acceleration relative to simpler ancestors, currently the fastest recorded acceleration of a resettable animal movement. We found that most morphological diversification occurred after evolution of latch-spring mechanisms, which evolved via minor realignments of mouthpart structures. This finding, whereby incremental changes in form lead to a change of function, followed by large morphological reorganization around the new function, provides a model for understanding the evolution of complex biomechanical traits, as well as insights into why such innovations often happen repeatedly.

An ancient, conserved gene regulatory network led to the rise of oral venom systems.

Oral venom systems evolved multiple times in numerous vertebrates enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on the toxins. However, using toxins present in modern day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the "metavenom network," a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with the toxins, and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon, evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.

Markerless tracking of an entire honey bee colony.

From cells in tissue, to bird flocks, to human crowds, living systems display a stunning variety of collective behaviors. Yet quantifying such phenomena first requires tracking a significant fraction of the group members in natural conditions, a substantial and ongoing challenge. We present a comprehensive, computational method for tracking an entire colony of the honey bee Apis mellifera using high-resolution video on a natural honeycomb background. We adapt a convolutional neural network (CNN) segmentation architecture to automatically identify bee and brood cell positions, body orientations and within-cell states. We achieve high accuracy (~10% body width error in position, ~10° error in orientation, and true positive rate > 90%) and demonstrate months-long monitoring of sociometric colony fluctuations. These fluctuations include ~24 h cycles in the counted detections, negative correlation between bee and brood, and nightly enhancement of bees inside comb cells. We combine detected positions with visual features of organism-centered images to track individuals over time and through challenging occluding events, recovering ~79% of bee trajectories from five observation hives over 5 min timespans. The trajectories reveal important individual behaviors, including waggle dances and crawling inside comb cells. Our results provide opportunities for the quantitative study of collective bee behavior and for advancing tracking techniques of crowded systems.

Intraspecific niche partition without speciation: individual level web polymorphism within a single island spider population.

Early in the process of adaptive radiation, allopatric disruption of gene flow followed by ecological specialization is key for speciation; but, do adaptive radiations occur on small islands without internal geographical barriers? Island populations sometimes harbour polymorphism in ecological specializations, but its significance remains unclear. On one hand, morphs may correspond to 'cryptic' species. Alternatively, they could result from population, developmental or behavioural plasticity. The spider Wendilgarda galapagensis (Araneae, Theridiosomatidae) is endemic to the small Isla del Coco and unique in spinning three different web types, each corresponding to a different microhabitat. We tested whether this variation is associated with 'cryptic' species or intraspecific behavioural plasticity. Despite analysing 36 803 loci across 142 individuals, we found no relationship between web type and population structure, which was only weakly geographically differentiated. The same pattern holds when looking within a sampling site or considering only Fst outliers. In line with genetic data, translocation experiments showed that web architecture is plastic within an individual. However, not all transitions between web types are equally probable, indicating the existence of individual preferences. Our data supports the idea that diversification on small islands might occur mainly at the behavioural level producing an intraspecific niche partition without speciation.

Transcription Profiling of Cultured Acropora digitifera Adult Cells Reveals the Existence of Ancestral Genome Regulatory Modules Underlying Pluripotency and Cell Differentiation in Cnidaria.

Due to their pluripotent nature and unlimited cell renewal, stem cells have been proposed as an ideal material for establishing long-term cnidarian cell cultures. However, the lack of unifying principles associated with "stemness" across the phylum complicates stem cells' identification and isolation. Here, we for the first time report gene expression profiles for cultured coral cells, focusing on regulatory gene networks underlying pluripotency and differentiation. Cultures were initiated from Acropora digitifera tip fragments, the fastest growing tissue in Acropora. Overall, in vitro transcription resembled early larvae, overexpressing orthologs of premetazoan and Hydra stem cell markers, and transcripts with roles in cell division, migration, and differentiation. Our results suggest the presence of pluripotent cell types in cultures and indicate the existence of ancestral genome regulatory modules underlying pluripotency and cell differentiation in cnidaria. Cultured cells appear to be synthesizing protein, differentiating, and proliferating.