A chromosome-level assembly of the Atlantic herring genome-detection of a supergene and other signals of selection

Copy number variation and elevated genetic diversity at immune trait loci in Atlantic and Pacific herring

BACKGROUND

Genome-wide comparisons of populations are widely used to explore the patterns of nucleotide diversity and sequence divergence to provide knowledge on how natural selection and genetic drift affect the genome. In this study we have compared whole-genome sequencing data from Atlantic and Pacific herring, two sister species that diverged about 2 million years ago, to explore the pattern of genetic differentiation between the two species.

RESULTS

The genome comparison of the two species revealed high genome-wide differentiation but with islands of remarkably low genetic differentiation, as measured by an FST analysis. However, the low FST observed in these islands is not caused by low interspecies sequence divergence (dxy) but rather by exceptionally high estimated intraspecies nucleotide diversity (π). These regions of low differentiation and elevated nucleotide diversity, termed high-diversity regions in this study, are not enriched for repeats but are highly enriched for immune-related genes. This enrichment includes genes from both the adaptive immune system, such as immunoglobulin, T-cell receptor and major histocompatibility complex genes, as well as a substantial number of genes with a role in the innate immune system, e.g. novel immune-type receptor, tripartite motif and tumor necrosis factor receptor genes. Analysis of long-read based assemblies from two Atlantic herring individuals revealed extensive copy number variation in these genomic regions, indicating that the elevated intraspecies nucleotide diversities were partially due to the cross-mapping of short reads.

CONCLUSIONS

This study demonstrates that copy number variation is a characteristic feature of immune trait loci in herring. Another important implication is that these loci are blind spots in classical genome-wide screens for genetic differentiation using short-read data, not only in herring, likely also in other species harboring qualitatively similar variation at immune trait loci. These loci stood out in this study because of the relatively high genome-wide baseline for FST values between Atlantic and Pacific herring.

Chromosome-level genome assembly and population genomic analysis reveal evolution and local adaptation in common hairfin anchovy (Setipinna tenuifilis)

Understanding the genetic structure and the factors associated with adaptive diversity has significant implications for the effective management of wild populations under threat from overfishing and climate change. The common hairfin anchovy (Setipinna tenuifilis) is an economically and ecologically important pelagic fish species, spanning a broad latitudinal gradient along marginal seas of the Northwest Pacific. In this study, we constructed the first reference genome of S. tenuifilis using PacBio long reads and high-resolution chromosome conformation capture (Hi-C) technology. The assembled genome was 798.38 Mb with a contig N50 of 1.43 Mb and a scaffold N50 of 32.42 Mb, which were anchored onto 24 pseudochromosomes. A total of 22,019 genes were functionally annotated, which accounted for 95.27% of the predicted protein-coding genes. Chromosomal collinearity analysis revealed chromosome fusion or fission events in Clupeiformes species. Three genetic groups of S. tenuifilis were revealed along the Chinese coast using restriction site-associated DNA sequencing (RADseq). We investigated the influence of four bioclimatic variables as potential drivers of adaptive divergence in S. tenuifilis, suggesting that these environmental variables, especially sea surface temperature, may play important roles as drivers of spatially varying selection for S. tenuifilis. We also identified candidate functional genes underlying adaptive mechanisms and ecological tradeoffs using redundancy analysis (RDA) and BayeScan analysis. In summary, this study sheds light on the evolution and spatial patterns of genetic variation of S. tenuifilis, providing a valuable genomic resource for further biological and genetic studies on this species and other closely related Clupeiformes.

Diverse pathways to speciation revealed by marine snails

Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.

Adaptation to seasonal reproduction and environment-associated factors drive temporal and spatial differentiation in northwest Atlantic herring despite gene flow

Understanding how marine organisms adapt to local environments is crucial for predicting how populations will respond to global climate change. The genomic basis, environmental factors and evolutionary processes involved in local adaptation are however not well understood. Here we use Atlantic herring, an abundant, migratory and widely distributed marine fish with substantial genomic resources, as a model organism to evaluate local adaptation. We examined genomic variation and its correlation with environmental variables across a broad environmental gradient, for 15 spawning aggregations in Atlantic Canada and the United States. We then compared our results with available genomic data of northeast Atlantic populations. We confirmed that population structure lies in a fraction of the genome including likely adaptive genetic variants of functional importance. We discovered 10 highly differentiated genomic regions distributed across four chromosomes. Nine regions show strong association with seasonal reproduction. One region, corresponding to a known inversion on chromosome 12, underlies a latitudinal pattern discriminating populations north and south of a biogeographic transition zone on the Scotian Shelf. Genome-environment associations indicate that winter seawater temperature best correlates with the latitudinal pattern of this inversion. The variation at two so-called 'islands of divergence' related to seasonal reproduction appear to be private to the northwest Atlantic. Populations in the northwest and northeast Atlantic share variation at four of these divergent regions, simultaneously displaying significant diversity in haplotype composition at another four regions, which includes an undescribed structural variant approximately 7.7 Mb long on chromosome 8. Our results suggest that the timing and geographic location of spawning and early development may be under diverse selective pressures related to allelic fitness across environments. Our study highlights the role of genomic architecture, ancestral haplotypes and selection in maintaining adaptive divergence in species with large population sizes and presumably high gene flow.

How Fish Population Genomics Can Promote Sustainable Fisheries: A Road Map

Maintenance of genetic diversity in marine fishes targeted by commercial fishing is a grand challenge for the future. Most of these species are abundant and therefore important for marine ecosystems and food security. Here, we present a road map of how population genomics can promote sustainable fisheries. In these species, the development of reference genomes and whole genome sequencing is key, because genetic differentiation at neutral loci is usually low due to large population sizes and gene flow. First, baseline allele frequencies representing genetically differentiated populations within species must be established. These can then be used to accurately determine the composition of mixed samples, forming the basis for population demographic analysis to inform sustainably set fish quotas. SNP-chip analysis is a cost-effective method for determining baseline allele frequencies and for population identification in mixed samples. Finally, we describe how genetic marker analysis can transform stock identification and management.

Divergence and gene flow history at two large chromosomal inversions underlying ecotype differentiation in the long-snouted seahorse

Chromosomal inversions can play an important role in divergence and reproductive isolation by building and maintaining distinct allelic combinations between evolutionary lineages. Alternatively, they can take the form of balanced polymorphisms that segregate within populations until one arrangement becomes fixed. Many questions remain about how inversion polymorphisms arise, how they are maintained over the long term, and ultimately, whether and how they contribute to speciation. The long-snouted seahorse (Hippocampus guttulatus) is genetically subdivided into geographic lineages and marine-lagoon ecotypes, with shared structural variation underlying lineage and ecotype divergence. Here, we aim to characterize structural variants and to reconstruct their history and suspected role in ecotype formation. We generated a near chromosome-level genome assembly and described genome-wide patterns of diversity and divergence through the analysis of 112 whole-genome sequences from Atlantic, Mediterranean, and Black Sea populations. By also analysing linked-read sequencing data, we found evidence for two chromosomal inversions that were several megabases in length and showed contrasting allele frequency patterns between lineages and ecotypes across the species range. We reveal that these inversions represent ancient intraspecific polymorphisms, one likely being maintained by divergent selection and the other by pseudo-overdominance. A possible selective coupling between the two inversions was further supported by the absence of specific haplotype combinations and a putative functional interaction between the two inversions in reproduction. Lastly, we detected gene flux eroding divergence between inverted alleles at varying levels for the two inversions, with a likely impact on their dynamics and contribution to divergence and speciation.

