Environmental DNA metabarcoding for monitoring metazoan biodiversity in Antarctic nearshore ecosystems.

Antarctic benthic ecosystems support high biodiversity but their characterization is limited to a few well-studied areas, due to the extreme environment and remoteness making access and sampling difficult. Our aim was to compare water and sediment as sources of environmental DNA (eDNA) to better characterise Antarctic benthic communities and further develop practical approaches for DNA-based biodiversity assessment in remote environments. We used a cytochrome c oxidase subunit I (COI) metabarcoding approach to characterise metazoan communities in 26 nearshore sites across 12 locations in the Vestfold Hills (East Antarctica) based on DNA extracted from either sediment cores or filtered seawater. We detected a total of 99 metazoan species from 12 phyla across 26 sites, with similar numbers of species detected in sediment and water eDNA samples. However, significantly different communities were detected in the two sample types at sites where both were collected (i.e., where paired samples were available). For example, nematodes and echinoderms were more likely to be detected exclusively in sediment and water eDNA samples, respectively. eDNA from water and sediment core samples are complementary sample types, with epifauna more likely to be detected in water column samples and infauna in sediment. More reference DNA sequences are needed for infauna/meiofauna to increase the proportion of sequences and number of taxa that can be identified. Developing a better understanding of the temporal and spatial dynamics of eDNA at low temperatures would also aid interpretation of eDNA signals from polar environments. Our results provide a preliminary scan of benthic metazoan communities in the Vestfold Hills, with additional markers required to provide a comprehensive biodiversity survey. However, our study demonstrates the choice of sample type for eDNA studies of benthic ecosystems (sediment, water or both) needs to be carefully considered in light of the research or monitoring question of interest.

Bacterial epibiont communities of panmictic Antarctic krill are spatially structured.

Antarctic krill (Euphausia superba) are amongst the most abundant animals on Earth, with a circumpolar distribution in the Southern Ocean. Genetic and genomic studies have failed to detect any population structure for the species, suggesting a single panmictic population. However, the hyper-abundance of krill slows the rate of genetic differentiation, masking potential underlying structure. Here we use high-throughput sequencing of bacterial 16S rRNA genes to show that krill bacterial epibiont communities exhibit spatial structuring, driven mainly by distance rather than environmental factors, especially for strongly krill-associated bacteria. Estimating the ecological processes driving bacterial community turnover indicated this was driven by bacterial dispersal limitation increasing with geographic distance. Furthermore, divergent epibiont communities generated from a single krill swarm split between aquarium tanks under near-identical conditions suggests physical isolation in itself can cause krill-associated bacterial communities to diverge. Our findings show that Antarctic krill-associated bacterial communities are geographically structured, in direct contrast with the lack of structure observed for krill genetic and genomic data.

Capturing open ocean biodiversity: Comparing environmental DNA metabarcoding to the continuous plankton recorder.

Environmental DNA (eDNA) metabarcoding is emerging as a novel, objective tool for monitoring marine metazoan biodiversity. Zooplankton biodiversity in the vast open ocean is currently monitored through continuous plankton recorder (CPR) surveys, using ship-based bulk plankton sampling and morphological identification. We assessed whether eDNA metabarcoding (2 L filtered seawater) could capture similar Southern Ocean zooplankton biodiversity as conventional CPR bulk sampling (~1,500 L filtered seawater per CPR sample). We directly compared eDNA metabarcoding with (a) conventional morphological CPR sampling and (b) bulk DNA metabarcoding of CPR collected plankton (two transects for each comparison, 40 and 44 paired samples, respectively). A metazoan-targeted cytochrome c oxidase I (COI) marker was used to characterize species-level diversity. In the 2 L seawater eDNA samples, this marker amplified large amounts of non-metazoan picoplanktonic algae, but eDNA metabarcoding still detected up to 1.6 times more zooplankton species than morphologically analysed bulk CPR samples. COI metabarcoding of bulk DNA samples mostly avoided nonmetazoan amplifications and recovered more zooplankton species than eDNA metabarcoding. However, eDNA metabarcoding detected roughly two thirds of metazoan species and identified similar taxa contributing to community differentiation across the subtropical front separating transects. We observed a diurnal pattern in eDNA data for copepods which perform diel vertical migrations, indicating a surprisingly short temporal eDNA signal. Compared to COI, a eukaryote-targeted 18S ribosomal RNA marker detected a higher proportion, but lower diversity, of metazoans in eDNA. With refinement and standardization of methodology, eDNA metabarcoding could become an efficient tool for monitoring open ocean biodiversity.

The genomic basis of adaptation to calcareous and siliceous soils in Arabidopsis lyrata.

