Approximate Bayesian Computation Untangles Signatures of Contemporary and Historical Hybridization between Two Endangered Species

Contemporary gene flow, when resumed after a period of isolation, can have crucial consequences for endangered species, as it can both increase the supply of adaptive alleles and erode local adaptation. Determining the history of gene flow and thus the importance of contemporary hybridization, however, is notoriously difficult. Here, we focus on two endangered plant species, Arabis nemorensis and A. sagittata, which hybridize naturally in a sympatric population located on the banks of the Rhine. Using reduced genome sequencing, we determined the phylogeography of the two taxa but report only a unique sympatric population. Molecular variation in chloroplast DNA indicated that A. sagittata is the principal receiver of gene flow. Applying classical D-statistics and its derivatives to whole-genome data of 35 accessions, we detect gene flow not only in the sympatric population but also among allopatric populations. Using an Approximate Bayesian computation approach, we identify the model that best describes the history of gene flow between these taxa. This model shows that low levels of gene flow have persisted long after speciation. Around 10 000 years ago, gene flow stopped and a period of complete isolation began. Eventually, a hotspot of contemporary hybridization was formed in the unique sympatric population. Occasional sympatry may have helped protect these lineages from extinction in spite of their extremely low diversity.     

A replicated climate change field experiment reveals rapid evolutionary response in an ecologically important soil invertebrate

Whether species can respond evolutionarily to current climate change is crucial for the persistence of many species. Yet, very few studies have examined genetic responses to climate change in manipulated experiments carried out in natural field conditions. We examined the evolutionary response to climate change in a common annelid worm using a controlled replicated experiment where climatic conditions were manipulated in a natural setting. Analyzing the transcribed genome of 15 local populations, we found that about 12% of the genetic polymorphisms exhibit differences in allele frequencies associated to changes in soil temperature and soil moisture. This shows an evolutionary response to realistic climate change happening over short‐time scale, and calls for incorporating evolution into models predicting future response of species to climate change. It also shows that designed climate change experiments coupled with genome sequencing offer great potential to test for the occurrence (or lack) of an evolutionary response.     

Inherent errors of the fixed-frame counting method for corneal endothelial cell density in eye banks

The aim of this work was to analyze the magnitude of inherent errors associated with the fixed-frame counting method for corneal endothelial cell density (ECD) measurements. This technique is common among most eye banks worldwide. Three types of mosaics were used: regular and irregular tessellated mosaics (eight increasing densities ranging from 800 to 3,600 cells/mm² by steps of 400 cells/mm²) generated by a computer, and real mosaics (four specimens) obtained from human corneal endothelium flat mounted and stained with Alizarin red. On the three mosaics, the fixed-frame counting method was applied using a computer program. The ECD was calculated for 3,000 successive random positions from calibrated grids which area ranged from 50 × 50 to 300 × 300 μm² (incremental steps of 25 μm). For each grid, the ECD was expressed either as a single count, a mean of five or a mean of 10 measures. The fixed-frame count was constantly associated with an inherent variability but repeatability increased with larger grid size and ECD. The mean calculated out of 10 measures was the most reliable, but still, we noted ±5 % of residual variability from the real ECD. The 100 × 100 μm² grid manual counts, performed in many eye banks, should be abandoned and upgraded to at least 200 × 200 μm² grid counts. Digital image analysis with a variable frame counting method would be the best alternative.     

Environmental metagenomics: An innovative resource for industrial biocatalysis

The current existing enzymes have been identified from cultivable micro-organisms, most frequently from bacteria. These bacterial biocatalytic capabilities have been widely used for biotransformations, resulting in the development of profitable industrial bioprocesses in the fields of feed and food processing, textiles, agro-chemistry, cosmetics, pharmaceuticals and fine chemistry. However, the originality of this bioresource is progressively drying up, while requests from industry for novel biocatalytic activities are increasing in the face of economic and environmental pressure. Metagenomics, through access to the huge reservoir of uncultivated bacteria which represents the majority of the present biodiversity, opens the door to new industrial sources of enzymes. Surmounting hurdles encountered with this technology (e.g. DNA extraction to obtain high quality DNA libraries with proper statistical representativity, setting up of relevant high throughput screenings assays, combining functional and genome-based identifications), gives unique opportunities to access novel biocatalysts that better fit with the required industrial specifications, thus providing new biocatalysis tool boxes.     

