Genic regions of a large salamander genome contain long introns and novel genes

Genome Assembly of the Dyeing Poison Frog Provides Insights into the Dynamics of Transposable Element and Genome-Size Evolution

Genome size varies greatly across the tree of life and transposable elements are an important contributor to this variation. Among vertebrates, amphibians display the greatest variation in genome size, making them ideal models to explore the causes and consequences of genome size variation. However, high-quality genome assemblies for amphibians have, until recently, been rare. Here, we generate a high-quality genome assembly for the dyeing poison frog, Dendrobates tinctorius. We compare this assembly to publicly available frog genomes and find evidence for both large-scale conserved synteny and widespread rearrangements between frog lineages. Comparing conserved orthologs annotated in these genomes revealed a strong correlation between genome size and gene size. To explore the cause of gene-size variation, we quantified the location of transposable elements relative to gene features and find that the accumulation of transposable elements in introns has played an important role in the evolution of gene size in D. tinctorius, while estimates of insertion times suggest that many insertion events are recent and species-specific. Finally, we carry out population-scale mobile-element sequencing and show that the diversity and abundance of transposable elements in poison frog genomes can complicate genotyping from repetitive element sequence anchors. Our results show that transposable elements have clearly played an important role in the evolution of large genome size in D. tinctorius. Future studies are needed to fully understand the dynamics of transposable element evolution and to optimize primer or bait design for cost-effective population-level genotyping in species with large, repetitive genomes.

The use of transgenics in the laboratory axolotl

The ability to generate transgenic animals sparked a wave of research committed to implementing such technology in a wide variety of model organisms. Building a solid base of ubiquitous and tissue-specific reporter lines has set the stage for later interrogations of individual cells or genetic elements. Compared to other widely used model organisms such as mice, zebrafish and fruit flies, there are only a few transgenic lines available in the laboratory axolotl (Ambystoma mexicanum), although their number is steadily expanding. In this review, we discuss a brief history of the transgenic methodologies in axolotl and their advantages and disadvantages. Next, we discuss available transgenic lines and insights we have been able to glean from them. Finally, we list challenges when developing transgenic axolotl, and where further work is needed in order to improve their standing as both a developmental and regenerative model.

The Axolotl's journey to the modern molecular era

The salamander Ambystoma mexicanum, commonly called "the axolotl" has a long, illustrious history as a model organism, perhaps with one of the longest track records as a laboratory-bred vertebrate, yet it also holds a prominent place among the emerging model organisms. Or rather it is by now an "emerged" model organism, boasting a full cohort molecular genetic tools that allows an expanding community of researchers in the field to explore the remarkable traits of this animal including regeneration, at cellular and molecular precision-which had been a dream for researchers over the years. This chapter describes the journey to this status, that could be helpful for those developing their respective animal or plant models.

Novel Expansion of Matrix Metalloproteases in the Laboratory Axolotl (Ambystoma mexicanum) and Other Salamander Species

Matrix metalloprotease (MMP) genes encode endopeptidases that cleave protein components of the extracellular matrix (ECM) as well as non-ECM proteins. Here we report the results of a comprehensive survey of MMPs in the laboratory axolotl and other representative salamanders. Surprisingly, 28 MMPs were identified in salamanders and 9 MMP paralogs were identified as unique to the axolotl and other salamander taxa, with several of these presenting atypical amino acid insertions not observed in other tetrapod vertebrates. Furthermore, as assessed by sequence information, all of the novel salamander MMPs are of the secreted type, rather than cell membrane anchored. This suggests that secreted type MMPs expanded uniquely within salamanders to presumably execute catalytic activities in the extracellular milieu. To facilitate future studies of salamander-specific MMPs, we annotated transcriptional information from published studies of limb and tail regeneration. Our analysis sets the stage for comparative studies to understand why MMPs expanded uniquely within salamanders.

Newt regeneration genes regulate Wingless signaling to restore patterning in Drosophila eye

Newts utilize their unique genes to restore missing parts by strategic regulation of conserved signaling pathways. Lack of genetic tools poses challenges to determine the function of such genes. Therefore, we used the Drosophila eye model to demonstrate the potential of 5 unique newt (Notophthalmus viridescens) gene(s), viropana1-viropana5 (vna1-vna5), which were ectopically expressed in L 2 mutant and GMR-hid, GMR-GAL4 eye. L 2 exhibits the loss of ventral half of early eye and head involution defective (hid) triggers cell-death during later eye development. Surprisingly, newt genes significantly restore missing photoreceptor cells both in L 2 and GMR>hid background by upregulating cell-proliferation and blocking cell-death, regulating evolutionarily conserved Wingless (Wg)/Wnt signaling pathway and exhibit non-cell-autonomous rescues. Further, Wg/Wnt signaling acts downstream of newt genes. Our data highlights that unique newt proteins can regulate conserved pathways to trigger a robust restoration of missing photoreceptor cells in Drosophila eye model with weak restoration capability.

The rax homeobox gene is mutated in the eyeless axolotl, Ambystoma mexicanum

BACKGROUND

Vertebrate eye formation requires coordinated inductive interactions between different embryonic tissue layers, first described in amphibians. A network of transcription factors and signaling molecules controls these steps, with mutations causing severe ocular, neuronal, and craniofacial defects. In eyeless mutant axolotls, eye morphogenesis arrests at the optic vesicle stage, before lens induction, and development of ventral forebrain structures is disrupted.

RESULTS

We identified a 5-bp deletion in the rax (retina and anterior neural fold homeobox) gene, which was tightly linked to the recessive eyeless (e) axolotl locus in an F2 cross. This frameshift mutation, in exon 2, truncates RAX protein within the homeodomain (P154fs35X). Quantitative RNA analysis shows that mutant and wild-type rax transcripts are equally abundant in E/e embryos. Translation appears to initiate from dual start codons, via leaky ribosome scanning, a conserved feature among gnathostome RAX proteins. Previous data show rax is expressed in the optic vesicle and diencephalon, deeply conserved among metazoans, and required for eye formation in other species.

CONCLUSION

The eyeless axolotl mutation is a null allele in the rax homeobox gene, with primary defects in neural ectoderm, including the retinal and hypothalamic primordia.

Characterization of axolotl lampbrush chromosomes by fluorescence in situ hybridization and immunostaining

