Turtle phylogeny: insights from a novel nuclear intron

Morphological Diversity of Turtle Hyoid Apparatus is Linked to Feeding Behavior

The hyoid apparatus of tetrapods is highly diverse in its morphology. It plays an important role in feeding, breathing, sound production, and various other behaviors. Among turtles, the diversity of the hyoid apparatus has been recurrently linked to their habitat. The ossification of the hyoid corpus is often the main trait used in correlations with "niche" occupancy, an ossified corpus being associated with aquatic environments and a cartilaginous corpus with terrestrial life. Most studies conducted so far have focused on species belonging to Testudinoidea, the clade that occupies the biggest diversity of habitats (i.e., terrestrial, semi-terrestrial, and aquatic animals), while other turtle lineages have been largely understudied. We assessed the adult anatomy of the hyoid apparatus of 92 turtle species from all "families", together with ossification sequences from embryological series of 11 species, some described for the first time here. Using nearly 40 different discrete anatomical characters, we discuss the evolutionary patterns and the biological significance of morphological transformations in the turtle hyoid elements. Morphological changes are strongly associated to feeding modes, with several instances of convergent evolution within and outside the Testudines clade, and are not as strongly connected to habitat as previously thought. Some of the hyoid character states we describe are diagnostic of specific turtle clades, thus providing phylogenetically relevant information.

Integrating morphological, molecular and cytogenetic data for F2 sea turtle hybrids diagnosis revealed balanced chromosomal sets

Hybridization could be considered part of the evolutionary history of many species. The hybridization among sea turtle species on the Brazilian coast is atypical and occurs where nesting areas and reproductive seasons overlap. Integrated analysis of morphology and genetics is still scarce, and there is no evidence of the parental chromosome set distribution in sea turtle interspecific hybrids. In this study, chromosome markers previously established for pure sea turtle species were combined with morphological and molecular analyses aiming to recognize genetic composition and chromosome sets in possible interspecific hybrids initially identified by mixed morphology. The data showed that one hybrid could be an F2 individual among Caretta caretta × Eretmochelys imbricata × Chelonia mydas, and another is resulting from backcross between C. caretta × Lepidochelys olivacea. Native alleles of different parental lineages were reported in the hybrids, and, despite this, it was verified that the hybrid chromosome sets were still balanced. Thus, how sea turtle hybridism can affect genetic features in the long term is a concern, as the implications of the crossing-over in hybrid chromosomal sets and the effects on genetic function are still unpredictable.

Phylogenomics reconciles molecular data with the rich fossil record on the origin of living turtles

Although a consensus exists that all living turtles fall within either Pleurodira or Cryptodira clades, estimating when these lineages split is still under debate. Most molecular studies date the split in the Triassic Period, whereas a Jurassic age is unanimous among morphological studies. Each hypothesis implies different paleobiogeographical scenarios to explain early turtle evolution. Here we explored the rich turtle fossil record with the Fossilized Birth-Death (FBD) and the traditional node dating (ND) methods using complete mitochondrial genomes (147 taxa) and a set of nuclear orthologs with over 10 million bp (25 taxa) to date the major splits in Testudines. Our results support an Early Jurassic split (191-182 Ma) for the crown Testudines with great consistency across different dating methods and datasets, with a narrow confidence interval. This result is independently supported by the oldest fossils of Testudines that postdate the Middle Jurassic (174 Ma), which were not used for calibration in this study. This age coincides with the Pangaea fragmentation and the formation of saltwater barriers such as the Atlantic Ocean and the Turgai Strait, supporting that diversification in Testudines was triggered by vicariance. Our ages of the splits in Pleurodira coincide with the geologic events of the Late Jurassic and Early Cretaceous. Conversely, the early Cryptodira radiation remained in Laurasia, and its diversification ensued as all its major lineages expanded their distribution into every continent during the Cenozoic. We provide the first detailed hypothesis of the evolution of Cryptodira in the Southern Hemisphere, in which our time estimates are correlated with each contact between landmasses derived from Gondwana and Laurasia. Although most South American Cryptodira arrived through the Great American Biotic Interchange, our results indicate that the Chelonoidis ancestor probably arrived from Africa through the chain islands of the South Atlantic during the Paleogene. Together, the presence of ancient turtle diversity and the vital role that turtles occupy in marine and terrestrial ecosystems underline South America as a chief area for conservation.

Genomic analyses of the scorpion mud turtle (Kinosternon scorpioides) (Linnaeus, 1766) in insular and continental Colombia: Evidence for multiple conservation and taxonomic units

The turtle genus Kinosternon is widespread with at least 25 species distributed from Mexico to northern Argentina. The taxonomy of this genus is controversial and requires a full revision using both morphological and molecular approaches. In this study, we did a genomic analysis on the species Kinosternon scorpioides distributed in insular and continental Colombia in order to define conservation units. Total DNA was extracted from 24 tissue samples and RADseq genotyping analysis was done. In addition, the intron R35 was amplified and sequenced for a subset of samples. A total of 35,507 SNPs combined with 1,047 bp of the intron were used for spatiotemporal colonization pattern reconstruction and phylogenetic analyses. In addition, SNPs were used for population structure inferences and allele frequency-based analyses. Reciprocal monophyly, significant differences in allele frequencies (Fst = 0.32 - 0.78), and evidence of reproductive isolation (no admixture/geneflow), indicate long-term divergence between groups (2-8 MYA), possibly due to geographical barriers. Four Evolutionarily Significant Units (ESUs) were defined within our samples. One ESU was represented by the insular subspecies K. scorpioides albogulare, found on San Andrés island, and three ESUs were defined for the subspecies K. s. scorpioides in continental Colombia: one trans-Andean, found in northwestern Colombia (Caribbean region) and two cis-Andean, found in eastern and southeastern Colombia in the Orinoco and Amazon regions, respectively. Colonization of this species occurred from an ancestral area on South of Central America region (~ 8.43 MYA), followed by an establishing of current populations on San Andrés Island and then, in the continent. First, in the Colombian Caribbean, next, in the Orinoco, and more recently, in the Amazon. We hypothesize that the emergence of the Panamá Isthmus, as well as the final uplift of the North Eastern Andes and Vaupes Arch, were key event leading to the differentiation of these ESUs. For management and conservation purposes, each of these ESUs should be considered as a separate management unit. A full revision of the taxonomy of the genus Kinosternon is warranted.

