E value cutoff and eukaryotic genome content phylogenetics

Comparative Transcriptomics to Identify RNA Writers and Erasers in Microalgae

Epitranscriptomics is considered as a new regulatory step in eukaryotes for developmental processes and stress responses. The aim of this study was, for the first time, to identify RNA methyltransferase (writers) and demethylase (erasers) in four investigated species, i.e., the dinoflagellates Alexandrium tamutum and Amphidinium carterae, the diatom Cylindrotheca closterium, and the green alga Tetraselmis suecica. As query sequences for the enzymatic classes of interest, we selected those ones that were previously detected in marine plants, evaluating their expression upon nutrient starvation stress exposure. The hypothesis was that upon stress exposure, the activation/deactivation of specific writers and erasers may occur. In microalgae, we found almost all plant writers and erasers (ALKBH9B, ALKBH10B, MTB, and FIP37), except for three writers (MTA, VIRILIZER, and HAKAI). A sequence similarity search by scanning the corresponding genomes confirmed their presence. Thus, we concluded that the three writer sequences were lacking from the studied transcriptomes probably because they were not expressed in those experimental conditions, rather than a real lack of these genes from their genomes. This study showed that some of them were expressed only in specific culturing conditions. We also investigated their expression in other culturing conditions (i.e., nitrogen depletion, phosphate depletion, and Zinc addition at two different concentrations) in A. carterae, giving new insights into their possible roles in regulating gene expression upon stress.

Placing human gene families into their evolutionary context

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

m6A RNA Methylation in Marine Plants: First Insights and Relevance for Biological Rhythms

Circadian regulations are essential for enabling organisms to synchronize physiology with environmental light-dark cycles. Post-transcriptional RNA modifications still represent an understudied level of gene expression regulation in plants, although they could play crucial roles in environmental adaptation. N⁶-methyl-adenosine (m⁶A) is the most prevalent mRNA modification, established by "writer" and "eraser" proteins. It influences the clockwork in several taxa, but only few studies have been conducted in plants and none in marine plants. Here, we provided a first inventory of m⁶A-related genes in seagrasses and investigated daily changes in the global RNA methylation and transcript levels of writers and erasers in Cymodocea nodosa and Zostera marina. Both species showed methylation peaks during the dark period under the same photoperiod, despite exhibiting asynchronous changes in the m⁶A profile and related gene expression during a 24-h cycle. At contrasting latitudes, Z. marina populations displayed overlapping daily patterns of the m⁶A level and related gene expression. The observed rhythms are characteristic for each species and similar in populations of the same species with different photoperiods, suggesting the existence of an endogenous circadian control. Globally, our results indicate that m⁶A RNA methylation could widely contribute to circadian regulation in seagrasses, potentially affecting the photo-biological behaviour of these plants.

Insights into the genome structure of four acetogenic bacteria with specific reference to the Wood-Ljungdahl pathway

Acetogenic bacteria are obligate anaerobes with the ability of converting carbon dioxide and other one‐carbon substrates into acetate through the Wood–Ljungdahl (WL) pathway. These substrates are becoming increasingly important feedstock in industrial microbiology. The main potential industrial application of acetogenic bacteria is the production of metabolites that constitute renewable energy sources (biofuel); such bacteria are of particular interest for this purpose thanks to their low energy requirements for large‐scale cultivation. Here, we report new genome sequences for four species, three of them are reported for the first time, namely Acetobacterium paludosum DSM 8237, Acetobacterium tundrae DSM 917, Acetobacterium bakii DSM 8239, and Alkalibaculum bacchi DSM 221123. We performed a comparative genomic analysis focused on the WL pathway's genes and their encoded proteins, using Acetobacterium woodii as a reference genome. The Average Nucleotide Identity (ANI) values ranged from 70% to 95% over an alignment length of 5.4–6.5 Mbp. The core genome consisted of 363 genes, whereas the number of unique genes in a single genome ranged from 486 in A. tundrae to 2360 in A.bacchi. No significant rearrangements were detected in the gene order for the Wood–Ljungdahl pathway however, two species showed variations in genes involved in formate metabolism: A. paludosum harbor two copies of fhs1, and A. bakii a truncated fdhF1. The analysis of protein networks highlighted the expansion of protein orthologues in A. woodii compared to A. bacchi, whereas protein networks involved in the WL pathway were more conserved. This study has increased our understanding on the evolution of the WL pathway in acetogenic bacteria.

Multilevel comparative bioinformatics to investigate evolutionary relationships and specificities in gene annotations: an example for tomato and grapevine

Background

“Omics” approaches may provide useful information for a deeper understanding of speciation events, diversification and function innovation. This can be achieved by investigating the molecular similarities at sequence level between species, allowing the definition of ortholog and paralog genes. However, the spreading of sequenced genome, often endowed with still preliminary annotations, requires suitable bioinformatics to be appropriately exploited in this framework.

Results

