UTRdb and UTRsite: specialized databases of sequences and functional elements of 5' and 3' untranslated regions of eukaryotic mRNAs. Update 2002

An in-silico study to determine susceptibility to cancer by evaluating the coding and non-coding non-synonymous single nucleotide variants in the SOCS3 gene

Single Nucleotide Variant (SNVs) affect gene expression as well as protein structure and activity, leading to reduced signaling capabilities and ultimately, increasing cancer risk. SOCS3 (suppressor of cytokine signaling 3), a critical tumor suppressor providing a substantial part in the feedback loop of the JAK/STAT pathway, is abnormally suppressed in various cancer. This study aims to screen non-coding and potentially deleterious coding SNVs in the SOCS3 gene. We performed six programs: PredictSNP1.0 (predicting Deleterious nsSNVs), ConSurf (analyzing sequence conservation), ModPred (analyzing SNVS in PTMs sites), I-Mutant and MUpro (to analyze SNVs effecting protein stability), and molecular docking and molecular dynamics (MD) (to assess the consequences of SOCS3 genetic variations on JAK interactions) for coding regions and three programs (UTRSite, SNP2TFBS, miRNA SNP) (to analyze SNVs effecting the gene expression) in non-coding regions, respectively. After screening 2786 SOCS3 SNVs, we found 10 SNVs, as well as 49 SNPs that change the function of non-coding areas. Out of 10 selected nsSNVs, 3 SNVs (W48R, R71C, N198S) predicted to be the most damaging by all the software programs, as well as one nsSNV (R194W) could be highly deleterious from Molecular Docking analysis combined with MD Simulations. Our findings propose a procedure for studying the structure-related consequences of SNVs on protein function in the future.Communicated by Ramaswamy H. Sarma.

Functional characterization of human genomic variation linked to polygenic diseases

The burden of human disease lies predominantly in polygenic diseases. Since the early 2000s, genome-wide association studies (GWAS) have identified genetic variants and loci associated with complex traits. These have ranged from variants in coding sequences to mutations in regulatory regions, such as promoters and enhancers, as well as mutations affecting mediators of mRNA stability and other downstream regulators, such as 5' and 3'-untranslated regions (UTRs), long noncoding RNA (lncRNA), and miRNA. Recent research advances in genetics have utilized a combination of computational techniques, high-throughput in vitro and in vivo screening modalities, and precise genome editing to impute the function of diverse classes of genetic variants identified through GWAS. In this review, we highlight the vastness of genomic variants associated with polygenic disease risk and address recent advances in how genetic tools can be used to functionally characterize them.

3'UTR Diversity: Expanding Repertoire of RNA Alterations in Human mRNAs

Genomic information stored in the DNA is transcribed to the mRNA and translated to proteins. The 3' untranslated regions (3'UTRs) of the mRNA serve pivotal roles in posttranscriptional gene expression, regulating mRNA stability, translation, and localization. Similar to DNA mutations producing aberrant proteins, RNA alterations expand the transcriptome landscape and change the cellular proteome. Recent global analyses reveal that many genes express various forms of altered RNAs, including 3'UTR length variants. Alternative polyadenylation and alternative splicing are involved in diversifying 3'UTRs, which could act as a hidden layer of eukaryotic gene expression control. In this review, we summarize the functions and regulations of 3'UTRs and elaborate on the generation and functional consequences of 3'UTR diversity. Given that dynamic 3'UTR length control contributes to phenotypic complexity, dysregulated 3'UTR diversity might be relevant to disease development, including cancers. Thus, 3'UTR diversity in cancer could open exciting new research areas and provide avenues for novel cancer theragnostics.

Association of rs13429458 and rs12478601 Single Nucleotide Polymorphisms of THADA Gene with Polycystic Ovary Syndrome

BACKGROUND

It is thought that genetic factors are influential in the etiology of polycystic ovarian syndrome (PCOS), the most frequent endocrinological disorder of females in their reproductive age. This study was carried out to elucidate the association of rs13429458 and rs12478601 single nucleotide polymorphisms (SNPs) of the THADA gene and the risk of the PCOS among a population of Iranian female patients.

MATERIALS AND METHODS

This case-control study contains 66 infertile women with PCOS (patient group) and 44 healthy women without PCOS (control group) that referred to the IVF Unit of the Infertility Research Center of the Academic Center for Education, Culture and Research (ACECR). The polymerase chain reaction (PCR) was utilized to amplify genome DNA as well as direct sequencing to determine SNPs. The THADA rs12478601 and rs13429458 genotypes were consequently examined with amplification refractory mutation system-PCR (ARMS-PCR).

RESULTS

In this study, we observed that rs13429458 polymorphism was not associated with PCOS risk in two groups (P=0.42). On the other hand, data analysis indicated that the rs12478601 genotype significantly increased the risk of PCOS in the case group (P=0.032) in compared with control group. We found that the "T" allele of rs12478601 in the THADA gene had a significant relation to PCOS in the case group (odds ratio [OR]: 2.574, 95% confidence interval [CI]: 1.439-4.604, P=0.001).

CONCLUSION

This study has presented further evidence that TT and CT genotype of THADA rs12478601 is associated with a high risk of PCOS.

Computational Detection of MicroRNA Targets

MicroRNAs (miRNAs) are small noncoding RNAs that are recognized as posttranscriptional regulators of gene expression. These molecules have been shown to play important roles in several cellular processes. MiRNAs act on their target by guiding the RISC complex and binding to the mRNA molecule. Thus, it is recognized that the function of a miRNA is determined by the function of its target (s). By using high-throughput methodologies, novel miRNAs are being identified, but their functions remain uncharted. Target validation is crucial to properly understand the specific role of a miRNA in a cellular pathway. However, molecular techniques for experimental validation of miRNA-target interaction are expensive, time-consuming, laborious, and can be not accurate in inferring true interactions. Thus, accurate miRNA target predictions are helpful to understand the functions of miRNAs. There are several algorithms proposed for target prediction and databases containing miRNA-target information. However, these available computational tools for prediction still generate a large number of false positives and fail to detect a considerable number of true targets, which indicates the necessity of highly confident approaches to identify bona fide miRNA-target interactions. This chapter focuses on tools and strategies used for miRNA target prediction, by providing practical insights and outlooks.

Peculiarities of Plasmodium falciparum Gene Regulation and Chromatin Structure

The highly complex life cycle of the human malaria parasite, Plasmodium falciparum, is based on an orchestrated and tightly regulated gene expression program. In general, eukaryotic transcription regulation is determined by a combination of sequence-specific transcription factors binding to regulatory DNA elements and the packaging of DNA into chromatin as an additional layer. The accessibility of regulatory DNA elements is controlled by the nucleosome occupancy and changes of their positions by an active process called nucleosome remodeling. These epigenetic mechanisms are poorly explored in P. falciparum. The parasite genome is characterized by an extraordinarily high AT-content and the distinct architecture of functional elements, and chromatin-related proteins also exhibit high sequence divergence compared to other eukaryotes. Together with the distinct biochemical properties of nucleosomes, these features suggest substantial differences in chromatin-dependent regulation. Here, we highlight the peculiarities of epigenetic mechanisms in P. falciparum, addressing chromatin structure and dynamics with respect to their impact on transcriptional control. We focus on the specialized chromatin remodeling enzymes and discuss their essential function in P. falciparum gene regulation.

RNAs as Regulators of Cellular Matchmaking

RNA molecules are increasingly being identified as facilitating or impeding the interaction of proteins and nucleic acids, serving as so-called scaffolds or decoys. Long non-coding RNAs have been commonly implicated in such roles, particularly in the regulation of nuclear processes including chromosome topology, regulation of chromatin state and gene transcription, and assembly of nuclear biomolecular condensates such as paraspeckles. Recently, an increased awareness of cytoplasmic RNA scaffolds and decoys has begun to emerge, including the identification of non-coding regions of mRNAs that can also function in a scaffold-like manner to regulate interactions of nascently translated proteins. Collectively, cytoplasmic RNA scaffolds and decoys are now implicated in processes such as mRNA translation, decay, protein localization, protein degradation and assembly of cytoplasmic biomolecular condensates such as P-bodies. Here, we review examples of RNA scaffolds and decoys in both the nucleus and cytoplasm, illustrating common themes, the suitability of RNA to such roles, and future challenges in identifying and better understanding RNA scaffolding and decoy functions.

