Junker: An Intergenic Explorer for Bacterial Genomes

In the past few decades, scientists from all over the world have taken a keen interest in novel functional units such as small regulatory RNAs, small open reading frames, pseudogenes, transposons, integrase binding attB/attP sites, repeat elements within the bacterial intergenic regions (IGRs) and in the analysis of those junk regions for ge- nomic complexity. Here we have developed a web server, named Junker, to facilitate the in-depth analysis of IGRs for examining their length distribution, four-quadrant...     

BBProF: An Asynchronous Application Server for Rapid Identification of Proteins Associated with Bacterial Bioleaching

Recovery of metal value, especially from low-grade ores and overburden minerals using acidophilic bacteria through the process of bioleaching is an environmentally benign and commercially scalable biotechnology. In recent years, while the “OMICS” landscape has been witnessing extensive application of computational tools to understand and interpret global biological sequence data, a dedicated bioinformatic server for analysis of bacterial information in the context of its bioleaching ability is not available. We have developed an on-line Bacterial Bioleaching Protein Finder (BBProF) System, which rapidly identifies novel proteins involved in a bacterial bioleaching process and also performs phylogenetic analysis of 16S rRNA genes. BBProF uses the features of Asynchronous Java Script and XML (AJAX) to provide an efficient and fast user experience with minimal requirement of network bandwidth. In the input module the server accepts any bacterial or archaeal complete genome sequence in RAW format and provides a list of proteins involved in the microbial leaching process. BBProF web server is integrated with the European Bioinformatics Institute (EBI) web services such as BLAST for homology search and InterProScan for functional characterization of output protein sequences. Studying evolutionary relationship of bacterial strains of interest using Muscle and ClustalW2 phylogeny web services from EBI is another key feature of our server, where 16S rRNA gene sequences are considered as input through a JQUERY interface along with the sequences present in the BBProF database library. Complete genome sequences of 24 bioleaching microorganism characterized by genomic and physiological study in the laboratory and their respective 16S rRNA gene sequences were stored in the database of the BBProF library. To our knowledge BBProF is the first integrated bioinformatic web server that demonstrates its utility in identifying potential bioleaching bacteria. We hope that the server facilitate ongoing comparative genomic studies on of bioleaching microorganisms and also assist in identification and design of novel microbial consortia that are optimally efficient bioleaching agents.     

ParIS Genome Rearrangement server

SUMMARY: ParIS Genome Rearrangement is a web server for a Bayesian analysis of unichromosomal genome pairs. The underlying model allows inversions, transpositions and inverted transpositions. The server generates a Markov chain using a Partial Importance Sampler technique, and samples trajectories of mutations from this chain. The user can specify several marginalizations to the posterior: the posterior distribution of number of mutations needed to transform one genome into another, length distribution of mutations, number of mutations that have occurred at a given site. Both text and graphical outputs are available. We provide a limited server, a downloadable unlimited server that can be installed locally on any linux/Unix operating system, and a database of mitochondrial gene orders.     

Complete nucleotide sequence of the sugarcane (Saccharum officinarum) chloroplast genome: a comparative analysis of four monocot chloroplast genomes.

The complete nucleotide sequence of the chloroplast genome of sugarcane (Saccharum officinarum) has been determined. It is a circular double-stranded DNA molecule, 141,182 bp in size, and is composed of a large single copy of 83,048 bp, a small single copy of 12,544 bp, and a pair of inverted repeat regions of 22,795 bp each. A comparative analysis among monocots showed that the sugarcane chloroplast genome was very similar to maize but not to rice or wheat. Between sugarcane and maize at the rps16-trnQ (UUG) region, however, a length polymorphism was identified. With regard to insertions/deletions equal to or longer than 5 bp, a total of 53 insertion and 31 deletion events were identified in the sugarcane chloroplast genome. Of the 84 loci identified, a pair of direct repeat sequences was located side by side in a tandem fashion in 47 loci (56.0%). A recombination event during plant evolution is discussed at two sites between the sugarcane and tobacco chloroplast genomes.     

Novel small RNA-encoding genes in the intergenic regions of Bacillus subtilis

