Population genomics of flat-tailed horned lizards (Phrynosoma mcallii) informs conservation and management across a fragmented Colorado Desert landscape.

Phrynosoma mcallii (flat-tailed horned lizards) is a species of conservation concern in the Colorado Desert of the United States and Mexico. We analysed ddRADseq data from 45 lizards to estimate population structure, infer phylogeny, identify migration barriers, map genetic diversity hotspots, and model demography. We identified the Colorado River as the main geographic feature contributing to population structure, with the populations west of this barrier further subdivided by the Salton Sea. Phylogenetic analysis confirms that northwestern populations are nested within southeastern populations. The best-fit demographic model indicates Pleistocene divergence across the Colorado River, with significant bidirectional gene flow, and a severe Holocene population bottleneck. These patterns suggest that management strategies should focus on maintaining genetic diversity on both sides of the Colorado River and the Salton Sea. We recommend additional lands in the United States and Mexico that should be considered for similar conservation goals as those in the Rangewide Management Strategy. We also recommend periodic rangewide genomic sampling to monitor ongoing attrition of diversity, hybridization, and changing structure due to habitat fragmentation, climate change, and other long-term impacts.

Lineage diversification of fringe-toed lizards (Phrynosomatidae: Uma notata complex) in the Colorado Desert: Delimiting species in the presence of gene flow.

Multi-locus nuclear DNA data were used to delimit species of fringe-toed lizards of the Uma notata complex, which are specialized for living in wind-blown sand habitats in the deserts of southwestern North America, and to infer whether Quaternary glacial cycles or Tertiary geological events were important in shaping the historical biogeography of this group. We analyzed ten nuclear loci collected using Sanger sequencing and genome-wide sequence/single-nucleotide polymorphism (SNP) data collected using restriction-associated DNA (RAD) sequencing. A combination of species discovery methods (concatenated phylogenies, parametric and non-parametric clustering algorithms) and species validation approaches (coalescent-based species tree/isolation-with-migration models) were used to delimit species, infer phylogenetic relationships, and to estimate effective population sizes, migration rates, and speciation times. Uma notata, U. inornata, U. cowlesi, and an undescribed species from Mohawk Dunes, Arizona (U. sp.) were supported as distinct in the concatenated analyses and by clustering algorithms, and all operational taxonomic units were decisively supported as distinct species by ranking hierarchical nested speciation models with Bayes factors based on coalescent-based species tree methods. However, significant unidirectional gene flow (2NM>1) from U. cowlesi and U. notata into U. rufopunctata was detected under the isolation-with-migration model. Therefore, we conservatively delimit four species-level lineages within this complex (U. inornata, U. notata, U. cowlesi, and U. sp.), treating U. rufopunctata as a hybrid population (U. notata×cowlesi). Both concatenated and coalescent-based estimates of speciation times support the hypotheses that speciation within the complex occurred during the late Pleistocene, and that the geological evolution of the Colorado River delta during this period was an important process shaping the observed phylogeographic patterns.

High resolution mapping of the tumor microenvironment using integrated single-cell, spatial and in situ analysis.

Single-cell and spatial technologies that profile gene expression across a whole tissue are revolutionizing the resolution of molecular states in clinical samples. Current commercially available technologies provide whole transcriptome single-cell, whole transcriptome spatial, or targeted in situ gene expression analysis. Here, we combine these technologies to explore tissue heterogeneity in large, FFPE human breast cancer sections. This integrative approach allowed us to explore molecular differences that exist between distinct tumor regions and to identify biomarkers involved in the progression towards invasive carcinoma. Further, we study cell neighborhoods and identify rare boundary cells that sit at the critical myoepithelial border confining the spread of malignant cells. Here, we demonstrate that each technology alone provides information about molecular signatures relevant to understanding cancer heterogeneity; however, it is the integration of these technologies that leads to deeper insights, ushering in discoveries that will progress oncology research and the development of diagnostics and therapeutics.

A new cryptic species of fringe-toed lizards from southwestern Arizona with a revised taxonomy of the Uma notata species complex (Squamata: Phrynosomatidae).

