Search

We investigated the controversial origin of domestic sheep (Ovis aries) using large samples of contemporary and ancient domestic individuals and their closest wild relatives: the Asiatic mouflon (Ovis gmelini), the urial (Ovis vignei) and the argali (Ovis ammon). A phylogeny based on mitochondrial DNA, including 213 new cytochrome-b sequences of wild Ovism confirmed that O. gmelini is the maternal ancestor of sheep and precluded mtDNA contributions from O. vignei (and O. gmelini × O. vignei hybrids) to domestic lineages. We also produced 54 new control region sequences showing shared haplogroups (A, B, C and E) between domestic sheep and wild O. gmelini which localized the domestication center in eastern Anatolia and central Zagros, excluding regions further east where exclusively wild haplogroups were found. This overlaps with the geographic distribution of O. gmelini gmelini, further suggesting that the maternal origin of domestic sheep derives from this subspecies. Additionally, we produced 57 new CR sequences of Neolithic sheep remains from a large area covering Anatolia to Europe, showing the early presence of at least three mitochondrial haplogroups (A, B and D) in Western colonization routes. This confirmed that sheep domestication was a large-scale process that captured diverse maternal lineages (haplogroups).

The way in which living organisms mobilize a combination of long-term adaptive mechanisms and short-term phenotypic plasticity to face environmental variations is still largely unknown. In the context of climate change, understanding the genetic and epigenetic bases for adaptation and plasticity is a major stake for preserving genomic resources and the resilience capacity of livestock populations. We characterized both epigenetic and genetic variations by contrasting 22 sheep and 21 goats from both sides of a climate gradient, focusing on free-ranging populations from Morocco. We produced for each individual Whole-Genome Sequence at 12X coverage and MeDIP-Seq data, to identify regions under selection and those differentially methylated. For both species, the analysis of genetic differences (FST) along the genome between animals from localities with high vs. low temperature annual variations detected candidate genes under selection in relation to environmental perception (5 genes), immunity (4 genes), reproduction (8 genes) and production (11 genes). Moreover, we found for each species one differentially methylated gene, namely AGPTA4 in goat and SLIT3 in sheep, which were both related, among other functions, to milk production and muscle development. In both sheep and goats, the comparison between genomic regions impacted by genetic and epigenetic variations suggests that climatic variations impacted similar biological pathways but different genes.

During domestication processes, changes in selective pressures induce multiple phenotypical, physiological, and behavioral changes in target species. The rise of next-generation sequencing has provided a chance to study the genetics bases of these changes, most of the time based on single nucleotide polymorphisms (SNPs). However, several studies have highlighted the impact of structural variations (SVs) on individual fitness, particularly in domestic species. We aimed at unraveling the role of SVs during the domestication and later improvement of small ruminants by analyzing whole-genome sequences of 40 domestic sheep and 11 of their close wild relatives (Ovis orientalis), and 40 goats and 18 of their close wild relatives (Capra aegagrus). Using a combination of detection tools, we called 45,796 SVs in Ovis and 15,047 SVs in Capra genomes, including insertions, deletions, inversions, copy number variations, and chromosomal translocations. Most of these SVs were previously unreported in small ruminants. 69 and 45 SVs in sheep and goats, respectively, were in genomic regions with neighboring SNPs highly differentiated between wilds and domestics (i.e., putatively related to domestication). Among them, 25 and 20 SVs were close to or overlapping with genes related to physiological and morpho-anatomical traits linked with productivity (e.g., size, meat or milk quality, wool color), reproduction, or immunity. Finally, several of the SVs differentiated between wilds and domestics would not have been detected by screening only the differentiation of SNPs surrounding them, highlighting the complementarity of SVs and SNPs based approaches to detect signatures of selection.

Whole genome sequences (WGS) greatly increase our ability to precisely infer population genetic parameters, demographic processes, and selection signatures. However, WGS may still be not affordable for a representative number of individuals/populations. In this context, our goal was to assess the efficiency of several SNP genotyping strategies by testing their ability to accurately estimate parameters describing neutral diversity and to detect signatures of selection. We analysed 110 WGS at 12× coverage for four different species, i.e., sheep, goats and their wild counterparts. From these data we generated 946 data sets corresponding to random panels of 1K to 5M variants, commercial SNP chips and exome capture, for sample sizes of five to 48 individuals. We also extracted low-coverage genome resequencing of 1×, 2× and 5× by randomly subsampling reads from the 12× resequencing data. Globally, 5K to 10K random variants were enough for an accurate estimation of genome diversity. Conversely, commercial panels and exome capture displayed strong ascertainment biases. Besides the characterization of neutral diversity, the detection of the signature of selection and the accurate estimation of linkage disequilibrium (LD) required high-density panels of at least 1M variants. Finally, genotype likelihoods increased the quality of variant calling from low coverage resequencing but proportions of incorrect genotypes remained substantial, especially for heterozygote sites. Whole genome resequencing coverage of at least 5× appeared to be necessary for accurate assessment of genomic variations. These results have implications for studies seeking to deploy low-density SNP collections or genome scans across genetically diverse populations/species showing similar genetic characteristics and patterns of LD decay for a wide variety of purposes.

