Phylogeny of members of the <Emphasis Type="Italic">Frankia</Emphasis> genus based on <Emphasis Type="Italic">gyr</Emphasis>B, <Emphasis Type="Italic">nif</Emphasis>H and <Emphasis Type="Italic">gln</Emphasis>II sequences

To construct an evolutionary hypothesis for the genus Frankia, gyrB (encoding gyrase B), nifH (encoding nitrogenase reductase) and glnII (encoding glutamine synthetase II) gene sequences were considered for 38 strains. The overall clustering pattern among Frankia strains based on the three analyzed sequences varied among themselves and with the previously established 16S rRNA gene phylogeny and they did not reliably reflect clear evolution of the four discerned Frankia clusters (1, 2, 3 and 4). Based on concatenated gyrB, nifH and glnII, robust phylogenetic trees were observed with the three treeing methods (Maximum Likelihood, Parsimony and Neighbor-Joining) and supported by strong bootstrap and posterior probability values (>75%) for overall branching. Cluster 4 (non-infective and/or non-nitrogen-fixing Frankia) was positioned at a deeper branch followed by cluster 3 (Rhamnaceae and Elaeagnaceae infective Frankia), while cluster 2 represents uncultured Frankia microsymbionts of the Coriariaceae, Datiscaceae, Rosaceae and of Ceanothus sp. (Rhamnaceae); Cluster 1 (Betulaceae, Myricaceae and Casuarinaceae infective Frankia) appears to have diverged more recently. The present study demonstrates the utility of phylogenetic analyses based upon concatenated gyrB, nifH and glnII sequences to help resolve previously unresolved or poorly resolved nodes and will aid in describing species among the genus Frankia.     

Lectin genes in the <Emphasis Type="Italic">Frankia alni</Emphasis> genome

Frankia alni strain ACN14a’s genome was scanned for the presence of determinants involved in interactions with its host plant, Alnus spp. One such determinant type is lectin, proteins that bind specifically to sugar motifs. The genome of F. alni was found to contain 7 such lectin-coding genes, five of which were of the ricinB-type. The proteins coded by these genes contain either only the lectin domain, or also a heat shock protein or a serine-threonine kinase domain upstream. These lectins were found to have several homologs in Streptomyces spp., and a few in other bacterial genomes among which none in Frankia EAN1pec and CcI3 and two in strain EUN1f. One of these F. alni genes, FRAAL0616, was cloned in E. coli, fused with a reporter gene yielding a fusion protein that was found to bind to both root hairs and to bacterial hyphae. This protein was also found to modify the dynamics of nodule formation in A. glutinosa, resulting in a higher number of nodules per root. Its role could thus be to permit binding of microbial cells to root hairs and help symbiosis to occur under conditions of low Frankia cell counts such as in pioneer situations.     

A non-sense mutation in the putative anti-mutator gene <Emphasis Type="Italic">ada/alkA</Emphasis> of <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis> and <Emphasis Type="Italic">M. bovis</Emphasis> isolates suggests convergent evolution

Background  

Previous studies have suggested that variations in DNA repair genes of W-Beijing strains may have led to transient mutator phenotypes which in turn may have contributed to host adaptation of this strain family. Single nucleotide polymorphism (SNP) in the DNA repair gene mutT1 was identified in MDR-prone strains from the Central African Republic. A Mycobacteriumtuberculosis H37Rv mutant inactivated in two DNA repair genes, namely ada/alkA and ogt, was shown to display a hypermutator phenotype. We then looked for polymorphisms in these genes in Central African Republic strains (CAR).     

The complete mitogenome and phylogenetic analysis of <Emphasis Type="Italic">Bombyx mandarina</Emphasis> strain Qingzhou