Fringe-toed lizards (Uma) are among the most specialized lizards in North America, adapted to insular windblown sand habitats in the hyper-arid southwestern deserts, with allopatric distributions, subtle morphological variation, and an unstable taxonomic history. We analyzed a morphological dataset of 40 characters for 65 specimens and a molecular dataset of 2,286 bases from three mitochondrial loci for 92 individuals and interpreted these data alongside published analyses of multi-locus genetic data with the goal of revising the taxonomy of the Uma notata (Baird 1858) species complex. We confirmed that fringe-toed lizards from the Mohawk Dunes in southwestern Arizona (U. sp.) constitute a cryptic species sister to the rest of the complex that can be diagnosed with DNA barcoding and geography, so we describe and name this species Uma thurmanae sp. nov. We also confirmed the evolutionary distinctiveness of U. inornata (Cope 1895), an endangered species endemic to Coachella Valley in southern California. We designate a lectotype for the taxon U. "rufopunctata", but we put its name in quotation marks to reflect its uncertain taxonomic status with respect to its neighboring species U. cowlesi and U. notata.

Zoogeography of the San Andreas Fault system: Great Pacific Fracture Zones correspond with spatially concordant phylogeographic boundaries in western North America.

The purpose of this article is to provide an ultimate tectonic explanation for several well-studied zoogeographic boundaries along the west coast of North America, specifically, along the boundary of the North American and Pacific plates (the San Andreas Fault system). By reviewing 177 references from the plate tectonics and zoogeography literature, I demonstrate that four Great Pacific Fracture Zones (GPFZs) in the Pacific plate correspond with distributional limits and spatially concordant phylogeographic breaks for a wide variety of marine and terrestrial animals, including invertebrates, fish, amphibians, reptiles, birds, and mammals. These boundaries are: (1) Cape Mendocino and the North Coast Divide, (2) Point Conception and the Transverse Ranges, (3) Punta Eugenia and the Vizcaíno Desert, and (4) Cabo Corrientes and the Sierra Transvolcanica. However, discussion of the GPFZs is mostly absent from the zoogeography and phylogeography literature likely due to a disconnect between biologists and geologists. I argue that the four zoogeographic boundaries reviewed here ultimately originated via the same geological process (triple junction evolution). Finally, I suggest how a comparative phylogeographic approach can be used to test the hypothesis presented here.

Phylogenomics of phrynosomatid lizards: conflicting signals from sequence capture versus restriction site associated DNA sequencing.

Sequence capture and restriction site associated DNA sequencing (RADseq) are popular methods for obtaining large numbers of loci for phylogenetic analysis. These methods are typically used to collect data at different evolutionary timescales; sequence capture is primarily used for obtaining conserved loci, whereas RADseq is designed for discovering single nucleotide polymorphisms (SNPs) suitable for population genetic or phylogeographic analyses. Phylogenetic questions that span both "recent" and "deep" timescales could benefit from either type of data, but studies that directly compare the two approaches are lacking. We compared phylogenies estimated from sequence capture and double digest RADseq (ddRADseq) data for North American phrynosomatid lizards, a species-rich and diverse group containing nine genera that began diversifying approximately 55 Ma. Sequence capture resulted in 584 loci that provided a consistent and strong phylogeny using concatenation and species tree inference. However, the phylogeny estimated from the ddRADseq data was sensitive to the bioinformatics steps used for determining homology, detecting paralogs, and filtering missing data. The topological conflicts among the SNP trees were not restricted to any particular timescale, but instead were associated with short internal branches. Species tree analysis of the largest SNP assembly, which also included the most missing data, supported a topology that matched the sequence capture tree. This preferred phylogeny provides strong support for the paraphyly of the earless lizard genera Holbrookia and Cophosaurus, suggesting that the earless morphology either evolved twice or evolved once and was subsequently lost in Callisaurus.

Speciation, population structure, and demographic history of the Mojave Fringe-toed Lizard (Uma scoparia), a species of conservation concern.

The North American deserts were impacted by both Neogene plate tectonics and Quaternary climatic fluctuations, yet it remains unclear how these events influenced speciation in this region. We tested published hypotheses regarding the timing and mode of speciation, population structure, and demographic history of the Mojave Fringe-toed Lizard (Uma scoparia), a sand dune specialist endemic to the Mojave Desert of California and Arizona. We sampled 109 individual lizards representing 22 insular dune localities, obtained DNA sequences for 14 nuclear loci, and found that U. scoparia has low genetic diversity relative to the U. notata species complex, comparable to that of chimpanzees and southern elephant seals. Analyses of genotypes using Bayesian clustering algorithms did not identify discrete populations within U. scoparia. Using isolation-with-migration (IM) models and a novel coalescent-based hypothesis testing approach, we estimated that U. scoparia diverged from U. notata in the Pleistocene epoch. The likelihood ratio test and the Akaike Information Criterion consistently rejected nested speciation models that included parameters for migration and population growth of U. scoparia. We reject the Neogene vicariance hypothesis for the speciation of U. scoparia and define this species as a single evolutionarily significant unit for conservation purposes.