ExRec: a python pipeline for generating recombination-filtered multi-locus datasets.

Summary: ExRec (Exclusion of Recombined DNA) is a dependency-free Python pipeline that implements the four-gamete test to automatically filter out recombined DNA blocks from thousands of DNA sequence loci. This procedure helps all loci better meet the "no intralocus recombination" assumption common to many coalescent-based analyses in population genomic, phylogeographic, and shallow-scale phylogenomic studies. The user-friendly pipeline contains five standalone applications-four file conversion scripts and one main script that performs the recombination filtering procedures. The pipeline outputs recombination-filtered data in a variety of common formats and a tab-delimited table that displays descriptive statistics for all loci and the analysis results. A novel feature of this software is that the user can select whether to output the longest nonrecombined sequence blocks from recombined loci (current best practice) or randomly select nonrecombined blocks from loci (a newer approach). We tested ExRec with six published phylogenomic datasets that ranged in size from 27 to 2237 loci and came in a variety of input file formats. In all trials the data could be easily analyzed in only seconds for the smaller datasets and <30 min for the largest using a simple laptop computer. Availability and implementation: ExRec was written in Python 3 under the MIT license. The program applications, user manual (including step-by-step tutorials), and sample data are freely available at https://github.com/Sammccarthypotter/ExRec.

Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles.

Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.

Evolutionary relationships among the snakelike pygopodid lizards: a review of phylogenetic studies of an enigmatic Australian adaptive radiation.

Here, I review phylogenetic studies of the lizard family Pygopodidae, a group of 47 extant species that diversified in Australia and New Guinea. The goal of this study was to examine published phylogenetic and phylogenomic hypotheses on pygopodids to identify the strengths and weaknesses in our understanding of their phylogeny. Many parts of the pygopodid family tree are well established by multiple independent tree inferences including: (1) all multispecies genera (i.e., Aprasia, Delma, Lialis, Pletholax, and Pygopus) are monophyletic groups; (2) the root of the pygopodid tree is located along the branch leading to the Delma clade, thus showing that Delma is the sister group to all other pygopodid genera; (3) the Aprasia repens group, Delma tincta group, and several other groups of closely related species are demonstrated to be monophyletic entities; and (4) the monotypic Paradelma orientalis is the sister lineage to the Pygopus clade. Based on accumulated phylogenetic evidence, two taxonomic recommendations are given: Paradelma merits generic status rather than being subsumed into Pygopus as some earlier studies had suggested, and the monotypic Aclys concinna should be recognized as a member of Delma (following current practice) until future studies clarify its placement inside or outside the Delma clade. One chronic problem with phylogenetic studies of pygopodids, which has limited the explanatory power of many tree hypotheses, concerns the undersampling of known species. Although the continual addition of newly described species, especially over the past two decades, has been a major reason for these taxon sampling gaps, deficits in species sampling for ingroups and/or outgroups in several studies of pygopodid species complexes has confounded the testing of some ingroup monophyly hypotheses. Ancient hybridization between non-sister lineages may also be confounding attempts to recover the relationships among pygopodids using molecular data. Indeed, such a phenomenon can explain at least five cases of mito-nuclear discordance and conflicts among trees based on nuclear DNA datasets. Another problem has been the lack of consensus on the relationships among most pygopodid genera, an issue that may stem from rapid diversification of these lineages early in the group's history. Despite current weaknesses in our understanding of pygopodid phylogeny, enough evidence exists to clarify many major and minor structural parts of their family tree. Accordingly, a composite tree for the Pygopodidae was able to be synthesized. This novel tree hypothesis contains all recognized pygopodid species and reveals that about half of the clades are corroborated by multiple independent tree hypotheses, while the remaining clades have less empirical support.

Barcoding the Neotropical freshwater fish fauna using a new pair of universal COI primers with a discussion of primer dimers and M13 primer tails.

