Taking advantage of reference-guided assembly in a slowly-evolving lineage: Application to Testudo graeca.

BACKGROUND: Obtaining de novo chromosome-level genome assemblies greatly enhances conservation and evolutionary biology studies. For many research teams, long-read sequencing technologies (that produce highly contiguous assemblies) remain unaffordable or unpractical. For the groups that display high synteny conservation, these limitations can be overcome by a reference-guided assembly using a close relative genome. Among chelonians, tortoises (Testudinidae) are considered one of the most endangered taxa, which calls for more genomic resources. Here we make the most of high synteny conservation in chelonians to produce the first chromosome-level genome assembly of the genus Testudo with one of the most iconic tortoise species in the Mediterranean basin: Testudo graeca. RESULTS: We used high-quality, paired-end Illumina sequences to build a reference-guided assembly with the chromosome-level reference of Gopherus evgoodei. We reconstructed a 2.29 Gb haploid genome with a scaffold N50 of 107.598 Mb and 5.37% gaps. We sequenced 25,998 protein-coding genes, and identified 41.2% of the assembly as repeats. Demographic history reconstruction based on the genome revealed two events (population decline and recovery) that were consistent with previously suggested phylogeographic patterns for the species. This outlines the value of such reference-guided assemblies for phylogeographic studies. CONCLUSIONS: Our results highlight the value of using close relatives to produce de novo draft assemblies in species where such resources are unavailable. Our annotated genome of T. graeca paves the way to delve deeper into the species' evolutionary history and provides a valuable resource to enhance direct conservation efforts on their threatened populations.

Genetic identification, morphology and distribution of Natrixhelvetica subspecies in southern and western Switzerland (Reptilia, Squamata, Serpentes).

Most of Switzerland is inhabited by the nominotypical subspecies of the barred grass snake (Natrixhelveticahelvetica), which is characterized by mitochondrial DNA lineage E. Only in the northeast of the country, the common grass snake (N.natrix) occurs and hybridizes with N.h.helvetica in a narrow contact zone. However, we discovered that in southern and western Switzerland barred grass snakes representing another mtDNA lineage (lineage C) are widely distributed. Lineage C is typical for Alpine populations of the southern subspecies N.h.sicula. Our microsatellite analyses of the Swiss samples revealed differences between the two subspecies and also a substructure with two clusters in each subspecies. Furthermore, we discovered a contact and hybrid zone of N.h.helvetica and N.h.sicula along the northern shore of Lake Geneva and also confirm that interbreeding with alien common grass snakes (N.n.moreotica, mtDNA lineage 7) occurs there. This finding is of concern for nature conservation and measures should be taken to prevent further genetic pollution. Using morphometrics, we found no differences between the two subspecies of N.helvetica, while N.natrix was slightly distinct from N.helvetica.

Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae).

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

Socioeconomic differences in the utilization of diagnostic imaging and non-pharmaceutical conservative therapies for spinal diseases.

BACKGROUND: A different utilization of health care services due to socioeconomic status on the same health plan contradicts the principle of equal treatment. We investigated the presence and magnitude of socioeconomic differences in utilization of diagnostic imaging and non-pharmaceutical conservative therapies for patients with spinal diseases. METHODS: The cohort study based on routine healthcare data from Germany with 11.7 million patient-years between 2012 and 2016 for patients with physician-confirmed spinal diseases (ICD-10: M40-M54), occupation and age 20 to 64 years. A Poisson model estimated the effects of the socioeconomic status (school education, professional education and occupational position) for the risk ratio of receiving diagnostic imaging (radiography, computed tomography, magnetic resonance imaging) and non-pharmaceutical conservative therapies (physical therapy including exercise therapy, manual therapy and massage, spinal manipulative therapy, acupuncture). RESULTS: Patients received diagnostic imaging in 26%, physical therapy in 32%, spinal manipulative therapy in 25%, and acupuncture in 4% of all patient-years. Similar to previous survey-based studies higher rates of utilization were associated with higher socioeconomic status. These differences were most pronounced for manual therapy, exercise therapy, and magnetic resonance imaging. CONCLUSIONS: The observed differences in health care utilization were highly related to socioeconomic status. Socioeconomic differences were higher for more expensive health services. Further research is necessary to identify barriers to equitable access to health services and to take appropriate action to decrease existing social disparities.

Genomics reveals broad hybridization in deeply divergent Palearctic grass and water snakes (Natrix spp.).

