No evidence for a role of MHC class II genotype in the chemical encoding of heterozygosity and relatedness in Antarctic fur seals

Despite decades of research, surprisingly little is known about the mechanism(s) by which an individual's genotype is encoded in odour. Many studies have focused on the role of the major histocompatibility complex (MHC) owing to its importance for survival and mate choice. However, the salience of MHC-mediated odours compared to chemicals influenced by the rest of the genome remains unclear, especially in wild populations where it is challenging to quantify and control for the effects of the genomic background. We addressed this issue in Antarctic fur seals by analysing skin swabs together with full-length MHC DQB II exon 2 sequences and data from 41 genome-wide distributed microsatellites. We did not find any effects of MHC relatedness on chemical similarity and there was also no relationship between MHC heterozygosity and chemical diversity. However, multilocus heterozygosity showed a significant positive association with chemical diversity, even after controlling for MHC heterozygosity. Our results appear to rule out a dominant role of the MHC in the chemical encoding of genetic information in a wild vertebrate population and highlight the need for genome-wide approaches to elucidate the mechanism(s) and specific genes underlying genotype-odour associations.

Diet in phenotypically divergent sympatric species of African weakly electric fish (genus: Campylomormyrus)-A hybrid capture/HTS metabarcoding approach

Ecological speciation within the mormyrid genus Campylomormyrus resulted in sympatric species exhibiting divergence in their feeding apparatus and electric organ discharge (EOD). This study documents the overall diet of the genus Campylomormyrus and examines the hypothesis that the Campylomormyrus radiation is caused by adaptation to different food sources. We performed diet assessment of five sympatric Campylomormyrus species (C. alces, C. compressirostris, C. curvirostris, C. tshokwe, C. numenius) and their sister taxon Gnathonemus petersii with markedly different snout morphologies and EODs using hybrid capture/HTS DNA metabarcoding of their stomach contents. Our approach allowed for high taxonomic resolution of prey items, including benthic invertebrates, allochthonous invertebrates and vegetation. Comparisons of the diet compositions using quantitative measures and diet overlap indices revealed that all species are able to exploit multiple food niches in their habitats, that is fauna at the bottom, the water surface and the water column. A major part of the diet is larvae of aquatic insects, such as dipterans, coleopterans and trichopterans, known to occur in holes and interstitial spaces of the substrate. The results indicate that different snout morphologies and the associated divergence in the EOD could translate into different prey spectra. This suggests that the diversification in EOD and/or morphology of the feeding apparatus could be under functional adaptation.

Gene and Allele-Specific Expression Underlying the Electric Signal Divergence in African Weakly Electric Fish

In the African weakly electric fish genus Campylomormyrus, electric organ discharge signals are strikingly different in shape and duration among closely related species, contribute to prezygotic isolation, and may have triggered an adaptive radiation. We performed mRNA sequencing on electric organs and skeletal muscles (from which the electric organs derive) from 3 species with short (0.4 ms), medium (5 ms), and long (40 ms) electric organ discharges and 2 different cross-species hybrids. We identified 1,444 upregulated genes in electric organ shared by all 5 species/hybrid cohorts, rendering them candidate genes for electric organ-specific properties in Campylomormyrus. We further identified several candidate genes, including KCNJ2 and KLF5, and their upregulation may contribute to increased electric organ discharge duration. Hybrids between a short (Campylomormyrus compressirostris) and a long (Campylomormyrus rhynchophorus) discharging species exhibit electric organ discharges of intermediate duration and showed imbalanced expression of KCNJ2 alleles, pointing toward a cis-regulatory difference at this locus, relative to electric organ discharge duration. KLF5 is a transcription factor potentially balancing potassium channel gene expression, a crucial process for the formation of an electric organ discharge. Unraveling the genetic basis of the species-specific modulation of the electric organ discharge in Campylomormyrus is crucial for understanding the adaptive radiation of this emerging model taxon of ecological (perhaps even sympatric) speciation.

DNA metabarcoding reveals impact of local recruitment, dispersal, and hydroperiod on assembly of a zooplankton metacommunity

Understanding the environmental impact on the assembly of local communities in relation to their spatial and temporal connectivity is still a challenge in metacommunity ecology. This study aims to unravel underlying metacommunity processes and environmental factors that result in observed zooplankton communities. Unlike most metacommunity studies, we jointly examine active and dormant zooplankton communities using a DNA metabarcoding approach to overcome limitations of morphological species identification. We applied two-fragment (COI and 18S) metabarcoding to monitor communities of 24 kettle holes over a two-year period to unravel (i) spatial and temporal connectivity of the communities, (ii) environmental factors influencing local communities, and (iii) dominant underlying metacommunity processes in this system. We found a strong separation of zooplankton communities from kettle holes of different hydroperiods (degree of permanency) throughout the season, while the community composition within single kettle holes did not differ between years. Species richness was primarily dependent on pH and permanency, while species diversity (Shannon Index) was influenced by kettle hole location. Community composition was impacted by kettle hole size and surrounding field crops. Environmental processes dominated temporal and spatial processes. Sediment communities showed a different composition compared to water samples but did not differ between ephemeral and permanent kettle holes. Our results suggest that communities are mainly structured by environmental filtering based on pH, kettle hole size, surrounding field crops, and permanency. Environmental filtering based on specific conditions in individual kettle holes seems to be the dominant process in community assembly in the studied zooplankton metacommunity.

Evolutionary history and seascape genomics of Harbour porpoises (Phocoena phocoena) across environmental gradients in the North Atlantic and adjacent waters

The Harbour porpoise (Phocoena phocoena) is a highly mobile cetacean species primarily occurring in coastal and shelf waters across the Northern hemisphere. It inhabits heterogeneous seascapes broadly varying in salinity and temperature. Here, we produced 74 whole genomes at intermediate coverage to study Harbour porpoise's evolutionary history and investigate the role of local adaptation in the diversification into subspecies and populations. We identified ~6 million high quality SNPs sampled at eight localities across the North Atlantic and adjacent waters, which we used for population structure, demographic and genotype-environment association analyses. Our results suggest a genetic differentiation between three subspecies (P.p. relicta, P.p. phocoena and P.p. meridionalis), and three distinct populations within P.p. phocoena: Atlantic, Belt Sea and Proper Baltic Sea. Effective population size and Tajima's D suggest population contraction in Black Sea and Iberian porpoises, but expansion in the P.p. phocoena populations. Phylogenetic trees indicate post-glacial colonization from a southern refugium. Genotype-environment association analysis identified salinity as major driver in genomic variation and we identified candidate genes putatively underlying adaptation to different salinity. Our study highlights the value of whole genome resequencing to unravel subtle population structure in highly mobile species, shows how strong environmental gradients and local adaptation may lead to population differentiation, and how neutral and adaptive markers can give different perspectives on population subdivision. The results have great conservation implications as we found inbreeding and low genetic diversity in the endangered Black Sea subspecies and identified the critically endangered Proper Baltic Sea porpoises as a separate population.

Recruitment and migration patterns reveal a key role for seed banks in the meta-population dynamics of an aquatic plant

