Analysis of distribution indicates diverse functions of simple sequence repeats in Mycoplasma genomes

Reevaluation of Parasynechococcus-like Strains and Genomic Analysis of Their Microsatellites and Compound Microsatellites

Morphologically similar to Synechococcus, a large number of Parasynechococcus strains were misclassified, resulting in extreme underestimation of their genetic diversity. In this study, 80 Synechococcus-like strains were reevaluated using a combination of 16S rRNA phylogeny and genomic approach, identifying 54 strains as Parasynechococcus-like strains and showing considerably intragenus genetic divergence among the subclades identified. Further, bioinformatics analysis disclosed diversified patterns of distribution, abundance, density, and diversity of microsatellites (SSRs) and compound microsatellites (CSSRs) in genomes of these Parasynechococcus-like strains. Variations of SSRs and CSSRs were observed amongst phylotypes and subclades. Both SSRs and CSSRs were in particular unequally distributed among genomes. Dinucleotide SSRs were the most widespread, while the genomes showed two patterns in the second most abundant repeat type (mononucleotide or trinucleotide SSRs). Both SSRs and CSSRs were predominantly observed in coding regions. These two types of microsatellites showed positive correlation with genome size (p < 0.01) but negative correlation with GC content (p < 0.05). Additionally, the motif (A)n, (AG)n and (AGC)n was a major one in the corresponding category. Meanwhile, distinctive motifs of CSSRs were found in 39 genomes. This study characterizes SSRs and CSSRs in genomes of Parasynechococcus-like strains and will be useful as a prerequisite for future studies regarding their distribution, function, and evolution. Moreover, the identified SSRs may facilitate fast acclimation of Parasynechococcus-like strains to fluctuating environments and contribute to the extensive distribution of Parasynechococcus species in global marine environments.

Methylation Motifs in Promoter Sequences May Contribute to the Maintenance of a Conserved m5C Methyltransferase in Helicobacter pylori

DNA methylomes of Helicobacter pylori strains are complex due to the large number of DNA methyltransferases (MTases) they possess. H. pylori J99 M.Hpy99III is a 5-methylcytosine (^m5C) MTase that converts GCGC motifs to G^m5CGC. Homologs of M.Hpy99III are found in essentially all H. pylori strains. Most of these homologs are orphan MTases that lack a cognate restriction endonuclease, and their retention in H. pylori strains suggest they have roles in gene regulation. To address this hypothesis, green fluorescent protein (GFP) reporter genes were constructed with six putative promoters that had a GCGC motif in the extended −10 region, and the expression of the reporter genes was compared in wild-type H. pylori G27 and a mutant lacking the M.Hpy99III homolog (M.HpyGIII). The expression of three of the GFP reporter genes was decreased significantly in the mutant lacking M.HpyGIII. In addition, the growth rate of the H. pylori G27 mutant lacking M.HpyGIII was reduced markedly compared to that of the wild type. These findings suggest that the methylation of the GCGC motif in many H. pylori GCGC-containing promoters is required for the robust expression of genes controlled by these promoters, which may account for the universal retention of M.Hpy99III homologs in H. pylori strains.

Genome-Wide Investigation and Analysis of Microsatellites and Compound Microsatellites in Leptolyngbya-like Species, Cyanobacteria

Microsatellites (simple sequence repeats, SSRs) are ubiquitously distributed in almost all known genomes. Here, the first investigation was designed to examine the SSRs and compound microsatellites (CSSRs) in genomes of Leptolyngbya-like strains. The results disclosed diversified patterns of distribution, abundance, density, and diversity of SSRs and CSSRs in genomes, indicating that they may be subject to rapid evolutionary change. The numbers of SSRs and CSSRs were extremely unevenly distributed among genomes, ranging from 11,086 to 24,000 and from 580 to 1865, respectively. Dinucleotide SSRs were the most abundant category in 31 genomes, while the other 15 genomes followed the pattern: mono- > di- > trinucleotide SSRs. The patterns related to SSRs and CSSRs showed differences among phylogenetic groups. Both SSRs and CSSRs were overwhelmingly distributed in coding regions. The numbers of SSRs and CSSRs were significantly positively correlated with genome size (p < 0.01) and negatively correlated with GC content (p < 0.05). Moreover, the motif (A/C)_n and (AG)_n was predominant in mononucleotide and dinucleotide SSRs, and unique motifs of CSSRs were identified in 39 genomes. This study provides the first insight into SSRs and CSSRs in genomes of Leptolyngbya-like strains and will be useful to understanding their distribution, predicting their function, and tracking their evolution. Additionally, the identified SSRs may provide an evolutionary advantage of fast adaptation to environmental changes and may play an important role in the cosmopolitan distribution of Leptolyngbya strains to globally diverse niches.

An empirical analysis of mtSSRs: could microsatellite distribution patterns explain the evolution of mitogenomes in plants?

Microsatellites (SSRs) are tandem repeat sequences in eukaryote genomes, including plant cytoplasmic genomes. The mitochondrial genome (mtDNA) has been shown to vary in size, number, and distribution of SSRs among different plant groups. Thus, SSRs contribute with genomic diversity in mtDNAs. However, the abundance, distribution, and evolutionary significance of SSRs in mtDNA from a wide range of algae and plants have not been explored. In this study, the mtDNAs of 204 plant and algal species were investigated related to the presence of SSRs. The number of SSRs was positively correlated with genome size. Its distribution is dependent on plant and algal groups analyzed, although the cluster analysis indicates the conservation of some common motifs in algal and terrestrial plants that reflect common ancestry of groups. Many SSRs in coding and non-coding regions can be useful for molecular markers. Moreover, mitochondrial SSRs are highly abundant, representing an important source for natural or induced genetic variation, i.e., for biotechnological approaches that can modulate mtDNA gene regulation. Thus, this comparative study increases the understanding of the plant and algal SSR evolution and brings perspectives for further studies.

The influence of regulatory elements on Mycoplasma hyopneumoniae 7448 transcriptional response during oxidative stress and heat shock

BACKGROUND

The comprehension of genome organization and gene modulation is essential for understanding pathogens' infection mechanisms. Mycoplasma hyopneumoniae 7448 genome is organized in transcriptional units (TUs), which are flanked by regulatory elements such as putative promoters, terminators and repetitive sequences. Yet the relationship between the presence of these elements and bacterial responses during stress conditions remains unclear. Thus, in this study, in silico and RT-qPCR analyses were associated to determine the effect of regulatory elements in gene expression regulation upon heat shock and oxidative stress conditions.

METHODS AND RESULTS

Thirteen TU's organizational profiles were found based on promoters and terminators distribution. Differential expression in genes sharing the same TUs was observed, suggesting the activity of internal regulatory elements. Moreover, 88.8% of tested genes were differentially expressed under oxidative stress in comparison to the control condition, being 81.3% of them surrounded by their own regulatory elements. Similarly, under heat shock, 44.4% of the genes showed regulation when compared to control condition, being 75.0% of them surrounded by their own regulatory elements.

CONCLUSIONS

Altogether, this data suggests the activity of internal regulatory elements in gene modulation of M. hyopneumoniae 7448 transcription.

Relatively semi-conservative replication and a folded slippage model for short tandem repeats

Background

