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

Third Generation Genome Sequencing Reveals That Endobacteria in Nematophagous Fungi Esteya vermicola Contain Multiple Genes Encoding for Nematicidal Proteins

Esteya vermicola is the first recorded endoparasitic nematophagous fungus with high infectivity capacity, attacking the pinewood nematode Bursaphelenchus xylophilus which causes pine wilt disease. Endosymbionts are found in the cytoplasm of E. vermicola from various geographical areas. We sequenced the genome of endobacteria residing in E. vermicola to discover possible biological functions of these widespread endobacteria. Multilocus phylogenetic analyses showed that the endobacteria form a previously unidentified lineage sister to Phyllobacterium myrsinacearum species. The number of genes in the endobacterium was 4542, with 87.8% of the proteins having a known function. It contained a high proportion of repetitive sequences, as well as more Acyl-CoA synthetase genes and genes encoding the electron transport chain, compared with compared with plant-associated P. zundukense Tri 48 and P. myrsinacearum DSM 5893. Thus, this symbiotic bacterium is likely to be more efficient in regulating gene expression and energy release. Furthermore, the endobacteria in nematophagous fungi Esteya vermicola contained multiple nematicidal subtilase/subtilisin encoding genes, so it is likely that endobacteria cooperate with the host to kill nematodes.

DNA tandem repeats contribute to the genetic diversity of Brevibacterium aurantiacum phages

This report presents the characterization of the first virulent phages infecting Brevibacterium aurantiacum, a bacterial species used during the manufacture of surface-ripened cheeses. These phages were also responsible for flavour and colour defects in surface-ripened cheeses. Sixteen phages (out of 62 isolates) were selected for genome sequencing and comparative analyses. These cos-type phages with a long non-contractile tail currently belong to the Siphoviridae family (Caudovirales order). Their genome sizes vary from 35,637 to 36,825 bp and, similar to their host, have a high GC content (~61%). Genes encoding for an immunity repressor, an excisionase and a truncated integrase were found, suggesting that these virulent phages may be derived from a prophage. Their genomic organization is highly conserved, with most of the diversity coming from the presence of long (198 bp) DNA tandem repeats (TRs) within an open reading frame coding for a protein of unknown function. We categorized these phages into seven genomic groups according to their number of TR, which ranged from two to eight. Moreover, we showed that TRs are widespread in phage genomes, found in more than 85% of the genomes available in public databases.

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.

A comparative survey of microsatellites among wild and domestic cat provides valuable resources for marker development

Information on the level and distribution of genetic variation is important for conservation plan of captive population of an endangered species such as tiger and cheetah. We assayed the frequency of microsatellites in the genomic and genic sequences of wild cats (Panthera tigris, Acinonyx jubatus) and compared it with the domestic cat (Felis catus). Frequency, relative abundance and density of microsatellites were highest in the domestic cat when compared with wild cats. The frequency of microsatellites was positively correlated with the G+C content of genomic and genic sequences. The maximum frequency of microsatellites among all three sequence sets was of di-nucleotide repeats (genomic-88.1%; genic-70.4%), whereas the hexa-nucleotide repeat represents < 0.5%. Motif conservation study among the genomic and genic sequences revealed conservation of 81.3% and 51.0% motif within the members of family Felidae. A total 40,233 primers from genic sequences were designed in order to enrich the members of family Felidae with genomic resources. The designed primers could be useful in determining the molecular genetics of population structure and individualization of a particular cat.

A Comparison of Microsatellites in Phytopathogenic Aspergillus Species in Order to Develop Markers for the Assessment of Genetic Diversity among Its Isolates

The occurrence of Microsatellites (SSRs) has been witnessed in most of the fungal genomes however its abundance varies across species. In the present study, we analyzed the frequency of SSRs in the whole genome and transcripts of two phyto-pathogenic (Aspergillus niger and Aspergillus terreus) and compared them with two non-pathogenic (Aspergillus nidulans and Aspergillus oryzae) Aspergillus. Higher relative abundance and relative density of SSRs were observed in the whole genome and transcript sequences of the pathogenic Aspergillus when compared to the non-pathogenic. The relative abundance and density of SSRs were positively correlated with the G+C content of transcripts. Among the different classes of SSR, the percentage of tetra-nucleotide SSRs were maximum in A. niger (36.7%) and A. oryzae (35.9%) whereas A. nidulans and A. terreus preferred tri-nucleotide SSRs (38.2 and 42.1%) in whole genome sequences. In transcripts, tri-nucleotide SSRs were the most abundant whereas di-nucleotide SSRs were the least favored. Motif conservation study among the transcripts revealed conservation of only 27% motif within Aspergillus species. Furthermore, a similar relationship among the Ascomycetes was obtained on the basis of motif conservation and conserved genes (rDNA). To analyze the diversity present within the Indian isolates of Aspergillus, primers were successfully designed for 692 motifs in A. niger and A. terreus of which 20 were selected for diversity analysis. Among all the markers amplified, 10 markers (83.3%) were polymorphic, whereas remaining two markers (16.6%) were monomorphic. Ten polymorphic markers acquired in this investigation showed the utility of recently created SSR markers in the assessment of genetic diversity among various isolates of Aspergillus.

