Limited genomic signatures of population collapse in the critically endangered black abalone (Haliotis cracherodii).

The black abalone, Haliotis cracherodii, is a large, long-lived marine mollusc that inhabits rocky intertidal habitats along the coast of California and Mexico. In 1985, populations were impacted by a bacterial disease known as withering syndrome (WS) that wiped out >90% of individuals, leading to the closure of all U.S. black abalone fisheries since 1993. Current conservation strategies include restoring diminished populations by translocating healthy individuals. However, population collapse on this scale may have dramatically lowered genetic diversity and strengthened geographic differentiation, making translocation-based recovery contentious. Additionally, the current prevalence of WS remains unknown. To address these uncertainties, we sequenced and analysed the genomes of 133 black abalone individuals from across their present range. We observed no spatial genetic structure among black abalone, with the exception of a single chromosomal inversion that increases in frequency with latitude. Outside the inversion, genetic differentiation between sites is minimal and does not scale with either geographic distance or environmental dissimilarity. Genetic diversity appears uniformly high across the range. Demographic inference does indicate a severe population bottleneck beginning just 15 generations in the past, but this decline is short lived, with present-day size far exceeding the pre-bottleneck status quo. Finally, we find the bacterial agent of WS is equally present across the sampled range, but only in 10% of individuals. The lack of population genetic structure, uniform diversity and prevalence of WS bacteria indicates that translocation could be a valid and low-risk means of population restoration for black abalone species' recovery.

A sex chromosome polymorphism maintains divergent plumage phenotypes between extensively hybridizing yellowhammers (Emberiza citrinella) and pine buntings (E. leucocephalos).

Under allopatric speciation, populations of a species become isolated by a geographic barrier and develop reproductive isolation through genetic differentiation. When populations meet in secondary contact, the strength of evolved reproductive barriers determines the extent of hybridization and whether the populations will continue to diverge or merge together. The yellowhammer (Emberiza citrinella) and pine bunting (E. leucocephalos) are avian sister species that diverged in allopatry on either side of Eurasia during the Pleistocene glaciations. Though they differ greatly in plumage and form distinct genetic clusters in allopatry, these taxa show negligible mitochondrial DNA differentiation and hybridize extensively where they overlap in central Siberia, lending uncertainty to the state of reproductive isolation in the system. To assess the strength of reproductive barriers between taxa, we examined genomic differentiation across the system. We found that extensive admixture has occurred in sympatry, indicating that reproductive barriers between taxa are weak. We also identified a putative Z chromosome inversion region that underlies plumage variation in the system, with the 'pine bunting' haplotype showing dominance over the 'yellowhammer' haplotype. Our results suggest that yellowhammers and pine buntings are currently at a crossroads and that evolutionary forces may push this system towards either continued differentiation or population merging. However, even if these taxa merge, recombination suppression between putative chromosome Z inversion haplotypes may maintain divergent plumage phenotypes within the system. In this way, our findings highlight the important role hybridization plays in increasing the genetic and phenotypic variation as well as the evolvability of a system.

The interplay of local adaptation and gene flow may lead to the formation of supergenes.

Supergenes are genetic architectures resulting in the segregation of alternative combinations of alleles underlying complex phenotypes. The co-segregation of alleles at linked loci is often facilitated by polymorphic chromosomal rearrangements suppressing recombination locally. Supergenes are involved in many complex polymorphisms, including sexual, colour or behavioural polymorphisms in numerous plants, fungi, mammals, fish, and insects. Despite a long history of empirical and theoretical research, the formation of supergenes remains poorly understood. Here, using a two-island population genetic model, we explore how gene flow and the evolution of overdominant chromosomal inversions may jointly lead to the formation of supergenes. We show that the evolution of inversions in differentiated populations, both under disruptive selection, leads to an increase in frequency of poorly adapted, immigrant haplotypes. Indeed, rare allelic combinations, such as immigrant haplotypes, are more frequently reshuffled by recombination than common allelic combinations, and therefore benefit from the recombination suppression generated by inversions. When an inversion capturing a locally adapted haplotype spreads but is associated with a fitness cost hampering its fixation (e.g. a recessive mutation load), the maintenance of a non-inverted haplotype in the population is enhanced; under certain conditions, the immigrant haplotype persists alongside the inverted local haplotype, while the standard local haplotype disappears. This establishes a stable, local polymorphism with two non-recombining haplotypes encoding alternative adaptive strategies, that is, a supergene. These results bring new light to the importance of local adaptation, overdominance, and gene flow in the formation of supergenes and inversion polymorphisms in general.

