Population transcriptogenomics highlights impaired metabolism and small population sizes in tree frogs living in the Chernobyl Exclusion Zone.

BACKGROUND: Individual functional modifications shape the ability of wildlife populations to cope with anthropogenic environmental changes. But instead of adaptive response, human-altered environments can generate a succession of deleterious functional changes leading to the extinction of the population. To study how persistent anthropogenic changes impacted local species' population status, we characterised population structure, genetic diversity and individual response of gene expression in the tree frog Hyla orientalis along a gradient of radioactive contamination around the Chernobyl nuclear power plant. RESULTS: We detected lower effective population size in populations most exposed to ionizing radiation in the Chernobyl Exclusion Zone that is not compensated by migrations from surrounding areas. We also highlighted a decreased body condition of frogs living in the most contaminated area, a distinctive transcriptomics signature and stop-gained mutations in genes involved in energy metabolism. While the association with dose will remain correlational until further experiments, a body of evidence suggests the direct or indirect involvement of radiation exposure in these changes. CONCLUSIONS: Despite ongoing migration and lower total dose rates absorbed than at the time of the accident, our results demonstrate that Hyla orientalis specimens living in the Chernobyl Exclusion Zone are still undergoing deleterious changes, emphasizing the long-term impacts of the nuclear disaster.

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.

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.

Evolutionary approach for pollution study: The case of ionizing radiation.

Estimating the consequences of environmental changes, specifically in a global change context, is essential for conservation issues. In the case of pollutants, the interest in using an evolutionary approach to investigate their consequences has been emphasized since the 2000s, but these studies remain rare compared to the characterization of direct effects on individual features. We focused on the study case of anthropogenic ionizing radiation because, despite its potential strong impact on evolution, the scarcity of evolutionary approaches to study the biological consequences of this stressor is particularly true. In this study, by investigating some particular features of the biological effects of this stressor, and by reviewing existing studies on evolution under ionizing radiation, we suggest that evolutionary approach may help provide an integrative view on the biological consequences of ionizing radiation. We focused on three topics: (i) the mutagenic properties of ionizing radiation and its disruption of evolutionary processes, (ii) exposures at different time scales, leading to an interaction between past and contemporary evolution, and (iii) the special features of contaminated areas called exclusion zones and how evolution could match field and laboratory observed effects. This approach can contribute to answering several key issues in radioecology: to explain species differences in the sensitivity to ionizing radiation, to improve our estimation of the impacts of ionizing radiation on populations, and to help identify the environmental features impacting organisms (e.g., interaction with other pollution, migration of populations, anthropogenic environmental changes). Evolutionary approach would benefit from being integrated to the ecological risk assessment process.

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.

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.

Evidence of en bloc duplication in vertebrate genomes.

It has been 30 years since it was first proposed that the vertebrate genome evolved through several rounds of genome-wide duplications (polyploidizations). Despite rapid advances in genetics, including sequencing of the complete genomes of several divergent species, this hypothesis has not been tested rigorously and is still a matter of debate. If polyploidizations occurred during chordate evolution, there should be a network of paralogous regions in the present-day jawed vertebrate (Gnathostomata) genomes. Here we present an investigation of the major histocompatibility complex (MHC) paralogous regions, which we accomplished by characterizing the corresponding region in amphioxus by identifying nine anchor genes and sequencing both the anchor genes and the regions that flank them (a total of 400 kb). Phylogenetic analysis of 31 genes (including the anchor genes) in these regions shows that duplications occurred after the divergence of cephalochordates and vertebrates but before the Gnathostomata radiation. The distribution of human and amphioxus orthologs in their respective genomes and the relationship between these distributions support the en bloc duplication events. Our analysis represents the first step towards demonstrating that the human ancestral genome has undergone polyploidization. Moreover, reconstruction of the pre-duplicated region indicates that one of the duplicated regions retains the ancestral organization.

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.

