Comparative analysis of transposable elements in the melanogaster subgroup sequenced genomes

Transposable elements (TEs) are indwelling components of genomes, and their dynamics have been a driving force in genome evolution. Although we now have more information concerning their amounts and characteristics in various organisms, we still have little data from overall comparisons of their sequences in very closely-related species. While the Drosophila melanogaster genome has been extensively studied, we have only limited knowledge regarding the precise TE sequences in the genomes of the related species Drosophila simulans, Drosophila sechellia and Drosophila yakuba. In this study we analyzed the number and structure of TE copies in the sequenced genomes of these four species. Our findings show that, unexpectedly, the number of TE insertions in D. simulans is greater than that in D. melanogaster, but that most of the copies in D. simulans are degraded and in small fragments, as in D. sechellia and D. yakuba. This suggests that all three species were invaded by numerous TEs a long time ago, but have since regulated their activity, as the present TE copies are degraded, with very few full-length elements. In contrast, in D. melanogaster, a recent activation of TEs has resulted in a large number of almost-identical TE copies. We have detected variants of some TEs in D. simulans and D. sechellia, that are almost identical to the reference TE sequences in D. melanogaster, suggesting that D. melanogaster has recently been invaded by active TE variants from the other species. Our results indicate that the three species D. simulans, D. sechellia, and D. yakuba seem to be at a different stage of their TE life cycle when compared to D. melanogaster. Moreover, we show that D. melanogaster has been invaded by active TE variants for several TE families likely to come from D. simulans or the ancestor of D. simulans and D. sechellia. The numerous horizontal transfer events implied to explain these results could indicate introgression events between these species.     

Comparative analysis of Acinetobacters: three genomes for three lifestyles

Acinetobacter baumannii is the source of numerous nosocomial infections in humans and therefore deserves close attention as multidrug or even pandrug resistant strains are increasingly being identified worldwide. Here we report the comparison of two newly sequenced genomes of A. baumannii. The human isolate A. baumannii AYE is multidrug resistant whereas strain SDF, which was isolated from body lice, is antibiotic susceptible. As reference for comparison in this analysis, the genome of the soil-living bacterium A. baylyi strain ADP1 was used. The most interesting dissimilarities we observed were that i) whereas strain AYE and A. baylyi genomes harbored very few Insertion Sequence elements which could promote expression of downstream genes, strain SDF sequence contains several hundred of them that have played a crucial role in its genome reduction (gene disruptions and simple DNA loss); ii) strain SDF has low catabolic capacities compared to strain AYE. Interestingly, the latter has even higher catabolic capacities than A. baylyi which has already been reported as a very nutritionally versatile organism. This metabolic performance could explain the persistence of A. baumannii nosocomial strains in environments where nutrients are scarce; iii) several processes known to play a key role during host infection (biofilm formation, iron uptake, quorum sensing, virulence factors) were either different or absent, the best example of which is iron uptake. Indeed, strain AYE and A. baylyi use siderophore-based systems to scavenge iron from the environment whereas strain SDF uses an alternate system similar to the Haem Acquisition System (HAS). Taken together, all these observations suggest that the genome contents of the 3 Acinetobacters compared are partly shaped by life in distinct ecological niches: human (and more largely hospital environment), louse, soil.     

A domain-centric analysis of oomycete plant pathogen genomes reveals unique protein organization