Progressive habitat fragmentation threatens plant species with narrow habitat requirements. While local environmental conditions define population growth rates and recruitment success at the patch level, dispersal is critical for population viability at the landscape scale. Identifying the dynamics of plant meta-populations is often confounded by the uncertainty about soil-stored population compartments. We combined a landscape-scale assessment of an amphibious plant's population structure with measurements of dispersal complexity in time to track dispersal and putative shifts in functional connectivity. Using 13 microsatellite markers, we analyzed the genetic structure of extant Oenanthe aquatica populations and their soil seed banks in a kettle hole system to uncover hidden connectivity among populations in time and space. Considerable spatial genetic structure and isolation-by-distance suggest limited gene flow between sites. Spatial isolation and patch size showed minor effects on genetic diversity. Genetic similarity found among extant populations and their seed banks suggests increased local recruitment, despite some evidence of migration and recent colonization. Results indicate stepping-stone dispersal across adjacent populations. Among permanent and ephemeral demes the resulting meta-population demography could be determined by source-sink dynamics. Overall, these spatiotemporal connectivity patterns support mainland-island dynamics in our system, highlighting the importance of persistent seed banks as enduring sources of genetic diversity.

A new genome assembly of an African weakly electric fish (Campylomormyrus compressirostris, Mormyridae) indicates rapid gene family evolution in Osteoglossomorpha

BACKGROUND

Teleost fishes comprise more than half of the vertebrate species. Within teleosts, most phylogenies consider the split between Osteoglossomorpha and Euteleosteomorpha/Otomorpha as basal, preceded only by the derivation of the most primitive group of teleosts, the Elopomorpha. While Osteoglossomorpha are generally species poor, the taxon contains the African weakly electric fish (Mormyroidei), which have radiated into numerous species. Within the mormyrids, the genus Campylomormyrus is mostly endemic to the Congo Basin. Campylomormyrus serves as a model to understand mechanisms of adaptive radiation and ecological speciation, especially with regard to its highly diverse species-specific electric organ discharges (EOD). Currently, there are few well-annotated genomes available for electric fish in general and mormyrids in particular. Our study aims at producing a high-quality genome assembly and to use this to examine genome evolution in relation to other teleosts. This will facilitate further understanding of the evolution of the osteoglossomorpha fish in general and of electric fish in particular.

RESULTS

A high-quality weakly electric fish (C. compressirostris) genome was produced from a single individual with a genome size of 862 Mb, consisting of 1,497 contigs with an N50 of 1,399 kb and a GC-content of 43.69%. Gene predictions identified 34,492 protein-coding genes, which is a higher number than in the two other available Osteoglossomorpha genomes of Paramormyrops kingsleyae and Scleropages formosus. A Computational Analysis of gene Family Evolution (CAFE5) comparing 33 teleost fish genomes suggests an overall faster gene family turnover rate in Osteoglossomorpha than in Otomorpha and Euteleosteomorpha. Moreover, the ratios of expanded/contracted gene family numbers in Osteoglossomorpha are significantly higher than in the other two taxa, except for species that had undergone an additional genome duplication (Cyprinus carpio and Oncorhynchus mykiss). As potassium channel proteins are hypothesized to play a key role in EOD diversity among species, we put a special focus on them, and manually curated 16 Kv1 genes. We identified a tandem duplication in the KCNA7a gene in the genome of C. compressirostris.

CONCLUSIONS

We present the fourth genome of an electric fish and the third well-annotated genome for Osteoglossomorpha, enabling us to compare gene family evolution among major teleost lineages. Osteoglossomorpha appear to exhibit rapid gene family evolution, with more gene family expansions than contractions. The curated Kv1 gene family showed seven gene clusters, which is more than in other analyzed fish genomes outside Osteoglossomorpha. The KCNA7a, encoding for a potassium channel central for EOD production and modulation, is tandemly duplicated which may related to the diverse EOD observed among Campylomormyrus species.

Variation in heat shock protein 40 kDa relates to divergence in thermotolerance among cryptic rotifer species

Genetic divergence and the frequency of hybridization are central for defining species delimitations, especially among cryptic species where morphological differences are merely absent. Rotifers are known for their high cryptic diversity and therefore are ideal model organisms to investigate such patterns. Here, we used the recently resolved Brachionus calyciflorus species complex to investigate whether previously observed between species differences in thermotolerance and gene expression are also reflected in their genomic footprint. We identified a Heat Shock Protein gene (HSP 40 kDa) which exhibits cross species pronounced sequence variation. This gene exhibits species-specific fixed sites, alleles, and sites putatively under positive selection. These sites are located in protein binding regions involved in chaperoning and may therefore reflect adaptive diversification. By comparing three genetic markers (ITS, COI, HSP 40 kDa), we revealed hybridization events between the cryptic species. The low frequency of introgressive haplotypes/alleles suggest a tight, but not fully impermeable boundary between the cryptic species.

Pleistocene climate variability in eastern Africa influenced hominin evolution

Despite more than half a century of hominin fossil discoveries in eastern Africa, the regional environmental context of hominin evolution and dispersal is not well established due to the lack of continuous palaeoenvironmental records from one of the proven habitats of early human populations, particularly for the Pleistocene epoch. Here we present a 620,000-year environmental record from Chew Bahir, southern Ethiopia, which is proximal to key fossil sites. Our record documents the potential influence of different episodes of climatic variability on hominin biological and cultural transformation. The appearance of high anatomical diversity in hominin groups coincides with long-lasting and relatively stable humid conditions from ~620,000 to 275,000 years bp (episodes 1-6), interrupted by several abrupt and extreme hydroclimate perturbations. A pattern of pronounced climatic cyclicity transformed habitats during episodes 7-9 (~275,000-60,000 years bp), a crucial phase encompassing the gradual transition from Acheulean to Middle Stone Age technologies, the emergence of Homo sapiens in eastern Africa and key human social and cultural innovations. Those accumulative innovations plus the alignment of humid pulses between northeastern Africa and the eastern Mediterranean during high-frequency climate oscillations of episodes 10-12 (~60,000-10,000 years bp) could have facilitated the global dispersal of H. sapiens.

Seed traits matter-Endozoochoric dispersal through a pervasive mobile linker

Although many plants are dispersed by wind and seeds can travel long distances across unsuitable matrix areas, a large proportion relies on co-evolved zoochorous seed dispersal to connect populations in isolated habitat islands. Particularly in agricultural landscapes, where remaining habitat patches are often very small and highly isolated, mobile linkers as zoochorous seed dispersers are critical for the population dynamics of numerous plant species. However, knowledge about the quali- or quantification of such mobile link processes, especially in agricultural landscapes, is still limited. In a controlled feeding experiment, we recorded the seed intake and germination success after complete digestion by the European brown hare (Lepus europaeus) and explored its mobile link potential as an endozoochoric seed disperser. Utilizing a suite of common, rare, and potentially invasive plant species, we disentangled the effects of seed morphological traits on germination success while controlling for phylogenetic relatedness. Further, we measured the landscape connectivity via hares in two contrasting agricultural landscapes (simple: few natural and semi-natural structures, large fields; complex: high amount of natural and semi-natural structures, small fields) using GPS-based movement data. With 34,710 seeds of 44 plant species fed, one of 200 seeds (0.51%) with seedlings of 33 species germinated from feces. Germination after complete digestion was positively related to denser seeds with comparatively small surface area and a relatively slender and elongated shape, suggesting that, for hares, the most critical seed characteristics for successful endozoochorous seed dispersal minimize exposure of the seed to the stomach and the associated digestive system. Furthermore, we could show that a hare's retention time is long enough to interconnect different habitats, especially grasslands and fields. Thus, besides other seed dispersal mechanisms, this most likely allows hares to act as effective mobile linkers contributing to ecosystem stability in times of agricultural intensification, not only in complex but also in simple landscapes.

A phylogeny of the genus Limia (Teleostei: Poeciliidae) suggests a single-lake radiation nested in a Caribbean-wide allopatric speciation scenario

OBJECTIVE

