Understanding Species Boundaries that Arise from Complex Histories: Gene Flow Across the Speciation Continuum in the Spotted Whiptail Lizards.

Gene flow between diverging lineages challenges the resolution of species boundaries and the understanding of evolutionary history in recent radiations. Here, we integrate phylogenetic and coalescent tools to resolve reticulate patterns of diversification and use a perspective focused on evolutionary mechanisms to distinguish interspecific and intraspecific taxonomic variation. We use this approach to resolve the systematics for one of the most intensively studied but difficult to understand groups of reptiles: the spotted whiptail lizards of the genus Aspidoscelis (A. gularis complex). Whiptails contain the largest number of unisexual species known within any vertebrate group and the spotted whiptail complex has played a key role in the generation of this diversity through hybrid speciation. Understanding lineage boundaries and the evolutionary history of divergence and reticulation within this group is therefore key to understanding the generation of unisexual diversity in whiptails. Despite this importance, long-standing confusion about their systematics has impeded understanding of which gonochoristic species have contributed to the formation of unisexual lineages. Using reduced representation genomic data, we resolve patterns of divergence and gene flow within the spotted whiptails and clarify patterns of hybrid speciation. We find evidence that biogeographically structured ecological and environmental variation has been important in morphological and genetic diversification, as well as the maintenance of species boundaries in this system. Our study elucidates how gene flow among lineages and the continuous nature of speciation can bias the practice of species delimitation and lead taxonomists operating under different frameworks to different conclusions (here we propose that a two species arrangement best reflects our current understanding). In doing so, this study provides conceptual and methodological insights into approaches to resolving diversification patterns and species boundaries in rapid radiations with complex histories, as well as long-standing taxonomic challenges in the field of systematic biology.

Investigating the impacts of artificial light via blackouts.

Natural experiments provide remarkable opportunities to test the large-scale effects of human activities. Widespread energy blackouts offer such an 'experiment' to test the impacts of artificial light at night (ALAN) on wildlife. We use the situation in South Africa, where regular scheduled blackouts are being implemented, to highlight this opportunity.

A primal-dual data-driven method for computational optical imaging with a photonic lantern.

Optical fibers aim to image in vivo biological processes. In this context, high spatial resolution and stability to fiber movements are key to enable decision-making processes (e.g. for microendoscopy). Recently, a single-pixel imaging technique based on a multicore fiber photonic lantern has been designed, named computational optical imaging using a lantern (COIL). A proximal algorithm based on a sparsity prior, dubbed SARA-COIL, has been further proposed to solve the associated inverse problem, to enable image reconstructions for high resolution COIL microendoscopy. In this work, we develop a data-driven approach for COIL. We replace the sparsity prior in the proximal algorithm by a learned denoiser, leading to a plug-and-play (PnP) algorithm. The resulting PnP method, based on a proximal primal-dual algorithm, enables to solve the Morozov formulation of the inverse problem. We use recent results in learning theory to train a network with desirable Lipschitz properties, and we show that the resulting primal-dual PnP algorithm converges to a solution to a monotone inclusion problem. Our simulations highlight that the proposed data-driven approach improves the reconstruction quality over variational SARA-COIL method on both simulated and real data.

The Expected Behaviors of Posterior Predictive Tests and Their Unexpected Interpretation.

Poor fit between models of sequence or trait evolution and empirical data is known to cause biases and lead to spurious conclusions about evolutionary patterns and processes. Bayesian posterior prediction is a flexible and intuitive approach for detecting such cases of poor fit. However, the expected behavior of posterior predictive tests has never been characterized for evolutionary models, which is critical for their proper interpretation. Here, we show that the expected distribution of posterior predictive P-values is generally not uniform, in contrast to frequentist P-values used for hypothesis testing, and extreme posterior predictive P-values often provide more evidence of poor fit than typically appreciated. Posterior prediction assesses model adequacy under highly favorable circumstances, because the model is fitted to the data, which leads to expected distributions that are often concentrated around intermediate values. Nonuniform expected distributions of P-values do not pose a problem for the application of these tests, however, and posterior predictive P-values can be interpreted as the posterior probability that the fitted model would predict a dataset with a test statistic value as extreme as the value calculated from the observed data.

