Automatic clustering of orthologs and inparalogs shared by multiple proteomes

Computational tools for plant genomics and breeding

Plant genomics and crop breeding are at the intersection of biotechnology and information technology. Driven by a combination of high-throughput sequencing, molecular biology and data science, great advances have been made in omics technologies at every step along the central dogma, especially in genome assembling, genome annotation, epigenomic profiling, and transcriptome profiling. These advances further revolutionized three directions of development. One is genetic dissection of complex traits in crops, along with genomic prediction and selection. The second is comparative genomics and evolution, which open up new opportunities to depict the evolutionary constraints of biological sequences for deleterious variant discovery. The third direction is the development of deep learning approaches for the rational design of biological sequences, especially proteins, for synthetic biology. All three directions of development serve as the foundation for a new era of crop breeding where agronomic traits are enhanced by genome design.

Comparison of gene clustering criteria reveals intrinsic uncertainty in pangenome analyses

BACKGROUND

A key step for comparative genomics is to group open reading frames into functionally and evolutionarily meaningful gene clusters. Gene clustering is complicated by intraspecific duplications and horizontal gene transfers that are frequent in prokaryotes. In consequence, gene clustering methods must deal with a trade-off between identifying vertically transmitted representatives of multicopy gene families, which are recognizable by synteny conservation, and retrieving complete sets of species-level orthologs. We studied the implications of adopting homology, orthology, or synteny conservation as formal criteria for gene clustering by performing comparative analyses of 125 prokaryotic pangenomes.

RESULTS

Clustering criteria affect pangenome functional characterization, core genome inference, and reconstruction of ancestral gene content to different extents. Species-wise estimates of pangenome and core genome sizes change by the same factor when using different clustering criteria, allowing robust cross-species comparisons regardless of the clustering criterion. However, cross-species comparisons of genome plasticity and functional profiles are substantially affected by inconsistencies among clustering criteria. Such inconsistencies are driven not only by mobile genetic elements, but also by genes involved in defense, secondary metabolism, and other accessory functions. In some pangenome features, the variability attributed to methodological inconsistencies can even exceed the effect sizes of ecological and phylogenetic variables.

CONCLUSIONS

Choosing an appropriate criterion for gene clustering is critical to conduct unbiased pangenome analyses. We provide practical guidelines to choose the right method depending on the research goals and the quality of genome assemblies, and a benchmarking dataset to assess the robustness and reproducibility of future comparative studies.

Genome-wide identification and analysis of a cotton secretome reveals its role in resistance against Verticillium dahliae

BACKGROUND

The extracellular space between the cell wall and plasma membrane is a battlefield in plant-pathogen interactions. Within this space, the pathogen employs its secretome to attack the host in a variety of ways, including immunity manipulation. However, the role of the plant secretome is rarely studied for its role in disease resistance.

RESULTS

Here, we examined the secretome of Verticillium wilt-resistant Gossypium hirsutum cultivar Zhongzhimian No.2 (ZZM2, encoding 95,327 predicted coding sequences) to determine its role in disease resistance against the wilt causal agent, Verticillium dahliae. Bioinformatics-driven analyses showed that the ZZM2 genome encodes 2085 secreted proteins and that these display disequilibrium in their distribution among the chromosomes. The cotton secretome displayed differences in the abundance of certain amino acid residues as compared to the remaining encoded proteins due to the localization of these putative proteins in the extracellular space. The secretome analysis revealed conservation for an allotetraploid genome, which nevertheless exhibited variation among orthologs and comparable unique genes between the two sub-genomes. Secretome annotation strongly suggested its involvement in extracellular stress responses (hydrolase activity, oxidoreductase activity, and extracellular region, etc.), thus contributing to resistance against the V. dahliae infection. Furthermore, the defense response genes (immunity marker NbHIN1, salicylic acid marker NbPR1, and jasmonic acid marker NbLOX4) were activated to varying degrees when Nicotina benthamiana leaves were agro-infiltrated with 28 randomly selected members, suggesting that the secretome plays an important role in the immunity response. Finally, gene silencing assays of 11 members from 13 selected candidates in ZZM2 displayed higher susceptibility to V. dahliae, suggesting that the secretome members confer the Verticillium wilt resistance in cotton.

CONCLUSIONS

Our data demonstrate that the cotton secretome plays an important role in Verticillium wilt resistance, facilitating the development of the resistance gene markers and increasing the understanding of the mechanisms regulating disease resistance.

NB-LRR Lineage-Specific Equipment Is Sorted Out by Sequence Pattern Adaptation and Domain Segment Shuffling

The nucleotide-binding and leucine-rich repeat (NB-LRR) genes, also known as resistance (R)-genes, play an important role in the activation of immune responses. In recent years, large-scale studies have been performed to highlight the diversification of plant NB-LRR repertories. It is well known that, to provide new functionalities, NB-LRR sequences are subject to duplication, domain fusions and acquisition and other kinds of mutations. Although some mechanisms that govern NB-LRR protein domain adaptations have been uncovered, to retrace the plant-lineage-specific evolution routes of R protein structure, a multi-genome comparative analysis was performed. This study allowed us to define groups of genes sharing homology relationships across different species. It is worth noting that the most populated groups contained well-characterized R proteins. The arsenal profile of such groups was investigated in five botanical families, including important crop species, to underline specific adaptation signatures. In addition, the dissection of 70 NB domains of well-characterized R-genes revealed the NB core motifs from which the three main R protein classes have been diversified. The structural remodeling of domain segments shaped the specific NB-LRR repertoires observed in each plant species. This analysis provided new evolutionary and functional insights on NB protein domain shuffling. Taken together, such findings improved our understanding of the molecular adaptive selection mechanisms occurring at plant R loci.

Genomic and transcriptomic analyses of the subterranean termite Reticulitermes speratus: Gene duplication facilitates social evolution

Gene duplication is a major source of evolutionary innovation and is associated with the increases in biological complexity and adaptive radiation. Termites are model social organisms characterized by a sophisticated caste system. We analyzed the genome of the Japanese subterranean termite, an ecologically and economically important insect acting as a destructive pest. The analyses revealed the significance of gene duplication in social evolution. Gene duplication associated with caste-biased gene expression was prevalent in the termite genome. Many of the duplicated genes were related to social functions, such as chemical communication, social immunity, and defense, and they were often expressed in caste-specific organs. We propose that gene duplication facilitates social evolution through regulatory diversification leading to caste-biased expression and functional specialization.