The Caribbean is an important global biodiversity hotspot. Adaptive radiations there lead to many speciation events within a limited period and hence are particularly prominent biodiversity generators. A prime example are freshwater fish of the genus Limia, endemic to the Greater Antilles. Within Hispaniola, nine species have been described from a single isolated site, Lake Miragoâne, pointing towards extraordinary sympatric speciation. This study examines the evolutionary history of the Limia species in Lake Miragoâne, relative to their congeners throughout the Caribbean.

RESULTS

For 12 Limia species, we obtained almost complete sequences of the mitochondrial cytochrome b gene, a well-established marker for lower-level taxonomic relationships. We included sequences of six further Limia species from GenBank (total N  = 18 species). Our phylogenies are in concordance with other published phylogenies of Limia. There is strong support that the species found in Lake Miragoâne in Haiti are monophyletic, confirming a recent local radiation. Within Lake Miragoâne, speciation is likely extremely recent, leading to incomplete lineage sorting in the mtDNA. Future studies using multiple unlinked genetic markers are needed to disentangle the relationships within the Lake Miragoâne clade.

Southern introgression increases adaptive immune gene variability in northern range margin populations of Fire-bellied toad

Northern range margin populations of the European fire-bellied toad (Bombina bombina) have rapidly declined during recent decades. Extensive agricultural land use has fragmented the landscape, leading to habitat disruption and loss, as well as eutrophication of ponds. In Northern Germany (Schleswig-Holstein) and Southern Sweden (Skåne), this population decline resulted in decreased gene flow from surrounding populations, low genetic diversity, and a putative reduction in adaptive potential, leaving populations vulnerable to future environmental and climatic changes. Previous studies using mitochondrial control region and nuclear transcriptome-wide SNP data detected introgressive hybridization in multiple northern B. bombina populations after unreported release of toads from Austria. Here, we determine the impact of this introgression by comparing the body conditions (proxy for fitness) of introgressed and nonintrogressed populations and the genetic consequences in two candidate genes for putative local adaptation (the MHC II gene as part of the adaptive immune system and the stress response gene HSP70 kDa). We detected regional differences in body condition and observed significantly elevated levels of within individual MHC allele counts in introgressed Swedish populations, associated with a tendency toward higher body weight, relative to regional nonintrogressed populations. These differences were not observed among introgressed and nonintrogressed German populations. Genetic diversity in both MHC and HSP was generally lower in northern than Austrian populations. Our study sheds light on the potential benefits of translocations of more distantly related conspecifics as a means to increase adaptive genetic variability and fitness of genetically depauperate range margin populations without distortion of local adaptation.

Genomic consequences of human-mediated translocations in margin populations of an endangered amphibian

Due to their isolated and often fragmented nature, range margin populations are especially vulnerable to rapid environmental change. To maintain genetic diversity and adaptive potential, gene flow from disjunct populations might therefore be crucial to their survival. Translocations are often proposed as a mitigation strategy to increase genetic diversity in threatened populations. However, this also includes the risk of losing locally adapted alleles through genetic swamping. Human-mediated translocations of southern lineage specimens into northern German populations of the endangered European fire-bellied toad (Bombina bombina) provide an unexpected experimental set-up to test the genetic consequences of an intraspecific introgression from central population individuals into populations at the species range margin. Here, we utilize complete mitochondrial genomes and transcriptome nuclear data to reveal the full genetic extent of this translocation and the consequences it may have for these populations. We uncover signs of introgression in four out of the five northern populations investigated, including a number of introgressed alleles ubiquitous in all recipient populations, suggesting a possible adaptive advantage. Introgressed alleles dominate at the MTCH2 locus, associated with obesity/fat tissue in humans, and the DSP locus, essential for the proper development of epidermal skin in amphibians. Furthermore, we found loci where local alleles were retained in the introgressed populations, suggesting their relevance for local adaptation. Finally, comparisons of genetic diversity between introgressed and nonintrogressed northern German populations revealed an increase in genetic diversity in all German individuals belonging to introgressed populations, supporting the idea of a beneficial transfer of genetic variation from Austria into North Germany.

Adaptive and nonadaptive plasticity in changing environments: Implications for sexual species with different life history strategies

Populations adapt to novel environmental conditions by genetic changes or phenotypic plasticity. Plastic responses are generally faster and can buffer fitness losses under variable conditions. Plasticity is typically modeled as random noise and linear reaction norms that assume simple one-to-one genotype-phenotype maps and no limits to the phenotypic response. Most studies on plasticity have focused on its effect on population viability. However, it is not clear, whether the advantage of plasticity depends solely on environmental fluctuations or also on the genetic and demographic properties (life histories) of populations. Here we present an individual-based model and study the relative importance of adaptive and nonadaptive plasticity for populations of sexual species with different life histories experiencing directional stochastic climate change. Environmental fluctuations were simulated using differentially autocorrelated climatic stochasticity or noise color, and scenarios of directional climate change. Nonadaptive plasticity was simulated as a random environmental effect on trait development, while adaptive plasticity as a linear, saturating, or sinusoidal reaction norm. The last two imposed limits to the plastic response and emphasized flexible interactions of the genotype with the environment. Interestingly, this assumption led to (a) smaller phenotypic than genotypic variance in the population (many-to-one genotype-phenotype map) and the coexistence of polymorphisms, and (b) the maintenance of higher genetic variation-compared to linear reaction norms and genetic determinism-even when the population was exposed to a constant environment for several generations. Limits to plasticity led to genetic accommodation, when costs were negligible, and to the appearance of cryptic variation when limits were exceeded. We found that adaptive plasticity promoted population persistence under red environmental noise and was particularly important for life histories with low fecundity. Populations producing more offspring could cope with environmental fluctuations solely by genetic changes or random plasticity, unless environmental change was too fast.

Ontogeny of the electric organ discharge and of the papillae of the electrocytes in the weakly electric fish Campylomormyrus rhynchophorus (Teleostei: Mormyridae)

The electric organ of the mormyrid weakly electric fish, Campylomormyrus rhynchophorus (Boulenger, 1898), undergoes changes in both the electric organ discharge (EOD) and the light and electron microscopic morphology as the fish mature from the juvenile to the adult form. Of particular interest was the appearance of papillae, surface specializations of the uninnervated anterior face of the electrocyte, which have been hypothesized to increase the duration of the EOD. In a 24.5 mm long juvenile the adult electric organ (EO) was not yet functional, and the electrocytes lacked papillae. A 40 mm long juvenile, which produced a short biphasic EOD of 1.3 ms duration, shows small papillae (average area 136 μm² ). In contrast, the EOD of a 79 mm long juvenile was triphasic. The large increase in duration of the EOD to 23.2 ms was accompanied by a small change in size of the papillae (average area 159 μm² ). Similarly, a 150 mm long adult produced a triphasic EOD of comparable duration to the younger stage (24.7 ms) but featured a prominent increase in size of the papillae (average area 402 μm² ). Thus, there was no linear correlation between EOD duration and papillary size. The most prominent ultrastructural change was at the level of the myofilaments, which regularly extended into the papillae, only in the oldest specimen-probably serving a supporting function. Physiological mechanisms, like gene expression levels, as demonstrated in some Campylomormyrus species, might be more important concerning the duration of the EOD.

Genetic Diversity and Connectivity in Plant Species Differing in Clonality and Dispersal Mechanisms in Wetland Island Habitats

