Selection-driven gene loss in bacteria

Gene loss by deletion is a common evolutionary process in bacteria, as exemplified by bacteria with small genomes that have evolved from bacteria with larger genomes by reductive processes. The driving force(s) for genome reduction remains unclear, and here we examined the hypothesis that gene loss is selected because carriage of superfluous genes confers a fitness cost to the bacterium. In the bacterium Salmonella enterica, we measured deletion rates at 11 chromosomal positions and the fitness effects of several spontaneous deletions. Deletion rates varied over 200-fold between different regions with the replication terminus region showing the highest rates. Approximately 25% of the examined deletions caused an increase in fitness under one or several growth conditions, and after serial passage of wild-type bacteria in rich medium for 1,000 generations we observed fixation of deletions that substantially increased bacterial fitness when reconstructed in a non-evolved bacterium. These results suggest that selection could be a significant driver of gene loss and reductive genome evolution.     

Evolution of ABC transporters by gene duplication and their role in human disease

The ATP-binding cassette (ABC) transporter genes represent the largest family of transporters and these genes are abundant in the genome of all vertebrates. Through analysis of the genome sequence databases we have characterized the full complement of ABC genes from several mammals and other vertebrates. Multiple gene duplication and deletion events were identified in ABC genes in different lineages indicating that the process of gene evolution is still ongoing. Gene duplication resulting in either gene birth or gene death plays a major role in the evolution of the vertebrate ABC genes. The understanding of this mechanism is important in the context of human health because these ABC genes are associated with human disease, involving nearly all organ systems of the body. In addition, ABC genes play an important role in the development of drug resistance in cancer cells. Future genetic, functional, and evolutionary studies of ABC transporters will provide important insight into human and animal biology.     

Metabolic networks of Sodalis glossinidius: a systems biology approach to reductive evolution

Background

Genome reduction is a common evolutionary process affecting bacterial lineages that establish symbiotic or pathogenic associations with eukaryotic hosts. Such associations yield highly reduced genomes with greatly streamlined metabolic abilities shaped by the type of ecological association with the host. Sodalis glossinidius, the secondary endosymbiont of tsetse flies, represents one of the few complete genomes available of a bacterium at the initial stages of this process. In the present study, genome reduction is studied from a systems biology perspective through the reconstruction and functional analysis of genome-scale metabolic networks of S. glossinidius.

Results

The functional profile of ancestral and extant metabolic networks sheds light on the evolutionary events underlying transition to a host-dependent lifestyle. Meanwhile, reductive evolution simulations on the extant metabolic network can predict possible future evolution of S. glossinidius in the context of genome reduction. Finally, knockout simulations in different metabolic systems reveal a gradual decrease in network robustness to different mutational events for bacterial endosymbionts at different stages of the symbiotic association.

Conclusions

Stoichiometric analysis reveals few gene inactivation events whose effects on the functionality of S. glossinidius metabolic systems are drastic enough to account for the ecological transition from a free-living to host-dependent lifestyle. The decrease in network robustness across different metabolic systems may be associated with the progressive integration in the more stable environment provided by the insect host. Finally, reductive evolution simulations reveal the strong influence that external conditions exert on the evolvability of metabolic systems.     

Reconstructing evolution: gene transfer from plastids to the nucleus

During evolution, the genomes of eukaryotic cells have undergone major restructuring to meet the new regulatory challenges associated with compartmentalization of the genetic material in the nucleus and the organelles acquired by endosymbiosis (mitochondria and plastids). Restructuring involved the loss of dispensable or redundant genes and the massive translocation of genes from the ancestral organelles to the nucleus. Genomics and bioinformatic data suggest that the process of DNA transfer from organelles to the nucleus still continues, providing raw material for evolutionary tinkering in the nuclear genome. Recent reconstruction of these events in the laboratory has provided a unique tool to observe genome evolution in real time and to study the molecular mechanisms by which plastid genes are converted into functional nuclear genes. Here, we summarize current knowledge about plastid-to-nuclear gene transfer in the context of genome evolution and discuss new insights gained from experiments that recapitulate endosymbiotic gene transfer in the laboratory.     

