Genetic innovations in animal-microbe symbioses.

Animal hosts have initiated myriad symbiotic associations with microorganisms and often have maintained these symbioses for millions of years, spanning drastic changes in ecological conditions and lifestyles. The establishment and persistence of these relationships require genetic innovations on the parts of both symbionts and hosts. The nature of symbiont innovations depends on their genetic population structure, categorized here as open, closed or mixed. These categories reflect modes of inter-host transmission that result in distinct genomic features, or genomic syndromes, in symbionts. Although less studied, hosts also innovate in order to preserve and control symbiotic partnerships. New capabilities to sequence host-associated microbial communities and to experimentally manipulate both hosts and symbionts are providing unprecedented insights into how genetic innovations arise under different symbiont population structures and how these innovations function to support symbiotic relationships.

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 sea spider Pycnogonum litorale overturns the paradigm of the absence of axial regeneration in molting animals.

Regenerative abilities and their evolution in the different animal lineages have fascinated generations of biologists. While some taxa are capable of restoring entire individuals from small body fragments, others can regrow only specific structures or lack structural regeneration completely. In contrast to many other protostomes, including the segmented annelids, molting animals (Ecdysozoa) are commonly considered incapable of primary body axis regeneration, which has been hypothesized to be linked to the evolution of their protective cuticular exoskeleton. This holds also for the extraordinarily diverse, segmented arthropods. Contradicting this long-standing paradigm, we here show that immatures of the sea spider Pycnogonum litorale reestablish the posterior body pole after transverse amputation and can regrow almost complete segments and the terminal body region, including the hindgut, anus, and musculature. Depending on the amputation level, normal phenotypes or hypomeric six-legged forms develop. Remarkably, also the hypomeric animals regain reproductive functionality by ectopic formation of gonoducts and gonopores. The discovery of such complex regenerative patterns in an extant arthropod challenges the hitherto widely assumed evolutionary loss of axial regeneration during ecdysozoan evolution. Rather, the branching of sea spiders at the base of Chelicerata and their likely ancestral anamorphic development suggests that the arthropod stem species may have featured similar regenerative capabilities. Accordingly, our results provide an incentive for renewed comparative regeneration studies across ecdysozoans, with the aim to resolve whether this trait was potentially even inherited from the protostome ancestor.

A thermal receptor for nonvisual sunlight detection in myriapods.

Organisms from cyanobacteria to humans have evolved a wide array of photoreceptive strategies to detect light. Sunlight avoidance behavior is common in animals without vision or known photosensory genes. While indirect light perception via photothermal conversion is a possible scenario, there is no experimental evidence for this hypothesis. Here, we show a nonvisual and extraocular sunlight detection mechanism by identifying the broad-range thermal receptor 1 (BRTNaC1, temperature range = 33 to 48 °C) in centipede antennae. BRTNaC1, a heat-activated cation-permeable ion channel, is structurally related to members of the epithelial sodium channel family. At the molecular level, heat activation of BRTNaC1 exhibits strong pH dependence controlled by two protonatable sites. Physiologically, temperature-dependent activation of BRTNaC1 upon sunlight exposure comes from a striking photothermal effect on the antennae, where a slightly acidic environment (pH 6.1) of the body fluid leads to the protonation of BRTNaC1 and switches on its high thermal sensitivity. Furthermore, testosterone potently inhibits heat activation of BRTNaC1 and the sunlight avoidance behavior of centipedes. Taken together, our study suggests a sophisticated strategy for nonvisual sunlight detection in myriapods.

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.

Identification of an arthropod molecular target for plant-derived natural repellents.

Arthropods maintain ecosystem balance while also contributing to the spread of disease. Plant-derived natural repellents represent an ecological method of pest control, but their direct molecular targets in arthropods remain to be further elucidated. Occupying a critical phylogenetic niche in arthropod evolution, scorpions retain an ancestral genetic profile. Here, using a behavior-guided screening of the Mesobuthus martensii genome, we identified a scorpion transient receptor potential (sTRP1) channel that senses Cymbopogon-derived natural repellents, while remaining insensitive to the synthetic chemical pesticide DEET. Scrutinizing orthologs of sTRP1 in Drosophila melanogaster, we further demonstrated dTRPγ ion channel as a chemosensory receptor of natural repellents to mediate avoidance behavior. This study sheds light on arthropod molecular targets of natural repellents, exemplifying the arthropod­plant adaptation. It should also help the rational design of insect control strategy and in conserving biodiversity.

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.

