Exploring the Effect of Whole-Genome Duplication on Salmonid LincRNA Repertoire
Long intergenic non-coding RNAs (lincRNAs) are key epigenetic regulators of genome function, yet their evolutionary dynamics following whole-genome duplication (WGD) events remain poorly understood. Salmonids, which underwent a lineage-specific autotetraploidization (salmonid-specific WGD, ~88–100 million years ago), provide an excellent model to investigate the retention, divergence, and functional potential of recently duplicated non-coding elements. LincRNA repertoires were compared across five genome-annotated salmonids (Oncorhynchus tshawytscha, O. kisutch, O. mykiss, Salmo salar, and S. trutta) and their closest non-duplicated relative, northern pike (Esox lucius). LincRNAs represented ~5–7% of annotated genes in all salmonids except S. salar (18%). Sequence conservation was low relative to coding genes, with only 11–68 highly similar (e-value < 1 × 10−30; similarity > 70% and alignments > 100 nucleotides) putative orthologues shared between salmonids and northern pike, and 161–338 among salmonids alone. Synteny conservation was modest in lincRNAs, with lower conservation in putative orthologues (8–16%) compared to putative ohnologues (8–33%). Secondary structure conservation was associated with sequence similarity (ρ = −0.45; p = 2.2 × 10−16), and the association was stronger among WGD ohnologues than orthologues. In S. salar and O. mykiss, lincRNA putative ohnologues showed weaker expression correlations than coding genes, suggesting widespread regulatory divergence, possibly through neo- and subfunctionalisation. Conserved salmonid lincRNAs showed enriched predicted interactions with miRNAs involved in tumour suppression, brain, bone, and muscle development (e.g., miR-455, miR-365, miR124, miR-133a, miR-140, and miR-9), a finding supported by limited transcriptomic data. Although salmonid WGD expanded lincRNA repertoires, lincRNAs have undergone rapid sequence and transcriptional divergence, with limited conservation across species based on sequence similarity, chromosomal position, synteny, and secondary structure. A subset of conserved lincRNAs retains structural features and regulatory signatures consistent with roles as miRNA sponges in brain, skeletal, and muscle development and tumour suppression, potentially acting within conserved regulatory networks. These findings provide new insights into lincRNA evolution following genome duplication and highlight the need for experimental validation of their regulatory functions.