Structure and Distribution of Endogenous Nonecotropic Murine Leukemia Viruses in Wild Mice

Virtually all of our present understanding of endogenous murine leukemia viruses (MLVs) is based on studies with inbred mice. To develop a better understanding of the interaction between endogenous retroviruses and their hosts, we have carried out a systematic investigation of endogenous nonecotropic MLVs in wild mice. Species studied included four major subspecies of Mus musculus (M. m. castaneus, M. m. musculus, M. m. molossinus, and M. m. domesticus) as well as four common inbred laboratory strains (AKR/J, HRS/J, C3H/HeJ, and C57BL/6J). We determined the detailed distribution of nonecotropic proviruses in the mice by using both env- and long terminal repeat (LTR)-derived oligonucleotide probes specific for the three different groups of endogenous MLVs. The analysis indicated that proviruses that react with all of the specific probes are present in most wild mouse DNAs tested, in numbers varying from 1 or 2 to more than 50. Although in common inbred laboratory strains the linkage of group-specific sequences in env and the LTR of the proviruses is strict, proviruses which combine env and the LTR sequences from different groups were commonly observed in the wild-mouse subspecies. The “recombinant” nonecotropic proviruses in the mouse genomes were amplified by PCR, and their genetic and recombinant natures were determined. These proviruses showed extended genetic variation and provide a valuable probe for study of the evolutionary relationship between MLVs and the murine hosts.     

Multiple sources of contamination in samples from patients reported to have XMRV infection

Xenotropic murine leukemia virus (MLV)-related retrovirus (XMRV) was reported to be associated with prostate cancer by Urisman, et al. in 2006 and chronic fatigue syndrome (CFS) by Lombardi, et al. in 2009. To investigate this association, we independently evaluated plasma samples from 4 patients with CFS reported by Lombardi, et al. to have XMRV infection and from 5 healthy controls reported to be XMRV uninfected. We also analyzed viral sequences obtained from supernatants of cell cultures found to contain XMRV after coculture with 9 clinical samples from 8 patients. A qPCR assay capable of distinguishing XMRV from endogenous MLVs showed that the viral sequences detected in the CFS patient plasma behaved like endogenous MLVs and not XMRV. Single-genome sequences (N = 89) from CFS patient plasma were indistinguishable from endogenous MLVs found in the mouse genome that are distinct from XMRV. By contrast, XMRV sequences were detected by qPCR in 2 of the 5 plasma samples from healthy controls (sequencing of the qPCR product confirmed XMRV not MLV). Single-genome sequences (N = 234) from the 9 culture supernatants reportedly positive for XMRV were indistinguishable from XMRV sequences obtained from 22Rv1 and XMRV-contaminated 293T cell-lines. These results indicate that MLV DNA detected in the plasma samples from CFS patients evaluated in this study was from contaminating mouse genomic DNA and that XMRV detected in plasma samples from healthy controls and in cultures of patient samples was due to cross-contamination with XMRV (virus or nucleic acid).     

Folate Deficiency Provides Protection against Colon Carcinogenesis in DNA Polymerase β Haploinsufficient Mice

Aging and DNA polymerase β deficiency (β-pol^+/−) interact to accelerate the development of malignant lymphomas and adenocarcinoma and increase tumor bearing load in mice. Folate deficiency (FD) has been shown to induce DNA damage repaired via the base excision repair (BER) pathway. We anticipated that FD and BER deficiency would interact to accelerate aberrant crypt foci (ACF) formation and tumor development in β-pol haploinsufficient animals. FD resulted in a significant increase in ACF formation in wild type (WT) animals exposed to 1,2-dimethylhydrazine, a known colon and liver carcinogen; however, FD reduced development of ACF in β-pol haploinsufficient mice. Prolonged feeding of the FD diet resulted in advanced ACF formation and liver tumors in wild type mice. However, FD attenuated onset and progression of ACF and prevented liver tumorigenesis in β-pol haploinsufficient mice, i.e. FD provided protection against tumorigenesis in a BER-deficient environment in all tissues where 1,2-dimethylhydrazine exerts its damage. Here we show a distinct down-regulation in DNA repair pathways, e.g. BER, nucleotide excision repair, and mismatch repair, and decline in cell proliferation, as well as an up-regulation in poly(ADP-ribose) polymerase, proapoptotic genes, and apoptosis in colons of FD β-pol haploinsufficient mice.     

