Thermal selection as a driver of marine ecological speciation

Population genetic structure and ecological differentiation in the bryozoan genus Reteporella across the Azores Archipelago (central North Atlantic)

The processes shaping population dynamics of benthic marine invertebrates with non-planktotrophic larvae are still poorly understood but have seen a renewed interest in applying integrative taxonomic approaches. We used mitochondrial and microsatellite (SSR-GBAS) data to estimate connectivity across islands and seamounts in the central North Atlantic Azores Archipelago in five species of the bryozoan genus Reteporella Busk, 1884. Discordant patterns were inferred between datasets, which might be due to methodological constraints related to the application of multilocus approaches based on amplification to multiple species or due to interspecific introgression in deep waters. A divergent cryptic ecotype of Reteporella atlantica (Busk, 1884) was found in shallow waters, likely resulting from ecologically-driven incipient speciation, posing new questions regarding the role of bathymetrical zonation as a promoter of differentiation. The occurrence of ecologically-driven differentiation and potential interspecific introgression in other bryozoans should be considered, both with potentially important evolutionary and biogeographical consequences. The discovery of incipient species, prompted by ecological factors, calls for the need to consider marine invertebrates when developing conservation strategies in oceanic insular ecosystems.

Rolling with the punches: Organism-environment interactions shape spatial pattern of adaptive differentiation in the widespread mantis shrimp Oratosquilla oratoria

Investigating spatial pattern of adaptive variation and its underlying processes can inform the adaptive potential distributed within species ranges, which is increasingly important in the context of a changing climate. A correct interpretation of adaptive variation pattern requires that population history and the ensuing population genetic structure are taken into account. Here we carried out such a study by integrating population genomic analyses, demographic model testing and species distribution modeling to investigate patterns and causes of adaptive differentiation in a widespread mantis shrimp, Oratosquilla oratoria, along a replicated, broad-scale temperature gradient in the northwestern Pacific (NWP). Our results supported a strong hierarchical ecogeographic structure dominated by habitat-linked divergence among O. oratoria populations accompanied with introgressive hybridization. A combined FST outlier and environmental correlation analyses revealed remarkable temperature-associated clines in allele frequency across paired North-South populations on Chinese and Japanese coasts, and identified a suite of loci associated with temperature adaptation. Further demographic model testing revealed the observed clinal variation derived partly from Pleistocene divergence followed by recent secondary contact. More importantly, the likelihood of hybridization is predicted to increase as climate change progresses, which would break barriers to gene flow and enable the spread of adaptive genetic variation. These results support that not only is temperature-driven adaptive differentiation occurs in O. oratoria but that such pattern is likely attributed to ancient adaptive variation, sustained by contemporary ocean conditions and a semi-permeable barrier to gene flow maintained by selection. They moreover provide genomic insights into the distribution of adaptive potential across O. oratoria' s species range. This work can serve as a case study to characterize adaptive diversity of marine species in the NWP by integrating environmental and genetic data at temporal and spatial scales in a population genomic framework, which would improve management and conservation actions under climate change.

Seascapes Shaped the Local Adaptation and Population Structure of South China Coast Yellowfin Seabream (Acanthopagrus latus)

Understanding the genetic composition and regional adaptation of marine species under environmental heterogeneity and fishing pressure is crucial for responsible management. In order to understand the genetic diversity and adaptability of yellowfin seabream (Acanthopagrus latus) along southern China coast, this study was conducted a seascape genome analysis on yellowfin seabream from the ecologically diverse coast, spanning over 1600 km. A total of 92 yellowfin seabream individuals from 15 sites were performed whole-genome resequencing, and 4,383,564 high-quality single nucleotide polymorphisms (SNPs) were called. By conducting a genotype-environment association analysis, 29,951 adaptive and 4,328,299 neutral SNPs were identified. The yellowfin seabream exhibited two distinct population structures, despite high gene flow between sites. The seascape genome analysis revealed that genetic structure was influenced by a variety of factors including salinity gradients, habitat distance, and ocean currents. The frequency of allelic variation at the candidate loci changed with the salinity gradient. Annotation of these loci revealed that most of the genes are associated with osmoregulation, such as kcnab2a, kcnk5a, and slc47a1. These genes are significantly enriched in pathways associated with ion transport including G protein-coupled receptor activity, transmembrane signaling receptor activity, and transporter activity. Overall, our findings provide insights into how seascape heterogeneity affects adaptive evolution, while providing important information for regional management in yellowfin seabream populations.

Local adaptation with gene flow in a highly dispersive shark

