The odds of duplicate gene persistence after polyploidization

Genomic insights into the chromosomal elongation in a family of Collembola

Collembola is a highly diverse and abundant group of soil arthropods with chromosome numbers ranging from 5 to 11. Previous karyotype studies indicated that the Tomoceridae family possesses an exceptionally long chromosome. To better understand chromosome size evolution in Collembola, we obtained a chromosome-level genome of Yoshiicerus persimilis with a size of 334.44 Mb and BUSCO completeness of 97.0% (n = 1013). Both genomes of Y. persimilis and Tomocerus qinae (recently published) have an exceptionally large chromosome (ElChr greater than 100 Mb), accounting for nearly one-third of the genome. Comparative genomic analyses suggest that chromosomal elongation occurred independently in the two species approximately 10 million years ago, rather than in the ancestor of the Tomoceridae family. The ElChr elongation was caused by large tandem and segmental duplications, as well as transposon proliferation, with genes in these regions experiencing weaker purifying selection (higher dN/dS) than conserved regions. Moreover, inter-genomic synteny analyses indicated that chromosomal fission/fusion events played a crucial role in the evolution of chromosome numbers (ranging from 5 to 7) within Entomobryomorpha. This study provides a valuable resource for investigating the chromosome evolution of Collembola.

Duplicate genes as sources for rapid adaptive evolution of sperm under environmental pollution in tree sparrow

Environmental pollution can result in poor sperm quality either directly or indirectly in birds. However, adaptive and compensatory sperm morphology changes and motility improvements have rapidly evolved in tree sparrows (Passer montanus) inhabiting polluted areas over the past 65 years. To identify the genetic underpinnings of the rapidly evolving sperm phenotype, we carried out population genomics and transcriptomics on tree sparrow populations in the two differently polluted places. We identified a gene encoding the serine/threonine protein kinase PIM1, which may drive rapid phenotypic evolution of sperm. An unprecedented and remarkable expansion of the PIM gene family, caused by tandem and segmental duplication of PIM1, was subsequently observed in the tree sparrow genome. Most PIM1 duplicates showed a testis-specific expression pattern, suggesting that their functions are related to male reproduction. Furthermore, the elevated expression level of PIM1 was consistent with our earlier findings of longer and faster swimming sperm in polluted sites, indicating an important role for duplicated PIM1 in facilitating the rapid evolution of sperm. Our results suggest that duplicated PIM1 provides sources of genetic variation that may enable the rapid evolution of sperm under environmental heavy metal pollution. The findings of this study indicated that duplicated genes can be targets of selection and predominant sources for rapid adaptation to environmental change and shed light on sperm evolution under pollution stress.

Enhanced Loss of Retinoic Acid Network Genes in Xenopus laevis Achieves a Tighter Signal Regulation

Retinoic acid (RA) is a major regulatory signal during embryogenesis produced from vitamin A (retinol) by an extensive, autoregulating metabolic and signaling network to prevent fluctuations that result in developmental malformations. Xenopus laevis is an allotetraploid hybrid frog species whose genome includes L (long) and S (short) chromosomes from the originating species. Evolutionarily, the X. laevis subgenomes have been losing either L or S homoeologs in about 43% of genes to generate singletons. In the RA network, out of the 47 genes, about 47% have lost one of the homoeologs, like the genome average. Interestingly, RA metabolism genes from storage (retinyl esters) to retinaldehyde production exhibit enhanced gene loss with 75% singletons out of 28 genes. The effect of this gene loss on RA signaling autoregulation was studied. Employing transient RA manipulations, homoeolog gene pairs were identified in which one homoeolog exhibits enhanced responses or looser regulation than the other, while in other pairs both homoeologs exhibit similar RA responses. CRISPR/Cas9 targeting of individual homoeologs to reduce their activity supports the hypothesis where the RA metabolic network gene loss results in tighter network regulation and more efficient RA robustness responses to overcome complex regulation conditions.

