Population genomics of rapidly invading lionfish in the Caribbean reveals signals of range expansion in the absence of spatial population structure

Invasion genomics of lionfish in the Mediterranean Sea

The rate of biological invasions is steadily increasing, with major ecological and economic impacts accounting for billions of dollars in damage as a result. One spectacular example is the western Atlantic invasion by lionfishes. In the Mediterranean Sea, invasions from the Red Sea via the Suez Canal (termed Lessepsian invasions) comprise more than 100 fish species, including a recent invasion by lionfish. In light of the devastating effects of lionfish in the Caribbean Sea, understanding the dynamics of Mediterranean lionfish invasion is crucial. The Lessepsian lionfish invasion started in 2012, and rapidly spread to the central Mediterranean. Here, we used thousands of RAD seq genomic markers to study the population dynamics of this invasion. While we did not find a reduction in genetic diversity between source (Red Sea) and invasive (Mediterranean) populations (i.e., bottleneck effects), we found evidence of population structure within the invasive range in the Mediterranean Sea. We found that loci that are potentially under selection may play an important role in invasion success (in particular, genes involved in osmoregulation and fin spine sizes). Genomic approaches proved powerful in examining the ecological and evolutionary patterns of successful invaders and may be used as tools to understand and potentially mitigate future invasions.

Computational reproducibility of Jupyter notebooks from biomedical publications

BACKGROUND

Jupyter notebooks facilitate the bundling of executable code with its documentation and output in one interactive environment, and they represent a popular mechanism to document and share computational workflows, including for research publications. The reproducibility of computational aspects of research is a key component of scientific reproducibility but has not yet been assessed at scale for Jupyter notebooks associated with biomedical publications.

APPROACH

We address computational reproducibility at 2 levels: (i) using fully automated workflows, we analyzed the computational reproducibility of Jupyter notebooks associated with publications indexed in the biomedical literature repository PubMed Central. We identified such notebooks by mining the article's full text, trying to locate them on GitHub, and attempting to rerun them in an environment as close to the original as possible. We documented reproduction success and exceptions and explored relationships between notebook reproducibility and variables related to the notebooks or publications. (ii) This study represents a reproducibility attempt in and of itself, using essentially the same methodology twice on PubMed Central over the course of 2 years, during which the corpus of Jupyter notebooks from articles indexed in PubMed Central has grown in a highly dynamic fashion.

RESULTS

Out of 27,271 Jupyter notebooks from 2,660 GitHub repositories associated with 3,467 publications, 22,578 notebooks were written in Python, including 15,817 that had their dependencies declared in standard requirement files and that we attempted to rerun automatically. For 10,388 of these, all declared dependencies could be installed successfully, and we reran them to assess reproducibility. Of these, 1,203 notebooks ran through without any errors, including 879 that produced results identical to those reported in the original notebook and 324 for which our results differed from the originally reported ones. Running the other notebooks resulted in exceptions.

CONCLUSIONS

We zoom in on common problems and practices, highlight trends, and discuss potential improvements to Jupyter-related workflows associated with biomedical publications.

Coral Reef Population Genomics in an Age of Global Change

Coral reefs are both exceptionally biodiverse and threatened by climate change and other human activities. Here, we review population genomic processes in coral reef taxa and their importance for understanding responses to global change. Many taxa on coral reefs are characterized by weak genetic drift, extensive gene flow, and strong selection from complex biotic and abiotic environments, which together present a fascinating test of microevolutionary theory. Selection, gene flow, and hybridization have played and will continue to play an important role in the adaptation or extinction of coral reef taxa in the face of rapid environmental change, but research remains exceptionally limited compared to the urgent needs. Critical areas for future investigation include understanding evolutionary potential and the mechanisms of local adaptation, developing historical baselines, and building greater research capacity in the countries where most reef diversity is concentrated.

Selection on growth rate and local adaptation drive genomic adaptation during experimental range expansions in the protist Tetrahymena thermophila

Populations that expand their range can undergo rapid evolutionary adaptation of life‐history traits, dispersal behaviour and adaptation to the local environment. Such adaptation may be aided or hindered by sexual reproduction, depending on the context.

