Simple sequence repeats in mycobacterial genomes

Investigating in vivo Mycobacterium avium subsp. paratuberculosis microevolution and mixed strain infections

Mycobacterium avium subsp. paratuberculosis (MAP) causes Johne's Disease (JD) in ruminants, which is responsible for significant economic loss to the global dairy industry. Mixed strain infection (MSI) refers to the concurrent infection of a susceptible host with genetically distinct strains of a pathogen, whereas within-host changes in an infecting strain leading to genetically distinguishable progeny is called microevolution. The two processes can influence host-pathogen dynamics, disease progression and outcomes, but not much is known about their prevalence and impact on JD. Therefore, we obtained up to 10 MAP isolates each from 14 high-shedding animals and subjected them to whole-genome sequencing. Twelve of the 14 animals examined showed evidence for the presence of MSIs and microevolution, while the genotypes of MAP isolates from the remaining two animals could be attributed solely to microevolution. All MAP isolates that were otherwise isogenic had differences in short sequence repeats (SSRs), of which SSR1 and SSR2 were the most diverse and homoplastic. Variations in SSR1 and SSR2, which are located in ORF1 and ORF2, respectively, affect the genetic reading frame, leading to protein products with altered sequences and computed structures. The ORF1 gene product is predicted to be a MAP surface protein with possible roles in host immune modulation, but nothing could be inferred regarding the function of ORF2. Both genes are conserved in Mycobacterium avium complex members, but SSR1-based modulation of ORF1 reading frames seems to only occur in MAP, which could have potential implications on the infectivity of this pathogen. IMPORTANCE Johne's disease (JD) is a major problem in dairy animals, and concerns have been raised regarding the association of Mycobacterium avium subsp. paratuberculosis (MAP) with Crohn's disease in humans. MAP is an extremely slow-growing bacterium with low genome evolutionary rates. Certain short sequence repeats (SSR1 and SSR2) in the MAP chromosome are highly variable and evolve at a faster rate than the rest of the chromosome. In the current study, multiple MAP isolates with genetic variations such as single-nucleotide polymorphisms, and more noticeably, diverse SSRs, could simultaneously infect animals. Variations in SSR1 and SSR2 affect the products of the respective genes containing them. Since multiple MAP isolates can infect the same animal and the possibility that the pathogen undergoes further changes within the host due to unstable SSRs, this could provide a compensative mechanism for an otherwise slow-evolving pathogen to increase phenotypic diversity for overcoming host responses.

Distinctive roles of translesion polymerases DinB1 and DnaE2 in diversification of the mycobacterial genome through substitution and frameshift mutagenesis

Antibiotic resistance of Mycobacterium tuberculosis is exclusively a consequence of chromosomal mutations. Translesion synthesis (TLS) is a widely conserved mechanism of DNA damage tolerance and mutagenesis, executed by translesion polymerases such as DinBs. In mycobacteria, DnaE2 is the only known agent of TLS and the role of DinB polymerases is unknown. Here we demonstrate that, when overexpressed, DinB1 promotes missense mutations conferring resistance to rifampicin, with a mutational signature distinct from that of DnaE2, and abets insertion and deletion frameshift mutagenesis in homo-oligonucleotide runs. DinB1 is the primary mediator of spontaneous −1 frameshift mutations in homo-oligonucleotide runs whereas DnaE2 and DinBs are redundant in DNA damage-induced −1 frameshift mutagenesis. These results highlight DinB1 and DnaE2 as drivers of mycobacterial genome diversification with relevance to antimicrobial resistance and host adaptation.

This manuscript elucidates new mechanisms of mutagenesis in mycobacteria by implicating two translesion DNA polymerases in genome diversification, including creating the mutations that underlie all antibiotic resistance in these global pathogens.

Identification of a Non-Pentapeptide Region Associated with Rapid Mycobacterial Evolution

