Characterization of species-specific repeats in 613 prokaryotic species

Prokaryotes are in general believed to possess small, compactly organized genomes, with repetitive sequences forming only a small part of them. Nonetheless, many prokaryotic genomes in fact contain species-specific repeats (>85 bp long genomic sequences with less than 60% identity to other species) as we have previously demonstrated. However, it is not known at present how frequent such species-specific repeats are and what their functional roles in bacterial genomes may be. Therefore, we have conducted a comprehensive survey of prokaryotic species-specific repeats and characterized them to examine as to whether there are functional classes among different repeats or not and how they are mutually related to each other. Of the 613 distinct prokaryotic species analyzed, 97% were found to contain at least one species-specific repeats. It seems interesting to note that the species-specific repeats thus identified appear to be functionally variable in different genomes: in some genomes, they are mostly associated with duplicated protein-coding genes, whereas in some other genomes with rRNA and tRNA genes. Contrary to what may be expected, only one-fourth of the species-specific repeats were found to be associated with mobile genetic elements.     

Screening for bacterial-fungal associations in a south-eastern US salt marsh using pre-established fungal monocultures

Both bacteria and fungi play critical roles in decomposition processes in many natural environments, yet only rarely have they been studied as an integrated community. We examined whether physical associations exist between individual bacterial and fungal species that co-occur on decaying smooth cordgrass, Spartina alterniflora, in a south-eastern US salt marsh. Fungal-pervaded decaying Spartina was used as "bait" for potential bacterial associates. The bundles (infiltrated with one of three dominant fungal members of the decomposer assemblage, or an autoclaved control) were placed in a salt marsh and collected biweekly for 6 weeks during the first experiment (late summer 2002), and weekly for 3 weeks during the second experiment (early summer 2003). Terminal-restriction fragment length polymorphism (T-RFLP) analysis of 16S rRNA genes was used to track colonization by bacterial taxa in association with the established fungal species. T-RFLP analysis of 18S-to-28S internal transcribed spacer (ITS) regions was used to monitor changes in fungal communities once bundles had been placed in the field. Results from both years were nearly identical, and showed that invasion by fungi other than the bait species was slow, resulting in a virtual fungal monoculture for several weeks into the experiments. Surprisingly, bacterial communities were unaffected by the identity of the fungal bait. Regardless of the fungal species, and even in the absence of prior fungal colonization, bacterial 16S rRNA profiles were remarkably similar. These results suggest that few species-specific associations, either positive or negative, exist between bacterial and fungal members of the Spartina decomposer community during initial colonization.     

A new experimental approach for studying bacterial genomic island evolution identifies island genes with bacterial host-specific expression patterns

Background  

Genomic islands are regions of bacterial genomes that have been acquired by horizontal transfer and often contain blocks of genes that function together for specific processes. Recently, it has become clear that the impact of genomic islands on the evolution of different bacterial species is significant and represents a major force in establishing bacterial genomic variation. However, the study of genomic island evolution has been mostly performed at the sequence level using computer software or hybridization analysis to compare different bacterial genomic sequences. We describe here a novel experimental approach to study the evolution of species-specific bacterial genomic islands that identifies island genes that have evolved in such a way that they are differentially-expressed depending on the bacterial host background into which they are transferred.     

Clues about the genetic basis of adaptation emerge from comparing the proteomes of two Ostreococcus ecotypes (Chlorophyta, Prasinophyceae)

We compared the proteomes of two picoplanktonic Ostreococcus unicellular green algal ecotypes to analyze the genetic basis of their adaptation with their ecological niches. We first investigated the function of the species-specific genes using Gene Ontology databases and similarity searches. Although most species-specific genes had no known function, we identified several species-specific functions involved in various cellular processes, which could be critical for environmental adaptations. Additionally, we investigated the rate of evolution of orthologous genes and its distribution across chromosomes. We show that faster evolving genes encode significantly more membrane or excreted proteins, consistent with the notion that selection acts on cell surface modifications that is driven by selection for resistance to viruses and grazers, keystone actors of phytoplankton evolution. The relationship between GC content and chromosome length also suggests that both strains have experienced recombination since their divergence and that lack of recombination on the two outlier chromosomes could explain part of their peculiar genomic features, including higher rates of evolution.     

