Microbiota Diversity Within and Between the Tissues of Two Wild Interbreeding Species.

Understanding the role of microbiota as reproductive barriers or sources of adaptive novelty in the fundamental biological phenomenon of speciation is an exciting new challenge necessitating exploration of microbiota variation in wild interbreeding species. We focused on two interbreeding cyprinid species, Chondrostoma nasus and Parachondrostoma toxostoma, which have geographic distributions characterized by a mosaic of hybrid zones. We described microbiota diversity and composition in the three main teleost mucosal tissues, the skin, gills and gut, in the parental parapatric populations. We found that tissue type was the principal determinant of bacterial community composition. In particular, there was strong microbiota differentiation between external and internal tissues, with secondary discrimination between the two species. These findings suggest that specific environmental and genetic filters associated with each species have shaped the bacterial communities, potentially reflecting deterministic assemblages of bacteria. We defined the core microbiota common to both Chondrostoma species for each tissue, highlighting the occurrence of microbe-host genome interactions at this critical level for studies of the functional consequences of hybridization. Further investigations will explore to what extend these specific tissue-associated microbiota signatures could be profoundly altered in hybrids, with functional consequences for post-mating reproductive isolation in relation to environmental constraints.

Challenges and advances for transcriptome assembly in non-model species.

Analyses of high-throughput transcriptome sequences of non-model organisms are based on two main approaches: de novo assembly and genome-guided assembly using mapping to assign reads prior to assembly. Given the limits of mapping reads to a reference when it is highly divergent, as is frequently the case for non-model species, we evaluate whether using blastn would outperform mapping methods for read assignment in such situations (>15% divergence). We demonstrate its high performance by using simulated reads of lengths corresponding to those generated by the most common sequencing platforms, and over a realistic range of genetic divergence (0% to 30% divergence). Here we focus on gene identification and not on resolving the whole set of transcripts (i.e. the complete transcriptome). For simulated datasets, the transcriptome-guided assembly based on blastn recovers 94.8% of genes irrespective of read length at 0% divergence; however, assignment rate of reads is negatively correlated with both increasing divergence level and reducing read lengths. Nevertheless, we still observe 92.6% of recovered genes at 30% divergence irrespective of read length. This analysis also produces a categorization of genes relative to their assignment, and suggests guidelines for data processing prior to analyses of comparative transcriptomics and gene expression to minimize potential inferential bias associated with incorrect transcript assignment. We also compare the performances of de novo assembly alone vs in combination with a transcriptome-guided assembly based on blastn both via simulation and empirically, using data from a cyprinid fish species and from an oak species. For any simulated scenario, the transcriptome-guided assembly using blastn outperforms the de novo approach alone, including when the divergence level is beyond the reach of traditional mapping methods. Combining de novo assembly and a related reference transcriptome for read assignment also addresses the bias/error in contigs caused by the dependence on a related reference alone. Empirical data corroborate these findings when assembling transcriptomes from the two non-model organisms: Parachondrostoma toxostoma (fish) and Quercus pubescens (plant). For the fish species, out of the 31,944 genes known from D. rerio, the guided and de novo assemblies recover respectively 20,605 and 20,032 genes but the performance of the guided assembly approach is much higher for both the contiguity and completeness metrics. For the oak, out of the 29,971 genes known from Vitis vinifera, the transcriptome-guided and de novo assemblies display similar performance, but the new guided approach detects 16,326 genes where the de novo assembly only detects 9,385 genes.

Multi-scale biomarker evaluation of the toxicity of a commercial azo dye (Disperse Red 1) in an animal model, the freshwater cnidarian Hydra attenuata.