Termites are model social organisms characterized by a polyphenic caste system. Subterranean termites (Rhinotermitidae) are ecologically and economically important species, including acting as destructive pests. Rhinotermitidae occupies an important evolutionary position within the clade representing a transitional taxon between the higher (Termitidae) and lower (other families) termites. Here, we report the genome, transcriptome, and methylome of the Japanese subterranean termite Reticulitermes speratus. Our analyses highlight the significance of gene duplication in social evolution in this termite. Gene duplication associated with caste-biased gene expression was prevalent in the R. speratus genome. The duplicated genes comprised diverse categories related to social functions, including lipocalins (chemical communication), cellulases (wood digestion and social interaction), lysozymes (social immunity), geranylgeranyl diphosphate synthase (social defense), and a novel class of termite lineage–specific genes with unknown functions. Paralogous genes were often observed in tandem in the genome, but their expression patterns were highly variable, exhibiting caste biases. Some of the assayed duplicated genes were expressed in caste-specific organs, such as the accessory glands of the queen ovary and the frontal glands of soldier heads. We propose that gene duplication facilitates social evolution through regulatory diversification, leading to caste-biased expression and subfunctionalization and/or neofunctionalization conferring caste-specialized functions.

Functional Genomics and Comparative Lineage-Specific Region Analyses Reveal Novel Insights into Race Divergence in Verticillium dahliae

Verticillium dahliae is a widespread soilborne fungus that causes Verticillium wilt on numerous economically important plant species. In tomato, until now, three races have been characterized based on the response of differential cultivars to V. dahliae, but the genetic basis of race divergence in V. dahliae remains undetermined. To investigate the genetic basis of race divergence, we sequenced the genomes of two race 2 strains and four race 3 strains for comparative analyses with two known race 1 genomes. The genetic basis of race divergence was described by the pathogenicity-related genes among the three races, orthologue analyses, and genomic structural variations. Global comparative genomics showed that chromosomal rearrangements are not the only source of race divergence and that race 3 should be split into two genotypes based on orthologue clustering. Lineage-specific regions (LSRs), frequently observed between genomes of the three races, encode several predicted secreted proteins that potentially function as suppressors of immunity triggered by known effectors. These likely contribute to the virulence of the three races. Two genes in particular that can act as markers for race 2 and race 3 (VdR2e and VdR3e, respectively) contribute to virulence on tomato, and the latter acts as an avirulence factor of race 3. We elucidated the genetic basis of race divergence through global comparative genomics and identified secreted proteins in LSRs that could potentially play critical roles in the differential virulence among the races in V. dahliae.

IMPORTANCE Deciphering the gene-for-gene relationships during host-pathogen interactions is the basis of modern plant resistance breeding. In the Verticillium dahliae-tomato pathosystem, two races (races 1 and 2) and their corresponding avirulence (Avr) genes have been identified, but strains that lack these two Avr genes exist in nature. In this system, race 3 has been described, but the corresponding Avr gene has not been identified. We de novo-sequenced genomes of six strains and identified secreted proteins within the lineage-specific regions (LSRs) distributed among the genomes of the three races that could potentially function as manipulators of host immunity. One of the LSR genes, VdR3e, was confirmed as the Avr gene for race 3. The results indicate that differences in transcriptional regulation may contribute to race differentiation. This is the first study to describe these differences and elucidate roles of secreted proteins in LSRs that play roles in race differentiation.

High-quality genome sequence assembly of R.A73 Enterococcus faecium isolated from freshwater fish mucus

BACKGROUND

Whole-genome sequencing using high throughput technologies has revolutionized and speeded up the scientific investigation of bacterial genetics, biochemistry, and molecular biology. Lactic acid bacteria (LABs) have been extensively used in fermentation and more recently as probiotics in food products that promote health. Genome sequencing and functional genomics investigations of LABs varieties provide rapid and important information about their diversity and their evolution, revealing a significant molecular basis. This study investigated the whole genome sequences of the Enterococcus faecium strain (HG937697), isolated from the mucus of freshwater fish in Tunisian dams. Genomic DNA was extracted using the Quick-GDNA kit and sequenced using the Illumina HiSeq2500 system. Sequences quality assessment was performed using FastQC software. The complete genome annotation was carried out with the Rapid Annotation using Subsystem Technology (RAST) web server then NCBI PGAAP.

RESULTS

The Enterococcus faecium R.A73 assembled in 28 contigs consisting of 2,935,283 bps. The genome annotation revealed 2884 genes in total including 2834 coding sequences and 50 RNAs containing 3 rRNAs (one rRNA 16 s, one rRNA 23 s and one rRNA 5 s) and 47 tRNAs. Twenty-two genes implicated in bacteriocin production are identified within the Enterococcus faecium R.A73 strain.

CONCLUSION

Data obtained provide insights to further investigate the effective strategy for testing this Enterococcus faecium R.A73 strain in the industrial manufacturing process. Studying their metabolism with bioinformatics tools represents the future challenge and contribution to improving the utilization of the multi-purpose bacteria in food.

Conserved Patterns in Developmental Processes and Phases, Rather than Genes, Unite the Highly Divergent Bilateria

Bilateria are the predominant clade of animals on Earth. Despite having evolved a wide variety of body plans and developmental modes, they are characterized by common morphological traits. By default, researchers have tried to link clade-specific genes to these traits, thus distinguishing bilaterians from non-bilaterians, by their gene content. Here we argue that it is rather biological processes that unite Bilateria and set them apart from their non-bilaterian sisters, with a less complex body morphology. To test this hypothesis, we compared proteomes of bilaterian and non-bilaterian species in an elaborate computational pipeline, aiming to search for a set of bilaterian-specific genes. Despite the limited confidence in their bilaterian specificity, we nevertheless detected Bilateria-specific functional and developmental patterns in the sub-set of genes conserved in distantly related Bilateria. Using a novel multi-species GO-enrichment method, we determined the functional repertoire of genes that are widely conserved among Bilateria. Analyzing expression profiles in three very distantly related model species-D. melanogaster, D. rerio and C. elegans-we find characteristic peaks at comparable stages of development and a delayed onset of expression in embryos. In particular, the expression of the conserved genes appears to peak at the phylotypic stage of different bilaterian phyla. In summary, our study illustrate how development connects distantly related Bilateria after millions of years of divergence, pointing to processes potentially separating them from non-bilaterians. We argue that evolutionary biologists should return from a purely gene-centric view of evolution and place more focus on analyzing and defining conserved developmental processes and periods.

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.

Identification and evolutionary analysis of the nucleolar proteome of Giardia lamblia

Background

The nucleoli, including their proteomes, of higher eukaryotes have been extensively studied, while few studies about the nucleoli of the lower eukaryotes – protists were reported. Giardia lamblia, a protist with the controversy of whether it is an extreme primitive eukaryote or just a highly evolved parasite, might be an interesting object for carrying out the nucleolar proteome study of protists and for further examining the controversy.

Results