Oomycetes comprise a diverse group of organisms that morphologically resemble fungi but belong to the stramenopile lineage within the supergroup of chromalveolates. Recent studies have shown that plant pathogenic oomycetes have expanded gene families that are possibly linked to their pathogenic lifestyle. We analyzed the protein domain organization of 67 eukaryotic species including four oomycete and five fungal plant pathogens. We detected 246 expanded domains in fungal and oomycete plant pathogens. The analysis of genes differentially expressed during infection revealed a significant enrichment of genes encoding expanded domains as well as signal peptides linking a substantial part of these genes to pathogenicity. Overrepresentation and clustering of domain abundance profiles revealed domains that might have important roles in host-pathogen interactions but, as yet, have not been linked to pathogenicity. The number of distinct domain combinations (bigrams) in oomycetes was significantly higher than in fungi. We identified 773 oomycete-specific bigrams, with the majority composed of domains common to eukaryotes. The analyses enabled us to link domain content to biological processes such as host-pathogen interaction, nutrient uptake, or suppression and elicitation of plant immune responses. Taken together, this study represents a comprehensive overview of the domain repertoire of fungal and oomycete plant pathogens and points to novel features like domain expansion and species-specific bigram types that could, at least partially, explain why oomycetes are such remarkable plant pathogens.     

Comparative analyses of three legume species reveals conserved and unique root extracellular proteins

Increasing evidence suggests that root extracellular proteins are involved in interactions between roots and their soil environment. In the present study, exudates released by 6‐day‐old roots of the three legume species white lupin (Lupinus albus), soybean (Glycine max), and cowpea (Vigna sinensis) were collected under axenic conditions, and their constitutively secreted proteomes were analyzed. Between 42 and 93 unique root extracellular proteins with 2 or more different peptide fragments per protein were identified by LC‐MS/MS. Functional annotation of these proteins classified them into 14–16 different functional categories. Among those 14 homologous proteins were identified in at least two legume species. Among the unique proteins, 58 in white lupin, 85 in soybean, and 31 in cowpea were specific for each plant species, and many of them were classified in the same functional categories. Interestingly, in contrast to soybean and cowpea, two protein bands of approximately 16 and 30 kDa were present on the SDS‐PAGE gel of white lupin. The identification of these bands revealed a class III chitinase and a thaumatin‐like protein. Both belong to the class of pathogenesis‐related proteins. The results imply that root extracellular proteins play important roles in the cross‐talk between plant roots and the rhizosphere.     

Comparative genomics of 274 Vibrio cholerae genomes reveals mobile functions structuring three niche dimensions

Background

Vibrio cholerae is a globally dispersed pathogen that has evolved with humans for centuries, but also includes non-pathogenic environmental strains. Here, we identify the genomic variability underlying this remarkable persistence across the three major niche dimensions space, time, and habitat.

Results

Taking an innovative approach of genome-wide association applicable to microbial genomes (GWAS-M), we classify 274 complete V. cholerae genomes by niche, including 39 newly sequenced for this study with the Ion Torrent DNA-sequencing platform. Niche metadata were collected for each strain and analyzed together with comprehensive annotations of genetic and genomic attributes, including point mutations (single-nucleotide polymorphisms, SNPs), protein families, functions and prophages.

Conclusions

Our analysis revealed that genomic variations, in particular mobile functions including phages, prophages, transposable elements, and plasmids underlie the metadata structuring in each of the three niche dimensions. This underscores the role of phages and mobile elements as the most rapidly evolving elements in bacterial genomes, creating local endemicity (space), leading to temporal divergence (time), and allowing the invasion of new habitats. Together, we take a data-driven approach for comparative functional genomics that exploits high-volume genome sequencing and annotation, in conjunction with novel statistical and machine learning analyses to identify connections between genotype and phenotype on a genome-wide scale.

Electronic supplementary material

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

Comparative analysis of the small heat shock proteins in three angiosperm genomes identifies new subfamilies and reveals diverse evolutionary patterns

