PhylomeDB v3.0: an expanding repository of genome-wide collections of trees, alignments and phylogeny-based orthology and paralogy predictions

Several secondary metabolite gene clusters in the genomes of ten Penicillium spp. raise the risk of multiple mycotoxin occurrence in chestnuts

Penicillium spp. produce a great variety of secondary metabolites, including several mycotoxins, on food substrates. Chestnuts represent a favorable substrate for Penicillium spp. development. In this study, the genomes of ten Penicillium species, virulent on chestnuts, were sequenced and annotated: P. bialowiezense. P. pancosmium, P. manginii, P. discolor, P. crustosum, P. palitans, P. viridicatum, P. glandicola, P. taurinense and P. terrarumae. Assembly size ranges from 27.5 to 36.8 Mb and the number of encoded genes ranges from 9,867 to 12,520. The total number of predicted biosynthetic gene clusters (BGCs) in the ten species is 551. The most represented families of BGCs are non ribosomal peptide synthase (191) and polyketide synthase (175), followed by terpene synthases (87). Genome-wide collections of gene phylogenies (phylomes) were reconstructed for each of the newly sequenced Penicillium species allowing for the prediction of orthologous relationships among our species, as well as other 20 annotated Penicillium species available in the public domain. We investigated in silico the presence of BGCs for 10 secondary metabolites, including 5 mycotoxins, whose production was validated in vivo through chemical analyses. Among the clusters present in this set of species we found andrastin A and its related cluster atlantinone A, mycophenolic acid, patulin, penitrem A and the cluster responsible for the synthesis of roquefortine C/glandicoline A/glandicoline B/meleagrin. We confirmed the presence of these clusters in several of the Penicillium species conforming our dataset and verified their capacity to synthesize them in a chestnut-based medium with chemical analysis. Interestingly, we identified mycotoxin clusters in some species for the first time, such as the andrastin A cluster in P. flavigenum and P. taurinense, and the roquefortine C cluster in P. nalgiovense and P. taurinense. Chestnuts proved to be an optimal substrate for species of Penicillium with different mycotoxigenic potential, opening the door to risks related to the occurrence of multiple mycotoxins in the same food matrix.

Genomic, transcriptomic, and metabolomic analysis of Oldenlandia corymbosa reveals the biosynthesis and mode of action of anti-cancer metabolites

Plants accumulate a vast array of secondary metabolites, which constitute a natural resource for pharmaceuticals. Oldenlandia corymbosa belongs to the Rubiaceae family, and has been used in traditional medicine to treat different diseases, including cancer. However, the active metabolites of the plant, their biosynthetic pathway and mode of action in cancer are unknown. To fill these gaps, we exposed this plant to eight different stress conditions and combined different omics data capturing gene expression, metabolic profiles, and anti-cancer activity. Our results show that O. corymbosa extracts are active against breast cancer cell lines and that ursolic acid is responsible for this activity. Moreover, we assembled a high-quality genome and uncovered two genes involved in the biosynthesis of ursolic acid. Finally, we also revealed that ursolic acid causes mitotic catastrophe in cancer cells and identified three high-confidence protein binding targets by Cellular Thermal Shift Assay (CETSA) and reverse docking. Altogether, these results constitute a valuable resource to further characterize the biosynthesis of active metabolites in the Oldenlandia group, while the mode of action of ursolic acid will allow us to further develop this valuable compound.

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

Evolutionary analyses of genes in Echinodermata offer insights towards the origin of metazoan phyla

Despite recent studies discussing the evolutionary impacts of gene duplications and losses among metazoans, the genomic basis for the evolution of phyla remains enigmatic. Here, we employ phylogenomic approaches to search for orthologous genes without known functions among echinoderms, and subsequently use them to guide the identification of their homologs across other metazoans. Our final set of 14 genes was obtained via a suite of homology prediction tools, gene expression data, gene ontology, and generating the Strongylocentrotus purpuratus phylome. The gene set was subjected to selection pressure analyses, which indicated that they are highly conserved and under negative selection. Their presence across broad taxonomic depths suggests that genes required to form a phylum are ancestral to that phylum. Therefore, rather than de novo gene genesis, we posit that evolutionary forces such as selection on existing genomic elements over large timescales may drive divergence and contribute to the emergence of phyla.

Candida glabrata Hst1-Rfm1-Sum1 complex evolved to control virulence-related genes

C. glabrata is an opportunistic fungal pathogen and the second most common cause of opportunistic fungal infections in humans, that has evolved virulence factors to become a successful pathogen: strong resistance to oxidative stress, capable to adhere and form biofilms in human epithelial cells as well as to abiotic surfaces and high resistance to xenobiotics. Hst1 (a NAD⁺-dependent histone deacetylase), Sum1 (putative DNA binding protein) and Rfm1 (connector protein) form a complex (HRS-C) and control the resistance to oxidative stress, to xenobiotics (the antifungal fluconazole), and adherence to epithelial cells. Hst1 is functionally conserved within the Saccharomycetaceae family, Rfm1 shows a close phylogenetic relation within the Saccharomycetaceae family while Sum1 displays a distant phylogenetic relation with members of the family and is not conserved functionally. CDR1 encodes for an ABC transporter (resistance to fluconazole) negatively controlled by HRS-C, for which its binding site is located within 223 bp upstream from the ATG of CDR1. The absence of Hst1 and Sum1 renders the cells hyper-adherent, possibly due to the overexpression of AED1, EPA1, EPA22 and EPA6, all encoding for adhesins. Finally, in a neutrophil survival assay, HST1 and SUM1, are not required for survival. We propose that Sum1 in the HRS-C diverged functionally to control a set of genes implicated in virulence: adherence, resistance to xenobiotics and oxidative stress.

Comparative transcriptomic analysis reveals conserved programmes underpinning organogenesis and reproduction in land plants

The appearance of plant organs mediated the explosive radiation of land plants, which shaped the biosphere and allowed the establishment of terrestrial animal life. The evolution of organs and immobile gametes required the coordinated acquisition of novel gene functions, the co-option of existing genes and the development of novel regulatory programmes. However, no large-scale analyses of genomic and transcriptomic data have been performed for land plants. To remedy this, we generated gene expression atlases for various organs and gametes of ten plant species comprising bryophytes, vascular plants, gymnosperms and flowering plants. A comparative analysis of the atlases identified hundreds of organ- and gamete-specific orthogroups and revealed that most of the specific transcriptomes are significantly conserved. Interestingly, our results suggest that co-option of existing genes is the main mechanism for evolving new organs. In contrast to female gametes, male gametes showed a high number and conservation of specific genes, which indicates that male reproduction is highly specialized. The expression atlas capturing pollen development revealed numerous transcription factors and kinases essential for pollen biogenesis and function.

Comparative transcriptomic analysis reveals conserved programmes underpinning organogenesis and reproduction in land plants

The appearance of plant organs mediated the explosive radiation of land plants, which shaped the biosphere and allowed the establishment of terrestrial animal life. The evolution of organs and immobile gametes required the coordinated acquisition of novel gene functions, the co-option of existing genes and the development of novel regulatory programmes. However, no large-scale analyses of genomic and transcriptomic data have been performed for land plants. To remedy this, we generated gene expression atlases for various organs and gametes of ten plant species comprising bryophytes, vascular plants, gymnosperms and flowering plants. A comparative analysis of the atlases identified hundreds of organ- and gamete-specific orthogroups and revealed that most of the specific transcriptomes are significantly conserved. Interestingly, our results suggest that co-option of existing genes is the main mechanism for evolving new organs. In contrast to female gametes, male gametes showed a high number and conservation of specific genes, which indicates that male reproduction is highly specialized. The expression atlas capturing pollen development revealed numerous transcription factors and kinases essential for pollen biogenesis and function.