The complete nucleotide sequence of the mitogenome of Bombyx mandarina strain Qingzhou was determined. The circular genome is 15,717 bp long and has the typical gene organization and order of lepidopteran mitogenomes. All protein-coding sequences are initiated with a typical ATN codon, except the COI gene, which has a 4-bp TTAG putative initiator codon. Eleven of the 13 protein-coding gene have a complete termination codon (all TAA), but the remaining two genes terminate with incomplete codons. All transfer RNAs (tRNAs) have a clover-leaf structure typical of the mitochondrial tRNAs, and some of them have a mismatch in the four-stem-and-loop structure. The length of the A + T rich region of B. mandarina strain Qingzhou is 495 bp, shorter than that of B. mandarina strain Tsukuba (747 bp) but similar to that of Bombyx mori. Phylogenetic analysis based on the whole mitochondrial genome sequences of the available sequenced species (B. mori strains C-108, Aojuku, Backokjam, and Xiafang, B. mandarina strains Tsukuba, Ankang, and Qingzhou, and Antheraea pernyi) shows the origin of the domesticated silkmoth B. mori to be the Chinese B. mandarina. Nuclear mitochondrial pseudogene sequences were detected in the nuclear genome of B. mori with the MEGA BLAST search program. A phylogenetic analysis of these nuclear mitochondrial pseudogene sequences suggests that B. mori was domesticated independently in different areas and periods.     

IS<Emphasis Type="Italic">Rm31</Emphasis>, a new insertion sequence of the IS<Emphasis Type="Italic">66</Emphasis> family in<Emphasis Type="Italic"> Sinorhizobium meliloti</Emphasis>

Sinorhizobium meliloti natural populations show a high level of genetic polymorphism possibly due to the presence of mobile genetic elements such as insertion sequences (IS), transposons, and bacterial mobile introns. The analysis of the DNA sequence polymorphism of the nod region of S. meliloti pSymA megaplasmid in an Italian isolate led to the discovery of a new insertion sequence, ISRm31. ISRm31 is 2,803 bp long and has 22-bp-long terminal inverted repeat sequences, 8-bp direct repeat sequences generated by transposition, and three ORFs (A, B, C) coding for proteins of 124, 115, and 541 amino acids, respectively. ORF A and ORF C are significantly similar to members of the transposase family. Amino acid and nucleotide sequences indicate that ISRm31 is a member of the IS66 family. ISRm31 sequences were found in 30.5% of the Italian strains analyzed, and were also present in several collection strains of the Rhizobiaceae family, including S. meliloti strain 1021. Alignment of targets sites in the genome of strains carrying ISRm31 suggested that ISRm31 inserts randomly into S. meliloti genomes. Moreover, analysis of ISRm31 insertion sites revealed DNA sequences not present in the recently sequenced S. meliloti strain 1021 genome. In fact, ISRm31 was in some cases linked to DNA fragments homologous to sequences found in other rhizobia species.     

Molecular diversity of <Emphasis Type="Italic">Frankia</Emphasis> from root nodules of <Emphasis Type="Italic">Hippophae salicifolia</Emphasis> D.Don found in Sikkim

Molecular diversity of Frankia was assessed directly from the root nodules of Hippophae salicifolia naturally occurring in North Sikkim. Amplicon restriction patterns (ARPs) were developed by digesting 16S-ITS-23S amplicons with RsaI. Three ARPs were detected, showing diversity among strains of Frankia that nodulate Hippophae. This was confirmed by sequencing one amplicon each for the three ARPs. Therefore, ARP can be used as a tool for screening amplicons for nucleotide sequencing.     

Evolution of nematode-resistant <Emphasis Type="Italic">Mi</Emphasis>-<Emphasis Type="Italic">1</Emphasis> gene homologs in three species of <Emphasis Type="Italic">Solanum</Emphasis>

