Tracing Genetic Lineages of Captive Desert Tortoises in Arizona

Abstract: We genotyped 180 captive desert tortoises (Gopherus agassizii) from Kingman (n = 45), Phoenix (n = 113), and Tucson (n = 22), Arizona, USA, to determine if the genetic lineage of captives is associated with that of wild tortoises in the local area (Sonoran Desert). We tested all samples for 16 short tandem repeats and sequenced 1,109 base pairs of mitochondrial DNA (mtDNA). To determine genetic origin, we performed assignment tests against a reference database of 997 desert tortoise samples collected throughout the Mojave and Sonoran Deserts. We found that >40% of our Arizona captive samples were genetically of Mojave Desert or hybrid origin, with the percentage of individuals exhibiting the Mojave genotype increasing as the sample locations approached the California, USA, border. In Phoenix, 11.5% were Sonoran–Mojave crosses, and 8.8% were hybrids between desert tortoise and Texas tortoise (G. berlandieri). Our findings present many potential implications for wild tortoises in the Sonoran Desert of Arizona. Escaped or released captive tortoises with Mojave or hybrid genotypes have the potential to affect the genetic composition of Sonoran wild populations. Genotyping captive desert tortoises could be used to inform the adoption process, and thereby provide additional protection to native desert-tortoise populations in Arizona.     

The evolution of different maternal investment strategies in two closely related desert vertebrates

We compared egg size phenotypes and tested several predictions from the optimal egg size (OES) and bet‐hedging theories in two North American desert‐dwelling sister tortoise taxa, Gopherus agassizii and G. morafkai, that inhabit different climate spaces: relatively unpredictable and more predictable climate spaces, respectively. Observed patterns in both species differed from the predictions of OES in several ways. Mean egg size increased with maternal body size in both species. Mean egg size was inversely related to clutch order in G. agassizii, a strategy more consistent with the within‐generation hypothesis arising out of bet‐hedging theory or a constraint in egg investment due to resource availability, and contrary to theories of density dependence, which posit that increasing hatchling competition from later season clutches should drive selection for larger eggs. We provide empirical evidence that one species, G. agassizii, employs a bet‐hedging strategy that is a combination of two different bet‐hedging hypotheses. Additionally, we found some evidence for G. morafkai employing a conservative bet‐hedging strategy. (e.g., lack of intra‐ and interclutch variation in egg size relative to body size). Our novel adaptive hypothesis suggests the possibility that natural selection favors smaller offspring in late‐season clutches because they experience a more benign environment or less energetically challenging environmental conditions (i.e., winter) than early clutch progeny, that emerge under harsher and more energetically challenging environmental conditions (i.e., summer). We also discuss alternative hypotheses of sexually antagonistic selection, which arise from the trade‐offs of son versus daughter production that might have different optima depending on clutch order and variation in temperature‐dependent sex determination (TSD) among clutches. Resolution of these hypotheses will require long‐term data on fitness of sons versus daughters as a function of incubation environment, data as yet unavailable for any species with TSD.     

Are captive tortoises a reservoir for conservation? An assessment of genealogical affiliation of captive Gopherus agassizii to local, wild populations

The conservation of tortoises poses a unique situation because several threatened species are commonly kept as pets within their native ranges. Thus, there is potential for captive populations to be a reservoir for repatriation efforts. We assess the utility of captive populations of the threatened Agassiz’s desert tortoise (Gopherus agassizii) for recovery efforts based on genetic affinity to local areas. We collected samples from 130 captive desert tortoises from three desert communities: two in California (Ridgecrest and Joshua Tree) and the Desert Tortoise Conservation Center (Las Vegas) in Nevada. We tested all samples for 25 short tandem repeats and sequenced 1,109 bp of the mitochondrial genome. We compared captive genotypes to a database of 1,258 Gopherus samples, including 657 wild caught G. agassizii spanning the full range of the species. We conducted population assignment tests to determine the genetic origins of the captive individuals. For our total sample set, only 44 % of captive individuals were assigned to local populations based on genetic units derived from the reference database. One individual from Joshua Tree, California, was identified as being a Morafka’s desert tortoise, G. morafkai, a cryptic species which is not native to the Mojave Desert. Our data suggest that captive desert tortoises kept within the native range of G. agassizii cannot be presumed to have a genealogical affiliation to wild tortoises in their geographic proximity. Precautions should be taken before considering the release of captive tortoises into the wild as a management tool for recovery.     