Comparative Genomics Used to Predict Virulence Factors and Metabolic Genes among Monilinia Species

Brown rot, caused by Monilinia spp., is among the most important diseases in stone fruits, and some pome fruits (mainly apples). This disease is responsible for significant yield losses, particularly in stone fruits, when weather conditions favorable for disease development appear. To achieve future sustainable strategies to control brown rot on fruit, one potential approach will be to characterize genomic variation among Monilinia spp. to define, among others, the capacity to infect fruit in this genus. In the present work, we performed genomic and phylogenomic comparisons of five Monilinia species and inferred differences in numbers of secreted proteins, including CAZy proteins and other proteins important for virulence. Duplications specific to Monilinia were sparse and, overall, more genes have been lost than gained. Among Monilinia spp., low variability in the CAZome was observed. Interestingly, we identified several secondary metabolism clusters based on similarity to known clusters, and among them was a cluster with homology to pyriculol that could be responsible for the synthesis of chloromonilicin. Furthermore, we compared sequences of all strains available from NCBI of these species to assess their MAT loci and heterokaryon compatibility systems. Our comparative analyses provide the basis for future studies into understanding how these genomic differences underlie common or differential abilities to interact with the host plant.

Predicting the Evolution of Syntenies—An Algorithmic Review

Syntenies are genomic segments of consecutive genes identified by a certain conservation in gene content and order. The notion of conservation may vary from one definition to another, the more constrained requiring identical gene contents and gene orders, while more relaxed definitions just require a certain similarity in gene content, and not necessarily in the same order. Regardless of the way they are identified, the goal is to characterize homologous genomic regions, i.e., regions deriving from a common ancestral region, reflecting a certain gene co-evolution that can enlighten important functional properties. In addition of being able to identify them, it is also necessary to infer the evolutionary history that has led from the ancestral segment to the extant ones. In this field, most algorithmic studies address the problem of inferring rearrangement scenarios explaining the disruption in gene order between segments with the same gene content, some of them extending the evolutionary model to gene insertion and deletion. However, syntenies also evolve through other events modifying their content in genes, such as duplications, losses or horizontal gene transfers, i.e., the movement of genes from one species to another. Although the reconciliation approach between a gene tree and a species tree addresses the problem of inferring such events for single-gene families, little effort has been dedicated to the generalization to segmental events and to syntenies. This paper reviews some of the main algorithmic methods for inferring ancestral syntenies and focus on those integrating both gene orders and gene trees.

Selection following Gene Duplication Shapes Recent Genome Evolution in the Pea Aphid Acyrthosiphon pisum

Ecology of insects is as wide as their diversity, which reflects their high capacity of adaptation in most of the environments of our planet. Aphids, with over 4,000 species, have developed a series of adaptations including a high phenotypic plasticity and the ability to feed on the phloem sap of plants, which is enriched in sugars derived from photosynthesis. Recent analyses of aphid genomes have indicated a high level of shared ancestral gene duplications that might represent a basis for genetic innovation and broad adaptations. In addition, there are a large number of recent, species-specific gene duplications whose role in adaptation remains poorly understood. Here, we tested whether duplicates specific to the pea aphid Acyrthosiphon pisum are related to genomic innovation by combining comparative genomics, transcriptomics, and chromatin accessibility analyses. Consistent with large levels of neofunctionalization, we found that most of the recent pairs of gene duplicates evolved asymmetrically, showing divergent patterns of positive selection and gene expression. Genes under selection involved a plethora of biological functions, suggesting that neofunctionalization and tissue specificity, among other evolutionary mechanisms, have orchestrated the evolution of recent paralogs in the pea aphid and may have facilitated host-symbiont cooperation. Our comprehensive phylogenomics analysis allowed us to tackle the history of duplicated genes to pave the road toward understanding the role of gene duplication in ecological adaptation.

MetaPhOrs 2.0: integrative, phylogeny-based inference of orthology and paralogy across the tree of life

Inferring homology relationships across genes in different species is a central task in comparative genomics. Therefore, a large number of resources and methods have been developed over the years. Some public databases include phylogenetic trees of homologous gene families which can be used to further differentiate homology relationships into orthology and paralogy. MetaPhOrs is a web server that integrates phylogenetic information from different sources to provide orthology and paralogy relationships based on a common phylogeny-based predictive algorithm and associated with a consistency-based confidence score. Here we describe the latest version of the web server which includes major new implementations and provides orthology and paralogy relationships derived from ∼8.2 million gene family trees-from 13 different source repositories across ∼4000 species with sequenced genomes. MetaPhOrs server is freely available, without registration, at http://orthology.phylomedb.org/.

An Overview of Duplicated Gene Detection Methods: Why the Duplication Mechanism Has to Be Accounted for in Their Choice

Gene duplication is an important evolutionary mechanism allowing to provide new genetic material and thus opportunities to acquire new gene functions for an organism, with major implications such as speciation events. Various processes are known to allow a gene to be duplicated and different models explain how duplicated genes can be maintained in genomes. Due to their particular importance, the identification of duplicated genes is essential when studying genome evolution but it can still be a challenge due to the various fates duplicated genes can encounter. In this review, we first describe the evolutionary processes allowing the formation of duplicated genes but also describe the various bioinformatic approaches that can be used to identify them in genome sequences. Indeed, these bioinformatic approaches differ according to the underlying duplication mechanism. Hence, understanding the specificity of the duplicated genes of interest is a great asset for tool selection and should be taken into account when exploring a biological question.

eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses

eggNOG is a public database of orthology relationships, gene evolutionary histories and functional annotations. Here, we present version 5.0, featuring a major update of the underlying genome sets, which have been expanded to 4445 representative bacteria and 168 archaea derived from 25 038 genomes, as well as 477 eukaryotic organisms and 2502 viral proteomes that were selected for diversity and filtered by genome quality. In total, 4.4M orthologous groups (OGs) distributed across 379 taxonomic levels were computed together with their associated sequence alignments, phylogenies, HMM models and functional descriptors. Precomputed evolutionary analysis provides fine-grained resolution of duplication/speciation events within each OG. Our benchmarks show that, despite doubling the amount of genomes, the quality of orthology assignments and functional annotations (80% coverage) has persisted without significant changes across this update. Finally, we improved eggNOG online services for fast functional annotation and orthology prediction of custom genomics or metagenomics datasets. All precomputed data are publicly available for downloading or via API queries at http://eggnog.embl.de

Soybean aphid biotype 1 genome: Insights into the invasive biology and adaptive evolution of a major agricultural pest

The soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae) is a serious pest of the soybean plant, Glycine max, a major world-wide agricultural crop. We assembled a de novo genome sequence of Ap. glycines Biotype 1, from a culture established shortly after this species invaded North America. 20.4% of the Ap. glycines proteome is duplicated. These in-paralogs are enriched with Gene Ontology (GO) categories mostly related to apoptosis, a possible adaptation to plant chemistry and other environmental stressors. Approximately one-third of these genes show parallel duplication in other aphids. But Ap. gossypii, its closest related species, has the lowest number of these duplicated genes. An Illumina GoldenGate assay of 2380 SNPs was used to determine the world-wide population structure of Ap. Glycines. China and South Korean aphids are the closest to those in North America. China is the likely origin of other Asian aphid populations. The most distantly related aphids to those in North America are from Australia. The diversity of Ap. glycines in North America has decreased over time since its arrival. The genetic diversity of Ap. glycines North American population sampled shortly after its first detection in 2001 up to 2012 does not appear to correlate with geography. However, aphids collected on soybean Rag experimental varieties in Minnesota (MN), Iowa (IA), and Wisconsin (WI), closer to high density Rhamnus cathartica stands, appear to have higher capacity to colonize resistant soybean plants than aphids sampled in Ohio (OH), North Dakota (ND), and South Dakota (SD). Samples from the former states have SNP alleles with high F_ST values and frequencies, that overlap with genes involved in iron metabolism, a crucial metabolic pathway that may be affected by the Rag-associated soybean plant response. The Ap. glycines Biotype 1 genome will provide needed information for future analyses of mechanisms of aphid virulence and pesticide resistance as well as facilitate comparative analyses between aphids with differing natural history and host plant range.

Phylogenomics Identifies an Ancestral Burst of Gene Duplications Predating the Diversification of Aphidomorpha

Aphids (Aphidoidea) are a diverse group of hemipteran insects that feed on plant phloem sap. A common finding in studies of aphid genomes is the presence of a large number of duplicated genes. However, when these duplications occurred remains unclear, partly due to the high relatedness of sequenced species. To better understand the origin of aphid duplications we sequenced and assembled the genome of Cinara cedri, an early branching lineage (Lachninae) of the Aphididae family. We performed a phylogenomic comparison of this genome with 20 other sequenced genomes, including the available genomes of five other aphids, along with the transcriptomes of two species belonging to Adelgidae (a closely related clade to the aphids) and Coccoidea. We found that gene duplication has been pervasive throughout the evolution of aphids, including many parallel waves of recent, species-specific duplications. Most notably, we identified a consistent set of very ancestral duplications, originating from a large-scale gene duplication predating the diversification of Aphidomorpha (comprising aphids, phylloxerids, and adelgids). Genes duplicated in this ancestral wave are enriched in functions related to traits shared by Aphidomorpha, such as association with endosymbionts, and adaptation to plant defenses and phloem-sap-based diet. The ancestral nature of this duplication wave (106-227 Ma) and the lack of sufficiently conserved synteny make it difficult to conclude whether it originated from a whole-genome duplication event or, alternatively, from a burst of large-scale segmental duplications. Genome sequencing of other aphid species belonging to different Aphidomorpha and related lineages may clarify these findings.

Comparative Genomics and Evolutionary Analysis of RNA-Binding Proteins of Burkholderia cenocepacia J2315 and Other Members of the B. cepacia Complex

RNA-binding proteins (RBPs) are important regulators of cellular functions, playing critical roles on the survival of bacteria and in the case of pathogens, on their interaction with the host. RBPs are involved in transcriptional, post-transcriptional, and translational processes. However, except for model organisms like Escherichia coli, there is little information about the identification or characterization of RBPs in other bacteria, namely in members of the Burkholderia cepacia complex (Bcc). Bcc is a group of bacterial species associated with a poor clinical prognosis in cystic fibrosis patients. These species have some of the largest bacterial genomes, and except for the presence of two-distinct Hfq-like proteins, their RBP repertoire has not been analyzed so far. Using in silico approaches, we identified 186 conventional putative RBPs in Burkholderia cenocepacia J2315, an epidemic and multidrug resistant pathogen of cystic fibrosis patients. Here we describe the comparative genomics and phylogenetic analysis of RBPs present in multiple copies and predicted to play a role in transcription, protein synthesis, and RNA decay in Bcc bacteria. In addition to the two different Hfq chaperones, five cold shock proteins phylogenetically close to E. coli CspD protein and three distinct RhlE-like helicases could be found in the B. cenocepacia J2315 genome. No RhlB, SrmB, or DeaD helicases could be found in the genomes of these bacteria. These results, together with the multiple copies of other proteins generally involved in RNA degradation, suggest the existence, in B. cenocepacia and in other Bcc bacteria, of some extra and unexplored functions for the mentioned RBPs, as well as of alternative mechanisms involved in RNA regulation and metabolism in these bacteria.

Transposons played a major role in the diversification between the closely related almond and peach genomes: results from the almond genome sequence

We sequenced the genome of the highly heterozygous almond Prunus dulcis cv. Texas combining short- and long-read sequencing. We obtained a genome assembly totaling 227.6 Mb of the estimated almond genome size of 238 Mb, of which 91% is anchored to eight pseudomolecules corresponding to its haploid chromosome complement, and annotated 27 969 protein-coding genes and 6747 non-coding transcripts. By phylogenomic comparison with the genomes of 16 additional close and distant species we estimated that almond and peach (Prunus persica) diverged around 5.88 million years ago. These two genomes are highly syntenic and show a high degree of sequence conservation (20 nucleotide substitutions per kb). However, they also exhibit a high number of presence/absence variants, many attributable to the movement of transposable elements (TEs). Transposable elements have generated an important number of presence/absence variants between almond and peach, and we show that the recent history of TE movement seems markedly different between them. Transposable elements may also be at the origin of important phenotypic differences between both species, and in particular for the sweet kernel phenotype, a key agronomic and domestication character for almond. Here we show that in sweet almond cultivars, highly methylated TE insertions surround a gene involved in the biosynthesis of amygdalin, whose reduced expression has been correlated with the sweet almond phenotype. Altogether, our results suggest a key role of TEs in the recent history and diversification of almond and its close relative peach.

Comparative Genomics and Transcriptomics To Analyze Fruiting Body Development in Filamentous Ascomycetes

Many filamentous ascomycetes develop three-dimensional fruiting bodies for production and dispersal of sexual spores. Fruiting bodies are among the most complex structures differentiated by ascomycetes; however, the molecular mechanisms underlying this process are insufficiently understood. Previous comparative transcriptomics analyses of fruiting body development in different ascomycetes suggested that there might be a core set of genes that are transcriptionally regulated in a similar manner across species. Conserved patterns of gene expression can be indicative of functional relevance, and therefore such a set of genes might constitute promising candidates for functional analyses. In this study, we have sequenced the genome of the Pezizomycete Ascodesmis nigricans, and performed comparative transcriptomics of developing fruiting bodies of this fungus, the Pezizomycete Pyronema confluens, and the Sordariomycete Sordaria macrospora With only 27 Mb, the A. nigricans genome is the smallest Pezizomycete genome sequenced to date. Comparative transcriptomics indicated that gene expression patterns in developing fruiting bodies of the three species are more similar to each other than to nonsexual hyphae of the same species. An analysis of 83 genes that are upregulated only during fruiting body development in all three species revealed 23 genes encoding proteins with predicted roles in vesicle transport, the endomembrane system, or transport across membranes, and 13 genes encoding proteins with predicted roles in chromatin organization or the regulation of gene expression. Among four genes chosen for functional analysis by deletion in S. macrospora, three were shown to be involved in fruiting body formation, including two predicted chromatin modifier genes.

Evolutionary and functional patterns of shared gene neighbourhood in fungi