The lampbrush chromosomes (LBCs) in oocytes of the Mexican axolotl (Ambystoma mexicanum) were identified some time ago by their relative lengths and predicted centromeres, but they have never been associated completely with the mitotic karyotype, linkage maps or genome assembly. We identified 9 of the axolotl LBCs using RNAseq to identify actively transcribed genes and 13 BAC (bacterial artificial clone) probes containing pieces of active genes. Using read coverage analysis to find candidate centromere sequences, we developed a centromere probe that localizes to all 14 centromeres. Measurements of relative LBC arm lengths and polymerase III localization patterns enabled us to identify all LBCs. This study presents a relatively simple and reliable way to identify each axolotl LBC cytologically and to anchor chromosome-length sequences (from the axolotl genome assembly) to the physical LBCs by immunostaining and fluorescence in situ hybridization. Our data will facilitate a more detailed transcription analysis of individual LBC loops.

Study of Natural Longlife Juvenility and Tissue Regeneration in Caudate Amphibians and Potential Application of Resulting Data in Biomedicine

The review considers the molecular, cellular, organismal, and ontogenetic properties of Urodela that exhibit the highest regenerative abilities among tetrapods. The genome specifics and the expression of genes associated with cell plasticity are analyzed. The simplification of tissue structure is shown using the examples of the sensory retina and brain in mature Urodela. Cells of these and some other tissues are ready to initiate proliferation and manifest the plasticity of their phenotype as well as the correct integration into the pre-existing or de novo forming tissue structure. Without excluding other factors that determine regeneration, the pedomorphosis and juvenile properties, identified on different levels of Urodele amphibians, are assumed to be the main explanation for their high regenerative abilities. These properties, being fundamental for tissue regeneration, have been lost by amniotes. Experiments aimed at mammalian cell rejuvenation currently use various approaches. They include, in particular, methods that use secretomes from regenerating tissues of caudate amphibians and fish for inducing regenerative responses of cells. Such an approach, along with those developed on the basis of knowledge about the molecular and genetic nature and age dependence of regeneration, may become one more step in the development of regenerative medicine.

Appendage regeneration in anamniotes utilizes genes active during larval-metamorphic stages that have been lost or altered in amniotes: The case for studying lizard tail regeneration

This review elaborates the idea that organ regeneration derives from specific evolutionary histories of vertebrates. Regenerative ability depends on genomic regulation of genes specific to the life-cycles that have differentially evolved in anamniotes and amniotes. In aquatic environments, where fish and amphibians live, one or multiple metamorphic transitions occur before the adult stage is reached. Each transition involves the destruction and remodeling of larval organs that are replaced with adult organs. After organ injury or loss in adult anamniotes, regeneration uses similar genes and developmental process than those operating during larval growth and metamorphosis. Therefore, the broad presence of regenerative capability across anamniotes is possible because generating new organs is included in their life history at metamorphic stages. Soft hyaluronate-rich regenerative blastemas grow in submersed or in hydrated environments, that is, essential conditions for regeneration, like during development. In adult anamniotes, the ability to regenerate different organs decreases in comparison to larval stages and becomes limited during aging. Comparisons of genes activated during metamorphosis and regeneration in anamniotes identify key genes unique to these processes, and include thyroid, wnt and non-coding RNAs developmental pathways. In the terrestrial environment, some genes or developmental pathways for metamorphic transitions were lost during amniote evolution, determining loss of regeneration. Among amniotes, the formation of soft and hydrated blastemas only occurs in lizards, a morphogenetic process that evolved favoring their survival through tail autotomy, leading to a massive although imperfect regeneration of the tail. Deciphering genes activity during lizard tail regeneration would address future attempts to recreate in other amniotes regenerative blastemas that grow into variably completed organs.

Exploring the role of microRNAs in axolotl regeneration

The axolotl, Ambystoma mexicanum, is used extensively for research in developmental biology, particularly for its ability to regenerate and restore lost organs, including in the nervous system, to full functionality. Regeneration in mammals typically depends on the healing process and scar formation with limited replacement of lost tissue. Other organisms, such as spiny mice (Acomys cahirinus), salamanders, and zebrafish, are able to regenerate some damaged body components. Blastema is a tissue that is formed after tissue injury in such organisms and is composed of progenitor cells or dedifferentiated cells that differentiate into various cell types during regeneration. Thus, identifying the molecules responsible for initiation of blastema formation is an important aspect for understanding regeneration. Introns, a major source of noncoding RNAs (ncRNAs), have characteristic sizes in the axolotl, particularly in genes associated with development. These ncRNAs, particularly microRNAs (miRNAs), exhibit dynamic regulation during regeneration. These miRNAs play an essential role in timing and control of gene expression to order and organize processes necessary for blastema creation. Master keys or molecules that underlie the remarkable regenerative abilities of the axolotl remain to be fully explored and exploited. Further and ongoing research on regeneration promises new knowledge that may allow improved repair and renewal of human tissues.

Integrative Analysis of Axolotl Gene Expression Data from Regenerative and Wound Healing Limb Tissues

Axolotl (Ambystoma mexicanum) is a urodele amphibian endowed with remarkable regenerative capacities manifested in scarless wound healing and restoration of amputated limbs, which makes it a powerful experimental model for regenerative biology and medicine. Previous studies have utilized microarrays and RNA-Seq technologies for detecting differentially expressed (DE) genes in different phases of the axolotl limb regeneration. However, sufficient consistency may be lacking due to statistical limitations arising from intra-laboratory analyses. This study aims to bridge such gaps by performing an integrative analysis of publicly available microarray and RNA-Seq data from axolotl limb samples having comparable study designs using the "merging" method. A total of 351 genes were found DE in regenerative samples compared to the control in data of both technologies, showing an adjusted p-value < 0.01 and log fold change magnitudes >1. Downstream analyses illustrated consistent correlations of the directionality of DE genes within and between data of both technologies, as well as concordance with the literature on regeneration related biological processes. qRT-PCR analysis validated the observed expression level differences of five of the top DE genes. Future studies may benefit from the utilized concept and approach for enhanced statistical power and robust discovery of biomarkers of regeneration.

Comparative iTRAQ proteomics revealed proteins associated with lobed fin regeneration in Bichirs

Background

Polypterus senegalus can fully regenerate its pectoral lobed fins, including a complex endoskeleton, with remarkable precision. However, despite the enormous potential of this species for use in medical research, its regeneration mechanisms remain largely unknown.

Methods

To identify the differentially expressed proteins (DEPs) during the early stages of lobed fin regeneration in P. senegalus, we performed a differential proteomic analysis using isobaric tag for relative and absolute quantitation (iTRAQ) approach based quantitative proteome from the pectoral lobed fins at 3 time points. Furthermore, we validated the changes in protein expression with multiple-reaction monitoring (MRM) analysis.

Results