Acute (24 h, 48 h, 72 h) and chronic (7 days) tests have been performed to evaluate the effects of the commercial azo dye Disperse Red 1 (DR1) using various biomarkers in the freshwater invertebrate Hydra attenuata. Morphological changes have been selected to calculate ecotoxicological thresholds for sublethal and lethal DR1 concentrations. A multinomial logistic model showed that the probability of each morphological stage occurrence was function of concentration, time and interaction between both. Results of oxidative balance parameter measurements (72 h and 7 days) suggest that polyps set up defense mechanisms to limit lipid peroxidation caused by DR1. DR1 exposure at hormetic concentrations induces increase of asexual reproductive rates. This result suggests (1) an impact on the fitness-related phenotypical traits and (2) trade-offs between reproduction and maintenance to allow the population to survive harsher conditions. Changes in serotonin immuno-labeling in polyps showing alterations in feeding behavior suggest that chronic DR1 exposure impaired neuronal processes related to ingesting behavior in H. attenuata. This ecotoxicity study sheds light on the possible serotonin function in Hydra model and reports for the first time that serotonin could play a significant role in feeding behavior. This study used a multi-scale biomarker approach investigating biochemical, morphological, reproductive and behavioral endpoints in Hydra attenuata. This organism is proposed for a pertinent animal model to assess ecotoxicological impact of pollutant mixtures in freshwater environment.

When Anthropogenic River Disturbance Decreases Hybridisation between Non-Native and Endemic Cyprinids and Drives an Ecomorphological Displacement towards Juvenile State in Both Species.

Understanding the impact of non-native species on native species is a major challenge in molecular ecology, particularly for genetically compatible fish species. Invasions are generally difficult to study because their effects may be confused with those of environmental or human disturbances. Colonized ecosystems are differently impacted by human activities, resulting in diverse responses and interactions between native and non-native species. We studied the dynamics between two Cyprinids species (invasive Chondrostoma nasus and endemic Parachondrostoma toxostoma) and their hybrids in 16 populations (from allopatric to sympatric situations and from little to highly fragmented areas) corresponding to 2,256 specimens. Each specimen was assigned to a particular species or to a hybrid pool using molecular identification (cytochrome b and 41 microsatellites). We carried out an ecomorphological analysis based on size, age, body shape, and diet (gut vacuity and molecular fecal contents). Our results contradicted our initial assumptions on the pattern of invasion and the rate of introgression. There was no sign of underperformance for the endemic species in areas where hybridisation occurred. In the unfragmented zone, the introduced species was found mostly downstream, with body shapes similar to those in allopatric populations while both species were found to be more insectivorous than the reference populations. However, high level of hybridisation was detected, suggesting interactions between the two species during spawning and/or the existence of hybrid swarm. In the disturbed zone, introgression was less frequent and slender body shape was associated with diatomivorous behaviour, smaller size (juvenile characteristics) and greater gut vacuity. Results suggested that habitat degradation induced similar ecomorphological trait changes in the two species and their hybrids (i.e. a transition towards a pedomorphic state) where the invasive species is more affected than the native species. Therefore, this study reveals a diversity of relationships between two genetically compatible species and emphasizes constraints on the invasion process in disturbed areas.

Evolutionary history of Chaetognatha inferred from molecular and morphological data: a case study for body plan simplification.

