BioConvert: a comprehensive format converter for life sciences.

Bioinformatics is a field known for the numerous standards and formats that have been developed over the years. This plethora of formats, sometimes complementary, and often redundant, poses many challenges to bioinformatics data analysts. They constantly need to find the best tool to convert their data into the suitable format, which is often a complex, technical and time consuming task. Moreover, these small yet important tasks are often difficult to make reproducible. To overcome these difficulties, we initiated BioConvert, a collaborative project to facilitate the conversion of life science data from one format to another. BioConvert aggregates existing software within a single framework and complemented them with original code when needed. It provides a common interface to make the user experience more streamlined instead of having to learn tens of them. Currently, BioConvert supports about 50 formats and 100 direct conversions in areas such as alignment, sequencing, phylogeny, and variant calling. In addition to being useful for end-users, BioConvert can also be utilized by developers as a universal benchmarking framework for evaluating and comparing numerous conversion tools. Additionally, we provide a web server implementing an online user-friendly interface to BioConvert, hence allowing direct use for the community.

Leishmania allelic selection during experimental sand fly infection correlates with mutational signatures of oxidative DNA damage.

Trypanosomatid pathogens are transmitted by blood-feeding insects, causing devastating human infections. These parasites show important phenotypic shifts that often impact parasite pathogenicity, tissue tropism, or drug susceptibility. The evolutionary mechanisms that allow for the selection of such adaptive phenotypes remain only poorly investigated. Here, we use Leishmania donovani as a trypanosomatid model pathogen to assess parasite evolutionary adaptation during experimental sand fly infection. Comparing the genome of the parasites before and after sand fly infection revealed a strong population bottleneck effect as judged by allele frequency analysis. Apart from random genetic drift caused by the bottleneck effect, our analyses revealed haplotype and allelic changes during sand fly infection that seem under natural selection given their convergence between independent biological replicates. Our analyses further uncovered signature mutations of oxidative DNA damage in the parasite genomes after sand fly infection, suggesting that Leishmania suffers from oxidative stress inside the insect digestive tract. Our results propose a model of Leishmania genomic adaptation during sand fly infection, with oxidative DNA damage and DNA repair processes likely driving haplotype and allelic selection. The experimental and computational framework presented here provides a useful blueprint to assess evolutionary adaptation of other eukaryotic pathogens inside their insect vectors, such as Plasmodium spp, Trypanosoma brucei, and Trypanosoma cruzi.

piRNAs initiate transcriptional silencing of spermatogenic genes during C. elegans germline development.

Eukaryotic genomes harbor invading transposable elements that are silenced by PIWI-interacting RNAs (piRNAs) to maintain genome integrity in animal germ cells. However, whether piRNAs also regulate endogenous gene expression programs remains unclear. Here, we show that C. elegans piRNAs trigger the transcriptional silencing of hundreds of spermatogenic genes during spermatogenesis, promoting sperm differentiation and function. This silencing signal requires piRNA-dependent small RNA biogenesis and loading into downstream nuclear effectors, which correlates with the dynamic reorganization of two distinct perinuclear biomolecular condensates present in germ cells. In addition, the silencing capacity of piRNAs is temporally counteracted by the Argonaute CSR-1, which targets and licenses spermatogenic gene transcription. The spatial and temporal overlap between these opposing small RNA pathways contributes to setting up the timing of the spermatogenic differentiation program. Thus, our work identifies a prominent role for piRNAs as direct regulators of endogenous transcriptional programs during germline development and gamete differentiation.

Translation and codon usage regulate Argonaute slicer activity to trigger small RNA biogenesis.