A large portion of the coding capacity of Mycobacterium tuberculosis is devoted to the production of proteins containing several copies of the pentapeptide-2 repeat, namely the PE/PPE_MPTR proteins. Protein domain repeats have a variety of binding properties and are involved in protein-protein interactions as well as binding to other ligands such as DNA and RNA. They are not as common in prokaryotes, compared to eukaryotes, but the enrichment of pentapeptide-2 repeats in Mycobacteria constitutes an exception to that rule. The genes encoding the PE/PPE_MPTR proteins have undergone many rearrangements and here we have identified the expansion patterns across the Mycobacteria. We have performed a reclassification of the PE/PPE_MPTR proteins using cohesive regions rather than sparse domain architectures. It is clear that these proteins have undergone large insertions of several pentapeptide-2 domains appearing adjacent to one another in a repetitive pattern. Further, we have identified a non-pentapeptide motif associated with rapid mycobacterial evolution. The sequence composition of this region suggests a different structure compared to pentapeptide-2 repeats. By studying the evolution of the PE/PPE_MPTR proteins, we have distinguished features pertaining to tuberculosis-inducing species. Further studies of the non-pentapeptide region associated with repeat expansions promises to shed light on the pathogenicity of Mycobacterium tuberculosis.

Revisiting an important component of plant genomes: microsatellites

Microsatellites are some of the most highly variable repetitive DNA tracts in genomes. Few studies focus on whether the characteristic instability of microsatellites is linked to phenotypic effects in plants. We summarise recent data to investigate how microsatellite variations affect gene expression and hence phenotype. We discuss how the basic characteristics of microsatellites may contribute to phenotypic effects. In summary, microsatellites in plants are universal and highly mutable, they coexist and coevolve with transposable elements, and are under selective pressure. The number of motif nucleotides, the type of motif and transposon activity all contribute to the nonrandom generation and decay of microsatellites, and to conservation and distribution biases. Although microsatellites are generated by accident, they mature through responses to environmental change before final decay. This process is mediated by organism adjustment mechanisms, which maintain a balance between birth versus death and growth versus decay in microsatellites. Close relationships also exist between the physical structure, variation and functionality of microsatellites: in most plant species, sequences containing microsatellites are associated with catalytic activity and binding functions, are expressed in the membrane and organelles, and participate in the developmental and metabolic processes. Microsatellites contribute to genome structure and functional plasticity, and may be considered to promote species evolution in plants in response to environmental changes. In conclusion, the generation, loss, functionality and evolution of microsatellites can be related to plant gene expression and functional alterations. The effect of microsatellites on phenotypic variation may be as significant in plants as it is in animals.

Direct and inverted repeats elicit genetic instability by both exploiting and eluding DNA double-strand break repair systems in mycobacteria

Repetitive DNA sequences with the potential to form alternative DNA conformations, such as slipped structures and cruciforms, can induce genetic instability by promoting replication errors and by serving as a substrate for DNA repair proteins, which may lead to DNA double-strand breaks (DSBs). However, the contribution of each of the DSB repair pathways, homologous recombination (HR), non-homologous end-joining (NHEJ) and single-strand annealing (SSA), to this sort of genetic instability is not fully understood. Herein, we assessed the genome-wide distribution of repetitive DNA sequences in the Mycobacterium smegmatis, Mycobacterium tuberculosis and Escherichia coli genomes, and determined the types and frequencies of genetic instability induced by direct and inverted repeats, both in the presence and in the absence of HR, NHEJ, and SSA. All three genomes are strongly enriched in direct repeats and modestly enriched in inverted repeats. When using chromosomally integrated constructs in M. smegmatis, direct repeats induced the perfect deletion of their intervening sequences ~1,000-fold above background. Absence of HR further enhanced these perfect deletions, whereas absence of NHEJ or SSA had no influence, suggesting compromised replication fidelity. In contrast, inverted repeats induced perfect deletions only in the absence of SSA. Both direct and inverted repeats stimulated excision of the constructs from the attB integration sites independently of HR, NHEJ, or SSA. With episomal constructs, direct and inverted repeats triggered DNA instability by activating nucleolytic activity, and absence of the DSB repair pathways (in the order NHEJ>HR>SSA) exacerbated this instability. Thus, direct and inverted repeats may elicit genetic instability in mycobacteria by 1) directly interfering with replication fidelity, 2) stimulating the three main DSB repair pathways, and 3) enticing L5 site-specific recombination.

Mining of EST-SSR markers of Musa and their transferability studies among the members of order the Zingiberales

