Comparative analysis of base correlations in 5' untranslated regions of various species

A simple shortcut to unsupervised alignment-free phylogenetic genome groupings, even from unassembled sequencing reads

We propose an extension to alignment-free approaches that can produce reasonably accurate phylogenetic groupings starting from unaligned genomes, for example, as fast as 1 min on a standard desktop computer for 25 bacterial genomes. A 6-fold speed-up and 11-fold reduction in memory requirements compared to previous alignment-free methods is achieved by reducing the comparison space to a representative sample of k-mers of optimal length and with specific tag motifs. This approach was applied to the test case of fitting the enterohemorrhagic O104:H4 E.coli strain from the 2011 outbreak in Germany into the phylogenetic network of previously known E.coli-related strains and extend the method to allow assigning any new strain to the correct phylogenetic group even directly from unassembled short sequence reads from next generation sequencing data. Hence, this approach is also useful to quickly identify the most suitable reference genome for subsequent assembly steps.

Selective constraint on the upstream open reading frames that overlap with coding sequences in animals

Upstream open reading frames (uORFs) are translational regulatory elements located in 5′ untranslated regions. They can significantly repress the translation of the downstream coding sequences (CDS), and participate in the spatio-temporal regulations of protein translation. Notwithstanding this biological significance, the selective constraint on uORFs remains underexplored. Particularly, the uORFs that partially overlap with CDS with a different reading frame (overlapping uORFs, or “VuORFs”) may lead to strong translational inhibition or N-terminal truncation of the peptides encoded by the affected CDS. By analyzing VuORF-containing transcripts (designated as “VuORF transcripts”) in human, mouse, and fruit fly, we demonstrate that VuORFs are in general slightly deleterious - the proportion of genes that encode at least one VuORF transcript is significantly smaller than expected in all of the three examined species. In addition, this proportion is significantly smaller in fruit fly than in mammals, indicating a higher efficiency of removing VuORFs in the former species because of its larger effective population size. Furthermore, the deleterious effect of a VuORF depends on the sequence context of its start codon (VuAUG). VuORFs with an optimal VuAUG context are more strongly disfavored than those with a suboptimal context in all of the three examined species. And the propensity to remove optimal-context VuAUGs is stronger in fruit fly than in mammals. Intriguingly, however, the currently observable optimal-context VuAUGs (but not suboptimal-context VuAUGs) are more conserved than expected. These observations suggest that the regulatory functions of VuORFs may have been gained fortuitously in organisms with a small effective population size because the slightly deleterious effect of these elements can be better tolerated in these organisms, thus allowing opportunities for the development of novel biological functions. Nevertheless, once the functions of VuORFs were established, they became subject to negative selection.

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.

Mapping of the transcription start site (TSS) and identification of SNPs in the bovine neuropeptide Y (NPY) gene

BACKGROUND

Neuropeptide Y is a key neurotransmitter of the central nervous system which plays a vital role in the feed energy homeostasis in mammals. Mutations in the regulatory and coding regions of the bovine NPY gene can potentially affect the neuronal regulation of appetite and feeding behaviour in cattle. The objectives of this experiment were to: a) fully characterize the bovine NPY gene transcript and b) identify the SNP diversity in both coding and non-coding regions of the NPY gene in a panel of Bos taurus and B. indicus cattle.

RESULTS

Bovine NPY gene consists of four exons (99, 188, 81 and 195 nucleotides) and three introns. The promoter region of the NPY gene consists of TATA and GC boxes which are separated from the transcription start site (TSS) by 29 and ~100 nt, respectively. Analyses of the tissue specific expression of the bovine NPY gene revealed the presence of highly abundant NPY gene transcripts in the arcuate nucleus, cerebral and cerebellar regions of the bovine brain. We identified a total of 59 SNPs in the 8.4 kb region of the bovine NPY gene. Seven out of nine total SNPs in the promoter region affect binding of putative transcription factors. A high level of nucleotide diversity was evident in the promoter regions (2.84 x 10(-3)) compared to the exonic (1.44 x 10(-3)), intronic (1.30 x 10(-3)) and 3' untranslated (1.26 x 10(-3)) regions.

CONCLUSION

The SNPs identified in different regions of bovine NPY gene may serve as a basis for understanding the regulation of the expression of the bovine NPY gene under a variety of physiological conditions and identification of genotypes with high feed energy conversion efficiency.

A mutation in the 5' untranslated region of the BRCA1 gene in sporadic breast cancer causes downregulation of translation efficiency

We screened 117 breast tumour samples in Chinese females for mutations in the breast cancer 1 (BRCA1) gene and identified a novel mutation in the 5' untranslated region (5' UTR) in two patients with grade III infiltrating ductal breast carcinoma. We examined whether this 5' UTR mutation affected the translational efficiency of BRCA1 protein. A vector was constructed containing the mutated 5' UTR up-stream of luciferase and we compared its translational efficiency with a wild-type 5' UTR. The expression of BRCA1 protein in breast tumour samples was evaluated using immunohistochemistry. The mutated 5' UTR of BRCA1 resulted in less luciferase activity compared with the wild-type 5' UTR, while there were no significant differences in luciferase mRNA levels. BRCA1 protein was much less expressed in breast tumour tissue from patients with the 5' UTR mutation than in samples from patients without the mutation. Our results show that a mutation in the 5' UTR of the BRCA1 gene downregulates translational efficiency of the protein.