In the Caenorhabditis elegans germline, thousands of mRNAs are concomitantly expressed with antisense 22G-RNAs, which are loaded into the Argonaute CSR-1. Despite their essential functions for animal fertility and embryonic development, how CSR-1 22G-RNAs are produced remains unknown. Here, we show that CSR-1 slicer activity is primarily involved in triggering the synthesis of small RNAs on the coding sequences of germline mRNAs and post-transcriptionally regulates a fraction of targets. CSR-1-cleaved mRNAs prime the RNA-dependent RNA polymerase, EGO-1, to synthesize 22G-RNAs in phase with translating ribosomes, in contrast to other 22G-RNAs mostly synthesized in germ granules. Moreover, codon optimality and efficient translation antagonize CSR-1 slicing and 22G-RNAs biogenesis. We propose that codon usage differences encoded into mRNA sequences might be a conserved strategy in eukaryotes to regulate small RNA biogenesis and Argonaute targeting.

Germline inherited small RNAs facilitate the clearance of untranslated maternal mRNAs in C. elegans embryos.

Inheritance and clearance of maternal mRNAs are two of the most critical events required for animal early embryonic development. However, the mechanisms regulating this process are still largely unknown. Here, we show that together with maternal mRNAs, C. elegans embryos inherit a complementary pool of small non-coding RNAs that facilitate the cleavage and removal of hundreds of maternal mRNAs. These antisense small RNAs are loaded into the maternal catalytically-active Argonaute CSR-1 and cleave complementary mRNAs no longer engaged in translation in somatic blastomeres. Induced depletion of CSR-1 specifically during embryonic development leads to embryonic lethality in a slicer-dependent manner and impairs the degradation of CSR-1 embryonic mRNA targets. Given the conservation of Argonaute catalytic activity, we propose that a similar mechanism operates to clear maternal mRNAs during the maternal-to-zygotic transition across species.

Small-RNA-mediated transgenerational silencing of histone genes impairs fertility in piRNA mutants.

PIWI-interacting RNAs (piRNAs) promote fertility in many animals. However, whether this is due to their conserved role in repressing repetitive elements (REs) remains unclear. Here, we show that the progressive loss of fertility in Caenorhabditis elegans lacking piRNAs is not caused by derepression of REs or other piRNA targets but, rather, is mediated by epigenetic silencing of all of the replicative histone genes. In the absence of piRNAs, downstream components of the piRNA pathway relocalize from germ granules and piRNA targets to histone mRNAs to synthesize antisense small RNAs (sRNAs) and induce transgenerational silencing. Removal of the downstream components of the piRNA pathway restores histone mRNA expression and fertility in piRNA mutants, and the inheritance of histone sRNAs in wild-type worms adversely affects their fertility for multiple generations. We conclude that sRNA-mediated silencing of histone genes impairs the fertility of piRNA mutants and may serve to maintain piRNAs across evolution.

The somatic piRNA pathway controls germline transposition over generations.

Transposable elements (TEs) are parasitic DNA sequences that threaten genome integrity by replicative transposition in host gonads. The Piwi-interacting RNAs (piRNAs) pathway is assumed to maintain Drosophila genome homeostasis by downregulating transcriptional and post-transcriptional TE expression in the ovary. However, the bursts of transposition that are expected to follow transposome derepression after piRNA pathway impairment have not yet been reported. Here, we show, at a genome-wide level, that piRNA loss in the ovarian somatic cells boosts several families of the endogenous retroviral subclass of TEs, at various steps of their replication cycle, from somatic transcription to germinal genome invasion. For some of these TEs, the derepression caused by the loss of piRNAs is backed up by another small RNA pathway (siRNAs) operating in somatic tissues at the post transcriptional level. Derepressed transposition during 70 successive generations of piRNA loss exponentially increases the genomic copy number by up to 10-fold.

Piwi Is Required during Drosophila Embryogenesis to License Dual-Strand piRNA Clusters for Transposon Repression in Adult Ovaries.