Designing primers for DNA barcoding is a significant challenge for the rich Neotropical fish fauna, which is comprised of ∼6000 species. Previously, researchers required multiple pairs of PCR primers or primer cocktails to obtain standard COI (i.e., mitochondrial cytochrome c oxidase subunit I) barcode sequences from assemblages of freshwater fish in this region. To simplify DNA barcoding and metabarcoding studies of Neotropical freshwater fish, we present a new pair of COI primers, which have yielded high quality barcodes across six teleost orders-Characiformes, Cichliformes, Cyprinodontiformes, Gymnotiformes, Siluriformes, and Synbranchiformes-native to South America. Following previous fish barcoding studies, we also tailed our primers with M13 forward and reverse primers to facilitate the DNA sequencing process. Although this practice generates primer dimers, we obtained complete and high quality COI barcode sequences for all samples. We discuss the problem of primer dimers and suggest strategies for neutralizing their influence on data quality.

The people living with HIV stigma survey UK 2015: HIV-related sexual rejection and other experiences of stigma and discrimination among gay and heterosexual men.

We aim to understand the difference in stigma and discrimination, in particular sexual rejection, experienced between gay and heterosexual men living with HIV in the UK. The People Living with HIV StigmaSurvey UK 2015 recruited a convenience sample of persons with HIV through over 120 cross sector community organisations and 46 HIV clinics to complete an online survey. 1162 men completed the survey, 969 (83%) gay men and 193 (17%) heterosexual men, 92% were on antiretroviral therapy. Compared to heterosexual men, gay men were significantly more likely to report worrying about workplace treatment in relation to their HIV (21% vs. 11%), worrying about HIV-related sexual rejection (42% vs 21%), avoiding sex because of their HIV status (37% vs. 23%), and experiencing HIV-related sexual rejection (27% vs. 9%) in the past 12 months. In a multivariate logistic regression controlling for other sociodemographic factors, being gay was a predictor of reporting HIV-related sexual rejection in the past 12 months (aOR 2.17, CI 1.16, 4.02). Both gay and heterosexual men living with HIV experienced stigma and discrimination in the past 12 months, and this was higher for gay men in terms of HIV-related sexual rejection. Due to the high proportion of men reporting sexual rejection, greater awareness and education of the low risk of transmission of HIV among people on effective treatment is needed to reduce stigma and sexual prejudice towards people living with HIV.

Experiences of stigma and discrimination in social and healthcare settings among trans people living with HIV in the UK.

The People Living with HIV StigmaSurvey UK 2015 was a community led national survey investigating experiences of people living with HIV in the UK in the past 12 months. Participants aged 18 and over were recruited through over 120 cross-sector community organisations and 46 HIV clinics to complete an anonymous online survey. Trans is an umbrella term which refers to individuals whose current gender identity is different to the gender they were assigned at birth. Trans participants self-identified via gender identity and gender at birth questions. Descriptive analyses of reported experiences in social and health care settings were conducted and multivariate logistic regression analyses were used to identify sociodemographic predictors of reporting being treated differently to non-HIV patients, and being delayed or refused healthcare treatment in the past 12 months. 31 out of 1576 participants (2%) identified as trans (19 trans women, 5 trans men, 2 gender queer/non-binary, 5 other). High levels of social stigma were reported for all participants, with trans participants significantly more likely to report worrying about verbal harassment (39% vs. 23%), and exclusion from family gatherings (23% vs. 9%) in the last 12 months, compared to cisgender participants. Furthermore, 10% of trans participants reported physical assault in the last 12 months, compared to 4% of cisgender participants. Identifying as trans was a predictor of reporting being treated differently to non-HIV patients (48% vs. 30%; aOR 2.61, CI 1.06, 6.42) and being delayed or refused healthcare (41% vs. 16%; aOR 4.58, CI 1.83, 11.44). Trans people living with HIV in the UK experience high levels of stigma and discrimination, including within healthcare settings, which is likely to impact upon health outcomes. Trans-specific education and awareness within healthcare settings could help to improve service provision for this demographic.

On the independent gene trees assumption in phylogenomic studies.

