Global mesozooplankton communities show lower connectivity in deep oceanic layers

Mesozooplankton is a key component of the ocean, regulating global processes such as the carbon pump, and ensuring energy transfer from lower to higher trophic levels. Yet, knowledge on mesozooplankton diversity, distribution and connectivity at global scale is still fragmented. To fill this gap, we applied DNA metabarcoding to mesozooplankton samples collected during the Malaspina-2010 circumnavigation expedition across the Atlantic, Indian and Pacific oceans from the surface to bathypelagic depths. We highlight the still scarce knowledge on global mesozooplankton diversity and identify the Indian Ocean and the deep sea as the oceanic regions with the highest proportion of hidden diversity. We report no consistent alpha-diversity patterns for mesozooplankton at a global scale, neither across vertical nor horizontal gradients. However, beta-diversity analysis suggests horizontal and vertical structuring of mesozooplankton communities mostly attributed to turnover and reveals an increase in mesozooplankton beta-diversity with depth, indicating reduced connectivity at deeper layers. Additionally, we identify a water mass type-mediated structuring of mesozooplankton bathypelagic communities instead of an oceanic basin-mediated as observed at upper layers. This suggests limited dispersal at deep ocean layers, most likely due to weaker currents and lower mixing of water mass types, thus reinforcing the importance of oceanic currents and barriers to dispersal in shaping global plankton communities.

Loss of species and genetic diversity during colonization: Insights from acanthocephalan parasites in northern European seals

Studies on host-parasite systems that have experienced distributional shifts, range fragmentation, and population declines in the past can provide information regarding how parasite community richness and genetic diversity will change as a result of anthropogenic environmental changes in the future. Here, we studied how sequential postglacial colonization, shifts in habitat, and reduced host population sizes have influenced species richness and genetic diversity of Corynosoma (Acanthocephala: Polymorphidae) parasites in northern European marine, brackish, and freshwater seal populations. We collected Corynosoma population samples from Arctic, Baltic, Ladoga, and Saimaa ringed seal subspecies and Baltic gray seals, and then applied COI barcoding and triple-enzyme restriction-site associated DNA (3RAD) sequencing to delimit species, clarify their distributions and community structures, and elucidate patterns of intraspecific gene flow and genetic diversity. Our results showed that Corynosoma species diversity reflected host colonization histories and population sizes, with four species being present in the Arctic, three in the Baltic Sea, two in Lake Ladoga, and only one in Lake Saimaa. We found statistically significant population-genetic differentiation within all three Corynosoma species that occur in more than one seal (sub)species. Genetic diversity tended to be high in Corynosoma populations originating from Arctic ringed seals and low in the landlocked populations. Our results indicate that acanthocephalan communities in landlocked seal populations are impoverished with respect to both species and intraspecific genetic diversity. Interestingly, the loss of genetic diversity within Corynosoma species seems to have been less drastic than in their seal hosts, possibly due to their large local effective population sizes resulting from high infection intensities and effective intra-host population mixing. Our study highlights the utility of genomic methods in investigations of community composition and genetic diversity of understudied parasites.

Genomic evidence for contrasting patterns of host-associated genetic differentiation across shared host-plant species in leaf- and bud-galling sawflies

Resource specialization and ecological speciation arising through host-associated genetic differentiation (HAD) are frequently invoked as an explanation for the high diversity of plant-feeding insects and other organisms with a parasitic lifestyle. While genetic studies have demonstrated numerous examples of HAD in insect herbivores, the rarity of comparative studies means that we still lack an understanding of how deterministic HAD is, and whether patterns of host shifts can be predicted over evolutionary timescales. We applied genome-wide single nucleotide polymorphism and mitochondrial DNA sequence data obtained through genome resequencing to define species limits and to compare host-plant use in population samples of leaf- and bud-galling sawflies (Hymenoptera: Tenthredinidae: Nematinae) collected from seven shared willow (Salicaceae: Salix) host species. To infer the repeatability of long-term cophylogenetic patterns, we also contrasted the phylogenies of the two galler groups with each other as well as with the phylogeny of their Salix hosts estimated based on RADseq data. We found clear evidence for host specialization and HAD in both of the focal galler groups, but also that leaf gallers are more specialized to single host species compared with most bud gallers. In contrast to bud gallers, leaf gallers also exhibited statistically significant cophylogenetic signal with their Salix hosts. The observed discordant patterns of resource specialization and host shifts in two related galler groups that have radiated in parallel across a shared resource base indicate a lack of evolutionary repeatability in the focal system, and suggest that short- and long-term host use and ecological diversification in plant-feeding insects are dominated by stochasticity and/or lineage-specific effects.