Genetic Variation and Population Structure of the Old World Bollworm Helicoverpa armigera (Hübner, 1808) (Lepidoptera: Noctuidae) in Ethiopia

Helicoverpa armigera is one of the most destructive insect pests of economically valuable crops in the world. Despite its economic importance, the population genetic structure of this insect remains unexplored in Ethiopia. To investigate the genetic diversity and population structure of H. armigera, we sampled 170 individuals from 15 populations throughout Ethiopia. We sequenced a fragment of the mitochondrial cytochrome b (cyt b) gene and five exon-primed intron-crossing (EPIC) markers. Twenty cyt b haplotypes with low-to-moderate haplotype diversity (mean Hd = 0.537) and high nucleotide diversity (mean Pi = 0.00339) were identified. The most frequently observed and widely distributed cyt b haplotype was designated as Hap_1 (67.058%), which is identical to sequences found across the globe. Tajima's D and Fu's F for the cyt b data were negative, supporting a model of population expansion. Within populations, a mean of 2.493 alleles/locus was recorded across the five EPIC loci, ranging from 1.200 to 3.600 alleles/locus. The highest mean effective number of alleles/population was 2.369 and the lowest was 1.178. The mean observed heterozygosity (HO) of the five loci (0-0.289; mean 0.104 ± 0.020) was lower than the expected heterozygosity (HE) (0.095-0.523; mean 0.258 ± 0.028). AMOVA detected significant genetic structure with 61% of the total molecular genetic variation of EPIC genotypes occurring between populations, suggesting a considerable degree of differentiation among populations. STRUCTURE analyses clustered the H. armigera populations into three distinct population groups but very low isolation by distance (R2 = 0.0132, P < 0.05).

A global phylogeny of turtles reveals a burst of climate-associated diversification on continental margins

Living turtles are characterized by extraordinarily low species diversity given their age. The clade's extensive fossil record indicates that climate and biogeography may have played important roles in determining their diversity. We investigated this hypothesis by collecting a molecular dataset for 591 individual turtles that, together, represent 80% of all turtle species, including representatives of all families and 98% of genera, and used it to jointly estimate phylogeny and divergence times. We found that the turtle tree is characterized by relatively constant diversification (speciation minus extinction) punctuated by a single threefold increase. We also found that this shift is temporally and geographically associated with newly emerged continental margins that appeared during the Eocene-Oligocene transition about 30 million years before present. In apparent contrast, the fossil record from this time period contains evidence for a major, but regional, extinction event. These seemingly discordant findings appear to be driven by a common global process: global cooling and drying at the time of the Eocene-Oligocene transition. This climatic shift led to aridification that drove extinctions in important fossil-bearing areas, while simultaneously exposing new continental margin habitat that subsequently allowed for a burst of speciation associated with these newly exploitable ecological opportunities.

The complete mitochondrial genome of the endangered Assam Roofed Turtle, Pangshura sylhetensis (Testudines: Geoemydidae): Genomic features and phylogeny

The Assam Roofed Turtle, Pangshura sylhetensis is an endangered and least studied species endemic to India and Bangladesh. The present study decodes the first complete mitochondrial genome of P. sylhetensis (16,568 bp) by using next-generation sequencing. The assembly encodes 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (CR). Most of the genes were encoded on the majority strand, except NADH dehydrogenase subunit 6 (nad6) and eight tRNAs. All PCGs start with an ATG initiation codon, except for Cytochrome oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 5 (nad5), which both start with GTG codon. The study also found the typical cloverleaf secondary structures in most of the predicted tRNA structures, except for serine (trnS1) which lacks of conventional DHU arm and loop. Both Bayesian and maximum-likelihood phylogenetic inference using 13 concatenated PCGs demonstrated strong support for the monophyly of all 52 Testudines species within their respective families and revealed Batagur trivittata as the nearest neighbor of P. sylhetensis. The mitogenomic phylogeny with other amniotes is congruent with previous research, supporting the sister relationship of Testudines and Archosaurians (birds and crocodilians). Additionally, the mitochondrial Gene Order (GO) analysis indicated plesiomorphy with the typical vertebrate GO in most of the Testudines species.

Genomic analyses reveal two species of the matamata (Testudines: Chelidae: Chelus spp.) and clarify their phylogeography