The experiment yielded a total of 3177 proteins and 15,091 unique peptides including 1006 non-redundant (nr) DEPs. Of these, 592 were upregulated while 349 were downregulated after lobed fin amputation when compared to the original tissue. Bioinformatics analyses showed that the DEPs were mainly associated with Ribosome and RNA transport, metabolic, ECM-receptor interaction, Golgi and endoplasmic reticulum, DNA replication, and Regulation of actin cytoskeleton.

Conclusions

To our knowledge, this is the first proteomic research to investigate alterations in protein levels and affected pathways in bichirs’ lobe-fin/limb regeneration. In addition, our study demonstrated a highly dynamic regulation during lobed fin regeneration in P. senegalus. These results not only provide a comprehensive dataset on differentially expressed proteins during the early stages of lobe-fin/limb regeneration but also advance our understanding of the molecular mechanisms underlying lobe-fin/limb regeneration.

Genomic Data Reveal Conserved Female Heterogamety in Giant Salamanders with Gigantic Nuclear Genomes

Systems of genetic sex determination and the homology of sex chromosomes in different taxa vary greatly across vertebrates. Much progress remains to be made in understanding systems of genetic sex determination in non-model organisms, especially those with homomorphic sex chromosomes and/or large genomes. We used reduced representation genome sequencing to investigate genetic sex determination systems in the salamander family Cryptobranchidae (genera Cryptobranchus and Andrias), which typifies both of these inherent difficulties. We tested hypotheses of male- or female-heterogamety by sequencing hundreds of thousands of anonymous genomic regions in a panel of known-sex cryptobranchids and characterized patterns of presence/absence, inferred zygosity, and depth of coverage to identify sex-linked regions of these 56 gigabase genomes. Our results strongly support the hypothesis that all cryptobranchid species possess homologous systems of female heterogamety, despite maintenance of homomorphic sex chromosomes over nearly 60 million years. Additionally, we report a robust, non-invasive genetic assay for sex diagnosis in Cryptobranchus and Andrias which may have great utility for conservation efforts with these endangered salamanders. Co-amplification of these W-linked markers in both cryptobranchid genera provides evidence for long-term sex chromosome stasis in one of the most divergent salamander lineages. These findings inform hypotheses about the ancestral mode of sex determination in salamanders, but suggest that comparative data from other salamander families are needed. Our results further demonstrate that massive genomes are not necessarily a barrier to effective genome-wide sequencing and that the resulting data can be highly informative about sex determination systems in taxa with homomorphic sex chromosomes.

Tail regeneration in Lepidosauria as an exception to the generalized lack of organ regeneration in amniotes

The present review hypothesizes that during the transition from water to land, amniotes lost part of the genetic program for metamorphosis utilized in larvae of their amphibian ancestors, a program that in extant fish and amphibians allows organ regeneration. The direct development of amniotes, with their growth from embryos to adults, occurred with the elimination of larval stages, increases the efficiency of immune responses and the complexity of nervous circuits. In amniotes, T-cells and macrophages likely eliminate embryonic-larval antigens that are replaced with the definitive antigens of adult organs. Among lepidosaurians numerous lizard families during the Permian and Triassic evolved the process of tail autotomy to escape predation, followed by tail regeneration. Autotomy limits inflammation allowing the formation of a regenerative blastema rich in the immunosuppressant and hygroscopic hyaluronic acid. Expression loss of developmental genes for metamorphosis and segmentation in addition to an effective immune system, determined an imperfect regeneration of the tail. Genes involved in somitogenesis were likely lost or are inactivated and the axial skeleton and muscles of the original tail are replaced with a nonsegmented cartilaginous tube and segmental myotomes. Lack of neural genes, negative influence of immune system, and isolation of the regenerating spinal cord within the cartilaginous tube impede the production of nerve and glial cells, and a stratified spinal cord with ganglia. Tissue and organ regeneration in other body regions of lizards and other reptiles is relatively limited, like in the other amniotes, although the cartilage shows a higher regenerative capability than in mammals.

A comprehensive reference transcriptome resource for the Iberian ribbed newt Pleurodeles waltl, an emerging model for developmental and regeneration biology

Urodele newts have unique biological properties, notably including prominent regeneration ability. The Iberian ribbed newt, Pleurodeles waltl, is a promising model amphibian distinguished by ease of breeding and efficient transgenic and genome editing methods. However, limited genetic information is available for P. waltl. We conducted an intensive transcriptome analysis of P. waltl using RNA-sequencing to build and annotate gene models. We generated 1.2 billion Illumina reads from a wide variety of samples across 12 different tissues/organs, unfertilized egg, and embryos at eight different developmental stages. These reads were assembled into 1,395,387 contigs, from which 202,788 non-redundant ORF models were constructed. The set is expected to cover a large fraction of P. waltl protein-coding genes, as confirmed by BUSCO analysis, where 98% of universal single-copy orthologs were identified. Ortholog analyses revealed the gene repertoire evolution of urodele amphibians. Using the gene set as a reference, gene network analysis identified regeneration-, developmental-stage-, and tissue-specific co-expressed gene modules. Our transcriptome resource is expected to enhance future research employing this emerging model animal for regeneration research as well as for investigations in other areas including developmental biology, stem cell biology, and cancer research. These data are available via our portal website, iNewt (http://www.nibb.ac.jp/imori/main/).

A chromosome-scale assembly of the axolotl genome

The axolotl (Ambystoma mexicanum) provides critical models for studying regeneration, evolution, and development. However, its large genome (∼32 Gb) presents a formidable barrier to genetic analyses. Recent efforts have yielded genome assemblies consisting of thousands of unordered scaffolds that resolve gene structures, but do not yet permit large-scale analyses of genome structure and function. We adapted an established mapping approach to leverage dense SNP typing information and for the first time assemble the axolotl genome into 14 chromosomes. Moreover, we used fluorescence in situ hybridization to verify the structure of these 14 scaffolds and assign each to its corresponding physical chromosome. This new assembly covers 27.3 Gb and encompasses 94% of annotated gene models on chromosomal scaffolds. We show the assembly's utility by resolving genome-wide orthologies between the axolotl and other vertebrates, identifying the footprints of historical introgression events that occurred during the development of axolotl genetic stocks, and precisely mapping several phenotypes including a large deletion underlying the cardiac mutant. This chromosome-scale assembly will greatly facilitate studies of the axolotl in biological research.

Using transcriptomics to enable a plethodontid salamander (Bolitoglossa ramosi) for limb regeneration research

Background

Tissue regeneration is widely distributed across the tree of life. Among vertebrates, salamanders possess an exceptional ability to regenerate amputated limbs and other complex structures. Thus far, molecular insights about limb regeneration have come from a relatively limited number of species from two closely related salamander families. To gain a broader perspective on the molecular basis of limb regeneration and enhance the molecular toolkit of an emerging plethodontid salamander (Bolitoglossa ramosi), we used RNA-Seq to generate a de novo reference transcriptome and identify differentially expressed genes during limb regeneration.

Results

Using paired-end Illumina sequencing technology and Trinity assembly, a total of 433,809 transcripts were recovered and we obtained functional annotation for 142,926 non-redundant transcripts of the B. ramosi de novo reference transcriptome. Among the annotated transcripts, 602 genes were identified as differentially expressed during limb regeneration. This list was further processed to identify a core set of genes that exhibit conserved expression changes between B. ramosi and the Mexican axolotl (Ambystoma mexicanum), and presumably their common ancestor from approximately 180 million years ago.

Conclusions

We identified genes from B. ramosi that are differentially expressed during limb regeneration, including multiple conserved protein-coding genes and possible putative species-specific genes. Comparative analyses reveal a subset of genes that show similar patterns of expression with ambystomatid species, which highlights the importance of developing comparative gene expression data for studies of limb regeneration among salamanders.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5076-0) contains supplementary material, which is available to authorized users.