Mountain- and brown hare genetic polymorphisms to survey local adaptations and conservation status of the heath hare (Lepus timidus sylvaticus, Nilsson 1831)

We provide the first whole genome sequences from three specimens of the mountain hare subspecies the heath hare (Lepus timidus sylvaticus), along with samples from two mountain hares (Lepus timidus timidus) and two brown hares (Lepus europaeus) from Sweden. The heath hare has a unique grey winter pelage as compared to other mountain hares (white) and brown hares (mostly brown), and face regional extinction, likely due to competitive exclusion from the non-native brown hare. Whole genome resequencing from the seven hare specimens were mapped to the Lepus timidus pseudoreference genome and used for detection of 11,363,883 polymorphic nucleotide positions. The data presented here could be useful for addressing local adaptations and conservation status of mountain hares and brown hares in Sweden, including unique subspecies.

Measurement(s)	whole genome sequencing
Technology Type(s)	Illumina HiSeq X
Factor Type(s)	Species
Sample Characteristic - Organism	Lepus timidus • Lepus europaeus • Lepus timidus sylvaticus
Sample Characteristic - Location	Sweden

DNA barcoding reveals different cestode helminth species in northern European marine and freshwater ringed seals

Three subspecies of the ringed seal (Pusa hispida) are found in northeastern Europe: P. h. botnica in the Baltic Sea, P. h saimensis in Lake Saimaa in Finland, and P. h. ladogensis in Lake Ladoga in Russia. We investigated the poorly-known cestode helminth communities of these closely related but ecologically divergent subspecies using COI barcode data. Our results show that, while cestodes from the Baltic Sea represent Schistocephalus solidus, all worms from the two lakes are identified as Ligula intestinalis, a species that has previously not been reported from seals. The observed shift in cestode communities appears to be driven by differential availability of intermediate fish host species in marine vs. freshwater environments. Both observed cestode species normally infect fish-eating birds, so further work is required to elucidate the health and conservation implications of cestode infections in European ringed seals, whether L. intestinalis occurs also in marine ringed seals, and whether the species is able to reproduce in seal hosts. In addition, a deep barcode divergence found within S. solidus suggests the presence of cryptic diversity under this species name.

•

COI barcoding reveals different cestodes in marine and freshwater ringed seals.

•

Ligula intestinalis is reported for the first time from seals.

•

A deep barcode divergence is found within Schistocephalus solidus in the Baltic Sea.

Microbiomes of willow-galling sawflies: effects of host plant, gall type, and phylogeny on community structure and function

While free-living herbivorous insects are thought to harbor microbial communities composed of transient bacteria derived from their diet, recent studies indicate that insects that induce galls on plants may be involved in more intimate host-microbe relationships. We used 16S rDNA metabarcoding to survey larval microbiomes of 20 nematine sawfly species that induce bud or leaf galls on 13 Salix species. The 391 amplicon sequence variants (ASVs) detected represented 69 bacterial genera in six phyla. Multi-variate statistical analyses showed that the structure of larval microbiomes is influenced by willow host species as well as by gall type. Nevertheless, a "core" microbiome composed of 58 ASVs is shared widely across the focal galler species. Within the core community, the presence of many abundant, related ASVs representing multiple distantly related bacterial taxa is reflected as a statistically significant effect of bacterial phylogeny on galler-microbe associations. Members of the core community have a variety of inferred functions, including degradation of phenolic compounds, nutrient supplementation, and production of plant hormones. Hence, our results support suggestions of intimate and diverse interactions between galling insects and microbes and add to a growing body of evidence that microbes may play a role in the induction of insect galls on plants.

Desert plant bacteria reveal host influence and beneficial plant growth properties

