Preferential subfunctionalization of slow-evolving genes after allopolyploidization inXenopus laevis

Abstract
As paleopolyploid genomes evolve, the expression profiles of retained gene pairs are expected to diverge. To examine this divergence process on a large scale in a vertebrate system, we compareXenopus laevis, which has retained ≈40% of loci in duplicate after a recent whole-genome duplication (WGD), with its unduplicated relativeSilurana (Xenopus)tropicalis. This comparison of ingroup pairs to an outgroup allows the direction of change in expression profiles to be inferred for a set of 1,300X. laevispairs, relative to their single orthologs inS. tropicalis, across 11 tissues. We identify 68 pairs in whichX. laevisis inferred to have undergone a significant reduction of expression in at least two tissues since WGD. Of these pairs, one-third show evidence of subfunctionalization, with decreases in the expression levels of different gene copies in two different tissues. Surprisingly, we find that genes with slow rates of evolution are particularly prone to subfunctionalization, even when the tendency for highly expressed genes to evolve slowly is controlled for. We interpret this result to be an effect of allopolyploidization. We then compare the outcomes of this WGD with an independent one that happened in the teleost fish lineage. We find that if a gene pair was retained in duplicate inX. laevis, the orthologous pair is more likely to have been retained in duplicate in zebrafish, suggesting that similar factors, among them subfunctionalization, determined which gene pairs survived in duplicate after the two WGDs.