Using bioinformatics methods, we reconstructed G. lamblia nucleolar proteome (GiNuP) and the common nucleolar proteome of the three representative higher eukaryotes (human, Arabidopsis, yeast) (HEBNuP). Comparisons of the two proteomes revealed that: 1) GiNuP is much smaller than HEBNuP, but 78.4% of its proteins have orthologs in the latter; 2) More than 68% of the GiNuP proteins are involved in the “Ribosome related” function, and the others participate in the other functions, and these two groups of proteins are much larger and much smaller than those in HEBNuP, respectively; 3) Both GiNuP and HEBNuP have their own specific proteins, but HEBNuP has a much higher proportion of such proteins to participate in more categories of nucleolar functions.

Conclusion

For the first time the nucleolar proteome of a protist - Giardia was reconstructed. The results of comparison of it with the common proteome of three representative higher eukaryotes -- HEBNuP indicated that the simplicity of GiNuP is most probably a reflection of primitiveness but not just parasitic reduction of Giardia, and simultaneously revealed some interesting evolutionary phenomena about the nucleolus and even the eukaryotic cell, compositionally and functionally.

Phylogeny and highland adaptation of Chinese species in Allium section Daghestanica (Amaryllidaceae) revealed by transcriptome sequencing

Allium L. is one of the largest monocotyledonous genera with extensive distribution in the Northern Hemisphere. The fundamental phylogenies of Allium have been investigated using many morphological and molecular characters. However, the morphological characters may not agree with the molecular results in some Allium groups or sections (such as the Chinese Allium section Daghestanica), which may result in ambiguous species relationships and hinder further evolutionary and adaptive researches. Here, transcriptome sequences of the six Chinese endemics from Allium section Daghestanica were collected, with their single-copy genes (SCGs) were extracted. The interspecies relationships were analyzed using concatenation and coalescent methods. The branch-site model (BSM) was conducted to detect the positively selected genes (PSGs) in five highland species of this section. Based on 1644, 1281 and 1580 SCGs in flowers, leaves, and flowers-leaves combination respectively, a robust consistent and well-resolved phylogeny was generated from the concatenation method. Strong conflicts among individual gene trees were detected in the coalescent method, and morphological characters were incongruent with molecular relationships to some degree. Many PSGs were involved in responses of various stresses and stimuli (e.g. hypoxia, low temperature, aridity), DNA repair, metabolism, nutrient or energy intake, photosynthesis, and signal transduction. Our study revealed a clear interspecies relationship of Chinese endemics in Allium section Daghestanica and suggested that the discordance between morphological characters and molecular relationships might result from that the former are more susceptible to convergence compared with the latter. PSGs detected in our study may provide some insights into highland adaptation in Allium species.

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.

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.

Deciphering the biological processes underlying tomato biomass production and composition

The huge amounts of biomass residues, remaining in the field after tomato fruits harvesting, can be utilized to produce bioenergy. A multiple level approach aimed to characterize two Solanum pennellii introgression lines (ILs), with contrasting phenotypes for plant architecture and biomass was carried out. The study of gene expression dynamics, microscopy cell traits and qualitative and quantitative cell wall chemical compounds variation enabled the discovery of key genes and cell processes involved biomass accumulation and composition. Enhanced biomass production observed in IL2-6 line is due to a more effective coordination of chloroplasts and mitochondria energy fluxes. Microscopy analysis revealed a higher number of cells and chloroplasts in leaf epidermis in the high biomass line whilst chemical measurements on the two lines pointed out striking differences in the cell wall composition and organization. Taken together, our findings shed light on the mechanisms underlying the tomato biomass production and processability.

Genome and transcriptome evolve separately in recently hybridized Trichosporon fungi

Genome hybridization is an important evolutionary event that gives rise to species with novel capabilities. However, the merging of distinct genomes also brings together incompatible regulatory networks that must be resolved during the course of evolution. Understanding of the early stages of post-hybridization evolution is particularly important because changes in these stages have long-term evolutionary consequences. Here, via comparative transcriptomic analyses of two closely related, recently hybridized Trichosporon fungi, T. coremiiforme and T. ovoides, and three extant relatives, we show that early post-hybridization evolutionary processes occur separately at the gene sequence and gene expression levels but together contribute to the stabilization of hybrid genome and transcriptome. Our findings also highlight lineage-specific consequences of genome hybridization, revealing that the transcriptional regulatory dynamics in these hybrids responded completely differently to gene loss events: one involving both subgenomes and another that is strictly subgenome-specific.

Evolutionary Rates of Bumblebee Genomes Are Faster at Lower Elevations

The importance of climate in determining biodiversity patterns has been well documented. However, the relationship between climate and rates of genetic evolution remains controversial. Latitude and elevation have been associated with rates of change in genetic markers such as cytochrome b. What is not known, however, is the strength of such associations and whether patterns found among these genes apply across entire genomes. Here, using bumblebee genetic data from seven subgenera of Bombus, we demonstrate that all species occupying warmer elevations have undergone faster genome-wide evolution than those in the same subgenera occupying cooler elevations. Our findings point to a critical biogeographic role in the relative rates of whole species evolution, potentially influencing global biodiversity patterns.

Pangloss: A Tool for Pan-Genome Analysis of Microbial Eukaryotes

Although the pan-genome concept originated in prokaryote genomics, an increasing number of eukaryote species pan-genomes have also been analysed. However, there is a relative lack of software intended for eukaryote pan-genome analysis compared to that available for prokaryotes. In a previous study, we analysed the pan-genomes of four model fungi with a computational pipeline that constructed pan-genomes using the synteny-dependent Pan-genome Ortholog Clustering Tool (PanOCT) approach. Here, we present a modified and improved version of that pipeline which we have called Pangloss. Pangloss can perform gene prediction for a set of genomes from a given species that the user provides, constructs and optionally refines a species pan-genome from that set using PanOCT, and can perform various functional characterisation and visualisation analyses of species pan-genome data. To demonstrate Pangloss's capabilities, we constructed and analysed a species pan-genome for the oleaginous yeast Yarrowialipolytica and also reconstructed a previously-published species pan-genome for the opportunistic respiratory pathogen Aspergillus fumigatus. Pangloss is implemented in Python, Perl and R and is freely available under an open source GPLv3 licence via GitHub.

Evolution in chronic cold: varied loss of cellular response to heat in Antarctic notothenioid fish

Background

Confined within the freezing Southern Ocean, the Antarctic notothenioids have evolved to become both cold adapted and cold specialized. A marked signature of cold specialization is an apparent loss of the cellular heat shock response (HSR). As the HSR has been examined in very few notothenioid species to-date, it remains unknown whether HSR loss pervades the Antarctic radiation, or whether the broader cellular responses to heat stress has sustained similar loss. Understanding the evolutionary status of these responses in this stenothermal taxon is crucial for evaluating its adaptive potential to ocean warming under climate change.

Results