A supergene affects androgen concentrations during early development in a bird with alternative reproductive morphs.

Species with alternative reproductive tactics typically show pronounced phenotypic variation between and within sexes. In some species, this variation culminates in discrete reproductive morphs that are genetically determined, facilitating studies on how genetic variation translates into phenotypic variation. In ruffs (Calidris pugnax), an autosomal inversion polymorphism underlies three reproductive morphs (Independents, Satellites and Faeders), which differ in circulating steroid concentrations in adults. Yet, it remains unknown whether morph differences in steroid concentrations already arise before adulthood. We examined variation in circulating testosterone, androstenedione and progesterone concentrations between morphs and sexes in ruff chicks and juveniles and compared the differences to those in adults. Since measured hormone concentrations only provide momentary states and show high within- and between-individual variation, we took repeated measurements to compare means, variances and skewness between groups. We found clear differences between morphs but not the sexes in early life. Between morphs, androgen concentrations in young ruffs differed in variance and skewness, but not in their means. For testosterone, Independents had a higher variance than Satellites/Faeders, whereas for androstenedione, we observed the opposite pattern. For progesterone, we did not detect clear differences between groups. Skewness values mirrored differences in morph variances. Compared to adults, premature ruffs had lower androgen concentrations. In both life stages, we detected morph-specific associations between androgen concentrations: androstenedione concentrations increased with testosterone concentrations more in Satellites/Faeders than in Independents. These observed morph differences during early life are consistent with a supergene-mediated regulation of androgen variation that underlies the diversification of adult behavioural phenotypes.

Targeted genomic translocations and inversions generated using a paired prime editing strategy.

A variety of cancers have been found to have chromosomal rearrangements, and the genomic abnormalities often induced expression of fusion oncogenes. To date, a pair of engineered nucleases including ZFNs, TALENs, and CRISPR-Cas9 nucleases have been used to generate chromosomal rearrangement in living cells and organisms for disease modeling. However, these methods induce unwanted indel mutations at the DNA break junctions, resulting in incomplete disease modeling. Here, we developed prime editor nuclease-mediated translocation and inversion (PETI), a method for programmable chromosomal translocation and inversion using prime editor 2 nuclease (PE2 nuclease) and paired pegRNA. Using PETI method, we successfully introduced DNA recombination in episomal fluorescence reporters as well as precise chromosomal translocations in human cells. We applied PETI to create cancer-associated translocations and inversions such as NPM1-ALK and EML4-ALK in human cells. Our findings show that PETI generated chromosomal translocation and inversion in a programmable manner with efficiencies comparable of Cas9. PETI methods, we believe, could be used to create disease models or for gene therapy.

Weak antagonistic fitness effects can maintain an inversion polymorphism.

