Head anatomy and phylogenomics show the Carboniferous giant Arthropleura belonged to a millipede-centipede group.

The Carboniferous myriapod Arthropleura is the largest arthropod of all time, but its fossils are usually incomplete, limiting the understanding of its anatomy, ecology, and relationships. Micro-computed tomography applied to exceptionally preserved specimens from the Carboniferous Montceau-les-Mines Lagerstätte (France) reveals unprecedented details of its functional anatomy, such as the head and mouthparts. Arthropleura shares features with both millipedes and centipedes. Total-evidence phylogeny combining morphological and transcriptomic data resolves Arthropleura alone as a stem group millipede, but the inclusion of the highly incomplete Siluro-Devonian Eoarthropleura draws it deeper into the myriapod stem. Arthropleura suggests transitional morphology between clades united primarily by molecular information and underscores the value of total-evidence phylogenetics to understanding evolutionary history.

Rapid evolution of mitochondrion-related genes in haplodiploid arthropods.

BACKGROUND: Mitochondrial genes and nuclear genes cooperate closely to maintain the functions of mitochondria, especially in the oxidative phosphorylation (OXPHOS) pathway. However, mitochondrial genes among arthropod lineages have dramatic evolutionary rate differences. Haplodiploid arthropods often show fast-evolving mitochondrial genes. One hypothesis predicts that the small effective population size of haplodiploid species could enhance the effect of genetic drift leading to higher substitution rates in mitochondrial and nuclear genes. Alternatively, positive selection or compensatory changes in nuclear OXPHOS genes could lead to the fast-evolving mitochondrial genes. However, due to the limited number of arthropod genomes, the rates of evolution for nuclear genes in haplodiploid species, besides hymenopterans, are largely unknown. To test these hypotheses, we used data from 76 arthropod genomes, including 5 independently evolved haplodiploid lineages, to estimate the evolutionary rates and patterns of gene family turnover of mitochondrial and nuclear genes. RESULTS: We show that five haplodiploid lineages tested here have fast-evolving mitochondrial genes and fast-evolving nuclear genes related to mitochondrial functions, while nuclear genes not related to mitochondrion showed no significant evolutionary rate differences. Among hymenopterans, bees and ants show faster rates of molecular evolution in mitochondrial genes and mitochondrion-related nuclear genes than sawflies and wasps. With genome data, we also find gene family expansions and contractions in mitochondrion-related genes of bees and ants. CONCLUSIONS: Our results reject the small population size hypothesis in haplodiploid species. A combination of positive selection and compensatory changes could lead to the observed patterns in haplodiploid species. The elevated evolutionary rates in OXPHOS complex 2 genes of bees and ants suggest a unique evolutionary history of social hymenopterans.

The Arthropoda-specific Tramtrack group BTB protein domains use previously unknown interface to form hexamers.

BTB (bric-a-brack, Tramtrack, and broad complex) is a diverse group of protein-protein interaction domains found within metazoan proteins. Transcription factors contain a dimerizing BTB subtype with a characteristic N-terminal extension. The Tramtrack group (TTK) is a distinct type of BTB domain, which can multimerize. Single-particle cryo-EM microscopy revealed that the TTK-type BTB domains assemble into a hexameric structure consisting of three canonical BTB dimers connected through a previously uncharacterized interface. We demonstrated that the TTK-type BTB domains are found only in Arthropods and have undergone lineage-specific expansion in modern insects. The Drosophila genome encodes 24 transcription factors with TTK-type BTB domains, whereas only four have non-TTK-type BTB domains. Yeast two-hybrid analysis revealed that the TTK-type BTB domains have an unusually broad potential for heteromeric associations presumably through a dimer-dimer interaction interface. Thus, the TTK-type BTB domains are a structurally and functionally distinct group of protein domains specific to Arthropodan transcription factors.

Revisiting the four Hexapoda classes: Protura as the sister group to all other hexapods.