In this study, we used an acute heat stress protocol followed by RNA-Seq analyses to study the evolution of cellular-wide transcriptional responses to heat stress across three select notothenioid lineages - the basal temperate and nearest non-Antarctic sister species Eleginops maclovinus serving as ancestral proxy, the cryopelagic Pagothenia borchgrevinki and the icefish Chionodraco rastrospinosus representing cold-adapted red-blooded and hemoglobinless Antarctic notothenioids respectively. E. maclovinus displayed robust cellular stress responses including the ER Unfolded Protein Response and the cytosolic HSR, cementing the HSR as a plesiomorphy that preceded Antarctic notothenioid radiation. While the transcriptional response to heat stress was minimal in P. borchgrevinki, C. rastrospinosus exhibited robust responses in the broader cellular networks especially in inflammatory responses despite lacking the classic HSR and UPR.

Conclusion

The disparate patterns observed in these two archetypal Antarctic species indicate the evolutionary status in cellular ability to mitigate acute heat stress varies even among Antarctic lineages, which may affect their adaptive potential in coping with a warming world.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1254-6) contains supplementary material, which is available to authorized users.

PGAweb: A Web Server for Bacterial Pan-Genome Analysis

An astronomical increase in microbial genome data in recent years has led to strong demand for bioinformatic tools for pan-genome analysis within and across species. Here, we present PGAweb, a user-friendly, web-based tool for bacterial pan-genome analysis, which is composed of two main pan-genome analysis modules, PGAP and PGAP-X. PGAweb provides key interactive and customizable functions that include orthologous clustering, pan-genome profiling, sequence variation and evolution analysis, and functional classification. PGAweb presents features of genomic structural dynamics and sequence diversity with different visualization methods that are helpful for intuitively understanding the dynamics and evolution of bacterial genomes. PGAweb has an intuitive interface with one-click setting of parameters and is freely available at http://PGAweb.vlcc.cn/.

De novo transcriptome assembly and positive selection analysis of an individual deep-sea fish

Background

High hydrostatic pressure and low temperatures make the deep sea a harsh environment for life forms. Actin organization and microtubules assembly, which are essential for intracellular transport and cell motility, can be disrupted by high hydrostatic pressure. High hydrostatic pressure can also damage DNA. Nucleic acids exposed to low temperatures can form secondary structures that hinder genetic information processing. To study how deep-sea creatures adapt to such a hostile environment, one of the most straightforward ways is to sequence and compare their genes with those of their shallow-water relatives.

Results

We captured an individual of the fish species Aldrovandia affinis, which is a typical deep-sea inhabitant, from the Okinawa Trough at a depth of 1550 m using a remotely operated vehicle (ROV). We sequenced its transcriptome and analyzed its molecular adaptation. We obtained 27,633 protein coding sequences using an Illumina platform and compared them with those of several shallow-water fish species. Analysis of 4918 single-copy orthologs identified 138 positively selected genes in A. affinis, including genes involved in microtubule regulation. Particularly, functional domains related to cold shock as well as DNA repair are exposed to positive selection pressure in both deep-sea fish and hadal amphipod.

Conclusions

Overall, we have identified a set of positively selected genes related to cytoskeleton structures, DNA repair and genetic information processing, which shed light on molecular adaptation to the deep sea. These results suggest that amino acid substitutions of these positively selected genes may contribute crucially to the adaptation of deep-sea animals. Additionally, we provide a high-quality transcriptome of a deep-sea fish for future deep-sea studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4720-z) contains supplementary material, which is available to authorized users.

Surveying alignment-free features for Ortholog detection in related yeast proteomes by using supervised big data classifiers

BACKGROUND

The development of new ortholog detection algorithms and the improvement of existing ones are of major importance in functional genomics. We have previously introduced a successful supervised pairwise ortholog classification approach implemented in a big data platform that considered several pairwise protein features and the low ortholog pair ratios found between two annotated proteomes (Galpert, D et al., BioMed Research International, 2015). The supervised models were built and tested using a Saccharomycete yeast benchmark dataset proposed by Salichos and Rokas (2011). Despite several pairwise protein features being combined in a supervised big data approach; they all, to some extent were alignment-based features and the proposed algorithms were evaluated on a unique test set. Here, we aim to evaluate the impact of alignment-free features on the performance of supervised models implemented in the Spark big data platform for pairwise ortholog detection in several related yeast proteomes.

RESULTS

The Spark Random Forest and Decision Trees with oversampling and undersampling techniques, and built with only alignment-based similarity measures or combined with several alignment-free pairwise protein features showed the highest classification performance for ortholog detection in three yeast proteome pairs. Although such supervised approaches outperformed traditional methods, there were no significant differences between the exclusive use of alignment-based similarity measures and their combination with alignment-free features, even within the twilight zone of the studied proteomes. Just when alignment-based and alignment-free features were combined in Spark Decision Trees with imbalance management, a higher success rate (98.71%) within the twilight zone could be achieved for a yeast proteome pair that underwent a whole genome duplication. The feature selection study showed that alignment-based features were top-ranked for the best classifiers while the runners-up were alignment-free features related to amino acid composition.

CONCLUSIONS

The incorporation of alignment-free features in supervised big data models did not significantly improve ortholog detection in yeast proteomes regarding the classification qualities achieved with just alignment-based similarity measures. However, the similarity of their classification performance to that of traditional ortholog detection methods encourages the evaluation of other alignment-free protein pair descriptors in future research.

Databases and Computational Tools for Evolutionary Analysis of Protein Phosphorylation

Advancements in MS-based phospho-proteomics techniques have helped uncover hundred thousands of protein phosphorylation sites in human and various model organisms. The majority of these sites are uncharacterized. The sheer number of uncharacterized sites necessitates systematic approaches to prioritize sites for more in-depth annotation. Analyzing the phosphorylation and sequence conservation of uncharacterized sites across species can help reveal a subset of the functionally important phosphorylation events. Here, we outline the workflow and provide an overview of publicly available computational resources for conservation analysis of novel phosphorylation sites.

A multi-omics study of the grapevine-downy mildew (Plasmopara viticola) pathosystem unveils a complex protein coding- and noncoding-based arms race during infection

Fungicides are applied intensively to prevent downy mildew infections of grapevines (Vitis vinifera) with high impact on the environment. In order to develop alternative strategies we sequenced the genome of the oomycete pathogen Plasmopara viticola causing this disease. We show that it derives from a Phytophthora-like ancestor that switched to obligate biotrophy by losing genes involved in nitrogen metabolism and γ-Aminobutyric acid catabolism. By combining multiple omics approaches we characterized the pathosystem and identified a RxLR effector that trigger an immune response in the wild species V. riparia. This effector is an ideal marker to screen novel grape resistant varieties. Our study reveals an unprecedented bidirectional noncoding RNA-based mechanism that, in one direction might be fundamental for P. viticola to proficiently infect its host, and in the other might reduce the effects of the infection on the plant.

Comparative genomics reveals cotton-specific virulence factors in flexible genomic regions in Verticillium dahliae and evidence of horizontal gene transfer from Fusarium