BACKGROUND: Chaetognatha are a phylum of marine carnivorous animals which includes more than 130 extant species. The internal systematics of this group have been intensively debated since it was discovered in the 18(th) century. While they can be traced back to the earlier Cambrian, they are an extraordinarily homogeneous phylum at the morphological level - a fascinating characteristic that puzzled many a scientist who has tried to clarify their taxonomy. Recent studies which have attempted to reconstruct a phylogeny using molecular data have relied on single gene analyses and a somewhat restricted taxon sampling. Here, we present the first large scale phylogenetic study of Chaetognatha based on a combined analysis of nearly the complete ribosomal RNA (rRNA) genes. We use this analysis to infer the evolution of some morphological characters. This work includes 36 extant species, mainly obtained from Tara Oceans Expedition 2009/2012, that represent 16 genera and 6 of the 9 extant families. RESULTS: Cladistic and phenetic analysis of morphological characters, geometric morphometrics and molecular small subunit (SSU rRNA) and large subunit (LSU rRNA) ribosomal genes phylogenies provided new insights into the relationships and the evolutionary history of Chaetognatha. We propose the following clade structure for the phylum: (((Sagittidae, Krohnittidae), Spadellidae), (Eukrohniidae, Heterokrohniidae)), with the Pterosagittidae included in the Sagittidae. The clade (Sagittidae, Krohnittidae) constitutes the monophyletic order of Aphragmophora. Molecular analyses showed that the Phragmophora are paraphyletic. The Ctenodontina/Flabellodontina and Syngonata/Chorismogonata hypotheses are invalidated on the basis of both morphological and molecular data. This new phylogeny also includes resurrected and modified genera within Sagittidae. CONCLUSIONS: The distribution of some morphological characters traditionally used in systematics and for species diagnosis suggests that the diversity in Chaetognatha was produced through a process of mosaic evolution. Moreover, chaetognaths have mostly evolved by simplification of their body plan and their history shows numerous convergent events of losses and reversions. The main morphological novelty observed is the acquisition of a second pair of lateral fins in Sagittidae, which represents an adaptation to the holoplanktonic niche.

QDD version 3.1: a user-friendly computer program for microsatellite selection and primer design revisited: experimental validation of variables determining genotyping success rate.

Microsatellite marker development has been greatly simplified by the use of high-throughput sequencing followed by in silico microsatellite detection and primer design. However, the selection of markers designed by the existing pipelines depends either on arbitrary criteria, or older studies on PCR success. Based on wet laboratory experiments, we have identified the following factors that are most likely to influence genotyping success rate: alignment score between the primers and the amplicon; the distance between primers and microsatellites; the length of the PCR product; target region complexity and the number of reads underlying the sequence. The QDD pipeline has been modified to include these most pertinent factors in the output to help the selection of markers. Furthermore, new features are also included in the present version: (i) not only raw sequencing reads are accepted as input, but also contigs, allowing the analysis of assembled high-coverage data; (ii) input data can be both in fasta and fastq format to facilitate the use of Illumina and IonTorrent reads; (iii) A comparison to known transposable elements allows their detection; (iv) A contamination check can be carried out by BLASTing potential markers against the nucleotide (nt) database of NCBI; (v) QDD3 is now also available imbedded into a virtual machine making installation easier and operating system independent. It can be used both on command-line version as well as integrated into a Galaxy server, providing a user-friendly interface, as well as the possibility to utilize a large variety of NGS tools.

Non-homogeneous combination of two porous genomes induces complex body shape trajectories in cyprinid hybrids.

INTRODUCTION: Hybridization is a common phenomenon in fish and is considered to be a major source of diversification. Deciphering the remoulding of genomic regions and phenotypes in zones where hybrid specimens occur is of particular interest to elucidate the emergence of evolutionary novelties. This approach is particularly challenging because the first step of hybridization seems to be the most important stage in the emergence of hybrid lineages. However, the signal can be significantly altered after only a few generations. RESULTS: We studied 41 microsatellites and partial cytochrome b gene sequences in 970 specimens belonging to two fish species (Chondrostoma nasus and Parachondrostoma toxostoma) in allopatric/parapatric zones, hybrids between them in a natural sympatric zone: the Ardèche basin. We showed that the genomic architecture in hybrids presented pattern heterogeneity of selection for the different loci. Indeed, the upstream part of the river (Rosières and Labeaume) presented an overdominant fitness of heterozygotes (12.20%) corresponding to a genomic compatibility, and underselection was observed for 4.88%-7.32% of the loci tested indicating a genomic incompatibility. Moreover the upstream station (Rosières) presented a positive selection of invasive C. nasus homozygotes (17.07% to 21.95%) indicating that hybridization may increase the fitness of admixed individuals.We showed that hybrid morphology (body shape based on 21 landmarks) correlated with genomic dilution indicating a species fingerprint. However, we demonstrated that the hybrid morphology was not a linear modification between the two parental species but a trade-off between several correlated traits. CONCLUSIONS: Hybrid specimens present a mosaic of genomic combination, showing regions with genomic compatibility and others with genomic incompatibility between the two species. Positive selection (invasive advantage ranging from 9.76% to 21.95% of the loci) was evidenced in the upstream part of the Ardèche indicating that environmental selection makes a substantial contribution. Although the presence of a dam is known to impose heterogeneous hybrid zones between these two species, we demonstrated in this study that a natural environment can also generate a hybrid zone with a large number (and diversity) of hybrids. The combination of the two genomes in the hybrids results in complex ontogenetic trajectories (with different morphological traits evolving at different rates) that correspond to novel developmental pathways.