Phylogenetic analyses of two "archaeal" genes in thermotoga maritima reveal multiple transfers between archaea and bacteria

The genome sequence of Thermotoga maritima revealed that 24% of its open reading frames (ORFs) showed the highest similarity scores to archaeal genes in BLAST analyses. Here we screened 16 strains from the genus Thermotoga and other related Thermotogales for the occurrence of two of these "archaeal" genes: the gene encoding the large subunit of glutamate synthase (gltB) and the myo-inositol 1P synthase gene (ino1). Both genes were restricted to the Thermotoga species within the Thermotogales. The distribution of the two genes, along with results from phylogenetic analyses, showed that they were acquired from Archaea during the divergence of the Thermotogales. Database searches revealed that three other bacteria-Dehalococcoides ethenogenes, Sinorhizobium meliloti, and Clostridium difficile-possess archaeal-type gltBs, and the phylogenetic analyses confirmed at least two lateral gene transfer (LGT) events between Bacteria and Archaea. These LGT events were also strongly supported by gene structure data, as the three domains in bacterial-type gltB are homologous to three independent ORFs in Archaea and Bacteria with archaeal-type gltBs. The ino1 gene has a scattered distribution among Bacteria, and apart from the Thermotoga strains it is found only in Aquifex aeolicus, D. ethenogenes, and some high-G+C Gram-positive bacteria. Phylogenetic analysis of the ino1 sequences revealed three highly supported prokaryotic clades, all containing a mixture of archaeal and bacterial sequences, and suggested that all bacterial ino1 genes had been recruited from archaeal donors. The Thermotoga strains and A. aeolicus acquired this gene independently from different archaeal species. Although transfer of genes from hyperthermophilic Archaea may have facilitated the evolution of bacterial hyperthermophily, between-domain transfers also affect mesophilic species. For hyperthermophiles, we hypothesize that LGT may be as much a consequence as the cause of adaptation to hyperthermophily.     

Horizontal gene transfer regulation in bacteria as a "spandrel" of DNA repair mechanisms

Horizontal gene transfer (HGT) is recognized as the major force for bacterial genome evolution. Yet, numerous questions remain about the transferred genes, their function, quantity and frequency. The extent to which genetic transformation by exogenous DNA has occurred over evolutionary time was initially addressed by an in silico approach using the complete genome sequence of the Ralstonia solanacearum GMI1000 strain. Methods based on phylogenetic reconstruction of prokaryote homologous genes families detected 151 genes (13.3%) of foreign origin in the R. solanacearum genome and tentatively identified their bacterial origin. These putative transfers were analyzed in comparison to experimental transformation tests involving 18 different genomic DNA positions in the genome as sites for homologous or homeologous recombination. Significant transformation frequency differences were observed among these positions tested regardless of the overall genomic divergence of the R. solanacearum strains tested as recipients. The genomic positions containing the putative exogenous DNA were not systematically transformed at the highest frequencies. The two genomic "hot spots", which contain recA and mutS genes, exhibited transformation frequencies from 2 to more than 4 orders of magnitude higher than positions associated with other genes depending on the recipient strain. These results support the notion that the bacterial cell is equipped with active mechanisms to modulate acquisition of new DNA in different genomic positions. Bio-informatics study correlated recombination "hot-spots" to the presence of Chi-like signature sequences with which recombination might be preferentially initiated. The fundamental role of HGT is certainly not limited to the critical impact that the very rare foreign genes acquired mainly by chance can have on the bacterial adaptation potential. The frequency to which HGT with homologous and homeologous DNA happens in the environment might have led the bacteria to hijack DNA repair mechanisms in order to generate genetic diversity without losing too much genomic stability.     

Genome Erosion in a Nitrogen-Fixing Vertically Transmitted Endosymbiotic Multicellular Cyanobacterium

Background

An ancient cyanobacterial incorporation into a eukaryotic organism led to the evolution of plastids (chloroplasts) and subsequently to the origin of the plant kingdom. The underlying mechanism and the identities of the partners in this monophyletic event remain elusive.

Methodology/Principal Findings