Lineage-Specific Genes and Family Expansions in Dictyostelid Genomes Display Expression Bias and Evolutionary Diversification during Development

Gene duplications generate new genes that can contribute to expression changes and the evolution of new functions. Genomes often consist of gene families that undergo expansions, some of which occur in specific lineages that reflect recent adaptive diversification. In this study, lineage-specific genes and gene family expansions were studied across five dictyostelid species to determine when and how they are expressed during multicellular development. Lineage-specific genes were found to be enriched among genes with biased expression (predominant expression in one developmental stage) in each species and at most developmental time points, suggesting independent functional innovations of new genes throughout the phylogeny. Biased duplicate genes had greater expression divergence than their orthologs and paralogs, consistent with subfunctionalization or neofunctionalization. Lineage-specific expansions in particular had biased genes with both molecular signals of positive selection and high expression, suggesting adaptive genetic and transcriptional diversification following duplication. Our results present insights into the potential contributions of lineage-specific genes and families in generating species-specific phenotypes during multicellular development in dictyostelids.

Copy number variations and young duplicate genes have high methylation levels in sticklebacks

Gene duplication is an important driver of genomic diversity that can promote adaptive evolution. However, like most mutations, a newly duplicated gene is often deleterious and removed from the genome by drift or natural selection. The early molecular changes that occur soon after duplication therefore may influence the long-term survival of gene duplicates, but relatively little empirical data exist on the events near the onset of duplication before mutations have time to accumulate. In this study, we contrast gene expression and DNA methylation levels of duplicate genes in the threespine stickleback, Gasterosteus aculeatus, including recently emerged duplications that segregate as copy number variations (CNVs). We find that younger duplicate genes have higher levels of promoter methylation than older genes, and that gene CNVs have higher promoter methylation than non-CNVs. These results suggest preferential duplication of highly methylated genes or rapid methylation changes soon after duplication. We also find a negative association between methylation and expression, providing a putative role for methylation in suppressing transcription that compensates for increases in gene copy numbers and promoting paralog retention. We propose that methylation contributes to the longevity of young duplicate genes, extending the window of opportunity for functional divergence via mutation.

Evolution of the Largest Mammalian Genome

The genome of the red vizcacha rat (Rodentia, Octodontidae, Tympanoctomys barrerae) is the largest of all mammals, and about double the size of their close relative, the mountain vizcacha rat Octomys mimax, even though the lineages that gave rise to these species diverged from each other only about 5 Ma. The mechanism for this rapid genome expansion is controversial, and hypothesized to be a consequence of whole genome duplication or accumulation of repetitive elements. To test these alternative but nonexclusive hypotheses, we gathered and evaluated evidence from whole transcriptome and whole genome sequences of T. barrerae and O. mimax. We recovered support for genome expansion due to accumulation of a diverse assemblage of repetitive elements, which represent about one half and one fifth of the genomes of T. barrerae and O. mimax, respectively, but we found no strong signal of whole genome duplication. In both species, repetitive sequences were rare in transcribed regions as compared with the rest of the genome, and mostly had no close match to annotated repetitive sequences from other rodents. These findings raise new questions about the genomic dynamics of these repetitive elements, their connection to widespread chromosomal fissions that occurred in the T. barrerae ancestor, and their fitness effects—including during the evolution of hypersaline dietary tolerance in T. barrerae.

Comparative transcriptomics of cyprinid minnows and carp in a common wild setting: a resource for ecological genomics in freshwater communities