Adaptive divergence in response to environmental clines are expected to be common in species occupying heterogeneous environments. Despite numerous advances in techniques appropriate for non-model species, gene-environment association studies in elasmobranchs are still scarce. The bronze whaler or copper shark (Carcharhinus brachyurus) is a large coastal shark with a wide distribution and one of the most exploited elasmobranchs in southern Africa. Here, we assessed the distribution of neutral and adaptive genomic diversity in C. brachyurus across a highly heterogeneous environment in southern Africa based on genome-wide SNPs obtained through a restriction site-associated DNA method (3RAD). A combination of differentiation-based genome-scan (outflank) and genotype-environment analyses (redundancy analysis, latent factor mixed models) identified a total of 234 differentiation-based outlier and candidate SNPs associated with bioclimatic variables. Analysis of 26,299 putatively neutral SNPs revealed moderate and evenly distributed levels of genomic diversity across sites from the east coast of South Africa to Angola. Multivariate and clustering analyses demonstrated a high degree of gene flow with no significant population structuring among or within ocean basins. In contrast, the putatively adaptive SNPs demonstrated the presence of two clusters and deep divergence between Angola and all other individuals from Namibia and South Africa. These results provide evidence for adaptive divergence in response to a heterogeneous seascape in a large, mobile shark despite high levels of gene flow. These results are expected to inform management strategies and policy at the national and regional level for conservation of C. brachyurus populations.

Insights into the Mitochondrial Genetic Makeup and Miocene Colonization of Primitive Flatfishes (Pleuronectiformes: Psettodidae) in the East Atlantic and Indo-West Pacific Ocean

The mitogenomic evolution of the Psettodes flatfishes is still poorly known from their range distribution in eastern Atlantic and Indo-West Pacific Oceans. The study delves into the matrilineal evolutionary pathway of these primitive flatfishes, with a specific focus on the complete mitogenome of the Psettodes belcheri species, as determined through next-generation sequencing. The mitogenome in question spans a length of 16,747 base pairs and comprises a total of 37 genes, including 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. Notably, the mitogenome of P. belcheri exhibits a bias towards AT base pairs, with a composition of 54.15%, mirroring a similar bias observed in its close relative, Psettodes erumei, which showcases percentages of 53.07% and 53.61%. Most of the protein-coding genes commence with an ATG initiation codon, except for Cytochrome c oxidase I (COI), which initiates with a GTG codon. Additionally, four protein-coding genes commence with a TAA termination codon, while seven others exhibit incomplete termination codons. Furthermore, two protein-coding genes, namely NAD1 and NAD6, terminate with AGG and TAG stop codons, respectively. In the mitogenome of P. belcheri, the majority of transfer RNAs demonstrate the classical cloverleaf secondary structures, except for tRNA-serine, which lacks a DHU stem. Comparative analysis of conserved blocks within the control regions of two Psettodidae species unveiled that the CSB-II block extended to a length of 51 base pairs, surpassing the other blocks and encompassing highly variable sites. A comprehensive phylogenetic analysis using mitochondrial genomes (13 concatenated PCGs) categorized various Pleuronectiformes species, highlighting the basal position of the Psettodidae family and showed monophyletic clustering of Psettodes species. The approximate divergence time (35-10 MYA) between P. belcheri and P. erumei was estimated, providing insights into their separation and colonization during the early Miocene. The TimeTree analysis also estimated the divergence of two suborders, Psettodoidei and Pleuronectoidei, during the late Paleocene to early Eocene (56.87 MYA). The distribution patterns of Psettodes flatfishes were influenced by ocean currents and environmental conditions, contributing to their ecological speciation. In the face of climate change and anthropogenic activities, the conservation implications of Psettodes flatfishes are emphasized, underscoring the need for regulated harvesting and adaptive management strategies to ensure their survival in changing marine ecosystems. Overall, this study contributes to understanding the evolutionary history, genetic diversity, and conservation needs of Psettodes flatfishes globally. However, the multifaceted exploration of mitogenome and larger-scale genomic data of Psettodes flatfish will provide invaluable insights into their genetic characterization, evolutionary history, environmental adaptation, and conservation in the eastern Atlantic and Indo-West Pacific Oceans.

Thermal niche partitioning and phenology of Nearctic and Palearctic flea (Siphonaptera) communities on rodents (Mammalia: Rodentia) from five ecoregions

Seasonality of fleas (Siphonaptera) may be due to species competition, prompting the idea that flea species partition temperature, along with correlated variables such as moisture (thermal-niche partitioning hypothesis). I compared the fleas of five rodent-flea communities described from the literature for thermal-niche optima by fitting non-linear LRF (Lobry-Rosso-Flandrois) curves to examine whether flea species in a community show distinct, partitioned thermal niches. LRF curves estimate physiological parameters of temperature minimum, optimum, maximum, and maximum abundance, and facilitate comparison between species by summarizing seasonal data. Flea-communities were on Nearctic Southern flying squirrel (Glaucomys volans volans), Richardson's ground-squirrel (Urocitellus richardsonii), North American deer-mouse (Peromyscus maniculatus), and Palearctic Midday jird (Meriones meridianus), and Wagner's gerbil (Dipodillus dasyurus). Flea communities appeared to show seasonality consistent with thermal-niche partitioning. Several flea families and genera had characteristic thermal niches: Ceratophyllidae had broad tolerance to extreme temperature, Leptopsyllidae (one species in this study) to cold, and Pulicidae to hot. In contrast, at the local, species level, climatic speciation could be significant in flea diversification. Non-competition hypotheses (environmental filtering, neutrality) require testing, too. Thermal-niche partitioning may increase flea species richness on hosts and could occur in other insect and plant communities. Implications for biodiversity conservation and disease ecology under global warming are wide-ranging.

Lack of a genetic cline and temporal genetic stability in an introduced barnacle along the Pacific coast of Japan