Expressed sequence tags (ESTs) databases of 11 Musa complementary DNA libraries were retrieved from National Center of Biotechnology Information and used for mining simple sequence repeats (SSRs). Out of 21,056 unique ESTs, SSR regions were found only in 5,158 ESTs. Among these SSR containing ESTs, the occurrence of trinucleotide repeats are the most abundant followed by mono-, di-, tetra-, hexa-, and pentanucleotides. Moreover, this study showed that the rate of class II SSRs (<20 nucleotides) was higher than the class I SSRs (<20 nucleotides), and proportion of class I and II SSRs as abundant for tri-repeats. As a representative sample, primers were synthesized for 24 ESTs, carrying >12 nucleotides of SSR region, and tested among the various genomic group of Musa accessions. The result showed that 88 % of primers were functional primers, and 43 % are showing polymorphism among the Musa accessions. Transferability studies of Musa EST-SSRs among the genera of the order Zingiberales exhibited 100 and 58 % transferability in Musaceae and Zingiberaceae, respectively. The sequence comparison of SSR regions among the different Musa accessions confirmed that polymorphism is mainly due to the variation in repeat length. High percentage of cross-species, cross-genera, and cross-family transferability also suggested that these Musa EST-SSR markers will be a valuable resource for the comparative mapping by developing COS markers, in evolutionary studies and in improvement of the members of Zingiberaceae and Musaceae.

A study on mutational dynamics of simple sequence repeats in relation to mismatch repair system in prokaryotic genomes

Mutational bias toward expansion or contraction of simple sequence repeats (SSRs) is referred to as directionality of SSR evolution. In this communication, we report the mutational bias exhibited by mononucleotide SSRs occurring in the non-coding regions of several prokaryotic genomes. Our investigations revealed that the strains or species lacking mismatch repair (MMR) system generally show higher number of polymorphic SSRs than those species/strains having MMR system. An exception to this observation was seen in the mycobacterial genomes that are MMR deficient where only a few SSR tracts were seen with mutations. This low incidence of SSR mutations even in the MMR-deficient background could be attributed to the high fidelity of the DNA polymerases as a consequence of high generation time of the mycobacteria. MMR system-deficient species generally did not show any bias toward mononucleotide SSR expansions or contractions indicating a neutral evolution of SSRs in these species. The MMR-proficient species in which the observed mutations correspond to secondary mutations showed bias toward contraction of polymononucleotide tracts, perhaps, indicating low efficiency of MMR system to repair SSR-induced slippage errors on template strands. This bias toward deletion in the mononucleotide SSR tracts might be a probable reason behind scarcity for long poly A|T and G|C tracts in prokaryotic systems which are mostly MMR proficient. In conclusion, our study clearly demonstrates mutational dynamics of SSRs in relation to the presence/absence of MMR system in the prokaryotic system.

Perspective on sequence evolution of microsatellite locus (CCG)n in Rv0050 gene from Mycobacterium tuberculosis

Background

The mycobacterial genome is inclined to polymerase slippage and a high mutation rate in microsatellite regions due to high GC content and absence of a mismatch repair system. However, the exact molecular mechanisms underlying microsatellite variation have not been fully elucidated. Here, we investigated mutation events in the hyper-variable trinucleotide microsatellite locus MML0050 located in the Rv0050 gene of W-Beijing and non-W-Beijing Mycobacterium tuberculosis strains in order to gain insight into the genomic structure and activity of repeated regions.

Results

Size analysis indicated the presence of five alleles that differed in length by three base pairs. Moreover, nucleotide gains occurred more frequently than loses in this trinucleotide microsatellite. Mutation frequency was not completely related with the total length, though the relative frequency in the longest allele was remarkably higher than that in the shortest. Sequence analysis was able to detect seven alleles and revealed that point mutations enhanced the level of locus variation. Introduction of an interruptive motif correlated with the total allele length and genetic lineage, rather than the length of the longest stretch of perfect repeats. Finally, the level of locus variation was drastically different between the two genetic lineages.

Conclusion

The Rv0050 locus encodes the bifunctional penicillin-binding protein ponA1 and is essential to mycobacterial survival. Our investigations of this particularly dynamic genomic region provide insights into the overall mode of microsatellite evolution. Specifically, replication slippage was implicated in the mutational process of this microsatellite and a sequence-based genetic analysis was necessary to determine that point mutation events acted to maintain microsatellite size integrity while providing genomic diversity.

How dormant is Mycobacterium tuberculosis during latency? A study integrating genomics and molecular epidemiology