The small heat shock proteins (sHSPs) are a diverse family of molecular chaperones. It is well established that these proteins are crucial components of the plant heat shock response. They also have important roles in other stress responses and in normal development. We have conducted a comparative sequence analysis of the sHSPs in three complete angiosperms genomes: Arabidopsis thaliana, Populus trichocarpa, and Oryza sativa. Our phylogenetic analysis has identified four additional plant sHSP subfamilies and thus has increased the number of plant sHSP subfamilies from 7 to 11. We have also identified a number of novel sHSP genes in each genome that lack close homologs in other genomes. Using publicly available gene expression data and predicted secondary structures, we have determined that the sHSPs in plants are far more diverse in sequence, expression profile, and in structure than had been previously known. Some of the newly identified subfamilies are not stress regulated, may not posses the highly conserved large oligomer structure, and may not even function as molecular chaperones. We found no consistent evolutionary patterns across the three species studied. For example, gene conversion was found among the sHSPs in O. sativa but not in A. thaliana or P. trichocarpa. Among the three species, P. trichocarpa had the most sHSPs. This was due to an expansion of the cytosolic I sHSPs that was not seen in the other two species. Our analysis indicates that the sHSPs are a dynamic protein family in angiosperms with unexpected levels of diversity. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comparative genome analysis of Bacillus cereus group genomes with Bacillus subtilis

Genome features of the Bacillus cereus group genomes (representative strains of Bacillus cereus, Bacillus anthracis and Bacillus thuringiensis sub spp. israelensis) were analyzed and compared with the Bacillus subtilis genome. A core set of 1381 protein families among the four Bacillus genomes, with an additional set of 933 families common to the B. cereus group, was identified. Differences in signal transduction pathways, membrane transporters, cell surface structures, cell wall, and S-layer proteins suggesting differences in their phenotype were identified. The B. cereus group has signal transduction systems including a tyrosine kinase related to two-component system histidine kinases from B. subtilis. A model for regulation of the stress responsive sigma factor sigmaB in the B. cereus group different from the well studied regulation in B. subtilis has been proposed. Despite a high degree of chromosomal synteny among these genomes, significant differences in cell wall and spore coat proteins that contribute to the survival and adaptation in specific hosts has been identified.     

Comparative study reveals unique features of the mycobiota in peat soils samples from Japan and Scotland

We isolated filamentous fungi from soil samples of peat layers in Aomori and Oita Prefectures in Japan and Perth and Kinross district in Scotland by a serial dilution plate technique. The mycobiota in each peat soil showed some common and characteristic features. The abundance of fungal isolates (CFU/g) from peat soil was low: about 1/3 to 1/30 compared with evergreen or coniferous forests or cultivated soil. Trichoderma or Mucorales species were scarcely observed; these fungi occupied only 3% of the total number of colonies. On the other hand, fungi such as Conioscypha and Tolypocladium that are normally isolated rather rarely were encountered at a comparatively high rate. Acremonium guillematii and Tolypocladium cylindrosporum were recorded for the first time in Japan. Sterile fungi occupied 50% of the total number of isolates. The low abundance of fast-growing fungi enabled us to pick slow-growing fungi up easily from the isolation medium. It is interesting that species not previously described in Japan, or scarcely reported, were isolated commonly from both Japanese and Scottish samples. A peat soil sample is therefore an attractive source of untapped microbial resources.     

Comparative analysis of 22 coronavirus HKU1 genomes reveals a novel genotype and evidence of natural recombination in coronavirus HKU1

