Vincent Castella

Reliable selfing rate estimates from imperfect population genetic data

Genotypic frequencies at codominant marker loci in population samples convey information on mating systems. A classical way to extract this information is to measure heterozygote deficiencies (FIS) and obtain the selfing rate s from FIS = s/(2 − s), assuming inbreeding equilibrium. A major drawback is that heterozygote deficiencies are often present without selfing, owing largely to technical artefacts such as null alleles or partial dominance. We show here that, in the absence of gametic disequilibrium, the multilocus structure can be used to derive estimates of s independent of FIS and free of technical biases. Their statistical power and precision are comparable to those of FIS, although they are sensitive to certain types of gametic disequilibria, a bias shared with progeny‐array methods but not FIS. We analyse four real data sets spanning a range of mating systems. In two examples, we obtain s = 0 despite positive FIS, strongly suggesting that the latter are artefactual. In the remaining examples, all estimates are consistent. All the computations have been implemented in a open‐access and user‐friendly software called rmes (robust multilocus estimate of selfing) available at http://ftp.cefe.cnrs.fr, and can be used on any multilocus data. Being able to extract the reliable information from imperfect data, our method opens the way to make use of the ever‐growing number of published population genetic studies, in addition to the more demanding progeny‐array approaches, to investigate selfing rates.

Is the Gibraltar Strait a barrier to gene flow for the bat Myotis myotis (Chiroptera: Vespertilionidae)?

Because of their role in limiting gene flow, geographical barriers like mountains or seas often coincide with intraspecific genetic discontinuities. Although the Strait of Gibraltar represents such a potential barrier for both plants and animals, few studies have been conducted on its impact on gene flow. Here we test this effect on a bat species (Myotis myotis) which is apparently distributed on both sides of the strait. Six colonies of 20 Myotis myotis each were sampled in southern Spain and northern Morocco along a linear transect of 1350 km. Results based on six nuclear microsatellite loci reveal no significant population structure within regions, but a complete isolation between bats sampled on each side of the strait. Variability at 600 bp of a mitochondrial gene (cytochrome b) confirms the existence of two genetically distinct and perfectly segregating clades, which diverged several million years ago. Despite the narrowness of the Gibraltar Strait (14 km), these molecular data suggest that neither males, nor females from either region have ever reproduced on the opposite side of the strait. Comparisons of molecular divergence with bats from a closely related species (M. blythii) suggest that the North African clade is possibly a distinct taxon warranting full species rank. We provisionally refer to it as Myotis cf punicus Felten 1977, but a definitive systematic understanding of the whole Mouse‐eared bat species complex awaits further genetic sampling, especially in the Eastern Mediterranean areas.

Contrasted patterns of mitochondrial and nuclear structure among nursery colonies of the bat Myotis myotis

Thirteen nursery colonies of Myotis myotis were sampled in central Europe to investigate the dispersal behaviour of this bat species. Mitochondrial DNA sequences of 260 bats reveal the occurrence of three evolutionary lineages that have probably originated in distinct glacial refugia and meet in a contact zone near the Alps. Moreover, the strong haplotypic segregation (ΦST=0.540) suggests that breeding females are philopatric. Contrastingly, the low population structure at 15 microsatellite loci (FST=0.022), suggests the homogenizing effect of nuclear gene flow. The different perspectives given by these two markers are consistent with strong male‐biased dispersal. As a result of female philopatry, the local haplotypic variability seems to be largely influenced by historical processes of colonization. Conversely, the homogeneity of nuclear variability within roosts that are located north of the Alps seems to mainly reflect contemporary gene flow. Finally, despite the fact that females are faithful to their natal colony, movements of both males and females occur outside the breeding period. Mitochondrial survey of individuals sampled exclusively in nurseries may thus poorly reflect the metapopulation dynamics of this species.

Genetic consequences of the ice ages on nurseries of the bat Myotis myotis: a mitochondrial and nuclear survey