Phylogenetic analysis of Thecosomata Blainville, 1824 (holoplanktonic opisthobranchia) using morphological and molecular data.

Thecosomata is a marine zooplankton group, which played an important role in the carbonate cycle in oceans due to their shell composition. So far, there is important discrepancy between the previous morphological-based taxonomies, and subsequently the evolutionary history of Thecosomata. In this study, the remarkable planktonic sampling of TARA Oceans expedition associated with a set of various other missions allowed us to assess the phylogenetic relationships of Thecosomata using morphological and molecular data (28 S and COI genes). The two gene trees showed incongruities (e.g. Hyalocylis, Cavolinia), and high congruence between morphological and 28S trees (e.g. monophyly of Euthecosomata). The monophyly of straight shell species led us to reviving the Orthoconcha, and the split of Limacinidae led us to the revival of Embolus inflata replacing Limacina inflata. The results also jeopardized the Euthecosomata families that are based on plesiomorphic character state as in the case for Creseidae which was not a monophyletic group. Divergence times were also estimated, and suggested that the evolutionary history of Thecosomata was characterized by four major diversifying events. By bringing the knowledge of palaeontology, we propose a new evolutionary scenario for which macro-evolution implying morphological innovations were rhythmed by climatic changes and associated species turn-over that spread from the Eocene to Miocene, and were shaped principally by predation and shell buoyancy.

A shot in the genome: how accurately do shotgun 454 sequences represent a genome?

BACKGROUND: Next generation sequencing (NGS) provides a valuable method to quickly obtain sequence information from non-model organisms at a genomic scale. In principle, if sequencing is not targeted for a genomic region or sequence type (e.g. coding region, microsatellites) NGS reads can be used as a genome snapshot and provide information on the different types of sequences in the genome. However, no study has ascertained if a typical 454 dataset of low coverage (1/4-1/8 of a PicoTiter plate leading to generally less than 0.1x of coverage) represents all parts of genomes equally. FINDINGS: Partial genome shotgun sequencing of total DNA (without enrichment) on a 454 NGS platform was used to obtain reads of Apis mellifera (454 reads hereafter). These 454 reads were compared to the assembled chromosomes of this species in three different aspects: (i) dimer and trimer compositions, (ii) the distribution of mapped 454 sequences along the chromosomes and (iii) the numbers of different classes of microsatellites. Highly significant chi-square tests for all three types of analyses indicated that the 454 data is not a perfect random sample of the genome. Only the number of 454 reads mapped to each of the 16 chromosomes and the number of microsatellites pooled by motif (repeat unit) length was not significantly different from the expected values. However, a very strong correlation (correlation coefficients greater than 0.97) was observed between most of the 454 variables (the number of different dimers and trimers, the number of 454 reads mapped to each chromosome fragments of one Mb, the number of 454 reads mapped to each chromosome, the number of microsatellites of each class) and their corresponding genomic variables. CONCLUSIONS: The results of chi square tests suggest that 454 shotgun reads cannot be regarded as a perfect representation of the genome especially if the comparison is done on a finer scale (e.g. chromosome fragments instead of whole chromosomes). However, the high correlation between 454 and genome variables tested indicate that a high proportion of the variability of 454 variables is explained by their genomic counterparts. Therefore, we conclude that using 454 data to obtain information on the genome is biologically meaningful.