Understanding speciation is one of the cornerstones of biological diversity research. Currently, speciation is often understood as a continuous process of divergence that continues until genetic or other incompatibilities minimize or prevent interbreeding. The Palearctic snake genus Natrix is an ideal group to study speciation, as it comprises taxa representing distinct stages of the speciation process, ranging from widely interbreeding parapatric taxa through parapatric species with very limited gene flow in narrow hybrid zones to widely sympatric species. To understand the evolution of reproductive isolation through time, we have sequenced the genomes of all five species within this genus and two additional subspecies. We used both long-read and short-read methods to sequence and de-novo-assemble two high-quality genomes (Natrix h. helvetica, Natrix n. natrix) to their 1.7 Gb length with a contig N50 of 4.6 Mbp and 1.5 Mbp, respectively, and used these as references to assemble the remaining short-read-based genomes. Our phylogenomic analyses yielded a well-supported dated phylogeny and evidence for a surprisingly complex history of interspecific gene flow, including between widely sympatric species. Furthermore, evidence for gene flow was also found for currently allopatric species pairs. Genetic exchange among these well-defined, distinct, and several million-year-old reptile species emphasizes that speciation and maintenance of species distinctness can occur despite continued genetic exchange.

Ancient DNA elucidates the lost world of western Indian Ocean giant tortoises and reveals a new extinct species from Madagascar.

Before humans arrived, giant tortoises occurred on many western Indian Ocean islands. We combined ancient DNA, phylogenetic, ancestral range, and molecular clock analyses with radiocarbon and paleogeographic evidence to decipher their diversity and biogeography. Using a mitogenomic time tree, we propose that the ancestor of the extinct Mascarene tortoises spread from Africa in the Eocene to now-sunken islands northeast of Madagascar. From these islands, the Mascarenes were repeatedly colonized. Another out-of-Africa dispersal (latest Eocene/Oligocene) produced on Madagascar giant, large, and small tortoise species. Two giant and one large species disappeared c. 1000 to 600 years ago, the latter described here as new to science using nuclear and mitochondrial DNA. From Madagascar, the Granitic Seychelles were colonized (Early Pliocene) and from there, repeatedly Aldabra (Late Pleistocene). The Granitic Seychelles populations were eradicated and later reintroduced from Aldabra. Our results underline that integrating ancient DNA data into a multi-evidence framework substantially enhances the knowledge of the past diversity of island faunas.

Skeletal repatterning enhances the protective capacity of the shell in African hinge-back tortoises (Kinixys).

Changes in the structural association of skeletal traits are crucial to the evolution of novel forms and functions. In vertebrates, such rearrangements often occur gradually and may precede or coincide with the functional activation of skeletal traits. To illustrate this process, we examined the ontogeny of African hinge-back tortoises (Kinixys spp.). Kinixys species feature a moveable "hinge" on the dorsal shell (carapace) that enables shell closure (kinesis) when the hind limbs are withdrawn. This hinge, however, is absent in juveniles. Herein, we describe how this unusual phenotype arises via alterations in the tissue configuration and shape of the carapace. The ontogenetic repatterning of osseous and keratinous tissue coincided with shifts in morphological integration and the establishment of anterior (static) and posterior (kinetic) carapacial modules. Based on ex vivo skeletal movement and raw anatomy, we propose that Kinixys employs a "sliding hinge" shell-closing system that overcomes thoracic rigidity and enhances the protective capacity of the carapace. Universal properties of the vertebrate skeleton, such as plasticity, modularity, and secondary maturation processes, contributed to adaptive evolutionary change in Kinixys. We discuss a hypothetical model to explain the delayed emergence of skeletal traits and its relevance to the origins of novel form-to-function relationships.

Dermacentor reticulatus - a tick on its way from glacial refugia to a panmictic Eurasian population.

The ornate dog tick (Dermacentor reticulatus) shows a recently expanding geographic distribution. Knowledge on its intraspecific variability, population structure, rate of genetic diversity and divergence, including its evolution and geographic distribution, is crucial to understand its dispersal capacity. All such information would help to evaluate the potential risk of future spread of associated pathogens of medical and veterinary concern. A set of 865 D. reticulatus ticks was collected from 65 localities across 21 countries, from Portugal in the west to Kazakhstan and southern Russia in the east. Cluster analyses of 16 microsatellite loci were combined with nuclear (ITS2, 18S) and mitochondrial (12S, 16S, COI) sequence data to uncover the ticks' population structures and geographical patterns. Approximate Bayesian computation was applied to model evolutionary relationships among the found clusters. Low variability and a weak phylogenetic signal showing an east-west cline were detected both for mitochondrial and nuclear sequence markers. Microsatellite analyses revealed three genetic clusters, where the eastern and western cluster gradient was supplemented by a third, northern cluster. Alternative scenarios could explain such a tripartite population structure by independent formation of clusters in separate refugia, limited gene flow connected with isolation by distance causing a "bipolar pattern", and the northern cluster deriving from admixture between the eastern and western populations. The best supported demographic scenario of this tick species indicates that the northern cluster derived from admixture between the eastern and western populations 441 (median) to 224 (mode) generations ago, suggesting a possible link with the end of the Little Ice Age in Europe.