Analyses of mitochondrial DNA (mtDNA) control region polymorphism and of variation at 10 nuclear microsatellite loci were used to investigate the mechanisms and genetic consequences of postglacial expansion of Myotis myotis in Europe. Initial sampling consisted of 480 bats genotyped in 24 nursery colonies arranged along a transect of ≈ 3000 km. The phylogeographical survey based on mtDNA sequences revealed the existence of major genetic subdivisions across this area, with several suture zones between haplogroups. Such zones of secondary contact were found in the Alps and Rhodopes, whereas other potential barriers to gene flow, like the Pyrenees, did not coincide with genetic discontinuities. Areas of population admixture increased locally the genetic diversity of colonies, which confounded the northward decrease in nucleotide diversity predicted using classical models of postglacial range expansion. However, when analyses were restricted to a subset of 15 nurseries originating from a single presumed glacial refugium, mtDNA polymorphism did indeed support a northwards decrease in diversity. Populations were also highly structured (ΦST = 0.384). Conversely, the same subset of colonies showed no significant latitudinal decrease in microsatellite diversity and much less population structure (FST = 0.010), but pairwise genetic differentiation at these nuclear markers was strongly correlated with increasing geographical distance. Together, this evidence suggests that alleles carried via male bats have maintained enough nuclear gene flow to counteract the effects of recurrent bottlenecks generally associated with recolonization processes. As females are highly philopatric, we argue that the maternally transmitted mtDNA marker better reflects the situation of past, historical gene flow, whereas current levels of gene flow are better reflected by microsatellite markers.

Differential DNA extraction of challenging simulated sexual-assault samples: a Swiss collaborative study

In sexual-assault cases, autosomal DNA analysis of gynecological swabs is a challenge, as the presence of a large quantity of female material may prevent detection of the male DNA. A solution to this problem is differential DNA extraction, but there is no established best practice for this. We decided to test the efficacy of a number of different protocols on simulated casework samples. Four difficult samples were sent to the nine Swiss laboratories active in forensic genetics. In each laboratory, staff used their routine protocols to separate the epithelial-cell fraction, enriched with the non-sperm DNA, from the sperm fraction. DNA extracts were then sent to the organizing laboratory for analysis. Estimates of male:female DNA ratio without differential DNA extraction ranged from 1:38 to 1:339, depending on the semen used to prepare the samples. After differential DNA extraction, most of the ratios ranged from 1:12 to 9:1, allowing detection of the male DNA. Compared with direct DNA extraction, cell separation resulted in losses of 94-98% of the male DNA. As expected, more male DNA was generally present in the sperm than in the epithelial-cell fraction. However, for about 30% of the samples, the reverse trend was seen. The recovery of male and female DNA was highly variable, depending on the laboratory involved. An experimental design similar to the one used in this study may be of assistance for local protocol testing and improvement.

Species identification in mammals from mixed biological samples based on mitochondrial DNA control region length polymorphism

The ability to identify the species origin of an unknown biological sample is relevant in the fields of human and wildlife forensics. However, the detection of several species mixed in the same sample still remains a challenge. We developed and tested a new approach for mammal DNA identification in mixtures of two or three species, based on the analysis of mitochondrial DNA control region interspecific length polymorphism followed by direct sequencing. Contrary to other published methods dealing with species mixtures, our protocol requires a single universal primer pair and is not based on a pre‐defined panel of species. Amplicons can be separated either on agarose gels or using CE. The advantages and limitations of the assay are discussed under different conditions, such as variable template concentration, amplicon sizes and size difference among the amplicons present in the mixture. For the first time, this protocol provides a simple, reliable and flexible method for simultaneous identification of multiple mammalian species from mixtures, without any prior knowledge of the species involved.

DIP–STR: Highly Sensitive Markers for the Analysis of Unbalanced Genomic Mixtures