Diffuse large B-cell lymphomas with CDKN2A deletion have a distinct gene expression signature and a poor prognosis under R-CHOP treatment: a GELA study.

Genomic alterations play a crucial role in the development and progression of diffuse large B-cell lymphomas (DLBCLs). We determined gene copy number alterations (GCNAs) of TP53, CDKN2A, CDKN1B, BCL2, MYC, REL, and RB1 with a single polymerase chain reaction (PCR) assay (quantitative multiplex PCR of short fragments [QMPSF]) in a cohort of 114 patients with DLBCL to assess their prognostic value and relationship with the gene expression profile. Losses of TP53 and CDKN2A, observed in 8% and 35% of patients, respectively, were significantly associated with a shorter survival after rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisone (R-CHOP) treatment, independently of the International Prognostic Index and of the cell of origin. Analysis of the 9p21 genomic region indicated that transcripts encoding p14ARF and p16INK4A were both disrupted in most patients with CDKN2A deletion. These patients predominantly had an activated B-cell profile and showed a specific gene expression signature, characterized by dysregulation of the RB/E2F pathway, activation of cellular metabolism, and decreased immune and inflammatory responses. These features may constitute the molecular basis sustaining the unfavorable outcome and chemoresistance of this DLBCL subgroup. Detection of TP53 and CDKN2A loss by QMPSF is a powerful tool that could be used for patient stratification in future clinical trials.

Unusual evolution of tree frog populations in the Chernobyl exclusion zone.

Despite the ubiquity of pollutants in the environment, their long-term ecological consequences are not always clear and still poorly studied. This is the case concerning the radioactive contamination of the environment following the major nuclear accident at the Chernobyl nuclear power plant. Notwithstanding the implications of evolutionary processes on the population status, few studies concern the evolution of organisms chronically exposed to ionizing radiation in the Chernobyl exclusion zone. Here, we examined genetic markers for 19 populations of Eastern tree frog (Hyla orientalis) sampled in the Chernobyl region about thirty years after the nuclear power plant accident to investigate microevolutionary processes ongoing in local populations. Genetic diversity estimated from nuclear and mitochondrial markers showed an absence of genetic erosion and higher mitochondrial diversity in tree frogs from the Chernobyl exclusion zone compared to other European populations. Moreover, the study of haplotype network permitted us to decipher the presence of an independent recent evolutionary history of Chernobyl exclusion zone's Eastern tree frogs caused by an elevated mutation rate compared to other European populations. By fitting to our data a model of haplotype network evolution, we suspected that Eastern tree frog populations in the Chernobyl exclusion zone have a high mitochondrial mutation rate and small effective population sizes. These data suggest that Eastern tree frog populations might offset the impact of deleterious mutations because of their large clutch size, but also question the long-term impact of ionizing radiation on the status of other species living in the Chernobyl exclusion zone.

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.

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.

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.

Lenalidomide/rituximab induces high molecular response in untreated follicular lymphoma: LYSA ancillary RELEVANCE study.

Complete molecular response (CMR) after first-line immunochemotherapy reflects treatment efficacy and may predict prognosis in patients with follicular lymphoma (FL). RELEVANCE is the first phase 3 trial comparing the chemotherapy-free regimen lenalidomide/rituximab (R2) vs rituximab/chemotherapy (R-Chemo) in previously untreated FL patients (ClinicalTrials.gov identifier: NCT01650701). The objective of the minimal residual disease (MRD) analysis was to determine the ability of a chemotherapy-free regimen to induce CMR. Of 440 French patients participating in the Lymphoma Study Association (LYSA) RELEVANCE MRD study, all 222 patients with a BIOMED-2-detectable BCL2-JH translocation at diagnosis were analyzed. MRD was quantified by droplet digital polymerase chain reaction with a sensitivity ≤10-4. At week 24 (end of induction treatment), 98% and 78% of patients achieved CMR in peripheral blood (PB) and bone marrow (BM), respectively. Achievement of CMR (in PB and/or BM) had a significant impact on progression-free survival (PFS), with 3-year PFS of 84% and 55% for patients with CMR and detectable MRD, respectively (P = .015). CMR at week 24 was reached more frequently in the R2 arm (105/117; 90%) than in the R-Chemo arm (70/90; 77%) (P = .022). The poor prognostic value in terms of PFS for the persistence of molecular disease was observed irrespective of treatment arm (interaction test, P = .31). In agreement with the clinical results of the RELEVANCE trial, our results show that R2 immunomodulatory treatment in first-line FL can achieve high rates of CMR.

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.