Plants have evolved several defense mechanisms, including resistance genes. Resistance to the root-knot nematode Meloidogyne incognita has been found in wild plant species. The molecular basis for this resistance has been best studied in the wild tomato Solanum peruvianum and it is based on a single dominant gene, Mi-1.2, which is found in a cluster of seven genes. This nematode attacks fiercely several crops, including potatoes. The genomic arrangement, number of copies, function and evolution of Mi-1 homologs in potatoes remain unknown. In this study, we analyzed partial genome sequences of the cultivated potato species S. tuberosum and S. phureja and identified 59 Mi-1 homologs. Mi-1 homologs in S. tuberosum seem to be arranged in clusters and located on chromosome 6 of the potato genome. Previous studies have suggested that Mi-1 genes in tomato evolved rapidly by frequent sequence exchanges among gene copies within the same cluster, losing orthologous relationships. In contrast, Mi-1 homologs from cultivated potato species (S. tuberosum and S. phureja) seem to have evolved by a birth-and-death process, in which genes evolve mostly by mutations and interallelic recombinations in addition to sequence exchanges.     

The <Emphasis Type="Italic">Salmonella enterica</Emphasis> Pan-genome

Salmonella enterica is divided into four subspecies containing a large number of different serovars, several of which are important zoonotic pathogens and some show a high degree of host specificity or host preference. We compare 45 sequenced S. enterica genomes that are publicly available (22 complete and 23 draft genome sequences). Of these, 35 were found to be of sufficiently good quality to allow a detailed analysis, along with two Escherichia coli strains (K-12 substr. DH10B and the avian pathogenic E. coli (APEC O1) strain). All genomes were subjected to standardized gene finding, and the core and pan-genome of Salmonella were estimated to be around 2,800 and 10,000 gene families, respectively. The constructed pan-genomic dendrograms suggest that gene content is often, but not uniformly correlated to serotype. Any given Salmonella strain has a large stable core, whilst there is an abundance of accessory genes, including the Salmonella pathogenicity islands (SPIs), transposable elements, phages, and plasmid DNA. We visualize conservation in the genomes in relation to chromosomal location and DNA structural features and find that variation in gene content is localized in a selection of variable genomic regions or islands. These include the SPIs but also encompass phage insertion sites and transposable elements. The islands were typically well conserved in several, but not all, isolates—a difference which may have implications in, e.g., host specificity.     

Analysis of insect toxin complex gene variability of <Emphasis Type="Italic">Yersinia pestis</Emphasis> and <Emphasis Type="Italic">Yersinia pseudotuberculosis</Emphasis> strains

The nucleotide sequences of the Tc’s insect toxin complex genes have been analyzed in 18 natural strains of the main and non-main subspecies of Yersinia pestis isolated in different natural foci in the Russian Federation, as well as neighboring and more remote countries, as compared to the data on Y. pestis and Y. pseudotuberculosis strains stored in the NCBI GenBank database. The nucleotide sequences of these genes in plague agent strains have been found to be highly conserved, in contrast to those of the pseudotuberculosis agent. The sequences of two genes, tcaC and tccC2, have been found to be almost identical in Y. pestis strains, whereas other three genes (tcaA, tcaB, and tccC1) contain a few mutations, which, however, are not common for all strains of the plague agent. Exceptions are only strains of the Y. pestis biovar orientalis, whose tcaB gene is in a nonfunctional state due to a nucleotide deletion. The results suggest that the formation of the species Y. pestis as an agent of a natural focal infection with a transmissive mechanism has not resulted in degradation of the Tc’s complex genes. Instead, these genes are likely to have been altered as the plague agent have been adapting to the new environment.     

Functional difference between <Emphasis Type="Italic">Sinorhizobium meliloti</Emphasis> NifA and <Emphasis Type="Italic">Enterobacter cloacae</Emphasis> NifA

The nifA gene is an important regulatory gene and its product, NifA protein, regulates the expression of many nif genes involved in the nitrogen fixation process. We introduced multiple copies of the constitutively expressed Sinorhizobium meliloti (Sm) or Enterobacter cloacae (Ec) nifA gene into both the nifA mutant strain SmY and the wild-type strain Sm1021. Root nodules produced by SmY containing a constitutively expressed Sm nifA gene were capable of fixing nitrogen, while nodules produced by SmY containing the Ec nifA gene remained unable to fix nitrogen, as is the case for SmY itself. However, transfer of an additional Sm nifA gene into Sm1021 improved the nitrogen-fixing efficiency of root nodules to a greater extent than that observed upon transfer of the Ec nifA gene into Sm1021. Comparative analysis of amino acid sequences between Sm NifA and Ec NifA showed that the N-terminal domain was the least similar, but this domain is indispensable for complementation of the Fix^- phenotype of SmY by Sm NifA. We conclude that more than one domain is involved in determining functional differences between Sm NifA and Ec NifA.     