Multiple quaternary refugia in the eastern Guiana shield revealed by comparative phylogeography of 12 frog species.

The Guiana Shield (GS) is one of the most pristine regions of Amazonia and biologically one of the richest areas on Earth. How and when this massive diversity arose remains the subject of considerable debate. The prevailing hypothesis of Quaternary glacial refugia suggests that a part of the eastern GS, among other areas in Amazonia, served as stable forested refugia during periods of aridity. However, the recently proposed disturbance-vicariance hypothesis proposes that fluctuations in temperature on orbital timescales, with some associated aridity, have driven Neotropical diversification. The expectations of the temporal and spatial organization of biodiversity differ between these two hypotheses. Here, we compare the genetic structure of 12 leaf-litter inhabiting frog species from the GS lowlands using a combination of mitochondrial and nuclear sequences in an integrative analytical approach that includes phylogenetic reconstructions, molecular dating, and Geographic Information System methods. This comparative and integrated approach overcomes the well-known limitations of phylogeographic inference based on single species and single loci. All of the focal species exhibit distinct phylogeographic patterns highlighting taxon-specific historical distributions, ecological tolerances to climatic disturbance, and dispersal abilities. Nevertheless, all but one species exhibit a history of fragmentation/isolation within the eastern GS during the Quaternary with spatial and temporal concordance among species. The signature of isolation in northern French Guiana (FG) during the early Pleistocene is particularly clear. Approximate Bayesian Computation supports the synchrony of the divergence between northern FG and other GS lineages. Substructure observed throughout the GS suggests further Quaternary fragmentation and a role for rivers. Our findings support fragmentation of moist tropical forest in the eastern GS during this period when the refuge hypothesis would have the region serving as a contiguous wet-forest refuge.

High-throughput microsatellite isolation through 454 GS-FLX Titanium pyrosequencing of enriched DNA libraries.

Microsatellites (or SSRs: simple sequence repeats) are among the most frequently used DNA markers in many areas of research. The use of microsatellite markers is limited by the difficulties involved in their de novo isolation from species for which no genomic resources are available. We describe here a high-throughput method for isolating microsatellite markers based on coupling multiplex microsatellite enrichment and next-generation sequencing on 454 GS-FLX Titanium platforms. The procedure was calibrated on a model species (Apis mellifera) and validated on 13 other species from various taxonomic groups (animals, plants and fungi), including taxa for which severe difficulties were previously encountered using traditional methods. We obtained from 11,497 to 34,483 sequences depending on the species and the number of detected microsatellite loci ranged from 199 to 5791. We thus demonstrated that this procedure can be readily and successfully applied to a large variety of taxonomic groups, at much lower cost than would have been possible with traditional protocols. This method is expected to speed up the acquisition of high-quality genetic markers for nonmodel organisms.

Accuracy and quality assessment of 454 GS-FLX Titanium pyrosequencing.

BACKGROUND: The rapid evolution of 454 GS-FLX sequencing technology has not been accompanied by a reassessment of the quality and accuracy of the sequences obtained. Current strategies for decision-making and error-correction are based on an initial analysis by Huse et al. in 2007, for the older GS20 system based on experimental sequences. We analyze here the quality of 454 sequencing data and identify factors playing a role in sequencing error, through the use of an extensive dataset for Roche control DNA fragments. RESULTS: We obtained a mean error rate for 454 sequences of 1.07%. More importantly, the error rate is not randomly distributed; it occasionally rose to more than 50% in certain positions, and its distribution was linked to several experimental variables. The main factors related to error are the presence of homopolymers, position in the sequence, size of the sequence and spatial localization in PT plates for insertion and deletion errors. These factors can be described by considering seven variables. No single variable can account for the error rate distribution, but most of the variation is explained by the combination of all seven variables. CONCLUSIONS: The pattern identified here calls for the use of internal controls and error-correcting base callers, to correct for errors, when available (e.g. when sequencing amplicons). For shotgun libraries, the use of both sequencing primers and deep coverage, combined with the use of random sequencing primer sites should partly compensate for even high error rates, although it may prove more difficult than previous thought to distinguish between low-frequency alleles and errors.

