Genomic divergence between the migratory and stationary ecotypes of Atlantic cod

Atlantic cod displays a range of phenotypic and genotypic variations, which includes the differentiation into coastal stationary and offshore migratory types of cod that co‐occur in several parts of its distribution range and are often sympatric on the spawning grounds. Differentiation of these ecotypes may involve both historical separation and adaptation to ecologically distinct environments, the genetic basis of which is now beginning to be unravelled. Genomic analyses based on recent sequencing advances are able to document genomic divergence in more detail and may facilitate the exploration of causes and consequences of genome‐wide patterns. We examined genomic divergence between the stationary and migratory types of cod in the Northeast Atlantic, using next‐generation sequencing of pooled DNA from each of two population samples. Sequence data was mapped to the published cod genome sequence, arranged in more than 6000 scaffolds (611 Mb). We identified 25 divergent scaffolds (26 Mb) with a higher than average gene density, against a backdrop of overall moderate genomic differentiation. Previous findings of localized genomic divergence in three linkage groups were confirmed, including a large (15 Mb) genomic region, which seems to be uniquely involved in the divergence of migratory and stationary cod. The results of the pooled sequencing approach support and extend recent findings based on single‐nucleotide polymorphism markers and suggest a high degree of reproductive isolation between stationary and migratory cod in the North‐east Atlantic.     

Behavioral responses of Atlantic cod to sea temperature changes

Understanding responses of marine species to temperature variability is essential to predict impacts of future climate change in the oceans. Most ectotherms are expected to adjust their behavior to avoid extreme temperatures and minimize acute changes in body temperature. However, measuring such behavioral plasticity in the wild is challenging. Combining 4 years of telemetry‐derived behavioral data on juvenile and adult (30–80 cm) Atlantic cod (Gadus morhua), and in situ ocean temperature measurements, we found a significant effect of sea temperature on cod depth use and activity level in coastal Skagerrak. During summer, cod were found in deeper waters when sea surface temperature increased. Further, this effect of temperature was stronger on larger cod. Diel vertical migration, which consists in a nighttime rise to shallow feeding habitats, was stronger among smaller cod. As surface temperature increased beyond ~15°C, their vertical migration was limited to deeper waters. In addition to larger diel vertical migrations, smaller cod were more active and travelled larger distances compared to larger specimens. Cold temperatures during winter tended, however, to reduce the magnitude of diel vertical migrations, as well as the activity level and distance moved by those smaller individuals. Our findings suggest that future and ongoing rises in sea surface temperature may increasingly deprive cod in this region from shallow feeding areas during summer, which may be detrimental for local populations of the species.     

Climate change and fishing: a century of shifting distribution in North Sea cod

Globally, spatial distributions of fish stocks are shifting but although the role of climate change in range shifts is increasingly appreciated, little remains known of the likely additional impact that high levels of fishing pressure might have on distribution. For North Sea cod, we show for the first time and in great spatial detail how the stock has shifted its distribution over the past 100 years. We digitized extensive historical fisheries data from paper charts in UK government archives and combined these with contemporary data to a time‐series spanning 1913–2012 (excluding both World Wars). New analysis of old data revealed that the current distribution pattern of cod – mostly in the deeper, northern‐ and north‐easternmost parts of the North Sea – is almost opposite to that during most of the Twentieth Century – mainly concentrated in the west, off England and Scotland. Statistical analysis revealed that the deepening, northward shift is likely attributable to warming; however, the eastward shift is best explained by fishing pressure, suggestive of significant depletion of the stock from its previous stronghold, off the coasts of England and Scotland. These spatial patterns were confirmed for the most recent 3½ decades by data from fisheries‐independent surveys, which go back to the 1970s. Our results demonstrate the fundamental importance of both climate change and fishing pressure for our understanding of changing distributions of commercially exploited fish.     

Divergent responses of Atlantic cod to ocean acidification and food limitation