Exchange of endogenous and heterogeneous yeast terminators in Pichia pastoris to tune mRNA stability and gene expression

In the yeast Saccharomyces cerevisiae, terminator sequences not only terminate transcription but also affect expression levels of the protein-encoded upstream of the terminator. The non-conventional yeast Pichia pastoris (syn. Komagataella phaffii) has frequently been used as a platform for metabolic engineering but knowledge regarding P. pastoris terminators is limited. To explore terminator sequences available to tune protein expression levels in P. pastoris, we created a 'terminator catalog' by testing 72 sequences, including terminators from S. cerevisiae or P. pastoris and synthetic terminators. Altogether, we found that the terminators have a tunable range of 17-fold. We also found that S. cerevisiae terminator sequences maintain function when transferred to P. pastoris. Successful tuning of protein expression levels was shown not only for the reporter gene used to define the catalog but also using betaxanthin production as an example application in pathway flux regulation. Moreover, we found experimental evidence that protein expression levels result from mRNA abundance and in silico evidence that levels reflect the stability of mRNA 3'-UTR secondary structure. In combination with promoter selection, the novel terminator catalog constitutes a basic toolbox for tuning protein expression levels in metabolic engineering and synthetic biology in P. pastoris.

Bacterial 3'UTRs: A Useful Resource in Post-transcriptional Regulation

Bacterial messenger RNAs (mRNAs) are composed of 5′ and 3′ untranslated regions (UTRs) that flank the coding sequences (CDSs). In eukaryotes, 3′UTRs play key roles in post-transcriptional regulatory mechanisms. Shortening or deregulation of these regions is associated with diseases such as cancer and metabolic disorders. Comparatively, little is known about the functions of 3′UTRs in bacteria. Over the past few years, 3′UTRs have emerged as important players in the regulation of relevant bacterial processes such as virulence, iron metabolism, and biofilm formation. This MiniReview is an update for the different 3′UTR-mediated mechanisms that regulate gene expression in bacteria. Some of these include 3′UTRs that interact with the 5′UTR of the same transcript to modulate translation, 3′UTRs that are targeted by specific ribonucleases, RNA-binding proteins and small RNAs (sRNAs), and 3′UTRs that act as reservoirs of trans-acting sRNAs, among others. In addition, recent findings regarding a differential evolution of bacterial 3′UTRs and its impact in the species-specific expression of orthologous genes are also discussed.

Regulation of ribosomal protein genes: An ordered anarchy

Ribosomal protein genes are among the most highly expressed genes in most cell types. Their products are generally essential for ribosome synthesis, which is the cornerstone for cell growth and proliferation. Many cellular resources are dedicated to producing ribosomal proteins and thus this process needs to be regulated in ways that carefully balance the supply of nascent ribosomal proteins with the demand for new ribosomes. Ribosomal protein genes have classically been viewed as a uniform interconnected regulon regulated in eukaryotic cells by target of rapamycin and protein kinase A pathway in response to changes in growth conditions and/or cellular status. However, recent literature depicts a more complex picture in which the amount of ribosomal proteins produced varies between genes in response to two overlapping regulatory circuits. The first includes the classical general ribosome‐producing program and the second is a gene‐specific feature responsible for fine‐tuning the amount of ribosomal proteins produced from each individual ribosomal gene. Unlike the general pathway that is mainly controlled at the level of transcription and translation, this specific regulation of ribosomal protein genes is largely achieved through changes in pre‐mRNA splicing efficiency and mRNA stability. By combining general and specific regulation, the cell can coordinate ribosome production, while allowing functional specialization and diversity. Here we review the many ways ribosomal protein genes are regulated, with special focus on the emerging role of posttranscriptional regulatory events in fine‐tuning the expression of ribosomal protein genes and its role in controlling the potential variation in ribosome functions.

This article is categorized under:

Translation > Ribosome Biogenesis

Translation > Ribosome Structure/Function

Translation > Translation Regulation

Comprehensive analysis of long noncoding RNA expression in dorsal root ganglion reveals cell-type specificity and dysregulation after nerve injury