Macroevolutionary shift in the size of amphibian genomes and the role of life history and climate

The evolution and great diversity of genome size has been of long-standing interest to biologists, but has seldom been investigated on a broad phylogenetic scale. Here we present a comparative quantitative analysis of factors shaping genome size evolution in amphibians, the extant class of vertebrates with the largest variation in genome size. We find that amphibian genomes have undergone saltations in size, although these are rare and the evolutionary history of genome size in amphibians has otherwise been one of gradual, time-dependent variation (that is, Brownian motion). This macroevolutionary homogeneity is remarkable given the evolutionary and ecological diversity of most other aspects of the natural history of amphibians. Contrary to previous claims, we find no evidence for associations between life cycle complexity and genome size despite the high diversity of reproductive modes and the multiple events of independent evolution of divergent life cycles in the group. Climate (temperature and humidity) affects genome size indirectly, at least in frogs, as a consequence of its effect on premetamorphic developmental period, although directionality of the relationship between developmental period and genome size is not unequivocal.

Forever young: Endocrinology of paedomorphosis in the Mexican axolotl (Ambystoma mexicanum)

The Mexican axolotl (Ambystoma mexicanum) is a salamander species that does not undergo metamorphosis, resulting in the retention of juvenile characteristics in the mature breeding stage (paedomorphosis). Here we review the endocrinological studies investigating the proximate cause of axolotl paedomorphosis with a focus on the hypothalamo-pituitary-thyroid (HPT) axis. It is well established that axolotl paedomorphosis is a consequence of low activity of the HPT axis. The pituitary hormone thyrotropin (TSH) is capable of inducing metamorphosis in the axolotl, which indicates that all processes and interactions in the HPT axis below the pituitary level are functional, but that TSH release is impaired. In metamorphosing species, TSH secretion is largely controlled by the hypothalamic neuropeptide corticotropin-releasing hormone (CRH), which seems to have lost its thyrotropic activity in the axolotl. However, preliminary experiments have not yet confirmed a role for faulty CRH signalling in axolotl paedomorphosis. Other hypothalamic factors and potential pituitary inhibitors need to be investigated to identify their roles in amphibian metamorphosis and axolotl paedomorphosis.

Causes and evolutionary consequences of primordial germ-cell specification mode in metazoans

In animals, primordial germ cells (PGCs) give rise to the germ lines, the cell lineages that produce sperm and eggs. PGCs form in embryogenesis, typically by one of two modes: a likely ancestral mode wherein germ cells are induced during embryogenesis by cell-cell signaling (induction) or a derived mechanism whereby germ cells are specified by using germ plasm-that is, maternally specified germ-line determinants (inheritance). The causes of the shift to germ plasm for PGC specification in some animal clades remain largely unknown, but its repeated convergent evolution raises the question of whether it may result from or confer an innate selective advantage. It has been hypothesized that the acquisition of germ plasm confers enhanced evolvability, resulting from the release of selective constraint on somatic gene networks in embryogenesis, thus leading to acceleration of an organism's protein-sequence evolution, particularly for genes expressed at early developmental stages, and resulting in high speciation rates in germ plasm-containing lineages (denoted herein as the "PGC-specification hypothesis"). Although that hypothesis, if supported, could have major implications for animal evolution, our recent large-scale coding-sequence analyses from vertebrates and invertebrates provided important examples of genera that do not support the hypothesis of liberated constraint under germ plasm. Here, we consider reasons why germ plasm might be neither a direct target of selection nor causally linked to accelerated animal evolution. We explore alternate scenarios that could explain the repeated evolution of germ plasm and propose potential consequences of the inheritance and induction modes to animal evolutionary biology.

Drug delivery and epimorphic salamander-type mouse regeneration: A full parts and labor plan

The capacity to regenerate entire body parts, tissues, and organs had generally been thought to be lost in evolution with very few exceptions (e.g. the liver) surviving in mammals. The discovery of the MRL mouse and the elucidation of the underlying molecular pathway centering around hypoxia inducible factor, HIF-1α, has allowed a drug and materials approach to regeneration in mice and hopefully humans. The HIF-1α pathway is ancient and permitted the transition from unicellular to multicellular organisms. Furthermore, HIF-1α and its regulation by PHDs, important oxygen sensors in the cell, provides a perfect drug target. We review the historical background of regeneration biology, the discovery of the MRL mouse, and its underlying biology, and novel approaches to drugs, targets, and delivery systems (see Fig. 1).

The axolotl genome and the evolution of key tissue formation regulators

Salamanders serve as important tetrapod models for developmental, regeneration and evolutionary studies. An extensive molecular toolkit makes the Mexican axolotl (Ambystoma mexicanum) a key representative salamander for molecular investigations. Here we report the sequencing and assembly of the 32-gigabase-pair axolotl genome using an approach that combined long-read sequencing, optical mapping and development of a new genome assembler (MARVEL). We observed a size expansion of introns and intergenic regions, largely attributable to multiplication of long terminal repeat retroelements. We provide evidence that intron size in developmental genes is under constraint and that species-restricted genes may contribute to limb regeneration. The axolotl genome assembly does not contain the essential developmental gene Pax3. However, mutation of the axolotl Pax3 paralogue Pax7 resulted in an axolotl phenotype that was similar to those seen in Pax3^-/- and Pax7^-/- mutant mice. The axolotl genome provides a rich biological resource for developmental and evolutionary studies.