In plants, long-distance dispersal is both attenuated and directed by specific movement vectors, including animals, wind, and/or water. Hence, movement vectors partly shape metapopulation genetic patterns that are, however, also influenced by other life-history traits such as clonal growth. We studied the relationship between area, isolation, plant-species richness, reproduction, and dispersal mechanisms with genetic diversity and divergence in 4 widespread wetland plant-species in a total of 20 island-like kettle-hole habitats surrounded by an intensive agricultural landscape. Our results showed that genetic parameters reflect the reproduction strategies with the highest genetic diversity being observed in the non-clonal, outcrossing Oenanthe aquatica compared to the clonal Lycopus europaeus, Typha latifolia, and Phragmites australis. Lycopus showed a positive relationship between genetic diversity and kettle-hole area, but a negative relationship with the number of neighboring kettle holes (less isolation). Genetic diversity increased with plant-species richness in the clonal species Phragmites and Lycopus; while it decreased in the non-clonal Oenanthe. Finally, genetic divergence and, therefore, connectivity differed between alternative dispersal strategies, where wind-dispersed Typha and Phragmites had a higher gene flow between the analyzed kettle holes compared with the insect-pollinated, hydrochorous Lycopus and Oenanthe. Our study provides information on genetic patterns related to reproduction and dispersal mechanisms of 4 common wetland species contributing to the understanding of the functioning of plant metacommunities occurring in kettle holes embedded in agricultural landscapes.

Transcriptome-wide single nucleotide polymorphisms related to electric organ discharge differentiation among African weakly electric fish species

African weakly electric fish of the mormyrid genus Campylomormyrus generate pulse-type electric organ discharges (EODs) for orientation and communication. Their pulse durations are species-specific and elongated EODs are a derived trait. So far, differential gene expression among tissue-specific transcriptomes across species with different pulses and point mutations in single ion channel genes indicate a relation of pulse duration and electrocyte geometry/excitability. However, a comprehensive assessment of expressed Single Nucleotide Polymorphisms (SNPs) throughout the entire transcriptome of African weakly electric fish, with the potential to identify further genes influencing EOD duration, is still lacking. This is of particular value, as discharge duration is likely based on multiple cellular mechanisms and various genes. Here we provide the first transcriptome-wide SNP analysis of African weakly electric fish species (genus Campylomormyrus) differing by EOD duration to identify candidate genes and cellular mechanisms potentially involved in the determination of an elongated discharge of C. tshokwe. Non-synonymous substitutions specific to C. tshokwe were found in 27 candidate genes with inferred positive selection among Campylomormyrus species. These candidate genes had mainly functions linked to transcriptional regulation, cell proliferation and cell differentiation. Further, by comparing gene annotations between C. compressirostris (ancestral short EOD) and C. tshokwe (derived elongated EOD), we identified 27 GO terms and 2 KEGG pathway categories for which C. tshokwe significantly more frequently exhibited a species-specific expressed substitution than C. compressirostris. The results indicate that transcriptional regulation as well cell proliferation and differentiation take part in the determination of elongated pulse durations in C. tshokwe. Those cellular processes are pivotal for tissue morphogenesis and might determine the shape of electric organs supporting the observed correlation between electrocyte geometry/tissue structure and discharge duration. The inferred expressed SNPs and their functional implications are a valuable resource for future investigations on EOD durations.

Indication of lethal interactions between a solitary bottlenose dolphin (Tursiops truncatus) and harbor porpoises (Phocoena phocoena) in the German Baltic Sea

Background

Aggressive interactions between bottlenose dolphins (Tursiops truncatus) and harbor porpoises (Phocoena phocoena) have been reported in different parts of the world since the late 1990s. In the Baltic Sea, harbor porpoises are the only native cetacean species, while bottlenose dolphins may appear there temporarily. In the fall of 2016, a solitary male photo-identified bottlenose dolphin stayed in the German Baltic Sea of Schleswig-Holstein for 3 months. During that time, the necropsies of the stranded harbor porpoises revealed types of trauma of varying degrees in six animals, which is unusual in this area. The purpose of this study was to determine if the appearance of the bottlenose dolphin could be linked to the trauma of the harbor porpoise carcasses.

Results

Pathological findings in these animals included subcutaneous, thoracic and abdominal hemorrhages, multiple, mainly bilateral, rib fractures, and one instance of lung laceration. These findings correspond with the previously reported dolphin-caused injuries in other regions. Moreover, public sighting reports showed a spatial and temporal correlation between the appearance of the dolphin and the stranding of fatally injured harbor porpoises.

Conclusion

Despite the fact that no attack has been witnessed in German waters to date, our findings indicate the first record of lethal interactions between a bottlenose dolphin and harbor porpoises in the German Baltic Sea. Furthermore, to our knowledge, this is the first report of porpoise aggression by a socially isolated bottlenose dolphin.

Intragenus F1-hybrids of African weakly electric fish (Mormyridae: Campylomormyrus tamandua ♂ × C. compressirostris ♀) are fertile

Hybridization is widespread in fish and constitutes an important mechanism in fish speciation. There is, however, little knowledge about hybridization in mormyrids. F1-interspecies hybrids between Campylomormyrus tamandua ♂ × C. compressirostris ♀ were investigated concerning: (1) fertility; (2) survival of F2-fish and (3) new gene combinations in the F2-generation concerning the structure of the electric organ and features of the electric organ discharge. These F1-hybrids achieved sexual maturity at about 12–13.5 cm total length. A breeding group comprising six males and 13 females spawned 28 times naturally proving these F1-fish to be fertile. On average 228 eggs were spawned, the average fertilization rate was 47.8%. Eggs started to hatch 70–72 h after fertilization, average hatching rate was 95.6%. Average mortality rate during embryonic development amounted to 2.3%. Average malformation rate during the free embryonic stage was 27.7%. Exogenous feeding started on day 11. In total, we raised 353 normally developed larvae all of which died consecutively, the oldest specimen reaching an age of 5 months. During survival, the activities of the larval and adult electric organs were recorded and the structure of the adult electric organ was investigated histologically.

Movement-mediated community assembly and coexistence

Organismal movement is ubiquitous and facilitates important ecological mechanisms that drive community and metacommunity composition and hence biodiversity. In most existing ecological theories and models in biodiversity research, movement is represented simplistically, ignoring the behavioural basis of movement and consequently the variation in behaviour at species and individual levels. However, as human endeavours modify climate and land use, the behavioural processes of organisms in response to these changes, including movement, become critical to understanding the resulting biodiversity loss. Here, we draw together research from different subdisciplines in ecology to understand the impact of individual-level movement processes on community-level patterns in species composition and coexistence. We join the movement ecology framework with the key concepts from metacommunity theory, community assembly and modern coexistence theory using the idea of micro-macro links, where various aspects of emergent movement behaviour scale up to local and regional patterns in species mobility and mobile-link-generated patterns in abiotic and biotic environmental conditions. These in turn influence both individual movement and, at ecological timescales, mechanisms such as dispersal limitation, environmental filtering, and niche partitioning. We conclude by highlighting challenges to and promising future avenues for data generation, data analysis and complementary modelling approaches and provide a brief outlook on how a new behaviour-based view on movement becomes important in understanding the responses of communities under ongoing environmental change.

Morphological differentiation in African weakly electric fish (genus Campylomormyrus) relates to substrate preferences

Under an ecological speciation scenario, the radiation of African weakly electric fish (genus Campylomormyrus) is caused by an adaptation to different food sources, associated with diversification of the electric organ discharge (EOD). This study experimentally investigates a phenotype–environment correlation to further support this scenario. Our behavioural experiments showed that three sympatric Campylomormyrus species with significantly divergent snout morphology differentially react to variation in substrate structure. While the short snout species (C. tamandua) exhibits preference to sandy substrate, the long snout species (C. rhynchophorus) significantly prefers a stone substrate for feeding. A third species with intermediate snout size (C. compressirostris) does not exhibit any substrate preference. This preference is matched with the observation that long-snouted specimens probe deeper into the stone substrate, presumably enabling them to reach prey more distant to the substrate surface. These findings suggest that the diverse feeding apparatus in the genus Campylomormyrus may have evolved in adaptation to specific microhabitats, i.e., substrate structures where these fish forage. Whether the parallel divergence in EOD is functionally related to this adaptation or solely serves as a prezygotic isolation mechanism remains to be elucidated.