The ubiquitous presence of short tandem repeats (STRs) in virtually all genomes implicates their functional relevance, while a widely-accepted definition of STR is yet to be established. Previous studies majorly focus on relatively longer STRs, while shorter repeats were generally excluded. Herein, we have adopted a more generous criteria to define shorter repeats, which has led to the definition of a much larger number of STRs that lack prior analysis. Using this definition, we analyzed the short repeats in 55 randomly selected segments in 55 randomly selected genomic sequences from a fairly wide range of species covering animals, plants, fungi, protozoa, bacteria, archaea and viruses.

Results

Our analysis reveals a high percentage of short repeats in all 55 randomly selected segments, indicating that the universal presence of high-content short repeats could be a common characteristic of genomes across all biological kingdoms. Therefore, it is reasonable to assume a mechanism for continuous production of repeats that can make the replicating process relatively semi-conservative. We have proposed a folded replication slippage model that considers the geometric space of nucleotides and hydrogen bond stability to explain the mechanism more explicitly, with improving the existing straight-line slippage model. The folded slippage model can explain the expansion and contraction of mono- to hexa- nucleotide repeats with proper folding angles. Analysis of external forces in the folding template strands also suggests that expansion exists more commonly than contraction in the short tandem repeats.

Conclusion

The folded replication slippage model provides a reasonable explanation for the continuous occurrences of simple sequence repeats in genomes. This model also contributes to the explanation of STR-to-genome evolution and is an alternative model that complements semi-conservative replication.

Structures and stability of simple DNA repeats from bacteria

DNA is a fundamentally important molecule for all cellular organisms due to its biological role as the store of hereditary, genetic information. On the one hand, genomic DNA is very stable, both in chemical and biological contexts, and this assists its genetic functions. On the other hand, it is also a dynamic molecule, and constant changes in its structure and sequence drive many biological processes, including adaptation and evolution of organisms. DNA genomes contain significant amounts of repetitive sequences, which have divergent functions in the complex processes that involve DNA, including replication, recombination, repair, and transcription. Through their involvement in these processes, repetitive DNA sequences influence the genetic instability and evolution of DNA molecules and they are located non-randomly in all genomes. Mechanisms that influence such genetic instability have been studied in many organisms, including within human genomes where they are linked to various human diseases. Here, we review our understanding of short, simple DNA repeats across a diverse range of bacteria, comparing the prevalence of repetitive DNA sequences in different genomes. We describe the range of DNA structures that have been observed in such repeats, focusing on their propensity to form local, non-B-DNA structures. Finally, we discuss the biological significance of such unusual DNA structures and relate this to studies where the impacts of DNA metabolism on genetic stability are linked to human diseases. Overall, we show that simple DNA repeats in bacteria serve as excellent and tractable experimental models for biochemical studies of their cellular functions and influences.

Mutational and Selective Processes Involved in Evolution during Bacterial Range Expansions

Bacterial populations have been shown to accumulate deleterious mutations during spatial expansions that overall decrease their fitness and ability to grow. However, it is unclear if and how they can respond to selection in face of this mutation load. We examine here if artificial selection can counteract the negative effects of range expansions. We examined the molecular evolution of 20 mutator lines selected for fast expansions (SEL) and compared them to 20 other mutator lines freely expanding without artificial selection (CONTROL). We find that the colony size of all 20 SEL lines have increased relative to the ancestral lines, unlike CONTROL lines, showing that enough beneficial mutations are produced during spatial expansions to counteract the negative effect of expansion load. Importantly, SEL and CONTROL lines have similar numbers of mutations indicating that they evolved for the same number of generations and that increased fitness is not due to a purging of deleterious mutations. We find that loss of function mutations better explain the increased colony size of SEL lines than nonsynonymous mutations or a combination of the two. Interestingly, most loss of function mutations are found in simple sequence repeats (SSRs) located in genes involved in gene regulation and gene expression. We postulate that such potentially reversible mutations could play a major role in the rapid adaptation of bacteria to changing environmental conditions by shutting down expensive genes and adjusting gene expression.

Comparative genomics in phytopathogenic prokaryotes reveals the higher relative abundance and density of long-SSRs in the smallest prokaryotic genome

Frequency and distribution of long-SSRs were studied in 18 phytopathogenic prokaryotes. Higher relative abundance of the long-SSRs was observed in phytopathogenic prokaryotes when compared to non-pathogenic control. The frequency of these SSRs was positively correlated with size and GC content of the genomes of phytopathogenic prokaryotes. Interestingly, phytopathogens with higher GC content in the genome were found to posses longer repeat motifs of SSRs, whereas those having lesser GC content were harbouring shorter repeat motifs. Higher abundance of tri- and hexa-nucleotide repeat motifs were the characteristic of actinomycetes, where as higher abundance of mono- and tetra-nucleotide repeats were the characteristic of the mollicutes. The maximum relative abundance and relative density of SSR were found in the smallest genome of host-adapted pathogen Aster yellow, however, length of microsatellite repeat units was the least. On the basis of presence of SSRs in the housekeeping genes, a phylogenetic relationship between these phytopathogenic prokaryotes was deduced and compared with the phylogeny developed based on 16S ribosomal RNA gene.

Microsatellite Diversity, Complexity, and Host Range of Mycobacteriophage Genomes of the Siphoviridae Family

The incidence, distribution, and variation of simple sequence repeats (SSRs) in viruses is instrumental in understanding the functional and evolutionary aspects of repeat sequences. Full-length genome sequences retrieved from NCBI were used for extraction and analysis of repeat sequences using IMEx software. We have also developed two MATLAB-based tools for extraction of gene locations from GenBank in tabular format and simulation of this data with SSR incidence data. Present study encompassing 147 Mycobacteriophage genomes revealed 25,284 SSRs and 1,127 compound SSRs (cSSRs) through IMEx. Mono- to hexa-nucleotide motifs were present. The SSR count per genome ranged from 78 (M100) to 342 (M58) while cSSRs incidence ranged from 1 (M138) to 17 (M28, M73). Though cSSRs were present in all the genomes, their frequency and SSR to cSSR conversion percentage varied from 1.08 (M138 with 93 SSRs) to 8.33 (M116 with 96 SSRs). In terms of localization, the SSRs were predominantly localized to coding regions (∼78%). Interestingly, genomes of around 50 kb contained a similar number of SSRs/cSSRs to that in a 110 kb genome, suggesting functional relevance for SSRs which was substantiated by variation in motif constitution between species with different host range. The three species with broad host range (M97, M100, M116) have around 90% of their mono-nucleotide repeat motifs composed of G or C and only M16 has both A and T mononucleotide motifs. Around 20% of the di-nucleotide repeat motifs in the genomes exhibiting a broad host range were CT/TC, which were either absent or represented to a much lesser extent in the other genomes.

In silico simulations of occurrence of transcription factor binding sites in bacterial genomes

Background

Interactions between transcription factors and their specific binding sites are a key component of regulation of gene expression. Until recently, it was generally assumed that most bacterial transcription factor binding sites are located at or near promoters. However, several recent works utilizing high-throughput technology to detect transcription factor binding sites in bacterial genomes found a large number of binding sites in unexpected locations, particularly inside genes, as opposed to known or expected promoter regions. While some of these intragenic binding sites likely have regulatory functions, an alternative scenario is that many of these binding sites arise by chance in the absence of selective constraints. The latter possibility was supported by in silico simulations for σ⁵⁴ binding sites in Salmonella.