SESAME (SEquence Sorter & AMplicon Explorer): genotyping based on high-throughput multiplex amplicon sequencing.

SUMMARY: Characterizing genetic diversity through genotyping short amplicons is central to evolutionary biology. Next-generation sequencing (NGS) technologies changed the scale at which these type of data are acquired. SESAME is a web application package that assists genotyping of multiplexed individuals for several markers based on NGS amplicon sequencing. It automatically assigns reads to loci and individuals, corrects reads if standard samples are available and provides an intuitive graphical user interface (GUI) for allele validation based on the sequences and associated decision-making tools. The aim of SESAME is to help allele identification among a large number of sequences. AVAILABILITY: SESAME and its documentation are freely available under the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported Licence for Windows and Linux from http://www1.montpellier.inra.fr/CBGP/NGS/ or http://tinyurl.com/ngs-sesame.

Representativeness of microsatellite distributions in genomes, as revealed by 454 GS-FLX titanium pyrosequencing.

BACKGROUND: Microsatellites are markers of choice in population genetics and genomics, as they provide useful insight into patterns and processes as diverse as genome evolutionary dynamics and demographic processes. The acquisition of microsatellites through multiplex-enriched libraries and 454 GS-FLX Titanium pyrosequencing is a promising new tool for the isolation of new markers in unknown genomes. This approach can also be used to evaluate the extent to which microsatellite-enriched libraries are representative of the genome from which they were isolated. In this study, we deciphered potential discrepancies in microsatellite content recovery for two reference genomes (Apis mellifera and Danio rerio), selected on the basis of their extreme heterogeneity in terms of the proportions and distributions of microsatellites on chromosomes. RESULTS: The A. mellifera genome, in particular, was found to be highly heterogeneous, due to extremely high rates of recombination, with hotspots, but the only bias consistently introduced into pyrosequenced multiplex-enriched libraries concerned sequence length, with the overrepresentation of sequences 160 to 320 bp in length. Other deviations from expected proportions or distributions of motifs on chromosomes were observed, but the significance and intensity of these deviations was mostly limited. Furthermore, no consistent adverse competition between multiplexed probes was observed during the motif enrichment phase. CONCLUSIONS: This approach therefore appears to be a promising strategy for improving the development of microsatellites, as it introduces no major bias in terms of the proportions and distribution of microsatellites.

A PCR-based method for diet analysis in freshwater organisms using 18S rDNA barcoding on faeces.

The development of DNA barcoding from faeces represents a promising method for animal diet analysis. However, current studies mainly rely on prior knowledge of prey diversity for a specific predator rather than on a range of its potential prey species. Considering that the feeding behaviour of teleosts may evolve with their environment, it could prove difficult to establish an exhaustive listing of their prey. In this article, we extend the DNA barcoding approach to diet analysis to allow the inclusion of a wide taxonomic range of potential prey items. Thirty-four ecological clade-specific primer sets were designed to cover a large proportion of prey species found in European river ecosystems. Selected primers sets were tested on isolated animal, algal or plant tissues and thereafter on fish faeces using nested PCR to increase DNA detection sensitivity. The PCR products were sequenced and analysed to confirm the identity of the taxa and to validate the method. The methodology developed here was applied to a diet analysis of three freshwater cyprinid species that are assumed to have similar feeding behaviour [Chondrostoma toxostoma toxostoma (Vallot 1837), Chondrostoma nasus nasus (Linnaeus, 1758) and Barbus barbus, (Linneaus 1758)]. These three species were sampled in four different hydrographic basins. Principal Component Analysis based on prey proportions identified distinct perilithon grazer and benthophagous behaviours. Furthermore, our results were consistent with the available literature on feeding behaviour in these fish. The simplicity of the PCR-based method and its potential generalization to other freshwater organisms may open new perspectives in food web ecology.

