Chance and predictability in evolution: The genomic basis of convergent dietary specializations in an adaptive radiation

Discovery of additional genomic segments reveals the fluidity of jingmenvirus genomic organization

Jingmenviruses are a distinct group of flavi-like viruses characterized by a genome consisting of four to five segments. Here, we report the discovery of three novel putative jingmenviruses, identified by mining publicly available metagenomics data from mosquito and arachnid samples. Strikingly, these novel jingmenvirus sequences contain up to six genomic segments, with pairs of homologous segments coding for putative structural proteins. Following this discovery, we found an additional homologous segment for two other jingmenvirus genomes, which had gone unnoticed in the initial publications. The presence of a single version of the segments coding for non-structural proteins suggests that we have indeed identified jingmenviruses with infectious units that contain up to six segments. We compared these novel jingmenvirus sequences to published sequences, in particular the segments with multiple open reading frames (ORFs), and we propose that the putative translation initiation mechanisms involved for these segments are ribosomal frameshift resulting in the fusion of ORFs and leaky scanning for overlapping ORFs. These putative mechanisms, conserved for all jingmenvirus sequences analysed, including in homologous segments, require biological confirmation. We also generated structural models of two putative structural proteins in the duplicated segments, and the corresponding alignments enabled us to confirm or identify the homologous relationship between sequences that shared limited nucleotide or amino acid identity. Altogether, these results highlight the fluid nature of jingmenviruses, which is a hallmark of multipartite viruses. Different combinations of segments packaged in different virus particles could facilitate the acquisition or loss of genomic segments and a segment duplication following genomic drift. Our data therefore contribute to the evidence of the multipartite nature of jingmenviruses and the evolutionary role this organization may play.

A targeted gene phylogenetic framework to investigate diversification in the highly diverse yet geographically restricted red devil spiders (Araneae, Dysderidae)

The family Dysderidae is a highly diverse group of nocturnal ground-dwelling and active-hunter spiders. Dysderids are mostly restricted to the Western Palearctic, and particularly rich and abundant around the Mediterranean region. Interestingly, the distribution of species richness among its 24 genera and three subfamilies is highly biased-80% of its 644 documented species belong to just two genera, Dysdera (326) and Harpactea (211). Dysderidae provides an excellent study case for evolutionary and ecological research. It includes cases of trophic specialization, which are uncommon among spiders, and exhibit other remarkable biological (e.g. holocentric chromosomes), behavioural (e.g. cryptic female choice), evolutionary (e.g. adaptive radiation) and ecological features (e.g. recurrent colonization of the subterranean environment). The lack of a quantitative hypothesis on its phylogenetic structure has hampered its potential as a testing ground for evolutionary, biogeographical and ecological hypotheses. Here, we present the results of a target, multi-locus phylogenetic analysis, using mitochondrial (cox1, 16s and 12s) and nuclear genes (h3, 28s and 18s), of the most exhaustive taxonomic sample within Dysderidae (104 spp.) to date and across related families (Synspermiata) (83 spp.). We estimate divergence times using a combination of fossil and biogeographic node calibrations and use this timeline to identify shifts in diversification rates. Our results support the monophyly of the Dysderidae subfamilies Rhodinae and Dysderinae but reject Harpacteinae as currently defined. Moreover, the clades recovered within Harpacteinae do not support its current taxonomy. The origin of the family most likely post-dated the break-up of Pangea, and cave colonization may be older than previously considered. After correcting for the taxonomic artefacts, we identified a significant shift in diversification rates at the base of the genus Dysdera. Although the unique coexistence of specialist and generalist diets within the lineage could be suggested as the potential driver for the rate acceleration, further quantitative analyses would be necessary to test this hypothesis.

Population Genomics of Adaptive Radiations: Exceptionally High Levels of Genetic Diversity and Recombination in an Endemic Spider From the Canary Islands

The spider genus Dysdera has undergone a remarkable diversification in the oceanic archipelago of the Canary Islands, with ~60 endemic species having originated during the 20 million years since the origin of the archipelago. This evolutionary radiation has been accompanied by substantial dietary shifts, often characterised by phenotypic modifications encompassing morphological, metabolic and behavioural changes. Hence, these endemic spiders represent an excellent model for understanding the evolutionary drivers and to pinpoint the genomic determinants underlying adaptive radiations. Recently, we achieved the first chromosome-level genome assembly of one of the endemic species, D. silvatica, providing a high-quality reference sequence for evolutionary genomics studies. Here, we conducted a low coverage-based resequencing study of a natural population of D. silvatica from La Gomera island. Taking advantage of the new high-quality genome, we characterised genome-wide levels of nucleotide polymorphism, divergence and linkage disequilibrium, and inferred the demographic history of this population. We also performed comprehensive genome-wide scans for recent positive selection. Our findings uncovered exceptionally high levels of nucleotide diversity and recombination in this geographically restricted endemic species, indicative of large historical effective population sizes. We also identified several candidate genomic regions that are potentially under positive selection, highlighting relevant biological processes, such as vision and nitrogen extraction as potential adaptation targets. These processes may ultimately drive species diversification in this genus. This pioneering study of spiders that are endemic to an oceanic archipelago lays the groundwork for broader population genomics analyses aimed at understanding the genetic mechanisms driving adaptive radiation in island ecosystems.

