Comparative genomics reveals insights into cyanobacterial evolution and habitat adaptation

Cyanobacteria are photosynthetic prokaryotes that inhabit diverse aquatic and terrestrial environments. However, the evolutionary mechanisms involved in the cyanobacterial habitat adaptation remain poorly understood. Here, based on phylogenetic and comparative genomic analyses of 650 cyanobacterial genomes, we investigated the genetic basis of cyanobacterial habitat adaptation (marine, freshwater, and terrestrial). We show: (1) the expansion of gene families is a common strategy whereby terrestrial cyanobacteria cope with fluctuating environments, whereas the genomes of many marine strains have undergone contraction to adapt to nutrient-poor conditions. (2) Hundreds of genes are strongly associated with specific habitats. Genes that are differentially abundant in genomes of marine, freshwater, and terrestrial cyanobacteria were found to be involved in light sensing and absorption, chemotaxis, nutrient transporters, responses to osmotic stress, etc., indicating the importance of these genes in the survival and adaptation of organisms in specific habitats. (3) A substantial fraction of genes that facilitate the adaptation of Cyanobacteria to specific habitats are contributed by horizontal gene transfer, and such genetic exchanges are more frequent in terrestrial cyanobacteria. Collectively, our results further our understandings of the adaptations of Cyanobacteria to different environments, highlighting the importance of ecological constraints imposed by the environment in shaping the evolution of Cyanobacteria.Subject terms: Phylogenetics, Microbial ecology     

Comparative genomic analysis reveals varying levels of mammalian adaptation to coronavirus infections

Severe acute respiratory coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is of zoonotic origin. Evolutionary analyses assessing whether coronaviruses similar to SARS-CoV-2 infected ancestral species of modern-day animal hosts could be useful in identifying additional reservoirs of potentially dangerous coronaviruses. We reasoned that if a clade of species has been repeatedly exposed to a virus, then their proteins relevant for viral entry may exhibit adaptations that affect host susceptibility or response. We perform comparative analyses across the mammalian phylogeny of angiotensin-converting enzyme 2 (ACE2), the cellular receptor for SARS-CoV-2, in order to uncover evidence for selection acting at its binding interface with the SARS-CoV-2 spike protein. We uncover that in rodents there is evidence for adaptive amino acid substitutions at positions comprising the ACE2-spike interaction interface, whereas the variation within ACE2 proteins in primates and some other mammalian clades is not consistent with evolutionary adaptations. We also analyze aminopeptidase N (APN), the receptor for the human coronavirus 229E, a virus that causes the common cold, and find evidence for adaptation in primates. Altogether, our results suggest that the rodent and primate lineages may have had ancient exposures to viruses similar to SARS-CoV-2 and HCoV-229E, respectively.     

Comparative genomics reveals novel biochemical pathways

How well do we understand which enzymes are involved in the primary metabolism of the cell? A recent study using comparative genomics and postgenomics approaches revealed a novel pathway in the most studied organism, Escherichia coli. The analysis of a new operon consisting of seven previously uncharacterized genes thought to be involved in the degradation of nucleic acid precursors shows the impact of comparative genomics on the discovery of novel pathways and enzymes.     

Comparative genomics of Riemerella anatipestifer reveals genetic diversity

Background

Riemerella anatipestifer is one of the most important pathogens of ducks. However, the molecular mechanisms of R. anatipestifer infection are poorly understood. In particular, the lack of genomic information from a variety of R. anatipestifer strains has proved severely limiting.

Results

In this study, we present the complete genomes of two R. anatipestifer strains, RA-CH-1 (2,309,519 bp, Genbank accession {"type":"entrez-nucleotide","attrs":{"text":"CP003787","term_id":"403311769","term_text":"CP003787"}}CP003787) and RA-CH-2 (2,166,321 bp, Genbank accession {"type":"entrez-nucleotide","attrs":{"text":"CP004020","term_id":"441482619","term_text":"CP004020"}}CP004020). Both strains are from isolates taken from two different sick ducks in the SiChuang province of China. A comparative genomics approach was used to identify similarities and key differences between RA-CH-1 and RA-CH-2 and the previously sequenced strain RA-GD, a clinical isolate from GuangDong, China, and ATCC11845.

Conclusion