Insects represent the most diverse animal group, yet previous phylogenetic analyses based on morphological and molecular data have failed to agree on the evolutionary relationships of early insects and their six-legged relatives (together constituting the clade Hexapoda). In particular, the phylogenetic positions of the three early-diverging hexapod lineages-the coneheads (Protura), springtails (Collembola), and two-pronged bristletails (Diplura)-have been debated for over a century, with alternative topologies implying drastically different scenarios of the evolution of the insect body plan and hexapod terrestrialization. We addressed this issue by sampling all hexapod orders and experimenting with a broad range of across-site compositional heterogeneous models designed to tackle ancient divergences. Our analyses support Protura as the earliest-diverging hexapod lineage ("Protura-sister") and Collembola as a sister group to Diplura, a clade corresponding to the original composition of Entognatha, and characterized by the shared possession of internal muscles in the antennal flagellum. The previously recognized 'Elliplura' hypothesis is recovered only under the site-homogeneous substitution models with partial supermatrices. Our cross-validation analysis shows that the site-heterogeneous CAT-GTR model, which recovers "Protura-sister," fits significantly better than homogeneous models. Furthermore, the morphologically unusual Protura are also supported as the earliest-diverging hexapod lineage by other lines of evidence, such as mitogenomes, comparative embryology, and sperm morphology, which produced results similar to those in this study. Our backbone phylogeny of hexapods will facilitate the exploration of the underpinnings of hexapod terrestrialization and megadiversity.

Arthropod Phylotranscriptomics With a Special Focus on the Basal Phylogeny of the Myriapoda.

Arthropoda represents the most diverse animal phylum, but clarifying the phylogenetic relationships among arthropod taxa remains challenging given the numerous arthropod lineages that diverged over a short period of time. In order to resolve the most controversial aspects of deep arthropod phylogeny, focusing on the Myriapoda, we conducted phylogenetic analyses based on ten super-matrices comprised of 751 to 1,233 orthologous genes across 64 representative arthropod species, including 28 transcriptomes that were newly generated in this study. Our findings provide unambiguous support for the monophyly of the higher arthropod taxa, Chelicerata, Mandibulata, Myriapoda, Pancrustacea, and Hexapoda, while the Crustacea are paraphyletic, with the class Remipedia supported as the lineage most closely related to hexapods. Within the Hexapoda, our results largely affirm previously proposed phylogenetic relationships among deep hexapod lineages, except that the Paraneoptera (Hemiptera, Thysanoptera, and Psocodea) was recovered as a monophyletic lineage in some analyses. The results corroborated the recently proposed phylogenetic framework of the four myriapod classes, wherein Symphyla and Pauropoda, as well as Chilopoda and Diplopoda, are each proposed to be sister taxa. The findings provide important insights into understanding the phylogeny and evolution of arthropods.

Adaptative responses of Neurospora crassa by histidine kinases upon the attack of the arthropod Sinella curviseta.

Histidine kinases (HKs) are important sensor proteins in fungi and play an essential role in environmental adaptation. However, the mechanisms by which fungi sense and respond to fungivores attack via HKs are not fully understood. In this study, we utilized Neurospora crassa to investigate the involvement of HKs in responding to fungivores attack. We found that the 11 HKs in N. crassa not only affected the growth and development, but also led to fluctuations in antioxidant production. Ten mutants in the genes encoding HKs (except ∆phy1) showed increased production of reactive oxygen species (ROS), especially upon Sinella curviseta attack. The ROS burst triggered changes in conidia and perithecial beaks formation, as well as accumulation of ß-glucan, ergothioneine, ergosterol, and carotenoids. ß-glucan was increased in ∆hk9, ∆os1, ∆hcp1, ∆nik2, ∆sln1, ∆phy1 and ∆phy2 mutants compared to the wild-type strain. In parallel, ergothioneine accumulation was improved in ∆phy1 and ∆hk16 mutants and further increased upon attack, except in ∆os1 and ∆hk16 mutants. Additionally, fungivores attack stimulated ergosterol and dehydroergosterol production in ∆hk9 and ∆os1 mutants. Furthermore, deletion of these genes altered carotenoid accumulation, with wild-type strain, ∆hk9, ∆os1, ∆hcp1, ∆sln1, ∆phy2, and ∆dcc1mutants showing an increase in carotenoids upon attack. Taken together, HKs are involved in regulating the production of conidia and antioxidants. Thus, HKs may act as sensors of fungivores attack and effectively improve the adaptive capacity of fungi to environmental stimuli.