Dorsal root ganglion (DRG) neurons provide connectivity between peripheral tissues and the spinal cord. Transcriptional plasticity within DRG sensory neurons after peripheral nerve injury contributes to nerve repair but also leads to maladaptive plasticity, including the development of neuropathic pain. This study presents tissue and neuron-specific expression profiling of both known and novel long noncoding RNAs (LncRNAs) in the rodent DRG after nerve injury. We have identified a large number of novel LncRNAs expressed within the rodent DRG, a minority of which were syntenically conserved between the mouse, rat, and human, and including, both intergenic and antisense LncRNAs. We have also identified neuron type-specific LncRNAs in the mouse DRG and LncRNAs that are expressed in human IPS cell-derived sensory neurons. We show significant plasticity in LncRNA expression after nerve injury, which in mice is strain and gender dependent. This resource is publicly available and will aid future studies of DRG neuron identity and the transcriptional landscape in both the naive and injured DRG. We present our work regarding novel antisense and intergenic LncRNAs as an online searchable database, accessible from PainNetworks (http://www.painnetworks.org/). We have also integrated all annotated gene expression data in PainNetworks, so they can be examined in the context of their protein interactions.

DCDE: An Efficient Deep Convolutional Divergence Encoding Method for Human Promoter Recognition

Efficient human promoter feature extraction is still a major challenge in genome analysis as it can better understand human gene regulation and will be useful for experimental guidance. Although many machine learning algorithms have been developed for eukaryotic gene recognition, performance on promoters is unsatisfactory due to the diverse nature. To extract discriminative features from human promoters, an efficient deep convolutional divergence encoding method (DCDE) is proposed based on statistical divergence (SD) and convolutional neural network (CNN). SD can help optimize kmer feature extraction for human promoters. CNN can also be used to automatically extract features in gene analysis. In DCDE, we first perform informative kmers settlement to encode original gene sequences. A series of SD methods can optimize the most discriminative kmers distributions while maintaining important positional information. Then, CNN is utilized to extract lower dimensional deep features by secondary encoding. Finally, we construct a hybrid recognition architecture with multiple support vector machines and a bilayer decision method. It is flexible to add new features or new models and can be extended to identify other genomic functional elements. The extensive experiments demonstrate that DCDE is effective in promoter encoding and can significantly improve the performance of promoter recognition.

In Silico Analysis of Coding/Noncoding SNPs of Human RETN Gene and Characterization of Their Impact on Resistin Stability and Structure

Resistin (RETN) is a gene coding for proinflammatory adipokine called resistin secreted by macrophages in humans. Single nucleotide polymorphisms (SNPs) in RETN are linked to obesity and insulin resistance in various populations. Using dbSNP, 78 nonsynonymous SNPs (nsSNPs) were retrieved and tested on a PredictSNP 1.0 megaserver. Among these, 15 nsSNPs were predicted as highly deleterious and thus subjected to further analyses, such as conservation, posttranscriptional modifications, and stability. The 3D structure of human resistin was generated by homology modeling using Swiss model. Root-mean-square deviation (RMSD), hydrogen bonds (h-bonds), and interactions were estimated. Furthermore, UTRscan served to identify UTR functional SNPs. Among the 15 most deleterious nsSNPs, 13 were predicted to be highly conserved including variants in posttranslational modification sites. Stability analysis predicted 9 nsSNPs (I32S, C51Y, G58E, G58R, C78S, G79C, W98C, C103G, and C104Y) which can decrease protein stability with at least three out of the four algorithms used in this study. These nsSNPs were chosen for structural analysis. Both variants C51Y and C104Y showed the highest RMS deviations (1.137 Å and 1.308 Å, respectively) which were confirmed by the important decrease in total h-bonds. The analysis of hydrophobic and hydrophilic interactions showed important differences between the native protein and the 9 mutants, particularly I32S, G79C, and C104Y. Six SNPs in the 3'UTR (rs920569876, rs74176247, rs1447199134, rs943234785, rs76346269, and rs78048640) were predicted to be implicated in polyadenylation signal. This study revealed 9 highly deleterious SNPs located in the human RETN gene coding region and 6 SNPs within the 3'UTR that may alter the protein structure. Interestingly, these SNPs are worth to be analyzed in functional studies to further elucidate their effect on metabolic phenotype occurrence.

Synthetic Promoters: Designing the cis Regulatory Modules for Controlled Gene Expression

Designing the expression cassettes with desired properties remains the most important consideration of gene engineering technology. One of the challenges for predictive gene expression is the modeling of synthetic gene switches to regulate one or more target genes which would directly respond to specific chemical, environmental, and physiological stimuli. Assessment of natural promoter, high-throughput sequencing, and modern biotech inventory aided in deciphering the structure of cis elements and molding the native cis elements into desired synthetic promoter. Synthetic promoters which are molded by rearrangement of cis motifs can greatly benefit plant biotechnology applications. This review gives a glimpse of the manual in vivo gene regulation through synthetic promoters. It summarizes the integrative design strategy of synthetic promoters and enumerates five approaches for constructing synthetic promoters. Insights into the pattern of cis regulatory elements in the pursuit of desirable "gene switches" to date has also been reevaluated. Joint strategies of bioinformatics modeling and randomized biochemical synthesis are addressed in an effort to construct synthetic promoters for intricate gene regulation.

Toxicogenomic responses of low level anticancer drug exposures in Daphnia magna

The use of anticancer drugs in chemotherapy is increasing, leading to growing environmental concentrations of imatinib mesylate (IMA), cisplatinum (CDDP), and etoposide (ETP) in aquatic systems. Previous studies have shown that these anticancer drugs cause DNA damage in the crustacean Daphnia magna at low, environmentally relevant concentrations. To explore the mechanism of action of these compounds and the downstream effects of DNA damage on D. magna growth and development at a sensitive life stage, we exposed neonates to low level concentrations equivalent to those that elicit DNA damage (IMA: 2000 ng/L, ETP: 300 ng/L, CDDP: 10 ng/L) and performed transcriptomic analysis using an RNA-seq approach. RNA sequencing generated 14 million reads per sample, which were aligned to the D. magna genome and assembled, producing approximately 23,000 transcripts per sample. Over 90% of the transcripts showed homology to proteins in GenBank, revealing a high quality transcriptome assembly, although functional annotation was much lower. RT-qPCR was used to identify robust biomarkers and confirmed the downregulation of an angiotensin converting enzyme-like gene (ance) involved in neuropeptide regulation across all three anticancer drugs and the down-regulation of DNA topoisomerase II by ETP. RNA-seq analysis also allowed for an in depth exploration of the differential splicing of transcripts revealing that regulation of different gene isoforms predicts potential impacts on translation and protein expression, providing a more meaningful assessment of transcriptomic data. Enrichment analysis and investigation of affected biological processes suggested that the DNA damage caused by ETP and IMA influences cell cycle regulation and GPCR signaling. This dysregulation is likely responsible for effects to neurological system processes and development, and overall growth and development. Our transcriptomic approach provided insight into the mechanisms that respond to DNA damage caused by anticancer drug exposure and generated novel hypotheses on how these chemicals may impact the growth and survival of this ecologically important zooplankton species.

De novo computational identification of stress-related sequence motifs and microRNA target sites in untranslated regions of a plant translatome

Gene regulation at the transcriptional and translational level leads to diversity in phenotypes and function in organisms. Regulatory DNA or RNA sequence motifs adjacent to the gene coding sequence act as binding sites for proteins that in turn enable or disable expression of the gene. Whereas the known DNA and RNA binding proteins range in the thousands, only a few motifs have been examined. In this study, we have predicted putative regulatory motifs in groups of untranslated regions from genes regulated at the translational level in Arabidopsis thaliana under normal and stressed conditions. The test group of sequences was divided into random subgroups and subjected to three de novo motif finding algorithms (Seeder, Weeder and MEME). In addition to identifying sequence motifs, using an in silico tool we have predicted microRNA target sites in the 3′ UTRs of the translationally regulated genes, as well as identified upstream open reading frames located in the 5′ UTRs. Our bioinformatics strategy and the knowledge generated contribute to understanding gene regulation during stress, and can be applied to disease and stress resistant plant development.

SD-MSAEs: Promoter recognition in human genome based on deep feature extraction

The prediction and recognition of promoter in human genome play an important role in DNA sequence analysis. Entropy, in Shannon sense, of information theory is a multiple utility in bioinformatic details analysis. The relative entropy estimator methods based on statistical divergence (SD) are used to extract meaningful features to distinguish different regions of DNA sequences. In this paper, we choose context feature and use a set of methods of SD to select the most effective n-mers distinguishing promoter regions from other DNA regions in human genome. Extracted from the total possible combinations of n-mers, we can get four sparse distributions based on promoter and non-promoters training samples. The informative n-mers are selected by optimizing the differentiating extents of these distributions. Specially, we combine the advantage of statistical divergence and multiple sparse auto-encoders (MSAEs) in deep learning to extract deep feature for promoter recognition. And then we apply multiple SVMs and a decision model to construct a human promoter recognition method called SD-MSAEs. Framework is flexible that it can integrate new feature extraction or new classification models freely. Experimental results show that our method has high sensitivity and specificity.

Computational Identification of Post Translational Modification Regulated RNA Binding Protein Motifs

RNA and its associated RNA binding proteins (RBPs) mitigate a diverse array of cellular functions and phenotypes. The interactions between RNA and RBPs are implicated in many roles of biochemical processing by the cell such as localization, protein translation, and RNA stability. Recent discoveries of novel mechanisms that are of significant evolutionary advantage between RBPs and RNA include the interaction of the RBP with the 3’ and 5’ untranslated region (UTR) of target mRNA. These mechanisms are shown to function through interaction of a trans-factor (RBP) and a cis-regulatory element (3’ or 5’ UTR) by the binding of a RBP to a regulatory-consensus nucleic acid motif region that is conserved throughout evolution. Through signal transduction, regulatory RBPs are able to temporarily dissociate from their target sites on mRNAs and induce translation, typically through a post-translational modification (PTM). These small, regulatory motifs located in the UTR of mRNAs are subject to a loss-of-function due to single polymorphisms or other mutations that disrupt the motif and inhibit the ability to associate into the complex with RBPs. The identification of a consensus motif for a given RBP is difficult, time consuming, and requires a significant degree of experimentation to identify each motif-containing gene on a genomic scale. We have developed a computational algorithm to analyze high-throughput genomic arrays that contain differential binding induced by a PTM for a RBP of interest–RBP-PTM Target Scan (RPTS). We demonstrate the ability of this application to accurately predict a PTM-specific binding motif to an RBP that has no antibody capable of distinguishing the PTM of interest, negating the use of in-vitro exonuclease digestion techniques.

Comparative transcriptomics reveals striking similarities between the bovine and feline isolates of Tritrichomonas foetus: consequences for in silico drug-target identification

BACKGROUND

Few, if any, protozoan parasites are reported to exhibit extreme organ tropism like the flagellate Tritrichomonas foetus. In cattle, T. foetus infects the reproductive system causing abortion, whereas the infection in cats results in chronic large bowel diarrhoea. In the absence of a T. foetus genome, we utilized a de novo approach to assemble the transcriptome of the bovine and feline genotype to identify host-specific adaptations and virulence factors specific to each genotype. Furthermore, a subset of orthologs was used to characterize putative druggable targets and expose complications of in silico drug target mining in species with indefinite host-ranges.

RESULTS

Illumina RNA-seq reads were assembled into two representative bovine and feline transcriptomes containing 42,363 and 36,559 contigs, respectively. Coding and non-coding regions of the genome libraries revealed striking similarities, with 24,620 shared homolog pairs reduced down to 7,547 coding orthologs between the two genotypes. The transcriptomes were near identical in functional category distribution; with no indication of selective pressure acting on orthologs despite differences in parasite origins/host. Orthologs formed a large proportion of highly expressed transcripts in both genotypes (bovine genotype: 76%, feline genotype: 56%). Mining the libraries for protease virulence factors revealed the cysteine proteases (CP) to be the most common. In total, 483 and 445 bovine and feline T. foetus transcripts were identified as putative proteases based on MEROPS database, with 9 hits to putative protease inhibitors. In bovine T. foetus, CP8 is the preferentially transcribed CP while in the feline genotype, transcription of CP7 showed higher abundance. In silico druggability analysis of the two genotypes revealed that when host sequences are taken into account, drug targets are genotype-specific.

CONCLUSION

Gene discovery analysis based on RNA-seq data analysis revealed prominent similarities between the bovine and feline T. foetus, suggesting recent adaptation to their respective host/niche. T. foetus represents a unique case of a mammalian protozoan expanding its parasitic grasp across distantly related host lineages. Consequences of the host-range for in silico drug targeting are exposed here, demonstrating that targets of the parasite in one host are not necessarily ideal for the same parasite in another host.

CFTR mRNA expression is regulated by an upstream open reading frame and RNA secondary structure in its 5' untranslated region

Post-transcriptional regulation of gene expression through 5' untranslated region (5'UTR)-encoded cis-acting elements is an important mechanism for the control of protein expression levels. Through controlling specific aspects of translation initiation, expression can be tightly regulated while remaining responsive to cellular requirements. With respect to cystic fibrosis (CF), the overexpression of cystic fibrosis transmembrane conductance regulator (CFTR) protein trafficking mutants, such as delta-F508, is of great biological and clinical interest. By understanding the post-transcriptional mechanisms that regulate CFTR expression, new procedures can be developed to enhance CFTR expression in homozygous delta-F508 CF patients. We have identified the key elements of a complex negative regulatory mechanism that is encoded within the human CFTR 5'UTR and show how these elements act in combination to restrict CFTR gene expression to a consistently low level in a transcript-specific manner. This study shows, for the first time, that endogenous human CFTR expression is post-transcriptionally regulated through a 5'UTR-mediated mechanism. We show that the very low levels of endogenous CFTR expression, compared with other low expression genes, are maintained through the co-operative inhibitory effects of an upstream open reading frame and a thermodynamically stable RNA secondary structure.

