Draft genome sequence of Sinorhizobium meliloti RU11/001, a model organism for flagellum structure,motility and chemotaxis

Sinorhizobium meliloti of the order Rhizobiales is a symbiotic nitrogen-fixing bacterium nodulating plants of the genera Medicago, Trigonella and Melilotus and hence is of great agricultural importance. In its free-living state it is motile and capable of modulating its movement patterns in response to chemical attractants. Here, the draft genome consisting of a circular chromosome, the megaplasmids pSymA and pSymB and three accessory plasmids of Sinorhizobium meliloti RU11/001, a model organism for flagellum structure, motility and chemotaxis, is reported.     

Different dynamics of genome content shuffling among host-specificity groups of the symbiotic actinobacterium Frankia

Background

Frankia is a genus of soil actinobacteria forming nitrogen-fixing root-nodule symbiotic relationships with non-leguminous woody plant species, collectively called actinorhizals, from eight dicotyledonous families. Frankia strains are classified into four host-specificity groups (HSGs), each of which exhibits a distinct host range. Genome sizes of representative strains of Alnus, Casuarina, and Elaeagnus HSGs are highly diverged and are positively correlated with the size of their host ranges.

Results

The content and size of 12 Frankia genomes were investigated by in silico comparative genome hybridization and pulsed-field gel electrophoresis, respectively. Data were collected from four query strains of each HSG and compared with those of reference strains possessing completely sequenced genomes. The degree of difference in genome content between query and reference strains varied depending on HSG. Elaeagnus query strains were missing the greatest number (22–32%) of genes compared with the corresponding reference genome; Casuarina query strains lacked the fewest (0–4%), with Alnus query strains intermediate (14–18%). In spite of the remarkable gene loss, genome sizes of Alnus and Elaeagnus query strains were larger than would be expected based on total length of the absent genes. In contrast, Casuarina query strains had smaller genomes than expected.

Conclusions

The positive correlation between genome size and host range held true across all investigated strains, supporting the hypothesis that size and genome content differences are responsible for observed diversity in host plants and host plant biogeography among Frankia strains. In addition, our results suggest that different dynamics of shuffling of genome content have contributed to these symbiotic and biogeographic adaptations. Elaeagnus strains, and to a lesser extent Alnus strains, have gained and lost many genes to adapt to a wide range of environments and host plants. Conversely, rather than acquiring new genes, Casuarina strains have discarded genes to reduce genome size, suggesting an evolutionary orientation towards existence as specialist symbionts.

Electronic supplementary material

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

Complete genome sequence of Yersinia pestis strain 91001, an isolate avirulent to humans.

Genomics provides an unprecedented opportunity to probe in minute detail into the genomes of the world's most deadly pathogenic bacteria- Yersinia pestis. Here we report the complete genome sequence of Y. pestis strain 91001, a human-avirulent strain isolated from the rodent Brandt's vole-Microtus brandti. The genome of strain 91001 consists of one chromosome and four plasmids (pPCP1, pCD1, pMT1 and pCRY). The 9609-bp pPCP1 plasmid of strain 91001 is almost identical to the counterparts from reference strains (CO92 and KIM). There are 98 genes in the 70,159-bp range of plasmid pCD1. The 106,642-bp plasmid pMT1 has slightly different architecture compared with the reference ones. pCRY is a novel plasmid discovered in this work. It is 21,742 bp long and harbors a cryptic type IV secretory system. The chromosome of 91001 is 4,595,065 bp in length. Among the 4037 predicted genes, 141 are possible pseudo-genes. Due to the rearrangements mediated by insertion elements, the structure of the 91001 chromosome shows dramatic differences compared with CO92 and KIM. Based on the analysis of plasmids and chromosome architectures, pseudogene distribution, nitrate reduction negative mechanism and gene comparison, we conclude that strain 91001 and other strains isolated from M. brandti might have evolved from ancestral Y. pestis in a different lineage. The large genome fragment deletions in the 91001 chromosome and some pseudogenes may contribute to its unique nonpathogenicity to humans and host-specificity.     