The matamata is one of the most charismatic turtles on earth, widely distributed in northern South America. Debates have occurred over whether or not there should be two subspecies or species recognized due to its geographic variation in morphology. Even though the matamata is universally known, its natural history, conservation status and biogeography are largely unexplored. In this study we examined the phylogeographic differentiation of the matamata based on three mitochondrial DNA fragments (2168 bp of the control region, cytochrome oxidase subunit I, and the cytochrome b gene), one nuclear genomic DNA fragment (1068 bp of the R35 intron) and 1661 Single Nucleotide Polymorphisms (SNPs). Our molecular and morphological analyses revealed the existence of two distinct, genetically deeply divergent evolutionary lineages of matamatas that separated in the late Miocene (approximately 12.7 million years ago), corresponding well to the time when the Orinoco Basin was established. As a result of our analyses, we describe the genetically and morphologically highly distinct matamata from the Orinoco and Río Negro Basins and the Essequibo drainage as a species new to science (Chelus orinocensis sp. nov.). Chelus fimbriata sensu stricto is distributed in the Amazon Basin and the Mahury drainage. Additionally, the analyses revealed that each species displays phylogeographic differentiation. For C. orinocensis, there is moderate mitochondrial differentiation between the Orinoco and the Río Negro. For C. fimbriata, there is more pronounced differentiation matching different river systems. One mitochondrial clade was identified from the Amazon, Ucayali, and Mahury Rivers, and another one from the Madeira and Jaci Paraná Rivers. The C. orinocensis in the Essequibo and Branco Rivers have haplotypes that constitute a third clade clustering with C. fimbriata. Phylogenetic analyses of the R35 intron and SNP data link the matamatas from the Essequibo and Branco with the new species, suggesting past gene flow and old mitochondrial introgression. Chelus orinocensis is collected for the pet trade in Colombia and Venezuela. However, neither the extent of the harvest nor its impact are known. Hence, it is crucial to gather more information and to assess its exploitation throughout its distribution range to obtain a better understanding of its conservation status and to design appropriate conservation and management procedures. RESUMEN: La matamata es una de las tortugas más carismáticas del mundo, ampliamente distribuida en el norte de Sudamérica. Debido a su variación morfológica geográfica, se debate sobre el reconocimiento de dos subespecies o especies. A pesar de que la matamata es universalmente conocida, su historia natural, estado de conservación y biogeografía han sido muy poco estudiados. En este estudio examinamos la diferenciación filogeográfica de las matamatas en base ​​a tres fragmentos de ADN mitocondrial (2168 pb de la región de control, la subunidad I del citocromo oxidasa y el gen del citocromo b), un fragmento de ADN genómico nuclear (1068 pb del intrón R35) y 1661 polimorfismos de nucleótido único (SNPs). Nuestros análisis moleculares y morfológicos revelaron la existencia de dos linajes evolutivos distintos de matamatas, genéticamente divergentes que se separaron en el Mioceno tardio (hace aproximadamente 12.7 millones de años), correspondiendo al tiempo en que se estableció la cuenca del Orinoco. Como resultado de nuestros análisis, describimos las genéticamente y morfológicamente distintas matamatas de las cuencas del Orinoco, Río Negro y Essequibo como una especie nueva para la ciencia (Chelus orinocensis sp. nov.). Chelus fimbriata sensu stricto se distribuye en la cuenca del Amazonas y en el drenaje del Mahury. Adicionalmente, los análisis revelaron que cada especie muestra diferenciación filogeográfica. Para C. orinocensis, hay una moderada diferenciación mitocondrial entre el Orinoco y el Río Negro. Para C. fimbriata, hay una diferenciación más pronunciada, concordando con los diferentes sistemas fluviales. Se identificó un clado de los ríos Amazonas, Ucayali y Mahury y otro de los ríos Madeira y Jaci Paraná. Las C. orinocensis de los ríos Essequibo y Branco tienen haplotipos que constituyen un tercer clado que se agrupa con C. fimbriata. Los análisis filogenéticos del intrón R35 y los datos de SNP asocian las matamatas de Essequibo y Branco con la nueva especie, sugiriendo flujo de genes pasado ​​e introgresión mitocondrial antigua. Chelus orinocensis se colecta para el comercio de mascotas en Colombia y Venezuela. Sin embargo, ni se conoce el alcance de las colectas ni su impacto. Por lo tanto, es crucial recopilar más información y evaluar su explotación en todo su rango de distribución, comprender mejor su estado de conservación y para diseñar acciones apropiadas de conservación y manejo.

Sex is determined by XX/XY sex chromosomes in Australasian side-necked turtles (Testudines: Chelidae)

Turtles demonstrate variability in sex determination and, hence, constitute an excellent model for the evolution of sex chromosomes. Notably, the sex determination of the freshwater turtles from the family Chelidae, a species-rich group with wide geographical distribution in the southern hemisphere, is still poorly explored. Here we documented the presence of an XX/XY sex determination system in seven species of the Australasian chelid genera Chelodina, Emydura, and Elseya by conventional (karyogram reconstruction, C-banding) and molecular cytogenetic methods (comparative genome hybridization, in situ hybridization with probes specific for GATA microsatellite motif, the rDNA loci, and the telomeric repeats). The sex chromosomes are microchromosomes in all examined species of the genus Chelodina. In contrast, the sex chromosomes are the 4^th largest pair of macrochromosomes in the genera Emydura and Elseya. Their X chromosomes are submetacentric, while their Y chromosomes are metacentric. The chelid Y chromosomes contain a substantial male-specific genomic region with an accumulation of the GATA microsatellite motif, and occasionally, of the rDNA loci and telomeric repeats. Despite morphological differences between sex chromosomes, we conclude that male heterogamety was likely already present in the common ancestor of Chelodina, Emydura and Elseya in the Mesozoic period.

Phylogeography of Ticks (Acari: Ixodida)

Improved understanding of tick phylogeny has allowed testing of some biogeographical patterns. On the basis of both literature data and a meta-analysis of available sequence data, there is strong support for a Gondwanan origin of Ixodidae, and probably Ixodida. A particularly strong pattern is observed for the genus Amblyomma, which appears to have originated in Antarctica/southern South America, with subsequent dispersal to Australia. The endemic Australian lineages of Ixodidae (no other continent has such a pattern) appear to result from separate dispersal events, probably from Antarctica. Minimum ages for a number of divergences are determined as part of an updated temporal framework for tick evolution. Alternative hypotheses for tick evolution, such as a very old Pangean group, a Northern hemisphere origin, or an Australian origin, fit less well with observed phylogeographic patterns.