A Plasmodium whole-genome synteny map: indels and synteny breakpoints as foci for species-specific genes

Whole-genome comparisons are highly informative regarding genome evolution and can reveal the conservation of genome organization and gene content, gene regulatory elements, and presence of species-specific genes. Initial comparative genome analyses of the human malaria parasite Plasmodium falciparum and rodent malaria parasites (RMPs) revealed a core set of 4,500 Plasmodium orthologs located in the highly syntenic central regions of the chromosomes that sharply defined the boundaries of the variable subtelomeric regions. We used composite RMP contigs, based on partial DNA sequences of three RMPs, to generate a whole-genome synteny map of P. falciparum and the RMPs. The core regions of the 14 chromosomes of P. falciparum and the RMPs are organized in 36 synteny blocks, representing groups of genes that have been stably inherited since these malaria species diverged, but whose relative organization has altered as a result of a predicted minimum of 15 recombination events. P. falciparum-specific genes and gene families are found in the variable subtelomeric regions (575 genes), at synteny breakpoints (42 genes), and as intrasyntenic indels (126 genes). Of the 168 non-subtelomeric P. falciparum genes, including two newly discovered gene families, 68% are predicted to be exported to the surface of the blood stage parasite or infected erythrocyte. Chromosomal rearrangements are implicated in the generation and dispersal of P. falciparum-specific gene families, including one encoding receptor-associated protein kinases. The data show that both synteny breakpoints and intrasyntenic indels can be foci for species-specific genes with a predicted role in host-parasite interactions and suggest that, besides rearrangements in the subtelomeric regions, chromosomal rearrangements may also be involved in the generation of species-specific gene families. A majority of these genes are expressed in blood stages, suggesting that the vertebrate host exerts a greater selective pressure than the mosquito vector, resulting in the acquisition of diversity.     

Codon Bias is a Major Factor Explaining Phage Evolution in Translationally Biased Hosts

The size and diversity of bacteriophage populations require methodologies to quantitatively study the landscape of phage differences. Statistical approaches are confronted with small genome sizes forbidding significant single-phage analysis, and comparative methods analyzing full phage genomes represent an alternative but they are of difficult interpretation due to lateral gene transfer, which creates a mosaic spectrum of related phage species. Based on a large-scale codon bias analysis of 116 DNA phages hosted by 11 translationally biased bacteria belonging to different phylogenetic families, we observe that phage genomes are almost always under codon selective pressure imposed by translationally biased hosts, and we propose a classification of phages with translationally biased hosts which is based on adaptation patterns. We introduce a computational method for comparing phages sharing homologous proteins, possibly accepted by different hosts. We observe that throughout phages, independently from the host, capsid genes appear to be the most affected by host translational bias. For coliphages, genes involved in virion morphogenesis, host interaction and ssDNA binding are also affected by adaptive pressure. Adaptation affects long and small phages in a significant way. We analyze in more detail the Microviridae phage space to illustrate the potentiality of the approach. The small number of directions in adaptation observed in phages grouped around ϕX174 is discussed in the light of functional bias. The adaptation analysis of the set of Microviridae phages defined around ϕMH2K shows that phage classification based on adaptation does not reflect bacterial phylogeny. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Human and simian immunodeficiency virus-mediated upregulation of the apoptotic factor TRAIL occurs in antigen-presenting cells from AIDS-susceptible but not from AIDS-resistant species

