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.

/pdf/genome-evolution-in-the-allotetraploid-frog-xenopus-laevis-5dt84by3ni.pdf

Genome evolution in the allotetraploid frog Xenopus laevis

There are two intriguing paradoxes in molecular biology-the inconsistent relationship between organismal complexity and (1) cellular DNA content and (2) the number of protein-coding genes-referred to as the C-value and G-value paradoxes, respectively. The C-value paradox may be largely explained by varying ploidy. The G-value paradox is more problematic, as the extent of protein coding sequence remains relatively static over a wide range of developmental complexity. We show by analysis of sequenced genomes that the relative amount of non-protein-coding sequence increases consistently with complexity. We also show that the distribution of introns in complex organisms is non-random. Genes composed of large amounts of intronic sequence are significantly overrepresented amongst genes that are highly expressed in the nervous system, and amongst genes downregulated in embryonic stem cells and cancers. We suggest that the informational paradox in complex organisms may be explained by the expansion of cis-acting regulatory elements and genes specifying trans-acting non-protein-coding RNAs.

The relationship between non-protein-coding DNA and eukaryotic complexity.

A central challenge in evolutionary biology is to identify genes underlying ecologically important traits and describe the fitness consequences of naturally occurring variation at these loci. To address this goal, several novel approaches have been developed, including 'population genomics,' where a large number of molecular markers are scored in individuals from different environments with the goal of identifying markers showing unusual patterns of variation, potentially due to selection at linked sites. Such approaches are appealing because of (1) the increasing ease of generating large numbers of genetic markers, (2) the ability to scan the genome without measuring phenotypes and (3) the simplicity of sampling individuals without knowledge of their breeding history. Although such approaches are inherently applicable to non-model systems, to date these studies have been limited in their ability to uncover functionally relevant genes. By contrast, quantitative genetics has a rich history, and more recently, quantitative trait locus (QTL) mapping has had some success in identifying genes underlying ecologically relevant variation even in novel systems. QTL mapping, however, requires (1) genetic markers that specifically differentiate parental forms, (2) a focus on a particular measurable phenotype and (3) controlled breeding and maintenance of large numbers of progeny. Here we present current advances and suggest future directions that take advantage of population genomics and quantitative genetic approaches - in both model and non-model systems. Specifically, we discuss advantages and limitations of each method and argue that a combination of the two provides a powerful approach to uncovering the molecular mechanisms responsible for adaptation.

/pdf/combining-population-genomics-and-quantitative-genetics-4yv2897y2l.pdf

Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits.

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

/pdf/the-african-coelacanth-genome-provides-insights-into-7lixs866y7.pdf

The African coelacanth genome provides insights into tetrapod evolution

Biological invasions dramatically affect the distribution, abundance and reproduction of many native species. Because of these ecological effects, exotic species can also influence the evolution of natives exposed to novel interactions with invaders. Evolutionary changes in natives in response to selection from exotics are usually overlooked, yet common responses include altered anti-predator defenses, changes in the spectrum of resources and habitats used, and other adaptations that allow native populations to persist in invaded areas. Whether a native population is capable of responding evolutionarily to selection from invaders will depend on the demographic impact of the invader, the genetic architecture and genetic variability of the native population and potentially the history of previous invasions. In some cases, natives will fail to evolve or otherwise adapt, and local or global extinction will result. In other cases, adaptive change in natives may diminish impacts of invaders and potentially promote coexistence between invaders and natives. Here, we review the evidence for evolutionary responses of native species to novel community members. We also discuss how the effects of introduced species may differ from those caused by natural range expansions of native species. Notably, introduced species may come from remote biotas with no previous evolutionary history with the native community. In addition, the rate of addition of introduced species into communities is much greater than all but the most extreme cases of historical biotic exchange. Understanding the evolutionary component of exotic/native species interactions is critical to recognizing the long-term impacts of biological invasions, and to understanding the role of evolutionary processes in the assembly and dynamics of natural communities.

/pdf/evolutionary-responses-of-natives-to-introduced-species-what-2zi6tua1o1.pdf

Evolutionary responses of natives to introduced species: what do introductions tell us about natural communities?