The genomes of RA-CH-2 and RA-GD were extremely similar, while RA-CH-1 was significantly different than ATCC11845. RA-CH-1 is 140,000 bp larger than the three other strains and has 16 unique gene families. Evolutionary analysis shows that RA-CH-1 and RA-CH-2 are closed and in a branch with ATCC11845, while RA-GD is located in another branch. Additionally, the detection of several iron/heme-transport related proteins and motility mechanisms will be useful in elucidating factors important in pathogenicity. This information will allow a better understanding of the phenotype of different R. anatipestifer strains and molecular mechanisms of infection.     

Comparative genomic analysis of inbred rat strains reveals the existence of ancestral polymorphisms

Mammalian Genome - In an alignment of closely related genomic sequences, the existence of discordant mutation sites, which do not reflect the phylogenetic relationship of the genomes, is often...     

Comparative genomic analysis reveals novel facts about Leptospirillum spp. cytochromes

Chemolithoautotrophic acidophilic bacteria, which belong to the genus Leptospirillum, can only grow with Fe(II) as electron donor and oxygen as an electron acceptor. Members of this genus play an important role in bioleaching sulfide ores. We used nearly complete genome sequences of Leptospirillum ferrooxidans (group I), Leptospirillum rubarum, Leptospirillum '5-way CG' (group II) and Leptospirillum ferrodiazotrophum (group III) to identify cytochromes that are likely involved in electron transfer chain(s). The results show the presence of genes encoding a number of c-type cytochromes (18-20 genes were identified in each species), as well as bd and cbb? oxidases. Genes encoding cbb? oxidase are clustered, with predicted genes involved in cbb? maturation proteins. Duplication of cbb? encoding genes (ccoNO) was detected in all four genomes. Interestingly, these micro-organisms also contain genes that potentially encode bc? and b?f-like complexes organized into two putative operon structures. To date, the Leptospirillum genus includes the only organisms reported to have genes coding for two different bc complexes. This study provides detailed insights into the components of electron transfer chains of Leptospirillum spp., revealing their conservation among leptospirilla groups and suggesting that there may be a single common pathway for electron transport between Fe(II) and oxygen.     

Comparative genomic analysis reveals evolutionary characteristics and patterns of microRNA clusters in vertebrates

MicroRNAs (miRNAs) are a class of small non-coding RNAs that can play important regulatory roles in many important biological processes. Although clustering patterns of miRNA clusters have been uncovered in animals, the origin and evolution of miRNA clusters in vertebrates are still poorly understood. Here, we performed comparative genomic analyses to construct 51 sets of orthologous miRNA clusters (SOMCs) across seven test vertebrate species, a collection of miRNA clusters from two or more species that are likely to have evolved from a common ancestral miRNA cluster, and used these to systematically examine the evolutionary characteristics and patterns of miRNA clusters in vertebrates. We found that miRNA clusters are continuously generated, and most of them tend to be conserved and maintained in vertebrate genomes, although some adaptive gains and losses of miRNA cluster have occurred during evolution. Furthermore, miRNA clusters appeared relatively early in the evolutionary history might suffer from more complicated adaptive gain-and-loss than those young miRNA clusters. Detailed analysis showed that genomic duplication events of ancestral miRNAs or miRNA clusters are likely to be major driving force and apparently contribute to origin and evolution of miRNA clusters. Comparison of conserved with lineage-specific miRNA clusters revealed that the contribution of duplication events for the formation of miRNA cluster appears to be more important for conserved miRNA clusters than lineage-specific. Our study provides novel insights for further exploring the origins and evolution of miRNA clusters in vertebrates at a genome scale.     

Comparative genomic analysis of Acinetobacter baumannii clinical isolates reveals extensive genomic variation and diverse antibiotic resistance determinants

Background

Acinetobacter baumannii is an important nosocomial pathogen that poses a serious health threat to immune-compromised patients. Due to its rapid ability to develop multidrug resistance (MDR), A. baumannii has increasingly become a focus of attention worldwide. To better understand the genetic variation and antibiotic resistance mechanisms of this bacterium at the genomic level, we reported high-quality draft genome sequences of 8 clinical isolates with various sequence types and drug susceptibility profiles.

Results

