Comparative genome analysis of four magnetotactic bacteria reveals a complex set of group-specific genes implicated in magnetosome biomineralization and function

Magnetotactic bacteria (MTB) are a heterogeneous group of aquatic prokaryotes with a unique intracellular organelle, the magnetosome, which orients the cell along magnetic field lines. Magnetotaxis is a complex phenotype, which depends on the coordinate synthesis of magnetosomes and the ability to swim and orient along the direction caused by the interaction with the Earth's magnetic field. Although a number of putative magnetotaxis genes were recently identified within a conserved genomic magnetosome island (MAI) of several MTB, their functions have remained mostly unknown, and it was speculated that additional genes located outside the MAI might be involved in magnetosome formation and magnetotaxis. In order to identify genes specifically associated with the magnetotactic phenotype, we conducted comparisons between four sequenced magnetotactic Alphaproteobacteria including the nearly complete genome of Magnetospirillum gryphiswaldense strain MSR-1, the complete genome of Magnetospirillum magneticum strain AMB-1, the complete genome of the magnetic coccus MC-1, and the comparative-ready preliminary genome assembly of Magnetospirillum magnetotacticum strain MS-1 against an in-house database comprising 426 complete bacterial and archaeal genome sequences. A magnetobacterial core genome of about 891 genes was found shared by all four MTB. In addition to a set of approximately 152 genus-specific genes shared by the three Magnetospirillum strains, we identified 28 genes as group specific, i.e., which occur in all four analyzed MTB but exhibit no (MTB-specific genes) or only remote (MTB-related genes) similarity to any genes from nonmagnetotactic organisms and which besides various novel genes include nearly all mam and mms genes previously shown to control magnetosome formation. The MTB-specific and MTB-related genes to a large extent display synteny, partially encode previously unrecognized magnetosome membrane proteins, and are either located within (18 genes) or outside (10 genes) the MAI of M. gryphiswaldense. These genes, which represent less than 1% of the 4,268 open reading frames of the MSR-1 genome, as yet are mostly of unknown functions but are likely to be specifically involved in magnetotaxis and, thus, represent prime targets for future experimental analysis.     

Cloning, characterization and taxonomic significance of genes for the 5S ribosomal RNA of Leptonema illini strain 3055.

The genes encoding the 5S ribosomal RNA (rRNA) for Leptonema illini strain 3055 were isolated and sequenced. The 5S RNA molecule encoded was 117 nucleotides long. The genome of strain 3055 contained two genes for 5S rRNA that were located close together. The nucleotide sequences of the Leptonema illini genes exhibited less similarity to the rRNA gene of Leptospira interrogans strain Moulton and also to those of typical eubacterial genes than did the rRNA genes of other leptospires. However, the overall secondary structure of the 5S rRNA encoded exhibited a strong similarity to that of typical eubacterial 5S rRNA. Southern hybridization of the 5S rRNA gene probe with the genomic DNA of strain 965, which is currently classified as Leptospira biflexa, showed the latter to have close similarity to that of strain 3055. The physical map of strain 965 was quite similar to that of strain 3055 and was greatly different from that of any other strains of L. biflexa. In the organization of 5S rRNA genes, strain 965 is sufficiently different from other members of the genus Leptospira to be regarded as a member of the genus Leptonema.     

Specifics of the coat protein gene in Russian strains of the cucumber green mottle mosaic virus

The primary structure of the coat protein (CP) gene was examined for pathogenic strain MS-1 and vaccine strain VIROG-43M of the cucumber green mottle mosaic virus (CGMMV). In CP amino acid composition, strains MS-1 and VIROG-43M are typical representatives of CGMMV: their CPs have 98–100% homology to CPs of other tobamoviruses of the group. The CP gene has the same nucleotide composition in pathogenic MS-1 and vaccine VIROG-43M, indicating that strain attenuation is not determined by this gene. The CP amino acid sequences of the two Russian strains are fully identical to the CP sequences of two Greek strains, GR-3 and GR-5. However, the nucleotide sequences of their genes differ in 13 bp, testifying to the difference between the Russian and Greek strains.     

Specifics of the coat protein gene in Russian strains of the cucumber green mottle mosaic virus

The primary structure of the coat protein (CP) gene was examined for pathogenic strain MS-1 and vaccine strain VIROG-43M of the cucumber green mottle mosaic virus (CGMMV). In CP amino acid composition, strains MS-1 and VIROG-43M are typical representatives of CGMMV: their CPs have 98-100% homology to CPs of other tobamoviruses of the group. The CP gene has the same nucleotide composition in pathogenic MS-1 and vaccine VIROG-43M, indicating that strain attenuation is not determined by this gene. The CP amino acid sequences of the two Russian strains are fully identical to the CP sequences of two Greek strains, GR-3 and GR-5. However, the nucleotide sequences of their genes differ in 13 bp, testifying to the difference between the Russian and Greek strains.     

