Assessing performance of orthology detection strategies applied to eukaryotic genomes

cascAGS: Comparative Analysis of SNP Calling Methods for Human Genome Data in the Absence of Gold Standard

The development of third-generation sequencing has accelerated the boom of single nucleotide polymorphism (SNP) calling methods, but evaluating accuracy remains challenging owing to the absence of the SNP gold standard. The definitions for without-gold-standard and performance metrics and their estimation are urgently needed. Additionally, the possible correlations between different SNP loci should also be further explored. To address these challenges, we first introduced the concept of a gold standard and imperfect gold standard under the consistency framework and gave the corresponding definitions of sensitivity and specificity. A latent class model (LCM) was established to estimate the sensitivity and specificity of callers. Furthermore, we incorporated different dependency structures into LCM to investigate their impact on sensitivity and specificity. The performance of LCM was illustrated by comparing the accuracy of BCFtools, DeepVariant, FreeBayes, and GATK on various datasets. Through estimations across multiple datasets, the results indicate that LCM is well-suitable for evaluating callers without the SNP gold standard, and accurate inclusion of the dependency between variations is crucial for better performance ranking. DeepVariant has a higher sum of sensitivity and specificity than other callers, followed by GATK and BCFtools. FreeBayes has low sensitivity but high specificity. Notably, appropriate sequencing coverage is another important factor for precise callers' evaluation. Most importantly, a web interface for assessing and comparing different callers was developed to simplify the evaluation process.

Genome-wide identification of HSP90 gene family in Rosa chinensis and its response to salt and drought stresses

Heat shock protein 90 (HSP90) is important for many organisms, including plants. Based on the whole genome information, the gene number, gene structure, evolutionary relationship, protein structure, and active site of the HSP90 gene family in Rosa chinensis and Rubus idaeus were determined, and the expression of the HSP90 gene under salt, and drought stresses in two rose varieties Wangxifeng and Sweet Avalanche were analyzed. Six and eight HSP90 genes were identified from R. chinensis and Ru. idaeus, respectively. Phylogenetic analysis revealed that the analyzed genes were divided into two Groups and four subgroups (Classes 1a, 1b, 2a, and 2b). Although members within the same classes displayed highly similar gene structures, while the gene structures and conserved domains of Group 1 (Class 1a and 1b) and the Group 2 (Class 2a and 2b) are different. Tandem and segmental duplication genes were found in Ru. idaeus, but not in R. chinensis, perhaps explaining the difference in HSP90 gene quantity in the two analyzed species. Analysis of cis-acting elements revealed abundant abiotic stress, photolight-response, and hormone-response elements in R. chinensis HSP90s. qRT-PCR analysis suggested that RcHSP90-1-1, RcHSP90-5-1 and RcHSP90-6-1 in Sweet Avalanche and Wangxifeng varieties played important regulatory roles under salt and drought stress. The analysis of protein structure and active sites indicate that the potential different roles of RcHSP90-1-1, RcHSP90-5-1, and RcHSP90-6-1 in salt and drought stresses may come from the differences of corresponding protein structures and activation sites. These data will provide information for the breeding of rose varieties with high stress resistance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04052-0.

OrthoRefine: automated enhancement of prior ortholog identification via synteny

BACKGROUND

Identifying orthologs continues to be an early and imperative step in genome analysis but remains a challenging problem. While synteny (conservation of gene order) has previously been used independently and in combination with other methods to identify orthologs, applying synteny in ortholog identification has yet to be automated in a user-friendly manner. This desire for automation and ease-of-use led us to develop OrthoRefine, a standalone program that uses synteny to refine ortholog identification.

RESULTS

We developed OrthoRefine to improve the detection of orthologous genes by implementing a look-around window approach to detect synteny. We tested OrthoRefine in tandem with OrthoFinder, one of the most used software for identification of orthologs in recent years. We evaluated improvements provided by OrthoRefine in several bacterial and a eukaryotic dataset. OrthoRefine efficiently eliminates paralogs from orthologous groups detected by OrthoFinder. Using synteny increased specificity and functional ortholog identification; additionally, analysis of BLAST e-value, phylogenetics, and operon occurrence further supported using synteny for ortholog identification. A comparison of several window sizes suggested that smaller window sizes (eight genes) were generally the most suitable for identifying orthologs via synteny. However, larger windows (30 genes) performed better in datasets containing less closely related genomes. A typical run of OrthoRefine with ~ 10 bacterial genomes can be completed in a few minutes on a regular desktop PC.

CONCLUSION

OrthoRefine is a simple-to-use, standalone tool that automates the application of synteny to improve ortholog detection. OrthoRefine is particularly efficient in eliminating paralogs from orthologous groups delineated by standard methods.

Identification and production of novel potential pathogen-specific biomarkers for diagnosis of histoplasmosis

IMPORTANCE

Histoplasmosis is considered one of the most important mycoses due to the increasing number of individuals susceptible to develop severe clinical forms, particularly those with HIV/AIDS or receiving immunosuppressive biological therapies, the high mortality rates reported when antifungal treatment is not initiated in a timely manner, and the limitations of conventional diagnostic methods. In this context, there is a clear need to improve the capacity of diagnostic tools to specifically detect the fungal pathogen, regardless of the patient's clinical condition or the presence of other co-infections. The proposed novel pathogen-specific biomarkers have the potential to be used in immunodiagnostic platforms and antifungal treatment monitoring in histoplasmosis. In addition, the bioinformatics strategy used in this study could be applied to identify potential diagnostic biomarkers in other models of fungal infection of public health importance.

Genome-wide identification and characterization of GRAS gene family in pigeonpea (Cajanus cajan (L.) Millspaugh)

The GRAS proteins are plant-specific transcription factors (TFs) that play a crucial role in various plant physiological processes, including tissue development and stress responses. To date, GRAS family has been comprehensively characterized in Arabidopsis, soybean, rice, chickpea and other plant species. To understand the structural and functional aspects of pigeonpea (C. cajan), we identified 60 putative GRAS (CcGRAS) genes from pigeonpea genome and further analysed their physicochemical properties, subcellular locations, evolutionary classification, exon-intron structures, conserved domains, gene duplication events and cis-promoter regions. Based on the sequence similarity, CcGRAS family was clustered into 9 subfamilies and the genes with a similar structure and motif distribution were clustered in the same group. The gene duplication studies revealed that these genes were derived from tandem and dispersed duplication events. The cis-promoter regulatory analysis of CcGRAS genes indicated the presence of three types of cis-acting elements including light-responsive, hormone-responsive and plant growth and development related. The expression profiling of CcGRAS genes revealed their tissue-specific functions and differential nature. Collectively, this study highlights relevant functional and regulatory elements of GRAS family in pigeonpea creating a significant resource for future functional studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03782-x.

Evolution of an Alternative Genetic Code in the Providencia Symbiont of the Hematophagous Leech Haementeria acuecueyetzin

Strict blood-feeding animals are confronted with a strong B-vitamin deficiency. Blood-feeding leeches from the Glossiphoniidae family, similarly to hematophagous insects, have evolved specialized organs called bacteriomes to harbor symbiotic bacteria. Leeches of the Haementeria genus have two pairs of globular bacteriomes attached to the esophagus which house intracellular "Candidatus Providencia siddallii" bacteria. Previous work analyzing a draft genome of the Providencia symbiont of the Mexican leech Haementeria officinalis showed that, in this species, the bacteria hold a reduced genome capable of synthesizing B vitamins. In this work, we aimed to expand our knowledge on the diversity and evolution of Providencia symbionts of Haementeria. For this purpose, we sequenced the symbiont genomes of three selected leech species. We found that all genomes are highly syntenic and have kept a stable genetic repertoire, mirroring ancient insect endosymbionts. Additionally, we found B-vitamin pathways to be conserved among these symbionts, pointing to a conserved symbiotic role. Lastly and most notably, we found that the symbiont of H. acuecueyetzin has evolved an alternative genetic code, affecting a portion of its proteome and showing evidence of a lineage-specific and likely intermediate stage of genetic code reassignment.

Major Revisions in Pancrustacean Phylogeny and Evidence of Sensitivity to Taxon Sampling