IN-MACA-MCC: Integrated Multiple Attractor Cellular Automata with Modified Clonal Classifier for Human Protein Coding and Promoter Prediction

Protein coding and promoter region predictions are very important challenges of bioinformatics (Attwood and Teresa, 2000). The identification of these regions plays a crucial role in understanding the genes. Many novel computational and mathematical methods are introduced as well as existing methods that are getting refined for predicting both of the regions separately; still there is a scope for improvement. We propose a classifier that is built with MACA (multiple attractor cellular automata) and MCC (modified clonal classifier) to predict both regions with a single classifier. The proposed classifier is trained and tested with Fickett and Tung (1992) datasets for protein coding region prediction for DNA sequences of lengths 54, 108, and 162. This classifier is trained and tested with MMCRI datasets for protein coding region prediction for DNA sequences of lengths 252 and 354. The proposed classifier is trained and tested with promoter sequences from DBTSS (Yamashita et al., 2006) dataset and nonpromoters from EID (Saxonov et al., 2000) and UTRdb (Pesole et al., 2002) datasets. The proposed model can predict both regions with an average accuracy of 90.5% for promoter and 89.6% for protein coding region predictions. The specificity and sensitivity values of promoter and protein coding region predictions are 0.89 and 0.92, respectively.

Potential roles of 5´ UTR and 3´ UTR regions in post-trans-criptional regulation of mouse Oct4 gene in BMSC and P19 cells

OBJECTIVES

OCT4 is a transcription factor required for pluripotency during early embryogenesis and the maintenance of identity of embryonic stem cells and pluripotent cells. Therefore, the effective expression regulation of this gene is highly critical. UTR regions are of great significance to gene regulation. In this study, we aimed to investigate the potential regulatory role played by 5´UTR and 3´UTR of the Oct4 gene in mouse BMSC and P19 cells.

MATERIALS AND METHODS

The Oct4 5´UTR and 3´UTR sequences were cloned into pGL3 luciferase plasmid which led to the generation of pGL3 5´-UTR, pGL3 5´&3´-UTRs and pGL3 3´-UTR vectors. The vectors were transfected into BMSC and P19 cells followed by luciferase assay.

RESULTS

The assay of luciferase expression exhibited a direct link between the presence of Oct4 3´- UTR and the decrease of luciferase count in both cell lines; whereas 5´UTR indicated diverse behaviors in two cells. This discrepancy could be explained in view of the difference of cellular contexts in which the Oct4 UTRs act.

CONCLUSION

This study sheds some light on the role of UTR regions of mouse Oct4 in regulating post-transcriptional gene expression in pluripotent cells. These data represent potential to be used for the development of novel therapeutic approaches for a variety of malignancies.

Nonsense-mediated mRNA decay immunity can help identify human polycistronic transcripts

Eukaryotic polycistronic transcription units are rare and only a few examples are known, mostly being the outcome of serendipitous discovery. We claim that nonsense-mediated mRNA decay (NMD) immune structure is a common characteristic of polycistronic transcripts, and that this immunity is an emergent property derived from all functional CDSs. The human RefSeq transcriptome was computationally screened for transcripts capable of eliciting NMD, and which contain an additional ORF(s) potentially capable of rescuing the transcript from NMD. Transcripts were further analyzed implementing domain-based strategies in order to estimate the potential of the candidate ORF to encode a functional protein. Consequently, we predict the existence of forty nine novel polycistronic transcripts.

Experimental verification was carried out utilizing two different types of analyses. First, five Gene Expression Omnibus (GEO) datasets from published NMD-inhibition studies were used, aiming to explore whether a given mRNA is indeed insensitive to NMD. All known bicistronic transcripts and eleven out of the twelve predicted genes that were analyzed, displayed NMD insensitivity using various NMD inhibitors. For three genes, a mixed expression pattern was observed presenting both NMD sensitivity and insensitivity in different cell types. Second, we used published global translation initiation sequencing data from HEK293 cells to verify the existence of translation initiation sites in our predicted polycistronic genes. In five of our genes, the predicted rescuing uORFs are indeed identified as translation initiation sites, and in two additional genes, one of two predicted rescuing uORF is verified. These results validate our computational analysis and reinforce the possibility that NMD-immune architecture is a parameter by which polycistronic genes can be identified. Moreover, we present evidence for NMD-mediated regulation controlling the production of one or more proteins encoded in the polycistronic transcript.

Base pairing interaction between 5'- and 3'-UTRs controls icaR mRNA translation in Staphylococcus aureus

The presence of regulatory sequences in the 3′ untranslated region (3′-UTR) of eukaryotic mRNAs controlling RNA stability and translation efficiency is widely recognized. In contrast, the relevance of 3′-UTRs in bacterial mRNA functionality has been disregarded. Here, we report evidences showing that around one-third of the mapped mRNAs of the major human pathogen Staphylococcus aureus carry 3′-UTRs longer than 100-nt and thus, potential regulatory functions. We selected the long 3′-UTR of icaR, which codes for the repressor of the main exopolysaccharidic compound of the S. aureus biofilm matrix, to evaluate the role that 3′-UTRs may play in controlling mRNA expression. We showed that base pairing between the 3′-UTR and the Shine-Dalgarno (SD) region of icaR mRNA interferes with the translation initiation complex and generates a double-stranded substrate for RNase III. Deletion or substitution of the motif (UCCCCUG) within icaR 3′-UTR was sufficient to abolish this interaction and resulted in the accumulation of IcaR repressor and inhibition of biofilm development. Our findings provide a singular example of a new potential post-transcriptional regulatory mechanism to modulate bacterial gene expression through the interaction of a 3′-UTR with the 5′-UTR of the same mRNA.