Background

Large numbers of exotic marine species have been introduced worldwide. Monitoring of introduced species is important to reveal mechanisms underlying their establishment and expansion. Balanus glandula is a common intertidal barnacle native to the northeastern Pacific. However, this species has been introduced to Japan, South America, South Africa, and Europe. While a latitudinal genetic cline is well known in its native range, it is unclear whether such a genetic cline occurs in introduced areas. Twenty years have passed since it was first identified in Japan and its distribution now ranges from temperate to subarctic regions.

Methods

In the present study, we examined genotypes of cytochrome oxidase subunit I (COI) of mitochondrial (mt)-DNA and elongation factor 1a (EF1) across the distribution of B. glandula in Japan at high and mid intertidal zones.

Results

At all sampling sites, native northern genotypes are abundant and I did not detect significant effects of latitude, tide levels, or their interaction on genotypic frequencies. Further, I did not detect any change of genotype composition between data collected during a study in 2004 and samples in the present study collected in 2019. Data from the present study offer an important baseline for future monitoring of this species and supply valuable insights into the mechanisms of establishment and expansion of introduced marine species generally.

Speciation genomics and the role of depth in the divergence of rockfishes (Sebastes) revealed through Pool-seq analysis of enriched sequences

Speciation in the marine environment is challenged by the wide geographic distribution of many taxa and potential for high rates of gene flow through larval dispersal mechanisms. Depth has recently been proposed as a potential driver of ecological divergence in fishes, and yet it is unclear how adaptation along these gradients' shapes genomic divergence. The genus Sebastes contains numerous species pairs that are depth-segregated and can provide a better understanding of the mode and tempo of genomic diversification. Here, we present exome data on two species pairs of rockfishes that are depth-segregated and have different degrees of divergence: S. chlorostictus-S. rosenblatti and S. crocotulus-S. miniatus. We were able to reliably identify "islands of divergence" in the species pair with more recent divergence (S. chlorostictus-S. rosenblatti) and discovered a number of genes associated with neurosensory function, suggesting a role for this pathway in the early speciation process. We also reconstructed demographic histories of divergence and found the best supported model was isolation followed by asymmetric secondary contact for both species pairs. These results suggest past ecological/geographic isolation followed by asymmetric secondary contact of deep to shallow species. Our results provide another example of using rockfish as a model for studying speciation and support the role of depth as an important mechanism for diversification in the marine environment.

Demographic, physiological, and genetic factors linked to the poleward range expansion of the snail Nerita yoldii along the shoreline of China

Species range shift is one of the most significant consequences of climate change in the Anthropocene. A comprehensive study, including demographic, physiological, and genetic factors linked to poleward range expansion, is crucial for understanding how the expanding population occupies the new habitat. In the present study, we investigated the demographic, physiological, and genetic features of the intertidal gastropod Nerita yoldii, which has extended its northern limit by ~200 km over the former biogeographic break of the Yangtze River Estuary during recent decades. The neutral SNPs data showed that the new marginal populations formed a distinct cluster established by a few founders. Demographic modelling analysis revealed that the new marginal populations experienced a strong genetic bottleneck followed by recent demographic expansion. Successful expansion that overcame the founder effect might be attributed to its high capacity of rapid population growth and multiple introductions. According to the non-neutral SNPs under diversifying selection, there were high levels of heterozygosity in the new marginal populations, which might be beneficial for adapting to the novel thermal conditions. The common garden experiment showed that the new marginal populations have evolved divergent transcriptomic and physiological responses to heat stress, allowing them to occupy and survive in the novel environment. Lower transcriptional plasticity was observed in the new marginal populations. These results suggest a new biogeographic pattern of N. yoldii has formed with the occurrence of demographic, physiologic, and genetic changes, and emphasize the roles of adaptation of marginal populations during range expansion.

Ecological speciation by sympatric host shifts in a clade of herbivorous sea slugs, with introgression and localized mitochondrial capture between species

Host shifting in insect-plant systems was historically important to the development of ecological speciation theory, yet surprisingly few studies have examined whether host shifting drives diversification of marine herbivores. When small-bodied consumers feed and also mate on a preferred host, disruptive selection can split a population into host races despite gene flow. Support for host shifts is notably lacking for invertebrates associated with macroalgae, where the scale of dispersal by planktonic larvae often far exceeds the grain of host patchiness, and adults are typically less specialized than terrestrial herbivores. Here, we present a candidate example of ecological speciation in a clade of sea slugs that primarily consume green algae in the genus Caulerpa, including highly invasive species. Ancestral character state reconstructions supported 'sea grapes' (C. racemosa, C. lentillifera) as the ancestral host for a tropical radiation of 12 Elysia spp., with one shift onto alternative Caulerpa spp. in the Indo-Pacific. A Caribbean radiation of three species included symaptric host shifts to Rhipocephalus brevicaulis in the ancestor of E. pratensis Ortea & Espinosa, 1996, and to C. prolifera in E. hamanni Krug, Vendetti & Valdes 2016, plus a niche expansion to a range of Caulerpa spp. in E. subornata Verrill, 1901. All three species are broadly sympatric across the Caribbean but are host-partitioned at a fine grain, and distinct by morphology and at nuclear loci. However, non-recombining mtDNA revealed a history of gene flow between E. pratensis and E. subornata: COI haplotypes from E. subornata were 10.4% divergent from E. pratensis haplotypes from four sites, but closely related to all E. pratensis haplotypes sampled from six Bahamian islands, indicating historical introgression and localized "mitochondrial capture." Disruptive selective likely fueled divergence and adaptation to distinct host environments, indicating ecological speciation may be an under-appreciated driver of diversification for marine herbivores as well as epibionts and other resource specialists.

