Extensive Morphological Convergence and Rapid Radiation in the Evolutionary History of the Family Geoemydidae (Old World Pond Turtles) Revealed by SINE Insertion Analysis

Matrilineal phylogeny and habitat suitability of the endangered spotted pond turtle (Geoclemys hamiltonii; Testudines: Geoemydidae): a two-dimensional approach to forecasting future conservation consequences

The spotted pond turtle (Geoclemys hamiltonii) is a threatened and less explored species endemic to Bangladesh, India, Nepal, and Pakistan. To infer structural variation and matrilineal phylogenetic interpretation, the present research decoded the mitogenome of G. hamiltonii (16,509 bp) using next-generation sequencing technology. The mitogenome comprises 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one AT-rich control region (CR) with similar strand symmetry in vertebrates. The ATG was identified as a start codon in most of the PCGs except Cytochrome oxidase subunit 1 (cox1), which started with the GTG codon. The non-coding CR of G. hamiltonii was determined to have a unique structure and variation in different domains and stem-loop secondary structure as compared with other Batagurinae species. The PCGs-based Bayesian phylogeny inferred strong monophyletic support for all Batagurinae species and confirmed the sister relationship of G. hamiltonii with Pangshura and Batagur taxa. We recommend generating more mitogenomic data for other Batagurinae species to confirm their population structure and evolutionary relationships. In addition, the present study aims to infer the habitat suitability and habitat quality of G. hamiltonii in its global distribution, both in the present and future climatic scenarios. We identify that only 58,542 km2 (7.16%) of the total range extent (817,341 km2) is suitable for this species, along with the fragmented habitats in both the eastern and western ranges. Comparative habitat quality assessment suggests the level of patch shape in the western range is higher (71.3%) compared to the eastern range. Our results suggest a massive decline of approximately 65.73% to 70.31% and 70.53% to 75.30% under ssp245 and ssp585 future scenarios, respectively, for the years between 2021-2040 and 2061-2080 compared with the current distribution. The present study indicates that proper conservation management requires greater attention to the causes and solutions to the fragmented distribution and safeguarding of this endangered species in the Indus, Ganges, and Brahmaputra (IGB) river basins.

SINEs as Credible Signs to Prove Common Ancestry in the Tree of Life: A Brief Review of Pioneering Case Studies in Retroposon Systematics

Currently, the insertions of SINEs (and other retrotransposed elements) are regarded as one of the most reliable synapomorphies in molecular systematics. The methodological mainstream of molecular systematics is the calculation of nucleotide (or amino acid) sequence divergences under a suitable substitution model. In contrast, SINE insertion analysis does not require any complex model because SINE insertions are unidirectional and irreversible. This straightforward methodology was named the “SINE method,” which resolved various taxonomic issues that could not be settled by sequence comparison alone. The SINE method has challenged several traditional hypotheses proposed based on the fossil record and anatomy, prompting constructive discussions in the Evo/Devo era. Here, we review our pioneering SINE studies on salmon, cichlids, cetaceans, Afrotherian mammals, and birds. We emphasize the power of the SINE method in detecting incomplete lineage sorting by tracing the genealogy of specific genomic loci with minimal noise. Finally, in the context of the whole-genome era, we discuss how the SINE method can be applied to further our understanding of the tree of life.

Geometric morphometrics and anatomical network analyses reveal ecospace partitioning among geoemydid turtles from the Uinta Formation, Utah

We present new fossil records of the geoemydid turtle Bridgeremys pusilla from the Uinta Formation of Utah. Turtles are abundant throughout the unit, and known taxa are similar to those from the older strata in the Upper Green River Basin in Wyoming from the Bridger and Washakie Formations. B. pusilla is known from Bridgerian deposits but was not previously known from after the Turtle Bluff Member of the Bridger Formation. The taxon was coveal with two species of the geoemydid Echmatemys (E. callopyge and E. wyomingensis), a common genus of extinct pond turtles known primarily from lacustrine and fluvial deposits in western North America, including the Uinta Basin. In addition to previously documented morphological differences, our geometric morphometric analyses revealed significant differences in epiplastral morphology between B. pusilla and the two coeval Echmatemys species. Bridgeremys pusilla shared several morphological characters with Testudinidae. However, our anatomical network analysis suggests that the carapace of B. pusilla distributed stress forces in a manner more similar to emydids (basal and derived) than to derived testudinoids (Testudinidae and Emydidae), including Echmatemys species. This finding changes our understanding of the ecology of the species and sheds light onto how geoemydid turtles of the Uinta Formation may have partitioned the available ecospace. These new Uintan records extend the geographic range of B. pusilla into the Uinta Basin and stratigraphically through the top of the Uinta Formation, extending the temporal range of the taxon by more than 4 million years through the Uintan North American Land Mammal Age to the base of the Duchesne River Formation.

