Search

Purpose: Methods for depth normalization have been assessed primarily with simulated data or cell-line-mixture data. There is a pressing need for benchmark data enabling a more realistic and objective assessment, especially in the context of small RNA sequencing., Methods: We collected a unique pair of microRNA sequencing data sets for the same set of tumor samples; one data set was collected with and the other without uniform handling and balanced design. The former provided a benchmark for evaluating evidence of differential expression and the latter served as a test bed for normalization. Next, we developed a data perturbation algorithm to simulate additional data set pairs. Last, we assembled a set of computational tools to visualize and quantify the assessment., Results: We validated the quality of the benchmark data and showed the need for normalization of the test data. For illustration, we applied the data and tools to assess the performance of 9 existing normalization methods. Among them, trimmed mean of M-values was a better scaling method, whereas the median and the upper quartiles were consistently the worst performers; one variation of remove unwanted variation had the best chance of capturing true positives but at the cost of increased false positives. In general, these methods were, at best, moderately helpful when the level of differential expression was extensive and asymmetric., Conclusion: Our study (1) provides the much-needed benchmark data and computational tools for assessing depth normalization, (2) shows the dependence of normalization performance on the underlying pattern of differential expression, and (3) calls for continued research efforts to develop more effective normalization methods.

Control of mRNA levels, a fundamental aspect in the regulation of gene expression, is achieved through a balance between mRNA synthesis and decay. E26-related gene (Erg) proteins are canonical transcription factors whose previously described functions are confined to the control of mRNA synthesis. Here, we report that ERG also regulates gene expression by affecting mRNA stability and identify the molecular mechanisms underlying this function in human cells. ERG is recruited to mRNAs via interaction with the RNA-binding protein RBPMS, and it promotes mRNA decay by binding CNOT2, a component of the CCR4-NOT deadenylation complex. Transcriptome-wide mRNA stability analysis revealed that ERG controls the degradation of a subset of mRNAs highly connected to Aurora signaling, whose decay during S phase is necessary for mitotic progression. Our data indicate that control of gene expression by mammalian transcription factors may follow a more complex scheme than previously anticipated, integrating mRNA synthesis and degradation.

RNA-binding protein with multiple splicing (designated RBPMS) is a higher vertebrate mRNA-binding protein containing a single RNA recognition motif (RRM). RBPMS has been shown to be involved in mRNA transport, localization and stability, with key roles in axon guidance, smooth muscle plasticity, as well as regulation of cancer cell proliferation and migration. We report on structure-function studies of the RRM domain of RBPMS bound to a CAC-containing single-stranded RNA. These results provide insights into potential topologies of complexes formed by the RBPMS RRM domain and the tandem CAC repeat binding sites as detected by photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation. These studies establish that the RRM domain of RBPMS forms a symmetrical dimer in the free state, with each monomer binding sequence-specifically to all three nucleotides of a CAC segment in the RNA bound state. Structure-guided mutations within the dimerization and RNA-binding interfaces of RBPMS RRM on RNA complex formation resulted in both disruption of dimerization and a decrease in RNA-binding affinity as observed by size exclusion chromatography and isothermal titration calorimetry. As anticipated from biochemical binding studies, over-expression of dimerization or RNA-binding mutants of Flag-HA-tagged RBPMS were no longer able to track with stress granules in HEK293 cells, thereby documenting the deleterious effects of such mutations in vivo.

