Behavior of restriction-modification systems as selfish mobile elements and their impact on genome evolution

Restriction–modification (RM) systems are composed of genes that encode a restriction enzyme and a modification methylase. RM systems sometimes behave as discrete units of life, like viruses and transposons. RM complexes attack invading DNA that has not been properly modified and thus may serve as a tool of defense for bacterial cells. However, any threat to their maintenance, such as a challenge by a competing genetic element (an incompatible plasmid or an allelic homologous stretch of DNA, for example) can lead to cell death through restriction breakage in the genome. This post-segregational or post-disturbance cell killing may provide the RM complexes (and any DNA linked with them) with a competitive advantage. There is evidence that they have undergone extensive horizontal transfer between genomes, as inferred from their sequence homology, codon usage bias and GC content difference. They are often linked with mobile genetic elements such as plasmids, viruses, transposons and integrons. The comparison of closely related bacterial genomes also suggests that, at times, RM genes themselves behave as mobile elements and cause genome rearrangements. Indeed some bacterial genomes that survived post-disturbance attack by an RM gene complex in the laboratory have experienced genome rearrangements. The avoidance of some restriction sites by bacterial genomes may result from selection by past restriction attacks. Both bacteriophages and bacteria also appear to use homologous recombination to cope with the selfish behavior of RM systems. RM systems compete with each other in several ways. One is competition for recognition sequences in post-segregational killing. Another is super-infection exclusion, that is, the killing of the cell carrying an RM system when it is infected with another RM system of the same regulatory specificity but of a different sequence specificity. The capacity of RM systems to act as selfish, mobile genetic elements may underlie the structure and function of RM enzymes.     

Complete Genome Analysis of Three Acinetobacter baumannii Clinical Isolates in China for Insight into the Diversification of Drug Resistance Elements

Background

The emergence and rapid spreading of multidrug-resistant Acinetobacter baumannii strains has become a major health threat worldwide. To better understand the genetic recombination related with the acquisition of drug-resistant elements during bacterial infection, we performed complete genome analysis on three newly isolated multidrug-resistant A. baumannii strains from Beijing using next-generation sequencing technology.

Methodologies/Principal Findings

Whole genome comparison revealed that all 3 strains share some common drug resistant elements including carbapenem-resistant bla _OXA-23 and tetracycline (tet) resistance islands, but the genome structures are diversified among strains. Various genomic islands intersperse on the genome with transposons and insertions, reflecting the recombination flexibility during the acquisition of the resistant elements. The blood-isolated BJAB07104 and ascites-isolated BJAB0868 exhibit high similarity on their genome structure with most of the global clone II strains, suggesting these two strains belong to the dominant outbreak strains prevalent worldwide. A large resistance island (RI) of about 121-kb, carrying a cluster of resistance-related genes, was inserted into the ATPase gene on BJAB07104 and BJAB0868 genomes. A 78-kb insertion element carrying tra-locus and bla _OXA-23 island, can be either inserted into one of the tniB gene in the 121-kb RI on the chromosome, or transformed to conjugative plasmid in the two BJAB strains. The third strains of this study, BJAB0715, which was isolated from spinal fluid, exhibit much more divergence compared with above two strains. It harbors multiple drug-resistance elements including a truncated AbaR-22-like RI on its genome. One of the unique features of this strain is that it carries both bla _OXA-23 and bla _OXA-58 genes on its genome. Besides, an Acinetobacter lwoffii adeABC efflux element was found inserted into the ATPase position in BJAB0715.

Conclusions

Our comparative analysis on currently completed Acinetobacter baumannii genomes revealed extensive and dynamic genome organizations, which may facilitate the bacteria to acquire drug-resistance elements into their genomes.     

Variation of the Virus-Related Elements within Syntenic Genomes of the Hyperthermophilic Archaeon Aeropyrum