The phylogenetic relationships of geoemydid turtles from the Eocene Messel Pit Quarry: a first assessment using methods for continuous and discrete characters

The geoemydid turtles of the Eocoene Messel Pit Quarry of Hesse, Germany, are part of a rich Western European fossil record of testudinoids. Originally referred to as "Ocadia" kehreri and "Ocadia" messeliana, their systematic relationships remain unclear. A previous study proposed that a majority of the Western European geoemydids, including the Messel geoemydids, are closely related to the Recent European representatives of the clade Mauremys. Another study hypothesised that the Western European geoemydid fauna is more phylogenetically diverse, and that the Messel geoemydids are closely related to the East Asian turtles Orlitia and Malayemys. Here we present the first quantitative analyses to date that investigate this question. We use continuous characters in the form of ratios to estimate the placement of the Messel geoemydids in a reference tree that was estimated from molecular data. We explore the placement error obtained from that data with maximum likelihood and Bayesian methods, as well as linear parsimony in combination with discrete characters. We find good overall performance with Bayesian and parsimony analyses. Parsimony performs even better when we also incorporated discrete characters. Yet, we cannot pin down the position of the Messel geoemydids with high confidence. Depending on how intraspecific variation of the ratio characters is treated, parsimony favours a placement of the Messel fossils sister to Orlitia borneensis or sister to Geoemyda spengleri, with weak bootstrap support. The latter placement is suspect because G. spengleri is a phylogenetically problematic species with molecular and morphological data. There is even less support for placements within the Mauremys clade.

Estimating the phylogeny of geoemydid turtles (Cryptodira) from landmark data: an assessment of different methods

BACKGROUND

In the last 20 years, a general picture of the evolutionary relationships between geoemydid turtles (ca. 70 species distributed over the Northern hemisphere) has emerged from the analysis of molecular data. However, there is a paucity of good traditional morphological characters that correlate with the phylogeny, which are essential for the robust integration of fossil and molecular data. Part of this problem might be due to intrinsic limitations of traditional discrete characters. Here, we explore the use of continuous data in the form of 3D coordinates of homologous landmarks on the turtle shell for phylogenetic inference and the phylogenetic placement of single species on a scaffold molecular tree. We focus on the performance yielded by sampling the carapace and/or plastral lobes and using various phylogenetic methods.

METHODS

We digitised the landmark coordinates of the carapace and plastron of 42 and 46 extant geoemydid species, respectively. The configurations were superimposed and we estimated the phylogenetic tree of geoemydids with landmark analysis under parsimony, traditional Farris parsimony, unweighted squared-change parsimony, maximum likelihood with a Brownian motion model, and neighbour-joining on a matrix of pairwise Procrustes distances. We assessed the performance of those analyses by comparing the trees against a reference phylogeny obtained from seven molecular markers. For comparisons between trees we used difference measures based on quartets and splits. We used the same reference tree to evaluate phylogenetic placement performance by a leave-one-out validation procedure.

RESULTS

Whatever method we used, similarity to the reference phylogeny was low. The carapace alone gave slightly better results than the plastron or the complete shell. Assessment of the potential for placement of single species on the reference tree with landmark data gave much better results, with similar accuracy and higher precision compared to the performance of discrete characters with parsimony.

Incomplete Lineage Sorting and Hybridization Statistics for Large-Scale Retroposon Insertion Data

