Genomic evidence for ameiotic evolution in the bdelloid rotifer Adineta vaga.

Loss of sexual reproduction is considered an evolutionary dead end for metazoans, but bdelloid rotifers challenge this view as they appear to have persisted asexually for millions of years. Neither male sex organs nor meiosis have ever been observed in these microscopic animals: oocytes are formed through mitotic divisions, with no reduction of chromosome number and no indication of chromosome pairing. However, current evidence does not exclude that they may engage in sex on rare, cryptic occasions. Here we report the genome of a bdelloid rotifer, Adineta vaga (Davis, 1873), and show that its structure is incompatible with conventional meiosis. At gene scale, the genome of A. vaga is tetraploid and comprises both anciently duplicated segments and less divergent allelic regions. However, in contrast to sexual species, the allelic regions are rearranged and sometimes even found on the same chromosome. Such structure does not allow meiotic pairing; instead, we find abundant evidence of gene conversion, which may limit the accumulation of deleterious mutations in the absence of meiosis. Gene families involved in resistance to oxidation, carbohydrate metabolism and defence against transposons are significantly expanded, which may explain why transposable elements cover only 3% of the assembled sequence. Furthermore, 8% of the genes are likely to be of non-metazoan origin and were probably acquired horizontally. This apparent convergence between bdelloids and prokaryotes sheds new light on the evolutionary significance of sex.

The genome of the white-rot fungus Pycnoporus cinnabarinus: a basidiomycete model with a versatile arsenal for lignocellulosic biomass breakdown.

BACKGROUND: Saprophytic filamentous fungi are ubiquitous micro-organisms that play an essential role in photosynthetic carbon recycling. The wood-decayer Pycnoporus cinnabarinus is a model fungus for the study of plant cell wall decomposition and is used for a number of applications in green and white biotechnology. RESULTS: The 33.6 megabase genome of P. cinnabarinus was sequenced and assembled, and the 10,442 predicted genes were functionally annotated using a phylogenomic procedure. In-depth analyses were carried out for the numerous enzyme families involved in lignocellulosic biomass breakdown, for protein secretion and glycosylation pathways, and for mating type. The P. cinnabarinus genome sequence revealed a consistent repertoire of genes shared with wood-decaying basidiomycetes. P. cinnabarinus is thus fully equipped with the classical families involved in cellulose and hemicellulose degradation, whereas its pectinolytic repertoire appears relatively limited. In addition, P. cinnabarinus possesses a complete versatile enzymatic arsenal for lignin breakdown. We identified several genes encoding members of the three ligninolytic peroxidase types, namely lignin peroxidase, manganese peroxidase and versatile peroxidase. Comparative genome analyses were performed in fungi displaying different nutritional strategies (white-rot and brown-rot modes of decay). P. cinnabarinus presents a typical distribution of all the specific families found in the white-rot life style. Growth profiling of P. cinnabarinus was performed on 35 carbon sources including simple and complex substrates to study substrate utilization and preferences. P. cinnabarinus grew faster on crude plant substrates than on pure, mono- or polysaccharide substrates. Finally, proteomic analyses were conducted from liquid and solid-state fermentation to analyze the composition of the secretomes corresponding to growth on different substrates. The distribution of lignocellulolytic enzymes in the secretomes was strongly dependent on growth conditions, especially for lytic polysaccharide mono-oxygenases. CONCLUSIONS: With its available genome sequence, P. cinnabarinus is now an outstanding model system for the study of the enzyme machinery involved in the degradation or transformation of lignocellulosic biomass.

An automated approach for the identification of horizontal gene transfers from complete genomes reveals the rhizome of Rickettsiales.