<Emphasis Type="Italic">Robinia pseudoacacia</Emphasis> in Poland and Japan is nodulated by <Emphasis Type="Italic">Mesorhizobium amorphae</Emphasis> strains

Robinia pseudoacacia microsymbionts from plants growing in Poland and Japan were evaluated for phylogeny and taxonomic position by genomic approach. Based on the comparative analyses of atpD (368 bp) and dnaK (573 bp) gene sequences as well as 16S rDNA restriction analysis (RFLP-16S rDNA), R. pseudoacacia microsymbionts were identified as Mesorhizobium strains. In dnaK and atpD gene phylograms R. pseudoacacia nodulators formed robust, monophyletic clusters with Mesorhizobium species with the nucleotide sequence similarity of 91–98% and 90–98%, respectively. The classification of R. pseudoacacia rhizobia to the genus Mesorhizobium was also supported by amplified 16S rDNA restriction analysis. The studied bacteria formed common clusters with Mesorhizobium species, and their DNA patterns were identical or nearly identical to Mesorhizobium genus strains. When DNA-DNA hybridization was performed, the total DNA of the representative R. pseudoacacia rhizobia exhibited 51–75% relatedness to DNA of Mesorhizobium amorphae ICMP15022 strain and below 41% to DNA of other Mesorhizobium species. These results showed that R. pseudoacacia and M. amorphae belong to the same genomospecies. The G+C content of DNA of R. pseudoacacia two microsymbionts was 59.7 and 60.6 mol% compared to 61–64 mol% across M. amorphae strains.     

Shared molecular characteristics of successfully transformed mitochondrial genomes in <Emphasis Type="Italic">Chlamydomonas</Emphasis><Emphasis Type="Italic">reinhardtii</Emphasis>

Three types of respiratory deficient mitochondrial strains have been reported in Chlamydomonas reinhardtii: a deficiency due to (i) two base substitutions causing an amino acid change in the apocytochrome b (COB) gene (i.e., strain named dum-15), (ii) one base deletion in the COXI gene (dum-19), or (iii) a large deletion extending from the left terminus of the genome to somewhere in the COB gene (dum-1, -14, and -16). We found that these respiratory deficient strains of C. reinhardtii can be divided into two groups: strains that are constantly transformable and those could not be transformed in our experiments. All transformable mitochondrial strains were limited to the type that has a large deletion in the left arm of the genome. For these mitochondria, transformation was successful not only with purified intact mitochondrial genomes but also with DNA-constructs containing the compensating regions. In comparison, mitochondria of all the non-transformable strains have both of their genome termini intact, leading us to speculate that mitochondria lacking their left genome terminus have unstable genomes and might have a higher potential for recombination. Analysis of mitochondrial gene organization in the resulting respiratory active transformants was performed by DNA sequencing and restriction enzyme digestion. Such analysis showed that homologous recombination occurred at various regions between the mitochondrial genome and the artificial DNA-constructs. Further analysis by Southern hybridization showed that the wild-type genome rapidly replaces the respiratory deficient monomer and dimer mitochondrial genomes, while the E. coli vector region of the artificial DNA-construct likely does not remain in the mitochondria.     

Polymorphism of lectin genes in <Emphasis Type="Italic">Lathyrus</Emphasis> plants