Resolving Marine-Freshwater Transitions by Diatoms Through a Fog of Gene Tree Discordance

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms.

Insights into the Genomics of Clownfish Adaptive Radiation: The Genomic Substrate of the Diversification

Clownfishes are an iconic group of coral reef fishes that evolved a mutualistic interaction with sea anemones, which triggered the rapid diversification of the group. Following the emergence of this mutualism, clownfishes diversified into different ecological niches and developed convergent phenotypes associated with their host use. The genetic basis of the initial acquisition of the mutualism with host anemones has been described, but the genomic architecture underlying clownfish diversification once the mutualism was established and the extent to which clownfish phenotypic convergence originated through shared genetic mechanisms are still unknown. Here, we investigated these questions by performing comparative genomic analyses on the available genomic data of five pairs of closely related but ecologically divergent clownfish species. We found that clownfish diversification was characterized by bursts of transposable elements, an overall accelerated coding evolution, incomplete lineage sorting, and ancestral hybridization events. Additionally, we detected a signature of positive selection in 5.4% of the clownfish genes. Among them, five presented functions associated with social behavior and ecology, and they represent candidate genes involved in the evolution of the size-based hierarchical social structure so particular to clownfishes. Finally, we found genes with patterns of either relaxation or intensification of purifying selection and signals of positive selection linked with clownfish ecological divergence, suggesting some level of parallel evolution during the diversification of the group. Altogether, this work provides the first insights into the genomic substrate of clownfish adaptive radiation and integrates the growing collection of studies investigating the genomic mechanisms governing species diversification.

Evolutionary genomics of oceanic island radiations

A recurring feature of oceanic archipelagos is the presence of adaptive radiations that generate endemic, species-rich clades that can offer outstanding insight into the links between ecology and evolution. Recent developments in evolutionary genomics have contributed towards solving long-standing questions at this interface. Using a comprehensive literature search, we identify studies spanning 19 oceanic archipelagos and 110 putative adaptive radiations, but find that most of these radiations have not yet been investigated from an evolutionary genomics perspective. Our review reveals different gaps in knowledge related to the lack of implementation of genomic approaches, as well as undersampled taxonomic and geographic areas. Filling those gaps with the required data will help to deepen our understanding of adaptation, speciation, and other evolutionary processes.

Topographically distinct adaptive landscapes for teeth, skeletons, and size explain the adaptive radiation of Carnivora (Mammalia)

Models of adaptive radiation were originally developed to explain the early, rapid appearance of distinct modes of life within diversifying clades. Phylogenetic tests of this hypothesis have yielded limited support for temporally declining rates of phenotypic evolution across diverse clades, but the concept of an adaptive landscape that links form to fitness, while also crucial to these models, has received more limited attention. Using methods that assess the temporal accumulation of morphological variation and estimate the topography of the underlying adaptive landscape, I found evidence of an early partitioning of mandibulo-dental morphological variation in Carnivora (Mammalia) that occurs on an adaptive landscape with multiple peaks, consistent with classic ideas about adaptive radiation. Although strong support for this mode of adaptive radiation is present in traits related to diet, its signal is not present in body mass data or for traits related to locomotor behavior and substrate use. These findings suggest that adaptive radiations may occur along some axes of ecomorphological variation without leaving a signal in others and that their dynamics are more complex than simple univariate tests might suggest.

Evolutionary transitions in diet influence the exceptional diversification of a lizard adaptive radiation

Background

Diet is a key component of a species ecological niche and plays critical roles in guiding the trajectories of evolutionary change. Previous studies suggest that dietary evolution can influence the rates and patterns of species diversification, with omnivorous (animal and plant, ‘generalist’) diets slowing down diversification compared to more restricted (‘specialist’) herbivorous and carnivorous diets. This hypothesis, here termed the “dietary macroevolutionary sink” hypothesis (DMS), predicts that transitions to omnivorous diets occur at higher rates than into any specialist diet, and omnivores are expected to have the lowest diversification rates, causing an evolutionary sink into a single type of diet. However, evidence for the DMS hypothesis remains conflicting. Here, we present the first test of the DMS hypothesis in a lineage of ectothermic tetrapods—the prolific Liolaemidae lizard radiation from South America.

