Programmed DNA elimination in Mesorhabditis nematodes.

Some species undergo programmed DNA elimination (PDE), whereby portions of the genome are systematically destroyed in somatic cells. PDE has emerged independently in several phyla, but its function is unknown. Although the mechanisms are partially solved in ciliates, PDE remains mysterious in metazoans because the study species were not yet amenable to functional approaches. We fortuitously discovered massive PDE in the free-living nematode genus Mesorhabditis, from the same family as C. elegans. As such, these species offer many experimental advantages to start elucidating the PDE mechanisms in an animal. Here, we used cytology to describe the dynamics of chromosome fragmentation and destruction in early embryos. Elimination occurs once in development, at the third embryonic cell division in the somatic blastomeres. Chromosomes are first fragmented during S phase. Next, some of the fragments fail to align on the mitotic spindle and remain outside the re-assembled nuclei after mitosis. These fragments are gradually lost after a few cell cycles. The retained fragments form new mini chromosomes, which are properly segregated in the subsequent cell divisions. With genomic approaches, we found that Mesorhabditis mainly eliminate repeated regions and also about a hundred genes. Importantly, none of the eliminated protein-coding genes are shared between closely related Mesorhabditis species. Our results strongly suggest PDE has not been selected for regulating genes with important biological functions in Mesorhabditis but rather mainly to irreversibly remove repeated sequences in the soma. We propose that PDE may target genes, provided their elimination in the soma is invisible to selection.

Bacterial diversity on an abandoned, industrial wasteland contaminated by polychlorinated biphenyls, dioxins, furans and trace metals.

Most former industrial sites are contaminated by mixtures of trace elements and organic pollutants. Levels of pollutants do not provide information regarding their biological impact, bioavailability and possible interactions between substances. There is genuine interest in combining chemical analyses with biological investigations. We studied a brownfield where several industrial activities were carried out starting in the 1970s, (incineration of pyralene transformers, recovery of copper by burning cables in the open air). Four representative plots showing different levels of polychlorobiphenyls were selected. Organic and trace metal levels were measured together with soil pedological characteristics. The bacterial community structure and functional diversity were assessed by 16S metagenomics with deep sequencing and community-level physiological profiling. Additionally, a vegetation survey was performed. Polychlorobiphenyls (8 mg.kg-1 to 1500 mg.kg-1) were from 2.4 × 103-fold to 6 × 105-fold higher than the European background level of 2.5 µg.kg-1. Polychlorinated dibenzo-p-dioxins and dibenzofurans ranged from 0.5 to 8.0 µg.kg-1. The soil was also contaminated with trace metals, i.e., Cu > 187, Zn > 217 and Pb > 372 mg.kg-1. Location within the study area, trace metal content and soil humidity were stronger determinants than organic pollutants of bacterial community structures and activities. Thus, the highest biological activity and the greatest bacteriological richness were observed in the plot that was less contaminated with trace metals, despite the high level of organic pollutants in the plot. Moreover, trace element pollution was associated with a relatively low presence of Actinobacteria and Rhizobia. The plot with the highest metal contamination was rich in metal-resistant bacteria such as Sphingomonadales, Geodermatophilaceae and KD4-96 (Chloroflexi phylum). Acidobacteria and Sphingomonadales, capable of resisting trace metals and degrading persistent organic pollutants, were dominant in the plots that had accumulated metal and organic contamination, but bacterial activity was lower in these plots than in the other plots.

Within-Generation Polygenic Selection Shapes Fitness-Related Traits across Environments in Juvenile Sea Bream.

Understanding the genetic underpinnings of fitness trade-offs across spatially variable environments remains a major challenge in evolutionary biology. In Mediterranean gilthead sea bream, first-year juveniles use various marine and brackish lagoon nursery habitats characterized by a trade-off between food availability and environmental disturbance. Phenotypic differences among juveniles foraging in different habitats rapidly appear after larval settlement, but the relative role of local selection and plasticity in phenotypic variation remains unclear. Here, we combine phenotypic and genetic data to address this question. We first report correlations of opposite signs between growth and condition depending on juvenile habitat type. Then, we use single nucleotide polymorphism (SNP) data obtained by Restriction Associated DNA (RAD) sequencing to search for allele frequency changes caused by a single generation of spatially varying selection between habitats. We found evidence for moderate selection operating at multiple loci showing subtle allele frequency shifts between groups of marine and brackish juveniles. We identified subsets of candidate outlier SNPs that, in interaction with habitat type, additively explain up to 3.8% of the variance in juvenile growth and 8.7% in juvenile condition; these SNPs also explained significant fraction of growth rate in an independent larval sample. Our results indicate that selective mortality across environments during early-life stages involves complex trade-offs between alternative growth strategies.

Detecting adaptive convergent amino acid evolution.

In evolutionary genomics, researchers have taken an interest in identifying substitutions that subtend convergent phenotypic adaptations. This is a difficult question that requires distinguishing foreground convergent substitutions that are involved in the convergent phenotype from background convergent substitutions. Those may be linked to other adaptations, may be neutral or may be the consequence of mutational biases. Furthermore, there is no generally accepted definition of convergent substitutions. Various methods that use different definitions have been proposed in the literature, resulting in different sets of candidate foreground convergent substitutions. In this article, we first describe the processes that can generate foreground convergent substitutions in coding sequences, separating adaptive from non-adaptive processes. Second, we review methods that have been proposed to detect foreground convergent substitutions in coding sequences and expose the assumptions that underlie them. Finally, we examine their power on simulations of convergent changes-including in the presence of a change in the efficacy of selection-and on empirical alignments. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Dynamic Interactions Between the Genome and an Endogenous Retrovirus: Tirant in Drosophila simulans Wild-Type Strains.

