W. Brian Simison

A phylogenomic analysis of turtles

Molecular analyses of turtle relationships have overturned prevailing morphological hypotheses and prompted the development of a new taxonomy. Here we provide the first genome-scale analysis of turtle phylogeny. We sequenced 2381 ultraconserved element (UCE) loci representing a total of 1,718,154bp of aligned sequence. Our sampling includes 32 turtle taxa representing all 14 recognized turtle families and an additional six outgroups. Maximum likelihood, Bayesian, and species tree methods produce a single resolved phylogeny. This robust phylogeny shows that proposed phylogenetic names correspond to well-supported clades, and this topology is more consistent with the temporal appearance of clades and paleobiogeography. Future studies of turtle phylogeny using fossil turtles should use this topology as a scaffold for their morphological phylogenetic analyses.

Genomic signatures of geographic isolation and natural selection in coral reef fishes

The drivers of speciation remain among the most controversial topics in evolutionary biology. Initially, Darwin emphasized natural selection as a primary mechanism of speciation, but the architects of the modern synthesis largely abandoned that view in favour of divergence by geographic isolation. The balance between selection and isolation is still at the forefront of the evolutionary debate, especially for the worlds tropical oceans where biodiversity is high, but isolating barriers are few. Here, we identify the drivers of speciation in Pacific reef fishes of the genus Acanthurus by comparative genome scans of two peripheral populations that split from a large Central‐West Pacific lineage at roughly the same time. Mitochondrial sequences indicate that populations in the Hawaiian Archipelago and the Marquesas Islands became isolated approximately 0.5 Ma. The Hawaiian lineage is morphologically indistinguishable from the widespread Pacific form, but the Marquesan form is recognized as a distinct species that occupies an unusual tropical ecosystem characterized by upwelling, turbidity, temperature fluctuations, algal blooms and little coral cover. An analysis of 3737 SNPs reveals a strong signal of selection at the Marquesas, with 59 loci under disruptive selection including an opsin Rh2 locus. While both the Hawaiian and Marquesan populations indicate signals of drift, the former shows a weak signal of selection that is comparable with populations in the Central‐West Pacific. This contrast between closely related lineages reveals one population diverging due primarily to geographic isolation and genetic drift, and the other achieving taxonomic species status under the influence of selection.

Global divergence of the human follicle mite Demodex folliculorum: Persistent associations between host ancestry and mite lineages

Significance Mites live in human hair follicles and have been implicated in medically important skin disorders, but we know surprisingly little about these residents of our skin. By analyzing the variation segregating among 241 mite sequences isolated from 70 human hosts, we showed that hosts with different regional ancestries harbor distinct lineages of mites and that these associations can persist despite generations spent in a new geographic region. These results suggest that some mite populations are better able to survive and reproduce on hosts from certain geographic regions. Improving our understanding of human follicle mites promises to shed light on human evolution and to provide important contextual information for their role in human health. Microscopic mites of the genus Demodex live within the hair follicles of mammals and are ubiquitous symbionts of humans, but little molecular work has been done to understand their genetic diversity or transmission. Here we sampled mite DNA from 70 human hosts of diverse geographic ancestries and analyzed 241 sequences from the mitochondrial genome of the species Demodex folliculorum. Phylogenetic analyses recovered multiple deep lineages including a globally distributed lineage common among hosts of European ancestry and three lineages that primarily include hosts of Asian, African, and Latin American ancestry. To a great extent, the ancestral geography of hosts predicted the lineages of mites found on them; 27% of the total molecular variance segregated according to the regional ancestries of hosts. We found that D. folliculorum populations are stable on an individual over the course of years and that some Asian and African American hosts maintain specific mite lineages over the course of years or generations outside their geographic region of birth or ancestry. D. folliculorum haplotypes were much more likely to be shared within families and between spouses than between unrelated individuals, indicating that transmission requires close contact. Dating analyses indicated that D. folliculorum origins may predate modern humans. Overall, D. folliculorum evolution reflects ancient human population divergences, is consistent with an out-of-Africa dispersal hypothesis, and presents an excellent model system for further understanding the history of human movement.

Introgression and selection shaped the evolutionary history of sympatric sister-species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.

Dating and biogeographical patterns in the sea slug genus Acanthodoris Gray, 1850 (Mollusca, Gastropoda, Nudibranchia).

Recent studies investigating vicariance and dispersal have been focused on correlating major geological events with instances of taxonomic expansion by incorporating the fossil record with molecular clock analyses. However, this approach becomes problematic for soft-bodied organisms that are poorly represented in the fossil record. Here, we estimate the phylogenetic relationships of the nudibranch genus Acanthodoris Gray, 1850 using three molecular markers (16S, COI, H3), and then test two alternative geologically calibrated molecular clock scenarios in BEAST and their effect on ancestral area reconstruction (AAR) estimates employed in LAGRANGE. The global temperate distribution of Acanthodoris spans multiple geological barriers, including the Bering Strait (∼5.32 Mya) and the Baja Peninsula (∼5.5 Mya), both of which are used in our dating estimates. The expansion of the Atlantic Ocean (∼95-105 Mya) is also used to calibrate the relationship between A. falklandica Eliot, 1905 and A. planca Fahey and Valdés, 2005, which are distributed in southern Chile and South Africa respectively. Phylogenetic analyses recovered strong biogeographical signal and recovered two major clades representing northern and southern hemispheric distributions of Acanthodoris. When all three geological events are applied to the calibration analyses, the age for Acanthodoris is estimated to be mid-Cretaceous. When the expansion of the Atlantic Ocean is excluded from our analyses, however, Acanthodoris is estimated to be much younger, with a divergence time estimate during the Miocene. Regardless of divergence estimates, our AAR suggests that Acanthodoris may have origins in the Atlantic Ocean with the Atlantic acting as a dispersal point to the northeastern Pacific. These results suggest that Acanthodoris exhibits a rare instance of western trans-arctic expansion. This study also shows that northeast Pacific specimens of A. pilosa should be regarded as A. atrogriseata and that A. serpentinotus should be regarded as a synonym of A. pina.