The Pre-NH2-Terminal Domain of the Herpes Simplex Virus 1 DNA Polymerase Catalytic Subunit Is Required for Efficient Viral Replication

The catalytic subunit of herpes simplex virus 1 DNA polymerase (HSV-1 Pol) has been extensively studied; however, its full complement of functional domains has yet to be characterized. A crystal structure has revealed a previously uncharacterized pre-NH₂-terminal domain (residues 1 to 140) within HSV-1 Pol. Due to the conservation of the pre-NH₂-terminal domain within the herpesvirus Pol family and its location in the crystal structure, we hypothesized that this domain provides an important function during viral replication in the infected cell distinct from 5′-3′ polymerase activity. We identified three pre-NH₂-terminal Pol mutants that exhibited 5′-3′ polymerase activity indistinguishable from that of wild-type Pol in vitro: deletion mutants PolΔN43 and PolΔN52 that lack the extreme N-terminal 42 and 51 residues, respectively, and mutant PolA₆, in which a conserved motif at residues 44 to 49 was replaced with alanines. We constructed the corresponding pol mutant viruses and found that the polΔN43 mutant displayed replication kinetics similar to those of wild-type virus, while polΔN52 and polA₆ mutant virus infection resulted in an 8-fold defect in viral yield compared to that achieved with wild type and their respective rescued derivative viruses. Additionally, both polΔN52 and polA₆ viruses exhibited defects in viral DNA synthesis that correlated with the observed reduction in viral yield. These results strongly indicate that the conserved motif within the pre-NH₂-terminal domain is important for viral DNA synthesis and production of infectious virus and indicate a functional role for this domain.     

The phylogeny of orthoretroviral long terminal repeats (LTRs)

LTRs are sequence elements in retroviruses and retrotransposons which are difficult to align due to their variability. One way of handling such cases is to use Hidden Markov Models (HMMs). In this work HMMs of LTRs were constructed for three groups of orthoretroviruses: the betaretroviruslike human MMTV-like (HML) endogenous retroviruses, the lentiviruses, including HIV, and gammaretroviruslike human endogenous retroviruses (HERVs). The HMM-generated LTR alignments and the phylogenetic trees constructed from them were compared with trees based on alignments of the pol gene at the nucleic acid level. The majority of branches in the LTR and pol based trees had the same order for the three retroviral genera, showing that HMM methods are successful in aligning and constructing phylogenies of LTRs. The HML LTR tree deviated somewhat from the pol tree for the groups HML3, HML7 and HML6. Among the gammaretroviruslike proviruses, the exogenous Mouse Leukemia Virus (MLV) was highly related to HERV-T in the pol based tree, but not in the LTR based tree. Aside from these differences, the similarity between the trees indicates that LTRs and pol coevolved in a largely monophyletic way.     

Developmental changes in histone macroH2A1-mediated gene regulation

macroH2A histone variants have been implicated to function in gene silencing by several studies, including ones showing a preferential association of macroH2A on the inactive X chromosome. To examine macroH2A function in vivo, we knocked out macroH2A1. macroH2A1 knockout mice are viable and fertile. A broad screen of liver gene expression showed no evidence of defects in X inactivation but did identify genes that have increased expression levels in macroH2A1 knockouts. macroH2A1-containing nucleosomes are enriched on the coding and/or upstream regions of these genes, suggesting that their increased expression levels are a direct effect of the absence of macroH2A1. The concentrations of macroH2A1 nucleosomes on these genes are low in the livers of newborn mice, and the macroH2A1 knockout had little effect on the expression levels of these genes in newborn liver. Our results indicate that an increase in liver macroH2A1 during the transition from newborn to young-adult status contributes to a decrease in the expression levels of these genes. These genes cluster in the area of lipid metabolism, and we observed metabolic effects in macroH2A1 knockouts. Our results indicate that the function of macroH2A1 histones is not restricted to gene silencing but also involves fine tuning the expression of specific genes.     

A Linkage Map of Endogenous Murine Leukemia Proviruses