At both sides of the protein-coding region, the mRNA molecule contains sequences that are not translated to protein. In eukaryotes, the untranslated 3′ region (3′-UTR), which comprises from the last codon used in translation to the 3′ end of the mRNA, controls mRNA stability, location and translation efficiency. In contrast, knowledge about the functions of 3′-UTRs in bacterial physiology is scarce. Here, we demonstrate that bacterial 3′-UTRs might play regulatory functions that might resemble those already described in eukaryotes. Transcriptome analysis of the human pathogen Staphylococcus aureus revealed that at least 30% of mRNAs contain long 3′-UTRs. Using the 3′-UTR of the mRNA encoding the main biofilm repressor IcaR as a model, we show that the 3′-UTR interferes with the translation initiation complex and promotes mRNA decay through base pairing with the ribosome binding site. This event contributes to adjusting IcaR level and modulating exopolysaccharide production and biofilm development in S. aureus. Our data illustrate that bacterial 3′-UTRs can provide strategies for fine-tuning control of gene expression.

Investigating the function of single nucleotide polymorphisms in the CTSB gene: a computational approach

Aim: Recent genome-wide association studies have revealed large numbers of single nucleotide polymorphisms (SNPs) related to Alzheimer’s disease. Here, we have investigated the gene CTSB, which plays a crucial role in encoding CTSB, a lysosomal cysteine proteinase protein. CTSB is also involved in the proteolytic processing of amyloid precursor protein (APP), which is believed to be a causative factor in Alzheimer’s disease. Materials & methods: Several bioinformatics algorithms such as, Sorting Intolerant from Tolerant (SIFT), Polymorphism Phenotyping (PolyPhen) and CUPSAT could identify the synonymous SNPs and nonsynonymous SNPs (nsSNPs), which are predicted to be deleterious and nondeleterious, respectively. Similar tools were used to predict the impact of single amino acid substitutions on CTSB protein activity. The FASTSNP server and UTRscan were used to predict the influence on splicing regulations. The stability and solvent-accessible surface area of modeled mutated proteins were analyzed using PBEQ solver and NetASA view. Furthermore, the DSP program was used to determine the secondary structures of the modeled protein. Results: A total of 999 SNPs in CTSB were retrieved from the SNP database; 55 nsSNPs, 35 synonymous SNPs, 165 mRNA were found in the 3´untranslated region SNPs, 12 SNPs were found in the 5´untranslated region in addition to 732 intronic SNPs. Potential functions of SNPs in the CTSB gene were identified using different web servers. For example, SIFT, PolyPhen and CUPSAT servers predicted ten nsSNPs to be intolerant, three nsSNPs to be damaging and eight nsSNPs to have the potential to destabilize protein structure. The FASTSNP server predicted 12 SNPs to influence splicing regulation, whereas two SNPs could predict a risk in the range of 3–4 (medium to high). Furthermore, mutant proteins were modeled and the total energy values were compared with the native CTSB protein. It was observed that on the surface of the protein, a mutation from threonine to serine at position 235 (rs17573) caused the greatest impact on stability. Conclusion: The genome-wide association studies database has already found rs7003814 of the CTSB gene reported against Alzheimer’s disease. Our study demonstrates the presence of other deleterious nsSNPs, which may play a crucial role in predicting Alzheimer’s disease risk.

Post-transcriptional regulation of connexin43 in H-Ras-transformed cells

Connexin43 (Cx43) expression is lost in cancer cells and many studies have reported that Cx43 is a tumor suppressor gene. Paradoxically, in a cellular NIH3T3 model, we have previously shown that Ha-Ras-mediated oncogenic transformation results in increased Cx43 expression. Although the examination of transcriptional regulation revealed essential regulatory elements, it could not solve this paradox. Here we studied post-transcriptional regulation of Cx43 expression in cancer using the same model in search of novel gene regulatory elements. Upon Ras transformation, both Cx43 mRNA stability and translation efficiency were increased. We investigated the role of Cx43 mRNA 3′ and 5′Untranslated regions (UTRs) and found an opposing effect; a 5′UTR-driven positive regulation is observed in Ras-transformed cells (NIH-3T3^Ras), while the 3′UTR is active only in normal NIH-3T3^Neo cells and completely silenced in NIH-3T3^Ras cells. Most importantly, we identified a previously unknown regulatory element within the 3′UTR, named S1516, which accounts for this 3′UTR-mediated regulation. We also examined the effect of other oncogenes and found that Ras- and Src-transformed cells show a different Cx43 UTRs post-transcriptional regulation than ErbB2-transformed cells, suggesting distinct regulatory pathways. Next, we detected different patterns of S1516 RNA-protein complexes in NIH-3T3^Neo compared to NIH-3T3^Ras cells. A proteomic approach identified most of the S1516-binding proteins as factors involved in post-transcriptional regulation. Building on our new findings, we propose a model to explain the discrepancy between the Cx43 expression in Ras-transformed NIH3T3 cells and the data in clinical specimens.

Evolutionary patterns of non-coding RNAs

A plethora of new functions of non-coding RNAs (ncRNAs) have been discovered in past few years. In fact, RNA is emerging as the central player in cellular regulation, taking on active roles in multiple regulatory layers from transcription, RNA maturation, and RNA modification to translational regulation. Nevertheless, very little is known about the evolution of this "Modern RNA World" and its components. In this contribution, we attempt to provide at least a cursory overview of the diversity of ncRNAs and functional RNA motifs in non-translated regions of regular messenger RNAs (mRNAs) with an emphasis on evolutionary questions. This survey is complemented by an in-depth analysis of examples from different classes of RNAs focusing mostly on their evolution in the vertebrate lineage. We present a survey of Y RNA genes in vertebrates and study the molecular evolution of the U7 snRNA, the snoRNAs E1/U17, E2, and E3, the Y RNA family, the let-7 microRNA (miRNA) family, and the mRNA-like evf-1 gene. We furthermore discuss the statistical distribution of miRNAs in metazoans, which suggests an explosive increase in the miRNA repertoire in vertebrates. The analysis of the transcription of ncRNAs suggests that small RNAs in general are genetically mobile in the sense that their association with a hostgene (e.g. when transcribed from introns of a mRNA) can change on evolutionary time scales. The let-7 family demonstrates, that even the mode of transcription (as intron or as exon) can change among paralogous ncRNA.

Computational identification of new structured cis-regulatory elements in the 3'-untranslated region of human protein coding genes

Messenger ribonucleic acids (RNAs) contain a large number of cis-regulatory RNA elements that function in many types of post-transcriptional regulation. These cis-regulatory elements are often characterized by conserved structures and/or sequences. Although some classes are well known, given the wide range of RNA-interacting proteins in eukaryotes, it is likely that many new classes of cis-regulatory elements are yet to be discovered. An approach to this is to use computational methods that have the advantage of analysing genomic data, particularly comparative data on a large scale. In this study, a set of structural discovery algorithms was applied followed by support vector machine (SVM) classification. We trained a new classification model (CisRNA-SVM) on a set of known structured cis-regulatory elements from 3'-untranslated regions (UTRs) and successfully distinguished these and groups of cis-regulatory elements not been strained on from control genomic and shuffled sequences. The new method outperformed previous methods in classification of cis-regulatory RNA elements. This model was then used to predict new elements from cross-species conserved regions of human 3'-UTRs. Clustering of these elements identified new classes of potential cis-regulatory elements. The model, training and testing sets and novel human predictions are available at: http://mRNA.otago.ac.nz/CisRNA-SVM.

Identification of 3'UTR sequence elements and a teloplasm localization motif sufficient for the localization of Hro-twist mRNA to the zygotic animal and vegetal poles

