Evaluating alternative explanations for an association of extinction risk and evolutionary uniqueness in multiple insular lineages

Genomic and ecological evidence shed light on the recent demographic history of two related invasive insects

Hypogeococcus pungens is a species complex native to southern South America that is composed of at least five putative species, each one specialized in the use of different host plants. Two of these undescribed species were registered as invasive in Central and North America: Hyp-C is a cactophagous mealybug that became an important pest that threatens endemic cactus species in Puerto Rico, and Hyp-AP feeds on Amaranthaceae and Portulacaceae hosts, but does not produce severe damage to the host plants. We quantified genomic variation and investigated the demographic history of both invasive species by means of coalescent-based simulations using high throughput sequencing data. We also evaluated the incidence of host plant infestation produced by both species and used an ecological niche modeling approach to assess potential distribution under current and future climatic scenarios. Our genetic survey evinced the footprints of strong effective population size reduction and signals of genetic differentiation among populations within each species. Incidence of plant attacks varied between species and among populations within species, with some host plant species preferred over others. Ecological niche modeling suggested that under future climatic scenarios both species would expand their distribution ranges in Puerto Rico. These results provide valuable information for the design of efficient management and control strategies of the Puerto Rican cactus pest and shed light on the evolutionary pathways of biological invasions.

gen3sis: A general engine for eco-evolutionary simulations of the processes that shape Earth's biodiversity

Understanding the origins of biodiversity has been an aspiration since the days of early naturalists. The immense complexity of ecological, evolutionary, and spatial processes, however, has made this goal elusive to this day. Computer models serve progress in many scientific fields, but in the fields of macroecology and macroevolution, eco-evolutionary models are comparatively less developed. We present a general, spatially explicit, eco-evolutionary engine with a modular implementation that enables the modeling of multiple macroecological and macroevolutionary processes and feedbacks across representative spatiotemporally dynamic landscapes. Modeled processes can include species' abiotic tolerances, biotic interactions, dispersal, speciation, and evolution of ecological traits. Commonly observed biodiversity patterns, such as α, β, and γ diversity, species ranges, ecological traits, and phylogenies, emerge as simulations proceed. As an illustration, we examine alternative hypotheses expected to have shaped the latitudinal diversity gradient (LDG) during the Earth's Cenozoic era. Our exploratory simulations simultaneously produce multiple realistic biodiversity patterns, such as the LDG, current species richness, and range size frequencies, as well as phylogenetic metrics. The model engine is open source and available as an R package, enabling future exploration of various landscapes and biological processes, while outputs can be linked with a variety of empirical biodiversity patterns. This work represents a key toward a numeric, interdisciplinary, and mechanistic understanding of the physical and biological processes that shape Earth's biodiversity.

Comparative transcriptomics of ice-crawlers demonstrates cold specialization constrains niche evolution in a relict lineage

Key changes in ecological niche space are often critical to understanding how lineages diversify during adaptive radiations. However, the converse, or understanding why some lineages are depauperate and relictual, is more challenging, as many factors may constrain niche evolution. In the case of the insect order Grylloblattodea, highly conserved thermal breadth is assumed to be closely tied to their relictual status, but has not been formerly tested. Here, we investigate whether evolutionary constraints in the physiological tolerance of temperature can help explain relictualism in this lineage. Using a comparative transcriptomics approach, we investigate gene expression following acute heat and cold stress across members of Grylloblattodea and their sister group, Mantophasmatodea. We additionally examine patterns of protein evolution, to identify candidate genes of positive selection. We demonstrate that cold specialization in Grylloblattodea has been accompanied by the loss of the inducible heat shock response under both acute heat and cold stress. Additionally, there is widespread evidence of selection on protein-coding genes consistent with evolutionary constraints due to cold specialization. This includes positive selection on genes involved in trehalose transport, metabolic function, mitochondrial function, oxygen reduction, oxidative stress, and protein synthesis. These patterns of molecular adaptation suggest that Grylloblattodea have undergone evolutionary trade-offs to survive in cold habitats and should be considered highly vulnerable to climate change. Finally, our transcriptomic data provide a robust backbone phylogeny for generic relationships within Grylloblattodea and Mantophasmatodea. Major phylogenetic splits in each group relate to arid conditions driving biogeographical patterns, with support for a sister-group relationship between North American Grylloblatta and Altai-Sayan Grylloblattella, and a range disjunction in Namibia splitting major clades within Mantophasmatodea.