Comparative mitogenomics of Clupeoid fish provides insights into the adaptive evolution of mitochondrial oxidative phosphorylation (OXPHOS) genes and codon usage in the heterogeneous habitats

Clupeoid fish can be considered excellent candidates to understand the role of mitochondrial DNA in adaptive evolution, as they have colonized different habitats (marine, brackish, freshwater, tropical and temperate regions) over millions of years. Here, we investigate patterns of tRNA location, codon usage bias, and lineage-specific diversifying selection signals to provide novel insights into how evolutionary improvements of mitochondrial metabolic efficiency have allowed clupeids to adapt to different habitats. Based on whole mitogenome data of 70 Clupeoids with a global distribution we find that purifying selection was the dominant force acting and that the mutational deamination pressure in mtDNA was stronger than the codon/amino acid constraints. The codon usage pattern appears evolved to achieve high translational efficiency (codon/amino acid-related constraints), as indicated by the complementarity of most codons to the GT-saturated tRNA anticodon sites (retained by deamination-induced pressure) and usage of the codons of the tRNA genes situated near to the control region (fixed by deamination pressure) where transcription efficiency was high. The observed shift in codon preference patterns between marine and euryhaline/freshwater Clupeoids indicates possible selection for improved translational efficiency in mitochondrial genes while adapting to low-salinity habitats. This mitogenomic plasticity and enhanced efficiency of the metabolic machinery may have contributed to the evolutionary success and abundance of Clupeoid fish.

Seascape genomics of coastal bottlenose dolphins along strong gradients of temperature and salinity

Heterogeneous seascapes and strong environmental gradients in coastal waters are expected to influence adaptive divergence, particularly in species with large population sizes where selection is expected to be highly efficient. However, these influences might also extend to species characterized by strong social structure, natal philopatry and small home ranges. We implemented a seascape genomic study to test this hypothesis in Indo-Pacific bottlenose dolphins (Tursiops aduncus) distributed along the environmentally heterogeneous coast of southern Australia. The datasets included oceanographic and environmental variables thought to be good predictors of local adaptation in dolphins and 8,081 filtered single nucleotide polymorphisms (SNPs) genotyped for individuals sampled from seven different bioregions. From a neutral perspective, population structure and connectivity of the dolphins were generally influenced by habitat type and social structuring. Genotype-environment association analysis identified 241 candidate adaptive loci and revealed that sea surface temperature and salinity gradients influenced adaptive divergence in these animals at both large- (1,000s km) and fine-scales (<100 km). Enrichment analysis and annotation of candidate genes revealed functions related to sodium-activated ion transport, kidney development, adipogenesis and thermogenesis. The findings of spatial adaptive divergence and inferences of putative physiological adaptations challenge previous suggestions that marine megafauna is most likely to be affected by environmental and climatic changes via indirect, trophic effects. Our work contributes to conservation management of coastal bottlenose dolphins subjected to anthropogenic disturbance and to efforts of clarifying how seascape heterogeneity influences adaptive diversity and evolution in small cetaceans.

Coastal dunefields maintain pre-Holocene genetic structure in a rocky shore red alga

Most intertidal algae have limited dispersal potential, and areas that lack hard substratum suitable for attachment are thus expected to isolate regional populations from each other. Here, we used nuclear and mitochondrial genetic data to compare genetic structure in two co-distributed intertidal red algae with different dispersal potential along the South African coastline. Gelidium pristoides is divided into a south-eastern and a south-western evolutionary lineage separated by extensive, continuous sandy shoreline habitat adjacent to coastal dunefields. In contrast, Hypnea spicifera is genetically homogeneous throughout its range. In G. pristoides, the genetic breaks are associated with contemporary coastal dunefields. The age of the divergence event suggests that this may reflect the effect of older dispersal barriers, and that genetic structure was subsequently maintained by the formation of contemporary coastal dunefields.

The sardine run in southeastern Africa is a mass migration into an ecological trap

The KwaZulu-Natal sardine run, popularly known as the “greatest shoal on Earth,” is a mass migration of South African sardines from their temperate core range into the subtropical Indian Ocean. It has been suggested that this represents the spawning migration of a distinct subtropical stock. Using genomic and transcriptomic data from sardines collected around the South African coast, we identified two stocks, one cool temperate (Atlantic) and the other warm temperate (Indian Ocean). Unexpectedly, we found that sardines participating in the sardine run are primarily of Atlantic origin and thus prefer colder water. These sardines separate from the warm-temperate stock and move into temporarily favorable Indian Ocean habitat during brief cold-water upwelling periods. Once the upwelling ends, they find themselves trapped in physiologically challenging subtropical habitat and subject to intense predation pressure. This makes the sardine run a rare example of a mass migration that has no apparent fitness benefits.