Results

In this work, we extend these simulations to more than forty transcription factors from E. coli and other bacteria. The results suggest that binding sites for all analyzed transcription factors are likely to arise throughout the genome by random genetic drift and many transcription factor binding sites found in genomes may not have specific regulatory functions. In addition, when comparing observed and expected patterns of occurrence of binding sites in genomes, we observed distinct differences among different transcription factors.

Conclusions

We speculate that transcription factor binding sites randomly occurring throughout the genome could be beneficial in promoting emergence of new regulatory interactions and thus facilitating evolution of gene regulatory networks.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1381-8) contains supplementary material, which is available to authorized users.

Comparative Analysis of Mycoplasma gallisepticum vlhA Promoters

Mycoplasma gallisepticum is an intracellular parasite affecting respiratory tract of poultry that belongs to class Mollicutes. M. gallisepticum features numerous variable lipoprotein hemagglutinin genes (vlhA) that play a role in immune escape. The vlhA promoters have a set of distinct properties in comparison to promoters of the other genes. The vlhA promoters carry a variable GAA repeats region at approximately 40 nts upstream of transcription start site. The promoters have been considered active only in the presence of exactly 12 GAA repeats. The mechanisms of vlhA expression regulation and GAA number variation are not described. Here we tried to understand these mechanisms using different computational methods. We conducted a comparative analysis among several M. gallisepticum strains. Nucleotide sequences analysis showed the presence of highly conserved regions flanking repeated trinucleotides that are not linked to GAA number variation. VlhA genes with 12 GAA repeats and their orthologs in 12 M. gallisepticum strains are more conserved than other vlhA genes and have narrower GAA number distribution. We conducted comparative analysis of physicochemical profiles of M. gallisepticum vlhA and sigma-70 promoters. Stress-induced duplex destabilization (SIDD) profiles showed that sigma-70 group is characterized by the common to prokaryotic promoters sharp maxima while vlhA promoters are hardly destabilized with the region between GAA repeats and transcription start site having zero opening probability. Electrostatic potential profiles of vlhA promoters indicate the presence of the distinct patterns that appear to govern initial stages of specific DNA-protein recognition. Open state dynamics profiles of vlhA demonstrate the pattern that might facilitate transcription bubble formation. Obtained data could be the basis for experimental identification of mechanisms of phase variation in M. gallisepticum.

Alternative transcriptional regulation in genome-reduced bacteria

Transcription is a core process of bacterial physiology, and as such it must be tightly controlled, so that bacterial cells maintain steady levels of each RNA molecule in homeostasis and modify them in response to perturbations. The major regulators of transcription in bacteria (and in eukaryotes) are transcription factors. However, in genome-reduced bacteria, the limited number of these proteins is insufficient to explain the variety of responses shown upon changes in their environment. Thus, alternative regulators may play a central role in orchestrating RNA levels in these microorganisms. These alternative mechanisms rely on intrinsic features within DNA and RNA molecules, suggesting they are ancestral mechanisms shared among bacteria that could have an increased relevance on transcriptional regulation in minimal cells. In this review, we summarize the alternative elements that can regulate transcript abundance in genome-reduced bacteria and how they contribute to the RNA homeostasis at different levels.

Characterization of the complete chloroplast genome of Arabis stellari and comparisons with related species

Arabis stellari var. japonica is an ornamental plant of the Brassicaceae family, and is widely distributed in South Korea. However, no information is available about its molecular biology and no genomic study has been performed on A. stellari. In this paper, the authors report the complete chloroplast genome sequence of A. stellari. The plastome of A. stellari was 153,683 bp in length with 36.4% GC and included a pair of inverted repeats (IRs) of 26,423 bp that separated a large single-copy (LSC) region of 82,807 bp and a small single-copy (SSC) region of 18,030 bp. It was also found to contain 113 unique genes, of which 79 were protein-coding genes, 30 were transfer RNAs, and four were ribosomal RNAs. The gene content and organization of the A. stellari chloroplast genome were similar to those of other Brassicaceae genomes except for the absence of the rps16 protein-coding gene. A total of 991 SSRs were identified in the genome. The chloroplast genome of A. stellari was compared with closely related species of the Brassicaceae family. Comparative analysis showed a minor divergence occurred in the protein-coding matK, ycf1, ccsA, accD and rpl22 genes and that the K_A/K_S nucleotide substitution ratio of the ndhA genes of A. stellari and A. hirsuta was 1.35135. The genes infA and rps16 were absent in the Arabis genus and phylogenetic evolutionary studies revealed that these genes evolved independently. However, phylogenetic analysis showed that the positions of Brassicaceae species are highly conserved. The present study provides A. stellari genomic information that may be found useful in conservation and molecular phylogenetic studies on Brassicaceae.

Novel DNA Binding and Regulatory Activities for σ54 (RpoN) in Salmonella enterica Serovar Typhimurium 14028s

The variable sigma (σ) subunit of the bacterial RNA polymerase (RNAP) holoenzyme, which is responsible for promoter specificity and open complex formation, plays a strategic role in the response to environmental changes. Salmonella enterica serovar Typhimurium utilizes the housekeeping σ⁷⁰ and five alternative sigma factors, including σ⁵⁴ The σ⁵⁴-RNAP differs from other σ-RNAP holoenzymes in that it forms a stable closed complex with the promoter and requires ATP hydrolysis by an activated cognate bacterial enhancer binding protein (bEBP) to transition to an open complex and initiate transcription. In S. Typhimurium, σ⁵⁴-dependent promoters normally respond to one of 13 different bEBPs, each of which is activated under a specific growth condition. Here, we utilized a constitutively active, promiscuous bEBP to perform a genome-wide identification of σ⁵⁴-RNAP DNA binding sites and the transcriptome of the σ⁵⁴ regulon of S. Typhimurium. The position and context of many of the identified σ⁵⁴ RNAP DNA binding sites suggest regulatory roles for σ⁵⁴-RNAP that connect the σ⁵⁴ regulon to regulons of other σ factors to provide a dynamic response to rapidly changing environmental conditions.IMPORTANCE The alternative sigma factor σ⁵⁴ (RpoN) is required for expression of genes involved in processes with significance in agriculture, bioenergy production, bioremediation, and host-microbe interactions. The characterization of the σ⁵⁴ regulon of the versatile pathogen S. Typhimurium has expanded our understanding of the scope of the σ⁵⁴ regulon and how it links to other σ regulons within the complex regulatory network for gene expression in bacteria.

Repetitive Elements in Mycoplasma hyopneumoniae Transcriptional Regulation