Human immunodeficiency virus (HIV) and simian immunodeficiency virus (SIV) infections lead to AIDS in humans and rhesus macaques (RM), while they are asymptomatic in species naturally infected with SIV, such as chimpanzees, sooty mangabeys (SM), and African green monkeys (AGM). Differential CD4(+) T-cell apoptosis may be responsible for these species-specific differences in susceptibility to disease. To identify factors that influence the different apoptotic responses of these species, we analyzed virus-infected human and nonhuman primate peripheral blood mononuclear cells (PBMC). We found that the apoptotic factor TRAIL was present at higher levels in human and RM PBMC cultures and was mediating, at least in part, CD4(+) T-cell apoptosis in these cultures. The species-specific increase in TRAIL and death receptor expression observed with cultures also occurred in vivo in SIV-infected RM but not in SIV-infected SM. In human and RM myeloid immature dendritic cells and macrophages, the virus-induced expression of TRAIL and other interferon-inducible genes, which did not occur in the same cells from chimpanzee, SM, and AGM, was Tat dependent. Our results link the differential induction of TRAIL in human and nonhuman primate cells to species-specific differences in disease susceptibility.     

Insertion of Horizontally Transferred Genes within Conserved Syntenic Regions of Yeast Genomes

Horizontal gene transfer has been occasionally mentioned in eukaryotic genomes, but such events appear much less numerous than in prokaryotes, where they play important functional and evolutionary roles. In yeasts, few independent cases have been described, some of which corresponding to major metabolic functions, but no systematic screening of horizontally transferred genes has been attempted so far. Taking advantage of the synteny conservation among five newly sequenced and annotated genomes of Saccharomycetaceae, we carried out a systematic search for HGT candidates amidst genes present in only one species within conserved synteny blocks. Out of 255 species-specific genes, we discovered 11 candidates for HGT, based on their similarity with bacterial proteins and on reconstructed phylogenies. This corresponds to a minimum of six transfer events because some horizontally acquired genes appear to rapidly duplicate in yeast genomes (e.g. YwqG genes in Kluyveromyces thermotolerans and serine recombinase genes of the IS607 family in Saccharomyces kluyveri). We show that the resulting copies are submitted to a strong functional selective pressure. The mechanisms of DNA transfer and integration are discussed, in relation with the generally small size of HGT candidates. Our results on a limited set of species expand by 50% the number of previously published HGT cases in hemiascomycetous yeasts, suggesting that this type of event is more frequent than usually thought. Our restrictive method does not exclude the possibility that additional HGT events exist. Actually, ancestral events common to several yeast species must have been overlooked, and the absence of homologs in present databases leaves open the question of the origin of the 244 remaining species-specific genes inserted within conserved synteny blocks.     

Genomic exploration of the hemiascomycetous yeasts: 21. Comparative functional classification of genes

We explored the biological diversity of hemiascomycetous yeasts using a set of 22000 newly identified genes in 13 species through BLASTX searches. Genes without clear homologue in Saccharomyces cerevisiae appeared to be conserved in several species, suggesting that they were recently lost by S. cerevisiae. They often identified well-known species-specific traits. Cases of gene acquisition through horizontal transfer appeared to occur very rarely if at all. All identified genes were ascribed to functional classes. Functional classes were differently represented among species. Species classification by functional clustering roughly paralleled rDNA phylogeny. Unequal distribution of rapidly evolving, ascomycete-specific, genes among species and functions was shown to contribute strongly to this clustering. A few cases of gene family amplification were documented, but no general correlation could be observed between functional differentiation of yeast species and variations of gene family sizes. Yeast biological diversity seems thus to result from limited species-specific gene losses or duplications, and for a large part from rapid evolution of genes and regulatory factors dedicated to specific functions.     

Gene gain and gene loss in streptococcus: is it driven by habitat?