Comparative transcriptomics can now be conducted on organisms in natural settings, which has greatly enhanced understanding of genome–environment interactions. Here, we demonstrate the utility and potential pitfalls of comparative transcriptomics of wild organisms, with an example from three cyprinid fish species (Teleostei:Cypriniformes). We present extensively filtered and annotated transcriptome assemblies that provide a valuable resource for studies of genome evolution (e.g. polyploidy), ecological and morphological diversification, speciation, and shared and unique responses to environmental variation in cyprinid fishes. Our results and analyses address the following points: (i) ‘essential developmental genes’ are shown to be ubiquitously expressed in a diverse suite of tissues across later ontogenetic stages (i.e. juveniles and adults), making these genes are useful for assessing the quality of transcriptome assemblies, (ii) the influence of microbiomes and other exogenous DNA, (iii) potentially novel, species-specific genes, and (iv) genomic rearrangements (e.g. whole genome duplication). The data we present provide a resource for future comparative work in cypriniform fishes and other taxa across a variety of sub-disciplines, including stress response, morphological diversification, community ecology, ecotoxicology, and climate change.

Homoeologous Recombination of the V1r1-V1r2 Gene Cluster of Pheromone Receptors in an Allotetraploid Lineage of Teleosts

In contrast to other olfactory receptor families that exhibit frequent lineage-specific expansions, the vomeronasal type 1 receptor (V1R) family exhibits a canonical six-member repertoire in teleosts. V1r1 and V1r2 are present in no more than one copy in all examined teleosts, including salmons, which are ancient polyploids, implying strict evolutionary constraints. However, recent polyploids have not been examined. Here, we identified a young allotetraploid lineage of weatherfishes and investigated their V1r1-V1r2 cluster. We found a novel pattern that the parental V1r1-V1r2 clusters had recombined in the tetraploid genome and that the recombinant was nearly fixed in the tetraploid population. Subsequent analyses suggested strong selective pressure, for both a new combination of paralogs and homogeneity among gene duplicates, acting on the V1r1-V1r2 pair.

Complex Genes Are Preferentially Retained After Whole-Genome Duplication in Teleost Fish

Gene duplication generates new genetic material which, if retained after duplication, may contribute to organismal evolution. A whole-genome duplication occurred in the ancestry of teleost fish and consequently there are many duplicated genes in teleost genomes. Indeed, it has been proposed that the evolutionary diversification of teleost fish may have been stimulated by the fish-specific genome duplication (FSGD). However, it is not clear which factors determine which genes are retained as duplicate copies and which return to a singleton state after duplication. In the present study, gene complexity, in terms of encoded protein length and functional domain number, is compared between duplicate and singleton genes for nine well-annotated teleost genomes. A total of 933 gene families with retained duplicates and 4590 singleton gene families are analysed. Genes with retained duplicates are found to be significantly longer (27.9-38.2%) and to have more functional domains (20.5-26.5%) than singleton genes in all the nine teleost genomes, suggesting that genes encoded longer proteins with and more functional domains were preferentially retained after whole-genome duplication in teleosts. This differential retention of duplicated genes will have increased the genomic complexity of teleost fish after FSGD which, together with differential duplicated gene retention as a lineage-splitting force, may have greatly contributed to the successful diversification of teleost fish.

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 demonstrate the allotetraploid origin of X. laevis by partitioning its genome into two homeologous subgenomes, marked by distinct families of “fossil” transposable elements. Based on the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged ~34 million years ago (Mya) and combined to form an allotetraploid ~17–18 Mya. 56% of all genes are retained in two homeologous copies. Protein function, gene expression, and the amount of flanking conserved 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.

Population Level Purifying Selection and Gene Expression Shape Subgenome Evolution in Maize

The maize ancestor experienced a recent whole-genome duplication (WGD) followed by gene erosion which generated two subgenomes, the dominant subgenome (maize1) experiencing fewer deletions than maize2. We take advantage of available extensive polymorphism and gene expression data in maize to study purifying selection and gene expression divergence between WGD retained paralog pairs. We first report a strong correlation in nucleotide diversity between duplicate pairs, except for upstream regions. We then show that maize1 genes are under stronger purifying selection than maize2. WGD retained genes have higher gene dosage and biased Gene Ontologies consistent with previous studies. The relative gene expression of paralogs across tissues demonstrates that 98% of duplicate pairs have either subfunctionalized in a tissuewise manner or have diverged consistently in their expression thereby preventing functional complementation. Tissuewise subfunctionalization seems to be a hallmark of transcription factors, whereas consistent repression occurs for macromolecular complexes. We show that dominant gene expression is a strong determinant of the strength of purifying selection, explaining the inferred stronger negative selection on maize1 genes. We propose a novel expression-based classification of duplicates which is more robust to explain observed polymorphism patterns than the subgenome location. Finally, upstream regions of repressed genes exhibit an enrichment in transposable elements which indicates a possible mechanism for expression divergence.