Reevaluating claims of ecological speciation in Halichoeres bivittatus

Allopatry has traditionally been viewed as the primary driver of speciation in marine taxa, but the geography of the marine environment and the larval dispersal capabilities of many marine organisms render this view somewhat questionable. In marine fishes, one of the earliest and most highly cited empirical examples of ecological speciation with gene flow is the slippery dick wrasse, Halichoeres bivittatus. Evidence for this cryptic or incipient speciation event was primarily in the form of a deep divergence in a single mitochondrial locus between the northern and southern Gulf of Mexico, combined with a finding that these two haplotypes were associated with different habitat types ("tropical" vs. "subtropical") in the Florida Keys and Bermuda, where they overlap. Here, we examine habitat assortment in the Florida Keys using a broader sampling of populations and habitat types than were available for the original study. We find no evidence to support the claim that haplotype frequencies differ between habitat types, and little evidence to support any differences between populations in the Keys. These results undermine claims of ecological speciation with gene flow in Halichoeres bivittatus. Future claims of this type should be supported by multiple lines of evidence that illuminate potential mechanisms and allow researchers to rule out alternative explanations for spatial patterns of genetic differences.

Population structure associated with bioregion and seasonal prey distribution for Indo-Pacific bottlenose dolphins (Tursiops aduncus) in South Africa

Many marine species exhibit fine-scale population structure despite high mobility and a lack of physical barriers to dispersal, but the evolutionary drivers of differentiation in these systems are generally poorly understood. Here we investigate the potential role of habitat transitions and seasonal prey distributions on the evolution of population structure in the Indo-Pacific bottlenose dolphin, Tursiops aduncus, off South Africa's coast, using double-digest restriction-site associated DNA sequencing. Population structure was identified between the eastern and southern coasts and correlated with the habitat transition between the temperate Agulhas (southern) and subtropical Natal (eastern) Bioregions, suggesting differentiation driven by resource specializations. Differentiation along the Natal coast was comparatively weak, but was evident in some analyses and varied depending on whether the samples were collected during or outside the seasonal sardine (Sardinops sagax) run. This local abundance of prey could influence the ranging patterns and apparent genetic structure of T. aduncus. These findings have significant and transferable management implications, most importantly in terms of differentiating populations inhabiting distinct bioregions and seasonal structural patterns within a region associated with the movement of prey resources.

Genomics-informed models reveal extensive stretches of coastline under threat by an ecologically dominant invasive species

Explaining why some species are widespread, while others are not, is fundamental to biogeography, ecology, and evolutionary biology. A unique way to study evolutionary and ecological mechanisms that either limit species' spread or facilitate range expansions is to conduct research on species that have restricted distributions. Nonindigenous species, particularly those that are highly invasive but have not yet spread beyond the introduced site, represent ideal systems to study range size changes. Here, we used species distribution modeling and genomic data to study the restricted range of a highly invasive Australian marine species, the ascidian Pyura praeputialis This species is an aggressive space occupier in its introduced range (Chile), where it has fundamentally altered the coastal community. We found high genomic diversity in Chile, indicating high adaptive potential. In addition, genomic data clearly showed that a single region from Australia was the only donor of genotypes to the introduced range. We identified over 3,500 km of suitable habitat adjacent to its current introduced range that has so far not been occupied, and importantly species distribution models were only accurate when genomic data were considered. Our results suggest that a slight change in currents, or a change in shipping routes, may lead to an expansion of the species' introduced range that will encompass a vast portion of the South American coast. Our study shows how the use of population genomics and species distribution modeling in combination can unravel mechanisms shaping range sizes and forecast future range shifts of invasive species.

Whole-genome re-sequencing data to infer historical demography and speciation processes in land snails: the study of two Candidula sister species

Despite the global biodiversity of terrestrial gastropods and their ecological and economic importance, the genomic basis of ecological adaptation and speciation in land snail taxa is still largely unknown. Here, we combined whole-genome re-sequencing with population genomics to evaluate the historical demography and the speciation process of two closely related species of land snails from western Europe, Candidula unifasciata and C. rugosiuscula. Historical demographic analysis indicated fluctuations in the size of ancestral populations, probably driven by Pleistocene climatic fluctuations. Although the current population distributions of both species do not overlap, our approximate Bayesian computation model selection approach on several speciation scenarios suggested that gene flow has occurred throughout the divergence process until recently. Positively selected genes diverging early in the process were associated with intragenomic and cyto-nuclear incompatibilities, respectively, potentially fostering reproductive isolation as well as ecological divergence. Our results suggested that the speciation between species entails complex processes involving both gene flow and ecological speciation, and that further research based on whole-genome data can provide valuable understanding on species divergence. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Mitochondrial genome announcements need to consider existing short sequences from closely related species to prevent taxonomic errors