Mycobacterium tuberculosis may survive for decades in the human body in a state termed latent tuberculosis infection (LTBI). We investigated the occurrence during LTBI of insertion/deletion events in a selected set of mononucleotide simple sequence repeats, DNA sequence changes in four M. tuberculosis genes, and large sequence variations in 4750 M. tuberculosis open reading frames. We studied 13 paired M. tuberculosis clinical isolates, with each pair representing a reactivation of LTBI more than three decades after primary infection. Absence of sequence variations between paired isolates in nearly all investigated loci suggests a low likelihood of bacterial replication during LTBI.

Mining of expressed sequence tag libraries of cacao for microsatellite markers using five computational tools

Expressed sequence tags (ESTs) provide researchers with a quick and inexpensive route for discovering new genes, data on gene expression and regulation, and also provide genic markers that help in constructing genome maps. Cacao is an important perennial crop of humid tropics. Cacao EST sequences, as available in the public domain, were downloaded and made into contigs. Microsatellites were located in these ESTs and contigs using five softwares (MISA, TRA, TROLL, SSRIT and SSR primer). MISA gave maximum coverage of SSRs in cacao ESTs and contigs, although TRA was able to detect higher order (5-mer) repeats. The frequency of SSRs was one per 26.9 kb in the known set of ESTs. One-third of the repeats in EST-contigs were found to be trimeric. A few rare repeats like 21-mer repeat were also located. A/T repeats were most abundant among the mononucleotide repeats and the AG/GA/TC/CT type was the most frequent among dimerics. Flanking primers were designed using Primer3 program and verified experimentally for PCR amplification. The results of the study are made available freely online database (http://riju.byethost31.com/cocoa/). Seven primer pairs amplified genomic DNA isolated from leaves were used to screen a representative set of 12 accessions of cacao.

Stabilization of the genome of the mismatch repair deficient Mycobacterium tuberculosis by context-dependent codon choice

BACKGROUND

The rate at which a stretch of DNA mutates is determined by the cellular systems for DNA replication and repair, and by the nucleotide sequence of the stretch itself. One sequence feature with a particularly strong influence on the mutation rate are nucleotide repeats. Some microbial pathogens use nucleotide repeats in their genome to stochastically vary phenotypic traits and thereby evade host defense. However, such unstable sequences also come at a cost, as mutations are often deleterious. Here, we analyzed how these opposing forces shaped genome stability in the human pathogen Mycobacterium tuberculosis. M. tuberculosis lacks a mismatch repair system, and this renders nucleotide repeats particularly unstable.

RESULTS

We found that proteins of M. tuberculosis are encoded by using codons in a context-dependent manner that prevents the emergence of nucleotide repeats. This context-dependent codon choice leads to a strong decrease in the estimated frame-shift mutation rate and thus to an increase in genome stability.

CONCLUSION

These results indicate that a context-specific codon choice can partially compensate for the lack of a mismatch repair system, and helps to maintain genome integrity in this pathogen.

IMEx: Imperfect Microsatellite Extractor

MOTIVATION

Microsatellites, also known as simple sequence repeats, are the tandem repeats of nucleotide motifs of the size 1-6 bp found in every genome known so far. Their importance in genomes is well known. Microsatellites are associated with various disease genes, have been used as molecular markers in linkage analysis and DNA fingerprinting studies, and also seem to play an important role in the genome evolution. Therefore, it is of importance to study distribution, enrichment and polymorphism of microsatellites in the genomes of interest. For this, the prerequisite is the availability of a computational tool for extraction of microsatellites (perfect as well as imperfect) and their related information from whole genome sequences. Examination of available tools revealed certain lacunae in them and prompted us to develop a new tool.

RESULTS

In order to efficiently screen genome sequences for microsatellites (perfect as well as imperfect), we developed a new tool called IMEx (Imperfect Microsatellite Extractor). IMEx uses simple string-matching algorithm with sliding window approach to screen DNA sequences for microsatellites and reports the motif, copy number, genomic location, nearby genes, mutational events and many other features useful for in-depth studies. IMEx is more sensitive, efficient and useful than the available widely used tools. IMEx is available in the form of a stand-alone program as well as in the form of a web-server.

AVAILABILITY

A World Wide Web server and the stand-alone program are available for free access at http://203.197.254.154/IMEX/ or http://www.cdfd.org.in/imex.