Low‐coverage whole‐genome sequencing reveals molecular markers for spawning season and sex identification in Gulf of Maine Atlantic cod (Gadus morhua,Linnaeus 1758)

Atlantic cod (Gadus morhua, Linnaeus 1758) in the western Gulf of Maine are managed as a single stock despite several lines of evidence supporting two spawning groups (spring and winter) that overlap spatially, while exhibiting seasonal spawning isolation. Low‐coverage whole‐genome sequencing was used to evaluate the genomic population structure of Atlantic cod spawning groups in the western Gulf of Maine and Georges Bank using 222 individuals collected over multiple years. Results indicated low total genomic differentiation, while also showing strong differentiation between spring and winter‐spawning groups at specific regions of the genome. Guided regularized random forest and ranked F _ST methods were used to select panels of single nucleotide polymorphisms (SNPs) that could reliably distinguish spring and winter‐spawning Atlantic cod (88.5% assignment rate), as well as males and females (95.0% assignment rate) collected in the western Gulf of Maine. These SNP panels represent a valuable tool for fisheries research and management of Atlantic cod in the western Gulf of Maine that will aid investigations of stock production and support accuracy of future assessments.     

Preliminary data on the variability of three microsatellite loci in Pacific <Emphasis Type="Italic">Gadus macrocephalus</Emphasis> and Atlantic <Emphasis Type="Italic">G. morhua</Emphasis> cod (Gadidae)

Variability of microsatellite DNA loci Gmo3, Gmo34, and Gmo35 is studied in samples of Pacific cod Gadus macrocephalus and Atlantic cod G. morhua. The results show high values of identity of the samples within the North Pacific basin (0.9766–0.9924) and within the Northeast Atlantic basin (0.9580). Based on the pairwise assessment of genetic differentiation, the F _ST values are significantly different in all variants between the samples of Pacific and Atlantic cod (F _ST = 0.5235–0.6719, p < 0.001). Within the basins, the significant differences in the frequencies of main alleles are revealed in the loci Gmo3 and Gmo34 for the samples from the Pacific and Atlantic oceans, respectively.     

Haplotype analysis of the mitochondrial DNA d‐loop region reveals the maternal origin and historical dynamics among the indigenous goat populations in east and west of the Democratic Republic of Congo

This study aimed at assessing haplotype diversity and population dynamics of three Congolese indigenous goat populations that included Kasai goat (KG), small goat (SG), and dwarf goat (DG) of the Democratic Republic of Congo (DRC). The 1169 bp d‐loop region of mitochondrial DNA (mtDNA) was sequenced for 339 Congolese indigenous goats. The total length of sequences was used to generate the haplotypes and evaluate their diversities, whereas the hypervariable region (HVI, 453 bp) was analyzed to define the maternal variation and the demographic dynamic. A total of 568 segregating sites that generated 192 haplotypes were observed from the entire d‐loop region (1169 bp d‐loop). Phylogenetic analyses using reference haplotypes from the six globally defined goat mtDNA haplogroups showed that all the three Congolese indigenous goat populations studied clustered into the dominant haplogroup A, as revealed by the neighbor‐joining (NJ) tree and median‐joining (MJ) network. Nine haplotypes were shared between the studied goats and goat populations from Pakistan (1 haplotype), Kenya, Ethiopia and Algeria (1 haplotype), Zimbabwe (1 haplotype), Cameroon (3 haplotypes), and Mozambique (3 haplotypes). The population pairwise analysis (F_ST ) indicated a weak differentiation between the Congolese indigenous goat populations. Negative and significant (p‐value <.05) values for Fu''s Fs (−20.418) and Tajima''s (−2.189) tests showed the expansion in the history of the three Congolese indigenous goat populations. These results suggest a weak differentiation and a single maternal origin for the studied goats. This information will contribute to the improvement of the management strategies and long‐term conservation of indigenous goats in DRC.     

Phylogeography and genetic diversity of the widespread katydid Ducetia japonica (Thunberg, 1815) across China

