Sequences within Both the 5′ UTR and Gag Are Required for Optimal In Vivo Packaging and Propagation of Mouse Mammary Tumor Virus (MMTV) Genomic RNA

Background

This study mapped regions of genomic RNA (gRNA) important for packaging and propagation of mouse mammary tumor virus (MMTV). MMTV is a type B betaretrovirus which preassembles intracellularly, a phenomenon distinct from retroviruses that assemble the progeny virion at cell surface just before budding such as the type C human and feline immunodeficiency viruses (HIV and FIV). Studies of FIV and Mason-Pfizer monkey virus (MPMV), a type D betaretrovirus with similar intracellular virion assembly processes as MMTV, have shown that the 5′ untranslated region (5′ UTR) and 5′ end of gag constitute important packaging determinants for gRNA.

Methodology

Three series of MMTV transfer vectors containing incremental amounts of gag or 5′ UTR sequences, or incremental amounts of 5′ UTR in the presence of 400 nucleotides (nt) of gag were constructed to delineate the extent of 5′ sequences that may be involved in MMTV gRNA packaging. Real time PCR measured the packaging efficiency of these vector RNAs into MMTV particles generated by co-transfection of MMTV Gag/Pol, vesicular stomatitis virus envelope glycoprotein (VSV-G Env), and individual transfer vectors into human 293T cells. Transfer vector RNA propagation was monitored by measuring transduction of target HeLaT4 cells following infection with viral particles containing a hygromycin resistance gene expression cassette on the packaged RNA.

Principal Findings

MMTV requires the entire 5′ UTR and a minimum of ∼120 nucleotide (nt) at the 5′ end of gag for not only efficient gRNA packaging but also propagation of MMTV-based transfer vector RNAs. Vector RNAs without the entire 5′ UTR were defective for both efficient packaging and propagation into target cells.

Conclusions/Significance

These results reveal that the 5′ end of MMTV genome is critical for both gRNA packaging and propagation, unlike the recently delineated FIV and MPMV packaging determinants that have been shown to be of bipartite nature.     

Identification of a New Target of miR-16, Vacuolar Protein Sorting 4a

Rationale

The rationale was to utilize a bioinformatics approach to identify miRNA binding sites in genes with single nucleotide mutations (SNPs) to discover pathways in heart failure (HF).

Objective

The objective was to focus on the genes containing miRNA binding sites with miRNAs that were significantly altered in end-stage HF and in response to a left ventricular assist device (LVAD).

Methods and Results

BEDTools v2.14.3 was used to discriminate SNPs within predicted 3′UTR miRNA binding sites. A member of the miR-15/107 family, miR-16, was decreased in the circulation of end-stage HF patients and increased in response to a LVAD (p<0.001). MiR-16 decreased Vacuolar Protein Sorting 4a (VPS4a) expression in HEK 293T cells (p<0.01). The SNP rs16958754 was identified in the miR-15/107 family binding site of VPS4a which abolished direct binding of miR-16 to the 3′UTR of VPS4a (p<0.05). VPS4a was increased in the circulation of end-stage HF patients (p<0.001), and led to a decrease in the number of HEK 293T cells in vitro (p<0.001).

Conclusions

We provide evidence that miR-16 decreases in the circulation of end-stage HF patients and increases with a LVAD. Modeling studies suggest that miR-16 binds to and decreases expression of VPS4a. Overexpression of VPS4a decreases cell number. Together, these experiments suggest that miR-16 and VPS4a expression are altered in end-stage HF and in response to unloading with a LVAD. This signaling pathway may lead to reduced circulating cell number in HF.     

The induction of microRNA targeting IRS-1 is involved in the development of insulin resistance under conditions of mitochondrial dysfunction in hepatocytes

Background

Mitochondrial dysfunction induces insulin resistance in myocytes via a reduction of insulin receptor substrate-1 (IRS-1) expression. However, the effect of mitochondrial dysfunction on insulin sensitivity is not understood well in hepatocytes. Although research has implicated the translational repression of target genes by endogenous non-coding microRNAs (miRNA) in the pathogenesis of various diseases, the identity and role of the miRNAs that are involved in the development of insulin resistance also remain largely unknown.