Evidence for concerted movement of nuclear and mitochondrial clines in a lizard hybrid zone

Moving hybrid zones provide compelling examples of evolution in action, yet long‐term studies that test the assumptions of hybrid zone stability are rare. Using replicated transect samples collected over a 10‐year interval from 2002 to 2012, we find evidence for concerted movement of genetic clines in a plateau fence lizard hybrid zone (Sceloporus tristichus) in Arizona. Cline‐fitting analyses of SNP and mtDNA data both provide evidence that the hybrid zone shifted northward by approximately 2 km during the 10‐year interval. For each sampling period, the mtDNA cline centre is displaced from the SNP cline centre and maintaining an introgression distance of approximately 3 km. The northward expansion of juniper trees into the Little Colorado River Basin in the early 1900s provides a plausible mechanism for hybrid zone formation and movement, and a broadscale quantification of recent land cover change provides support for increased woody species encroachment at the southern end of the hybrid zone. However, population processes can also contribute to hybrid zone movement, and the current stability of the ecotone habitats in the centre of the hybrid zone suggests that movement could decelerate in the future.     

Chronic disease in the Mojave desert tortoise: Host physiology and recrudescence obscure patterns of pathogen transmission

A seminatural, factorial‐design experiment was used to quantify dynamics of the pathogen Mycoplasma agassizii and upper respiratory tract disease in the Mojave desert tortoise (Gopherus agassizii) over 2 years. Groups of initially healthy animals were separated into serologically positive (seropositive), seronegative, and artificially infected groups and paired into 23 pens. We found no evidence of long‐term immune protection to M. agassizii or of immunological memory. Initially seronegative, healthy tortoises experienced an equal amount of disease when paired with other seronegative groups as when paired with seropositive and artificially infected groups—suggesting that recrudescence is as significant as transmission in introducing disease in individuals in this host–pathogen system. Artificially infected groups of tortoises showed reduced levels of morbidity when paired with initially seronegative animals—suggesting either a dilution effect or a strong effect of pathogen load in this system. Physiological dynamics within the host appear to be instrumental in producing morbidity, recrudescence, and infectiousness, and thus of population‐level dynamics. We suggest new avenues for studying diseases in long‐lived ectothermic vertebrates and a shift in modeling such diseases.     

Genomic replacement of native Cobitis lutheri with introduced C. tetralineata through a hybrid swarm following the artificial connection of river systems

River connections via artificial canals will bring about secondary contacts between previously isolated fish species. Here, we present a genetic consequence of such a secondary contact between Cobitis fish species, C. lutheri in the Dongjin River, and C. tetralineata in the Seomjin River in Korea. The construction of water canals about 80 years ago has unidirectionally introduced C. tetralineata into the native habitat of C. lutheri, and then these species have hybridized in the main stream section of the Dongjin River. According to the divergence population genetic analyses of DNA sequence data, the two species diverged about 3.3 million years ago, which is interestingly coincident with the unprecedented paleoceanographic change that caused isolations of the paleo‐river systems in northeast Asia due to sea‐level changes around the late Pliocene. Multilocus genotypic data of nine microsatellites and three nuclear loci revealed an extensively admixed structure in the hybrid zone with a high proportion of various post‐F1 hybrids. Surprisingly, pure native C. lutheri was absent in the hybrid zone in contrast to the 7% of pure C. tetralineata. Such a biased proportion must have resulted from the dominant influence of continually introducing C. tetralineata on the native C. lutheri which has no supply of natives from other tributaries to the hybrid zone due to numerous low‐head dams. In addition, mating experiments indicated that there is no discernible reproductive isolation between them. All the results suggest that the gene pool of native C. lutheri is being rapidly replaced by that of continually introducing C. tetralineata through a hybrid swarm for the last 80 years, which will ultimately lead to the genomic extinction of natives in this hybrid zone.     