Genomic insights into the chromosomal elongation in a family of Collembola.

Collembola is a highly diverse and abundant group of soil arthropods with chromosome numbers ranging from 5 to 11. Previous karyotype studies indicated that the Tomoceridae family possesses an exceptionally long chromosome. To better understand chromosome size evolution in Collembola, we obtained a chromosome-level genome of Yoshiicerus persimilis with a size of 334.44 Mb and BUSCO completeness of 97.0% (n = 1013). Both genomes of Y. persimilis and Tomocerus qinae (recently published) have an exceptionally large chromosome (ElChr greater than 100 Mb), accounting for nearly one-third of the genome. Comparative genomic analyses suggest that chromosomal elongation occurred independently in the two species approximately 10 million years ago, rather than in the ancestor of the Tomoceridae family. The ElChr elongation was caused by large tandem and segmental duplications, as well as transposon proliferation, with genes in these regions experiencing weaker purifying selection (higher dN/dS) than conserved regions. Moreover, inter-genomic synteny analyses indicated that chromosomal fission/fusion events played a crucial role in the evolution of chromosome numbers (ranging from 5 to 7) within Entomobryomorpha. This study provides a valuable resource for investigating the chromosome evolution of Collembola.

Differential viral RNA methylation contributes to pathogen blocking in Wolbachia-colonized arthropods.

Arthropod endosymbiont Wolbachia pipientis is part of a global biocontrol strategy to reduce the replication of mosquito-borne RNA viruses such as alphaviruses. We previously demonstrated the importance of a host cytosine methyltransferase, DNMT2, in Drosophila and viral RNA as a cellular target during pathogen-blocking. Here we report a role for DNMT2 in Wolbachia-induced alphavirus inhibition in Aedes species. Expression of DNMT2 in mosquito tissues, including the salivary glands, is elevated upon virus infection. Notably, this is suppressed in Wolbachia-colonized animals, coincident with reduced virus replication and decreased infectivity of progeny virus. Ectopic expression of DNMT2 in cultured Aedes cells is proviral, increasing progeny virus infectivity, and this effect of DNMT2 on virus replication and infectivity is dependent on its methyltransferase activity. Finally, examining the effects of Wolbachia on modifications of viral RNA by LC-MS show a decrease in the amount of 5-methylcytosine modification consistent with the down-regulation of DNMT2 in Wolbachia colonized mosquito cells and animals. Collectively, our findings support the conclusion that disruption of 5-methylcytosine modification of viral RNA is a vital mechanism operative in pathogen blocking. These data also emphasize the essential role of epitranscriptomic modifications in regulating fundamental alphavirus replication and transmission processes.

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.

Comparative genomic analysis of chemosensory-related gene families in gastropods.

Chemoreception is critical for the survival and reproduction of animals. Except for a reduced group of insects and chelicerates, the molecular identity of chemosensory proteins is poorly understood in invertebrates. Gastropoda is the extant mollusk class with the greatest species richness, including marine, freshwater, and terrestrial lineages, and likely, highly diverse chemoreception systems. Here, we performed a comprehensive comparative genome analysis taking advantage of the chromosome-level information of two Gastropoda species, one of which belongs to a lineage that underwent a whole genome duplication event. We identified thousands of previously uncharacterized chemosensory-related genes, the majority of them encoding G protein-coupled receptors (GPCR), mostly organized into clusters distributed across all chromosomes. We also detected gene families encoding degenerin epithelial sodium channels (DEG-ENaC), ionotropic receptors (IR), sensory neuron membrane proteins (SNMP), Niemann-Pick type C2 (NPC2) proteins, and lipocalins, although with a lower number of members. Our phylogenetic analysis of the GPCR gene family across protostomes revealed: (i) remarkable gene family expansions in Gastropoda; (ii) clades including members from all protostomes; and (iii) species-specific clades with a substantial number of receptors. For the first time, we provide new and valuable knowledge into the evolution of the chemosensory gene families in invertebrates other than arthropods.