Abundant Intergenic TAACTGA Direct Repeats and Putative Alternate RNA Polymerase β' Subunits in Marine Beggiatoaceae Genomes: Possible Regulatory Roles and Origins

The genome sequences of several giant marine sulfur-oxidizing bacteria present evidence of a possible post-transcriptional regulatory network that may have been transmitted to or from two distantly related bacteria lineages. The draft genome of a Cand. “Maribeggiatoa” filament from the Guaymas Basin (Gulf of California, Mexico) seafloor contains 169 sets of TAACTGA direct repeats and one indirect repeat, with two to six copies per set. Related heptamers are rarely or never found as direct repeats. TAACTGA direct repeats are also found in some other Beggiatoaceae, Thiocystis violascens, a range of Cyanobacteria, and five Bacteroidetes. This phylogenetic distribution suggests they may have been transmitted horizontally, but no mechanism is evident. There is no correlation between total TAACTGA occurrences and repeats per genome. In most species the repeat units are relatively short, but longer arrays of up to 43 copies are found in several Bacteroidetes and Cyanobacteria. The majority of TAACTGA repeats in the Cand. “Maribeggiatoa” Orange Guaymas (BOGUAY) genome are within several nucleotides upstream of a putative start codon, suggesting they may be binding sites for a post-transcriptional regulator. Candidates include members of the ribosomal protein S1, Csp (cold shock protein), and Csr (carbon storage regulator) families. No pattern was evident in the predicted functions of the open reading frames (ORFs) downstream of repeats, but some encode presumably essential products such as ribosomal proteins. Among these is an ORF encoding a possible alternate or modified RNA polymerase beta prime subunit, predicted to have the expected subunit interaction domains but lacking most catalytic residues. A similar ORF was found in the Thioploca ingrica draft genome, but in no others. In both species they are immediately upstream of putative sensor kinase genes with nearly identical domain structures. In the marine Beggiatoaceae, a role for the TAACTGA repeats in translational regulation is suggested. More speculatively, the putative alternate RNA polymerase subunit could be a negative transcriptional regulator.

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.

The discovery, function and development of the variable number tandem repeats in different Mycobacterium species

The method of genotyping by variable number tandem repeats (VNTRs) facilitates the epidemiological studies of different Mycobacterium species worldwide. Until now, the VNTR method is not fully understood, for example, its discovery, function and classification. The inconsistent nomenclature and terminology of VNTR is especially confusing. In this review, we first describe in detail the VNTRs in Mycobacterium tuberculosis (M. tuberculosis), as this pathogen resulted in more deaths than any other microbial pathogen as well as for which extensive studies of VNTRs were carried out, and then we outline the recent progress of the VNTR-related epidemiological research in several other Mycobacterium species, such as M. abscessus, M. africanum, M. avium, M. bovis, M. canettii, M. caprae, M. intracellulare, M. leprae, M. marinum, M. microti, M. pinnipedii and M. ulcerans from different countries and regions. This article is aimed mainly at the practical notes of VNTR to help the scientists in better understanding and performing this method.

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.

Spatial genetic structure in natural populations of Phragmites australis in a mosaic of saline habitats in the Yellow River Delta, China

Determination of spatial genetic structure (SGS) in natural populations is important for both theoretical aspects of evolutionary genetics and their application in species conservation and ecological restoration. In this study, we examined genetic diversity within and among the natural populations of a cosmopolitan grass Phragmites australis (common reed) in the Yellow River Delta (YRD), China, where a mosaic of habitat patches varying in soil salinity was detected. We demonstrated that, despite their close geographic proximity, the common reed populations in the YRD significantly diverged at six microsatellite loci, exhibiting a strong association of genetic variation with habitat heterogeneity. Genetic distances among populations were best explained as a function of environmental difference, rather than geographical distance. Although the level of genetic divergence among populations was relatively low (F’_ST = 0.073), weak but significant genetic differentiation, as well as the concordance between ecological and genetic landscapes, suggests spatial structuring of genotypes in relation to patchy habitats. These findings not only provided insights into the population dynamics of common reed in changing environments, but also demonstrated the feasibility of using habitat patches in a mosaic landscape as test systems to identify appropriate genetic sources for ecological restoration.

A longitudinal study of the diversity and dynamics of Mycoplasma hyopneumoniae infections in pig herds