In order to understand the effect of global change on marine fishes, it is imperative to quantify the effects on fundamental parameters such as survival and growth. Larval survival and recruitment of the Atlantic cod (Gadus morhua) were found to be heavily impaired by end‐of‐century levels of ocean acidification. Here, we analysed larval growth among 35–36 days old surviving larvae, along with organ development and ossification of the skeleton. We combined CO₂ treatments (ambient: 503 µatm, elevated: 1,179 µatm) with food availability in order to evaluate the effect of energy limitation in addition to the ocean acidification stressor. As expected, larval size (as a proxy for growth) and skeletogenesis were positively affected by high food availability. We found significant interactions between acidification and food availability. Larvae fed ad libitum showed little difference in growth and skeletogenesis due to the CO₂ treatment. Larvae under energy limitation were significantly larger and had further developed skeletal structures in the elevated CO₂ treatment compared to the ambient CO₂ treatment. However, the elevated CO₂ group revealed impairments in critically important organs, such as the liver, and had comparatively smaller functional gills indicating a mismatch between size and function. It is therefore likely that individual larvae that had survived acidification treatments will suffer from impairments later during ontogeny. Our study highlights important allocation trade‐off between growth and organ development, which is critically important to interpret acidification effects on early life stages of fish.     

Forecasting future recruitment success for Atlantic cod in the warming and acidifying Barents Sea

Productivity of marine fish stocks is known to be affected by environmental and ecological drivers, and global climate change is anticipated to alter recruitment success of many stocks. While the direct effects of environmental drivers on fish early life stage survival can be quantified experimentally, indirect effects in marine ecosystems and the role of adaptation are still highly uncertain. We developed an integrative model for the effects of ocean warming and acidification on the early life stages of Atlantic cod in the Barents Sea, termed SCREI (Simulator of Cod Recruitment under Environmental Influences). Experimental results on temperature and CO₂ effects on egg fertilization, egg and larval survival and development times are incorporated. Calibration using empirical time series of egg production, temperature, food and predator abundance reproduces age‐0 recruitment over three decades. We project trajectories of recruitment success under different scenarios and quantify confidence limits based on variation in experiments. A publicly accessible web version of the SCREI model can be run under www.oceanchange.uni-bremen.de/ ;SCREI. Severe reductions in average age‐0 recruitment success of Barents Sea cod are projected under uncompensated warming and acidification toward the middle to end of this century. Although high population stochasticity was found, considerable rates of evolutionary adaptation to acidification and shifts in organismal thermal windows would be needed to buffer impacts on recruitment. While increases in food availability may mitigate short‐term impacts, an increase in egg production achieved by stock management could provide more long‐term safety for cod recruitment success. The SCREI model provides a novel integration of multiple driver effects in different life stages and enables an estimation of uncertainty associated with interindividual and ecological variation. The model thus helps to advance toward an improved empirical foundation for quantifying climate change impacts on marine fish recruitment, relevant for ecosystem‐based assessments of marine systems under climate change.     

Analysis of ROH patterns in the Noriker horse breed reveals signatures of selection for coat color and body size

Overlapping runs of homozygosity (ROH islands) shared by the majority of a population are hypothesized to be the result of selection around a target locus. In this study we investigated the impact of selection for coat color within the Noriker horse on autozygosity and ROH patterns. We analyzed overlapping homozygous regions (ROH islands) for gene content in fragments shared by more than 50% of horses. Long‐term assortative mating of chestnut horses and the small effective population size of leopard spotted and tobiano horses resulted in higher mean genome‐wide ROH coverage (S_ROH) within the range of 237.4–284.2 Mb, whereas for bay, black and roan horses, where rotation mating is commonly applied, lower autozygosity (S_ROH from 176.4–180.0 Mb) was determined. We identified seven common ROH islands considering all Noriker horses from our dataset. Specific islands were documented for chestnut, leopard spotted, roan and bay horses. The ROH islands contained, among others, genes associated with body size (ZFAT, LASP1 and LCORL/NCAPG), coat color (MC1R in chestnut and the factor PATN1 in leopard spotted horses) and morphogenesis (HOXB cluster in all color strains except leopard spotted horses). This study demonstrates that within a closed population sharing the same founders and ancestors, selection on a single phenotypic trait, in this case coat color, can result in genetic fragmentation affecting levels of autozygosity and distribution of ROH islands and enclosed gene content.     

Metagenomic analysis of historical herbarium specimens reveals a postmortem microbial community