Toward a genome sequence for every animal: Where are we now?

In less than 25 y, the field of animal genome science has transformed from a discipline seeking its first glimpses into genome sequences across the Tree of Life to a global enterprise with ambitions to sequence genomes for all of Earth's eukaryotic diversity [H. A. Lewin et al., Proc. Natl. Acad. Sci. U.S.A. 115, 4325-4333 (2018)]. As the field rapidly moves forward, it is important to take stock of the progress that has been made to best inform the discipline's future. In this Perspective, we provide a contemporary, quantitative overview of animal genome sequencing. We identified the best available genome assemblies in GenBank, the world's most extensive genetic database, for 3,278 unique animal species across 24 phyla. We assessed taxonomic representation, assembly quality, and annotation status for major clades. We show that while tremendous taxonomic progress has occurred, stark disparities in genomic representation exist, highlighted by a systemic overrepresentation of vertebrates and underrepresentation of arthropods. In terms of assembly quality, long-read sequencing has dramatically improved contiguity, whereas gene annotations are available for just 34.3% of taxa. Furthermore, we show that animal genome science has diversified in recent years with an ever-expanding pool of researchers participating. However, the field still appears to be dominated by institutions in the Global North, which have been listed as the submitting institution for 77% of all assemblies. We conclude by offering recommendations for improving genomic resource availability and research value while also broadening global representation.

Genomic and transcriptomic analyses support a silk gland origin of spider venom glands.

BACKGROUND: Spiders comprise a hyperdiverse lineage of predators with venom systems, yet the origin of functionally novel spider venom glands remains unclear. Previous studies have hypothesized that spider venom glands originated from salivary glands or evolved from silk-producing glands present in early chelicerates. However, there is insufficient molecular evidence to indicate similarity among them. Here, we provide comparative analyses of genome and transcriptome data from various lineages of spiders and other arthropods to advance our understanding of spider venom gland evolution. RESULTS: We generated a chromosome-level genome assembly of a model spider species, the common house spider (Parasteatoda tepidariorum). Module preservation, GO semantic similarity, and differentially upregulated gene similarity analyses demonstrated a lower similarity in gene expressions between the venom glands and salivary glands compared to the silk glands, which questions the validity of the salivary gland origin hypothesis but unexpectedly prefers to support the ancestral silk gland origin hypothesis. The conserved core network in the venom and silk glands was mainly correlated with transcription regulation, protein modification, transport, and signal transduction pathways. At the genetic level, we found that many genes in the venom gland-specific transcription modules show positive selection and upregulated expressions, suggesting that genetic variation plays an important role in the evolution of venom glands. CONCLUSIONS: This research implies the unique origin and evolutionary path of spider venom glands and provides a basis for understanding the diverse molecular characteristics of venom systems.

Expression of netrin and its receptors uncoordinated-5 and frazzled in arthropods and onychophorans suggests conserved and diverged functions in neuronal pathfinding and synaptogenesis.

BACKGROUND: Development of the nervous system and the correct connection of nerve cells require coordinated axonal pathfinding through an extracellular matrix. Outgrowing axons exhibit directional growth toward or away from external guidance cues such as Netrin. Guidance cues can be detected by growth cones that are located at the end of growing axons through membrane-bound receptors such as Uncoordianted-5 and Frazzled. Binding of Netrin causes reformation of the cytoskeleton and growth of the axon toward (or away from) the source of Netrin production. RESULTS: Here, we investigate the embryonic mRNA expression patterns of netrin genes and their potential receptors, uncoordinated-5 and frazzled in arthropod species that cover all main branches of Arthropoda, that is, Pancrustacea, Myriapoda, and Chelicerata. We also studied the expression patterns in a closely related outgroup species, the onychophoran Euperipatoides kanangrensis, and provide data on expression profiles of these genes in larval tissues of the fly Drosophila melanogaster including the brain and the imaginal disks. CONCLUSION: Our data reveal conserved and diverged aspects of neuronal guidance in Drosophila with respect to the other investigated species and suggest a conserved function in nervous system patterning of the developing appendages.