Gene clusters comprise genomically co-localized and potentially co-regulated genes that tend to be conserved across species. In eukaryotes, multiple examples of metabolic gene clusters are known, particularly among fungi and plants. However, little is known about how gene clustering patterns vary among taxa or with respect to functional roles. Furthermore, mechanisms of the formation, maintenance and evolution of gene clusters remain unknown. We surveyed 341 fungal genomes to discover gene clusters shared by different species, independently of their functions. We inferred 12,120 cluster families, which comprised roughly one third of the gene space and were enriched in genes associated with diverse cellular functions. Additionally, most clusters did not encode transcription factors, suggesting that they are regulated distally. We used phylogenomics to characterize the evolutionary history of these clusters. We found that most clusters originated once and were transmitted vertically, coupled to differential loss. However, convergent evolution-that is, independent appearance of the same cluster-was more prevalent than anticipated. Finally, horizontal gene transfer of entire clusters was somewhat restricted, with the exception of those associated with secondary metabolism. Altogether, our results provide insights on the evolution of gene clustering as well as a broad catalogue of evolutionarily conserved gene clusters whose function remains to be elucidated.

Adaptation of S. cerevisiae to Fermented Food Environments Reveals Remarkable Genome Plasticity and the Footprints of Domestication

The budding yeast Saccharomyces cerevisiae can be found in the wild and is also frequently associated with human activities. Despite recent insights into the phylogeny of this species, much is still unknown about how evolutionary processes related to anthropogenic niches have shaped the genomes and phenotypes of S. cerevisiae. To address this question, we performed population-level sequencing of 82 S. cerevisiae strains from wine, flor, rum, dairy products, bakeries, and the natural environment (oak trees). These genomic data enabled us to delineate specific genetic groups corresponding to the different ecological niches and revealed high genome content variation across the groups. Most of these strains, compared with the reference genome, possessed additional genetic elements acquired by introgression or horizontal transfer, several of which were population-specific. In addition, several genomic regions in each population showed evidence of nonneutral evolution, as shown by high differentiation, or of selective sweeps including genes with key functions in these environments (e.g., amino acid transport for wine yeast). Linking genetics to lifestyle differences and metabolite traits has enabled us to elucidate the genetic basis of several niche-specific population traits, such as growth on galactose for cheese strains. These data indicate that yeast has been subjected to various divergent selective pressures depending on its niche, requiring the development of customized genomes for better survival in these environments. These striking genome dynamics associated with local adaptation and domestication reveal the remarkable plasticity of the S. cerevisiae genome, revealing this species to be an amazing complex of specialized populations.

Genome Sequence of the Brown Rot Fungal Pathogen Monilinia laxa

Monilinia laxa (phylum Ascomycota) is a plant pathogen responsible for the brown rot blossom blight disease in stone fruit trees of the Rosaceae family, such as apricots. We report here the genome sequence of strain 8L of this species, which was assembled into 618 scaffolds, having a total size of 40.799 Mb and encoding 9,567 unique protein-coding genes.

Probing the Phylogenomics and Putative Pathogenicity Genes of Pythium insidiosum by Oomycete Genome Analyses

Pythium insidiosum is a human-pathogenic oomycete. Many patients infected with it lose organs or die. Toward the goal of developing improved treatment options, we want to understand how Py. insidiosum has evolved to become a successful human pathogen. Our approach here involved the use of comparative genomic and other analyses to identify genes with possible functions in the pathogenicity of Py. insidiosum. We generated an Oomycete Gene Table and used it to explore the genome contents and phylogenomic relationships of Py. insidiosum and 19 other oomycetes. Initial sequence analyses showed that Py. insidiosum is closely related to Pythium species that are not pathogenic to humans. Our analyses also indicated that the organism harbours secreted and adhesin-like proteins, which are absent from related species. Putative virulence proteins were identified by comparison to a set of known virulence genes. Among them is the urease Ure1, which is absent from humans and thus a potential diagnostic and therapeutic target. We used mass spectrometric data to successfully validate the expression of 30% of 14,962 predicted proteins and identify 15 body temperature (37 °C)-dependent proteins of Py. insidiosum. This work begins to unravel the determinants of pathogenicity of Py. insidiosum.

Practical considerations for plant phylogenomics

The past decade has seen a major breakthrough in our ability to easily and inexpensively sequence genome-scale data from diverse lineages. The development of high-throughput sequencing and long-read technologies has ushered in the era of phylogenomics, where hundreds to thousands of nuclear genes and whole organellar genomes are routinely used to reconstruct evolutionary relationships. As a result, understanding which options are best suited for a particular set of questions can be difficult, especially for those just starting in the field. Here, we review the most recent advances in plant phylogenomic methods and make recommendations for project-dependent best practices and considerations. We focus on the costs and benefits of different approaches in regard to the information they provide researchers and the questions they can address. We also highlight unique challenges and opportunities in plant systems, such as polyploidy, reticulate evolution, and the use of herbarium materials, identifying optimal methodologies for each. Finally, we draw attention to lingering challenges in the field of plant phylogenomics, such as reusability of data sets, and look at some up-and-coming technologies that may help propel the field even further.

Phylogenomics of the olive tree (Olea europaea) reveals the relative contribution of ancient allo- and autopolyploidization events

BACKGROUND

Polyploidization is one of the major evolutionary processes that shape eukaryotic genomes, being particularly common in plants. Polyploids can arise through direct genome doubling within a species (autopolyploidization) or through the merging of genomes from distinct species after hybridization (allopolyploidization). The relative contribution of both mechanisms in plant evolution is debated. Here we used phylogenomics to dissect the tempo and mode of duplications in the genome of the olive tree (Olea europaea), one of the first domesticated Mediterranean fruit trees.

RESULTS

Our results depict a complex scenario involving at least three past polyploidization events, of which two-at the bases of the family Oleaceae and the tribe Oleeae, respectively-are likely to be the result of ancient allopolyploidization. A more recent polyploidization involves specifically the olive tree and relatives.

CONCLUSION

Our results show the power of phylogenomics to distinguish between allo- and auto polyploidization events and clarify the contributions of duplications in the evolutionary history of the olive tree.

Improved transcriptome sampling pinpoints 26 ancient and more recent polyploidy events in Caryophyllales, including two allopolyploidy events

Studies of the macroevolutionary legacy of polyploidy are limited by an incomplete sampling of these events across the tree of life. To better locate and understand these events, we need comprehensive taxonomic sampling as well as homology inference methods that accurately reconstruct the frequency and location of gene duplications. We assembled a data set of transcriptomes and genomes from 168 species in Caryophyllales, of which 43 transcriptomes were newly generated for this study, representing one of the most densely sampled genomic-scale data sets available. We carried out phylogenomic analyses using a modified phylome strategy to reconstruct the species tree. We mapped the phylogenetic distribution of polyploidy events by both tree-based and distance-based methods, and explicitly tested scenarios for allopolyploidy. We identified 26 ancient and more recent polyploidy events distributed throughout Caryophyllales. Two of these events were inferred to be allopolyploidy. Through dense phylogenomic sampling, we show the propensity of polyploidy throughout the evolutionary history of Caryophyllales. We also provide a framework for utilizing transcriptome data to detect allopolyploidy, which is important as it may have different macroevolutionary implications compared with autopolyploidy.

Regulatory Mechanisms of a Highly Pectinolytic Mutant of Penicillium occitanis and Functional Analysis of a Candidate Gene in the Plant Pathogen Fusarium oxysporum