Transcriptional regulation, a multiple-step process, is still poorly understood in the important pig pathogen Mycoplasma hyopneumoniae. Basic motifs like promoters and terminators have already been described, but no other cis-regulatory elements have been found. DNA repeat sequences have been shown to be an interesting potential source of cis-regulatory elements. In this work, a genome-wide search for tandem and palindromic repetitive elements was performed in the intergenic regions of all coding sequences from M. hyopneumoniae strain 7448. Computational analysis demonstrated the presence of 144 tandem repeats and 1,171 palindromic elements. The DNA repeat sequences were distributed within the 5' upstream regions of 86% of transcriptional units of M. hyopneumoniae strain 7448. Comparative analysis between distinct repetitive sequences found in related mycoplasma genomes demonstrated different percentages of conservation among pathogenic and nonpathogenic strains. qPCR assays revealed differential expression among genes showing variable numbers of repetitive elements. In addition, repeats found in 206 genes already described to be differentially regulated under different culture conditions of M. hyopneumoniae strain 232 showed almost 80% conservation in relation to M. hyopneumoniae strain 7448 repeats. Altogether, these findings suggest a potential regulatory role of tandem and palindromic DNA repeats in the M. hyopneumoniae transcriptional profile.

Inter- and intra-strain variability of tandem repeats in Mycoplasma pneumoniae based on next-generation sequencing data

AIM

To characterize inter- and intra-strain variability of variable-number tandem repeats (VNTRs) in Mycoplasma pneumoniae to determine the optimal multilocus VNTR analysis scheme for improved strain typing.

METHODS

Whole genome assemblies and next-generation sequencing data from diverse M. pneumoniae isolates were used to characterize VNTRs and their variability, and to compare the strain discriminability of new VNTR and existing markers.

RESULTS

We identified 13 VNTRs including five reported previously. These VNTRs displayed different levels of inter- and intra-strain copy number variations. All new markers showed similar or higher discriminability compared with existing VNTR markers and the P1 typing system.

CONCLUSION

Our study provides novel insights into VNTR variations and potential new multilocus VNTR analysis schemes for improved genotyping of M. pneumoniae.

Defying Muller's Ratchet: Ancient Heritable Endobacteria Escape Extinction through Retention of Recombination and Genome Plasticity

Heritable endobacteria, which are transmitted from one host generation to the next, are subjected to evolutionary forces that are different from those experienced by free-living bacteria. In particular, they suffer consequences of Muller’s ratchet, a mechanism that leads to extinction of small asexual populations due to fixation of slightly deleterious mutations combined with the random loss of the most-fit genotypes, which cannot be recreated without recombination. Mycoplasma-related endobacteria (MRE) are heritable symbionts of fungi from two ancient lineages, Glomeromycota (arbuscular mycorrhizal fungi) and Mucoromycotina. Previous studies revealed that MRE maintain unusually diverse populations inside their hosts and may have been associated with fungi already in the early Paleozoic. Here we show that MRE are vulnerable to genomic degeneration and propose that they defy Muller’s ratchet thanks to retention of recombination and genome plasticity. We suggest that other endobacteria may be capable of raising similar defenses against Muller’s ratchet.

Survey and analysis of simple sequence repeats (SSRs) in three genomes of Candida species

Simple sequence repeats (SSRs) or microsatellites, which composed of tandem repeated short units of 1-6 bp, have been paying attention continuously. Here, the distribution, composition and polymorphism of microsatellites and compound microsatellites were analyzed in three available genomes of Candida species (Candida dubliniensis, Candida glabrata and Candida orthopsilosis). The results show that there were 118,047, 66,259 and 61,119 microsatellites in genomes of C. dubliniensis, C. glabrata and C. orthopsilosis, respectively. The SSRs covered more than 1/3 length of genomes in the three species. The microsatellites, which just consist of bases A and (or) T, such as (A)n, (T)n, (AT)n, (TA)n, (AAT)n, (TAA)n, (TTA)n, (ATA)n, (ATT)n and (TAT)n, were predominant in the three genomes. The length of microsatellites was focused on 6 bp and 9 bp either in the three genomes or in its coding sequences. What's more, the relative abundance (19.89/kbp) and relative density (167.87 bp/kbp) of SSRs in sequence of mitochondrion of C. glabrata were significantly great than that in any one of genomes or chromosomes of the three species. In addition, the distance between any two adjacent microsatellites was an important factor to influence the formation of compound microsatellites. The analysis may be helpful for further studying the roles of microsatellites in genomes' origination, organization and evolution of Candida species.

Investigation of a Quadruplex-Forming Repeat Sequence Highly Enriched in Xanthomonas and Nostoc sp

In prokaryotes simple sequence repeats (SSRs) with unit sizes of 1–5 nucleotides (nt) are causative for phase and antigenic variation. Although an increased abundance of heptameric repeats was noticed in bacteria, reports about SSRs of 6–9 nt are rare. In particular G-rich repeat sequences with the propensity to fold into G-quadruplex (G4) structures have received little attention. In silico analysis of prokaryotic genomes show putative G4 forming sequences to be abundant. This report focuses on a surprisingly enriched G-rich repeat of the type GGGNATC in Xanthomonas and cyanobacteria such as Nostoc. We studied in detail the genomes of Xanthomonas campestris pv. campestris ATCC 33913 (Xcc), Xanthomonas axonopodis pv. citri str. 306 (Xac), and Nostoc sp. strain PCC7120 (Ana). In all three organisms repeats are spread all over the genome with an over-representation in non-coding regions. Extensive variation of the number of repetitive units was observed with repeat numbers ranging from two up to 26 units. However a clear preference for four units was detected. The strong bias for four units coincides with the requirement of four consecutive G-tracts for G4 formation. Evidence for G4 formation of the consensus repeat sequences was found in biophysical studies utilizing CD spectroscopy. The G-rich repeats are preferably located between aligned open reading frames (ORFs) and are under-represented in coding regions or between divergent ORFs. The G-rich repeats are preferentially located within a distance of 50 bp upstream of an ORF on the anti-sense strand or within 50 bp from the stop codon on the sense strand. Analysis of whole transcriptome sequence data showed that the majority of repeat sequences are transcribed. The genetic loci in the vicinity of repeat regions show increased genomic stability. In conclusion, we introduce and characterize a special class of highly abundant and wide-spread quadruplex-forming repeat sequences in bacteria.

Intrastrand triplex DNA repeats in bacteria: a source of genomic instability

Repetitive nucleic acid sequences are often prone to form secondary structures distinct from B-DNA. Prominent examples of such structures are DNA triplexes. We observed that certain intrastrand triplex motifs are highly conserved and abundant in prokaryotic genomes. A systematic search of 5246 different prokaryotic plasmids and genomes for intrastrand triplex motifs was conducted and the results summarized in the ITxF database available online at http://bioinformatics.uni-konstanz.de/utils/ITxF/. Next we investigated biophysical and biochemical properties of a particular G/C-rich triplex motif (TM) that occurs in many copies in more than 260 bacterial genomes by CD and nuclear magnetic resonance spectroscopy as well as in vivo footprinting techniques. A characterization of putative properties and functions of these unusually frequent nucleic acid motifs demonstrated that the occurrence of the TM is associated with a high degree of genomic instability. TM-containing genomic loci are significantly more rearranged among closely related Escherichia coli strains compared to control sites. In addition, we found very high frequencies of TM motifs in certain Enterobacteria and Cyanobacteria that were previously described as genetically highly diverse. In conclusion we link intrastrand triplex motifs with the induction of genomic instability. We speculate that the observed instability might be an adaptive feature of these genomes that creates variation for natural selection to act upon.

Occurrence, distribution and possible functional roles of simple sequence repeats in phytoplasma genomes