The identification, adaptive evolutionary analyses and mRNA expression levels of homeobox (hox) genes in the Chinese mitten crab Eriocheir sinensis.

BACKGROUND: Arthropods are the largest group in the animal kingdom and are morphologically characterized by heterorhythmic segments. Brachyuran decapod crustaceans undergo brachyurization metamorphosis in the early developmental process, characterized by a reduced abdomen that is folded beneath the cephalothorax and inserted between the pereiopods or in a special cavity. As the main cause of major alterations in the evolution of animal body plans, Hox genes encode transcription factors and are involved in bilaterian anterior-posterior axis patterning. RESULTS: We found eight Hox genes (labial, proboscipedia, Deformed, zerknüllt, Sex combs reduced, Antennapedia, Ultrabithorax, fushi tarazu, abdominal-A and Abdominal-B) in Eriocheir sinensis. The phylogenetic topology of 13 arthropod Hox genes was closely related to traditional taxonomic groupings. Genome collinearity analysis was performed using genomic data and chromosomal location data of E. sinensis and Portunus trituratus. We found that their chromosomes were highly collinear, and there was a corresponding collinear relationship between the three Hox genes (lab, ftz and Abd-B). The mRNA expression levels of Scr and Antp fluctuated significantly in different developmental stages of E. sinensis, especially in the brachyurization stages. Evolutionary analysis indicated the presence of positively selected sites in Ubx. CONCLUSIONS: In this study, we used genome-wide analysis to identify and analyze all members of the Hox genes in E. sinensis. Our data will contribute to a better understanding of Hox genes in E. sinensis and provide useful molecular evolutionary information for further investigation on their roles in the brachyurization of crabs.

Expression of posterior Hox genes and opisthosomal appendage development in a mygalomorph spider.

Spiders represent an evolutionary successful group of chelicerate arthropods. The body of spiders is subdivided into two regions (tagmata). The anterior tagma, the prosoma, bears the head appendages and four pairs of walking legs. The segments of the posterior tagma, the opisthosoma, either lost their appendages during the course of evolution or their appendages were substantially modified to fulfill new tasks such as reproduction, gas exchange, and silk production. Previous work has shown that the homeotic Hox genes are involved in shaping the posterior appendages of spiders. In this paper, we investigate the expression of the posterior Hox genes in a tarantula that possesses some key differences of posterior appendages compared to true spiders, such as the lack of the anterior pair of spinnerets and a second set of book lungs instead of trachea. Based on the observed differences in posterior Hox gene expression in true spiders and tarantulas, we argue that subtle changes in the Hox gene expression of the Hox genes abdA and AbdB are possibly responsible for at least some of the morphological differences seen in true spiders versus tarantulas.

Advances in genome sequencing reveal changes in gene content that contribute to arthropod macroevolution.

Current sequencing technology allows for the relatively affordable generation of highly contiguous genomes. Technological advances have made it possible for researchers to investigate the consequences of diverse sorts of genomic variants, such as gene gain and loss. With the extraordinary number of high-quality genomes now available, we take stock of how these genomic variants impact phenotypic evolution. We take care to point out that the identification of genomic variants of interest is only the first step in understanding their impact. Painstaking lab or fieldwork is still required to establish causal relationships between genomic variants and phenotypic evolution. We focus mostly on arthropod research, as this phylum has an impressive degree of phenotypic diversity and is also the subject of much evolutionary genetics research. This article is intended to both highlight recent advances in the field and also to be a primer for learning about evolutionary genetics and genomics.

Virome of Bat-Infesting Arthropods: Highly Divergent Viruses in Different Vectors.