Penicillium occitanis is a model system for enzymatic regulation. A mutant strain exhibiting constitutive overproduction of different pectinolytic enzymes both under inducing (pectin) or repressing conditions (glucose) was previously isolated after chemical mutagenesis. In order to identify the molecular basis of this regulatory mechanism, the genomes of the wild type and the derived mutant strain were sequenced and compared, providing the first reference genome for this species. We used a phylogenomic approach to compare P. occitanis with other pectinolytic fungi and to trace expansions of gene families involved in carbohydrate degradation. Genome comparison between wild type and mutant identified seven mutations associated with predicted proteins. The most likely candidate was a mutation in a highly conserved serine residue of a conserved fungal protein containing a GAL4-like Zn2Cys6 binuclear cluster DNA-binding domain and a fungus-specific transcription factor regulatory middle homology region. To functionally characterize the role of this candidate gene, the mutation was recapitulated in the predicted orthologue Fusarium oxysporum, a vascular wilt pathogen which secretes a wide array of plant cell wall degrading enzymes, including polygalacturonases, pectate lyases, xylanases and proteases, all of which contribute to infection. However, neither the null mutant nor a mutant carrying the analogous point mutation exhibited a deregulation of pectinolytic enzymes. The availability, annotation and phylogenomic analysis of the P. occitanis genome sequence represents an important resource for understanding the evolution and biology of this species, and sets the basis for the discovery of new genes of biotechnological interest for the degradation of complex polysaccharides.

Three, two, one yeast fatty acid desaturases: regulation and function

Fatty acid composition of biological membranes functionally adapts to environmental conditions by changing its composition through the activity of lipid biosynthetic enzymes, including the fatty acid desaturases. Three major desaturases are present in yeasts, responsible for the generation of double bonds in position C9-C10, C12-C13 and C15-C16 of the carbon backbone. In this review, we will report data addressed to define the functional role of basidiomycete and ascomycete yeast desaturase enzymes in response to various external signals and the regulation of the expression of their corresponding genes. Many yeast species have the complete set of three desaturases; however, only the Δ9 desaturase seems to be necessary and sufficient to ensure yeast viability. The evolutionary issue of this observation will be discussed.

Widespread Inter- and Intra-Domain Horizontal Gene Transfer of d-Amino Acid Metabolism Enzymes in Eukaryotes

Analysis of the growing number of available fully-sequenced genomes has shown that Horizontal Gene Transfer (HGT) in eukaryotes is more common than previously thought. It has been proposed that genes with certain functions may be more prone to HGT than others, but we still have a very poor understanding of the selective forces driving eukaryotic HGT. Recent work uncovered that d-amino acid racemases have been commonly transferred from bacteria to fungi, but their role in the receiving organisms is currently unknown. Here, we set out to assess whether d-amino acid racemases are commonly transferred to and between eukaryotic groups. For this we performed a global survey that used a novel automated phylogeny-based HGT-detection algorithm (Abaccus). Our results revealed that at least 7.0% of the total eukaryotic racemase repertoire is the result of inter- or intra-domain HGT. These transfers are significantly enriched in plant-associated fungi. For these, we hypothesize a possible role for the acquired racemases allowing to exploit minoritary nitrogen sources in plant biomass, a nitrogen-poor environment. Finally, we performed experiments on a transferred aspartate-glutamate racemase in the fungal human pathogen Candida glabrata, which however revealed no obvious biological role.

Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx

BACKGROUND

Genomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction.

RESULTS

We generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species.

CONCLUSIONS

The genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.

High-Throughput Proteomics Reveals the Unicellular Roots of Animal Phosphosignaling and Cell Differentiation

Cell-specific regulation of protein levels and activity is essential for the distribution of functions among multiple cell types in animals. The finding that many genes involved in these regulatory processes have a premetazoan origin raises the intriguing possibility that the mechanisms required for spatially regulated cell differentiation evolved prior to the appearance of animals. Here, we use high-throughput proteomics in Capsaspora owczarzaki, a close unicellular relative of animals, to characterize the dynamic proteome and phosphoproteome profiles of three temporally distinct cell types in this premetazoan species. We show that life-cycle transitions are linked to extensive proteome and phosphoproteome remodeling and that they affect key genes involved in animal multicellularity, such as transcription factors and tyrosine kinases. The observation of shared features between Capsaspora and metazoans indicates that elaborate and conserved phosphosignaling and proteome regulation supported temporal cell-type differentiation in the unicellular ancestor of animals.

SICLE: a high-throughput tool for extracting evolutionary relationships from phylogenetic trees

We present the phylogeny analysis software SICLE (Sister Clade Extractor), an easy-to-use, high-throughput tool to describe the nearest neighbors to a node of interest in a phylogenetic tree as well as the support value for the relationship. The application is a command line utility that can be embedded into a phylogenetic analysis pipeline or can be used as a subroutine within another C++ program. As a test case, we applied this new tool to the published phylome of Salinibacter ruber, a species of halophilic Bacteriodetes, identifying 13 unique sister relationships to S. ruber across the 4,589 gene phylogenies. S. ruber grouped with bacteria, most often other Bacteriodetes, in the majority of phylogenies, but 91 phylogenies showed a branch-supported sister association between S. ruber and Archaea, an evolutionarily intriguing relationship indicative of horizontal gene transfer. This test case demonstrates how SICLE makes it possible to summarize the phylogenetic information produced by automated phylogenetic pipelines to rapidly identify and quantify the possible evolutionary relationships that merit further investigation. SICLE is available for free for noncommercial use at http://eebweb.arizona.edu/sicle/.

Efficient Gene Tree Correction Guided by Genome Evolution

MOTIVATIONS

Gene trees inferred solely from multiple alignments of homologous sequences often contain weakly supported and uncertain branches. Information for their full resolution may lie in the dependency between gene families and their genomic context. Integrative methods, using species tree information in addition to sequence information, often rely on a computationally intensive tree space search which forecloses an application to large genomic databases.

RESULTS

We propose a new method, called ProfileNJ, that takes a gene tree with statistical supports on its branches, and corrects its weakly supported parts by using a combination of information from a species tree and a distance matrix. Its low running time enabled us to use it on the whole Ensembl Compara database, for which we propose an alternative, arguably more plausible set of gene trees. This allowed us to perform a genome-wide analysis of duplication and loss patterns on the history of 63 eukaryote species, and predict ancestral gene content and order for all ancestors along the phylogeny.

AVAILABILITY

A web interface called RefineTree, including ProfileNJ as well as a other gene tree correction methods, which we also test on the Ensembl gene families, is available at: http://www-ens.iro.umontreal.ca/~adbit/polytomysolver.html. The code of ProfileNJ as well as the set of gene trees corrected by ProfileNJ from Ensembl Compara version 73 families are also made available.

Emergence and Evolution of Hominidae-Specific Coding and Noncoding Genomic Sequences

Family Hominidae, which includes humans and great apes, is recognized for unique complex social behavior and intellectual abilities. Despite the increasing genome data, however, the genomic origin of its phenotypic uniqueness has remained elusive. Clade-specific genes and highly conserved noncoding sequences (HCNSs) are among the high-potential evolutionary candidates involved in driving clade-specific characters and phenotypes. On this premise, we analyzed whole genome sequences along with gene orthology data retrieved from major DNA databases to find Hominidae-specific (HS) genes and HCNSs. We discovered that Down syndrome critical region 4 (DSCR4) is the only experimentally verified gene uniquely present in Hominidae. DSCR4 has no structural homology to any known protein and was inferred to have emerged in several steps through LTR/ERV1, LTR/ERVL retrotransposition, and transversion. Using the genomic distance as neutral evolution threshold, we identified 1,658 HS HCNSs. Polymorphism coverage and derived allele frequency analysis of HS HCNSs showed that these HCNSs are under purifying selection, indicating that they may harbor important functions. They are overrepresented in promoters/untranslated regions, in close proximity of genes involved in sensory perception of sound and developmental process, and also showed a significantly lower nucleosome occupancy probability. Interestingly, many ancestral sequences of the HS HCNSs showed very high evolutionary rates. This suggests that new functions emerged through some kind of positive selection, and then purifying selection started to operate to keep these functions.