The study of chromosomal inversion polymorphisms has received much recent attention, particularly in cases where inversions have drastic effects on phenotypes and fitness (e.g. lethality of homozygotes). Less attention has been paid to the question of the maintenance of inversion polymorphisms that show only weak effects. Here, we study the maintenance of such an inversion polymorphism that links 250 genes on chromosome Tgu11 in the zebra finch (Taeniopygia guttata). Based on data from over 6000 captive birds, we estimated the effects of this inversion on a wide range of fitness-related traits. We found that, compared with the ancestral allele A, the inverted allele D had small additive beneficial effects on male siring success and on female fecundity. These fitness-enhancing effects may explain the initial spread of the derived D allele (allele frequency 53%). However, individuals that were homozygous for D had a slightly lower survival rate, which may explain why the D allele has not spread to fixation. We used individual-based simulations to examine how an inversion polymorphism with such antagonistic fitness effects behaves over time. Our results indicate that polymorphisms become stabilized at an intermediate allele frequency if the inversion links an additively beneficial allele of small effect size to a recessive weakly deleterious mutation, overall resulting in weak net heterosis. Importantly, this conclusion remains valid over a wide range of selection coefficients against the homozygous DD (up to lethality), suggesting that the conditions needed to maintain the polymorphism may frequently be met. However, the simulations also suggest that in our zebra finch populations, the estimated recessive deleterious effect of the D allele (on survival in captivity) is not quite large enough to prevent fixation of the D allele in the long run. Estimates of fitness effects from free-living populations are needed to validate these results.

Genomic architecture underlying morphological and physiological adaptation to high elevation in a songbird.

Organisms often acquire physiological and morphological modifications to conquer ecological challenges when colonizing new environments which lead to their adaptive evolution. However, deciphering the genomic mechanism of ecological adaptation is difficult because ecological environments are often too complex for straightforward interpretation. Thus, we examined the adaptation of a widespread songbird-the rufous-capped babbler (Cyanoderma ruficeps)-to a relatively simple system: distinct environments across elevational gradients on the mountainous island of Taiwan. We focused on the genomic sequences of 43 birds from five populations to show that the Taiwan group split from its sister group in mainland China around 1-2 million years ago (Ma) and colonized the montane habitats of Taiwan at least twice around 0.03-0.22 Ma. The montane and lowland Taiwan populations diverged with gene flow between them, suggesting strong selection associated with different elevations. We found that the montane babblers had smaller beaks than the lowland ones, consistent with Allen's rule, and identified candidate genes-COL9A1 and SOX11-underlying the beak size changes. We also found that altitudinally divergent mutations were mostly located in noncoding regions and tended to accumulate in chromosomal inversions and autosomes. The altitudinally divergent mutations might regulate genes related to haematopoietic, metabolic, immune, auditory and vision functions, as well as cerebrum morphology and plumage development. The results reveal the genomic bases of morphological and physiological adaptation in this species to the low temperature, hypoxia and high UV light environment at high elevation. These findings improve our understanding of how ecological adaptation drives population divergence from the perspective of genomic architecture.

Examination of Large Chromosomal Inversions in the Genome of Erwinia amylovora Strains Reveals Worldwide Distribution and North America-Specific Types.

Erwinia amylovora is a relatively homogeneous species with low genetic diversity at the nucleotide level. However, phenotypic differences and genomic structural variations among E. amylovora strains have been documented. In this study, we identified 10 large chromosomal inversion (LCI) types in the Spiraeoideae-infecting (SI) E. amylovora strains by combining whole genome sequencing and PCR-based molecular markers. It was found that LCIs were mainly caused by homologous recombination events among seven rRNA operons (rrns) in SI E. amylovora strains. Although ribotyping results identified inter- and intra-variations in the internal transcribed spacer (ITS1 and ITS2) regions among rrns, LCIs tend to occur between rrns transcribed in the opposite directions and with the same tRNA content (tRNA-Glu or tRNA-Ile/Ala) in ITS1. Based on the LCI types, physical/estimated replichore imbalance (PRI/ERI) was examined and calculated. Among the 117 SI strains evaluated, the LCI types of Ea1189, CFBP1430, and Ea273 were the most common, with ERI values at 1.31, 7.87, and 4.47°, respectively. These three LCI types had worldwide distribution, whereas the remaining seven LCI types were restricted to North America (or certain regions of the United States). Our results indicated ongoing chromosomal recombination events in the SI E. amylovora population and showed that LCI events are mostly symmetrical, keeping the ERI less than 15°. These findings provide initial evidence about the prevalence of certain LCI types in E. amylovora strains, how LCI occurs, and its potential evolutionary advantage and history, which might help track the movement of the pathogen.