Bats are reservoirs of important zoonotic viruses like Nipah and SARS viruses. However, whether the blood-sucking arthropods on the body surface of bats also carry these viruses and the relationship between viruses carried by the blood-sucking arthropods and viruses carried by bats have not been reported. This study collected 686 blood-sucking arthropods on the body surface of bats from Yunnan Province, China, between 2012 and 2015, and they included wingless bat flies, bat flies, ticks, mites, and fleas. The viruses carried by these arthropods were analyzed using a meta-transcriptomic approach, and 144 highly diverse positive-sense single-stranded RNA, negative-sense single-stranded RNA, and double-stranded RNA viruses were found, of which 138 were potentially new viruses. These viruses were classified into 14 different virus families or orders, including Bunyavirales, Mononegavirales, Reoviridae, and Picornavirales. Further analyses found that Bunyavirales were the most abundant virus group (84% of total virus RNA) in ticks, whereas narnaviruses were the most abundant (52 to 92%) in the bat flies and wingless bat flies libraries, followed by solemoviruses (1 to 29%) and reoviruses (0 to 43%). These viruses were highly structured based on the arthropod types. It is worth noting that no bat-borne zoonotic viruses were found in the virome of bat-infesting arthropod, seemingly not supporting that bat surface arthropods are vectors of zoonotic viruses carried by bats. IMPORTANCE Bats are reservoirs of many important viral pathogens. To evaluate whether bat-parasitic blood-sucking arthropods participate in the circulation of these important viruses, it is necessary to conduct unbiased virome studies on these arthropods. We evaluated five types of blood-sucking parasitic arthropods on the surface of bats in Yunnan, China, and identified a variety of viruses, some of which had high prevalence and abundance levels, although there is limited overlap in virome between distant arthropods. While most of the virome discovered here is potentially arthropod-specific viruses, we identified three possible arboviruses, including one orthobunyavirus and two vesiculoviruses (family Rhabdoviridae), suggesting bat-parasitic arthropods carry viruses with risk of spillage, which warrants further study.

Genomic characterization of 99 viruses from the bunyavirus families Nairoviridae, Peribunyaviridae, and Phenuiviridae, including 35 previously unsequenced viruses.

Bunyaviruses (Negarnaviricota: Bunyavirales) are a large and diverse group of viruses that include important human, veterinary, and plant pathogens. The rapid characterization of known and new emerging pathogens depends on the availability of comprehensive reference sequence databases that can be used to match unknowns, infer evolutionary relationships and pathogenic potential, and make response decisions in an evidence-based manner. In this study, we determined the coding-complete genome sequences of 99 bunyaviruses in the Centers for Disease Control and Prevention's Arbovirus Reference Collection, focusing on orthonairoviruses (family Nairoviridae), orthobunyaviruses (Peribunyaviridae), and phleboviruses (Phenuiviridae) that either completely or partially lacked genome sequences. These viruses had been collected over 66 years from 27 countries from vertebrates and arthropods representing 37 genera. Many of the viruses had been characterized serologically and through experimental infection of animals but were isolated in the pre-sequencing era. We took advantage of our unusually large sample size to systematically evaluate genomic characteristics of these viruses, including reassortment, and co-infection. We corroborated our findings using several independent molecular and virologic approaches, including Sanger sequencing of 197 genome segments, and plaque isolation of viruses from putative co-infected virus stocks. This study contributes to the described genetic diversity of bunyaviruses and will enhance the capacity to characterize emerging human pathogenic bunyaviruses.

Ecological network structure in response to community assembly processes over evolutionary time.

The dynamic structure of ecological communities results from interactions among taxa that change with shifts in species composition in space and time. However, our ability to study the interplay of ecological and evolutionary processes on community assembly remains relatively unexplored due to the difficulty of measuring community structure over long temporal scales. Here, we made use of a geological chronosequence across the Hawaiian Islands, representing 50 years to 4.15 million years of ecosystem development, to sample 11 communities of arthropods and their associated plant taxa using semiquantitative DNA metabarcoding. We then examined how ecological communities changed with community age by calculating quantitative network statistics for bipartite networks of arthropod-plant associations. The average number of interactions per species (linkage density), ratio of plant to arthropod species (vulnerability) and uniformity of energy flow (interaction evenness) increased significantly in concert with community age. The index of specialization H 2 ' has a curvilinear relationship with community age. Our analyses suggest that younger communities are characterized by fewer but stronger interactions, while biotic associations become more even and diverse as communities mature. These shifts in structure became especially prominent on East Maui (~0.5 million years old) and older volcanos, after enough time had elapsed for adaptation and specialization to act on populations in situ. Such natural progression of specialization during community assembly is probably impeded by the rapid infiltration of non-native species, with special risk to younger or more recently disturbed communities that are composed of fewer specialized relationships.