Molecular detection of hemogregarines and haemosporidians in Brazilian free-living testudines

Morphological and molecular techniques were used to investigate the presence of hemogregarines and haemosporidians in biological samples of free-living Geoffroy's side-necked turtles (Phrynops geoffroanus) and Giant Amazon turtles (Podocnemis expansa) from Brazil. No evolutionary form of haemosporidians or hemogregarines were observed in the blood smears of 83 P. geoffroanus samples, and there were no meronts in the histological sections of 31 necropsied P. geoffroanus samples. All DNA samples extracted from P. geoffroanus tissues and blood aliquots were negative in haemosporidian PCR assays (based on the mitochondrial cytochrome b gene) and hemogregarine PCR assays (based on the 18S rRNA gene). In the analysis of blood smears of all seven Podocnemis expansa evaluated, gametocytes of hemogregarines were observed. The seven P. expansa were negative in the haemosporidian PCR assays. Moreover, hemogregarine DNA was detected in blood samples from all of the sampled P. expansa. The phylogenetic maximum likelihood inference and probabilistic Bayesian inference revealed five closely related genotypes that formed a monophyletic group. There was also a sister group to the lineage that consisted of Haemogregarina spp. of freshwater turtles from Canada, Italy, Mozambique, Kenya, Gabon, Vietnam, and China. The findings suggest that free-living P. expansa were parasitized by a new genotype or even a possible new species of the genus Haemogregarina. Haemosporidians and hemogregarines are not frequently found in P. geoffroanus in the studied region under the local conditions of that period.

Molecular detection of hemogregarines and haemosporidians in Brazilian free-living testudines

Morphological and molecular techniques were used to investigate the presence of hemogregarines and haemosporidians in biological samples of free-living Geoffroy's side-necked turtles (Phrynops geoffroanus) and Giant Amazon turtles (Podocnemis expansa) from Brazil. No evolutionary form of haemosporidians or hemogregarines were observed in the blood smears of 83 P. geoffroanus samples, and there were no meronts in the histological sections of 31 necropsied P. geoffroanus samples. All DNA samples extracted from P. geoffroanus tissues and blood aliquots were negative in haemosporidian PCR assays (based on the mitochondrial cytochrome b gene) and hemogregarine PCR assays (based on the 18S rRNA gene). In the analysis of blood smears of all seven Podocnemis expansa evaluated, gametocytes of hemogregarines were observed. The seven P. expansa were negative in the haemosporidian PCR assays. Moreover, hemogregarine DNA was detected in blood samples from all of the sampled P. expansa. The phylogenetic maximum likelihood inference and probabilistic Bayesian inference revealed five closely related genotypes that formed a monophyletic group. There was also a sister group to the lineage that consisted of Haemogregarina spp. of freshwater turtles from Canada, Italy, Mozambique, Kenya, Gabon, Vietnam, and China. The findings suggest that free-living P. expansa were parasitized by a new genotype or even a possible new species of the genus Haemogregarina. Haemosporidians and hemogregarines are not frequently found in P. geoffroanus in the studied region under the local conditions of that period.

•

Hemogregarines was detected in blood samples of free-living Brazilian testudines.

•

Gametocytes of hemogregarines were observed in Podocnemis expansa blood smears.

•

Hemogregarines DNA fragments based on the 18S rRNA gene were detected in P. expansa.

•

We propose that P. expansa were parasitized by a new genotype of Haemogregarina.

•

Haemosporidians was not observed in either P. expansa or Phrynops geoffroanus samples.

Phylogenetic relationships and divergence dates of softshell turtles (Testudines: Trionychidae) inferred from complete mitochondrial genomes

The softshell turtles (Trionychidae) are one of the most widely distributed reptile groups in the world, and fossils have been found on all continents except Antarctica. The phylogenetic relationships among members of this group have been previously studied; however, disagreements regarding its taxonomy, its phylogeography and divergence times are still poorly understood as well. Here, we present a comprehensive mitogenomic study of softshell turtles. We sequenced the complete mitochondrial genomes of 10 softshell turtles, in addition to the GenBank sequence of Dogania subplana, Lissemys punctata, Trionyx triunguis, which cover all extant genera within Trionychidae except for Cyclanorbis and Cycloderma. These data were combined with other mitogenomes of turtles for phylogenetic analyses. Divergence time calibration and ancestral reconstruction were calculated using BEAST and RASP software, respectively. Our phylogenetic analyses indicate that Trionychidae is the sister taxon of Carettochelyidae, and support the monophyly of Trionychinae and Cyclanorbinae, which is consistent with morphological data and molecular analysis. Our phylogenetic analyses have established a sister taxon relationship between the Asian Rafetus and the Asian Palea + Pelodiscus + Dogania + Nilssonia + Amyda, whereas a previous study grouped the Asian Rafetus with the American Apalone. The results of divergence time estimates and area ancestral reconstruction show that extant Trionychidae originated in Asia at around 108 million years ago (MA), and radiations mainly occurred during two warm periods, namely Late Cretaceous-Early Eocene and Oligocene. By combining the estimated divergence time and the reconstructed ancestral area of softshell turtles, we determined that the dispersal of softshell turtles out of Asia may have taken three routes. Furthermore, the times of dispersal seem to be in agreement with the time of the India-Asia collision and opening of the Bering Strait, which provide evidence for the accuracy of our estimation of divergence time. Overall, the mitogenomes of this group were used to explore the origin and dispersal route of Trionychidae and have provided new insights on the evolution of this group.