The carbohydrate-binding sequences of the lectin genes from spring vetchling Lathyrus vernus (L.) Bernh., marsh vetchling L. palustris (L.), and Gmelin’s vetchling L. gmelinii (Fitsch) (Fabaceae) were determined. Computer-aided analysis revealed substantial differences between nucleotide and predicted amino acid sequences of the lectin gene regions examined in each of the three vetchling species tested. In the phylogenetic trees based on sequence similarity of carbohydrate-biding regions of legume lectins, the sequences examined formed a compact cluster with the lectin genes of the plants belonging to the tribe Fabeae. In each plant, L. vernus, L. palustris, and L. gmelinii, three different lectin-encoding genes were detected. Most of the substitutions were identified within the gene sequence responsible for coding the carbohydrate-binding protein regions. This finding may explain different affinity of these lectins to different carbohydrates, and as a consequence, can affect the plant host specificity upon development of symbiosis with rhizobium bacteria.     

Discrimination among <Emphasis Type="Italic">Bacillus thuringiensis</Emphasis> H serotypes,serovars and strains based on 16S rRNA, <Emphasis Type="Italic">gyrB</Emphasis> and <Emphasis Type="Italic">aroE</Emphasis> gene sequence analyses

Our aim was to investigate the capability of each of three genes, 16S rRNA, gyrB and aroE, to discriminate, first, among Bacillus thuringiensis H serotypes; second, among B. thuringiensis serovars from the same H serotype; and third, among B. thuringiensis strains from the same serovar. The 16S rRNA, gyrB and aroE genes were amplified from 21 B. thuringiensis H serotypes and their nucleotide sequences determined. Additional strains from four B. cereus sensu lato species were included for comparison purposes. These sequences were pair-wise compared and phylogenetic relationships were revealed. Each of the three genes under study could discriminate among B. thuringiensis H serotypes. The gyrB and aroE genes showed a discriminatory power among B. thuringiensis H serotypes up to nine fold greater than that of the 16S rRNA gene. The gyrB gene was retained for subsequent analyses to discriminate B. thuringiensis serovars from the same H serotype and to discriminate strains from same serovar. A total of 42 B. thuringiensis strains, which encompassed 25 serovars from 12 H serotypes, were analyzed. The gyrB gene nucleotide sequences were different enough as to be sufficient to discriminate among B. thuringiensis serovars from the same H serotype and among B. thuringiensis strains from the same serovar. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

BAHD superfamily of acyl-CoA dependent acyltransferases in <Emphasis Type="Italic">Populus</Emphasis> and <Emphasis Type="Italic">Arabidopsis</Emphasis>: bioinformatics and gene expression

Plant acyl-CoA dependent acyltransferases constitute a large specific protein superfamily, named BAHD. Using the conserved sequence motifs of BAHD members, we searched the genome sequences of Populus and Arabidopsis, and identified, respectively, 94- and 61-putative genes. Subsequently, we analyzed the phylogeny, gene structure, and chromosomal distribution of BAHD members of both species; then, we profiled expression patterns of BAHD genes by “in silico” northern- and microarray-analyses based on public databases, and by RT-PCR. While our genomic- and bioinformatic- analyses provided full sets of BAHD superfamily genes, and cleaned up a few existing annotation errors, importantly it led to our recognizing several unique Arabidopsis BAHD genes that inversely overlapped with their neighboring genes on the genome, and disclosing a potential natural anti-sense regulation for gene expressions. Systemic gene-expression profiling of BAHD members revealed distinct tissue-specific/preferential expression patterns, indicating their diverse biological functions. Our study affords a strong knowledge base for understanding BAHD members’ evolutionary relationships and gene functions implicated in plant growth, development and metabolism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Nucleotide substitutions in <Emphasis Type="Italic">rolC</Emphasis> and <Emphasis Type="Italic">nptII</Emphasis> gene sequences during long-term cultivation of <Emphasis Type="Italic">Panax ginseng</Emphasis> cell cultures