Non-contiguous finished genome sequence and description of Anaerococcus senegalensis sp. nov

Anaerococcus senegalensis strain JC48^T sp. nov. is the type strain of A. senegalensis sp. nov. a new species within the genus Anaerococcus. This strain whose genome is described here was isolated from the fecal flora of a healthy patient. A. senegalensis is an obligate anaerobic coccus. Here we describe the features of this organism together with the complete genome sequence and annotation. The 1,790,835 bp long genome (1 chromosome but no plasmid) contains 1,721 protein-coding and 53 RNA genes including 5 rRNA genes     

Genome mapping of Arcobacter butzleri

A genome map of the potential food-poisoning pathogen Arcobacter butzleri, strain NCTC 12481, has been constructed. The map was ordered by a combination of one- and two-dimensional pulsed field gel electrophoresis, using the restriction enzymes EagI, SacII, SalI and XhoI, and Southern hybridization with digoxygenin-labelled probes. The size of the genome was estimated at 2.57+/-0.01 Mb; this is 0.9 Mb larger than the genome of the closely related species Campylobacter jejuni. We deduce that there are five copies of the 23S and 16S rRNA genes present in the A. butzleri NCTC 12481 genome and two copies of the gene glyA. A gene homologous to pglF, involved in bacterial glycosylation in C. jejuni, was also located on the genome map. In C. jejuni and Wolinella succinogenes, most glycosylation genes are clustered together and the gene order is similar. In A. butzleri NCTC 12481, no analogous organization was found in the region sequenced.     

一株高产苯乳酸的古墓土壤葡糖醋杆菌FBFS97的全基因组测序与分析

【背景】苯乳酸(phenyllactic acid,PLA)是一种应用潜力巨大的天然广谱抑菌物质。本课题组前期分离得到一株高产PLA的醋酸菌(acetic acid bacteria,AAB)——葡糖醋杆菌(Gluconacetobacter sp.)FBFS97,但尚未鉴定到种,而且其产PLA的分子机理尚不清楚。【目的】确定FBFS97的种属关系,解析FBFS97的遗传信息,特别是与PLA产生相关的基因。【方法】采用光学显微镜和扫描电镜对FBFS97的菌体形态进行表征,通过16S rRNA基因序列分析对FBFS97进行分类鉴定,并以高效液相色谱分析苯丙氨酸对其产PLA的影响。在此基础上,对FBFS97进行全基因组测序、拼接和基因预测,并进行GO/COG聚类、KEGG代谢通路和VFDB毒力等分析,以及PLA生物合成途径的预测。【结果】根据16S rRNA基因序列的比对结果,结合形态学分析,该菌被鉴定为古墓土壤葡糖醋杆菌(Gluconacetobacter tumulisoli)。将1 000 mg/L苯丙氨酸添加到FBFS97液体培养基中,发酵液中PLA最高浓度可达400 mg/L,为对照组的8倍。该菌的基因组大小为3 988 308 bp,(G+C)mol%含量为66.62%,编码基因3 500个;KEGG代谢通路分析表明,该菌基因组中存在经莽草酸途径合成PLA的所有基因;VFDB毒力预测结果显示,该菌基因组中不存在产生毒素的相关基因。【结论】首次报道了一株高产PLA的AAB——古墓土壤葡糖醋杆菌FBFS97的全基因组序列信息,并发现该菌株的基因组中含有合成PLA的所有相关基因,为后续进一步研究FBFS97产生PLA的生物合成途径提供了理论依据。     

Use of disomic strains to study the arrangement of ribosomal cistrons in saccharomyces of ribosomal cistrons in Saccharomyces cerevisiae

Summary Disomic strains ofSaccharomyces cerevisiae were studied by DNA-rRNA hybridization to examine the arrangement of rRNA cistrons on yeast chromosomes as well as to identify a disomic strain which was enriched for rRNA cistrons. Four of the five disomic strains tested showed a per cent hybridization lower than wild type. Two of these strains were found to be disomic for more than one chromosome. A slight increase in the per cent hybridization was observed with DNA isolated from one disomic strain. It was concluded that some chromosomes inSaccharomyces cerevisiae had few if any rRNA cistrons suggesting that the rRNA cistrons are non randomly distributed over the genome. From DNA-tRNA hybridization experiments, evidence for the presence of tyrosine tRNA genes on chromosomes VI was obtained.     