Assessing models of speciation under different biogeographic scenarios; an empirical study using multi-locus and RNA-seq analyses

Evolutionary biology often seeks to decipher the drivers of speciation, and much debate persists over the relative importance of isolation and gene flow in the formation of new species. Genetic studies of closely related species can assess if gene flow was present during speciation, because signatures of past introgression often persist in the genome. We test hypotheses on which mechanisms of speciation drove diversity among three distinct lineages of desert tortoise in the genus Gopherus. These lineages offer a powerful system to study speciation, because different biogeographic patterns (physical vs. ecological segregation) are observed at opposing ends of their distributions. We use 82 samples collected from 38 sites, representing the entire species' distribution and generate sequence data for mtDNA and four nuclear loci. A multilocus phylogenetic analysis in *BEAST estimates the species tree. RNA‐seq data yield 20,126 synonymous variants from 7665 contigs from two individuals of each of the three lineages. Analyses of these data using the demographic inference package ∂a∂i serve to test the null hypothesis of no gene flow during divergence. The best‐fit demographic model for the three taxa is concordant with the *BEAST species tree, and the ∂a∂i analysis does not indicate gene flow among any of the three lineages during their divergence. These analyses suggest that divergence among the lineages occurred in the absence of gene flow and in this scenario the genetic signature of ecological isolation (parapatric model) cannot be differentiated from geographic isolation (allopatric model).

The differentiation of tuna (family: Scombridae) products through the PCR-based analysis of the cytochrome b gene and parvalbumin introns

BACKGROUND

In spite of the many studies performed over the years, there are still problems in the authentication of closely related tuna species, not only for canned fish but also for raw products. With the aim of providing screening methods to identify different tuna species and related scombrids, segments of mitochondrial cytochrome b (cyt b) and nuclear parvalbumin genes were amplified and sequenced or subjected to single-strand conformation polymorphism (SSCP) and restriction fragment length polymorphism (RFLP) analyses.

RESULTS

The nucleotide diagnostic sites in the cyt b gene of five tuna species from Indonesia were determined in this study and used to construct a phylogenetic tree. In addition, the suitability of the nuclear gene that encodes parvalbumin for the differentiation of tuna species was determined by SSCP and RFLP analyses of an intron segment. RFLP differentiated Thunnus albacares and from T. obesus, and fish species in the Thunnus genus could be distinguished from bullet tuna (Auxis rochei) by SSCP.

CONCLUSIONS

Parvalbumin-based polymerase chain reaction systems could serve as an additional tool in the detection and identification of tuna and other Scombridae fish species for routine seafood control. This reaction can be performed in addition to the cyt b analysis as previously described.

Complete mitochondrial genome of Camponotus atrox (Hymenoptera: Formicidae): a new tRNA arrangement in Hymenoptera

We sequenced the complete mitochondrial (mt) genome of Camponotus atrox (Hymenoptera: Formicidae), which is only distributed in Korea. The genome was 16 540 bp in size and contained typical sets of genes (13 protein-coding genes, 22 tRNAs, and 2 rRNAs). The C. atrox A+T-rich region, at 1402 bp, was the longest of all sequenced ant genomes and was composed of an identical tandem repeat consisting of six 100-bp copies and one 96-bp copy. A total of 315 bp of intergenic spacer sequence was spread over 23 regions. An alignment of the spacer sequences in ants was largely feasible among congeneric species, and there was substantial sequence divergence, indicating their potential use as molecular markers for congeneric species. The A/T contents at the first and second codon positions of protein-coding genes (PCGs) were similar for ant species, including C. atrox (73.9% vs. 72.3%, on average). With increased taxon sampling among hymenopteran superfamilies, differences in the divergence rates (i.e., the non-synonymous substitution rates) between the suborders Symphyta and Apocrita were detected, consistent with previous results. The C. atrox mt genome had a unique gene arrangement, trnI-trnM-trnQ, at the A+T-rich region and ND2 junction (underline indicates inverted gene). This may have originated from a tandem duplication of trnM-trnI, resulting in trnM-trnI-trnM-trnI-trnQ, and the subsequent loss of the first trnM and second trnI, resulting in trnI-trnM-trnQ.

Portlandemys gracilis n. sp., a New Coastal Marine Turtle from the Late Jurassic of Porrentruy (Switzerland) and a Reconsideration of Plesiochelyid Cranial Anatomy

BACKGROUND

Several groups of stem cryptodires became adapted to coastal marine environments as early as the Late Jurassic, 40 million years before the Pan-Chelonioidea. The Plesiochelyidae are a major component of this first radiation of crown-group turtles into marine habitats. They are abundant in many European localities, but their systematics is still greatly confused. Only three species are represented by cranial material: Plesiochelys etalloni, Plesiochelys planiceps, and Portlandemys mcdowelli.

METHODOLOGY/PRINCIPAL FINDINGS

In the present study, we describe a cranium and a mandible from the Kimmeridgian of Porrentruy (Switzerland), which we refer to a new species, Portlandemys gracilis n. sp. This new taxon differs from Portlandemys mcdowelli in several aspects of the cranium and mandible, notably in being generally more gracile, but the two species share a narrow skull, a more acute angle between the labial ridges on the mandible, and a unique configuration of the anterodorsal part of the basicranium. The cranial anatomy of plesiochelyid turtles is discussed in details based primarily on these new specimens and new cranial material of Plesiochelys etalloni from Solothurn, Switzerland.

CONCLUSIONS/SIGNIFICANCE