Sheep, the Jaagsiekte sheep retrovirus (JSRV) and its endogenous forms (enJSRVs) are a good model to study long-time relationships between retroviruses and their hosts. Taking advantage of 76 whole genome resequencing data of wild and domestic Ovis, we investigated the evolution of this relationship. An innovative analysis of re-sequencing data allowed characterizing 462 enJSRVs insertion sites (including 435 newly described insertions) in the Ovis genus. We focused our study on endogenous copies inserted in the q13 locus of chromosome 6 (6q13). Those copies are known to confer resistance against exogenous JSRV thanks to alleles bearing a mutation in the gag gene. We characterized (i) the distribution of protective and non-protective alleles across Ovis species and (ii) the copy number variation of the 6q13 locus. Our results challenged the previous hypothesis of fixation and amplification of the protective copies in relation with domestication, and allowed building a new model for the evolution of the 6q13 locus. JSRV would have integrated the 6q13 locus after the Ovis-Capra divergence (5-11 MYA) and before the Ovis diversification (2.4-5 MYA). The protective mutation in the enJSRV 6q13 copy appeared shortly after its insertion and was followed by genomic amplifications, after the divergence between Pachyform lineage on one side and the Argaliform and moufloniform lineages on the other (2.4-5 MYA). Considering the potential selective advantage of the protective mutation, its fixation in both sheep and its closest wild relative Ovis orientalis may be due to natural selection before domestication from O. orientalis populations.