Copy number variation and genetic diversity of MHC Class IIb alleles in an alien population of Xenopus laevis

Xenopus laevis (the African clawed frog), which originated through hybridisation and whole genome duplication, has been used as a model for genetics and development for many years, but surprisingly little is known about immune gene variation in natural populations. The purpose of this study was to use an isolated population of X. laevis that was introduced to Wales, UK in the past 50 years to investigate how variation at the MHC compares to that at other loci, following a severe population bottleneck. Among 18 individuals, we found nine alleles based on exon 2 sequences of the Class IIb region (which includes the peptide binding region). Individuals carried from one to three of the loci identified from previous laboratory studies. Genetic variation was an order of magnitude higher at the MHC compared with three single-copy nuclear genes, but all loci showed high levels of heterozygosity and nucleotide diversity and there was not an excess of homozygosity or decrease in diversity over time that would suggest extensive inbreeding in the introduced population. Tajima’s D was positive for all loci, which is consistent with a bottleneck. Moreover, comparison with published sequences identified the source of the introduced population as the Western Cape region of South Africa, where most commercial suppliers have obtained their stocks. These factors suggest that despite founding by potentially already inbred individuals, the alien population in Wales has maintained substantial genetic variation at both adaptively important and neutral genes.

The Rise and Fall of TRP-N, an Ancient Family of Mechanogated Ion Channels, in Metazoa

Mechanoreception, the sensing of mechanical forces, is an ancient means of orientation and communication and tightly linked to the evolution of motile animals. In flies, the transient-receptor-potential N protein (TRP-N) was found to be a cilia-associated mechanoreceptor. TRP-N belongs to a large and diverse family of ion channels. Its unusually long N-terminal repeat of 28 ankyrin domains presumably acts as the gating spring by which mechanical energy induces channel gating. We analyzed the evolutionary origins and possible diversification of TRP-N. Using a custom-made set of highly discriminative sequence profiles we scanned a representative set of metazoan genomes and subsequently corrected several gene models. We find that, contrary to other ion channel families, TRP-N is remarkably conserved in its domain arrangements and copy number (1) in all Bilateria except for amniotes, even in the wake of several whole-genome duplications. TRP-N is absent in Porifera but present in Ctenophora and Placozoa. Exceptional multiplications of TRP-N occurred in Cnidaria, independently along the Hydra and the Nematostella lineage. Molecular signals of subfunctionalization can be attributed to different mechanisms of activation of the gating spring. In Hydra this is further supported by in situ hybridization and immune staining, suggesting that at least three paralogs adapted to nematocyte discharge, which is key for predation and defense. We propose that these new candidate proteins help explain the sensory complexity of Cnidaria which has been previously observed but so far has lacked a molecular underpinning. Also, the ancient appearance of TRP-N supports a common origin of important components of the nervous systems in Ctenophores, Cnidaria, and Bilateria.

Benchmarking Transcriptome Quantification Methods for Duplicated Genes in Xenopus laevis

