XB-ART-36320BMC Biol. July 25, 2007; 5 31.
Accelerated gene evolution and subfunctionalization in the pseudotetraploid frog Xenopus laevis.
Ancient whole genome duplications have been implicated in the vertebrate and teleost radiations, and in the emergence of diverse angiosperm lineages, but the evolutionary response to such a perturbation is still poorly understood. The African clawed frog Xenopus laevis experienced a relatively recent tetraploidization ~40 million years ago. Analysis of the considerable amount of EST sequence available for this species together with the genome sequence of the related diploid Xenopus tropicalis provides a unique opportunity to study the genomic response to whole genome duplication. We identified 2218 gene triplets in which a single gene in X. tropicalis corresponds to precisely two co-orthologous genes in X. laevis--the largest such collection published from any duplication event in animals. Analysis of these triplets reveals accelerated evolution or relaxation of constraint in the peptides of the X. laevis pairs compared with the orthologous sequences in X. tropicalis and other vertebrates. In contrast, single-copy X. laevis genes do not show this acceleration. Duplicated genes can differ substantially in expression levels and patterns. We find no significant difference in gene content in the duplicated set, versus the single-copy set based on molecular and biological function ontologies. These results support a scenario in which duplicate genes are retained through a process of subfunctionalization and/or relaxation of constraint on both copies of an ancestral gene.
PubMed ID: 17651506
PMC ID: PMC1949811
Article link: BMC Biol.
Grant support: GM66684 NIGMS NIH HHS , K08-HD42550 NICHD NIH HHS , K08 HD042550 NICHD NIH HHS , R01 GM066684 NIGMS NIH HHS , GM66684 NIGMS NIH HHS , K08-HD42550 NICHD NIH HHS
Genes referenced: clock foxa1 skp1 sri
Article Images: [+] show captions
|Figure 2. Symmetric evolution of paralogs. Scatter plot of relative evolution between X. tropicalis peptides and their co-orthologous sequences in X. laevis. A total of 578 gene triples with 16 or more highly-conserved positions are shown (see text for details). L1 and L2 refer to co-orthologous genes 1 and 2 in X. laevis. The diagonal line represents a null model assuming symmetric evolution of L1 and L2. Black boxes are L1�L2 pairs inconsistent with this model at P < 0.01.|
|Figure 3. Normalized peptide to nucleotide evolutionary rates show an accelerated divergence of duplicated X. laevis peptides. The chart shows the ratio of peptide evolution (P-distance) to synonymous transversion rates (4 DTv), normalized by the human-mouse P-distance/4 DTV value of 0.242 � 0.004, for three sets of multiple alignments corresponding to genes found in single copy in each of human, mouse, rat, and X. tropicalis, and two copies in X. laevis (sextuplets); pentuplets obtained by randomly selecting one X. laevis paralog from each sextuplet (5 A), and pentuplets in which only a single X. laevis sequence is known (5B).|
|Figure 4. Expression of specific X. laevis paralogs and their X. tropicalis ortholog. Panels depict the expression of skp1a (a-c), isocitrate dehydrogenase (isoD) (d-f), foxA1 (g-i), sorcin (j-l). X. laevis paralogs were arbitrarily assigned as a (a,d,g,j) or b (b,e,h,k) and are compared to the X. tropicalis ortholog (c, f, i, l). All views are lateral with anterior to the left. Embryos (a-f, j-l) are at stages 31 while embryos (g-i) are at stage 37�38. The arrowhead in (b) indicates kidney expression of skp1a in X. laevis paralog b that is not seen in the a paralog. Insets in (d) and (e) magnify somite expression revealing the differential expression between X. laevis paralogs ((d) with narrow expression, (e) with broad expression). The arrow in (g) highlights posterior expression of foxA1 seen in paralog a but absent in paralog b. The arrow in (k) indicates weak lateral expression of sorcin in X. laevis paralog b that is not seen in paralog a. X. tropicalis embryos are shown at a higher magnification than X. laevis embryos, reflecting their smaller size.|
|Figure 5. 4DS distances identify genome duplication event. Histograms of the 4DS distances for the 9905 mutual highest scoring L-T pairs (blue line) as well as for the 3358 unambigous L-L pairs (red bars). The L-L pairs with 0.05 < 4DS < 0.25, peaking around 0.16, are selected as originating from the genome duplication event.|