All genomes contain repeated sequences that are known as transposable elements (TEs). Among these are endogenous retroviruses (ERVs), which are sequences similar to retroviruses and are transmitted across generations from parent to progeny. These sequences are controlled in genomes through epigenetic mechanisms. At the center of the epigenetic control of TEs are small interfering RNAs of the piRNA class, which trigger heterochromatinization of TE sequences. The tirant ERV of Drosophila simulans displays intra-specific variability in copy numbers, insertion sites, and transcription levels, providing us with a well-suited model to study the dynamic relationship between a TE family and the host genome through epigenetic mechanisms. We show that tirant transcript amounts and piRNA amounts are positively correlated in ovaries in normal conditions, unlike what was previously described following divergent crosses. In addition, we describe tirant insertion polymorphism in the genomes of three D. simulans wild-type strains, which reveals a limited number of insertions that may be associated with gene transcript level changes through heterochromatin spreading and have phenotypic impacts. Taken together, our results participate in the understanding of the equilibrium between the host genome and its TEs.

Developmental and comparative transcriptomic identification of iridophore contribution to white barring in clownfish.

Actinopterygian fishes harbor at least eight distinct pigment cell types, leading to a fascinating diversity of colors. Among this diversity, the cellular origin of the white color appears to be linked to several pigment cell types such as iridophores or leucophores. We used the clownfish Amphiprion ocellaris, which has a color pattern consisting of white bars over a darker body, to characterize the pigment cells that underlie the white hue. We observe by electron microscopy that cells in white bars are similar to iridophores. In addition, the transcriptomic signature of clownfish white bars exhibits similarities with that of zebrafish iridophores. We further show by pharmacological treatments that these cells are necessary for the white color. Among the top differentially expressed genes in white skin, we identified several genes (fhl2a, fhl2b, saiyan, gpnmb, and apoD1a) and show that three of them are expressed in iridophores. Finally, we show by CRISPR/Cas9 mutagenesis that these genes are critical for iridophore development in zebrafish. Our analyses provide clues to the genomic underpinning of color diversity and allow identification of new iridophore genes in fish.

CAARS: comparative assembly and annotation of RNA-Seq data.

MOTIVATION: RNA sequencing (RNA-Seq) is a widely used approach to obtain transcript sequences in non-model organisms, notably for performing comparative analyses. However, current bioinformatic pipelines do not take full advantage of pre-existing reference data in related species for improving RNA-Seq assembly, annotation and gene family reconstruction. RESULTS: We built an automated pipeline named CAARS to combine novel data from RNA-Seq experiments with existing multi-species gene family alignments. RNA-Seq reads are assembled into transcripts by both de novo and assisted assemblies. Then, CAARS incorporates transcripts into gene families, builds gene alignments and trees and uses phylogenetic information to classify the genes as orthologs and paralogs of existing genes. We used CAARS to assemble and annotate RNA-Seq data in rodents and fishes using distantly related genomes as reference, a difficult case for this kind of analysis. We showed CAARS assemblies are more complete and accurate than those assembled by a standard pipeline consisting of de novo assembly coupled with annotation by sequence similarity on a guide species. In addition to annotated transcripts, CAARS provides gene family alignments and trees, annotated with orthology relationships, directly usable for downstream comparative analyses. AVAILABILITY AND IMPLEMENTATION: CAARS is implemented in Python and Ocaml and is freely available at https://github.com/carinerey/caars. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Accurate Detection of Convergent Amino-Acid Evolution with PCOC.

In the history of life, some phenotypes have been acquired several times independently, through convergent evolution. Recently, lots of genome-scale studies have been devoted to identify nucleotides or amino acids that changed in a convergent manner when the convergent phenotypes evolved. These efforts have had mixed results, probably because of differences in the detection methods, and because of conceptual differences about the definition of a convergent substitution. Some methods contend that substitutions are convergent only if they occur on all branches where the phenotype changed toward the exact same state at a given nucleotide or amino acid position. Others are much looser in their requirements and define a convergent substitution as one that leads the site at which they occur to prefer a phylogeny in which species with the convergent phenotype group together. Here, we suggest to look for convergent shifts in amino acid preferences instead of convergent substitutions to the exact same amino acid. We define as convergent shifts substitutions that occur on all branches where the phenotype changed and such that they correspond to a change in the type of amino acid preferred at this position. We implement the corresponding model into a method named PCOC. We show on simulations that PCOC better recovers convergent shifts than existing methods in terms of sensitivity and specificity. We test it on a plant protein alignment where convergent evolution has been studied in detail and find that our method recovers several previously identified convergent substitutions and proposes credible new candidates.

RiboDB Database: A Comprehensive Resource for Prokaryotic Systematics.

Ribosomal proteins (r-proteins) are increasingly used as an alternative to ribosomal rRNA for prokaryotic systematics. However, their routine use is difficult because r-proteins are often not or wrongly annotated in complete genome sequences, and there is currently no dedicated exhaustive database of r-proteins. RiboDB aims at fulfilling this gap. This weekly updated comprehensive database allows the fast and easy retrieval of r-protein sequences from publicly available complete prokaryotic genome sequences. The current version of RiboDB contains 90 r-proteins from 3,750 prokaryotic complete genomes encompassing 38 phyla/major classes and 1,759 different species. RiboDB is accessible at http://ribodb.univ-lyon1.fr and through ACNUC interfaces.