Ancient retroposon insertions can be used as virtually homoplasy-free markers to reconstruct the phylogenetic history of species. Inherited, orthologous insertions in related species offer reliable signals of a common origin of the given species. One prerequisite for such a phylogenetically informative insertion is that the inserted element was fixed in the ancestral population before speciation; if not, polymorphically inserted elements may lead to random distributions of presence/absence states during speciation and possibly to apparently conflicting reconstructions of their ancestry. Fortunately, such misleading fixed cases are relatively rare but nevertheless, need to be considered. Here, we present novel, comprehensive statistical models applicable for (1) analyzing any pattern of rare genomic changes, (2) testing and differentiating conflicting phylogenetic reconstructions based on rare genomic changes caused by incomplete lineage sorting or/and ancestral hybridization, and (3) differentiating between search strategies involving genome information from one or several lineages. When the new statistics are applied, in non-conflicting cases a minimum of three elements present in both of two species and absent in a third group are considered significant support (p<0.05) for the branching of the third from the other two, if all three of the given species are screened equally for genome or experimental data. Five elements are necessary for significant support (p<0.05) if a diagnostic locus derived from only one of three species is screened, and no conflicting markers are detected. Most potentially conflicting patterns can be evaluated for their significance and ancestral hybridization can be distinguished from incomplete lineage sorting by considering symmetric or asymmetric distribution of rare genomic changes among possible tree configurations. Additionally, we provide an R-application to make the new KKSC insertion significance test available for the scientific community at http://retrogenomics.uni-muenster.de:3838/KKSC_significance_test/.

The presence/absence patterns of transposed elements, so called jumping genes, provide invaluable information about evolution. Unfortunately, there is still no clear all-encompassing analysis of the statistical significance of insertion patterns, and the single existing model of insertion data is no longer sufficient for the emerging genomic era. Here, we have provided a comprehensive statistical framework for testing the significance of support for phylogenetic hypotheses derived from genome-level presence/absence data such as retroposon insertions and for evaluating such data for different evolutionary scenarios, including polytomy, incomplete lineage sorting, and ancestral hybridization. This statistical framework is especially important for high-throughput applications of current and upcoming genome projects due to its treatment of unlimited numbers of testable markers, and is embedded in a user-friendly R-application available to the scientific community online. Finally, a reliable, adaptable calculation for the significance of support for phylogenetic trees derived from retroposon presence/absence data is now available.

Feeding kinematics of freshwater turtles: what advantage do invasive species possess?

The red-eared slider (Trachemys scripta) is an invasive turtle species that is displacing the populations of native freshwater turtles in many countries. However, the mechanism that makes red-eared sliders superior competitors has been less well studied. In this study, we compare the feeding kinematics of the red-eared slider with those of Reeves' pond turtle (Mauremys reevesii), a turtle native to East Asia, and offer an explanation as to why red-eared sliders are superior in food competition. Reeves' pond turtles sympatric to red-eared sliders have been reported to have a mainly durophagous diet in contrast to the preference for a soft diet in areas of allopatry to red-eared sliders, indicating the dietary shift resulted from food competition. Maximum neck extension and retraction speed are considered to be indicators of striking ability, but were not found to be superior in red-eared sliders. In fact, maximum neck extension speed was significantly higher in Reeves' pond turtles, though this advantage may be counteracted by the longer neck of Reeves' pond turtles, resulting in similar neck extension times. On the other hand, red-eared sliders had a significantly shorter gape cycle time and neck retraction time, indicating that they can complete feeding in a short time. Therefore, red-eared sliders are suggested to be superior in food competition not due to their striking ability, but due to exploiting preferable food in a shorter time.

The population history of endogenous retroviruses in mule deer (Odocoileus hemionus)

Mobile elements are powerful agents of genomic evolution and can be exceptionally informative markers for investigating species and population-level evolutionary history. While several studies have utilized retrotransposon-based insertional polymorphisms to resolve phylogenies, few population studies exist outside of humans. Endogenous retroviruses are LTR-retrotransposons derived from retroviruses that have become stably integrated in the host genome during past infections and transmitted vertically to subsequent generations. They offer valuable insight into host-virus co-evolution and a unique perspective on host evolutionary history because they integrate into the genome at a discrete point in time. We examined the evolutionary history of a cervid endogenous gammaretrovirus (CrERVγ) in mule deer (Odocoileus hemionus). We sequenced 14 CrERV proviruses (CrERV-in1 to -in14), and examined the prevalence and distribution of 13 proviruses in 262 deer among 15 populations from Montana, Wyoming, and Utah. CrERV absence in white-tailed deer (O. virginianus), identical 5' and 3' long terminal repeat (LTR) sequences, insertional polymorphism, and CrERV divergence time estimates indicated that most endogenization events occurred within the last 200000 years. Population structure inferred from CrERVs (F ST = 0.008) and microsatellites (θ = 0.01) was low, but significant, with Utah, northwestern Montana, and a Helena herd being particularly differentiated. Clustering analyses indicated regional structuring, and non-contiguous clustering could often be explained by known translocations. Cluster ensemble results indicated spatial localization of viruses, specifically in deer from northeastern and western Montana. This study demonstrates the utility of endogenous retroviruses to elucidate and provide novel insight into both ERV evolutionary history and the history of contemporary host populations.