Prevalence of pSmeSM11a-like plasmids in indigenous Sinorhizobium meliloti strains isolated in the course of a field release experiment with genetically modified S. meliloti strains

Plasmid pSmeSM11a, residing in the indigenous Sinorhizobium meliloti strain SM11 originating from a field in Strassmoos (Bavaria, Germany), was analysed previously at the genomic level. Thirty-seven indigenous S. meliloti strains, originating from two different locations in Germany, were screened for genes identified previously on pSmeSM11a. Seven of these strains harbour accessory plasmids that are very similar to pSmeSM11a. The identified pSmeSM11a-like plasmids are c. 130-150 kb in size and possess nearly identical restriction profiles. Up to 30 genes identified previously on pSmeSM11a could be detected on these plasmids by hybridisation experiments, e.g., the nodulation genes nodP and nodQ, the ethylene level modulation gene acdS and the taurine metabolism gene tauD. A few pSmeSM11a genes were also detected on other plasmids. The reference plasmid pSmeSM11a contains a region that is similar to a segment of S. meliloti strain Rm1021 pSymA. Regions with similarity to pSymA were also detected on the aforementioned seven pSmeSM11a-like plasmids. The specifications of these regions are nearly identical to the one on pSmeSM11a and differ from Rm1021 pSymA as determined by nucleotide sequence analysis. Two further plasmids similar to pSmeSM11a completely lack the pSymA-region. Those strains carrying accessory plasmids that contain the acdS gene encoding 1-aminocyclopropane-1-carboxylate deaminase are able to grow on 1-aminocyclopropane-1-carboxylate as the sole source of nitrogen, demonstrating functionality of the acdS gene product. About 36% of the analysed plasmids, including three pSmeSM11a-like plasmids, could be transferred to another S. meliloti recipient strain, allowing for their dissemination in S. meliloti populations.     

Comparative genomic hybridisation and ultrafast pyrosequencing revealed remarkable differences between the Sinorhizobium meliloti genomes of the model strain Rm1021 and the field isolate SM11

Genomic variation between the Sinorhizobium meliloti model strain Rm1021 and the field isolate SM11 was assessed by using the genome-wide S. meliloti Rm1021 Sm6k-oligonucleotide microarray in a comparative genomic hybridisation experiment. Several gene clusters present in the Rm1021 genome are missing in the SM11 genome. In detail, three missing gene clusters were identified for the chromosome, five for megaplasmid pSymA and two for megaplasmid pSymB. To confirm these hybridisation results, the draft genome sequence of the S. meliloti field isolate SM11 was established by 454-pyrosequencing. Three sequencing runs on the ultrafast Genome Sequencer 20 System yielded 112.5 million bases. These could be assembled into 905 larger contigs resulting in a nearly 15-fold coverage of the 7.1Mb SM11 genome. The missing gene regions identified by comparative genomic hybridisation could be confirmed by the results of the 454-sequencing project. An in-depth analysis of these gene regions resulted in the following findings: (i) a complete type I restriction/modification system encoded by a composite transposon is absent in the chromosome of strain SM11. (ii) Most of the Rm1021 denitrification genes and the complete siderophore biosynthesis operon were found to be missing on SM11 megaplasmid pSymA. (iii) S. meliloti SM11 megaplasmid pSymB lacks a complete cell surface carbohydrate synthesis gene cluster. (iv) Several genes that are absent in the SM11 genome could be assigned to insertion sequences and transposons.     

Complete genome sequencing of Agrobacterium sp. H13-3, the former Rhizobium lupini H13-3, reveals a tripartite genome consisting of a circular and a linear chromosome and an accessory plasmid but lacking a tumor-inducing Ti-plasmid