Methodology

To determine whether mitochondrial dysfunction induced by genetic or metabolic inhibition causes insulin resistance in hepatocytes, we analyzed the expression and insulin-stimulated phosphorylation of insulin signaling intermediates in SK-Hep1 hepatocytes. We used qRT-PCR to measure cellular levels of selected miRNAs that are thought to target IRS-1 3′ untranslated regions (3′UTR). Using overexpression of miR-126, we determined whether IRS-1-targeting miRNA causes insulin resistance in hepatocytes.

Principal Findings

Mitochondrial dysfunction resulting from genetic (mitochondrial DNA depletion) or metabolic inhibition (Rotenone or Antimycin A) induced insulin resistance in hepatocytes via a reduction in the expression of IRS-1 protein. In addition, we observed a significant up-regulation of several miRNAs presumed to target IRS-1 3′UTR in hepatocytes with mitochondrial dysfunction. Using reporter gene assay we confirmed that miR-126 directly targeted to IRS-1 3′UTR. Furthermore, the overexpression of miR-126 in hepatocytes caused a substantial reduction in IRS-1 protein expression, and a consequent impairment in insulin signaling.

Conclusions/Significance

We demonstrated that miR-126 was actively involved in the development of insulin resistance induced by mitochondrial dysfunction. These data provide novel insights into the molecular basis of insulin resistance, and implicate miRNA in the development of metabolic disease.     

MicroRNA-18a Attenuates DNA Damage Repair through Suppressing the Expression of Ataxia Telangiectasia Mutated in Colorectal Cancer

Background

miR-18a is one of the most up-regulated miRNAs in colorectal cancers (CRC) based on miRNA profiling. In this study, we examined the functional significance of miR-18a in CRC.

Methods

Expression of miR-18a was investigated in 45 CRC patients. Potential target genes of miR-18a were predicted by in silico search and confirmed by luciferase activity assay and Western blot. DNA damage was measured by comet assay. Gene function was measured by cell viability, colony formation and apoptosis assays.

Results

The up-regulation of miR-18a was validated and confirmed in 45 primary CRC tumors compared with adjacent normal tissues (p<0.0001). Through in silico search, the 3′UTR of Ataxia telangiectasia mutated (ATM) contains a conserved miR-18a binding site. Expression of ATM was down-regulated in CRC tumors (p<0.0001) and inversely correlated with miR-18a expression (r = -0.4562, p<0.01). Over-expression of miR-18a in colon cancer cells significantly reduced the luciferase activity of the construct with wild-type ATM 3′UTR but not that with mutant ATM 3′UTR, inferring a direct interaction of miR-18a with ATM 3′UTR. This was further confirmed by the down-regulation of ATM protein by miR-18a. As ATM is a key enzyme in DNA damage repair, we evaluated the effect of miR-18a on DNA double-strand breaks. Ectopic expression of miR-18a significantly inhibited the repair of DNA damage induced by etoposide (p<0.001), leading to accumulation of DNA damage, increase in cell apoptosis and poor clonogenic survival.

Conclusion

miR-18a attenuates cellular repair of DNA double-strand breaks by directly suppressing ATM, a key enzyme in DNA damage repair.     

Role of Polypyrimidine Tract Binding Protein in Mediating Internal Initiation of Translation of Interferon Regulatory Factor 2 RNA

Background

Earlier we have reported translational control of interferon regulatory factor 2 (IRF2) by internal initiation (Dhar et al, Nucleic Acids Res, 2007). The results implied possible role of IRF2 in controlling the intricate balance of cellular gene expression under stress conditions in general. Here we have investigated the secondary structure of the Internal Ribosome Entry Site of IRF2 RNA and demonstrated the role of PTB protein in ribosome assembly to facilitate internal initiation.

Methodology/Principal Findings