Deserts, such as those found in Saudi Arabia, are one of the most hostile places for plant growth. However, desert plants are able to impact their surrounding microbial community and select beneficial microbes that promote their growth under these extreme conditions. In this study, we examined the soil, rhizosphere and endosphere bacterial communities of four native desert plants Tribulus terrestris, Zygophyllum simplex, Panicum turgidum and Euphorbia granulata from the Southwest (Jizan region), two of which were also found in the Midwest (Al Wahbah area) of Saudi Arabia. While the rhizosphere bacterial community mostly resembled that of the highly different surrounding soils, the endosphere composition was strongly correlated with its host plant phylogeny. In order to assess whether any of the native bacterial endophytes might have a role in plant growth under extreme conditions, we analyzed the properties of 116 cultured bacterial isolates that represent members of the phyla Proteobacteria, Bacteroidetes, Actinobacteria and Firmicutes. Our analysis shows that different strains have highly different biochemical properties with respect to nutrient acquisition, hormone production and growth under stress conditions. More importantly, eleven of the isolated strains could confer salinity stress tolerance to the experimental model plant Arabidopsis thaliana suggesting some of these plant-associated bacteria might be useful for improving crop desert agriculture.

The Genome Sequence of the Wild Tomato Solanum pimpinellifolium Provides Insights Into Salinity Tolerance

Solanum pimpinellifolium, a wild relative of cultivated tomato, offers a wealth of breeding potential for desirable traits such as tolerance to abiotic and biotic stresses. Here, we report the genome assembly and annotation of S. pimpinellifolium 'LA0480.' Moreover, we present phenotypic data from one field experiment that demonstrate a greater salinity tolerance for fruit- and yield-related traits in S. pimpinellifolium compared with cultivated tomato. The 'LA0480' genome assembly size (811 Mb) and the number of annotated genes (25,970) are within the range observed for other sequenced tomato species. We developed and utilized the Dragon Eukaryotic Analyses Platform (DEAP) to functionally annotate the 'LA0480' protein-coding genes. Additionally, we used DEAP to compare protein function between S. pimpinellifolium and cultivated tomato. Our data suggest enrichment in genes involved in biotic and abiotic stress responses. To understand the genomic basis for these differences in S. pimpinellifolium and S. lycopersicum, we analyzed 15 genes that have previously been shown to mediate salinity tolerance in plants. We show that S. pimpinellifolium has a higher copy number of the inositol-3-phosphate synthase and phosphatase genes, which are both key enzymes in the production of inositol and its derivatives. Moreover, our analysis indicates that changes occurring in the inositol phosphate pathway may contribute to the observed higher salinity tolerance in 'LA0480.' Altogether, our work provides essential resources to understand and unlock the genetic and breeding potential of S. pimpinellifolium, and to discover the genomic basis underlying its environmental robustness.

Finding Nemo's Genes: A chromosome-scale reference assembly of the genome of the orange clownfish Amphiprion percula

The iconic orange clownfish, Amphiprion percula, is a model organism for studying the ecology and evolution of reef fishes, including patterns of population connectivity, sex change, social organization, habitat selection and adaptation to climate change. Notably, the orange clownfish is the only reef fish for which a complete larval dispersal kernel has been established and was the first fish species for which it was demonstrated that antipredator responses of reef fishes could be impaired by ocean acidification. Despite its importance, molecular resources for this species remain scarce and until now it lacked a reference genome assembly. Here, we present a de novo chromosome-scale assembly of the genome of the orange clownfish Amphiprion percula. We utilized single-molecule real-time sequencing technology from Pacific Biosciences to produce an initial polished assembly comprised of 1,414 contigs, with a contig N50 length of 1.86 Mb. Using Hi-C-based chromatin contact maps, 98% of the genome assembly were placed into 24 chromosomes, resulting in a final assembly of 908.8 Mb in length with contig and scaffold N50s of 3.12 and 38.4 Mb, respectively. This makes it one of the most contiguous and complete fish genome assemblies currently available. The genome was annotated with 26,597 protein-coding genes and contains 96% of the core set of conserved actinopterygian orthologs. The availability of this reference genome assembly as a community resource will further strengthen the role of the orange clownfish as a model species for research on the ecology and evolution of reef fishes.