Samples containing highly unbalanced DNA mixtures from two individuals commonly occur both in forensic mixed stains and in peripheral blood DNA microchimerism induced by pregnancy or following organ transplant. Because of PCR amplification bias, the genetic identification of a DNA that contributes trace amounts to a mixed sample represents a tremendous challenge. This means that standard genetic markers, namely microsatellites, also referred as short tandem repeats (STR), and single‐nucleotide polymorphism (SNP) have limited power in addressing common questions of forensic and medical genetics. To address this issue, we developed a molecular marker, named DIP–STR that relies on pairing deletion–insertion polymorphisms (DIP) with STR. This novel analytical approach allows for the unambiguous genotyping of a minor component in the presence of a major component, where DIP–STR genotypes of the minor were successfully procured at ratios up to 1:1,000. The compound nature of this marker generates a high level of polymorphism that is suitable for identity testing. Here, we demonstrate the power of the DIP–STR approach on an initial set of nine markers surveyed in a Swiss population. Finally, we discuss the limitations and potential applications of our new system including preliminary tests on clinical samples and estimates of their performance on simulated DNA mixtures.

Genetic identification of decomposed cadavers using nails as DNA source

Blood or muscle can be used as a DNA source for the genetic identification of recently deceased persons. If the post mortem interval increases, bones and teeth are used. In this case, collection and DNA isolation will be more difficult and time consuming. The aim of this study was to evaluate the use of nails as an alternative DNA source for the genetic identification of decomposed cadavers. DNA extraction from 5mg of fingernails from 7 volunteers using 1h cell lysis in a standard buffer and a DNA purification on QIAamp DNA mini kit columns allowed to acquire a mean quantity of 100 ng DNA/mg nail. This was unexpected, as blood and muscle contain comparable amounts of DNA. Our protocol allowed to obtain full PowerPlex 16 DNA profiles from 10 cadavers characterized by post mortem intervals ranging from 5 days to more than 6 months. The good quality of these profiles indicated that DNA from nail is well preserved. In conclusion, nails are very easy to collect and contain large amounts of good quality DNA that can be extracted within a few hours. They may therefore represent an attractive DNA source not only for routine, but also for urgent genetic identification of decomposed cadavers.

A novel set of DIP-STR markers for improved analysis of challenging DNA mixtures

Currently available molecular biology tools allow forensic scientists to characterize DNA evidence found at crime scenes for a large variety of samples, including those of limited quantity and quality, and achieve high levels of individualization. Yet, standard forensic markers provide limited or no results when applied to mixed DNA samples where the contributors are present in very different proportions (unbalanced DNA mixtures). This becomes an issue mostly for the analysis of trace samples collected on the victim or from touched objects. To this end, we recently proposed an innovative type of genetic marker, named DIP-STR that relies on pairing deletion/insertion polymorphisms (DIP) with standard short tandem repeats (STR). This novel compound marker allows detection of the minor DNA contributor in a DNA mixture of any gender and cellular origin with unprecedented resolution (beyond a DNA ratio of 1:1000). To provide a novel analytical tool useful in practice to common forensic laboratories, this article describes the first set of 10 DIP-STR markers selected according to forensic technical standards. The novel DIP-STR regions are short (between 146 and 271 bp), include only highly polymorphic tri-, tetra- and pentanucleotide tandem repeats and are located on different chromosomes or chromosomal arms to provide statistically independent results. This novel set of DIP-STR can target the amplification of 0.03-0.1 ng of DNA when mixed with a 1000-fold excess of major DNA. DIP-STR relative allele frequencies are estimated based on a survey of 103 Swiss individuals. Finally, this study provides an estimate of the occurrence of informative alleles and a calculation of the corresponding random match probability of the detected minor DIP-STR genotype assessed across 10,506 pairwise conceptual mixtures.

Forensic identification of urine samples: a comparison between nuclear and mitochondrial DNA markers

Urine samples from 20 male volunteers of European Caucasian origin were stored at 4°C over a 4-month period in order to compare the identification potential of nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) markers. The amount of nDNA recovered from urines dramatically declined over time. Consequently, nDNA likelihood ratios (LRs) greater than 1,000 were obtained for 100, 70 and 55% of the urines analysed after 6, 60 and 120 days, respectively. For the mtDNA, HVI and HVII sequences were obtained for all samples tested, whatever the period considered. Nevertheless, the highest mtDNA LR of 435 was relatively low compared to its nDNA equivalent. Indeed, LRs obtained with only three nDNA loci could easily exceed this value and are quite easier to obtain. Overall, the joint use of nDNA and mtDNA markers enabled the 20 urine samples to be identified, even after the 4-month period.