The increasing number of genome sequences of archaea and bacteria show their adaptation to different environmental conditions at the genomic level. Aeropyrum spp. are aerobic and hyperthermophilic archaea. Aeropyrum camini was isolated from a deep-sea hydrothermal vent, and Aeropyrum pernix was isolated from a coastal solfataric vent. To investigate the adaptation strategy in each habitat, we compared the genomes of the two species. Shared genome features were a small genome size, a high GC content, and a large portion of orthologous genes (86 to 88%). The genomes also showed high synteny. These shared features may have been derived from the small number of mobile genetic elements and the lack of a RecBCD system, a recombinational enzyme complex. In addition, the specialized physiology (aerobic and hyperthermophilic) of Aeropyrum spp. may also contribute to the entire-genome similarity. Despite having stable genomes, interference of synteny occurred with two proviruses, A. pernix spindle-shaped virus 1 (APSV1) and A. pernix ovoid virus 1 (APOV1), and clustered regularly interspaced short palindromic repeat (CRISPR) elements. Spacer sequences derived from the A. camini CRISPR showed significant matches with protospacers of the two proviruses infecting A. pernix, indicating that A. camini interacted with viruses closely related to APSV1 and APOV1. Furthermore, a significant fraction of the nonorthologous genes (41 to 45%) were proviral genes or ORFans probably originating from viruses. Although the genomes of A. camini and A. pernix were conserved, we observed nonsynteny that was attributed primarily to virus-related elements. Our findings indicated that the genomic diversification of Aeropyrum spp. is substantially caused by viruses.     

Positive Selection and Horizontal Gene Transfer in the Genome of a Male-Killing Wolbachia

Wolbachia are a genus of widespread bacterial endosymbionts in which some strains can hijack or manipulate arthropod host reproduction. Male killing is one such manipulation in which these maternally transmitted bacteria benefit surviving daughters in part by removing competition with the sons for scarce resources. Despite previous findings of interesting genome features of microbial sex ratio distorters, the population genomics of male-killers remain largely uncharacterized. Here, we uncover several unique features of the genome and population genomics of four Arizonan populations of a male-killing Wolbachia strain, wInn, that infects mushroom-feeding Drosophila innubila. We first compared the wInn genome with other closely related Wolbachia genomes of Drosophila hosts in terms of genome content and confirm that the wInn genome is largely similar in overall gene content to the wMel strain infecting D. melanogaster. However, it also contains many unique genes and repetitive genetic elements that indicate lateral gene transfers between wInn and non-Drosophila eukaryotes. We also find that, in line with literature precedent, genes in the Wolbachia prophage and Octomom regions are under positive selection. Of all the genes under positive selection, many also show evidence of recent horizontal transfer among Wolbachia symbiont genomes. These dynamics of selection and horizontal gene transfer across the genomes of several Wolbachia strains and diverse host species may be important underlying factors in Wolbachia’s success as a male-killer of divergent host species.     

Characterization of Imcrop, a Mutator-like MITE family in the rice genome

Sequence comparisons of ammonium transporter 1?C2 genes (OsAMT1-2) in different rice accessions revealed a MITE insertion in the upstream region of the gene. The 391-bp MITE, classified as a Mutator superfamily member and named Imcrop, included terminal inverted repeat (TIR) and 9-bp target site duplication (TSD) sequences. We identified 151 Imcrop elements dispersed on 12 chromosomes of the japonica reference genome. Of these, 12.6% were found in genic regions and 33.1% were located within 1.5 kb of annotated rice genes. We constructed comparative insertion maps with 111 and 102 intact Imcrop elements in the japonica and indica reference genomes, respectively. The Imcrop elements showed relatively even distribution across all chromosomes although their frequency was higher on chromosomes 1, 3, and 4 in both genomes. Seventy seven Imcrop elements were detected in both subspecies, whereas 34 and 25 insertions were found only in the japonica or indica genome, respectively. We compared insertion polymorphisms of 19 Imcrop elements found inside genes in 48 Korean rice cultivars, consisting of 42 japonica and six Tongil-types (indica-japonica cross). Thirteen insertions were common to all cultivars indicating these elements were present before indica-japonica divergence. The six other elements showed insertion polymorphisms among accessions, showing their recent insertion history or no critical positive effect of their insertion on the rice genome.     