BACKGROUND: Horizontal gene transfer (HGT) is considered to be a major force driving the evolutionary history of prokaryotes. HGT is widespread in prokaryotes, contributing to the genomic repertoire of prokaryotic organisms, and is particularly apparent in Rickettsiales genomes. Gene gains from both distantly and closely related organisms play crucial roles in the evolution of bacterial genomes. In this work, we focus on genes transferred from distantly related species into Rickettsiales species. RESULTS: We developed an automated approach for the detection of HGT from other organisms (excluding alphaproteobacteria) into Rickettsiales genomes. Our systematic approach consisted of several specialized features including the application of a parsimony method for inferring phyletic patterns followed by blast filter, automated phylogenetic reconstruction and the application of patterns for HGT detection. We identified 42 instances of HGT in 31 complete Rickettsiales genomes, of which 38 were previously unidentified instances of HGT from Anaplasma, Wolbachia, Candidatus Pelagibacter ubique and Rickettsia genomes. Additionally, putative cases with no phylogenetic support were assigned gene ontology terms. Overall, these transfers could be characterized as "rhizome-like". CONCLUSIONS: Our analysis provides a comprehensive, systematic approach for the automated detection of HGTs from several complete proteome sequences that can be applied to detect instances of HGT within other genomes of interest.

The butyrophilin (BTN) gene family: from milk fat to the regulation of the immune response.

Butyrophilins (BTN) belong to the immunoglobulin (Ig) superfamily of transmembrane proteins. These molecules are of increasing interest to immunologists, as they share a structural homology with B7 family members at the extracellular domain level. Moreover, a role of these molecules has been suggested in the negative regulation of lymphocyte activation for almost all the BTN that have been studied. In addition, the expression of some BTN family members has been reported to be associated with autoimmune diseases. Over the last few years, the number of BTN and BTN-like members has greatly increased. In this study, the butyrophilin family in mammals has been revisited, using phylogenetic analysis to identify all the family members and the phylogenetic relations among them, and to establish a standard nomenclature. Fourteen BTN groups were identified that are not all conserved between mammalian species. In addition, an overview of expression profiles and functional BTN data demonstrates that these molecules represent a new area of investigation for the design of future strategies in the modulation of the immune system.

Peptidoglycan: a post-genomic analysis.

BACKGROUND: To derive post-genomic, neutral insight into the peptidoglycan (PG) distribution among organisms, we mined 1,644 genomes listed in the Carbohydrate-Active Enzymes database for the presence of a minimal 3-gene set that is necessary for PG metabolism. This gene set consists of one gene from the glycosyltransferase family GT28, one from family GT51 and at least one gene belonging to one of five glycoside hydrolase families (GH23, GH73, GH102, GH103 and GH104). RESULTS: None of the 103 Viruses or 101 Archaea examined possessed the minimal 3-gene set, but this set was detected in 1/42 of the Eukarya members (Micromonas sp., coding for GT28, GT51 and GH103) and in 1,260/1,398 (90.1%) of Bacteria, with a 100% positive predictive value for the presence of PG. Pearson correlation test showed that GT51 family genes were significantly associated with PG with a value of 0.963 and a p value less than 10(-3). This result was confirmed by a phylogenetic comparative analysis showing that the GT51-encoding gene was significantly associated with PG with a Pagel's score of 60 and 51 (percentage of error close to 0%). Phylogenetic analysis indicated that the GT51 gene history comprised eight loss and one gain events, and suggested a dynamic on-going process. CONCLUSIONS: Genome analysis is a neutral approach to explore prospectively the presence of PG in uncultured, sequenced organisms with high predictive values.

Validation of new HLA-J alleles assigned by next generation sequencing.

We describe nine novel HLA-J alleles validated by in silico and experimental data.

Identification of 11 novel HLA-A, -B, -C, -DRB1, and -DQB1 alleles in the 1000 Genomes Project panel.

We isolated 11 novel HLA alleles in the 1000 Genomes Project using PolyPheMe software.

Blood group typing from whole-genome sequencing data.