Thirty endogenous proviruses belonging to the modified polytropic (Mpmv) class of murine leukemia virus (MLV) were identified by proviral-cellular DNA junction fragment segregation in several sets of recombinant inbred mice. Twenty-six Mpmv loci were mapped to chromosomal regions by matching proviral strain distribution patterns to those of previously assigned genes. Like other endogenous nonecotropic MLVs, Mpmv loci were present on several chromosomes in all strains examined. We pooled recombinant inbred strain linkage data from 110 MLV loci and selected marker genes in order to construct a chromosomal linkage map. Every mouse chromosome was found to harbor at least one proviral insertion, and several regions contained multiple integrations. However, the overall distribution of the 110 mapped proviruses did not deviate significantly from a random distribution. Because of their polymorphism in inbred strains of mice, and the ability to score as many as 57 proviruses per strain using only three hybridization probes, the nonecotropic MLVs mapped in common strains of mice offer a significant advantage over older methods (e.g., biochemical or individual restriction fragment polymorphisms) as genetic markers. These endogenous insertion elements should also be useful for assessing strain purity, and for studying the relatedness of common and not-so-common inbred strains.     

Endogenous murine leukemia retroviral variation across wild European and inbred strains of house mouse

Background

Endogenous murine leukemia retroviruses (MLVs) are high copy number proviral elements difficult to comprehensively characterize using standard low throughput sequencing approaches. However, high throughput approaches generate data that is challenging to process, interpret and present.

Results

Next generation sequencing (NGS) data was generated for MLVs from two wild caught Mus musculus domesticus (from mainland France and Corsica) and for inbred laboratory mouse strains C3H, LP/J and SJL. Sequence reads were grouped using a novel sequence clustering approach as applied to retroviral sequences. A Markov cluster algorithm was employed, and the sequence reads were queried for matches to specific xenotropic (Xmv), polytropic (Pmv) and modified polytropic (Mpmv) viral reference sequences.

Conclusions

Various MLV subtypes were more widespread than expected among the mice, which may be due to the higher coverage of NGS, or to the presence of similar sequence across many different proviral loci. The results did not correlate with variation in the major MLV receptor Xpr1, which can restrict exogenous MLVs, suggesting that endogenous MLV distribution may reflect gene flow more than past resistance to infection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1766-z) contains supplementary material, which is available to authorized users.     

Mice without MacroH2A Histone Variants

MacroH2A core histone variants have a unique structure that includes a C-terminal nonhistone domain. They are highly conserved in vertebrates and are thought to regulate gene expression. However, the nature of genes regulated by macroH2As and their biological significance remain unclear. Here, we examine macroH2A function in vivo by knocking out both macroH2A1 and macroH2A2 in the mouse. While macroH2As are not required for early development, the absence of macroH2As impairs prenatal and postnatal growth and can significantly reduce reproductive efficiency. The distributions of macroH2A.1- and macroH2A.2-containing nucleosomes show substantial overlap, as do their effects on gene expression. Our studies in fetal and adult liver indicate that macroH2As can exert large positive or negative effects on gene expression, with macroH2A.1 and macroH2A.2 acting synergistically on the expression of some genes and apparently having opposing effects on others. These effects are very specific and in the adult liver preferentially involve genes related to lipid metabolism, including the leptin receptor. MacroH2A-dependent gene regulation changes substantially in postnatal development and can be strongly affected by fasting. We propose that macroH2As produce adaptive changes to gene expression, which in the liver focus on metabolism.     

siRNAs containing 2′-fluorinated Northern-methanocarbacyclic (2′-F-NMC) nucleotides: in vitro and in vivo RNAi activity and inability of mitochondrial polymerases to incorporate 2′-F-NMC NTPs

We recently reported the synthesis of 2′-fluorinated Northern-methanocarbacyclic (2′-F-NMC) nucleotides, which are based on a bicyclo[3.1.0]hexane scaffold. Here, we analyzed RNAi-mediated gene silencing activity in cell culture and demonstrated that a single incorporation of 2′-F-NMC within the guide or passenger strand of the tri-N-acetylgalactosamine-conjugated siRNA targeting mouse Ttr was generally well tolerated. Exceptions were incorporation of 2′-F-NMC into the guide strand at positions 1 and 2, which resulted in a loss of the in vitro activity. Activity at position 1 was recovered when the guide strand was modified with a 5′ phosphate, suggesting that the 2′-F-NMC is a poor substrate for 5′ kinases. In mice, the 2′-F-NMC-modified siRNAs had comparable RNAi potencies to the parent siRNA. 2′-F-NMC residues in the guide seed region position 7 and at positions 10, 11 and 12 were well tolerated. Surprisingly, when the 5′-phosphate mimic 5′-(E)-vinylphosphonate was attached to the 2′-F-NMC at the position 1 of the guide strand, activity was considerably reduced. The steric constraints of the bicyclic 2′-F-NMC may impair formation of hydrogen-bonding interactions between the vinylphosphonate and the MID domain of Ago2. Molecular modeling studies explain the position- and conformation-dependent RNAi-mediated gene silencing activity of 2′-F-NMC. Finally, the 5′-triphosphate of 2′-F-NMC is not a substrate for mitochondrial RNA and DNA polymerases, indicating that metabolites should not be toxic.     