Virtual Genome Walking across the 32 Gb Ambystoma mexicanum genome; assembling gene models and intronic sequence

Large repeat rich genomes present challenges for assembly using short read technologies. The 32 Gb axolotl genome is estimated to contain ~19 Gb of repetitive DNA making an assembly from short reads alone effectively impossible. Indeed, this model species has been sequenced to 20× coverage but the reads could not be conventionally assembled. Using an alternative strategy, we have assembled subsets of these reads into scaffolds describing over 19,000 gene models. We call this method Virtual Genome Walking as it locally assembles whole genome reads based on a reference transcriptome, identifying exons and iteratively extending them into surrounding genomic sequence. These assemblies are then linked and refined to generate gene models including upstream and downstream genomic, and intronic, sequence. Our assemblies are validated by comparison with previously published axolotl bacterial artificial chromosome (BAC) sequences. Our analyses of axolotl intron length, intron-exon structure, repeat content and synteny provide novel insights into the genic structure of this model species. This resource will enable new experimental approaches in axolotl, such as ChIP-Seq and CRISPR and aid in future whole genome sequencing efforts. The assembled sequences and annotations presented here are freely available for download from https://tinyurl.com/y8gydc6n . The software pipeline is available from https://github.com/LooseLab/iterassemble .

A Tissue-Mapped Axolotl De Novo Transcriptome Enables Identification of Limb Regeneration Factors

Mammals have extremely limited regenerative capabilities; however, axolotls are profoundly regenerative and can replace entire limbs. The mechanisms underlying limb regeneration remain poorly understood, partly because the enormous and incompletely sequenced genomes of axolotls have hindered the study of genes facilitating regeneration. We assembled and annotated a de novo transcriptome using RNA-sequencing profiles for a broad spectrum of tissues that is estimated to have near-complete sequence information for 88% of axolotl genes. We devised expression analyses that identified the axolotl orthologs of cirbp and kazald1 as highly expressed and enriched in blastemas. Using morpholino anti-sense oligonucleotides, we find evidence that cirbp plays a cytoprotective role during limb regeneration whereas manipulation of kazald1 expression disrupts regeneration. Our transcriptome and annotation resources greatly complement previous transcriptomic studies and will be a valuable resource for future research in regenerative biology.

Amphibian and Avian Karyotype Evolution: Insights from Lampbrush Chromosome Studies

Amphibian and bird karyotypes typically have a complex organization, which makes them difficult for standard cytogenetic analysis. That is, amphibian chromosomes are generally large, enriched with repetitive elements, and characterized by the absence of informative banding patterns. The majority of avian karyotypes comprise a small number of relatively large macrochromosomes and numerous tiny morphologically undistinguishable microchromosomes. A good progress in investigation of amphibian and avian chromosome evolution became possible with the usage of giant lampbrush chromosomes typical for growing oocytes. Due to the giant size, peculiarities of organization and enrichment with cytological markers, lampbrush chromosomes can serve as an opportune model for comprehensive high-resolution cytogenetic and cytological investigations. Here, we review the main findings on chromosome evolution in amphibians and birds that were obtained using lampbrush chromosomes. In particular, we discuss the data on evolutionary chromosomal rearrangements, accumulation of polymorphisms, evolution of sex chromosomes as well as chromosomal changes during clonal reproduction of interspecies hybrids.

Advances in Decoding Axolotl Limb Regeneration

Humans and other mammals are limited in their natural abilities to regenerate lost body parts. By contrast, many salamanders are highly regenerative and can spontaneously replace lost limbs even as adults. Because salamander limbs are anatomically similar to human limbs, knowing how they regenerate should provide important clues for regenerative medicine. Although interest in understanding the mechanics of this process has never wavered, until recently researchers have been vexed by seemingly impenetrable logistics of working with these creatures at a molecular level. Chief among the problems has been the very large size of salamander genomes, and not a single salamander genome has been fully sequenced to date. Recently the enormous gap in sequence information has been bridged by approaches that leverage mRNA as the starting point. Together with functional experimentation, these data are rapidly enabling researchers to finally uncover the molecular mechanisms underpinning the astonishing biological process of limb regeneration.

A tongue for all seasons: extreme phenotypic flexibility in salamandrid newts

Many organisms faced with seasonally fluctuating abiotic and biotic conditions respond by altering their phenotype to account for the demands of environmental changes. Here we discovered that newts, which switch seasonally between an aquatic and terrestrial lifestyle, grow a complex adhesive system on their tongue pad consisting of slender lingual papillae and mucus-producing cells to increase the efficiency of prey capture as they move from water onto land. The adhesive system is reduced again as newts switch back to their aquatic stage, where they use suction to capture prey. As suction performance is also enhanced seasonally by reshaping of the mouth due to the growth of labial lobes, our results show that newts are exceptional in exhibiting phenotypic flexibility in two alternating components (i.e. tongue pad and labial lobes) within a single functional system, and suggest that this form of phenotypic flexibility demands complex genetic regulation.

Analysis of embryonic development in the unsequenced axolotl: Waves of transcriptomic upheaval and stability

The axolotl (Ambystoma mexicanum) has long been the subject of biological research, primarily owing to its outstanding regenerative capabilities. However, the gene expression programs governing its embryonic development are particularly underexplored, especially when compared to other amphibian model species. Therefore, we performed whole transcriptome polyA+ RNA sequencing experiments on 17 stages of embryonic development. As the axolotl genome is unsequenced and its gene annotation is incomplete, we built de novo transcriptome assemblies for each stage and garnered functional annotation by comparing expressed contigs with known genes in other organisms. In evaluating the number of differentially expressed genes over time, we identify three waves of substantial transcriptome upheaval each followed by a period of relative transcriptome stability. The first wave of upheaval is between the one and two cell stage. We show that the number of differentially expressed genes per unit time is higher between the one and two cell stage than it is across the mid-blastula transition (MBT), the period of zygotic genome activation. We use total RNA sequencing to demonstrate that the vast majority of genes with increasing polyA+ signal between the one and two cell stage result from polyadenylation rather than de novo transcription. The first stable phase begins after the two cell stage and continues until the mid-blastula transition, corresponding with the pre-MBT phase of transcriptional quiescence in amphibian development. Following this is a peak of differential gene expression corresponding with the activation of the zygotic genome and a phase of transcriptomic stability from stages 9-11. We observe a third wave of transcriptomic change between stages 11 and 14, followed by a final stable period. The last two stable phases have not been documented in amphibians previously and correspond to times of major morphogenic change in the axolotl embryo: gastrulation and neurulation. These results yield new insights into global gene expression during early stages of amphibian embryogenesis and will help to further develop the axolotl as a model species for developmental and regenerative biology.