The relationship between neutral genetic diversity and performance in wild arthropod populations.

Larger effective populations (Ne) are characterized by higher genetic diversity, which is expected to predict population performance (average individual performance that influences fitness). Empirical studies of the relationship between neutral diversity and performance mostly represent species with small Ne, while there is limited data from the species-rich and ecologically important arthropods that are assumed to have large Ne but are threatened by massive declines. We performed a systematic literature search and used meta-analytical models to test the prediction of a positive association between neutral genetic diversity and performance in wild arthropods. From 14 relevant studies of 286 populations, we detected a weak (r = 0.15) but nonsignificant positive association both in the full data set (121 effect sizes) and a reduced data set accounting for dependency (14 effect sizes). Theory predicts that traits closely associated with fitness show a relatively stronger correlation with neutral diversity; this relationship was upheld for longevity and marginally for reproduction. Our analyses point to major knowledge gaps in our understanding of relationships between neutral diversity and performance. Future studies using genome-wide data sets across populations could guide more powerful designs to evaluate relationships between adaptive, deleterious and neutral diversity and performance.

Recovering plant-associated arthropod communities by eDNA metabarcoding historical herbarium specimens.

Natural history collections are a priceless resource for understanding patterns and processes of biodiversity change in the Anthropocene.1 Herbaria, which house millions of historical plant records from all over the globe, are particularly valuable to study population genetics of the plants themselves and to understand the assembly of plant-associated microbial communities.2 Here we test if herbaria can serve yet another essential purpose, namely to provide information on the historical assembly of plant-arthropod interactions. The specificity and temporal stability of these associations are poorly known.3 Considering their pivotal role in the assembly of terrestrial food webs,4 this knowledge is paramount to understanding the consequences of global change. We use environmental DNA (eDNA) metabarcoding to characterize communities of plant-associated arthropods from archived herbarium specimens of different ages and origins. The herbarium specimens yield arthropod DNA across various ecological guilds and trophic levels over multiple decades. In an experiment, we also show that the typical dry storage of plants in herbaria does not alter the recovered arthropod diversity and community composition. By analyzing a time series of leaf samples from a forest monitoring project, we then characterize changes in arthropod biodiversity over two decades, showing that archived plants can also provide the time series data that are urgently needed to understand arthropod declines.5 This use of herbaria and plant archives promises unprecedented insights into plant-arthropod interactions and revolutionizes our ability to monitor spatiotemporal changes in interaction diversity.

Dissection of the synthesis of polyunsaturated fatty acids in nematodes and Collembola of the soil fauna.

It is becoming increasingly clear that not only unicellular, photoautotrophic eukaryotes, plants, and fungi, but also invertebrates are capable of synthesizing ω3 long-chain polyunsaturated fatty acids (LC-PUFA) de novo. However, the distribution of this anabolic capacity among different invertebrate groups and its implementation at the gene and protein level are often still unknown. This study investigated the PUFA pathways in common soil fauna, i.e. two nematode and two Collembola species. Of these, one species each (Panagrellus redivivus, Folsomia candida) was assumed to produce ω3 LC-PUFA de novo, while the others (Acrobeloides bodenheimeri, Isotoma caerulea) were supposed to be unable to do so. A highly labeled oleic acid (99 % 13C) was supplemented and the isotopic signal was used to trace its metabolic path. All species followed the main pathway of lipid biosynthesis. However, in A. bodenheimeri this terminated at arachidonic acid (ω6 PUFA), whereas the other three species continued the pathway to eicosapentaenoic acid (ω3 PUFA), including I. caerulea. For the nematode P. redivivus the identification and functional characterization of four new fatty acid desaturase (FAD) genes was performed. These genes encode the FAD activities Δ9, Δ6, and Δ5, respectively. Additionally, the Δ12 desaturase was analyzed, yet the observed activity of an ω3 FAD could not be attributed to a coding gene. In the Collembola F. candida, 11 potential first desaturases (Δ9) and 13 front-end desaturases (Δ6 or Δ5 FADs) have been found. Further sequence analysis indicates the presence of omega FADs, specifically Δ12, which are likely derived from Δ9 FADs.

