The structure of bacteria-fungi bipartite networks along elevational gradients in contrasting climates

Climate change is altering species distribution and modifying interactions in microbial communities. Understanding microbial community structure and their interactions is crucial to interpreting ecosystem responses to climate change. Here, we examined the assemblages of stream bacteria and fungi, and the associations between the two groups along elevational gradients in two regions with contrasting precipitation and temperature, that is the Galong and Qilian mountains of the Tibetan Plateau. In the wetter and warmer region, the species richness significantly increased and decreased with elevation for bacteria and fungi, respectively, while were nonsignificant in the drier and colder region. Their bipartite network structure was also different by showing significant increases in connectance and nestedness towards higher elevations only in the wetter and warmer region. In addition, these correlation network structure generally exhibited similar positive association with species richness in the wetter and warmer region and the drier and colder region. In the wetter and warmer region, climatic change along elevation was more important in determining connectance and nestedness, whereas microbial species richness exerted a stronger influence on network structure and robustness in the drier and colder region. These findings indicate substantial forthcoming changes in microbial diversity and network structure in warming climates, especially in wetter and warmer regions on Earth, advancing the understanding of microbial bipartite interactions' response to climate change.

Generalism in species interactions is more the consequence than the cause of ecological success

Generalism in resource use is commonly considered a critical driver of population success, species distribution and extinction risk. This idea can be questioned as generalism may be a result rather than the cause of species abundance and range size. We tested these contrasting causal hypotheses focusing on host use in three databases encompassing approximately 44,000 mutualistic (hummingbird-plant), commensalistic (lichen-plant) and parasitic (flea-mammal) interactions in 617 ecological communities across the Americas and Eurasia. Across all interaction types, our analyses indicated that range size and abundance influence the probability of encountering hosts and set the arena for species to express generalism potentials or adapt to new hosts. Hence, our findings support the hypothesis that generalism is a consequence of species ecological success. This highlights the importance of ecological opportunity in driving species characteristics considered key for their survival and conservation.

Structure of compound and component communities of fleas parasitic on small mammals in six different regions as revealed by environmental-based co-occurrence geometry analyses

We inferred the patterns of co-occurrence of flea species in compound (across all host species) and component (across conspecific hosts) communities from six regions of the world (Mongolia, Northwest Argentina, Argentinian Patagonia, West Siberia, Slovakia, and South Africa) using the novel eigenvector ellipsoid method. This method allows us to infer structural community patterns by comparing species' environmental requirements with the pattern of their co-occurrences. We asked whether: (a) communities are characterized by species segregation, nestedness, or modularity; (b) patterns detected by the novel method conform to the patterns identified by traditional methods that search for non-randomness in community structure; and (c) the pattern of flea species co-occurrences in component communities is associated with host species traits. The results of the application of the eigenvector ellipsoid method suggested that the co-occurrence of flea species was random in all compound communities except in South Africa, where this community demonstrated a tendency to be nested. Flea species co-occurrences were random in many component communities. Species segregation was detected in the flea community of one host, whereas the flea communities of 14 hosts from different regions appeared to be nested. No indication of a modular structure in any community was found. The nestedness of flea component communities was mainly characteristic of hosts with a low relative brain mass. We concluded that the application of this novel method that combines data on species distribution and their environmental requirements allows better identification of the community structural patterns and produces more reliable results as compared with traditional methods.

Environment and traits affect parasite and host species positions but not roles in flea-mammal networks