Results

Ancestral reconstructions suggest that the stem ancestor was probably insectivorous. The best supported trait model is a diet-dependent speciation rate, with independent extinction rates. Herbivory has the highest net diversification rate, omnivory ranks second, and insectivory has the lowest. The extinction rate is the same for all three diet types and is much lower than the speciation rates. The highest transition rate was from omnivory to insectivory, and the lowest transition rates were between insectivory and herbivory.

Conclusions

Our findings challenge the core prediction of the DMS hypothesis that generalist diets represent an ‘evolutionary sink’. Interestingly, liolaemid lizards have rapidly and successfully proliferated across some of the world’s coldest climates (at high elevations and latitudes), where species have evolved mixed arthropod-plant (omnivore) or predominantly herbivore diets. This longstanding observation is consistent with the higher net diversification rates found in both herbivory and omnivory. Collectively, just like the evolution of viviparity has been regarded as a ‘key adaptation’ during the liolaemid radiation across cold climates, our findings suggest that transitions from insectivory to herbivory (bridged by omnivory) are likely to have played a role as an additional key adaptation underlying the exceptional diversification of these reptiles across extreme climates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02028-3.

The Tetragnatha kauaiensis Genome Sheds Light on the Origins of Genomic Novelty in Spiders

Spiders (Araneae) have a diverse spectrum of morphologies, behaviors, and physiologies. Attempts to understand the genomic-basis of this diversity are often hindered by their large, heterozygous, and AT-rich genomes with high repeat content resulting in highly fragmented, poor-quality assemblies. As a result, the key attributes of spider genomes, including gene family evolution, repeat content, and gene function, remain poorly understood. Here, we used Illumina and Dovetail Chicago technologies to sequence the genome of the long-jawed spider Tetragnatha kauaiensis, producing an assembly distributed along 3,925 scaffolds with an N50 of ∼2 Mb. Using comparative genomics tools, we explore genome evolution across available spider assemblies. Our findings suggest that the previously reported and vast genome size variation in spiders is linked to the different representation and number of transposable elements. Using statistical tools to uncover gene-family level evolution, we find expansions associated with the sensory perception of taste, immunity, and metabolism. In addition, we report strikingly different histories of chemosensory, venom, and silk gene families, with the first two evolving much earlier, affected by the ancestral whole genome duplication in Arachnopulmonata (∼450 Ma) and exhibiting higher numbers. Together, our findings reveal that spider genomes are highly variable and that genomic novelty may have been driven by the burst of an ancient whole genome duplication, followed by gene family and transposable element expansion.

Integrated miRNA and transcriptome profiling to explore the molecular determinism of convergent adaptation to corn in two lepidopteran pests of agriculture

Background

The degree to which adaptation to same environment is determined by similar molecular mechanisms, is a topic of broad interest in evolutionary biology, as an indicator of evolutionary predictability. We wished to address if adaptation to the same host plant in phytophagous insects involved related gene expression patterns. We compared sRNA-Seq and RNA-Seq data between two pairs of taxa of Ostrinia and Spodoptera frugiperda sharing maize as host-plant. For the latter, we had previously carried out a reciprocal transplant experiment by feeding of the larvae of the Corn strain (Sf-C) and the Rice strain (Sf-R) on corn versus rice and characterized the mRNA and miRNA responses.

Results

First, we predicted the genes encoding miRNA in Ostrinia nubilalis (On) and O. scapulalis (Os). Respectively 67 and 65 known miRNA genes, as well as 196 and 190 novel ones were predicted with Os genome using sncRNAs extracted from whole larvae feeding on corn or mugwort. In On, a read counts analysis showed that 37 (55.22%) known miRNAs and 19 (9.84%) novel miRNAs were differentially expressed (DE) on mugwort compared to corn (in Os, 25 known miRs (38.46%) and 8 novel ones (4.34%)). Between species on corn, 8 (12.5%) known miRNAs and 8 (6.83%) novel ones were DE while only one novel miRNA showed expression variation between species on mugwort. Gene target prediction led to the identification of 2953 unique target genes in On and 2719 in Os, among which 11.6% (344) were DE when comparing species on corn. 1.8% (54) of On miR targets showed expression variation upon a change of host-plant.

