Disruption of Metapopulation Structure Reduces Tasmanian Devil Facial Tumour Disease Spread at the Expense of Abundance and Genetic Diversity

Metapopulation structure plays a fundamental role in the persistence of wildlife populations. It can also drive the spread of infectious diseases and transmissible cancers such as the Tasmanian devil facial tumour disease (DFTD). While disrupting this structure can reduce disease spread, it can also impair host resilience by disrupting gene flow and colonisation dynamics. Using an individual-based metapopulation model we investigated the synergistic effects of host dispersal, disease transmission rate and inter-individual contact distance for transmission, on the spread and persistence of DFTD from local to regional scales. Disease spread, and the ensuing population declines, are synergistically determined by individuals’ dispersal, disease transmission rate and within-population mixing. Transmission rates can be magnified by high dispersal and inter-individual transmission distance. The isolation of local populations effectively reduced metapopulation-level disease prevalence but caused severe declines in metapopulation size and genetic diversity. The relative position of managed (i.e., isolated) local populations had a significant effect on disease prevalence, highlighting the importance of considering metapopulation structure when implementing metapopulation-scale disease control measures. Our findings suggest that population isolation is not an ideal management method for preventing disease spread in species inhabiting already fragmented landscapes, where genetic diversity and extinction risk are already a concern.

Disruption of metapopulation structure effectively reduces Tasmanian devil facial tumour disease spread at the expense of abundance and genetic diversity

Metapopulation structure (i.e. the spatial arrangement of local populations and corridors between them) plays a fundamental role in the persistence of wildlife populations, but can also drive the spread of infectious diseases. While the disruption of metapopulation connectivity can reduce disease spread, it can also impair host resilience by disrupting gene flow and colonisation dynamics. Thus, a pressing challenge for many wildlife populations is to elucidate whether the benefits of disease management methods that reduce metapopulation connectivity outweigh the associated risks. Directly transmissible cancers are clonal malignant cell lines capable to spread through host populations without immune recognition, when susceptible and infected hosts become in close contact. Using an individual-based metapopulation model we investigate the effects of the interplay between host dispersal, disease transmission rate and inter-individual contact distance for transmission (determining within-population mixing) on the spread and persistence of a transmissible cancer, Tasmanian devil facial tumour disease (DFTD), from local to regional scales. Further, we explore population isolation scenarios to devise management strategies to mitigate disease spread. Disease spread, and the ensuing population declines, are synergistically determined by individuals' dispersal, disease transmission rate and within-population mixing. Low to intermediate transmission rates can be magnified by high dispersal and inter-individual transmission distance. Once disease transmission rate is high, dispersal and inter-individual contact distance do not impact the outcome of the disease transmission dynamics. Isolation of local populations effectively reduced metapopulation-level disease prevalence but caused severe declines in metapopulation size and genetic diversity. The relative position of managed (i.e. isolated) populations within the metapopulation had a significant effect on disease prevalence, highlighting the importance of considering metapopulation structure when implementing metapopulation-scale disease control measures. Our findings suggests that population isolation is not an ideal management method for preventing disease spread in species inhabiting already fragmented landscapes, where genetic diversity and extinction risk are already a concern, such as the Tasmanian devil.

Dongsha Atoll is an important stepping-stone that promotes regional genetic connectivity in the South China Sea

Background Understanding region-wide patterns of larval connectivity and gene flow is crucial for managing and conserving marine biodiversity. Dongsha Atoll National Park (DANP), located in the northern South China Sea (SCS), was established in 2007 to study and conserve this diverse and remote coral atoll. However, the role of Dongsha Atoll in connectivity throughout the SCS is seldom studied. In this study, we aim to evaluate the role of DANP in conserving regional marine biodiversity. Methods In total, 206 samples across nine marine species were collected and sequenced from Dongsha Atoll, and these data were combined with available sequence data from each of these nine species archived in the Genomic Observatories Metadatabase (GEOME). Together, these data provide the most extensive population genetic analysis of a single marine protected area. We evaluate metapopulation structure for each species by using a coalescent sampler, selecting among panmixia, stepping-stone, and island models of connectivity in a likelihood-based framework. We then completed a heuristic graph theoretical analysis based on maximum dispersal distance to get a sense of Dongsha’s centrality within the SCS. Results Our dataset yielded 111 unique haplotypes across all taxa at DANP, 58% of which were not sampled elsewhere. Analysis of metapopulation structure showed that five out of nine species have strong regional connectivity across the SCS such that their gene pools are effectively panmictic (mean pelagic larval duration (PLD) = 78 days, sd = 60 days); while four species have stepping-stone metapopulation structure, indicating that larvae are exchanged primarily between nearby populations (mean PLD = 37 days, sd = 15 days). For all but one species, Dongsha was ranked within the top 15 out of 115 large reefs in the South China Sea for betweenness centrality. Thus, for most species, Dongsha Atoll provides an essential link for maintaining stepping-stone gene flow across the SCS. Conclusions This multispecies study provides the most comprehensive examination of the role of Dongsha Atoll in marine connectivity in the South China Sea to date. Combining new and existing population genetic data for nine coral reef species in the region with a graph theoretical analysis, this study provides evidence that Dongsha Atoll is an important hub for sustaining connectivity for the majority of coral-reef species in the region.