Has recombination changed during the recent evolution of the guppy Y chromosome?

Genome sequencing and genetic mapping of molecular markers have demonstrated nearly complete Y-linkage across much of the guppy (Poecilia reticulata) XY chromosome pair. Predominant Y-linkage of factors controlling visible male-specific coloration traits also suggested that these polymorphisms are sexually antagonistic (SA). However, occasional exchanges with the X are detected, and recombination patterns also appear to differ between natural guppy populations, suggesting ongoing evolution of recombination suppression under selection created by partially sex-linked SA polymorphisms. We used molecular markers to directly estimate genetic maps in sires from 4 guppy populations. The maps are very similar, suggesting that their crossover patterns have not recently changed. Our maps are consistent with population genomic results showing that variants within the terminal 5 Mb of the 26.5 Mb sex chromosome, chromosome 12, are most clearly associated with the maleness factor, albeit incompletely. We also confirmed occasional crossovers proximal to the male-determining region, defining a second, rarely recombining, pseudo-autosomal region, PAR2. This fish species may therefore have no completely male-specific region (MSY) more extensive than the male-determining factor. The positions of the few crossover events suggest a location for the male-determining factor within a physically small repetitive region. A sex-reversed XX male had few crossovers in PAR2, suggesting that this region's low crossover rate depends on the phenotypic, not the genetic, sex. Thus, rare individuals whose phenotypic and genetic sexes differ, and/or occasional PAR2 crossovers in males can explain the failure to detect fully Y-linked variants.

Epigenetic patterns in Atlantic herring (Clupea harengus): Temperature and photoperiod as environmental stressors during larval development

Understanding the molecular mechanisms underlying individual responses to environmental changes is crucial for species conservation and management. Pelagic fishes including Atlantic herring (Clupea harengus) are of particular interest because of their key ecological and economic roles and their susceptibility to a changing ocean from global warming. Temperature and photoperiod have been linked with spawning time and location in adult herring, but no study has thus far investigated the role of environmental factors on gene regulation during the vulnerable early developmental stages. Here, we examine DNA methylation patterns of larval herring bred under two temperatures (11°C and 13°C) and photoperiod (6 and 12 h) regimes in a 2 × 2 factorial design. We found consistently high levels of global methylation across all individuals and a decline in global methylation with increased developmental stage that was more pronounced at 13°C (p ≤ 0.007) than at 11°C (p ≥ 0.21). Most of the differentially methylated sites were in exon and promoter regions for genes linked to metabolism and development, some of which were hypermethylated at higher temperature. These results demonstrate the important role of DNA methylation during larval development and suggest that this molecular mechanism might be key in regulating early-stage responses to environmental stressors in Atlantic herring.

Non-synonymous variation and protein structure of candidate genes associated with selection in farm and wild populations of turbot (Scophthalmus maximus)

Non-synonymous variation (NSV) of protein coding genes represents raw material for selection to improve adaptation to the diverse environmental scenarios in wild and livestock populations. Many aquatic species face variations in temperature, salinity and biological factors throughout their distribution range that is reflected by the presence of allelic clines or local adaptation. The turbot (Scophthalmus maximus) is a flatfish of great commercial value with a flourishing aquaculture which has promoted the development of genomic resources. In this study, we developed the first atlas of NSVs in the turbot genome by resequencing 10 individuals from Northeast Atlantic Ocean. More than 50,000 NSVs where detected in the ~ 21,500 coding genes of the turbot genome, and we selected 18 NSVs to be genotyped using a single Mass ARRAY multiplex on 13 wild populations and three turbot farms. We detected signals of divergent selection on several genes related to growth, circadian rhythms, osmoregulation and oxygen binding in the different scenarios evaluated. Furthermore, we explored the impact of NSVs identified on the 3D structure and functional relationship of the correspondent proteins. In summary, our study provides a strategy to identify NSVs in species with consistently annotated and assembled genomes to ascertain their role in adaptation.

Independent Evolution of Sex Chromosomes and Male Pregnancy-Related Genes in Two Seahorse Species

Unlike birds and mammals, many teleosts have homomorphic sex chromosomes, and changes in the chromosome carrying the sex-determining locus, termed "turnovers", are common. Recent turnovers allow studies of several interesting questions. One question is whether the new sex-determining regions evolve to become completely non-recombining, and if so, how and why. Another is whether (as predicted) evolutionary changes that benefit one sex accumulate in the newly sex-linked region. To study these questions, we analyzed the genome sequences of two seahorse species of the Syngnathidae, a fish group in which many species evolved a unique structure, the male brood pouch. We find that both seahorse species have XY sex chromosome systems, but their sex chromosome pairs are not homologs, implying that at least one turnover event has occurred. The Y-linked regions occupy 63.9% and 95.1% of the entire sex chromosome of the two species and do not exhibit extensive sequence divergence with their X-linked homologs. We find evidence for occasional recombination between the extant sex chromosomes that may account for their homomorphism. We argue that these Y-linked regions did not evolve by recombination suppression after the turnover, but by the ancestral nature of the low crossover rates in these chromosome regions. With such an ancestral crossover landscape, a turnover can instantly create an extensive Y-linked region. Finally, we test for adaptive evolution of male pouch-related genes after they became Y-linked in the seahorse.

The role of genomic signatures of directional selection and demographic history in the population structure of a marine teleost with high gene flow

Recent studies have uncovered patterns of genomic divergence in marine teleosts where panmixia due to high gene flow has been the general paradigm. These signatures of divergent selection are often impacted by structural variants, acting as "supergenes" facilitating local adaptation. The highly dispersing European plaice (Pleuronectes platessa)-in which putative structural variants (i.e., inversions) have been identified-has successfully colonized the brackish water ecosystem of the Baltic Sea. Thus, the species represents an ideal opportunity to investigate how the interplay of gene flow, structural variants, natural selection, past demographic history, and gene flow impacts on population (sub)structuring in marine systems. Here, we report on the generation of an annotated draft plaice genome assembly in combination with population sequencing data-following the salinity gradient from the Baltic Sea into the North Sea together with samples from Icelandic waters-to illuminate genome-wide patterns of divergence. Neutral markers pointed at large-scale panmixia across the European continental shelf associated with high gene flow and a common postglacial colonization history of shelf populations. However, based on genome-wide outlier loci, we uncovered signatures of population substructuring among the European continental shelf populations, i.e., suggesting signs of ongoing selection. Genome-wide selection analyses (xp-EHH) and the identification of genes within genomic regions of recent selective sweeps-overlapping with the outlier loci-suggest that these represent the signs of divergent selection. Our findings provide support for genomic divergence driven by local adaptation in the face of high gene flow and elucidate the relative importance of demographic history versus adaptive divergence in shaping the contemporary population genetic structure of a marine teleost. The role of the putative inversion(s) in the substructuring-and potentially ongoing adaptation-was seemingly not substantial.