The clade Pancrustacea, comprising crustaceans and hexapods, is the most diverse group of animals on earth, containing over 80% of animal species and half of animal biomass. It has been the subject of several recent phylogenomic analyses, yet relationships within Pancrustacea show a notable lack of stability. Here, the phylogeny is estimated with expanded taxon sampling, particularly of malacostracans. We show small changes in taxon sampling have large impacts on phylogenetic estimation. By analyzing identical orthologs between two slightly different taxon sets, we show that the differences in the resulting topologies are due primarily to the effects of taxon sampling on the phylogenetic reconstruction method. We compare trees resulting from our phylogenomic analyses with those from the literature to explore the large tree space of pancrustacean phylogenetic hypotheses and find that statistical topology tests reject the previously published trees in favor of the maximum likelihood trees produced here. Our results reject several clades including Caridoida, Eucarida, Multicrustacea, Vericrustacea, and Syncarida. Notably, we find Copepoda nested within Allotriocarida with high support and recover a novel relationship between decapods, euphausiids, and syncarids that we refer to as the Syneucarida. With denser taxon sampling, we find Stomatopoda sister to this latter clade, which we collectively name Stomatocarida, dividing Malacostraca into three clades: Leptostraca, Peracarida, and Stomatocarida. A new Bayesian divergence time estimation is conducted using 13 vetted fossils. We review our results in the context of other pancrustacean phylogenetic hypotheses and highlight 15 key taxa to sample in future studies.

An Overview of PRR- and NLR-Mediated Immunities: Conserved Signaling Components across the Plant Kingdom That Communicate Both Pathways

Cell-surface-localized pattern recognition receptors (PRRs) and intracellular nucleotide-binding domain and leucine-rich repeat receptors (NLRs) are plant immune proteins that trigger an orchestrated downstream signaling in response to molecules of microbial origin or host plant origin. Historically, PRRs have been associated with pattern-triggered immunity (PTI), whereas NLRs have been involved with effector-triggered immunity (ETI). However, recent studies reveal that such binary distinction is far from being applicable to the real world. Although the perception of plant pathogens and the final mounting response are achieved by different means, central hubs involved in signaling are shared between PTI and ETI, blurring the zig-zag model of plant immunity. In this review, we not only summarize our current understanding of PRR- and NLR-mediated immunities in plants, but also highlight those signaling components that are evolutionarily conserved across the plant kingdom. Altogether, we attempt to offer an overview of how plants mediate and integrate the induction of the defense responses that comprise PTI and ETI, emphasizing the need for more evolutionary molecular plant-microbe interactions (EvoMPMI) studies that will pave the way to a better understanding of the emergence of the core molecular machinery involved in the so-called evolutionary arms race between plants and microbes.

Age-Associated Different Transcriptome Profiling in Zebrafish and Rats: an Insight into the Diversity of Vertebrate Aging

Most mammals, including humans, show obvious aging phenotypes, for example, loss of tissue plasticity and sarcopenia. In this regard, fish can be attractive models to study senescence because of their unique aging characteristics. The lifespan of fish varies widely, and several species can live for over 200 years. Moreover, some fish show anti-aging features and indeterminate growth throughout their life. Therefore, exploring the aging mechanism in fish could provide new insights into vertebrate aging. To this end, we conducted RNA sequencing (RNA-seq) assays for various organs and growth stages of zebrafish and compared the data with previously published RNA-seq data of rats. Age-associated differentially expressed genes (DEGs) for all zebrafish tissue samples reveal the upregulation of circadian genes and downregulation of hmgb3a. On one hand, a comparative analysis of DEG profiles associated with aging between zebrafish and rats identifies upregulation of circadian genes and downregulation of collagen genes as conserved transcriptome changes. On the other hand, in zebrafish, upregulation of autophagy-related genes in muscles and AP-1 transcription factor genes in various tissues is observed, which may imply fish-specific anti-aging characteristics. Consistent with our knowledge of mammalian aging, DEG profiles related to tissue senescence are observed in rats. We also detect age-associated downregulation of muscle homeostasis and differentiation-related genes in zebrafish gills, indicating a fish-specific senescence phenotype. Our results indicate both common and different aging profiles between fish and mammals, which could be used for future translational research.

Genome-wide identification and expression characterization of the DoG gene family of moso bamboo (Phyllostachys edulis)

BACKGROUND

The DoG (Delay of Germination1) family plays a key regulatory role in seed dormancy and germination. However, to date, there is no complete genomic overview of the DoG gene family of any economically valuable crop, including moso bamboo (Phyllostachys edulis), and no studies have been conducted to characterize its expression profile. To identify the DoG gene members of moso bamboo (PeDoG) and to investigate their family structural features and tissue expression profile characteristics, a study was conducted. Based on the whole genome and differential transcriptome data, in this investigation, we have scrutinized the physicochemical properties, gene structure, cis-acting elements, phylogenetic relationships, conserved structural (CS) domains, CS motifs and expression patterns of the PeDoG1 family of moso bamboo.

RESULTS

The DoG family genes of moso bamboo were found distributed across 16 chromosomal scaffolds with 24 members. All members were found to carry DoG1 structural domains, while 23 members additionally possessed basic leucine zipper (bZIP) structural domains. We could divide the PeDoG genes into three subfamilies based on phylogenetic relationships. Covariance analysis revealed that tandem duplication was the main driver of amplification of the PeDoG genes. The upstream promoter of these genes containing several cis-acting elements indicates a plausible role in abiotic stress and hormone induction. Gene expression pattern according to transcriptome data revealed participation of the PeDoG genes in tissue and organ development. Analysis using Short Time-series Expression Miner (STEM) tool revealed that the PeDoG gene family is also associated with rapid early shoot growth. Gene ontology (GO) and KEGG analyses showed a dual role of the PeDoG genes. We found that PeDoGs has a possible role as bZIP transcription factors by regulating Polar like1 (PL1) gene expression, and thereby playing a disease response role in moso bamboo. Quantitative gene expression of the PeDoG genes revealed that they were abundantly expressed in roots and leaves, and could be induced in response to gibberellin (GA).

CONCLUSION

In this study, we found that the PeDoG genes are involved in a wide range of activities such as growth and development, stress response and transcription. This forms the first report of PeDoG genes and their potential roles in moso bamboo.

Tracing Eukaryotic Ribosome Biogenesis Factors Into the Archaeal Domain Sheds Light on the Evolution of Functional Complexity

Ribosome assembly is an essential and carefully choreographed cellular process. In eukaryotes, several 100 proteins, distributed across the nucleolus, nucleus, and cytoplasm, co-ordinate the step-wise assembly of four ribosomal RNAs (rRNAs) and approximately 80 ribosomal proteins (RPs) into the mature ribosomal subunits. Due to the inherent complexity of the assembly process, functional studies identifying ribosome biogenesis factors and, more importantly, their precise functions and interplay are confined to a few and very well-established model organisms. Although best characterized in yeast (Saccharomyces cerevisiae), emerging links to disease and the discovery of additional layers of regulation have recently encouraged deeper analysis of the pathway in human cells. In archaea, ribosome biogenesis is less well-understood. However, their simpler sub-cellular structure should allow a less elaborated assembly procedure, potentially providing insights into the functional essentials of ribosome biogenesis that evolved long before the diversification of archaea and eukaryotes. Here, we use a comprehensive phylogenetic profiling setup, integrating targeted ortholog searches with automated scoring of protein domain architecture similarities and an assessment of when search sensitivity becomes limiting, to trace 301 curated eukaryotic ribosome biogenesis factors across 982 taxa spanning the tree of life and including 727 archaea. We show that both factor loss and lineage-specific modifications of factor function modulate ribosome biogenesis, and we highlight that limited sensitivity of the ortholog search can confound evolutionary conclusions. Projecting into the archaeal domain, we find that only few factors are consistently present across the analyzed taxa, and lineage-specific loss is common. While members of the Asgard group are not special with respect to their inventory of ribosome biogenesis factors (RBFs), they unite the highest number of orthologs to eukaryotic RBFs in one taxon. Using large ribosomal subunit maturation as an example, we demonstrate that archaea pursue a simplified version of the corresponding steps in eukaryotes. Much of the complexity of this process evolved on the eukaryotic lineage by the duplication of ribosomal proteins and their subsequent functional diversification into ribosome biogenesis factors. This highlights that studying ribosome biogenesis in archaea provides fundamental information also for understanding the process in eukaryotes.