DNA methylation regulates transcriptional homeostasis of algal endosymbiosis in the coral model Aiptasia

The symbiotic relationship between cnidarians and dinoflagellates is the cornerstone of coral reef ecosystems. Although research has focused on the molecular mechanisms underlying this symbiosis, the role of epigenetic mechanisms, that is, the study of heritable changes that do not involve changes in the DNA sequence, is unknown. To assess the role of DNA methylation in the cnidarian-dinoflagellate symbiosis, we analyzed genome-wide CpG methylation, histone associations, and transcriptomic states of symbiotic and aposymbiotic anemones in the model system Aiptasia. We found that methylated genes are marked by histone 3 lysine 36 trimethylation (H3K36me3) and show significant reduction of spurious transcription and transcriptional noise, revealing a role of DNA methylation in the maintenance of transcriptional homeostasis. Changes in DNA methylation and expression show enrichment for symbiosis-related processes, such as immunity, apoptosis, phagocytosis recognition, and phagosome formation, and reveal intricate interactions between the underlying pathways. Our results demonstrate that DNA methylation provides an epigenetic mechanism of transcriptional homeostasis that responds to symbiosis.

Epigenome-associated phenotypic acclimatization to ocean acidification in a reef-building coral

There are increasing concerns that the current rate of climate change might outpace the ability of reef-building corals to adapt to future conditions. Work on model systems has shown that environmentally induced alterations in DNA methylation can lead to phenotypic acclimatization. While DNA methylation has been reported in corals and is thought to associate with phenotypic plasticity, potential mechanisms linked to changes in whole-genome methylation have yet to be elucidated. We show that DNA methylation significantly reduces spurious transcription in the coral Stylophora pistillata. Furthermore, we find that DNA methylation also reduces transcriptional noise by fine-tuning the expression of highly expressed genes. Analysis of DNA methylation patterns of corals subjected to long-term pH stress showed widespread changes in pathways regulating cell cycle and body size. Correspondingly, we found significant increases in cell and polyp sizes that resulted in more porous skeletons, supporting the hypothesis that linear extension rates are maintained under conditions of reduced calcification. These findings suggest an epigenetic component in phenotypic acclimatization that provides corals with an additional mechanism to cope with environmental change.

Pushing the limits of photoreception in twilight conditions: The rod-like cone retina of the deep-sea pearlsides

Most vertebrates have a duplex retina comprising two photoreceptor types, rods for dim-light (scotopic) vision and cones for bright-light (photopic) and color vision. However, deep-sea fishes are only active in dim-light conditions; hence, most species have lost their cones in favor of a simplex retina composed exclusively of rods. Although the pearlsides, Maurolicus spp., have such a pure rod retina, their behavior is at odds with this simplex visual system. Contrary to other deep-sea fishes, pearlsides are mostly active during dusk and dawn close to the surface, where light levels are intermediate (twilight or mesopic) and require the use of both rod and cone photoreceptors. This study elucidates this paradox by demonstrating that the pearlside retina does not have rod photoreceptors only; instead, it is composed almost exclusively of transmuted cone photoreceptors. These transmuted cells combine the morphological characteristics of a rod photoreceptor with a cone opsin and a cone phototransduction cascade to form a unique photoreceptor type, a rod-like cone, specifically tuned to the light conditions of the pearlsides' habitat (blue-shifted light at mesopic intensities). Combining properties of both rods and cones into a single cell type, instead of using two photoreceptor types that do not function at their full potential under mesopic conditions, is likely to be the most efficient and economical solution to optimize visual performance. These results challenge the standing paradigm of the function and evolution of the vertebrate duplex retina and emphasize the need for a more comprehensive evaluation of visual systems in general.

The genome of Chenopodium quinoa