We sequenced 7 MDR and 1 drug-sensitive clinical A. baumannii isolates and performed comparative genomic analysis of these draft genomes with 16 A. baumannii complete genomes from GenBank. We found a high degree of variation in A. baumannii, including single nucleotide polymorphisms (SNPs) and large DNA fragment variations in the AbaR-like resistance island (RI) regions, the prophage and the type VI secretion system (T6SS). In addition, we found several new AbaR-like RI regions with highly variable structures in our MDR strains. Interestingly, we found a novel genomic island (designated as GI_BJ4) in the drug-sensitive strain BJ4 carrying metal resistance genes instead of antibiotic resistance genes inserted into the position where AbaR-like RIs commonly reside in other A. baumannii strains. Furthermore, we showed that diverse antibiotic resistance determinants are present outside the RIs in A. baumannii, including antibiotic resistance-gene bearing integrons, the bla_OXA-23-containing transposon Tn2009, and chromosomal intrinsic antibiotic resistance genes.

Conclusions

Our comparative genomic analysis revealed that extensive genomic variation exists in the A. baumannii genome. Transposons, genomic islands and point mutations are the main contributors to the plasticity of the A. baumannii genome and play critical roles in facilitating the development of antibiotic resistance in the clinical isolates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1163) contains supplementary material, which is available to authorized users.     

Comparative genomic analysis of Ralstonia solanacearum reveals candidate genes for host specificity

Background

Ralstonia solanacearum is a vascular soil-borne plant pathogen with an unusually broad host range. This economically destructive and globally distributed bacterium has thousands of distinct lineages within a heterogeneous and taxonomically disputed species complex. Some lineages include highly host-adapted strains (ecotypes), such as the banana Moko disease-causing strains, the cold-tolerant potato brown rot strains (also known as R3bv2) and the recently emerged Not Pathogenic to Banana (NPB) strains.

Results

These distinct ecotypes offer a robust model to study host adaptation and the emergence of ecotypes because the polyphyletic Moko strains include lineages that are phylogenetically close to the monophyletic brown rot and NPB strains. Draft genomes of eight new strains belonging to these three model ecotypes were produced to complement the eleven publicly available R. solanacearum genomes. Using a suite of bioinformatics methods, we searched for genetic and evolutionary features that distinguish ecotypes and propose specific hypotheses concerning mechanisms of host adaptation in the R. solanacearum species complex. Genome-wide, few differences were identified, but gene loss events, non-synonymous polymorphisms, and horizontal gene transfer were identified among type III effectors and were associated with host range differences.

Conclusions

This extensive comparative genomics analysis uncovered relatively few divergent features among closely related strains with contrasting biological characteristics; however, several virulence factors were associated with the emergence of Moko, NPB and brown rot and could explain host adaptation.

Electronic supplementary material

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

Comparative genomics reveals molecular features unique to the songbird lineage

Background

Songbirds (oscine Passeriformes) are among the most diverse and successful vertebrate groups, comprising almost half of all known bird species. Identifying the genomic innovations that might be associated with this success, as well as with characteristic songbird traits such as vocal learning and the brain circuits that underlie this behavior, has proven difficult, in part due to the small number of avian genomes available until recently. Here we performed a comparative analysis of 48 avian genomes to identify genomic features that are unique to songbirds, as well as an initial assessment of function by investigating their tissue distribution and predicted protein domain structure.

Results

Using BLAT alignments and gene synteny analysis, we curated a large set of Ensembl gene models that were annotated as novel or duplicated in the most commonly studied songbird, the Zebra finch (Taeniopygia guttata), and then extended this analysis to 47 additional avian and 4 non-avian genomes. We identified 10 novel genes uniquely present in songbird genomes. A refined map of chromosomal synteny disruptions in the Zebra finch genome revealed that the majority of these novel genes localized to regions of genomic instability associated with apparent chromosomal breakpoints. Analyses of in situ hybridization and RNA-seq data revealed that a subset of songbird-unique genes is expressed in the brain and/or other tissues, and that 2 of these (YTHDC2L1 and TMRA) are highly differentially expressed in vocal learning-associated nuclei relative to the rest of the brain.

Conclusions

Our study reveals novel genes unique to songbirds, including some that may subserve their unique vocal control system, substantially improves the quality of Zebra finch genome annotations, and contributes to a better understanding of how genomic features may have evolved in conjunction with the emergence of the songbird lineage.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1082) contains supplementary material, which is available to authorized users.     

Comparative genomics reveals differences in algal galactan biosynthesis and related pathways in early and late diverging red algae