Several characters (e.g., the contribution of the parietal to the foramen nervi trigemini, the configuration of the dorsum sellae and sella turcica, the presence of an infolding ridge on the posterior surface of the quadrate) appear as potential candidates to help elucidate plesiochelyid relationships. Some of these characters are included in a previously published phylogenetic dataset and help to stabilize the relationships of plesiochelyid turtles and closely related taxa. For the first time, our results suggest that plesiochelyids, 'Thalassemys' moseri, and Solnhofia parsonsi (representing the Eurysternidae) form a clade at the base of Eucryptodira.

A phylogenetic analysis of the genus Fragaria (strawberry) using intron-containing sequence from the ADH-1 gene

The genus Fragaria encompasses species at ploidy levels ranging from diploid to decaploid. The cultivated strawberry, Fragaria×ananassa, and its two immediate progenitors, F. chiloensis and F. virginiana, are octoploids. To elucidate the ancestries of these octoploid species, we performed a phylogenetic analysis using intron-containing sequences of the nuclear ADH-1 gene from 39 germplasm accessions representing nineteen Fragaria species and one outgroup species, Dasiphora fruticosa. All trees from Maximum Parsimony and Maximum Likelihood analyses showed two major clades, Clade A and Clade B. Each of the sampled octoploids contributed alleles to both major clades. All octoploid-derived alleles in Clade A clustered with alleles of diploid F. vesca, with the exception of one octoploid allele that clustered with the alleles of diploid F. mandshurica. All octoploid-derived alleles in clade B clustered with the alleles of only one diploid species, F. iinumae. When gaps encoded as binary characters were included in the Maximum Parsimony analysis, tree resolution was improved with the addition of six nodes, and the bootstrap support was generally higher, rising above the 50% threshold for an additional nine branches. These results, coupled with the congruence of the sequence data and the coded gap data, validate and encourage the employment of sequence sets containing gaps for phylogenetic analysis. Our phylogenetic conclusions, based upon sequence data from the ADH-1 gene located on F. vesca linkage group II, complement and generally agree with those obtained from analyses of protein-encoding genes GBSSI-2 and DHAR located on F. vesca linkage groups V and VII, respectively, but differ from a previous study that utilized rDNA sequences and did not detect the ancestral role of F. iinumae.

The advantages of going large: genome-wide SNPs clarify the complex population history and systematics of the threatened western pond turtle

As the field of phylogeography has matured, it has become clear that analyses of one or a few genes may reveal more about the history of those genes than the populations and species that are the targets of study. To alleviate these concerns, the discipline has moved towards larger analyses of more individuals and more genes, although little attention has been paid to the qualitative or quantitative gains that such increases in scale and scope may yield. Here, we increase the number of individuals and markers by an order of magnitude over previously published work to comprehensively assess the phylogeographical history of a well-studied declining species, the western pond turtle (Emys marmorata). We present a new analysis of 89 independent nuclear SNP markers and one mitochondrial gene sequence scored for rangewide sampling of >900 individuals, and compare these to smaller-scale, rangewide genetic and morphological analyses. Our enlarged SNP data fundamentally revise our understanding of evolutionary history for this lineage. Our results indicate that the gains from greatly increasing both the number of markers and individuals are substantial and worth the effort, particularly for species of high conservation concern such as the pond turtle, where accurate assessments of population history are a prerequisite for effective management.

Cryptic population structuring and the role of the Isthmus of Tehuantepec as a gene flow barrier in the critically endangered Central American River Turtle

The critically endangered Central American River Turtle (Dermatemys mawii) is the only remaining member of the Dermatemydidae family, yet little is known about its population structuring. In a previous study of mitochondrial (mt) DNA in the species, three main lineages were described. One lineage (Central) was dominant across most of the range, while two other lineages were restricted to Papaloapan (PAP; isolated by the Isthmus of Tehuantepec and the Sierra de Santa Marta) or the south-eastern part of the range (1D). Here we provide data from seven polymorphic microsatellite loci and the R35 intron to re-evaluate these findings using DNA from the nuclear genome. Based on a slightly expanded data set of a total of 253 samples from the same localities, we find that mtDNA and nuclear DNA markers yield a highly congruent picture of the evolutionary history and population structuring of D. mawii. While resolution provided by the R35 intron (sequenced for a subset of the samples) was very limited, the microsatellite data revealed pronounced population structuring. Within the Grijalva-Usumacinta drainage basin, however, many populations separated by more than 300 kilometers showed signals of high gene flow. Across the entire range, neither mitochondrial nor nuclear DNA show a significant isolation-by-distance pattern, but both genomes highlight that the D. mawii population in the Papaloapan basin is genetically distinctive. Further, both marker systems detect unique genomic signals in four individuals with mtDNA clade 1D sampled on the southeast edge of the Grijalva-Usumacinta basin. These individuals may represent a separate cryptic taxon that is likely impacted by recent admixture.

Resolving the phylogenetic history of the short-necked turtles, genera Elseya and Myuchelys (Testudines: Chelidae) from Australia and New Guinea

Phylogenetic relationships and taxonomy of the short-necked turtles of the genera Elseya, Myuchelys, and Emydura in Australia and New Guinea have long been debated as a result of conflicting hypotheses supported by different data sets and phylogenetic analyses. To resolve this contentious issue, we analyzed sequences from two mitochondrial genes (cytochrome b and ND4) and one nuclear intron gene (R35) from all species of the genera Elseya, Myuchelys, Emydura, and their relatives. Phylogenetic analyses using three methods (maximum parsimony, maximum likelihood, and Bayesian inference) produce a single, well resolved, and strongly corroborated hypothesis, which provides support for the three genera, with the exception that the genus Myuchelys is paraphyletic - Myuchelys purvisi is the sister taxon to the remaining Elseya, Myuchelys and Emydura. A new genus is proposed for the species Myuchelys purvisi to address this paraphyletic relationship. Time-calibration analysis suggests that diversification of the group in Australia coincides with periods of aridification in the late Eocene and between the mid-Miocene and early Pliocene. Other speciation events occurred during the faunal exchange between Australia and the island of New Guinea during the late Miocene and early Pliocene. Lineages distributed in New Guinea are likely influenced by the complex geologic history of the island, and include cryptic species diversity.