Central nervous system (CNS) is an increasingly common site of isolated metastasis for patients with Stage 4 neuroblastoma. To explore the microRNA (miRNA) profile of this metastatic process, miRNA sequencing was performed to identify miRNA sequence families with differential expression between tumor pairs (pre-CNS primary and CNS metastasis) from 13 patients with Stage 4 neuroblastoma. Seven miRNA sequence families had distinct expression in CNS metastases when compared with their corresponding pre-CNS primaries. MiR-7 was upregulated (3.75-fold), and miR-21, miR-22, miR-29a, miR-143, miR-199a-1-3p, and miR-199a-1-5p were downregulated (3.5-6.1-fold), all confirmed by quantitative reverse transcription-PCR. MiR-29a, previously shown to be downregulated in a broad spectrum of solid tumors including neuroblastoma, had the most significant decrease in all 13 CNS metastases (P = 0.001). Its known onco-targets CDC6, CDK6, and DNMT3A, as well as B7-H3, an inhibitory ligand for T cells, and natural killer cells, were found to have higher differential expression in these 13 CNS metastases when compared with their paired primaries. Additionally, miR-29a expression in primary tumors was significantly lower among patients who eventually relapsed in the CNS. Irrespective of the amplification status of MYCN, which is known to be associated with metastasis, pre-CNS primaries, and CNS metastases had significantly lower miR-29a expression than non-CNS primary tumors. Among MYCN amplified cell lines, those from CNS relapse also had lower miR-29a expression than non-CNS relapse. These findings raised the hypothesis that miR-29a could be a biomarker for neuroblastoma CNS metastasis, and its downregulation may play a pivotal role in CNS progression., (© 2014 Wiley Periodicals, Inc.)