The early localization of mRNA transcripts is critical in sorting cell fate determinants in the developing embryo. In the glossiphoniid leech, Helobdella robusta, maternal mRNAs, such as Hro-twist, localize to the zygotic teloplasm. Ten seven nucleotide repeat elements (AAUAAUA) called ARE2 and a predicted secondary structural motif, called teloplasm localization motif (TLM), are present in the 3'UTR of Hro-twist mRNA. We used site-directed mutagenesis, deletions, and microinjection of labeled, exogenous transcripts to determine if ARE2 elements, and the TLM, play a role in Hro-twist mRNA localization. Deleting the poly-A tail and the cytoplasmic polyadenylation element (CPE) had no effect on Hro-twist mRNA localization. Site-directed mutagenesis of nucleotides that altered ARE2 element sequences or the TLM suggest that the ARE2 elements and the TLM are important for Hro-twist mRNA localization to the teloplasm of pre-cleavage zygotes. Hro-Twist protein expression data suggest that the localization of Hro-twist transcripts in zygotes and stage two embryos is not involved in ensuring mesoderm specification, as Hro-Twist protein is expressed uniformly in most cells before gastrulation. Our data may support a shared molecular mechanism for leech transcripts that localize to the teloplasm.

Yellow lupin (Lupinus luteus L.) transcriptome sequencing: molecular marker development and comparative studies

Background

Yellow lupin (Lupinus luteus L.) is a minor legume crop characterized by its high seed protein content. Although grown in several temperate countries, its orphan condition has limited the generation of genomic tools to aid breeding efforts to improve yield and nutritional quality. In this study, we report the construction of 454-expresed sequence tag (EST) libraries, carried out comparative studies between L. luteus and model legume species, developed a comprehensive set of EST-simple sequence repeat (SSR) markers, and validated their utility on diversity studies and transferability to related species.

Results

Two runs of 454 pyrosequencing yielded 205 Mb and 530 Mb of sequence data for L1 (young leaves, buds and flowers) and L2 (immature seeds) EST- libraries. A combined assembly (L1L2) yielded 71,655 contigs with an average contig length of 632 nucleotides. L1L2 contigs were clustered into 55,309 isotigs. 38,200 isotigs translated into proteins and 8,741 of them were full length. Around 57% of L. luteus sequences had significant similarity with at least one sequence of Medicago, Lotus, Arabidopsis, or Glycine, and 40.17% showed positive matches with all of these species. L. luteus isotigs were also screened for the presence of SSR sequences. A total of 2,572 isotigs contained at least one EST-SSR, with a frequency of one SSR per 17.75 kbp. Empirical evaluation of the EST-SSR candidate markers resulted in 222 polymorphic EST-SSRs. Two hundred and fifty four (65.7%) and 113 (30%) SSR primer pairs were able to amplify fragments from L. hispanicus and L. mutabilis DNA, respectively. Fifty polymorphic EST-SSRs were used to genotype a sample of 64 L. luteus accessions. Neighbor-joining distance analysis detected the existence of several clusters among L. luteus accessions, strongly suggesting the existence of population subdivisions. However, no clear clustering patterns followed the accession’s origin.

Conclusion

L. luteus deep transcriptome sequencing will facilitate the further development of genomic tools and lupin germplasm. Massive sequencing of cDNA libraries will continue to produce raw materials for gene discovery, identification of polymorphisms (SNPs, EST-SSRs, INDELs, etc.) for marker development, anchoring sequences for genome comparisons and putative gene candidates for QTL detection.

Application of serial analysis of gene expression to the study of the gene expression profile of Leishmania infantum chagasi promastigote

This study describes the application of the LongSAGE methodology to study the gene expression profile in promastigotes of Leishmania infantum chagasi. A tag library was created using the LongSAGE method and consisted of 14,208 tags of 17 bases. Of these, 8,427 (59.3%) were distinct. BLAST research of the 1,645 most abundant tags showed that 12.8% of them identified the coding sequences of genes, while 82% (1,349/1,645) identified one or more genomic sequences that did not correspond with open reading frames. Only 5.2% (84/1,645) of the tags were not aligned to any position in the L. infantum genome. The UTR size of Leishmania and the lack of CATG sites in some transcripts were decisive for the generation of tags in these regions. Additional analysis will allow a better understanding of the expression profile and discovering the key genes in this life cycle.

Genetic variation in RTN4 3'-UTR and susceptibility to cervical squamous cell carcinoma

Recent studies have suggested that RTN4 is a multifunctional gene, including inhibition of axonal regeneration, vascular remodeling, apoptosis, and tumor suppression. The TATC and CAA insertion/deletion polymorphisms of RTN4 3'-UTR have been linked to schizophrenia, depression, and dilated cardiomyopathy. To test whether these two polymorphisms are associated with cervical squamous cell carcinoma (CSCC), in this research, by using polymerase chain reaction-polyacrylamide gel electrophoresis, we determined the genotypes of the TATC and CAA polymorphisms in 336 CSCC patients and 450 unrelated control subjects. Allele frequencies of TATC and CAA polymorphisms were not significantly different between CSCC patients and control subjects (odds ratio [OR]=1.22, 95% confidence interval [CI]=0.98-1.50 for TATC; OR=0.95, 95% CI=0.76-1.18 for CAA). Decreased CSCC risk was associated with TATC polymorphism in a recessive model (OR=0.49, 95% CI=0.30-0.77), while no significant association was observed between CAA polymorphism and CSCC in different genetic models. Results of stratified analysis revealed that both TATC and CAA polymorphisms were associated with high clinical stage, and CAA polymorphism was also associated with positive parametrial invasion (OR=0.69, 95% CI=0.48-0.98). The present study provides evidence that TATC and CAA insertion/deletion polymorphisms are associated with CSCC, indicating that genetic variation in RTN4 3'-UTR contributes to the susceptibility to CSCC. It is necessary to confirm these findings in ethnically different populations and with a larger sample.

RNAi-mediated targeting of noncoding and coding sequences in DNA repair gene messages efficiently radiosensitizes human tumor cells

Human tumor cell death during radiotherapy is caused mainly by ionizing radiation (IR)-induced DNA double-strand breaks (DSB), which are repaired by either homologous recombination repair (HRR) or nonhomologous end-joining (NHEJ). Although siRNA-mediated knockdown of DNA DSB repair genes can sensitize tumor cells to IR, this approach is limited by inefficiencies of gene silencing. In this study, we show that combining an artificial miRNA (amiR) engineered to target 3'-untranslated regions of XRCC2 (an HRR factor) or XRCC4 (an NHEJ factor) along with an siRNA to target the gene coding region can improve silencing efficiencies to achieve more robust radiosensitization than a single approach alone. Mechanistically, the combinatorial knockdown decreased targeted gene expression through both a reduction in mRNA stability and a blockade to mRNA translation. Together, our findings establish a general method of gene silencing that is more efficient and particularly suited for suppressing genes that are difficult to downregulate by amiR- or siRNA-based methods alone.

Control of messenger RNA fate by RNA-binding proteins: an emphasis on mammalian spermatogenesis

Posttranscriptional status of messenger RNAs (mRNA) can be affected by many factors, most of which are RNA-binding proteins (RBP) that either bind mRNA in a nonspecific manner or through specific motifs, usually located in the 3' untranslated regions. RBPs can also be recruited by small noncoding RNAs (sncRNA), which have been shown to be involved in posttranscriptional regulations and transposon repression (eg, microRNAs or P-element-induced wimpy testis-interacting RNA) as components of the sncRNA effector complex. Non-sncRNA-binding RBPs have much more diverse effects on their target mRNAs. Some can cause degradation of their target transcripts and/or repression of translation, whereas others can stabilize and/or activate translation. The splicing and exportation of transcripts from the nucleus to the cytoplasm are often mediated by sequence-specific RBPs. The mechanisms by which RBPs regulate mRNA transcripts involve manipulating the 3' poly(A) tail, targeting the transcript to polysomes or to other ribonuclear protein particles, recruiting regulatory proteins, or competing with other RBPs. Here, we briefly review the known mechanisms of posttranscriptional regulation mediated by RBPs, with an emphasis on how these mechanisms might control spermatogenesis in general.