Agrobacterium sp. H13-3, formerly known as Rhizobium lupini H13-3, is a soil bacterium that was isolated from the rhizosphere of Lupinus luteus. The isolate has been established as a model system for studying novel features of flagellum structure, motility and chemotaxis within the family Rhizobiaceae. The complete genome sequence of Agrobacterium sp. H13-3 has been established and the genome structure and phylogenetic assignment of the organism was analysed. For de novo sequencing of the Agrobacterium sp. H13-3 genome, a combined strategy comprising 454-pyrosequencing on the Genome Sequencer FLX platform and PCR-based amplicon sequencing for gap closure was applied. The finished genome consists of three replicons and comprises 5,573,770 bases. Based on phylogenetic analyses, the isolate could be assigned to the genus Agrobacterium biovar I and represents a genomic species G1 strain within this biovariety. The highly conserved circular chromosome (2.82 Mb) of Agrobacterium sp. H13-3 mainly encodes housekeeping functions characteristic for an aerobic, heterotrophic bacterium. Agrobacterium sp. H13-3 is a motile bacterium driven by the rotation of several complex flagella. Its behaviour towards external stimuli is regulated by a large chemotaxis regulon and a total of 17 chemoreceptors. Comparable to the genome of Agrobacterium tumefaciens C58, Agrobacterium sp. H13-3 possesses a linear chromosome (2.15 Mb) that is related to its reference replicon and features chromosomal and plasmid-like properties. The accessory plasmid pAspH13-3a (0.6 Mb) is only distantly related to the plasmid pAtC58 of A. tumefaciens C58 and shows a mosaic structure. A tumor-inducing Ti-plasmid is missing in the sequenced strain H13-3 indicating that it is a non-virulent isolate.     

The regulatory status of the fixL- and fixJ-like genes in Bradyrhizobium japonicum may be different from that in Rhizobium meliloti

Summary The cloning, sequencing and mutational analysis of the Bradyrhizobium japonicum symbiotic nitrogen fixation genes fixL and fixJ are reported here. The two genes were adjacent and probably formed an operon, fixLJ. The predicted FixL and FixJ proteins, members of the two-component sensor/regulator family, were homologous over almost their entire lengths to the corresponding Rhizobium meliloti proteins (approx. 50% identity). Downstream of the B. japonicum fixJ gene was found an open reading frame with 138 codons (ORF138) whose product shared 36% homology with the N-terminal part of FixJ. Deletion and insertion mutations within fixL and fixJ led to a loss of approximately 90% wildtype symbiotic nitrogen fixation (Fix) activity, whereas an ORF138 mutant was Fix⁺. In fixL, fixJ and ORF138 mutant backgrounds, the aerobic expression of the fixR-nifA operon was not affected. NifA itself did not regulate the expression of the fixJ gene. Thus, the B. japonicum FixL and FixJ proteins were neither involved in the regulation of aerobic nifA gene expression nor in the anaerobic NifA-dependent autoregulation of the fixRnifA operon; rather they appeared to control symbiotically important genes other than those whose expression was dependent on the NifA protein. The fixL and fixJ mutant strains were unable to grow anaerobically with nitrate as the terminal electron acceptor. Therefore, some of the FixJ-dependent genes in B. japonicum may be concerned with anaerobic respiration.     

Complete sequence and analysis of the plastid genome of the unicellular red alga Cyanidioschyzon merolae.

The complete nucleotide sequence of the plastid genome of the unicellular primitive red alga Cyanidioschyzon merolae 10D (Cyanidiophyceae) was determined. The genome is a circular DNA composed of 149,987 bp with no inverted repeats. The G + C content of this plastid genome is 37.6%. The C. merolae plastid genome contains 243 genes, which are distributed on both strands and consist of 36 RNA genes (3 rRNAs, 31 tRNAs, tmRNA, and a ribonuclease P RNA component) and 207 protein genes, including unidentified open reading frames. The striking feature of this genome is the high degree of gene compaction; it has very short intergenic distances (approximately 40% of the protein genes were overlapped) and no genes have introns. This genome encodes several genes that are rarely found in other plastid genomes. A gene encoding a subunit of sulfate transporter (cysW) is the first to be identified in a plastid genome. The cysT and cysW genes are located in the C. merolae plastid genome in series, and they probably function together with other nuclear-encoded components of the sulfate transport system. Our phylogenetic results suggest that the Cyanidiophyceae, including C. merolae, are a basal clade within the red lineage plastids.     