A linkage map for the Newt Notophthalmus viridescens: Insights in vertebrate genome and chromosome evolution

Genetic linkage maps are fundamental resources that enable diverse genetic and genomic approaches, including quantitative trait locus (QTL) analyses and comparative studies of genome evolution. It is straightforward to build linkage maps for species that are amenable to laboratory culture and genetic crossing designs, and that have relatively small genomes and few chromosomes. It is more difficult to generate linkage maps for species that do not meet these criteria. Here, we introduce a method to rapidly build linkage maps for salamanders, which are known for their enormous genome sizes. As proof of principle, we developed a linkage map with thousands of molecular markers (N=2349) for the Eastern newt (Notophthalmus viridescens). The map contains 12 linkage groups (152.3-934.7cM), only one more than the number of chromosome pairs. Importantly, this map was generated using RNA isolated from a single wild caught female and her 28 offspring. We used the map to reveal chromosome-scale conservation of synteny among N. viridescens, A. mexicanum (Urodela), and chicken (Amniota), and to identify large conserved segments between N. viridescens and Xenopus tropicalis (Anura). We also show that met1, a major effect QTL that regulates the expression of alternate metamorphic and paedomorphic modes of development in Ambystoma, associates with a chromosomal fusion that is not found in the N. viridescens map. Our results shed new light on the ancestral amphibian karyotype and reveal specific fusion and translocation events that shaped the genomes of three amphibian model taxa. The ability to rapidly build linkage maps for large salamander genomes will enable genetic and genomic analyses within this important vertebrate group, and more generally, empower comparative studies of vertebrate biology and evolution.

The relationship between growth and pattern formation

Successful development depends on the creation of spatial gradients of transcription factors within developing fields, and images of graded distributions of gene products populate the pages of developmental biology journals. Therefore the challenge is to understand how the graded levels of intracellular transcription factors are generated across fields of cells. We propose that transcription factor gradients are generated as a result of an underlying gradient of cell cycle lengths. Very long cell cycles will permit accumulation of a high level of a gene product encoded by a large transcription unit, whereas shorter cell cycles will permit progressively fewer transcripts to be completed due to gating of transcription by the cell cycle. We also propose that the gradients of cell cycle lengths are generated by gradients of extracellular morphogens/growth factors. The model of cell cycle gated transcriptional regulation brings focus back to the functional role of morphogens as cell cycle regulators, and proposes a specific and testable mechanism by which morphogens, in their roles as growth factors (how they were originally discovered), also determine cell fate.

Transposons, Genome Size, and Evolutionary Insights in Animals

The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.

Initial characterization of the large genome of the salamander Ambystoma mexicanum using shotgun and laser capture chromosome sequencing

Vertebrates exhibit substantial diversity in genome size, and some of the largest genomes exist in species that uniquely inform diverse areas of basic and biomedical research. For example, the salamander Ambystoma mexicanum (the Mexican axolotl) is a model organism for studies of regeneration, development and genome evolution, yet its genome is ~10× larger than the human genome. As part of a hierarchical approach toward improving genome resources for the species, we generated 600 Gb of shotgun sequence data and developed methods for sequencing individual laser-captured chromosomes. Based on these data, we estimate that the A. mexicanum genome is ~32 Gb. Notably, as much as 19 Gb of the A. mexicanum genome can potentially be considered single copy, which presumably reflects the evolutionary diversification of mobile elements that accumulated during an ancient episode of genome expansion. Chromosome-targeted sequencing permitted the development of assemblies within the constraints of modern computational platforms, allowed us to place 2062 genes on the two smallest A. mexicanum chromosomes and resolves key events in the history of vertebrate genome evolution. Our analyses show that the capture and sequencing of individual chromosomes is likely to provide valuable information for the systematic sequencing, assembly and scaffolding of large genomes.

What's in a genome? The C-value enigma and the evolution of eukaryotic genome content

Some notable exceptions aside, eukaryotic genomes are distinguished from those of Bacteria and Archaea in a number of ways, including chromosome structure and number, repetitive DNA content, and the presence of introns in protein-coding regions. One of the most notable differences between eukaryotic and prokaryotic genomes is in size. Unlike their prokaryotic counterparts, eukaryotes exhibit enormous (more than 60,000-fold) variability in genome size which is not explained by differences in gene number. Genome size is known to correlate with cell size and division rate, and by extension with numerous organism-level traits such as metabolism, developmental rate or body size. Less well described are the relationships between genome size and other properties of the genome, such as gene content, transposable element content, base pair composition and related features. The rapid expansion of 'complete' genome sequencing projects has, for the first time, made it possible to examine these relationships across a wide range of eukaryotes in order to shed new light on the causes and correlates of genome size diversity. This study presents the results of phylogenetically informed comparisons of genome data for more than 500 species of eukaryotes. Several relationships are described between genome size and other genomic parameters, and some recommendations are presented for how these insights can be extended even more broadly in the future.

Co-option and evolution of non-olfactory proteinaceous pheromones in a terrestrial lungless salamander

Gene co-option is a major force in the evolution of novel biological functions. In plethodontid salamanders, males deliver proteinaceous courtship pheromones to the female olfactory system or transdermally to the bloodstream. Molecular studies identified three families of highly duplicated, rapidly evolving pheromones (PRF, PMF, and SPF). Analyses for Plethodon salamanders revealed pheromone mixtures of primarily PRF and PMF. The current study demonstrates that in Desmognathus ocoee--a plesiomorphic species with transdermal delivery--SPF is the major pheromone component representing >30% of total protein. Chromatographic profiles of D. ocoee pheromones were consistent from May through October. LC/MS-MS analysis suggested uniform SPF isoform expression between individual male D. ocoee. A gene ancestry for SPF with the Three-Finger Protein superfamily was supported by intron-exon boundaries, but not by the disulfide bonding pattern. Further analysis of the pheromone mixture revealed paralogs to peptide hormones that contained mutations in receptor binding regions, such that these novel molecules may alter female physiology by acting as hormone agonists/antagonists. Cumulatively, gene co-option, duplication, and neofunctionalization have permitted recruitment of additional gene families for pheromone activity. Such independent co-option events may be playing a key role in salamander speciation by altering male traits that influence reproductive success.

Hox cluster characterization of Banna caecilian (Ichthyophis bannanicus) provides hints for slow evolution of its genome

BACKGROUND