Urbanization correlates with the prevalence and richness of blood parasites in Eurasian Blackbirds (Turdus merula).

Urbanization is increasing worldwide, producing severe environmental impacts. Biodiversity is affected by the expansion of cities, with many species being unable to cope with the different human-induced stressors present in these landscapes. However, this knowledge is mainly based on research from taxa such as plants or vertebrates, while other organisms like protozoa have been less studied in this context. The impact of urbanization on the transmission of vector-borne pathogens in wildlife is still unclear despite its relevance for animal and human health. Here, we investigated whether cities are associated with changes in the prevalence and richness of lineages of three vector-borne protozoans (Plasmodium, Haemoproteus and Leucocytozoon) in Eurasian blackbirds (Turdus merula) from multiple urban and forest areas in Europe. Our results show important species-specific differences between these two habitat types. We found a significant lower prevalence of Leucocytozoon in urban birds compared to forest birds, but no differences for Plasmodium and Haemoproteus. Furthermore, the richness of parasite lineages in European cities was higher for Plasmodium but lower for Leucocytozoon than in forests. We also found one Plasmodium lineage exclusively from cities while another of Leucocytozoon was only found in forests suggesting a certain level of habitat specialization for these protozoan vectors. Overall, our findings show that cities provide contrasting opportunities for the transmission of different vector-borne pathogens and generate new scenarios for the interactions between hosts, vectors and parasites.

Fiber-connectorized ultrafast-laser-inscribed K-band integrated optics beam combiner for the CHARA telescope array.

A fiber-connectorized K-band integrated-optics two-telescope beam combiner was developed for long-baseline interferometry at the CHARA telescope array utilizing the ultrafast laser inscription (ULI) technique. Single-mode waveguide insertion losses were measured to be ∼1.1d B over the 2-2.3 µm window. The development of asymmetric directional couplers enabled the construction of a beam combiner that includes a 50:50 coupler for interferometric combination and two ∼75:25 couplers for photometric calibration. The visibility of the bare beam combiner was measured at 87% and then at 82% after fiber-connectorization by optimizing the input polarization. These results indicate that ULI technique can fabricate efficient fiber-connectorized K-band beam combiners for astronomical purposes.

Axi-Stack: a method for manufacturing freeform air-silica optical fibre.

We present a method with potential for fabricating freeform air-silica optical fibre preforms which is free from the stacking constraints associated with conventional stack-and-draw. The method, termed Axi-Stack, is enabled by the precision machining of short cross-sectional preform discs by ultrafast laser assisted etching; a laser-based microfabrication technique which facilitates near arbitrary shaping of the preform structure. Several preform discs are stacked axially and fused together via ultrafast laser welding to construct the preform, which can be drawn to fibre using conventional methods. To illustrate the Axi-Stack process, we detail the fabrication of a 30 cm long solid-core photonic crystal fibre preform with a square lattice of cladding holes and characterise fibre drawn from it.

Testing concordance and conflict in spatial replication of landscape genetics inferences.

The degree to which landscape genetics findings can be extrapolated to different areas of a species range is poorly understood. Here, we used a broadly distributed ectothermic lizard (Sceloporus occidentalis, Western Fence lizard) as a model species to evaluate the full role of topography, climate, vegetation, and roads on dispersal and genetic differentiation. We conducted landscape genetics analyses with a total of 119 individuals in five areas within the Sierra Nevada mountain range. Genetic distances calculated from thousands of ddRAD markers were used to optimize landscape resistance surfaces and infer the effects of landscape and topographic features on genetic connectivity. Across study areas, we found a great deal of consistency in the primary environmental gradients impacting genetic connectivity, along with some site-specific differences, and a range in the proportion of genetic variance explained by environmental factors across study sites. High-elevation colder areas were consistently found to be barriers to gene flow, as were areas of high ruggedness and slope. High temperature seasonality and high precipitation during the winter wet season also presented a substantial barrier to gene flow in a majority of study areas. The effect of other landscape variables on genetic differentiation was more idiosyncratic and depended on specific attributes at each site. Across study areas, canyon valleys were always implicated as facilitators to dispersal and key features linking populations and maintaining genetic connectivity, though the relative importance varied in different areas. We emphasize that spatial data layers are complex and multidimensional, and careful consideration of spatial data correlation structure and robust analytic frameworks will be critical to our continued understanding of spatial genetics processes.