Xenopus is an important model organism for the study of genome duplication in vertebrates. With the full genome sequence of diploid Xenopus tropicalis available, and that of allotetraploid X. laevis close to being finished, we will be able to expand our understanding of how duplicated genes have evolved. One of the key features in the study of the functional consequence of gene duplication is how their expression patterns vary across different conditions, and RNA-seq seems to have enough resolution to discriminate the expression of highly similar duplicated genes. However, most of the current RNA-seq analysis methods were not designed to study samples with duplicate genes such as in X. laevis. Here, various computational methods to quantify gene expression in RNA-seq data were evaluated, using 2 independent X. laevis egg RNA-seq datasets and 2 reference databases for duplicated genes. The fact that RNA-seq can measure expression levels of similar duplicated genes was confirmed, but long paired-end reads are more informative than short single-end reads to discriminate duplicated genes. Also, it was found that bowtie, one of the most popular mappers in RNA-seq analysis, reports significantly smaller numbers of unique hits according to a mapping quality score compared to other mappers tested (BWA, GSNAP, STAR). Calculated from unique hits based on a mapping quality score, both expression levels and the expression ratio of duplicated genes can be estimated consistently among biological replicates, demonstrating that this method can successfully discriminate the expression of each copy of a duplicated gene pair. This comprehensive evaluation will be a useful guideline for studying gene expression of organisms with genome duplication using RNA-seq in the future.

Molecular Polymorphism and Divergence of Duplicated Genes in Tetraploid African Clawed Frogs (Xenopus)

Genome duplication creates redundancy in proteins and their interaction networks, and subsequent smaller-scale gene duplication can further amplify genetic redundancy. Mutations then lead to the loss, maintenance or functional divergence of duplicated genes. Genome duplication occurred many times in African clawed frogs (genus Xenopus), and almost all extant species in this group evolved from a polyploid ancestor. To better understand the nature of selective constraints in a polyploid genome, we examined molecular polymorphism and divergence of duplicates and single-copy genes in 2 tetraploid African clawed frog species, Xenopus laevis and X. victorianus. We found that molecular polymorphism in the coding regions of putative duplicated genes was higher than in singletons, but not significantly so. Our findings also suggest that transcriptome evolution in polyploids is influenced by variation in the genome-wide mutation rate, and do not reject the hypothesis that gene dosage balance is also important.

Evidence for widespread subfunctionalization of splice forms in vertebrate genomes

Gene duplication and alternative splicing are important sources of proteomic diversity. Despite research indicating that gene duplication and alternative splicing are negatively correlated, the evolutionary relationship between the two remains unclear. One manner in which alternative splicing and gene duplication may be related is through the process of subfunctionalization, in which an alternatively spliced gene upon duplication divides distinct splice isoforms among the newly generated daughter genes, in this way reducing the number of alternatively spliced transcripts duplicate genes produce. Previously, it has been shown that splice form subfunctionalization will result in duplicate pairs with divergent exon structure when distinct isoforms become fixed in each paralog. However, the effects of exon structure divergence between paralogs have never before been studied on a genome-wide scale. Here, using genomic data from human, mouse, and zebrafish, we demonstrate that gene duplication followed by exon structure divergence between paralogs results in a significant reduction in levels of alternative splicing. In addition, by comparing the exon structure of zebrafish duplicates to the co-orthologous human gene, we have demonstrated that a considerable fraction of exon divergent duplicates maintain the structural signature of splice form subfunctionalization. Furthermore, we find that paralogs with divergent exon structure demonstrate reduced breadth of expression in a variety of tissues when compared to paralogs with identical exon structures and singletons. Taken together, our results are consistent with subfunctionalization partitioning alternatively spliced isoforms among duplicate genes and as such highlight the relationship between gene duplication and alternative splicing.

Pan-African phylogeography of a model organism, the African clawed frog 'Xenopus laevis'

The African clawed frog Xenopus laevis has a large native distribution over much of sub-Saharan Africa and is a model organism for research, a proposed disease vector, and an invasive species. Despite its prominent role in research and abundance in nature, surprisingly little is known about the phylogeography and evolutionary history of this group. Here, we report an analysis of molecular variation of this clade based on 17 loci (one mitochondrial, 16 nuclear) in up to 159 individuals sampled throughout its native distribution. Phylogenetic relationships among mitochondrial DNA haplotypes were incongruent with those among alleles of the putatively female-specific sex-determining gene DM-W, in contrast to the expectation of strict matrilineal inheritance of both loci. Population structure and evolutionarily diverged lineages were evidenced by analyses of molecular variation in these data. These results further contextualize the chronology, and evolutionary relationships within this group, support the recognition of X. laevis sensu stricto, X. petersii, X. victorianus and herein revalidated X. poweri as separate species. We also propose that portions of the currently recognized distributions of X. laevis (north of the Congo Basin) and X. petersii (south of the Congo Basin) be reassigned to X. poweri.

