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.