Advances in DNA extraction and next‐generation sequencing have made a vast number of historical herbarium specimens available for genomic investigation. These specimens contain not only genomic information from the individual plants themselves, but also from associated microorganisms such as bacteria and fungi. These microorganisms may have colonized the living plant (e.g., pathogens or host‐associated commensal taxa) or may result from postmortem colonization that may include decomposition processes or contamination during sample handling. Here we characterize the metagenomic profile from shotgun sequencing data from herbarium specimens of two widespread plant species (Ambrosia artemisiifolia and Arabidopsis thaliana) collected up to 180 years ago. We used blast searching in combination with megan and were able to infer the metagenomic community even from the oldest herbarium sample. Through comparison with contemporary plant collections, we identify three microbial species that are nearly exclusive to herbarium specimens, including the fungus Alternaria alternata, which can comprise up to 7% of the total sequencing reads. This species probably colonizes the herbarium specimens during preparation for mounting or during storage. By removing the probable contaminating taxa, we observe a temporal shift in the metagenomic composition of the invasive weed Am. artemisiifolia. Our findings demonstrate that it is generally possible to use herbarium specimens for metagenomic analyses, but that the results should be treated with caution, as some of the identified species may be herbarium contaminants rather than representing the natural metagenomic community of the host plant.     

Long‐term change in a behavioural trait: truncated spawning distribution and demography in Northeast Arctic cod

Harvesting may be a potent driver of demographic change and contemporary evolution, which both may have great impacts on animal populations. Research has focused on changes in phenotypic traits that are easily quantifiable and for which time series exist, such as size, age, sex, or gonad size, whereas potential changes in behavioural traits have been under‐studied. Here, we analyse potential drivers of long‐term changes in a behavioural trait for the Northeast Arctic stock of Atlantic cod Gadus morhua, namely choice of spawning location. For 104 years (1866–1969), commercial catches were recorded annually and reported by county along the Norwegian coast. During this time period, spawning ground distribution has fluctuated with a trend towards more northerly spawning. Spawning location is analysed against a suite of explanatory factors including climate, fishing pressure, density dependence, and demography. We find that demography (age or age at maturation) had the highest explanatory power for variation in spawning location, while climate had a limited effect below statistical significance. As to potential mechanisms, some effects of climate may act through demography, and explanatory variables for demography may also have absorbed direct evolutionary change in migration distance for which proxies were unavailable. Despite these caveats, we argue that fishing mortality, either through demographic or evolutionary change, has served as an effective driver for changing spawning locations in cod, and that additional explanatory factors related to climate add no significant information.     

Genetic signatures through space,time and multiple disturbances in a ubiquitous brooding coral

The predominance of self‐recruitment in many reef‐building corals has fundamental and complex consequences for their genetic diversity, population persistence and responses to climate change. Knowledge of genetic structure over local scales needs to be placed within a broad spatial context, and also integrated with genetic monitoring through time to disentangle these consequences. Here, we examined patterns of genetic diversity over multiple spatio‐temporal scales across tropical Australia in the ubiquitous brooding coral, Seriatopora hystrix. We also analysed complimentary environmental and demographic data to elucidate the seascape drivers of these patterns. Large genetic differences were detected between the east vs. west coasts of Australia. In northwest Australia, geographic differentiation dominated genetic structure over multiple scales. However, three sympatric lineages were detected at the largest offshore reef system (Scott Reef). Similar to the differences observed among putative species in eastern Australia, these lineages were associated with different levels of wave exposure. Local genetic structure within the Scott Reef system was relatively stable over 10 years, but temporal differences were observed that reflected small but important genetic changes over a few generations during recovery after severe bleaching. These results highlight the importance of self‐recruitment together with occasional longer distance connectivity for the persistence of a metapopulation across spatially and temporally variable environments. Our multidimensional research provides a foundation for further long‐term genetic monitoring to inform conservation strategies and highlights that sampling scales, ecological effects and cryptic diversity are important considerations to develop realistic understanding of the evolutionary resilience of corals.     

Detecting selection along environmental gradients: analysis of eight methods and their effectiveness for outbreeding and selfing populations