Transcriptome-wide analysis of SAMe superfamily to novelty phosphoethanolamine N-methyltransferase copy in Lonicera japonica

The S-adenosyl-L-methionine-dependent methyltransferase superfamily plays important roles in plant development. The buds of Lonicera japonica are used as Chinese medical material and foods; chinese people began domesticating L. japonica thousands of years ago. Compared to the wild species, L. japonica var. chinensis, L. japonica gives a higher yield of buds, a fact closely related to positive selection over the long cultivation period of the species. Genome duplications, which are always detected in the domestic species, are the source of the multifaceted roles of the functional gene. In this paper, we investigated the evolution of the SAMe genes in L. japonica and L. japonica var. chinensis and further analyzed the roles of the duplicated genes among special groups. The SAMe protein sequences were subdivided into three clusters and several subgroups. The difference in transcriptional levels of the duplicated genes showed that seven SAMe genes could be related to the differences between the wild and the domesticated varieties. The sequence diversity of seven SAMe genes was also analyzed, and the results showed that different gene expression levels between the varieties could not be related to amino acid variation. The transcriptional level of duplicated PEAMT could be regulated through the SAM-SAH cycle.

Whole-genome duplication in teleost fishes and its evolutionary consequences

Whole-genome duplication (WGD) events have shaped the history of many evolutionary lineages. One such duplication has been implicated in the evolution of teleost fishes, by far the most species-rich vertebrate clade. After initial controversy, there is now solid evidence that such event took place in the common ancestor of all extant teleosts. It is termed teleost-specific (TS) WGD. After WGD, duplicate genes have different fates. The most likely outcome is non-functionalization of one duplicate gene due to the lack of selective constraint on preserving both. Mechanisms that act on preservation of duplicates are subfunctionalization (partitioning of ancestral gene functions on the duplicates), neofunctionalization (assigning a novel function to one of the duplicates) and dosage selection (preserving genes to maintain dosage balance between interconnected components). Since the frequency of these mechanisms is influenced by the genes' properties, there are over-retained classes of genes, such as highly expressed ones and genes involved in neural function. The consequences of the TS-WGD, especially its impact on the massive radiation of teleosts, have been matter of controversial debate. It is evident that gene duplications are crucial for generating complexity and that WGDs provide large amounts of raw material for evolutionary adaptation and innovation. However, it is less clear whether the TS-WGD is directly linked to the evolutionary success of teleosts and their radiation. Recent studies let us conclude that TS-WGD has been important in generating teleost complexity, but that more recent ecological adaptations only marginally related to TS-WGD might have even contributed more to diversification. It is likely, however, that TS-WGD provided teleosts with diversification potential that can become effective much later, such as during phases of environmental change.

Faster evolving Drosophila paralogs lose expression rate and ubiquity and accumulate more non-synonymous SNPs

Background

Duplicated genes can indefinately persist in genomes if either both copies retain the original function due to dosage benefit (gene conservation), or one of the copies assumes a novel function (neofunctionalization), or both copies become required to perform the function previously accomplished by a single copy (subfunctionalization), or through a combination of these mechanisms. Different models of duplication retention imply different predictions about substitution rates in the coding portion of paralogs and about asymmetry of these rates.

Results

