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).

Genetic assessments and parentage analysis of captive Bolson tortoises (Gopherus flavomarginatus) inform their "rewilding" in New Mexico

The Bolson tortoise (Gopherus flavomarginatus) is the first species of extirpated megafauna to be repatriated into the United States. In September 2006, 30 individuals were translocated from Arizona to New Mexico with the long-term objective of restoring wild populations via captive propagation. We evaluated mtDNA sequences and allelic diversity among 11 microsatellite loci from the captive population and archived samples collected from wild individuals in Durango, Mexico (n = 28). Both populations exhibited very low genetic diversity and the captive population captured roughly 97.5% of the total wild diversity, making it a promising founder population. Genetic screening of other captive animals (n = 26) potentially suitable for reintroduction uncovered multiple hybrid G. flavomarginatus×G. polyphemus, which were ineligible for repatriation; only three of these individuals were verified as purebred G. flavomarginatus. We used these genetic data to inform mate pairing, reduce the potential for inbreeding and to monitor the maintenance of genetic diversity in the captive population. After six years of successful propagation, we analyzed the parentage of 241 hatchlings to assess the maintenance of genetic diversity. Not all adults contributed equally to successive generations. Most yearly cohorts of hatchlings failed to capture the diversity of the parental population. However, overlapping generations of tortoises helped to alleviate genetic loss because the entire six-year cohort of hatchlings contained the allelic diversity of the parental population. Polyandry and sperm storage occurred in the captives and future management strategies must consider such events.

The genetic crisis of the Mexican Bolson Tortoise (Gopherus flavomarginatus: Testudinidae)

Gopherusflavomarginatus(Testudinidae) is endemic to a series of discontinuous, isolated basins collectively known as the Bolson de Mapimí in the Chihuahuan Desert. Its numbers declined after catastrophic levels of exploitation during the mid-20th century. However currently, the Bolson Tortoise appears to be on a path to recovery owing to intensive, sustained conservation efforts. We sequenced an 842-bp-long fragment of the D-loop from 76 individuals distributed throughout the species’ range. The results revealed only two haplotypes. An AMOVA showed that 95% of the variance occurred among populations, whereas the remaining 5% was explained by genetic differences within populations. Tectonic processes together with ecological transformation during the Pleistocene and Holocene may be responsible for the reduction in this species’ genetic variation. A bottleneck during which a significant percentage of the haplotype diversity was lost would result in genetic homogeneity. Although there is demographic growth, the lack of genetic diversity is indicative of the potential crisis that the Bolson Tortoise is facing, and awareness must be brought to this situation.