3D Evaluation of the acetabular coverage assessed by biplanar X-rays or single anteroposterior X-ray compared with CT-scan

Considering the increasing development of three dimensional (3D) imaging, the 3D assessment of the acetabular coverage is to become the most interesting tool for the detection of acetabular pathologies. Biplanar X-rays based methods allow a 3D reconstruction of the hip with a reduced radiation dose. This study proposes a 3D assessment method of the acetabular coverage from biplanar X-rays or from an anteroposterior X-ray (conventional clinical imaging). An in vitro evaluation of the method was performed on six hip joints in comparison with computed tomography. The global coverage, the local coverage and the acetabular rim orientation were estimated in 3D. The mean global acetabular coverage was 40% with an estimated mean accuracy of 1.3% for the biplanar X-rays based method. This study evaluated a 3D assessment method of the acetabular coverage from biplanar X-rays or anteroposterior X-ray and open the way for clinical in vivo applications.     

Competing risks regression for stratified data

For competing risks data, the Fine-Gray proportional hazards model for subdistribution has gained popularity for its convenience in directly assessing the effect of covariates on the cumulative incidence function. However, in many important applications, proportional hazards may not be satisfied, including multicenter clinical trials, where the baseline subdistribution hazards may not be common due to varying patient populations. In this article, we consider a stratified competing risks regression, to allow the baseline hazard to vary across levels of the stratification covariate. According to the relative size of the number of strata and strata sizes, two stratification regimes are considered. Using partial likelihood and weighting techniques, we obtain consistent estimators of regression parameters. The corresponding asymptotic properties and resulting inferences are provided for the two regimes separately. Data from a breast cancer clinical trial and from a bone marrow transplantation registry illustrate the potential utility of the stratified Fine-Gray model.     

Seqcrawler: biological data indexing and browsing platform

ABSTRACT: BACKGROUND: Seqcrawler takes its roots in software like SRS or Lucegene. It provides an indexing platform to ease the search of data and meta-data in biological banks and it can scale to face the current flow of data. While many biological bank search tools are available on the Internet, mainly provided by large organizations to search in their data, there is a lack of free and open source solution to browse one own set of data with a flexible query system and able to scale from single computer to a cloud system. A personal index platform will help labs and bioinformaticians to search in their meta-data but also to build a larger information system with custom subsets of data. RESULTS: The software is scalable from a single computer to a cloud-based infrastructure. It has been successfully tested in a private cloud with 3 index shards (piece of index) hosting ~400 millions of sequence information (whole GenBank, UniProt, PDB and others) for a total size of 600 GB in a fault tolerant architecture (high-availability). It has also been successfully integrated with software to add extra meta-data from blast results to enhance user's result analysis. CONCLUSIONS: Seqcrawler provides a complete open source search and store solution for labs or platforms needing to manage large amount of data/meta-data with a flexible and customizable web interface. All components (search engine, visualization and data storage), though independent, share a common and coherent data system that can be queried with a simple HTTP interface. The solution scales easily and can also provide a high availability infrastructure.     

A model for the evolution of crèching behaviour in gulls

A crèche is an aggregation of chicks outside nesting territories, within chicks continue to be fed only by their own parents. Several adaptive functions of crèching have been proposed, the most frequent being a reduction in predator pressure. Using an evolutionary stable strategy approach based on the computation of individuals' fecundity, we examined which regime of aerial and terrestrial predation is likely to favour the evolution and stability of the crèching strategy (CS) in gulls. Our results confirm the hypothesis that habitat instability associated with high levels of terrestrial predation favours the evolution and maintenance of crèching behaviour. Moreover, our results suggest that a low aggressiveness against predators may be a pre-adaptation to a CS. In contrast, the high synchronisation often observed in crèching species does not favour the evolution of a crèching behaviour and is thus probably under selection pressures different from those modelled here.