Edaphic conditions are important determinants of plant fitness. While much has been learnt in recent years about plant adaptation to heavy metal contaminated soils, the genomic basis underlying adaptation to calcareous and siliceous substrates remains largely unknown. We performed a reciprocal germination experiment and whole-genome resequencing in natural calcareous and siliceous populations of diploid Arabidopsis lyrata to test for edaphic adaptation and detect signatures of selection at loci associated with soil-mediated divergence. In parallel, genome scans on respective diploid ecotypes from the Arabidopsis arenosa species complex were undertaken, to search for shared patterns of adaptive genetic divergence. Soil ecotypes of A. lyrata display significant genotype-by-treatment responses for seed germination. Sequence (SNPs) and copy-number variants (CNVs) point towards loci involved in ion transport as the main targets of adaptive genetic divergence. Two genes exhibiting high differentiation among soil types in A. lyrata further share trans-specific single nucleotide polymorphisms with A. arenosa. This work applies experimental and genomic approaches to study edaphic adaptation in A. lyrata and suggests that physiological response to elemental toxicity and deficiency underlies the evolution of calcareous and siliceous ecotypes. The discovery of shared adaptive variation between sister species indicates that ancient polymorphisms contribute to adaptive evolution.

Antarctic Krill Are Reservoirs for Distinct Southern Ocean Microbial Communities.

Host-associated bacterial communities have received limited attention in polar habitats, but are likely to represent distinct nutrient-rich niches compared to the surrounding environment. Antarctic krill (Euphausia superba) are a super-abundant species with a circumpolar distribution, and the krill microbiome may make a substantial contribution to marine bacterial diversity in the Southern Ocean. We used high-throughput sequencing of the bacterial 16S ribosomal RNA gene to characterize bacterial diversity in seawater and krill tissue samples from four locations south of the Kerguelen Plateau, one of the most productive regions in the Indian Sector of the Southern Ocean. Krill-associated bacterial communities were distinct from those of the surrounding seawater, with different communities inhabiting the moults, digestive tract and faecal pellets, including several phyla not detected in the surrounding seawater. Digestive tissues from many individuals contained a potential gut symbiont (order: Mycoplasmoidales) shown to improve survival on a low quality diet in other crustaceans. Antarctic krill swarms thus influence Southern Ocean microbial communities not only through top-down grazing of eukaryotic cells and release of nutrients into the water column, but also by transporting distinct microbial assemblages horizontally via migration and vertically via sinking faecal pellets and moulted exuviae. Changes to Antarctic krill demographics or distribution through fishing pressure or climate-induced range shifts will also influence the composition and dispersal of Southern Ocean microbial communities.

Gene regulatory variation mediates flowering responses to vernalization along an altitudinal gradient in Arabidopsis.

Steep environmental gradients provide ideal settings for studies of potentially adaptive phenotypic and genetic variation in plants. The accurate timing of flowering is crucial for reproductive success and is regulated by several pathways, including the vernalization pathway. Among the numerous genes known to enable flowering in response to vernalization, the most prominent is FLOWERING LOCUS C (FLC). FLC and other genes of the vernalization pathway vary extensively among natural populations and are thus candidates for the adaptation of flowering time to environmental gradients such as altitude. We used 15 natural Arabidopsis (Arabidopsis thaliana) genotypes originating from an altitudinal gradient (800-2,700 m above sea level) in the Swiss Alps to test whether flowering time correlated with altitude under different vernalization scenarios. Additionally, we measured the expression of 12 genes of the vernalization pathway and its downstream targets. Flowering time correlated with altitude in a nonlinear manner for vernalized plants. Flowering time could be explained by the expression and regulation of the vernalization pathway, most notably by AGAMOUS LIKE19 (AGL19), FLOWERING LOCUS T (FT), and FLC. The expression of AGL19, FT, and VERNALIZATION INSENSITIVE3 was associated with altitude, and the regulation of MADS AFFECTING FLOWERING2 (MAF2) and MAF3 differed between low- and high-altitude genotypes. In conclusion, we found clinal variation across an altitudinal gradient both in flowering time and the expression and regulation of genes in the flowering time control network, often independent of FLC, suggesting that the timing of flowering may contribute to altitudinal adaptation.

Inheritance and reproductive consequences of floral anthocyanin deficiency in Silene dioica (Caryophyllaceae).