Caecilians, with a discrete lifestyle, are the least explored group of amphibians. Though with distinct traits, many aspects of their biology are poorly investigated. Obtaining the caecilian genomic sequences will offer new perspectives and aid the fundamental studies in caecilian biology. The caecilian genomic sequences are also important and practical in the comparative genomics of amphibians. Currently, however, only sparse genomic sequences of caecilians are available. Hox genes, an old family of transcription factors playing central roles in the establishment of metazoan body plan. Understanding their structure and genomic organization may provide insights into the animal's genome, which is valuable for animals without a sequenced genome.

RESULTS

We sequenced and characterized the Hox clusters of Banna caecilian (Ichthyophis bannanicus) with a strategy combining long range PCR and genome walking. We obtained the majority of the four caecilian Hox clusters and identified 39 Hox genes, 5 microRNA genes and 1 pseudogene (ψHoxD12). There remained seven intergenic gaps we were unable to fill. From the obtained sequences, the caecilian Hox clusters contained less repetitive sequences and more conserved noncoding elements (CNEs) than the frog counterparts. We found that caecilian and coelacanth shared many more CNEs than frog and coelacanth did. Relative rate of sequence evolution showed that caecilian Hox genes evolved significantly more slowly than the other tetrapod species used in this study and were comparable to the slowly evolving coelacanth Hox genes. Phylogenetic tree of the four Hox clusters also revealed shorter branch length especially for the caecilian HoxA, HoxB and HoxD clusters. These features of the caecilian Hox clusters suggested a slowly evolving genome, which was supported by further analysis of a large orthologous protein dataset.

CONCLUSIONS

Our analyses greatly extended the knowledge about the caecilian Hox clusters from previous PCR surveys. From the obtained Hox sequences and the orthologous protein dataset, the caecilian Hox loci and its genome appear evolving comparatively slowly. As the basal lineage of amphibians and land vertebrate, this characteristic of the caecilian genome is valuable in the study concerning the genome biology and evolution of amphibians and early tetrapods.

Hellbender genome sequences shed light on genomic expansion at the base of crown salamanders

Among animals, genome sizes range from 20 Mb to 130 Gb, with 380-fold variation across vertebrates. Most of the largest vertebrate genomes are found in salamanders, an amphibian clade of 660 species. Thus, salamanders are an important system for studying causes and consequences of genomic gigantism. Previously, we showed that plethodontid salamander genomes accumulate higher levels of long terminal repeat (LTR) retrotransposons than do other vertebrates, although the evolutionary origins of such sequences remained unexplored. We also showed that some salamanders in the family Plethodontidae have relatively slow rates of DNA loss through small insertions and deletions. Here, we present new data from Cryptobranchus alleganiensis, the hellbender. Cryptobranchus and Plethodontidae span the basal phylogenetic split within salamanders; thus, analyses incorporating these taxa can shed light on the genome of the ancestral crown salamander lineage, which underwent expansion. We show that high levels of LTR retrotransposons likely characterize all crown salamanders, suggesting that disproportionate expansion of this transposable element (TE) class contributed to genomic expansion. Phylogenetic and age distribution analyses of salamander LTR retrotransposons indicate that salamanders' high TE levels reflect persistence and diversification of ancestral TEs rather than horizontal transfer events. Finally, we show that relatively slow DNA loss rates through small indels likely characterize all crown salamanders, suggesting that a decreased DNA loss rate contributed to genomic expansion at the clade's base. Our identification of shared genomic features across phylogenetically distant salamanders is a first step toward identifying the evolutionary processes underlying accumulation and persistence of high levels of repetitive sequence in salamander genomes.

Long bone histology of the stem salamander Kokartus honorarius (Amphibia: Caudata) from the Middle Jurassic of Kyrgyzstan