Transcription factors interact with RNA to regulate genes.

Transcription factors (TFs) orchestrate the gene expression programs that define each cell's identity. The canonical TF accomplishes this with two domains, one that binds specific DNA sequences and the other that binds protein coactivators or corepressors. We find that at least half of TFs also bind RNA, doing so through a previously unrecognized domain with sequence and functional features analogous to the arginine-rich motif of the HIV transcriptional activator Tat. RNA binding contributes to TF function by promoting the dynamic association between DNA, RNA, and TF on chromatin. TF-RNA interactions are a conserved feature important for vertebrate development and disrupted in disease. We propose that the ability to bind DNA, RNA, and protein is a general property of many TFs and is fundamental to their gene regulatory function.

Computational Fluorescence Suppression in Shifted Excitation Raman Spectroscopy.

Fiber-based Raman spectroscopy in the context of in vivo biomedical application suffers from the presence of background fluorescence from the surrounding tissue that might mask the crucial but inherently weak Raman signatures. One method that has shown potential for suppressing the background to reveal the Raman spectra is shifted excitation Raman spectroscopy (SER). SER collects multiple emission spectra by shifting the excitation by small amounts and uses these spectra to computationally suppress the fluorescence background based on the principle that Raman spectrum shifts with excitation while fluorescence spectrum does not. We introduce a method that utilizes the spectral characteristics of the Raman and fluorescence spectra to estimate them more effectively, and compare this approach against existing methods on real world datasets.

Temperate Zone Isolation by Climate: An Extension of Janzen's 1967 Hypothesis.

AbstractOne of the most stunning patterns of the distribution of life on Earth is the latitudinal biodiversity gradient. In an influential article, Janzen (1967) predicted that tropical mountains are more effective migration barriers than temperate mountains of the same elevation, because annual temperature variation in the tropics is lower. A great deal of research has demonstrated that the mechanism envisioned by Janzen operates at broad latitudinal scales. However, the extent that the mechanism mediates biodiversity generally, and at smaller scales, is far less understood. We investigated whether climate overlap is associated with genetic similarity between populations within temperate regions using lizards in the Sierra Nevada mountain range of California as a study system. By comparing genetic differentiation between high- and low-elevation populations, we found that in addition to the expected strong pattern of isolation by distance, high climate overlap was negatively associated with genetic differentiation, indicating that population pairs that inhabit climatically similar environments are less genetically differentiated. Moreover, while climate overlap between high- and low-elevation sites is predicted to increase from the equator to temperate regions, we find that in adjacent mountain ranges at the same latitude in temperate regions, climate overlap values can vary widely. This study suggests that in addition to the well-studied main effect of latitude on climate overlap and population differentiation, local climate factors within bioclimatic regions can also influence genetic differentiation between populations and do so by the same general mechanism that operates at larger geographic scales.

Biogeographic inferences across spatial and evolutionary scales.

The field of biogeography unites landscape genetics and phylogeography under a common conceptual framework. Landscape genetics traditionally focuses on recent-time, population-based, spatial genetics processes at small geographical scales, while phylogeography typically investigates deep past, lineage- and species-based processes at large geographical scales. Here, we evaluate the link between landscape genetics and phylogeographical methods using the western fence lizard (Sceloporus occidentalis) as a model species. First, we conducted replicated landscape genetics studies across several geographical scales to investigate how population genetics inferences change depending on the spatial extent of the study area. Then, we carried out a phylogeographical study of population structure at two evolutionary scales informed by inferences derived from landscape genetics results to identify concordance and conflict between these sets of methods. We found significant concordance in landscape genetics processes at all but the largest geographical scale. Phylogeographical results indicate major clades are restricted to distinct river drainages or distinct hydrological regions. At a more recent timescale, we find minor clades are restricted to single river canyons in the majority of cases, while the remainder of river canyons include samples from at most two clades. Overall, the broad-scale pattern implicating stream and river valleys as key features linking populations in the landscape genetics results, and high degree of clade specificity within major topographic subdivisions in the phylogeographical results, is consistent. As landscape genetics and phylogeography share many of the same objectives, synthesizing theory, models and methods between these fields will help bring about a better understanding of ecological and evolutionary processes structuring genetic variation across space and time.