Phylogeography of the European brook lamprey (Lampetra planeri) and the European river lamprey (Lampetra fluviatilis) species pair based on mitochondrial data

The European river lamprey Lampetra fluviatilis and the European brook lamprey Lampetra planeri (Block 1784) are classified as a paired species, characterized by notably different life histories but morphological similarities. Previous work has further shown limited genetic differentiation between these two species at the mitochondrial DNA level. Here, we expand on this previous work, which focused on lamprey species from the Iberian Peninsula in the south and mainland Europe in the north, by sequencing three mitochondrial marker regions of Lampetra individuals from five river systems in Ireland and five in southern Italy. Our results corroborate the previously identified pattern of genetic diversity for the species pair. We also show significant genetic differentiation between Irish and mainland European lamprey populations, suggesting another ichthyogeographic district distinct from those previously defined. Finally, our results stress the importance of southern Italian L. planeri populations, which maintain several private alleles and notable genetic diversity.

A simple dynamic model explains the diversity of island birds worldwide

Colonization, speciation and extinction are dynamic processes that influence global patterns of species richness^1-6. Island biogeography theory predicts that the contribution of these processes to the accumulation of species diversity depends on the area and isolation of the island^7,8. Notably, there has been no robust global test of this prediction for islands where speciation cannot be ignored⁹, because neither the appropriate data nor the analytical tools have been available. Here we address both deficiencies to reveal, for island birds, the empirical shape of the general relationships that determine how colonization, extinction and speciation rates co-vary with the area and isolation of islands. We compiled a global molecular phylogenetic dataset of birds on islands, based on the terrestrial avifaunas of 41 oceanic archipelagos worldwide (including 596 avian taxa), and applied a new analysis method to estimate the sensitivity of island-specific rates of colonization, speciation and extinction to island features (area and isolation). Our model predicts-with high explanatory power-several global relationships. We found a decline in colonization with isolation, a decline in extinction with area and an increase in speciation with area and isolation. Combining the theoretical foundations of island biogeography^7,8 with the temporal information contained in molecular phylogenies¹⁰ proves a powerful approach to reveal the fundamental relationships that govern variation in biodiversity across the planet.

Within species expressed genetic variability and gene expression response to different temperatures in the rotifer Brachionus calyciflorus sensu stricto

Genetic divergence is impacted by many factors, including phylogenetic history, gene flow, genetic drift, and divergent selection. Rotifers are an important component of aquatic ecosystems, and genetic variation is essential to their ongoing adaptive diversification and local adaptation. In addition to coding sequence divergence, variation in gene expression may relate to variable heat tolerance, and can impose ecological barriers within species. Temperature plays a significant role in aquatic ecosystems by affecting species abundance, spatio-temporal distribution, and habitat colonization. Recently described (formerly cryptic) species of the Brachionus calyciflorus complex exhibit different temperature tolerance both in natural and in laboratory studies, and show that B. calyciflorus sensu stricto (s.s.) is a thermotolerant species. Even within B. calyciflorus s.s., there is a tendency for further temperature specializations. Comparison of expressed genes allows us to assess the impact of stressors on both expression and sequence divergence among disparate populations within a single species. Here, we have used RNA-seq to explore expressed genetic diversity in B. calyciflorus s.s. in two mitochondrial DNA lineages with different phylogenetic histories and differences in thermotolerance. We identify a suite of candidate genes that may underlie local adaptation, with a particular focus on the response to sustained high or low temperatures. We do not find adaptive divergence in established candidate genes for thermal adaptation. Rather, we detect divergent selection among our two lineages in genes related to metabolism (lipid metabolism, metabolism of xenobiotics).

Elevated mutation rates are unlikely to evolve in sexual species, not even under rapid environmental change

BACKGROUND

Organisms are expected to respond to changing environmental conditions through local adaptation, range shift or local extinction. The process of local adaptation can occur by genetic changes or phenotypic plasticity, and becomes especially relevant when dispersal abilities or possibilities are somehow constrained. For genetic changes to occur, mutations are the ultimate source of variation and the mutation rate in terms of a mutator locus can be subject to evolutionary change. Recent findings suggest that the evolution of the mutation rate in a sexual species can advance invasion speed and promote adaptation to novel environmental conditions. Following this idea, this work uses an individual-based model approach to investigate if the mutation rate can also evolve in a sexual species experiencing different conditions of directional climate change, under different scenarios of colored stochastic environmental noise, probability of recombination and of beneficial mutations. The color of the noise mimicked investigating the evolutionary dynamics of the mutation rate in different habitats.

RESULTS

The results suggest that the mutation rate in a sexual species experiencing directional climate change scenarios can evolve and reach relatively high values mainly under conditions of complete linkage of the mutator locus and the adaptation locus. In contrast, when they are unlinked, the mutation rate can slightly increase only under scenarios where at least 50% of arising mutations are beneficial and the rate of environmental change is relatively fast. This result is robust under different scenarios of stochastic environmental noise, which supports the observation of no systematic variation in the mutation rate among organisms experiencing different habitats.

CONCLUSIONS

Given that 50% beneficial mutations may be an unrealistic assumption, and that recombination is ubiquitous in sexual species, the evolution of an elevated mutation rate in a sexual species experiencing directional climate change might be rather unlikely. Furthermore, when the percentage of beneficial mutations and the population size are small, sexual species (especially multicellular ones) producing few offspring may be expected to react to changing environments not by adaptive genetic change, but mainly through plasticity. Without the ability for a plastic response, such species may become - at least locally - extinct.

Enhanced North Pacific deep-ocean stratification by stronger intermediate water formation during Heinrich Stadial 1

The deglacial history of CO₂ release from the deep North Pacific remains unresolved. This is due to conflicting indications about subarctic Pacific ventilation changes based on various marine proxies, especially for Heinrich Stadial 1 (HS-1) when a rapid atmospheric CO₂ rise occurs. Here, we use a complex Earth System Model to investigate the deglacial North Pacific overturning and its control on ocean stratification. Our results show an enhanced intermediate-to-deep ocean stratification coeval with intensified North Pacific Intermediate Water (NPIW) formation during HS-1, compared to the Last Glacial Maximum. The stronger NPIW formation causes lower salinities and higher temperatures at intermediate depths. By lowering NPIW densities, this enlarges vertical density gradient and thus enhances intermediate-to-deep ocean stratification during HS-1. Physically, this process prevents the North Pacific deep waters from a better communication with the upper oceans, thus prolongs the existing isolation of glacial Pacific abyssal carbons during HS-1.

Supporting evidence for PCB pollution threatening global killer whale population