To shed light on this evolutionary process, we sequenced the genome of a cyanobacterium residing extracellularly in an endosymbiosis with a plant, the water-fern Azolla filiculoides Lam. This symbiosis was selected as it has characters which make it unique among extant cyanobacterial plant symbioses: the cyanobacterium lacks autonomous growth and is vertically transmitted between plant generations. Our results reveal features of evolutionary significance. The genome is in an eroding state, evidenced by a large proportion of pseudogenes (31.2%) and a high frequency of transposable elements (∼600) scattered throughout the genome. Pseudogenization is found in genes such as the replication initiator dnaA and DNA repair genes, considered essential to free-living cyanobacteria. For some functional categories of genes pseudogenes are more prevalent than functional genes. Loss of function is apparent even within the ‘core’ gene categories of bacteria, such as genes involved in glycolysis and nutrient uptake. In contrast, serving as a critical source of nitrogen for the host, genes related to metabolic processes such as cell differentiation and nitrogen-fixation are well preserved.

Conclusions/Significance

This is the first finding of genome degradation in a plant symbiont and phenotypically complex cyanobacterium and one of only a few extracellular endosymbionts described showing signs of reductive genome evolution. Our findings suggest an ongoing selective streamlining of this cyanobacterial genome which has resulted in an organism devoted to nitrogen fixation and devoid of autonomous growth. The cyanobacterial symbiont of Azolla can thus be considered at the initial phase of a transition from free-living organism to a nitrogen-fixing plant entity, a transition process which may mimic what drove the evolution of chloroplasts from a cyanobacterial ancestor.     

Massive presence of insertion sequences in the genome of SOPE, the primary endosymbiont of the rice weevil Sitophilus oryzae

Bacteria that establish an obligate intracellular relationship with eukaryotic hosts undergo an evolutionary genomic reductive process. Recent studies have shown an increase in the number of mobile elements in the first stage of the adaptive process towards intracellular life, although these elements are absent in ancient endosymbionts. Here, the genome of SOPE, the obligate mutualistic endosymbiont of rice weevils, was used as a model to analyze the initial events that occur after symbiotic integration. During the first phases of the SOPE genome project, four different types of insertion sequence (IS) elements, belonging to well-characterized IS families from gamma-proteobacteria, were identified. In the present study, these elements, which may represent more than 20% of the complete genome, were completely characterized; their relevance as a source of gene inactivation, chromosomal rearrangements, and as participants in the genome reductive process are discussed herein.     

Introns in,introns out in plant gene families: a genomic approach of the dynamics of gene structure

Gene duplication is considered to be a source of genetic information for the creation of new functions. The Arabidopsis thaliana genome sequence revealed that a majority of plant genes belong to gene families. Regarding the problem of genes involved in the genesis of novel organs or functions during evolution, the reconstitution of the evolutionary history of gene families is of critical importance. A comparison of the intron/exon gene structure may provide clues for the understanding of the evolutionary mechanisms underlying the genesis of gene families. An extensive study of A. thaliana genome showed that families of duplicated genes may be organized according to the number and/or density of intron and the diversity in gene structure. In this paper, we propose a genomic classification of several A. thaliana gene families based on introns in an evolutionary perspective. abbreviations BGAL, -galactosidases; PCMP, plant combinatorial and modular protein     

Phylogenomic analysis reveals dynamic evolutionary history of the Drosophila heterochromatin protein 1 (HP1) gene family