Assessing DNA barcoding as a tool for species identification and data quality control

In recent years, the number of sequences of diverse species submitted to GenBank has grown explosively and not infrequently the data contain errors. This problem is extensively recognized but not for invalid or incorrectly identified species, sample mixed-up, and contamination. DNA barcoding is a powerful tool for identifying and confirming species and one very important application involves forensics. In this study, we use DNA barcoding to detect erroneous sequences in GenBank by evaluating deep intraspecific and shallow interspecific divergences to discover possible taxonomic problems and other sources of error. We use the mitochondrial DNA gene encoding cytochrome b (Cytb) from turtles to test the utility of barcoding for pinpointing potential errors. This gene is widely used in phylogenetic studies of the speciose group. Intraspecific variation is usually less than 2.0% and in most cases it is less than 1.0%. In comparison, most species differ by more than 10.0% in our dataset. Overlapping intra- and interspecific percentages of variation mainly involve problematic identifications of species and outdated taxonomies. Further, we detect identical problems in Cytb from Insectivora and Chiroptera. Upon applying this strategy to 47,524 mammalian CoxI sequences, we resolve a suite of potentially problematic sequences. Our study reveals that erroneous sequences are not rare in GenBank and that the DNA barcoding can serve to confirm sequencing accuracy and discover problems such as misidentified species, inaccurate taxonomies, contamination, and potential errors in sequencing.

Genetic introgression and hybridization in Antillean freshwater turtles (Trachemys) revealed by coalescent analyses of mitochondrial and cloned nuclear markers

Determining whether a conflict between gene trees and species trees represents incomplete lineage sorting (ILS) or hybridization involving native and/or invasive species has implications for reconstructing evolutionary relationships and guiding conservation decisions. Among vertebrates, turtles represent an exceptional case for exploring these issues because of the propensity for even distantly related lineages to hybridize. In this study we investigate a group of freshwater turtles (Trachemys) from a part of its range (the Greater Antilles) where it is purported to have undergone reticulation events from both natural and anthropogenic processes. We sequenced mtDNA for 83 samples, sequenced three nuDNA markers for 45 samples, and cloned 29 polymorphic sequences, to identify species boundaries, hybridization, and intergrade zones for Antillean Trachemys and nearby mainland populations. Initial coalescent analyses of phased nuclear alleles (using (*)BEAST) recovered a Bayesian species tree that strongly conflicted with the mtDNA phylogeny and traditional taxonomy, and appeared to be confounded by hybridization. Therefore, we undertook exploratory phylogenetic analyses of mismatched alleles from the "coestimated" gene trees (Heled and Drummond, 2010) in order to identify potential hybrid origins. The geography, morphology, and sampling context of most samples with potential introgressed alleles suggest hybridization over ILS. We identify contact zones between different species on Jamaica (T. decussata × T. terrapen), on Hispaniola (T. decorata × T. stejnegeri), and in Central America (T. emolli × T. venusta). We are unable to determine whether the distribution of T. decussata on Jamaica is natural or the result of prehistoric introduction by Native Americans. This uncertainty means that the conservation status of the Jamaican T. decussata populations and contact zone with T. terrapen are unresolved. Human-mediated dispersal events were more conclusively implicated for the prehistoric translocation of T. stejnegeri between Puerto Rico and Hispaniola, as well as the more recent genetic pollution of native species by an invasive pet turtle native to the USA (T. scripta elegans). Finally, we test the impact of introgressed alleles using the multispecies coalescent in a Bayesian framework and show that studies that do not phase heterozygote sequences of hybrid individuals may recover the correct species tree, but overall support for clades that include hybrid individuals may be reduced.

Forelimb kinematics during swimming in the pig-nosed turtle, Carettochelys insculpta, compared with other turtle taxa: rowing versus flapping, convergence versus intermediacy

Animals that swim using appendages do so by way of rowing and/or flapping motions. Often considered discrete categories, rowing and flapping are more appropriately regarded as points along a continuum. The pig-nosed turtle, Carettochelys insculpta, is unusual in that it is the only freshwater turtle to have limbs modified into flippers and swim via synchronous forelimb motions that resemble dorsoventral flapping, traits that evolved independently from their presence in sea turtles. We used high-speed videography to quantify forelimb kinematics in C. insculpta and a closely related, highly aquatic rower (Apalone ferox). Comparisons of our new data with those previously collected for a generalized freshwater rower (Trachemys scripta) and a flapping sea turtle (Caretta caretta) allow us to: (1) more precisely quantify and characterize the range of limb motions used by flappers versus rowers, and (2) assess whether the synchronous forelimb motions of C. insculpta can be classified as flapping (i.e. whether they exhibit forelimb kinematics and angles of attack more similar to closely related rowing species or more distantly related flapping sea turtles). We found that the forelimb kinematics of previously recognized rowers (T. scripta and A. ferox) were most similar to each other, but that those of C. insculpta were more similar to rowers than to flapping C. caretta. Nevertheless, of the three freshwater species, C. insculpta was most similar to flapping C. caretta. 'Flapping' in C. insculpta is achieved through humeral kinematics very different from those in C. caretta, with C. insculpta exhibiting significantly more anteroposterior humeral motion and protraction, and significantly less dorsoventral humeral motion and depression. Based on several intermediate kinematic parameters and angle of attack data, C. insculpta may in fact represent a synchronous rower or hybrid rower-flapper, suggesting that traditional views of C. insculpta as a flapper should be revised.

