Ylenia Chiari

Comparative performance of the 16S rRNA gene in DNA barcoding of amphibians

BackgroundIdentifying species of organisms by short sequences of DNA has been in the center of ongoing discussions under the terms DNA barcoding or DNA taxonomy. A C-terminal fragment of the mitochondrial gene for cytochrome oxidase subunit I (COI) has been proposed as universal marker for this purpose among animals.ResultsHerein we present experimental evidence that the mitochondrial 16S rRNA gene fulfills the requirements for a universal DNA barcoding marker in amphibians. In terms of universality of priming sites and identification of major vertebrate clades the studied 16S fragment is superior to COI. Amplification success was 100% for 16S in a subset of fresh and well-preserved samples of Madagascan frogs, while various combination of COI primers had lower success rates.COI priming sites showed high variability among amphibians both at the level of groups and closely related species, whereas 16S priming sites were highly conserved among vertebrates. Interspecific pairwise 16S divergences in a test group of Madagascan frogs were at a level suitable for assignment of larval stages to species (1–17%), with low degrees of pairwise haplotype divergence within populations (0–1%).ConclusionWe strongly advocate the use of 16S rRNA as standard DNA barcoding marker for vertebrates to complement COI, especially if samples a priori could belong to various phylogenetically distant taxa and false negatives would constitute a major problem.

Comparative population genomics in animals uncovers the determinants of genetic diversity

Genetic diversity is the amount of variation observed between DNA sequences from distinct individuals of a given species. This pivotal concept of population genetics has implications for species health, domestication, management and conservation. Levels of genetic diversity seem to vary greatly in natural populations and species, but the determinants of this variation, and particularly the relative influences of species biology and ecology versus population history, are still largely mysterious. Here we show that the diversity of a species is predictable, and is determined in the first place by its ecological strategy. We investigated the genome-wide diversity of 76 non-model animal species by sequencing the transcriptome of two to ten individuals in each species. The distribution of genetic diversity between species revealed no detectable influence of geographic range or invasive status but was accurately predicted by key species traits related to parental investment: long-lived or low-fecundity species with brooding ability were genetically less diverse than short-lived or highly fecund ones. Our analysis demonstrates the influence of long-term life-history strategies on species response to short-term environmental perturbations, a result with immediate implications for conservation policies.

Phylogenomic analyses support the position of turtles as the sister group of birds and crocodiles (Archosauria).

BackgroundThe morphological peculiarities of turtles have, for a long time, impeded their accurate placement in the phylogeny of amniotes. Molecular data used to address this major evolutionary question have so far been limited to a handful of markers and/or taxa. These studies have supported conflicting topologies, positioning turtles as either the sister group to all other reptiles, to lepidosaurs (tuatara, lizards and snakes), to archosaurs (birds and crocodiles), or to crocodilians. Genome-scale data have been shown to be useful in resolving other debated phylogenies, but no such adequate dataset is yet available for amniotes.ResultsIn this study, we used next-generation sequencing to obtain seven new transcriptomes from the blood, liver, or jaws of four turtles, a caiman, a lizard, and a lungfish. We used a phylogenomic dataset based on 248 nuclear genes (187,026 nucleotide sites) for 16 vertebrate taxa to resolve the origins of turtles. Maximum likelihood and Bayesian concatenation analyses and species tree approaches performed under the most realistic models of the nucleotide and amino acid substitution processes unambiguously support turtles as a sister group to birds and crocodiles. The use of more simplistic models of nucleotide substitution for both concatenation and species tree reconstruction methods leads to the artefactual grouping of turtles and crocodiles, most likely because of substitution saturation at third codon positions. Relaxed molecular clock methods estimate the divergence between turtles and archosaurs around 255 million years ago. The most recent common ancestor of living turtles, corresponding to the split between Pleurodira and Cryptodira, is estimated to have occurred around 157 million years ago, in the Upper Jurassic period. This is a more recent estimate than previously reported, and questions the interpretation of controversial Lower Jurassic fossils as being part of the extant turtles radiation.ConclusionsThese results provide a phylogenetic framework and timescale with which to interpret the evolution of the peculiar morphological, developmental, and molecular features of turtles within the amniotes.

Evidence for recent gene flow between north-eastern and south-eastern Madagascan poison frogs from a phylogeography of the Mantella cowani group