Plant and habitat use by Black‐breasted Pufflegs (Eriocnemis nigrivestis), a critically endangered hummingbird

Black‐breasted Pufflegs (Eriocnemis nigrivestis) are hummingbirds endemic to Ecuador and are considered critically endangered because of their limited distribution, a population estimated at fewer than 1000 individuals, and ongoing habitat degradation. From November 2013 to June 2016, we examined the foraging preferences of these hummingbirds using a combination of direct observations, time‐lapse cameras, and motion‐detection software. We first identified 21 species of ornithophilous plants distributed among five sites in the northwestern flanks of the Pichincha volcano in northwest Ecuador. We then monitored these plant species using time‐lapse cameras and recorded 144 visits by Black‐breasted Pufflegs to seven of the 21 species. Most visits (128 of 144 visits, 89%) were to just two species of plants, Macleania rupestris and Palicourea fuchsioides, the latter of which is also an endemic and threatened species. In addition, Black‐breasted Pufflegs were only observed in the most pristine habitats. Given the potential negative effects of climate change for species in the tropical Andes plus the possible loss and degradation of habitat resulting from human activities, efforts are needed to conserve habitats currently used by Black‐breasted Pufflegs, recover degraded habitats, and connect isolated patches of suitable habitat. Our results concerning species of flowering plants used most by Black‐breasted Pufflegs (P. fuchsioides and M. rupestris) should help guide any habitat restoration initiatives.     

Reproductive isolation and cryptic introgression in a sky island enclave of Appalachian birds

Reproductive isolation is central to the speciation process, and cases where the strength of reproductive isolation varies geographically can inform our understanding of speciation mechanisms. Although generally treated as separate species, Black‐capped chickadees (Poecile atricapillus) and Carolina chickadees (P. carolinensis) hybridize and undergo genetic introgression in many areas where they come into contact across the eastern United States and in the northern Appalachian Mountains. The Great Smoky Mountains harbor the last large breeding population of atricapillus in the southern Appalachians, isolated from the species’ main range by nearly 200 km. This population is believed to be reproductively isolated from local carolinensis due to an unusual, behaviorally mediated elevational range gap, which forms during the breeding season and may function as an incipient reproductive isolating mechanism. We examined the effectiveness of this putative isolating mechanism by looking for genetic introgression from carolinensis in Great Smoky Mountain atricapillus. We characterized this population and parental controls genetically using hundreds of amplified fragment length polymorphism (AFLP) loci as well as mitochondrial DNA (mtDNA) sequence data from cytochrome‐b. Great Smoky Mountain atricapillus have experienced nuclear genetic introgression from carolinensis, but at much lower levels than other populations near the hybrid zone to the north. No mitochondrial introgression was detected, in contrast to northern contact areas. Thus, the seasonal elevational range gap appears to have been effective in reducing gene flow between these closely related taxa.     

Phylogeography and population genetics of pine butterflies: Sky islands increase genetic divergence

The sky islands of southeastern Arizona (AZ) mark a major transition zone between tropical and temperate biota and are considered a neglected biodiversity hotspot. Dispersal ability and host plant specificity are thought to impact the history and diversity of insect populations across the sky islands. We aimed to investigate the population structure and phylogeography of two pine‐feeding pierid butterflies, the pine white (Neophasia menapia) and the Mexican pine white (Neophasia terlooii), restricted to these “islands” at this transition zone. Given their dependence on pines as the larval hosts, we hypothesized that habitat connectivity affects population structure and is at least in part responsible for their allopatry. We sampled DNA from freshly collected butterflies from 17 sites in the sky islands and adjacent high‐elevation habitats and sequenced these samples using ddRADSeq. Up to 15,399 SNPs were discovered and analyzed in population genetic and phylogenetic contexts with Stacks and pyRAD pipelines. Low genetic differentiation in N. menapia suggests that it is panmictic. Conversely, there is strong evidence for population structure within N. terlooii. Each sky island likely contains a population of N. terlooii, and clustering is hierarchical, with populations on proximal mountains being more related to each other. The N. menapia habitat, which is largely contiguous, facilitates panmixia, while the N. terlooii habitat, restricted to the higher elevations on each sky island, creates distinct population structure. Phylogenetic results corroborate those from population genetic analyses. The historical climate‐driven fluxes in forest habitat connectivity have implications for understanding the biodiversity of fragmented habitats.     