Verticillium dahliae isolates are most virulent on the host from which they were originally isolated. Mechanisms underlying these dominant host adaptations are currently unknown. We sequenced the genome of V. dahliae Vd991, which is highly virulent on its original host, cotton, and performed comparisons with the reference genomes of JR2 (from tomato) and VdLs.17 (from lettuce). Pathogenicity-related factor prediction, orthology and multigene family classification, transcriptome analyses, phylogenetic analyses, and pathogenicity experiments were performed. The Vd991 genome harbored several exclusive, lineage-specific (LS) genes within LS regions (LSRs). Deletion mutants of the seven genes within one LSR (G-LSR2) in Vd991 were less virulent only on cotton. Integration of G-LSR2 genes individually into JR2 and VdLs.17 resulted in significantly enhanced virulence on cotton but did not affect virulence on tomato or lettuce. Transcription levels of the seven LS genes in Vd991 were higher during the early stages of cotton infection, as compared with other hosts. Phylogenetic analyses suggested that G-LSR2 was acquired from Fusarium oxysporum f. sp. vasinfectum through horizontal gene transfer. Our results provide evidence that horizontal gene transfer from Fusarium to Vd991 contributed significantly to its adaptation to cotton and may represent a significant mechanism in the evolution of an asexual plant pathogen.

MM, Prakash Mishra G, Singh D, Ye W, Tyler BM, Tripathy S. EumicrobeDBLite: a lightweight genomic resource and analytic platform for draft oomycete genomes

We have developed EumicrobeDBLite-a lightweight comprehensive genome resource and sequence analysis platform for oomycete organisms. EumicrobeDBLite is a successor of the VBI Microbial Database (VMD) that was built using the Genome Unified Schema (GUS). In this version, GUS has been greatly simplified with the removal of many obsolete modules and the redesign of others to incorporate contemporary data. Several dependences, such as perl object layers used for data loading in VMD, have been replaced with independent lightweight scripts. EumicrobeDBLite now runs on a powerful annotation engine developed at our laboratory, called 'Genome Annotator Lite'. Currently, this database has 26 publicly available genomes and 10 expressed sequence tag (EST) datasets of oomycete organisms. The browser page has dynamic tracks presenting comparative genomics analyses, coding and non-coding data, tRNA genes, repeats and EST alignments. In addition, we have defined 44 777 core conserved proteins from 12 oomycete organisms which form 2974 clusters. Synteny viewing is enabled by the incorporation of the Genome Synteny Viewer (GSV) tool. The user interface has undergone major changes for ease of browsing. Queryable comparative genomics information, conserved orthologous genes and pathways are among the new key features updated in this database. The browser has been upgraded to enable user upload of GFF files for quick view of genome annotation comparisons. The toolkit page integrates the EMBOSS package and has a gene prediction tool. Annotations for the organisms are updated once every 6 months to ensure quality. The database resource is available at www.eumicrobedb.org.

Pan-Genome Storage and Analysis Techniques

Computational pan-genome analysis has emerged from the rapid increase of available genome sequencing data. Starting from a microbial pan-genome, the concept has spread to a variety of species, such as plants or viruses. Characterizing a pan-genome provides insights into intra-species evolution, functions, and diversity. However, researchers face challenges such as processing and maintaining large datasets while providing accurate and efficient analysis approaches. Comparative genomics methods are required for detecting conserved and unique regions between a set of genomes. This chapter gives an overview of tools available for indexing pan-genomes, identifying the sub-regions of a pan-genome and offering a variety of downstream analysis methods. These tools are categorized into two groups, gene-based and sequence-based, according to the pan-genome identification method. We highlight the differences, advantages, and disadvantages between the tools, and provide information about the general workflow, methodology of pan-genome identification, covered functionalities, usability and availability of the tools.

A Trichosporonales genome tree based on 27 haploid and three evolutionarily conserved 'natural' hybrid genomes

To construct a backbone tree consisting of basidiomycetous yeasts, draft genome sequences from 25 species of Trichosporonales (Tremellomycetes, Basidiomycota) were generated. In addition to the hybrid genomes of Trichosporon coremiiforme and Trichosporon ovoides that we described previously, we identified an interspecies hybrid genome in Cutaneotrichosporon mucoides (formerly Trichosporon mucoides). This hybrid genome had a gene retention rate of ~55%, and its closest haploid relative was Cutaneotrichosporon dermatis. After constructing the C. mucoides subgenomes, we generated a phylogenetic tree using genome data from the 27 haploid species and the subgenome data from the three hybrid genome species. It was a high-quality tree with 100% bootstrap support for all of the branches. The genome-based tree provided superior resolution compared with previous multi-gene analyses. Although our backbone tree does not include all Trichosporonales genera (e.g. Cryptotrichosporon), it will be valuable for future analyses of genome data. Interest in interspecies hybrid fungal genomes has recently increased because they may provide a basis for new technologies. The three Trichosporonales hybrid genomes described in this study are different from well-characterized hybrid genomes (e.g. those of Saccharomyces pastorianus and Saccharomyces bayanus) because these hybridization events probably occurred in the distant evolutionary past. Hence, they will be useful for studying genome stability following hybridization and speciation events. Copyright © 2017 John Wiley & Sons, Ltd.

Evolutionary recruitment of flexible Esrp-dependent splicing programs into diverse embryonic morphogenetic processes

Epithelial-mesenchymal interactions are crucial for the development of numerous animal structures. Thus, unraveling how molecular tools are recruited in different lineages to control interplays between these tissues is key to understanding morphogenetic evolution. Here, we study Esrp genes, which regulate extensive splicing programs and are essential for mammalian organogenesis. We find that Esrp homologs have been independently recruited for the development of multiple structures across deuterostomes. Although Esrp is involved in a wide variety of ontogenetic processes, our results suggest ancient roles in non-neural ectoderm and regulating specific mesenchymal-to-epithelial transitions in deuterostome ancestors. However, consistent with the extensive rewiring of Esrp-dependent splicing programs between phyla, most developmental defects observed in vertebrate mutants are related to other types of morphogenetic processes. This is likely connected to the origin of an event in Fgfr, which was recruited as an Esrp target in stem chordates and subsequently co-opted into the development of many novel traits in vertebrates.

Popmarker: Identifying Phylogenetic Markers at the Population Level

As phylogenomic approach becomes a common practice for constructing true bacterial phylogenies, it has become apparent that single molecular markers such as 16S ribosomal DNA often lead to misclassification of species. In this study, we present a program called Popmarker that uses the true species phylogeny and identifies a minimum set of molecular markers reflecting the bacterial evolution history and phylogenetic relationship at the resolution of populations. Popmarker ranks the proteome according to the correlation of whole species tree or subtree branch length against orthologous sequence distances. We demonstrate that 5 proteins of 2 top ranks achieve the same resolution as concatenation of 2203 single-copy orthologous genes and the right species classification as well as correct split of the 2 groups of Vibrio campbellii . The top-ranking genes selected by Popmarker are candidates that lead to speciation and are useful in distinguishing close related species in microbiome study.

Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination

Background

Root and tuber crops are a major food source in tropical Africa. Among these crops are several species in the monocotyledonous genus Dioscorea collectively known as yam, a staple tuber crop that contributes enormously to the subsistence and socio-cultural lives of millions of people, principally in West and Central Africa. Yam cultivation is constrained by several factors, and yam can be considered a neglected “orphan” crop that would benefit from crop improvement efforts. However, the lack of genetic and genomic tools has impeded the improvement of this staple crop.

Results

To accelerate marker-assisted breeding of yam, we performed genome analysis of white Guinea yam (Dioscorea rotundata) and assembled a 594-Mb genome, 76.4% of which was distributed among 21 linkage groups. In total, we predicted 26,198 genes. Phylogenetic analyses with 2381 conserved genes revealed that Dioscorea is a unique lineage of monocotyledons distinct from the Poales (rice), Arecales (palm), and Zingiberales (banana). The entire Dioscorea genus is characterized by the occurrence of separate male and female plants (dioecy), a feature that has limited efficient yam breeding. To infer the genetics of sex determination, we performed whole-genome resequencing of bulked segregants (quantitative trait locus sequencing [QTL-seq]) in F1 progeny segregating for male and female plants and identified a genomic region associated with female heterogametic (male = ZZ, female = ZW) sex determination. We further delineated the W locus and used it to develop a molecular marker for sex identification of Guinea yam plants at the seedling stage.

Conclusions

Guinea yam belongs to a unique and highly differentiated clade of monocotyledons. The genome analyses and sex-linked marker development performed in this study should greatly accelerate marker-assisted breeding of Guinea yam. In addition, our QTL-seq approach can be utilized in genetic studies of other outcrossing crops and organisms with highly heterozygous genomes. Genomic analysis of orphan crops such as yam promotes efforts to improve food security and the sustainability of tropical agriculture.

Electronic supplementary material

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

Transcriptomes and expression profiling of deep-sea corals from the Red Sea provide insight into the biology of azooxanthellate corals

Despite the importance of deep-sea corals, our current understanding of their ecology and evolution is limited due to difficulties in sampling and studying deep-sea environments. Moreover, a recent re-evaluation of habitat limitations has been suggested after characterization of deep-sea corals in the Red Sea, where they live at temperatures of above 20 °C at low oxygen concentrations. To gain further insight into the biology of deep-sea corals, we produced reference transcriptomes and studied gene expression of three deep-sea coral species from the Red Sea, i.e. Dendrophyllia sp., Eguchipsammia fistula, and Rhizotrochus typus. Our analyses suggest that deep-sea coral employ mitochondrial hypometabolism and anaerobic glycolysis to manage low oxygen conditions present in the Red Sea. Notably, we found expression of genes related to surface cilia motion that presumably enhance small particle transport rates in the oligotrophic deep-sea environment. This is the first study to characterize transcriptomes and in situ gene expression for deep-sea corals. Our work offers several mechanisms by which deep-sea corals might cope with the distinct environmental conditions present in the Red Sea As such, our data provide direction for future research and further insight to organismal response of deep-sea coral to environmental change and ocean warming.

Investigation of orthologous pathogen recognition gene-rich regions in solanaceous species

Pathogen receptor proteins such as receptor-like protein (RLP), receptor-like kinase (RLK), and nucleotide-binding leucine-rich repeat (NLR) play a leading role in plant immunity activation. The genome architecture of such genes has been extensively investigated in several plant species. However, we still know little about their elaborate reorganization that arose during the plant speciation process. Using recently released pepper and eggplant genome sequences, we were able to identify 1097 pathogen recognition genes (PRGs) in the cultivated pepper Zunla-1 and 775 in the eggplant line Nakate-Shinkuro. The retrieved genes were analysed for their tendency to cluster, using different methods to infer the means of grouping. Orthologous relationships among clustering loci were found, and interesting reshuffling within given loci was observed for each analysed species. The information obtained was integrated into a comparative map to highlight the evolutionary dynamics in which the PRG loci were involved. Diversification of 14 selected PRG-rich regions was also explored using a DNA target-enrichment approach. A large number of gene variants were found as well as rearrangements of sequences encoding single protein domain and changes in chromosome gene order among species. Gene duplication and transposition activity have clearly influenced plant genome R-gene architecture and diversification. Our findings contribute to addressing several biological questions concerning the parallel evolution that occurred between genomes of the family Solanaceae. Moreover, the integration of different methods proved a powerful approach to reconstruct the evolutionary history in plant families and to transfer important biology findings among plant genomes.

Comparative analysis of the genomes of clinical isolates of Mycobacterium avium subsp. hominissuis regarding virulence-related genes

PURPOSE

Mycobacterium avium subsp. hominissuis is a member of the M. avium complex, a heterogeneous group of bacteria that cause lung infection in immunocompetent patients or disseminated infection in patients with immunosuppression. The bacteria belonging to this complex have variable virulence, depending on the strain considered, and therefore a representative of the most common clinical phenotype was analysed.

METHODOLOGY

The genomic sequences of four M. avium subsp. hominissuis isolates obtained from clinical specimens were completed. Mav101, Mav100 and MavA5 were isolated from the blood of patients with AIDS. MavA5 was disseminated from the lung, while Mav3388 was isolated from the lungs of a patient with chronic lung disease. The sequences were annotated using the published Mav104 genome as a blueprint. Functional and virulence analyses of the sequences were carried out. Mice studies comparing the virulence of the strains were performed.

RESULTS

Findings showed that while Mav101 was very similar to Mav104, there were numerous differences between Mav104 and the remaining strains at nucleotide and predicted protein levels. The presence of genes associated with biofilm formation and several known virulence-related genes were sometimes differentially present among the isolates, suggesting overlapping functions by different genetic determinants.

CONCLUSIONS

The sequences provided important information about M. avium heterogenicity and evolution as a pathogen. The limitation is the lack of understanding on possible overlapping functions of genes/proteins.

Transcriptome analysis explores genes related to shikonin biosynthesis in Lithospermeae plants and provides insights into Boraginales' evolutionary history

Shikonin and its derivatives extracted from Lithospermeae plants' red roots have current applications in food and pharmaceutical industries. Previous studies have cloned some genes related to shikonin biosynthesis. However, most genes related to shikonin biosynthesis remain unclear, because the lack of the genome/transcriptome of the Lithospermeae plants. Therefore, in order to provide a new understanding of shikonin biosynthesis, we obtained transcriptome data and unigenes expression profiles in three shikonin-producing Lithospermeae plants, i.e., Lithospermum erythrorhizon, Arnebia euchroma and Echium plantagineum. As a result, two unigenes (i.e., G10H and 12OPR) that are involved in "shikonin downstream biosynthesis" and "methyl jasmonate biosynthesis" were deemed to relate to shikonin biosynthesis in this study. Furthermore, we conducted a Lamiids phylogenetic model and identified orthologous unigenes under positive selection in above three Lithospermeae plants. The results indicated Boraginales was more relative to Solanales/Gentianales than to Lamiales.