Phytoplasmas are unculturable, cell-wall-less bacteria that parasitize plants and insects. This transkingdom life cycle requires rapid responses to vastly different environments, including transitions from plant phloem sieve elements to various insect tissues and alternations among diverse plant hosts. Features that enable such flexibility in other microbes include simple sequence repeats (SSRs) — mutation-prone, phase-variable short DNA tracts that function as ‘evolutionary rheostats’ and enhance rapid adaptations. To gain insights into the occurrence, distribution and potentially functional roles of SSRs in phytoplasmas, we performed computational analysis on the genomes of five completely sequenced phytoplasma strains, ‘Candidatus Phytoplasma asteris’-related strains OYM and AYWB, ‘Candidatus Phytoplasma australiense’-related strains CBWB and SLY and ‘Candidatus Phytoplasma mali’-related strain AP-AT. The overall density of SSRs in phytoplasma genomes was higher than in representative strains of other prokaryotes. While mono- and trinucleotide SSRs were significantly overrepresented in the phytoplasma genomes, dinucleotide SSRs and other higher-order SSRs were underrepresented. The occurrence and distribution of long SSRs in the prophage islands and phytoplasma-unique genetic loci indicated that SSRs played a role in compounding the complexity of sequence mosaics in individual genomes and in increasing allelic diversity among genomes. Findings from computational analyses were further complemented by an examination of SSRs in varied additional phytoplasma strains, with a focus on potential contingency genes. Some SSRs were located in regions that could profoundly alter the regulation of transcription and translation of affected genes and/or the composition of protein products.

The Cipher Code of Simple Sequence Repeats in "Vampire Pathogens"

Blood inside mammals is a forbidden area for the majority of prokaryotic microbes; however, red blood cells tropism microbes, like “vampire pathogens” (VP), succeed in matching scarce nutrients and surviving strong immunity reactions. Here, we found VP of Mycoplasma, Rhizobiales, and Rickettsiales showed significantly higher counts of (AG)_n dimeric simple sequence repeats (Di-SSRs) in the genomes, coding and non-coding regions than non Vampire Pathogens (N_VP). Regression analysis indicated a significant correlation between GC content and the span of (AG)_n-Di-SSR variation. Gene Ontology (GO) terms with abundance of (AG)₃-Di-SSRs shared by the VP strains were associated with purine nucleotide metabolism (FDR < 0.01), indicating an adaptation to the limited availability of purine and nucleotide precursors in blood. Di-amino acids coded by (AG)_n-Di-SSRs included all three six-fold code amino acids (Arg, Leu and Ser) and significantly higher counts of Di-amino acids coded by (AG)₃, (GA)₃, and (TC)₃ in VP than N_VP. Furthermore, significant differences (P < 0.001) on the numbers of triplexes formed from (AG)_n-Di-SSRs between VP and N_VP in Mycoplasma suggested the potential role of (AG)_n-Di-SSRs in gene regulation.

Specificity and Strain-Typing Capabilities of Nanorod Array-Surface Enhanced Raman Spectroscopy for Mycoplasma pneumoniae Detection

Mycoplasma pneumoniae is a cell wall-less bacterial pathogen of the human respiratory tract that accounts for > 20% of all community-acquired pneumonia (CAP). At present the most effective means for detection and strain-typing is quantitative polymerase chain reaction (qPCR), which can exhibit excellent sensitivity and specificity but requires separate tests for detection and genotyping, lacks standardization between available tests and between labs, and has limited practicality for widespread, point-of-care use. We have developed and previously described a silver nanorod array-surface enhanced Raman Spectroscopy (NA-SERS) biosensing platform capable of detecting M. pneumoniae with statistically significant specificity and sensitivity in simulated and true clinical throat swab samples, and the ability to distinguish between reference strains of the two main genotypes of M. pneumoniae. Furthermore, we have established a qualitative lower endpoint of detection for NA-SERS of < 1 genome equivalent (cell/μl) and a quantitative multivariate detection limit of 5.3 ± 1 cells/μl. Here we demonstrate using partial least squares- discriminatory analysis (PLS-DA) of sample spectra that NA-SERS correctly identified M. pneumoniae clinical isolates from globally diverse origins and distinguished these from a panel of 12 other human commensal and pathogenic mycoplasma species with 100% cross-validated statistical accuracy. Furthermore, PLS-DA correctly classified by strain type all 30 clinical isolates with 96% cross-validated accuracy for type 1 strains, 98% cross-validated accuracy for type 2 strains, and 90% cross-validated accuracy for type 2V strains.

Comparative analysis of microsatellites in chloroplast genomes of lower and higher plants

Microsatellites, or simple sequence repeats (SSRs), contain repetitive DNA sequence where tandem repeats of one to six base pairs are present number of times. Chloroplast genome sequences have been  shown to possess extensive variations in the length, number and distribution of SSRs. However, a comparative analysis of chloroplast microsatellites is not available. Considering their potential importance in generating genomic diversity, we have systematically analysed the abundance and distribution of simple and compound microsatellites in 164 sequenced chloroplast genomes from wide range of plants. The key findings of these studies are (1) a large number of mononucleotide repeats as compared to SSR(2-6)(di-, tri-, tetra-, penta-, hexanucleotide repeats) are present in all chloroplast genomes investigated, (2) lower plants such as algae show wide variation in relative abundance, density and distribution of microsatellite repeats as compared to flowering plants, (3) longer SSRs are excluded from coding regions of most chloroplast genomes, (4) GC content has a weak influence on number, relative abundance and relative density of mononucleotide as well as SSR(2-6). However, GC content strongly showed negative correlation with relative density (R (2) = 0.5, P < 0.05) and relative abundance (R (2) = 0.6, P < 0.05) of cSSRs. In summary, our comparative studies of chloroplast genomes illustrate the variable distribution of microsatellites and revealed that chloroplast genome of smaller plants possesses relatively more genomic diversity compared to higher plants.

Minimal genomes of mycoplasma-related endobacteria are plastic and contain host-derived genes for sustained life within Glomeromycota

Arbuscular mycorrhizal fungi (AMF, Glomeromycota) colonize roots of the majority of terrestrial plants. They provide essential minerals to their plant hosts and receive photosynthates in return. All major lineages of AMF harbor endobacteria classified as Mollicutes, and known as mycoplasma-related endobacteria (MRE). Except for their substantial intrahost genetic diversity and ability to transmit vertically, virtually nothing is known about the life history of these endobacteria. To understand MRE biology, we sequenced metagenomes of three MRE populations, each associated with divergent AMF hosts. We found that each AMF species harbored a genetically distinct group of MRE. Despite vertical transmission, all MRE populations showed extensive chromosomal rearrangements, which we attributed to genetic recombination, activity of mobile elements, and a history of plectroviral invasion. The MRE genomes are characterized by a highly reduced gene content, indicating metabolic dependence on the fungal host, with the mechanism of energy production remaining unclear. Several MRE genes encode proteins with domains involved in protein-protein interactions with eukaryotic hosts. In addition, the MRE genomes harbor genes horizontally acquired from AMF. Some of these genes encode small ubiquitin-like modifier (SUMO) proteases specific to the SUMOylation systems of eukaryotes, which MRE likely use to manipulate their fungal host. The extent of MRE genome plasticity and reduction, along with the large number of horizontally acquired host genes, suggests a high degree of adaptation to the fungal host. These features, together with the ubiquity of the MRE-Glomeromycota associations, emphasize the significance of MRE in the biology of Glomeromycota.