Thanks to genome‐scale diversity data, present‐day studies can provide a detailed view of how natural and cultivated species adapt to their environment and particularly to environmental gradients. However, due to their sensitivity, up‐to‐date studies might be more sensitive to undocumented demographic effects such as the pattern of migration and the reproduction regime. In this study, we provide guidelines for the use of popular or recently developed statistical methods to detect footprints of selection. We simulated 100 populations along a selective gradient and explored different migration models, sampling schemes and rates of self‐fertilization. We investigated the power and robustness of eight methods to detect loci potentially under selection: three designed to detect genotype–environment correlations and five designed to detect adaptive differentiation (based on F_ST or similar measures). We show that genotype–environment correlation methods have substantially more power to detect selection than differentiation‐based methods but that they generally suffer from high rates of false positives. This effect is exacerbated whenever allele frequencies are correlated, either between populations or within populations. Our results suggest that, when the underlying genetic structure of the data is unknown, a number of robust methods are preferable. Moreover, in the simulated scenario we used, sampling many populations led to better results than sampling many individuals per population. Finally, care should be taken when using methods to identify genotype–environment correlations without correcting for allele frequency autocorrelation because of the risk of spurious signals due to allele frequency correlations between populations.     

Comparative mitochondrial genome analysis reveals the evolutionary rearrangement mechanism in Brassica

The genus Brassica has many species that are important for oil, vegetable and other food products. Three mitochondrial genome types (mitotype) originated from its common ancestor. In this paper, a B. nigra mitochondrial main circle genome with 232,407 bp was generated through de novo assembly. Synteny analysis showed that the mitochondrial genomes of B. rapa and B. oleracea had a better syntenic relationship than B. nigra. Principal components analysis and development of a phylogenetic tree indicated maternal ancestors of three allotetraploid species in Us triangle of Brassica. Diversified mitotypes were found in allotetraploid B. napus, in which napus‐type B. napus was derived from B. oleracea, while polima‐type B. napus was inherited from B. rapa. In addition, the mitochondrial genome of napus‐type B. napus was closer to botrytis‐type than capitata‐type B. oleracea. The sub‐stoichiometric shifting of several mitochondrial genes suggested that mitochondrial genome rearrangement underwent evolutionary selection during domestication and/or plant breeding. Our findings clarify the role of diploid species in the maternal origin of allotetraploid species in Brassica and suggest the possibility of breeding selection of the mitochondrial genome.     

Genome‐wide detection of signatures of selection in Korean Hanwoo cattle

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26–0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise F_ST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome‐wide smoothed F_ST, 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty‐four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.     

Population structure of Atlantic mackerel inferred from RAD‐seq‐derived SNP markers: effects of sequence clustering parameters and hierarchical SNP selection

Restriction‐site‐associated DNA sequencing (RAD‐seq) and related methods are revolutionizing the field of population genomics in nonmodel organisms as they allow generating an unprecedented number of single nucleotide polymorphisms (SNPs) even when no genomic information is available. Yet, RAD‐seq data analyses rely on assumptions on nature and number of nucleotide variants present in a single locus, the choice of which may lead to an under‐ or overestimated number of SNPs and/or to incorrectly called genotypes. Using the Atlantic mackerel (Scomber scombrus L.) and a close relative, the Atlantic chub mackerel (Scomber colias), as case study, here we explore the sensitivity of population structure inferences to two crucial aspects in RAD‐seq data analysis: the maximum number of mismatches allowed to merge reads into a locus and the relatedness of the individuals used for genotype calling and SNP selection. Our study resolves the population structure of the Atlantic mackerel, but, most importantly, provides insights into the effects of alternative RAD‐seq data analysis strategies on population structure inferences that are directly applicable to other species.     

Stocking activities for the Arctic charr in Lake Geneva: Genetic effects in space and time