Tri-mode optical biopsy probe with fluorescence endomicroscopy, Raman spectroscopy, and time-resolved fluorescence spectroscopy.

We present an endoscopic probe that combines three distinct optical fibre technologies including: A high-resolution imaging fibre for optical endomicroscopy, a multimode fibre for time-resolved fluorescence spectroscopy, and a hollow-core fibre with multimode signal collection cores for Raman spectroscopy. The three fibers are all enclosed within a 1.2 mm diameter clinical grade catheter with a 1.4 mm end cap. To demonstrate the probe's flexibility we provide data acquired with it in loops of radii down to 2 cm. We then use the probe in an anatomically accurate model of adult human airways, showing that it can be navigated to any part of the distal lung using a commercial bronchoscope. Finally, we present data acquired from fresh ex vivo human lung tissue. Our experiments show that this minimally invasive probe can deliver real-time optical biopsies from within the distal lung - simultaneously acquiring co-located high-resolution endomicroscopy and biochemical spectra.

Miniaturised gap sensor using fibre optic Fabry-Pérot interferometry for structural health monitoring.

A miniaturised structural health monitoring device has been developed capable of measuring the absolute distance between close parallel surfaces using Fabry-Pérot interferometry with nm-scale sensitivity. This is achieved by fabricating turning mirrors on two opposite cores of a multi-core fibre to produce a probe with dimensions limited only be the fibre diameter. Two fabrication processes have been investigated: Focused ion beam milling, which has resulted in a sensor measurement accuracy, sensitivity and range of ±0.056 µm, ±0.006 µm and ∼16000  µm respectively; and ultrafast laser assisted etching of the cleaved fibre end, where a sensor measurement accuracy, sensitivity and range of ±0.065 µm, ±0.006 µm and ∼7500 µm have been demonstrated.

Finite element human body models with active reflexive muscles suitable for sex based whiplash injury prediction.

Previous research has not produced a satisfactory resource to study reflexive muscle activity for investigating potentially injurious whiplash motions. Various experimental and computational studies are available, but none provided a comprehensive biomechanical representation of human response during rear impacts. Three objectives were addressed in the current study to develop female and male finite element human body models with active reflexive neck muscles: 1) eliminate the buckling in the lower cervical spine of the model observed in earlier active muscle controller implementations, 2) evaluate and quantify the influence of the individual features of muscle activity, and 3) evaluate and select the best model configuration that can be used for whiplash injury predictions. The current study used an open-source finite element model of the human body for injury assessment representing an average 50th percentile female anthropometry, together with the derivative 50th percentile male morphed model. Based on the head-neck kinematics and CORelation and Analyis (CORA) tool for evaluation, models with active muscle controller and parallel damping elements showed improved head-neck kinematics agreement with the volunteers over the passive models. It was concluded that this model configuration would be the most suitable for gender-based whiplash injury prediction when different impact severities are to be studied.

Investigating focus elongation using a spatial light modulator for high-throughput ultrafast-laser-induced selective etching in fused silica.

Ultrafast-laser-induced selective chemical etching is an enabling microfabrication technology compatible with optical materials such as fused silica. The technique offers unparalleled three-dimensional manufacturing freedom and feature resolution but can be limited by long laser inscription times and widely varying etching selectivity depending on the laser irradiation parameters used. In this paper, we aim to overcome these limitations by employing beam shaping via a spatial light modulator to generate a vortex laser focus with controllable depth-of-focus (DOF), from diffraction limited to several hundreds of microns. We present the results of a thorough parameter-space investigation of laser irradiation parameters, documenting the observed influence on etching selectivity and focus elongation in the polarization-insensitive writing regime, and show that etching selectivity greater than 800 is maintained irrespective of the DOF. To demonstrate high-throughput laser writing with an elongated DOF, geometric shapes are fabricated with a 12-fold reduction in writing time compared to writing with a phase-unmodulated Gaussian focus.

Breeding near heterospecifics as a defence against brood parasites: can redstarts lower probability of cuckoo parasitism using neighbours?