We analyse sequence evolution asymmetry in paralogs present in 12 Drosophila genomes using the nearest non-duplicated orthologous outgroup as a reference. Those paralogs present in D. melanogaster are analysed in conjunction with the asymmetry of expression rate and ubiquity and of segregating non-synonymous polymorphisms in the same paralogs. Paralogs accumulate substitutions, on average, faster than their nearest singleton orthologs. The distribution of paralogs’ substitution rate asymmetry is overdispersed relative to that of orthologous clades, containing disproportionally more unusually symmetric and unusually asymmetric clades. We show that paralogs are more asymmetric in: a) clades orthologous to highly constrained singleton genes; b) genes with high expression level; c) genes with ubiquitous expression and d) non-tandem duplications. We further demonstrate that, in each asymmetrically evolving pair of paralogs, the faster evolving member of the pair tends to have lower average expression rate, lower expression uniformity and higher frequency of non-synonymous SNPs than its slower evolving counterpart.

Conclusions

Our findings are consistent with the hypothesis that many duplications in Drosophila are retained despite stabilising selection being more relaxed in one of the paralogs than in the other, suggesting a widespread unfinished pseudogenization. This phenomenon is likely to make detection of neo- and subfunctionalization signatures difficult, as these models of duplication retention also predict asymmetries in substitution rates and expression profiles.

Reviewers

This article has been reviewed by Dr. Jia Zeng (nominated by Dr. I. King Jordan), Dr. Fyodor Kondrashov and Dr. Yuri Wolf.

A large pseudoautosomal region on the sex chromosomes of the frog Silurana tropicalis

Sex chromosome divergence has been documented across phylogenetically diverse species, with amphibians typically having cytologically nondiverged (“homomorphic”) sex chromosomes. With an aim of further characterizing sex chromosome divergence of an amphibian, we used “RAD-tags” and Sanger sequencing to examine sex specificity and heterozygosity in the Western clawed frog Silurana tropicalis (also known as Xenopus tropicalis). Our findings based on approximately 20 million genotype calls and approximately 200 polymerase chain reaction-amplified regions across multiple male and female genomes failed to identify a substantially sized genomic region with genotypic hallmarks of sex chromosome divergence, including in regions known to be tightly linked to the sex-determining region. We also found that expression and molecular evolution of genes linked to the sex-determining region did not differ substantially from genes in other parts of the genome. This suggests that the pseudoautosomal region, where recombination occurs, comprises a large portion of the sex chromosomes of S. tropicalis. These results may in part explain why African clawed frogs have such a high incidence of polyploidization, shed light on why amphibians have a high rate of sex chromosome turnover, and raise questions about why homomorphic sex chromosomes are so prevalent in amphibians.

Ancient eudicot hexaploidy meets ancestral eurosid gene order

BACKGROUND

A hexaploidization event over 125 Mya underlies the evolutionary lineage of the majority of flowering plants, including very many species of agricultural importance. Half of these belong to the rosid subgrouping, containing severals whose genome sequences have been published. Although most duplicate and triplicate genes have been lost in all descendants, clear traces of the original chromosome triples can be discerned, their internal contiguity highly conserved in some genomes and very fragmented in others. To understand the particular evolutionary patterns of plant genomes, there is a need to systematically survey the fate of the subgenomes of polyploids, including the retention of a small proportion of the duplicate and triplicate genes and the reconstruction of putative ancestral intermediates between the original hexaploid and modern species, in this case the ancestor of the eurosid clade.

RESULTS

We quantitatively trace the fate of gene triples originating in the hexaploidy across seven core eudicot flowering plants, and fit this to a two-stage model, pre- and post-radiation. We also measure the simultaneous dynamics of duplicate orthologous gene loss in three rosids, as influenced by biological functional class. We propose a new protocol for reconstructing ancestral gene order using only gene adjacency data from pairwise genomic analyses, based on repeating MAXIMUM WEIGHT MATCHING at two levels of resolution, an approach designed to transcend limitations on reconstructed contig size, while still avoiding the ambiguities of a multiplicity of solutions. Applied to three high-quality rosid genomes without subsequent polyploidy events, our automated procedure reconstructs the ancestor of the eurosid clade.