PanWeb: A web interface for pan-genomic analysis

With increased production of genomic data since the advent of next-generation sequencing (NGS), there has been a need to develop new bioinformatics tools and areas, such as comparative genomics. In comparative genomics, the genetic material of an organism is directly compared to that of another organism to better understand biological species. Moreover, the exponentially growing number of deposited prokaryote genomes has enabled the investigation of several genomic characteristics that are intrinsic to certain species. Thus, a new approach to comparative genomics, termed pan-genomics, was developed. In pan-genomics, various organisms of the same species or genus are compared. Currently, there are many tools that can perform pan-genomic analyses, such as PGAP (Pan-Genome Analysis Pipeline), Panseq (Pan-Genome Sequence Analysis Program) and PGAT (Prokaryotic Genome Analysis Tool). Among these software tools, PGAP was developed in the Perl scripting language and its reliance on UNIX platform terminals and its requirement for an extensive parameterized command line can become a problem for users without previous computational knowledge. Thus, the aim of this study was to develop a web application, known as PanWeb, that serves as a graphical interface for PGAP. In addition, using the output files of the PGAP pipeline, the application generates graphics using custom-developed scripts in the R programming language. PanWeb is freely available at http://www.computationalbiology.ufpa.br/panweb.

Molecular adaptation in the world's deepest-living animal: Insights from transcriptome sequencing of the hadal amphipod Hirondellea gigas

The Challenger Deep in the Mariana Trench is the deepest point in the oceans of our planet. Understanding how animals adapt to this harsh environment characterized by high hydrostatic pressure, food limitation, dark and cold is of great scientific interest. Of the animals dwelling in the Challenger Deep, amphipods have been captured using baited traps. In this study, we sequenced the transcriptome of the amphipod Hirondellea gigas collected at a depth of 10,929 m from the East Pond of the Challenger Deep. Assembly of these sequences resulted in 133,041 contigs and 22,046 translated proteins. Functional annotation of these contigs was made using the go and kegg databases. Comparison of these translated proteins with those of four shallow-water amphipods revealed 10,731 gene families, of which 5659 were single-copy orthologs. Base substitution analysis on these single-copy orthologs showed that 62 genes are positively selected in H. gigas, including genes related to β-alanine biosynthesis, energy metabolism and genetic information processing. For multiple-copy orthologous genes, gene family expansion analysis revealed that cold-inducible proteins (i.e., transcription factors II A and transcription elongation factor 1) as well as zinc finger domains are expanded in H. gigas. Overall, our results indicate that genetic adaptation to the hadal environment by H. gigas may be mediated by both gene family expansion and amino acid substitutions of specific proteins.

Genome-wide identification of genes probably relevant to the adaptation of schizothoracins (Teleostei: Cypriniformes) to the uplift of the Qinghai-Tibet Plateau

BACKGROUND

Molecular adaptation to the severe environments present during the uplift of the Qinghai-Tibet Plateau has attracted the attention of researchers. The divergence of the three specialization groups of schizothoracins (Primitive, Specialized and Highly Specialized) may correspond to the three phases of plateau uplift. Based on the transcripts of representative species of the three specialized groups and an outgroup, genes in schizothoracins that may have played important roles during the adaptation to new environments were investigated.

RESULTS

The contigs of Gymnodiptychus dybowskii and Schizothorax pseudaksaiensis were compared with those of Gymnocypris przewalskii ganzihonensis and the outgroup Sinocyclocheilus angustiporus, and 5,894 ortholog groups with an alignment length longer than 90 nt after deleting gaps were retained. Evolutionary analyses indicated that the average evolutionary rate of the branch leading to the Specialized group was faster than that of the branch leading to the Highly Specialized group. Moreover, the numbers of gene categories in which more than half of the genes evolved faster than the average values of the genome were 117 and 15 along the branches leading to the Specialized and Highly Specialized groups, respectively. A total of 40, 36, and 55 genes were likely subject to positive selection along the branches leading to the Primitive, Specialized and Highly Specialized groups, respectively, and many of these genes are likely relevant to adaptation to the cold temperatures, low oxygen concentrations, and strong ultraviolet radiation that result from elevation.

CONCLUSIONS

By selecting representative species of the three groups of schizothoracins and applying next-generation sequencing technology, several candidate genes corresponding to adaptation to the three phases of plateau uplift were identified. Some of the genes identified in this report that were likely subject to positive selection are good candidates for subsequent evolutionary and functional analyses of adaptation to high altitude.

Genomes of Fasciola hepatica from the Americas Reveal Colonization with Neorickettsia Endobacteria Related to the Agents of Potomac Horse and Human Sennetsu Fevers

Food borne trematodes (FBTs) are an assemblage of platyhelminth parasites transmitted through the food chain, four of which are recognized as neglected tropical diseases (NTDs). Fascioliasis stands out among the other NTDs due to its broad and significant impact on both human and animal health, as Fasciola sp., are also considered major pathogens of domesticated ruminants. Here we present a reference genome sequence of the common liver fluke, Fasciola hepatica isolated from sheep, complementing previously reported isolate from cattle. A total of 14,642 genes were predicted from the 1.14 GB genome of the liver fluke. Comparative genomics indicated that F. hepatica Oregon and related food-borne trematodes are metabolically less constrained than schistosomes and cestodes, taking advantage of the richer millieux offered by the hepatobiliary organs. Protease families differentially expanded between diverse trematodes may facilitate migration and survival within the heterogeneous environments and niches within the mammalian host. Surprisingly, the sequencing of Oregon and Uruguay F. hepatica isolates led to the first discovery of an endobacteria in this species. Two contigs from the F. hepatica Oregon assembly were joined to complete the 859,205 bp genome of a novel Neorickettsia endobacterium (nFh) closely related to the etiological agents of human Sennetsu and Potomac horse fevers. Immunohistochemical studies targeting a Neorickettsia surface protein found nFh in specific organs and tissues of the adult trematode including the female reproductive tract, eggs, the Mehlis' gland, seminal vesicle, and oral suckers, suggesting putative routes for fluke-to-fluke and fluke-to-host transmission. The genomes of F. hepatica and nFh will serve as a resource for further exploration of the biology of F. hepatica, and specifically its newly discovered trans-kingdom interaction with nFh and the impact of both species on disease in ruminants and humans.

Whole-genome Sequencing for Surveillance of Invasive Pneumococcal Diseases in Ontario, Canada: Rapid Prediction of Genotype, Antibiotic Resistance and Characterization of Emerging Serotype 22F