Multilocus coalescent methods for inferring species trees or historical demographic parameters typically require the assumption that gene trees for sampled SNPs or DNA sequence loci are conditionally independent given their species tree. In practice, researchers have used different criteria to delimit "independent loci." One criterion identifies sampled loci as being independent of each other if they undergo Mendelian independent assortment (IA criterion). O'Neill et al. (2013, Molecular Ecology, 22, 111-129) used this approach in their phylogeographic study of North American tiger salamander species complex. In two other studies, researchers developed a pair of related methods that employ an independent genealogies criterion (IG criterion), which considers the effects of population-level recombination on correlations between the gene trees of intrachromosomal loci. Here, I explain these three methods, illustrate their use with example data, and evaluate their efficacies. I show that the IA approach is more conservative, is simpler to use and requires fewer assumptions than the IG approaches. However, IG approaches can identify much larger numbers of independent loci than the IA method, which, in turn, allows researchers to obtain more precise and accurate estimates of species trees and historical demographic parameters. A disadvantage of the IG methods is that they require an estimate of the population recombination rate. Despite their drawbacks, IA and IG approaches provide molecular ecologists with promising a priori methods for selecting SNPs or DNA sequence loci that likely meet the independence assumption in coalescent-based phylogenomic studies.

Stereodynamics and edge-to-face CH-π aromatic interactions in imino compounds containing heterocyclic rings.

By comparison with close contact interactions between benzene rings there is a paucity of experimental data available for attractive interactions involving aromatic heterocyclic rings, especially for small molecules in solution. Herein we describe aromatic heterocyclic and carbocyclic edge-to face interactions and conformational stereodynamics of N-1,2-diphenylethyl imines bearing a phenyl group and either a 2-pyridyl, 3-pyridyl, 2-thiophene or 2-furanyl moiety on the imino carbon. X-ray crystal structures have been determined for two compounds. Slow rotation about the phenyl-imino bond in the E-isomers and around the heterocycle-imino bond in the Z-isomers of the pyridyl compounds was observed at low temperatures by NMR. Abnormally large shielding of one ortho hydrogen indicates that both the imino phenyl and heterocycle rings can engage in an edge-to-face interaction with the N-terminal phenyl moiety in the appropriate isomer. Some rotational barriers around the phenyl-imino and heterocycle-imino bonds were measured.

In silico phylogenomics using complete genomes: a case study on the evolution of hominoids.

The increasing availability of complete genome data is facilitating the acquisition of phylogenomic data sets, but the process of obtaining orthologous sequences from other genomes and assembling multiple sequence alignments remains piecemeal and arduous. We designed software that performs these tasks and outputs anonymous loci (AL) or anchored enrichment/ultraconserved element loci (AE/UCE) data sets in ready-to-analyze formats. We demonstrate our program by applying it to the hominoids. Starting with human, chimpanzee, gorilla, and orangutan genomes, our software generated an exhaustive data set of 292 ALs (∼1 kb each) in ∼3 h. Not only did analyses of our AL data set validate the program by yielding a portrait of hominoid evolution in agreement with previous studies, but the accuracy and precision of our estimated ancestral effective population sizes and speciation times represent improvements. We also used our program with a published set of 512 vertebrate-wide AE "probe" sequences to generate data sets consisting of 171 and 242 independent loci (∼1 kb each) in 11 and 13 min, respectively. The former data set consisted of flanking sequences 500 bp from adjacent AEs, while the latter contained sequences bordering AEs. Although our AE data sets produced the expected hominoid species tree, coalescent-based estimates of ancestral population sizes and speciation times based on these data were considerably lower than estimates from our AL data set and previous studies. Accordingly, we suggest that loci subjected to direct or indirect selection may not be appropriate for coalescent-based methods. Complete in silico approaches, combined with the burgeoning genome databases, will accelerate the pace of phylogenomics.

Continuity of cannabis use and violent offending over the life course.