Small, non-coding RNAs (ncRNAs) perform diverse functions in a variety of organisms, but few ncRNAs have been identified in Bacillus subtilis. To search the B. subtilis genome for genes encoding ncRNAs, we focused on 123 intergenic regions (IGRs) over 500 bp in length and analyzed expression from these regions. Seven IGRs termed bsrC, bsrD, bsrE, bsrF, bsrG, bsrH and bsrI expressed RNAs smaller than 380 nt. All small RNAs except BsrD RNA were expressed in transformed Escherichia coli cells harboring a plasmid with PCR-amplified IGRs of B. subtilis, indicating that their own promoters independently express small RNAs. Under the non-stressed condition, depletion of the genes for the small RNAs did not affect growth. Although their functions are unknown, gene expression profiles at several time points showed that most of the genes except for bsrD were expressed during the vegetative phase (4-6 h), but undetectable during the stationary phase (8 h). Mapping the 5' ends of the 6 small RNAs revealed that the genes for BsrE, BsrF, BsrG, BsrH, and BsrI RNAs are preceded by a recognition site for RNA polymerase sigma factor sigma(A). These small RNAs might lack an SD sequence and exert their actions as ncRNAs.     

Identification of prophages in bacterial genomes by dinucleotide relative abundance difference

Background

Prophages are integrated viral forms in bacterial genomes that have been found to contribute to interstrain genetic variability. Many virulence-associated genes are reported to be prophage encoded. Present computational methods to detect prophages are either by identifying possible essential proteins such as integrases or by an extension of this technique, which involves identifying a region containing proteins similar to those occurring in prophages. These methods suffer due to the problem of low sequence similarity at the protein level, which suggests that a nucleotide based approach could be useful.

Methodology

Earlier dinucleotide relative abundance (DRA) have been used to identify regions, which deviate from the neighborhood areas, in genomes. We have used the difference in the dinucleotide relative abundance (DRAD) between the bacterial and prophage DNA to aid location of DNA stretches that could be of prophage origin in bacterial genomes. Prophage sequences which deviate from bacterial regions in their dinucleotide frequencies are detected by scanning bacterial genome sequences. The method was validated using a subset of genomes with prophage data from literature reports. A web interface for prophage scan based on this method is available at http://bicmku.in:8082/prophagedb/dra.html. Two hundred bacterial genomes which do not have annotated prophages have been scanned for prophage regions using this method.

Conclusions

The relative dinucleotide distribution difference helps detect prophage regions in genome sequences. The usefulness of this method is seen in the identification of 461 highly probable loci pertaining to prophages which have not been annotated so earlier. This work emphasizes the need to extend the efforts to detect and annotate prophage elements in genome sequences.     

MRD: a microsatellite repeats database for prokaryotic and eukaryotic genomes

MRD is a database system to access the microsatellite repeats information of genomes such as archea, eubacteria, and other eukaryotic genomes whose sequence information is available in public domains. MRD stores information about simple tandemly repeated k-mer sequences where k= 1 to 6, i.e. monomer to hexamer. The web interface allows the users to search for the repeat of their interest and to know about the association of the repeat with genes and genomic regions in the specific organism. The data contains the abundance and distribution of microsatellites in the coding and non-coding regions of the genome. The exact location of repeats with respect to genomic regions of interest (such as UTR, exon, intron or intergenic regions) whichever is applicable to organism is highlighted. MRD is available on the World Wide Web at and/or . The database is designed as an open-ended system to accommodate the microsatellite repeats information of other genomes whose complete sequences will be available in future through public domain.     

Tandem repeat regions within the Burkholderia pseudomallei genome and their application for high resolution genotyping

Background

The facultative, intracellular bacterium Burkholderia pseudomallei is the causative agent of melioidosis, a serious infectious disease of humans and animals. We identified and categorized tandem repeat arrays and their distribution throughout the genome of B. pseudomallei strain K96243 in order to develop a genetic typing method for B. pseudomallei. We then screened 104 of the potentially polymorphic loci across a diverse panel of 31 isolates including B. pseudomallei, B. mallei and B. thailandensis in order to identify loci with varying degrees of polymorphism. A subset of these tandem repeat arrays were subsequently developed into a multiple-locus VNTR analysis to examine 66 B. pseudomallei and 21 B. mallei isolates from around the world, as well as 95 lineages from a serial transfer experiment encompassing ~18,000 generations.

Results