We sequenced and compared the complete genomes of 22 strains of coronavirus HKU1 (CoV HKU1) obtained from nasopharyngeal aspirates of patients with respiratory tract infections over a 2-year period. Phylogenetic analysis of 24 putative proteins and polypeptides showed that the 22 CoV HKU1 strains fell into three clusters (genotype A, 13 strains; genotype B, 3 strains and genotype C, 6 strains). However, different phylogenetic relationships among the three clusters were observed in different regions of their genomes. From nsp4 to nsp6, the genotype A strains were clustered with the genotype B strains. For nsp7 and nsp8 and from nsp10 to nsp16, the genotype A strains were clustered with the genotype C strains. From hemagglutinin esterase (HE) to nucleocapsid (N), the genotype B strains were clustered closely with the genotype C strains. Bootscan analysis showed possible recombination between genotypes B and C from nucleotide positions 11,500 to 13,000, corresponding to the nsp6-nsp7 junction, giving rise to genotype A, and between genotypes A and B from nucleotide positions 21,500 to 22,500, corresponding to the nsp16-HE junction, giving rise to genotype C. Multiple alignments further narrowed the sites of crossover to a 143-bp region between nucleotide positions 11,750 and 11,892 and a 29-bp region between nucleotide positions 21,502 and 21,530. Genome analysis also revealed various numbers of tandem copies of a perfect 30-base acidic tandem repeat (ATR) which encodes NDDEDVVTGD and various numbers and sequences of imperfect repeats in the N terminus of nsp3 inside the acidic domain upstream of papain-like protease 1 among the 22 genomes. All 10 CoV HKU1 strains with incomplete imperfect repeats (1.4 and 4.4) belonged to genotype A. The present study represents the first evidence for natural recombination in coronavirus associated with human infection. Analysis of a single gene is not sufficient for the genotyping of CoV HKU1 strains but requires amplification and sequencing of at least two gene loci, one from nsp10 to nsp16 (e.g., pol or helicase) and another from HE to N (e.g., spike or N). Further studies will delineate whether the ATR is useful for the molecular typing of CoV HKU1.     

Comparative analysis reveals distinctive epigenetic features of the human cerebellum

Identifying the molecular underpinnings of the neural specializations that underlie human cognitive and behavioral traits has long been of considerable interest. Much research on human-specific changes in gene expression and epigenetic marks has focused on the prefrontal cortex, a brain structure distinguished by its role in executive functions. The cerebellum shows expansion in great apes and is gaining increasing attention for its role in motor skills and cognitive processing, including language. However, relatively few molecular studies of the cerebellum in a comparative evolutionary context have been conducted. Here, we identify human-specific methylation in the lateral cerebellum relative to the dorsolateral prefrontal cortex, in a comparative study with chimpanzees (Pan troglodytes) and rhesus macaques (Macaca mulatta). Specifically, we profiled genome-wide methylation levels in the three species for each of the two brain structures and identified human-specific differentially methylated genomic regions unique to each structure. We further identified which differentially methylated regions (DMRs) overlap likely regulatory elements and determined whether associated genes show corresponding species differences in gene expression. We found greater human-specific methylation in the cerebellum than the dorsolateral prefrontal cortex, with differentially methylated regions overlapping genes involved in several conditions or processes relevant to human neurobiology, including synaptic plasticity, lipid metabolism, neuroinflammation and neurodegeneration, and neurodevelopment, including developmental disorders. Moreover, our results show some overlap with those of previous studies focused on the neocortex, indicating that such results may be common to multiple brain structures. These findings further our understanding of the cerebellum in human brain evolution.     

Comparative analysis of human and bovine protein kinases reveals unique relationship and functional diversity

Reversible protein phosphorylation by protein kinases and phosphatases is a common event in various cellular processes. The eukaryotic protein kinase superfamily, which is one of the largest superfamilies of eukaryotic proteins, plays several roles in cell signaling and diseases. We identified 482 eukaryotic protein kinases and 39 atypical protein kinases in the bovine genome, by searching publicly accessible genetic-sequence databases. Bovines have 512 putative protein kinases, each orthologous to a human kinase. Whereas orthologous kinase pairs are, on an average, 90.6% identical, orthologous kinase catalytic domain pairs are, on an average, 95.9% identical at the amino acid level. This bioinformatic study of bovine protein kinases provides a suitable framework for further characterization of their functional and structural properties.     

Comparative multivariate analysis of codon and amino acid usage in three Leishmania genomes