Mitochondrial genomes of the genus Ephydatia Lamouroux, 1816: can palindromic elements be used in species-level studies?

Poriferan mitochondrial DNA (mtDNA), especially large intergenic regions, is a target for the insertion of repetitive hairpin-forming elements. These elements are responsible for the large mt genome size differences observed even among closely related sponge taxa. In this study, we present the new, nearly complete, mt genome sequence of Ephydatia fluviatilis and compare it with previously published mt genomes of freshwater sponges. Special emphasis was placed on comparison with the closely related species Ephydatia muelleri, thereby comparing the only two species of the genus Ephydatia on the western Balkan Peninsula. In particular, we analyzed repetitive palindromic elements within the mitochondrial intergenic regions. The genomic distribution of these repetitive elements was analyzed and their potential role in the evolution of mt genomes discussed. We show here that palindromic elements are widespread through the whole mt genome, including the protein coding genes, thus introducing genetic variability into mt genomes.     

MIPS Arabidopsis thaliana Database (MAtDB): an integrated biological knowledge resource based on the first complete plant genome

Arabidopsis thaliana is the first plant for which the complete genome has been sequenced and published. Annotation of complex eukaryotic genomes requires more than the assignment of genetic elements to the sequence. Besides completing the list of genes, we need to discover their cellular roles, their regulation and their interactions in order to understand the workings of the whole plant. The MIPS Arabidopsis thaliana Database (MAtDB; http://mips.gsf.de/proj/thal/db) started out as a repository for genome sequence data in the European Scientists Sequencing Arabidopsis (ESSA) project and the Arabidopsis Genome Initiative. Our aim is to transform MAtDB into an integrated biological knowledge resource by integrating diverse data, tools, query and visualization capabilities and by creating a comprehensive resource for Arabidopsis as a reference model for other species, including crop plants.     

Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

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

Ma-LMM01 infecting toxic Microcystis aeruginosa illuminates diverse cyanophage genome strategies

Cyanobacteria and their phages are significant microbial components of the freshwater and marine environments. We identified a lytic phage, Ma-LMM01, infecting Microcystis aeruginosa, a cyanobacterium that forms toxic blooms on the surfaces of freshwater lakes. Here, we describe the first sequenced freshwater cyanomyovirus genome of Ma-LMM01. The linear, circularly permuted, and terminally redundant genome has 162,109 bp and contains 184 predicted protein-coding genes and two tRNA genes. The genome exhibits no colinearity with previously sequenced genomes of cyanomyoviruses or other Myoviridae. The majority of the predicted genes have no detectable homologues in the databases. These findings indicate that Ma-LMM01 is a member of a new lineage of the Myoviridae family. The genome lacks homologues for the photosynthetic genes that are prevalent in marine cyanophages. However, it has a homologue of nblA, which is essential for the degradation of the major cyanobacteria light-harvesting complex, the phycobilisomes. The genome codes for a site-specific recombinase and two prophage antirepressors, suggesting that it has the capacity to integrate into the host genome. Ma-LMM01 possesses six genes, including three coding for transposases, that are highly similar to homologues found in cyanobacteria, suggesting that recent gene transfers have occurred between Ma-LMM01 and its host. We propose that the Ma-LMM01 NblA homologue possibly reduces the absorption of excess light energy and confers benefits to the phage living in surface waters. This phage genome study suggests that light is central in the phage-cyanobacterium relationships where the viruses use diverse genetic strategies to control their host's photosynthesis.     

Genome-Wide Analysis of Repetitive Elements in Papaya

Papaya (Carica papaya L.) is an important fruit crop cultivated in tropical and subtropical regions worldwide. A first draft of its genome sequence has been recently released. Together with Arabidopsis, rice, poplar, grapevine and other genomes in the pipeline, it represents a good opportunity to gain insight into the organization of plant genomes. Here we report a detailed analysis of repetitive elements in the papaya genome, including transposable elements (TEs), tandemly-arrayed sequences, and high copy number genes. These repetitive sequences account for ～56% of the papaya genome with TEs being the most abundant at 52%, tandem repeats at 1.3% and high copy number genes at 3%. Most common types of TEs are represented in the papaya genome with retrotransposons being the dominant class, accounting for 40% of the genome. The most prevalent retrotransposons are Ty3-gypsy (27.8%) and Ty1-copia (5.5%). Among the tandem repeats, microsatellites are the most abundant in number, but represent only 0.19% of the genome. Minisatellites and satellites are less abundant, but represent 0.68% and 0.43% of the genome, respectively, due to greater repeat length. Despite an overall smaller gene repertoire in papaya than many other angiosperms, a significant fraction of genes (>2%) are present in large gene families with copy number greater than 20. This repeat database clarified a major part of the papaya genome organization and partly explained the lower gene repertoire in papaya than in Arabidopsis.     

Genome Reduction Promotes Increase in Protein Functional Complexity in Bacteria

Obligate pathogenic and endosymbiotic bacteria typically experience gene loss due to functional redundancy, asexuality, and genetic drift. We hypothesize that reduced genomes increase their functional complexity through protein multitasking, in which many genes adopt new roles to counteract gene loss. Comparisons of interaction networks among six bacteria that have varied genome sizes (Mycoplasma pneumoniae, Treponema pallidum, Helicobacter pylori, Campylobacter jejuni, Synechocystis sp., and Mycobacterium tuberculosis) reveal that proteins in small genomes interact with proteins from a wider range of functions than do their orthologs in larger genomes. This suggests that surviving proteins form increasingly complex functional relationships to compensate for genes that are lost.     

What It Takes to Be a Pseudomonas aeruginosa? The Core Genome of the Opportunistic Pathogen Updated

Pseudomonas aeruginosa is an opportunistic bacterial pathogen able to thrive in highly diverse ecological niches and to infect compromised patients. Its genome exhibits a mosaic structure composed of a core genome into which accessory genes are inserted en bloc at specific sites. The size and the content of the core genome are open for debate as their estimation depends on the set of genomes considered and the pipeline of gene detection and clustering. Here, we redefined the size and the content of the core genome of P. aeruginosa from fully re-analyzed genomes of 17 reference strains. After the optimization of gene detection and clustering parameters, the core genome was defined at 5,233 orthologs, which represented ~ 88% of the average genome. Extrapolation indicated that our panel was suitable to estimate the core genome that will remain constant even if new genomes are added. The core genome contained resistance determinants to the major antibiotic families as well as most metabolic, respiratory, and virulence genes. Although some virulence genes were accessory, they often related to conserved biological functions. Long-standing prophage elements were subjected to a genetic drift to eventually display a G+C content as higher as that of the core genome. This contrasts with the low G+C content of highly conserved ribosomal genes. The conservation of metabolic and respiratory genes could guarantee the ability of the species to thrive on a variety of carbon sources for energy in aerobiosis and anaerobiosis. Virtually all the strains, of environmental or clinical origin, have the complete toolkit to become resistant to the major antipseudomonal compounds and possess basic pathogenic mechanisms to infect humans. The knowledge of the genes shared by the majority of the P. aeruginosa isolates is a prerequisite for designing effective therapeutics to combat the wide variety of human infections.     

Definition of the minimal cis-acting sequences necessary for genome maturation of the herpesvirus murine cytomegalovirus

Herpesvirus DNA replication proceeds via concatemeric replicative intermediates that are comprised of head-to-tail-linked genomes. Genome maturation is carried out by the terminase, a protein complex that mediates both insertion of concatemer DNA into capsids and its subsequent cleavage to release genomes within these capsids. This cleavage is sequence specific, but the governing cis-acting DNA sequences are only partially characterized. Two highly conserved motifs called pac1 and pac2 lie near the ends of herpesvirus genomes and are known to be critical for genome maturation. However, the potential importance of other sequences has not been fully investigated. We have undertaken to define all of the sequences necessary for efficient genome maturation for a herpesvirus by inserting ectopic cleavage sites into the murine cytomegalovirus genome and assessing their ability to mediate genome maturation. A combination of deletion and substitution mutations revealed that the minimal cleavage site is large (~180 bp) and complex. Sequences distal of pac1 (relative to the point of cleavage) were dispensable, suggesting that pac1 may be the sole cis-acting element on this side of the cleavage site. In contrast, a region distal to pac2 up to 150 bp from the point of cleavage was essential. Scanning substitutions revealed that the pac2 side of the cleavage site is complex and may contain multiple cis-acting sequence elements in addition to pac2. These results should facilitate the identification of trans-acting factors that bind to these elements and the elucidation of their functions. Such information will be critical for understanding the molecular basis of this complex process.     

Functional annotation of the genome unravels probiotic potential of Bacillus coagulans HS243

Spore forming Bacillus species are widely used as probiotics for human dietary supplements and in animal feeds. However, information on genetic basis of their probiotic action is obscure. Therefore, the present investigation was undertaken to elucidate probiotic traits of B. coagulans HS243 through its genome analysis. Genome mining revealed the presence of an arsenal of marker genes attributed to genuine probiotic traits. In silico analysis of HS243 genome revealed the presence of multi subunit ATPases, ADI pathway genes, chologlycine hydrolase, adhesion proteins for surviving and colonizing harsh gastric transit. HS243 genome harbored vitamin and essential amino acid biosynthetic genes, suggesting the use of HS243 as a nutrient supplement. Bacteriocin producing genes highlighted the disease preventing potential of HS243. Thus, this work established that HS243 possessed the genetic repertoire required for surviving harsh gastric transit and conferring health benefits to the host which were further validated by wet lab evidences.     

Passport, a native Tc1 transposon from flatfish, is functionally active in vertebrate cells

The Tc1/mariner family of DNA transposons is widespread across fungal, plant and animal kingdoms, and thought to contribute to the evolution of their host genomes. To date, an active Tc1 transposon has not been identified within the native genome of a vertebrate. We demonstrate that Passport, a native transposon isolated from a fish (Pleuronectes platessa), is active in a variety of vertebrate cells. In transposition assays, we found that the Passport transposon system improved stable cellular transgenesis by 40-fold, has an apparent preference for insertion into genes, and is subject to overproduction inhibition like other Tc1 elements. Passport represents the first vertebrate Tc1 element described as both natively intact and functionally active, and given its restricted phylogenetic distribution, may be contemporaneously active. The Passport transposon system thus complements the available genetic tools for the manipulation of vertebrate genomes, and may provide a unique system for studying the infiltration of vertebrate genomes by Tc1 elements.     

Massive colonization of protein-coding exons by selfish genetic elements in Paramecium germline genomes

Ciliates are unicellular eukaryotes with both a germline genome and a somatic genome in the same cytoplasm. The somatic macronucleus (MAC), responsible for gene expression, is not sexually transmitted but develops from a copy of the germline micronucleus (MIC) at each sexual generation. In the MIC genome of Paramecium tetraurelia, genes are interrupted by tens of thousands of unique intervening sequences called internal eliminated sequences (IESs), which have to be precisely excised during the development of the new MAC to restore functional genes. To understand the evolutionary origin of this peculiar genomic architecture, we sequenced the MIC genomes of 9 Paramecium species (from approximately 100 Mb in Paramecium aurelia species to >1.5 Gb in Paramecium caudatum). We detected several waves of IES gains, both in ancestral and in more recent lineages. While the vast majority of IESs are single copy in present-day genomes, we identified several families of mobile IESs, including nonautonomous elements acquired via horizontal transfer, which generated tens to thousands of new copies. These observations provide the first direct evidence that transposable elements can account for the massive proliferation of IESs in Paramecium. The comparison of IESs of different evolutionary ages indicates that, over time, IESs shorten and diverge rapidly in sequence while they acquire features that allow them to be more efficiently excised. We nevertheless identified rare cases of IESs that are under strong purifying selection across the aurelia clade. The cases examined contain or overlap cellular genes that are inactivated by excision during development, suggesting conserved regulatory mechanisms. Similar to the evolution of introns in eukaryotes, the evolution of Paramecium IESs highlights the major role played by selfish genetic elements in shaping the complexity of genome architecture and gene expression.

A comparative genomics study of nine Paramecium species reveals successful invasion of genes by transposable elements in their germline genomes, showing that the internal eliminated sequences (IESs) followed an evolutionary trajectory remarkably similar to that of spliceosomal introns.     

Genome-Scale Genotype-Phenotype Matching of Two Lactococcus lactis Isolates from Plants Identifies Mechanisms of Adaptation to the Plant Niche

Lactococcus lactis is a primary constituent of many starter cultures used for the manufacturing of fermented dairy products, but the species also occurs in various nondairy niches such as (fermented) plant material. Three genome sequences of L. lactis dairy strains (IL-1403, SK11, and MG1363) are publicly available. An extensive molecular and phenotypic diversity analysis was now performed on two L. lactis plant isolates. Diagnostic sequencing of their genomes resulted in over 2.5 Mb of sequence for each strain. A high synteny was found with the genome of L. lactis IL-1403, which was used as a template for contig mapping and locating deletions and insertions in the plant L. lactis genomes. Numerous genes were identified that do not have homologs in the published genome sequences of dairy L. lactis strains. Adaptation to growth on substrates derived from plant cell walls is evident from the presence of gene sets for the degradation of complex plant polymers such as xylan, arabinan, glucans, and fructans but also for the uptake and conversion of typical plant cell wall degradation products such as α-galactosides, β-glucosides, arabinose, xylose, galacturonate, glucuronate, and gluconate. Further niche-specific differences are found in genes for defense (nisin biosynthesis), stress response (nonribosomal peptide synthesis and various transporters), and exopolysaccharide biosynthesis, as well as the expected differences in various mobile elements such as prophages, plasmids, restriction-modification systems, and insertion sequence elements. Many of these genes were identified for the first time in Lactococcus lactis. In most cases good correspondence was found with the phenotypic characteristics of these two strains.     

Whole genome comparisons of serotype 4b and 1/2a strains of the food-borne pathogen Listeria monocytogenes reveal new insights into the core genome components of this species

The genomes of three strains of Listeria monocytogenes that have been associated with food-borne illness in the USA were subjected to whole genome comparative analysis. A total of 51, 97 and 69 strain-specific genes were identified in L.monocytogenes strains F2365 (serotype 4b, cheese isolate), F6854 (serotype 1/2a, frankfurter isolate) and H7858 (serotype 4b, meat isolate), respectively. Eighty-three genes were restricted to serotype 1/2a and 51 to serotype 4b strains. These strain- and serotype-specific genes probably contribute to observed differences in pathogenicity, and the ability of the organisms to survive and grow in their respective environmental niches. The serotype 1/2a-specific genes include an operon that encodes the rhamnose biosynthetic pathway that is associated with teichoic acid biosynthesis, as well as operons for five glycosyl transferases and an adenine-specific DNA methyltransferase. A total of 8603 and 105 050 high quality single nucleotide polymorphisms (SNPs) were found on the draft genome sequences of strain H7858 and strain F6854, respectively, when compared with strain F2365. Whole genome comparative analyses revealed that the L.monocytogenes genomes are essentially syntenic, with the majority of genomic differences consisting of phage insertions, transposable elements and SNPs.