Habitat fragmentation can lower migration rates and genetic connectivity among remaining populations of native species. Ducetia japonica is one of the most widespread katydids in China, but little is known about its genetic structure and phylogeographic distribution. We combined the five‐prime region of cytochrome c oxidase subunit I (COI‐5P), 11 newly developed microsatellite loci coupled with an ecological niche model (ENM) to examine the genetic diversity and population structure of D. japonica in China and beyond to Laos and Singapore. Both Bayesian inference (BI) and haplotype network methods revealed six mitochondrial COI‐5P lineages. The distribution of COI‐5P haplotypes may not demonstrate significant phylogeographic structure (N _ST > G _ST, p > .05). The STRUCTURE analysis based on microsatellite data also revealed six genetic clusters, but discordant with those obtained from COI‐5P haplotypes. For both COI‐5P and microsatellite data, Mantel tests revealed a significant positive correlation between geographic and genetic distances in mainland China. Bayesian skyline plot (BSP) analyses indicated that the population size of D. japonica''s three major mitochondrial COI‐5P lineages were seemingly not affected by last glacial maximum (LGM, 0.015–0.025 Mya). The ecological niche models showed that the current distribution of D. japonica was similar to the species’ distribution during the LGM period and only slightly extended in northern China. Further phylogeographic studies based on more extensive sampling are needed to identify specific locations of glacial refugia in northern China.     

Resource‐driven colonization by cod in a high Arctic food web

1. Climate change is commonly associated with many species redistributions and the influence of other factors may be marginalized, especially in the rapidly warming Arctic.
2. The Barents Sea, a high latitude large marine ecosystem in the Northeast Atlantic has experienced above‐average temperatures since the mid‐2000s with divergent bottom temperature trends at subregional scales.
3. Concurrently, the Barents Sea stock of Atlantic cod Gadus morhua, one of the most important commercial fish stocks in the world, increased following a large reduction in fishing pressure and expanded north of 80°N.
4. We examined the influence of food availability and temperature on cod expansion using a comprehensive data set on cod stomach fullness stratified by subregions characterized by divergent temperature trends. We then tested whether food availability, as indexed by cod stomach fullness, played a role in cod expansion in subregions that were warming, cooling, or showed no trend.
5. The greatest increase in cod occupancy occurred in three northern subregions with contrasting temperature trends. Cod apparently benefited from initial high food availability in these regions that previously had few large‐bodied fish predators.
6. The stomach fullness in the northern subregions declined rapidly after a few years of high cod abundance, suggesting that the arrival of cod caused a top‐down effect on the prey base. Prolonged cod residency in the northern Barents Sea is, therefore, not a certainty.

     

Population Structure of Pacific Cod <Emphasis Type="Italic">Gadus macrocephalus</Emphasis> in the Southern Part of the Range Based on the Microsatellite Analyses

The population structure of Pacific cod Gadus macrocephalus in the southern part of the range and adjacent regions is studied on the basis of the results of microsatellite analyses. Collected data indicate heterogeneity of this species population within the studied area. According to the obtained F_ST values, Pacific cod from the waters of the Republic of Korea (Yellow Sea side) and northwestern part of the Sea of Okhotsk significantly differ from all other studied regions (Table 4). Significant differentiation was also revealed between samples from the waters of the Tatar Strait and all other regions except for South Kurils Pacific cod (both Sea of Okhotsk and Pacific Ocean sides). These two latter sample collections were similar to each other as well. A low level of differentiation was also shown for the Peter the Great Bay and the East Sea/Sea of Japan waters of the Republic of Korea.     

Population structure of chum salmon and selection on the markers collected for stock identification

Genetic stock identification (GSI) is a major management tool of Pacific salmon (Oncorhynchus Spp.) that has provided rich genetic baseline data of allozymes, microsatellites, and single‐nucleotide polymorphisms (SNPs) across the Pacific Rim. Here, we analyzed published data sets for adult chum salmon (Oncorhynchus keta), namely 10 microsatellites, 53 SNPs, and a mitochondrial DNA locus (mtDNA3, control region, and NADH‐3 combined) in samples from 495 locations in the same distribution range (n = 61,813). TreeMix analysis of the microsatellite loci identified the greatest convergence toward Japanese/Korean populations and suggested two admixture events from Japan/Korea to Russia and the Alaskan Peninsula. The SNPs had been purposively collected from rapidly evolving genes to increase the power of GSI. The largest expected heterozygosity was observed in Japanese/Korean populations for microsatellites, whereas it was largest in Western Alaskan populations for SNPs, reflecting the SNP discovery process. A regression of SNP population structures on those of microsatellites indicated the selection of the SNP loci according to deviations from the predicted structures. Specifically, we matched the sampling locations of the SNPs with those of the microsatellites and performed regression analyses of SNP allele frequencies on a 2‐dimensional scaling (MDS) of matched locations obtained from microsatellite pairwise F _ST values. The MDS first axis indicated a latitudinal cline in American and Russian populations, whereas the second axis showed differentiation of Japanese/Korean populations. The top five outlier SNPs included mtDNA3, U502241 (unknown), GnRH373, ras1362, and TCP178, which were identified by principal component analysis. We summarized the functions of 53 nuclear genes surrounding SNPs and the mtDNA3 locus by referring to a gene database system and propose how they may influence the fitness of chum salmon.     