Single-plant,Sterile Microcosms for Nodulation and Growth of the Legume Plant Medicago truncatula with the Rhizobial Symbiont Sinorhizobium meliloti

Rhizobial bacteria form symbiotic, nitrogen-fixing nodules on the roots of compatible host legume plants. One of the most well-developed model systems for studying these interactions is the plant Medicago truncatula cv. Jemalong A17 and the rhizobial bacterium Sinorhizobium meliloti 1021. Repeated imaging of plant roots and scoring of symbiotic phenotypes requires methods that are non-destructive to either plants or bacteria. The symbiotic phenotypes of some plant and bacterial mutants become apparent after relatively short periods of growth, and do not require long-term observation of the host/symbiont interaction. However, subtle differences in symbiotic efficiency and nodule senescence phenotypes that are not apparent in the early stages of the nodulation process require relatively long growth periods before they can be scored. Several methods have been developed for long-term growth and observation of this host/symbiont pair. However, many of these methods require repeated watering, which increases the possibility of contamination by other microbes. Other methods require a relatively large space for growth of large numbers of plants. The method described here, symbiotic growth of M. truncatula/S. meliloti in sterile, single-plant microcosms, has several advantages. Plants in these microcosms have sufficient moisture and nutrients to ensure that watering is not required for up to 9 weeks, preventing cross-contamination during watering. This allows phenotypes to be quantified that might be missed in short-term growth systems, such as subtle delays in nodule development and early nodule senescence. Also, the roots and nodules in the microcosm are easily viewed through the plate lid, so up-rooting of the plants for observation is not required.     

The complete mitochondrial genome sequence of the little egret (Egretta garzetta)

Many phylogenetic questions in the Ciconiiformes remain unresolved and complete mitogenome data are urgently needed for further molecular investigation. In this work, we determined the complete mitogenome sequence of the little egret (Egretta garzetta). The genome was 17,361 bp in length and the gene organization was typical of other avian mtDNA. In protein-coding genes (PCGs), a C insertion was found in ND3, and COIII and ND4 terminated with incomplete stop codons (T). tRNA-Val and tRNA-Ser (AGY) were unable to fold into canonical cloverleaf secondary structures because they had lost the DHU arms. Long repetitive sequences consisting of five types of tandem repeats were found at the 3′ end of Domain III in the control region. A phylogenetic analysis of 11 species of Ciconiiformes was done using complete mitogenome data and 12 PCGs. The tree topologies obtained with these two strategies were identical, which strongly confirmed the monophyly of Ardeidae, Threskiorothidae and Ciconiidae. The phylogenetic analysis also revealed that Egretta was more closely related to Ardea than to Nycticorax in the Ardeidae, and Platalea was more closely related to Threskiornis than to Nipponia in the Threskiornithidae. These findings contribute to our understanding of the phylogenetic relationships of Ciconiiformes based on complete mitogenome data.     

The complete mitochondrial genome sequence of the black-capped capuchin (Cebus apella)

The phylogenetic relationships of primates have been extensively investigated, but key issues remain unresolved. Complete mitochondrial genome (mitogenome) data have many advantages in phylogenetic analyses, but such data are available for only 46 primate species. In this work, we determined the complete mitogenome sequence of the black-capped capuchin (Cebus apella). The genome was 16,538 bp in size and consisted of 13 protein-coding genes, 22 tRNAs, two rRNAs and a control region. The genome organization, nucleotide composition and codon usage did not differ significantly from those of other primates. The control region contained several distinct repeat motifs, including a putative termination-associated sequence (TAS) and several conserved sequence blocks (CSB-F, E, D, C, B and 1). Among the protein-coding genes, the COII gene had lower nonsynonymous and synonymous substitutions rates while the ATP8 and ND4 genes had higher rates. A phylogenetic analysis using Maximum likelihood and Bayesian methods and the complete mitogenome data for platyrrhine species confirmed the basal position of the Callicebinae and the sister relationship between Atelinae and Cebidae, as well as the sister relationship between Aotinae (Aotus) and Cebinae (Cebus/Saimiri) in Cebidae. These conclusions agreed with the most recent molecular phylogenetic investigations on primates. This work provides a framework for the use of complete mitogenome information in phylogenetic analyses of the Platyrrhini and primates in general.     