Population genetics and sympatric divergence of the freshwater gudgeon,Gobiobotia filifer,in the Yangtze River inferred from mitochondrial DNA

The ecosystem and Pleistocene glaciations play important roles in population demography. The freshwater gudgeon, Gobiobotia filifer, is an endemic benthic fish in the Yangtze River and is a good model for ecological and evolutionary studies. This study aimed to decode the population structure of G. filifer in the Yangtze River and reveal whether divergence occurred before or after population radiation. A total of 292 specimens from eight locations in the upper and middle reaches of the Yangtze River were collected from 2014 to 2016 and analyzed via mitochondrial DNA Cyt b gene sequencing. A moderately high level of genetic diversity was found without structures among the population. However, phylogenetic and network topology showed two distinct haplotype groups, and each group contained a similar proportion of individuals from all sampled sites. This suggested the existence of two genetically divergent source populations in G. filifer. We deduced that a secondary contact of distinct glacial refugia was the main factor creating sympatric populations of G. filifer, and climate improvement promoted population expansion and colonization.     

Mito‐nuclear discordance across a recent contact zone for California voles

To examine the processes that maintain genetic diversity among closely related taxa, we investigated the dynamics of introgression across a contact zone between two lineages of California voles (Microtus californicus). We tested the prediction that introgression of nuclear loci would be greater than that for mitochondrial loci, assuming ongoing gene flow across the contact zone. We also predicted that genomic markers would show a mosaic pattern of differentiation across this zone, consistent with genomes that are semi‐permeable. Using mitochondrial cytochrome b sequences and genome‐wide loci developed via ddRAD‐seq, we analyzed genetic variation for 10 vole populations distributed along the central California coast; this transect included populations from within the distributions of both parental lineages as well as the putative contact zone. Our analyses revealed that (1) the two lineages examined are relatively young, having diverged ca. 8.5–54 kya, (2) voles from the contact zone in Santa Barbara County did not include F1 or early generation backcrossed individuals, and (3) there appeared to be little to no recurrent gene flow across the contact zone. Introgression patterns for mitochondrial and nuclear markers were not concordant; only mitochondrial markers revealed evidence of introgression, putatively due to historical hybridization. These differences in genetic signatures are intriguing given that the contact zone occurs in a region of continuous vole habitat, with no evidence of past or present physical barriers. Future studies that examine specific isolating mechanisms, such as microhabitat use and mate choice, will facilitate our understanding of how genetic boundaries are maintained in this system.     

Implications of Anthropogenic Landscape Change on Inter-Population Movements of the Desert Tortoise (Gopherus agassizii)

In the Sonoran Desert of North America, populations of the desert tortoise (Gopherus agassizii) occur in rocky foothills throughout southwestern Arizona and northwestern Mexico. Although tortoise populations appear to be isolated from each other by low desert valleys, individuals occasionally move long distances between populations. Increasingly, these movements are hindered by habitat fragmentation due to anthropogenic landscape changes. We used molecular techniques and radiotelemetry to examine movement patterns of desert tortoises in southern Arizona. We collected blood samples from 170 individuals in nine mountain ranges and analyzed variability in seven microsatellite loci to determine genetic differentiation among populations. Gene flow estimates between populations indicate that populations exchanged individuals historically at a rate greater than one migrant per generation, and positive correlation between genetic and geographic distance of population pairs suggests that the limiting factor for gene flow among populations is isolation by distance. Life history traits of the desert tortoise, a long-lived species with delayed sexual maturity, may severely constrain the ability of small populations to respond to disturbances that increase adult mortality. Historic gene flow estimates among populations suggests that recovery of declining populations may rely heavily on the immigration of new individuals from adjacent mountain ranges. Management strategies compatible with the evolutionary history of gene flow among disjunct populations will help ensure the long-term persistence of Sonoran desert tortoise populations.     