C-T variant in a miRNA target site of BCL2 is associated with increased risk of human papilloma virus related cervical cancer--an in silico approach

MicroRNAs control gene expression at the posttranscriptional level by base-pairing to the 3'-UTR of their target mRNAs, thus leading to mRNA degradation of protein fabrication. We hypothesize, SNPs within miRNAs and their targets could be of significance to an individual's risk of developing cancer. We analyzed in silico SNP information on cervical cancer associated aberrant alleles and further investigated this in a case-control study by examining eleven SNPs from different genes. It was observed that a C to T polymorphism in putative miRNA target site of BCL2 was significantly conspicuous for the aberrant SNP allele in cancer tissues as compared to controls. This study provides evidence that SNPs in miRNA-binding sites may play an important role in increasing risk of cancer. The results also paves way for future studies to validate these results in other well-characterized populations as well as to explore the biological significance of these particular SNPs.

The 3'-untranslated region length and AU-rich RNA location modulate RNA-protein interaction and translational control of β2-adrenergic receptor mRNA

Posttranscriptional controls play a major role in β(2)-adrenergic receptor (β(2)-AR) expression. We recently reported that β(2)-AR mRNA translation is suppressed by elements in its 3'-untranslated region (UTR). We also identified T-cell-restricted intracellular antigen-related protein (TIAR) and HuR as prominent AU-rich (ARE) RNA-binding proteins that associate with β(2)-AR mRNA 3'-UTR. In this study, we identified a poly(U) region at the distal end of the 3'-UTR as critical for TIAR binding to β(2)-AR mRNA and for translational suppression. Here, we also report that the locations of the poly(U) and ARE sequences within the 3'-UTR are important determinants that control the translation of β(2)-AR mRNA. Consistent with this finding, a 20-nucleotide ARE RNA from the proximal 3'-UTR that did not inhibit mRNA translation in its native position was able to suppress translation when re-located to the distal 3'-UTR of the receptor mRNA. Immunoprecipitation and polysome profile analysis demonstrated the importance of 3'-UTR length and the ARE RNA location within the 3'-UTR, as key determinants of RNA/protein interactions and translational control of β(2)-AR mRNA. Further, the importance of 3'-UTR length and ARE location in TIAR and HuR association with mRNA and translational suppression was demonstrated using a chimeric luciferase reporter gene.

Association of the human CYP11B2 gene and essential hypertension in southwest Han Chinese population: a haplotype-based case-control study

Aldosterone synthase produces aldosterone, which regulates electrolytes and thereby blood pressure (BP). The aldosterone-synthase gene (CYP11B2) has been regarded as a candidate gene for essential hypertension. To address this issue, we carried out a haplotype-based, case-control study to explore the association between a human CYP11B2 gene and essential hypertension (EH) in the southwest Han population of China (n = 1020 individuals). Four tag single-nucleotide polymorphisms (SNPs) (rs4536, rs4545, rs3097, and rs3802230) and the C-344T polymorphism, as well as the K173R polymorphism in the CYP11B2 gene, were genotyped using the PCR-RFLP method. Single-locus analysis showed that the C allele of rs3802230 was significantly more prevalent in the EH subjects as compared to control subjects, adjusted for covariates. Haplotype analysis showed that the haplotype AAGC constructed by the tag SNPs (rs4536, rs4545, rs3097, and rs3802230), which carried the susceptible rs3802230 C allele, significantly increased the risk of essential hypertension with an odds ratios equal to 3.56 (P = 0.0001). The present results indicated that the rs3802230 C allele might be a risk marker for essential hypertension and haplotype AAGC might confer high genetic susceptibility to essential hypertension in a southwest Han Chinese population.

MiRenSVM: towards better prediction of microRNA precursors using an ensemble SVM classifier with multi-loop features

Background

MicroRNAs (simply miRNAs) are derived from larger hairpin RNA precursors and play essential regular roles in both animals and plants. A number of computational methods for miRNA genes finding have been proposed in the past decade, yet the problem is far from being tackled, especially when considering the imbalance issue of known miRNAs and unidentified miRNAs, and the pre-miRNAs with multi-loops or higher minimum free energy (MFE). This paper presents a new computational approach, miRenSVM, for finding miRNA genes. Aiming at better prediction performance, an ensemble support vector machine (SVM) classifier is established to deal with the imbalance issue, and multi-loop features are included for identifying those pre-miRNAs with multi-loops.

Results

We collected a representative dataset, which contains 697 real miRNA precursors identified by experimental procedure and other computational methods, and 5428 pseudo ones from several datasets. Experiments showed that our miRenSVM achieved a 96.5% specificity and a 93.05% sensitivity on the dataset. Compared with the state-of-the-art approaches, miRenSVM obtained better prediction results. We also applied our method to predict 14 Homo sapiens pre-miRNAs and 13 Anopheles gambiae pre-miRNAs that first appeared in miRBase13.0, MiRenSVM got a 100% prediction rate. Furthermore, performance evaluation was conducted over 27 additional species in miRBase13.0, and 92.84% (4863/5238) animal pre-miRNAs were correctly identified by miRenSVM.

Conclusion

MiRenSVM is an ensemble support vector machine (SVM) classification system for better detecting miRNA genes, especially those with multi-loop secondary structure.

Computational identification and structural analysis of deleterious functional SNPs in MLL gene causing acute leukemia

A promising application of the huge amounts of data from the Human Genome Project currently available offers new opportunities for identifying the genetic predisposition and developing a better understanding of complex diseases such as cancers. The main focus of cancer genetics is the study of mutations that are causally implicated in tumorigenesis. The identification of such causal mutations does not only provide insight into cancer biology but also presents anticancer therapeutic targets and diagnostic markers. In this study, we evaluated the Single Nucleotide Polymorphisms (SNPs) that can alter the expression and the function in MLL gene through computational methods. We applied an evolutionary perspective to screen the SNPs using a sequence homologybased SIFT tool, suggested that 10 non-synonymous SNPs (nsSNPs) (50%) were found to be deleterious. Structure based approach PolyPhen server suggested that 5 nsSNPS (25%) may disrupt protein function and structure. PupaSuite tool predicted the phenotypic effect of SNPs on the structure and function of the affected protein. Structure analysis was carried out with the major mutations that occurred in the native protein coded by MLL gene is at amino acid positions Q1198P and K1203Q. The solvent accessibility results showed that 7 residues changed from exposed state in the native type protein to buried state in Q1198P mutant protein and remained unchanged in the case of K1203Q. From the overall results obtained, nsSNP with id (rs1784246) at the amino acid position Q1198P could be considered as deleterious mutation in the acute leukemia caused by MLL gene.

SCS: signal, context, and structure features for genome-wide human promoter recognition

This paper integrates the signal, context, and structure features for genome-wide human promoter recognition, which is important in improving genome annotation and analyzing transcriptional regulation without experimental supports of ESTs, cDNAs, or mRNAs. First, CpG islands are salient biological signals associated with approximately 50 percent of mammalian promoters. Second, the genomic context of promoters may have biological significance, which is based on n-mers (sequences of n bases long) and their statistics estimated from training samples. Third, sequence-dependent DNA flexibility originates from DNA 3D structures and plays an important role in guiding transcription factors to the target site in promoters. Employing decision trees, we combine above signal, context, and structure features to build a hierarchical promoter recognition system called SCS. Experimental results on controlled data sets and the entire human genome demonstrate that SCS is significantly superior in terms of sensitivity and specificity as compared to other state-of-the-art methods. The SCS promoter recognition system is available online as supplemental materials for academic use and can be found on the Computer Society Digital Library at http://doi.ieeecomputersociety.org/10.1109/TCBB.2008.95.

Analysis of bovine foamy virus btas mRNA transcripts during persistent infection