Most piRNAs in the Drosophila female germline are transcribed from heterochromatic regions called dual-strand piRNA clusters. Histone 3 lysine 9 trimethylation (H3K9me3) is required for licensing piRNA production by these clusters. However, it is unclear when and how they acquire this permissive heterochromatic state. Here, we show that transient Piwi depletion in Drosophila embryos results in H3K9me3 decrease at piRNA clusters in ovaries. This is accompanied by impaired biogenesis of ovarian piRNAs, accumulation of transposable element transcripts, and female sterility. Conversely, Piwi depletion at later developmental stages does not disturb piRNA cluster licensing. These results indicate that the identity of piRNA clusters is epigenetically acquired in a Piwi-dependent manner during embryonic development, which is reminiscent of the widespread genome reprogramming occurring during early mammalian zygotic development.

microRNA target prediction programs predict many false positives.

According to the current view, each microRNA regulates hundreds of genes. Computational tools aim at identifying microRNA targets, usually selecting evolutionarily conserved microRNA binding sites. While the false positive rates have been evaluated for some prediction programs, that information is rarely put forward in studies making use of their predictions. Here, we provide evidence that such predictions are often biologically irrelevant. Focusing on miR-223-guided repression, we observed that it is often smaller than inter-individual variability in gene expression among wild-type mice, suggesting that most predicted targets are functionally insensitive to that microRNA. Furthermore, we found that human haplo-insufficient genes tend to bear the most highly conserved microRNA binding sites. It thus appears that biological functionality of microRNA binding sites depends on the dose-sensitivity of their host gene and that, conversely, it is unlikely that every predicted microRNA target is dose-sensitive enough to be functionally regulated by microRNAs. We also observed that some mRNAs can efficiently titrate microRNAs, providing a reason for microRNA binding site conservation for inefficiently repressed targets. Finally, many conserved microRNA binding sites are conserved in a microRNA-independent fashion: Sequence elements may be conserved for other reasons, while being fortuitously complementary to microRNAs. Collectively, our data suggest that the role of microRNAs in normal and pathological conditions has been overestimated due to the frequent overlooking of false positive rates.

MicroRNA-Dependent Transcriptional Silencing of Transposable Elements in Drosophila Follicle Cells.

RNA interference-related silencing mechanisms concern very diverse and distinct biological processes, from gene regulation (via the microRNA pathway) to defense against molecular parasites (through the small interfering RNA and the Piwi-interacting RNA pathways). Small non-coding RNAs serve as specificity factors that guide effector proteins to ribonucleic acid targets via base-pairing interactions, to achieve transcriptional or post-transcriptional regulation. Because of the small sequence complementarity required for microRNA-dependent post-transcriptional regulation, thousands of microRNA (miRNA) putative targets have been annotated in Drosophila. In Drosophila somatic ovarian cells, genomic parasites, such as transposable elements (TEs), are transcriptionally repressed by chromatin changes induced by Piwi-interacting RNAs (piRNAs) that prevent them from invading the germinal genome. Here we show, for the first time, that a functional miRNA pathway is required for the piRNA-mediated transcriptional silencing of TEs in this tissue. Global miRNA depletion, caused by tissue- and stage-specific knock down of drosha (involved in miRNA biogenesis), AGO1 or gawky (both responsible for miRNA activity), resulted in loss of TE-derived piRNAs and chromatin-mediated transcriptional de-silencing of TEs. This specific TE de-repression was also observed upon individual titration (by expression of the complementary miRNA sponge) of two miRNAs (miR-14 and miR-34) as well as in a miR-14 loss-of-function mutant background. Interestingly, the miRNA defects differentially affected TE- and 3' UTR-derived piRNAs. To our knowledge, this is the first indication of possible differences in the biogenesis or stability of TE- and 3' UTR-derived piRNAs. This work is one of the examples of detectable phenotypes caused by loss of individual miRNAs in Drosophila and the first genetic evidence that miRNAs have a role in the maintenance of genome stability via piRNA-mediated TE repression.

Compositional biases among synonymous substitutions cause conflict between gene and protein trees for plastid origins.