It has been shown previously that the rolC gene from Agrobacterium tumefaciens gene was stably and highly expressed in 15-year-old Panax ginseng transgenic cell cultures. In the present report, we analyze in detail the nucleotide composition of the rolC and nptII (neomycin phosphotransferase) genes, which is the selective marker used for transgenic cell cultures of P. ginseng. It has been established that the nucleotide sequences of the rolC and nptII genes underwent mutagenesis during cultivation. Particularly, 1–4 nucleotide substitutions were found per sequence in the 540 and 798 bp segments of the complete rolC and nptII genes, respectively. Approximately half of these nucleotide substitutions caused changes in the structure of the predicted gene product. In addition, we attempted to determine the rate of accumulation of these changes by comparison of DNA extracted from P. ginseng cell cultures from 1995 to 2007. It was observed that the frequency of nucleotide substitutions for the rolC and nptII genes in 1995 was 1.21 ± 0.02 per 1,000 nucleotides analyzed, while in 2007, the nucleotide substitutions significantly increased (1.37 ± 0.07 per 1,000 nucleotides analyzed). Analyzing the nucleotide substitutions, we found that substitution to G or to C nucleotides significantly increased (in 1.9 times) in the rolC and nptII genes compared with P. ginseng actin gene. Finally, the level of nucleotide substitutions in the rolC gene was 1.1-fold higher when compared with the nptII gene. Thus, for the first time, we have experimentally demonstrated the level of nucleotide substitutions in transferred genes in transgenic plant cell cultures.     

The <Emphasis Type="Italic">Caenorhabditis</Emphasis><Emphasis Type="Italic">briggsae</Emphasis> genome contains active <Emphasis Type="Italic">CbmaT1</Emphasis> and <Emphasis Type="Italic">Tcb1</Emphasis> transposons

The maT clade of transposons is a group of transposable elements intermediate in sequence and predicted protein structure to mariner and Tc transposons, with a distribution thus far limited to a few invertebrate species. We present evidence, based on searches of publicly available databases, that the nematode Caenorhabditis briggsae has several maT-like transposons, which we have designated as CbmaT elements, dispersed throughout its genome. We also describe two additional transposon sequences that probably share their evolutionary history with the CbmaT transposons. One resembles a fold back variant of a CbmaT element, with long (380-bp) inverted terminal repeats (ITRs) that show a high degree (71%) of identity to CbmaT1. The other, which shares only the 26-bp ITR sequences with one of the CbmaT variants, is present in eight nearly identical copies, but does not have a transposase gene and may therefore be cross mobilised by a CbmaT transposase. Using PCR-based mobility assays, we show that CbmaT1 transposons are capable of excising from the C. briggsae genome. CbmaT1 excised approximately 500 times less frequently than Tcb1 in the reference strain AF16, but both CbmaT1 and Tcb1 excised at extremely high frequencies in the HK105 strain. The HK105 strain also exhibited a high frequency of spontaneous induction of unc-22 mutants, suggesting that it may be a mutator strain of C. briggsae.     

Application of <Emphasis Type="Italic">gyrB</Emphasis> and <Emphasis Type="Italic">parE</Emphasis> sequence similarity analyses for differentiation of species within the genus <Emphasis Type="Italic">Geobacillus</Emphasis>

The primary structures of the genes encoding the β-subunits of a type II topoisomerase (gyrase, gyrB) and a type IV topoisomerase (parE) were determined for 15 strains of thermophilic bacteria of the genus Geobacillus. The obtained sequences were used for analysis of the phylogenetic similarity between members of this genus. Comparison of the phylogenetic trees of geobacilli constructed on the basis of the 16S rRNA, gyrB, and parE gene sequences demonstrated that the level of genetic distance between the sequences of the genes encoding the β-subunits of type II topoisomerases significantly exceeded the values obtained by comparative analysis of the 16S rRNA gene sequences of Geobacillus strains. It was shown that, unlike the 16S rRNA gene analysis, comparative analysis of the gyrB and parE gene sequences provided a more precise determination of the phylogenetic position of bacteria at the species level. The data obtained suggest the possibility of using the genes encoding the β-subunits of type II topoisomerases as phylogenetic markers for determination of the species structure of geobacilli.