A recent Science report predicted the global killer whale population to collapse due to PCB pollution. Here we present empirical evidence, which supports and extends the reports' statement. In 2016, a neonate male killer whale stranded on the German island of Sylt. Neonatal attributes indicated an age of at least 3 days. The stomach contained ∼20 mL milk residue and no pathologies explaining the cause of death could be detected. Blubber samples presenting low lipid concentrations were analysed for persistent organic pollutants. Skin samples were collected for genotyping of the mitochondrial control region. The blubber PCB concentrations were very high [SPCBs, 225 mg/kg lipid weight (lw)], largely exceeding the PCB toxicity thresholds reported for the onset of immunosuppression [9 mg/kg lw ∑PCB] and for severe reproductive impairment [41 mg/kg lw ∑PCB] reported for marine mammals. Additionally, this individual showed equally high concentrations in p,p'-DDE [226 mg/kg lw], PBDEs [5 mg/kg lw] and liver mercury levels [1.1 μg/g dry weight dw]. These results suggest a high placental transfer of pollutants from mother to foetus. Consequently, blubber and plasma PCB concentrations and calf mortality rates are both high in primiparous females. With such high pollutant levels, this neonate had poor prerequisites for survival. The neonate belonged to Ecotype I (generalist feeder) and carried the mitochondrial haplotype 35 present in about 16% of the North Atlantic killer whale from or close to the North Sea. The relevance of this data becomes apparent in the UK West Coast Community, the UK's only residentorca population, which is currently composed of only eight individuals (each four males and females) and no calves have been reported over the last 19 years.Despite worldwide regulations, PCBs persist in the environment and remain a severe concern for killer whale populations, placing calves at high risk due to the mother-offspring PCB-transfer resulting in a high toxicological burden of the neonates.

Temporal and spatial patterns of mitochondrial haplotype and species distributions in Siberian larches inferred from ancient environmental DNA and modeling

Changes in species' distributions are classically projected based on their climate envelopes. For Siberian forests, which have a tremendous significance for vegetation-climate feedbacks, this implies future shifts of each of the forest-forming larch (Larix) species to the north-east. However, in addition to abiotic factors, reliable projections must assess the role of historical biogeography and biotic interactions. Here, we use sedimentary ancient DNA and individual-based modelling to investigate the distribution of larch species and mitochondrial haplotypes through space and time across the treeline ecotone on the southern Taymyr peninsula, which at the same time presents a boundary area of two larch species. We find spatial and temporal patterns, which suggest that forest density is the most influential driver determining the precise distribution of species and mitochondrial haplotypes. This suggests a strong influence of competition on the species' range shifts. These findings imply possible climate change outcomes that are directly opposed to projections based purely on climate envelopes. Investigations of such fine-scale processes of biodiversity change through time are possible using paleoenvironmental DNA, which is available much more readily than visible fossils and can provide information at a level of resolution that is not reached in classical palaeoecology.

Differential response to heat stress among evolutionary lineages of an aquatic invertebrate species complex

Under global warming scenarios, rising temperatures can constitute heat stress to which species may respond differentially. Within a described species, knowledge on cryptic diversity is of further relevance, as different lineages/cryptic species may respond differentially to environmental change. The Brachionus calyciflorus species complex (Rotifera), which was recently described using integrative taxonomy, is an essential component of aquatic ecosystems. Here, we tested the hypothesis that these (formerly cryptic) species differ in their heat tolerance. We assigned 47 clones with nuclear ITS1 (nuITS1) and mitochondrial COI (mtCOI) markers to evolutionary lineages, now named B. calyciflorus sensu stricto (s.s.) and B. fernandoi We selected 15 representative clones and assessed their heat tolerance as a bi-dimensional phenotypic trait affected by both the intensity and duration of heat stress. We found two distinct groups, with B. calyciflorus s.s. clones having higher heat tolerance than the novel species B. fernandoi This apparent temperature specialization among former cryptic species underscores the necessity of a sound species delimitation and assignment, when organismal responses to environmental changes are investigated.

High-quality whole-genome sequence of an abundant Holarctic odontocete, the harbour porpoise (Phocoena phocoena)

The harbour porpoise (Phocoena phocoena) is a highly mobile cetacean found across the Northern hemisphere. It occurs in coastal waters and inhabits basins that vary broadly in salinity, temperature and food availability. These diverse habitats could drive subtle differentiation among populations, but examination of this would be best conducted with a robust reference genome. Here, we report the first harbour porpoise genome, assembled de novo from an individual originating in the Kattegat Sea (Sweden). The genome is one of the most complete cetacean genomes currently available, with a total size of 2.39 Gb and 50% of the total length found in just 34 scaffolds. Using 122 of the longest scaffolds, we were able to show high levels of synteny with the genome of the domestic cattle (Bos taurus). Our draft annotation comprises 22,154 predicted genes, which we further annotated through matches to the NCBI nucleotide database, GO categorization and motif prediction. Within the predicted genes, we have confirmed the presence of >20 genes or gene families that have been associated with adaptive evolution in other cetaceans. Overall, this genome assembly and draft annotation represent a crucial addition to the genomic resources currently available for the study of porpoises and Phocoenidae evolution, phylogeny and conservation.

Inter-individual differences in contamination profiles as tracer of social group association in stranded sperm whales

Ecological and physiological factors lead to different contamination patterns in individual marine mammals. The objective of the present study was to assess whether variations in contamination profiles are indicative of social structures of young male sperm whales as they might reflect a variation in feeding preferences and/or in utilized feeding grounds. We used a total of 61 variables associated with organic compounds and trace element concentrations measured in muscle, liver, kidney and blubber gained from 24 sperm whales that stranded in the North Sea in January and February 2016. Combining contaminant and genetic data, there is evidence for at least two cohorts with different origin among these stranded sperm whales; one from the Canary Island region and one from the northern part of the Atlantic. While genetic data unravel relatedness and kinship, contamination data integrate over areas, where animals occured during their lifetime. Especially in long-lived animals with a large migratory potential, as sperm whales, contamination data may carry highly relevant information about aggregation through time and space.

Electric pulse characteristics can enable species recognition in African weakly electric fish species

Communication is key to a wide variety of animal behaviours and multiple modalities are often involved in this exchange of information from sender to receiver. The communication of African weakly electric fish, however, is thought to be predominantly unimodal and is mediated by their electric sense, in which species-specific electric organ discharges (EODs) are generated in a context-dependent and thus variable sequence of pulse intervals (SPI). While the primary function of the electric sense is considered to be electrolocation, both of its components likely carry information regarding identity of the sender. However, a clear understanding of their contribution to species recognition is incomplete. We therefore analysed these two electrocommunication components (EOD waveform and SPI statistics) in two sympatric mormyrid Campylomormyrus species. In a set of five playback conditions, we further investigated which components may drive interspecific recognition and discrimination. While we found that both electrocommunication components are species-specific, the cues necessary for species recognition differ between the two species studied. While the EOD waveform and SPI were both necessary and sufficient for species recognition in C. compressirostris males, C. tamandua males apparently utilize other, non-electric modalities. Mapped onto a recent phylogeny, our results suggest that discrimination by electric cues alone may be an apomorphic trait evolved during a recent radiation in this taxon.

North Pacific freshwater events linked to changes in glacial ocean circulation

There is compelling evidence that episodic deposition of large volumes of freshwater into the oceans strongly influenced global ocean circulation and climate variability during glacial periods^1,2. In the North Atlantic region, episodes of massive freshwater discharge to the North Atlantic Ocean were related to distinct cold periods known as Heinrich Stadials^1-3. By contrast, the freshwater history of the North Pacific region remains unclear, giving rise to persistent debates about the existence and possible magnitude of climate links between the North Pacific and North Atlantic oceans during Heinrich Stadials^4,5. Here we find that there was a strong connection between changes in North Atlantic circulation during Heinrich Stadials and injections of freshwater from the North American Cordilleran Ice Sheet to the northeastern North Pacific. Our record of diatom δ¹⁸O (a measure of the ratio of the stable oxygen isotopes ¹⁸O and ¹⁶O) over the past 50,000 years shows a decrease in surface seawater δ¹⁸O of two to three per thousand, corresponding to a decline in salinity of roughly two to four practical salinity units. This coincided with enhanced deposition of ice-rafted debris and a slight cooling of the sea surface in the northeastern North Pacific during Heinrich Stadials 1 and 4, but not during Heinrich Stadial 3. Furthermore, results from our isotope-enabled model⁶ suggest that warming of the eastern Equatorial Pacific during Heinrich Stadials was crucial for transmitting the North Atlantic signal to the northeastern North Pacific, where the associated subsurface warming resulted in a discernible freshwater discharge from the Cordilleran Ice Sheet during Heinrich Stadials 1 and 4. However, enhanced background cooling across the northern high latitudes during Heinrich Stadial 3-the coldest period in the past 50,000 years⁷-prevented subsurface warming of the northeastern North Pacific and thus increased freshwater discharge from the Cordilleran Ice Sheet. In combination, our results show that nonlinear ocean-atmosphere background interactions played a complex role in the dynamics linking the freshwater discharge responses of the North Atlantic and North Pacific during glacial periods.