We have probed the putative secondary structure of the IRF2 5′UTR RNA using various enzymatic and chemical modification agents to constrain the secondary structure predicted from RNA folding algorithm Mfold. The IRES activity was found to be influenced by the interaction of trans-acting factor, polypyrimidine tract binding protein (PTB). Deletion of 25 nts from the 3′terminus of the 5′untranslated region resulted in reduced binding with PTB protein and also showed significant decrease in IRES activity compared to the wild type. We have also demonstrated putative contact points of PTB on the IRF2–5′UTR using primer extension inhibition assay. Majority of the PTB toe-prints were found to be restricted to the 3′end of the IRES. Additionally, Circular Dichroism (CD) spectra analysis suggested change in the conformation of the RNA upon PTB binding. Further, binding studies using S10 extract from HeLa cells, partially silenced for PTB gene expression, resulted in reduced binding by other trans-acting factors. Finally, we have demonstrated that addition of recombinant PTB enhances ribosome assembly on IRF2 IRES suggesting possible role of PTB in mediating internal initiation of translation of IRF2 RNA.

Conclusion/Significance

It appears that PTB binding to multiple sites within IRF2 5′UTR leads to a conformational change in the RNA that facilitate binding of other trans-acting factors to mediate internal initiation of translation.     

Genetic Variants in RKIP Are Associated with Clear Cell Renal Cell Carcinoma Risk in a Chinese Population

Background

Raf-1 kinase inhibitor protein (RKIP) plays a critical role in tumor development by regulating cell functions such as invasion, apoptosis and differentiation. Down-regulation of RKIP expression has been implicated in the development and progression of renal cell carcinoma (RCC). Herein, we hypothesized that genetic polymorphisms in RKIP might be associated with susceptibility and progression of RCC.

Methods

A total of 5 tagging single-nucleotide polymorphisms (tSNPs) in RKIP were selected and genotyped by SNapShot method in a case-control study of 859 RCC patients and 1004 controls. The logistic regression was used to evaluate the genetic association with occurrence and progression of RCC. The functionality of the important SNP was preliminary examined by qRT-PCR.

Result

We found that the rs17512051 in the promoter region of RKIP was significantly associated with decreased clear cell RCC (ccRCC) risk (TA/AA vs. TT: P = 0.039, OR = 0.78, 95%CI = 0.62–0.99). Another SNP (rs1051470) in the 3′UTR region of RKIP was marginally associated with increased ccRCC risk (TT vs. CC+CT: OR = 1.45, 95%CI = 1.01–2.09). In the stratified analysis, the protective effect of rs17512051 was more predominant in the subgroups of male, non-smokers, non-drinkers as well as subjects without history of diabetes. Furthermore, we observed higher RKIP mRNA levels in the presence of the rs17512051A allele in normal renal tissues.

Conclusion

Our results suggest that the potentially functional RKIP rs17512051 polymorphism may affect ccRCC susceptibility through altering the endogenous RKIP expression level. Risk effects and the functional impact of this polymorphism need further validation.     

Alternative Processing of Primary microRNA Transcripts by Drosha Generates 5′ End Variation of Mature microRNA

Background

It is generally believed that the miRNA processing machinery ensures the generation of a mature miRNA with a fixed sequence, particularly at its 5′ end. However, we and others have recently noted that the ends of a given mature miRNA are not absolutely fixed, but subject to variation. Neither the significance nor the mechanism behind the generation of such miRNA polymorphism is understood. miR-142 is an abundantly expressed miRNA in hematopoietic cells and exhibits a high frequency of 5′ end polymorphism.

Methodology/Principal Findings

Here we show that a shift in the Drosha processing of pri-miRNA generates multiple forms of miR-142s in vivo with differing 5′ ends that might target different genes. Sequence analysis of several pre-miRNA ends cloned from T cells reveals that unlike many other pri-miRNAs that are processed into a single pre-miRNA, pri-miR-142 is processed into 3 distinct pre-miR-142s. Dicer processing studies suggest that each of the 3 pre-miR-142s is processed into a distinct double-stranded miRNA, giving rise to 4 mature miRNA variants that might regulate different target gene pools.

Conclusions/Significance

Thus, alternative Drosha processing might be a novel mechanism for diversification of the miRNA target gene pool.     