Many questions can be explored thanks to whole-genome data. The aim of this study was to overcome their main limits, software availability and database accuracy, and estimate the feasibility of red blood cell (RBC) antigen typing from whole-genome sequencing (WGS) data. We analyzed whole-genome data from 79 individuals for HLA-DRB1 and 9 RBC antigens. Whole-genome sequencing data was analyzed with software allowing phasing of variable positions to define alleles or haplotypes and validated for HLA typing from next-generation sequencing data. A dedicated database was set up with 1648 variable positions analyzed in KEL (KEL), ACKR1 (FY), SLC14A1 (JK), ACHE (YT), ART4 (DO), AQP1 (CO), CD44 (IN), SLC4A1 (DI) and ICAM4 (LW). Whole-genome sequencing typing was compared to that previously obtained by amplicon-based monoallelic sequencing and by SNaPshot analysis. Whole-genome sequencing data were also explored for other alleles. Our results showed 93% of concordance for blood group polymorphisms and 91% for HLA-DRB1. Incorrect typing and unresolved results confirm that WGS should be considered reliable with read depths strictly above 15x. Our results supported that RBC antigen typing from WGS is feasible but requires improvements in read depth for SNV polymorphisms typing accuracy. We also showed the potential for WGS in screening donors with rare blood antigens, such as weak JK alleles. The development of WGS analysis in immunogenetics laboratories would offer personalized care in the management of RBC disorders.

PhyloPattern: regular expressions to identify complex patterns in phylogenetic trees.

BACKGROUND: To effectively apply evolutionary concepts in genome-scale studies, large numbers of phylogenetic trees have to be automatically analysed, at a level approaching human expertise. Complex architectures must be recognized within the trees, so that associated information can be extracted. RESULTS: Here, we present a new software library, PhyloPattern, for automating tree manipulations and analysis. PhyloPattern includes three main modules, which address essential tasks in high-throughput phylogenetic tree analysis: node annotation, pattern matching, and tree comparison. PhyloPattern thus allows the programmer to focus on: i) the use of predefined or user defined annotation functions to perform immediate or deferred evaluation of node properties, ii) the search for user-defined patterns in large phylogenetic trees, iii) the pairwise comparison of trees by dynamically generating patterns from one tree and applying them to the other. CONCLUSION: PhyloPattern greatly simplifies and accelerates the work of the computer scientist in the evolutionary biology field. The library has been used to automatically identify phylogenetic evidence for domain shuffling or gene loss events in the evolutionary histories of protein sequences. However any workflow that relies on phylogenetic tree analysis, could be automated with PhyloPattern.

CASSIOPE: an expert system for conserved regions searches.

BACKGROUND: Understanding genome evolution provides insight into biological mechanisms. For many years comparative genomics and analysis of conserved chromosomal regions have helped to unravel the mechanisms involved in genome evolution and their implications for the study of biological systems. Detection of conserved regions (descending from a common ancestor) not only helps clarify genome evolution but also makes it possible to identify quantitative trait loci (QTLs) and investigate gene function.The identification and comparison of conserved regions on a genome scale is computationally intensive, making process automation essential. Three key requirements are necessary: consideration of phylogeny to identify orthologs between multiple species, frequent updating of the annotation and panel of compared genomes and computation of statistical tests to assess the significance of identified conserved gene clusters. RESULTS: We developed a modular system superimposed on a multi-agent framework, called CASSIOPE (Clever Agent System for Synteny Inheritance and Other Phenomena in Evolution). CASSIOPE automatically identifies statistically significant conserved regions between multiple genomes based on automated phylogenies and statistical testing. Conserved regions were searched for in 19 species and 1,561 hits were found. To our knowledge, CASSIOPE is the first system to date that integrates evolutionary biology-based concepts and fulfills all three key requirements stated above. All results are available at http://194.57.197.245/cassiopeWeb/displayCluster?clusterId=1 CONCLUSION: CASSIOPE makes it possible to study conserved regions from a chosen query genetic region and to infer conserved gene clusters based on phylogenies and statistical tests assessing the significance of these conserved regions.Source code is freely available, please contact: Pierre.pontarotti@univ-provence.fr.

HLAIb worldwide genetic diversity: New HLA-H alleles and haplotype structure description.