A supergene-linked estrogen receptor drives alternative phenotypes in a polymorphic songbird.

Behavioral evolution relies on genetic changes, yet few behaviors can be traced to specific genetic sequences in vertebrates. Here we provide experimental evidence showing that differentiation of a single gene has contributed to the evolution of divergent behavioral phenotypes in the white-throated sparrow, a common backyard songbird. In this species, a series of chromosomal inversions has formed a supergene that segregates with an aggressive phenotype. The supergene has captured ESR1, the gene that encodes estrogen receptor α (ERα); as a result, this gene is accumulating changes that now distinguish the supergene allele from the standard allele. Our results show that in birds of the more aggressive phenotype, ERα knockdown caused a phenotypic change to that of the less aggressive phenotype. We next showed that in a free-living population, aggression is predicted by allelic imbalance favoring the supergene allele. Finally, we identified cis-regulatory features, both genetic and epigenetic, that explain the allelic imbalance. This work provides a rare illustration of how genotypic divergence has led to behavioral phenotypic divergence in a vertebrate.

A large chromosomal inversion affects antimicrobial sensitivity of Escherichia coli to sodium deoxycholate.

Resistance to antimicrobials is normally caused by mutations in the drug targets or genes involved in antimicrobial activation or expulsion. Here we show that an Escherichia coli strain, named DOC14, selected for increased resistance to the bile salt sodium deoxycholate, has no mutations in any ORF, but instead has a 2.1 Mb chromosomal inversion. The breakpoints of the inversion are two inverted copies of an IS5 element. Besides lowering deoxycholate susceptibility, the IS5-mediated chromosomal inversion in the DOC14 mutant was found to increase bacterial survival upon exposure to ampicillin and vancomycin, and sensitize the cell to ciprofloxacin and meropenem, but does not affect bacterial growth or cell morphology in a rich medium in the absence of antibacterial molecules. Overall, our findings support the notion that a large chromosomal inversion can benefit bacterial cells under certain conditions, contributing to genetic variability available for selection during evolution. The DOC14 mutant paired with its isogenic parental strain form a useful model as bacterial ancestors in evolution experiments to study how a large chromosomal inversion influences the evolutionary trajectory in response to various environmental stressors.

Temperature-associated selection linked to putative chromosomal inversions in king scallop (Pecten maximus).

The genomic landscape of divergence-the distribution of differences among populations or species across the genome-is increasingly characterized to understand the role that microevolutionary forces such as natural selection and recombination play in causing and maintaining genetic divergence. This line of inquiry has also revealed chromosome structure variation to be an important factor shaping the landscape of adaptive genetic variation. Owing to a high prevalence of chromosome structure variation and the strong pressure for local adaptation necessitated by their sessile nature, bivalve molluscs are an ideal taxon for exploring the relationship between chromosome structure variation and local adaptation. Here, we report a population genomic survey of king scallop (Pecten maximus) across its natural range in the northeastern Atlantic Ocean, using a recent chromosome-level genome assembly. We report the presence of at least three large (12-22 Mb), putative chromosomal inversions associated with sea surface temperature and whose frequencies are in contrast to neutral population structure. These results highlight a potentially large role for recombination-suppressing chromosomal inversions in local adaptation and suggest a hypothesis to explain the maintenance of differences in reproductive timing found at relatively small spatial scales across king scallop populations.

Neo-sex chromosome evolution and phenotypic differentiation across an elevational gradient in horned larks (Eremophila alpestris).