QDD: a user-friendly program to select microsatellite markers and design primers from large sequencing projects.

SUMMARY: QDD is an open access program providing a user-friendly tool for microsatellite detection and primer design from large sets of DNA sequences. The program is designed to deal with all steps of treatment of raw sequences obtained from pyrosequencing of enriched DNA libraries, but it is also applicable to data obtained through other sequencing methods, using FASTA files as input. The following tasks are completed by QDD: tag sorting, adapter/vector removal, elimination of redundant sequences, detection of possible genomic multicopies (duplicated loci or transposable elements), stringent selection of target microsatellites and customizable primer design. It can treat up to one million sequences of a few hundred base pairs in the tag-sorting step, and up to 50,000 sequences in a single input file for the steps involving estimation of sequence similarity. AVAILABILITY: QDD is freely available under the GPL licence for Windows and Linux from the following web site: http://www.univ-provence.fr/gsite/Local/egee/dir/meglecz/QDD.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Trade-off between morphological convergence and opportunistic diet behavior in fish hybrid zone.

BACKGROUND: The invasive Chondrostoma nasus nasus has colonized part of the distribution area of the protected endemic species Chondrostoma toxostoma toxostoma. This hybrid zone is a complex system where multiple effects such as inter-species competition, bi-directional introgression, strong environmental pressure and so on are combined. Why do sympatric Chondrostoma fish present a unidirectional change in body shape? Is this the result of inter-species interactions and/or a response to environmental effects or the result of trade-offs? Studies focusing on the understanding of a trade-off between multiple parameters are still rare. Although this has previously been done for Cichlid species flock and for Darwin finches, where mouth or beak morphology were coupled to diet and genetic identification, no similar studies have been done for a fish hybrid zone in a river. We tested the correlation between morphology (body and mouth morphology), diet (stable carbon and nitrogen isotopes) and genomic combinations in different allopatric and sympatric populations for a global data set of 1330 specimens. To separate the species interaction effect from the environmental effect in sympatry, we distinguished two data sets: the first one was obtained from a highly regulated part of the river and the second was obtained from specimens coming from the less regulated part. RESULTS: The distribution of the hybrid combinations was different in the two part of the sympatric zone, whereas all the specimens presented similar overall changes in body shape and in mouth morphology. Sympatric specimens were also characterized by a larger diet behavior variance than reference populations, characteristic of an opportunistic diet. No correlation was established between the body shape (or mouth deformation) and the stable isotope signature. CONCLUSION: The Durance River is an untamed Mediterranean river despite the presence of numerous dams that split the river from upstream to downstream. The sympatric effect on morphology and the large diet behavior range can be explained by a tendency toward an opportunistic behavior of the sympatric specimens. Indeed, the similar response of the two species and their hybrids implied an adaptation that could be defined as an alternative trade-off that underline the importance of epigenetics mechanisms for potential success in a novel environment.

Inter- and extra-Indian admixture and genetic diversity in reunion island revealed by analysis of mitochondrial DNA.

Reunion Island is a French territory located in the western Indian Ocean. The genetic pattern of the Reunionese population has been shaped by contributions from highly contrasting regions of the world. Over the last 350 years, several migration waves and cultural and socio-economic factors have led to the emergence of six main ethnic groups in Reunion. India is one of the principal regions that contributed to the setting up of the Reunionese population. Diversity, demographic and admixture analyses were performed on mtDNA variation of the Reunionese of Indian ancestry, including the Malbar and Zarab ethnic groups, in order to question their history. Using a phylogeographical approach, we generated and analysed quantitative data on the contribution of the Indian parental populations. Furthermore, we showed that the settlement of Reunion Island by Indians did not involve a founder effect, except in the very beginning of the Reunionese settlement (at the end of the 17(th) century). The accuracy of our results was optimised by a re-evaluation of the classification of the Southern Asian mtDNA haplogroups. Finally, by comparing our results to a previous study dealing with the Reunionese population, we highlighted how ethno-historical data are critical for reconstructing the complex history of multiethnic populations.