Impact of cyclones on hard coral and metapopulation structure, connectivity and genetic diversity of coral reef fish

Cyclones have one of the greatest effects on the biodiversity of coral reefs and the associated species. But it is unknown how stochastic alterations in habitat structure influence metapopulation structure, connectivity and genetic diversity. From 1993 to 2018, the reefs of the Capricorn Bunker Reef group in the southern part of the Great Barrier Reef were impacted by three tropical cyclones including cyclone Hamish (2009, category 5). This resulted in substantial loss of live habitat-forming coral and coral reef fish communities. Within 6–8 years after cyclones had devastated, live hard corals recovered by 50–60%. We show the relationship between hard coral cover and the abundance of the neon damselfish (Pomacentrus coelestis), the first fish colonizing destroyed reefs. We present the first long-term (2008–2015 years corresponding to 16–24 generations of P. coelestis) population genetic study to understand the impact of cyclones on the meta-population structure, connectivity and genetic diversity of the neon damselfish. After the cyclone, we observed the largest change in the genetic structure at reef populations compared to other years. Simultaneously, allelic richness of genetic microsatellite markers dropped indicating a great loss of genetic diversity, which increased again in subsequent years. Over years, metapopulation dynamics were characterized by high connectivity among fish populations associated with the Capricorn Bunker reefs (2200 km2); however, despite high exchange, genetic patchiness was observed with annual strong genetic divergence between populations among reefs. Some broad similarities in the genetic structure in 2015 could be explained by dispersal from a source reef and the related expansion of local populations. This study has shown that alternating cyclone-driven changes and subsequent recovery phases of coral habitat can greatly influence patterns of reef fish connectivity. The frequency of disturbances determines abundance of fish and genetic diversity within species.

Metapopulation Patterns of Iberian Butterflies Revealed by Fuzzy Logic

Metapopulation theory considers that the populations of many species are fragmented into patches connected by the migration of individuals through an interterritorial matrix. We applied fuzzy set theory and environmental favorability (F) functions to reveal the metapopulational structure of the 222 butterfly species in the Iberian Peninsula. We used the sets of contiguous grid cells with high favorability (F ≥ 0.8), to identify the favorable patches for each species. We superimposed the known occurrence data to reveal the occupied and empty favorable patches, as unoccupied patches are functional in a metapopulation dynamics analysis. We analyzed the connectivity between patches of each metapopulation by focusing on the territory of intermediate and low favorability for the species (F < 0.8). The friction that each cell opposes to the passage of individuals was computed as 1-F. We used the r.cost function of QGIS to calculate the cost of reaching each cell from a favorable patch. The inverse of the cost was computed as connectivity. Only 126 species can be considered to have a metapopulation structure. These metapopulation structures are part of the dark biodiversity of butterflies because their identification is not evident from the observation of the occurrence data but was revealed using favorability functions.

Metapopulation Structure of Diatom-associated Marine Bacteria

Marine bacteria-phytoplankton interaction ultimately shapes ecosystem productivity. The biochemical mechanisms underlying their interactions become increasingly known, yet how these ubiquitous interactions drive bacterial evolution has not been illustrated. Here, we sequenced genomes of 294 bacterial isolates associated with 19 coexisting diatom cells. These bacteria constitute eight genetically monomorphic populations of the globally abundant Roseobacter group. Six of these populations are members of Sulfitobacter, arguably the most prevalent bacteria associated with marine diatoms. A key finding is that populations varying at the intra-specific level have been differentiated and each are either associated with a single diatom host or with multiple hosts not overlapping with those of other populations. These closely related populations further show functional differentiation; they differ in motility phenotype and they harbor distinct types of secretion systems with implication for mediating organismal interactions. This interesting host-dependent population structure is even evident for demes within a genetically monomorphic population but each associated with a distinct diatom cell, as shown by a greater similarity in genome content between isolates from the same host compared to those from different hosts. Importantly, the intra- and inter-population differentiation pattern remains when the analyses are restricted to isolates from intra-specific diatom hosts, ruling out distinct selective pressures and instead suggesting coexisting microalgal cells as physical barriers of bacterial gene flow. Taken together, microalgae-associated bacteria display a unique microscale metapopulation structure, which consists of numerous small populations whose evolution is driven by random genetic drift.

