XB-ART-52612Nature January 1, 2016; 538 (7625): 336-343.
Genome evolution in the allotetraploid frog Xenopus laevis.
To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of ''fossil'' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.
PubMed ID: 27762356
PMC ID: PMC5313049
Article link: Nature
Grant support: HD065705 , GM086321 , R01 HD069344 NICHD NIH HHS , P41 HD064556 NICHD NIH HHS , HD080708 , R35 GM118183 NIGMS NIH HHS , R01 GM086627 NIGMS NIH HHS , R21 GM119021 NIGMS NIH HHS , R01 HL117164 NHLBI NIH HHS , T32 HG000047 NHGRI NIH HHS , R01 GM042341 NIGMS NIH HHS , R01 GM086321 NIGMS NIH HHS , R01 GM104853 NIGMS NIH HHS , S10 RR027303 NCRR NIH HHS , R01 OD010549 NIH HHS , R21 HD084072 NICHD NIH HHS , R01 HD080708 NICHD NIH HHS
Genes referenced: six6 actr3 ap3s1 aqp7 arrb1 atp6ap1.2 avd bix1.1 bix1.2 bix1.3 brd2 capza1 ccna1 ccna2 ccnb1.2 ccnb3 ccnd1 ccnd2 ccndx ccne1 ccne2 ccnh ccny cdk1 cdk2 cdk4 cdk6 cdk7 cers1 chrd.1 coil cox7c crb2 crb3 crebbp ctnnb1 cttnbp2 dnase1l1 eda evx1 fam50a foxa1 fubp1 gdf1 gdf3 gdi1 gnao1 gsk3a hhat hhip hoxb2 hoxb4 htt inhba kcnd3 kcnh1 klf2 lats2 lifr lrp5 lypd1 map4k4 mob1a mob1b mst1 myh6 nf2 nodal1 nodal2 nodal3.1 nodal5 nodal6 numbl ptch1 rab11fip4 rcor3 rhoc rnf4 rpl10 rpl13a rps11 sav1 scpep1 scrib sf3b4 shh skp1 slco2b1 smo spcs3 stau1 suz12 sv2a tacc1 tap2 tapbp taz tcf3 tcf4 tcf7 tdrd6 tgfb1 tmem161b trim25 ubxn1 utp25 utp6 vdac1 vps41 wnt1 wnt11b wnt2b yap1
Article Images: [+] show captions
|Figure 1. Chromosome evolution in Xenopusa. Comparative cytogenetic map of XLA and XTR chromosomes. Magenta lines show relationships of chromosomal locations of 198 homeologous gene pairs between XLAL and S chromosomes, identified by FISH mapping using BAC clones (Supplemental Table 1, see Supplementary Note 3.1). Blue lines show relationships of chromosomal locations of orthologous genes between XTR chromosomes and (i) both XLA L and S chromosomes (solid line) (lines between XLAL and S are omitted), (ii) only XLAL (dashed), or (iii) only XLAS (dotted), which were taken from our previous studies14,15. Light blue lines indicate positional relationships of actr3 and lypd1 on XTR9q and rpl13a and rps11 on XTR10q with those on XLA9_10LS chromosomes (see Supplementary Note 6.2). Double-headed arrows on the right of XLAS chromosomes indicate the chromosomal regions in which inversions occurred. Ideograms of XTR and XLA chromosomes were taken from our previous reports15,16.b. Distribution of homeologous genes (purple), singletons (grey), and subgenome-specific repeats across XLA1L (top) and XLA1S (bottom). Xl-TpL_Harb is red, Xl-TpS_Harb is blue, and Xl-TpS_Mar is green. Purple lines mark homeologous genes present in both L and S chromosomes, the black line marks the approximate centromere location on each chromosome. The homeologous gene pairs, from left to right: rnf4, spcs3, intsl2, foxa1, sds, ap3s1, lifr, aqp7. Each bin is 3 MB in size, with 0.5 MB overlap with the previous bin.c. Chromosomal localization of the Xl-TpS_Mar sequence. Hybridization signals were only observed on the S chromosomes. Scale bar shows 10 μm.|
|Figure 2. Molecular evolution and allotetraploidya. (top) The distribution of pseudogene ages, as described in Supplemental Note 9. (bottom) Phylogenetic tree illustrating the different epochs in Xenopus, with times based on protein-coding gene phylogeny of pipids, including Xenopus, Pipa carvalhoi, Hymenochirus boettgeri, and Rana pipiens (only Xenopus depicted). We date the speciation of X. tropicalis and the X. laevis ancestor at 48 Mya, the L and S polyploid progenitors at 34 Mya, and the divergence of the polyploid Xenopus radiation at 17 Mya. Using these times as calibration points, we estimate bursts of transposon activity at 18 Mya (mariner, blue star) and 33–34 Mya (harbinger, red star). The purple star is the time of hybridization, around 17–18 Mya.b. Phylogenetic tree based on protein-coding genes of tetrapods, rooted by elephant shark (not shown). Alignments were done by MACSE, the maximum-likelihood tree was built by PhyML. Branch length scale shown on bottom. The difference in branch length between Xenopus laevis-L and Xenopus laevis-S is similar to that seen between mouse and rat. Both subgenomes of X. laevis have longer branch lengths than X. tropicalis.|
|Figure 3. Structural response to allotetraploidya. Distributions of consecutive retentions (left) and deletions (right) in the L (red) and S (blue) subgenomes. The distributions were fit using the equation y= a*(ebx) + c*(edx). The y-axis is shown on a log scale. Significant differences were seen between L and S subgenomes in both distributions (Student’s t-test, retention p=3.6E-22, deletion p=4.5E-84).b. Evolutionary conservation of the Xenopus MHC and differential MHC silencing on the two X. laevis subgenomes. Selected gene names shown above. The ‘Adaptive MHC’ encodes tightly-linked essential genes involved in antigen presentation to T cells; this group of genes is the primordial linkage group and has been preserved in most non-mammalian vertebrates, including Xenopus. Differential gene silencing is particularly pronounced as four genes around the class I gene are functional on S chromosome but absent (dma, dmb) or pseudogenes (ring3, lmp2) on L chromosome. The gene map is not to scale; pseudogenes (p) are noted as indicated. HSA, XLA, GGA: human, Xenopus, and chicken MHC, respectively. Refer to the Supplemental Table 8 for a more detailed MHC map.c. Hox gene clusters. X. laevis retains eight Hox clusters, consisting of pairs of HoxA, B, C and D clusters, on L and S chromosomes. even-skipped genes (evx1 or evx2) are positioned flanking Hoxa and Hoxd clusters. hox genes are classified into four, labial, proboscipedia, central, and posterior groups. Note that hoxb2.L (black) is a pseudogene.d. Syntenies around the mix gene family. Abbreviations for species and chromosome numbers: human (H. sapiens; HSA1), chicken (G. gallus; GGA3), X. tropicalis (XTR5), X. laevis (XLA5L and XLA5S), zebrafish (D. rerio; DRE20). Each Xenopus (sub)genome experienced its own independent expansion of the family (see Extended Data Fig. 5 for details).|
|Figure 4. Retention and functional differentiationa. (left) Comparison of L and S gene loss by KEGG categories. X-axis is fraction loss of L genes, and the y-axis is fraction loss of S genes. Blue line is expected L/S loss based on genome-wide average (56.4%). Red points are functional categories that show a high degree of loss (χ2 test p<0.01). Magenta points are functional categories that show a high degree of retention (p<0.01). (right) Similar scatterplot for tissue WGCNA categories. See Supplemental Note 10.1 for a more detailed discussion.b. Boxplot of log10(Ltpm/Stpm) for homeologous gene pairs, zoomed in to show medians. Ovary and maternally-controlled developmental time points are on the left (light blue and dark blue bars respectively), zygotically-controlled developmental time points and adult tissues are on the right (red and green bars respectively). The red line shows the equal ratio log10(1). On average maternal data sets express the L gene of a homeologous pair 12% higher than as S (median = 0%), while the zygotic tissues and time points express the L gene of a homeologous pair 25% higher than S (median = 1.8%). The difference between the mean and medians is explained by many genes with large differences between homeologs, illustrated by the full distribution in Extended Data Fig. 8c. Here, to illustrate the difference in median of zygotic expression, we zoom in on the center of the boxplot.c. (left) Developmental expression plot and (right) epigenetic landscape surrounding hoxb4. L expression is red, S expression is blue, tropicalis expression is shown in black. The right panel shows the genomic profiles of H3K4me3 (green) and p300 (yellow) ChIP-seq tracks, as well as DNA methylation levels determined by whole-genome bisulfite sequencing (grey). The gene annotation track shows the hoxb4 gene on L (top) and S. The conservation between the L and S genomic sequence is shown in grey between the gene annotation tracks.d. (left) Developmental expression plot and (right) epigenetic landscape surrounding numbl. L expression is red, S expression is blue, tropicalis expression is shown in black. The small amount of expression seen in maternal numbl and numbl.L is consistent between replicates. In addition to the tracks described for c), the right panel shows RNA Polymerase II (RNAPII; purple) and H3K36me3 (blue) ChIP-seq profiles.e. Representative embryos with GFP expression driven by either six6.L-CNE or six6.S-CNE linked to a basal promoter-GFP cassette (six6.L-CNE:GFP and six6.S-CNE:GFP, respectively). GFP expression was detected by in situ hybridization. Semi-quantitative image analysis revealed a statistically significant difference in their average expression level (p < 0.01); the expression driven by six6.S-CNE (n = 27) was 0.6-fold weaker than that by six6.L-CNE in the eye region (n = 32). Given eye-specific patterns of their endogenous expression, the six6 genes likely have additional silencers for restricting enhancer activity of the CNEs in the eye.|
|a–e, Scenarios for allotetraploid formation from distinct ancestral diploid species A and B. Horizontal single lines indicate normal gametes, horizontal double lines indicate unreduced gametes; black square represents fertilization; vertical double lines indicate spontaneous (somatic) genome doubling. a, (i) Fusion of unreduced gametes from species A and B. (ii) Interspecific hybridization followed by spontaneous doubling. (iii) Fusion of unreduced gametes produced by interspecific hybrids. (iv) Interspecific hybrids produce unreduced gametes, which fuse with normal gametes from species A. The resulting triploid again produces unreduced gametes, which fuse with normal gametes from species B. (v) Unreduced gamete from species A fuses with normal gamete from species B. The resulting AAB triploid produces unreduced gametes that are fertilized by normal gametes species B. See Supplementary Note 1.1 for a more detailed discussion. b, History of the J strain. See Supplementary Note 2.1 for details. The years of events and generation numbers (such as frog transfer to another institute, establishment of homozygosity, construction of materials) are indicated in the scheme. Generation numbers are estimates due to loss of old breeding records. c, The nucleotide distance of orthologues (green), homoeologues (red) and alleles (blue) is discussed in Supplementary Note 8.7. The distances are shown on a log scale to differentiate between the distributions. d, Frequency histogram showing the number of 51-mers with specified count in the shotgun dataset. The prominent peak implies that each genomic locus is sampled 29× in 51-mers. Note the absence of a feature at twice this depth, indicating that homoeologous features with high identity are rare. e, Cumulative proportion of 51-mers as a function of relative depth (that is, depth/29). Relative depth provides an estimate of genomic copy number. The rapid rise at relative depth 1 implies that 70–75% of the X. laevis genome is a single copy with respect to 51-mers. The remainder of the genome is primarily concentrated in repetitive sequences with copy number > 100. Note logarithmic scale. f, The contact map of 85,260 TCC read pairs for JGIv72.000090484.chr4S. Read pairs were binned at 10-kb intervals. For each read pair, the forward and reverse reads map with a map quality score of at least 20. g, The contact map of 85,260 Chicago read pairs for JGIv72.000090484.chr4S, a 3.1-Mb scaffold in the XENLA_JGI_v72 assembly. h, The insert distribution of TCC and Chicago read pairs that map to the same scaffold of XENLA_JGI_v72 with a map quality score of at least 20. The x axis is the read pair separation distance. The y axis is the counts for that bin divided by the total number of reads. The bins are 1 kb.|
|a, Structure of the sex chromosome of X. laevis (XLA2L) and comparison with XLA2S and XTR2. The W version of XLA2L harbours a W-specific sequence containing the female sex-determining gene dmw (red) while Z has a different Z-specific sequence (blue). Pentagon arrows and black triangles indicate genes and olfactory receptor genes, respectively. Their tips correspond to their 3′-ends. b, Alignment of the q-terminal regions of XTR9 and 10 with corresponding regions of XLA9_10L and XLA9_10S. Genes near the q-terminal regions of XTR 9 and XTR10 were missing in the X. tropicalis genome assembly v9, but rps11, rpl13a, lypd1 and actr3 were expected to be located there based on the synteny with human chromosomes, and then verified by cDNA FISH (upper panels). Small triangles on XLA9_10L and S indicate the distribution of gene models showing both identity and coverage greater than 30%, against the human and chicken peptide sequences from Ensembl, in the region ±2 Mb from the prospective 9/10 junction. HSA, human chromosome; GGA, chicken chromosome. The magnified view represents syntenic genes to scale with colours corresponding to human genes. c, The orders of orthologous genes across XTR9, XTR10, XLA9_10L and XLA9_10S. Green arrowheads: positions of centromeres in XTR9 and 10 predicted by examination of the cytogenetic chromosome length ratio of p versus q arms15. Blue arrowheads: positions of centromere repeats, frog centromeric repeat-1 (ref. 55), in XLA9_10L and S. Magenta and yellow ellipses, chromosomal locations of snrpn (magenta) and stau1 (yellow) from X. tropicalis v9 and X. laevis v9.1 assemblies. Red ellipses, chromosomal locations of four genes, rps11, rpl13a, lypd1 and actr3. XTR9 is inverted to facilitate comparison. Blue bidirectional arrows indicate the homologous regions where pericentric inversions may have occurred on proto-chromosomes (see Extended Data Fig. 2d). d, Schematic representation for the two hypothetical processes of chromosomal rearrangements (fusion and inversion) that occurred between the hypothetical proto-XTR9 and 10 to produce proto-XLA9_10, and eventually XLA9_10L and S. The process of chromosome rearrangements is explained parsimoniously in two different ways (left and right panels), starting from proto-XTR9 and 10. Actual and hypothetical ancestral chromosomal locations of snrpn and stau1 are shown by magenta and yellow circles, respectively. Note that the chromosomal locations of these genes on the proto-XTR10 differ between the two models. Chromosome segments homologous to XTR9 and XTR10 are shown in red and blue, respectively. XTR9 is inverted to facilitate comparison. Bidirectional arrows indicate the regions where pericentric inversions may have occurred. Black arrows indicate the direction of chromosomal evolution.|
|EDF 3. TRANSPOSONS Density of the subgenome specific transposons on each chromosome (coverage length of transposable element [bp]/chromosome length [Mbp]). The coverage lengths of transposons were calculated from the results of BLASTN search (E-value cutoff 1E-5) using the consensus sequences as queries.Jukes-Cantor distances across non-CpG sites, corrected as in Supplemental Note 7.5. Distances between X. tropicalis and X. laevis transposons consensus sequences are shown. The X. laevis-specific transposon differences are each individual transposon sequence against the consensus sequence for that subfamily.Phylogenetic tree of Xl-TpS_Mar transposon expansions in the X. laevis genome, built using Jukes-Cantor corrected distances (Supplemental Note 7.5). Sub-clusters with enough members to determine accurate timings are highlighted. The scale bar represents the corrected Jukes-Cantor distance.|
|EDF 4. PHYLOGENY Phylogenetic tree of pan-vertebrate conserved non-coding elements (pvCNEs), rooted by elephant shark. Alignments were done by MUSCLE, and the maximum-likelihood tree was built by PhyML. Branch length scale shown on bottom. The difference in branch lengths of tetrapods follows the same topology as the protein-coding tree (Fig. 2b).Complete phylogenetic tree from Fig. 2a, with divergence times computed by r8s.Distribution of Ks and Ka on specific subgenomes during the time between L and S speciation, before X. laevis and X. borealis speciation. We find accelerated mutations rates between T2 and T3 in Ks and Ka (p=1.4e-5 (left), 8.6e-3 (right)).Distribution of Ks and Ka on specific subgenomes during the time after laevis and borealis speciation. We do not find significantly accelerated substitution rates. (p= 0.10 (left) and 0.03 (right)).Table showing the number of homeologs and singletons identified as homeologs from the ancient vertebrate duplication (or ohnologs as they are historically called)63. 79.9% of ohnologs retain both copies in X. laevis today, significantly more than the 54.3% of the rest of the genome after excluding ohnologs (χ2 test p-value= 4.44E-69).Table showing the branch lengths of bootstrapped maximum likelihood trees described in Supplemental Note 12.5. The columns refer to the X. tropicalis (XTR), L chromosome of X. laevis (XLA.L), S chromosome of X. laevis (XLA.S), and XLA.L/XLA.S branch lengths respectively. The first row is triplets where all genes show expression, the second row is triplets where L is a thanagene, and the third row is triplets where S is a thanagene. The L branch length is significantly smaller when all genes are expressed, or when S is a thanagene (Wilcoxon p-value=1.7E-216 and 6.4E-212 respectively). The S branch length is smaller when L is a thanagene (p=2.4E-223). The ratio of branch lengths (L/S) is significantly different for either L or S thanagene datasets compared to when all genes are expressed (p=3.55E-214 and 7.48E-220 respectively). The ratio is different between the two thanagene datasets as well (p=1.79E-217).|
|EDF 5. STRUCTURAL EVOLUTION Chromosomal locations of the 45S pre-ribosomal RNA gene (rna45s), which encodes a precursor RNA for 18S, 5.8S, and 28S rRNAs, was determined using pHr21Ab (5.8-kb for the 5′ portion) and pHr14E3 (7.3-kb for the 3′ portion) fragments as FISH probes. DNA fragments used for the probes were provided by National Institutes of Biomedical Innovation, Health and Nutrition, Osaka, and labeled with biotin-16-dUTP (Roche Diagnostics) by nick translation. After hybridization, the slides were incubated with FITC-avidin (Vector Laboratories). Hybridization signals (arrows) were detected to the short arm of XLA3L, but not XLA3S. Scale bar represents 5 μm.A large deletion including an olfactory receptor gene (or) cluster. Schematic structures of or gene clusters and adjacent genes on the 8th chromosomes of X. tropicalis (XTR8) and X. laevis (XLA8L and XLA8S). Chromosomal locations: XTR8:107,524,547-108,927,581; XLA8L:105,062,063–106,610,199; XLA8S:91,630,596–92,060,451. Horizontal bars, genomic DNA sequences; triangles, genes. Outside of or gene cluster, only representative genes are shown. The length of triangle is to scale. The orientation of triangles indicates 5′ to 3′ direction of genes. Thin lines connect orthologous/homeologous genes. Magenta triangles, or genes; green triangles, pseudogenes (point-mutated or truncated or genes). The number of or genes is shown underneath gene clusters. Dotted lines, a deleted region in XLA8S in comparison with XLA8L. The centromere is located on the left side and the telomere is on the right.The relative frequency (left panel) and size (right panel) of genomic regions deleted in the S (blue) and L (green) chromosomes respectively. Both subgenomes experienced sequence loss through deletions, however, the deletions on the S subgenome are larger and have been more frequent. Deletions were called based on the progressive Cactus sequence alignment between the X.laevis L and S subgenomes and the X. tropicalis genome. Chromosome 9_10 of laevis was split into 9 and 10 on basis of alignment with the X. tropicalis chromosomes. Sequences from L that were not present on S, but could at least partially be identified in X. tropicalis, and consisted of gaps for no more than 25% of their length were called as deleted regions in S. The same procedure was followed for deleted regions in L.Identification of triplet loci is described in Supplemental Note 8.1. Loci were classified into groups based on the presence of gene 2 in both X. laevis subgenomes (homeolog retained), versus those that had a pseudogene in the middle (pseudogene) or no remnant of the middle gene as assessed by Exonerate (deletion). To normalize the intergenic lengths we divided the nucleotide distance between genes 1 and 3 in either X. laevis subgenome by the orthologous distance in X. tropicalis. The median of the normalized ratio distribution is plotted on the bar chart. On average S deletions appear to be larger than L deletions (52.9% length vs 80.2% the size of the orthologous X. tropicalis region respectively).The number of RNA-seq reads aligning +/− 1kb of precursor miRNA loci (red) was compared to the read count for 10,000 random unannotated 2.1 kb regions of the genome (blue). All 83 homeologous, intergenic miRNA pairs showed alignment within their regions, as opposed to 4,127/10,000 (41.27%) of the randomly chosen intergenic sequences. The putative primary-miRNA loci have a higher read count than the expressed randomly chosen regions as well (Wilcoxon p=1.4E-38).The CACTUS alignment was parsed to identify flanking CNE around each X. tropicalis gene. The number of CNEs > 50bp in length for singletons is shown in red, homeologs in blue. Komologrov-Smirnov test p-value is 1E-11.The average distance to the nearest gene was computed for each chromosomal locus in X. tropicalis. The average intergenic distance for those with a single X. laevis gene is shown in red, those with two shown in blue. Wilcoxon p-value= 9.8E-24.The distribution of gene retention by genomic footprint of the X. tropicalis ortholog. We define genomic footprint as the genomic distance from the start signal of the CDS to the stop signal, including introns. The x axis shows log10(genomic footprint), the y-axis is the retention rate of each bin. The error bars are the standard deviation of the total divided by the number of genes in each bin. We tested for significant differences in length between homeologs and singletons by a Wilcoxon test (p-value = 2.4E-96).The distribution of gene retention by CDS length of the X. tropicalis ortholog. The x axis shows log10(CDS length), the y-axis is the retention rate of each bin. The error bars are the standard deviation of the total divided by the number of genes in each bin. We tested for significant differences in length between homeologs and singletons by a Wilcoxon test (p-value= 1.7E-21).The distribution of gene retention by exon number of the X. tropicalis ortholog. The x axis shows number of exons; the y-axis is the retention rate of each bin. The error bars are the standard deviation of the total divided by the number of genes in each bin. We tested for significant differences in length between homeologs and singletons by a Wilcoxon test (p-value= 3.2E-8).|
|EDF 6. PSEUDOGENES Illustration of htt.S pseudogene alignment to X. tropicalis htt, and the extant X. laevis htt.L, translated to amino acids. The amino acid position is shown at the beginning of each line. Missing codons are marked by ‘−’. Frameshifts and premature stops are marked by ‘X’ and ‘*’ respectively (and pointed to with red arrows). (top) the first exon of the pseudogene is completely missing from the S chromosome. The characteristic poly-Q region is maintained by both htt and htt.L. (bottom) An exon with conservation in the pseudogene, illustrating that despite many frameshifts, premature stops, and the lack of a proper start, and insertions of new sequence, we identify many codons in the pseudogene that occur in large conserved blocks.Illustration of our model to compute pseudogene ages. The star represents the point of nonfunctionalization for a currently pseudogenized locus. We assume the expected rate of nonsynonymous changes can be estimated by the Ka of the extant gene and X. tropicalis. We then compare the Ks and Ka of the pseudogene sequence to estimate the time of nonfunctionalization. See Supplemental Note 9 for a more detailed discussion.Estimated epochs of pseudogenization for 430 genes are indistinguishable from a burst of pseudogenization > 10 Mya (Ks > 0.03). See Supplemental Note 9 for a more detailed discussion.Correlation of pseudogene expression with its extant homeolog. The little expression seen in pseudogenes tends to be uncorrelated with the extant homeolog.Histogram of pseudogene expression values across all 28 tissues and developmental stages (red) compared to all extant genes (blue). The pseudogenes are rarely expressed, and tend to be expressed at lower levels than extant protein-coding genes.Histograms of expression variance of pseudogenes (red) compared to extant genes (blue). The small amount of pseudogene expression observed does not tend to vary across tissues and developmental stages in the same way that extant genes do.|
|EDF 7. TANDEM DUPLICATIONS Phylogenetic trees of the mix/bix cluster. Nucleotide sequences were aligned using MUSCLE, and a phylogenetic diagram was generated by the ML method with 1,000 bootstraps (MEGA6). Circles with different colors represent X. laevis L genes (magenta), X. laevis S genes (blue), and X. tropicalis genes (green). The table shows the correspondence of bix gene names proposed in this study and previously used (synonyms).FISH analysis showing XLA3S-specific deletion of the nodal5 gene cluster. One unit of the nodal5 gene region, including exons, introns, and an intergenic region was used as a probe for FISH (counterstained with Hoechst). Arrows indicate the hybridization signals of nodal5s. Scale bar indicates 5 um.Comparison of the nodal5 gene cluster. Genome sequencing revealed that nodal5.e1.L~.e5.L (in pink) and nodal6.L are clustered. Amplification of nodal5 gene in XLA3L and loss of this cluster in XLA3S were confirmed. Pseudogenes (nodal5p1.L~p4.L and nodal5p1.S) are indicated in black. The nodal5 cluster of X. tropicails does not contain any pseudogene.X. laevis L chromosome has four complete copies of nodal3 (nodal3.e1.L~.e4.L), whereas the gene cluster is lost from the X. laevis S chromosome. A truncated nodal3 gene (nodal3p1.L) is likely to be a pseudogene, and highly degenerate pseudogenes (nodal3p2.L and nodal3p3.L) also exist on the L chromosome.Like nodal3, vg1 is lost from the S chromosome although there is a pseudogene (vg1p.S). vg1 is specifically amplified on the X. laevis L chromosome (vg1.e1.L~.e3.L) in comparison with X. tropicalis. An amino acid change (Ser20 to Pro20) in Vg1 protein has been shown to result in functional differences (Supplementary Note 13.9). vg1 and derrière are orthologous to mammalian gdf1.Fraction of all genes duplicated and retained to present epoch per 1 expected 4DTV(four-fold degenerate transversions) at different epochs (semi-log scale). Shown also are linear fits, which would be consistent with constant birth- and death rate models (first epoch is omitted from both fitted data sets, as is second epoch from X. laevis). See Supplemental Note 11 for a more detailed discussion.Same, but for “short genes” (CDS < 600 bp) and “long genes” (CDS > 1200 bp) separately. The loss rate of new duplicates appears to be similar. If the extra copy of a newly duplicated gene were lost when the first 100% disabling mutation occurred, we would expect, on average, the longer genes to be lost.|
|EDF 8. GENE EXPRESSION ANALYSIS Pairwise Pearson correlation distributions between homeologous genes (red) and all genes (blue). The left histogram is for stage data; right is for adult data. The x-axis is the correlation; the y-axis is the percent of data. The homeologous genes have a correlation distribution closer to one due to their being the same locus recently. X. laevis TPM values 0.5 were lowered to 0. Any gene with no TPM > 0 was removed from analysis. We then added 0.1 to all TPM values and log transformed (log10).Scatterplot comparing binned genes by their median X. tropicalis expression64 to the retention rate of their X. laevis (co)-orthologs. Error bars are the standard deviation for the whole data set divided by the square root of the number of genes analyzed in a bin. We assessed significance by a Wilcoxon test of the homeologous and singleton distributions, p-value = 6.31E-113.Complete boxplot shown in Fig. 4c. The difference between subgenomes is difficult to see at this magnification, illustrating that many loci deviate from the whole genome median of preferring the L homeolog. There are some L outliers expressed 104 as much as their S homeologs, whereas no S genes shows such a strong trend. These differences are discussed in more detail in Supplemental Note 12.Boxplot of 4DTv (four-fold degenerate transversions) by homeolog class defined in Supplemental Note 12.4. Significant differences are marked by a red asterisk (Wilcoxon p<1E-5). HCSE group shows lower sequence change than others (p=3.7E-12) and the NCDE group shows high rates of sequence change (p=5.6E-14).Boxplot of CDS length difference between X. laevis homeologs by homeolog class defined in Supplemental Note 12.4. Significant differences are marked by a red asterisk (Wilcoxon p<1E-5). HCSE group shows smaller CDS length differences than others (p=2.4E-13) and the NCDE group shows large differences in homeolog CDS length (p=2.1E-32).Boxplot of Ka/Ks between X. laevis homeologs by homeolog class defined in Supplemental Note 12.4. Significant differences are marked by a red asterisk (t-test p<1E-5). HCSE group shows lower non-synonymous sequence change than others (p=8.2E-19) and the NCDE and NCSE groups shows higher rates of non-synonymous sequence change (p=2.0E-12 and p=7.0E-9 respectively).RNA-seq analysis of six6.L (red) and six6.S (blue) during X. laevis development (left panel) and in the adult tissues (right panel). Expression levels of six6.S were lower than those of six6.L at most developmental stages and in adult tissues.Diagram of Homo sapiens, X. tropicalis and X. laevis six6 loci (upper panel). Magenta and black boxes indicate CNEs and exons, respectively. The phylogenetic tree analyses of H. sapiens, X. tropicalis and X. laevis six6 CNEs (lower left panel), and Six6 proteins (lower right panel). Notably, six6.S is more diverged from X. tropicalis six6 than six6.L, both in the encoded protein sequences and in conserved non-coding elements (CNEs) within 3 kb from the transcription start sites. Materials, methods and the CNE locations on genome assemblies are described in Supplementary Materials (Supplementary Note 13.1).On the basis of chromatin state properties, a Random Forest machine-learning algorithm can accurately predict L versus S expression bias. The classification is based on all genes with greater than 3-fold expression difference at NF stage 10.5 (a set of 1,129 genes). The mean (dotted black line) of the ROC area under the curve is 0.778 (10-fold cross-validation). Features were selected using Linear Support Vector Classification and are shown in Extended Data Fig. 8j.Relative importance (based on Gini impurity) of selected features used in the Random Forest classification. All features used in the classification are shown. Among various variables, the ratios of H3K4me3 and DNA methylation at the promoter contributed most to the decision tree model. A difference in p300 binding in the genomic region surrounding the gene also contributed to the Random Forest classification, as did the presence or absence of a number of specific transcription factor motifs in the promoter.|
|EDF 9. EXAMPLES OF PATHWAY RESPONSES Wnt pathway. Left panel: Several key components of the canonical Wnt pathway in the X. laevis genome. The numbers in brackets show the number of paralogs. Components that have homeologous pair of genes or singleton are shown in blue and red, respectively. Each (Wnt:21 genes, LRP:2 genes, Fzd:10 genes, Dvl:3 genes, Frat(GBP):1 gene, GSK3:2 genes, Axin:2 genes, bcatenin: 1 gene, APC: 2 genes, TCF/LEF:4 genes) were classified into 4 groups according to subcellular localization, and the number of singleton and homeolog retained genes is shown by pie charts. Right panel: Syntenies around four singleton genes.Cell cycle. Upper right panel: Diagram of the cell cycle and regulatory proteins critical to each phase. Cyclin H (CcnH) and Cdk7 constitute Cdk-activating kinase (CAK), a key factor required for activation of all Cdks. Genes encoding Cyclin H and Cdk7 (red), but not other regulators (blue), became singletons. Upper left panel: Pie charts show the numbers of homeologous pairs (blue) and singletons (red) in each functional category as indicated. Lower left panel: Syntenies of ccnh and cdk7 loci in X. tropicalis and X. laevis. Abbreviations for species and chromosome numbers: X. tropicalis (XTR1), X. laevis (XLA1L and XLA1S). Lower right table: Individual genes used for drawing the pie charts are shown in the table.Hippo pathway. Upper panel: Hippo pathway components and retention of their homeologous gene pairs. All genes for Hippo pathway components as indicated were identified in the whole genome of X. laevis. Blue icons indicate that both of the homeologous genes are expressed in normal development and adult organs. The red icon, Taz, indicates a singleton. Yap is interchangeable with Taz in most cases, but TAZ, but not YAP, serves as a mediator of Wnt signaling (broken line). Pie charts show the numbers of homeolog pairs (blue) and singleton (red) in each category of Hippo pathway components classified according to subcellular localization. Lower panel: Comparative analysis of syntenies around the taz gene. X. tropicalis scaffold247 is not incorporated into the chromosome-scale assembly (v9) and hence its chromosomal location is not known yet. The p arm termini of XLA8L and XLA8S are on the left.|
|EDF 10. PATHWAYS PART 2 TGF-beta pathway. Pie charts indicate the ratio of differentially expressed homeologous pairs (orange) and singleton (red). A large portion of the extracellular regulatory factors is either differentially regulated or became singleton. Genes for a type I receptor, co-receptors, an inhibitory Smad are also differentially regulated. Multicopy genes like nodal3, nodal5, and vg1 are not counted as singletons, though those S genes are deleted. Instead, these and duplicated chordin genes are categorized into differentially regulated genes.Hedgehog pathway. Upper panel: The simplified Hedgehog pathway known in Shh signalling is schematically shown. Most signalling components are encoded by both homeologous genes, whereas Hhat (shown in red) is encoded by a singleton gene. In case paralogs exist, the numbers of paralogs are shown in parenthesis. In the left cell, the Shh precursor (Hh precursor) is matured through the process involving Hhat and Hhatl, and secreted. In the right cell, the binding of Shh (Hh) to Ptch1 (Ptch) receptor inhibits Ptch1-mediated repression of Smo, leading to Smo activation and subsequent inhibition of PKA; otherwise PKA converts Gli activators to truncated repressors. As a consequence, Gli proteins activate target genes, such as Ptch1 and Hhip. The transmembrane protein Hhip binds Shh and suppresses Shh activity. Lower panel: Schematic comparison of syntenies around hhat genes of X. tropicalis chromosome 5 (top) and X. laevis 5L chromosome (middle), and the corresponding region of X. laevis 5S chromosome (bottom). The diagram is not drawn to scale.Deletions rates on L (x-axis), vs S (y-axis) for different Pfam groups. For Pfam groups we computed the number of X. laevis single-copy genes (singletons) vs homeolog pairs and computes a fraction retained. The line is expected L/S loss based on genome-wide average (56.4%). Red points show groups with high or low rates of loss (p<.01). See Supplemental Table 5 for more information.Deletions rates on L (x-axis), vs S (y-axis) for different stage WGCNA groups (visualized as a heatmap in Fig. 4a). For stage WGCNA groups we computed the number of X. laevis single-copy genes (singletons) vs homeolog pairs and computes a fraction retained. The line is expected L/S loss based on genome-wide average (56.4%). Red points show groups with high or low rates of loss (p<.01).Deletion rates on L (x-axis), vs S (y-axis) for different GO groups. For GO groups we computed the number of X. laevis single-copy genes (singletons) vs homeolog pairs and computes a fraction retained. The line is expected L/S loss based on genome-wide average (56.4%). Red points show groups with high or low rates of loss (p<0.01). See Supplemental Table 5 for more information.|
|EDF 2. CHROMOSOME STRUCTURE Structure of the sex chromosome of X. laevis (XLA2L) and comparison with XLA2S and XTR2. The W version of XLA2L harbors W-specific sequence containing the female sex-determining gene, dmw (red), while Z has a different Z-specific sequence (blue). Pentagon arrows and black triangles indicate genes and olfactory receptor genes, respectively. Their tips correspond to their 3’-ends.Alignment of the q-terminal regions of XTR9 and 10 with corresponding regions of XLA9_10L and XLA9_10S. Genes near the q terminal regions of XTR 9 and XTR10 were missing in the X. tropicalis genome assembly v9, but rps11, rpl13a, lypd1, and actr3 were expected to be located there based on the synteny with human chromosomes, and then verified by cDNA FISH (upper panels). Small triangles on XLA9_10L and S indicate the distribution of gene models showing both identity and coverage greater than 30%, against the human and chicken peptide sequences from Ensembl, in the region between ±2 Mb from the prospective 9/10 junction. HSA: human chromosome. GGA: chicken chromosome. The magnified view represents syntenic genes to scale with colors corresponding to human genes.The orders of orthologous genes across XTR9, XTR10, XLA9_10L and XLA9_10S. Green arrowheads: positions of centromeres in XTR9 and 10 predicted by examination of the cytogenetic chromosome length ratio of p versus q arms15. Blue arrowheads: positions of centromere repeats, frog centromeric repeat-162, in XLA9_10L and S. Magenta and yellow ellipses: chromosomal locations of snrpn (magenta) and stau1 (yellow) from X. tropicalis v9 and X. laevis v9.1 assemblies. Red ellipses: chromosomal locations of four genes, rps11, rpl13a, lypd1, and actr3. XTR9 is flipped to facilitate comparison. Blue bidirectional arrows indicate the homologous regions where pericentric inversions may have occurred on proto-chromosomes (see Extended Data Fig. 2d).Schematic representation for the two hypothetical processes of chromosomal rearrangements (fusion and inversion) that occurred between the hypothetical proto-XTR9 and 10 to produce proto-XLA9_10, and eventually XLA9_10L and S. The process of chromosome rearrangements is explained parsimoniously in two different ways (left and right panels), starting from proto-XTR9 and 10. Actual and hypothetical ancestral chromosomal locations of snrpn and stau1 are shown by magenta and yellow circles, respectively. Note that the chromosomal locations of these genes on the proto-XTR10 differ between the two models. Chromosome segments homologous to XTR9 and XTR10 are shown in red and blue, respectively. XTR9 is inverted to facilitate comparison. Bidirectional arrows indicate the regions where pericentric inversions may have occurred. Black arrows indicate the direction of chromosomal evolution.|