Male size, not female preferences influence female reproductive success in a poeciliid fish (Poecilia latipinna): a combined behavioural/genetic approach

Objective

We investigated the potential role of indirect benefits for female mate preferences in a highly promiscuous species of live-bearing fishes, the sailfin molly Poecilia latipinna using an integrative approach that combines methods from animal behavior, life-history evolution, and genetics. Males of this species solely contribute sperm for reproduction, and consequently females do not receive any direct benefits. Despite this, females typically show clear mate preferences. It has been suggested that females can increase their reproductive success through indirect benefits from choosing males of higher quality.

Results

Although preferences for large body size have been recorded as an honest signal for genetic quality, this particular study resulted in female preference being unaffected by male body size. Nonetheless, larger males did sire more offspring, but with no effect on offspring quality. This study presents a methodical innovation by combining preference testing with life history measurements—such as the determination of the dry weight of fish embryos—and paternity analyses on single fish embryos.

Electronic supplementary material

The online version of this article (10.1186/s13104-018-3487-2) contains supplementary material, which is available to authorized users.

Male-mediated species recognition among African weakly electric fishes

Effective communication among sympatric species is often instrumental for behavioural isolation, where the failure to successfully discriminate between potential mates could lead to less fit hybrid offspring. Discrimination between con- and heterospecifics tends to occur more often in the sex that invests more in offspring production, i.e. females, but males may also mediate reproductive isolation. In this study, we show that among two Campylomormyrus African weakly electric fish species, males preferentially associate with conspecific females during choice tests using live fish as stimuli, i.e. when all sensory modalities potentially used for communication were present. We then conducted playback experiments to determine whether the species-specific electric organ discharge (EOD) used for electrocommunication serves as the cue for this conspecific association preference. Interestingly, only C. compressirostris males associated significantly more with the conspecific EOD waveform when playback stimuli were provided, while no such association preference was observed in C. tamandua males. Given our results, the EOD appears to serve, in part, as a male-mediated pre-zygotic isolation mechanism among sympatric species. However, the failure of C. tamandua males to discriminate between con- and heterospecific playback discharges suggests that multiple modalities may be necessary for species recognition in some African weakly electric fish species.

The gonadal transcriptome of the unisexual Amazon molly Poecilia formosa in comparison to its sexual ancestors, Poecilia mexicana and Poecilia latipinna

Background

The unisexual Amazon molly (Poecilia formosa) originated from a hybridization between two sexual species, the sailfin molly (Poecilia latipinna) and the Atlantic molly (Poecilia mexicana). The Amazon molly reproduces clonally via sperm-dependent parthenogenesis (gynogenesis), in which the sperm of closely related species triggers embryogenesis of the apomictic oocytes, but typically does not contribute genetic material to the next generation. We compare for the first time the gonadal transcriptome of the Amazon molly to those of both ancestral species, P. mexicana and P. latipinna.

Results

We sequenced the gonadal transcriptomes of the P. formosa and its parental species P. mexicana and P. latipinna using Illumina RNA-sequencing techniques (paired-end, 100 bp). De novo assembly of about 50 million raw read pairs for each species was performed using Trinity, yielding 106,922 transcripts for P. formosa, 115,175 for P. latipinna, and 133,025 for P. mexicana after eliminating contaminations. On the basis of sequence similarity comparisons to other teleost species and the UniProt databases, functional annotation, and differential expression analysis, we demonstrate the similarity of the transcriptomes among the three species. More than 40% of the transcripts for each species were functionally annotated and about 70% were assigned to orthologous genes of a closely related species. Differential expression analysis between the sexual and unisexual species uncovered 2035 up-regulated and 564 down-regulated genes in P. formosa. This was exemplary validated for six genes by qRT-PCR.

Conclusions

We identified more than 130 genes related to meiosis and reproduction within the apomictically reproducing P. formosa. Overall expression of these genes seems to be down-regulated in the P. formosa transcriptome compared to both ancestral species (i.e., 106 genes down-regulated, 29 up-regulated). A further 35 meiosis and reproduction related genes were not found in the P. formosa transcriptome, but were only expressed in the sexual species. Our data support the hypothesis of general down-regulation of meiosis-related genes in the apomictic Amazon molly. Furthermore, the obtained dataset and identified gene catalog will serve as a resource for future research on the molecular mechanisms behind the reproductive mode of this unisexual species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4382-2) contains supplementary material, which is available to authorized users.

Equilibrium Bird Species Diversity in Atlantic Islands

Half a century ago, MacArthur and Wilson proposed that the number of species on islands tends toward a dynamic equilibrium diversity around which species richness fluctuates [1]. The current prevailing view in island biogeography accepts the fundamentals of MacArthur and Wilson's theory [2] but questions whether their prediction of equilibrium can be fulfilled over evolutionary timescales, given the unpredictable and ever-changing nature of island geological and biotic features [3-7]. Here we conduct a complete molecular phylogenetic survey of the terrestrial bird species from four oceanic archipelagos that make up the diverse Macaronesian bioregion-the Azores, the Canary Islands, Cape Verde, and Madeira [8, 9]. We estimate the times at which birds colonized and speciated in the four archipelagos, including many previously unsampled endemic and non-endemic taxa and their closest continental relatives. We develop and fit a new multi-archipelago dynamic stochastic model to these data, explicitly incorporating information from 91 taxa, both extant and extinct. Remarkably, we find that all four archipelagos have independently achieved and maintained a dynamic equilibrium over millions of years. Biogeographical rates are homogeneous across archipelagos, except for the Canary Islands, which exhibit higher speciation and colonization. Our finding that the avian communities of the four Macaronesian archipelagos display an equilibrium diversity pattern indicates that a diversity plateau may be rapidly achieved on islands where rates of in situ radiation are low and extinction is high. This study reveals that equilibrium processes may be more prevalent than recently proposed, supporting MacArthur and Wilson's 50-year-old theory.

Intragenus (Campylomormyrus) and intergenus hybrids in mormyrid fish: Physiological and histological investigations of the electric organ ontogeny

African weakly electric mormyrid fish show a high diversity of their electric organ discharge (EOD) both across and within genera. Thanks to a recently developed technique of artificial reproduction in mormyrid fish, we were able to perform hybridizations between different genera and within one genus (Campylomormyrus). The hybrids of intergenus hybridizations exhibited different degrees of reduced survival related to the phylogenetic distance of the parent species: hybrids of the crosses between C. rhynchophorus and its sister genus Gnathonemus survived and developed normally. Hybrids between C. rhynchophorus and a Mormyrus species (a more basal clade compared to Campylomormyrus s) survived up to 42days and developed many malformations, e.g., at the level of the unpaired fins. Hybrids between C. numenius and Hippopotamyrus pictus (a derived clade, only distantly related to Campylomormyrus) only survived for two days during embryological development. Eight different hybrid combinations among five Campylomormyrus species (C. tamandua, C. compressirostris, C. tshokwe, C. rhynchophorus, C. numenius) were performed. The aim of the hybridizations was to combine species with (1) either caudal or rostral position of the main stalk innervating the electrocytes in the electric organ and (2) short, median or long duration of their EOD. The hybrids, though they are still juveniles, show very interesting features concerning electrocyte geometry as well as EOD form and duration: the caudal position of the stalk is prevailing over the rostral position, and the penetration of the stalk is dominant over the non-penetrating feature (in the Campylomormyrus hybrids); in the hybrid between C. rhynchophorus and Gnathonemus petersii it is the opposite. When crossing species with long and short EODs, it is always the long duration EOD that is expressed in the hybrids. The F1-Hybrids of the cross C. tamandua×C. compressirostris are fertile: viable F2-fish could be obtained with artificial reproduction.