We presented here a multilevel comparative approach to investigate on genome evolutionary relationships and peculiarities of two fleshy fruit species of relevant agronomic interest, Solanum lycopersicum (tomato) and Vitis vinifera (grapevine). We defined 17,823 orthology relationships between tomato and grapevine reference gene annotations. The resulting orthologs are associated with the detected paralogs in each species, permitting the definition of gene networks, useful to investigate the different relationships. The reconciliation of the compared collections in terms of an updating of the functional descriptions was also exploited. All the results were made accessible in ComParaLogs, a dedicated bioinformatics platform available at http://biosrv.cab.unina.it/comparalogs/gene/search.

Conclusions

The aim of the work was to suggest a reliable approach to detect all similarities of gene loci between two species based on the integration of results from different levels of information, such as the gene, the transcript and the protein sequences, overcoming possible limits due to exclusive protein versus protein comparisons. This to define reliable ortholog and paralog genes, as well as species specific gene loci in the two species, overcoming limits due to the possible draft nature of preliminary gene annotations. Moreover, reconciled functional descriptions, as well as common or peculiar enzymatic classes and protein domains from tomato and grapevine, together with the definition of species-specific gene sets after the pairwise comparisons, contributed a comprehensive set of information useful to comparatively exploit the two species gene annotations and investigate on differences between species with climacteric and non-climacteric fruits. In addition, the definition of networks of ortholog genes and of associated paralogs, and the organization of web-based interfaces for the exploration of the results, defined a friendly computational bench-work in support of comparative analyses between two species.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2420-y) contains supplementary material, which is available to authorized users.

New Insights Into Beclin-1: Evolution and Pan-Malignancy Inhibitor Activity

Autophagy is a functionally conserved self-degradation process that facilitates the survival of eukaryotic life via the management of cellular bioenergetics and maintenance of the fidelity of genomic DNA. The first known autophagy inducer was Beclin-1. Beclin-1 is expressed in multicellular eukaryotes ranging throughout plants to animals, comprising a nonmonophyllic group, as shown in this report via aggressive BLAST searches. In humans, Beclin-1 is a haploinsuffient tumor suppressor as biallelic deletions have not been observed in patient tumors clinically. Therefore, Beclin-1 fails the Knudson hypothesis, implicating expression of at least one Beclin-1 allele is essential for cancer cell survival. However, Beclin-1 is frequently monoallelically deleted in advanced human cancers and the expression of two Beclin-1 allelles is associated with greater anticancer effects. Overall, experimental evidence suggests that Beclin-1 inhibits tumor formation, angiogenesis, and metastasis alone and in cooperation with the tumor suppressive molecules UVRAG, Bif-1, Ambra1, and MDA-7/IL-24 via diverse mechanisms of action. Conversely, Beclin-1 is upregulated in cancer stem cells (CSCs), portending a role in cancer recurrence, and highlighting this molecule as an intriguing molecular target for the treatment of CSCs. Many aspects of Beclin-1's biological effects remain to be studied. The consequences of these BLAST searches on the molecular evolution of Beclin-1, and the eukaryotic branches of the tree of life, are discussed here in greater detail with future inquiry focused upon protist taxa. Also in this review, the effects of Beclin-1 on tumor suppression and cancer malignancy are discussed. Beclin-1 holds significant promise for the development of novel targeted cancer therapeutics and is anticipated to lead to a many advances in our understanding of eukaryotic evolution, multicellularity, and even the treatment of CSCs in the coming decades.

Akaryotes and Eukaryotes are independent descendants of a universal common ancestor

We reconstructed a global tree of life (ToL) with non-reversible and non-stationary models of genome evolution that root trees intrinsically. We implemented Bayesian model selection tests and compared the statistical support for four conflicting ToL hypotheses. We show that reconstructions obtained with a Bayesian implementation (Klopfstein et al., 2015) are consistent with reconstructions obtained with an empirical Sankoff parsimony (ESP) implementation (Harish et al., 2013). Both are based on the genome contents of coding sequences for protein domains (superfamilies) from hundreds of genomes. Thus, we conclude that the independent descent of Eukaryotes and Akaryotes (archaea and bacteria) from the universal common ancestor (UCA) is the most probable as well as the most parsimonious hypothesis for the evolutionary origins of extant genomes. Reconstructions of ancestral proteomes by both Bayesian and ESP methods suggest that at least 70% of unique domain-superfamilies known in extant species were present in the UCA. In addition, identification of a vast majority (96%) of the mitochondrial superfamilies in the UCA proteome precludes a symbiotic hypothesis for the origin of eukaryotes. Accordingly, neither the archaeal origin of eukaryotes nor the bacterial origin of mitochondria is supported by the data. The proteomic complexity of the UCA suggests that the evolution of cellular phenotypes in the two primordial lineages, Akaryotes and Eukaryotes, was driven largely by duplication of common superfamilies as well as by loss of unique superfamilies. Finally, innovation of novel superfamilies has played a surprisingly small role in the evolution of Akaryotes and only a marginal role in the evolution of Eukaryotes.