Novelties in hybrid zones: crossroads between population genomic and ecological approaches.

BACKGROUND: Interspecific hybridization is widespread, occurring in a taxonomically diverse array of species. The Cyprinidae family, which displays more than 30% hybridization, is a good candidate for studies of processes underlying isolation and speciation, such as genetic exchange between previously isolated lineages. This is particularly relevant in the case of recent hybridization between an invasive species, Chondrostoma nasus nasus (from Eastern Europe), and C. toxostoma toxostoma (a threatened species endemic to southern France), in which bidirectional introgressive hybridization has been demonstrated. METHODOLOGY/PRINCIPAL FINDINGS: We studied 128 specimens from reference populations and 1495 hybrid zone specimens (two years of sampling and four stations), using five molecular markers (one mitochondrial gene, four nuclear introns), morphology (meristic and plastic characters) and life history traits (weight, size, coefficient of condition, sex, age, shoaling). We identified 65 hybrid combinations and visualized spatial and temporal changes in composition. The direction of mitochondrial introgression was density-dependent in favor of the rarer species and we demonstrate that the sexual selection hypothesis is a preponderant explanation in the asymmetry of introgression. Despite genomic evolution in the hybrid zone, convergence was observed for body shape and coefficient of condition, indicating changes in foraging behavior with respect to reference populations, reflecting strong environmental pressure. CONCLUSIONS/SIGNIFICANCE: The complex rules of hybrid zone dynamics are established very early in the contact zone. We propose "inheritance from the rare species" as a new evolutionary hypothesis for animal models. The endemic species was not assimilated by the invasive species. Survival rates for this species were highest in the middle of the river (the warmest part) due to a trade-off between food availability and fecundity. The environment-independent hybrid combination may result from nuclear-mitochondrial interactions involving the Tpi1b gene or a gene linked to this gene (Chromosome 16). This genomic region is also responsible for shoaling behavior in Danio rerio and is a promising zone for studies of changes in population dynamics and advances in integrated studies of hybrid zones.

A new scale for measuring dynamic patterns of sexual partnership and concurrency: application to three French Caribbean regions.

BACKGROUND: Quantitative assessment of multiple sexual partnerships and concurrency may help to elucidate the large observed differences in the prevalence of AIDS among population subgroups and countries. GOAL: The goals of the study were (1) to develop a global scale of dynamic patterns of sexual partnerships, including concurrency with new partners and stable concurrency; (2) to apply this scale to three Caribbean regions characterized by different cumulative rates of incidence of AIDS; and (3) to compare the concurrency rates given by this scale with those of other published methods. STUDY DESIGN: We defined an individual scale based on 6 patterns of sexual behavior over the previous 12-month period, by using a simple algorithm to combine 7 variables. We then applied this scale to cross-sectional data collected from men living in three French Caribbean regions: Guadeloupe, Martinique, and Guyana. RESULTS: We found that all adults of all age classes in the three regions studied frequently had multiple (>2) and concurrent partnerships. The patterns of sexual behavior in the three regions were consistent with the respective cumulative incidence rates of AIDS, and a lower rate of concurrency with new partners and a higher rate of stable partnership concurrency were noted in Martinique, especially among 45- to 59-year-olds. The rate of concurrent partnerships was found to depend on the criteria used to define them and on the observation period (a given moment, or a defined period). Our definition gave a higher rate of concurrency than previously published indicators. CONCLUSION: The proposed scale can be applied to easy-to-collect data in cross-sectional population surveys and takes into account a wide variety of behaviors, including different types of concurrency.