BackgroundThe genus Mantella, endemic poison frogs of Madagascar with 16 described species, are known in the field of international pet trade and entered under the CITES control for the last four years. The phylogeny and phylogeography of this genus have been recently subject of study for conservation purposes. Here we report on the studies of the phylogeography of the Mantella cowani group using a fragment of 453 bp of the mitochondrial cytochrome b gene from 195 individuals from 21 localities. This group is represented by five forms: M. cowani, a critically endangered species, a vulnerable species, M. haraldmeieri, and the non-threatened M. baroni, M. aff. baroni, and M. nigricans.ResultsThe Bayesian phylogenetic and haplotype network analyses revealed the presence of three separated haplotype clades: (1) M. baroni, M. aff. baroni, M. nigricans, and putative hybrids of M. cowani and M. baroni, (2) M. cowani and putative hybrids of M. cowani and M. baroni, and (3) M. haraldmeieri. The putative hybrids were collected from sites where M. cowani and M. baroni live in sympatry.ConclusionThese results suggest (a) a probable hybridization between M. cowani and M. baroni, (b) a lack of genetic differentiation between M. baroni/M. aff. baroni and M. nigricans, (c) evidence of recent gene-flow between the northern (M. nigricans), eastern (M. baroni), and south-eastern (M. aff. baroni) forms of distinct coloration, and (d) the existence of at least three units for conservation in the Mantella cowani group.

Reference-free transcriptome assembly in non-model animals from next-generation sequencing data.

Next‐generation sequencing (NGS) technologies offer the opportunity for population genomic study of non‐model organisms sampled in the wild. The transcriptome is a convenient and popular target for such purposes. However, designing genetic markers from NGS transcriptome data requires assembling gene‐coding sequences out of short reads. This is a complex task owing to gene duplications, genetic polymorphism, alternative splicing and transcription noise. Typical assembling programmes return thousands of predicted contigs, whose connection to the species true gene content is unclear, and from which SNP definition is uneasy. Here, the transcriptomes of five diverse non‐model animal species (hare, turtle, ant, oyster and tunicate) were assembled from newly generated 454 and Illumina sequence reads. In two species for which a reference genome is available, a new procedure was introduced to annotate each predicted contig as either a full‐length cDNA, fragment, chimera, allele, paralogue, genomic sequence or other, based on the number of, and overlap between, blast hits to the appropriate reference. Analyses showed that (i) the highest quality assemblies are obtained when 454 and Illumina data are combined, (ii) typical de novo assemblies include a majority of irrelevant cDNA predictions and (iii) assemblies can be appropriately cleaned by filtering contigs based on length and coverage. We conclude that robust, reference‐free assembly of thousands of genes from transcriptomic NGS data is possible, opening promising perspectives for transcriptome‐based population genomics in animals. A Galaxy pipeline implementing our best‐performing assembling strategy is provided.

Reference-free transcriptome assembly in non-model animals from next generation sequencing data

Next‐generation sequencing (NGS) technologies offer the opportunity for population genomic study of non‐model organisms sampled in the wild. The transcriptome is a convenient and popular target for such purposes. However, designing genetic markers from NGS transcriptome data requires assembling gene‐coding sequences out of short reads. This is a complex task owing to gene duplications, genetic polymorphism, alternative splicing and transcription noise. Typical assembling programmes return thousands of predicted contigs, whose connection to the species true gene content is unclear, and from which SNP definition is uneasy. Here, the transcriptomes of five diverse non‐model animal species (hare, turtle, ant, oyster and tunicate) were assembled from newly generated 454 and Illumina sequence reads. In two species for which a reference genome is available, a new procedure was introduced to annotate each predicted contig as either a full‐length cDNA, fragment, chimera, allele, paralogue, genomic sequence or other, based on the number of, and overlap between, blast hits to the appropriate reference. Analyses showed that (i) the highest quality assemblies are obtained when 454 and Illumina data are combined, (ii) typical de novo assemblies include a majority of irrelevant cDNA predictions and (iii) assemblies can be appropriately cleaned by filtering contigs based on length and coverage. We conclude that robust, reference‐free assembly of thousands of genes from transcriptomic NGS data is possible, opening promising perspectives for transcriptome‐based population genomics in animals. A Galaxy pipeline implementing our best‐performing assembling strategy is provided.

New evidence for parallel evolution of colour patterns in Malagasy poison frogs ( Mantella )