Foamy virus (FV) is an unconventional retrovirus that possesses a complex genome and a special mechanism for gene expression regulation. The genome encodes transcriptional protein Tas which is found to regulate both the internal promoter (IP) and the long terminal repeat promoter (LTR). However, the detailed mechanism of Tas-mediated gene expression remains unknown. In this study, we provided the first evidence for the temporal production and utilization of four different bovine foamy virus (BFV) btas mRNAs during persistent infection. These four forms of btas mRNA transcripts initiated either at BFV LTR or IP and spliced or unspliced have a differential ability to activate BFV promoters. Furthermore, by developing an MS2 translational operator/coat protein combined system to track mRNA exportation from the nucleus and distribution throughout the cytoplasm, we observed that the IP spliced transcript could be exported into the cytoplasm more efficiently than unspliced transcripts. These findings provide evidence for the hypothesis that the functional interplay of both promoters contributes to the temporal pattern of BFV transcription and suggest that a post-transcriptional regulation exist in BFV replication.

PAP1 [poly(A) polymerase 1] homozygosity and hyperadenylation are major determinants of increased mRNA stability of CDR1 in azole-resistant clinical isolates of Candida albicans

Using genetically matched azole-susceptible (AS) and azole-resistant (AR) clinical isolates of Candida albicans, we recently demonstrated that CDR1 overexpression in AR isolates is due to its enhanced transcriptional activation and mRNA stability. This study examines the molecular mechanisms underlying enhanced CDR1 mRNA stability in AR isolates. Mapping of the 3' untranslated region (3' UTR) of CDR1 revealed that it was rich in adenylate/uridylate (AU) elements, possessed heterogeneous polyadenylation sites, and had putative consensus sequences for RNA-binding proteins. Swapping of heterologous and chimeric lacZ-CDR1 3' UTR transcriptional reporter fusion constructs did not alter the reporter activity in AS and AR isolates, indicating that cis-acting sequences within the CDR1 3' UTR itself are not sufficient to confer the observed differential mRNA decay. Interestingly, the poly(A) tail of the CDR1 mRNA of AR isolates was approximately 35-50 % hyperadenylated as compared with AS isolates. C. albicans poly(A) polymerase (PAP1), responsible for mRNA adenylation, resides on chromosome 5 in close proximity to the mating type-like (MTL) locus. Two different PAP1 alleles, PAP1-a/PAP1-alpha, were recovered from AS (MTL-a/MTL-alpha), while a single type of PAP1 allele (PAP1-alpha) was recovered from AR isolates (MTL-alpha/MTL-alpha). Among the heterozygous deletions of PAP1-a (Deltapap1-a/PAP1-alpha) and PAP1-alpha (PAP1-a/Deltapap1-alpha), only the former led to relatively enhanced drug resistance, to polyadenylation and to transcript stability of CDR1 in the AS isolate. This suggests a dominant negative role of PAP1-a in CDR1 transcript polyadenylation and stability. Taken together, our study provides the first evidence, to our knowledge, that loss of heterozygosity at the PAP1 locus is linked to hyperadenylation and subsequent increased stability of CDR1 transcripts, thus contributing to enhanced drug resistance.

UTRdb and UTRsite (RELEASE 2010): a collection of sequences and regulatory motifs of the untranslated regions of eukaryotic mRNAs

The 5' and 3' untranslated regions of eukaryotic mRNAs (UTRs) play crucial roles in the post-transcriptional regulation of gene expression through the modulation of nucleo-cytoplasmic mRNA transport, translation efficiency, subcellular localization and message stability. UTRdb is a curated database of 5' and 3' untranslated sequences of eukaryotic mRNAs, derived from several sources of primary data. Experimentally validated functional motifs are annotated and also collated as the UTRsite database where more specific information on the functional motifs and cross-links to interacting regulatory protein are provided. In the current update, the UTR entries have been organized in a gene-centric structure to better visualize and retrieve 5' and 3'UTR variants generated by alternative initiation and termination of transcription and alternative splicing. Experimentally validated miRNA targets and conserved sequence elements are also annotated. The integration of UTRdb with genomic data has allowed the implementation of an efficient annotation system and a powerful retrieval resource for the selection and extraction of specific UTR subsets. All internet resources implemented for retrieval and functional analysis of 5' and 3' untranslated regions of eukaryotic mRNAs are accessible at http://utrdb.ba.itb.cnr.it/.

Phytochrome B and histone deacetylase 6 control light-induced chromatin compaction in Arabidopsis thaliana

Natural genetic variation in Arabidopsis thaliana exists for many traits and often reflects acclimation to local environments. Studying natural variation has proven valuable in the characterization of phenotypic traits and, in particular, in identifying genetic factors controlling these traits. It has been previously shown that chromatin compaction changes during development and biotic stress. To gain more insight into the genetic control of chromatin compaction, we investigated the nuclear phenotype of 21 selected Arabidopsis accessions from different geographic origins and habitats. We show natural variation in chromatin compaction and demonstrate a positive correlation with latitude of geographic origin. The level of compaction appeared to be dependent on light intensity. A novel approach, combining Quantitative Trait Locus (QTL) mapping and microscopic examination, pointed at PHYTOCHROME-B (PHYB) and HISTONE DEACETYLASE-6 (HDA6) as positive regulators of light-controlled chromatin compaction. Indeed, mutant analyses demonstrate that both factors affect global chromatin organization. HDA6, in addition, strongly promotes the light-mediated compaction of the Nucleolar Organizing Regions (NORs). The accession Cape Verde Islands-0 (Cvi-0), which shows sequence polymorphism in the PHYB gene and in the HDA6 promotor, resembles the hda6 mutant in having reduced chromatin compaction and decreased methylation levels of DNA and histone H3K9 at the NORs. We provide evidence that chromatin organization is controlled by light intensity. We propose that chromatin plasticity is associated with acclimation of Arabidopsis to its environment. The polymorphic alleles such as PHYB and HDA6 control this process.

The habitat of the plant model species Arabidopsis thaliana can be found throughout the Northern hemisphere. As a consequence, individual populations have acclimated to a great diversity of environmental conditions. This is reflected by a wealth of natural genetic variation in many phenotypic traits. We utilized this natural variation via a novel approach, combining microscopic examination, quantitative genetics, and analysis of environmental parameters, to understand the regulation of nuclear chromatin compaction in leaf mesophyll cells. We show that the level of chromatin compaction among natural Arabidopsis thaliana accessions correlates with latitude of origin and depends on local light intensity. Our study provides evidence that the photoreceptor PHYTOCHROME-B (PHYB) and the histone modifier HISTONE DEACETYLASE 6 (HDA6) are positive regulators of global chromatin organization in a light-dependent manner. In addition, HDA6 specifically controls light-mediated chromatin compaction of the Nucleolar Organizing Regions (NORs). We propose that the observed light-controlled plasticity of chromatin plays a role in acclimation and survival of plants in their natural environment.

Regulatory element identification in subsets of transcripts: comparison and integration of current computational methods

Regulatory elements in mRNA play an often pivotal role in post-transcriptional regulation of gene expression. However, a systematic approach to efficiently identify putative regulatory elements from sets of post-transcriptionally coregulated genes is lacking, hampering studies of coregulation mechanisms. Although there are several analytical methods that can be used to detect conserved mRNA regulatory elements in a set of transcripts, there has been no systematic study of how well any of these methods perform individually or as a group. We therefore compared how well three algorithms, each based on a different principle (enumeration, optimization, or structure/sequence profiles), can identify elements in unaligned untranslated sequence regions. Two algorithms were originally designed to detect transcription factor binding sites, Weeder and BioProspector; and one was designed to detect RNA elements conserved in structure, RNAProfile. Three types of elements were examined: (1) elements conserved in both primary sequence and secondary structure; (2) elements conserved only in primary sequence; and (3) microRNA targets. Our results indicate that all methods can uniquely identify certain known RNA elements, and therefore, integrating the output from all algorithms leads to the most complete identification of elements. We therefore developed an approach to integrate results and guide selection of candidate elements from several algorithms presented as a web service (https://dbw.msi.umn.edu:8443/recit). These findings together with the approach for integration can be used to identify candidate elements from genome-wide post-transcriptional profiling data sets.