Complete genome structure of the nitrogen-fixing symbiotic bacterium Mesorhizobium loti.

The complete nucleotide sequence of the genome of a symbiotic bacterium Mesorhizobium loti strain MAFF303099 was determined. The genome of M. loti consisted of a single chromosome (7,036,071 bp) and two plasmids, designated as pMLa (351,911 bp) and pMLb (208, 315 bp). The chromosome comprises 6752 potential protein-coding genes, two sets of rRNA genes and 50 tRNA genes representing 47 tRNA species. Fifty-four percent of the potential protein genes showed sequence similarity to genes of known function, 21% to hypothetical genes, and the remaining 25% had no apparent similarity to reported genes. A 611-kb DNA segment, a highly probable candidate of a symbiotic island, was identified, and 30 genes for nitrogen fixation and 24 genes for nodulation were assigned in this region. Codon usage analysis suggested that the symbiotic island as well as the plasmids originated and were transmitted from other genetic systems. The genomes of two plasmids, pMLa and pMLb, contained 320 and 209 potential protein-coding genes, respectively, for a variety of biological functions. These include genes for the ABC-transporter system, phosphate assimilation, two-component system, DNA replication and conjugation, but only one gene for nodulation was identified.     

The (in)complete organelle genome: exploring the use and nonuse of available technologies for characterizing mitochondrial and plastid chromosomes

Not long ago, scientists paid dearly in time, money and skill for every nucleotide that they sequenced. Today, DNA sequencing technologies epitomize the slogan ‘faster, easier, cheaper and more’, and in many ways, sequencing an entire genome has become routine, even for the smallest laboratory groups. This is especially true for mitochondrial and plastid genomes. Given their relatively small sizes and high copy numbers per cell, organelle DNAs are currently among the most highly sequenced kind of chromosome. But accurately characterizing an organelle genome and the information it encodes can require much more than DNA sequencing and bioinformatics analyses. Organelle genomes can be surprisingly complex and can exhibit convoluted and unconventional modes of gene expression. Unravelling this complexity can demand a wide assortment of experiments, from pulsed‐field gel electrophoresis to Southern and Northern blots to RNA analyses. Here, we show that it is exactly these types of ‘complementary’ analyses that are often lacking from contemporary organelle genome papers, particularly short ‘genome announcement’ articles. Consequently, crucial and interesting features of organelle chromosomes are going undescribed, which could ultimately lead to a poor understanding and even a misrepresentation of these genomes and the genes they express. High‐throughput sequencing and bioinformatics have made it easy to sequence and assemble entire chromosomes, but they should not be used as a substitute for or at the expense of other types of genomic characterization methods.     

The complete chloroplast genome sequence of Pelargonium x hortorum: organization and evolution of the largest and most highly rearranged chloroplast genome of land plants

The chloroplast genome of Pelargonium x hortorum has been completely sequenced. It maps as a circular molecule of 217,942 bp and is both the largest and most rearranged land plant chloroplast genome yet sequenced. It features 2 copies of a greatly expanded inverted repeat (IR) of 75,741 bp each and, consequently, diminished single-copy regions of 59,710 and 6,750 bp. Despite the increase in size and complexity of the genome, the gene content is similar to that of other angiosperms, with the exceptions of a large number of pseudogenes, the recognition of 2 open reading frames (ORF56 and ORF42) in the trnA intron with similarities to previously identified mitochondrial products (ACRS and pvs-trnA), the losses of accD and trnT-ggu and, in particular, the presence of a highly divergent set of rpoA-like ORFs rather than a single, easily recognized gene for rpoA. The 3-fold expansion of the IR (relative to most angiosperms) accounts for most of the size increase of the genome, but an additional 10% of the size increase is related to the large number of repeats found. The Pelargonium genome contains 35 times as many 31 bp or larger repeats than the unrearranged genome of Spinacia. Most of these repeats occur near the rearrangement hotspots, and 2 different associations of repeats are localized in these regions. These associations are characterized by full or partial duplications of several genes, most of which appear to be nonfunctional copies or pseudogenes. These duplications may also be linked to the disruption of at least 1 but possibly 2 or 3 operons. We propose simple models that account for the major rearrangements with a minimum of 8 IR boundary changes and 12 inversions in addition to several insertions of duplicated sequence.     