Island songbirds as windows into evolution in small populations

Due to their limited ranges and inherent isolation, island species have long been recognized as crucial systems for tackling a range of evolutionary questions, including in the early study of speciation.^1,2 Such species have been less studied in the understanding of the evolutionary forces driving DNA sequence evolution. Island species usually have lower census population sizes (N) than continental species and, supposedly, lower effective population sizes (Ne). Given that both the rates of change caused by genetic drift and by selection are dependent upon Ne, island species are theoretically expected to exhibit (1) lower genetic diversity, (2) less effective natural selection against slightly deleterious mutations,^3,4 and (3) a lower rate of adaptive evolution.^5-8 Here, we have used a large set of newly sequenced and published whole-genome sequences of Passerida species (14 insular and 11 continental) to test these predictions. We confirm that island species exhibit lower census size and Ne, supporting the hypothesis that the smaller area available on islands constrains the upper bound of Ne. In the insular species, we find lower nucleotide diversity in coding regions, higher ratios of non-synonymous to synonymous polymorphisms, and lower adaptive substitution rates. Our results provide robust evidence that the lower Ne experienced by island species has affected both the ability of natural selection to efficiently remove weakly deleterious mutations and also the adaptive potential of island species, therefore providing considerable empirical support for the nearly neutral theory. We discuss the implications for both evolutionary and conservation biology.

Geo-Climatic Changes and Apomixis as Major Drivers of Diversification in the Mediterranean Sea Lavenders (Limonium Mill.)

The Mediterranean realm, comprising the Mediterranean and Macaronesian regions, has long been recognized as one of the world's biodiversity hotspots, owing to its remarkable species richness and endemism. Several hypotheses on biotic and abiotic drivers of species diversification in the region have been often proposed but rarely tested in an explicit phylogenetic framework. Here, we investigate the impact of both species-intrinsic and -extrinsic factors on diversification in the species-rich, cosmopolitan Limonium, an angiosperm genus with center of diversity in the Mediterranean. First, we infer and time-calibrate the largest Limonium phylogeny to date. We then estimate ancestral ranges and diversification dynamics at both global and regional scales. At the global scale, we test whether the identified shifts in diversification rates are linked to specific geological and/or climatic events in the Mediterranean area and/or asexual reproduction (apomixis). Our results support a late Paleogene origin in the proto-Mediterranean area for Limonium, followed by extensive in situ diversification in the Mediterranean region during the late Miocene, Pliocene, and Pleistocene. We found significant increases of diversification rates in the "Mediterranean lineage" associated with the Messinian Salinity Crisis, onset of Mediterranean climate, Plio-Pleistocene sea-level fluctuations, and apomixis. Additionally, the Euro-Mediterranean area acted as the major source of species dispersals to the surrounding areas. At the regional scale, we infer the biogeographic origins of insular endemics in the oceanic archipelagos of Macaronesia, and test whether woodiness in the Canarian Nobiles clade is a derived trait linked to insular life and a biotic driver of diversification. We find that Limonium species diversity on the Canary Islands and Cape Verde archipelagos is the product of multiple colonization events followed by in situ diversification, and that woodiness of the Canarian endemics is indeed a derived trait but is not associated with a significant shift to higher diversification rates. Our study expands knowledge on how the interaction between abiotic and biotic drivers shape the uneven distribution of species diversity across taxonomic and geographical scales.