KinOrtho: a method for mapping human kinase orthologs across the tree of life and illuminating understudied kinases

BACKGROUND

Protein kinases are among the largest druggable family of signaling proteins, involved in various human diseases, including cancers and neurodegenerative disorders. Despite their clinical relevance, nearly 30% of the 545 human protein kinases remain highly understudied. Comparative genomics is a powerful approach for predicting and investigating the functions of understudied kinases. However, an incomplete knowledge of kinase orthologs across fully sequenced kinomes severely limits the application of comparative genomics approaches for illuminating understudied kinases. Here, we introduce KinOrtho, a query- and graph-based orthology inference method that combines full-length and domain-based approaches to map one-to-one kinase orthologs across 17 thousand species.

RESULTS

Using multiple metrics, we show that KinOrtho performed better than existing methods in identifying kinase orthologs across evolutionarily divergent species and eliminated potential false positives by flagging sequences without a proper kinase domain for further evaluation. We demonstrate the advantage of using domain-based approaches for identifying domain fusion events, highlighting a case between an understudied serine/threonine kinase TAOK1 and a metabolic kinase PIK3C2A with high co-expression in human cells. We also identify evolutionary fission events involving the understudied OBSCN kinase domains, further highlighting the value of domain-based orthology inference approaches. Using KinOrtho-defined orthologs, Gene Ontology annotations, and machine learning, we propose putative biological functions of several understudied kinases, including the role of TP53RK in cell cycle checkpoint(s), the involvement of TSSK3 and TSSK6 in acrosomal vesicle localization, and potential functions for the ULK4 pseudokinase in neuronal development.

CONCLUSIONS

In sum, KinOrtho presents a novel query-based tool to identify one-to-one orthologous relationships across thousands of proteomes that can be applied to any protein family of interest. We exploit KinOrtho here to identify kinase orthologs and show that its well-curated kinome ortholog set can serve as a valuable resource for illuminating understudied kinases, and the KinOrtho framework can be extended to any protein-family of interest.

TIAMMAt: Leveraging biodiversity to revise protein domain models, evidence from innate immunity

Sequence annotation is fundamental for studying the evolution of protein families, particularly when working with nonmodel species. Given the rapid, ever-increasing number of species receiving high-quality genome sequencing, accurate domain modeling that is representative of species diversity is crucial for understanding protein family sequence evolution and their inferred function(s). Here, we describe a bioinformatic tool called Taxon-Informed Adjustment of Markov Model Attributes (TIAMMAt) which revises domain profile hidden Markov models (HMMs) by incorporating homologous domain sequences from underrepresented and nonmodel species. Using innate immunity pathways as a case study, we show that revising profile HMM parameters to directly account for variation in homologs among underrepresented species provides valuable insight into the evolution of protein families. Following adjustment by TIAMMAt, domain profile HMMs exhibit changes in their per-site amino acid state emission probabilities and insertion/deletion probabilities while maintaining the overall structure of the consensus sequence. Our results show that domain revision can heavily impact evolutionary interpretations for some families (i.e., NLR’s NACHT domain), whereas impact on other domains (e.g., rel homology domain and interferon regulatory factor domains) is minimal due to high levels of sequence conservation across the sampled phylogenetic depth (i.e., Metazoa). Importantly, TIAMMAt revises target domain models to reflect homologous sequence variation using the taxonomic distribution under consideration by the user. TIAMMAt’s flexibility to revise any subset of the Pfam database using a user-defined taxonomic pool will make it a valuable tool for future protein evolution studies, particularly when incorporating (or focusing) on nonmodel species.

CoreCruncher: Fast and Robust Construction of Core Genomes in Large Prokaryotic Data Sets

The core genome represents the set of genes shared by all, or nearly all, strains of a given population or species of prokaryotes. Inferring the core genome is integral to many genomic analyses, however, most methods rely on the comparison of all the pairs of genomes; a step that is becoming increasingly difficult given the massive accumulation of genomic data. Here, we present CoreCruncher; a program that robustly and rapidly constructs core genomes across hundreds or thousands of genomes. CoreCruncher does not compute all pairwise genome comparisons and uses a heuristic based on the distributions of identity scores to classify sequences as orthologs or paralogs/xenologs. Although it is much faster than current methods, our results indicate that our approach is more conservative than other tools and less sensitive to the presence of paralogs and xenologs. CoreCruncher is freely available from: https://github.com/lbobay/CoreCruncher. CoreCruncher is written in Python 3.7 and can also run on Python 2.7 without modification. It requires the python library Numpy and either Usearch or Blast. Certain options require the programs muscle or mafft.

Identification and Expression Analysis of miR160 and Their Target Genes in Cucumber

miR160 plays a crucial role in various biological processes by regulating their target gene auxin response factor (ARF) in plants. However, little is known about miR160 and ARF in cucumber fruit expansion. Here, 4 Csa-MIR160 family members and 17 CsARFs were identified through a genome-wide search. Csa-miR160 showed a closer relationship with those in melon. Phylogenetic analysis revealed that CsARFs were divided into four classes and most of CsARFs presented a closer evolutionary relationship with those from tomato. Putative cis-elements analysis predicted that Csa-MIR160 and CsARFs were involved in light, phytohormone and stress response, which proved that they might take part in light, phytohormone and stress signaling pathway by the miR160-ARF module. In addition, CsARF5, CsARF11, CsARF13 and CsARF14 were predicted as the target genes of Csa-miR160. qRT-PCR revealed that Csa-miR160 and their target gene CsARFs were differentially expressed in differential cucumber tissues and developmental stages. Csa-miR160d was only expressed in the expanded cucumber fruit. CsARF5, CsARF11 and CsARF13 exhibited the lower expression in the expanded fruit than those in the ovary, while, CsARF14 showed the reverse trend. Our results suggested that Csa-miR160d might play a crucial role in cucumber fruit expansion by negatively targeting CsARF5, CsARF11 and CsARF13. This is the first genome-wide analysis of miR160 in cucumber. These findings provide useful information and resources for further studying the role of miR160 and ARF in cucumber fruit expansion.

Cytoplasmic incompatibility between Old and New World populations of a tramp ant

Reproductive manipulation by endosymbiotic Wolbachia can cause unequal inheritance, allowing the manipulator to spread and potentially impacting evolutionary dynamics in infected hosts. Tramp and invasive species are excellent models to study the dynamics of host-Wolbachia associations because introduced populations often diverge in their microbiomes after colonizing new habitats, resulting in infection polymorphisms between native and introduced populations. Ants are the most abundant group of insects on earth, and numerous ant species are classified as highly invasive. However, little is known about the role of Wolbachia in these ecologically dominant insects. Here, we provide the first description of reproductive manipulation by Wolbachia in an ant. We show that Old and New World populations of the cosmotropic tramp ant Cardiocondyla obscurior harbor distinct Wolbachia strains, and that only the Old World strain manipulates host reproduction by causing cytoplasmic incompatibility (CI) in hybrid crosses. By uncovering a symbiont-induced mechanism of reproductive isolation in a social insect, our study provides a novel perspective on the biology of tramp ants and introduces a new system for studying the evolutionary consequences of CI.

A natural symbiotic bacterium drives mosquito refractoriness to Plasmodium infection via secretion of an antimalarial lipase

The stalling global progress in the fight against malaria prompts the urgent need to develop new intervention strategies. Whilst engineered symbiotic bacteria have been shown to confer mosquito resistance to parasite infection, a major challenge for field implementation is to address regulatory concerns. Here, we report the identification of a Plasmodium-blocking symbiotic bacterium, Serratia ureilytica Su_YN1, isolated from the midgut of wild Anopheles sinensis in China that inhibits malaria parasites via secretion of an antimalarial lipase. Analysis of Plasmodium vivax epidemic data indicates that local malaria cases in Tengchong (Yunnan province, China) are significantly lower than imported cases and importantly, that the local vector A. sinensis is more resistant to infection by P. vivax than A. sinensis from other regions. Analysis of the gut symbiotic bacteria of mosquitoes from Yunnan province led to the identification of S. ureilytica Su_YN1. This bacterium renders mosquitoes resistant to infection by the human parasite Plasmodium falciparum or the rodent parasite Plasmodium berghei via secretion of a lipase that selectively kills parasites at various stages. Importantly, Su_YN1 rapidly disseminates through mosquito populations by vertical and horizontal transmission, providing a potential tool for blocking malaria transmission in the field.