The classical HLA class I genes (HLA Ia) were extensively studied because of their implication in clinical fields and anthropology. Less is known about worldwide genetic diversity and linkage disequilibrium for non-classical HLA class I genes (HLA Ib) and HLA pseudogenes. Notably, HLA-H, which is deleted in a fraction of the population, remains scarcely explored. The aims of this study were 1/ to get further insight into HLA-H genetic diversity and into how this variability potentially affects its expression and 2/ to define HLA Ib worldwide allelic diversity and linkage. Exome sequence data from the 1000 Genomes Project were used to define second field HLA-A, -E, -F, -G and -H typing using PolyPheMe software. Allelic and two-loci haplotype frequencies were estimated using Gene[Rate] software both at worldwide and continental levels. Eleven novel HLA-H alleles identified in exome data were validated by NGS performed on 25 genomic DNA samples from the same cohort. Phylogenetic analysis and frequency distribution of HLA-H alleles revealed three clades, each predominantly represented in Admixed American, European and East Asian populations, African populations and South Asian populations. Among these eleven novel alleles, two potentially encode complete transmembrane HLA proteins. We confirm the high LD between HLA-H and -A, and between HLA-H and -G, and show the three genes have distinct worldwide allelic distribution. Conversely, HLA-E and HLA-F both showed little LD, displayed restricted allelic diversity and practically no difference in their distribution across the planet. Our work thus reveals an unexpectedly high HLA-H genetic diversity, with alleles highly represented in Asia possibly encoding a functional HLA protein. Functional implication of these results remains to be explored, both in physiological and pathological contexts.

Immune diversity sheds light on missing variation in worldwide genetic diversity panels.

Defining worldwide human genetic variation is a critical step to reveal how genome plasticity contributes to disease. Yet, there is currently no metric to assess the representativeness and completeness of current and widely used data on genetic variation. We show here that Human Leukocyte Antigen (HLA) genes can serve as such metric as they are both the most polymorphic and the most studied genetic system. As a test case, we investigated the 1,000 Genomes Project panel. Using high-accuracy in silico HLA typing, we find that over 20% of the common HLA variants and over 70% of the rare HLA variants are missing in this reference panel for worldwide genetic variation, due to undersampling and incomplete geographical coverage, in particular in Oceania and West Asia. Because common and rare variants both contribute to disease, this study thus illustrates how HLA diversity can detect and help fix incomplete sampling and hence accelerate efforts to draw a comprehensive overview of the genetic variation that is relevant to health and disease.

HLA-G Haplotypes Are Differentially Associated with Asthmatic Features.

Human leukocyte antigen (HLA)-G, a HLA class Ib molecule, interacts with receptors on lymphocytes such as T cells, B cells, and natural killer cells to influence immune responses. Unlike classical HLA molecules, HLA-G expression is not found on all somatic cells, but restricted to tissue sites, including human bronchial epithelium cells (HBEC). Individual variation in HLA-G expression is linked to its genetic polymorphism and has been associated with many pathological situations such as asthma, which is characterized by epithelium abnormalities and inflammatory cell activation. Studies reported both higher and equivalent soluble HLA-G (sHLA-G) expression in different cohorts of asthmatic patients. In particular, we recently described impaired local expression of HLA-G and abnormal profiles for alternatively spliced isoforms in HBEC from asthmatic patients. sHLA-G dosage is challenging because of its many levels of polymorphism (dimerization, association with ß2-microglobulin, and alternative splicing), thus many clinical studies focused on HLA-G single-nucleotide polymorphisms as predictive biomarkers, but few analyzed HLA-G haplotypes. Here, we aimed to characterize HLA-G haplotypes and describe their association with asthmatic clinical features and sHLA-G peripheral expression and to describe variations in transcription factor (TF) binding sites and alternative splicing sites. HLA-G haplotypes were differentially distributed in 330 healthy and 580 asthmatic individuals. Furthermore, HLA-G haplotypes were associated with asthmatic clinical features showed. However, we did not confirm an association between sHLA-G and genetic, biological, or clinical parameters. HLA-G haplotypes were phylogenetically split into distinct groups, with each group displaying particular variations in TF binding or RNA splicing sites that could reflect differential HLA-G qualitative or quantitative expression, with tissue-dependent specificities. Our results, based on a multicenter cohort, thus support the pertinence of HLA-G haplotypes as predictive genetic markers for asthma.