Population dynamics of Baltic herring since the Viking Age revealed by ancient DNA and genomics

The world's oceans are currently facing major stressors in the form of overexploitation and anthropogenic climate change. The Baltic Sea was home to the first "industrial" fishery ∼800 y ago targeting the Baltic herring, a species that is still economically and culturally important today. Yet, the early origins of marine industries and the long-term ecological consequences of historical and contemporary fisheries remain debated. Here, we study long-term population dynamics of Baltic herring to evaluate the past impacts of humans on the marine environment. We combine modern whole-genome data with ancient DNA (aDNA) to identify the earliest-known long-distance herring trade in the region, illustrating that extensive fish trade began during the Viking Age. We further resolve population structure within the Baltic and observe demographic independence for four local herring stocks over at least 200 generations. It has been suggested that overfishing at Øresund in the 16th century resulted in a demographic shift from autumn-spawning to spring-spawning herring dominance in the Baltic. We show that while the Øresund fishery had a negative impact on the western Baltic herring stock, the demographic shift to spring-spawning dominance did not occur until the 20th century. Instead, demographic reconstructions reveal population trajectories consistent with expected impacts of environmental change and historical reports on shifting fishing targets over time. This study illustrates the joint impact of climate change and human exploitation on marine species as well as the role historical ecology can play in conservation and management policies.

Ensembl 2023

