Comparative analysis of core promoter region: information content from mono and dinucleotide substitution matrices

Information theoretic perspective on genome clustering

Shannon’s information theoretic perspective of communication helps one to understand the storage and processing of information in one-dimensional sequences. An information theoretic analysis of 937 available completely sequenced prokaryotic genomes and 238 eukaryotic chromosomes is presented. Information content (Id) values were used to cluster these chromosomes. Chargaff’s second parity rule i.e compositional self-complementarity, an empirical fact is observed in all the genomes, except for the proteobacteria Candidatus Hodgkinia cicadicola. High information content, arising out of biased base composition in all the 14 chromosomes of Plasmodium falciparum is found among two other genomes of prokaryotes viz. Buchnera aphidicola str. Cc (Cinara cedri) and Candidatus Carsonella ruddii PV. Despite size and compositional variations, both prokaryotic and eukaryotic genomes do not deviate significantly from an equiprobable and random situation. Eukaryotic chromosomes of an organism tend to have similar informational restraints as seen when a simple distance based method is used to cluster them. In eukaryotes, in certain cases, Id values are also similar for the two arms (p and q arm) of the chromosomes. The results of this current study confirm that the information content can provide insights into the clustering of genomes and the evolution of messaging strategies of the genomes. An efficient and robust Perl CGI standalone tool is created based on this information theory algorithm for the analysis of the whole genomes and is made available at https://github.com/AlagurajVeluchamy/InformationTheory.

In silico analysis of high affinity potassium transporter (HKT) isoforms in different plants

BACKGROUND

High affinity potassium transporters (HKTs) are located in the plasma membrane of the vessels and have significant influence on salt tolerance in some plants. They exclude Na(+) from the parenchyma cells to reduce Na(+) concentration. Despite many studies, the underlying regulatory mechanisms and the exact functions of HKTs within different genomic backgrounds are relatively unknown. In this study, various bioinformatics techniques, including promoter analysis, identification of HKT-surrounding genes, and construction of gene networks, were applied to investigate the HKT regulatory mechanism.

RESULTS

Promoter analysis showed that rice HKTs carry ABA response elements. Additionally, jasmonic acid response elements were detected on promoter region of TmHKT1;5. In silico synteny highlighted several unknown and new loci near rice, Arabidopsis thaliana and Physcomitrella patent HKTs, which may play a significant role in salt stress tolerance in concert with HKTs. Gene network prediction unravelled that crosstalk between jasmonate and ethylene reduces AtHKT1;1 expression. Furthermore, antiporter and transferase proteins were found in AtHKT1;1 gene network. Interestingly, regulatory elements on the promoter region of HKT in wild genotype (TmHKT1;5) were more frequent and variable than the ones in cultivated wheat (TaHKT1;5) which provides the possibility of rapid response and better understanding of environmental conditions for wild genotype.

CONCLUSION

Detecting ABA and jasmonic acid response elements on promoter regions of HKTs provide valuable clues on underlying regulatory mechanisms of HKTs. In silico synteny and pathway discovery indicated several candidates which act in concert with HKTs in stress condition. We highlighted different arrangement of regulatory elements on promoter region of wild wheat (TmHKT1;5) compared to bread wheat (TaHKT1;5) in this study.

Identifying statistical dependence in genomic sequences via mutual information estimates

Questions of understanding and quantifying the representation and amount of information in organisms have become a central part of biological research, as they potentially hold the key to fundamental advances. In this paper, we demonstrate the use of information-theoretic tools for the task of identifying segments of biomolecules (DNA or RNA) that are statistically correlated. We develop a precise and reliable methodology, based on the notion of mutual information, for finding and extracting statistical as well as structural dependencies. A simple threshold function is defined, and its use in quantifying the level of significance of dependencies between biological segments is explored. These tools are used in two specific applications. First, they are used for the identification of correlations between different parts of the maize zmSRp32 gene. There, we find significant dependencies between the 5' untranslated region in zmSRp32 and its alternatively spliced exons. This observation may indicate the presence of as-yet unknown alternative splicing mechanisms or structural scaffolds. Second, using data from the FBI's combined DNA index system (CODIS), we demonstrate that our approach is particularly well suited for the problem of discovering short tandem repeats-an application of importance in genetic profiling.

Comparative analysis of splice site regions by information content

We have applied concepts from information theory for a comparative analysis of donor (gt) and acceptor (ag) splice site regions in the genes of five different organisms by calculating their mutual information content (relative entropy) over a selected block of nucleotides. A similar pattern that the information content decreases as the block size increases was observed for both regions in all the organisms studied. This result suggests that the information required for splicing might be contained in the consensus of ~6–8 nt at both regions. We assume from our study that even though the nucleotides are showing some degrees of conservation in the flanking regions of the splice sites, certain level of variability is still tolerated, which leads the splicing process to occur normally even if the extent of base pairing is not fully satisfied. We also suggest that this variability can be compensated by recognizing different splice sites with different spliceosomal factors.

Comparative analysis of transcription start sites using mutual information

The transcription start site (TSS) region shows greater variability compared with other promoter elements. We are interested to search for its variability by using information content as a measure. We note in this study that the variability is significant in the block of 5 nucleotides (nt) surrounding the TSS region compared with the block of 15 nt. This suggests that the actual region that may be involved is in the range of 5-10 nt in size. For Escherichia coli, we note that the information content from dinucleotide substitution matrices clearly shows a better discrimination, suggesting the presence of some correlations. However, for human this effect is much less, and for mouse it is practically absent. We can conclude that the presence of short-range correlations within the TSS region is species-dependent and is not universal. We further observe that there are other variable regions in the mitochondrial control element apart from TSS. It is also noted that effective comparisons can only be made on blocks, while single nucleotide comparisons do not give us any detectable signals.