We found molecular changes matching convergent phenotype, i.e., a set of nine miRNAs that are regulated either according to the host-plant both in On and Sf-C or between them on the same plant, corn. Among DE miR target genes between taxa, 13.7% shared exactly the same annotation between the two pairs of taxa and had function related to insect host-plant interaction.

Conclusion

There is some similarity in underlying genetic mechanisms of convergent evolution of two distant Lepidopteran species having adopted corn in their host range, highlighting possible adaptation genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07905-7.

The chromosome-scale assembly of the Canary Islands endemic spider Dysdera silvatica (Arachnida, Araneae) sheds light on the origin and genome structure of chemoreceptor gene families in chelicerates

Here, we present the chromosome-level genome assembly of Dysdera silvatica Schmidt, 1981, a nocturnal ground-dwelling spider endemic from the Canary Islands. The genus Dysdera has undergone a remarkable diversification in this archipelago mostly associated with shifts in the level of trophic specialization, becoming an excellent model to study the genomic drivers of adaptive radiations. The new assembly (1.37 Gb; scaffold N50 of 174.2 Mb), was performed using the chromosome conformation capture scaffolding technique, represents a continuity improvement of more than 4500 times with respect to the previous version. The seven largest scaffolds or pseudochromosomes, which cover 87% of the total assembly size, probably correspond with the seven chromosomes of the karyotype of this species, including a characteristic large X chromosome. To illustrate the value of this new resource we performed a comprehensive analysis of the two major arthropod chemoreceptor gene families (i.e., gustatory and ionotropic receptors). We identified 545 chemoreceptor sequences distributed across all pseudochromosomes, with a notable underrepresentation in the X chromosome. At least 54% of them localize in 83 genomic clusters with a significantly lower evolutionary distances between them than the average of the family, suggesting a recent origin of many of them. This chromosome-level assembly is the first high-quality genome representative of the Synspermiata clade, and just the third among spiders, representing a new valuable resource to gain insights into the structure and organization of chelicerate genomes, including the role that structural variants, repetitive elements and large gene families played in the extraordinary biology of spiders.

Prey acceptance and metabolic specialisations in some Canarian Dysdera spiders

Spiders of the genus Dysdera are peculiar for their preying on terrestrial isopods though the preference for this prey type varies between species. We tested prey acceptance of two isopod and two non-isopod prey in 8 species endemic to the Canary Islands, and analyzed growth and metabolic parameters (growth efficiency; dry mass, lipid and N extraction efficiency; lipid:protein consumption ratio) of 6 of these species when fed either house flies (Musca domestica) or isopods (Porcellio scaber) in the laboratory. The species represented four morphological types (unmodified chelicerae, slightly elongated chelicerae, concave chelicerae, flattened fang), supposedly reflecting different specializations to isopod prey. The results showed reduced relative acceptance of non-isopod and increased acceptance of isopod prey in species groups with specialised morphologies compared to the unmodified species group. All species had similar or lower growth and growth efficiency when feeding on isopods than on flies. Extraction efficiency of dry mass and lipid were higher for flies than for isopods, while extraction efficiency of protein was higher for isopods than for flies. All species also utilized isopod protein equally well, but protein utilization of flies was lower in the presumed specialist compared to generalist species, indicating a possible metabolic trade-off from isopod specialisation. Thus, morphological adaptations were associated with increased behavioural preferences for isopods and reduced metabolic ability to handle 'generalist' prey.

Mosaic patterns of homoplasy accompany the parallel evolution of suspensory adaptations in the forelimb of tree sloths (Folivora: Xenarthra)

We examine how derived functional signal and phylogenetic inheritance interact in the forelimb of tree sloths, to understand the relative contribution of each in the evolution of a novel morphobehavioural suite. Molecular and craniodental data demonstrate that extant tree sloths evolved suspensory behaviours and associated morphologies from a non-suspensory ancestor independently of one another, making them a useful model system. We find that convergence in univariate traits is expressed mosaically, although the signal is largely functional. Three-dimensional analyses suggest there is greater conservatism of gross morphology in more proximal bones than in more distal elements. Convergence in some univariate scapular traits is independent of the gross morphology of the scapula itself, demonstrating that functionally relevant morphologies were mapped on to a more conserved scapular shape. Our results suggest that morphological homoplasy is expressed in a mosaic manner. The relationship between homoplasy and trait integration may be more nuanced than previously thought, even within a single adaptive system.

Evolutionary History of Major Chemosensory Gene Families across Panarthropoda