The reconstruction of complete mitochondrial genomes (mitogenomes) has considerable potential to clarify species relationships in cases where morphological analysis and DNA sequencing of individual genes are inconclusive. However, the trend to use only mitogenomes for the phylogenies presented in mitogenome announcements carries the inherent risk that the study species’ taxonomy is incorrect because no mitogenomes have yet been reconstructed for its sister species. Here, I illustrate this problem using the mitogenomes of two seahorses, Hippocampus capensis and H. queenslandicus. Both specimens used for mitogenome reconstruction originated from traditional Chinese medicine markets rather than native habitats. Although mitogenome phylogenies placed these specimens correctly among the seahorses from which mitogenomes were available at the time, incorporating single-marker sequence from closely related species into the phylogenies revealed that both mitogenomes are problematic. The mitogenome of the endemic South African H. capensis did not cluster among single-marker DNA sequences of seahorses from the species’ native habitat, but among sequences submitted under the names H. casscsio, H. fuscus and H. kuda that originated from all over the Indo-Pacific, including China. Phylogenetic placement of the mitogenome of H. queenslandicus within a cluster of seahorses that also included H. spinosissimus confirms an earlier finding that H. queenslandicus is a synonym of H. spinosissimus, a widespread Indo-Pacific species that also occurs in China. It is recommended that mitogenome announcements incorporate available single-marker sequences of closely related species, not only mitogenomes. The reconstruction of mitogenomes can exacerbate taxonomic confusion if existing information is ignored.

Genomic divergence and differential gene expression between crustacean ecotypes across a marine thermal gradient

Environmental gradients between marine biogeographical provinces separate distinct faunal communities. However, the absence of absolute dispersal barriers allows numerous species to occur on both sides of such boundaries. While the regional populations of such widespread species are often morphologically indistinguishable from each other, genetic evidence suggests that they represent unique ecotypes, and likely even cryptic species, that may be uniquely adapted to their local environment. Here, we explored genomic divergence in four sympatric southern African decapod crustaceans whose ranges span the boundary between the cool-temperate west coast (south-eastern Atlantic) and the warm-temperate south coast (south-western Indian Ocean) near the southern tip of the African continent. Using genome-wide data, we found that all four species comprise distinct west coast and south coast ecotypes, with molecular dating suggesting divergence during the Pleistocene. Transcriptomic data from the hepatopancreas of twelve specimens of one of these species, the mudprawn Upogebia africana, which were exposed to either 10 °C or 20 °C, showed a clear difference in gene expression profiles between the west- and south coast ecotypes. This difference was particularly clear at 10 °C, where individuals from the south coast experienced a 'transcriptomic shock'. This low temperature is more typical of the west coast during upwelling events, and the physiological stress experienced by the south coast ecotype under such conditions may explain its absence from that coastline. Our results shed new light on the processes involved in driving genomic divergence and incipient speciation along coastlines with porous dispersal barriers.

Phylogeography of Mytilisepta virgata (Mytilidae: Bivalvia) in the northwestern Pacific: Cryptic mitochondrial lineages and mito-nuclear discordance

The purplish bifurcate mussel Mytilisepta virgata is widely distributed and represents one of the major components of the intertidal community in the northwestern Pacific (NWP). Here, we characterized population genetic structure of NWP populations throughout nearly their whole distribution range using both mitochondrial (mtDNA cox1) and nuclear (ITS1) markers. Population genetic analyses for mtDNA cox 1 sequences revealed two monophyletic lineages (i.e., southern and northern lineages) geographically distributed according to the two different surface water temperature zones in the NWP. The timing of the lineage split is estimated at the Pliocene- mid-Pleistocene (5.49-1.61 Mya), which is consistent with the timing of the historical isolation of the East Sea/Sea of Japan from the South and East China Seas due to sea level decline during glacial cycles. Historical sea level fluctuation during the Pliocene-Pleistocene and subsequent adaptation of mussels to different surface water temperature zones may have contributed to shaping the contemporary genetic diversity and deep divergence of the two mitochondrial lineages. In contrast to mtDNA sequences, a clear lineage split between the two mitochondrial lineages was not found in ITS1 sequences, which showed a star-like structure composed of a mixture of southern and northern mitochondrial lineages. Possible reasons for this type of mito-nuclear discordance include stochastic divergence in the coalescent processes of the two molecular markers, or balancing selection under different marine environments. Cryptic speciation cannot be ruled out from these results, and future work using genomic analyses is required to address whether the thermal physiology of these mussels corresponds to the deep divergence of their mitochondrial genes and to test for the existence of morphologically indistinguishable but genetically separate cryptic species.

Multi-model seascape genomics identifies distinct environmental drivers of selection among sympatric marine species

BACKGROUND

As global change and anthropogenic pressures continue to increase, conservation and management increasingly needs to consider species' potential to adapt to novel environmental conditions. Therefore, it is imperative to characterise the main selective forces acting on ecosystems, and how these may influence the evolutionary potential of populations and species. Using a multi-model seascape genomics approach, we compare putative environmental drivers of selection in three sympatric southern African marine invertebrates with contrasting ecology and life histories: Cape urchin (Parechinus angulosus), Common shore crab (Cyclograpsus punctatus), and Granular limpet (Scutellastra granularis).

RESULTS