Chenopodium quinoa (quinoa) is a highly nutritious grain identified as an important crop to improve world food security. Unfortunately, few resources are available to facilitate its genetic improvement. Here we report the assembly of a high-quality, chromosome-scale reference genome sequence for quinoa, which was produced using single-molecule real-time sequencing in combination with optical, chromosome-contact and genetic maps. We also report the sequencing of two diploids from the ancestral gene pools of quinoa, which enables the identification of sub-genomes in quinoa, and reduced-coverage genome sequences for 22 other samples of the allotetraploid goosefoot complex. The genome sequence facilitated the identification of the transcription factor likely to control the production of anti-nutritional triterpenoid saponins found in quinoa seeds, including a mutation that appears to cause alternative splicing and a premature stop codon in sweet quinoa strains. These genomic resources are an important first step towards the genetic improvement of quinoa.

A bacterial community-based index to assess the ecological status of estuarine and coastal environments

Biotic indices for monitoring marine ecosystems are mostly based on the analysis of benthic macroinvertebrate communities. Due to their high sensitivity to pollution and fast response to environmental changes, bacterial assemblages could complement the information provided by benthic metazoan communities as indicators of human-induced impacts, but so far, this biological component has not been well explored for this purpose. Here we performed 16S rRNA gene amplicon sequencing to analyze the bacterial assemblage composition of 51 estuarine and coastal stations characterized by different environmental conditions and human-derived pressures. Using the relative abundance of putative indicator bacterial taxa, we developed a biotic index that is significantly correlated with a sediment quality index calculated on the basis of organic and inorganic compound concentrations. This new index based on bacterial assemblage composition can be a sensitive tool for providing a fast environmental assessment and allow a more comprehensive integrative ecosystem approach for environmental management.

Differential specificity between closely related corals and abundant Endozoicomonas endosymbionts across global scales

Reef-building corals are well regarded not only for their obligate association with endosymbiotic algae, but also with prokaryotic symbionts, the specificity of which remains elusive. To identify the central microbial symbionts of corals, their specificity across species and conservation over geographic regions, we sequenced partial SSU ribosomal RNA genes of Bacteria and Archaea from the common corals Stylophora pistillata and Pocillopora verrucosa across 28 reefs within seven major geographical regions. We demonstrate that both corals harbor Endozoicomonas bacteria as their prevalent symbiont. Importantly, catalyzed reporter deposition–fluorescence in situ hybridization (CARD–FISH) with Endozoicomonas-specific probes confirmed their residence as large aggregations deep within coral tissues. Using fine-scale genotyping techniques and single-cell genomics, we demonstrate that P. verrucosa harbors the same Endozoicomonas, whereas S. pistillata associates with geographically distinct genotypes. This specificity may be shaped by the different reproductive strategies of the hosts, potentially uncovering a pattern of symbiont selection that is linked to life history. Spawning corals such as P. verrucosa acquire prokaryotes from the environment. In contrast, brooding corals such as S. pistillata release symbiont-packed planula larvae, which may explain a strong regional signature in their microbiome. Our work contributes to the factors underlying microbiome specificity and adds detail to coral holobiont functioning.

Gene Expression Variation Resolves Species and Individual Strains among Coral-Associated Dinoflagellates within the Genus Symbiodinium

Reef-building corals depend on symbiotic mutualisms with photosynthetic dinoflagellates in the genus Symbiodinium. This large microalgal group comprises many highly divergent lineages (“Clades A–I”) and hundreds of undescribed species. Given their ecological importance, efforts have turned to genomic approaches to characterize the functional ecology of Symbiodinium. To date, investigators have only compared gene expression between representatives from separate clades—the equivalent of contrasting genera or families in other dinoflagellate groups—making it impossible to distinguish between clade-level and species-level functional differences. Here, we examined the transcriptomes of four species within one Symbiodinium clade (Clade B) at ∼20,000 orthologous genes, as well as multiple isoclonal cell lines within species (i.e., cultured strains). These species span two major adaptive radiations within Clade B, each encompassing both host-specialized and ecologically cryptic taxa. Species-specific expression differences were consistently enriched for photosynthesis-related genes, likely reflecting selection pressures driving niche diversification. Transcriptional variation among strains involved fatty acid metabolism and biosynthesis pathways. Such differences among individuals are potentially a major source of physiological variation, contributing to the functional diversity of coral holobionts composed of unique host–symbiont genotype pairings. Our findings expand the genomic resources available for this important symbiont group and emphasize the power of comparative transcriptomics as a method for studying speciation processes and interindividual variation in nonmodel organisms.