B. pseudomallei contains a preponderance of tandem repeat loci throughout its genome, many of which are duplicated elsewhere in the genome. The majority of these loci are composed of repeat motif lengths of 6 to 9 bp with 4 to 10 repeat units and are predominately located in intergenic regions of the genome. Across geographically diverse B. pseudomallei and B.mallei isolates, the 32 VNTR loci displayed between 7 and 28 alleles, with Nei's diversity values ranging from 0.47 and 0.94. Mutation rates for these loci are comparable (>10^-5 per locus per generation) to that of the most diverse tandemly repeated regions found in other less diverse bacteria.

Conclusion

The frequency, location and duplicate nature of tandemly repeated regions within the B. pseudomallei genome indicate that these tandem repeat regions may play a role in generating and maintaining adaptive genomic variation. Multiple-locus VNTR analysis revealed extensive diversity within the global isolate set containing B. pseudomallei and B. mallei, and it detected genotypic differences within clonal lineages of both species that were identical using previous typing methods. Given the health threat to humans and livestock and the potential for B. pseudomallei to be released intentionally, MLVA could prove to be an important tool for fine-scale epidemiological or forensic tracking of this increasingly important environmental pathogen.     

Relationships Among Powered Flight,Metabolic Rate,Body Mass,Genome Size,and the Retrotransposon Complement of Volant Birds

Avian genomes are of interest because the rapid metabolic rate associated with powered flight requires small cells which constrain genome size. Consequently, flying birds tend to have small genomes relative to other vertebrates such as mammals. It thus stands to reason that flying birds should have smaller genomes than ground-dwelling birds with lower metabolic rates. Small genomes could be condensed but uncompromised in a number of ways, including smaller intergenic intervals, shorter introns, and/or a reduced transposable element (TE) complement. We evaluated genome size in light of the orthologous TE complement among 41 flying (FY) and seven ground-dwelling (GD) bird species to determine if a preponderance of deletions in orthologous TEs might explain the compact genomes of flying birds with high metabolic rates. We measured, across multiple loci in all 48 species, the lengths of 50 contemporary orthologous chicken repeat 1 (CR1, a non-LTR retrotransposon) copies relative to inferred ancestral CR1 sequences. We found genome sizes in GD birds were not different than those in FY birds, but the mean lengths of orthologous CR1 loci were significantly shorter in FY birds than in GD birds. Moreover, we observed a negative correlation between basal metabolic rate and length of orthologous CR1 loci. Finally, we observed positive correlations between body mass and both genome sizes as well as length of orthologous CR1 loci, which we expected given that body mass correlates negatively with metabolic rates. Our results support the contention that metabolism helps shape the avian TE complement and thus indirectly contributes to the compact genomes of birds.     

Genomic heterogeneity in the density of noncoding single-nucleotide and microsatellite polymorphisms in Plasmodium falciparum

The density and distribution of single-nucleotide polymorphisms (SNPs) across the genome has important implications for linkage disequilibrium mapping and association studies, and the level of simple-sequence microsatellite polymorphisms has important implications for the use of oligonucleotide hybridization methods to genotype SNPs. To assess the density of these types of polymorphisms in P. falciparum, we sampled introns and noncoding DNA upstream and downstream of coding regions among a variety of geographically diverse parasites. Across 36,229 base pairs of noncoding sequence representing 41 genetic loci, a total of 307 polymorphisms including 248 polymorphic microsatellites and 39 SNPs were identified. We found a significant excess of microsatellite polymorphisms having a repeat unit length of one or two, compared to those with longer repeat lengths, as well as a nonrandom distribution of SNP polymorphisms. Almost half of the SNPs localized to only three of the 41 genetic loci sampled. Furthermore, we find significant differences in the frequency of polymorphisms across the two chromosomes (2 and 3) examined most extensively, with an excess of SNPs and a surplus of polymorphic microsatellites on chromosome 3 as compared to chromosome 2 (P=0.0001). Furthermore, at some individual genetic loci we also find a nonrandom distribution of polymorphisms between coding and flanking noncoding sequences, where completely monomorphic regions may flank highly polymorphic genes. These data, combined with our previous findings of nonrandom distribution of SNPs across chromosome 2, suggest that the Plasmodium falciparum genome may be a mosaic with regard to genetic diversity, containing chromosomal regions that are highly polymorphic interspersed with regions that are much less polymorphic.     

Native homing endonucleases can target conserved genes in humans and in animal models