Challenges in benchmarking metagenomic profilers

Accurate microbial identification and abundance estimation are crucial for metagenomics analysis. Various methods for classification of metagenomic data and estimation of taxonomic profiles, broadly referred to as metagenomic profilers, have been developed. Nevertheless, benchmarking of metagenomic profilers remains challenging because some tools are designed to report relative sequence abundance while others report relative taxonomic abundance. Here we show how misleading conclusions can be drawn by neglecting this distinction between relative abundance types when benchmarking metagenomic profilers. Moreover, we show compelling evidence that interchanging sequence abundance and taxonomic abundance will influence both per-sample summary statistics and cross-sample comparisons. We suggest that the microbiome research community pay attention to potentially misleading biological conclusions arising from this issue when benchmarking metagenomic profilers, by carefully considering the type of abundance data that were analyzed and interpreted and clearly stating the strategy used for metagenomic profiling.

Ten Years of Collaborative Progress in the Quest for Orthologs

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

Total Ortholog Median Matrix as an alternative unsupervised approach for phylogenomics based on evolutionary distance between protein coding genes

The increasing number of available genomic data allowed the development of phylogenomic analytical tools. Current methods compile information from single gene phylogenies, whether based on topologies or multiple sequence alignments. Generally, phylogenomic analyses elect gene families or genomic regions to construct phylogenomic trees. Here, we presented an alternative approach for Phylogenomics, named TOMM (Total Ortholog Median Matrix), to construct a representative phylogram composed by amino acid distance measures of all pairwise ortholog protein sequence pairs from desired species inside a group of organisms. The procedure is divided two main steps, (1) ortholog detection and (2) creation of a matrix with the median amino acid distance measures of all pairwise orthologous sequences. We tested this approach within three different group of organisms: Kinetoplastida protozoa, hematophagous Diptera vectors and Primates. Our approach was robust and efficacious to reconstruct the phylogenetic relationships for the three groups. Moreover, novel branch topologies could be achieved, providing insights about some phylogenetic relationships between some taxa.

Comparative Analysis of the Secretome and Interactome of Trypanosoma cruzi and Trypanosoma rangeli Reveals Species Specific Immune Response Modulating Proteins

Chagas disease, a zoonosis caused by the flagellate protozoan Trypanosoma cruzi, is a chronic and systemic parasitic infection that affects ~5–7 million people worldwide, mainly in Latin America. Chagas disease is an emerging public health problem due to the lack of vaccines and effective treatments. According to recent studies, several T. cruzi secreted proteins interact with the human host during cell invasion. Moreover, some comparative studies with T. rangeli, which is non-pathogenic in humans, have been performed to identify proteins directly involved in the pathogenesis of the disease. In this study, we present an integrated analysis of canonical putative secreted proteins (PSPs) from both species. Additionally, we propose an interactome with human host and gene family clusters, and a phylogenetic inference of a selected protein. In total, we identified 322 exclusively PSPs in T. cruzi and 202 in T. rangeli. Among the PSPs identified in T. cruzi, we found several trans-sialidases, mucins, MASPs, proteins with phospholipase 2 domains (PLA2-like), and proteins with Hsp70 domains (Hsp70-like) which have been previously characterized and demonstrated to be related to T. cruzi virulence. PSPs found in T. rangeli were related to protozoan metabolism, specifically carboxylases and phosphatases. Furthermore, we also identified PSPs that may interact with the human immune system, including heat shock and MASP proteins, but in a lower number compared to T. cruzi. Interestingly, we describe a hypothetical hybrid interactome of PSPs which reveals that T. cruzi secreted molecules may be down-regulating IL-17 whilst T. rangeli may enhance the production of IL-15. These results will pave the way for a better understanding of the pathophysiology of Chagas disease and may ultimately lead to the identification of molecular targets, such as key PSPs, that could be used to minimize the health outcomes of Chagas disease by modulating the immune response triggered by T. cruzi infection.

The genetic factors of bilaterian evolution

The Cambrian explosion was a unique animal radiation ~540 million years ago that produced the full range of body plans across bilaterians. The genetic mechanisms underlying these events are unknown, leaving a fundamental question in evolutionary biology unanswered. Using large-scale comparative genomics and advanced orthology evaluation techniques, we identified 157 bilaterian-specific genes. They include the entire Nodal pathway, a key regulator of mesoderm development and left-right axis specification; components for nervous system development, including a suite of G-protein-coupled receptors that control physiology and behaviour, the Robo-Slit midline repulsion system, and the neurotrophin signalling system; a high number of zinc finger transcription factors; and novel factors that previously escaped attention. Contradicting the current view, our study reveals that genes with bilaterian origin are robustly associated with key features in extant bilaterians, suggesting a causal relationship.

The Protector within: Comparative Genomics of APSE Phages across Aphids Reveals Rampant Recombination and Diverse Toxin Arsenals

Phages can fundamentally alter the physiology and metabolism of their hosts. Although these phages are ubiquitous in the bacterial world, they have seldom been described among endosymbiotic bacteria. One notable exception is the APSE phage that is found associated with the gammaproteobacterial Hamiltonella defensa, hosted by several insect species. This secondary facultative endosymbiont is not necessary for the survival of its hosts but can infect certain individuals or even whole populations. Its infection in aphids is often associated with protection against parasitoid wasps. This protective phenotype has actually been linked to the infection of the symbiont strain with an APSE, which carries a toxin cassette that varies among so-called "types." In the present work, we seek to expand our understanding of the diversity of APSE phages as well as the relations of their Hamiltonella hosts. For this, we assembled and annotated the full genomes of 16 APSE phages infecting Hamiltonella symbionts across ten insect species. Molecular and phylogenetic analyses suggest that recombination has occurred repeatedly among lineages. Comparative genomics of the phage genomes revealed two variable regions that are useful for phage typing. Additionally, we find that mobile elements could play a role in the acquisition of new genes in the toxin cassette. Altogether, we provide an unprecedented view of APSE diversity and their genome evolution across aphids. This genomic investigation will provide a valuable resource for the design and interpretation of experiments aiming at understanding the protective phenotype these phages confer to their insect hosts.

Genome-Wide Characterization of GRAS Family and Their Potential Roles in Cold Tolerance of Cucumber (Cucumis sativus L.)

Cucumber (Cucumis sativus L.) is one of the most important cucurbit vegetables but is often subjected to stress during cultivation. GRAS (gibberellic acid insensitive, repressor of GAI, and scarecrow) genes encode a family of transcriptional factors that regulate plant growth and development. In the model plant Arabidopsis thaliana, GRAS family genes function in formation of axillary meristem and root radial structure, phytohormone (gibberellin) signal transduction, light signal transduction and abiotic/biological stress. In this study, a gene family was comprehensively analyzed from the aspects of evolutionary tree, gene structure, chromosome location, evolutionary and expression pattern by means of bioinformatics; 37 GRAS gene family members have been screened from cucumber. We reconstructed an evolutionary tree based on multiple sequence alignment of the typical GRAS domain and conserved motif sequences with those of other species (A. thaliana and Solanum lycopersicum). Cucumber GRAS family was divided into 10 groups according to the classification of Arabidopsis and tomato genes. We conclude that tandem and segmental duplication have played important roles in the expansion and evolution of the cucumber GRAS (CsaGRAS) family. Expression patterns of CsaGRAS genes in different tissues and under cold treatment, combined with gene ontology annotation and interaction network analysis, revealed potentially different functions for CsaGRAS genes in response to cold tolerance, with members of the SHR, SCR and DELLA subfamilies likely playing important roles. In conclusion, this study provides valuable information and candidate genes for improving cucumber tolerance to cold stress.

Whole genome sequence analysis of the mosquitocidal Bacillus thuringiensis LLP29