BACKGROUND: Although the association between cannabis use and violence has been reported in the literature, the precise nature of this relationship, especially the directionality of the association, is unclear. METHOD: Young males from the Cambridge Study of Delinquent Development (n = 411) were followed up between the ages of 8 and 56 years to prospectively investigate the association between cannabis use and violence. A multi-wave (eight assessments, T1-T8) follow-up design was employed that allowed temporal sequencing of the variables of interest and the analysis of violent outcome measures obtained from two sources: (i) criminal records (violent conviction); and (ii) self-reports. A combination of analytic approaches allowing inferences as to the directionality of associations was employed, including multivariate logistic regression analysis, fixed-effects analysis and cross-lagged modelling. RESULTS: Multivariable logistic regression revealed that compared with never-users, continued exposure to cannabis (use at age 18, 32 and 48 years) was associated with a higher risk of subsequent violent behaviour, as indexed by convictions [odds ratio (OR) 7.1, 95% confidence interval (CI) 2.19-23.59] or self-reports (OR 8.9, 95% CI 2.37-46.21). This effect persisted after controlling for other putative risk factors for violence. In predicting violence, fixed-effects analysis and cross-lagged modelling further indicated that this effect could not be explained by other unobserved time-invariant factors. Furthermore, these analyses uncovered a bi-directional relationship between cannabis use and violence. CONCLUSIONS: Together, these results provide strong indication that cannabis use predicts subsequent violent offending, suggesting a possible causal effect, and provide empirical evidence that may have implications for public policy.

The complete mitochondrial genome of the ruby-topaz hummingbird Chrysolampis mosquitus through Illumina sequencing.

The complete mitochondrial genome of the Ruby-Topaz Hummingbird, Chrysolampis mosquitus, was determined using 1/11 of an Illumina Hi-seq lane ran with a Nextera kit. We assembled the mitogenome in a two-step approach using both (i) de novo (SOAPdenovo-Trans) and (ii) reference-based (MITObim) genome assembly software. A circular molecule containing 17,767 bp was assembled. As expected for most vertebrates, the C. mosquitus mitogenome contained 13 protein-coding genes, 22 transfer RNA, 2 ribosomal RNA genes, and 1 non-coding control region. We assembled the whole mitogenome using 0.45% of the total amount of reads and produced a high-coverage mitochondrial genome (>1000×). We deposited the assembled mitogenome into GenBank (accession number KJ619585).

Complete mitochondrial genome of the versicoloured emerald hummingbird Amazilia versicolor, a polymorphic species.

The genome of the versicoloured emerald hummingbird (Amazilia versicolor) was partially sequenced in one-sixth of an Illumina HiSeq lane. The mitochondrial genome was assembled using MIRA and MITObim software, yielding a circular molecule of 16,861 bp in length and deposited in GenBank under the accession number KF624601. The mitogenome contained 13 protein-coding genes, 22 transfer tRNAs, 2 ribosomal RNAs and 1 non-coding control region. The molecule was assembled using 21,927 sequencing reads of 100 bp each, resulting in ∼130 × coverage of uniformly distributed reads along the genome. This is the forth mitochondrial genome described for this highly diverse family of birds and may benefit further phylogenetic, phylogeographic, population genetic and species delimitation studies of hummingbirds.

DNA barcoding reveals species level divergence between populations of the microhylid frog genus Arcovomer (Anura: Microhylidae) in the Atlantic Rainforest of southeastern Brazil.

The microhylid frogs belonging to the genus Arcovomer have been reported from lowland Atlantic Rainforest in the Brazilian states of Espírito Santo, Rio de Janeiro, and São Paulo. Here, we use DNA barcoding to assess levels of genetic divergence between apparently isolated populations in Espírito Santo and Rio de Janeiro. Our mtDNA data consisting of cytochrome oxidase subunit I (COI) nucleotide sequences reveals 13.2% uncorrected and 30.4% TIM2 + I + Γ corrected genetic divergences between these two populations. This level of divergence exceeds the suggested 10% uncorrected divergence threshold for elevating amphibian populations to candidate species using this marker, which implies that the Espírito Santo population is a species distinct from Arcovomer passarellii. Calibration of our model-corrected sequence divergence estimates suggests that the time of population divergence falls between 12 and 29 million years ago.

Desert tortoises (Gopherus agassizii) are selective herbivores that track the flowering phenology of their preferred food plants.