Amphibian genomes differ greatly in DNA content and chromosome size, morphology, and number. Investigations of this diversity are needed to identify mechanisms that have shaped the evolution of vertebrate genomes. We used comparative mapping to investigate the organization of genes in the Mexican axolotl (Ambystoma mexicanum), a species that presents relatively few chromosomes (n = 14) and a gigantic genome (>20 pg/N). We show extensive conservation of synteny between Ambystoma, chicken, and human, and a positive correlation between the length of conserved segments and genome size. Ambystoma segments are estimated to be four to 51 times longer than homologous human and chicken segments. Strikingly, genes demarking the structures of 28 chicken chromosomes are ordered among linkage groups defining the Ambystoma genome, and we show that these same chromosomal segments are also conserved in a distantly related anuran amphibian (Xenopus tropicalis). Using linkage relationships from the amphibian maps, we predict that three chicken chromosomes originated by fusion, nine to 14 originated by fission, and 12-17 evolved directly from ancestral tetrapod chromosomes. We further show that some ancestral segments were fused prior to the divergence of salamanders and anurans, while others fused independently and randomly as chromosome numbers were reduced in lineages leading to Ambystoma and Xenopus. The maintenance of gene order relationships between chromosomal segments that have greatly expanded and contracted in salamander and chicken genomes, respectively, suggests selection to maintain synteny relationships and/or extremely low rates of chromosomal rearrangement. Overall, the results demonstrate the value of data from diverse, amphibian genomes in studies of vertebrate genome evolution.

/pdf/origin-of-amphibian-and-avian-chromosomes-by-fission-fusion-3xec82jyfa.pdf

Origin of amphibian and avian chromosomes by fission, fusion, and retention of ancestral chromosomes

The evolution of alternate modes of development may occur through genetic changes in metamorphic timing. This hypothesis was examined by crossing salamanders that express alternate developmental modes: metamorphosis vs. paedomorphosis. Three strains were used in the crossing design: Ambystoma tigrinum tigrinum (Att; metamorph), wild-caught A. mexicanum (Am; paedomorph), and laboratory Am (paedomorph). Att/Am hybrids were created for each Am strain and then backcrossed to their respective Am line. Previous studies have shown that a dominant allele from Att (metAtt) and a recessive allele from lab Am (metlab) results in metamorphosis in Att/Am hybrids, and metAtt/metlab and metlab/metlab backcross genotypes are strongly associated with metamorphosis and paedomorphosis, respectively. We typed a molecular marker (contig325) linked to met and found that metAtt/metlab and metAtt/metwild were associated with metamorphosis in 99% of the cases examined. However, the frequency of paedomorphosis was 4.5 times higher for metlab/metlab than for metwild/metwild. We also found that metAtt/metwild and metwild/metwild genotypes discriminated distributions of early and late metamorphosing individuals. Two forms of phenotypic variation are contributed by met: continuous variation of metamorphic age and expression of discrete, alternate morphs. We suggest that the evolution of paedomorphosis is associated with genetic changes that delay metamorphic timing in biphasic life cycles.

https://www.genetics.org/content/genetics/170/1/275.full.pdf

Evolution of salamander life cycles : A major-effect quantitative trait locus contributes to discrete and continuous variation for metamorphic timing

Expressed sequence tag (EST) markers were developed for Ambystoma tigrinum tigrinum (Eastern tiger salamander) and for A. mexicanum (Mexican axolotl) to generate the first comprehensive linkage map for these model amphibians. We identified 14 large linkage groups (125.5–836.7 cM) that presumably correspond to the 14 haploid chromosomes in the Ambystoma genome. The extent of genome coverage for these linkage groups is apparently high because the total map size (5251 cM) falls within the range of theoretical estimates and is consistent with independent empirical estimates. Unlike most vertebrate species, linkage map size in Ambystoma is not strongly correlated with chromosome arm number. Presumably, the large physical genome size (∼30 Gbp) is a major determinant of map size in Ambystoma. To demonstrate the utility of this resource, we mapped the position of two historically significant A. mexicanum mutants, white and melanoid, and also met, a quantitative trait locus (QTL) that contributes to variation in metamorphic timing. This new collection of EST-based PCR markers will better enable the Ambystoma system by facilitating development of new molecular probes, and the linkage map will allow comparative studies of this important vertebrate group.