Transposable elements: from DNA parasites to architects of metazoan evolution

One of the most unexpected insights that followed from the completion of the human genome a decade ago was that more than half of our DNA is derived from transposable elements (TEs). Due to advances in high throughput sequencing technologies it is now clear that TEs comprise the largest molecular class within most metazoan genomes. TEs, once categorised as "junk DNA", are now known to influence genomic structure and function by increasing the coding and non-coding genetic repertoire of the host. In this way TEs are key elements that stimulate the evolution of metazoan genomes. This review highlights several lines of TE research including the horizontal transfer of TEs through host-parasite interactions, the vertical maintenance of TEs over long periods of evolutionary time, and the direct role that TEs have played in generating morphological novelty.

Universal artifacts affect the branching of phylogenetic trees, not universal scaling laws

BACKGROUND

The superficial resemblance of phylogenetic trees to other branching structures allows searching for macroevolutionary patterns. However, such trees are just statistical inferences of particular historical events. Recent meta-analyses report finding regularities in the branching pattern of phylogenetic trees. But is this supported by evidence, or are such regularities just methodological artifacts? If so, is there any signal in a phylogeny?

METHODOLOGY

In order to evaluate the impact of polytomies and imbalance on tree shape, the distribution of all binary and polytomic trees of up to 7 taxa was assessed in tree-shape space. The relationship between the proportion of outgroups and the amount of imbalance introduced with them was assessed applying four different tree-building methods to 100 combinations from a set of 10 ingroup and 9 outgroup species, and performing covariance analyses. The relevance of this analysis was explored taking 61 published phylogenies, based on nucleic acid sequences and involving various taxa, taxonomic levels, and tree-building methods.

PRINCIPAL FINDINGS

All methods of phylogenetic inference are quite sensitive to the artifacts introduced by outgroups. However, published phylogenies appear to be subject to a rather effective, albeit rather intuitive control against such artifacts. The data and methods used to build phylogenetic trees are varied, so any meta-analysis is subject to pitfalls due to their uneven intrinsic merits, which translate into artifacts in tree shape. The binary branching pattern is an imposition of methods, and seldom reflects true relationships in intraspecific analyses, yielding artifactual polytomies in short trees. Above the species level, the departure of real trees from simplistic random models is caused at least by two natural factors--uneven speciation and extinction rates; and artifacts such as choice of taxa included in the analysis, and imbalance introduced by outgroups and basal paraphyletic taxa. This artifactual imbalance accounts for tree shape convergence of large trees.

SIGNIFICANCE

There is no evidence for any universal scaling in the tree of life. Instead, there is a need for improved methods of tree analysis that can be used to discriminate the noise due to outgroups from the phylogenetic signal within the taxon of interest, and to evaluate realistic models of evolution, correcting the retrospective perspective and explicitly recognizing extinction as a driving force. Artifacts are pervasive, and can only be overcome through understanding the structure and biological meaning of phylogenetic trees. Catalan Abstract in Translation S1.

Developing markers for multilocus phylogenetics in non-model organisms: A test case with turtles

We present a strategy for phylogenetic marker development in non-model systems. Rather than using the traditional approach of comparing distantly related taxa to develop conserved primers for unknown species, we explore an alternative strategy that builds primers directly from a single, relatively well characterized species and applies those primers to increasingly distantly related taxa. We develop and test our protocol with turtles. Using a single BAC end-sequence library consisting of 3461 sequences totaling 2.43 million base pairs of data, we outline a procedure to flag repeat elements, followed by a BLAST approach to categorize sequences into high, low, and no similarity compartments compared to GenBank sequences. We developed and tested a panel of 96 primer pairs with a set of turtle tissues that forms a series of increasingly distantly related taxa with respect to the BAC reference species. Finally, we sequenced 11 of these newly discovered markers across a diverse set of 18 turtle species that spans the 210 million years of chelonian crown-group history and that includes representatives of most of the major clades of extant turtles. Our results indicate that large numbers of new, phylogenetically informative markers can be developed quickly and inexpensively from a single BAC, EST, or similar genomic resource, and that those markers provide reliable phylogenetic information across both shallow and deep levels of phylogenetic history. Our results also highlight the importance of screening for and managing repetitive elements found in randomly sequenced DNA fragments. We presume that our strategy should work well across any similarly divergent clade, suggesting that many-marker datasets can be developed quickly and efficiently for phylogenetic analysis.