Bacterial genomes can evolve either by gene gain, gene loss, mutating existing genes, and/or by duplication of existing genes. Recent studies have clearly demonstrated that the acquisition of new genes by lateral gene transfer (LGT) is a predominant force in bacterial evolution. To better understand the significance of LGT, we employed a comparative genomics approach to model species-specific and intraspecies gene insertions/deletions (ins/del among 12 sequenced streptococcal genomes using a maximum likelihood method. This study indicates that the rate of gene ins/del is higher on the external branches and varies dramatically for each species. We have analyzed here some of the experimentally characterized species-specific genes that have been acquired by LGT and conclude that at least a portion of these genes have a role in adaptation.     

Single gene mutations in Drosophila: what can they tell us about the evolution of sexual behaviour?

The molecular analysis of specific mutant genes that affect the courtship behaviours of Drosophila melanogaster males and females is discussed in the light of the possibility that they may contribute to mate choice. There is clear evidence that some genes can act as a reservoir of species-specific behaviour, particularly for the male actions during courtship. However, to date there has not been a single genetic locus that has been isolated at the molecular level and shown to be associated with a change in female preference. There are some promising avenues of exploration, in that recent genetic analyses suggest that a small number of genes may make major contributions to female preferences. Finally a candidate gene approach is advocated in which orthologous genes from other species of Drosophila are used as natural mutations, and transformed into D. melanogaster hosts to investigate whether they carry species-specific mating information of the donor.     

Detection of bacterial pathogens in environmental samples using DNA microarrays

Polymerase chain reaction (PCR) is an important tool for pathogen detection, but historically, it has not been possible to accurately identify PCR products without sequencing, Southern blots, or dot-blots. Microarrays can be coupled with PCR where they serve as a set of parallel dot-blots to enhance product detection and identification. Microarrays are composed of many discretely located probes on a solid substrate such as glass. Each probe is composed of a sequence that is complimentary to a pathogen-specific gene sequence. PCR is used to amplify one or more genes and the products are then hybridized to the array to identify species-specific polymorphism within one or more genes. We illustrate this type of array using 16S rDNA probes suitable for distinguishing between several salmonid pathogens. We also describe the use of microarrays for direct detection of either RNA or DNA without the aid of PCR, although the sensitivity of these systems currently limits their application for pathogen detection. Finally, microarrays can also be used to "fingerprint" bacterial isolates and they can be used to identify diagnostic markers suitable for developing new PCR-based detection assays. We illustrate this type of array for subtyping an important food-borne pathogen, Listeria monocytogenes.     

Seven GC-rich microbial genomes adopt similar codon usage patterns regardless of their phylogenetic lineages

Seven GC-rich (group I) and three AT-rich (group II) microbial genomes are analyzed in this paper. The seven microbes in group I belong to different phylogenetic lineages, even different domains of life. The common feature is that they are highly GC-rich organisms, with more than 60% genomic GC content. Group II includes three bacteria, which belong to the same subdivision as Pseudomonas aeruginosa in group I. The genomic GC content of the three bacteria is in the range of 26-50%. It is shown that although the phylogenetic lineages of the organisms in group I are remote, the common feature of highly genomic GC content forces them to adopt similar codon usage patterns, which constitutes the basis of an algorithm using a set of universal parameters to recognize known genes in the seven genomes. The common codon usage pattern of function known genes in the seven genomes is GGS type, where G, G, and S are the bases of G, non-G, and G/C, respectively. On the contrary, although the phylogenetic lineages of the three bacteria in group II are quite close, the codon usage patterns of function known genes in these genomes are obviously distinct. There are no universal parameters to identify known genes in the three genomes in group II. It can be deduced that the genomic GC content is more important than phylogenetic lineage in gene recognition programs. We hope that the work might be useful for understanding the common characteristics in the organization of microbial genomes.     

The Drosophila subobscura Adh genomic region contains valuable evolutionary markers.

We have sequenced 4 kb of the genomic region comprising the Adh (Alcohol dehydrogenase) gene of Drosophila subobscura. In agreement with other species which belong to the same subgenus, two structural genes, Adh and Adh-dup, are contained in this region. The main features of these two genes of D. subobscura have been inferred from the sequence data and compared with the homologous region of D. ambigua and D. pseudoobscura. Drosophila subobscura Adh and Adh-dup differ from those of D. ambigua at a corrected estimation of 10.1% and 12.5%, respectively, while from those of D. pseudoobscura they differ by 9.5% and 8.1%, respectively. Our data suggest that Adh and Adh-dup are evolving independently, showing a species-specific pattern. Moreover, particular features of some regions of these genes make them valuable evolutionary hallmarks. For instance, replacement substitutions in the third exon of Adh may indicate the branching of the melanogaster-obscura groups, whereas replacement substitutions in the third exon of the Adh-dup could be used to assess speciation within the obscura group.