The genome of <Emphasis Type="Italic">Streptomyces rimosus</Emphasis> subsp. <Emphasis Type="Italic">rimosus</Emphasis> shows a novel structure compared to other <Emphasis Type="Italic">Streptomyces</Emphasis> using DNA/DNA microarray analysis

DNA/DNA genome microarray analysis together with genome sequencing suggests that the genome of members of the genus Streptomyces would seem to have a common structure including a linear genomic structure, a core of common syntenous Actinomycete genes, the presence of species specific terminal regions and two intermediate group of syntenous genes that seem to be genus specific. We analyzed Streptomyces species using DNA/DNA microarray comparative genome analysis. Only Streptomyces rimosus failed to give a congruent genome pattern for the genes found in Streptomyces coelicolor. We expanded the analysis to include a number of strains related to the type strain of S. rimosus and obtained a similar divergence from the main body of Streptomyces species. These strains showed very close identity to the original strain with no gene deletion or duplication detected. The 16S rRNA sequences of these S. rimosus strains were confirmed as very similar to the S. rimosus sequences available from the Ribosomal Database Project. When the SSU ribosomal RNA phylogeny of S. rimosus is analyzed, the species is positioned at the edge of the Streptomyces clade. We conclude that S. rimosus represents a distinct evolutionary lineage making the species a worthy possibility for genome sequencing.     

The complete mitochondrial genome of the gall-forming fly, Fergusonina taylori Nelson and Yeates (Diptera: Fergusoninidae)

The monogeneric family Fergusoninidae consists of gall-forming flies that, together with Fergusobia (Tylenchida: Neotylenchidae) nematodes, form the only known mutualistic association between insects and nematodes. In this study, the entire 16,000 bp mitochondrial genome of Fergusonina taylori Nelson and Yeates was sequenced. The circular genome contains one encoding region including 27 genes and one non-coding A+T-rich region. The arrangement of the protein-coding, ribosomal RNA (rRNA) and transfer RNA (tRNA) genes was the same as that found in the ancestral insect. Nucleotide composition is highly A+T biased. All of the protein initiation codons are ATN, except for nad1 which begins with TTT. All 22 tRNA anticodons of F. taylori match those observed in Drosophila yakuba, and all form the typical cloverleaf structure except for tRNA-Ser((AGN)) which lacks a dihydrouridine (DHU) arm. Secondary structural features of the rRNA genes of Fergusonina are similar to those proposed for other insects, with minor modifications. The mitochondrial genome of Fergusonina presented here may prove valuable for resolving the sister group to the Fergusoninidae, and expands the available mtDNA data sources for acalyptrates overall.     

Genetic relationship of two highly studied Synechococcus strains designated Anacystis nidulans.

The cyanobacteria Synechococcus sp. strain PCC 7942 and Synechococcus sp. strain PCC 6301 are very closely related and both have been designated by the binomial Anacystis nidulans. The only established difference between the two strains is the superior transformation properties of strain PCC 7942. Significant homology between the rRNA genes of these strains was demonstrated by the ability of an rRNA operon from strain PCC 6301, interrupted by a spectinomycin and streptomycin resistance marker, to transform strain PCC 7942 by recombining with and replacing an endogenous rRNA operon. Restriction fragment length polymorphism data indicated that the chromosomes of the two strains were conserved around the three psbA loci, the two rRNA operons, and the psbDI locus. However, multiple polymorphisms were detected downstream of the psbDII locus, identifying a DNA rearrangement such as an inversion, insertion, or deletion within the chromosome. Analysis of genome structure by pulsed-field gel electrophoresis of large NotI restriction fragments showed only two bands that were visibly shifted between the chromosomes of the two strains. These data support their very close genetic relationship and the feasibility of studying genes derived from strain PCC 6301 in the highly transformable PCC 7942 strain.     

Organization of the ribosomal RNA genes in Treponema phagedenis and Treponema pallidum.

The genomic DNA fragment which contains ribosomal RNA (rRNA) genes for Treponema phagedenis was cloned into bacteriophage vector lambda EMBL3. A restriction map of the fragment was constructed and the organization of the rRNA genes was determined. The fragment contained at least one copy of the 16S, 23S and 5S sequences and the genes are arranged in the order 16S-23S-5S. Southern hybridization using radiolabeled rRNA gene probes to genomic DNA from T. phagedenis strain Reiter and T. pallidum strain Nichols showed that these organisms have two radioactive fragments which hybridize to the probes in their genome. These results suggest that both pathogenic and non-pathogenic strains of Treponema may carry at least two sets of rRNA genes on their chromosomes.     