Porównawcza ocena programów analizy żywotności populacji (PVA) w rankingu scenariuszy przekształceń krajobrazu = A comparative assessment of PVA software packages applied to rank the landscape management scenarios

Because of the scale and speed of species extinctions conservationists require methods that facilitate decision making. Therefore, a wide range of habitat and population viability analysis (PVA) software has been developed. Given the diversity of available programs it is currently challenging to decide which program is the most appropriate for a particular problem and what has to be considered when interpreting and comparing results from different approaches. Previous comparisons of PVA software addressed more generic questions such as data requirements, assumptions and predictive accuracy. In contract, we focus on a more applied problem that is still unresolved: how do simple habitat models and PVA software packages affect the ranking of alternative management scenarios? We addressed this problem by comparing different packages (LARCH, META-X, VORTEX and RAMAS GIS). As a test case, we studied the impact of alternative landscape development scenarios (river regulation, grassland restoration, reforestation and renaturalisation) for the Vistula valley, Poland, on the natterjack toad (Bufo calamita). In this context we also aimed to assess whether the use of at least two different PVA packages can enable users to better understand the differences in model predictions, which would imply a greater awareness and critical use of the packages. Our model selection represents different approaches to population viability analysis, including habitat, local population and stochastic patch occupancy models. The models can be evaluated in regard to the complexity of parameters and to the way the landscape is handled. We used RAMAS GIS to create a habitat model (RAMASh) and a detailed spatially explicit stochastic metapopulation model (RAMASp) which combined served as a complete “virtual” dataset for parameterisation of other programs. As an example of a stochastic patch occupancy model, we selected the META-X software. For a more independent comparison we added VORTEX – another package that includes explicit population dynamics, similar to RAMAS. Additionally, we included the habitat model LARCH because this type of model is often used by policy makers. We compared the metapopulation structure produced by RAMASh and LARCH. Scenario ranking according to the predicted carrying capacity in both programs was exactly the same, because the quantitative results for each scenario were almost identical in both programs. However, the metapopulation structure showed big differences between the programs, especially in the number of small populations. The analyses of results of different PVA programs (RAMASp, VORTEX and META-X) showed that absolute values of viability measures partly differed among these programs. Slight differences in population growth rate in RAMASp and VORTEX were amplified by stochasticity and resulted in visibly lower values of final abundance in VORTEX than in RAMASp. Also the absolute values of intrinsic mean time to extinction showed some discrepancies in VORTEX and META-X. These results are in agreement with findings of previous PVA comparisons, which emphasizes that absolute values of viability measures produced by any single model should be treated with caution. Nevertheless, despite these differences the rankings of the scenarios were the same in all three programs. However the order of the scenarios was different than in habitat models. In addition, these rankings were robust to the choice of viability measure. Taken together, these results emphasize that scenario ranking delivered by PVA software is robust and thus very useful for conservation management. Furthermore, we recommend using at least two PVA software packages in parallel, as this forces user to scrutinize the simplifying assumptions of the underlying models and of the viability metrics used.

Comparative genomics revealed adaptive admixture in Cryptosporidium hominis in Africa

Cryptosporidiosis is a major cause of diarrhoeal illness among African children, and is associated with childhood mortality, malnutrition, cognitive development and growth retardation. Cryptosporidium hominis is the dominant pathogen in Africa, and genotyping at the glycoprotein 60 (gp60) gene has revealed a complex distribution of different subtypes across this continent. However, a comprehensive exploration of the metapopulation structure and evolution based on whole-genome data has yet to be performed. Here, we sequenced and analysed the genomes of 26 C. hominis isolates, representing different gp60 subtypes, collected at rural sites in Gabon, Ghana, Madagascar and Tanzania. Phylogenetic and cluster analyses based on single-nucleotide polymorphisms showed that isolates predominantly clustered by their country of origin, irrespective of their gp60 subtype. We found a significant isolation-by-distance signature that shows the importance of local transmission, but we also detected evidence of hybridization between isolates of different geographical regions. We identified 37 outlier genes with exceptionally high nucleotide diversity, and this group is significantly enriched for genes encoding extracellular proteins and signal peptides. Furthermore, these genes are found more often than expected in recombinant regions, and they show a distinct signature of positive or balancing selection. We conclude that: (1) the metapopulation structure of C. hominis can only be accurately captured by whole-genome analyses; (2) local anthroponotic transmission underpins the spread of this pathogen in Africa; (3) hybridization occurs between distinct geographical lineages; and (4) genetic introgression provides novel substrate for positive or balancing selection in genes involved in host–parasite coevolution.