Genome content analysis yields new insights into the relationship between the human malaria parasite Plasmodium falciparum and its anopheline vectors

Background

The persistent and growing gap between the availability of sequenced genomes and the ability to assign functions to sequenced genes led us to explore ways to maximize the information content of automated annotation for studies of anopheline mosquitos. Specifically, we use genome content analysis of a large number of previously sequenced anopheline mosquitos to follow the loss and gain of protein families over the evolutionary history of this group.

The importance of this endeavor lies in the potential for comparative genomic studies between Anopheles and closely related non-vector species to reveal ancestral genome content dynamics involved in vector competence. In addition, comparisons within Anopheles could identify genome content changes responsible for variation in the vectorial capacity of this family of important parasite vectors.

Results

The competence and capacity of P. falciparum vectors do not appear to be phylogenetically constrained within the Anophelinae. Instead, using ancestral reconstruction methods, we suggest that a previously unexamined component of vector biology, anopheline nucleotide metabolism, may contribute to the unique status of anophelines as P. falciparum vectors. While the fitness effects of nucleotide co-option by P. falciparum parasites on their anopheline hosts are not yet known, our results suggest that anopheline genome content may be responding to selection pressure from P. falciparum. Whether this response is defensive, in an attempt to redress improper nucleotide balance resulting from P. falciparum infection, or perhaps symbiotic, resulting from an as-yet-unknown mutualism between anophelines and P. falciparum, is an open question that deserves further study.

Conclusions

Clearly, there is a wealth of functional information to be gained from detailed manual genome annotation, yet the rapid increase in the number of available sequences means that most researchers will not have the time or resources to manually annotate all the sequence data they generate. We believe that efforts to maximize the amount of information obtained from automated annotation can help address the functional annotation deficit that most evolutionary biologists now face, and here demonstrate the value of such an approach.

Electronic supplementary material

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

A whole genome gene content phylogenetic analysis of anopheline mosquitoes

Construction of stringent gene content matrices was accomplished for 21 Anopheline mosquito species and strains and four outgroups species. The presence absence matrix using e^-75 as a cutoff in single linkage clustering had over 17,000 ortholog groups. We used the gene content matrix to generate a phylogenetic hypothesis that is in general agreement with gene sequence based phylogenies. In addition to establishing a congruent gene content phylogeny we examined the consistency of three methods for analyzing presence absence data - unweighted parsimony, dollo parsimonly and maximum likelihood using a BINGAMMA model. An examination of the chromosomal location of the gains and losses in the presence absence matrix revealed a low frequency of gains and losses at centromeres and tips of chromosomes.

The phylogenetic position of ctenophores and the origin(s) of nervous systems

Ctenophores have traditionally been treated as eumetazoans, but some recent whole genome studies have revived the idea that they are, rather, the sister group to all other metazoans. This deep branching position implies either that nervous systems have evolved twice, in Ctenophora and in Eumetazoa, or that an ancestral metazoan nervous system has been lost in sponges and placozoans. We caution, however, that phylogenetic-tree construction artifacts may have placed ctenophores too deep in the metazoan tree. We discuss nervous system origins under these alternative phylogenies and in light of comparative data of ctenophore and eumetazoan nervous systems. We argue that characters like neuropeptide signaling, ciliary photoreceptors, gap junctions and presynaptic molecules are consistent with a shared ancestry of nervous systems. However, if ctenophores are the sister group to all other metazoans, this ancestral nervous system was likely very simple. Further studies are needed to resolve the deep phylogeny of metazoans and to have a better understanding of the early steps of nervous system evolution.

Exploiting a Reference Genome in Terms of Duplications: The Network of Paralogs and Single Copy Genes in Arabidopsis thaliana

Arabidopsis thaliana became the model organism for plant studies because of its small diploid genome, rapid lifecycle and short adult size. Its genome was the first among plants to be sequenced, becoming the reference in plant genomics. However, the Arabidopsis genome is characterized by an inherently complex organization, since it has undergone ancient whole genome duplications, followed by gene reduction, diploidization events and extended rearrangements, which relocated and split up the retained portions. These events, together with probable chromosome reductions, dramatically increased the genome complexity, limiting its role as a reference. The identification of paralogs and single copy genes within a highly duplicated genome is a prerequisite to understand its organization and evolution and to improve its exploitation in comparative genomics. This is still controversial, even in the widely studied Arabidopsis genome. This is also due to the lack of a reference bioinformatics pipeline that could exhaustively identify paralogs and singleton genes. We describe here a complete computational strategy to detect both duplicated and single copy genes in a genome, discussing all the methodological issues that may strongly affect the results, their quality and their reliability. This approach was used to analyze the organization of Arabidopsis nuclear protein coding genes, and besides classifying computationally defined paralogs into networks and single copy genes into different classes, it unraveled further intriguing aspects concerning the genome annotation and the gene relationships in this reference plant species. Since our results may be useful for comparative genomics and genome functional analyses, we organized a dedicated web interface to make them accessible to the scientific community