Bacillus thuringiensis (Bt) is efficient, strongly specific, and avirulent to humans, making it one of the most popular biopesticides in the world. Bt LLP29 is a mosquitocidal strain that was first isolated from Magnolia denudata. To understand its molecular mechanism against mosquitoes, the genome of Bt LLP29 was sequenced and annotated in this study. The LLP29 genome was found to have a total length of 5.99 Mb, with an average G + C content of 35.21%. A total of 6107 coding sequences were also detected, together with 42 rRNAs and 124 tRNAs and 135 other RNAs. With the help of annotation databases, including GO, COG, KEGG, Nr and Swiss-Prot, most unigene functions were identified. At the same time, a collinear analysis was performed on the genome of LLP29. There were also some virulence genes detected, including cry, chitinase, zwittermicin and vip.

Microbial genome analysis: the COG approach

For the past 20 years, the Clusters of Orthologous Genes (COG) database had been a popular tool for microbial genome annotation and comparative genomics. Initially created for the purpose of evolutionary classification of protein families, the COG have been used, apart from straightforward functional annotation of sequenced genomes, for such tasks as (i) unification of genome annotation in groups of related organisms; (ii) identification of missing and/or undetected genes in complete microbial genomes; (iii) analysis of genomic neighborhoods, in many cases allowing prediction of novel functional systems; (iv) analysis of metabolic pathways and prediction of alternative forms of enzymes; (v) comparison of organisms by COG functional categories; and (vi) prioritization of targets for structural and functional characterization. Here we review the principles of the COG approach and discuss its key advantages and drawbacks in microbial genome analysis.

Comparative Genomics of Pediococcus pentosaceus Isolated From Different Niches Reveals Genetic Diversity in Carbohydrate Metabolism and Immune System

Pediococcus pentosaceus isolated from fermented food and the gastrointestinal tracts of humans and animals have been widely identified, and some strains have been reported to reduce inflammation, encephalopathy, obesity and fatty liver in animals. In this study, the genomes of 65 P. pentosaceus strains isolated from human and animal feces and different fermented food were sequenced and comparative genomics analysis was performed on all strains along with nine sequenced representative strains to preliminarily reveal the lifestyle of P. pentosaceus, and investigate the genomic diversity within this species. The results reveal that P. pentosaceus is not host-specific, and shares core genes encoding proteins related to translation, ribosomal structure and biogenesis and signal transduction mechanisms, while its genetic diversity relates mainly to carbohydrate metabolism, and horizontally transferred DNA, especially prophages and bacteriocins encoded on plasmids. Additionally, this is the first report of a type IIA CRISPR/Cas system in P. pentosaceus. This work provides expanded resources of P. pentosaceus genomes, and offers a framework for understanding the biotechnological potential of this species.

Serial horizontal transfer of vitamin-biosynthetic genes enables the establishment of new nutritional symbionts in aphids' di-symbiotic systems

Many insects depend on obligate mutualistic bacteria to provide essential nutrients lacking from their diet. Most aphids, whose diet consists of phloem, rely on the bacterial endosymbiont Buchnera aphidicola to supply essential amino acids and B vitamins. However, in some aphid species, provision of these nutrients is partitioned between Buchnera and a younger bacterial partner, whose identity varies across aphid lineages. Little is known about the origin and the evolutionary stability of these di-symbiotic systems. It is also unclear whether the novel symbionts merely compensate for losses in Buchnera or carry new nutritional functions. Using whole-genome endosymbiont sequences of nine Cinara aphids that harbour an Erwinia-related symbiont to complement Buchnera, we show that the Erwinia association arose from a single event of symbiont lifestyle shift, from a free-living to an obligate intracellular one. This event resulted in drastic genome reduction, long-term genome stasis, and co-divergence with aphids. Fluorescence in situ hybridisation reveals that Erwinia inhabits its own bacteriocytes near Buchnera's. Altogether these results depict a scenario for the establishment of Erwinia as an obligate symbiont that mirrors Buchnera's. Additionally, we found that the Erwinia vitamin-biosynthetic genes not only compensate for Buchnera's deficiencies, but also provide a new nutritional function; whose genes have been horizontally acquired from a Sodalis-related bacterium. A subset of these genes have been subsequently transferred to a new Hamiltonella co-obligate symbiont in one specific Cinara lineage. These results show that the establishment and dynamics of multi-partner endosymbioses can be mediated by lateral gene transfers between co-ocurring symbionts.

Identification and expression analysis of GRAS transcription factors in the wild relative of sweet potato Ipomoea trifida

Background

GRAS gene is an important transcription factor gene family that plays a crucial role in plant growth, development, adaptation to adverse environmental condition. Sweet potato is an important food, vegetable, industrial raw material, and biofuel crop in the world, which plays an essential role in food security in China. However, the function of sweet potato GRAS genes remains unknown.

Results

In this study, we identified and characterised 70 GRAS members from Ipomoea trifida, which is the progenitor of sweet potato. The chromosome distribution, phylogenetic tree, exon-intron structure and expression profiles were analysed. The distribution map showed that GRAS genes were randomly located in 15 chromosomes. In combination with phylogenetic analysis and previous reports in Arabidopsis and rice, the GRAS proteins from I. trifida were divided into 11 subfamilies. Gene structure showed that most of the GRAS genes in I. trifida lacked introns. The tissue-specific expression patterns and the patterns under abiotic stresses of ItfGRAS genes were investigated via RNA-seq and further tested by RT-qPCR. Results indicated the potential functions of ItfGRAS during plant development and stress responses.

Conclusions

Our findings will further facilitate the functional study of GRAS gene and molecular breeding of sweet potato.

Genome-wide identification and characterization of the MADS-box gene family in Salix suchowensis

MADS-box genes encode transcription factors that participate in various plant growth and development processes, particularly floral organogenesis. To date, MADS-box genes have been reported in many species, the completion of the sequence of the willow genome provides us with the opportunity to conduct a comprehensive analysis of the willow MADS-box gene family. Here, we identified 60 willow MADS-box genes using bioinformatics-based methods and classified them into 22 M-type (11 Mα, seven Mβ and four Mγ) and 38 MIKC-type (32 MIKCc and six MIKC*) genes based on a phylogenetic analysis. Fifty-six of the 60 SsMADS genes were randomly distributed on 19 putative willow chromosomes. By combining gene structure analysis with evolutionary analysis, we found that the MIKC-type genes were more conserved and played a more important role in willow growth. Further study showed that the MIKC* type was a transition between the M-type and MIKC-type. Additionally, the number of MADS-box genes in gymnosperms was notably lower than that in angiosperms. Finally, the expression profiles of these willow MADS-box genes were analysed in five different tissues (root, stem, leave, bud and bark) and validated by RT-qPCR experiments. This study is the first genome-wide analysis of the willow MADS-box gene family, and the results establish a basis for further functional studies of willow MADS-box genes and serve as a reference for related studies of other woody plants.

SonicParanoid: fast, accurate and easy orthology inference

Motivation

Orthology inference constitutes a common base of many genome-based studies, as a pre-requisite for annotating new genomes, finding target genes for biotechnological applications and revealing the evolutionary history of life. Although its importance keeps rising with the ever-growing number of sequenced genomes, existing tools are computationally demanding and difficult to employ.

Results

Here, we present SonicParanoid, which is faster than, but comparably accurate to, the well-established tools with a balanced precision-recall trade-off. Furthermore, SonicParanoid substantially relieves the difficulties of orthology inference for those who need to construct and maintain their own genomic datasets.

Availability and implementation

SonicParanoid is available with a GNU GPLv3 license on the Python Package Index and BitBucket. Documentation is available at http://iwasakilab.bs.s.u-tokyo.ac.jp/sonicparanoid.

Supplementary information

Supplementary data are available at Bioinformatics online.

A Co-Expression Network in Hexaploid Wheat Reveals Mostly Balanced Expression and Lack of Significant Gene Loss of Homeologous Meiotic Genes Upon Polyploidization