Whole-genome-based phylogenetic analyses provide new insights into the evolution of springtails (Hexapoda: Collembola).

Springtails (Collembola) stand as one of the most abundant, widespread, and ancient terrestrial arthropods on earth. However, their evolutionary history and deep phylogenetic relationships remain elusive. In this study, we employed phylogenomic approaches to elucidate the basal relationships among Collembola. We sampled whole-genome data representing all major collembolan lineages in proportion to their known diversity. To account for potential phylogenomic biases, we implemented various data extraction, locus sampling, and signal filtering strategies to generate matrices. Subsequently, we applied a diverse array of tree-searching and rate-modelling methods to reconstruct the phylogeny. Our analyses, utilizing different matrices and methods, converged on the same unrooted relationships among collembolan ingroups, supporting the current ordinal classification and challenging the monophyly of Arthropleona and Symphypleona s.l. However, discrepancies across analyses existed in the root of Collembola. Among various root positions, those based on more informative matrices and biologically realistic models, favoring a basal topology of Entomobryomorpha + (Symphypleona s.s. + (Neelipleona + Poduromorpha)), were supported by subsequent methodological assessment, topology tests, and rooting analyses. This optimal topology suggests multiple independent reduction of the pronotum in non-poduromorph orders and aligns with the plesiomorphic status of neuroendocrine organs and epicuticular structure of Entomobryomorpha. Fossil-calibrated dating analyses based on the optimal topology indicated late-Paleozoic to mid-Mesozoic origins of the crown Collembola and four orders. In addition, our results questioned the monophyly of Isotomidae and Neanuridae, underscoring the need for further attention to the systematics of these families. Overall, this study provides novel insights into the phylogenetic backbone of Collembola, which will inform future studies on the systematics, ecology, and evolution of this significant arthropod lineage.

Gene editing in agricultural, health, and veterinary pest arthropods: recent advances.

Pest arthropods cause significant crop damage or are vectors of pathogens for both plants and animals. The current standard of pest management prevents against crop losses and protects human and animal health, but shortcomings exist, such as insecticide resistance and environmental damage to nontarget organisms. New management methods are therefore needed. The development of new tools, such as site-specific gene editing, has accelerated the study of gene function and phenotype in nonmodel arthropod species and may enable the development of new strategies for pathogen and arthropod control. Here, the most recent developments in gene editing in arthropod pests are briefly reviewed. Additionally, technological advances that could be applicable to new species or enhance the success rates of gene editing in species with already established protocols are highlighted.

Selection Across the Three-Dimensional Structure of Venom Proteins from North American Scolopendromorph Centipedes.

Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.

Assessing arthropod biodiversity with DNA barcoding in Jinnah Garden, Lahore, Pakistan.

