Dataset Information

NASP: a parallel program for identifying evolutionarily conserved nucleic acid secondary structures from nucleotide sequence alignments.

ABSTRACT:

Summary

Many natural nucleic acid sequences have evolutionarily conserved secondary structures with diverse biological functions. A reliable computational tool for identifying such structures would be very useful in guiding experimental analyses of their biological functions. NASP (Nucleic Acid Structure Predictor) is a program that takes into account thermodynamic stability, Boltzmann base pair probabilities, alignment uncertainty, covarying sites and evolutionary conservation to identify biologically relevant secondary structures within multiple sequence alignments. Unique to NASP is the consideration of all this information together with a recursive permutation-based approach to progressively identify and list the most conserved probable secondary structures that are likely to have the greatest biological relevance. By focusing on identifying only evolutionarily conserved structures, NASP forgoes the prediction of complete nucleotide folds but outperforms various other secondary structure prediction methods in its ability to selectively identify actual base pairings.

Availability

Downloable and web-based versions of NASP are freely available at http://web.cbio.uct.ac.za/~yves/nasp_portal.php

Contact

yves@cbio.uct.ac.za

Supplementary information

Supplementary data are available at Bioinformatics online.

SUBMITTER: Semegni JY

PROVIDER: S-EPMC3157929 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Significance statementMutations in hepatocyte nuclear factor-1 β ( HNF1B ) are the most common monogenic causes of congenital renal malformations. HNF1B is necessary to directly reprogram fibroblasts to induced renal tubule epithelial cells (iRECs) and, as we demonstrate, can induce ectopic pronephric tissue in Xenopus ectodermal organoids. Using these two systems, we analyzed the effect of HNF1B mutations found in patients with cystic dysplastic kidney disease. We found cross-species conserved targets of HNF1B, identified transcripts that are differentially regulated by the patient-specific mutant protein, and functionally validated novel HNF1B targets in vivo . These results highlight evolutionarily conserved transcriptional mechanisms and provide insights into the genetic circuitry of nephrogenesis.BackgroundHepatocyte nuclear factor-1 β (HNF1B) is an essential transcription factor during embryogenesis. Mutations in HNF1B are the most common monogenic causes of congenital cystic dysplastic renal malformations. The direct functional consequences of mutations in HNF1B on its transcriptional activity are unknown.MethodsDirect reprogramming of mouse fibroblasts to induced renal tubular epithelial cells was conducted both with wild-type HNF1B and with patient mutations. HNF1B was expressed in Xenopus ectodermal explants. Transcriptomic analysis by bulk RNA-Seq identified conserved targets with differentially regulated expression by the wild-type or R295C mutant. CRISPR/Cas9 genome editing in Xenopus embryos evaluated transcriptional targets in vivo .ResultsHNF1B is essential for reprogramming mouse fibroblasts to induced renal tubular epithelial cells and induces development of ectopic renal organoids from pluripotent Xenopus cells. The mutation R295C retains reprogramming and inductive capacity but alters the expression of specific sets of downstream target genes instead of diminishing overall transcriptional activity of HNF1B. Surprisingly, targets associated with polycystic kidney disease were less affected than genes affected in congenital renal anomalies. Cross-species-conserved transcriptional targets were dysregulated in hnf1b CRISPR-depleted Xenopus embryos, confirming their dependence on hnf1b .ConclusionsHNF1B activates an evolutionarily conserved program of target genes that disease-causing mutations selectively disrupt. These findings provide insights into the renal transcriptional network that controls nephrogenesis.

Project description:The molecular basis for function of the mammalian H19 as a tumor suppressor is poorly understood. Large, conserved open reading frames (ORFs) are absent from both the human and mouse cDNAs, suggesting that it may act as an RNA. Contradicting earlier reports, however, recent studies have shown that the H19 transcript exists in polysomal form and is likely translated. To distinguish between possible functional roles for the gene product, we have characterized the sequence requirements for H19-mediated in vitro suppression of tumor cell clonogenicity and analyzed the sequence of the gene cloned from a range of mammals. A cDNA version of the human gene, lacking the unusually short introns characteristic of imprinted genes, is as effective as a genomic copy in blocking anchorage-independent growth by G401 cells. The first 710 nucleotides of the gene can be deleted with no effect on in vitro activity. Further truncations from either the 5'- or 3'-end, however, cause a loss of suppression of clonogenicity. Using conserved sequences within the H19 gene as PCR primers, genomic DNA fragments were amplified from a range of mammalian species that span the functional domain defined by deletion analysis. Sequences from cat, lynx, elephant, gopher and orangutan complement the previous database of sequences from human, mouse, rat and rabbit. Hypothetical translation of the resulting sequences shows an absence of conserved ORFs of any size. Free energy and covariational analysis of the RNA sequences was used to identify potential helical pairings within the H19 transcript. A set of 16 helices are supported by covariation (i.e. conservation of base pairing potential in the absence of primary sequence conservation). The predicted RNA pairings consist largely of local hairpins but also include several long range interactions that bridge the 5'- and 3'-ends of the functional domain. Given the evolutionary conservation of structure at the RNA level and the absence of conservation at the protein level, we presume that the functional product of the H19 gene is a structured RNA.

Project description:BACKGROUND:NASP is an essential protein in mammals that functions in histone transport pathways and maintenance of a soluble reservoir of histones H3/H4. NASP has been studied exclusively in Opisthokonta lineages where some functional diversity has been reported. In humans, growing evidence implicates NASP miss-regulation in the development of a variety of cancers. Although a comprehensive phylogenetic analysis is lacking, NASP-family proteins that possess four TPR motifs are thought to be widely distributed across eukaryotes. RESULTS:We characterize the molecular evolution of NASP by systematically identifying putative NASP orthologs across diverse eukaryotic lineages ranging from excavata to those of the crown group. We detect extensive silent divergence at the nucleotide level suggesting the presence of strong purifying selection acting at the protein level. We also observe a selection bias for high frequencies of acidic residues which we hypothesize is a consequence of their critical function(s), further indicating the role of functional constraints operating on NASP evolution. Our data indicate that TPR1 and TPR4 constitute the most rapidly evolving functional units of NASP and may account for the functional diversity observed among well characterized family members. We also show that NASP paralogs in ray-finned fish have different genomic environments with clear differences in their GC content and have undergone significant changes at the protein level suggesting functional diversification. CONCLUSION:We draw four main conclusions from this study. First, wide distribution of NASP throughout eukaryotes suggests that it was likely present in the last eukaryotic common ancestor (LECA) possibly as an important innovation in the transport of H3/H4. Second, strong purifying selection operating at the protein level has influenced the nucleotide composition of NASP genes. Further, we show that selection has acted to maintain a high frequency of functionally relevant acidic amino acids in the region that interrupts TPR2. Third, functional diversity reported among several well characterized NASP family members can be explained in terms of quickly evolving TPR1 and TPR4 motifs. Fourth, NASP fish specific paralogs have significantly diverged at the protein level with NASP2 acquiring a NNR domain.