Genotyping the phenotypic diversity in Aegean Natrixnatrixmoreotica (Bedriaga, 1882) (Reptilia, Serpentes, Natricidae).

We examined the mitochondrial identity of Aegean Natrixnatrixmoreotica representing different morphotypes, with a focus on new material from Milos and Skyros. We found no correlation between distinct morphotypes and mitochondrial identity. Our results support that grass snake populations are polyphenetic and that southern subspecies, including island populations, show a higher variability than northern ones.

Indication for spinal surgery: associated factors and regional differences in Germany.

BACKGROUND: Rising surgery rates have raised questions about the indications for spinal surgery. The study investigated patient-level and regional factors associated with spinal surgery for patients with spinal diseases. METHODS: We undertook a cohort study based on routine healthcare data from Germany of 18.4 million patients within 60.9 million episodes of two patient-years before a possible spinal surgery in the time period 2008 to 2016. Using a Poisson model, the effects of a broad range of patient-related (sociodemographic, morbidity, social status), disease- and healthcare-related (physicians' specialty, conservative treatments) and regional variables were analyzed. RESULTS: There was substantial regional heterogeneity in the occurrence of spinal surgery which decreased by only one quarter when controlling for the various determinants assessed. Previous musculoskeletal and mental health disorders as well as physical therapy were associated with a lower probability of surgery in the fully-adjusted model. Prescriptions for pain medication and consultations of specialists were associated with a higher probability of surgery. However, the specific severity of the vertebral diseases could not be taken into account in the analysis. Furthermore, a substantial proportion of patients with surgery did not receive a consultation with an outpatient specialist (29.5%), preoperative diagnostics (37.0%) or physical therapy (48.3%) before hospital admission. CONCLUSION: This large study on spinal diseases in Germany highlights important patterns in medical care of spinal diseases and their association with the probability of spinal surgery. However, only a relatively small proportion of the regional heterogeneity in spinal surgery could be explained by the extensive consideration of confounders, which suggests the relevance of other unmeasured factors like physicians' preferences.

Disentangling the Pelodiscus axenaria complex, with the description of a new Chinese species and neotype designation for P. axenaria (Zhou, Zhang amp; Fang, 1991).

We describe a new species from the Pelodiscus axenaria complex from Hunan and Jiangxi Provinces, China. Also, the application of the name P. axenaria (Zhou, Zhang Fang, 1991) is clarified by designating a neotype for this species. Besides its genetic divergence, the new species differs from all other Pelodiscus species, including the two other taxa constituting the complex (P. axenaria, P. huangshanensis), in the following combination of morphological traits: (1) small adult size, <15 cm carapace length; (2) carapace distinctly keeled, more or less strongly tuberculated, usually olive clay-coloured and adorned with greenish black marbling; (3) plastron yellowish white, typically immaculate except for a blurred-edged blotch behind each axilla that does not extend to the entoplastron and some slight black suffusion along its anterior border; (4) underside of the leathery margin of the carapace with varying amounts of dark pigmentation; (5) head olive clay-coloured with numerous black splotches; (6) chin grey brown with pale stipples, throat dark grey, finely spotted with black; (7) neck with a wide yellow lateral band stretching from the tympanum posteriorly, which tends to fade with age; (8) entoplastron boomerang-shaped, the amount of bending of the transverse bar between the two posteriolaterally directed rami >90.

"Ancient DNA" reveals that the scientific name for an extinct tortoise from Cape Verde refers to an extant South American species.

We examined the type material of the extinct tortoise species Geochelone atlantica López-Jurado, Mateo and García-Márquez, 1998 from Sal Island, Cape Verde, using aDNA approaches and AMS radiocarbon dating. High-quality mitochondrial genomes obtained from the three type specimens support that all type material belongs to the same individual. In phylogenetic analyses using mitochondrial genomes of all species groups and genera of extant and some recently extinct tortoises, the type material clusters within the extant South American red-footed tortoise Chelonoidis carbonarius (Spix, 1824). Our radiocarbon date indicates that the tortoise from which the type series of G. atlantica originates was still alive during 1962 and 1974. These results provide firm evidence that the type material of G. atlantica does not belong to the Quaternary tortoise bones excavated on Sal Island in the 1930s, as originally thought. Thus, the extinct tortoise species remains unstudied and lacks a scientific name, and the name G. atlantica has to be relegated into the synonymy of C. carbonarius. The circumstances that led to this confusion currently cannot be disentangled.

Population structure and gene flow of the syntopic turtles Emys and Mauremys from coastal and inland regions of Anatolia (Turkey): results from mitochondrial and microsatellite data.

Revealing the genetic basis of the existence of different species living together in different geographic regions provides clarification of this phylogeographic differentiation. In this study, we investigated the population genetics and evaluated the level of genetic variation of inland and coastal populations of Mauremys and Emys in Turkey. Tissue samples of 196 terrapins were studied which were collected from syntopic coastal (Gölbent-Söke/Aydin; M. rivulata and E. orbicularis) and inland populations (Bahçesaray/Aksaray; M. caspica and E. orbicularis). DNA was isolated using the InnuPREP DNA Mini Kit. Mitochondrial DNA sequences and allelic variation at 13 microsatellite loci for Mauremys and 12 microsatellite loci for Emys were examined.  Three haplotypes were found for Emys orbicularis (Im, Ip and Iw) collected from the coastal region and two haplotypes for Emys orbicularis (Ig and Im) collected from inland. Two haplotypes were identified for M. caspica (Cmt8 and Cmt9) and three haplotypes were identified for M. rivulata (Rmt3, Rmt24 and Rmt26). Using microsatellites and the software STRUCTURE the most probable value for K was revealed as two 2 for both species. The FST value between M. rivulata and M. caspica was 0.39, and between the coastal and inland populations of E. orbicularis 0.09. It can be concluded that Emys populations tend to evolve by somehow preserving the allelic richness they have and Mauremys populations continue to differentiate so that new species emerge in the evolutionary process to reach the ideal allelic structure.

Biogeography of Italy revisited: genetic lineages confirm major phylogeographic patterns and a pre-Pleistocene origin of its biota.

Warm-adapted species survived the cold stages of the past glacial-interglacial cycles in southern European peninsulas and recolonized major parts of Central and Northern Europe in the wake of postglacial warming. However, many of the genetic lineages which differentiated within these refugia predate even the Pleistocene. One of these refugia is the Italian Peninsula with its surrounding islands. In this meta-analysis, we compiled phylogeographic patterns of multiple species across this region. We transformed existing geographic information on 78 animal and plant species (with a total of 471 genetic lineages) within 17 predefined areas into presence/absence matrices. We elaborated three matrices: (i) only old lineages exclusively distinguished by deep splits, (ii) only young lineages distinguished by shallow (i.e. recent) splits, and (iii) presence/absence of the respective species. To infer biogeographic relationships between the predefined areas, we performed bootstrapped neighbour joining cluster analyses on these three matrices. In addition, we reviewed the geological history of Italy to identify causes of the observed biogeographic patterns. We found Sardinia and Corsica to be biogeographically closely linked with each other, and that they diverge strongly from all other regions. Sicily also diverges strongly from all other regions, while the intra-island differentiation was comparatively low. On the Italian mainland, Calabria exhibited the most pronounced biogeographic differentiation, often with several lineages present, resulting from old vicariance events within the region. Furthermore, southern Apulia and the Po Plain with adjoining areas of northern peninsular Italy displayed considerable distinctiveness. Admixture prevailed in the areas between these three regions. The ancient isolation of Sicily, as well as Sardinia plus Corsica, resulted in endemic lineages with only moderate recent exchange with adjacent mainland regions. Pronounced diversification occurs within the Italian Peninsula. The complex tectonic activities, such as shifting (micro)plates, submergence of major parts of peninsular Italy with the genesis of numerous Pliocene islands, in combination with the climatic cycles during the Pleistocene have most likely generated the current biogeographic pattern of species. Immigrations from the Balkan Peninsula into northern Italy partly accounted for the distinctiveness of this region.

Evolutionary history of mental glands in turtles reveals a single origin in an aquatic ancestor and recurrent losses independent of macrohabitat.

Despite the relevance of chemical communication in vertebrates, comparative examinations of macroevolutionary trends in chemical signaling systems are scarce. Many turtle and tortoise species are reliant on chemical signals to communicate in aquatic and terrestrial macrohabitats, and many of these species possess specialized integumentary organs, termed mental glands (MGs), involved in the production of chemosignals. We inferred the evolutionary history of MGs and tested the impact of macrohabitat on their evolution. Inference of ancestral states along a time-calibrated phylogeny revealed a single origin in the ancestor of the subclade Testudinoidea. Thus, MGs represent homologous structures in all descending lineages. We also inferred multiple independent losses of MGs in both terrestrial and aquatic clades. Although MGs first appeared in an aquatic turtle (the testudinoid ancestor), macrohabitat seems to have had little effect on MG presence or absence in descendants. Instead, we find clade-specific evolutionary trends, with some clades showing increased gland size and morphological complexity, whereas others exhibiting reduction or MG loss. In sister clades inhabiting similar ecological niches, contrasting patterns (loss vs. maintenance) may occur. We conclude that the multiple losses of MGs in turtle clades have not been influenced by macrohabitat and that other factors have affected MG evolution.

Ancient mitogenomics elucidates diversity of extinct West Indian tortoises.

We present 10 nearly complete mitochondrial genomes of the extinct tortoise Chelonoidis alburyorum from the Bahamas. While our samples represent morphologically distinct populations from six islands, their genetic divergences were shallow and resembled those among Galápagos tortoises. Our molecular clock estimates revealed that divergence among Bahamian tortoises began ~ 1.5 mya, whereas divergence among the Galápagos tortoises (C. niger complex) began ~ 2 mya. The inter-island divergences of tortoises from within the Bahamas and within the Galápagos Islands are much younger (0.09-0.59 mya, and 0.08-1.43 mya, respectively) than the genetic differentiation between any other congeneric pair of tortoise species. The shallow mitochondrial divergences of the two radiations on the Bahamas and the Galápagos Islands suggest that each archipelago sustained only one species of tortoise, and that the taxa currently regarded as distinct species in the Galápagos should be returned to subspecies status. The extinct tortoises from the Bahamas have two well-supported clades: the first includes one sample from Great Abaco and two from Crooked Island; the second clade includes tortoises from Great Abaco, Eleuthera, Crooked Island, Mayaguana, Middle Caicos, and Grand Turk. Tortoises belonging to both clades on Great Abaco and Crooked Island suggest late Holocene inter-island transport by prehistoric humans.

Cytogenetic Analysis of the Asian Box Turtles of the Genus Cuora (Testudines, Geoemydidae).

The Asian box turtle genus Cuora currently comprises 13 species with a wide distribution in Southeast Asia, including China and the islands of Indonesia and Philippines. The populations of these species are rapidly declining due to human pressure, including pollution, habitat loss, and harvesting for food consumption. Notably, the IUCN Red List identifies almost all species of the genus Cuora as Endangered (EN) or Critically Endangered (CR). In this study, we explore the karyotypes of 10 Cuora species with conventional (Giemsa staining, C-banding, karyogram reconstruction) and molecular cytogenetic methods (in situ hybridization with probes for rDNA loci and telomeric repeats). Our study reveals a diploid chromosome number of 2n = 52 chromosomes in all studied species, with karyotypes of similar chromosomal morphology. In all examined species, rDNA loci are detected at a single medium-sized chromosome pair and the telomeric repeats are restricted to the expected terminal position across all chromosomes. In contrast to a previous report, sex chromosomes are neither detected in Cuoragalbinifrons nor in any other species. Therefore, we assume that these turtles have either environmental sex determination or genotypic sex determination with poorly differentiated sex chromosomes. The conservation of genome organization could explain the numerous observed cases of interspecific hybridization both within the genus Cuora and across geoemydid turtles.

Turtles and Tortoises Are in Trouble.

Turtles and tortoises (chelonians) have been integral components of global ecosystems for about 220 million years and have played important roles in human culture for at least 400,000 years. The chelonian shell is a remarkable evolutionary adaptation, facilitating success in terrestrial, freshwater and marine ecosystems. Today, more than half of the 360 living species and 482 total taxa (species and subspecies combined) are threatened with extinction. This places chelonians among the groups with the highest extinction risk of any sizeable vertebrate group. Turtle populations are declining rapidly due to habitat loss, consumption by humans for food and traditional medicines and collection for the international pet trade. Many taxa could become extinct in this century. Here, we examine survival threats to turtles and tortoises and discuss the interventions that will be needed to prevent widespread extinction in this group in coming decades.