Chemosensory perception is a fundamental biological process of particular relevance in basic and applied arthropod research. However, apart from insects, there is little knowledge of specific molecules involved in this system, which is restricted to a few taxa with uneven phylogenetic sampling across lineages. From an evolutionary perspective, onychophorans (velvet worms) and tardigrades (water bears) are of special interest since they represent the closest living relatives of arthropods, altogether comprising the Panarthropoda. To get insights into the evolutionary origin and diversification of the chemosensory gene repertoire in panarthropods, we sequenced the antenna- and head-specific transcriptomes of the velvet worm Euperipatoides rowelli and analyzed members of all major chemosensory families in representative genomes of onychophorans, tardigrades, and arthropods. Our results suggest that the NPC2 gene family was the only family encoding soluble proteins in the panarthropod ancestor and that onychophorans might have lost many arthropod-like chemoreceptors, including the highly conserved IR25a receptor of protostomes. On the other hand, the eutardigrade genomes lack genes encoding the DEG-ENaC and CD36-sensory neuron membrane proteins, the chemosensory members of which have been retained in arthropods; these losses might be related to lineage-specific adaptive strategies of tardigrades to survive extreme environmental conditions. Although the results of this study need to be further substantiated by an increased taxon sampling, our findings shed light on the diversification of chemosensory gene families in Panarthropoda and contribute to a better understanding of the evolution of animal chemical senses.

Woodlice change the habitat use of spiders in a different food chain

Background

In old field systems, the common woodlouse may have an indirect effect on a nursery web spider. Woodlice and nursery web spiders feed in different food chains, yet previous work demonstrated that the presence of woodlice is correlated with higher predation success by nursery web spiders upon their grasshopper prey. This finding suggested a new hypothesis which links two seemingly disparate food chains: when woodlice are present, the spider predator or the grasshopper prey changes their location in the vegetative canopy in a way that increases their spatial overlap and therefore predation rate. However, warming temperatures may complicate this phenomenon. The spider cannot tolerate thermal stress, meaning warming temperatures may cause the spider to move downwards in the vegetative canopy or otherwise alter its response to woodlice. Therefore, we would expect warming and woodlice presence to have an interactive effect on predation rate.

Methods

We conducted behavioral experiments in 2015, 2017, and 2018 to track habitat domains-the use of the vegetative canopy space by grasshoppers and spiders-in experimental cages. Then, we used three models of spider movement to try to explain the response of spiders to woodlice: expected net energy gain, signal detection theory, and individual-based modelling.

Results

Habitat domain observations revealed that spiders shift upward in the canopy when woodlice are present, but the corresponding effect on grasshopper prey survival was variable over the different years of study. Under warming conditions, spiders remained lower in the canopy regardless of the presence of woodlice, suggesting that thermal stress is more important than the effect of woodlice. Our modelling results suggest that spiders do not need to move away from woodlice to maximize net energy gain (expected net energy gain and signal detection theory models). Instead spider behavior is consistent with the null hypothesis that they move away from unsuccessful encounters with woodlice (individual-based simulation). We conclude that mapping how predator behavior changes across biotic (e.g. woodlouse presence) and abiotic conditions (e.g. temperature) may be critical to anticipate changes in ecosystem dynamics.

The draft genome sequence of the spider Dysdera silvatica (Araneae, Dysderidae): A valuable resource for functional and evolutionary genomic studies in chelicerates

BACKGROUND

We present the draft genome sequence of Dysdera silvatica, a nocturnal ground-dwelling spider from a genus that has undergone a remarkable adaptive radiation in the Canary Islands.

RESULTS

The draft assembly was obtained using short (Illumina) and long (PaciBio and Nanopore) sequencing reads. Our de novo assembly (1.36 Gb), which represents 80% of the genome size estimated by flow cytometry (1.7 Gb), is constituted by a high fraction of interspersed repetitive elements (53.8%). The assembly completeness, using BUSCO and core eukaryotic genes, ranges from 90% to 96%. Functional annotations based on both ab initio and evidence-based information (including D. silvatica RNA sequencing) yielded a total of 48,619 protein-coding sequences, of which 36,398 (74.9%) have the molecular hallmark of known protein domains, or sequence similarity with Swiss-Prot sequences. The D. silvatica assembly is the first representative of the superfamily Dysderoidea, and just the second available genome of Synspermiata, one of the major evolutionary lineages of the "true spiders" (Araneomorphae).

CONCLUSIONS

Dysderoids, which are known for their numerous instances of adaptation to underground environments, include some of the few examples of trophic specialization within spiders and are excellent models for the study of cryptic female choice. This resource will be therefore useful as a starting point to study fundamental evolutionary and functional questions, including the molecular bases of the adaptation to extreme environments and ecological shifts, as well of the origin and evolution of relevant spider traits, such as the venom and silk.