Distribution of 5S ribosomal RNA genes in somatic and germ cells of the house cricket,Acheta domesticus

In the house cricket,Acheta domesticus, the 110 genes per haploid genome encoding 18S and 28S rRNA are contained within rDNA repeats which are amplified during oogenesis. The 5S rRNA coding sequences of this cricket are found in two sizes of 5S DNA repeating units (measuring 2.1 and 3.0 kb). The 3.0 kb repeats account for more than 90% of the totalAcheta 5S DNA. We have determined the number of cricket 5S rRNA genes by RNA-DNA hybridization analysis: 310 5S DNA repeats/haploid genome clearly approximates the number of 18S and 28S rRNA genes. Because of the relatively low copy number of 5S rRNA genes the possibility of 5S DNA amplification in oocytes ofA. domesticus was also examined. Although amplification of rDNA is readily detectable, amplification of 5S DNA is not observed in oocytes ofA. domesticus. Unlike the genes coding for 18S and 28S rRNA which are localized at a single chromosomal site in the genome ofA. domesticus, the 5S rRNA genes occupy numerous sites distributed along the length of most chromosomes.     

一株羊源A型多杀性巴氏杆菌的全基因组测序及生物学特性分析

[背景] 多杀性巴氏杆菌（Pasteurella multocida,Pm）是一种革兰氏阴性菌,可引起动物和人类的呼吸道疾病和败血症等。本实验室前期分离鉴定一株A型Pm HN02菌株。[目的] 通过对HN02菌株的全基因组测序及生物信息学分析,扩充多杀性巴氏杆菌的基因组数据库信息;通过毒力基因鉴定和系统进化树分析,明确该菌株含有的毒力基因和遗传进化关系,为临床预防和诊断提供理论依据。[方法] 使用单分子实时测序（Single Molecule Real Time Sequencing,SMRT）技术对Pm HN02菌株进行全基因组测序,利用Illumina测序校正后进行基因功能注释和生物信息学分析。使用PCR鉴定菌株毒力基因,并构建进化树进行分析。[结果] Pm HN02菌株全基因组大小为2 333 292 bp,GC含量为40.15mol%,预测到的编码基因有2 389个,包含19个rRNA （6个23S rRNA、6个16S rRNA、7个5S rRNA）、62个tRNA基因、5个sRNA;含84个串联重复序列、66个小卫星DNA、2个微卫星DNA、9个基因岛、9个前噬菌体;分别有1 648、2 190和1 917个基因注释在GO、KEGG和COG数据库中,而且大部分富集于Pm的代谢过程;还有85个III型分泌系统效应蛋白、191个表型突变基因、165个毒力因子相关基因。根据分析结果绘制该菌株的全基因组圈图,并将基因组信息提交至NCBI后获得登录号cp037865。PCR鉴定发现该菌株含有fimA、toxA等14个毒力基因,缺失了tadD等毒力基因。系统进化树分析发现该菌株同北京的Pm3菌株（MH150895.1）进化关系最接近。[结论] 研究完成了A型Pm HN02株的全基因组测序和生物学特性鉴定,揭示了其同国内外Pm分离株的进化关系,为预防Pm疾病流行和探索Pm致病机制提供了参考。     

枯草芽孢杆菌N2-10的基因组测序和比较基因组分析