Variability of microsatellite loci of Greenland cod <Emphasis Type="Italic">Gadus ogac</Emphasis> Richardson 1836: Comparison with other species of <Emphasis Type="Italic">Gadus</Emphasis> genus (Gadidae)

Comparative analysis of variability of seven microsatellite loci—Gmo3, Gmo-G12, Gmo-G18, Gmo19, Gmo34, Gmo35 and Pgmo32—was performed for the Greenland cod Gadus ogac, Pacific cod G. macrocephalus, Atlantic cod G. morhua, and White Sea cod G. morhua marisalbi. High genetic identity was observed between the Greenland cod and Pacific cod (I = 0.9520). Pair analysis of genetic differentiation was performed on the studied microsatellite loci according to θ (analogue of F _ST). The Greenland cod differed significantly from the Pacific, Atlantic, and the White Sea cod; however, the differentiation level varied. The lowest value was observed for the pair Greenland cod-Pacific cod (0.123), and the highest levels were registered for the pairs Greenland cod-Atlantic cod (0.605) and Greenland cod-White Sea cod (0.535).     

Worldwide Distribution of the MYH9 Kidney Disease Susceptibility Alleles and Haplotypes: Evidence of Historical Selection in Africa

MYH9 was recently identified as renal susceptibility gene (OR 3–8, p<10⁻⁸) for major forms of kidney disease disproportionately affecting individuals of African descent. The risk haplotype (E-1) occurs at much higher frequencies in African Americans (≥60%) than in European Americans (<4%), revealing a genetic basis for a major health disparity. The population distributions of MYH9 risk alleles and the E-1 risk haplotype and the demographic and selective forces acting on the MYH9 region are not well explored. We reconstructed MYH9 haplotypes from 4 tagging single nucleotide polymorphisms (SNPs) spanning introns 12–23 using available data from HapMap Phase II, and by genotyping 938 DNAs from the Human Genome Diversity Panel (HGDP). The E-1 risk haplotype followed a cline, being most frequent within sub-Saharan African populations (range 50–80%), less frequent in populations from the Middle East (9–27%) and Europe (0–9%), and rare or absent in Asia, the Americas, and Oceania. The fixation indexes (F_ST) for pairwise comparisons between the risk haplotypes for continental populations were calculated for MYH9 haplotypes; F_ST ranged from 0.27–0.40 for Africa compared to other continental populations, possibly due to selection. Uniquely in Africa, the Yoruba population showed high frequency extended haplotype length around the core risk allele (C) compared to the alternative allele (T) at the same locus (rs4821481, iHs = 2.67), as well as high population differentiation (F_{ST(CEU vs. YRI)} = 0.51) in HapMap Phase II data, also observable only in the Yoruba population from HGDP (F_ST = 0.49), pointing to an instance of recent selection in the genomic region. The population-specific divergence in MYH9 risk allele frequencies among the world''s populations may prove important in risk assessment and public health policies to mitigate the burden of kidney disease in vulnerable populations.     

Mitochondrial DNA polymorphism of atlantic cod of the Barents and White seas

High level of intraspecific polymorphism of the mtDNA cytochrome b gene in the Atlantic cod Gadus morhua L. from the northeastern part of the species range (Barents and White seas) is detected. The absence of genetic differentiation between the samples from different areas of the Norwegian and Barents seas is shown, and the fact that these samples belong to the same population of Northeast Arctic cod is confirmed. Significant differences in the frequences of mitochondrial DNA haplotypes between the cod from the Barents and White seas, as well as from the northwestern part of the range, are found. The cod Gadus morhua kildinensis Derjugin, 1920 from Mogilnoye Lake and Gadus morhua marisalbi Derjugin, 1920 from the White Sea were revealed to have common haplotypes with the cod from the Barents Sea, confirming the recent origin of these forms from the Atlantic cod G. morhua.     