Population genetic structure,migration, and polyploidy origin of a medicinal species Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum, a member of family Cucurbitaceae, is a perennial creeping herb used as a traditional medicinal plant in China. In this study, six polymorphic nSSR and four EST‐SSR primers were used to genotype 1,020 individuals in 72 wild populations of G. pentaphyllum. The genetic diversity and population structure were investigated, and ecological niche modeling was performed to reveal the evolution and demographic history of its natural populations. The results show that G. pentaphyllum has a low level of genetic diversity and high level of variation among populations because of pervasive asexual propagation, genetic drift, and long‐term habitat isolation. Results of the Mantel test demonstrate that the genetic distance and geographic distance are significantly correlated among G. pentaphyllum natural populations. The populations can be divided into two clusters on the basis of genetic structure. Asymmetrical patterns of historical gene flow were observed among the clusters. For the contemporary, almost all the bidirectional gene flow of the related pairs was symmetrical with slight differences. Recent bottlenecks were experienced by 34.72% of the studied populations. The geographic range of G. pentaphyllum continues to expand northward and eastward from Hengduan Mountains. The present distribution of G. pentaphyllum is a consequence of its complex evolution. Polyploidy in G. pentaphyllum is inferred to be polygenetic. Finally, G. pentaphyllum is a species in need of protection, so in situ and ex situ measures should be considered in the future.     

Formation of a recent hybrid zone offers insight into the geographic puzzle and maintenance of species boundaries in musk turtles

Speciation is the result of an accumulation of reproductive barriers between populations, but pinpointing the factors involved is often difficult. However, hybrid zones can form when these barriers are not complete, especially when lineages come into contact in intermediate or modified habitats. We examine a hybrid zone between two closely related riverine turtle species, Sternotherus depressus and S. peltifer, and use dual‐digest RAD sequencing to understand how this hybrid zone formed and elucidate genomic patterns of reproductive isolation. First, the geographical extent and timing of formation of the hybrid zone is established to provide context for understanding the role of extrinsic and intrinsic reproductive isolating mechanisms in this system. The strength of selection on taxon‐specific contributions to maintenance of the hybrid zone is then inferred using a Bayesian genomic cline model. These analyses identify a role for selection inhibiting introgression in some genomic regions at one end of the hybrid zone and promoting introgression in many loci at the other. When selective pressures necessary to generate outliers to the genomic cline are considered with the geographical and temporal context of this hybrid zone, we conclude that habitat‐specific selection probably limits introgression from S. depressus to S. peltifer in the direction of river flow. However, selection is mediating rapid, unidirectional introgression from S. peltifer to S. depressus, which is probably facilitated by anthropogenic habitat alteration. These findings indicate a potentially imminent threat of population‐level genomic extinction for an already imperiled species due to ongoing human‐caused habitat alteration.     

Estimating the probability of illegal desert tortoise collection in the Sonoran Desert

The expansion of road networks in desert tortoise (Gopherus agassizii) habitat in the Sonoran Desert has raised questions concerning appropriate mitigation to reduce impacts at the population level. Although some effects, namely road-kill and habitat loss, have been well documented, illegal tortoise collection has been insufficiently addressed. It has become increasingly important for wildlife and land-use managers to understand the cumulative impacts of roads on tortoises and the effect that those impacts have on population persistence. We estimated the probability of desert tortoise detection and collection along 2-lane paved, maintained gravel, and non-maintained gravel roads to evaluate whether collection probabilities were related to road type. Although collection probability did not vary by road type, the probability of desert tortoise detection by passing motorists was greatest on maintained gravel roads and fewest on non-maintained gravel and paved roads. These results have implications for effectively mitigating the impacts of roads on desert tortoises. Published 2011. This article is a U.S. Government work and is in the public domain in the USA.     

Geographically widespread honeybee‐gut symbiont subgroups show locally distinct antibiotic‐resistant patterns