Previous studies of desert tortoise foraging ecology in the western Mojave Desert suggest that these animals are selective herbivores, which alter their diet according to the temporal availability of preferred food plants. These studies, however, did not estimate availability of potential food plants by taking into account the spatial and temporal variability in ephemeral plant abundance that occurs within the spring season. In this study, we observed 18 free-ranging adult tortoises take 35,388 bites during the spring foraging season. We also estimated the relative abundance of potential food plants by stratifying our sampling across different phenological periods of the 3-month long spring season and by different habitats and microhabitats. This methodology allowed us to conduct statistical tests comparing tortoise diet against plant abundance. Our results show that tortoises choose food plants non-randomly throughout the foraging season, a finding that corroborates the hypothesis that desert tortoises rely on key plants during different phenological periods of spring. Moreover, tortoises only consumed plants in a succulent state until the last few weeks of spring, at which time most annuals and herbaceous perennials had dried and most tortoises had ceased foraging. Many species of food plants--including several frequently eaten species--were not detected in our plant surveys, yet tortoises located these rare plants in their home ranges. Over 50% of bites consumed were in the group of undetected species. Interestingly, tortoises focused heavily on several leguminous species, which could be nutritious foods owing to their presumably high nitrogen contents. We suggest that herbaceous perennials, which were rare on our study area but represented ~30% of tortoise diet, may be important in sustaining tortoise populations during droughts when native annuals are absent. These findings highlight the vulnerability of desert tortoises to climate change if such changes alter the availability of their preferred food plants.

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.

An evaluation of substituent effects on aromatic edge-to-face interactions and CF-π versus CH-π interactions using an imino torsion balance model.

A selection of imines derived from phenyl t-butyl ketones and substituted 2-phenylethylamines or phenylalanine exhibit slow rotation around the aryl­imino bond at ambient temperature, resulting in a large non-equivalence of the ortho hydrogens in the 1H NMR spectra. This facilitates assessment of aryl substituent effects on the face tilted-T CH­π interaction between a phenyl ring (A) on the imino carbon proximate to the terminal phenyl ring (B). Analysis of the marked temperature dependence of the chemical shift of the interacting ortho hydrogen affords estimates of the opposing enthalpic and entropic factors involved in the rapid equilibrium between the closed edge-to-face conformation and alternative open conformations devoid of a CH­π interaction while in solution. Above ca. 80 °C the entropy term (TΔS) cancels out the enthalpy (ΔH) favouring the closed conformation and open conformations are preferred. Accordingly, commonly reported binding free energies may not be a good measure of the energetic strength of intramolecular aromatic interactions. Investigation of an ortho fluoro substituted compound indicates that a CF­π interaction is at least 1.0 kcal mol−1 weaker in enthalpy than the CH­π interaction. Several X-ray crystal structures depicting an intramolecular edge-to-face interaction are presented.

Defining evolutionary boundaries across parapatric ecomorphs of Black Salamanders (Aneides flavipunctatus) with conservation implications.

The accurate delimitation of evolutionary population units represents an important component in phylogeographic and conservation genetic studies. Here, we used a combined population assignment and historical demographic approach to study a complex of ecomorphologically distinctive populations of Black Salamanders (Aneides flavipunctatus) that are parapatrically distributed and meet at a three-way contact zone in north-western California. We used mitochondrial tree-based and multilocus clustering methods to evaluate a priori two- (Northern and Southern) and three (Northern, Coast and Inland) population hypotheses derived from previous studies. Mitochondrial results were consistent with the two- and three-population hypotheses, while the nDNA clustering results supported only the two-population hypothesis. Historical demographic analyses and mtDNA gene divergence estimates revealed that the Northern and Southern populations split during the Pliocene (2-5 Ma). Subdivision of the Southern population into Coast and Inland populations was estimated to be late Pleistocene (0.24 Ma), although our mtDNA results suggested a Pliocene divergence. Effective gene flow estimates (2N(e)m) suggest that either the two- or three-population hypotheses remain valid. However, our results unexpectedly revealed that the Northern population might instead represent two parapatric populations that separated nearly 4 Ma. These results are surprising because the Pliocene divergence between these ecomorphologically conservative forms is similar or older than for the ecomorphologically divergent Coast and Inland sister populations. We conclude that Black Salamanders in north-western California belong to at least three or four populations or species, and these all meet criteria for being Evolutionary Significant Units or 'ESUs' and therefore warrant conservation consideration.

Evaluation of the Bruker SMART X2S: crystallography for the nonspecialist?

An evaluation of the Bruker SMART X2S for the collection of crystallographic diffraction data, structure solution and refinement is carried out with a variety of materials with different electron densities, presenting some of the successes and challenges of automation in chemical crystallography.