Multivariate analysis of codon and amino acid usage was performed for three Leishmania species, including L. donovani, L. infantum and L. major. It was revealed that all three species are under mutational bias and translational selection. Lower GC 12 and higher GC 3S in all three parasites suggests that the ancestral highly expressed genes (HEGs), compared to lowly expressed genes (LEGs), might have been rich in AT-content. This also suggests that there must have been a faster rate of evolution under GC-bias in LEGs. It was observed from the estimation of synonymous/non-synonymous substitutions in HEGs that the HEG dataset of L. donovani is much closer to L. major evolutionarily. This is also supported by the higher d N value as compared to d S between L. donovani and L. major, suggesting the conservation of synonymous codon positions between these two species and the role of translational selection in shaping the composition of protein-coding genes.     

N49 phospholipase A2, a unique subgroup of snake venom group II phospholipase A2

A novel phospholipase A2 (PLA2) with Asn at its site 49 was purified from the snake venom of Protobothrops mucrosquamatus by using SP-Sephadex C25, Superdex 75, Heparin-Sepharose (FF) and HPLC reverse-phage C18 chromatography and designated as TM-N49. It showed a molecular mass of 13.875 kDa on MALDI-TOF. TM-N49 does not possess enzymatic, hemolytic and hemorrhagic activities. It fails to induce platelet aggregation by itself, and does not inhibit the platelet aggregation induced by ADP. However, it exhibits potent myotoxic activity causing inflammatory cell infiltration, severe myoedema, myonecrosis and myolysis in the gastrocnemius muscles of BALB/c mice. Phylogenetic analysis found that that TM-N49 combined with two phospholipase A2s from Trimeresurus stejnegeri, TsR6 and CTs-R6 cluster into one group. Structural and functional analysis indicated that these phospholipase A2s are distinct from the other subgroups (D49 PLA2, S49 PLA2 and K49 PLA2) and represent a unique subgroup of snake venom group II PLA2, named N49 PLA2 subgroup.     

Comparative profiling of Pseudomonas aeruginosa strains reveals differential expression of novel unique and conserved small RNAs

Pseudomonas aeruginosa is a highly adaptable bacterium that thrives in a broad range of ecological niches and can infect multiple hosts as diverse as plants, nematodes and mammals. In humans, it is an important opportunistic pathogen. This wide adaptability correlates with its broad genetic diversity. In this study, we used a deep-sequencing approach to explore the complement of small RNAs (sRNAs) in P. aeruginosa as the number of such regulatory molecules previously identified in this organism is relatively low, considering its genome size, phenotypic diversity and adaptability. We have performed a comparative analysis of PAO1 and PA14 strains which share the same host range but differ in virulence, PA14 being considerably more virulent in several model organisms. Altogether, we have identified more than 150 novel candidate sRNAs and validated a third of them by Northern blotting. Interestingly, a number of these novel sRNAs are strain-specific or showed strain-specific expression, strongly suggesting that they could be involved in determining specific phenotypic traits.     

Comparative computational analysis of SARS-CoV-2 nucleocapsid protein epitopes in taxonomically related coronaviruses

Several research lines are currently ongoing to address the multitude of facets of the pandemic COVID-19. In line with the One-Health concept, extending the target of the studies to the animals which humans are continuously interacting with may favor a better understanding of the SARS-CoV-2 biology and pathogenetic mechanisms; thus, helping to adopt the most suitable containment measures. The last two decades have already faced severe manifestations of the coronavirus infection in both humans and animals, thus, circulating epitopes from previous outbreaks might confer partial protection from SARS-CoV-2 infections. In the present study, we provide an in-silico survey of the major nucleocapsid protein epitopes and compare them with the homologues of taxonomically-related coronaviruses with tropism for animal species that are closely inter-related with the human beings population all over the world. Protein sequence alignment provides evidence of high sequence homology for some of the investigated proteins. Moreover, structural epitope mapping by homology modelling revealed a potential immunogenic value also for specific sequences scoring a lower identity with SARS-CoV-2 nucleocapsid proteins. These evidence provide a molecular structural rationale for a potential role in conferring protection from SARS-CoV-2 infection and identifying potential candidates for the development of diagnostic tools and prophylactic-oriented strategies.