Using pooled (Pool-seq), restriction-site associated DNA sequencing (RAD-seq), and seven outlier detection methods, we characterise genomic variation between populations along a strong biogeographical gradient. Of the three species, only S. granularis showed significant isolation-by-distance, and isolation-by-environment driven by sea surface temperatures (SST). In contrast, sea surface salinity (SSS) and range in air temperature correlated more strongly with genomic variation in C. punctatus and P. angulosus. Differences were also found in genomic structuring between the three species, with outlier loci contributing to two clusters in the East and West Coasts for S. granularis and P. angulosus, but not for C. punctatus.

CONCLUSION

The findings illustrate distinct evolutionary potential across species, suggesting that species-specific habitat requirements and responses to environmental stresses may be better predictors of evolutionary patterns than the strong environmental gradients within the region. We also found large discrepancies between outlier detection methodologies, and thus offer a novel multi-model approach to identifying the principal environmental selection forces acting on species. Overall, this work highlights how adding a comparative approach to seascape genomics (both with multiple models and species) can elucidate the intricate evolutionary responses of ecosystems to global change.

Signals of selection in the mitogenome provide insights into adaptation mechanisms in heterogeneous habitats in a widely distributed pelagic fish

Oceans are vast, dynamic, and complex ecosystems characterized by fluctuations in environmental parameters like sea surface temperature (SST), salinity, oxygen availability, and productivity. Environmental variability acts as the driver of organismal evolution and speciation as organisms strive to cope with the challenges. We investigated the evolutionary consequences of heterogeneous environmental conditions on the mitogenome of a widely distributed small pelagic fish of Indian ocean, Indian oil sardine, Sardinella longiceps. Sardines were collected from different eco-regions of the Indian Ocean and selection patterns analyzed in coding and non-coding regions. Signals of diversifying selection were observed in key functional regions involved in OXPHOS indicating OXPHOS gene regulation as the critical factor to meet enhanced energetic demands. A characteristic control region with 38–40 bp tandem repeat units under strong selective pressure as evidenced by sequence conservation and low free energy values was also observed. These changes were prevalent in fishes from the South Eastern Arabian Sea (SEAS) followed by the Northern Arabian Sea (NAS) and rare in Bay of Bengal (BoB) populations. Fishes belonging to SEAS exhibited accelerated substitution rate mainly due to the selective pressures to survive in a highly variable oceanic environment characterized by seasonal hypoxia, variable SST, and food availability.

Trans-Atlantic Distribution and Introgression as Inferred from Single Nucleotide Polymorphism: Mussels Mytilus and Environmental Factors

Large-scale climate changes influence the geographic distribution of biodiversity. Many taxa have been reported to extend or reduce their geographic range, move poleward or displace other species. However, for closely related species that can hybridize in the natural environment, displacement is not the only effect of changes of environmental variables. Another option is subtler, hidden expansion, which can be found using genetic methods only. The marine blue mussels Mytilus are known to change their geographic distribution despite being sessile animals. In addition to natural dissemination at larval phase—enhanced by intentional or accidental introductions and rafting—they can spread through hybridization and introgression with local congeners, which can create mixed populations sustaining in environmental conditions that are marginal for pure taxa. The Mytilus species have a wide distribution in coastal regions of the Northern and Southern Hemisphere. In this study, we investigated the inter-regional genetic differentiation of the Mytilus species complex at 53 locations in the North Atlantic and adjacent Arctic waters and linked this genetic variability to key local environmental drivers. Of seventy-nine candidate single nucleotide polymorphisms (SNPs), all samples were successfully genotyped with a subset of 54 SNPs. There was a clear interregional separation of Mytilus species. However, all three Mytilus species hybridized in the contact area and created hybrid zones with mixed populations. Boosted regression trees (BRT) models showed that inter-regional variability was important in many allele models but did not prevail over variability in local environmental factors. Local environmental variables described over 40% of variability in about 30% of the allele frequencies of Mytilus spp. For the 30% of alleles, variability in their frequencies was only weakly coupled with local environmental conditions. For most studied alleles the linkages between environmental drivers and the genetic variability of Mytilus spp. were random in respect to “coding” and “non-coding” regions. An analysis of the subset of data involving functional genes only showed that two SNPs at Hsp70 and ATPase genes correlated with environmental variables. Total predictive ability of the highest performing models (r² between 0.550 and 0.801) were for alleles that discriminated most effectively M. trossulus from M. edulis and M. galloprovincialis, whereas the best performing allele model (BM101A) did the best at discriminating M. galloprovincialis from M. edulis and M. trossulus. Among the local environmental variables, salinity, water temperature, ice cover and chlorophyll a concentration were by far the greatest predictors, but their predictive performance varied among different allele models. In most cases changes in the allele frequencies along these environmental gradients were abrupt and occurred at a very narrow range of environmental variables. In general, regions of change in allele frequencies for M. trossulus occurred at 8–11 psu, 0–10 °C, 60%–70% of ice cover and 0–2 mg m⁻³ of chlorophyll a, M. edulis at 8–11 and 30–35 psu, 10–14 °C and 60%–70% of ice cover and for M. galloprovincialis at 30–35 psu, 14–20 °C.

Fine-Scale Population Genetic Structure and Parapatric Cryptic Species of Kuruma Shrimp (Marsupenaeus japonicus), Along the Northwestern Pacific Coast of China