Kokartus honorarius from the Middle Jurassic (Bathonian) of Kyrgyzstan is one of the oldest salamanders in the fossil record, characterized by a mixture of plesiomorphic morphological features and characters shared with crown-group salamanders. Here we present a detailed histological analysis of its long bones. The analysis of a growth series demonstrates a significant histological maturation during ontogeny, expressed by the progressive appearance of longitudinally oriented primary vascular canals, primary osteons, growth marks, remodelling features in primary bone tissues, as well as progressive resorption of the calcified cartilage, formation of endochondral bone and development of cartilaginous to bony trabeculae in the epiphyses. Apart from the presence of secondary osteons, the long bone histology of Kokartus is very similar to that of miniaturized temnospondyls, other Jurassic stem salamanders, miniaturized seymouriamorphs and modern crown-group salamanders. We propose that the presence of secondary osteons in Kokartus honorarius is a plesiomorphic feature, and the loss of secondary osteons in the long bones of crown-group salamanders as well as in those of miniaturized temnospondyls is the result of miniaturization processes. Hitherto, all stem salamander long bong histology (Kokartus, Marmorerpeton and 'salamander A') has been generally described as having paedomorphic features (i.e. the presence of Katschenko's Line and a layer of calcified cartilage), these taxa were thus most likely neotenic forms. The absence of clear lines of arrested growth and annuli in long bones of Kokartus honorarius suggests that the animals lived in an environment with stable local conditions.

A whole-genome shotgun approach for assembling and anchoring the hexaploid bread wheat genome

Polyploid species have long been thought to be recalcitrant to whole-genome assembly. By combining high-throughput sequencing, recent developments in parallel computing, and genetic mapping, we derive, de novo, a sequence assembly representing 9.1 Gbp of the highly repetitive 16 Gbp genome of hexaploid wheat, Triticum aestivum, and assign 7.1 Gb of this assembly to chromosomal locations. The genome representation and accuracy of our assembly is comparable or even exceeds that of a chromosome-by-chromosome shotgun assembly. Our assembly and mapping strategy uses only short read sequencing technology and is applicable to any species where it is possible to construct a mapping population.

A whole-genome shotgun approach for assembling and anchoring the hexaploid bread wheat genome

Polyploid species have long been thought to be recalcitrant to whole-genome assembly. By combining high-throughput sequencing, recent developments in parallel computing, and genetic mapping, we derive, de novo, a sequence assembly representing 9.1 Gbp of the highly repetitive 16 Gbp genome of hexaploid wheat, Triticum aestivum, and assign 7.1 Gb of this assembly to chromosomal locations. The genome representation and accuracy of our assembly is comparable or even exceeds that of a chromosome-by-chromosome shotgun assembly. Our assembly and mapping strategy uses only short read sequencing technology and is applicable to any species where it is possible to construct a mapping population.

Low levels of LTR retrotransposon deletion by ectopic recombination in the gigantic genomes of salamanders

Across the tree of life, species vary dramatically in nuclear genome size. Mutations that add or remove sequences from genomes-insertions or deletions, or indels-are the ultimate source of this variation. Differences in the tempo and mode of insertion and deletion across taxa have been proposed to contribute to evolutionary diversity in genome size. Among vertebrates, most of the largest genomes are found within the salamanders, an amphibian clade with genome sizes ranging from ~14 to ~120 Gb. Salamander genomes have been shown to experience slower rates of DNA loss through small (i.e., <30 bp) deletions than do other vertebrate genomes. However, no studies have addressed DNA loss from salamander genomes resulting from larger deletions. Here, we focus on one type of large deletion-ectopic-recombination-mediated removal of LTR retrotransposon sequences. In ectopic recombination, double-strand breaks are repaired using a "wrong" (i.e., ectopic, or non-allelic) template sequence-typically another locus of similar sequence. When breaks occur within the LTR portions of LTR retrotransposons, ectopic-recombination-mediated repair can produce deletions that remove the internal transposon sequence and the equivalent of one of the two LTR sequences. These deletions leave a signature in the genome-a solo LTR sequence. We compared levels of solo LTRs in the genomes of four salamander species with levels present in five vertebrates with smaller genomes. Our results demonstrate that salamanders have low levels of solo LTRs, suggesting that ectopic-recombination-mediated deletion of LTR retrotransposons occurs more slowly than in other vertebrates with smaller genomes.

Sal-Site: research resources for the Mexican axolotl

Sal-Site serves axolotl research efforts by providing Web access to genomic data and information, and living stocks that are reared and made available by the Ambystoma Genetic Stock Center (AGSC). In this chapter, we detail how investigators can search for genes of interest among Sal-Site resources to identify orthologous nucleotide and protein-coding sequences, determine genome positions within the Ambystoma meiotic map, and obtain estimates of gene expression. In the near future, additional genomic resources will be made available for the axolotl, including a listing of genes that are partially or wholly contained within Bacterial Artificial Chromosome (BAC) vectors, a prioritized collection of deeply sequenced BAC clones, chromosome-specific assemblies of genomic DNA, and transgenic axolotls that are engineered using TALENs and CRISPRs. Also, services provided by the AGSC will be expanded to include microinjection of user constructs into single cell embryos and distribution of axolotl tissues, DNA, and RNA. In conclusion, Sal-Site is a useful resource that generates, shares, and evolves Ambystoma associated information and databases to serve research and education.

Generating and identifying axolotls with targeted mutations using Cas9 RNA-guided nuclease

The CRISPR/Cas9 RNA-guided nuclease now enables a reverse genetics approach to investigate the function of genes of interest during regeneration in the axolotl. The process of generating the constructs necessary for targeting a gene of interest is considerably less labor intensive than for other methods of targeted mutagenesis such as Zinc finger nucleases or Transcription activator-like effector nucleases. Here, we describe the identification of targetable sequences in the gene of interest, the construction of unique guide RNAs, the microinjection of these RNAs with Cas9-encoding mRNA, the selection of well-injected animals, and an inexpensive, PCR-based method for identifying highly mutagenized animals.

Highly efficient targeted mutagenesis in axolotl using Cas9 RNA-guided nuclease

Among tetrapods, only urodele salamanders, such as the axolotl Ambystoma mexicanum, can completely regenerate limbs as adults. The mystery of why salamanders, but not other animals, possess this ability has for generations captivated scientists seeking to induce this phenomenon in other vertebrates. Although many recent advances in molecular biology have allowed limb regeneration and tissue repair in the axolotl to be investigated in increasing detail, the molecular toolkit for the study of this process has been limited. Here, we report that the CRISPR-Cas9 RNA-guided nuclease system can efficiently create mutations at targeted sites within the axolotl genome. We identify individual animals treated with RNA-guided nucleases that have mutation frequencies close to 100% at targeted sites. We employ this technique to completely functionally ablate EGFP expression in transgenic animals and recapitulate developmental phenotypes produced by loss of the conserved gene brachyury. Thus, this advance allows a reverse genetic approach in the axolotl and will undoubtedly provide invaluable insight into the mechanisms of salamanders' unique regenerative ability.

Development, validation and high-throughput analysis of sequence markers in nonmodel species

DNA sequences derived from multiple regions of the nuclear genome are essential for historical inferences in the fields of phylogeography and phylogenetics. The appropriate markers should be single-copy, variable, easy to amplify from multiple samples and easy to sequence using high-throughput technologies. This may be difficult to achieve for species lacking sequenced genomes and particularly challenging for species possessing large genomes, which consist mostly of repetitive sequences. Here, we present a cost-effective, broadly applicable framework for designing, validating and high-throughput sequencing of multiple markers in nonmodel species without sequenced genomes. We demonstrate its utility in two closely related species of newts, representatives of urodeles, a vertebrate group characterized by large genomes. We show that over 80 markers, c. 600 bp each, developed mainly from 3' untranslated transcript regions (3'UTR) may be effectively multiplexed and sequenced. Data are further processed using standard, freely available bioinformatic tools, producing phase-resolved sequences. The approach does not require barcoded PCR primers, and the cost of library preparation is independent of the number of markers investigated. We hope that this approach will be of broad interest for researchers working at the interface of population genetics and phylogenetics, exploring deep intraspecific genetic structure, species boundaries and phylogeographies of closely related species.

Microsatellite analyses across three diverse vertebrate transcriptomes (Acipenser fulvescens, Ambystoma tigrinum, and Dipodomys spectabilis)

Historically, many population genetics studies have utilized microsatellite markers sampled at random from the genome and presumed to be selectively neutral. Recent studies, however, have shown that microsatellites can occur in transcribed regions, where they are more likely to be under selection. In this study, we mined microsatellites from transcriptomes generated by 454-pyrosequencing for three vertebrate species: lake sturgeon (Acipenser fulvescens), tiger salamander (Ambystoma tigrinum), and kangaroo rat (Dipodomys spectabilis). We evaluated (i) the occurrence of microsatellites across species; (ii) whether particular gene ontology terms were over-represented in genes that contained microsatellites; (iii) whether repeat motifs were located in untranslated regions or coding sequences of genes; and (iv) in silico polymorphism. Microsatellites were less common in tiger salamanders than in either lake sturgeon or kangaroo rats. Across libraries, trinucleotides were found more frequently than any other motif type, presumably because they do not cause frameshift mutations. By evaluating variation across reads assembled to a given contig, we were able to identify repeat motifs likely to be polymorphic. Our study represents one of the first comparative data sets on the distribution of vertebrate microsatellites within expressed genes. Our results reinforce the idea that microsatellites do not always occur in noncoding DNA, but commonly occur in expressed genes.