Genome wide survey and analysis of small repetitive sequences in caulimoviruses

Microsatellites are known to exhibit ubiquitous presence across all kingdoms of life including viruses. Members of the Caulimoviridae family severely affect growth of vegetable and fruit plants and reduce economic yield in diverse cropping systems worldwide. Here, we analyzed the nature and distribution of both simple and complex microsatellites present in complete genome of 44 species of Caulimoviridae. Our results showed, in all analyzed genomes, genome size and GC content had a weak influence on number, relative abundance and relative density of microsatellites, respectively. For each genome, mono- and dinucleotide repeats were found to be highly predominant and are overrepresented in genome of majority of caulimoviruses. AT/TA and GAA/AAG/AGA was the most abundant di- and trinucleotide repeat motif, respectively. Repeats larger than trinucleotide were rarely found in these genomes. Comparative study of occurrence, abundance and density of microsatellite among available RNA and DNA viral genomes indicated that simple repeats were least abundant in genomes of caulimoviruses. Polymorphic repeats even though rare were observed in the large intergenic region of the genome, indicating strand slippage and/or unequal recombination processes do occur in caulimoviruses. To our knowledge, this is the first analysis of microsatellites occurring in any dsDNA viral genome. Characterization of such variations in repeat sequences would be important in deciphering the origin, mutational processes, and role of repeat sequences in viral genomes.

Assessing diversity of DNA structure-related sequence features in prokaryotic genomes

Prokaryotic genomes are diverse in terms of their nucleotide and oligonucleotide composition as well as presence of various sequence features that can affect physical properties of the DNA molecule. We present a survey of local sequence patterns which have a potential to promote non-canonical DNA conformations (i.e. different from standard B-DNA double helix) and interpret the results in terms of relationships with organisms' habitats, phylogenetic classifications, and other characteristics. Our present work differs from earlier similar surveys not only by investigating a wider range of sequence patterns in a large number of genomes but also by using a more realistic null model to assess significant deviations. Our results show that simple sequence repeats and Z-DNA-promoting patterns are generally suppressed in prokaryotic genomes, whereas palindromes and inverted repeats are over-represented. Representation of patterns that promote Z-DNA and intrinsic DNA curvature increases with increasing optimal growth temperature (OGT), and decreases with increasing oxygen requirement. Additionally, representations of close direct repeats, palindromes and inverted repeats exhibit clear negative trends with increasing OGT. The observed relationships with environmental characteristics, particularly OGT, suggest possible evolutionary scenarios of structural adaptation of DNA to particular environmental niches.

Simple sequence repeats together with mismatch repair deficiency can bias mutagenic pathways in Pseudomonas aeruginosa during chronic lung infection

Pseudomonas aeruginosa is an opportunistic pathogen that chronically infects the airways of cystic fibrosis (CF) patients and undergoes a process of genetic adaptation based on mutagenesis. We evaluated the role of mononucleotide G:C and A:T simple sequence repeats (SSRs) in this adaptive process. An in silico survey of the genome sequences of 7 P. aeruginosa strains showed that mononucleotide G:C SSRs but not A:T SSRs were greatly under-represented in coding regions, suggesting a strong counterselection process for G:C SSRs with lengths >5 bp but not for A:T SSRs. A meta-analysis of published whole genome sequence data for a P. aeruginosa strain from a CF patient with chronic airway infection showed that G:C SSRs but not A:T SSRs were frequently mutated during the infection process through the insertion or deletion of one or more SSR subunits. The mutation tendency of G:C SSRs was length-dependent and increased exponentially as a function of SSR length. When this strain naturally became a stable Mismatch Repair System (MRS)-deficient mutator, the degree of increase of G:C SSRs mutations (5-fold) was much higher than that of other types of mutation (2.2-fold or less). Sequence analysis of several mutated genes reported for two different collections, both containing mutator and non-mutator strains of P. aeruginosa from CF chronic infections, showed that the proportion of G:C SSR mutations was significantly higher in mutators than in non-mutators, whereas no such difference was observed for A:T SSR mutations. Our findings, taken together, provide genome-scale evidences that under a MRS-deficient background, long G:C SSRs are able to stochastically bias mutagenic pathways by making the genes in which they are harbored more prone to mutation. The combination of MRS deficiency and virulence-related genes that contain long G:C SSRs is therefore a matter of concern in P. aeruginosa CF chronic infection.

A genome-wide analysis of simple sequence repeats in maize and the development of polymorphism markers from next-generation sequence data

Background

Maize (Zea mays ssp. mays L.), as the most important plant for staple food of several million people, animal feed and bioenergy productions, is widely cultivated around the world. Simple sequence repeats (SSRs) are widely used as molecular markers in maize genetics and breeding, but only two thousands pairs of SSRs have been published currently, which hardly satisfies for the increasing needs of geneticists and breeders. Furthermore, the increasing studies have revealed that SSRs also play a vital role in functional regulation and evolution. It is fortunate that the development of sequencing technology and bio-software provides the basis for characterization and development of SSRs in maize.

Results

In this study, MISA was applied to identify overall 179,681 SSRs in maize reference genome B73, with an average distance of 11.46 Kbp. Their distributions within the genome in different regions were non-random, and the density followed in a descending order of UTR, promotor, intron, intergenic and CDS. Meanwhile, 82,694 (46.02%) SSRs with unique flanking sequences were selected, and then applied to analyze the polymorphism of next-generation sequencing data from 345 maize inbred lines and data from maize reference genome B73. There were 58,946 SSRs with length information results in ten or more than ten genomes, accounting for 71.28% of SSRs with unique flanking sequences, while 55,621 SSRs had polymorphism, with an average PIC value of 0.498. 250 pairs of SSR primers in different genomic regions covering all maize chromosomes were randomly chosen for the experimental validation, with an average PIC value of 0.63 in 11 elite maize inbred lines.

Conclusions

Our work provided insight into the non-random distribution spatterns and compositions of SSRs in different regions of maize genome, and also developed more polymorphic SSR markers using next-generation sequencing reads. The genome-wide SSRs polymorphism markers could be useful for genetic analysis and marker-assisted selection in breeding practice, and it was also proved to be high efficient for molecular marker development via next-generation sequencing reads.

The role of variable DNA tandem repeats in bacterial adaptation

DNA tandem repeats (TRs), also designated as satellite DNA, are inter- or intragenic nucleotide sequences that are repeated two or more times in a head-to-tail manner. Because TR tracts are prone to strand-slippage replication and recombination events that cause the TR copy number to increase or decrease, loci containing TRs are hypermutable. An increasing number of examples illustrate that bacteria can exploit this instability of TRs to reversibly shut down or modulate the function of specific genes, allowing them to adapt to changing environments on short evolutionary time scales without an increased overall mutation rate. In this review, we discuss the prevalence and distribution of inter- and intragenic TRs in bacteria and the mechanisms of their instability. In addition, we review evidence demonstrating a role of TR variations in bacterial adaptation strategies, ranging from immune evasion and tissue tropism to the modulation of environmental stress tolerance. Nevertheless, while bioinformatic analysis reveals that most bacterial genomes contain a few up to several dozens of intra- and intergenic TRs, only a small fraction of these have been functionally studied to date.