Recent studies implicated the RNA-binding protein with multiple splicing (RBPMS) family of proteins in oocyte, retinal ganglion cell, heart, and gastrointestinal smooth muscle development. These RNA-binding proteins contain a single RNA recognition motif (RRM), and their targets and molecular function have not yet been identified. We defined transcriptome-wide RNA targets using photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) in HEK293 cells, revealing exonic mature and intronic pre-mRNA binding sites, in agreement with the nuclear and cytoplasmic localization of the proteins. Computational and biochemical approaches defined the RNA recognition element (RRE) as a tandem CAC trinucleotide motif separated by a variable spacer region. Similar to other mRNA-binding proteins, RBPMS family of proteins relocalized to cytoplasmic stress granules under oxidative stress conditions suggestive of a support function for mRNA localization in large and/or multinucleated cells where it is preferentially expressed., (© 2014 Farazi et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

CLP1 is a RNA kinase involved in tRNA splicing. Recently, CLP1 kinase-dead mice were shown to display a neuromuscular disorder with loss of motor neurons and muscle paralysis. Human genome analyses now identified a CLP1 homozygous missense mutation (p.R140H) in five unrelated families, leading to a loss of CLP1 interaction with the tRNA splicing endonuclease (TSEN) complex, largely reduced pre-tRNA cleavage activity, and accumulation of linear tRNA introns. The affected individuals develop severe motor-sensory defects, cortical dysgenesis, and microcephaly. Mice carrying kinase-dead CLP1 also displayed microcephaly and reduced cortical brain volume due to the enhanced cell death of neuronal progenitors that is associated with reduced numbers of cortical neurons. Our data elucidate a neurological syndrome defined by CLP1 mutations that impair tRNA splicing. Reduction of a founder mutation to homozygosity illustrates the importance of rare variations in disease and supports the clan genomics hypothesis., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Background: Various microRNAs (miRNAs) are up- or downregulated in tumors. However, the repression of cognate miRNA targets responsible for the phenotypic effects of this dysregulation in patients remains largely unexplored. To define miRNA targets and associated pathways, together with their relationship to outcome in breast cancer, we integrated patient-paired miRNA-mRNA expression data with a set of validated miRNA targets and pathway inference., Results: To generate a biochemically-validated set of miRNA-binding sites, we performed argonaute-2 photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation (AGO2-PAR-CLIP) in MCF7 cells. We then defined putative miRNA-target interactions using a computational model, which ranked and selected additional TargetScan-predicted interactions based on features of our AGO2-PAR-CLIP binding-site data. We subselected modeled interactions according to the abundance of their constituent miRNA and mRNA transcripts in tumors, and we took advantage of the variability of miRNA expression within molecular subtypes to detect miRNA repression. Interestingly, our data suggest that miRNA families control subtype-specific pathways; for example, miR-17, miR-19a, miR-25, and miR-200b show high miRNA regulatory activity in the triple-negative, basal-like subtype, whereas miR-22 and miR-24 do so in the HER2 subtype. An independent dataset validated our findings for miR-17 and miR-25, and showed a correlation between the expression levels of miR-182 targets and overall patient survival. Pathway analysis associated miR-17, miR-19a, and miR-200b with leukocyte transendothelial migration., Conclusions: We combined PAR-CLIP data with patient expression data to predict regulatory miRNAs, revealing potential therapeutic targets and prognostic markers in breast cancer.

Characteristic small RNA biogenesis processing patterns are used for the discovery of novel microRNAs (miRNAs) from next-generation sequencing data. Here, we highlight and discuss key criteria for mammalian - specifically human - miRNA database curation based on small RNA sequencing data. Sequence reads obtained from small RNA cDNA libraries are aligned to reference genomic regions, and miRNA genes are revealed by their distinct read length and bimodal read frequency distribution, the predicted secondary structure of the deduced miRNA stem-loop precursor molecule, and, to a lesser degree, based on evolutionary conservation of small RNAs from other vertebrates. Properly curated miRNA databases are an important resource for investigators interested in miRNA biology, diagnostics, and therapeutics.

Mature microRNAs (miRNAs) are single-stranded RNA molecules of 20-23-nucleotide (nt) length that control gene expression in many cellular processes. These molecules typically reduce the translation and stability of mRNAs, including those of genes that mediate processes in tumorigenesis, such as inflammation, cell cycle regulation, stress response, differentiation, apoptosis, and invasion. miRNA targeting is initiated through specific base-pairing interactions between the 5' end ("seed" region) of the miRNA and sites within coding and untranslated regions (UTRs) of mRNAs; target sites in the 3' UTR lead to more effective mRNA destabilization. Since miRNAs frequently target hundreds of mRNAs, miRNA regulatory pathways are complex. To provide a critical overview of miRNA dysregulation in cancer, we first discuss the methods currently available for studying the role of miRNAs in cancer and then review miRNA genomic organization, biogenesis, and mechanism of target recognition, examining how these processes are altered in tumorigenesis. Given the critical role miRNAs play in tumorigenesis processes and their disease specific expression, they hold potential as therapeutic targets and novel biomarkers.

The characterization of post-transcriptional gene regulation by small regulatory (20-30 nt) RNAs, particularly miRNAs and piRNAs, has become a major focus of research in recent years. A prerequisite for characterizing small RNAs is their identification and quantification across different developmental stages, and in normal and disease tissues, as well as model cell lines. Here we present a step-by-step protocol for generating barcoded small RNA cDNA libraries compatible with Illumina HiSeq sequencing, thereby facilitating miRNA and other small RNA profiling of large sample collections., (Copyright © 2012 Elsevier Inc. All rights reserved.)

The characterization of post-transcriptional gene regulation by small regulatory RNAs of 20-30 nt length, particularly miRNAs and piRNAs, has become a major focus of research in recent years. A prerequisite for the characterization of small RNAs is their identification and quantification across different developmental stages, normal and diseased tissues, as well as model cell lines. Here we present a step-by-step protocol for the bioinformatic analysis of barcoded cDNA libraries for small RNA profiling generated by Illumina sequencing, thereby facilitating miRNA and other small RNA profiling of large sample collections., (Copyright © 2012 Elsevier Inc. All rights reserved.)

FUS, EWSR1 and TAF15, constituting the FET protein family, are abundant, highly conserved RNA-binding proteins with important roles in oncogenesis and neuronal disease, yet their RNA targets and recognition elements are unknown. Using PAR-CLIP, we defined global RNA targets for all human FET proteins and two ALS-causing human FUS mutants. FET members showed similar binding profiles, whereas FUS mutants showed a drastically altered binding pattern, consistent with changes in subcellular localization.

MicroRNAs (miRNA) regulate many genes critical for tumorigenesis. We profiled miRNAs from 11 normal breast tissues, 17 noninvasive, 151 invasive breast carcinomas, and 6 cell lines by in-house-developed barcoded Solexa sequencing. miRNAs were organized in genomic clusters representing promoter-controlled miRNA expression and sequence families representing seed sequence-dependent miRNA target regulation. Unsupervised clustering of samples by miRNA sequence families best reflected the clustering based on mRNA expression available for this sample set. Clustering and comparative analysis of miRNA read frequencies showed that normal breast samples were separated from most noninvasive ductal carcinoma in situ and invasive carcinomas by increased miR-21 (the most abundant miRNA in carcinomas) and multiple decreased miRNA families (including miR-98/let-7), with most miRNA changes apparent already in the noninvasive carcinomas. In addition, patients that went on to develop metastasis showed increased expression of mir-423, and triple-negative breast carcinomas were most distinct from other tumor subtypes due to upregulation of the mir~17-92 cluster. However, absolute miRNA levels between normal breast and carcinomas did not reveal any significant differences. We also discovered two polymorphic nucleotide variations among the more abundant miRNAs miR-181a (T19G) and miR-185 (T16G), but we did not identify nucleotide variations expected for classical tumor suppressor function associated with miRNAs. The differentiation of tumor subtypes and prediction of metastasis based on miRNA levels is statistically possible but is not driven by deregulation of abundant miRNAs, implicating far fewer miRNAs in tumorigenic processes than previously suggested., (©2011 AACR.)

Mature microRNAs (miRNAs) are single-stranded RNA molecules of 20-23 nucleotide (nt) length that control gene expression in many cellular processes. These molecules typically reduce the stability of mRNAs, including those of genes that mediate processes in tumorigenesis, such as inflammation, cell cycle regulation, stress response, differentiation, apoptosis and invasion. miRNA targeting is mostly achieved through specific base-pairing interactions between the 5' end ('seed' region) of the miRNA and sites within coding and untranslated regions (UTRs) of mRNAs; target sites in the 3' UTR lead to more effective mRNA destabilization. Since miRNAs frequently target hundreds of mRNAs, miRNA regulatory pathways are complex. To provide a critical overview of miRNA dysregulation in cancer, we first discuss the methods currently available for studying the role of miRNAs in cancer and then review miRNA genomic organization, biogenesis and mechanism of target recognition, examining how these processes are altered in tumorigenesis. Given the critical role miRNAs play in tumorigenesis processes and their disease-specific expression, they hold potential as therapeutic targets and novel biomarkers., (Copyright © 2010 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.)

Several distinct classes of small RNAs, some newly identified, have been discovered to play important regulatory roles in diverse cellular processes. These classes include siRNAs, miRNAs, rasiRNAs and piRNAs. Each class binds to distinct members of the Argonaute/Piwi protein family to form ribonucleoprotein complexes that recognize partially, or nearly perfect, complementary nucleic acid targets, and that mediate a variety of regulatory processes, including transcriptional and post-transcriptional gene silencing. Based on the known relationship of Argonaute/Piwi proteins with distinct classes of small RNAs, we can now predict how many new classes of small RNAs or silencing processes remain to be discovered.

MyristoylCoA:protein N-myristoyltransferase (Nmt, EC 2.3.1.97), a member of the GCN5 acetyltransferase (GNAT) superfamily, is an essential eukaryotic enzyme that catalyzes covalent attachment of myristate (C14:0) to the N-terminal Gly of proteins involved in myriad cellular functions. The 2.5 A resolution structure of a ternary complex of Saccharomyces cerevisiae Nmt1p with a bound substrate peptide (GLYASKLA) and nonhydrolyzable myristoylCoA analogue [Farazi, T. A., et al. (2001) Biochemistry 40, 6335] was used as the basis for a series of mutagenesis experiments designed to define the enzyme's catalytic mechanism. The kinetic properties of an F170A/L171A Nmt mutant are consistent with the proposal that their main chain amides, located in a beta-bulge structure conserved among GNATs, function as an oxyanion hole to polarize the thioester carbonyl of bound myristoylCoA prior to subsequent nucleophilic attack. Removal of the two C-terminal residues (M454 and L455) produces a 300--400-fold reduction in the chemical transformation rate and converts the rate-limiting step from a step after the transformation to the transformation event itself. This finding is consistent with the main chain C-terminal carboxylate of L455 functioning as a catalytic base that abstracts a proton from the N-terminal Gly ammonium of the bound peptide to generate the nucleophilic amine. Mutating N169 and T205 in concert reduces the rate of the chemical transformation, supporting their role as components of an H-bonding network that facilitates attack of the Gly1 amine and stabilizes the tetrahedral intermediate.

MyristoylCoA:protein N-myristoyltransferase (Nmt) attaches myristate to the N-terminal Gly residue of proteins involved in a variety of signal transduction cascades, and other critical cellular functions. To gain insight about the structural basis of substrate recognition and catalysis, we determined the structures of a binary complex of Saccharomyces cerevisiae Nmt1p with myristoylCoA to 2.2 A resolution and of a ternary complex of Nmt1p with a nonhydrolyzable myristoylCoA analogue [S-(2-oxo)pentadecylCoA] and an octapeptide substrate (GLYASKLA) to 2.5 A resolution. The binary complex reveals how myristoylCoA alters the conformation of the enzyme to promote binding of both myristoylCoA and peptide and identifies the backbone amides of F170 and L171 as an oxyanion hole which polarizes the reactive thioester carbonyl. The ternary complex structure reveals details of the enzyme's peptide binding specificity and illuminates its mechanism of acyl transfer. The N-terminal Gly ammonium is positioned in close proximity to the C-terminal carboxylate of the protein, where it is poised to undergo the required deprotonation to an amine. In this conformation, the nucleophile is 6.3 A away from the thioester carbonyl. A catalytic mechanism is proposed whereby, once deprotonation is initiated, the N-terminal Gly amine can approximate the thioester carbonyl by rotating along Psi. This motion is facilitated by a H-bond network and leads to reaction between the glycine nitrogen nucleophile and the carbonyl. Loss of CoA from the tetrahedral intermediate may be facilitated by intramolecular H-bonding of the sulfur to the adenylamine of CoA. This affords a compact leaving group and lends a role for the observed bends in the CoA structure. The absolute requirement for Gly at the N-terminus of substrates is explained by the requirement for flexible rotation of its amine.

MyristoylCoA:protein N-myristoyltransferase is a member of the superfamily of GCN5-related N-acetyltransferases and catalyzes the covalent attachment of myristate to the N-terminal Gly residue of proteins with diverse functions. Saccharomyces cerevisiae Nmt1p is a monomeric protein with an ordered bi-bi reaction mechanism: myristoylCoA is bound prior to peptide substrate; after catalysis, CoA is released followed by myristoylpeptide. Analysis of the X-ray structure of Nmt1p with bound substrate analogues indicates that the active site contains an oxyanion hole and a catalytic base and that catalysis proceeds through the nucleophilic addition-elimination mechanism. To determine the rate-limiting step in the enzyme reaction, pre-steady-state kinetic analyses were performed using a new, sensitive nonradioactive assay that detects CoA. Multiple turnover quenched flow studies disclosed that a step after the chemical transformation limits the overall rate of the reaction. Multiple and single turnover analyses revealed that the rate for the chemical transformation step is 13.8+/-0.6 s(-1) while the slower steady-state phase is 0.10+/-0.01 s(-1). Stopped flow kinetic studies of substrate acquisition indicated that binding of myristoylCoA to the apo-enzyme occurs through at least a two-step process, with a fast phase rate of 3.2 x 10(8) M(-1) s(-1) and a slow phase rate of 23+/-2 s(-1) (defined at 5 degrees C). Binding of an octapeptide substrate, representing the N-terminal sequence of a known yeast N-myristoylprotein (Cnb1p), to a binary complex composed of Nmt1p and a nonhydrolyzable myristoylCoA analogue (S-(2-oxo)pentadecylCoA) has a second-order rate constant of 2.1+/-0.3 x 10(6) M(-1) s(-1) and a dissociation rate of 26+/-15 s(-1) (defined at 10 degrees C). These results are interpreted in light of the X-ray structures of this enzyme.

N-myristoyltransferase (Nmt) attaches myristate to the N-terminal glycine of many important eukaryotic and viral proteins. It is a target for anti-fungal and anti-viral therapy. We have determined the structure, to 2.9 A resolution, of a ternary complex of Saccharomyces cerevisiae Nmt1p with bound myristoylCoA and peptide substrate analogs. The model reveals structural features that define the enzyme's substrate specificities and regulate the ordered binding and release of substrates and products. A novel catalytic mechanism is proposed involving deprotonation of the N-terminal ammonium of a peptide substrate by the enzyme's C-terminal backbone carboxylate.

Saccharomyces cerevisiae contains four known acyl-CoA synthetases (fatty acid activation proteins, Faaps). Faa1p and Faa4p activate exogenously derived fatty acids. Acyl-CoA metabolism plays a critical role in regulating protein N-myristoylation by the essential enzyme, myristoyl-CoA:protein N-myristoyltransferase (Nmt1p). In this report, we have examined whether Faa1p and Faa4p have distinct roles in affecting protein N-myristoylation as cells transition from growth in rich media to a growth-arrested state during nutrient deprivation (stationary phase). The colony-forming potential of 10 isogenic strains was defined as a function of time spent in stationary phase. These strains contained either a wild type or mutant NMT1 allele, and wild type or null alleles of each FAA. Only the combination of the Nmt mutant (nmt451Dp; reduced affinity for myristoyl-CoA) and loss of Faa4p produced a dramatic loss of colony-forming units (CFU). The progressive millionfold reduction in CFU was associated with a deficiency in protein N-myristoylation that first appeared during logarithmic growth, worsened through the post-diauxic phase, and became extreme in stationary phase. Northern and Western blot analyses plus N-myristoyltransferase assays showed that Nmt is normally present only during the log and diauxic/post-diauxic periods, indicating that N-myristoylproteins present in stationary phase are "inherited" from these earlier phases. Moreover, FAA4 is the only FAA induced during the critical diauxic/early post-diauxic transition. Although substitution of nmt1-451D for NMT1 results in deficiencies in protein N-myristoylation, these deficiencies are modest and limited by compensatory responses that include augmented expression of nmt1-451D and precocious induction of FAA4 in log phase. Loss of Faa4p from nmt1-451D cells severely compromises their capacity to adequately myristoylate Nmt substrates prior to entry into stationary phase since none of the other Faaps are able to functionally compensate for its absence. To identify Nmt1p substrates that may affect maintenance of proliferative potential during stationary phase, we searched the yeast genome for known and putative N-myristoylproteins. Of the 64 genes found, 48 were successfully deleted in NMT1 cells. Removal of any one of the following nine substrates produced a loss of CFU similar to that observed in nmt1-451Dfaa4Delta cells: Arf1p, Arf2p, Sip2p, Van1p, Ptc2p, YBL049W (homology to Snf7p), YJR114W, YKR007W, and YMR077C. These proteins provide opportunities to further define the molecular mechanisms that regulate survival during stationary phase.