https://www.genetics.org/content/genetics/early/2005/08/03/genetics.105.046433.full.pdf

A Comprehensive Expressed Sequence Tag Linkage Map for Tiger Salamander and Mexican Axolotl: Enabling Gene Mapping and Comparative Genomics in Ambystoma

Thyroid hormones (TH) induce gene expression programs that orchestrate amphibian metamorphosis. In contrast to anurans, many salamanders do not undergo metamorphosis in nature. However, they can be induced to undergo metamorphosis via exposure to thyroxine (T4). We induced metamorphosis in juvenile Mexican axolotls (Ambystoma mexicanum) using 5 and 50 nM T4, collected epidermal tissue from the head at four time points (Days 0, 2, 12, 28), and used microarray analysis to quantify mRNA abundances. Individuals reared in the higher T4 concentration initiated morphological and transcriptional changes earlier and completed metamorphosis by Day 28. In contrast, initiation of metamorphosis was delayed in the lower T4 concentration and none of the individuals completed metamorphosis by Day 28. We identified 402 genes that were statistically differentially expressed by ≥ two-fold between T4 treatments at one or more non-Day 0 sampling times. To complement this analysis, we used linear and quadratic regression to identify 542 and 709 genes that were differentially expressed by ≥ two-fold in the 5 and 50 nM T4 treatments, respectively. We found that T4 concentration affected the timing of gene expression and the shape of temporal gene expression profiles. However, essentially all of the identified genes were similarly affected by 5 and 50 nM T4. We discuss genes and biological processes that appear to be common to salamander and anuran metamorphosis, and also highlight clear transcriptional differences. Our results show that gene expression in axolotls is diverse and precise, and that axolotls provide new insights about amphibian metamorphosis.

/pdf/effect-of-thyroid-hormone-concentration-on-the-2f9ykbr7fn.pdf

Effect of thyroid hormone concentration on the transcriptional response underlying induced metamorphosis in the Mexican axolotl (Ambystoma)

The transformation of ancestral phenotypes into novel traits is poorly understood for many examples of evolutionary novelty. Ancestrally, salamanders have a biphasic life cycle with an aquatic larval stage, a brief and pronounced metamorphosis, followed by a terrestrial adult stage. Repeatedly during evolution, metamorphic timing has been delayed to exploit growth-permissive environments, resulting in paedomorphic salamanders that retain larval traits as adults. We used thyroid hormone (TH) to rescue metamorphic phenotypes in paedomorphic salamanders and then identified quantitative trait loci (QTL) for life history traits that are associated with amphibian life cycle evolution: metamorphic timing and adult body size. We demonstrate that paedomorphic tiger salamanders (Ambystoma tigrinum complex) carry alleles at three moderate effect QTL (met1–3) that vary in responsiveness to TH and additively affect metamorphic timing. Salamanders that delay metamorphosis attain significantly larger body sizes as adults and met2 explains a significant portion of this variation. Thus, substitution of alleles at TH-responsive loci suggests an adaptive pleiotropic basis for two key life-history traits in amphibians: body size and metamorphic timing. Our study demonstrates a likely pathway for the evolution of novel paedomorphic species from metamorphic ancestors via selection of TH-response alleles that delay metamorphic timing and increase adult body size.

/pdf/thyroid-hormone-responsive-qtl-and-the-evolution-of-1lr06jzida.pdf

S. R. Voss

Papers

Origin of amphibian and avian chromosomes by fission, fusion, and retention of ancestral chromosomes

Evolution of salamander life cycles : A major-effect quantitative trait locus contributes to discrete and continuous variation for metamorphic timing

A Comprehensive Expressed Sequence Tag Linkage Map for Tiger Salamander and Mexican Axolotl: Enabling Gene Mapping and Comparative Genomics in Ambystoma

Effect of thyroid hormone concentration on the transcriptional response underlying induced metamorphosis in the Mexican axolotl (Ambystoma)

Thyroid hormone responsive QTL and the evolution of paedomorphic salamanders