A Type III restriction-modification system in Mycoplasma mycoides subsp. capri

The sequenced genome of Mycoplasma mycoides subsp. capri revealed the presence of a Type III restriction–modification system (MmyCI). The methyltransferase (modification) subunit of MmyCI (M.MmyCI) was shown to recognize the sequence 5′-TGAG-3′ and methylate the adenine. The coding region of the methyltransferase gene contains 12 consecutive AG dinucleotide repeats that result in a translational termination at a TAA codon immediately beyond the repeat region. This strain does not have MmyCI activity. A clone was found with 10 AG repeats such that the gene is in frame, and this strain has MmyCI activity, suggesting that the expression of the MmyCI methyltransferase may be phase variable.

Coevolution between simple sequence repeats (SSRs) and virus genome size

Background

Relationship between the level of repetitiveness in genomic sequence and genome size has been investigated by making use of complete prokaryotic and eukaryotic genomes, but relevant studies have been rarely made in virus genomes.

Results

In this study, a total of 257 viruses were examined, which cover 90% of genera. The results showed that simple sequence repeats (SSRs) is strongly, positively and significantly correlated with genome size. Certain repeat class is distributed in a certain range of genome sequence length. Mono-, di- and tri- repeats are widely distributed in all virus genomes, tetra- SSRs as a common component consist in genomes which more than 100 kb in size; in the range of genome < 100 kb, genomes containing penta- and hexa- SSRs are not more than 50%. Principal components analysis (PCA) indicated that dinucleotide repeat affects the differences of SSRs most strongly among virus genomes. Results showed that SSRs tend to accumulate in larger virus genomes; and the longer genome sequence, the longer repeat units.

Conclusions

We conducted this research standing on the height of the whole virus. We concluded that genome size is an important factor in affecting the occurrence of SSRs; hosts are also responsible for the variances of SSRs content to a certain degree.

Differences in DNA curvature-related sequence periodicity between prokaryotic chromosomes and phages, and relationship to chromosomal prophage content

Background

Periodic spacing of A-tracts (short runs of A or T) with the DNA helical period of ~10–11 bp is characteristic of intrinsically bent DNA. In eukaryotes, the DNA bending is related to chromatin structure and nucleosome positioning. However, the physiological role of strong sequence periodicity detected in many prokaryotic genomes is not clear.

Results

We developed measures of intensity and persistency of DNA curvature-related sequence periodicity and applied them to prokaryotic chromosomes and phages. The results indicate that strong periodic signals present in chromosomes are generally absent in phage genomes. Moreover, chromosomes containing prophages are less likely to possess a persistent periodic signal than chromosomes with no prophages.

Conclusions

Absence of DNA curvature-related sequence periodicity in phages could arise from constraints associated with DNA packaging in the viral capsid. Lack of prophages in chromosomes with persistent periodic signal suggests that the sequence periodicity and concomitant DNA curvature could play a role in protecting the chromosomes from integration of phage DNA.

Identification and characterization of a lipopolysaccharide α,2,3-sialyltransferase from the human pathogen Helicobacter bizzozeronii

Terminal sialic acid in the lipopolysaccharides (LPSs) of mucosal pathogens is an important virulence factor. Here we report the characterization of a Helicobacter sialyltransferase involved in the biosynthesis of sialylated LPS in Helicobacter bizzozeronii, the only non-pylori gastric Helicobacter species isolated from humans thus far. Starting from the genome sequences of canine and human strains, we identified potential sialyltransferases downstream of three genes involved in the biosynthesis of N-acetylneuraminic acid. One of these candidates showed monofunctional α,2,3-sialyltransferase activity with a preference for N-acetyllactosamine as a substrate. The LPSs from different strains were shown by SDS-PAGE and high-performance anion-exchange chromatography with pulsed amperometric detection (HPAEC-PAD) to contain sialic acid after neuraminidase treatment. The expression of this sialyltransferase and sialyl-LPS appeared to be a phase-variable characteristic common to both human and canine H. bizzozeronii strains. The sialylation site of the LPSs of two H. bizzozeronii strains was determined to be NeuAc-Hex-HexNAc, suggesting terminal 3'-sialyl-LacNAc. Moreover, serological typing revealed the possible presence of sialyl-Lewis X in two additional strains, indicating that H. bizzozeronii could also mimic the surface glycans of mammalian cells. The expression of sialyl-glycans may influence the adaptation process of H. bizzozeronii during the host jump from dogs to humans.

Comparative genomics of Helicobacter pylori and the human-derived Helicobacter bizzozeronii CIII-1 strain reveal the molecular basis of the zoonotic nature of non-pylori gastric Helicobacter infections in humans

BACKGROUND

The canine Gram-negative Helicobacter bizzozeronii is one of seven species in Helicobacter heilmannii sensu lato that are detected in 0.17-2.3% of the gastric biopsies of human patients with gastric symptoms. At the present, H. bizzozeronii is the only non-pylori gastric Helicobacter sp. cultivated from human patients and is therefore a good alternative model of human gastric Helicobacter disease. We recently sequenced the genome of the H. bizzozeronii human strain CIII-1, isolated in 2008 from a 47-year old Finnish woman suffering from severe dyspeptic symptoms. In this study, we performed a detailed comparative genome analysis with H. pylori, providing new insights into non-pylori Helicobacter infections and the mechanisms of transmission between the primary animal host and humans.

RESULTS

H. bizzozeronii possesses all the genes necessary for its specialised life in the stomach. However, H. bizzozeronii differs from H. pylori by having a wider metabolic flexibility in terms of its energy sources and electron transport chain. Moreover, H. bizzozeronii harbours a higher number of methyl-accepting chemotaxis proteins, allowing it to respond to a wider spectrum of environmental signals. In this study, H. bizzozeronii has been shown to have high level of genome plasticity. We were able to identify a total of 43 contingency genes, 5 insertion sequences (ISs), 22 mini-IS elements, 1 genomic island and a putative prophage. Although H. bizzozeronii lacks homologues of some of the major H. pylori virulence genes, other candidate virulence factors are present. In particular, we identified a polysaccharide lyase (HBZC1_15820) as a potential new virulence factor of H. bizzozeronii.

CONCLUSIONS

The comparative genome analysis performed in this study increased the knowledge of the biology of gastric Helicobacter species. In particular, we propose the hypothesis that the high metabolic versatility and the ability to react to a range of environmental signals, factors which differentiate H. bizzozeronii as well as H. felis and H. suis from H. pylori, are the molecular basis of the of the zoonotic nature of H. heilmannii sensu lato infection in humans.

Sequence analysis of the p1 adhesin gene of Mycoplasma pneumoniae in clinical isolates collected in Beijing in 2008 to 2009

The p1 genes of 60 Mycoplasma pneumoniae clinical isolates were sequenced and compared to previously reported p1 gene sequences. An AGT trinucleotide variable-number tandem repeat was identified that ranged in copy number from 5 to 14 among the isolates. In addition, a novel p1 gene variant named 2c was identified in 6 of the isolates.

In silico comparative analysis of SSR markers in plants

BACKGROUND

