The functional divergence of duplicate genes (ohnologues) retained from whole genome duplication (WGD) is thought to promote evolutionary diversification. However, species radiation and phenotypic diversification are often temporally separated from WGD. Salmonid fish, whose ancestor underwent WGD by autotetraploidization ~95 million years ago, fit such a ‘time-lag’ model of post-WGD radiation, which occurred alongside a major delay in the rediploidization process. Here we propose a model, ‘lineage-specific ohnologue resolution’ (LORe), to address the consequences of delayed rediploidization. Under LORe, speciation precedes rediploidization, allowing independent ohnologue divergence in sister lineages sharing an ancestral WGD event. Using cross-species sequence capture, phylogenomics and genome-wide analyses of ohnologue expression divergence, we demonstrate the major impact of LORe on salmonid evolution. One-quarter of each salmonid genome, harbouring at least 4550 ohnologues, has evolved under LORe, with rediploidization and functional divergence occurring on multiple independent occasions >50 million years post-WGD. We demonstrate the existence and regulatory divergence of many LORe ohnologues with functions in lineage-specific physiological adaptations that potentially facilitated salmonid species radiation. We show that LORe ohnologues are enriched for different functions than ‘older’ ohnologues that began diverging in the salmonid ancestor. LORe has unappreciated significance as a nested component of post-WGD divergence that impacts the functional properties of genes, whilst providing ohnologues available solely for lineage-specific adaptation. Under LORe, which is predicted following many WGD events, the functional outcomes of WGD need not appear ‘explosively’, but can arise gradually over tens of millions of years, promoting lineage-specific diversification regimes under prevailing ecological pressures.

Background: Whole genome duplication (WGD) events have played a major role in eukaryotic genome evolution, but the consequence of these extreme events in adaptive genome evolution is still not well understood. To address this knowledge gap, we used a comparative phylogenetic model and transcriptomic data from seven species to infer selection on gene expression in duplicated genes (ohnologs) following the salmonid WGD 80–100 million years ago.Results: We find rare cases of tissue-specific expression evolution but pervasive expression evolution affecting many tissues, reflecting strong selection on maintenance of genome stability following genome doubling. Ohnolog expression levels have evolved mostly asymmetrically, by diverting one ohnolog copy down a path towards lower expression and possible pseudogenization. Loss of expression in one ohnolog is significantly associated with transposable element insertions in promoters and likely driven by selection on gene dosage including selection on stoichiometric balance. We also find symmetric expression shifts, and these are associated with genes under strong evolutionary constraints such as ribosome subunit genes. This possibly reflects selection operating to achieve a gene dose reduction while avoiding accumulation of “toxic mutations”. Mechanistically, ohnolog regulatory divergence is dictated by the number of bound transcription factors in promoters, with transposable elements being one likely source of novel binding sites driving tissue-specific gains in expression.Conclusions: Our results imply pervasive adaptive expression evolution following WGD to overcome the immediate challenges posed by genome doubling and to exploit the long-term genetic opportunities for novel phenotype evolution.

The European flat oyster (Ostrea edulis L.) is a bivalve naturally distributed across Europe, which was an integral part of human diets for centuries, until anthropogenic activities and disease outbreaks severely reduced wild populations. Despite a growing interest in genetic applications to support population management and aquaculture, a reference genome for this species is lacking to date. Here, we report a chromosome‐level assembly and annotation for the European Flat oyster genome, generated using Oxford Nanopore, Illumina, Dovetail OmniC™ proximity ligation and RNA sequencing. A contig assembly (N50: 2.38 Mb) was scaffolded into the expected karyotype of 10 pseudochromosomes. The final assembly is 935.13 Mb, with a scaffold‐N50 of 95.56 Mb, with a predicted repeat landscape dominated by unclassified elements specific to O. edulis. The assembly was verified for accuracy and completeness using multiple approaches, including a novel linkage map built with ddRAD‐Seq technology, comprising 4016 SNPs from four full‐sib families (eight parents and 163 F1 offspring). Annotation of the genome integrating multitissue transcriptome data, comparative protein evidence and ab‐initio gene prediction identified 35,699 protein‐coding genes. Chromosome‐level synteny was demonstrated against multiple high‐quality bivalve genome assemblies, including an O. edulis genome generated independently for a French O. edulis individual. Comparative genomics was used to characterize gene family expansions during Ostrea evolution that potentially facilitated adaptation. This new reference genome for European flat oyster will enable high‐resolution genomics in support of conservation and aquaculture initiatives, and improves our understanding of bivalve genome evolution.

European flat oyster (Ostrea edulis) is an ecologically and economically important marine bivalve, that has been severely affected by the intracellular parasite Bonamia ostreae. In this study, a flat oyster SNP array (~14,000 SNPs) was used to validate previously reported outlier loci for divergent selection associated with B. ostreae exposure in the Northeast Atlantic Area. A total of 134 wild and hatchery individuals from the North Sea, collected in naïve (NV) and long‐term affected (LTA) areas, were analysed. Genetic diversity and differentiation were related to the sampling origin (wild vs. hatchery) when using neutral markers, and to bonamiosis status (NV vs. LTA) when using outlier loci for divergent selection. Two genetic clusters appeared intermingled in all sampling locations when using outlier loci, and their frequency was associated with their bonamiosis status. When both clusters were compared, outlier data sets showed high genetic divergence (F_ST > 0.25) unlike neutral loci (F_ST not ≠ 0). Moreover, the cluster associated with LTA samples showed much higher genetic diversity and significant heterozygote excess with outlier loci, but not with neutral data. Most outliers mapped on chromosome 8 (OE‐C8) of the flat oyster genome, supporting a main genomic region underlying resilience to bonamiosis. Furthermore, differentially expressed genes previously reported between NV and LTA strains showed higher mapping density on OE‐C8. A range of relevant immune functions were specifically enriched among genes annotated on OE‐C8, providing hypotheses for resilience mechanisms to an intracellular parasite. The results suggest that marker‐assisted selection could be applied to breed resilient strains of O. edulis to bonamiosis, if lower parasite load and/or higher viability of the LTA genetic cluster following B. ostreae infection is demonstrated.

Standardised analysis pipelines are an important part of FAIR bioinformatics research. Over the last decade, there has been a notable shift from point-and-click pipeline solutions such as Galaxy towards command-line solutions such as Nextflow and Snakemake. We report on recent developments in the nf-core and Nextflow frameworks that have led to widespread adoption across many scientific communities. We describe how adopting nf-core standards enables faster development, improved interoperability, and collaboration with the >8,000 members of the nf-core community. The recent development of Nextflow Domain-Specific Language 2 (DSL2) allows pipeline components to be shared and combined across projects. The nf-core community has harnessed this with a library of modules and subworkflows that can be integrated into any Nextflow pipeline, enabling research communities to progressively transition to nf-core best practices. We present a case study of nf-core adoption by six European research consortia, grouped under the EuroFAANG umbrella and dedicated to farmed animal genomics. We believe that the process outlined in this report can inspire many large consortia to seek harmonisation of their data analysis procedures.