Malagasy poison frogs of the genus Mantella are diurnal and toxic amphibians of highly variable and largely aposematic coloration. Previous studies provided evidence for several instances of homoplastic colour evolution in this genus but were unable to sufficiently resolve relationships among major species groups or to clarify the phylogenetic position of several crucial taxa. Here, we provide cytochrome b data for 143 individuals of three species in the Mantella madagascariensis group, including four newly discovered populations. Three of these new populations are characterized by highly variable coloration and patterns but showed no conspicuous increase of haplotype diversity which would be expected under a scenario of secondary hybridization or admixture of chromatically uniform populations. Several populations of these variable forms and of M. crocea were geographically interspersed between the distribution areas of Mantella aurantiaca and Mantella milotympanum. This provides further support for the hypothesis that the largely similar uniformly orange colour of the last two species evolved in parallel. Phylogenies based on over 2000 bp of two nuclear genes (Rag‐1 and Rag‐2) identified reliably a clade of the Mantella betsileo and Mantella laevigata groups as sister lineage to the M. madagascariensis group, but did not support species within the latter group as monophyletic. The evolutionary history of these frogs might have been characterized by fast and recurrent evolution of colour patterns, possibly triggered by strong selection pressures and mimicry effects, being too complex to be represented by simple bifurcating models of phylogenetic reconstruction.

Next-generation sequencing of transcriptomes: a guide to RNA isolation in nonmodel animals

Next Generation Sequencing technologies (NGS) are rapidly invading many evolutionary and ecological fields, such as phylogenomics, molecular evolution, population genomics and molecular ecology. Among the potential targets of NGS is transcriptome sequencing, a fast and relatively cheap way to generate massive amounts of coding sequence data, offering promising perspectives for the analysis of molecular diversity in the wild. A number of molecular ecology research groups therefore may switch from DNA‐based to RNA‐based typing in the near future. Sample preparation from natural populations, however, requires specific care and protocols when RNA is the target. Furthermore, NGS sequencing of transcriptome requires high amount of good‐quality RNA. Here we present the results of RNA extraction experiments from various samples of 39 animal species caught in the wild. We compared tissue preparation and storage conditions, evaluated and improved standard RNA extraction protocols, and achieved RNA yield and quality suitable for NGS in all cases. We derive general guidelines for the production of ready‐to‐sequence RNA in nonmodel animals sampled in the field.

Dating cryptodiran nodes: Origin and diversification of the turtle superfamily Testudinoidea

The superfamily Testudinoidea is the most diverse and widely distributed clade of extant turtles. Surprisingly, despite an extensive fossil record, and increasing amount of molecular data available, the temporal origin of this group is still largely unknown. To address this issue, we used a comprehensive molecular dataset to perform phylogenetic and molecular dating analyses, as well as seven fossil constraints to calibrate the ages of the nodes in the phylogeny. The molecular dataset includes the complete mitochondrial genomes of 37 turtle species, including newly sequenced mitochondrial genomes of Phrynops hilarii, Emys orbicularis, Rhinoclemmys punctularia, and Chelonoidis nigra, and four nuclear markers. Our results revealed that the earliest divergences within crown testudinoids occurred around 95.0 Mya, in the early Late Cretaceous, earlier than previously reported, raising new questions about the historical biogeography of this group.

High mitochondrial diversity within and among populations of Malagasy poison frogs.

The diurnal, brightly colored, and toxic frogs of the genus Mantella are among the most prominent representatives of the endemic anuran fauna of Madagascar. Especially three closely related species, M. aurantiaca, M. crocea, and M. milotympanum, are intensively collected for the pet trade, although basic data on their natural history and genetic diversity are still lacking. Our phylogenetic analyses based on 2.8 kbp of partial 16S rRNA, 12S rRNA, cytochrome b, and rhodopsin DNA sequences confirmed that these species belong to one of the five major clades in Mantella, the M. madagascariensis group. A haplotype network constructed using 830 bp of cytochrome b in 49 individuals from seven populations revealed that M. milotympanum and M. crocea have largely similar haplotypes sharing, confirming doubts about the species validity of M. milotympanum and indicating independent evolution of bright orange pattern in M. milotympanum and M. aurantiaca. Further, clustering of four individuals of M. aurantiaca from Andranomena with M. crocea suggests incomplete lineage sorting or introgression resulting from secondary contact of refugial populations. AMOVA confirmed significant intrapopulation nucleotide diversity (>20%). These diversity patterns and our field observations indicate relatively large population sizes. Hence, overcollecting is probably a minor problem and conservation efforts should rather focus on saving some large populations from habitat destruction through logging and forest fires.