Genetic structure and phenotypic variation among populations are affected by both geographic distance and environmental variation across species' distributions. Understanding the relative contributions of isolation by distance (IBD) and isolation by environment (IBE) is important for elucidating population dynamics across habitats and ecological gradients. In this study, we compared phenotypic and genetic variation among Horned Lark (Eremophila alpestris) populations from 10 sites encompassing an elevational gradient from low-elevation desert scrub in Death Valley (285 a.s.l.) to high-elevation meadows in the White Mountains of the Sierra Nevada of California (greater than 3000 m a.s.l.). Using a ddRAD data set of 28,474 SNPs aligned to a high-quality reference genome, we compared genetic structure with elevational, environmental, and spatial distance to quantify how different aspects of the landscape drive genomic and phenotypic differentiation in Horned Larks. We found larger-bodied birds were associated with sites that had less seasonality and higher annual precipitation, and longer spurs occurred in soils with more clay and silt content, less sand, and finer fragments. Larks have large neo-sex chromosomes, and we found that associations with elevation and environmental variation were much stronger among neo-sex chromosomes compared to autosomes. Furthermore, we found that putative chromosomal translocations, fusions, and inversions were associated with elevation and may underlie local adaptation across an elevational gradient in Horned Larks. Our results suggest that genetic variation in Horned Larks is affected more by IBD than IBE, but specific phenotypes and genomic regions-particually on neo-sex chromosomes-bear stronger associations with the environment.

Do chromosomal inversion carriers really need preimplantation genetic testing?

PURPOSE: This study aimed to evaluate the rates of euploidy, aneuploidy, and mosaicism in preimplantation genetic testing for structural rearrangements (PGT-SR) cycles from chromosomal inversion carriers. In addition, this work also focused on assessing the impact of some contributors on the incidence of parental originating aneuploidy and mosaicism. METHODS: This retrospective review enrolled chromosomal inversion carrier couples of whom the females were under 38 years old undergoing PGT-SR at a single academic reproductive center. Subgroups were divided according to the gender of carriers, the inversion type, and the semen parameters of male carriers (male factor infertility (MF) or non-MF). Patient demographics, cycle characteristics, and PGT-SR outcomes were compared among subgroups. RESULTS: A total of 71 PGT-SR cycles from 57 inversion carrier couples were included for analysis. Among the 283 blastocysts, 48.4% were identified as euploidy, 27.9% as aneuploidy, and the remaining 23.7% as mosaicism. Only 32.9% of aneuploid embryos and 1.5% of mosaic embryos involved the parental inversion chromosomes. Notably, the female inversion carriers seemed to produce more parental originating aneuploid embryos than male inversion carriers (45.5% vs 23.9%, p = 0.044). CONCLUSIONS: The type of inversion and sperm parameters of male chromosomal inversion carriers did not affect the ploidy status of embryos. The incidence of parental originating aneuploidy in inversion carrier couples is lower than expected. For male chromosomal inversion carriers with normal sperm condition whose female partners are under 38 years old, natural conception combined with prenatal diagnosis could be provided as an option during fertility counseling.

New Wolbachia pipientis Genotype Increasing Heat Stress Resistance of Drosophila melanogaster Host Is Characterized by a Large Chromosomal Inversion.

The maternally transmitted endocellular bacteria Wolbachia is a well-known symbiont of insects, demonstrating both negative and positive effects on host fitness. The previously found Wolbachia strain wMelPlus is characterized by a positive effect on the stress-resistance of its host Drosophila melanogaster, under heat stress conditions. This investigation is dedicated to studying the genomic underpinnings of such an effect. We sequenced two closely related Wolbachia strains, wMelPlus and wMelCS112, assembled their complete genomes, and performed comparative genomic analysis engaging available Wolbachia genomes from the wMel and wMelCS groups. Despite the two strains under study sharing very close gene-composition, we discovered a large (>1/6 of total genome) chromosomal inversion in wMelPlus, spanning through the region that includes the area of the inversion earlier found in the wMel group of Wolbachia genotypes. A number of genes in unique inversion blocks of wMelPlus were identified that might be involved in the induction of a stress-resistant phenotype in the host. We hypothesize that such an inversion could rearrange established genetic regulatory-networks, causing the observed effects of such a complex fly phenotype as a modulation of heat stress resistance. Based on our findings, we propose that wMelPlus be distinguished as a separate genotype of the wMelCS group, named wMelCS3.