Spatially Explicit Analysis of Genome-Wide SNPs Detects Subtle Population Structure in a Mobile Marine Mammal, the Harbor Porpoise

The population structure of the highly mobile marine mammal, the harbor porpoise (Phocoena phocoena), in the Atlantic shelf waters follows a pattern of significant isolation-by-distance. The population structure of harbor porpoises from the Baltic Sea, which is connected with the North Sea through a series of basins separated by shallow underwater ridges, however, is more complex. Here, we investigated the population differentiation of harbor porpoises in European Seas with a special focus on the Baltic Sea and adjacent waters, using a population genomics approach. We used 2872 single nucleotide polymorphisms (SNPs), derived from double digest restriction-site associated DNA sequencing (ddRAD-seq), as well as 13 microsatellite loci and mitochondrial haplotypes for the same set of individuals. Spatial principal components analysis (sPCA), and Bayesian clustering on a subset of SNPs suggest three main groupings at the level of all studied regions: the Black Sea, the North Atlantic, and the Baltic Sea. Furthermore, we observed a distinct separation of the North Sea harbor porpoises from the Baltic Sea populations, and identified splits between porpoise populations within the Baltic Sea. We observed a notable distinction between the Belt Sea and the Inner Baltic Sea sub-regions. Improved delineation of harbor porpoise population assignments for the Baltic based on genomic evidence is important for conservation management of this endangered cetacean in threatened habitats, particularly in the Baltic Sea proper. In addition, we show that SNPs outperform microsatellite markers and demonstrate the utility of RAD-tags from a relatively small, opportunistically sampled cetacean sample set for population diversity and divergence analysis.

Kin Recognition in a Clonal Fish, Poecilia formosa

Relatedness strongly influences social behaviors in a wide variety of species. For most species, the highest typical degree of relatedness is between full siblings with 50% shared genes. However, this is poorly understood in species with unusually high relatedness between individuals: clonal organisms. Although there has been some investigation into clonal invertebrates and yeast, nothing is known about kin selection in clonal vertebrates. We show that a clonal fish, the Amazon molly (Poecilia formosa), can distinguish between different clonal lineages, associating with genetically identical, sister clones, and use multiple sensory modalities. Also, they scale their aggressive behaviors according to the relatedness to other females: they are more aggressive to non-related clones. Our results demonstrate that even in species with very small genetic differences between individuals, kin recognition can be adaptive. Their discriminatory abilities and regulation of costly behaviors provides a powerful example of natural selection in species with limited genetic diversity.

Sequence Evolution and Expression of the Androgen Receptor and Other Pathway-Related Genes in a Unisexual Fish, the Amazon Molly, Poecilia formosa, and Its Bisexual Ancestors

The all-female Amazon molly (Poecilia formosa) originated from a single hybridization of two bisexual ancestors, Atlantic molly (Poecilia mexicana) and sailfin molly (Poecilia latipinna). As a gynogenetic species, the Amazon molly needs to copulate with a heterospecific male, but the genetic information of the sperm-donor does not contribute to the next generation, as the sperm only acts as the trigger for the diploid eggs' embryogenesis. Here, we study the sequence evolution and gene expression of the duplicated genes coding for androgen receptors (ars) and other pathway-related genes, i.e., the estrogen receptors (ers) and cytochrome P450, family19, subfamily A, aromatase genes (cyp19as), in the Amazon molly, in comparison to its bisexual ancestors. Mollies possess-as most other teleost fish-two copies of the ar, er, and cyp19a genes, i.e., arα/arβ, erα/erβ1, and cyp19a1 (also referred as cyp19a1a)/cyp19a2 (also referred to as cyp19a1b), respectively. Non-synonymous single nucleotide polymorphisms (SNPs) among the ancestral bisexual species were generally predicted not to alter protein function. Some derived substitutions in the P. mexicana and one in P. formosa are predicted to impact protein function. We also describe the gene expression pattern of the ars and pathway-related genes in various tissues (i.e., brain, gill, and ovary) and provide SNP markers for allele specific expression research. As a general tendency, the levels of gene expression were lowest in gill and highest in ovarian tissues, while expression levels in the brain were intermediate in most cases. Expression levels in P. formosa were conserved where expression did not differ between the two bisexual ancestors. In those cases where gene expression levels significantly differed between the bisexual species, P. formosa expression was always comparable to the higher expression level among the two ancestors. Interestingly, erβ1 was expressed neither in brain nor in gill in the analyzed three molly species, which implies a more important role of erα in the estradiol synthesis pathway in these tissues. Furthermore, our data suggest that interactions of steroid-signaling pathway genes differ across tissues, in particular the interactions of ars and cyp19as.

Radiocarbon constraints on the extent and evolution of the South Pacific glacial carbon pool

During the last deglaciation, the opposing patterns of atmospheric CO2 and radiocarbon activities (Δ(14)C) suggest the release of (14)C-depleted CO2 from old carbon reservoirs. Although evidences point to the deep Pacific as a major reservoir of this (14)C-depleted carbon, its extent and evolution still need to be constrained. Here we use sediment cores retrieved along a South Pacific transect to reconstruct the spatio-temporal evolution of Δ(14)C over the last 30,000 years. In ∼2,500-3,600 m water depth, we find (14)C-depleted deep waters with a maximum glacial offset to atmospheric (14)C (ΔΔ(14)C=-1,000‰). Using a box model, we test the hypothesis that these low values might have been caused by an interaction of aging and hydrothermal CO2 influx. We observe a rejuvenation of circumpolar deep waters synchronous and potentially contributing to the initial deglacial rise in atmospheric CO2. These findings constrain parts of the glacial carbon pool to the deep South Pacific.

Genetic data from algae sedimentary DNA reflect the influence of environment over geography

Genetic investigations on eukaryotic plankton confirmed the existence of modern biogeographic patterns, but analyses of palaeoecological data exploring the temporal variability of these patterns have rarely been presented. Ancient sedimentary DNA proved suitable for investigations of past assemblage turnover in the course of environmental change, but genetic relatedness of the identified lineages has not yet been undertaken. Here, we investigate the relatedness of diatom lineages in Siberian lakes along environmental gradients (i.e. across treeline transects), over geographic distance and through time (i.e. the last 7000 years) using modern and ancient sedimentary DNA. Our results indicate that closely-related Staurosira lineages occur in similar environments and less-related lineages in dissimilar environments, in our case different vegetation and co-varying climatic and limnic variables across treeline transects. Thus our study reveals that environmental conditions rather than geographic distance is reflected by diatom-relatedness patterns in space and time. We tentatively speculate that the detected relatedness pattern in Staurosira across the treeline could be a result of adaptation to diverse environmental conditions across the arctic boreal treeline, however, a geographically-driven divergence and subsequent repopulation of ecologically different habitats might also be a potential explanation for the observed pattern.