Detecting short spatial scale local adaptation and epistatic selection in climate‐related candidate genes in European beech (Fagus sylvatica) populations

Detecting signatures of selection in tree populations threatened by climate change is currently a major research priority. Here, we investigated the signature of local adaptation over a short spatial scale using 96 European beech (Fagus sylvatica L.) individuals originating from two pairs of populations on the northern and southern slopes of Mont Ventoux (south‐eastern France). We performed both single and multilocus analysis of selection based on 53 climate‐related candidate genes containing 546 SNPs. F_ST outlier methods at the SNP level revealed a weak signal of selection, with three marginally significant outliers in the northern populations. At the gene level, considering haplotypes as alleles, two additional marginally significant outliers were detected, one on each slope. To account for the uncertainty of haplotype inference, we averaged the Bayes factors over many possible phase reconstructions. Epistatic selection offers a realistic multilocus model of selection in natural populations. Here, we used a test suggested by Ohta based on the decomposition of the variance of linkage disequilibrium. Overall populations, 0.23% of the SNP pairs (haplotypes) showed evidence of epistatic selection, with nearly 80% of them being within genes. One of the between gene epistatic selection signals arose between an F_ST outlier and a nonsynonymous mutation in a drought response gene. Additionally, we identified haplotypes containing selectively advantageous allele combinations which were unique to high or low elevations and northern or southern populations. Several haplotypes contained nonsynonymous mutations situated in genes with known functional importance for adaptation to climatic factors.     

Genetic structure and distribution of Gelidium elegans (Gelidiales,Rhodophyta) in Korea based on mitochondrial cox1 sequence data

Gelidium elegans Kützing is commonly found in Korea, China, and Japan, and is the economically most important agarophyte in the northwest Pacific. To assess the genetic structure of the Korean species, we analyzed 1200 base pairs of the mitochondrial cox1 gene from 272 individuals collected from 36 locations. A total of 34 haplotypes were found, most of which were unique, including 27 (79%) ‘private’ haplotypes. The nucleotide and haplotype diversities of cox1 within G. elegans were 0.711 ± 0.028 (H) and 0.00736 ± 0.00038 (π), respectively. The distribution of cox1 haplotypes, pairwise F_ST values, results of neutrality tests, AMOVA, and mismatch distribution revealed the existence of a deep genetic break between central Pacific Japan and all the other locations, corresponding to the surface seawater current patterns as well as the genetic signature of potential demographic expansion.     

Environmental confirmation of multiple ice age refugia for Pacific cod, Gadus macrocephalus

The concept of species surviving through quaternary climatic extremes by retreating to glacial refugia, and then evolving genetically during re-population movements of the following interglacial, has been in the literature for over 40 years. Recently, advances in genetic analysis have enabled this concept to be validated and theories regarding population expansions and contractions to be built. For the major Northern Hemisphere species of cod, Gadus morhua (Atlantic cod) and Gadus macrocephalus (Pacific cod), genetic analysis has suggested retreat to separate refugia on both sides of their respective ocean basins during the last glacial period. Ecosystem niche modelling has previously confirmed that environmental conditions during the last glacial were compatible with the existence of these separate refugia for Atlantic cod. Here it is shown that such modelling also confirms a reduced core glacial distribution for G. macrocephalus, but probable refugia on either side of the Pacific. Existing mitochondrial DNA analyses suggest two separate glacial populations in the northwest Pacific, which modelling confirms, with predicted separate marine refugia in the land-locked Sea of Japan basin and the Sea of Okhotsk. Existing mitochondrial DNA for the northeast Pacific populations is less conclusive regarding whether there were one or two separate refugia off this coast, and their location. Using environmental niche models this study shows the glacial NE Pacific environment could support two marine refugia, one centred in the Aleutians/Gulf of Alaska and the other off British Columbia. The intervening Cordilleran Ice Sheet, and the glacially sub-aerial and ice-free Queen Charlotte Islands shelf, is hypothesised to have constrained exchange between glacial stocks either side of the Islands. The postulated southern marine refugium is off-shore from an established terrestrial refugium, suggesting greater dependence of species and ecosystems during environmental change. An earth system approach to evolutionary change could enhance understanding of past and future ecosystems.     

