A holistic approach for the selection of forensic DNA swabs.

In the present study, we compared the performance of five different ISO 18385 certified forensic swabs for DNA sampling in practice over a time period of five months. Comparisons were made for DNA profiling success rates, measured as the percentage of CODIS (Combined DNA Index System) suitable profiles as well as for practical suitability during sampling at the scene, measured through a survey among collaborators. More than forty members of our crime scene investigation (CSI) unit took part in the test series and provided structured feedback concerning different aspects of swab handling. A total number of 1094 "touch" DNA samples have been subjected to DNA analysis. Swabs performed significantly different in terms of DNA profiling success rates. We also observed significant differences in DNA extraction efficiency between swabs. The evaluation by the collaborators of various aspects of handling differed significantly between swabs. We can assume that a more convenient handling decreases the risk of contamination or sample mislabelling and increases sampling efficiency and staff satisfaction. Our results demonstrate that the selection of disposable sampling devices such as forensic swabs for DNA sampling should be made based on a holistic approach. To be able to select the best performing swab for a given combination of CSI and DNA laboratory procedures, it might not be sufficient to only perform DNA extraction comparisons and trace sampling under controlled laboratory conditions.

Human background DNA on stones in an urban environment.

Stones are frequently used as tools in criminal acts. In our department, around 5 % of all analysed crime scene related trace samples are contact or touch DNA traces swabbed from stones. These samples are primarily related to cases of damage to property and burglary. In court, questions can arise about DNA transfer and the persistence of background DNA not related to the respective crime. To shed some light on the question of how likely it is to detect human DNA as background DNA on stones from an urban environment, the surfaces of 108 stones sampled throughout the city of Bern, the Swiss capital, were swabbed. We detected a median quantity of 33 pg on the sampled stones. STR-profiles suitable for a CODIS (Combined DNA Index System) registration in the Swiss DNA database were established from 6.5 % of all sampled stone surfaces. For comparison, retrospective casework data analysis from routine crime scene samples demonstrates a success rate of 20.6 % for the establishment of CODIS-suitable DNA profiles from stones sampled for touch DNA. We further investigated how climatic conditions, location and properties of the stones affected the quantity and quality of the recovered DNA. In this study, we show that the quantity of the measurable DNA decreases significantly with increasing temperature. Furthermore, less DNA could be recovered from porous stones, compared to smooth ones.

Population genetic analysis of 12 X-chromosomal STRs in a Swiss sample.

X-chromosomal STRs are a powerful tool to assess a broad variety of complex kinship scenarios. We introduce herewith the first Swiss X-STR dataset based on 1198 individuals (592 female, 606 male), characterized with the Qiagen Investigator® Argus X-12 QS multiplex kit. Anomalous allele patterns, allele and haplotype frequencies, and forensic and population genetic parameters are presented. We detected linkage disequilibrium within three out of the four designated linkage groups and no apparent intra-national population substructure. We compared the dataset to a global panel of X-STR datasets and it fits well in the European context, as expected.

Expanding the Swiss autosomal marker set to 32 STRs.

By genotyping 1198 individuals with the Qiagen Investigator® HDplex Kit, we expand the Swiss autosomal STR dataset to 32 loci, providing additional resources for complex kinship cases. We present the first high-quality allele frequency dataset for loci D2S1360, D5S2500, D7S1517, and D10S2325 that will be accessible through the ENFSI reference database STRidER. For loci D3S1744, D4S2366, D6S474, D8S1132, and D21S2055, we provide a first European STRidER dataset.

fastsimcoal2: demographic inference under complex evolutionary scenarios.

MOTIVATION: fastsimcoal2 extends fastsimcoal, a continuous time coalescent-based genetic simulation program, by enabling the estimation of demographic parameters under very complex scenarios from the site frequency spectrum under a maximum-likelihood framework. RESULTS: Other improvements include multi-threading, handling of population inbreeding, extended input file syntax facilitating the description of complex demographic scenarios, and more efficient simulations of sparsely structured populations and of large chromosomes. AVAILABILITY AND IMPLEMENTATION: fastsimcoal2 is freely available on http://cmpg.unibe.ch/software/fastsimcoal2/. It includes console versions for Linux, Windows and MacOS, additional scripts for the analysis and visualization of simulated and estimated scenarios, as well as a detailed documentation and ready-to-use examples.

Recent Evolutionary History of Tigers Highlights Contrasting Roles of Genetic Drift and Selection.

Species conservation can be improved by knowledge of evolutionary and genetic history. Tigers are among the most charismatic of endangered species and garner significant conservation attention. However, their evolutionary history and genomic variation remain poorly known, especially for Indian tigers. With 70% of the world's wild tigers living in India, such knowledge is critical. We re-sequenced 65 individual tiger genomes representing most extant subspecies with a specific focus on tigers from India. As suggested by earlier studies, we found strong genetic differentiation between the putative tiger subspecies. Despite high total genomic diversity in India, individual tigers host longer runs of homozygosity, potentially suggesting recent inbreeding or founding events, possibly due to small and fragmented protected areas. We suggest the impacts of ongoing connectivity loss on inbreeding and persistence of Indian tigers be closely monitored. Surprisingly, demographic models suggest recent divergence (within the last 20,000 years) between subspecies and strong population bottlenecks. Amur tiger genomes revealed the strongest signals of selection related to metabolic adaptation to cold, whereas Sumatran tigers show evidence of weak selection for genes involved in body size regulation. We recommend detailed investigation of local adaptation in Amur and Sumatran tigers prior to initiating genetic rescue.