Heterochromatin is the gene-poor, satellite-rich eukaryotic genome compartment that supports many essential cellular processes. The functional diversity of proteins that bind and often epigenetically define heterochromatic DNA sequence reflects the diverse functions supported by this enigmatic genome compartment. Moreover, heterogeneous signatures of selection at chromosomal proteins often mirror the heterogeneity of evolutionary forces that act on heterochromatic DNA. To identify new such surrogates for dissecting heterochromatin function and evolution, we conducted a comprehensive phylogenomic analysis of the Heterochromatin Protein 1 gene family across 40 million years of Drosophila evolution. Our study expands this gene family from 5 genes to at least 26 genes, including several uncharacterized genes in Drosophila melanogaster. The 21 newly defined HP1s introduce unprecedented structural diversity, lineage-restriction, and germline-biased expression patterns into the HP1 family. We find little evidence of positive selection at these HP1 genes in both population genetic and molecular evolution analyses. Instead, we find that dynamic evolution occurs via prolific gene gains and losses. Despite this dynamic gene turnover, the number of HP1 genes is relatively constant across species. We propose that karyotype evolution drives at least some HP1 gene turnover. For example, the loss of the male germline-restricted HP1E in the obscura group coincides with one episode of dramatic karyotypic evolution, including the gain of a neo-Y in this lineage. This expanded compendium of ovary- and testis-restricted HP1 genes revealed by our study, together with correlated gain/loss dynamics and chromosome fission/fusion events, will guide functional analyses of novel roles supported by germline chromatin.     

Synonymous codon usage and gene function are strongly related in Oryza sativa

The relationship between codon usage and gene function was investigated while considering a dataset of 2106 nuclear genes of Oryza sativa. The results of standard chi(2) test and F-statistic showed that for every 59 synonymous codons, a strongly significant association with gene functional categories existed in rice, indicating that codon usage was generally coordinated with gene function whether it was at the level of individual amino acids or at the level of nucleotides. However, it could not be directly said that the use of every codons differed significantly between any two functional categories. Notably, there existed large difference both in selection for biased codons or selection intensity among functional categories. Therefore, we identified at least two classes of genes: one group of genes, mainly belonging to the "METABOLISM" category, was tended to use G- and/or C-ending codons while the other was more biased to choose codons ending with A and/or U. The latter group contained genes of various functions, especially those genes classified into the "Nuclear Structure" category. These observations will be more important for molecular genetic engineering and genome functional annotation.     

Indications for acquisition of reductive dehalogenase genes through horizontal gene transfer by Dehalococcoides ethenogenes strain 195

The genome of Dehalococcoides ethenogenes strain 195, an anaerobic dehalorespiring bacterium, contains 18 copies of putative reductive dehalogenase genes, including the well-characterized tceA gene, whose gene product functions as the key enzyme in the environmentally important dehalorespiration process. The genome of D. ethenogenes was analyzed using a bioinformatic tool based on the frequency of oligonucleotides. The results in the form of a genomic signature revealed several local disruptions of the host signature along the genome sequence. These fractures represent DNA segments of potentially foreign origin, so-called atypical regions, which may have been acquired by an ancestor through horizontal gene transfer. Most interestingly, 15 of the 18 reductive dehalogenase genes, including the tceA gene, were found to be located in these regions, strongly indicating the foreign nature of the dehalorespiration activity. The GC content and the presence of recombinase genes within some of these regions corroborate this hypothesis. A hierarchical classification of the atypical regions containing the reductive dehalogenase genes indicated that these regions were probably acquired by several gene transfer events.     

Opsin gene duplication and divergence in ray-finned fish

Opsin gene sequences were first reported in the 1980s. The goal of that research was to test the hypothesis that human opsins were members of a single gene family and that variation in human color vision was mediated by mutations in these genes. While the new data supported both hypotheses, the greatest contribution of this work was, arguably, that it provided the data necessary for PCR-based surveys in a diversity of other species. Such studies, and recent whole genome sequencing projects, have uncovered exceptionally large opsin gene repertoires in ray-finned fishes (taxon, Actinopterygii). Guppies and zebrafish, for example, have 10 visual opsin genes each. Here we review the duplication and divergence events that have generated these gene collections. Phylogenetic analyses revealed that large opsin gene repertories in fish have been generated by gene duplication and divergence events that span the age of the ray-finned fishes. Data from whole genome sequencing projects and from large-insert clones show that tandem duplication is the primary mode of opsin gene family expansion in fishes. In some instances gene conversion between tandem duplicates has obscured evolutionary relationships among genes and generated unique key-site haplotypes. We mapped amino acid substitutions at so-called key-sites onto phylogenies and this exposed many examples of convergence. We found that dN/dS values were higher on the branches of our trees that followed gene duplication than on branches that followed speciation events, suggesting that duplication relaxes constraints on opsin sequence evolution. Though the focus of the review is opsin sequence evolution, we also note that there are few clear connections between opsin gene repertoires and variation in spectral environment, morphological traits, or life history traits.     