Reassessing the conservation status of the shrew Crocidura thomensis , endemic to São Tomé Island

The shrew Crocidura thomensis is a little-known species endemic to Sao Tome Island. We review its distribution, ecology and conservation status based on nine published and 23 new records. The species has a wide distribution across Sao Tome, preferring rugged forested areas with high rainfall. The location of new records coincides with that of historical records but the proportion of records in plantations has declined, possibly as a result of agricultural intensification, increased use of pesticides, and presence of exotic species. The shrew is restricted to a single island, its extent of occurrence is 2 and its habitat is declining in extent and quality, and thus its categorization as Endangered on the IUCN Red List remains appropriate. It is important to gain a better knowledge of population trends, ecological preferences and sensitivity to potential threats, but the effective protection of Sao Tome Obo Natural Park and surrounding forests is the most important measure to ensure the long-term survival of this mammal.

Isolation and characterization of microsatellite markers for identifying hybridization and genetic pollution associated with red-eared slider turtles (Trachemys scripta elegans)

Fourteen microsatellite markers were isolated and characterized from a microsatellite enriched genomic library obtained for the invasive red-eared slider turtle Trachemys scripta elegans (TSE). All but one locus conformed to Hardy–Weinberg equilibrium and none displayed linkage disequilibrium. The number of alleles ranged from 5 to 17 with observed heterozygosity values ranging from 0.26 to 0.96 for 23 individuals collected from its native range in the Mississippi River drainage. These markers will be useful in characterizing TSE populations and invaded populations of other Trachemys and could aid in the development of conservation management policies for invasive TSE.

Whole-genome assembly of the coral reef Pearlscale Pygmy Angelfish ( Centropyge vrolikii)

The diversity of DNA sequencing methods and algorithms for genome assemblies presents scientists with a bewildering array of choices. Here, we construct and compare eight candidate assemblies combining overlapping shotgun read data, mate-pair and Chicago libraries and four different genome assemblers to produce a high-quality draft genome of the iconic coral reef Pearlscale Pygmy Angelfish, Centropyge vrolikii (family Pomacanthidae). The best candidate assembly combined all four data types and had a scaffold N50 127.5 times higher than the candidate assembly obtained from shotgun data only. Our best candidate assembly had a scaffold N50 of 8.97 Mb, contig N50 of 189,827, and 97.4% complete for BUSCO v2 (Actinopterygii set) and 95.6% complete for CEGMA matches. These contiguity and accuracy scores are higher than those of any other fish assembly released to date that did not apply linkage map information, including those based on more expensive long-read sequencing data. Our analysis of how different data types improve assembly quality will help others choose the most appropriate de novo genome sequencing strategy based on resources and target applications. Furthermore, the draft genome of the Pearlscale Pygmy angelfish will play an important role in future studies of coral reef fish evolution, diversity and conservation.

Sequence-Based Analysis of Thermal Adaptation and Protein Energy Landscapes in an Invasive Blue Mussel (Mytilus galloprovincialis)

Abstract Adaptive responses to thermal stress in poikilotherms plays an important role in determining competitive ability and species distributions. Amino acid substitutions that affect protein stability and modify the thermal optima of orthologous proteins may be particularly important in this context. Here, we examine a set of 2,770 protein-coding genes to determine if proteins in a highly invasive heat tolerant blue mussel (Mytilus galloprovincialis) contain signals of adaptive increases in protein stability relative to orthologs in a more cold tolerant M. trossulus. Such thermal adaptations might help to explain, mechanistically, the success with which the invasive marine mussel M. galloprovincialis has displaced native species in contact zones in the eastern (California) and western (Japan) Pacific. We tested for stabilizing amino acid substitutions in warm tolerant M. galloprovincialis relative to cold tolerant M. trossulus with a generalized linear model that compares in silico estimates of recent changes in protein stability among closely related congeners. Fixed substitutions in M. galloprovincialis were 3,180.0 calories per mol per substitution more stabilizing at genes with both elevated dN/dS ratios and transcriptional responses to heat stress, and 705.8 calories per mol per substitution more stabilizing across all 2,770 loci investigated. Amino acid substitutions concentrated in a small number of genes were more stabilizing in M. galloprovincialis compared with cold tolerant M. trossulus. We also tested for, but did not find, enrichment of a priori GO terms in genes with elevated dN/dS ratios in M. galloprovincialis. This might indicate that selection for thermodynamic stability is generic across all lineages, and suggests that the high change in estimated protein stability that we observed in M. galloprovincialis is driven by selection for extra stabilizing substitutions, rather than by higher incidence of selection in a greater number of genes in this lineage. Nonetheless, our finding of more stabilizing amino acid changes in the warm adapted lineage is important because it suggests that adaption for thermal stability has contributed to M. galloprovincialis’ superior tolerance to heat stress, and that pairing tests for positive selection and tests for transcriptional response to heat stress can identify candidates of protein stability adaptation.