A euryhaline fish, lost in the desert: The unexpected metapopulation structure of Coptodon guineensis (Günther, 1862) in the Sebkha of Imlili.

Euryhaline Cichlid fish of the species Coptodon guineensis are present in different water holes situated in a dried depression in the desert in the extreme South of Morocco, the Sebkha of Imlili. A genetic survey of this population, using complete sequences of the ND2 gene (mtDNA) and sixteen microsatellite loci, revealed that the fish in the sebkha did not form a single population, but rather a metapopulation. This metapopulational structure may be regarded as good news from the point of view of the conservation of fish in the sebkha. Although small individual populations may have short, finite life spans, the metapopulation as a whole is more stable, because immigrants from one population are likely to re-colonize the habitat, left open by the extinction of another.

Clinical application of targeted next-generation sequencing for colorectal cancer patients: a multicentric Belgian experience.

International guidelines made RAS (KRAS and NRAS) status a prerequisite for the use of anti-EGFR agents for metastatic colorectal cancer (CRC) patients. Daily, new data emerges on the theranostic and prognostic role of molecular biomarkers; this is a strong incentive for a validated, sensitive, and broadly available molecular screening test. Next-generation sequencing (NGS) has begun to supplant other technologies for genomic profiling. We report here our 2 years of clinical practice using NGS results to guide therapeutic decisions. The Ion Torrent AmpliSeq colon/lung cancer panel, which allows mutation detection in 22 cancer-related genes, was prospectively used in clinical practice (BELAC ISO 15189 accredited method). The DNA of 741 formalin-fixed paraffin-embedded CRC tissues, including primary tumors and metastasis, was obtained from 14 different Belgian institutions and subjected to targeted NGS. Of the tumors tested, 98% (727) were successfully sequenced and 89% (650) harbored at least one mutation. KRAS, BRAF and NRAS mutations were found in 335 (46%), 78 (11%) and 32 (4%) samples, respectively. These mutation frequencies were consistent with those reported in public databases. Moreover, mutations and amplifications in potentially actionable genes were identified in 464 samples (64%), including mutations in PIK3CA (14%), ERBB2 (0.4%), AKT1 (0.6%), and MAP2K1 (0.1%), as well as amplifications of ERBB2 (0.3%) and EGFR (0.3%). The median turnaround time between reception of the sample in the laboratory and report release was 8 calendar days. Overall, the AmpliSeq colon/lung cancer panel was successfully applied in daily practice and provided reliable clinically relevant information for CRC patients.

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.

[Urbilateria, a complex organism?]. / Urbilateria, un être évolué?

In order to establish the portrait of Urbilateria, the common ancestor of triblastic metazoan, this paper focuses on the antero-posterior segmentation frequently considered as characterising the bilaterian bauplan. The synthesis presented here describes the morphological, anatomical and functional aspects of this organisation. Furthermore it analyses the conditions of its emergence during the ontogenesis of Annelids, Arthropods and Chordates and identifies its genetic bases. The provided data exhibit the unitary character of the segmentation modalities among protostomian and deuterostomian organisms. This process occurs in two phases, involving a posterior proliferative zone after the gastrulation. It shows the similarity of the expression patterns of orthologous genes, the implication of comparable signalisation and regulation pathways. The congruence of the results obtained at both structural and molecular levels reinforce the segmental organisation conception of the common ancestor of Bilaterians.