Polyploidization has played an important role in plant evolution. However, upon polyploidization, the process of meiosis must adapt to ensure the proper segregation of increased numbers of chromosomes to produce balanced gametes. It has been suggested that meiotic gene (MG) duplicates return to a single copy following whole genome duplication to stabilize the polyploid genome. Therefore, upon the polyploidization of wheat, a hexaploid species with three related (homeologous) genomes, the stabilization process may have involved rapid changes in content and expression of MGs on homeologous chromosomes (homeologs). To examine this hypothesis, sets of candidate MGs were identified in wheat using co-expression network analysis and orthology informed approaches. In total, 130 RNA-Seq samples from a range of tissues including wheat meiotic anthers were used to define co-expressed modules of genes. Three modules were significantly correlated with meiotic tissue samples but not with other tissue types. These modules were enriched for GO terms related to cell cycle, DNA replication, and chromatin modification and contained orthologs of known MGs. Overall, 74.4% of genes within these meiosis-related modules had three homeologous copies which was similar to other tissue-related modules. Amongst wheat MGs identified by orthology, rather than co-expression, the majority (93.7%) were either retained in hexaploid wheat at the same number of copies (78.4%) or increased in copy number (15.3%) compared to ancestral wheat species. Furthermore, genes within meiosis-related modules showed more balanced expression levels between homeologs than genes in non-meiosis-related modules. Taken together, our results do not support extensive gene loss nor changes in homeolog expression of MGs upon wheat polyploidization. The construction of the MG co-expression network allowed identification of hub genes and provided key targets for future studies.

SwiftOrtho: A fast, memory-efficient, multiple genome orthology classifier

BACKGROUND

Gene homology type classification is required for many types of genome analyses, including comparative genomics, phylogenetics, and protein function annotation. Consequently, a large variety of tools have been developed to perform homology classification across genomes of different species. However, when applied to large genomic data sets, these tools require high memory and CPU usage, typically available only in computational clusters.

FINDINGS

Here we present a new graph-based orthology analysis tool, SwiftOrtho, which is optimized for speed and memory usage when applied to large-scale data. SwiftOrtho uses long k-mers to speed up homology search, while using a reduced amino acid alphabet and spaced seeds to compensate for the loss of sensitivity due to long k-mers. In addition, it uses an affinity propagation algorithm to reduce the memory usage when clustering large-scale orthology relationships into orthologous groups. In our tests, SwiftOrtho was the only tool that completed orthology analysis of proteins from 1,760 bacterial genomes on a computer with only 4 GB RAM. Using various standard orthology data sets, we also show that SwiftOrtho has a high accuracy.

CONCLUSIONS

SwiftOrtho enables the accurate comparative genomic analyses of thousands of genomes using low-memory computers. SwiftOrtho is available at https://github.com/Rinoahu/SwiftOrtho.

Natural selection and repeated patterns of molecular evolution following allopatric divergence

Although geographic isolation is a leading driver of speciation, the tempo and pattern of divergence at the genomic level remain unclear. We examine genome-wide divergence of putatively single-copy orthologous genes (POGs) in 20 allopatric species/variety pairs from diverse angiosperm clades, with 16 pairs reflecting the classic eastern Asia-eastern North America floristic disjunction. In each pair, >90% of POGs are under purifying selection, and <10% are under positive selection. A set of POGs are under strong positive selection, 14 of which are shared by 10-15 pairs, and one shared by all pairs; 15 POGs are annotated to biological processes responding to various stimuli. The relative abundance of POGs under different selective forces exhibits a repeated pattern among pairs despite an ~10 million-year difference in divergence time. Species divergence times are positively correlated with abundance of POGs under moderate purifying selection, but negatively correlated with abundance of POGs under strong purifying selection.

Paternally Expressed Imprinted Genes under Positive Darwinian Selection in Arabidopsis thaliana

Genomic imprinting is an epigenetic phenomenon where autosomal genes display uniparental expression depending on whether they are maternally or paternally inherited. Genomic imprinting can arise from parental conflicts over resource allocation to the offspring, which could drive imprinted loci to evolve by positive selection. We investigate whether positive selection is associated with genomic imprinting in the inbreeding species Arabidopsis thaliana. Our analysis of 140 genes regulated by genomic imprinting in the A. thaliana seed endosperm demonstrates they are evolving more rapidly than expected. To investigate whether positive selection drives this evolutionary acceleration, we identified orthologs of each imprinted gene across 34 plant species and elucidated their evolutionary trajectories. Increased positive selection was sought by comparing its incidence among imprinted genes with nonimprinted controls. Strikingly, we find a statistically significant enrichment of imprinted paternally expressed genes (iPEGs) evolving under positive selection, 50.6% of the total, but no such enrichment for positive selection among imprinted maternally expressed genes (iMEGs). This suggests that maternally- and paternally expressed imprinted genes are subject to different selective pressures. Almost all positively selected amino acids were fixed across 80 sequenced A. thaliana accessions, suggestive of selective sweeps in the A. thaliana lineage. The imprinted genes under positive selection are involved in processes important for seed development including auxin biosynthesis and epigenetic regulation. Our findings support a genomic imprinting model for plants where positive selection can affect paternally expressed genes due to continued conflict with maternal sporophyte tissues, even when parental conflict is reduced in predominantly inbreeding species.

CAARS: comparative assembly and annotation of RNA-Seq data

MOTIVATION

RNA sequencing (RNA-Seq) is a widely used approach to obtain transcript sequences in non-model organisms, notably for performing comparative analyses. However, current bioinformatic pipelines do not take full advantage of pre-existing reference data in related species for improving RNA-Seq assembly, annotation and gene family reconstruction.

RESULTS

We built an automated pipeline named CAARS to combine novel data from RNA-Seq experiments with existing multi-species gene family alignments. RNA-Seq reads are assembled into transcripts by both de novo and assisted assemblies. Then, CAARS incorporates transcripts into gene families, builds gene alignments and trees and uses phylogenetic information to classify the genes as orthologs and paralogs of existing genes. We used CAARS to assemble and annotate RNA-Seq data in rodents and fishes using distantly related genomes as reference, a difficult case for this kind of analysis. We showed CAARS assemblies are more complete and accurate than those assembled by a standard pipeline consisting of de novo assembly coupled with annotation by sequence similarity on a guide species. In addition to annotated transcripts, CAARS provides gene family alignments and trees, annotated with orthology relationships, directly usable for downstream comparative analyses.

AVAILABILITY AND IMPLEMENTATION

CAARS is implemented in Python and Ocaml and is freely available at https://github.com/carinerey/caars.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Hidden genomic evolution in a morphospecies-The landscape of rapidly evolving genes in Tetrahymena

A morphospecies is defined as a taxonomic species based wholly on morphology, but often morphospecies consist of clusters of cryptic species that can be identified genetically or molecularly. The nature of the evolutionary novelty that accompanies speciation in a morphospecies is an intriguing question. Morphospecies are particularly common among ciliates, a group of unicellular eukaryotes that separates 2 kinds of nuclei—the silenced germline nucleus (micronucleus [MIC]) and the actively expressed somatic nucleus (macronucleus [MAC])—within a common cytoplasm. Because of their very similar morphologies, members of the Tetrahymena genus are considered a morphospecies. We explored the hidden genomic evolution within this genus by performing a comprehensive comparative analysis of the somatic genomes of 10 species and the germline genomes of 2 species of Tetrahymena. These species show high genetic divergence; phylogenomic analysis suggests that the genus originated about 300 million years ago (Mya). Seven universal protein domains are preferentially included among the species-specific (i.e., the youngest) Tetrahymena genes. In particular, leucine-rich repeat (LRR) genes make the largest contribution to the high level of genome divergence of the 10 species. LRR genes can be sorted into 3 different age groups. Parallel evolutionary trajectories have independently occurred among LRR genes in the different Tetrahymena species. Thousands of young LRR genes contain tandem arrays of exactly 90-bp exons. The introns separating these exons show a unique, extreme phase 2 bias, suggesting a clonal origin and successive expansions of 90-bp–exon LRR genes. Identifying LRR gene age groups allowed us to document a Tetrahymena intron length cycle. The youngest 90-bp exon LRR genes in T. thermophila are concentrated in pericentromeric and subtelomeric regions of the 5 micronuclear chromosomes, suggesting that these regions act as genome innovation centers. Copies of a Tetrahymena Long interspersed element (LINE)-like retrotransposon are very frequently found physically adjacent to 90-bp exon/intron repeat units of the youngest LRR genes. We propose that Tetrahymena species have used a massive exon-shuffling mechanism, involving unequal crossing over possibly in concert with retrotransposition, to create the unique 90-bp exon array LRR genes.

Genomic comparison of ten morphologically very similar species of ciliate from the genus Tetrahymena reveals how parallel microevolutionary processes have shaped their genomes and created unique genes through retrotransposition.

