A versatile and highly efficient toolkit including 102 nuclear markers for vertebrate phylogenomics, tested by resolving the higher level relationships of the caudata

The Influence of Genome and Cell Size on Brain Morphology in Amphibians

In amphibians, nerve cell size is highly correlated with genome size, and increases in genome and cell size cause a retardation of the rate of development of nervous (as well as nonnervous) tissue leading to secondary simplification. This yields an inverse relationship between genome and cell size on the one hand and morphological complexity of the tectum mesencephali as the main visual center, the size of the torus semicircularis as the main auditory center, the size of the amphibian papilla as an important peripheral auditory structure, and the size of the cerebellum as a major sensorimotor center. Nervous structures developing later (e.g., torus and cerebellum) are more affected by secondary simplification than those that develop earlier (e.g., the tectum). This effect is more prominent in salamanders and caecilians than in frogs owing to larger genome and cells sizes in the former two taxa. We hypothesize that because of intragenomic evolutionary processes, important differences in brain morphology can arise independently of specific environmental selection.

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.

Efficient Detection of Novel Nuclear Markers for Brassicaceae by Transcriptome Sequencing

The lack of DNA sequence information for most non-model organisms impairs the design of primers that are universally applicable for the study of molecular polymorphisms in nuclear markers. Next-generation sequencing (NGS) techniques nowadays provide a powerful approach to overcome this limitation. We present a flexible and inexpensive method to identify large numbers of nuclear primer pairs that amplify in most Brassicaceae species. We first obtained and mapped NGS transcriptome sequencing reads from two of the distantly related Brassicaceae species, Cardamine hirsuta and Arabis alpina, onto the Arabidopsis thaliana reference genome, and then identified short conserved sequence motifs among the three species bioinformatically. From these, primer pairs to amplify coding regions (nuclear protein coding loci, NPCL) and exon-primed intron-crossing sequences (EPIC) were developed. We identified 2,334 universally applicable primer pairs, targeting 1,164 genes, which provide a large pool of markers as readily usable genomic resource that will help addressing novel questions in the Brassicaceae family. Testing a subset of the newly designed nuclear primer pairs revealed that a great majority yielded a single amplicon in all of the 30 investigated Brassicaceae taxa. Sequence analysis and phylogenetic reconstruction with a subset of these markers on different levels of phylogenetic divergence in the mustard family were compared with previous studies. The results corroborate the usefulness of the newly developed primer pairs, e.g., for phylogenetic analyses or population genetic studies. Thus, our method provides a cost-effective approach for designing nuclear loci across a broad range of taxa and is compatible with current NGS technologies.

Parallel tagged amplicon sequencing of relatively long PCR products using the Illumina HiSeq platform and transcriptome assembly

In phylogenetics and population genetics, a large number of loci are often needed to accurately resolve species relationships. Normally, loci are enriched by PCR and sequenced by Sanger sequencing, which is expensive when the number of amplicons is large. Next-generation sequencing (NGS) techniques are increasingly used for parallel amplicon sequencing, which reduces sequencing costs tremendously, but has not reduced preparation costs very much. Moreover, for most current NGS methods, amplicons need to be purified and quantified before sequencing and their lengths are also restricted (normally <700 bp). Here, we describe an approach to sequence pooled amplicons of any length using the Illumina platform. Using this method, amplicons are pooled at equal volume rather than at equal concentration, thus eliminating the laborious purification and quantification steps. We then shear the pooled amplicons, repair the ends, add sample identifying linkers and pool multiple samples prior to Illumina library preparation. Data are then assembled using the transcriptome assembly program trinity, which is optimized to deal with templates of highly varying quantities. We demonstrated the utility of our approach by recovering 93.5% of the target amplicons (size up to 1650 bp) in full length for a 16 taxa × 101 loci project, using ~2.0 GB of Illumina HiSeq paired-end 90-bp data. Overall, we validate a rapid, cost-effective and scalable approach to sequence a large number of targeted loci from a large number of samples that is particularly suitable for both phylogenetics and population genetics studies that require a modest scale of data.