A hapless mathematical contribution to biology : Chromosome inversions in Drosophila, 1937-1941.

This is the story, told in the light of a new analysis of historical data, of a mathematical biology problem that was explored in the 1930s in Thomas Morgan's laboratory at the California Institute of Technology. It is one of the early developments of evolutionary genetics and quantitative phylogeny, and deals with the identification and counting of chromosomal inversions in Drosophila species from comparisons of genetic maps. A re-analysis of the data produced in the 1930s using current mathematics and computational technologies reveals how a team of biologists, with the help of a renowned mathematician and against their first intuition, came to an erroneous conclusion regarding the presence of phylogenetic signals in gene arrangements. This example illustrates two different aspects of a same piece: (1) the appearance of a mathematical in biology problem solved with the development of a combinatorial algorithm, which was unusual at the time, and (2) the role of errors in scientific activity. Also underlying is the possible influence of computational complexity in understanding the directions of research in biology.

A 1.7-Mb chromosomal inversion downstream of a PpOFP1 gene is responsible for flat fruit shape in peach.

Flat peaches have become popular worldwide due to their novelty and convenience. The peach flat fruit trait is genetically controlled by a single gene at the S locus, but its genetic basis remains unclear. Here, we report a 1.7-Mb chromosomal inversion downstream of a candidate gene encoding OVATE Family Protein, designated PpOFP1, as the causal mutation for the peach flat fruit trait. Genotyping of 727 peach cultivars revealed an occurrence of this large inversion in flat peaches, but absent in round peaches. Ectopic overexpression of PpOFP1 resulted in oval-shaped leaves and shortened siliques in Arabidopsis, suggesting its role in repressing cell elongation. Transcriptional activation of PpOFP1 by the chromosomal inversion may repress vertical elongation in flat-shaped fruits at early stages of development, resulting in the flat fruit shape. Moreover, PpOFP1 can interact with fruit elongation activator PpTRM17, suggesting a regulatory network controlling fruit shape in peach. Additionally, screening of peach wild relatives revealed an exclusive presence of the chromosomal inversion in P. ferganensis, supporting that this species is the ancestor of the domesticated peach. This study provides new insights into mechanisms underlying fruit shape evolution and molecular tools for genetic improvement of fruit shape trait in peach breeding programmes.

Inversion breakpoints and the evolution of supergenes.

The coexistence of discrete morphs that differ in multiple traits is common within natural populations of many taxa. Such morphs are often associated with chromosomal inversions, presumably because the recombination suppressing effects of inversions help maintain alternate adaptive combinations of alleles across the multiple loci affecting these traits. However, inversions can also harbour selected mutations at their breakpoints, leading to their rise in frequency in addition to (or independent from) their role in recombination suppression. In this review, we first describe the different ways that breakpoints can create mutations. We then critically examine the evidence for the breakpoint-mutation and recombination suppression hypotheses for explaining the existence of discrete morphs associated with chromosomal inversions. We find that the evidence that inversions are favoured due to recombination suppression is often indirect. The evidence that breakpoints harbour mutations that are adaptive is also largely indirect, with the characterization of inversion breakpoints at the sequence level being incomplete in most systems. Direct tests of the role of suppressed recombination and breakpoint mutations in inversion evolution are thus needed. Finally, we emphasize how the two hypotheses of recombination suppression and breakpoint mutation can act in conjunction, with implications for understanding the emergence of supergenes and their evolutionary dynamics. We conclude by discussing how breakpoint characterization could improve our understanding of complex, discrete phenotypic forms in nature.