In recent years, both homing endonucleases (HEases) and zinc-finger nucleases (ZFNs) have been engineered and selected for the targeting of desired human loci for gene therapy. However, enzyme engineering is lengthy and expensive and the off-target effect of the manufactured endonucleases is difficult to predict. Moreover, enzymes selected to cleave a human DNA locus may not cleave the homologous locus in the genome of animal models because of sequence divergence, thus hampering attempts to assess the in vivo efficacy and safety of any engineered enzyme prior to its application in human trials. Here, we show that naturally occurring HEases can be found, that cleave desirable human targets. Some of these enzymes are also shown to cleave the homologous sequence in the genome of animal models. In addition, the distribution of off-target effects may be more predictable for native HEases. Based on our experimental observations, we present the HomeBase algorithm, database and web server that allow a high-throughput computational search and assignment of HEases for the targeting of specific loci in the human and other genomes. We validate experimentally the predicted target specificity of candidate fungal, bacterial and archaeal HEases using cell free, yeast and archaeal assays.     

A comparative map of bovine chromosome 19 based on a combination of mapping on a bacterial artificial chromosome scaffold map, a whole genome radiation hybrid panel and the human draft sequence

We have constructed a medium density physical map of bovine chromosome 19 using a combination of mapping loci on both a bovine bacterial artificial chromosome (BAC) scaffold map and a whole genome radiation hybrid (WGRH) panel. The resulting map contains 70 loci spanning the length of bovine chromosome 19. Three contiguous groups of BACs were identified on the basis of multiple loci mapping to individual BAC clones. Bovine chromosome 19 was found in this study to be comprised almost entirely from regions of human chromosome 17, with a small region putatively assigned to human chromosome 10. Fourteen breakpoints between the bovine and human chromosomes were detected, with a possibility of five more based on ordering of the WGRH map.     

Expansion of GAA trinucleotide repeats in mammals

We have previously shown that GAA trinucleotide repeats have undergone significant expansion in the human genome. Here we present the analysis of the length distribution of all 10 nonredundant trinucleotide repeat motifs in 20 complete eukaryotic genomes (6 mammalian, 2 nonmammalian vertebrates, 4 arthropods, 4 fungi, and 1 each of nematode, amoebozoa, alveolate, and plant), which showed that the abundance of large expansions of GAA trinucleotide repeats is specific to mammals. Analysis of human-chimpanzee-gorilla orthologs revealed that loci with large expansions are species-specific and have occurred after divergence from the common ancestor. PCR analysis of human controls revealed large expansions at multiple human (GAA)(30+) loci; nine loci showed expanded alleles containing >65 triplets, analogous to disease-causing expansions in Friedreich ataxia, including two that are in introns of genes of unknown function. The abundance of long GAA trinucleotide repeat tracts in mammalian genomes represents a significant mutation potential and source of interindividual variability.     

Intrinsic polymorphism of variable number tandem repeat loci in the human genome.

In the human genome, short tandem repetitive (STR) DNA sequences often show restriction fragment length polymorphisms (RFLPs) due to variation in the number of copies of the repeat unit. For a subset of these sequences known as minisatellites or variable number tandem repeat loci (VNTR), it has been proposed that a homologous "core" sequence of 10-12 nucleotides is involved in the mechanism(s) generating the polymorphism. In our present study we have prepared oligonucleotide probes complementary to one or two repeat units of several VNTR loci. Under stringent hybridization and wash conditions these probes hybridize locus specifically thus allowing the evaluation of the intrinsic polymorphism of individual loci. Our results indicate that not all of the loci having STR DNA sequences are polymorphic despite the fact that they share the "core" sequence. This suggests that more than the DNA sequence of the locus is involved in the mechanism(s) generating the polymorphism.     

A hypervariable microsatellite revealed by in vitro amplification of a dinucleotide repeat within the cardiac muscle actin gene.

The human genome contains approximately 50,000 copies of an interspersed repeat with the sequence (dT-dG)n, where n = approximately 10-60. In humans, (TG)n repeats have been found in several sequenced regions. Since minisatellite regions with larger repeat elements often display extensive length polymorphisms, we suspected that (TG)n repeats ("microsatellites") might also be polymorphic. Using the polymerase chain reaction to amplify a (TG)n microsatellite in the human cardiac actin gene, we detected 12 different allelic fragments in 37 unrelated individuals, 32 of whom were heterozygous. Codominant Mendelian inheritance of fragments was observed in three families with a total of 24 children. Because of the widespread distribution of (TG)n microsatellites, polymorphisms of this type may be generally abundant and present in regions where minisatellites are rare, making such microsatellite loci very useful for linkage studies in humans.