Comparative Transcriptome Analysis in Eggplant Reveals Selection Trends during Eggplant Domestication

Eggplant (Solanum melongena L.) is an economically and nutritionally important fruit crop of the Solanaceae family, which was domesticated in India and southern China. However, the genome regions subjected to selective sweeps in eggplant remain unknown. In the present study, we performed comparative transcriptome analysis of cultivated and wild eggplant species with emphasis on the selection pattern during domestication. In total, 44,073 (S. sisymbriifolium) to 58,677 (S. melongena cultivar S58) unigenes were generated for the six eggplant accessions with total lengths of 36.6-46 Mb. The orthologous genes were assessed using the ratio of nonsynonymous (K a) to synonymous (K s) nucleotide substitutions to characterize selective patterns during eggplant domestication. We identified 19 genes under positive selection across the phylogeny that were classified into four groups. The gene (OG12205) under positive selection was possibly associated with fruit-related traits in eggplant, which may have resulted from human manipulation. Eight positive selected genes were potentially involved in stress tolerance or disease resistance, suggesting that environmental changes and biotic stresses were important selective pressures in eggplant domestication. Taken together, our results shed light on the effects of artificial and natural selection on the transcriptomes of eggplant and its wild relatives. Identification of the selected genes will facilitate the understanding of genetic architecture of domesticated-related traits and provide resources for resistant breeding in eggplant.

The Evolutionary Traceability of a Protein

Orthologs document the evolution of genes and metabolic capacities encoded in extant and ancient genomes. However, the similarity between orthologs decays with time, and ultimately it becomes insufficient to infer common ancestry. This leaves ancient gene set reconstructions incomplete and distorted to an unknown extent. Here we introduce the “evolutionary traceability” as a measure that quantifies, for each protein, the evolutionary distance beyond which the sensitivity of the ortholog search becomes limiting. Using yeast, we show that genes that were thought to date back to the last universal common ancestor are of high traceability. Their functions mostly involve catalysis, ion transport, and ribonucleoprotein complex assembly. In turn, the fraction of yeast genes whose traceability is not sufficient to infer their presence in last universal common ancestor is enriched for regulatory functions. Computing the traceabilities of genes that have been experimentally characterized as being essential for a self-replicating cell reveals that many of the genes that lack orthologs outside bacteria have low traceability. This leaves open whether their orthologs in the eukaryotic and archaeal domains have been overlooked. Looking at the example of REC8, a protein essential for chromosome cohesion, we demonstrate how a traceability-informed adjustment of the search sensitivity identifies hitherto missed orthologs in the fast-evolving microsporidia. Taken together, the evolutionary traceability helps to differentiate between true absence and nondetection of orthologs, and thus improves our understanding about the evolutionary conservation of functional protein networks. “protTrace,” a software tool for computing evolutionary traceability, is freely available at https://github.com/BIONF/protTrace.git; last accessed February 10, 2019.

Evolutionary Patterns of Non-Coding RNA in Cardiovascular Biology

Cardiovascular diseases (CVDs) affect the heart and the vascular system with a high prevalence and place a huge burden on society as well as the healthcare system. These complex diseases are often the result of multiple genetic and environmental risk factors and pose a great challenge to understanding their etiology and consequences. With the advent of next generation sequencing, many non-coding RNA transcripts, especially long non-coding RNAs (lncRNAs), have been linked to the pathogenesis of CVD. Despite increasing evidence, the proper functional characterization of most of these molecules is still lacking. The exploration of conservation of sequences across related species has been used to functionally annotate protein coding genes. In contrast, the rapid evolutionary turnover and weak sequence conservation of lncRNAs make it difficult to characterize functional homologs for these sequences. Recent studies have tried to explore other dimensions of interspecies conservation to elucidate the functional role of these novel transcripts. In this review, we summarize various methodologies adopted to explore the evolutionary conservation of cardiovascular non-coding RNAs at sequence, secondary structure, syntenic, and expression level.

Relation of the pdxB-usg-truA-dedA Operon and the truA Gene to the Intracellular Survival of Salmonella enterica Serovar Typhimurium

Salmonella is the genus of Gram-negative, facultative intracellular pathogens that have the ability to infect large numbers of animal or human hosts. The S. enterica usg gene is associated with intracellular survival based on ortholog screening and identification. In this study, the λ-Red recombination system was used to construct gene deletion strains and to investigate whether the identified operon was related to intracellular survival. The pdxB-usg-truA-dedA operon enhanced the intracellular survival of S. enterica by resisting the oxidative environment and the usg and truA gene expression was induced by H₂O₂. Moreover, the genes in this operon (except for dedA) contributed to virulence in mice. These findings indicate that the pdxB-usg-truA-dedA operon functions in resistance to oxidative environments during intracellular survival and is required for in vivo S. enterica virulence. This study provides insight toward a better understand of the characteristics of intracellular pathogens and explores the gene modules involved in their intracellular survival.

Inferring Orthology and Paralogy

The distinction between orthologs and paralogs, genes that started diverging by speciation versus duplication, is relevant in a wide range of contexts, most notably phylogenetic tree inference and protein function annotation. In this chapter, we provide an overview of the methods used to infer orthology and paralogy. We survey both graph-based approaches (and their various grouping strategies) and tree-based approaches, which solve the more general problem of gene/species tree reconciliation. We discuss conceptual differences among the various orthology inference methods and databases and examine the difficult issue of verifying and benchmarking orthology predictions. Finally, we review typical applications of orthologous genes, groups, and reconciled trees and conclude with thoughts on future methodological developments.

Molecular docking and in silico studies of the physicochemical properties of potential inhibitors for the phosphotransferase system of Streptococcus mutans

This study identified potential inhibitory compounds of the phosphoenolpyruvate-sugar. Phosphotransferase system of S. mutans, specifically enzyme II mannose transporter (EII^Man) in its subunits IIA, IIB and IIC by means of a selection protocol and in silico molecular analysis. Intervening the phosphotransferase system would compromise the physiological behavior and the pathogenic expression of S. mutans, and possibly other acidogenic bacteria that use phosphotransferases in their metabolism-making the phosphotransferase system a therapeutic target for the selective control of acidogenic microorganisms in caries control. Several computational techniques were used to evaluate molecular, physicochemical, and toxicological aspects of various compounds. Molecular docking was used to calculate the binding potential (ΔG) between receptor protein subunits and more than 836,000 different chemical compounds from the ZINC database. Physicochemical parameters related to the compounds' pharmacokinetic and pharmacodynamic indicators were evaluated, including absorption, distribution, metabolism, excretion, and toxicity (ADMET), and chemical analysis characterized the compounds structures. Thirteen compounds with EII binding potential of the phosphotransferase system of S. mutans and favorable ADMET properties were identified. Six spirooxindoles and three pyrrolidones stand out from the found compounds; unique structural characteristics of spirooxindoles and pyrrolidones associated with various reported biological activities like anti-microbial, antiinflammatory, anticancer, nootropic, neuroprotective and antiepileptic effects, among other pharmacological effects with surprising differences in terms of mechanisms of action. Following studies will provide more evidence of the action of these compounds on the phosphotransferase system of S. mutans, and its possible applications.

Comparative analysis of selection mode reveals different evolutionary rate and expression pattern in Arachis duranensis and Arachis ipaënsis duplicated genes

Our results reveal that K_s is a determining factor affecting selective pressure and different evolution and expression patterns are detected between PSGs and NSGs in wild Arachis duplicates. Selective pressure, including purifying (negative) and positive selection, can be detected in organisms. However, studies on comparative evolutionary rates, gene expression patterns and gene features between negatively selected genes (NSGs) and positively selected genes (PSGs) are lagging in paralogs of plants. Arachis duranensis and Arachis ipaënsis are ancestors of the cultivated peanut, an important oil and protein crop. Here, we carried out a series of systematic analyses, comparing NSG and PSG in paralogs, using genome sequences and transcriptome datasets in A. duranensis and A. ipaënsis. We found that synonymous substitution rate (K_s) is a determining factor affecting selective pressure in A. duranensis and A. ipaënsis duplicated genes. Lower expression level, lower gene expression breadth, higher codon bias and shorter polypeptide length were found in PSGs and not in NSGs. The correlation analyses showed that gene expression breadth was positively correlated with polypeptide length and GC content at the first codon site (GC1) in PSGs and NSGs, respectively. There was a negative correlation between expression level and polypeptide length in PSGs. In NSGs, the K_s was positively correlated with expression level, gene expression breadth, GC1, and GC content at the third codon site (GC3), but selective pressure was negatively correlated with expression level, gene expression breadth, polypeptide length, GC1, and GC3 content. The function of most duplicated gene pairs was divergent under drought and nematode stress. Taken together, our results show that different evolution and expression patterns occur between PSGs and NSGs in paralogs of two wild Arachis species.