Effect of β2-adrenergic receptor gene (ADRB2) 3′ untranslated region polymorphisms on inhaled corticosteroid/long-acting β2-adrenergic agonist response

Background

Evidence suggests that variation in the length of the poly-C repeat in the 3′ untranslated region (3′UTR) of the β₂-adrenergic receptor gene (ADRB2) may contribute to interindividual variation in β-agonist response. However, methodology in previous studies limited the assessment of the effect of sequence variation in the context of poly-C repeat length. The objectives of this study were to design a novel genotyping method to fully characterize sequence variation in the ADRB2 3′UTR poly-C repeat in asthma patients treated with inhaled corticosteroid and long-acting β₂-adrenergic agonist (ICS/LABA) combination therapy, and to analyze the effect of the poly-C repeat polymorphism on clinical response.

Methods

In 2,250 asthma patients randomized to treatment with budesonide/formoterol or fluticasone/salmeterol in a six-month study (AstraZeneca study code: SD-039-0735), sequence diversity in the ADRB2 poly-C repeat region was determined using a novel sequencing-based genotyping method. The relationship between the poly-C repeat polymorphism and the incidence of severe asthma exacerbations, and changes in pulmonary function and asthma symptoms from baseline to the average during the treatment period, were analyzed.

Results

Poly-C repeat genotypes were assigned in 97% (2,192/2,250) of patients. Of the 13 different poly-C repeat alleles identified, six alleles occurred at a frequency of >5% in one or more population in this study. The repeat length of these six common alleles ranged from 10 to 14 nucleotides. Twelve poly-C repeat genotypes were observed at a frequency of >1%. No evidence of an association between poly-C repeat genotype and the incidence of severe asthma exacerbations was observed. Patients’ pulmonary function measurements improved and asthma symptoms declined when treated with ICS/LABA combination therapy regardless of poly-C repeat genotype.

Conclusions

The extensive sequence diversity present in the poly-C repeat region of the ADRB2 3′UTR did not predict therapeutic response to ICS/LABA therapy.     

Rare BRCA1 haplotypes including 3′UTR SNPs associated with breast cancer risk

Genetic markers identifying women at an increased risk of developing breast cancer exist, yet the majority of inherited risk remains elusive. While numerous BRCA1 coding sequence mutations are associated with breast cancer risk, BRCA1 mutations account for less then 5% of breast cancer risk. Since 3′ untranslated region (3′UTR) polymorphisms disrupting microRNA (miRNA) binding can be functional and can act as genetic markers of cancer risk, we tested the hypothesis that such polymorphisms in the 3′UTR of BRCA1 and haplotypes containing these functional polymorphisms may be associated with breast cancer risk. We sequenced the BRCA1 3′UTR from breast cancer patients to identify miRNA disrupting polymorphisms. We further evaluated haplotypes of this region including the identified 3′UTR variants in a large population of controls and breast cancer patients (n = 221) with known breast cancer subtypes and ethnicities. We identified three 3′UTR variants in BRCA1 that are polymorphic in breast cancer populations, and haplotype analysis including these variants revealed that breast cancer patients harbor five rare haplotypes not generally found among controls (9.50% for breast cancer chromosomes, 0.11% for control chromosomes, p = 0.0001). Three of these rare haplotypes contain the rs8176318 BRCA1 3′UTR functional variant. These haplotypes are not biomarkers for BRCA1 coding mutations, as they are found rarely in BRCA1 mutant breast cancer patients (1/129 patients = 0.78%). These rare BRCA1 haplotypes and 3′UTR SNPs may represent new genetic markers of breast cancer risk.Key words: BRCA1, haplotype, microRNA, SNP, 3′UTR, breast cancer, triple negative breast cancer     

DAZL Relieves miRNA-Mediated Repression of Germline mRNAs by Controlling Poly(A) Tail Length in Zebrafish

Background

During zebrafish embryogenesis, microRNA (miRNA) miR-430 contributes to restrict Nanos1 and TDRD7 to primordial germ cells (PGCs) by inducing mRNA deadenylation, mRNA degradation, and translational repression of nanos1 and tdrd7 mRNAs in somatic cells. The nanos1 and tdrd7 3′UTRs include cis-acting elements that allow activity in PGCs even in the presence of miRNA-mediated repression.