How long‐term antibiotic treatment affects host bacterial associations is still largely unknown. The honeybee‐gut microbiota has a simple composition, so we used this gut community to investigate how long‐term antibiotic treatment affects host‐associated microbiota. We investigated the phylogenetic relatedness, genomic content (GC percentage, genome size, number of genes and CRISPR) and antibiotic‐resistant genes (ARG) for strains from two abundant members of the honeybee core gut microbiota (Gilliamella apicola and Snodgrassella alvi). Domesticated honeybees are subjected to geographically different management policies, so we used two research apiaries, representing different antibiotic treatment regimens in their apiculture: low antibiotic usage (Norway) and high antibiotic usage (Arizona, USA). We applied whole‐genome shotgun sequencing on 48 G. apicola and 22 S. alvi. We identified three predominating subgroups of G. apicola in honeybees from both Norway and Arizona. For G. apicola, genetic content substantially varied between subgroups and distance similarity calculations showed similarity discrepancy between subgroups. Functional differences between subgroups, such as pectin‐degrading enzymes (G. apicola), were also identified. In addition, we identified horizontal gene transfer (HGT) of transposon (Tn10)‐associated tetracycline resistance (Tet B) across the G. apicola subgroups in the Arizonan honeybees, using interspace polymorphisms in the Tet B determinant. Our results support that honeybee‐gut symbiont subgroups can resist long‐term antibiotic treatment and maintain functionality through acquisition of geographically distinct antibiotic‐resistant genes by HGT.     

Secondary contact between Lycaeides idas and L. melissa in the Rocky Mountains: extensive admixture and a patchy hybrid zone

Studies of hybridization have increased our understanding of the nature of species boundaries, the process of speciation, and the effects of hybridization on the evolution of populations and species. In the present study we use genetic and morphological data to determine the outcome and consequences of secondary contact and hybridization between the butterfly species Lycaeides idas and L. melissa in the Rocky Mountains. Admixture proportions estimated from structure and geographical cline analysis indicate L. idas and L. melissa have hybridized extensively in the Rocky Mountains and that reproductive isolation was insufficient to prevent introgression for much of the genome. Geographical patterns of admixture suggest that hybridization between L. idas and L. melissa has led to the formation of a hybrid zone. The hybrid zone is relatively wide, given estimates of dispersal for Lycaeides butterflies, and does not show strong evidence of cline concordance among characters. We believe the structure of the Lycaeides hybrid zone might be best explained by the patchy distribution of Lycaeides, local extinction and colonization of habitat patches, environmental variation and weak overall selection against hybrids. We found no evidence that hybridization in the Rocky Mountains has resulted in the formation of independent hybrid species, in contrast to the outcome of hybridization between L. idas and L. melissa in the Sierra Nevada. Finally, our results suggest that differences in male morphology between L. idas and L. melissa might contribute to isolation, or perhaps even that selection has favoured the spread of L. melissa male genitalia alleles.     

Interploidy hybridization in sympatric zones: the formation of Epidendrum fulgens × E. puniceoluteum hybrids (Epidendroideae,Orchidaceae)

Interspecific hybridization is a primary cause of extensive morphological and chromosomal variation and plays an important role in plant species diversification. However, the role of interploidal hybridization in the formation of hybrid swarms is less clear. Epidendrum encompasses wide variation in chromosome number and lacks strong premating barriers, making the genus a good model for clarifying the role of chromosomes in postzygotic barriers in interploidal hybrids. In this sense, hybrids from the interploidal sympatric zone between E. fulgens (2n = 2x = 24) and E. puniceoluteum (2n = 4x = 56) were analyzed using cytogenetic techniques to elucidate the formation and establishment of interploidal hybrids. Hybrids were not a uniform group: two chromosome numbers were observed, with the variation being a consequence of severe hybrid meiotic abnormalities and backcrossing with E. puniceoluteum. The hybrids were triploids (2n = 3x = 38 and 40) and despite the occurrence of enormous meiotic problems associated with triploidy, the hybrids were able to backcross, producing successful hybrid individuals with broad ecological distributions. In spite of the nonpolyploidization of the hybrid, its formation is a long‐term evolutionary process rather than a product of a recent disturbance, and considering other sympatric zones in Epidendrum, these events could be recurrent.