Standing variation rather than recent adaptive introgression probably underlies differentiation of the texanus subspecies of Helianthus annuus.

The origins of geographic races in wide-ranging species are poorly understood. In Texas, the texanus subspecies of Helianthus annuus has long been thought to have acquired its defining phenotypic traits via introgression from a local congener, H. debilis, but previous tests of this hypothesis were inconclusive. Here, we explore the origins of H. a. texanus using whole genome sequencing data from across the entire range of H. annuus and possible donor species, as well as phenotypic data from a common garden study. We found that although it is morphologically convergent with H. debilis, H. a. texanus has conflicting signals of introgression. Genome wide tests (Patterson's D and TreeMix) only found evidence of introgression from H. argophyllus (sister species to H. annuus and also sympatric), but not H. debilis, with the exception of one individual of 109 analysed. We further scanned the genome for localized signals of introgression using PCAdmix and found minimal but nonzero introgression from H. debilis and significant introgression from H. argophyllus in some populations. Given the paucity of introgression from H. debilis, we argue that the morphological convergence observed in Texas is probably from standing genetic variation. We also found that genomic differentiation in H. a. texanus is mostly driven by large segregating inversions, several of which have signatures of natural selection based on haplotype frequencies.

Systems genetic analysis of inversion polymorphisms in the malaria mosquito Anopheles gambiae.

Inversion polymorphisms in the African malaria vector Anopheles gambiae segregate along climatic gradients of aridity. Despite indirect evidence of their adaptive significance, little is known of the phenotypic targets of selection or the underlying genetic mechanisms. Here we adopt a systems genetics approach to explore the interaction of two inversions on opposite arms of chromosome 2 with gender, climatic conditions, and one another. We measure organismal traits and transcriptional profiles in 8-d-old adults of both sexes and four alternative homokaryotypic classes reared under two alternative climatic regimes. We show that karyotype strongly influences both organismal traits and transcriptional profiles but that the strength and direction of the effects depend upon complex interactions with gender and environmental conditions and between inversions on independent arms. Our data support the suppressed recombination model for the role of inversions in local adaptation, and-supported by transcriptional and physiological measurements following perturbation with the drug rapamycin-suggest that one mechanism underlying their adaptive role may be the maintenance of energy homeostasis.

Life-stage-dependent supergene haplotype frequencies and metapopulation neutral genetic patterns of Atlantic cod, Gadus morhua, from Canada's Northern cod stock region and adjacent areas.

Among highly migratory fish species, nursery areas occupied by juveniles often differ from adult habitats. To better understand the spatial dynamics of Canada's Northern cod stock, juveniles caught off the east coast of Newfoundland and Labrador were compared to adults from the same region as well as individuals from other areas in Atlantic Canada using double-digest restriction site-associated DNA sequencing-derived single nucleotide polymorphisms. A reduced proportion of homozygotes with a chromosomal inversion located in linkage group 1 (LG1) was detected between juvenile and adult samples in the Northern cod stock region, potentially indicating age-dependent habitat use or ontogenetic selection for attributes associated with the many genes located in LG1. No selectively neutral genetic differences were found between samples from the Northern cod stock; nevertheless, significant differences were found between some of these samples and cod collected from St. Pierre Bank, Bay of Fundy, Browns Bank and the southern Scotian Shelf. Clustering analysis of variants at neutral loci provided evidence for three major genetic units: (a) the Newfoundland Atlantic Coast, (b) eastern and southern Gulf of St. Lawrence and Burgeo Bank and (c) the Bay of Fundy, Browns Bank and southern Scotian Shelf. Both adaptive and neutral population structure within the Northern cod stock should be considered by managers to promote demographic rebuilding of the stock.