Archaeplastida (=Kingdom Plantae) are primary plastid-bearing organisms that evolved via the endosymbiotic association of a heterotrophic eukaryote host cell and a cyanobacterial endosymbiont approximately 1,400 Ma. Here, we present analyses of cyanobacterial and plastid genomes that show strongly conflicting phylogenies based on 75 plastid (or nuclear plastid-targeted) protein-coding genes and their direct translations to proteins. The conflict between genes and proteins is largely robust to the use of sophisticated data- and tree-heterogeneous composition models. However, by using nucleotide ambiguity codes to eliminate synonymous substitutions due to codon-degeneracy, we identify a composition bias, and dependent codon-usage bias, resulting from synonymous substitutions at all third codon positions and first codon positions of leucine and arginine, as the main cause for the conflicting phylogenetic signals. We argue that the protein-coding gene data analyses are likely misleading due to artifacts induced by convergent composition biases at first codon positions of leucine and arginine and at all third codon positions. Our analyses corroborate previous studies based on gene sequence analysis that suggest Cyanobacteria evolved by the early paraphyletic splitting of Gloeobacter and a specific Synechococcus strain (JA33Ab), with all other remaining cyanobacterial groups, including both unicellular and filamentous species, forming the sister-group to the Archaeplastida lineage. In addition, our analyses using better-fitting models suggest (but without statistically strong support) an early divergence of Glaucophyta within Archaeplastida, with the Rhodophyta (red algae), and Viridiplantae (green algae and land plants) forming a separate lineage.

Unravelling the distinct strains of Tharu ancestry.

The northern region of the Indian subcontinent is a vast landscape interlaced by diverse ecologies, for example, the Gangetic Plain and the Himalayas. A great number of ethnic groups are found there, displaying a multitude of languages and cultures. The Tharu is one of the largest and most linguistically diverse of such groups, scattered across the Tarai region of Nepal and bordering Indian states. Their origins are uncertain. Hypotheses have been advanced postulating shared ancestry with Austroasiatic, or Tibeto-Burman-speaking populations as well as aboriginal roots in the Tarai. Several Tharu groups speak a variety of Indo-Aryan languages, but have traditionally been described by ethnographers as representing East Asian phenotype. Their ancestry and intra-population diversity has previously been tested only for haploid (mitochondrial DNA and Y-chromosome) markers in a small portion of the population. This study presents the first systematic genetic survey of the Tharu from both Nepal and two Indian states of Uttarakhand and Uttar Pradesh, using genome-wide SNPs and haploid markers. We show that the Tharu have dual genetic ancestry as up to one-half of their gene pool is of East Asian origin. Within the South Asian proportion of the Tharu genetic ancestry, we see vestiges of their common origin in the north of the South Asian Subcontinent manifested by mitochondrial DNA haplogroup M43.

rely.py, a python script to detect reliable clades.

rely.py is a program implementing the method to detect independently repeated clades by comparing phylogenies as described in Li and Lecointre (2009) and adapted to incompletely overlapping datasets in Li et al. (2009). The comparison can be performed on trees obtained by any inference method (maximum parsimony, Bayesian inference, maximum likelihood). The program computes repetition indices, provides greedy summary trees for each validity domain and a nexus matrix representation of the clades weighted by their repetition indices. The additional script concatnexus.py assists the user in preparing the primary analyses, but it can also be used separately to concatenate nexus datasets.

SuperTRI: A new approach based on branch support analyses of multiple independent data sets for assessing reliability of phylogenetic inferences.