Genomic divergence in sympatry indicates strong reproductive barriers and cryptic species within Eucalyptus salubris

Genetic studies are increasingly detecting cryptic taxa that likely represent a significant component of global biodiversity. However, cryptic taxa are often criticized because they are typically detected serendipitously and may not receive the follow‐up study required to verify their geographic or evolutionary limits. Here, we follow‐up a study of Eucalyptus salubris that unexpectedly detected two divergent lineages but was not sampled sufficiently to make clear interpretations. We undertook comprehensive sampling for an independent genomic analysis (3,605 SNPs) to investigate whether the two purported lineages remain discrete genetic entities or if they intergrade throughout the species’ range. We also assessed morphological and ecological traits, and sequenced chloroplast DNA. SNP results showed strong genome‐wide divergence (F _ST = 0.252) between two discrete lineages: one dominated the north and one the southern regions of the species’ range. Within lineages, gene flow was high, with low differentiation (mean F _ST = 0.056) spanning hundreds of kilometers. In the central region, the lineages were interspersed but maintained their genomic distinctiveness: an indirect demonstration of reproductive isolation. Populations of the southern lineage exhibited significantly lower specific leaf area and occurred on soils with lower phosphorus relative to the northern lineage. Finally, two major chloroplast haplotypes were associated with each lineage but were shared between lineages in the central distribution. Together, these results suggest that these lineages have non‐contemporary origins and that ecotypic adaptive processes strengthened their divergence more recently. We conclude that these lineages warrant taxonomic recognition as separate species and provide fascinating insight into eucalypt speciation.     

Population differences in Chinook salmon (Oncorhynchus tshawytscha) DNA methylation: Genetic drift and environmental factors

Local adaptation and phenotypic differences among populations have been reported in many species, though most studies focus on either neutral or adaptive genetic differentiation. With the discovery of DNA methylation, questions have arisen about its contribution to individual variation in and among natural populations. Previous studies have identified differences in methylation among populations of organisms, although most to date have been in plants and model animal species. Here we obtained eyed eggs from eight populations of Chinook salmon (Oncorhynchus tshawytscha) and assayed DNA methylation at 23 genes involved in development, immune function, stress response, and metabolism using a gene‐targeted PCR‐based assay for next‐generation sequencing. Evidence for population differences in methylation was found at eight out of 23 gene loci after controlling for developmental timing in each individual. However, we found no correlation between freshwater environmental parameters and methylation variation among populations at those eight genes. A weak correlation was identified between pairwise DNA methylation dissimilarity among populations and pairwise F _ST based on 15 microsatellite loci, indicating weak effects of genetic drift or geographic distance on methylation. The weak correlation was primarily driven by two genes, GTIIBS and Nkef. However, single‐gene Mantel tests comparing methylation and pairwise F _ST were not significant after Bonferroni correction. Thus, population differences in DNA methylation are more likely related to unmeasured oceanic environmental conditions, local adaptation, and/or genetic drift. DNA methylation is an additional mechanism that contributes to among population variation, with potential influences on organism phenotype, adaptive potential, and population resilience.     

Evidence for transoceanic migrations by loggerhead sea turtles in the southern Pacific Ocean

Post-hatchling loggerhead turtles (Caretta caretta) in the northern Pacific and northern Atlantic Oceans undertake transoceanic developmental migrations. Similar migratory behaviour is hypothesized in the South Pacific Ocean as post-hatchling loggerhead turtles are observed in Peruvian fisheries, yet no loggerhead rookeries occur along the coast of South America. This hypothesis was supported by analyses of the size-class distribution of 123 post-hatchling turtles in the South Pacific and genetic analysis of mtDNA haplotypes of 103 nesting females in the southwest Pacific, 19 post-hatchlings stranded on the southeastern Australian beaches and 22 post-hatchlings caught by Peruvian longline fisheries. Only two haplotypes (CCP1 93% and CCP5 7%) were observed across all samples, and there were no significant differences in haplotype frequencies between the southwest Pacific rookeries and the post-hatchlings. By contrast, the predominant CCP1 haplotype is rarely observed in North Pacific rookeries and haplotype frequencies were strongly differentiated between the two regions (F_st=0.82; p=<0.00001). These results suggest that post-hatchling loggerhead turtles emerging from the southwest Pacific rookeries are undertaking transoceanic migrations to the southeastern Pacific Ocean, thus emphasizing the need for a broader focus on juvenile mortality throughout the South Pacific to develop effective conservation strategies.     