Breeding habitat choice based on the attraction to other species can provide valuable social information and protection benefits. In birds, species with overlapping resources can be a cue of good quality habitats; species with shared predators and/or brood parasites can increase joint vigilance or cooperative mobbing, while raptors may provide a protective umbrella against these threats. We tested whether the migratory common redstart (Phoenicurus phoenicurus) is attracted to breed near active nests of the great tit (Parus major), a keystone-information source for migrant passerine birds, or a top predator, the northern goshawk (Accipiter gentilis). This system is unique to test these questions because the redstart is a regular host for the common cuckoo (Cuculus canorus). Therefore, we also evaluated other possible benefits coming from the heterospecific attraction, especially in terms of reducing brood parasitism risk. We monitored redstart occupancy rates, onset of breeding, reproductive investment, and followed nest outcomes in terms of brood parasitism, nest predation risk and overall reproductive success. Redstarts avoided breeding near goshawks, but showed neither attraction nor avoidance to breed next to great tits. Both neighbours neither reduced brood parasitism risk nor affected overall nesting success in redstarts. Redstarts may not use heterospecific attraction for settlement decisions, as associations with other species can only exist when some benefits are gained. Thus, environmental cues may be more important than social information for redstarts when breeding habitat choice. Other front-line defence strategies may have a better impact reducing breeding negative interactions, such brood parasitism.

Analysis of control strategies for VIVA OpenHBM with active reflexive neck muscles.

Modeling muscle activity in the neck muscles of a finite element (FE) human body model can be based on two biological reflex systems. One approach is to approximate the Vestibulocollic reflex (VCR) function, which maintains the head orientation relative to a fixed reference in space. The second system tries to maintain the head posture relative to the torso, similar to the Cervicocolic reflex (CCR). Strategies to combine these two neck muscle controller approaches in a single head-neck FE model were tested, optimized, and compared to rear-impact volunteer data. The first approach, Combined-Control, assumed that both controllers simultaneously controlled all neck muscle activations. In the second approach, Distributed-Control, one controller was used to regulate activation of the superficial muscles while a different controller acted on deep neck muscles. The results showed that any muscle controller that combined the two approaches was less effective than only using one of VCR- or CCR-based systems on its own. A passive model had the best objective rating for cervical spine kinematics, but the addition of a single active controller provided the best response for both head and cervical spine kinematics. The present study demonstrates the difficulty in completely capturing representative head and cervical spine responses to rear-impact loading and identified a controller capturing the VCR reflex as the best candidate to investigate whiplash injury mechanisms through FE modeling.

The evolutionary network of whiptail lizards reveals predictable outcomes of hybridization.

Hybridization between diverging lineages is associated with the generation and loss of species diversity, introgression, adaptation, and changes in reproductive mode, but it is unknown when and why it results in these divergent outcomes. We estimate a comprehensive evolutionary network for the largest group of unisexual vertebrates and use it to understand the evolutionary outcomes of hybridization. Our results show that rates of introgression between species decrease with time since divergence and suggest that species must attain a threshold of evolutionary divergence before hybridization results in transitions to unisexuality. Rates of hybridization also predict genome-wide patterns of genetic diversity in whiptail lizards. These results distinguish among models for hybridization that have not previously been tested and suggest that the evolutionary outcomes can be predictable.

A LAMP sequencing approach for high-throughput co-detection of SARS-CoV-2 and influenza virus in human saliva.

The COVID-19 pandemic has created an urgent need for rapid, effective, and low-cost SARS-CoV-2 diagnostic testing. Here, we describe COV-ID, an approach that combines RT-LAMP with deep sequencing to detect SARS-CoV-2 in unprocessed human saliva with a low limit of detection (5-10 virions). Based on a multi-dimensional barcoding strategy, COV-ID can be used to test thousands of samples overnight in a single sequencing run with limited labor and laboratory equipment. The sequencing-based readout allows COV-ID to detect multiple amplicons simultaneously, including key controls such as host transcripts and artificial spike-ins, as well as multiple pathogens. Here, we demonstrate this flexibility by simultaneous detection of 4 amplicons in contrived saliva samples: SARS-CoV-2, influenza A, human STATHERIN, and an artificial SARS calibration standard. The approach was validated on clinical saliva samples, where it showed excellent agreement with RT-qPCR. COV-ID can also be performed directly on saliva absorbed on filter paper, simplifying collection logistics and sample handling.