Global gene expression profile of Orientia tsutsugamushi

Orientia tsutsugamushi, an obligate intracellular bacterium, is the causative agent of Scrub typhus. The control mechanisms for bacterial gene expression are largely unknown. Here, the global gene expression of O. tsutsugamushi within eukaryotic cells was examined using a microarray and proteomic approaches for the first time. These approaches identified 643 genes, corresponding to approximately 30% of the genes encoded in the genome. The majority of expressed genes belonged to several functional categories including protein translation, protein processing/secretion, and replication/repair. We also searched the conserved sequence blocks (CSBs) in the O. tsutsugamushi genome which is unique in that up to 40% of its genome consists of dispersed repeated sequences. Although extensive shuffling of genomic sequences was observed between two different strains, 204 CSBs, covering 48% of the genome, were identified. When combining the data of CSBs and global gene expression, the CSBs correlates well with the location of expressed genes, suggesting the functional conservation between gene expression and genomic location. Finally, we compared the gene expression of the bacteria‐infected fibroblasts and macrophages using microarray analysis. Some major changes were the downregulation of genes involved in translation, protein processing and secretion, which correlated with the reduction in bacterial translation rates and growth within macrophages.     

Characterization of genes with tissue-specific differential expression patterns in Populus

Like many plants, Populus has an evolutionary history in which several, both recent and more ancient, genome duplication events have occurred and, therefore, constitutes an excellent model system for studying the functional evolution of genes. In the present study, we have focused on the properties of genes with tissue-specific differential expression patterns in poplar. We identified the genes by analyzing digital expression profiles derived by mapping 90,000+ expressed sequence tags (ESTs) from 18 sources to the predicted genes of Populus. Our sequence analysis suggests that tissue-specific differentially expressed genes have less diverged paralogs than average, indicating that gene duplication events is an important event in the pathway leading to this type of expression pattern. The functional analysis showed that genes coding for proteins involved in processes of functional importance for the specific tissue(s) in which they are expressed and genes coding for regulatory or responsive proteins are most common among the differentially expressed genes, demonstrating that the expression differentiation process is under strong selective pressure. Thus, our data supports a model where gene duplication followed by gene specialization or expansion of the regulatory and responsive networks leads to tissue-specific differential expression patterns. We have also searched for clustering of genes with similar expression pattern into gene-expression neighborhoods within the Populus genome. However, we could not detect any major clustering among the analyzed genes with highly specific expression patterns. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Correlations Between Bacterial Ecology and Mobile DNA

Several factors can affect the density of mobile DNA in bacterial genomes including rates of exposure to novel gene pools, recombination, and reductive evolution. These traits are difficult to measure across a broad range of bacterial species, but the ecological niches occupied by an organism provide some indication of the relative magnitude of these forces. Here, by analyzing 384 bacterial genomes assigned to three ecological categories (obligate intracellular, facultative intracellular, and extracellular), we address two, related questions: How does the density of mobile DNA vary across the Bacteria? And is there a statistically supported relationship between ecological niche and mobile element gene density? We report three findings. First, the fraction of mobile element genes in bacterial genomes ranges from 0 to 21% and decreases significantly: facultative intracellular > extracellular > obligate intracellular bacteria. Results further show that the obligate intracellular bacteria that host switch have a higher mobile DNA gene density than the obligate intracellular bacteria that are vertically transmitted. Second, while bacteria from the three ecological niches differ in their average mobile DNA contents, the ranges of mobile DNA found in each category overlap a surprising extent, suggesting bacteria with different lifestyles can tolerate similar amounts of mobile DNA. Third, mobile DNA gene densities increase with genome size across the entire dataset, and the significance of this correlation is dependent on the obligate intracellular bacteria. Further, mobile DNA gene densities do not correlate with evolutionary relationships in a 16S rDNA phylogeny. These findings statistically support a compelling link between mobile element evolution and bacterial ecology.