In silico characterization and expression profiling of the diacylglycerol acyltransferase gene family (DGAT1, DGAT2, DGAT3 and WS/DGAT) from oil palm, Elaeis guineensis

The diacylglycerol acyltransferases (DGAT) (diacylglycerol:acyl-CoA acyltransferase, EC 2.3.1.20) are a key group of enzymes that catalyse the final and usually the most important rate-limiting step of triacylglycerol biosynthesis in plants and other organisms. Genes encoding four distinct functional families of DGAT enzymes have been characterised in the genome of the African oil palm, Elaeis guineensis. The contrasting features of the various isoforms within the four families of DGAT genes, namely DGAT1, DGAT2, DGAT3 and WS/DGAT are presented both in the oil palm itself and, for comparative purposes, in 12 other oil crop or model/related plants, namely Arabidopsis thaliana, Brachypodium distachyon, Brassica napus, Elaeis oleifera, Glycine max, Gossypium hirsutum, Helianthus annuus, Musa acuminata, Oryza sativa, Phoenix dactylifera, Sorghum bicolor, and Zea mays. The oil palm genome contains respectively three, two, two and two distinctly expressed functional copies of the DGAT1, DGAT2, DGAT3 and WS/DGAT genes. Phylogenetic analyses of the four DGAT families showed that the E. guineensis genes tend to cluster with sequences from P. dactylifera and M. acuminata rather than with other members of the Commelinid monocots group, such as the Poales which include the major cereal crops such as rice and maize. Comparison of the predicted DGAT protein sequences with other animal and plant DGATs was consistent with the E. guineensis DGAT1 being ER located with its active site facing the lumen while DGAT2, although also ER located, had a predicted cytosol-facing active site. In contrast, DGAT3 and some (but not all) WS/DGAT in E. guineensis are predicted to be soluble, cytosolic enzymes. Evaluation of E. guineensis DGAT gene expression in different tissues and developmental stages suggests that the four DGAT groups have distinctive physiological roles and are particularly prominent in developmental processes relating to reproduction, such as flowering, and in fruit/seed formation especially in the mesocarp and endosperm tissues.

A Freeloader? The Highly Eroded Yet Large Genome of the Serratia symbiotica Symbiont of Cinara strobi

Genome reduction is pervasive among maternally inherited bacterial endosymbionts. This genome reduction can eventually lead to serious deterioration of essential metabolic pathways, thus rendering an obligate endosymbiont unable to provide essential nutrients to its host. This loss of essential pathways can lead to either symbiont complementation (sharing of the nutrient production with a novel co-obligate symbiont) or symbiont replacement (complete takeover of nutrient production by the novel symbiont). However, the process by which these two evolutionary events happen remains somewhat enigmatic by the lack of examples of intermediate stages of this process. Cinara aphids (Hemiptera: Aphididae) typically harbor two obligate bacterial symbionts: Buchnera and Serratia symbiotica. However, the latter has been replaced by different bacterial taxa in specific lineages, and thus species within this aphid lineage could provide important clues into the process of symbiont replacement. In the present study, using 16S rRNA high-throughput amplicon sequencing, we determined that the aphid Cinara strobi harbors not two, but three fixed bacterial symbionts: Buchnera aphidicola, a Sodalis sp., and S. symbiotica. Through genome assembly and genome-based metabolic inference, we have found that only the first two symbionts (Buchnera and Sodalis) actually contribute to the hosts' supply of essential nutrients while S. symbiotica has become unable to contribute towards this task. We found that S. symbiotica has a rather large and highly eroded genome which codes only for a few proteins and displays extensive pseudogenization. Thus, we propose an ongoing symbiont replacement within C. strobi, in which a once "competent" S. symbiotica does no longer contribute towards the beneficial association. These results suggest that in dual symbiotic systems, when a substitute cosymbiont is available, genome deterioration can precede genome reduction and a symbiont can be maintained despite the apparent lack of benefit to its host.

OrthoList 2: A New Comparative Genomic Analysis of Human and Caenorhabditis elegans Genes

OrthoList, a compendium of Caenorhabditis elegans genes with human orthologs compiled in 2011 by a meta-analysis of four orthology-prediction methods, has been a popular tool for identifying conserved genes for research into biological and disease mechanisms. However, the efficacy of orthology prediction depends on the accuracy of gene-model predictions, an ongoing process, and orthology-prediction algorithms have also been updated over time. Here we present OrthoList 2 (OL2), a new comparative genomic analysis between C. elegans and humans, and the first assessment of how changes over time affect the landscape of predicted orthologs between two species. Although we find that updates to the orthology-prediction methods significantly changed the landscape of C. elegans-human orthologs predicted by individual programs and-unexpectedly-reduced agreement among them, we also show that our meta-analysis approach "buffered" against changes in gene content. We show that adding results from more programs did not lead to many additions to the list and discuss reasons to avoid assigning "scores" based on support by individual orthology-prediction programs; the treatment of "legacy" genes no longer predicted by these programs; and the practical difficulties of updating due to encountering deprecated, changed, or retired gene identifiers. In addition, we consider what other criteria may support claims of orthology and alternative approaches to find potential orthologs that elude identification by these programs. Finally, we created a new web-based tool that allows for rapid searches of OL2 by gene identifiers, protein domains [InterPro and SMART (Simple Modular Architecture Research Tool], or human disease associations ([OMIM (Online Mendelian Inheritence in Man], and also includes available RNA-interference resources to facilitate potential translational cross-species studies.

Refining a steroidogenic model: an analysis of RNA-seq datasets from insect prothoracic glands

BACKGROUND

The prothoracic gland (PG), the principal steroidogenic organ of insects, has been proposed as a model for steroid hormone biosynthesis and regulation.

RESULTS

To validate the robustness of the model, we present an analysis of accumulated transcriptomic data from PGs of two model species, Drosophila melanogaster and Bombyx mori. We identify that the common core components of the model in both species are encoded by nine genes. Five of these are Halloween genes whose expression differs substantially between the PGs of these species.

CONCLUSIONS

We conclude that the PGs can be a model for steroid hormone synthesis and regulation within the context of mitochondrial cholesterol transport and steroid biosynthesis but beyond these core mechanisms, gene expression in insect PGs is too diverse to fit in a context-specific model and should be analysed within a species-specific framework.

HiMAP: Robust phylogenomics from highly multiplexed amplicon sequencing

High-throughput sequencing has fundamentally changed how molecular phylogenetic data sets are assembled, and phylogenomic data sets commonly contain 50- to 100-fold more loci than those generated using traditional Sanger sequencing-based approaches. Here, we demonstrate a new approach for building phylogenomic data sets using single-tube, highly multiplexed amplicon sequencing, which we name HiMAP (highly multiplexed amplicon-based phylogenomics) and present bioinformatic pipelines for locus selection based on genomic and transcriptomic data resources and postsequencing consensus calling and alignment. This method is inexpensive and amenable to sequencing a large number (hundreds) of taxa simultaneously and requires minimal hands-on time at the bench (<1/2 day), and data analysis can be accomplished without the need for read mapping or assembly. We demonstrate this approach by sequencing 878 amplicons in single reactions for 82 species of tephritid fruit flies across seven genera (384 individuals), including some of the most economically important agricultural insect pests. The resulting filtered data set (>150,000-bp concatenated alignment, ~20% missing character sites across all individuals and amplicons) contained >40,000 phylogenetically informative characters, and although some discordance was observed between analyses, it provided unparalleled resolution of many phylogenetic relationships in this group. Most notably, we found high support for the generic status of Zeugodacus and the sister relationship between Dacus and Zeugodacus. We discuss HiMAP, with regard to its molecular and bioinformatic strengths, and the insight the resulting data set provides into relationships of this diverse insect group.