Polygenic Patterns of Adaptive Introgression in Modern Humans Are Mainly Shaped by Response to Pathogens.

Anatomically modern humans carry many introgressed variants from other hominins in their genomes. Some of them affect their phenotype and can thus be negatively or positively selected. Several individual genes have been proposed to be the subject of adaptive introgression, but the possibility of polygenic adaptive introgression has not been extensively investigated yet. In this study, we analyze archaic introgression maps with refined functional enrichment methods to find signals of polygenic adaptation of introgressed variants. We first apply a method to detect sets of connected genes (subnetworks) within biological pathways that present higher-than-expected levels of archaic introgression. We then introduce and apply a new statistical test to distinguish between epistatic and independent selection in gene sets of present-day humans. We identify several known targets of adaptive introgression, and we show that they belong to larger networks of introgressed genes. After correction for genetic linkage, we find that signals of polygenic adaptation are mostly explained by independent and potentially sequential selection episodes. However, we also find some gene sets where introgressed variants present significant signals of epistatic selection. Our results confirm that archaic introgression has facilitated local adaptation, especially in immunity related and metabolic functions and highlight its involvement in a coordinated response to pathogens out of Africa.

Relaxed Selection During a Recent Human Expansion.

Humans have colonized the planet through a series of range expansions, which deeply impacted genetic diversity in newly settled areas and potentially increased the frequency of deleterious mutations on expanding wave fronts. To test this prediction, we studied the genomic diversity of French Canadians who colonized Quebec in the 17th century. We used historical information and records from ∼4000 ascending genealogies to select individuals whose ancestors lived mostly on the colonizing wave front and individuals whose ancestors remained in the core of the settlement. Comparison of exomic diversity reveals that: (i) both new and low-frequency variants are significantly more deleterious in front than in core individuals, (ii) equally deleterious mutations are at higher frequencies in front individuals, and (iii) front individuals are two times more likely to be homozygous for rare very deleterious mutations present in Europeans. These differences have emerged in the past six to nine generations and cannot be explained by differential inbreeding, but are consistent with relaxed selection mainly due to higher rates of genetic drift on the wave front. Demographic inference and modeling of the evolution of rare variants suggest lower effective size on the front, and lead to an estimation of selection coefficients that increase with conservation scores. Even though range expansions have had a relatively limited impact on the overall fitness of French Canadians, they could explain the higher prevalence of recessive genetic diseases in recently settled regions of Quebec.

Detecting gene subnetworks under selection in biological pathways.

Advances in high throughput sequencing technologies have created a gap between data production and functional data analysis. Indeed, phenotypes result from interactions between numerous genes, but traditional methods treat loci independently, missing important knowledge brought by network-level emerging properties. Therefore, detecting selection acting on multiple genes affecting the evolution of complex traits remains challenging. In this context, gene network analysis provides a powerful framework to study the evolution of adaptive traits and facilitates the interpretation of genome-wide data. We developed a method to analyse gene networks that is suitable to evidence polygenic selection. The general idea is to search biological pathways for subnetworks of genes that directly interact with each other and that present unusual evolutionary features. Subnetwork search is a typical combinatorial optimization problem that we solve using a simulated annealing approach. We have applied our methodology to find signals of adaptation to high-altitude in human populations. We show that this adaptation has a clear polygenic basis and is influenced by many genetic components. Our approach, implemented in the R package signet, improves on gene-level classical tests for selection by identifying both new candidate genes and new biological processes involved in adaptation to altitude.

STRAF-A convenient online tool for STR data evaluation in forensic genetics.

Population data in forensic genetics has to be checked for a variety of statistical parameters before it can be employed for case work. A lot of very powerful statistical tools are available for this task, most of them developed by labs having their research focus on population genetics or evolution. However, most of these programs require a substantial amount of experience. In addition, to our knowledge, none of the freely available programs calculates all the common parameters for a population study in forensic genetics at once, based on a single input file. We present here a convenient online tool that fills this gap. STRAF (STR Analysis for Forensics) provides an intuitive interface and input file format and computes all the relevant parameters for a classical population study based on autosomal STR data at once and in a convenient way. In addition, STRAF includes a PCA module that can be used for population substructure detection or quality control. The results generated by the program were verified by recalculating parameters from an already published population study.

The summary-likelihood method and its implementation in the Infusion package.

In recent years, simulation methods such as approximate Bayesian computation have extensively been used to infer parameters of population genetic models where the likelihood is intractable. We describe an alternative approach, summary likelihood, that provides a likelihood-based analysis of the information retained in the summary statistics whose distribution is simulated. We provide an automated implementation as a standard R package, Infusion, and we test the method, in particular for a scenario of inference of population-size change from genetic data. We show that the method provides confidence intervals with controlled coverage independently of a prior distribution on parameters, in contrast to approximate Bayesian computation. We expect the method to be applicable for at least six-parameter models and discuss possible modifications for higher-dimensional inference problems.