HIV-1 Suppressive Sequences Are Modulated by Rev Transport of Unspliced RNA and Are Required for Efficient HIV-1 Production

The unspliced human immunodeficiency virus type 1 (HIV-1) RNAs are translated as Gag and Gag-Pol polyproteins or packaged as genomes into viral particles. Efficient translation is necessary before the transition to produce infective virions. The viral protein Rev exports all intron-containing viral RNAs; however, it also appears to enhance translation. Cellular microRNAs target cellular and viral mRNAs to silence their translation and enrich them at discrete cytoplasmic loci that overlap with the putative interim site of Gag and the genome. Here, we analyzed how Rev-mediated transport and the splicing status of the mRNA influenced the silencing status imposed by microRNA. Through identification and mutational analysis of the silencing sites in the HIV-1 genome, we elucidated the effect of silencing on virus production. Renilla luciferase mRNA, which contains a let-7 targeting site in its 3′ untranslated region, was mediated when it was transported by Rev and not spliced, but it was either not mediated when it was spliced even in a partial way or it was Rev-independent. The silencing sites in the pol and env-nef regions of the HIV-1 genome, which were repressed in T cells and other cell lines, were Drosha-dependent and could also be modulated by Rev in an unspliced state. Mutant viruses that contained genomic mutations that reflect alterations to show more derepressive effects in the 3′ untranslated region of the Renilla luciferase gene replicated more slowly than wild-type virus. These findings yield insights into the HIV-1 silencing sites that might allow the genome to avoid translational machinery and that might be utilized in coordinating virus production during initial virus replication. However, the function of Rev to modulate the silencing sites of unspliced RNAs would be advantageous for the efficient translation that is required to support protein production prior to viral packaging and particle production.     

Analysis of acyclic nucleoside modifications in siRNAs finds sensitivity at position 1 that is restored by 5′-terminal phosphorylation both in vitro and in vivo

Small interfering RNAs (siRNAs) are short, double-stranded RNAs that use the endogenous RNAi pathway to mediate gene silencing. Phosphorylation facilitates loading of a siRNA into the Ago2 complex and subsequent cleavage of the target mRNA. In this study, 2′, 3′ seco nucleoside modifications, which contain an acylic ribose ring and are commonly called unlocked nucleic acids (UNAs), were evaluated at all positions along the guide strand of a siRNA targeting apolipoprotein B (ApoB). UNA modifications at positions 1, 2 and 3 were detrimental to siRNA activity. UNAs at positions 1 and 2 prevented phosphorylation by Clp1 kinase, abrogated binding to Ago2, and impaired Ago2-mediated cleavage of the mRNA target. The addition of a 5′-terminal phosphate to siRNA containing a position 1 UNA restored ApoB mRNA silencing, Ago2 binding, and Ago2 mediated cleavage activity. Position 1 UNA modified siRNA containing a 5′-terminal phosphate exhibited a partial restoration of siRNA silencing activity in vivo. These data reveal the complexity of interpreting the effects of chemical modification on siRNA activity, and exemplify the importance of using multiple biochemical, cell-based and in vivo assays to rationally design chemically modified siRNA destined for therapeutic use.     

siRNAs compete with miRNAs for methylation by HEN1 in Arabidopsis

Plant microRNAs (miRNAs) and small interfering RNAs (siRNAs) bear a 2′-O-methyl group on the 3′-terminal nucleotide. This methyl group is post-synthetically added by the methyltransferase protein HEN1 and protects small RNAs from enzymatic activities that target the 3′-OH. A mutagenesis screen for suppressors of the partial loss-of-function hen1-2 allele in Arabidopsis identified second-site mutations that restore miRNA methylation. These mutations affect two subunits of the DNA-dependent RNA polymerase IV (Pol IV), which is essential for the biogenesis of 24 nt endogenous siRNAs. A mutation in RNA-dependent RNA polymerase 2, another essential gene for the biogenesis of endogenous 24-nt siRNAs, also rescued the defects in miRNA methylation of hen1-2, revealing a previously unsuspected, negative influence of siRNAs on HEN1-mediated miRNA methylation. In addition, our findings imply the existence of a negative modifier of HEN1 activity in the Columbia genetic background.     