Supermatrix and supertree are two methods for constructing a phylogenetic tree by using multiple data sets. However, these methods are not a panacea, as conflicting signals between data sets can lead to misinterpret the evolutionary history of taxa. In particular, the supermatrix approach is expected to be misleading if the species-tree signal is not dominant after the combination of the data sets. Moreover, most current supertree methods suffer from two limitations: (i) they ignore or misinterpret secondary (non-dominant) phylogenetic signals of the different data sets; and (ii) the logical basis of node robustness measures is unclear. To overcome these limitations, we propose a new approach, called SuperTRI, which is based on the branch support analyses of the independent data sets, and where the reliability of the nodes is assessed using three measures: the supertree Bootstrap percentage and two other values calculated from the separate analyses: the mean branch support (mean Bootstrap percentage or mean posterior probability) and the reproducibility index. The SuperTRI approach is tested on a data matrix including seven genes for 82 taxa of the family Bovidae (Mammalia, Ruminantia), and the results are compared to those found with the supermatrix approach. The phylogenetic analyses of the supermatrix and independent data sets were done using four methods of tree reconstruction: Bayesian inference, maximum likelihood, and unweighted and weighted maximum parsimony. The results indicate, firstly, that the SuperTRI approach shows less sensitivity to the four phylogenetic methods, secondly, that it is more accurate to interpret the relationships among taxa, and thirdly, that interesting conclusions on introgression and radiation can be drawn from the comparisons between SuperTRI and supermatrix analyses.

RNF213, a new nuclear marker for acanthomorph phylogeny.

We show that RNF213 is a nuclear gene suitable for investigating large scale acanthomorph teleosteans interrelationships. The gene recovers many clades already found by several independent studies of acanthomorph molecular phylogenetics and considered as reliable. Moreover, we performed phylogenetic analyses of three other independent nuclear markers, first separately and then of all possible combinations (Dettaï, A., Lecointre, G., 2004. In search of nothothenioid (Teleostei) relatives. Antarct. Sci. 16 (1), 71-85. URL http://dx.doi.org/10.1017/S0954102004) of the four genes. This was coupled with an assessment of the reliability of clades using the repetition index of Li and Lecointre (Li, B., Lecointre, G., 2008. Formalizing reliability in the taxonomic congruence approach. Article accepted by Zoologica Scripta. URL http://dx.doi.org/10.1111/j.1463-6409.2008.00361.x). This index was improved here to handle the incomplete taxonomic overlap among datasets. The results lead to the identification of new reliable clades within the 'acanthomorph bush'. Within a clade containing the Atherinomorpha, the Mugiloidei, the Plesiopidae, the Blennioidei, the Gobiesocoidei, the Cichlidae and the Pomacentridae, the Plesiopidae is the sister-group of the Mugiloidei. The Apogonidae are closely related to the Gobioidei. A clade named 'H' grouping a number of families close to stromateids and scombrids (Stromateidae, Scombridae, Trichiuridae, Chiasmodontidae, Nomeidae, Bramidae, Centrolophidae) is related to another clade named 'E' (Aulostomidae, Macrorhamphosidae, Dactylopteridae). The Sciaenidae is closely related to the Haemulidae. Within clade 'X' (Dettaï, A., Lecointre, G., 2004. In search of nothothenioid (Teleostei) relatives. Antarct. Sci. 16 (1), 71-85. URL http://dx.doi.org/10.1017/S0954102004), the Cottoidei, the Zoarcoidei, the Gasterosteidae, the Triglidae, the Scorpaenidae, the Sebastidae, the Synanceiidae, and the Congiopodidae form a clade. Within clade 'L' (Chen, W.-J., Bonillo, C., Lecointre, G., 2003. Repeatability of clades as a criterion of reliability: a case study for molecular phylogeny of Acanthomorpha (Teleostei) with larger number of taxa. Mol. Phylogenet. Evol. 26, 262-288; Dettaï, A., Lecointre, G., 2004. In search of nothothenioid (Teleostei) relatives. Antarct. Sci. 16 (1), 71-85. URL http://dx.doi.org/10.1017/S0954102004) grouping carangoids with flatfishes and other families (Centropomidae, Menidae, Sphyraenidae, Polynemidae, Echeneidae, Toxotidae, Xiphiidae), carangids are the stem-group of echeneids and coryphaenids, and sphyraenids are the sister-group to the Carangoidei. The Howellidae, the Epigonidae and the Lateolabracidae are closely related. We propose names for most of the clades repeatedly found in acanthomorph phylogenetic studies of various teams of the past decade.