Identification of the orphan gene Prod 1 in basal and other salamander families

BACKGROUND

The urodele amphibians (salamanders) are the only adult tetrapods able to regenerate the limb. It is unclear if this is an ancestral property that is retained in salamanders but lost in other tetrapods or if it evolved in salamanders. The three-finger protein Prod 1 is implicated in the mechanism of newt limb regeneration, and no orthologs have been found in other vertebrates, thus providing evidence for the second viewpoint. It has also been suggested that this protein could play a role in salamander-specific aspects of limb development. There are ten families of extant salamanders, and Prod 1 has only been identified in two of them to date. It is important to determine if it is present in other families and, particularly, the basal group of two families which diverged approximately 200 MYA.

FINDINGS

We have used polymerase chain reaction (PCR) to identify Prod 1 in a Chinese hynobiid species Batrachuperus longdongensis. We obtained an intestinal transcriptome of the plethodontid Aneides lugubris and, from this, identified a primer which allowed PCR of two Prod 1 genes from this species. All known Prod 1 sequences from nine species in four families have been aligned, and a phylogenetic tree has been derived.

CONCLUSIONS

Prod 1 is found in basal salamanders of the family Hynobiidae, and in at least three other families, so it may be present in all extant salamanders. It remains a plausible candidate to have been involved in the origins of limb regeneration, as well as the apomorphic aspects of limb development.

Variation in salamanders: an essay on genomes, development, and evolution

Regeneration is studied in a few model species of salamanders, but the ten families of salamanders show considerable variation, and this has implications for our understanding of salamander biology. The most recent classification of the families identifies the cryptobranchoidea as the basal group which diverged in the early Jurassic. Variation in the sizes of genomes is particularly obvious, and reflects a major contribution from transposable elements which is already present in the basal group.Limb development has been a focus for evodevo studies, in part because of the variable property of pre-axial dominance which distinguishes salamanders from other tetrapods. This is thought to reflect the selective pressures that operate on a free-living aquatic larva, and might also be relevant for the evolution of limb regeneration. Recent fossil evidence suggests that both pre-axial dominance and limb regeneration were present 300 million years ago in larval temnospondyl amphibians that lived in mountain lakes. A satisfying account of regeneration in salamanders may need to address all these different aspects in the future.

Comparative and phylogenetic perspectives of the cleavage process in tailed amphibians

The order Caudata includes about 660 species and displays a variety of important developmental traits such as cleavage pattern and egg size. However, the cleavage process of tailed amphibians has never been analyzed within a phylogenetic framework. We use published data on the embryos of 36 species concerning the character of the third cleavage furrow (latitudinal, longitudinal or variable) and the magnitude of synchronous cleavage period (up to 3–4 synchronous cell divisions in the animal hemisphere or a considerably longer series of synchronous divisions followed by midblastula transition). Several species from basal caudate families Cryptobranchidae (Andrias davidianus and Cryptobranchus alleganiensis) and Hynobiidae (Onychodactylus japonicus) as well as several representatives from derived families Plethodontidae (Desmognathus fuscus and Ensatina eschscholtzii) and Proteidae (Necturus maculosus) are characterized by longitudinal furrows of the third cleavage and the loss of synchrony as early as the 8-cell stage. By contrast, many representatives of derived families Ambystomatidae and Salamandridae have latitudinal furrows of the third cleavage and extensive period of synchronous divisions. Our analysis of these ontogenetic characters mapped onto a phylogenetic tree shows that the cleavage pattern of large, yolky eggs with short series of synchronous divisions is an ancestral trait for the tailed amphibians, while the data on the orientation of third cleavage furrows seem to be ambiguous with respect to phylogeny. Nevertheless, the midblastula transition, which is characteristic of the model species Ambystoma mexicanum (Caudata) and Xenopus laevis (Anura), might have evolved convergently in these two amphibian orders.