International audience; North Africa represents a rich and early reservoir of goat genetic diversity, from which the main African breeds have been derived. In this study, the genetic diversity of four indigenous Algerian goat breeds (i.e., Arabia, Makatia, M'Zabite and Kabyle, with n = 12 for each breed) has been investigated for the first time by genome-wide SNP genotyping; moreover in a broader context, genetic structuration of Algerian and Moroccan goats was explored (via FST, MDS, STRUCTURE, FineSTRUCTURE, BAPS, sPCA and DAPC analyses). At national level, the study revealed high level of genetic diversity and a significant phenomenon of admixture affecting all the Algerian breeds. At broader scale, clear global genetic homogeneity appeared considering both Algerian and Moroccan stocks. Indeed, genetic structuration was almost nonexistent among Arabia (from Algeria), Draa, Black and Nord (from Morocco), while the ancestral Kabyle and M'Zabite breeds, reared by Berber peoples, showed genetic distinctness. The study highlighted the threat to the Maghrebin stock, probably induced by unsupervised cross-breeding practices which have intensified in recent centuries. Moreover, it underlined the necessity to deepen our understanding of the genetic resources represented by the resilient North-African goat stock.

DNA barcoding approaches are used to describe biodiversity by analysing specimens or environmental samples in taxonomic, phylogenetic and ecological studies. While sharing data among these disciplines would be highly valuable, this remains difficult because of contradictory requirements. The properties making a DNA barcode efficient for specimen identification or species delimitation are hardly reconcilable with those required for a powerful analysis of degraded DNA from environmental samples. The use of next generation sequencing methods open up the way towards the development of new markers (e.g., multilocus barcodes) that would overcome such limitations. However, several challenges should be taken up for coordinating actions at the interface between taxonomy, ecology, molecular biology and bioinformatics in order to develop methods and protocols compatible with both taxonomic and ecological studies.

The Podismini are melanopline grasshoppers with a Holarctic distribution and well represented in the Eurasian fauna. To investigate their controversial taxonomy and evolutionary history, we studied 86%, 78% and 33% respectively of the Eurasian, European and Asian Palaearctic genera (Otte, 1995; Eades et al., 2013). We reconstructed parsimony, maximum likelihood and Bayesian phylogenies using fragments of four genes (ITS1, 16S, 12S, CO2). We applied a Bayesian molecular clock to estimate the times of species divergence, and the event-based parsimony method to depict the biogeographic framework of the diversification. Our results suggest that the selected Eurasian Podismini constitute a monophyletic group inside the Melanoplinae, provided it includes the North American genus Phaulotettix. The clades proposed by the present study inside the Podismini do not fit the older morphological or cytological classifications, but are in agreement with more recent proposals. Furthermore, our results can be explained by a plausible biogeographic history in which the present geographical distribution of the Eurasian Podismini resulted from known changes, to the Cenozoic climate and vegetation, induced by major geological events including the genesis of high mountain chains (e.g., Himalayas, Altay, Alps) and large deserts (e.g., Gobi, Karakoum, Taklamakan), and the opening of marginal seas (e.g., Bering, Japanese and Yellow Seas).

Ecosystems across the globe are threatened by climate change and human activities. New rapid survey approaches for monitoring biodiversity would greatly advance assessment and understanding of these threats. Taking advantage of next-generation DNA sequencing, we tested an approach we call metabarcoding: high-throughput and simultaneous taxa identification based on a very short (usually

Virtually all empirical ecological studies require species identification during data collection. DNA metabarcoding refers to the automated identification of multiple species from a single bulk sample containing entire organisms or from a single environmental sample containing degraded DNA (soil, water, faeces, etc.). It can be implemented for both modern and ancient environmental samples. The availability of next-generation sequencing platforms and the ecologists' need for high-throughput taxon identification have facilitated the emergence of DNA metabarcoding. The potential power of DNA metabarcoding as it is implemented today is limited mainly by its dependency on PCR and by the considerable investment needed to build comprehensive taxonomic reference libraries. Further developments associated with the impressive progress in DNA sequencing will eliminate the currently required DNA amplification step, and comprehensive taxonomic reference libraries composed of whole organellar genomes and repetitive ribosomal nuclear DNA can be built based on the well-curated DNA extract collections maintained by standardized barcoding initiatives. The near-term future of DNA metabarcoding has an enormous potential to boost data acquisition in biodiversity research.

DNA metabarcoding refers to the DNA-based identification of multiple species from a single complex and degraded environmental sample. We developed new sampling and extraction protocols suitable for DNA metabarcoding analyses targeting soil extracellular DNA. The proposed sampling protocol has been designed to reduce, as much as possible, the influence of local heterogeneity by processing a large amount of soil resulting from the mixing of many different cores. The DNA extraction is based on the use of saturated phosphate buffer. The sampling and extraction protocols were validated first by analysing plant DNA from a set of 12 plots corresponding to four plant communities in alpine meadows, and, second, by conducting pilot experiments on fungi and earthworms. The results of the validation experiments clearly demonstrated that sound biological information can be retrieved when following these sampling and extraction procedures. Such a protocol can be implemented at any time of the year without any preliminary knowledge of specific types of organisms during the sampling. It offers the opportunity to analyse all groups of organisms using a single sampling/extraction procedure and opens the possibility to fully standardize biodiversity surveys.

Background The current extensive use of the domestic goat (Capra hircus) is the result of its medium size and high adaptability as multiple breeds. The extent to which its genetic variability was influenced by early domestication practices is largely unknown. A common standard by which to analyze maternally-inherited variability of livestock species is through complete sequencing of the entire mitogenome (mitochondrial DNA, mtDNA). Results We present the first extensive survey of goat mitogenomic variability based on 84 complete sequences selected from an initial collection of 758 samples that represent 60 different breeds of C. hircus, as well as its wild sister species, bezoar (Capra aegagrus) from Iran. Our phylogenetic analyses dated the most recent common ancestor of C. hircus to ~460,000 years (ka) ago and identified five distinctive domestic haplogroups (A, B1, C1a, D1 and G). More than 90 % of goats examined were in haplogroup A. These domestic lineages are predominantly nested within C. aegagrus branches, diverged concomitantly at the interface between the Epipaleolithic and early Neolithic periods, and underwent a dramatic expansion starting from ~12–10 ka ago. Conclusions Domestic goat mitogenomes descended from a small number of founding haplotypes that underwent domestication after surviving the last glacial maximum in the Near Eastern refuges. All modern haplotypes A probably descended from a single (or at most a few closely related) female C. aegagrus. Zooarchaelogical data indicate that domestication first occurred in Southeastern Anatolia. Goats accompanying the first Neolithic migration waves into the Mediterranean were already characterized by two ancestral A and C variants. The ancient separation of the C branch (~130 ka ago) suggests a genetically distinct population that could have been involved in a second event of domestication. The novel diagnostic mutational motifs defined here, which distinguish wild and domestic haplogroups, could be used to understand phylogenetic relationships among modern breeds and ancient remains and to evaluate whether selection differentially affected mitochondrial genome variants during the development of economically important breeds. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2342-2) contains supplementary material, which is available to authorized users.

The analysis of food webs and their dynamics facilitates understanding of the mechanistic processes behind community ecology and ecosystem functions. Having accurate techniques for determining dietary ranges and components is critical for this endeavour. While visual analyses and early molecular approaches are highly labour intensive and often lack resolution, recent DNA-based approaches potentially provide more accurate methods for dietary studies. A suite of approaches have been used based on the identification of consumed species by characterization of DNA present in gut or faecal samples. In one approach, a standardized DNA region (DNA barcode) is PCR amplified, amplicons are sequenced and then compared to a reference database for identification. Initially, this involved sequencing clones from PCR products, and studies were limited in scale because of the costs and effort required. The recent development of next generation sequencing (NGS) has made this approach much more powerful, by allowing the direct characterization of dozens of samples with several thousand sequences per PCR product, and has the potential to reveal many consumed species simultaneously (DNA metabarcoding). Continual improvement of NGS technologies, on-going decreases in costs and current massive expansion of reference databases make this approach promising. Here we review the power and pitfalls of NGS diet methods. We present the critical factors to take into account when choosing or designing a suitable barcode. Then, we consider both technical and analytical aspects of NGS diet studies. Finally, we discuss the validation of data accuracy including the viability of producing quantitative data.

The present study characterizes the mitochondrial DNA diversity of Moroccan goats. 150 goats of different phenotypic entities were sampled over four geographic regions covering most of the Moroccan territory and the HVI segment of their mitochondrial DNA (mtDNA) control region was sequenced. The 150 Moroccan goats represented 97 haplotypes for this mtDNA fragment. Most of this diversity was present within phenotypic entities and within geographic regions. This weak genetic structure may result from the fact that all haplotypes were already mixed in the populations that colonized Morocco and/or the existence of recurrent gene flows from Mediterranean routes. Comparing the Moroccan haplotype diversity to that of 21 haplotypes representative of the worldwide diversity showed that all the Moroccan goats studied belonged to the A haplogroup that is preponderant in the world. The haplotypes of the Northern region appeared to be less diverse, what would probably reflect a stronger founder effect in this region.

New insights for the systematic and evolution of the wild sheep are provided by molecular phylogenies inferred from Maximum parsimony, Bayesian, Maximum likelihood, and Neighbor-Joining methods. The phylogeny of the wild sheep was based on cytochrome b sequences of 290 samples representative of most of the sub-species described in the genus Ovis. The result was confirmed by a combined tree based on cytochrome b and nuclear sequences for 79 Ovis samples representative of the robust clades established with mitochondrial data. Urial and mouflon, which are either considered as a single or two separate species, form two monophyletic groups (O. orientalis and O. vignei). Their hybrids appear in one or the other group, independently from their geographic origin. The European mouflon O. musimon is clearly in the O. orientalis clade. The others species, O. dalli, O. canadensis, O. nivicola, and O. ammon are monophyletic. The results support an Asiatic origin of the genus Ovis, followed by a migration to North America through North-Eastern Asia and the Bering Strait and a diversification of the genus in Eurasia less than 3 million years ago. Our results show that the evolution of the genus Ovis is a striking example of successive speciation events occurring along the migration routes propagating from the ancestral area.

The development of DNA barcoding (species identification using a standardized DNA sequence), and the availability of recent DNA sequencing techniques offer new possibilities in diet analysis. DNA fragments shorter than 100-150 bp remain in a much higher proportion in degraded DNA samples and can be recovered from faeces. As a consequence, by using universal primers that amplify a very short but informative DNA fragment, it is possible to reliably identify the plant taxon that has been eaten. According to our experience and using this identification system, about 50% of the taxa can be identified to species using the trnL approach, that is, using the P6 loop of the chloroplast trnL (UAA) intron. We demonstrated that this new method is fast, simple to implement, and very robust. It can be applied for diet analyses of a wide range of phytophagous species at large scales. We also demonstrated that our approach is efficient for mammals, birds, insects and molluscs. This method opens new perspectives in ecology, not only by allowing large-scale studies on diet, but also by enhancing studies on resource partitioning among competing species, and describing food webs in ecosystems.

International audience; For about 10 000 years, farmers have been managing cattle, sheep, and goats in a sustainable way, leading to animals that are well adapted to the local conditions. About 200 years ago, the situation started to change dramatically, with the rise of the concept of breed. All animals from the same breed began to be selected for the same phenotypic characteristics, and reproduction among breeds was seriously reduced. This corresponded to a strong fragmentation of the initial populations. A few decades ago, the selection pressures were increased again in order to further improve productivity, without enough emphasis on the preservation of the overall genetic diversity. The efficiency of modern selection methods successfully increased the production, but with a dramatic loss of genetic variability. Many industrial breeds now suffer from inbreeding, with effective population sizes falling below 50. With the development of these industrial breeds came economic pressure on farmers to abandon their traditional breeds, and many of these have recently become extinct as a result. This means that genetic resources in cattle, sheep, and goats are highly endangered, particularly in developed countries. It is therefore important to take measures that promote a sustainable management of these genetic resources; first, by in situ preservation of endangered breeds; second, by using selection programmes to restore the genetic diversity of industrial breeds; and finally, by protecting the wild relatives that might provide useful genetic resources.

Genotyping errors occur when the genotype determined after molecular analysis does not correspond to the real genotype of the individual under consideration. Virtually every genetic data set includes some erroneous genotypes, but genotyping errors remain a taboo subject in population genetics, even though they might greatly bias the final conclusions, especially for studies based on individual identification. Here, we consider four case studies representing a large variety of population genetics investigations differing in their sampling strategies (noninvasive or traditional), in the type of organism studied (plant or animal) and the molecular markers used [microsatellites or amplified fragment length polymorphisms (AFLPs)]. In these data sets, the estimated genotyping error rate ranges from 0.8% for microsatellite loci from bear tissues to 2.6% for AFLP loci from dwarf birch leaves. Main sources of errors were allelic dropouts for microsatellites and differences in peak intensities for AFLPs, but in both cases human factors were non-negligible error generators. Therefore, tracking genotyping errors and identifying their causes are necessary to clean up the data sets and validate the final results according to the precision required. In addition, we propose the outline of a protocol designed to limit and quantify genotyping errors at each step of the genotyping process. In particular, we recommend (i) several efficient precautions to prevent contaminations and technical artefacts; (ii) systematic use of blind samples and automation; (iii) experience and rigor for laboratory work and scoring; and (iv) systematic reporting of the error rate in population genetics studies.

Environmental DNA (eDNA) metabarcoding is increasingly used to study the present and past biodiversity. eDNA analyses often rely on amplification of very small quantities or degraded DNA. To avoid missing detection of taxa that are actually present (false negatives), multiple extractions and amplifications of the same samples are often performed. However, the level of replication needed for reliable estimates of the presence/absence patterns remains an unaddressed topic. Furthermore, degraded DNA and PCR/sequencing errors might produce false positives. We used simulations and empirical data to evaluate the level of replication required for accurate detection of targeted taxa in different contexts and to assess the performance of methods used to reduce the risk of false detections. Furthermore, we evaluated whether statistical approaches developed to estimate occupancy in the presence of observational errors can successfully estimate true prevalence, detection probability and false-positive rates. Replications reduced the rate of false negatives; the optimal level of replication was strongly dependent on the detection probability of taxa. Occupancy models successfully estimated true prevalence, detection probability and false-positive rates, but their performance increased with the number of replicates. At least eight PCR replicates should be performed if detection probability is not high, such as in ancient DNA studies. Multiple DNA extractions from the same sample yielded consistent results; in some cases, collecting multiple samples from the same locality allowed detecting more species. The optimal level of replication for accurate species detection strongly varies among studies and could be explicitly estimated to improve the reliability of results.

International audience; Roads are of major concern in conservation biology, as they are known to restrict animal movements through landscape fragmentation, and may therefore impact genetic patterns in native terrestrial organisms. We assessed the effect of two large-scale transportation infrastructures (LTIs), a 40-year-old highway and a 30 year-old high-speed railway, on the spatial genetic structure of the alpine newt Ichthyosaura alpestris, a highly nomadic amphibian. Genetic data were gathered following a targeted individual-based sampling scheme and analysed using both overlay and correlative methods. While simulations suggested that the highway may be old enough for a significant barrier effect to be detected, LTIs were never detected as barriers to gene flow: inferred genetic boundaries rather coincided with transition zones between major landscape entities. Furthermore, spatial principal component analysis, a method designed to reveal cryptic genetic spatial patterns in high gene flow species, counter-intuitively suggested that the highway may act as a potential dispersal corridor in low-quality habitats, thus challenging traditional hypotheses on road impacts in amphibians. Our study showed that considering local interactions between species, infrastructures and landscape-specific characteristics is essential for better understanding the potential impacts of roads on movement patterns in terrestrial organisms.

International audience; Genetic data are increasingly used in landscape ecology for the indirect assessment of functional connectivity, that is, the permeability of landscape to movements of organ- isms. Among available tools, matrix correlation analyses (e.g. Mantel tests or mixed models) are commonly used to test for the relationship between pairwise genetic dis- tances and movement costs incurred by dispersing individuals. When organisms are spatially clustered, a population-based sampling scheme (PSS) is usually performed, so that a large number of genotypes can be used to compute pairwise genetic distances on the basis of allelic frequencies. Because of financial constraints, this kind of sam- pling scheme implies a drastic reduction in the number of sampled aggregates, thereby reducing sampling coverage at the landscape level. We used matrix correlation analy- ses on simulated and empirical genetic data sets to investigate the efficiency of an individual-based sampling scheme (ISS) in detecting isolation-by-distance and isola- tion-by-barrier patterns. Provided that pseudo-replication issues are taken into account (e.g. through restricted permutations in Mantel tests), we showed that the use of inter- individual measures of genotypic dissimilarity may efficiently replace interpopulation measures of genetic differentiation: the sampling of only three or four individuals per aggregate may be sufficient to efficiently detect specific genetic patterns in most situa- tions. The ISS proved to be a promising methodological alternative to the more con- ventional PSS, offering much flexibility in the spatial design of sampling schemes and ensuring an optimal representativeness of landscape heterogeneity in data, with few aggregates left unsampled. Each strategy offering specific advantages, a combined use of both sampling schemes is discussed.

International audience; The goat (Capra hircus) is one of the earliest domesticated species ca. 10,500 years ago in the Middle-East where its wild ancestor, the bezoar (Capra aegagrus), still occurs. During the Neolithic dispersal, the domestic goat was then introduced in Europe, including the main Mediterranean islands. Islands are interesting models as they maintain traces of ancient colonization, historical exchanges or of peculiar systems of husbandry. Here, we compare the mitochondrial genetic diversity of both medieval and extant goats in the Island of Corsica that presents an original and ancient model of breeding with free-ranging animals. We amplified a fragment of the Control Region for 21 medieval and 28 current goats. Most of them belonged to the A haplogroup, the most worldwide spread and frequent today, but the C haplogroup is also detected at low frequency in the current population. Present Corsican goats appeared more similar to medieval goats than to other European goat populations. Moreover, 16 out of the 26 haplotypes observed were endemic to Corsica and the inferred demographic history suggests that the population has remained constant since the Middle Ages. Implications of these results on management and conservation of endangered Corsican goats currently decimated by a disease are addressed.

International audience; This study introduces a novel DNA sampling method in amphibians using skin swabs. We assessed the relevancy of skin swabs relevancy for genetic studies by amplifying a set of 17 microsatellite markers in the alpine newt Ichthyosaura alpestris, including 14 new polymorphic loci, and a set of 11 microsatellite markers in Hyla arborea, from DNA collected with buccal swabs (the standard swab method), dorsal skin swabs and ventral skin swabs. We tested for quality and quantity of collected DNA with each method by comparing electrophoresis migration patterns. The consistency between genotypes obtained from skin swabs and buccal swabs was assessed. Dorsal swabs performed better than ventral swabs in both species, possibly due to differences in skin structure. Skin swabbing proved to be a useful alternative to buccal swabbing for small or vulnerable animals: by drastically limiting handling, this method may improve the trade-off between the scientific value of collected data, individual welfare and species conservation. In addition, the 14 new polymorphic microsatellites for the alpine newt will increase the power of genetic studies in this species. In four populations from France (n = 19-25), the number of alleles per locus varied from 2 to 16 and expected heterozygosities ranged from 0.04 to 0.91. Presence of null alleles was detected in two markers and two pairs displayed gametic disequilibrium. No locus appeared to be sex-linked.

International audience; Using non-conventional markers, DNA metabarcoding allows biodiversity assessment from complex substrates. In this article, we present ecoPrimers, a software for identifying new barcode markers and their associated PCR primers. ecoPrimers scans whole genomes to find such markers without a priori knowledge. ecoPrimers optimizes two quality indices measuring taxonomical range and discrimination to select the most efficient markers from a set of reference sequences, according to specific experimental constraints such as marker length or specifically targeted taxa. The key step of the algorithm is the identification of conserved regions among reference sequences for anchoring primers. We propose an efficient algorithm based on data mining, that allows the analysis of huge sets of sequences. We evaluate the efficiency of ecoPrimers by running it on three different sequence sets: mitochondrial, chloroplast and bacterial genomes. Identified barcode markers correspond either to barcode regions already in use for plants or animals, or to new potential barcodes. Results from empirical experiments carried out on a promising new barcode for analyzing vertebrate diversity fully agree with expectations based on bioinformatics analysis. These tests demonstrate the efficiency of ecoPrimers for inferring new barcodes fitting with diverse experimental contexts. ecoPrimers is available as an open source project at: http://www.grenoble.prabi.fr/trac/ecoPrimers.

The Melanoplinae is one of the largest subfamilies of the Acrididae grasshoppers, with a Holarctic–Neotropical distribution. To investigate its controversial taxonomy and evolutionary history, we studied 19 species representative of its main tribes, and 7 Acridoidea outgroup species. More than 1650 base pairs of three regions of nuclear rDNA (18S, ITS1, 28S) and one mitochondrial rDNA (12S) were combined and used to construct parsimony, maximum likelihood and Bayesian phylogenies. Our results correspond with the present geographical distribution of the taxa rather than the existing taxonomy based on morphological characters. The morphologically unclassified and atypical taxa group with the Neotropical Melanoplinae. The currently recognized Melanoplinae appear to be polyphyletic due to inclusion of the Mexican genus Netrosoma. The distribution of the American and Eurasiatic Melanoplinae fauna can be explained by climatic and geological events, such as the Andean uplift, that would have affected the diversification and migration of Neotropical taxa.

The emergence of farming during the Neolithic transition, including the domestication of livestock, was a critical point in the evolution of human kind. The goat ( Capra hircus ) was one of the first domesticated ungulates. In this study, we compared the genetic diversity of domestic goats to that of the modern representatives of their wild ancestor, the bezoar, by analyzing 473 samples collected over the whole distribution range of the latter species. This partly confirms and significantly clarifies the goat domestication scenario already proposed by archaeological evidence. All of the mitochondrial DNA haplogroups found in current domestic goats have also been found in the bezoar. The geographic distribution of these haplogroups in the wild ancestor allowed the localization of the main domestication centers. We found no haplotype that could have been domesticated in the eastern half of the Iranian Plateau, nor further to the east. A signature of population expansion in bezoars of the C haplogroup suggests an early domestication center on the Central Iranian Plateau (Yazd and Kerman Provinces) and in the Southern Zagros (Fars Province), possibly corresponding to the management of wild flocks. However, the contribution of this center to the current domestic goat population is rather low (1.4%). We also found a second domestication center covering a large area in Eastern Anatolia, and possibly in Northern and Central Zagros. This last domestication center is the likely origin of almost all domestic goats today. This finding is consistent with archaeological data identifying Eastern Anatolia as an important domestication center.

Times Cited: 4; International audience; DNA barcoding should provide rapid, accurate and automatable species identifications by using a standardized DNA region as a tag. Based on sequences available in GenBank and sequences produced for this study, we evaluated the resolution power of the whole chloroplast trnL (UAA) intron (254-767 bp) and of a shorter fragment of this intron (the P6 loop, 10-143 bp) amplified with highly conserved primers. The main limitation of the whole trnL intron for DNA barcoding remains its relatively low resolution (67.3% of the species from GenBank unambiguously identified). The resolution of the P6 loop is lower (19.5% identified) but remains higher than those of existing alternative systems. The resolution is much higher in specific contexts such as species originating from a single ecosystem, or commonly eaten plants. Despite the relatively low resolution, the whole trnL intron and its P6 loop have many advantages: the primers are highly conserved, and the amplification system is very robust. The P6 loop can even be amplified when using highly degraded DNA from processed food or from permafrost samples, and has the potential to be extensively used in food industry, in forensic science, in diet analyses based on feces and in ancient DNA studies.

Times Cited: 24; International audience; Today, with the rapid development of population genomics, the genetic basis of adaptation can be unraveled directly at the genome level, without any prerequisites about the selectively advantageous genes or traits. For nonmodel species, it is now possible to screen many markers randomly scattered across the genome and to distinguish between the neutral genetic background and outlier loci displaying an atypical behavior (e.g., a higher differentiation between populations). This study investigated the genetic frame of adaptation to a gradient of altitude in the common frog (Rana temporaria) by means of a genome scan based on 392 amplified fragment length polymorphism markers. Using two outlier detection methods never applied to dominant data so far, we sought for loci with a genetic differentiation diverging from neutral expectations when comparing populations from different altitudes. All the detected loci were sorted out according to their most probable cause for outlier behavior and classified as false positives, outliers due to local effects, or outliers associated with altitude. Altogether, eight good candidate loci were identified as potentially involved in adaptation to altitude because they were picked out in several independent interaltitude comparisons. This result illustrated the potential of genome-wide surveys to reveal selection signatures along selection gradients, where the association between environmental variables and fitness-related traits may be complex and/or cryptic. In this article, we also underlined the need for confirmation of the selection footprints for the outlier loci. Finally, we provided some preliminary insights into the genetic basis of adaptation along an altitudinal cline in the common frog.

Times Cited: 66; International audience; Although genotyping errors affect most data and can markedly influence the biological conclusions of a study, they are too often neglected. Errors have various causes, but their occurrence and effect can be limited by considering these causes in the production and analysis of the data. Procedures that have been developed for dealing with errors in linkage studies, forensic analyses and non-invasive genotyping should be applied more broadly to any genetic study. We propose a protocol for estimating error rates and recommend that these measures be systemically reported to attest the reliability of published genotyping studies.

Times Cited: 4; International audience; The amplified fragment length polymorphism (AFLP) technique is one of the most informative and cost-effective fingerprinting methods. It produces polymerase chain reaction (PCR)-based multi-locus genotypes helpful in many areas of population genetics. This chapter focuses on technical laboratory information to successfully develop the AFLP technique for vertebrates. Several AFLP protocols are described, as well as recommendations about important factors of the procedure such as the choice of enzyme and primer combinations, the choice and scoring of markers, the influence of the genome size on the AFLP procedure, and the control and estimation of genotyping errors. Finally, this chapter proposes a troubleshooting guide to help resolve the main technical difficulties encountered during the AFLP procedure.

International audience; Phylogenetic relationships among the Chiastocheta species (Diptera, Anthomyiidae) were investigated using a 1320-bp mitochondrial DNA fragment including parts of cytochrome oxidase subunits I and II. Larvae of the Chiastocheta genus feed exclusively on Trollius (Ranunculaceae) seeds. Six Chiastocheta species (and two vicariant taxa) coexist on Trollius europaeus in Europe, where they were shown to be the sole pollinator of their host plant. In Asia, several Trollius species are parasitized by Chiastocheta spp., and several other Trollius species are free of Chiastocheta. The plant-insect association is thus obligate in Europe but facultative in Asia. This system therefore provides a unique opportunity to study the evolution of an obligate mutualism between a plant and its seed parasite pollinators and its consequences on insect diversification. Most parsimonious, distance, and maximum likelihood analyses of 16 haplotypes from Europe and 14 from Asia show that European species do not form a monophyletic group. The species relationship indicated by egg morphology is only partly supported by molecular data. Moreover, a strong discrepancy between nucleotide variation and male genitalia morphology variation was found, particularly within and among European species coexisting on a single host plant. Molecular divergence across species ranges from 0 to 4.6%, suggesting a recent origin of the genus and several events of diversification, involving both host shifts and within-host-plant radiation. We discuss alternative evolutionary scenarios that are consistent with molecular data. (C) 2002 Elsevier Science (USA).

DNA metabarcoding enables efficient characterization of species composition in environmental DNA or bulk biodiversity samples, and this approach is making significant and unique contributions in the field of ecology. In metabarcoding of animals, the cytochrome c oxidase subunit I (COI) gene is frequently used as the marker of choice because no other genetic region can be found in taxonomically verified databases with sequences covering so many taxa. However, the accuracy of metabarcoding datasets is dependent on recovery of the targeted taxa using conserved amplification primers. We argue that COI does not contain suitably conserved regions for most amplicon-based metabarcoding applications. Marker selection deserves increased scrutiny and available marker choices should be broadened in order to maximize potential in this exciting field of research.

Among apicomplexan parasites, ticks are known vector of several species belongs to three protozoan genera (Babesia, Theileria and Hepatozoon). During their life cycle, tick-transmitted apicomplexan parasites alternate between asexual (in vertebrate host) and sexual (in ticks) developmental stages. The major constraint for the proper diagnostics of those pathogens is the high possibility of mix infection, both in ticks and vertebrate hosts, with several species or genotypes. The aim of this study was to apply the Next Generation Sequencing (NGS) as a method of choice for simultaneous determination of the full spectrum of apicomplexan pathogens in ticks and the mammalian hosts. Therefore, A pair of universal primers were designed to flank a 167 bp barcode region of the 18s rRNA gene of all Babesia, Theileria and Hepatozoon species. The new protocol was evaluated on DNA samples isolated from 195 dogs and 144 ticks (Rhipicephalus armatus and R. pulchellus) collected from Northern Kenya. In total 301 sample (89%) were positive for apicomplexan infections; ranging from single to multiple infection with one species or several species and/or genotypes in a single sample. The most abundant apicomplexan pathogens were Hepatozoon followed by Babesia and Theileria, respectively. Further, the result shows that the barcode region entails enough variability that allows identifying the pathogens up to the subspecies and genotypes level. The exact methodological and results detailed will be presented later. This work was supported by the project OP VK CZ.1.07/2.3.00/30.0014.

International audience; Since the time of their domestication, goats (Capra hircus) have evolved in a large variety of locally adapted populations in response to different human and environmental pressures. In the present era, many indigenous populations are threatened with extinction due to their substitution by cosmopolitan breeds, while they might represent highly valuable genomic resources. It is thus crucial to characterize the neutral and adaptive genetic diversity of indigenous populations. A fine characterization of whole genome variation in farm animals is now possible by using new sequencing technologies. We sequenced the complete genome at 12× coverage of 44 goats geographically representative of the three phenotypically distinct indigenous populations in Morocco. The study of mitochondrial genomes showed a high diversity exclusively restricted to the haplogroup A. The 44 nuclear genomes showed a very high diversity (24 million variants) associated with low linkage disequilibrium. The overall genetic diversity was weakly structured according to geography and phenotypes. When looking for signals of positive selection in each population we identified many candidate genes, several of which gave insights into the metabolic pathways or biological processes involved in the adaptation to local conditions (e.g., panting in warm/desert conditions). This study highlights the interest of WGS data to characterize livestock genomic diversity. It illustrates the valuable genetic richness present in indigenous populations that have to be sustainably managed and may represent valuable genetic resources for the long-term preservation of the species.

Background DNA barcoding is a key tool for assessing biodiversity in both taxonomic and environmental studies. Essential features of barcodes include their applicability to a wide spectrum of taxa and their ability to identify even closely related species. Several DNA regions have been proposed as barcodes and the region selected strongly influences the output of a study. However, formal comparisons between barcodes remained limited until now. Here we present a standard method for evaluating barcode quality, based on the use of a new bioinformatic tool that performs in silico PCR over large databases. We illustrate this approach by comparing the taxonomic coverage and the resolution of several DNA regions already proposed for the barcoding of vertebrates. To assess the relationship between in silico and in vitro PCR, we also developed specific primers amplifying different species of Felidae, and we tested them using both kinds of PCR Results Tests on specific primers confirmed the correspondence between in silico and in vitro PCR. Nevertheless, results of in silico and in vitro PCRs can be somehow different, also because tuning PCR conditions can increase the performance of primers with limited taxonomic coverage. The in silico evaluation of DNA barcodes showed a strong variation of taxonomic coverage (i.e., universality): barcodes based on highly degenerated primers and those corresponding to the conserved region of the Cyt-b showed the highest coverage. As expected, longer barcodes had a better resolution than shorter ones, which are however more convenient for ecological studies analysing environmental samples. Conclusions In silico PCR could be used to improve the performance of a study, by allowing the preliminary comparison of several DNA regions in order to identify the most appropriate barcode depending on the study aims.