Multi-Gene Analysis Reveals a Lack of Genetic Divergence between Calanus agulhensis and C. sinicus (Copepoda; Calanoida)

The discrimination and taxonomic identification of marine species continues to pose a challenge despite the growing number of diagnostic metrics and approaches. This study examined the genetic relationship between two sibling species of the genus Calanus (Crustacea; Copepoda; Calanidae), C. agulhensis and C. sinicus, using a multi-gene analysis. DNA sequences were determined for portions of the mitochondrial cytochrome c oxidase I (mtCOI); nuclear citrate synthase (CS), and large subunit (28S) rRNA genes for specimens collected from the Sea of Japan and North East (NE) Pacific Ocean for C. sinicus and from the Benguela Current and Agulhas Bank, off South Africa, for C. agulhensis. For mtCOI, C. sinicus and C. agulhensis showed similar levels of haplotype diversity (H_d = 0.695 and 0.660, respectively) and nucleotide diversity (π = 0.003 and 0.002, respectively). Pairwise F_ST distances for mtCOI were significant only between C. agulhensis collected from the Agulhas and two C. sinicus populations: the Sea of Japan (F_ST = 0.152, p<0.01) and NE Pacific (F_ST = 0.228, p<0.005). Between the species, F_ST distances were low for both mtCOI (F_ST = 0.083, p = 0.003) and CS (F_ST = 0.050, p = 0.021). Large subunit (28S) rRNA showed no variation between the species. Our results provide evidence of the lack of genetic distinction of C. sinicus and C. agulhensis, raise questions of whether C. agulhensis warrants status as a distinct species, and indicate the clear need for more intensive and extensive ecological and genetic analysis.     

Genetic variation and local differences in Pacific cod Gadus macrocephalus around Japan

The population structure of the Pacific cod Gadus macrocephalus was examined using 15 microsatellite loci and mitochondrial DNA (ND2 region). In total, 274 individuals were sampled from 16 locations around Japan to estimate the level of genetic differentiation and effective population size (N_e). Pairwise F_ST, analysis of molecular variance and Bayesian clustering analysis suggested the presence of two genetically distinct groups in waters around Japan, with a higher N_e value in the eastern group than in the western group. A possible factor that restricts gene flow between groups may be related to the water temperature differences in the south‐western part of the Sea of Japan, where the Tsushima Warm Current flows around the area inhabited by the western group, which may limit migration between the west and east.     

Genetic differentiation and connectivity of morphological types of the broadcast‐spawning coral Galaxea fascicularis in the Nansei Islands,Japan

Population connectivity resulting from larval dispersal is essential for the maintenance or recovery of populations in marine ecosystems, including coral reefs. Studies of species diversity and genetic connectivity within species are essential for the conservation of corals and coral reef ecosystems. We analyzed mitochondrial DNA sequence types and microsatellite genotypes of the broadcast‐spawning coral, Galaxea fascicularis, from four regions in the subtropical Nansei Islands in the northwestern Pacific Ocean. Two types (soft and hard types) of nematocyst morphology are known in G. fascicularis and are significantly correlated with the length of a mitochondrial DNA noncoding sequence (soft type: mt‐L; hard type: mt‐S type). Using microsatellites, significant genetic differentiation was detected between the mitochondrial DNA sequence types in all regions. We also found a third genetic cluster (mt‐L+), and this unexpected type may be a cryptic species of Galaxea. High clonal diversity was detected in both mt‐L and mt‐S types. Significant genetic differentiation, which was found among regions within a given type (F _ST = 0.009–0.024, all Ps ≤ 0.005 in mt‐L; 0.009–0.032, all Ps ≤ 0.01 in mt‐S), may result from the shorter larval development than in other broadcast‐spawning corals, such as the genus Acropora. Nevertheless, intraspecific genetic diversity and connectivity have been maintained, and with both sexual and asexual reproduction, this species appears to have a potential for the recovery of populations after disturbance.