Large-scale Comparative Genomics Research Articles

Severe Bottleneck Impacted the Genomic Structure of Egg-Eating Cichlids in Lake Victoria.

Within 15,000 years, the explosive adaptive radiation of haplochromine cichlids in Lake Victoria, East Africa, generated 500 endemic species. In the 1980s, the upsurge of Nile perch, a carnivorous fish artificially introduced to the lake, drove the extinction of more than 200 endemic cichlids. The Nile perch predation particularly harmed piscivorous cichlids, including paedophages, cichlids eat eggs and fries, which is an example of the unique trophic adaptation seen in African cichlids. Here, aiming to investigate past demographic events possibly triggered by the invasion of Nile perch and the subsequent impacts on the genetic structure of cichlids, we conducted large-scale comparative genomics. We discovered evidence of recent bottleneck events in 4 species, including 2 paedophages, which began during the 1970s to 1980s, and population size rebounded during the 1990s to 2000s. The timing of the bottleneck corresponded to the historical records of endemic haplochromines" disappearance and later resurgence, which is likely associated with the introduction of Nile perch by commercial demand to Lake Victoria in the 1950s. Interestingly, among the 4 species that likely experienced bottleneck, Haplochromis sp. "matumbi hunter," a paedophagous cichlid, showed the most severe bottleneck signatures. The components of shared ancestry inferred by ADMIXTURE suggested a high genetic differentiation between matumbi hunter and other species. In contrast, our phylogenetic analyses highly supported the monophyly of the 5 paedophages, consistent with the results of previous studies. We conclude that high genetic differentiation of matumbi hunter occurred due to the loss of shared genetic components among haplochromines in Lake Victoria caused by the recent severe bottleneck.

FEVER: an interactive web-based resource for evolutionary transcriptomics across fishes.

Teleost fish represent one of the largest and most diverse clades of vertebrates, which makes them great models in various research areas such as ecology and evolution. Recent sequencing endeavors provided high-quality genomes for species covering the main fish evolutionary lineages, opening up large-scale comparative genomics studies. However, transcriptomic data across fish species and organs are heterogenous and have not been integrated with newly sequenced genomes making gene expression quantification and comparative analyses particularly challenging. Thus, resources integrating genomic and transcriptomic data across fish species and organs are still lacking. Here, we present FEVER, a web-based resource allowing evolutionary transcriptomics across species and tissues. First, based on query genes FEVER reconstructs gene trees providing orthologous and paralogous relationships as well as their evolutionary dynamics across 13 species covering the major fish lineages, and 4 model species as evolutionary outgroups. Second, it provides unbiased gene expression across 11 tissues using up-to-date fish genomes. Finally, genomic and transcriptomic data are combined together allowing the exploration of gene expression evolution following speciation and duplication events. FEVER is freely accessible at https://fever.sk8.inrae.fr/.

BacLIFE: a user-friendly computational workflow for genome analysis and prediction of lifestyle-associated genes in bacteria