Ensembl (https://www.ensembl.org) has produced high-quality genomic resources for vertebrates and model organisms for more than twenty years. During that time, our resources, services and tools have continually evolved in line with both the publicly available genome data and the downstream research and applications that utilise the Ensembl platform. In recent years we have witnessed a dramatic shift in the genomic landscape. There has been a large increase in the number of high-quality reference genomes through global biodiversity initiatives. In parallel, there have been major advances towards pangenome representations of higher species, where many alternative genome assemblies representing different breeds, cultivars, strains and haplotypes are now available. In order to support these efforts and accelerate downstream research, it is our goal at Ensembl to create high-quality annotations, tools and services for species across the tree of life. Here, we report our resources for popular reference genomes, the dramatic growth of our annotations (including haplotypes from the first human pangenome graphs), updates to the Ensembl Variant Effect Predictor (VEP), interactive protein structure predictions from AlphaFold DB, and the beta release of our new website.

A baseline for the genetic stock identification of Atlantic herring, Clupea harengus, in ICES Divisions 6.a, 7.b-c

Atlantic herring in International Council for Exploration of the Sea (ICES) Divisions 6.a, 7.b-c comprises at least three populations, distinguished by temporal and spatial differences in spawning, which have until recently been managed as two stocks defined by geographical delineators. Outside of spawning the populations form mixed aggregations, which are the subject of acoustic surveys. The inability to distinguish the populations has prevented the development of separate survey indices and separate stock assessments. A panel of 45 single-nucleotide polymorphisms, derived from whole-genome sequencing, were used to genotype 3480 baseline spawning samples (2014-2021). A temporally stable baseline comprising 2316 herring from populations known to inhabit Division 6.a was used to develop a genetic assignment method, with a self-assignment accuracy greater than 90%. The long-term temporal stability of the assignment model was validated by assigning archive (2003-2004) baseline samples (270 individuals) with a high level of accuracy. Assignment of non-baseline samples (1514 individuals) from Divisions 6.a, 7.b-c indicated previously unrecognized levels of mixing of populations outside of the spawning season. The genetic markers and assignment models presented constitute a 'toolbox' that can be used for the assignment of herring caught in mixed survey and commercial catches in Division 6.a into their population of origin with a high level of accuracy.

Transcriptomics of a Greenlandic Snailfish Reveals Exceptionally High Expression of Antifreeze Protein Transcripts

Polar fishes have evolved antifreeze proteins (AFPs) that allow them to survive in subzero temperatures. We performed deep transcriptomic sequencing on a postlarval/juvenile variegated snailfish, Liparis gibbus (Actinopterygii: Scorpaeniformes: Cottoidei: Liparidae), living in an iceberg habitat (-2°C) in Eastern Greenland and report detection of highly expressed transcripts that code for putative AFPs from 2 gene families, Type I and LS-12-like proteins (putative Type IV AFPs). The transcripts encoding both proteins have expression levels among the top <1% of expressed genes in the fish. The Type I AFP sequence is different from a reported Type I AFP from the same species, possibly expressed from a different genetic locus. While prior findings from related adult sculpins suggest that LS-12-like/Type IV AFPs may not have a role in antifreeze protection, our finding of very high relative gene expression of the LS-12-like gene suggests that highly active transcription of the gene is important to the fish in the iceberg habitat and raises the possibility that weak or combinatorial antifreeze activity could be beneficial. These findings highlight the physiological importance of antifreeze proteins to the survival of fishes living in polar habitats.

A chromosome-level genome assembly of Paracymoriza distinctalis (Lepidoptera: Crambidae: Acentropinae)

Paracymoriza distinctalis is a semiaquatic lepidopteran insect, which is of great value for studying the differentiation of the Pyraloidea superfamily. However, the understanding of heredity, evolution, and functional genomics of P. distinctalis are limited by few genome-wide resources. Here, we applied PacBio sequencing and the chromosome capture technique to assemble the first P. distinctalis genome from a single female individual. The genome size is 1.2 Gb with 32 chromosomes and the N50 is 38.91 Mb. Approximately 576.37 Mb, accounting for 48.93% of the genome, was identified as repeats. The genome comprises 39,003 protein-coding genes, 66.56% of which were functionally annotated. Comparative genomics analysis suggested that the common ancestor of P. distinctalis and Chilo suppressalis lived ~83.5 million years ago. This chromosome-level genome assembly work is not only conducive to the understanding of P. distinctalis, but also may promote the study of the genomes of other lepidopteran species.

Reflections on antifreeze proteins and their evolution

The discovery of radically different antifreeze proteins (AFPs) in fishes during the 1970s and 1980s suggested that these proteins had recently and independently evolved to protect teleosts from freezing in icy seawater. Early forays into the isolation and characterization of AFP genes in these fish showed they were massively amplified, often in long tandem repeats. The work of many labs in the 1980s onward led to the discovery and characterization of AFPs in other kingdoms, such as insects, plants, and many different microorganisms. The distinct ice-binding property that these ice-binding proteins (IBPs) share has facilitated their purification through adsorption to ice, and the ability to produce recombinant versions of IBPs has enabled their structural characterization and the mapping of their ice-binding sites (IBSs) using site-directed mutagenesis. One hypothesis for their ice affinity is that the IBS organizes surface waters into an ice-like pattern that freezes the protein onto ice. With access now to a rapidly expanding database of genomic sequences, it has been possible to trace the origins of some fish AFPs through the process of gene duplication and divergence, and to even show the horizontal transfer of an AFP gene from one species to another.

The Genetic Population Structure of Lake Tanganyika's Lates Species Flock, an Endemic Radiation of Pelagic Top Predators

Understanding genetic connectivity plays a crucial role in species conservation decisions, and genetic connectivity is an important component of modern fisheries management. In this study, we investigated the population genetics of four endemic Lates species of Lake Tanganyika (Lates stappersii, L. microlepis, L. mariae, and L. angustifrons) using reduced-representation genomic sequencing methods. We find the four species to be strongly differentiated from one another (mean interspecific FST = 0.665), with no evidence for contemporary admixture. We also find evidence for strong genetic structure within L. mariae, with the majority of individuals from the most southern sampling site forming a genetic group that is distinct from the individuals at other sampling sites. We find evidence for much weaker structure within the other three species (L. stappersii, L. microlepis, and L. angustifrons). Our ability to detect this weak structure despite small and unbalanced sample sizes and imprecise geographic sampling locations suggests the possibility for further structure undetected in our study. We call for further research into the origins of the genetic differentiation in these four species-particularly that of L. mariae-which may be important for conservation and management of this culturally and economically important clade of fishes.

Atlantic herring (Clupea harengus) population structure in the Northeast Atlantic Ocean

The Atlantic herring Clupea harengus L has a vast geographical distribution and a complex population structure with a few very large migratory units and many small local populations. Each population has its own spawning ground and/or time, thereby maintaining their genetic integrity. Several herring populations migrate between common feeding grounds and over-wintering areas resulting in frequent mixing of populations. Thus, many herring fisheries are based on mixed populations of different demographic status. In order to avoid over-exploitation of weak populations and to conserve biodiversity, understanding the population structure and population mixing is important for maintaining biologically sustainable herring fisheries. The aim of this study was to investigate the genetic population structure of herring in the Faroese and surrounding waters, and to develop genetic markers for distinguishing between four herring management units (often called stocks), namely the Norwegian spring-spawning herring (NSSH), Icelandic summer-spawning herring (ISSH), North Sea autumn-spawning herring (NSAH), and Faroese autumn-spawning herring (FASH). Herring from the four stocks were sequenced at low coverage, and single nucleotide polymorphisms (SNPs) were called and used for population structure analysis and individual assignment. An ancestry-informative SNP panel with 118 SNPs was developed and tested on 240 individuals. The results showed that all four stocks appeared to be genetically differentiated populations, but at lower levels of differentiation between FASH and ISSH than the other two populations. Overall assignment rate with the SNP panel was 80.7%, and agreement between the genetic and traditional visual assignment was 75.5%. The NSAH and NSSH samples had the highest assignment rate (100% and 98.3%, respectively) and highest agreement between traditional and genetic assignment methods (96.6% and 94.9%, respectively). The FASH and ISSH samples had substantially lower assignment rates (72.9% and 51.7%, respectively) and agreement between traditional and genetic methods (39.5% and 48.4%, respectively).

Executive function deficit in bipolar offspring: A neurocognitive endophenotype?

BACKGROUND

Recent studies in bipolar offspring (BO) showed that a low cognitive performance, especially executive function deficit, could be an early marker of bipolar disorder (BD). Nevertheless, these findings have not been replicated (specifically attentional control, flexibility, and working memory). In addition, most studies have focused on children and adolescents, but few studies analyze the executive function performance in BO adults.

OBJECTIVE

Our goal was to compare the neurocognitive performance of BO with control parent-offspring (CO) in a sample that included various age groups.

METHOD

We conducted a cohort study, including subjects between six to 30 years old. We evaluated 129 BO and 113 CO subjects using validated psychiatric diagnostic interviews and an extensive neuropsychological battery.

RESULTS

Compared to the CO group, the BO group presented a lower performance in several executive functioning domains, mainly in tasks of attentional control, flexibility, and working memory. All age groups exhibited these findings.

CONCLUSIONS

BO group presents executive function deficits, regardless of the age group: children, adolescents, and adults. This neurocognitive deficit should be accountable as a neurocognitive endophenotype candidate in BD.

A highly contiguous nuclear genome assembly of the mandarinfish Synchiropus splendidus (Syngnathiformes: Callionymidae)

The fish order Syngnathiformes has been referred to as a collection of misfit fishes, comprising commercially important fish such as red mullets as well as the highly diverse seahorses, pipefishes, and seadragons-the well-known family Syngnathidae, with their unique adaptations including male pregnancy. Another ornate member of this order is the species mandarinfish. No less than two types of chromatophores have been discovered in the spectacularly colored mandarinfish: the cyanophore (producing blue color) and the dichromatic cyano-erythrophore (producing blue and red). The phylogenetic position of mandarinfish in Syngnathiformes, and their promise of additional genetic discoveries beyond the chromatophores, made mandarinfish an appealing target for whole-genome sequencing. We used linked sequences to create synthetic long reads, producing a highly contiguous genome assembly for the mandarinfish. The genome assembly comprises 483 Mbp (longest scaffold 29 Mbp), has an N50 of 12 Mbp, and an L50 of 14 scaffolds. The assembly completeness is also high, with 92.6% complete, 4.4% fragmented, and 2.9% missing out of 4584 BUSCO genes found in ray-finned fishes. Outside the family Syngnathidae, the mandarinfish represents one of the most contiguous syngnathiform genome assemblies to date. The mandarinfish genomic resource will likely serve as a high-quality outgroup to syngnathid fish, and furthermore for research on the genomic underpinnings of the evolution of novel pigmentation.

An accurate assignment test for extremely low-coverage whole-genome sequence data

Genomic assignment tests can provide important diagnostic biological characteristics, such as population of origin or ecotype. Yet, assignment tests often rely on moderate- to high-coverage sequence data that can be difficult to obtain for fields such as molecular ecology and ancient DNA. We have developed a novel approach that efficiently assigns biologically relevant information (i.e., population identity or structural variants such as inversions) in extremely low-coverage sequence data. First, we generate databases from existing reference data using a subset of diagnostic single nucleotide polymorphisms (SNPs) associated with a biological characteristic. Low-coverage alignment files are subsequently compared to these databases to ascertain allelic state, yielding a joint probability for each association. To assess the efficacy of this approach, we assigned haplotypes and population identity in Heliconius butterflies, Atlantic herring, and Atlantic cod using chromosomal inversion sites and whole-genome data. We scored both modern and ancient specimens, including the first whole-genome sequence data recovered from ancient Atlantic herring bones. The method accurately assigns biological characteristics, including population membership, using extremely low-coverage data (as low as 0.0001x) based on genome-wide SNPs. This approach will therefore increase the number of samples in evolutionary, ecological and archaeological research for which relevant biological information can be obtained.

Genetic structure of marine and lake forms of Pacific herring Clupea pallasii

The Pacific herring (Clupea pallasii) is one of the most important species in the commercial fisheries distributed in the North Pacific Ocean and the northeastern European seas. This teleost has marine and lake ecological forms a long its distribution in the Holarctic. However, the level of genetic differentiation between these two forms is not well known. In the present study, we used ddRAD-sequencing to genotype 54 specimens from twelve wild Pacific herring populations from the Kara Sea and the Russian part of the northwestern Pacific Ocean for unveiling the genetic structure of Pacific herring. We found that the Kara Sea population is significantly distinct from Pacific Ocean populations. It was demonstrated that lake populations of Pacific herring differ from one another as well as from marine specimens. Our results show that fresh and brackish water Pacific herring, which inhabit lakes, can be distinguished as a separate lake ecological form. Moreover, we demonstrate that each observed lake Pacific herring population has its own and unique genetic legacy.

Genome structural variation in human evolution

Structural variation (SV) is a large difference (typically >100 bp) in the genomic structure of two genomes and includes both copy number variation and variation that does not change copy number of a genomic region, such as an inversion. Improved reference genomes, combined with widespread genome sequencing using short-read sequencing technology, and increasingly using long-read sequencing, have reignited interest in SV. Recent large-scale studies and functional focused analyses have highlighted the role of SV in human evolution. In this review, we highlight human-specific SVs involved in changes in the brain, population-specific SVs that affect response to the environment, including adaptation to diet and infectious diseases, and summarise the contribution of archaic hominin admixture to present-day human SV.

Genomic Architecture of Rapid Parallel Adaptation to Fresh Water in a Wild Fish

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. However, the genetic architecture underlying these adaptive changes is still poorly understood. Using population genomic approaches, we investigated the genomic architecture that underlies rapid parallel adaptation of Coilia nasus to fresh water by comparing four freshwater-resident populations with their ancestral anadromous population. Linkage disequilibrium network analysis and population genetic analyses revealed two putative large chromosome inversions on LG6 and LG22, which were enriched for outlier loci and exhibited parallel association with freshwater adaptation. Drastic frequency shifts and elevated genetic differentiation were observed for the two chromosome inversions among populations, suggesting that both inversions would undergo divergent selection between anadromous and resident ecotypes. Enrichment analysis of genes within chromosome inversions showed significant enrichment of genes involved in metabolic process, immunoregulation, growth, maturation, osmoregulation, and so forth, which probably underlay differences in morphology, physiology and behavior between the anadromous and freshwater-resident forms. The availability of beneficial standing genetic variation, large optimum shift between marine and freshwater habitats, and high efficiency of selection with large population size could lead to the observed rapid parallel adaptive genomic change. We propose that chromosomal inversions might have played an important role during the evolution of rapid parallel ecological divergence in the face of environmental heterogeneity in C. nasus. Our study provides insights into the genomic basis of rapid adaptation of complex traits in novel habitats and highlights the importance of structural genomic variants in analyses of ecological adaptation.

Evolution of master sex determiners: TGF-β signalling pathways at regulatory crossroads

To date, more than 20 different vertebrate master sex-determining genes have been identified on different sex chromosomes of mammals, birds, frogs and fish. Interestingly, six of these genes are transcription factors (Dmrt1- or Sox3- related) and 13 others belong to the TGF-β signalling pathway (Amh, Amhr2, Bmpr1b, Gsdf and Gdf6). This pattern suggests that only a limited group of factors/signalling pathways are prone to become top regulators again and again. Although being clearly a subordinate member of the sex-regulatory network in mammals, the TGF-β signalling pathway made it to the top recurrently and independently. Facing this rolling wave of TGF-β signalling pathways, this review will decipher how the TGF-β signalling pathways cope with the canonical sex gene regulatory network and challenge the current evolutionary concepts accounting for the diversity of sex-determining mechanisms. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Distinct genetic clustering in the weakly differentiated polar cod, Boreogadus saida Lepechin, 1774 from East Siberian Sea to Svalbard

The cold-adapted polar cod Boreogadus saida, a key species in Arctic ecosystems, is vulnerable to global warming and ice retreat. In this study, 1257 individuals sampled in 17 locations within the latitudinal range of 75–81°N from Svalbard to East Siberian Sea were genotyped with a dedicated suite of 116 single-nucleotide polymorphic loci (SNP). The overall pattern of isolation by distance (IBD) found was driven by the two easternmost samples (East Siberian Sea and Laptev Sea), whereas no differentiation was registered in the area between the Kara Sea and Svalbard. Eleven SNP under strong linkage disequilibrium, nine of which could be annotated to chromosome 2 in Atlantic cod, defined two genetic groups of distinct size, with the major cluster containing seven-fold larger number of individuals than the minor. No underlying geographic basis was evident, as both clusters were detected throughout all sampling sites in relatively similar proportions (i.e. individuals in the minor cluster ranging between 4 and 19% on the location basis). Similarly, females and males were also evenly distributed between clusters and age groups. A differentiation was, however, found regarding size at age: individuals belonging to the major cluster were significantly longer in the second year. This study contributes to increasing the population genetic knowledge of this species and suggests that an appropriate management should be ensured to safeguard its diversity.

Functional differences between TSHR alleles associate with variation in spawning season in Atlantic herring

The underlying molecular mechanisms that determine long day versus short day breeders remain unknown in any organism. Atlantic herring provides a unique opportunity to examine the molecular mechanisms involved in reproduction timing, because both spring and autumn spawners exist within the same species. Although our previous whole genome comparisons revealed a strong association of TSHR alleles with spawning seasons, the functional consequences of these variants remain unknown. Here we examined the functional significance of six candidate TSHR mutations strongly associated with herring reproductive seasonality. We show that the L471M missense mutation in the spring-allele causes enhanced cAMP signaling. The best candidate non-coding mutation is a 5.2 kb retrotransposon insertion upstream of the TSHR transcription start site, near an open chromatin region, which is likely to affect TSHR expression. The insertion occurred prior to the split between Pacific and Atlantic herring and was lost in the autumn-allele. Our study shows that strongly associated coding and non-coding variants at the TSHR locus may both contribute to the regulation of seasonal reproduction in herring.

Junfeng Chen et al. examine potential functional consequences of reproduction timing-associated TSHR alleles segregating in Atlantic herring. By comparing fish that spawn during the spring to those that spawn in the autumn, they find that the spring-allele is correlated with enhanced cAMP signaling and that both coding and non-coding variants in the TSHR locus contribute to seasonal reproduction.

Chromosome-Level Assembly of the Atlantic Silverside Genome Reveals Extreme Levels of Sequence Diversity and Structural Genetic Variation

The levels and distribution of standing genetic variation in a genome can provide a wealth of insights about the adaptive potential, demographic history, and genome structure of a population or species. As structural variants are increasingly associated with traits important for adaptation and speciation, investigating both sequence and structural variation is essential for wholly tapping this potential. Using a combination of shotgun sequencing, 10x Genomics linked reads and proximity-ligation data (Chicago and Hi-C), we produced and annotated a chromosome-level genome assembly for the Atlantic silverside (Menidia menidia)-an established ecological model for studying the phenotypic effects of natural and artificial selection-and examined patterns of genomic variation across two individuals sampled from different populations with divergent local adaptations. Levels of diversity varied substantially across each chromosome, consistently being highly elevated near the ends (presumably near telomeric regions) and dipping to near zero around putative centromeres. Overall, our estimate of the genome-wide average heterozygosity in the Atlantic silverside is among the highest reported for a fish, or any vertebrate (1.32-1.76% depending on inference method and sample). Furthermore, we also found extreme levels of structural variation, affecting ∼23% of the total genome sequence, including multiple large inversions (> 1 Mb and up to 12.6 Mb) associated with previously identified haploblocks showing strong differentiation between locally adapted populations. These extreme levels of standing genetic variation are likely associated with large effective population sizes and may help explain the remarkable adaptive divergence among populations of the Atlantic silverside.

Accessing Livestock Resources in Ensembl

Genome assembly is cheaper, more accurate and more automated than it has ever been. This is due to a combination of more cost-efficient chemistries, new sequencing technologies and better algorithms. The livestock community has been at the forefront of this new wave of genome assembly, generating some of the highest quality vertebrate genome sequences. Ensembl's goal is to add functional and comparative annotation to these genomes, through our gene annotation, genomic alignments, gene trees, regulatory, and variation data. We run computationally complex analyses in a high throughput and consistent manner to help accelerate downstream science. Our livestock resources are continuously growing in both breadth and depth. We annotate reference genome assemblies for newly sequenced species and regularly update annotation for existing genomes. We are the only major resource to support the annotation of breeds and other non-reference assemblies. We currently provide resources for 13 pig breeds, maternal and paternal haplotypes for hybrid cattle and various other non-reference or wild type assemblies for livestock species. Here, we describe the livestock data present in Ensembl and provide protocols for how to view data in our genome browser, download via it our FTP site, manipulate it via our tools and interact with it programmatically via our REST API.

Evolution at two time frames: ancient structural variants involved in post-glacial divergence of the European plaice (Pleuronectes platessa)

Changing environmental conditions can lead to population diversification through differential selection on standing genetic variation. Structural variant (SV) polymorphisms provide examples of ancient alleles that in time become associated with novel environmental gradients. The European plaice (Pleuronectes platessa) is a marine flatfish showing large allele-frequency differences at two putative SVs associated with environmental variation. In this study, we explored the contribution of these SVs to population structure across the North East Atlantic. We compared genome-wide population structure using sets of RAD-sequencing SNPs with the spatial structure of the SVs. We found that in contrast to the rest of the genome, the SVs were only weakly associated with an isolation-by-distance pattern. Indeed, both SVs showed important variation in haplogroup frequencies, with the same haplogroup increasing both along the salinity gradient of the Baltic Sea, and found in high frequency in the northern-range margin of the Atlantic. Phylogenetic analyses suggested that the SV alleles are much older than the age of the Baltic Sea itself. These results suggest that the SVs are older than the age of the environmental gradients with which they currently co-vary. Altogether, our results suggest that the plaice SVs were shaped by evolutionary processes occurring at two time frames, firstly following their origin, ancient spread and maintenance in the ancestral populations, and secondly related to their current association with more recently formed environmental gradients such as those found in the North Sea-Baltic Sea transition zone.

The Genomic Architecture and Evolutionary Fates of Supergenes

Supergenes are genomic regions containing sets of tightly linked loci that control multi-trait phenotypic polymorphisms under balancing selection. Recent advances in genomics have uncovered significant variation in both the genomic architecture as well as the mode of origin of supergenes across diverse organismal systems. Although the role of genomic architecture for the origin of supergenes has been much discussed, differences in the genomic architecture also subsequently affect the evolutionary trajectory of supergenes and the rate of degeneration of supergene haplotypes. In this review, we synthesize recent genomic work and historical models of supergene evolution, highlighting how the genomic architecture of supergenes affects their evolutionary fate. We discuss how recent findings on classic supergenes involved in governing ant colony social form, mimicry in butterflies, and heterostyly in flowering plants relate to theoretical expectations. Furthermore, we use forward simulations to demonstrate that differences in genomic architecture affect the degeneration of supergenes. Finally, we discuss implications of the evolution of supergene haplotypes for the long-term fate of balanced polymorphisms governed by supergenes.

Horizontal Gene Transfer in Vertebrates: A Fishy Tale

The recent assembly of the herring genome suggests this fish acquired its antifreeze protein gene by horizontal transfer and then passed a copy on to the smelt. The direction of gene transfer is confirmed by some accompanying transposable elements and by the breakage of gene synteny.

A New Zebrafish Model for Pseudoxanthoma Elasticum

Calcification of various tissues is a significant health issue associated with aging, cancer and autoimmune diseases. There are both environmental and genetic factors behind this phenomenon and understanding them is essential for the development of efficient therapeutic approaches. Pseudoxanthoma elasticum (PXE) is a rare genetic disease, a prototype for calcification disorders, resulting from the dysfunction of ABCC6, a transport protein found in the membranes of cells. It is identified by excess calcification in a variety of tissues (e.g., eyes, skin, arteries) and currently it has no cure, known treatments target the symptoms only. Preclinical studies of PXE have been successful in mice, proving the usefulness of animal models for the study of the disease. Here, we present a new zebrafish (Danio rerio) model for PXE. By resolving some ambiguous assemblies in the zebrafish genome, we show that there are two functional and one non-functional paralogs for ABCC6 in zebrafish (abcc6a, abcc6b.1, and abcc6b.2, respectively). We created single and double mutants for the functional paralogs and characterized their calcification defects with a combination of techniques. Zebrafish deficient in abcc6a show defects in their vertebral calcification and also display ectopic calcification foci in their soft tissues. Our results also suggest that the impairment of abcc6b.1 does not affect this biological process.

Functional genetic diversity in an exploited marine species and its relevance to fisheries management

The timing of reproduction influences key evolutionary and ecological processes in wild populations. Variation in reproductive timing may be an especially important evolutionary driver in the marine environment, where the high mobility of many species and few physical barriers to migration provide limited opportunities for spatial divergence to arise. Using genomic data collected from spawning aggregations of Pacific herring (Clupea pallasii) across 1600 km of coastline, we show that reproductive timing drives population structure in these pelagic fish. Within a specific spawning season, we observed isolation by distance, indicating that gene flow is also geographically limited over our study area. These results emphasize the importance of considering both seasonal and spatial variation in spawning when delineating management units for herring. On several chromosomes, we detected linkage disequilibrium extending over multiple Mb, suggesting the presence of chromosomal rearrangements. Spawning phenology was highly correlated with polymorphisms in several genes, in particular SYNE2, which influences the development of retinal photoreceptors in vertebrates. SYNE2 is probably within a chromosomal rearrangement in Pacific herring and is also associated with spawn timing in Atlantic herring (Clupea harengus). The observed genetic diversity probably underlies resource waves provided by spawning herring. Given the ecological, economic and cultural significance of herring, our results support that conserving intraspecific genetic diversity is important for maintaining current and future ecosystem processes.

Ecological adaptation in European eels is based on phenotypic plasticity

The relative role of genetic adaptation and phenotypic plasticity is of fundamental importance in evolutionary ecology [M. J. West-Eberhard, Proc. Natl. Acad. Sci. U.S.A. 102 (suppl. 1), 6543-6549 (2005)]. European eels have a complex life cycle, including transitions between life stages across ecological conditions in the Sargasso Sea, where spawning occurs, and those in brackish and freshwater bodies from northern Europe to northern Africa. Whether continental eel populations consist of locally adapted and genetically distinct populations or comprise a single panmictic population has received conflicting support. Here we use whole-genome sequencing and show that European eels belong to one panmictic population. A complete lack of geographical genetic differentiation is demonstrated. We postulate that this is possible because the most critical life stages-spawning and embryonic development-take place under near-identical conditions in the Sargasso Sea. We further show that within-generation selection, which has recently been proposed as a mechanism for genetic adaptation in eels, can only marginally change allele frequencies between cohorts of eels from different geographic regions. Our results strongly indicate plasticity as the predominant mechanism for how eels respond to diverse environmental conditions during postlarval stages, ultimately solving a long-standing question for a classically enigmatic species.

Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird-of-paradise

Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat-rich and GC-rich regions (genomic "dark matter") limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long-read, linked-read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC-rich microchromosomes and the repeat-rich W chromosome. Telomere-to-telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes.

Ecological adaptation in Atlantic herring is associated with large shifts in allele frequencies at hundreds of loci

Atlantic herring is widespread in North Atlantic and adjacent waters and is one of the most abundant vertebrates on earth. This species is well suited to explore genetic adaptation due to minute genetic differentiation at selectively neutral loci. Here, we report hundreds of loci underlying ecological adaptation to different geographic areas and spawning conditions. Four of these represent megabase inversions confirmed by long read sequencing. The genetic architecture underlying ecological adaptation in herring deviates from expectation under a classical infinitesimal model for complex traits because of large shifts in allele frequencies at hundreds of loci under selection.

Structure and Characterization of Phosphoglucomutase 5 from Atlantic and Baltic Herring-An Inactive Enzyme with Intact Substrate Binding

Phosphoglucomutase 5 (PGM5) in humans is known as a structural muscle protein without enzymatic activity, but detailed understanding of its function is lacking. PGM5 belongs to the alpha-D-phosphohexomutase family and is closely related to the enzymatically active metabolic enzyme PGM1. In the Atlantic herring, Clupea harengus, PGM5 is one of the genes strongly associated with ecological adaptation to the brackish Baltic Sea. We here present the first crystal structures of PGM5, from the Atlantic and Baltic herring, differing by a single substitution Ala330Val. The structure of PGM5 is overall highly similar to structures of PGM1. The structure of the Baltic herring PGM5 in complex with the substrate glucose-1-phosphate shows conserved substrate binding and active site compared to human PGM1, but both PGM5 variants lack phosphoglucomutase activity under the tested conditions. Structure comparison and sequence analysis of PGM5 and PGM1 from fish and mammals suggest that the lacking enzymatic activity of PGM5 is related to differences in active-site loops that are important for flipping of the reaction intermediate. The Ala330Val substitution does not alter structure or biophysical properties of PGM5 but, due to its surface-exposed location, could affect interactions with protein-binding partners.

Strong population differentiation in lingcod (Ophiodon elongatus) is driven by a small portion of the genome

Delimiting intraspecific genetic variation in harvested species is crucial to the assessment of population status for natural resource management and conservation purposes. Here, we evaluated genetic population structure in lingcod (Ophiodon elongatus), a commercially and recreationally important fishery species along the west coast of North America. We used 16,749 restriction site-associated DNA sequencing (RADseq) markers, in 611 individuals collected from across the bulk of the species range from Southeast Alaska to Baja California, Mexico. In contrast to previous population genetic work on this species, we found strong evidence for two distinct genetic clusters. These groups separated latitudinally with a break near Point Reyes off Northern California, and there was a high frequency of admixed individuals in close proximity to the break. F-statistics corroborate this genetic break between northern and southern sampling sites, although most loci are characterized by low FST values, suggesting high gene flow throughout most of the genome. Outlier analyses identified 182 loci putatively under divergent selection, most of which mapped to a single genomic region. When individuals were grouped by cluster assignment (northern, southern, and admixed), 71 loci were fixed between the northern and southern cluster, all of which were identified in the outlier scans. All individuals identified as admixed exhibited near 50:50 assignment to northern and southern clusters and were heterozygous for most fixed loci. Alignments of RADseq loci to a draft lingcod genome assembly and three other teleost genomes with chromosome-level assemblies suggest that outlier and fixed loci are concentrated on a single chromosome. Similar genomic patterns have been attributed to chromosomal inversions in diverse taxonomic groups. Regardless of the evolutionary mechanism, these results represent novel observations of genetic structure in lingcod and designate clear evolutionary units that could be used to inform fisheries management.

Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal, poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab-adapted ecotypes that exhibit partial reproductive isolation from wave-adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co-association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes-consistent with models of adaptation with gene flow-but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co-association of loci during ecological speciation with ongoing gene flow.

Vertebrate Chromosome Evolution

The study of chromosome evolution is undergoing a resurgence of interest owing to advances in DNA sequencing technology that facilitate the production of chromosome-scale whole-genome assemblies de novo. This review focuses on the history, methods, discoveries, and current challenges facing the field, with an emphasis on vertebrate genomes. A detailed examination of the literature on the biology of chromosome rearrangements is presented, specifically the relationship between chromosome rearrangements and phenotypic evolution, adaptation, and speciation. A critical review of the methods for identifying, characterizing, and visualizing chromosome rearrangements and computationally reconstructing ancestral karyotypes is presented. We conclude by looking to the future, identifying the enormous technical and scientific challenges presented by the accumulation of hundreds and eventually thousands of chromosome-scale assemblies.

Twelve quick steps for genome assembly and annotation in the classroom

Eukaryotic genome sequencing and de novo assembly, once the exclusive domain of well-funded international consortia, have become increasingly affordable, thus fitting the budgets of individual research groups. Third-generation long-read DNA sequencing technologies are increasingly used, providing extensive genomic toolkits that were once reserved for a few select model organisms. Generating high-quality genome assemblies and annotations for many aquatic species still presents significant challenges due to their large genome sizes, complexity, and high chromosome numbers. Indeed, selecting the most appropriate sequencing and software platforms and annotation pipelines for a new genome project can be daunting because tools often only work in limited contexts. In genomics, generating a high-quality genome assembly/annotation has become an indispensable tool for better understanding the biology of any species. Herein, we state 12 steps to help researchers get started in genome projects by presenting guidelines that are broadly applicable (to any species), sustainable over time, and cover all aspects of genome assembly and annotation projects from start to finish. We review some commonly used approaches, including practical methods to extract high-quality DNA and choices for the best sequencing platforms and library preparations. In addition, we discuss the range of potential bioinformatics pipelines, including structural and functional annotations (e.g., transposable elements and repetitive sequences). This paper also includes information on how to build a wide community for a genome project, the importance of data management, and how to make the data and results Findable, Accessible, Interoperable, and Reusable (FAIR) by submitting them to a public repository and sharing them with the research community.

Chromosome-level genome assembly for the largemouth bass Micropterus salmoides provides insights into adaptation to fresh and brackish water

Largemouth bass (LMB; Micropterus salmoides) has been an economically important fish in North America, Europe, and China. This study obtained a chromosome-level genome assembly of LMB using PacBio and Hi-C sequencing. The final assembled genome is 964 Mb, with contig N50 and scaffold N50 values of 1.23 Mb and 36.48 Mb, respectively. Combining with RNA sequencing data, we annotated a total of 23,701 genes. Chromosomal assembly and syntenic analysis proved that, unlike most Perciformes with the popular haploid chromosome number of 24, LMB has only 23 chromosomes (Chr), among which the Chr1 seems to be resulted from a chromosomal fusion event. LMB is phylogenetically closely related to European seabass and spotted seabass, diverging 64.1 million years ago (mya) from the two seabass species. Eight gene families comprising 294 genes associated with ionic regulation were identified through positive selection, transcriptome and genome comparisons. These genes involved in iron facilitated diffusion (such as claudin, aquaporins, sodium channel protein and so on) and others related to ion active transport (such as sodium/potassium-transporting ATPase and sodium/calcium exchanger). The claudin gene family, which is critical for regulating cell tight junctions and osmotic homeostasis, showed a significant expansion in LMB with 27 family members and 68 copies for salinity adaptation. In summary, we reported the first high-quality LMB genome, and provided insights into the molecular mechanisms of LMB adaptation to fresh and brackish water. The chromosome-level LMB genome will also be a valuable genomic resource for in-depth biological and evolutionary studies, germplasm conservation and genetic breeding of LMB.

Transcriptome characterization of BPG axis and expression profiles of ovarian steroidogenesis-related genes in the Japanese sardine

Background

The clupeoid fishes are ecologically and commercially important fish species worldwide that exhibit a high level of population fluctuation, accompanied by alteration of reproductive traits. However, knowledge about their reproductive physiology in order to understand mechanisms underlying such population dynamics is limited. The endocrine system along with the brain–pituitary–gonadal (BPG) axis is critical for regulating reproduction. The aims of this study were to provide transcript data and genes related to the BPG axis, and to characterize the expression profiles of ovarian steroidogenesis-related genes in the Japanese sardine (Sardinops melanostictus, Clupeidae).

Results

RNA sequencing was performed using the sardine brain, pituitary, and gonad in both sexes. A total of 290,119 contigs were obtained and 115,173 non-redundant ORFs were annotated. The genes differentially expressed between ovary and testis were strongly associated with GO terms related to gamete production. The tissue-specific profile of the abundance of transcripts was characterized for the major regulators in the BPG axis, such as gonadotropin-releasing hormone, gonadotropin, and steroidogenic enzyme. By comparing between ovary and testis, out of eight different 17β-hydroxysteroid dehydrogenase (Hsd17b) genes identified, higher hsd17b7 expression was found in testis, whereas higher expression of hsd17b8, hsd17b10, hsd17b12a, and hsd17b12b was found in ovary. The cDNAs encoding key endocrine factors in the ovarian steroidogenic pathway were cloned, sequenced, and quantitatively assayed. In the pituitary, follicle-stimulating hormone beta peaked during vitellogenesis, while luteinizing hormone beta peaked at the completion of vitellogenesis. In the ovary, follicle-stimulating hormone receptor and luteinizing hormone receptor were upregulated from mid- to late phase of vitellogenesis. Furthermore, three steroidogenic enzyme genes (cyp11a1, cyp17a1, and cyp19a1a) gradually increased their expression during ovarian development, accompanying a rise in serum estradiol-17β, while 3β-hydroxysteroid dehydrogenase and steroidogenic acute regulatory protein did not change significantly.

Conclusions

This is the first report of deep RNA sequencing analysis of Japanese sardine, in which many key genes involved in the BPG axis were identified. Expression profiles of ovarian steroidogenesis-related genes provide a molecular basis of the physiological processes underlying ovarian development in the sardine. Our study will be a valuable resource for clarifying the molecular biology of clupeoid fishes.

Reconstruction of the birth of a male sex chromosome present in Atlantic herring

Understanding the evolution of sex determination mechanisms and sex chromosomes is of fundamental importance in biology. Here we have reconstructed the evolution of the sex-determining region in the Atlantic herring. The region is small and contains only three intact genes. The candidate sex-determining factor BMPR1BBY is an evolutionary innovation in the herring lineage. It encodes a truncated form of a BMP type I receptor, which originated by gene duplication and underwent rapid protein evolution. The receptor has maintained its kinase activity and has the potential to induce development of testis. The other two genes in the sex-determining region, CATSPERG and CATSPER3Y, are male beneficial genes because they encode proteins predicted to be essential for sperm to fertilize the egg.

The mechanisms underlying sex determination are astonishingly plastic. Particularly the triggers for the molecular machinery, which recalls either the male or female developmental program, are highly variable and have evolved independently and repeatedly. Fish show a huge variety of sex determination systems, including both genetic and environmental triggers. The advent of sex chromosomes is assumed to stabilize genetic sex determination. However, because sex chromosomes are notoriously cluttered with repetitive DNA and pseudogenes, the study of their evolution is hampered. Here we reconstruct the birth of a Y chromosome present in the Atlantic herring. The region is tiny (230 kb) and contains only three intact genes. The candidate male-determining gene BMPR1BBY encodes a truncated form of a BMP1B receptor, which originated by gene duplication and translocation and underwent rapid protein evolution. BMPR1BBY phosphorylates SMADs in the absence of ligand and thus has the potential to induce testis formation. The Y region also contains two genes encoding subunits of the sperm-specific Ca²⁺ channel CatSper required for male fertility. The herring Y chromosome conforms with a characteristic feature of many sex chromosomes, namely, suppressed recombination between a sex-determining factor and genes that are beneficial for the given sex. However, the herring Y differs from other sex chromosomes in that suppression of recombination is restricted to an ∼500-kb region harboring the male-specific and sex-associated regions. As a consequence, any degeneration on the herring Y chromosome is restricted to those genes located in the small region affected by suppressed recombination.

A Nanopore Based Chromosome-Level Assembly Representing Atlantic Cod from the Celtic Sea

Currently available genome assemblies for Atlantic cod (Gadus morhua) have been constructed from fish belonging to the Northeast Arctic Cod (NEAC) population; a migratory population feeding in the Barents Sea. These assemblies have been crucial for the development of genetic markers which have been used to study population differentiation and adaptive evolution in Atlantic cod, pinpointing four discrete islands of genomic divergence located on linkage groups 1, 2, 7 and 12. In this paper, we present a high-quality reference genome from a male Atlantic cod representing a southern population inhabiting the Celtic sea. The genome assembly (gadMor_Celtic) was produced from long-read nanopore data and has a combined contig length of 686 Mb with an N50 of 10 Mb. Integrating contigs with genetic linkage mapping information enabled us to construct 23 chromosome sequences which mapped with high confidence to the latest NEAC population assembly (gadMor3) and allowed us to characterize, to an extent not previously reported large chromosomal inversions on linkage groups 1, 2, 7 and 12. In most cases, inversion breakpoints could be located within single nanopore contigs. Our results suggest the presence of inversions in Celtic cod on linkage groups 6, 11 and 21, although these remain to be confirmed. Further, we identified a specific repetitive element that is relatively enriched at predicted centromeric regions. Our gadMor_Celtic assembly provides a resource representing a 'southern' cod population which is complementary to the existing 'northern' population based genome assemblies and represents the first step toward developing pan-genomic resources for Atlantic cod.