However, few empirical and experimental studies have investigated the genetic basis of adaptive evolution during range expansions. Even less attention has been given to the question how sexual reproduction may modulate such adaptive evolution during range expansions.

We here studied genomic adaptation during experimental range expansions of the protist Tetrahymena thermophila in landscapes with a uniform environment or a pH gradient. Specifically, we investigated two aspects of genomic adaptation during range expansion. First, we investigated adaptive genetic change in terms of the underlying numbers of allele frequency changes from standing genetic variation and de novo variants. We focused on how sexual reproduction may alter this adaptive genetic change. Second, we identified genes subject to selection caused by the expanding range itself, and directional selection due to the presence or absence of the pH gradient. We focused this analysis on alleles with large frequency changes that occurred in parallel in more than one population to identify the most likely candidate targets of selection.

We found that sexual reproduction altered adaptive genetic change both in terms of de novo variants and standing genetic variation. However, sexual reproduction affected allele frequency changes in standing genetic variation only in the absence of long‐distance gene flow. Adaptation to the range expansion affected genes involved in cell divisions and DNA repair, whereas adaptation to the pH gradient additionally affected genes involved in ion balance and oxidoreductase reactions. These genetic changes may result from selection on growth and adaptation to low pH.

In the absence of gene flow, sexual reproduction may have aided genetic adaptation. Gene flow may have swamped expanding populations with maladapted alleles, thus reducing the extent of evolutionary adaptation during range expansion. Sexual reproduction also altered the genetic basis of adaptation in our evolving populations via de novo variants, possibly by purging deleterious mutations or by revealing fitness benefits of rare genetic variants.

Balanced polymorphism fuels rapid selection in an invasive crab despite high gene flow and low genetic diversity

Adaptation across environmental gradients has been demonstrated in numerous systems with extensive dispersal, despite high gene flow and consequently low genetic structure. The speed and mechanisms by which such adaptation occurs remain poorly resolved, but are critical to understanding species spread and persistence in a changing world. Here, we investigate these mechanisms in the European green crab Carcinus maenas, a globally distributed invader. We focus on a northwestern Pacific population that spread across >12 degrees of latitude in 10 years from a single source, following its introduction <35 years ago. Using six locations spanning >1500 km, we examine genetic structure using 9376 single nucleotide polymorphisms (SNPs). We find high connectivity among five locations, with significant structure between these locations and an enclosed lagoon with limited connectivity to the coast. Among the five highly connected locations, the only structure observed was a cline driven by a handful of SNPs strongly associated with latitude and winter temperature. These SNPs are almost exclusively found in a large cluster of genes in strong linkage disequilibrium that was previously identified as a candidate for cold tolerance adaptation in this species. This region may represent a balanced polymorphism that evolved to promote rapid adaptation in variable environments despite high gene flow, and which now contributes to successful invasion and spread in a novel environment. This research suggests an answer to the paradox of genetically depauperate yet successful invaders: populations may be able to adapt via a few variants of large effect despite low overall diversity.

Protect the Natives to Combat the Aliens: Could Octopus vulgaris Cuvier, 1797 Be a Natural Agent for the Control of the Lionfish Invasion in the Mediterranean Sea?

Biological invasions constitute a major threat to native ecosystems and to global biodiversity [...]

Factors influencing taxonomic unevenness in scientific research: a mixed-methods case study of non-human primate genomic sequence data generation