Lokiarchaeota Marks the Transition between the Archaeal and Eukaryotic Selenocysteine Encoding Systems

Selenocysteine (Sec) is the 21st amino acid in the genetic code, inserted in response to UGA codons with the help of RNA structures, the SEC Insertion Sequence (SECIS) elements. The three domains of life feature distinct strategies for Sec insertion in proteins and its utilization. While bacteria and archaea possess similar sets of selenoproteins, Sec biosynthesis is more similar among archaea and eukaryotes. However, SECIS elements are completely different in the three domains of life. Here, we analyze the archaeon Lokiarchaeota that resolves the relationships among Sec insertion systems. This organism has selenoproteins representing five protein families, three of which have multiple Sec residues. Remarkably, these archaeal selenoprotein genes possess conserved RNA structures that strongly resemble the eukaryotic SECIS element, including key eukaryotic protein-binding sites. These structures also share similarity with the SECIS element in archaeal selenoprotein VhuD, suggesting a relation of direct descent. These results identify Lokiarchaeota as an intermediate form between the archaeal and eukaryotic Sec-encoding systems and clarify the evolution of the Sec insertion system.

Computational Identification of the Paralogs and Orthologs of Human Cytochrome P450 Superfamily and the Implication in Drug Discovery

The human cytochrome P450 (CYP) superfamily consisting of 57 functional genes is the most important group of Phase I drug metabolizing enzymes that oxidize a large number of xenobiotics and endogenous compounds, including therapeutic drugs and environmental toxicants. The CYP superfamily has been shown to expand itself through gene duplication, and some of them become pseudogenes due to gene mutations. Orthologs and paralogs are homologous genes resulting from speciation or duplication, respectively. To explore the evolutionary and functional relationships of human CYPs, we conducted this bioinformatic study to identify their corresponding paralogs, homologs, and orthologs. The functional implications and implications in drug discovery and evolutionary biology were then discussed. GeneCards and Ensembl were used to identify the paralogs of human CYPs. We have used a panel of online databases to identify the orthologs of human CYP genes: NCBI, Ensembl Compara, GeneCards, OMA ("Orthologous MAtrix") Browser, PATHER, TreeFam, EggNOG, and Roundup. The results show that each human CYP has various numbers of paralogs and orthologs using GeneCards and Ensembl. For example, the paralogs of CYP2A6 include CYP2A7, 2A13, 2B6, 2C8, 2C9, 2C18, 2C19, 2D6, 2E1, 2F1, 2J2, 2R1, 2S1, 2U1, and 2W1; CYP11A1 has 6 paralogs including CYP11B1, 11B2, 24A1, 27A1, 27B1, and 27C1; CYP51A1 has only three paralogs: CYP26A1, 26B1, and 26C1; while CYP20A1 has no paralog. The majority of human CYPs are well conserved from plants, amphibians, fishes, or mammals to humans due to their important functions in physiology and xenobiotic disposition. The data from different approaches are also cross-validated and validated when experimental data are available. These findings facilitate our understanding of the evolutionary relationships and functional implications of the human CYP superfamily in drug discovery.

The link between orthology relations and gene trees: a correction perspective

Background

While tree-oriented methods for inferring orthology and paralogy relations between genes are based on reconciling a gene tree with a species tree, many tree-free methods are also available (usually based on sequence similarity). Recently, the link between orthology relations and gene trees has been formally considered from the perspective of reconstructing phylogenies from orthology relations. In this paper, we consider this link from a correction point of view. Indeed, a gene tree induces a set of relations, but the converse is not always true: a set of relations is not necessarily in agreement with any gene tree. A natural question is thus how to minimally correct an infeasible set of relations. Another natural question, given a gene tree and a set of relations, is how to minimally correct a gene tree so that the resulting gene tree fits the set of relations.

Results

We consider four variants of relation and gene tree correction problems, and provide hardness results for all of them. More specifically, we show that it is NP-Hard to edit a minimum of set of relations to make them consistent with a given species tree. We also show that the problem of finding a maximum subset of genes that share consistent relations is hard to approximate. We then demonstrate that editing a gene tree to satisfy a given set of relations in a minimum way is NP-Hard, where “minimum” refers either to the number of modified relations depicted by the gene tree or the number of clades that are lost. We also discuss some of the algorithmic perspectives given these hardness results.

Whole genome sequencing of turbot (Scophthalmus maximus; Pleuronectiformes): a fish adapted to demersal life

The turbot is a flatfish (Pleuronectiformes) with increasing commercial value, which has prompted active genomic research aimed at more efficient selection. Here we present the sequence and annotation of the turbot genome, which represents a milestone for both boosting breeding programmes and ascertaining the origin and diversification of flatfish. We compare the turbot genome with model fish genomes to investigate teleost chromosome evolution. We observe a conserved macrosyntenic pattern within Percomorpha and identify large syntenic blocks within the turbot genome related to the teleost genome duplication. We identify gene family expansions and positive selection of genes associated with vision and metabolism of membrane lipids, which suggests adaptation to demersal lifestyle and to cold temperatures, respectively. Our data indicate a quick evolution and diversification of flatfish to adapt to benthic life and provide clues for understanding their controversial origin. Moreover, we investigate the genomic architecture of growth, sex determination and disease resistance, key traits for understanding local adaptation and boosting turbot production, by mapping candidate genes and previously reported quantitative trait loci. The genomic architecture of these productive traits has allowed the identification of candidate genes and enriched pathways that may represent useful information for future marker-assisted selection in turbot.

Late acquisition of mitochondria by a host with chimaeric prokaryotic ancestry

The origin of eukaryotes stands as a major conundrum in biology. Current evidence indicates that the last eukaryotic common ancestor already possessed many eukaryotic hallmarks, including a complex subcellular organization. In addition, the lack of evolutionary intermediates challenges the elucidation of the relative order of emergence of eukaryotic traits. Mitochondria are ubiquitous organelles derived from an alphaproteobacterial endosymbiont. Different hypotheses disagree on whether mitochondria were acquired early or late during eukaryogenesis. Similarly, the nature and complexity of the receiving host are debated, with models ranging from a simple prokaryotic host to an already complex proto-eukaryote. Most competing scenarios can be roughly grouped into either mito-early, which consider the driving force of eukaryogenesis to be mitochondrial endosymbiosis into a simple host, or mito-late, which postulate that a significant complexity predated mitochondrial endosymbiosis. Here we provide evidence for late mitochondrial endosymbiosis. We use phylogenomics to directly test whether proto-mitochondrial proteins were acquired earlier or later than other proteins of the last eukaryotic common ancestor. We find that last eukaryotic common ancestor protein families of alphaproteobacterial ancestry and of mitochondrial localization show the shortest phylogenetic distances to their closest prokaryotic relatives, compared with proteins of different prokaryotic origin or cellular localization. Altogether, our results shed new light on a long-standing question and provide compelling support for the late acquisition of mitochondria into a host that already had a proteome of chimaeric phylogenetic origin. We argue that mitochondrial endosymbiosis was one of the ultimate steps in eukaryogenesis and that it provided the definitive selective advantage to mitochondria-bearing eukaryotes over less complex forms.

Genome and transcriptome analysis of the Mesoamerican common bean and the role of gene duplications in establishing tissue and temporal specialization of genes

BACKGROUND

Legumes are the third largest family of angiosperms and the second most important crop class. Legume genomes have been shaped by extensive large-scale gene duplications, including an approximately 58 million year old whole genome duplication shared by most crop legumes.