Previous difficulties in arthropod taxonomy (such as limitations in conventional morphological approaches, the possibility of cryptic species and a shortage of knowledgeable taxonomists) has been overcome by the powerful tool of DNA barcoding. This study presents a thorough analysis of DNA barcoding in regards to Pakistani arthropods, which were collected from Lahore's Jinnah Garden. The 88 % (9,451) of the 10,792 specimens that were examined were able to generate DNA barcodes and 83% (8,974) of specimens were assigned 1,361 barcode index numbers (BINs). However, the success rate differed significantly between the orders of arthropods, from 77% for Thysanoptera to an astounding 93% for Diptera. Through morphological exams, DNA barcoding, and cross-referencing with the Barcode of Life Data system (BOLD), the Barcode Index Numbers (BINs) were assigned with a high degree of accuracy, both at the order (100%) and family (98%) levels. Though, identifications at the genus (37%) and species (15%) levels showed room for improvement. This underscores the ongoing need for enhancing and expanding the DNA barcode reference library. This study identified 324 genera and 191 species, underscoring the advantages of DNA barcoding over traditional morphological identification methods. Among the 17 arthropod orders identified, Coleoptera, Diptera, Hemiptera, Hymenoptera, and Lepidoptera from the class Insecta dominated, collectively constituting 94% of BINs. Expected malaise trap Arthropod fauna in Jinnah Garden could contain approximately 2,785 BINs according to Preston log-normal species distribution, yet the Chao-1 Index predicts 2,389.74 BINs. The Simpson Index of Diversity (1-D) is 0.989, signaling high species diversity, while the Shannon Index is 5.77, indicating significant species richness and evenness. These results demonstrated that in Pakistani arthropods, DNA barcoding and BOLD are an invaluable tool for improving taxonomic understanding and biodiversity assessment, opening the door for further eDNA and metabarcoding research.

Lower Ordovician synziphosurine reveals early euchelicerate diversity and evolution.

Euchelicerata is a clade of arthropods comprising horseshoe crabs, scorpions, spiders, mites and ticks, as well as the extinct eurypterids (sea scorpions) and chasmataspidids. The understanding of the ground plans and relationships between these crown-group euchelicerates has benefited from the discovery of numerous fossils. However, little is known regarding the origin and early evolution of the euchelicerate body plan because the relationships between their Cambrian sister taxa and synziphosurines, a group of Silurian to Carboniferous stem euchelicerates with chelicerae and an unfused opisthosoma, remain poorly understood owing to the scarce fossil record of appendages. Here we describe a synziphosurine from the Lower Ordovician (ca. 478 Ma) Fezouata Shale of Morocco. This species possesses five biramous appendages with stenopodous exopods bearing setae in the prosoma and a fully expressed first tergite in the opisthosoma illuminating the ancestral anatomy of the group. Phylogenetic analyses recover this fossil as a member of the stem euchelicerate family Offacolidae, which is characterized by biramous prosomal appendages. Moreover, it also shares anatomical features with the Cambrian euarthropod Habelia optata, filling the anatomical gap between euchelicerates and Cambrian stem taxa, while also contributing to our understanding of the evolution of euchelicerate uniramous prosomal appendages and tagmosis.

Annotation and visualization of parasite, fungi and arthropod genomes with Companion.

As sequencing genomes has become increasingly popular, the need for annotation of the resulting assemblies is growing. Structural and functional annotation is still challenging as it includes finding the correct gene sequences, annotating other elements such as RNA and being able to submit those data to databases to share it with the community. Compared to de novo assembly where contiguous chromosomes are a sign of high quality, it is difficult to visualize and assess the quality of annotation. We developed the Companion web server to allow non-experts to annotate their genome using a reference-based method, enabling them to assess the output before submitting to public databases. In this update paper, we describe how we have included novel methods for gene finding and made the Companion server more efficient for annotation of genomes of up to 1 Gb in size. The reference set was increased to include genomes of interest for human and animal health from the fungi and arthropod kingdoms. We show that Companion outperforms existing comparable tools where closely related references are available.

The role of ATP-binding cassette transporters in arthropod pesticide toxicity and resistance.

Pesticide resistance in arthropods threatens agricultural productivity and the control of vector-borne diseases. The ATP-binding cassette (ABC) transporters have emerged as important factors in the toxicity of synthetic pesticides, as well as for Bacillus thuringiensis insecticidal Cry protein binding. Depending on the localization of expression, both higher and lower expression of ABCs have been linked with pesticide resistance. The recent development of genetic-based approaches such as RNAi and CRISPR/Cas9 gene editing in nonmodel species, has greatly contributed to unveil their functional importance in pesticide toxicity and resistance. Using these tools, we are now poised to further unravel the molecular genetic mechanisms of gene regulation uncovering more elusive regulatory resistance genes.