Artificial stocking practices are widely used by resource managers worldwide, in order to sustain fish populations exploited by both recreational and commercial activities, but their benefits are controversial. Former practices involved exotic strains, although current programs rather consider artificial breeding of local fishes (supportive breeding). Understanding the complex genetic effects of these management strategies is an important challenge with economic and conservation implications, especially in the context of population declines. In this study, we focus on the declining Arctic charr (Salvelinus alpinus) population from Lake Geneva (Switzerland and France), which has initially been restocked with allochtonous fishes in the early eighties, followed by supportive breeding. In this context, we conducted a genetic survey to document the evolution of the genetic diversity and structure throughout the last 50 years, before and after the initiation of hatchery supplementation, using contemporary and historical samples. We show that the introduction of exotic fishes was associated with a genetic bottleneck in the 1980–1990s, a break of Hardy–Weinberg Equilibrium (HWE), a reduction in genetic diversity, an increase in genetic structure among spawning sites, and a change in their genetic composition. Together with better environmental conditions, three decades of subsequent supportive breeding using local fishes allowed to re‐establish HWE and the initial levels of genetic variation. However, current spawning sites have not fully recovered their original genetic composition and were extensively homogenized across the lake. Our study demonstrates the drastic genetic consequences of different restocking tactics in a comprehensive spatiotemporal framework and suggests that genetic alteration by nonlocal stocking may be partly reversible through supportive breeding. We recommend that conservation‐based programs consider local diversity and implement adequate protocols to limit the genetic homogenization of this Arctic charr population.     

Mechanistic insights into the effects of climate change on larval cod

Understanding the biophysical mechanisms that shape variability in fisheries recruitment is critical for estimating the effects of climate change on fisheries. In this study, we used an Earth System Model (ESM) and a mechanistic individual‐based model (IBM) for larval fish to analyze how climate change may impact the growth and survival of larval cod in the North Atlantic. We focused our analysis on five regions that span the current geographical range of cod and are known to contain important spawning populations. Under the SRES A2 (high emissions) scenario, the ESM‐projected surface ocean temperatures are expected to increase by >1 °C for 3 of the 5 regions, and stratification is expected to increase at all sites between 1950–1999 and 2050–2099. This enhanced stratification is projected to decrease large (>5 μm ESD) phytoplankton productivity and mesozooplankton biomass at all 5 sites. Higher temperatures are projected to increase larval metabolic costs, which combined with decreased food resources will reduce larval weight, increase the probability of larvae dying from starvation and increase larval exposure to visual and invertebrate predators at most sites. If current concentrations of piscivore and invertebrate predators are maintained, larval survival is projected to decrease at all five sites by 2050–2099. In contrast to past observed responses to climate variability in which warm anomalies led to better recruitment in cold‐water stocks, our simulations indicated that reduced prey availability under climate change may cause a reduction in larval survival despite higher temperatures in these regions. In the lower prey environment projected under climate change, higher metabolic costs due to higher temperatures outweigh the advantages of higher growth potential, leading to negative effects on northern cod stocks. Our results provide an important first large‐scale assessment of the impacts of climate change on larval cod in the North Atlantic.     

Signatures of selection in the Iberian honey bee (Apis mellifera iberiensis) revealed by a genome scan analysis of single nucleotide polymorphisms

Understanding the genetic mechanisms of adaptive population divergence is one of the most fundamental endeavours in evolutionary biology and is becoming increasingly important as it will allow predictions about how organisms will respond to global environmental crisis. This is particularly important for the honey bee, a species of unquestionable ecological and economical importance that has been exposed to increasing human‐mediated selection pressures. Here, we conducted a single nucleotide polymorphism (SNP)‐based genome scan in honey bees collected across an environmental gradient in Iberia and used four F_ST‐based outlier tests to identify genomic regions exhibiting signatures of selection. Additionally, we analysed associations between genetic and environmental data for the identification of factors that might be correlated or act as selective pressures. With these approaches, 4.4% (17 of 383) of outlier loci were cross‐validated by four F_ST‐based methods, and 8.9% (34 of 383) were cross‐validated by at least three methods. Of the 34 outliers, 15 were found to be strongly associated with one or more environmental variables. Further support for selection, provided by functional genomic information, was particularly compelling for SNP outliers mapped to different genes putatively involved in the same function such as vision, xenobiotic detoxification and innate immune response. This study enabled a more rigorous consideration of selection as the underlying cause of diversity patterns in Iberian honey bees, representing an important first step towards the identification of polymorphisms implicated in local adaptation and possibly in response to recent human‐mediated environmental changes.