Genome-wide distribution of macroH2A1 histone variants in mouse liver chromatin

Studies of macroH2A histone variants indicate that they have a role in regulating gene expression. To identify direct targets of the macroH2A1 variants, we produced a genome-wide map of the distribution of macroH2A1 nucleosomes in mouse liver chromatin using high-throughput DNA sequencing. Although macroH2A1 nucleosomes are widely distributed across the genome, their local concentration varies over a range of 100-fold or more. The transcribed regions of most active genes are depleted of macroH2A1, often in sharply localized domains that show depletion of 4-fold or more relative to bulk mouse liver chromatin. We used macroH2A1 enrichment to help identify genes that appear to be directly regulated by macroH2A1 in mouse liver. These genes functionally cluster in the area of lipid metabolism. All but one of these genes has increased expression in macroH2A1 knockout mice, indicating that macroH2A1 functions primarily as a repressor in adult liver. This repressor activity is further supported by the substantial and relatively uniform macroH2A1 enrichment along the inactive X chromosome, which averages 4-fold. Genes that escape X inactivation stand out as domains of macroH2A1 depletion. The rarity of such genes indicates that few genes escape X inactivation in mouse liver, in contrast to what has been observed in human cells.     

2'-O-[2-[(N,N-dimethylamino)oxy]ethyl]-modified oligonucleotides inhibit expression of mRNA in vitro and in vivo

Synthesis and antisense activity of oligonucleotides modified with 2′-O-[2-[(N,N-dimethylamino)oxy] ethyl] (2′-O-DMAOE) are described. The 2′-O-DMAOE-modified oligonucleotides showed superior metabolic stability in mice. The phosphorothioate oligonucleotide ‘gapmers’, with 2′-O-DMAOE- modified nucleoside residues at the ends and 2′-deoxy nucleosides residues in the central region, showed dose-dependent inhibition of mRNA expression in cell culture for two targets. ‘Gapmer’ oligonucleotides have one or two 2′-O-modified regions and a 2′-deoxyoligonucleotide phosphorothioate region that allows RNase H digestion of target mRNA. To determine the in vivo potency and efficacy, BalbC mice were treated with 2′-O-DMAOE gapmers and a dose-dependent reduction in the targeted C-raf mRNA expression was observed. Oligonucleotides with 2′-O-DMAOE modifications throughout the sequences reduced the intercellular adhesion molecule-1 (ICAM-1) protein expression very efficiently in HUVEC cells with an IC₅₀ of 1.8 nM. The inhibition of ICAM-1 protein expression by these uniformly modified 2′-O-DMAOE oligonucleotides may be due to selective interference with the formation of the translational initiation complex. These results demonstrate that 2′-O-DMAOE- modified oligonucleotides are useful for antisense-based therapeutics when either RNase H-dependent or RNase H-independent target reduction mechanisms are employed.     

Regulation of small RNA stability: methylation and beyond

As central components of RNA silencing, small RNAs play diverse and important roles in many biological processes in eukaryotes. Aberrant reduction or elevation in the levels of small RNAs is associated with many developmental and physiological defects. The in vivo levels of small RNAs are precisely regulated through modulating the rates of their biogenesis and turnover. 2′-O-methylation on the 3′ terminal ribose is a major mechanism that increases the stability of small RNAs. The small RNA methyltransferase HUA ENHANCER1 (HEN1) and its homologs methylate microRNAs and small interfering RNAs (siRNAs) in plants, Piwi-interacting RNAs (piRNAs) in animals, and siRNAs in Drosophila. 3′ nucleotide addition, especially uridylation, and 3′-5′ exonucleolytic degradation are major mechanisms that turnover small RNAs. Other mechanisms impacting small RNA stability include complementary RNAs, cis-elements in small RNA sequences and RNA-binding proteins. Investigations are ongoing to further understand how small RNA stability impacts their accumulation in vivo in order to improve the utilization of RNA silencing in biotechnology and therapeutic applications.