UNLABELLED: • PREMISE OF THE STUDY: Flower color is one of the most important traits for pollinator attraction. However, natural plant populations often harbor rare flower color variants resulting from mutations in biochemical pathways for floral pigment production. It is unclear how such mutations can persist because they can affect not only pollinator visitation but also plant fertility and performance.• METHODS: We collected rare white-flowered (anthocyanin-deficient) and common pink-flowered morphs of Silene dioica from natural populations in Switzerland. First- and second-generation hybrids between pink and white morphs, as well as backcrosses toward white morphs were produced, and the proportion of white-flowered offspring was determined. We compared seed siring ability and seed production between morphs using hand pollination experiments. Moreover, we scored the transfer of pollen analogues (fluorescent dyes) in a 50:50 array of the two morphs.• KEY RESULTS: The proportions of white-flowered plants in the offspring of our crosses were consistent with more than one recessive mutation as the cause of floral anthocyanin deficiency and further suggested a role of maternal effects for flower color. Seed siring ability and seed set did not differ significantly between pink and white morphs. Pollen transfer occurred preferentially within morphs.• CONCLUSIONS: Overall, our results imply that the white morph of S. dioica likely is caused by recessive mutations that do not impair reproduction. However, as this flower color polymorphism led to assortative mating in our experiment, it may represent standing genetic variation with the potential to contribute to evolutionary divergence.

Phenotypic effects of salt and heat stress over three generations in Arabidopsis thaliana.

Current and predicted environmental change will force many organisms to adapt to novel conditions, especially sessile organisms such as plants. It is therefore important to better understand how plants react to environmental stress and to what extent genotypes differ in such responses. It has been proposed that adaptation to novel conditions could be facilitated by heritable epigenetic changes induced by environmental stress, independent of genetic variation. Here we assessed phenotypic effects of heat and salt stress within and across three generations using four highly inbred Arabidopsis thaliana genotypes (Col, Cvi, Ler and Sha). Salt stress generally decreased fitness, but genotypes were differently affected, suggesting that susceptibility of A. thaliana to salt stress varies among genotypes. Heat stress at an early rosette stage had less detrimental effects but accelerated flowering in three out of four accessions. Additionally, we found three different modes of transgenerational effects on phenotypes, all harboring the potential of being adaptive: heat stress in previous generations induced faster rosette growth in Sha, both under heat and control conditions, resembling a tracking response, while in Cvi, the phenotypic variance of several traits increased, resembling diversified bet-hedging. Salt stress experienced in earlier generations altered plant architecture of Sha under salt but not control conditions, similar to transgenerational phenotypic plasticity. However, transgenerational phenotypic effects depended on the type of stress as well as on genotype, suggesting that such effects may not be a general response leading to adaptation to novel environmental conditions in A. thaliana.

Environmental heat and salt stress induce transgenerational phenotypic changes in Arabidopsis thaliana.

Plants that can adapt their phenotype may be more likely to survive changing environmental conditions. Heritable epigenetic variation could provide a way to rapidly adapt to such changes. Here we tested whether environmental stress induces heritable, potentially adaptive phenotypic changes independent of genetic variation over few generations in Arabidopsis thaliana. We grew two accessions (Col-0, Sha-0) of A. thaliana for three generations under salt, heat and control conditions and tested for induced heritable phenotypic changes in the fourth generation (G4) and in reciprocal F1 hybrids generated in generation three. Using these crosses we further tested whether phenotypic changes were maternally or paternally transmitted. In generation five (G5), we assessed whether phenotypic effects persisted over two generations in the absence of stress. We found that exposure to heat stress in previous generations accelerated flowering under G4 control conditions in Sha-0, but heritable effects disappeared in G5 after two generations without stress exposure. Previous exposure to salt stress increased salt tolerance in one of two reciprocal F1 hybrids. Transgenerational effects were maternally and paternally inherited. Lacking genetic variability, maternal and paternal inheritance and reversibility of transgenerational effects together indicate that stress can induce heritable, potentially adaptive phenotypic changes, probably through epigenetic mechanisms. These effects were strongly dependent on plant genotype and may not be a general response to stress in A. thaliana.

Admixture facilitates adaptation from standing variation in the European aspen (Populus tremula L.), a widespread forest tree.

Adaptation to new environments can start from new mutations or from standing variation already present in natural populations. Whether admixture constrains or facilitates adaptation from standing variation is largely unknown, especially in ecological keystone or foundation species. We examined patterns of neutral and adaptive population divergence in Populus tremula L., a widespread forest tree, using mapped molecular genetic markers. We detected the genetic signature of postglacial admixture between a Western and an Eastern lineage of P. tremula in Scandinavia, an area suspected to represent a zone of postglacial contact for many species of animals and plants. Stringent divergence-based neutrality tests provided clear indications for locally varying selection at the European scale. Six of 12 polymorphisms under selection were located less than 1 kb away from the nearest gene predicted by the Populus trichocarpa genome sequence. Few of these loci exhibited a signature of 'selective sweeps' in diversity-based tests, which is to be expected if adaptation occurs primarily from standing variation. In Scandinavia, admixture explained genomic patterns of ancestry and the nature of clinal variation and strength of selection for bud set, a phenological trait of great adaptive significance in temperate trees, measured in a common garden trial. Our data provide a hitherto missing direct link between past range shifts because of climatic oscillations, and levels of standing variation currently available for selection and adaptation in a terrestrial foundation species.