【背景】枯草芽孢杆菌N2-10是一株具有较强抑菌能力且能产纤维素酶等多种水解酶的革兰氏阳性菌,在发酵饲料中具有较大的应用潜力。【目的】通过获得枯草芽孢杆菌N2-10的全基因组序列信息,进一步解析菌株次级代谢产物合成基因信息,并通过比较基因组学分析菌株N2-10与模式菌株的差异性,为阐明N2-10抑菌和益生机制提供理论基础。【方法】通过二代Illumina NovaSeq联合三代PacBio Sequel测序平台,对菌株N2-10进行全基因组测序,将测序数据进行基因组组装、基因预测与功能注释,并利用比较基因组学分析N2-10与其他菌株的差异。【结果】菌株N2-10基因组大小为4 036 899 bp,GC含量为43.88%;共编码4 163个编码基因,所有编码基因总长度为3594369bp,编码区总长度占基因组总长度的89.04%;含有85个tRNA、10个5S rRNA、10个16S rRNA、10个23S rRNA,以及2个CRISPR-Cas、1个前噬菌体和6个基因岛;在GO (gene ontolog)、COG (clusters of orthologous groups of...     

A tandem repeat gene in a picornavirus.

Three closely related genes for the small genome-linked protein (VPg) of picornaviruses have been identified by sequence analysis as a tandem repeat in the genome of Foot and Mouth Disease Virus (FMDV), strain O1K. This unusual structure was also found in the genome of strain C1O, belonging to a different FMDV serotype. Predicted biochemical properties of the three VPg gene products are in excellent agreement with the data from protein analysis of a heterogeneous VPg population from a third FMDV serotype, strain A10 (1). Taken together, these data indicate that the VPgs from all three genes function equally well in vivo. This is the first report of a tandem repeat gene in a viral genome.     

Intragenomic diversity of the V1 regions of 16S rRNA genes in high-alkaline protease-producing Bacillus clausii spp

Alkaliphilic Bacillus sp. strain KSM-K16, which produces high-alkaline M-protease, was characterized phenotypically, biochemically and genetically. This strain was identified as Bacillus clausii based on the results of taxonomic studies, including sequencing of the 16S rRNA gene and DNA-DNA hybridization. Seven rRNA operons in the genome were identified by pulsed-field gel electrophoresis. Sequencing of cloned 16S rRNA genes revealed two distinct types of variable region V1. Moreover, some cloned 16S rRNA genes in some of the reference strains of B. clausii had a V1 region of yet another type. The B. clausii strains could clearly be divided into at least two subgroups based on the frequencies of the types of cloned V1 sequence. Bacillus sp. strain KSM-K16 was found to be in a different phylogenetic position from other high-alkaline protease-producing strains of B. clausii.     

Highly conserved gene arrangement of the mitochondrial genomes of 23 Plasmodium species

Mitochondrial (mt) genomes from diverse phylogenetic groups vary considerably in size, structure and organization. The genus Plasmodium, the causative agent of malaria, has the smallest mt genome in the form of a tandemly repeated, linear element of 6 kb. The Plasmodium mt genome encodes only three protein genes (cox1, cox3 and cob) and large- and small-subunit ribosomal RNA (rRNA) genes, which are highly fragmented with 19 identified rRNA pieces. The complete mt genome sequences of 21 Plasmodium species have been published but a thorough investigation of the arrangement of rRNA gene fragments has been undertaken for only Plasmodium falciparum, the human malaria parasite. In this study, we determined the arrangement of mt rRNA gene fragments in 23 Plasmodium species, including two newly determined mt genome sequences from P. gallinaceum and P. vinckei vinckei, as well as Leucocytozoon caulleryi, an outgroup of Plasmodium. Comparative analysis reveals complete conservation of the arrangement of rRNA gene fragments in the mt genomes of all the 23 Plasmodium species and L. caulleryi. Surveys for a new rRNA gene fragment using hidden Markov models enriched with recent mt genome sequences led us to suggest the mtR-26 sequence as a novel candidate LSU rRNA fragment in the mt genomes of the 24 species. Additionally, we found 22-25 bp-inverted repeat sequences, which may be involved in the generation of lineage-specific mt genome arrangements after divergence from a common ancestor of the genera Eimeria and Plasmodium/Leucocytozoon.     

Genome sequence of Diplorickettsia massiliensis, an emerging Ixodes ricinus-associated human pathogen

Diplorickettsia massiliensis is a gammaproteobacterium in the order Legionellales and an agent of tick-borne infection. We sequenced the genome from strain 20B, isolated from an Ixodes ricinus tick. The genome consists of a 1,727,973-bp chromosome but no plasmid and includes 2,269 protein-coding genes and 42 RNA genes, including 3 rRNA genes.     

The gene coding for small ribosomal subunit RNA in the basidiomycete Ustilago maydis contains a group I intron.

The nucleotide sequence of the gene coding for small ribosomal subunit RNA in the basidiomycete Ustilago maydis was determined. It revealed the presence of a group I intron with a length of 411 nucleotides. This is the third occurrence of such an intron discovered in a small subunit rRNA gene encoded by a eukaryotic nuclear genome. The other two occurrences are in Pneumocystis carinii, a fungus of uncertain taxonomic status, and Ankistrodesmus stipitatus, a green alga. The nucleotides of the conserved core structure of 101 group I intron sequences present in different genes and genome types were aligned and their evolutionary relatedness was examined. This revealed a cluster including all group I introns hitherto found in eukaryotic nuclear genes coding for small and large subunit rRNAs. A secondary structure model was designed for the area of the Ustilago maydis small ribosomal subunit RNA precursor where the intron is situated. It shows that the internal guide sequence pairing with the intron boundaries fits between two helices of the small subunit rRNA, and that minimal rearrangement of base pairs suffices to achieve the definitive secondary structure of the 18S rRNA upon splicing.