The adverse environmental conditions impose extreme limitation to growth and plant development, restricting the genetic potential and reflecting on plant yield losses. The progress obtained by classic plant breeding methods aiming at increasing abiotic stress tolerances have not been enough to cope with increasing food demands. New target genes need to be identified to reach this goal, which requires extensive studies of the related biological mechanisms. Comparative analyses in ancestral plant groups can help to elucidate yet unclear biological processes.

RESULTS

In this study, we surveyed the occurrence patterns of expressed sequence tag-derived microsatellite markers for model plants. A total of 13,133 SSR markers were discovered using the SSRLocator software in non-redundant EST databases made for all eleven species chosen for this study. The dimer motifs are more frequent in lower plant species, such as green algae and mosses, and the trimer motifs are more frequent for the majority of higher plant groups, such as monocots and dicots. With this in silico study we confirm several microsatellite plant survey results made with available bioinformatics tools.

CONCLUSIONS

The comparative studies of EST-SSR markers among all plant lineages is well suited for plant evolution studies as well as for future studies of transferability of molecular markers.

Genome nucleotide composition shapes variation in simple sequence repeats

Simple sequence repeats (SSRs) or microsatellites are a common component of genomes but vary greatly across species in their abundance. We tested the hypothesis that this variation is due in part to AT/GC content of genomes, with genomes biased toward either high AT or high CG generating more short random repeats that are long enough to enhance expansion through slippage during replication. To test this hypothesis, we identified repeats with perfect tandem iterations of 1-6 bp from 25 protists with complete or near-complete genome sequences. As expected, the density and the frequency are highly related to genome AT content, with excellent fits to quadratic regressions with minima near a 50% AT content and rising toward both extremes. Within species, the same trends hold, except the limited variation in AT content within each species places each mainly on the descending (GC rich), middle, or ascending (AT rich) part of the curve. The base usages of repeat motifs are also significantly correlated with genome nucleotide compositions: Percentages of AT-rich motifs rise with the increase of genome AT content but vice versa for GC-rich subgroups. Amino acid homopolymer repeats also show the expected quadratic relationship, with higher abundance in species with AT content biased in either direction. Our results show that genome nucleotide composition explains up to half of the variance in the abundance and motif constitution of SSRs.

Comparative analysis of sequence periodicity among prokaryotic genomes points to differences in nucleoid structure and a relationship to gene expression

Regular spacing of short runs of A or T nucleotides in DNA sequences with a period close to the helical period of the DNA double helix has been associated with intrinsic DNA bending and nucleosome positioning in eukaryotes. Analogous periodic signals were also observed in prokaryotic genomes. While the exact role of this periodicity in prokaryotes is not known, it has been proposed to facilitate the DNA packaging in the prokaryotic nucleoid and/or to promote negative or positive supercoiling. We developed a methodology for assessments of intragenomic heterogeneity of these periodic patterns and applied it in analysis of 1,025 prokaryotic chromosomes. This technique allows more detailed analysis of sequence periodicity than previous methods where sequence periodicity was assessed in an integral form across the whole chromosome. We found that most genomes have the periodic signal confined to several chromosomal segments while most of the chromosome lacks a strong sequence periodicity. Moreover, there are significant differences among different prokaryotes in both the intensity and persistency of sequence periodicity related to DNA curvature. We proffer that the prokaryotic nucleoid consists of relatively rigid sections stabilized by short intrinsically bent DNA segments and characterized by locally strong periodic patterns alternating with regions featuring a weak periodic signal, which presumably permits higher structural flexibility. This model applies to most bacteria and archaea. In genomes with an exceptionally persistent periodic signal, highly expressed genes tend to concentrate in aperiodic sections, suggesting that structural heterogeneity of the nucleoid is related to local differences in transcriptional activity.

Finding sequence motifs in prokaryotic genomes--a brief practical guide for a microbiologist

Finding significant nucleotide sequence motifs in prokaryotic genomes can be divided into three types of tasks: (1) supervised motif finding, where a sample of motif sequences is used to find other similar sequences in genomes; (2) unsupervised motif finding, which typically relates to the task of finding regulatory motifs and protein binding sites and (3) exploratory motif finding, which aims to identify potential functionally significant sequence motifs as those that are unusual in some statistical sense. This article provides a conceptual overview for each type of task, a brief description of basic algorithms used in their solution, and a review of selected relevant software available online.

Mhp493 (P216) is a proteolytically processed, cilium and heparin binding protein of Mycoplasma hyopneumoniae

Mycoplasma hyopneumoniae induces respiratory disease in swine by colonizing cilia causing ciliostasis, cilial loss and epithelial cell death. Heparin binds to M. hyopneumoniae cells in a dose-dependent manner and blocks its ability to adhere to porcine cilia. We show here that Mhp493 (P216), a paralogue of the cilium adhesin P97 (Mhp183), is cleaved between amino acids 1040 and 1089 generating surface-accessible, heparin-binding proteins P120 and P85. Antiphosphoserine antibodies recognized P85 in 2-D immunoblotting studies and TiO(2) chromatography of trypsin digests of P85 isolated a single peptide with an m/z of 917.3. A phosphoserine residue in the tryptic peptide (90)VSELpSFR(96) (position 94 in P85) was identified by MALDI-MS/MS. Polyhistidine fusion proteins (F1(P216), F2(P216), F3(P216)) spanning Mhp493 bound heparin with biologically significant Kd values, and heparin, fucoidan and mucin inhibited this interaction. Latex beads coated with F1(P216), F2(P216) and F3(P216) adhered to and entered porcine kidney epithelial-like (PK15) cell monolayers. Microtitre plate-based assays showed that sequences within P120 and P85 bind to porcine cilia and are recognized by serum antibodies elicited during infection by M. hyopneumoniae. Mhp493 contributes significantly to the surface architecture of M. hyopneumoniae and is the first cilium adhesin to be described that lacks an R1 cilium-binding domain.

Long simple sequence repeats in host-adapted pathogens localize near genes encoding antigens, housekeeping genes, and pseudogenes

Simple sequence repeats (SSRs) in DNA sequences are tandem iterations of a single nucleotide or a short oligonucleotide. SSRs are subject to slipped-strand mutations and a common source of phase variation in bacteria and antigenic variation in pathogens. Significantly long SSRs are generally rare in prokaryotic genomes, and long SSRs composed of iterations of mono-, di-, tri-, and tetranucleotides are mostly restricted to host-adapted pathogens. We present new results concerning associations between long SSRs and genes related to different cellular functions in genomes of host-adapted pathogens. We found that in the majority of the analyzed genomes, at least some of the genes associated with SSRs encode potential antigens, which is expected if the primary function of SSRs is their contribution to antigenic variation. However, we also found a number of long SSRs associated with housekeeping genes, including rRNA and tRNA genes, genes encoding ribosomal proteins, amino acyl-tRNA synthetases, chaperones, and important metabolic enzymes. Many of these genes are probably essential and it is unlikely that they are phase-variable. Few statistically significant associations between SSRs and gene functional classifications were detected, suggesting that most long SSRs are not related to a particular cellular function or process. Long SSRs in Mycobacterium leprae are mostly associated with pseudogenes and may be contributing to gene loss following the adaptation to an obligate pathogenic lifestyle. We speculate that LSSRs may have played a similar role in genome reduction of other host-adapted pathogens.