Species boundaries and phylogenetic relationships in the critically endangered Asian box turtle genus Cuora

Turtles are currently the most endangered major clade of vertebrates on earth, and Asian box turtles (Cuora) are in catastrophic decline. Effective management of this diverse turtle clade has been hampered by human-mediated, and perhaps natural hybridization, resulting in discordance between mitochondrial and nuclear markers and confusion regarding species boundaries and phylogenetic relationships among hypothesized species of Cuora. Here, we present analyses of mitochondrial and nuclear DNA data for all 12 currently hypothesized species to resolve both species boundaries and phylogenetic relationships. Our 15-gene, 40-individual nuclear data set was frequently in conflict with our mitochondrial data set; based on its general concordance with published morphological analyses and the strength of 15 independent estimates of evolutionary history, we interpret the nuclear data as representing the most reliable estimate of species boundaries and phylogeny of Cuora. Our results strongly reiterate the necessity of using multiple nuclear markers for phylogeny and species delimitation in these animals, including any form of DNA "barcoding", and point to Cuora as an important case study where reliance on mitochondrial DNA can lead to incorrect species identification.

A novel primer set for multilocus phylogenetic inference in East African cichlid fishes

The cichlid fishes in the East African Great Lakes are a prime model system for the study of adaptive radiation. Therefore, the availability of an elaborate phylogenetic framework is an important prerequisite. Previous phylogenetic hypotheses on East African cichlids are mainly based on mitochondrial and/or fragment-based markers, and, to date, no taxon-rich phylogeny exists that is based on multilocus DNA sequence data. Here, we present the design of an extensive new primer set (24 nuclear makers) for East African cichlids that will be used for multilocus phylogenetic analyses in the future. The primers are designed to work for both Sanger sequencing and next-generation sequencing with the 454 technology. As a proof of principle, we validate these primers in a phylogenetically representative set of 16 cichlid species from Lake Tanganyika and main river systems in the area and provide a basic evaluation of the markers with respect to marker length and diversity indices.

Weak genetic divergence between the two South American toad-headed turtles Mesoclemmys dahli and M. zuliae (Testudines: Pleurodira: Chelidae)

Mesoclemmys dahli and M. zuliae are two endangered, little-known toad-headed turtles with small distribution ranges in Colombia and Venezuela, respectively. Using the mitochondrial cytochrome b gene as a marker, we investigate their phylogeographic differentiation. Furthermore, based on 2341 bp of mtDNA and 2109 bp of nDNA of M. dahli, M. zuliae and allied chelid turtles, we infer their divergence time using a fossil-calibrated relaxed molecular clock approach. Mesoclemmys dahli and M. zuliae are closely related species, with an estimated mean divergence time of 10.6 million years. This estimate correlates with the uplift of the Serranía de Perijá, an Andean mountain chain separating their distribution ranges, suggesting that this event could have caused the evolution of the two species. Haplotype and nucleotide diversities of M. dahli are markedly higher than in Podocnemis lewyana, another endemic turtle species of Colombia. This pronounced dissimilarity may reflect differences in the phylogeographies and demographic histories of the two species, but also different habitat preferences.

Comparing and combining distance-based and character-based approaches for barcoding turtles

Molecular barcoding can serve as a powerful tool in wildlife forensics and may prove to be a vital aid in conserving organisms that are threatened by illegal wildlife trade, such as turtles (Order Testudines). We produced cytochrome oxidase subunit one (COI) sequences (650 bp) for 174 turtle species and combined these with publicly available sequences for 50 species to produce a data set representative of the breadth of the order. Variability within the barcode region was assessed, and the utility of both distance-based and character-based methods for species identification was evaluated. For species in which genetic material from more than one individual was available (n = 69), intraspecific divergences were 1.3% on average, although divergences greater than the customary 2% barcode threshold occurred within 15 species. High intraspecific divergences could indicate species with a high degree of internal genetic structure or possibly even cryptic species, although introgression is also probable in some of these taxa. Divergences between species of the same genus were 6.4% on average; however, 49 species were <2% divergent from congeners. Low levels of interspecific divergence could be caused by recent evolutionary radiations coupled with the low rates of mtDNA evolution previously observed in turtles. Complementing distance-based barcoding with character-based methods for identifying diagnostic sets of nucleotides provided better resolution in several cases where distance-based methods failed to distinguish species. An online identification engine was created to provide character-based identifications. This study constitutes the first comprehensive barcoding effort for this seriously threatened order.

Paleogenomic data suggest mammal-like genome size in the ancestral amniote and derived large genome size in amphibians

An unsolved question in evolutionary genomics is whether amniote genomes have been expanding or contracting since the common ancestor of this diverse group. Here, we report on the polarity of amniote genome size evolution using genome size estimates for 14 extinct tetrapod genera from the Paleozoic and early Mesozoic Eras using osteocyte lacunae size as a correlate. We find substantial support for a phylogenetically controlled regression model relating genome size to osteocyte lacunae size (P of slopes <0.01, r²=0.65, phylogenetic signal λ=0.83). Genome size appears to have been homogeneous across Paleozoic crown-tetrapod lineages (average haploid genome size 2.9-3.7 pg) with values similar to those of extant mammals. The differentiation in genome size and underlying architecture among extant tetrapod lineages likely evolved in the Mesozoic and Cenozoic Eras, with expansion in amphibians, contractions along the diapsid lineage, and no directional change within the synapsid lineage leading to mammals.