We studied spatial variation in the effects of environment and network size on species positions and roles in multiple flea-mammal networks from four biogeographic realms. We asked whether species positions (measured as species strength [SS], the degree of interaction specialization [d'], and the eigenvector centrality [C]) or the roles of fleas and their hosts in the interaction networks: (a) are repeatable/conserved within a flea or a host species; (b) vary in dependence on environmental variables and/or network size; and (c) the effects of environment and network size on species positions or roles in the networks depend on species traits. The repeatability analysis of species position indices for 441 flea and 429 host species, occurring in at least two networks, demonstrated that the repeatability of SS, d', and C within a species was significant, although not especially high, suggesting that the indices' values were affected by local factors. The majority of flea and host species in the majority of networks demonstrated a peripheral role. A value of at least one index of species position was significantly affected by environmental variables or network size in 41 and 36, respectively, of the 52 flea and 52 host species that occurred in multiple networks. In both fleas and hosts, the occurrence of the significant effect of environment or network size on at least one index of species position, but not on a species' role in a network, was associated with some species traits.

Active degradation-nitrification microbial assemblages in the hypoxic zone in a subtropical estuary

In 2017 summer, we observed widespread bottom hypoxia at the lower estuary of the Pearl River estuary (PRE). Our previous study noticed that AOA and bacteria were highly abundant and clustered within the hypoxia zone. Moreover, nitrification and respiration rates were also evidently higher in these hypoxic waters. These observations prompt us to investigate whether these two oxygen-consuming microorganisms have symbiotic relationships and whether specific groups consistently coexist and form ecological-meaningful associations. In this study, we use network analysis to investigate the presence and active communities (DNA-RNA) based on bacterial and AOA communities sequencing (inferred from the 16S rRNA and amoA gene, respectively) to gain more insight into ecological-meaningful associations. We observed a highly diverse and active bacterial community in the hypoxia zone. The RNA networks were more modulized than the corresponding DNA networks, indicating that the active communities were better parsed into functional microbial assemblages. The network topology revealed that Gammaproteobacteria, Bacteroidetes (Flavobacteriales), Alphaproteobacteria (Rhodobacterales and Rhodospirillales), Marinimicrobia, Cyanobacteria (Synechococcales), and AOA sublineages were module hubs and connectors, indicating that they were the keystone taxa of the microbial communities. The hub-subnetwork further showed robust co-occurrence between Gammaproteobacteria, Bacteroidetes (Flavobacteriales), Alphaproteobacteria (Rhodobacterales and Rhodospirillales), Marinimicrobia with AOA sublineages, and Nitrospinae (presumably NOB) reflecting the formation of Degradation-Nitrification (sequential oxidation of Organic matter degradation to ammonia, then nitrate) microbial assemblage in the hypoxia zone. The subnetworks revealed AOA ecotype-specific modularization and niche partitioning of different AOA sublineages. Interestingly, the recurring co-occurrence of nitrifiers assemblage in the RNA subnetworks (SCM1-like-II (AOA) and Nitrospinae OTUs (NOB) suggests an active interaction via nitrite exchange. The Degradation-Nitrification microbial assemblage may contribute substantially to the oxygen consumption in the hypoxia formation in PRE. Our results provide new insight into the functional microbial assemblages, which is worth further investigation on their ecological implication in estuarine waters.

Parasite-host network analysis provides insights into the evolution of two mistletoe lineages (Loranthaceae and Santalaceae)

Mistletoes are ecologically important parasitic plants, with > 1600 species from five lineages worldwide. Mistletoe lineages exhibit distinct patterns of species diversification and host specificity, however, the mechanisms underlying these differences are poorly understood. In this study, we analysed a comprehensive parasite-host network, including 280 host species from 60 families and 22 mistletoe species from two lineages (Santalaceae and Loranthaceae) in Xishuangbanna, located in a biodiversity hotspot of tropical Asia. We identified the factors that predict the infection strength of mistletoes. We also detected host specificity and the phylogenetic signal of mistletoes and their hosts. We found that this interaction network could be largely explained by a model based on the relative abundance of species. Host infection was positively correlated with diameter at breast height and tree coverage, but negatively correlated with wood density. Overall, closely related mistletoe species tend to interact more often with similar hosts. However, the two lineages showed a significantly different network pattern. Rates of host generality were higher in Loranthaceae than in Santalaceae, although neither lineage showed phylogenetic signal for host generality. This study demonstrates that the neutral interaction hypothesis provides suitable predictions of the mistletoe-host interaction network, and mistletoe species show significant phylogenetic signals for their hosts. Our findings also indicate that high species diversification in Loranthaceae may be explained by high rates of host generality and the evolutionary history shared by Loranthaceae species with diverse host plants in the tropics.

Effect of altitude and spatial heterogeneity on the host-parasite relationship in anurans from a remnant humid forest in the brazilian semiarid

In the present study, we investigated the effect of habitat heterogeneity, elevation gradient, and phylogenetic distance of host species on the abundance and richness of anuran endoparasites, assuming that parasites follow the distribution of their hosts independently of environmental variation. We collected 192 anurans distributed in three altitude ranges: 100-200 m, 400-500 m, and 700-800 m. We performed discriminant principal component analysis to analyze the interrelationships between environmental heterogeneity and the distribution of parasite and host species in the formation of species groups in each altitude range. We estimated the niche width and parasite overlap, using host species as a variable, and assessed whether parasite abundance is more influenced by historical (distance host phylogeny) or ecological effects in each altitude category and overall. Finally, we use network analyses to understand how interactions between parasites and hosts are formed along the altitude gradient. We found 22 parasite species, and the overall prevalence of infection was 74%. In our study, we did not identify environmental (altitude gradients and heterogeneity) or phylogenetic effects acting on the parasite species diversity. Overall, our results suggest that the parasites are distributed following the dispersal of their hosts and are dispersed among most anuran species.

Host specificity of flea parasites of mammals from the Andean Biogeographic Region

Host specificity of fleas affects their biodiversity that plays a major role in determining the potential transmission routes by pathogens through vertebrate hosts, including humans. In the Biogeographic Andean region, numerous systematic and ecological studies have been conducted, revealing a high diversity of flea taxa of mammals and the presence of pathogenic organisms transmitted by fleas; however, the degree of preference with which each flea species associates with a mammal host remains poorly understood in this region. Herein, host specificity in mammal fleas from the Andean region was analysed. We employed the number of host species for each flea species and the index of host specificity STD *. Following the literature, 144 species and 13 subspecies of fleas (31 genera and 10 families) have been described in the Andean biogeographic region; 76 taxa are endemic to this region. To carry out the analyses of host specificity, we considered 1759 records of fleas collected from 124 species and 59 genera of wild and domestic mammals, mostly rodent species (85.9%). Our results indicate that typical Andean fleas are genus or family host specific (mostly STD * less than 3.0). More diverse mammal hosts are parasitized by more diverse flea genera and families and these hosts are phylogenetically related. Otherwise, these hosts are associated with different flea lineages, suggesting the interaction of ecological and evolutionary mechanisms (host-switching, ecological adaptations and co-evolutionary alternation). The fields of disease ecology and One Health are considering the host specificity of arthropod vectors as an important point to understand the mechanisms of emergence and re-emergence of diseases. Our results allow us to estimate the risk of diseases involving fleas in the Andean region.

A network-based approach to identifying correlations between phylogeny, morphological traits and occurrence of fish species in US river basins

The complex network framework has been successfully used to model interactions between entities in Complex Systems in the Biological Sciences such as Proteomics, Genomics, Neuroscience, and Ecology. Networks of organisms at different spatial scales and in different ecosystems have provided insights into community assembly patterns and emergent properties of ecological systems. In the present work, we investigate two questions pertaining to fish species assembly rules in US river basins, a) if morphologically similar fish species also tend to be phylogenetically closer, and b) to what extent are co-occurring species that are phylogenetically close also morphologically similar? For the first question, we construct a network of Hydrologic Unit Code 8 (HUC8) regions as nodes with interaction strengths (edges) governed by the number of common species. For each of the modules of this network, which are found to be geographically separated, there is differential yet significant evidence that phylogenetic distance predicts morphological distance. For the second question, we construct and analyze nearest neighbor directed networks of species based on their morphological distances and phylogenetic distances. Through module detection on these networks and comparing the module-level mean phylogenetic distance and mean morphological distance with the number of basins of common occurrence of species in modules, we find that both phylogeny and morphology of species have significant roles in governing species co-occurrence, i.e. phylogenetically and morphologically distant species tend to co-exist more. In addition, between the two quantities (morphological distance and phylogentic distance), we find that morphological distance is a stronger determinant of species co-occurrences.

Modularity in host-parasite mixed networks: interaction configuration shifts based on human perturbation and parasitism form

Parasitism is an association based on host individual traits and environmental factors. The complexity of this type of interaction is often lost when studying species-by-species interaction networks. Here we analyze changes in modularity - a metric describing groups of nodes interacting much more frequently among themselves than they do with nodes of other modules, considering the host individual variation and the different forms of parasitism: ecto- and endo-parasitism. For this, we studied mixed networks: bipartite networks comprising host individuals and parasite species as two sets of nodes interacting with each other. We used a fish-parasite mixed network from a highly perturbed coastal river to understand how an anthropogenic perturbation gradient influences the modular structure of host-parasite networks. In addition, we tested how host individual traits drove module configuration within host-parasite mixed networks. Our results showed that different forms of parasitism respond differently to the environment: modularity in fish-ectoparasite networks increased with human perturbation, but modularity was not related to human perturbation in fish-endoparasite networks. In addition, mixed network modules were intrinsically related to individual variation, with host intensity of infection being the most important trait, regardless of the parasite's life form. The effect of total abundance over network structure indicates signs of changes in community equilibrium, with an increase in species with opportunistic behaviors. Module composition was also related to host fitness and body size, which were most predictive in more preserved and diverse river sections. Overall, our results indicate that host-parasite networks are sensitive to ecological gradients marked by human perturbation and that host individual fitness helps to determine network structure.

Bat-parasite interaction networks in urban green areas in northeastern Brazil

Interaction networks can provide detailed information regarding ecological systems, helping us understand how communities are organized and species are connected. The goals of this study were to identify the pattern of interaction between bats and ectoparasites in urban green areas of Grande Aracaju, Sergipe, and calculate connectance, specialization, nesting, modularity and centrality metrics. Bats were captured using 10 mist nets inside and on the edges of the fragments, and the collected ectoparasites were stored in 70% alcohol. All analyses were performed using R software. The interaction network consisted of 10 species of bats and 13 ectoparasites. Connectivity was considered low (0.12). The specialization indices for ectoparasites ranged from 0.50 to 1.00, and the value obtained for the network was 0.96, which is high. The observed nesting metric was low (wNODF = 1.47), whereas the modularity was high (wQ = 0.74), indicating that the studied network had a modular topology. All centrality metrics had low values. The observed modularity may have been caused by the evolutionary history of the bats and ectoparasites involved and the high specificity index of the interactions. The low centrality values may be associated with low connectivity and a high degree of specialization. This study provides relevant information on bat–parasite interactions in an urban environment, highlighting the need for further studies to improve our understanding of host–parasite interaction networks.

Sensitivity of bipartite network analyses to incomplete sampling and taxonomic uncertainty

Bipartite network analysis is a powerful tool to study the processes structuring interactions in ecological communities. In applying the method, it is assumed that the sampled interactions provide an accurate representation of the actual community. However, acquiring a representative sample may be difficult as not all species are equally abundant or easily identifiable. Two potential sampling issues can compromise the conclusions of bipartite network analyses: failure to capture the full range of interactions (sampling completeness) and use of a taxonomic level higher than species to evaluate the network (taxonomic resolution). We asked how commonly used descriptors of bipartite antagonistic communities (modularity, nestedness, connectance, and specialization [H₂ ']) are affected by reduced host sampling completeness, parasite taxonomic resolution, and their crossed effect, as they are likely to co-occur. We used a quantitative niche model to generate weighted bipartite networks that resembled natural host-parasite communities. The descriptors were more sensitive to uncertainty in parasite taxonomic resolution than to host sampling completeness. When only 10% of parasite taxonomic resolution was retained, modularity and specialization decreased by ~76% and ~12%, respectively, and nestedness and connectance increased by ~114% and ~345% respectively. The loss of taxonomic resolution led to a wide range of possible communities, which made it difficult to predict its effects on a given network. With regards to host sampling completeness, standardized nestedness, connectance, and specialization were robust, whereas modularity was sensitive (~30% decrease). The combination of both sampling issues had an additive effect on modularity. In communities with low effort for both sampling issues (50%-10% of sampling completeness and taxonomic resolution), estimators of modularity, and nestedness could not be distinguished from those of random assemblages. Thus, the categorical description of communities with low sampling effort (e.g., if a community is modular or not) should be done with caution. We recommend evaluating both sampling completeness and taxonomic certainty when conducting bipartite network analyses. Care should also be exercised when using nonrobust descriptors (the four descriptors for parasite taxonomic resolution; modularity for host sampling completeness) when sampling issues are likely to affect a dataset.

A 25-Year Retrospective Review of Mortality in Chimpanzees (Pan troglodytes) in Accredited U.S. Zoos from a Management and Welfare Perspective

Understanding causes of death allows adjustment of health management strategies for animals in managed care. From 224 documented chimpanzee deaths occurring from 1995 to 2019 in 42 accredited U.S. zoos, post-mortem records and necropsy reports were analyzed for the primary cause of death, which were available for 214 individuals. In total, 37 cases of stillbirth and neonatal deaths were assessed (16.5%); however, the focus was otherwise placed on the remaining 177 cases in which the death occurred in individuals aged greater than 1 month. There were no sex-related differences in etiology; however, age variation in the cause of death was statistically significant (p < 0.001). Elderly (35 years and older) chimpanzees tended to die of intrinsic, often degenerative, etiologies, whereas infants, juveniles, and adolescents (less than 15 years) were more likely to be involved in fatal trauma. Overall, there were 27 deaths (15.3% of all post-neonatal deaths) related to trauma and 13 of these were directly or indirectly related to conspecific aggression. Understanding causes of mortality and the interrelation with management can benefit managed populations of chimpanzees.

Phylogenetic structure of specialization: A new approach that integrates partner availability and phylogenetic diversity to quantify biotic specialization in ecological networks

Biotic specialization holds information about the assembly, evolution, and stability of biological communities. Partner availabilities can play an important role in enabling species interactions, where uneven partner availabilities can bias estimates of biotic specialization when using phylogenetic diversity indices. It is therefore important to account for partner availability when characterizing biotic specialization using phylogenies. We developed an index, phylogenetic structure of specialization (PSS), that avoids bias from uneven partner availabilities by uncoupling the null models for interaction frequency and phylogenetic distance. We incorporate the deviation between observed and random interaction frequencies as weights into the calculation of partner phylogenetic α‐diversity. To calculate the PSS index, we then compare observed partner phylogenetic α‐diversity to a null distribution generated by randomizing phylogenetic distances among the same number of partners. PSS quantifies the phylogenetic structure (i.e., clustered, overdispersed, or random) of the partners of a focal species. We show with simulations that the PSS index is not correlated with network properties, which allows comparisons across multiple systems. We also implemented PSS on empirical networks of host–parasite, avian seed‐dispersal, lichenized fungi–cyanobacteria, and hummingbird pollination interactions. Across these systems, a large proportion of taxa interact with phylogenetically random partners according to PSS, sometimes to a larger extent than detected with an existing method that does not account for partner availability. We also found that many taxa interact with phylogenetically clustered partners, while taxa with overdispersed partners were rare. We argue that species with phylogenetically overdispersed partners have often been misinterpreted as generalists when they should be considered specialists. Our results highlight the important role of randomness in shaping interaction networks, even in highly intimate symbioses, and provide a much‐needed quantitative framework to assess the role that evolutionary history and symbiotic specialization play in shaping patterns of biodiversity. PSS is available as an R package at https://github.com/cjpardodelahoz/pss.

Scenes From Tick Physiology: Proteins of Sialome Talk About Their Biological Processes

Ticks are blood-sucking parasites with different strategies of feeding depending on the tick family. The major families are Ixodidae or Argasidae, being slow or fast feeders, respectively. In the recent years, the advances in molecular sequencing techniques have enabled to gain knowledge about the proteome of the tick's salivary glands. But an holistic view of the biological processes underlying the expression of the sialome has been neglected. In this study we propose the use of standard biological processes as a tool to draw the physiology of the tick's salivary glands. We used published data on the sialome of Rhipicephalus sanguineus s.l. (Ixodidae) and Ornithodoros rostratus (Argasidae). A partial set of proteins obtained by these studies were used to define the biological process(es) in which proteins are involved. We used a directed network construction in which the nodes are proteins (source) and biological processes (target), separately for the low-level processes ("children") and the top-level ones ("parents"). We applied the method to feeding R. sanguineus at different time slices, and to different organs of O. rostratus. The network connects the proteins and the processes with a strength directly proportional to the transcript per millions of each protein. We used PageRank as a measure of the importance of each biological process. As suggested in previous studies, the sialome of unfed R. sanguineus express about 30% less biological processes than feeding ticks. Another decrease (25%) is noticed at the middle of the feeding and before detachment. However, top-level processes are deeply affected only at the onset of feeding, demonstrating a redundancy in the feeding. When ixodid-argasid are compared, large differences were observed: they do not share 91% of proteins, but share 90% of the biological processes. However, caution must be observed when examining these results. The hypothesis of different proteins linked to similar biological process(es) in both ticks is an extreme not confirmed in this study. Considering the limitations of this study, carried out with a selected set of proteins, we propose the networks of proteins of sialome linked to their biological processes as a tool aimed to explain the biological processes behind families of proteins.

The structure of tropical bat-plant interaction networks during an extreme El Niño-Southern Oscillation event

Interaction network structure reflects the ecological mechanisms acting within biological communities, which are affected by environmental conditions. In tropical forests, higher precipitation usually increases fruit production, which may lead frugivores to increase specialization, resulting in more modular and less nested animal–plant networks. In these ecosystems, El Niño is a major driver of precipitation, but we still lack knowledge of how species interactions change under this influence. To understand bat–plant network structure during an extreme El Niño‐Southern Oscillation event, we determined the links between plantivorous bat species and the plants they consume by DNA barcoding seeds and pulp in bat faeces. These interactions were recorded in the dry forest and rainforest of Costa Rica, during the dry and the wet seasons of an extreme El Niño year. From these we constructed seasonal and whole‐year bat–plant networks and analysed their structures and dissimilarities. In general, networks had low nestedness, had high modularity, and were dominated by one large compartment which included most species and interactions. Contrary to our expectations, networks were less nested and more modular in drier conditions, both in the comparison between forest types and between seasons. We suggest that increased competition, when resources are scarce during drier seasons and habitats, lead to higher resource partitioning among bats and thus higher modularity. Moreover, we have found similar network structures between dry and rainforests during El Niño and non‐El Niño years. Finally, most interaction dissimilarity among networks occurred due to interaction rewiring among species, potentially driven by seasonal changes in resource availability.

Endoparasite community structure of an anuran assemblage in the Caatinga, Northeastern Neotropical Region

Amphibians are a widespread Chordata taxon and are important for maintaining the balance of both terrestrial and aquatic ecosystems. Brazil has a rich amphibian fauna; however, little is known about the role of their ecology and phylogenetic relationships during the assembly processes of associated endoparasite communities. Herein, we describe an endoparasite community in an anuran assemblage in the Caatinga, a unique biome of dry forests in north-eastern Brazil. We studied endoparasite diversity, as well as the effects of body length, body mass, body volume and sex on parasite abundance. We also investigated the influence of ecological and historical factors and anuran microhabitat use on endoparasite composition. We analysed individuals from 13 anuran species distributed across five families: Odontophrynidae (Proceratophrys cristiceps); Leptodactylidae (Leptodactylus fuscus, Leptodactylus vastus, Leptodactylus macrosternum, Leptodactylus troglodytes and Physalaemus cuvieri); Hylidae (Pithecopus gonzagai, Scinax x-signatus, Boana raniceps and Dendropsophus nanus); Bufonidae (Rhinella diptycha and Rhinella granulosa); and Microhylidae (Dermatonotus muelleri). We found nine species of endoparasites, including seven nematodes (Aplectana membranosa, Cosmocerca sp., Oswaldocruzia mazzai, Raillietnema spectans, Rhabdias fuelleborni, Schrankiana sp. and Physaloptera sp.), one species of Trematoda (Glypthelmins pseudium) and one non-identified cestode. There was no significant relationship between endoparasite abundance and host body length, body mass, body volume and sex. A phylogenetic principal component analysis showed that ecological factors had a greater influence on endoparasite assemblage than historical factors. Similarly, our results showed that ecological factors had a greater influence on anuran microhabitat use compared to historical factors, which contributed to the generalist characteristics presented by most of the sampled endoparasite species.

Identifying co-phylogenetic hotspots for zoonotic disease

The incidence of zoonotic diseases is increasing worldwide, which makes identifying parasites likely to become zoonotic and hosts likely to harbour zoonotic parasites a critical concern. Prior work indicates that there is a higher risk of zoonotic spillover accruing from closely related hosts and from hosts that are infected with a high phylogenetic diversity of parasites. This suggests that host and parasite evolutionary history may be important drivers of spillover, but identifying whether host-parasite associations are more strongly structured by the host, parasite or both requires co-phylogenetic analyses that combine host-parasite association data with host and parasite phylogenies. Here, we use host-parasite datasets containing associations between helminth taxa and free-range mammals in combination with phylogenetic models to explore whether host, parasite, or both host and parasite evolutionary history influences host-parasite associations. We find that host phylogenetic history is most important for driving patterns of helminth-mammal association, indicating that zoonoses are most likely to come from a host's close relatives. More broadly, our results suggest that co-phylogenetic analyses across broad taxonomic scales can provide a novel perspective for surveying potential emerging infectious diseases. This article is part of the theme issue 'Infectious disease macroecology: parasite diversity and dynamics across the globe'.

Modularity and specialization in bat-fly interaction networks are remarkably consistent across patches within urbanized landscapes and spatial scales

Patterns of specialization and the structure of interactions between bats and ectoparasitic flies have been studied mostly on non-urban environments and at local scales. Thus, how anthropogenic disturbances influence species interactions and network structure in this system remain poorly understood. Here, we investigated patterns of interaction between Phyllostomidae bats and ectoparasitic Streblidae flies, and variations in network specialization and structure across Cerrado patches within urbanized landscapes in Brazil and between local and regional scales. We found high similarity in the richness and composition of bat and fly species across communities, associated with low turnover of interactions between networks. The high specialization of bat–streblid interactions resulted in little connected and modular networks, with the emergence of modules containing subsets of species that interact exclusively or primarily with each other. Such similarities in species and interaction composition and network structure across communities and scales suggest that bat–fly interactions within Cerrado patches are little affected by the degree of human modification in the surrounding matrix. This remarkable consistency is likely promoted by specific behaviors, the tolerance of Phyllostomidae bats to surrounding urbanized landscapes as well as by the specificity of the streblid–bat interactions shaped over evolutionary time.

Host selection pattern and flavivirus screening of mosquitoes in a disturbed Colombian rainforest

Studies on the feeding behavior of hematophagous insects, particularly those of medical importance, are relevant for tracking possible pathogen transmission routes and identifying biases in the choice of vertebrates. We evaluated host selection of blood-feeding mosquitoes in a disturbed forest in the Magdalena Medio valley in Colombia from March 2017 to April 2018, after the introduction of Zika virus to the Americas from the 2015-2016 outbreak. We estimated vertebrate diversity and collected blood-engorged female mosquitoes. Genomic DNA/RNA was extracted from the mosquito's abdomen for vertebrate host identification and pathogen detection. We performed conventional PCR and sequencing, using universal primers targeting vertebrate regions of the eukaryotic mitochondrial genome to determine bloodmeal host. Additionally, we tested for the presence of flaviviruses in all mosquito samples with RT-PCR. Based on the identity and quantity of detected bloodmeals, we performed mosquito-vertebrate interaction network analysis and estimated topology metrics. In total, we collected 292 engorged female mosquitoes representing 20 different species. Bloodmeal analyses identified 26 vertebrate species, the majority of which were mammals (N = 16; 61.5%). No flaviviruses of medical importance were detected from the samples. Although feeding patterns varied, network analyses showed a high degree of specialization by mosquitoes and revealed ecological and phylogenetic relationships among the host community. We conclude that host selection or preference by mosquitoes is species specific.

The roles of anurans in antagonistic networks are explained by life-habit and body-size

Interactions among living beings are the structuring basis of ecosystems, and studies of networks allow us to identify the patterns and consistency of such interactions. Antagonistic networks reflect the energy flow of communities, and identifying network structure and the biological aspects that influence its stability is crucial to understanding ecosystem functioning. We used antagonistic anuran interactions-predator-prey and host-parasite-to assess structural patterns and to identify the key anuran species structuring these networks. We tested whether anuran body-size and life-habit are related to their roles in these networks. We collected individuals of 9 species of anurans from an area of the Atlantic Forest in Brazil and identified their prey and helminth parasites. We used network (modularity, specialization, and nestedness) and centrality metrics (degree, closeness, and betweenness) to identify the role of anuran species in both networks. We then evaluated whether anuran body-size or life-habit were related to anuran centrality using generalized linear mixed models. The networks formed specialized interactions in compartments composed by key species from different habits. In our networks, anurans with rheophilic and cryptozoic habit are central in predator-prey networks, and those with larger body size and arboreal and cryptozoic habit in the host-parasite network. This study represents a step towards a better understanding of the influential factors that affect the structure of anuran antagonist networks, as well as to recognize the functioning roles of anuran species.

Haemosporidian taxonomic composition, network centrality and partner fidelity between resident and migratory avian hosts

Migration can modify interaction dynamics between parasites and their hosts with migrant hosts able to disperse parasites and impact local community transmission. Thus, studying the relationships among migratory hosts and their parasites is fundamental to elucidate how migration shapes host-parasite interactions. Avian haemosporidians are some of the most prevalent and diverse group of wildlife parasites and are also widely studied as models in ecological and evolutionary research. Here, we contrast partner fidelity, network centrality and parasite taxonomic composition among resident and non-resident avian hosts using presence/absence data on haemosporidians parasitic in South American birds as study model. We ran multilevel Bayesian models to assess the role of migration in determining partner fidelity (i.e., normalized degree) and centrality (i.e., weighted closeness) in host-parasite networks of avian hosts and their respective haemosporidian parasites. In addition, to evaluate parasite taxonomic composition, we performed permutational multivariate analyses of variance to quantify dissimilarity in haemosporidian lineages infecting different host migratory categories. We observed similar partner fidelity and parasite taxonomic composition among resident and migratory hosts. Conversely, we demonstrate that migratory hosts play a more central role in host-parasite networks than residents. However, when evaluating partially and fully migratory hosts separately, we observed that only partially migratory species presented higher network centrality when compared to resident birds. Therefore, migration does not lead to differences in both partner fidelity and parasite taxonomic composition. However, migratory behavior is positively associated with network centrality, indicating migratory hosts play more important roles in shaping host-parasite interactions and influence local transmission.

Disentangling the beta-diversity in anuran parasite communities

There is great heterogeneity in parasite communities among hosts, understanding the nature and drivers of such variations is still a great scientific quest. Here, we analyse the variation in parasite communities by addressing the following questions: (i) What is the beta-diversity component (nestedness or turnover) that most contributes to beta diversity in parasite communities among anuran species? (ii) Does the beta diversity of parasite communities follow a non-random pattern? (iii) Is the dissimilarity in composition of parasite communities related to the phylogenetic or functional dissimilarity among hosts? We found that turnover in parasite assemblages was the main component of beta diversity, but the variation observed both in the total beta diversity and in its components did not differ from the respective null models. The dissimilarity among parasite communities was not related to the phylogenetic species variability or functional dissimilarity among anuran species for most localities. In short, our findings may indicate a process of resource tracking by the parasite species, in which the resource may not necessarily be conserved phylogenetically in their hosts.

Network Analysis: Ten Years Shining Light on Host-Parasite Interactions

Biological interactions are key drivers of ecological and evolutionary processes. The complexity of such interactions hinders our understanding of ecological systems and our ability to make effective predictions in changing environments. However, network analysis allows us to better tackle the complexity of ecosystems because it extracts the properties of an ecological system according to the number and distribution of links among interacting entities. The number of studies using network analysis to solve ecological and evolutionary questions in parasitology has increased over the past decade. Here, we synthesise the contribution of network analysis toward disentangling host-parasite processes. Furthermore, we identify current trends in mainstream ecology and novel applications of network analysis that present opportunities for research on host-parasite interactions.

Nematodes of Proceratophrys ararype (Anura: Odontophrynidae), an endemic frog from the Araripe Plateau, northeastern Brazil

Abstract: Parasites are an important component of the global biomass, having significant roles in several regulatory mechanisms in the ecosystem. Parasitism is one of the most common ecological interactions on the planet. Studies have shown that the helminth fauna of only 8% of amphibian species in Brazil have been studied, and this percentage is lower for the Odontophrynidae family, with only four of the 50 species kown to occur in Brazil having been investigated. Here, we present the helminth fauna of Proceratophrys ararype, an anuran endemic to the “Brejo de Altitude” Chapada do Araripe (Araripe Plateau), northeastern Brazil. The infection parameters analyzed were prevalence, mean intensity of infection and mean abundance of parasites. We used the Pearson’s linear correlation coefficient to check the correlations between the abundance of the parasites with the snout-vent length (SVL) of hosts. To verify the degree of aggregation of parasites in hosts, we used the dispersion index. Out of 40 specimens examined, 19 specimens were infected with at least one parasite. The overall prevalence was 47.5% with a mean infection intensity of 18.93 ± 10.77. The endoparasitic community associated with P. ararype consisted of six species of parasites, with Falcaustra mascula having the highest prevalence (25%). Most parasites had a uniform dispersion index in the hosts and their abundance was unrelated to host size. Here, we also present a compilation of all parasites associated with host species of the genus Proceratophrys from South America. Overall, we found 23 species of parasites associated with five host species (P. ararype, P. appendiculata, P. boiei, P. cristiceps, P. mantiqueira). Of these, two species of parasites (Oswaldocruzia mazzai and Strongyloides sp.) represent new records for the genus Proceratophrys. Our results demonstrate the lack of studies on amphibian helminth fauna and fill an important knowledge gap on the diversity of parasites of Proceratophrys ararype, an endemic frog from the Araripe Plateau, northeastern Brazil.

Nonadaptive host-use specificity in tropical armored scale insects

Most herbivorous insects are diet specialists in spite of the apparent advantages of being a generalist. This conundrum might be explained by fitness trade-offs on alternative host plants, yet the evidence of such trade-offs has been elusive. Another hypothesis is that specialization is nonadaptive, evolving through neutral population-genetic processes and within the bounds of historical constraints. Here, we report on a striking lack of evidence for the adaptiveness of specificity in tropical canopy communities of armored scale insects. We find evidence of pervasive diet specialization, and find that host use is phylogenetically conservative, but also find that more-specialized species occur on fewer of their potential hosts than do less-specialized species, and are no more abundant where they do occur. Of course local communities might not reflect regional diversity patterns. But based on our samples, comprising hundreds of species of hosts and armored scale insects at two widely separated sites, more-specialized species do not appear to outperform more generalist species.

The Structure of Ecological Networks Across Levels of Organization

Interactions connect the units of ecological systems, forming networks. Individual-based networks characterize variation in niches among individuals within populations. These individual-based networks merge with each other, forming species-based networks and food webs that describe the architecture of ecological communities. Networks at broader spatiotemporal scales portray the structure of ecological interactions across landscapes and over macroevolutionary time. Here, I review the patterns observed in ecological networks across multiple levels of biological organization. A fundamental challenge is to understand the amount of interdependence as we move from individual-based networks to species-based networks and beyond. Despite the uneven distribution of studies, regularities in network structure emerge across scales due to the fundamental architectural patterns shared by complex networks and the interplay between traits and numerical effects. I illustrate the integration of these organizational scales by exploring the consequences of the emergence of highly connected species for network structures across scales.

The network structure and eco-evolutionary dynamics of CRISPR-induced immune diversification

As a heritable sequence-specific adaptive immune system, CRISPR-Cas is a powerful force shaping strain diversity in host-virus systems. While the diversity of CRISPR alleles has been explored, the associated structure and dynamics of host-virus interactions have not. We explore the role of CRISPR in mediating the interplay between host-virus interaction structure and eco-evolutionary dynamics in a computational model and compare the results with three empirical datasets from natural systems. We show that the structure of the networks describing who infects whom and the degree to which strains are immune, are respectively modular (containing groups of hosts and viruses that interact strongly) and weighted-nested (specialist hosts are more susceptible to subsets of viruses that in turn also infect the more generalist hosts with many spacers matching many viruses). The dynamic interplay between these networks influences transitions between dynamical regimes of virus diversification and host control. The three empirical systems exhibit weighted-nested immunity networks, a pattern our theory shows is indicative of hosts able to suppress virus diversification. Previously missing from studies of microbial host-pathogen systems, the immunity network plays a key role in the coevolutionary dynamics.

High fidelity defines the temporal consistency of host-parasite interactions in a tropical coastal ecosystem

Host-parasite interactions represent a selective force that may reduce hosts’ lifespan, their reproductive success and survival. Environmental conditions can affect host-parasite communities, leading to distinct patterns of interactions with divergent ecological and evolutionary consequences for their persistence. Here, we tested whether climatic oscillation shapes the temporal dynamics of bird-haemosporidian associations, assessing the main mechanisms involved in the temporal dissimilarity of their interactions’ networks. For two years, we monthly sampled birds in a tropical coastal ecosystem to avian malaria molecular diagnosis. The studied networks exhibited high specialization, medium modularity, with low niche overlap among parasites lineages. Moreover, alpha and β-diversity of hosts, parasites and their interactions, as well as the structure of their networks were temporally consistent, i.e., stable under fluctuations in temperature or precipitation over seasons. The structure and temporal consistency of the studied antagonistic networks suggest a high fidelity between partners, which is likely relevant for their evolutionary persistence.

An individual-based model for the eco-evolutionary emergence of bipartite interaction networks

How ecological interaction networks emerge on evolutionary time scales remains unclear. Here we build an individual-based eco-evolutionary model for the emergence of mutualistic, antagonistic and neutral bipartite interaction networks. Exploring networks evolved under these scenarios, we find three main results. First, antagonistic interactions tend to foster species and trait diversity, while mutualistic interactions reduce diversity. Second, antagonistic interactors evolve higher specialisation, which results in networks that are often more modular than neutral ones; resource species in these networks often display phylogenetic conservatism in interaction partners. Third, mutualistic interactions lead to networks that are more nested than neutral ones, with low phylogenetic conservatism in interaction partners. These results tend to match overall empirical trends, demonstrating that structures of empirical networks that have most often been explained by ecological processes can result from an evolutionary emergence. Our model contributes to the ongoing effort of better integrating ecological interactions and macroevolution.

The role of ecological opportunity in shaping host–parasite networks

Despite the great interest to quantify the structure of host–parasite interaction networks, the real influence of some factors such as taxonomy, host body size and ecological opportunity remains poorly understood. In this paper, we investigate the general patterns of organization and structure of interactions in two anuran–parasite networks in the Brazilian Pantanal (seasonally flooded environment) and Atlantic Forest (non-flooded forest). We present theoretical models to test whether the structures of these host–parasite interaction networks are influenced by neutrality, host taxonomy and host body size. Subsequently, we calculated metrics of connectance, nestedness and modularity to characterize the network structure. We demonstrated the structure networks were influenced mainly by body size and taxonomy of the host. Moreover, our results showed that the seasonally flooded environment present networks with higher connectance/nestedness and lower modularity compared to the other environment. The results also suggest that seasonal floods may promote ecological opportunities for new species associations.

Partner Fidelity and Asymmetric Specialization in Ecological Networks

AbstractSpecies are embedded in complex networks of interdependencies that may change across geographic locations. Yet most approaches to investigate the architecture of this entangled web of life have considered exclusively local communities. To quantify to what extent species interactions change at a biogeographic scale, we need to shed light on how among-community variation affects the occurrence of species interactions. Here we quantify the probability for two partners to interact wherever they co-occur (i.e., partner fidelity) by analyzing the most extensive database on species interaction networks worldwide. We found that mutualistic species show more fidelity in their interactions than antagonistic ones when there is asymmetric specialization (i.e., when specialist species interact with generalist partners). Moreover, resources (e.g., plants in plant-pollinator mutualisms or hosts in host-parasite interactions) show a higher partner fidelity in mutualistic interactions than in antagonistic interactions, which can be explained neither by sampling effort nor by phylogenetic constraints developed during their evolutionary histories. In spite of the general belief that mutualistic interactions among free-living species are labile, asymmetric specialization is very much conserved across large geographic areas.

Generalists are more specialized in low-resource habitats, increasing stability of ecological network structure

Linking mechanistic processes to the stability of ecological networks is a key frontier in ecology. In trophic networks, "modules"-groups of species that interact more with each other than with other members of the community-confer stability, mitigating effects of species loss or perturbation. Modularity, in turn, is shaped by the interplay between species' diet breadth traits and environmental influences, which together dictate interaction structure. Despite the importance of network modularity, variation in this emergent property is poorly understood in complex natural systems. Using two years of field data, we quantified interactions between a rich community of lepidopteran herbivores and their host plants across a mosaic of low-resource serpentine and high-resource nonserpentine soils. We used literature and our own observations to categorize herbivore species as generalists (feeding on more than one plant family) or specialists (feeding on one plant family). In both years, the plant-herbivore network was more modular on serpentine than on nonserpentine soils-despite large differences in herbivore assemblage size across years. This structural outcome was primarily driven by reduction in the breadth of host plant use by generalist species, rather than by changes in the composition of species with different fundamental diet breadths. Greater modularity-and thus greater stability-reflects environmental conditions and plastic responses by generalist herbivores to low host plant quality. By considering the dual roles of species traits and ecological processes, we provide a deeper mechanistic understanding of network modularity, and suggest a role for resource availability in shaping network persistence.

A new model explaining the origin of different topologies in interaction networks

Nestedness and modularity have been recurrently observed in species interaction networks. Some studies argue that those topologies result from selection against unstable networks, and others propose that they likely emerge from processes driving the interactions between pairs of species. Here we present a model that simulates the evolution of consumer species using resource species following simple rules derived from the integrative hypothesis of specialization (IHS). Without any selection on stability, our model reproduced all commonly observed network topologies. Our simulations demonstrate that resource heterogeneity drives network topology. On the one hand, systems containing only homogeneous resources form generalized nested networks, in which generalist consumers have higher performance on each resource than specialists. On the other hand, heterogeneous systems tend to have a compound topology: modular with internally nested modules, in which generalists that divide their interactions between modules have low performance. Our results demonstrate that all real-world topologies likely emerge through processes driving interactions between pairs of species. Additionally, our simulations suggest that networks containing similar species differ from heterogeneous networks and that modules may not present the topology of entire networks.