Eleven ancestral gene families lost in mammals and vertebrates while otherwise universally conserved in animals.

BACKGROUND: Gene losses played a role which may have been as important as gene and genome duplications and rearrangements, in modelling today species' genomes from a common ancestral set of genes. The set and diversity of protein-coding genes in a species has direct output at the functional level. While gene losses have been reported in all the major lineages of the metazoan tree of life, none have proposed a focus on specific losses in the vertebrates and mammals lineages. In contrast, genes lost in protostomes (i.e. arthropods and nematodes) but still present in vertebrates have been reported and extensively detailed. This probable over-anthropocentric way of comparing genomes does not consider as an important phenomena, gene losses in species that are usually described as "higher". However reporting universally conserved genes throughout evolution that have recently been lost in vertebrates and mammals could reveal interesting features about the evolution of our genome, particularly if these losses can be related to losses of capability. RESULTS: We report 11 gene families conserved throughout eukaryotes from yeasts (such as Saccharomyces cerevisiae) to bilaterian animals (such as Drosophila melanogaster or Caenorhabditis elegans). This evolutionarily wide conservation suggests they were present in the last common ancestors of fungi and metazoan animals. None of these 11 gene families are found in human nor mouse genomes, and their absence generally extends to all vertebrates. A total of 8 out of these 11 gene families have orthologs in plants, suggesting they were present in the Last Eukaryotic Common Ancestor (LECA). We investigated known functional information for these 11 gene families. This allowed us to correlate some of the lost gene families to loss of capabilities. CONCLUSION: Mammalian and vertebrate genomes lost evolutionary conserved ancestral genes that are probably otherwise not dispensable in eukaryotes. Hence, the human genome, which is generally viewed as being the result of increased complexity and gene-content, has also evolved through simplification and gene losses. This acknowledgement confirms, as already suggested, that the genome of our far ancestor was probably more complex than ever considered.

Tracking the connection between evolutionary and functional shifts using the fungal lipase/feruloyl esterase A family.

BACKGROUND: There have been many claims of adaptive molecular evolution, but what role does positive selection play in functional divergence? The aim of this study was to test the relationship between evolutionary and functional shifts with special emphasis on the role of the environment. For this purpose, we studied the fungal lipase/feruloyl esterase A family, whose functional diversification makes it a very promising candidate. RESULTS: The results suggested functional shift following a duplication event where neofunctionalisation of feruloyl esterase A had occurred with conservation of the ancestral lipase function. Evolutionary shift was detected using the branch-site model for testing positive selection on individual codons along specific lineages. Positively selected amino acids were detected. Furthermore, biological data obtained from site-directed mutagenesis experiments clearly demonstrated that certain amino acids under positive selection were involved in the functional shift. We reassessed evolutionary history in terms of environmental response, and hypothesized that environmental changes such as colonisation by terrestrial plants might have driven adaptation by functional diversification in Euascomycetes (Aspergilli), thus conferring a selective advantage on this group. CONCLUSION: The results reported here illustrate a rare example of connection between fundamental events in molecular evolution. We demonstrated an unequivocal connection between evolutionary and functional shifts, which led us to conclude that these events were probably linked to environmental change.

FIGENIX: intelligent automation of genomic annotation: expertise integration in a new software platform.