RESULTS

We report the genome and the transcription atlas of coding and non-coding genes of a Mesoamerican genotype of common bean (Phaseolus vulgaris L., BAT93). Using a comprehensive phylogenomics analysis, we assessed the past and recent evolution of common bean, and traced the diversification of patterns of gene expression following duplication. We find that successive rounds of gene duplications in legumes have shaped tissue and developmental expression, leading to increased levels of specialization in larger gene families. We also find that many long non-coding RNAs are preferentially expressed in germ-line-related tissues (pods and seeds), suggesting that they play a significant role in fruit development. Our results also suggest that most bean-specific gene family expansions, including resistance gene clusters, predate the split of the Mesoamerican and Andean gene pools.

CONCLUSIONS

The genome and transcriptome data herein generated for a Mesoamerican genotype represent a counterpart to the genomic resources already available for the Andean gene pool. Altogether, this information will allow the genetic dissection of the characters involved in the domestication and adaptation of the crop, and their further implementation in breeding strategies for this important crop.

Genome sequencing and secondary metabolism of the postharvest pathogen Penicillium griseofulvum

BACKGROUND

Penicillium griseofulvum is associated in stored apples with blue mould, the most important postharvest disease of pome fruit. This pathogen can simultaneously produce both detrimental and beneficial secondary metabolites (SM). In order to gain insight into SM synthesis in P. griseofulvum in vitro and during disease development on apple, we sequenced the genome of P. griseofulvum strain PG3 and analysed important SM clusters.

RESULTS

PG3 genome sequence (29.3 Mb) shows that P. griseofulvum branched off after the divergence of P. oxalicum but before the divergence of P. chrysogenum. Genome-wide analysis of P. griseofulvum revealed putative gene clusters for patulin, griseofulvin and roquefortine C biosynthesis. Furthermore, we quantified the SM production in vitro and on apples during the course of infection. The expression kinetics of key genes of SM produced in infected apple were examined. We found additional SM clusters, including those potentially responsible for the synthesis of penicillin, yanuthone D, cyclopiazonic acid and we predicted a cluster putatively responsible for the synthesis of chanoclavine I.

CONCLUSIONS

These findings provide relevant information to understand the molecular basis of SM biosynthesis in P. griseofulvum, to allow further research directed to the overexpression or blocking the synthesis of specific SM.

Contrasting Genomic Diversity in Two Closely Related Postharvest Pathogens: Penicillium digitatum and Penicillium expansum

Penicillium digitatum and Penicillium expansum are two closely related fungal plant pathogens causing green and blue mold in harvested fruit, respectively. The two species differ in their host specificity, being P. digitatum restricted to citrus fruits and P. expansum able to infect a wide range of fruits after harvest. Although host-specific Penicillium species have been found to have a smaller gene content, it is so far unclear whether these different host specificities impact genome variation at the intraspecific level. Here we assessed genome variation across four P. digitatum and seven P. expansum isolates from geographically distant regions. Our results show very high similarity (average 0.06 SNPs [single nucleotide polymorphism] per kb) between globally distributed isolates of P. digitatum pointing to a recent expansion of a single lineage. This low level of genetic variation found in our samples contrasts with the higher genetic variability observed in the similarly distributed P. expansum isolates (2.44 SNPs per kb). Patterns of polymorphism in P. expansum indicate that recombination exists between genetically diverged strains. Consistent with the existence of sexual recombination and heterothallism, which was unknown for this species, we identified the two alternative mating types in different P. expansum isolates. Patterns of polymorphism in P. digitatum indicate a recent clonal population expansion of a single lineage that has reached worldwide distribution. We suggest that the contrasting patterns of genomic variation between the two species reflect underlying differences in population dynamics related with host specificities and related agricultural practices. It should be noted, however, that this results should be confirmed with a larger sampling of strains, as new strains may broaden the diversity so far found in P. digitatum.

Horizontal acquisition of toxic alkaloid synthesis in a clade of plant associated fungi

•

Loline and ergot alkaloids are fundamental for Clavicipitaceae.

•

These secondary metabolism clusters are found in other, far related, species.

•

Phylogenetics establish a close relationship between clusters of different species.

•

Detection of large amounts of gene loss supports the HGT events.

Clavicipitaceae is a fungal group that comprises species that closely interact with plants as pathogens, parasites or symbionts. A key factor in these interactions is the ability of these fungi to synthesize toxic alkaloid compounds that contribute to the protection of the plant host against herbivores. Some of these compounds such as ergot alkaloids are toxic to humans and have caused important epidemics throughout history. The gene clusters encoding the proteins responsible for the synthesis of ergot alkaloids and lolines in Clavicipitaceae have been elucidated. Notably, homologs to these gene clusters can be found in distantly related species such as Aspergillus fumigatus and Penicillium expansum, which diverged from Clavicipitaceae more than 400 million years ago. We here use a phylogenetic approach to analyze the evolution of these gene clusters. We found that the gene clusters conferring the ability to synthesize ergot alkaloids and loline emerged first in Eurotiomycetes and were then likely transferred horizontally to Clavicipitaceae. Horizontal gene transfer is known to play a role in shaping the distribution of secondary metabolism clusters across distantly related fungal species. We propose that HGT events have played an important role in the capability of Clavicipitaceae to produce two key secondary metabolites that have enhanced the ability of these species to protect their plant hosts, therefore favoring their interactions.

Genome of Rhodnius prolixus, an insect vector of Chagas disease, reveals unique adaptations to hematophagy and parasite infection

Rhodnius prolixus not only has served as a model organism for the study of insect physiology, but also is a major vector of Chagas disease, an illness that affects approximately seven million people worldwide. We sequenced the genome of R. prolixus, generated assembled sequences covering 95% of the genome (∼ 702 Mb), including 15,456 putative protein-coding genes, and completed comprehensive genomic analyses of this obligate blood-feeding insect. Although immune-deficiency (IMD)-mediated immune responses were observed, R. prolixus putatively lacks key components of the IMD pathway, suggesting a reorganization of the canonical immune signaling network. Although both Toll and IMD effectors controlled intestinal microbiota, neither affected Trypanosoma cruzi, the causal agent of Chagas disease, implying the existence of evasion or tolerance mechanisms. R. prolixus has experienced an extensive loss of selenoprotein genes, with its repertoire reduced to only two proteins, one of which is a selenocysteine-based glutathione peroxidase, the first found in insects. The genome contained actively transcribed, horizontally transferred genes from Wolbachia sp., which showed evidence of codon use evolution toward the insect use pattern. Comparative protein analyses revealed many lineage-specific expansions and putative gene absences in R. prolixus, including tandem expansions of genes related to chemoreception, feeding, and digestion that possibly contributed to the evolution of a blood-feeding lifestyle. The genome assembly and these associated analyses provide critical information on the physiology and evolution of this important vector species and should be instrumental for the development of innovative disease control methods.

The Leishmania metaphylome: a comprehensive survey of Leishmania protein phylogenetic relationships

Background

Leishmaniasis is a neglected parasitic disease with diverse clinical manifestations and a complex epidemiology. It has been shown that its parasite-related traits vary between species and that they modulate infectivity, pathogenicity, and virulence. However, understanding of the species-specific adaptations responsible for these features and their evolutionary background is limited. To improve our knowledge regarding the parasite biology and adaptation mechanisms of different Leishmania species, we conducted a proteome-wide phylogenomic analysis to gain insights into Leishmania evolution.

Results