The Patterns of Codon Usage between Chordates and Arthropods are Different but Co-evolving with Mutational Biases.

Different frequencies amongst codons that encode the same amino acid (i.e. synonymous codons) have been observed in multiple species. Studies focused on uncovering the forces that drive such codon usage showed that a combined effect of mutational biases and translational selection works to produce different frequencies of synonymous codons. However, only few have been able to measure and distinguish between these forces that may leave similar traces on the coding regions. Here, we have developed a codon model that allows the disentangling of mutation, selection on amino acids and synonymous codons, and GC-biased gene conversion (gBGC) which we employed on an extensive dataset of 415 chordates and 191 arthropods. We found that chordates need 15 more synonymous codon categories than arthropods to explain the empirical codon frequencies, which suggests that the extent of codon usage can vary greatly between animal phyla. Moreover, methylation at CpG sites seems to partially explain these patterns of codon usage in chordates but not in arthropods. Despite the differences between the two phyla, our findings demonstrate that in both, GC-rich codons are disfavored when mutations are GC-biased, and the opposite is true when mutations are AT-biased. This indicates that selection on the genomic coding regions might act primarily to stabilize its GC/AT content on a genome-wide level. Our study shows that the degree of synonymous codon usage varies considerably among animals, but is likely governed by a common underlying dynamic.

Characterization of the complete mitochondrial genome of Spirobolus grahami (Diplopoda: Spirobolidae) with phylogenetic analysis.

Diplopoda is one of the most diverse and important groups of soil arthropods, but little research has been done on their phylogenetic relationship and evolution. Here, we sequenced and annotated the complete mitochondrial genomes of Spirobolus grahami. The total mitogenome of S. grahami was typical circular, double-stranded molecules, with 14,875 bp in length, including 13 protein-coding genes, 22 tRNAs, two rRNAs, and one control region. Base composition analysis suggested that the mitochondrial sequences were biased toward A and T, with A + T content of 58.68%. The mitogenomes of S. grahami exhibited negative AT and positive GC skews. Most of the 13 PCGs had ATN as the start codon, except COX1 start with CGA, and most PCGs ended with the T stop codon. The dN/dS values for most PCGs were lower than 1, suggesting that purifying selection was likely the main driver of mitochondrial PCG evolution. Phylogenetic analyses based on 13 PCGs using BI and ML methods support the classification of genus Spirobolus and Tropostreptus. Glomeridesmus spelaeus is distantly related to the other Diplopoda species.

A long-read draft assembly of the Chinese mantis (Mantodea: Mantidae: Tenodera sinensis) genome reveals patterns of ion channel gain and loss across Arthropoda.

Praying mantids (Mantodea: Mantidae) are iconic insects that have captivated biologists for decades, especially the species with cannibalistic copulatory behavior. This behavior has been cited as evidence that insects lack nociceptive capacities and cannot feel pain; however, this behaviorally driven hypothesis has never been rigorously tested at the genetic or functional level. To enable future studies of nociceptive capabilities in mantids, we sequenced and assembled a draft genome of the Chinese praying mantis (Tenodera sinensis) and identified multiple classes of nociceptive ion channels by comparison to orthologous gene families in Arthropoda. Our assembly-produced using PacBio HiFi reads-is fragmented (total size = 3.03 Gb; N50 = 1.8 Mb; 4,966 contigs), but is highly complete with respect to gene content (BUSCO complete = 98.7% [odb10_insecta]). The size of our assembly is substantially larger than that of most other insects, but is consistent with the size of other mantid genomes. We found that most families of nociceptive ion channels are present in the T. sinensis genome; that they are most closely related to those found in the damp-wood termite (Zootermopsis nevadensis); and that some families have expanded in T. sinensis while others have contracted relative to nearby lineages. Our findings suggest that mantids are likely to possess nociceptive capabilities and provide a foundation for future experimentation regarding ion channel functions and their consequences for insect behavior.

Genetic underpinnings of arthropod community distributions in Populus trichocarpa.

Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. We surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.