A longitudinal study was carried out to investigate the diversity and persistence of Mycoplasma hyopneumoniae (M. hyopneumoniae) strains in four infected pig herds. In each herd, 20 pigs were randomly selected and blood and/or bronchoalveolar lavage (BAL) fluid was collected at 6, 10, 14 and 26 weeks of age. In the BAL fluid, quantitative PCR and MLVA (multiple-locus variable number of tandem repeats (VNTR) analysis) testing were performed for detection and typing of M. hyopneumoniae strains, respectively. At 26 weeks of age, the prevalence and severity of lung lesions were recorded at slaughter (minimum 50 pigs belonging to the same batch as the investigated pigs). The percentage of pigs testing positive on qPCR increased from 35% at 6 weeks to 96% at 26 weeks of age. With MLVA testing, positive pigs were found from 14 weeks onwards. Within each herd, only one distinct strain was detected, although clonal variants were identified in two herds. In three of the herds, the strain remained present until slaughter age. The percentage of pigs with Mycoplasma-like lesions ranged from 38% to 98%, and the average pneumonia score ranged from 1.7 to 11.9, respectively. The present field study documented that within a herd, mainly one distinct M. hyopneumoniae strain was present that persisted in the same animals for at least 12 weeks. This implies that the immune response of the animals following infection is not able to rapidly clear the infection from the respiratory tract.

The pseudogenes of Mycobacterium leprae reveal the functional relevance of gene order within operons

Almost 50 years following the discovery of the prokaryotic operon, the functional relevance of gene order within operons remains unclear. In this work, we take advantage of the eroded genome of Mycobacterium leprae to add evidence supporting the notion that functionally less important genes have a tendency to be located at the end of its operons. M. leprae’s genome includes 1133 pseudogenes and 1614 protein-coding genes and can be compared with the close genome of M. tuberculosis. Assuming M. leprae’s pseudogenes to represent dispensable genes, we have studied the position of these pseudogenes in the operons of M. leprae and of their orthologs in M. tuberculosis. We observed that both tend to be located in the 3′ (downstream) half of the operon (P-values of 0.03 and 0.18, respectively). Analysis of pseudogenes in all available prokaryotic genomes confirms this trend (P-value of 7.1 × 10⁻⁷). In a complementary analysis, we found a significant tendency for essential genes to be located at the 5′ (upstream) half of the operon (P-value of 0.006). Our work provides an indication that, in prokarya, functionally less important genes have a tendency to be located at the end of operons, while more relevant genes tend to be located toward operon starts.

Comparative genomic analysis of simple sequence repeats in three Plasmodium species

Simple sequence repeats (SSRs) are known to be responsible for genetic complexities and play major roles in gene and genome evolution. To this respect, malaria parasites are known to have rapidly evolving and complex genomes with complicated and differential pathogenic behaviors. Hence, by studying the whole genome comparative SSRs patterns, one can understand genomic complexities and differential evolutionary patterns of these species. We herein utilized the whole genome sequence information of three Plasmodium species, Plasmodium falciparum, Plasmodium vivax, and Plasmodium knowlesi, to comparatively analyze genome-wide distribution of SSRs. The study revealed that despite having the smallest genome size, P. falciparum bears the highest SSR content among the three Plasmodium species. Furthermore, distribution patterns of different SSRs types (e.g., mono, di, tri, tetra, penta, and hexa) in term of relative abundance and relative density provide evidences for greater accumulation of di-repeats and marked decrease of mono-repeats in P. falciparum in comparison to other two species. Overall, the types and distribution of SSRs in P. falciparum genome was found to be different than that of P. vivax and P. knowlesi. The latter two species have quite similar SSR organizations in many aspects of the data. The results were discussed in terms of comparative SSR patterns among the three Plasmodium species, uniqueness of P. falciparum in SSR organization and general pattern of evolution of SSRs in Plasmodium.

Repeat performance: how do genome packaging and regulation depend on simple sequence repeats?

Non-coding DNA has consistently increased during evolution of higher eukaryotes. Since the number of genes has remained relatively static during the evolution of complex organisms, it is believed that increased degree of sophisticated regulation of genes has contributed to the increased complexity. A higher proportion of non-coding DNA, including repeats, is likely to provide more complex regulatory potential. Here, we propose that repeats play a regulatory role by contributing to the packaging of the genome during cellular differentiation. Repeats, and in particular the simple sequence repeats, are proposed to serve as landmarks that can target regulatory mechanisms to a large number of genomic sites with the help of very few factors and regulate the linked loci in a coordinated manner. Repeats may, therefore, function as common target sites for regulatory mechanisms involved in the packaging and dynamic compartmentalization of the chromatin into active and inactive regions during cellular differentiation.

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.