Scholars have noted major disparities in the extent of scientific research conducted among taxonomic groups. Such trends may cascade if future scientists gravitate towards study species with more data and resources already available. As new technologies emerge, do research studies employing these technologies continue these disparities? Here, using non-human primates as a case study, we identified disparities in massively parallel genomic sequencing data and conducted interviews with scientists who produced these data to learn their motivations when selecting study species. We tested whether variables including publication history and conservation status were significantly correlated with publicly available sequence data in the NCBI Sequence Read Archive (SRA). Of the 179.6 terabases (Tb) of sequence data in SRA for 519 non-human primate species, 135 Tb (approx. 75%) were from only five species: rhesus macaques, olive baboons, green monkeys, chimpanzees and crab-eating macaques. The strongest predictors of the amount of genomic data were the total number of non-medical publications (linear regression; r ² = 0.37; p = 6.15 × 10^-12) and number of medical publications (r ² = 0.27; p = 9.27 × 10^-9). In a generalized linear model, the number of non-medical publications (p = 0.00064) and closer phylogenetic distance to humans (p = 0.024) were the most predictive of the amount of genomic sequence data. We interviewed 33 authors of genomic data-producing publications and analysed their responses using grounded theory. Consistent with our quantitative results, authors mentioned their choice of species was motivated by sample accessibility, prior published work and relevance to human medicine. Our mixed-methods approach helped identify and contextualize some of the driving factors behind species-uneven patterns of scientific research, which can now be considered by funding agencies, scientific societies and research teams aiming to align their broader goals with future data generation efforts.

Genetic connectivity of lionfish (Pterois volitans) in marine protected areas of the Gulf of Mexico and Caribbean Sea

Lionfish (Pterois volitans) have rapidly invaded the tropical Atlantic and spread across the wider Caribbean in a relatively short period of time. Because of its high invasion capacity, we used it as a model to identify the connectivity among nine marine protected areas (MPAs) situated in four countries in the Gulf of Mexico and the Caribbean Sea. This study provides evidence of local genetic differentiation of P. volitans in the Gulf of Mexico and the Caribbean Sea. A total of 475 lionfish samples were characterized with 12 microsatellites, with 6-20 alleles per locus. Departures from Hardy-Weinberg equilibrium (HWE) were found in 10 of the 12 loci, all caused by heterozygous excess. Moderate genetic differentiation was observed between Chiriviche, Venezuela and Xcalak, México localities (F ST = 0.012), and between the Los Roques and the Veracruz (F ST = 0.074) sites. STRUCTURE analysis found that four genetic entities best fit our data. A unique genetic group in the Gulf of Mexico may imply that the lionfish invasion unfolded both in a counterclockwise manner in the Gulf of Mexico. In spite of the notable dispersion of P. volitans, our results show some genetic structure, as do other noninvasive Caribbean fish species, suggesting that the connectivity in some MPAs analyzed in the Caribbean is limited and caused by only a few source individuals with subsequent genetic drift leading to local genetic differentiation. This indicates that P. volitans dispersion could be caused by mesoscale phenomena, which produce stochastic connectivity pulses. Due to the isolation of some MPAs from others, these findings may hold a promise for local short-term control of by means of intensive fishing, even in MPAs, and may have regional long-term effects.

Lionfish (Pterois volitans) as biomonitoring species for oil pollution effects in coral reef ecosystems

With oil spills, and other sources of aromatic hydrocarbons, being a continuous threat to coral reef systems, and most reef fish species being protected or difficult to collect, the use of the invasive lionfish (Pterois volitans) might be a good model species to monitor biomarkers in potentially exposed fish in the Caribbean and western Atlantic. The rapid expansion of lionfish in the Caribbean and western Atlantic, and the unregulated fishing for this species, would make the lionfish a suitable candidate as biomonitoring species for oil pollution effects. However, to date little has been published about the responses of lionfish to environmental pollutants. For this study lionfish were collected in the Florida Keys a few weeks after Hurricane Irma, which sank numerous boats resulting in leaks of oil and fuel, and during the winter and early spring after that. Several biomarkers indicative of exposure to PAHs (bile fluorescence, cytochrome P450-1A induction, glutathione S-transferase activity) were measured. To establish if these biomarkers are inducible in PAH exposed lionfish, dosing experiments with different concentrations of High Energy Water Accommodated Fraction of crude oil were performed. The results revealed no significant effects in the biomarkers in the field collected fish, while the exposure experiments demonstrated that lionfish did show strong effects in the measured biomarkers, even at the lowest concentration tested (0.3% HEWAF, or 25 μg/l ƩPAH50). Based on its widespread distribution, relative ease of collection, and significant biomarker responses in the controlled dosing experiment, it is concluded that lionfish has good potential to be used as a standardized biomonitoring species for oil pollution in its neotropical realm.