Sequence analysis of the 181-kb accessory plasmid pSmeSM11b, isolated from a dominant Sinorhizobium meliloti strain identified during a long-term field release experiment

The 181 251 bp accessory plasmid pSmeSM11b of Sinorhizobium meliloti strain SM11, belonging to a dominant indigenous S. meliloti subpopulation identified during a long-term field release experiment, was sequenced. This plasmid has 166 coding sequences (CDSs), 42% of which encode proteins with homology to proteins of known function. Plasmid pSmeSM11b is a member of the repABC replicon family and contains a large gene region coding for a conjugation system similar to that of other self-transmissible plasmids in Rhizobium and Agrobacterium. Another pSmeSM11b gene region, possibly involved in sugar metabolism and polysaccharide catabolism, resembled a region of S. meliloti 1021 megaplasmid pSymB and in the genome of Sinorhizobium medicae WSM419. Another module of plasmid pSmeSM11b encodes proteins similar to those of the nitrogen-fixing actinomycete Frankia CcI3, and which are likely to be involved in the synthesis of a secondary metabolite. Several ORFs of pSmeSM11b were predicted to play a role in nonribosomal peptide synthesis. Plasmid pSmeSM11b has many mobile genetic elements, which contribute to the mosaic composition of the plasmid.     

The chloroplast genome sequence of Chara vulgaris sheds new light into the closest green algal relatives of land plants

The phylum Streptophyta comprises all land plants and six monophyletic groups of charophycean green algae (Mesostigmatales, Chlorokybales, Klebsormidiales, Zygnematales, Coleochaetales, and Charales). Phylogenetic analyses of four genes encoded in three cellular compartments suggest that the Charales are sister to land plants and that charophycean green algae evolved progressively toward an increasing cellular complexity. To validate this phylogenetic hypothesis and to understand how and when the highly conservative pattern displayed by land plant chloroplast DNAs (cpDNAs) originated in the Streptophyta, we have determined the complete chloroplast genome sequence (184,933 bp) of a representative of the Charales, Chara vulgaris, and compared this genome to those of Mesostigma (Mesostigmatales), Chlorokybus (Chlorokybales), Staurastrum and Zygnema (Zygnematales), Chaetosphaeridium (Coleochaetales), and selected land plants. The phylogenies we inferred from 76 cpDNA-encoded proteins and genes using various methods favor the hypothesis that the Charales diverged before the Coleochaetales and Zygnematales. The Zygnematales were identified as sister to land plants in the best tree topology (T1), whereas Chaetosphaeridium (T2) or a clade uniting the Zygnematales and Chaetosphaeridium (T3) occupied this position in alternative topologies. Chara remained at the same basal position in trees including more land plant taxa and inferred from 56 proteins/genes. Phylogenetic inference from gene order data yielded two most parsimonious trees displaying the T1 and T3 topologies. Analyses of additional structural cpDNA features (gene order, gene content, intron content, and indels in coding regions) provided better support for T1 than for the topology of the above-mentioned four-gene tree. Our structural analyses also revealed that many of the features conserved in land plant cpDNAs were inherited from their green algal ancestors. The intron content data predicted that at least 15 of the 21 land plant group II introns were gained early during the evolution of streptophytes and that a single intron was acquired during the transition from charophycean green algae to land plants. Analyses of genome rearrangements based on inversions predicted no alteration in gene order during the transition from charophycean green algae to land plants.