Bacteria have an extensive adaptive ability to live in close association with eukaryotic hosts, exhibiting detrimental, neutral or beneficial effects on host growth and health. However, the genes involved in niche adaptation are mostly unknown and their functions poorly characterized. Here, we present bacLIFE (https://github.com/Carrion-lab/bacLIFE) a streamlined computational workflow for genome annotation, large-scale comparative genomics, and prediction of lifestyle-associated genes (LAGs). As a proof of concept, we analyzed 16,846 genomes from the Burkholderia/Paraburkholderia and Pseudomonas genera, which led to the identification of hundreds of genes potentially associated with a plant pathogenic lifestyle. Site-directed mutagenesis of 14 of these predicted LAGs of unknown function, followed by plant bioassays, showed that 6 predicted LAGs are indeed involved in the phytopathogenic lifestyle of Burkholderia plantarii and Pseudomonas syringae pv. phaseolicola. These 6 LAGs encompassed a glycosyltransferase, extracellular binding proteins, homoserine dehydrogenases and hypothetical proteins. Collectively, our results highlight bacLIFE as an effective computational tool for prediction of LAGs and the generation of hypotheses for a better understanding of bacteria-host interactions.

Brief Communication: The Predictable Network Topology of Evolutionary Genomic Constraint.

Large-scale comparative genomics studies offer valuable resources for understanding both functional and evolutionary rate constraints. It is suggested that constraint aligns with the topology of genomic networks, increasing toward the center, with intermediate nodes combining relaxed constraint with higher contributions to the phenotype due to pleiotropy. However, this pattern has yet to be demonstrated in vertebrates. This study shows that constraint intensifies toward the network's center in placental mammals. Genes with rate changes associated with emergence of hibernation cluster mostly toward intermediate positions, with higher constraint in faster-evolving genes, which is indicative of a "sweet spot" for adaptation. If this trend holds universally, network node metrics could predict high-constraint regions even in clades lacking empirical constraint data.

A large-scale comparative genomics study reveals niche-driven and within-sample intra-species functional diversification in Lacticaseibacillus rhamnosus

Lacticaseibacillus rhamnosus (L. rhamnosus) is widely recognized as a probiotic species, and it exists in a variety of environments including host gut and dairy products. This work aimed at conducting a large-scale comparative genomics analysis of 384 L. rhamnosus genomes (257 whole-sequence or metagenomic-assembled genomes from gut-associated isolates [122 and 135 retrieved from the UHGG and NCBI databases, respectively] and 127 genomes from dairy isolates [34 from the NCBI database; 93 isolated from a cheese sample and sequenced here]). Our results showed that L. rhamnosus had a large and open pan-genome (15,253 pan-genes identified from all 384 genomes; 15,028 pan-genes if the 93 cheese-originated isolates were excluded). The core-gene phylogenetic tree constructed from the 384 L. rhamnosus genomes comprised five phylogenetic branches, with a random distribution of dairy and gut-associated isolates/genomes across the tree. No significant difference was identified in the overall profile of metabolism-related genes between dairy and gut-associated genomes; however, notably, the gut-associated strains/isolates contained more genes coding for specific metabolic pathways and carbohydrate-active enzymes, e.g., lacto-N-biosidase (EC 3.2.1.140; GT20) and lacto-N-biose phosphorylase/galacto-N-biose phosphorylase (EC 2.4.1.211; GH112). Further, we found that there was obvious intra-species diversification of the 93 cheese-originated L. rhamnosus isolates, forming three clades (Clades A, B, and C) in the reconstructed core-gene phylogenetic tree. There were numerous single nucleotide variations (over 10,000) across the three clades. Moreover, significant differences were observed in the content of metabolism-related genes across clades (p < 0.05, Adonis test), characterized by the enrichment in glycoside hydrolases in Clade C and the possession of unique metabolic pathways in each clade. These results implicated genomics/functional diversification of L. rhamnosus in a single food matrix and niche-driven adaptive evolution of isolates from dairy and host gut-associated origins. Our study shed insights into the selection of candidate strains for food industry applications.

Minmers are a generalization of minimizers that enable unbiased local jaccard estimation.

The Jaccard similarity on k-mer sets has shown to be a convenient proxy for sequence identity. By avoiding expensive base-level alignments and comparing reduced sequence representations, tools such as MashMap can scale to massive numbers of pairwise comparisons while still providing useful similarity estimates. However, due to their reliance on minimizer winnowing, previous versions of MashMap were shown to be biased and inconsistent estimators of Jaccard similarity. This directly impacts downstream tools that rely on the accuracy of these estimates. To address this, we propose the minmer winnowing scheme, which generalizes the minimizer scheme by use of a rolling minhash with multiple sampled k-mers per window. We show both theoretically and empirically that minmers yield an unbiased estimator of local Jaccard similarity, and we implement this scheme in an updated version of MashMap. The minmer-based implementation is over 10 times faster than the minimizer-based version under the default ANI threshold, making it well-suited for large-scale comparative genomics applications. MashMap3 is available at https://github.com/marbl/MashMap.

Large-scale comparative genomics to refine the organization of the global Salmonella enterica population structure.

The White-Kauffmann-Le Minor (WKL) scheme is the most widely used Salmonella typing scheme for reporting the disease prevalence of the enteric pathogen. With the advent of whole-genome sequencing (WGS), in silico methods have increasingly replaced traditional serotyping due to reproducibility, speed and coverage. However, despite integrating genomic-based typing by in silico serotyping tools such as SISTR, in silico serotyping in certain contexts remains ambiguous and insufficiently informative. Specifically, in silico serotyping does not attempt to resolve polyphyly. Furthermore, in spite of the widespread acknowledgement of polyphyly from genomic studies, the prevalence of polyphyletic serovars is not well characterized. Here, we applied a genomics approach to acquire the necessary resolution to classify genetically discordant serovars and propose an alternative typing scheme that consistently reflect natural Salmonella populations. By accessing the unprecedented volume of bacterial genomic data publicly available in GenomeTrakr and PubMLST databases (>180 000 genomes representing 723 serovars), we characterized the global Salmonella population structure and systematically identified putative non-monophyletic serovars. The proportion of putative non-monophyletic serovars was estimated higher than previous reports, reinforcing the inability of antigenic determinants to depict the complexity of Salmonella evolutionary history. We explored the extent of genetic diversity masked by serotyping labels and found significant intra-serovar molecular differences across many clinically important serovars. To avoid false discovery due to incorrect in silico serotyping calls, we cross-referenced reported serovar labels and concluded a low error rate in in silico serotyping. The combined application of clustering statistics and genome-wide association methods demonstrated effective characterization of stable bacterial populations and explained functional differences. The collective methods adopted in our study have practical values in establishing genomic-based typing nomenclatures for an entire microbial species or closely related subpopulations. Ultimately, we foresee an improved typing scheme to be a hybrid that integrates both genomic and antigenic information such that the resolution from WGS is leveraged to improve the precision of subpopulation classification while preserving the common names defined by the WKL scheme.

Twenty-five years of Genomes OnLine Database (GOLD): data updates and new features in v.9.

The Genomes OnLine Database (GOLD) (https://gold.jgi.doe.gov/) at the Department of Energy Joint Genome Institute (DOE-JGI) continues to maintain its role as one of the flagship genomic metadata repositories of the world. The ever-increasing number of projects and metadata are freely available to the user community world-wide. GOLD's metadata is consumed by scientists and remains an important source for large-scale comparative genomics analysis initiatives. Encouraged by this active user engagement and growth, GOLD has continued to add new components and capabilities. The new features such as a public Application Programming Interface (API) and Ecosystem landing page as well as the growth of different entities in this current GOLD v.9 edition are described in detail in this manuscript.

VIVID: A Web Application for Variant Interpretation and Visualization in Multi-dimensional Analyses.

Large-scale comparative genomics- and population genetic studies generate enormous amounts of polymorphism data in the form of DNA variants. Ultimately, the goal of many of these studies is to associate genetic variants to phenotypes or fitness. We introduce VIVID, an interactive, user-friendly web application that integrates a wide range of approaches for encoding genotypic to phenotypic information in any organism or disease, from an individual or population, in three-dimensional (3D) space. It allows mutation mapping and annotation, calculation of interactions and conservation scores, prediction of harmful effects, analysis of diversity and selection, and 3D visualization of genotypic information encoded in Variant Call Format on AlphaFold2 protein models. VIVID enables the rapid assessment of genes of interest in the study of adaptive evolution and the genetic load, and it helps prioritizing targets for experimental validation. We demonstrate the utility of VIVID by exploring the evolutionary genetics of the parasitic protist Plasmodium falciparum, revealing geographic variation in the signature of balancing selection in potential targets of functional antibodies.

Thousands of small, novel genes predicted in global phage genomes.

Small genes (<150 nucleotides) have been systematically overlooked in phage genomes. We employ a large-scale comparative genomics approach to predict >40,000 small-gene families in ∼2.3 million phage genome contigs. We find that small genes in phage genomes are approximately 3-fold more prevalent than in host prokaryotic genomes. Our approach enriches for small genes that are translated in microbiomes, suggesting the small genes identified are coding. More than 9,000 families encode potentially secreted or transmembrane proteins, more than 5,000 families encode predicted anti-CRISPR proteins, and more than 500 families encode predicted antimicrobial proteins. By combining homology and genomic-neighborhood analyses, we reveal substantial novelty and diversity within phage biology, including small phage genes found in multiple host phyla, small genes encoding proteins that play essential roles in host infection, and small genes that share genomic neighborhoods and whose encoded proteins may share related functions.

A Large-Scale Genome-Based Survey of Acidophilic Bacteria Suggests That Genome Streamlining Is an Adaption for Life at Low pH.

The genome streamlining theory suggests that reduction of microbial genome size optimizes energy utilization in stressful environments. Although this hypothesis has been explored in several cases of low-nutrient (oligotrophic) and high-temperature environments, little work has been carried out on microorganisms from low-pH environments, and what has been reported is inconclusive. In this study, we performed a large-scale comparative genomics investigation of more than 260 bacterial high-quality genome sequences of acidophiles, together with genomes of their closest phylogenetic relatives that live at circum-neutral pH. A statistically supported correlation is reported between reduction of genome size and decreasing pH that we demonstrate is due to gene loss and reduced gene sizes. This trend is independent from other genome size constraints such as temperature and G + C content. Genome streamlining in the evolution of acidophilic bacteria is thus supported by our results. The analyses of predicted Clusters of Orthologous Genes (COG) categories and subcellular location predictions indicate that acidophiles have a lower representation of genes encoding extracellular proteins, signal transduction mechanisms, and proteins with unknown function but are enriched in inner membrane proteins, chaperones, basic metabolism, and core cellular functions. Contrary to other reports for genome streamlining, there was no significant change in paralog frequencies across pH. However, a detailed analysis of COG categories revealed a higher proportion of genes in acidophiles in the following categories: “replication and repair,” “amino acid transport,” and “intracellular trafficking”. This study brings increasing clarity regarding the genomic adaptations of acidophiles to life at low pH while putting elements, such as the reduction of average gene size, under the spotlight of streamlining theory.

Genome Informatics and Machine Learning-Based Identification of Antimicrobial Resistance-Encoding Features and Virulence Attributes in Escherichia coli Genomes Representing Globally Prevalent Lineages, Including High-Risk Clonal Complexes.

ABSTRACTEscherichia coli, a ubiquitous commensal/pathogenic member from the Enterobacteriaceae family, accounts for high infection burden, morbidity, and mortality throughout the world. With emerging multidrug resistance (MDR) on a massive scale, E. coli has been listed as one of the Global Antimicrobial Resistance and Use Surveillance System (GLASS) priority pathogens. Understanding the resistance mechanisms and underlying genomic features appears to be of utmost importance to tackle further spread of these multidrug-resistant superbugs. While a few of the globally prevalent sequence types (STs) of E. coli, such as ST131, ST69, ST405, and ST648, have been previously reported to be highly virulent and harboring MDR, there is no clarity if certain ST lineages have a greater propensity to acquire MDR. In this study, large-scale comparative genomics of a total of 5,653 E. coli genomes from 19 ST lineages revealed ST-wide prevalence patterns of genomic features, such as antimicrobial resistance (AMR)-encoding genes/mutations, virulence genes, integrons, and transposons. Interpretation of the importance of these features using a Random Forest Classifier trained with 11,988 genomic features from whole-genome sequence data identified ST-specific or phylogroup-specific signature proteins mostly belonging to different protein superfamilies, including the toxin-antitoxin systems. Our study provides a comprehensive understanding of a myriad of genomic features, ST-specific proteins, and resistance mechanisms entailing different lineages of E. coli at the level of genomes; this could be of significant downstream importance in understanding the mechanisms of AMR, in clinical discovery, in epidemiology, and in devising control strategies.

A chromosome-scale Gastrodia elata genome and large-scale comparative genomic analysis indicate convergent evolution by gene loss in mycoheterotrophic and parasitic plants.

Mycoheterotrophic and parasitic plants are heterotrophic and parasitize on fungi and plants, respectively, to obtain nutrients. Large-scale comparative genomics analysis has not been conducted in mycoheterotrophic or parasitic plants or between these two groups of parasites. We assembled a chromosome-level genome of the fully mycoheterotrophic plant Gastrodia elata (Orchidaceae) and performed comparative genomic analyses on the genomes of G. elata and four orchids (initial mycoheterotrophs), three parasitic plants (Cuscuta australis, Striga asiatica, and Sapria himalayana), and 36 autotrophs from various angiosperm lineages. It was found that while in the hemiparasite S. asiatica and initial mycoheterotrophic orchids, approximately 4-5% of the conserved orthogroups were lost, the fully heterotrophic G. elata and C. australis both lost approximately 10% of the conserved orthogroups, indicating that increased heterotrophy is positively associated with gene loss. Importantly, many genes that are essential for autotrophs, including those involved in photosynthesis, the circadian clock, flowering time regulation, immunity, nutrient uptake, and root and leaf development, were convergently lost in both G. elata and C. australis. The high-quality genome of G. elata will facilitate future studies on the physiology, ecology, and evolution of mycoheterotrophic plants, and our findings highlight the critical role of gene loss in the evolution of plants with heterotrophic lifestyles.

Identification of new serovar-specific detection targets against salmonella B serogroup using large-scale comparative genomics

Salmonella serotyping is a critical step for the identification of pathological isolates in humans and animals. Whole-genome analysis is a proven method for Salmonella serovars identification. In this study, we analysed 1668 genomic sequences of 23 serogroup B serovars and 129 genomic sequences of 80 non-target serovars to detect serovar-specific molecular targets for serogroup B (O:4 epitopes). Serogroup B is the most common serogroup found in China, whose presence is increasing in recent years. We identified five, eight, one, and seven new serovar-specific genes for Indiana, Derby, Typhimurium, and Agona serovars, respectively. We then designed specific PCR primer sets for these genes to validate our results. All target Salmonella serovars could be detected successfully by PCR in different types of food samples. The analytical sensitivity of PCR for 19 out of 21 pairs of primers was 101–103 CFU/mL and 102–105 CFU/g in pure cultures of the bacteria and artificially contaminated fresh pork samples, respectively. Furthermore, we performed a functional analysis of the newly-identified serovar-specific genes. We found that the function was unknown for most of the genes, although for some of them we were able to identify their generic functional categories (membrane protein genes and secreted proteins genes). Our findings revealed that whole-genome analysis is a valuable tool to identify new serovar-specific targets, which could be used to develop rapid and accurate molecular detection methods and provide additional information about different Salmonella serogroup B serovars.

The genus bifidobacterium: From genomics to functionality of an important component of the mammalian gut microbiota running title: Bifidobacterial adaptation to and interaction with the host.

Members of the genus Bifidobacterium are dominant and symbiotic inhabitants of the mammalian gastrointestinal tract. Being vertically transmitted, bifidobacterial host colonization commences immediately after birth and leads to a phase of host infancy during which bifidobacteria are highly prevalent and abundant to then transit to a reduced, yet stable abundance phase during host adulthood. However, in order to reach and stably colonize their elective niche, i.e. the large intestine, bifidobacteria have to cope with a multitude of oxidative, osmotic and bile salt/acid stress challenges that occur along the gastrointestinal tract (GIT). Concurrently, bifidobacteria not only have to compete with the myriad of other gut commensals for nutrient acquisition, but they also require protection against bacterial viruses. In this context, Next-Generation Sequencing (NGS) techniques, allowing large-scale comparative and functional genome analyses have helped to identify the genetic strategies that bifidobacteria have developed in order to colonize, survive and adopt to the highly competitive mammalian gastrointestinal environment. The current review is aimed at providing a comprehensive overview concerning the molecular strategies on which bifidobacteria rely to stably and successfully colonize the mammalian gut.

Myanmar Burkholderia pseudomallei strains are genetically diverse and originate from Asia with phylogenetic evidence of reintroductions from neighbouring countries

Melioidosis was first identified in Myanmar in 1911 but for the last century it has remained largely unreported there. Burkholderia pseudomallei was first isolated from the environment of Myanmar in 2016, confirming continuing endemicity. Recent genomic studies showed that B. pseudomallei originated in Australia and spread to Asia, with phylogenetic evidence of repeated reintroduction of B. pseudomallei across countries bordered by the Mekong River and the Malay Peninsula. We present the first whole-genome sequences of B. pseudomallei isolates from Myanmar: nine clinical and seven environmental isolates. We used large-scale comparative genomics to assess the genetic diversity, phylogeography and potential origins of B. pseudomallei in Myanmar. Global phylogenetics demonstrated that Myanmar isolates group in two distantly related clades that reside in a more ancestral Asian clade with high amounts of genetic diversity. The diversity of B. pseudomallei from Myanmar and divergence within our global phylogeny suggest that the original introduction of B. pseudomallei to Myanmar was not a recent event. Our study provides new insights into global patterns of B. pseudomallei dissemination, most notably the dynamic nature of movement of B. pseudomallei within densely populated Southeast Asia. The role of anthropogenic influences in both ancient and more recent dissemination of B. pseudomallei to Myanmar and elsewhere in Southeast Asia and globally requires further study.

A Chromosome-Scale Genome Assembly for the Fusarium oxysporum Strain Fo5176 To Establish a Model Arabidopsis-Fungal Pathosystem.

Plant pathogens cause widespread yield losses in agriculture. Understanding the drivers of plant-pathogen interactions requires decoding the molecular dialog leading to either resistance or disease. However, progress in deciphering pathogenicity genes has been severely hampered by suitable model systems and incomplete fungal genome assemblies. Here, we report a significant improvement of the assembly and annotation of the genome of the Fusarium oxysporum (Fo) strain Fo5176. Fo comprises a large number of serious plant pathogens on dozens of plant species with largely unresolved pathogenicity factors. The strain Fo5176 infects Arabidopsis thaliana and, hence, constitutes a highly promising model system. We use high-coverage Pacific Biosciences Sequel long-read and Hi-C sequencing data to assemble the genome into 19 chromosomes and a total genome size of 67.98 Mb. The genome has a N50 of 4 Mb and a 99.1% complete BUSCO score. Phylogenomic analyses based on single-copy orthologs clearly place the Fo5176 strain in the Fo f sp. conglutinans clade as expected. We generated RNAseq data from culture medium and plant infections to train gene predictions and identified ∼18,000 genes including ten effector genes known from other Fo clades. We show that Fo5176 is able to infect cabbage and Brussel sprouts of the Brassica oleracea, expanding the usefulness of the Fo5176 model pathosystem. Finally, we performed large-scale comparative genomics analyses comparing the Fo5176 to 103 additional Fo genomes to define core and accessory genomic regions. In conjunction with the molecular tool sets available for A. thaliana, the Fo5176 genome and annotation provides a crucial step toward the establishment of a highly promising pathosystem.

Phylogenetic and genomic analysis reveals high genomic openness and genetic diversity of Clostridium perfringens.

Clostridium perfringens is associated with a variety of diseases in both humans and animals. Recent advances in genomic sequencing make it timely to re-visit this important pathogen. Although the genome sequence of C. perfringens was first determined in 2002, large-scale comparative genomics with isolates of different origins is still lacking. In this study, we used whole-genome sequencing of 45 C . perfringens isolates with isolation time spanning an 80‐year period and performed comparative analysis of 173 genomes from worldwide strains. We also conducted phylogenetic lineage analysis and introduced an openness index (OI) to evaluate the openness of bacterial genomes. We classified all these genomes into five lineages and hypothesized that the origin of C. perfringens dates back to ~80 000 years ago. We showed that the pangenome of the 173 C . perfringens strains contained a total of 26 954 genes, while the core genome comprised 1020 genes, accounting for about a third of the genome of each isolate. We demonstrated that C. perfringens had the highest OI compared with 51 other bacterial species. Intact prophage sequences were found in nearly 70.0 % of C. perfringens genomes, while CRISPR sequences were found only in ~40.0 %. Plasmids were prevalent in C. perfringens isolates, and half of the virulence genes and antibiotic resistance genes (ARGs) identified in all the isolates could be found in plasmids. ARG-sharing network analysis showed that C. perfringens shared its 11 ARGs with 55 different bacterial species, and a high frequency of ARG transfer may have occurred between C. perfringens and species in the genera Streptococcus and Staphylococcus . Correlation analysis showed that the ARG number in C. perfringens strains increased with time, while the virulence gene number was relative stable. Our results, taken together with previous studies, revealed the high genome openness and genetic diversity of C. perfringens and provide a comprehensive view of the phylogeny, genomic features, virulence gene and ARG profiles of worldwide strains.

Evol2Circos: A web-based tool for genome synteny and collinearity analysis and its visualization in fishes.

The advent of high throughput next generation sequencing technologies and improved assembly algorithms have ensued in accumulation of voluminous genomic data in public domains. It has opened up entries for large scale comparative genome studies, especially the identification of conserved syntenic blocks among the species, facilitating the evolutionary importance of the conservation and variation in genomic organization. Synteny construction and visualization requires computational and bioinformatics skills to prepare input file for the synteny analysis pipeline. The syntenic information in fishes is still in juvenile stage and are scattered in different research domains. Here, we present a web-based tool 'Evol2Circos' to provide a user-friendly GUI- and web-based tool to analyse user specific data for synteny construction and visualization, and to facilitate the browsing of syntenic information of different fishes using the circos, bar, dual and dot plots. The information generated from the tool can also be used for further downstream analyses. Evol2Circos software tool is tested under Ubuntu Linux. The web-browser, source code, documentation, user manual, example dataset and scripts are available online at: 203.190.147.148/evole2circos/.

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.