Methodology/Principal Findings

Using a GFP reporter mRNA that was fused with tdrd7 3′UTR, we show that a germline-specific RNA-binding protein DAZ-like (DAZL) can relieve the miR-430-mediated repression of tdrd7 mRNA by inducing poly(A) tail elongation (polyadenylation) in zebrafish. We also show that DAZL enhances protein synthesis via the 3′UTR of dazl mRNA, another germline mRNA targeted by miR-430.

Conclusions/Significance

Our present study indicated that DAZL acts as an “anti-miRNA factor” during vertebrate germ cell development. Our data also suggested that miRNA-mediated regulation can be modulated on specific target mRNAs through the poly(A) tail control.     

A Transposon in Comt Generates mRNA Variants and Causes Widespread Expression and Behavioral Differences among Mice

Background

Catechol-O-methyltransferase (COMT) is a key enzyme responsible for the degradation of dopamine and norepinephrine. COMT activity influences cognitive and emotional states in humans and aggression and drug responses in mice. This study identifies the key sequence variant that leads to differences in Comt mRNA and protein levels among mice, and that modulates synaptic function and pharmacological and behavioral traits.

Methodology/Principal Findings

We examined Comt expression in multiple tissues in over 100 diverse strains and several genetic crosses. Differences in expression map back to Comt and are generated by a 230 nt insertion of a B2 short interspersed element (B2 SINE) in the proximal 3′ UTR of Comt in C57BL/6J. This transposon introduces a premature polyadenylation signal and creates a short 3′ UTR isoform. The B2 SINE is shared by a subset of strains, including C57BL/6J, A/J, BALB/cByJ, and AKR/J, but is absent in others, including DBA/2J, FVB/NJ, SJL/J, and wild subspecies. The short isoform is associated with increased protein expression in prefrontal cortex and hippocampus relative to the longer ancestral isoform. The Comt variant causes downstream differences in the expression of genes involved in synaptic function, and also modulates phenotypes such as dopamine D1 and D2 receptor binding and pharmacological responses to haloperidol.

Conclusions/Significance

We have precisely defined the B2 SINE as the source of variation in Comt and demonstrated that a transposon in a 3′ UTR can alter mRNA isoform use and modulate behavior. The recent fixation of the variant in a subset of strains may have contributed to the rapid divergence of inbred strains.     

Impact of GC content on gene expression pattern in chicken

Background

GC content varies greatly between different genomic regions in many eukaryotes. In order to determine whether this organization named isochore organization influences gene expression patterns, the relationship between GC content and gene expression has been investigated in man and mouse. However, to date, this question is still a matter for debate. Among the avian species, chicken (Gallus gallus) is the best studied representative with a complete genome sequence. The distinctive features and organization of its sequence make it a good model to explore important issues in genome structure and evolution.

Methods

Only nuclear genes with complete information on protein-coding sequence with no evidence of multiple-splicing forms were included in this study. Chicken protein coding sequences, complete mRNA sequences (or full length cDNA sequences), and 5^′ untranslated region sequences (5^′ UTR) were downloaded from Ensembl and chicken expression data originated from a previous work. Three indices i.e. expression level, expression breadth and maximum expression level were used to measure the expression pattern of a given gene. CpG islands were identified using hgTables of the UCSC Genome Browser. Correlation analysis between variables was performed by SAS Proprietary Software Release 8.1.

Results

In chicken, the GC content of 5^′ UTR is significantly and positively correlated with expression level, expression breadth, and maximum expression level, whereas that of coding sequences and introns and at the third coding position are negatively correlated with expression level and expression breadth, and not correlated with maximum expression level. These significant trends are independent of recombination rate, chromosome size and gene density. Furthermore, multiple linear regression analysis indicated that GC content in genes could explain approximately 10% of the variation in gene expression.

Conclusions

GC content is significantly associated with gene expression pattern and could be one of the important regulation factors in the chicken genome.