CONCLUSIONS

The gene loss analysis and the ancestor reconstruction present complementary assessments of post-hexaploidization evolution, the first at the level of individual gene families within and across sister genomes and the second at the chromosome level. Despite the loss of more than 95% of gene duplicates and triplicates, and despite major structural rearrangement, our reconstructed eurosid ancestor clearly identifies the three regions corresponding to each of the seven original chromosomes of the earlier pre-hexaploid ancestor. Functional analysis confirmed trends reported for more recent plant polyploidy events: genes involved with regulation and responses were retained in multiple copies, while genes involved with metabolic processes were lost.

Genetic and genomic interactions of animals with different ploidy levels

Polyploid animals have independently evolved from diploids in diverse taxa across the tree of life. We review a few polyploid animal species or biotypes where recently developed molecular and cytogenetic methods have significantly improved our understanding of their genetics, reproduction and evolution. Mitochondrial sequences that target the maternal ancestor of a polyploid show that polyploids may have single (e.g. unisexual salamanders in the genus Ambystoma) or multiple (e.g. parthenogenetic polyploid lizards in the genus Aspidoscelis) origins. Microsatellites are nuclear markers that can be used to analyze genetic recombinations, reproductive modes (e.g. Ambystoma) and recombination events (e.g. polyploid frogs such as Pelophylax esculentus). Hom(e)ologous chromosomes and rare intergenomic exchanges in allopolyploids have been distinguished by applying genome-specific fluorescent probes to chromosome spreads. Polyploids arise, and are maintained, through perturbations of the 'normal' meiotic program that would include pre-meiotic chromosome replication and genomic integrity of homologs. When possible, asexual, unisexual and bisexual polyploid species or biotypes interact with diploid relatives, and genes are passed from diploid to polyploid gene pools, which increase genetic diversity and ultimately evolutionary flexibility in the polyploid. When diploid relatives do not exist, polyploids can interact with another polyploid (e.g. species of African Clawed Frogs in the genus Xenopus). Some polyploid fish (e.g. salmonids) and frogs (Xenopus) represent independent lineages whose ancestors experienced whole genome duplication events. Some tetraploid frogs (P. esculentus) and fish (Squaliusalburnoides) may be in the process of becoming independent species, but diploid and triploid forms of these 'species' continue to genetically interact with the comparatively few tetraploid populations. Genetic and genomic interaction between polyploids and diploids is a complex and dynamic process that likely plays a crucial role for the evolution and persistence of polyploid animals. See also other articles in this themed issue.

Genomic dynamics of transposable elements in the western clawed frog (Silurana tropicalis)

Transposable elements (TEs) are repetitive DNA sequences that can make new copies of themselves that are inserted elsewhere in a host genome. The abundance and distributions of TEs vary considerably among phylogenetically diverse hosts. With the aim of exploring the basis of this variation, we evaluated correlations between several genomic variables and the presence of TEs and non-TE repeats in the complete genome sequence of the Western clawed frog (Silurana tropicalis). This analysis reveals patterns of TE insertion consistent with gene disruption but not with the insertional preference model. Analysis of non-TE repeats recovered unique features of their genome-wide distribution when compared with TE repeats, including no strong correlation with exons and a particularly strong negative correlation with GC content. We also collected polymorphism data from 25 TE insertion sites in 19 wild-caught S. tropicalis individuals. DNA transposon insertions were fixed at eight of nine sites and at a high frequency at one of nine, whereas insertions of long terminal repeat (LTR) and non-LTR retrotransposons were fixed at only 4 of 16 sites and at low frequency at 12 of 16. A maximum likelihood model failed to attribute these differences in insertion frequencies to variation in selection pressure on different classes of TE, opening the possibility that other phenomena such as variation in rates of replication or duration of residence in the genome could play a role. Taken together, these results identify factors that sculpt heterogeneity in TE distribution in S. tropicalis and illustrate that genomic dynamics differ markedly among TE classes and between TE and non-TE repeats.