Red algae are a major source of marine sulfated galactans. In this study, orthologs and inparalogs from seven red algae were analyzed and compared with the aim to discover differences in algal galactan biosynthesis and related pathways of these algae. Red algal orthologs for putative carbohydrate sulfotransferases were found to be prevalent in Porphyridium purpureum, Florideophytes and Bangiophytes, while red algal orthologs for putative chondroitin sulfate synthases, sulfurylases, and porphyranases /carrageenases were found exclusively in Florideophytes and Bangiophytes. The acquirement of these genes could have happened after the divergence from Cyanidiales red algae. Cyanidiales red algae were found to have more number and types of putative sulfate permeases, suggesting that these genes could have been acquired in adaptation to the environmental stresses and biogeochemistry of respective habitats. The findings of this study shed lights on the evolution of different homeostasis mechanisms by the early and late diverging red algal orders.     

Comparative genomics for reliable protein-function prediction from genomic data

Genomic data provide invaluable, yet unreliable information about protein function. However, if the overlap in information among various genomic datasets is taken into account, one observes an increase in the reliability of the protein-function predictions that can be made. Recently published approaches achieved this either by comparing the same type of data from multiple species (horizontal comparative genomics) or by using subtle, Bayesian methods to compare different types of genomic data from a single species (vertical comparative genomics). In this article, we discuss these methods, illustrating horizontal comparative genomics by comparing yeast two-hybrid (Y2H) data from Saccharomyces cerevisiae with Y2H data from Drosophila melanogaster, and illustrating vertical comparative genomics by comparing RNA expression data with proteomic data from Plasmodium falciparum.     

Comparative genomics analysis of completely sequenced microbial genomes reveals the ubiquity of N-linked glycosylation in prokaryotes

Glycosylation of proteins in prokaryotes has been known for the last few decades. Glycan structures and/or the glycosylation pathways have been experimentally characterized in only a small number of prokaryotes. Even this has become possible only during the last decade or so, primarily due to technological and methodological developments. Glycosylated proteins are diverse in their function and localization. Glycosylation has been shown to be associated with a wide range of biological phenomena. Characterization of the various types of glycans and the glycosylation machinery is critical to understand such processes. Such studies can help in the identification of novel targets for designing drugs, diagnostics, and engineering of therapeutic proteins. In view of this, the experimentally characterized pgl system of Campylobacter jejuni, responsible for N-linked glycosylation, has been used in this study to identify glycosylation loci in 865 prokaryotes whose genomes have been completely sequenced. Results from the present study show that only a small number of organisms have homologs for all the pgl enzymes and a few others have homologs for none of the pgl enzymes. Most of the organisms have homologs for only a subset of the pgl enzymes. There is no specific pattern for the presence or absence of pgl homologs vis-à-vis the 16S rRNA sequence-based phylogenetic tree. This may be due to differences in the glycan structures, high sequence divergence, horizontal gene transfer or non-orthologous gene displacement. Overall, the presence of homologs for pgl enzymes in a large number of organisms irrespective of their habitat, pathogenicity, energy generation mechanism, etc., hints towards the ubiquity of N-linked glycosylation in prokaryotes.     

Comparative genomics reveals unusually long motifs in mammalian genomes

MOTIVATION: The recent discovery of the first small modulatory RNA (smRNA) presents the challenge of finding other molecules of similar length and conservation level. Unlike short interfering RNA (siRNA) and micro-RNA (miRNA), effective computational and experimental screening methods are not currently known for this species of RNA molecule, and the discovery of the one known example was partly fortuitous because it happened to be complementary to a well-studied DNA binding motif (the Neuron Restrictive Silencer Element). RESULTS: The existing comparative genomics approaches (e.g., phylogenetic footprinting) rely on alignments of orthologous regions across multiple genomes. This approach, while extremely valuable, is not suitable for finding motifs with highly diverged "non-alignable" flanking regions. Here we show that several unusually long and well conserved motifs can be discovered de novo through a comparative genomics approach that does not require an alignment of orthologous upstream regions. These motifs, including Neuron Restrictive Silencer Element, were missed in recent comparative genomics studies that rely on phylogenetic footprinting. While the functions of these motifs remain unknown, we argue that some may represent biologically important sites. AVAILABILITY: Our comparative genomics software, a web-accessible database of our results and a compilation of experimentally validated binding sites for NRSE can be found at http://www.cse.ucsd.edu/groups/bioinformatics.