Background: Molecular typing is essential for inferring genetic relatedness between bacterial pathogens. In this study, we applied whole genome sequencing (WGS) for rapid prediction of sequence type and antibiotic resistance for invasive pneumococcal isolates. Methods: 240 isolates from adults (≥50 years old) in Ontario, Canada during 2009 to 2013 were subjected to WGS. Sequence type, antibiotic susceptibility and resistance were predicted directly from short reads. Emerging non-vaccine serotype 22F was further characterized by WGS. Results: Sequence type was successfully determined for 98.3% of isolates. The overall sensitivity and specificity for antibiotic resistance prediction were 95 and 100% respectively, compared to standard susceptibility testing methods. WGS-based phylogeny divided emerging 22F (ST433) strains into two distinct clades: clade A harboring a 23 kb-prophage and anti-phage PhD/Doc system and clade B with virulence-related proteases. Five isolates in clade A developed macrolide resistance via 5.1 kb mega element recombination (encoding mefE and msrD), while one isolate in clade B displayed quinolone resistance via a gyrA mutation. Conclusions: WGS is valuable for routine surveillance of pneumococcal clinical isolates and facilitates prediction of genotype and antibiotic resistance. The emergence of 22F in Ontario in the post-vaccine era and evidence of evolution and divergence of the 22F population warrants heightened pneumococcal molecular surveillance.

Transcriptomic comparison of invasive bigheaded carps (Hypophthalmichthys nobilis and Hypophthalmichthys molitrix) and their hybrids

Bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix), collectively called bigheaded carps, are invasive species in the Mississippi River Basin (MRB). Interspecific hybridization between bigheaded carps has been considered rare within their native rivers in China; however, it is prevalent in the MRB. We conducted de novo transcriptome analysis of pure and hybrid bigheaded carps and obtained 40,759 to 51,706 transcripts for pure, F₁ hybrid, and backcross bigheaded carps. The search against protein databases resulted in 20,336–28,133 annotated transcripts (over 50% of the transcriptome) with over 13,000 transcripts mapped to 23 Gene Ontology biological processes and 127 KEGG metabolic pathways. More transcripts were detected in silver carp than in bighead carp; however, comparable numbers of transcripts were annotated. Transcriptomic variation detected between two F₁ hybrids may indicate a potential loss of fitness in hybrids. The neighbor‐joining distance tree constructed using over 2,500 one‐to‐one orthologous sequences suggests transcriptomes could be used to infer the history of introgression and hybridization. Moreover, we detected 24,792 candidate SNPs that can be used to identify different species. The transcriptomes, orthologous sequences, and candidate SNPs obtained in this study should provide further knowledge of interspecific hybridization and introgression.

Analysis of the Core Genome and Pan-Genome of Autotrophic Acetogenic Bacteria

Acetogens are obligate anaerobic bacteria capable of reducing carbon dioxide (CO₂) to multicarbon compounds coupled to the oxidation of inorganic substrates, such as hydrogen (H₂) or carbon monoxide (CO), via the Wood-Ljungdahl pathway. Owing to the metabolic capability of CO₂ fixation, much attention has been focused on understanding the unique pathways associated with acetogens, particularly their metabolic coupling of CO₂ fixation to energy conservation. Most known acetogens are phylogenetically and metabolically diverse bacteria present in 23 different bacterial genera. With the increased volume of available genome information, acetogenic bacterial genomes can be analyzed by comparative genome analysis. Even with the genetic diversity that exists among acetogens, the Wood-Ljungdahl pathway, a central metabolic pathway, and cofactor biosynthetic pathways are highly conserved for autotrophic growth. Additionally, comparative genome analysis revealed that most genes in the acetogen-specific core genome were associated with the Wood-Ljungdahl pathway. The conserved enzymes and those predicted as missing can provide insight into biological differences between acetogens and allow for the discovery of promising candidates for industrial applications.

Protein domain architectures provide a fast, efficient and scalable alternative to sequence-based methods for comparative functional genomics

A functional comparative genome analysis is essential to understand the mechanisms underlying bacterial evolution and adaptation. Detection of functional orthologs using standard global sequence similarity methods faces several problems; the need for defining arbitrary acceptance thresholds for similarity and alignment length, lateral gene acquisition and the high computational cost for finding bi-directional best matches at a large scale. We investigated the use of protein domain architectures for large scale functional comparative analysis as an alternative method. The performance of both approaches was assessed through functional comparison of 446 bacterial genomes sampled at different taxonomic levels. We show that protein domain architectures provide a fast and efficient alternative to methods based on sequence similarity to identify groups of functionally equivalent proteins within and across taxonomic boundaries, and it is suitable for large scale comparative analysis. Running both methods in parallel pinpoints potential functional adaptations that may add to bacterial fitness.

Global deceleration of gene evolution following recent genome hybridizations in fungi

Polyploidization events such as whole-genome duplication and inter-species hybridization are major evolutionary forces that shape genomes. Although long-term effects of polyploidization have been well-characterized, early molecular evolutionary consequences of polyploidization remain largely unexplored. Here, we report the discovery of two recent and independent genome hybridizations within a single clade of a fungal genus, Trichosporon. Comparative genomic analyses revealed that redundant genes are experiencing decelerations, not accelerations, of evolutionary rates. We identified a relationship between gene conversion and decelerated evolution suggesting that gene conversion may improve the genome stability of young hybrids by restricting gene functional divergences. Furthermore, we detected large-scale gene losses from transcriptional and translational machineries that indicate a global compensatory mechanism against increased gene dosages. Overall, our findings illustrate counteracting mechanisms during an early phase of post-genome hybridization and fill a critical gap in existing theories on genome evolution.

Mechanical cell competition kills cells via induction of lethal p53 levels

Cell competition is a quality control mechanism that eliminates unfit cells. How cells compete is poorly understood, but it is generally accepted that molecular exchange between cells signals elimination of unfit cells. Here we report an orthogonal mechanism of cell competition, whereby cells compete through mechanical insults. We show that MDCK cells silenced for the polarity gene scribble (scrib^KD) are hypersensitive to compaction, that interaction with wild-type cells causes their compaction and that crowding is sufficient for scrib^KD cell elimination. Importantly, we show that elevation of the tumour suppressor p53 is necessary and sufficient for crowding hypersensitivity. Compaction, via activation of Rho-associated kinase (ROCK) and the stress kinase p38, leads to further p53 elevation, causing cell death. Thus, in addition to molecules, cells use mechanical means to compete. Given the involvement of p53, compaction hypersensitivity may be widespread among damaged cells and offers an additional route to eliminate unfit cells.

Cell competition is a quality control mechanism to eliminate unfit cells. Here the authors show that physical compaction of less fit cells surrounded by healthy neighbours leads to increased expression of tumour suppressor p53 in the compacted cells, causing cell death.