The analysis of the reconstructed phylomes (totaling 45,918 phylogenies) allowed us to detect genes that are shared in pathogenic Leishmania species, such as calpain-like cysteine peptidases and 3'a2rel-related proteins, or genes that could be associated with visceral or cutaneous development. This analysis also established the phylogenetic relationship of several hypothetical proteins whose roles remain to be characterized. Our findings demonstrated that gene duplication constitutes an important evolutionary force in Leishmania, acting on protein families that mediate host-parasite interactions, such as amastins, GP63 metallopeptidases, cathepsin L-like proteases, and our methods permitted a deeper analysis of their phylogenetic relationships.

Conclusions

Our results highlight the importance of proteome wide phylogenetic analyses to detect adaptation and evolutionary processes in different organisms and underscore the need to characterize the role of expanded and species-specific proteins in the context of Leishmania evolution by providing a framework for the phylogenetic relationships of Leishmania proteins.

Phylogenomic data are publicly available for use through PhylomeDB (http://www.phylomedb.org).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2091-2) contains supplementary material, which is available to authorized users.

An Effective Big Data Supervised Imbalanced Classification Approach for Ortholog Detection in Related Yeast Species

Orthology detection requires more effective scaling algorithms. In this paper, a set of gene pair features based on similarity measures (alignment scores, sequence length, gene membership to conserved regions, and physicochemical profiles) are combined in a supervised pairwise ortholog detection approach to improve effectiveness considering low ortholog ratios in relation to the possible pairwise comparison between two genomes. In this scenario, big data supervised classifiers managing imbalance between ortholog and nonortholog pair classes allow for an effective scaling solution built from two genomes and extended to other genome pairs. The supervised approach was compared with RBH, RSD, and OMA algorithms by using the following yeast genome pairs: Saccharomyces cerevisiae-Kluyveromyces lactis, Saccharomyces cerevisiae-Candida glabrata, and Saccharomyces cerevisiae-Schizosaccharomyces pombe as benchmark datasets. Because of the large amount of imbalanced data, the building and testing of the supervised model were only possible by using big data supervised classifiers managing imbalance. Evaluation metrics taking low ortholog ratios into account were applied. From the effectiveness perspective, MapReduce Random Oversampling combined with Spark SVM outperformed RBH, RSD, and OMA, probably because of the consideration of gene pair features beyond alignment similarities combined with the advances in big data supervised classification.

Reconstructing a SuperGeneTree minimizing reconciliation

Combining a set of trees on partial datasets into a single tree is a classical method for inferring large phylogenetic trees. Ideally, the combined tree should display each input partial tree, which is only possible if input trees do not contain contradictory phylogenetic information. The simplest version of the supertree problem is thus to state whether a set of trees is compatible, and if so, construct a tree displaying them all. Classically, supertree methods have been applied to the reconstruction of species trees. Here we rather consider reconstructing a super gene tree in light of a known species tree S. We define the supergenetree problem as finding, among all supertrees displaying a set of input gene trees, one supertree minimizing a reconciliation distance with S. We first show how classical exact methods to the supertree problem can be extended to the supergenetree problem. As all these methods are highly exponential, we also exhibit a natural greedy heuristic for the duplication cost, based on minimizing the set of duplications preceding the first speciation event. We then show that both the supergenetree problem and its restriction to minimizing duplications preceding the first speciation are NP-hard to approximate within a n^1-ϵ factor, for any 0 < ϵ < 1. Finally, we show that a restriction of this problem to uniquely labeled speciation gene trees, which is relevant to many biological applications, is also NP-hard. Therefore, we introduce new avenues in the field of supertrees, and set the theoretical basis for the exploration of various algorithmic aspects of the problems.

The Genomic Aftermath of Hybridization in the Opportunistic Pathogen Candida metapsilosis

Candida metapsilosis is a rarely-isolated, opportunistic pathogen that belongs to a clade of pathogenic yeasts known as the C. parapsilosis sensu lato species complex. To gain insight into the recent evolution of C. metapsilosis and the genetic basis of its virulence, we sequenced the genome of 11 clinical isolates from various locations, which we compared to each other and to the available genomes of the two remaining members of the complex: C. orthopsilosis and C. parapsilosis. Unexpectedly, we found compelling genomic evidence that C. metapsilosis is a highly heterozygous hybrid species, with all sequenced clinical strains resulting from the same past hybridization event involving two parental lineages that were approximately 4.5% divergent in sequence. This result indicates that the parental species are non-pathogenic, but that hybridization between them formed a new opportunistic pathogen, C. metapsilosis, that has achieved a worldwide distribution. We show that these hybrids are diploid and we identified strains carrying loci for both alternative mating types, which supports mating as the initial mechanism for hybrid formation. We trace the aftermath of this hybridization at the genomic level, and reconstruct the evolutionary relationships among the different strains. Recombination and introgression -resulting in loss of heterozygosis- between the two subgenomes have been rampant, and includes the partial overwriting of the MTLa mating locus in all strains. Collectively, our results shed light on the recent genomic evolution within the C. parapsilosis sensu lato complex, and argue for a re-definition of species within this clade, with at least five distinct homozygous lineages, some of which having the ability to form hybrids.

C5orf30 is a negative regulator of tissue damage in rheumatoid arthritis

The variant rs26232, in the first intron of the chromosome 5 open reading frame 30 (C5orf30) locus, has recently been associated with both risk of developing rheumatoid arthritis (RA) and severity of tissue damage. The biological activities of human C5orf30 are unknown, and neither the gene nor protein show significant homology to any other characterized human sequences. The C5orf30 gene is present only in vertebrate genomes with a high degree of conservation, implying a central function in these organisms. Here, we report that C5orf30 is highly expressed in the synovium of RA patients compared with control synovial tissue, and that it is predominately expressed by synovial fibroblast (RASF) and macrophages in the lining and sublining layer of the tissue. These cells play a central role in the initiation and perpetuation of RA and are implicated in cartilage destruction. RASFs lacking C5orf30 exhibit increased cell migration and invasion in vitro, and gene profiling following C5orf30 inhibition confirmed up-regulation of genes involved in cell migration, adhesion, angiogenesis, and immune and inflammatory pathways. Importantly, loss of C5orf30 contributes to the pathology of inflammatory arthritis in vivo, because inhibition of C5orf30 in the collagen-induced arthritis model markedly accentuated joint inflammation and tissue damage. Our study reveal C5orf30 to be a previously unidentified negative regulator of tissue damage in RA, and this protein may act by modulating the autoaggressive phenotype that is characteristic of RASFs.

Role of Cln1 during melanization of Cryptococcus neoformans

Cryptococcus neoformans is an opportunistic fungal pathogen that has several well-described virulence determinants. A polysaccharide capsule and the ability to produce melanin are among the most important. Melanization occurs both in vitro, in the presence of catecholamine and indole compounds, and in vivo during the infection. Despite the importance of melanin production for cryptococcal virulence, the component and mechanisms involved in its synthesis have not been fully elucidated. In this work, we describe the role of a G1/S cyclin (Cln1) in the melanization process. Cln1 has evolved specifically with proteins present only in other basidiomycetes. We found that Cln1 is required for the cell wall stability and production of melanin in C. neoformans. Absence of melanization correlated with a defect in the expression of the LAC1 gene. The relation between cell cycle elements and melanization was confirmed by the effect of drugs that cause cell cycle arrest at a specific phase, such as rapamycin. The cln1 mutant was consistently more susceptible to oxidative damage in a medium that induces melanization. Our results strongly suggest a novel and hitherto unrecognized role for C. neoformans Cln1 in the expression of virulence traits.