The kuruma shrimp (Marsupenaeus japonicus) includes two cryptic species, which are distributed mostly allopatrically but co-occur in the northern South China Sea (from Huilai to Beihai). To obtain a better understanding of the fine-scale genetic structure and parapatric diversification of these two varieties in the northwestern Pacific region, we used a genotyping-by-sequencing (GBS) and comparative transcriptomics approach to establish their phylogenetic relationships. Using the GBS technique, we genotyped 28891 SNPs in 160 individuals in the Northwest Pacific. The results supported two highly diverged evolutionary lineages of kuruma shrimp (var. I and II). The ND and XM populations showed complex genetic patterns, which might be affected by the complex environment of the Taiwan Strait. In addition, the migration rates and inbreeding coefficients of XM and BH were much lower than those of the other populations, which might be related to the land-sea changes and complex ocean currents in the Taiwan Strait and Qiongzhou Strait. Based on the synonymous substitution rates (ds) of 2,491 candidate orthologs, we estimated that the divergence time between the two varieties was 0.26~0.69 Mya. Choice and no-choice interbreeding experiments provided support for the biological species concept, by showing the existence of reproductive isolation or incompatibility. In view of these differences between the two Marsupenaeus species, we believe that it is essential and urgent to establish a genetic database for each and reevaluate their ecological suitable conditions in order to improve species-specific culturing techniques. Moreover, this research can serve as a case study for future research on speciation and hybridization.

Rejection of the genetic implications of the "Abundant Centre Hypothesis" in marine mussels

The ‘Abundant-Centre Hypothesis’ is a well-established but controversial hypothesis stating that the abundance of a species is highest at the centre of its range and decreases towards the edges, where conditions are unfavourable. As genetic diversity depends on population size, edge populations are expected to show lower intra-population genetic diversity than core populations, while showing high inter-population genetic divergence. Here, the genetic implications of the Abundant-Centre Hypothesis were tested on two coastal mussels from South Africa that disperse by means of planktonic larvae, the native Perna perna and the invasive Mytilus galloprovincialis. Genetic structure was found within P. perna, which, together with evidence from Lagrangian particle simulations, points to significant reductions in gene flow between sites. Despite this, the expected diversity pattern between centre and edge populations was not found for either species. We conclude that the genetic predictions of the Abundant-Centre Hypothesis are unlikely to be met by high-dispersal species with large population sizes, and may only become evident in species with much lower levels of connectivity.

Ecological divergence and speciation in common bottlenose dolphins in the western South Atlantic

Coastal and offshore ecotypes of common bottlenose dolphins have been recognized in the western South Atlantic, and it is possible that trophic niche divergence associated with social interactions is leading them to genetic and phenotypic differentiation. The significant morphological differentiation observed between these ecotypes suggests they represent two different subspecies. However, there is still a need to investigate whether there is congruence between morphological and genetic data to rule out the possibility of ecophenotypic variation accompanied by gene flow. Mitochondrial DNA (mtDNA) control region sequence data and 10 microsatellite loci collected from stranded and biopsied dolphins sampled in coastal and offshore waters of Brazil as well as 106 skulls for morphological analyses were used to determine whether the morphological differentiation was supported by genetic differentiation. There was congruence among the data sets, reinforcing the presence of two distinct ecotypes. The divergence may be relatively recent, however, given the moderate values of mtDNA nucleotide divergence (dA = 0.008), presence of one shared mtDNA haplotype and possibly low levels of gene flow (around 1% of migrants per generation). Results suggest the ecotypes may be in the process of speciation and reinforce they are best described as two different subspecies until the degree of nuclear genetic divergence is thoroughly evaluated: Tursiops truncatus gephyreus (coastal ecotype) and T. t. truncatus (offshore ecotype). The endemic distribution of T. t. gephyreus in the western South Atlantic and number of anthropogenic threats in the area reinforces the importance of protecting this ecotype and its habitat.

Genomic footprints of a biological invasion: Introduction from Asia and dispersal in Europe of the topmouth gudgeon (Pseudorasbora parva)

Facilitated by the intensification of global trading, the introduction and dispersal of species to areas in which they are historically non-native is nowadays common. From an evolutionary standpoint, invasions are paradoxical: not only non-native environments could be different from native ones for which introduced individuals would be ill-adapted, but also small founding population size should be associated with reduced adaptive potential. As such, biological invasions are considered valuable real-time evolutionary experiments. Here, we investigated the population structure and adaptive potential of the highly invasive topmouth gudgeon (Pseudorasbora parva) across Europe and East Asia. We RAD-sequenced 301 specimens from sixteen populations and three distinct within-catchment invaded regions as well as two locations in the native range. With 13,785 single nucleotide polymorphisms, we provide conclusive evidence for a genome-wide signature of two distinct invasion events, in Slovakia and Turkey, each originating from a specific area in the native range. A third invaded area, in France, appears to be the result of dispersal within the invasive range. Few loci showed signs of selection, the vast majority of which being identified in the Slovakian region. Functional annotation suggests that faster early stage development, resistance to pollution and immunocompetence contribute to the invasion success of the local habitats. By showing that populations in the invasive range have different evolutionary histories, our study reinforces the idea that populations, rather than species, are the units to consider in invasion biology.