BACKGROUND: Two of the main objectives of the genomic and post-genomic era are to structurally and functionally annotate genomes which consists of detecting genes' position and structure, and inferring their function (as well as of other features of genomes). Structural and functional annotation both require the complex chaining of numerous different software, algorithms and methods under the supervision of a biologist. The automation of these pipelines is necessary to manage huge amounts of data released by sequencing projects. Several pipelines already automate some of these complex chaining but still necessitate an important contribution of biologists for supervising and controlling the results at various steps. RESULTS: Here we propose an innovative automated platform, FIGENIX, which includes an expert system capable to substitute to human expertise at several key steps. FIGENIX currently automates complex pipelines of structural and functional annotation under the supervision of the expert system (which allows for example to make key decisions, check intermediate results or refine the dataset). The quality of the results produced by FIGENIX is comparable to those obtained by expert biologists with a drastic gain in terms of time costs and avoidance of errors due to the human manipulation of data. CONCLUSION: The core engine and expert system of the FIGENIX platform currently handle complex annotation processes of broad interest for the genomic community. They could be easily adapted to new, or more specialized pipelines, such as for example the annotation of miRNAs, the classification of complex multigenic families, annotation of regulatory elements and other genomic features of interest.

Mycobacteriophage-drived diversification of Mycobacterium abscessus.

BACKGROUND: Mycobacterium abscessus is an emerging opportunistic pathogen which diversity was acknowledged by the recent description of two subspecies accommodating M. abscessus, Mycobacterium bolletii and Mycobacterium massiliense isolates. RESULTS: Here, genome analysis found 1-8 prophage regions in 47/48 M. abscessus genomes ranging from small prophage-like elements to complete prophages. A total of 20,304 viral and phage proteins clustered into 853 orthologous groups. Phylogenomic and phylogenetic analyses based on prophage region homology found three main clusters corresponding to M. abscessus, M. bolletii and M. massiliense. Analysing 135 annotated Tape Measure Proteins found thirteen clusters and four singletons, suggesting that at least 17 mycobacteriophages had infected M. abscessus during its evolution. The evolutionary history of phages differed from that of their mycobacterial hosts. In particular, 33 phage-related proteins have been horizontally transferred within M. abscessus genomes. They comprise of an integrase, specific mycobacteriophage proteins, hypothetical proteins and DNA replication and metabolism proteins. Gene exchanges, loss and gains which occurred in M. abscessus genomes have been driven by several mycobacteriophages. CONCLUSIONS: This analysis of phage-mycobacterium co-evolution suggests that mycobacteriophages are playing a key-role in the on-going diversification of M. abscessus. REVIEWERS: This article was reviewed by Eric Bapteste, Patrick Forterre and Eugene Koonin.

Reliable Phylogenetic Trees Building: A New Web Interface for FIGENIX.

The community needed a reliable and user friendly tool to quickly produce robust phylogenetic trees which are crucial in evolutionary studies and genomes' functional annotation. FIGENIX is software dedicated to this and was published in 2005. Several laboratories around the world use it in their research, but it was difficult to use for non-expert users, thus we developed a new graphical user interface for the benefit of all biologists.

GLADX: an automated approach to analyze the lineage-specific loss and pseudogenization of genes.

A well-established ancestral gene can usually be found, in one or multiple copies, in different descendant species. Sometimes during the course of evolution, all the representatives of a well-established ancestral gene disappear in specific lineages; such gene losses may occur in the genome by deletion of a DNA fragment or by pseudogenization. The loss of an entire gene family in a given lineage may reflect an important phenomenon, and could be due either to adaptation, or to a relaxation of selection that leads to neutral evolution. Therefore, the lineage-specific gene loss analyses are important to improve the understanding of the evolutionary history of genes and genomes. In order to perform this kind of study from the increasing number of complete genome sequences available, we developed a unique new software module called GLADX in the DAGOBAH framework, based on a comparative genomic approach. The software is able to automatically detect, for all the species of a phylum, the presence/absence of a representative of a well-established ancestral gene, and by systematic steps of re-annotation, confirm losses, detect and analyze pseudogenes and find novel genes. The approach is based on the use of highly reliable gene phylogenies, of protein predictions and on the analysis of genomic mutations. All the evidence associated to evolutionary approach provides accurate information for building an overall view of the evolution of a given gene in a selected phylum. The reliability of GLADX has been successfully tested on a benchmark analysis of 14 reported cases. It is the first tool that is able to fully automatically study the lineage-specific losses and pseudogenizations. GLADX is available at http://ioda.univ-provence.fr/IodaSite/gladx/.