A Viral (Arc)hive for Metazoan Memory.

Arc, a master regulator of synaptic plasticity, contains sequence elements that are evolutionarily related to retrotransposon Gag genes. Two related papers in this issue of Cell show that Arc retains retroviral-like capsid-forming ability and can transmit mRNA between cells in the nervous system, a process that may be important for synaptic function.

A hominoid-specific endogenous retrovirus may have rewired the gene regulatory network shared between primordial germ cells and naïve pluripotent cells.

Mammalian germ cells stem from primordial germ cells (PGCs). Although the gene regulatory network controlling the development of germ cells such as PGCs is critical for ensuring gamete integrity, substantial differences exist in this network among mammalian species, suggesting that this network has been modified during mammalian evolution. Here, we show that a hominoid-specific group of endogenous retroviruses, LTR5_Hs, discloses enhancer-like signatures in human in vitro-induced PGCs, PGC-like cells (PGCLCs). Human PGCLCs exhibit a transcriptome signature similar to that of naïve-state pluripotent cells. LTR5_Hs are epigenetically activated in both PGCLCs and naïve pluripotent cells, and the expression of genes in the vicinity of LTR5_Hs is coordinately upregulated in these cell types, contributing to the establishment of the transcriptome similarity between these cell types. LTR5_Hs are preferentially bound by transcription factors that are highly expressed in both PGCLCs and naïve pluripotent cells (KLF4, TFAP2C, NANOG, and CBFA2T2), suggesting that these transcription factors contribute to the epigenetic activation of LTR5_Hs in these cells. Comparative transcriptome analysis between humans and macaques suggests that the expression of many genes in PGCLCs and naïve pluripotent cells is upregulated by LTR5_Hs insertions in the hominoid lineage. Together, this study suggests that LTR5_Hs insertions may have finetuned the gene regulatory network shared between PGCLCs and naïve pluripotent cells and coordinately altered the gene expression in these cells during hominoid evolution.

Virus-derived variation in diverse human genomes.

Acquisition of genetic material from viruses by their hosts can generate inter-host structural genome variation. We developed computational tools enabling us to study virus-derived structural variants (SVs) in population-scale whole genome sequencing (WGS) datasets and applied them to 3,332 humans. Although SVs had already been cataloged in these subjects, we found previously-overlooked virus-derived SVs. We detected non-germline SVs derived from squirrel monkey retrovirus (SMRV), human immunodeficiency virus 1 (HIV-1), and human T lymphotropic virus (HTLV-1); these variants are attributable to infection of the sequenced lymphoblastoid cell lines (LCLs) or their progenitor cells and may impact gene expression results and the biosafety of experiments using these cells. In addition, we detected new heritable SVs derived from human herpesvirus 6 (HHV-6) and human endogenous retrovirus-K (HERV-K). We report the first solo-direct repeat (DR) HHV-6 likely to reflect DR rearrangement of a known full-length endogenous HHV-6. We used linkage disequilibrium between single nucleotide variants (SNVs) and variants in reads that align to HERV-K, which often cannot be mapped uniquely using conventional short-read sequencing analysis methods, to locate previously-unknown polymorphic HERV-K loci. Some of these loci are tightly linked to trait-associated SNVs, some are in complex genome regions inaccessible by prior methods, and some contain novel HERV-K haplotypes likely derived from gene conversion from an unknown source or introgression. These tools and results broaden our perspective on the coevolution between viruses and humans, including ongoing virus-to-human gene transfer contributing to genetic variation between humans.

Virus-like insertions with sequence signatures similar to those of endogenous nonretroviral RNA viruses in the human genome.

Understanding the genetics and taxonomy of ancient viruses will give us great insights into not only the origin and evolution of viruses but also how viral infections played roles in our evolution. Endogenous viruses are remnants of ancient viral infections and are thought to retain the genetic characteristics of viruses from ancient times. In this study, we used machine learning of endogenous RNA virus sequence signatures to identify viruses in the human genome that have not been detected or are already extinct. Here, we show that the k-mer occurrence of ancient RNA viral sequences remains similar to that of extant RNA viral sequences and can be differentiated from that of other human genome sequences. Furthermore, using this characteristic, we screened RNA viral insertions in the human reference genome and found virus-like insertions with phylogenetic and evolutionary features indicative of an exogenous origin but lacking homology to previously identified sequences. Our analysis indicates that animal genomes still contain unknown virus-derived sequences and provides a glimpse into the diversity of the ancient virosphere.

The regulation of persistent Borna disease virus infection by RNA silencing factors in human cells.

Viral infection induces diverse cellular immune responses. Some viruses induce the production of antiviral cytokines, alterations of endogenous gene expression, and apoptosis; however, other viruses replicate without inducing such responses, enabling them to persistently infect cells. Infection by Borna disease virus type 1 (BoDV-1) can result in fatal immune-mediated encephalitis, including in humans, yet infection of cells in vitro is generally persistent. The regulatory mechanisms underlying this persistent infection remain unclear. Here, we show that an enhancer of RNA-silencing, TRBP, positively regulates BoDV RNA level in human cells. Knockdown of TRBP decreased BoDV RNA levels in persistently-infected cells, whereas overexpression of TRBP increased BoDV RNA levels. To investigate the mechanism underlying this phenomenon, we performed immunoprecipitation assays and found that TRBP interacts with BoDV RNA. Furthermore, we performed cell fractionation, which revealed that persistent infection with BoDV does not alter the localization of TRBP and other RNA silencing factors in cells. Our results showed the regulation of persistent BoDV infection by RNA-silencing factors in human cells.

Mammalian antiviral systems directed by small RNA.

There are strong incentives for human populations to develop antiviral systems. Similarly, genomes that encode antiviral systems have had strong selective advantages. Protein-guided immune systems, which have been well studied in mammals, are necessary for survival in our virus-laden environments. Small RNA-directed antiviral immune systems suppress invasion of cells by non-self genetic material via complementary base pairing with target sequences. These RNA silencing-dependent systems operate in diverse organisms. In mammals, there is strong evidence that microRNAs (miRNAs) regulate endogenous genes important for antiviral immunity, and emerging evidence that virus-derived nucleic acids can be directly targeted by small interfering RNAs (siRNAs), PIWI-interacting RNAs (piRNAs), and transfer RNAs (tRNAs) for protection in some contexts. In this review, we summarize current knowledge of the antiviral functions of each of these small RNA types and consider their conceptual and mechanistic overlap with innate and adaptive protein-guided immunity, including mammalian antiviral cytokines, as well as the prokaryotic RNA-guided immune system, CRISPR. In light of recent successes in delivery of RNA for antiviral purposes, most notably for vaccination, we discuss the potential for development of small noncoding RNA-directed antiviral therapeutics and prophylactics.

Endogenization and excision of human herpesvirus 6 in human genomes.

Sequences homologous to human herpesvirus 6 (HHV-6) are integrated within the nuclear genome of about 1% of humans, but it is not clear how this came about. It is also uncertain whether integrated HHV-6 can reactivate into an infectious virus. HHV-6 integrates into telomeres, and this has recently been associated with polymorphisms affecting MOV10L1. MOV10L1 is located on the subtelomere of chromosome 22q (chr22q) and is required to make PIWI-interacting RNAs (piRNAs). As piRNAs block germline integration of transposons, piRNA-mediated repression of HHV-6 integration has been proposed to explain this association. In vitro, recombination of the HHV-6 genome along its terminal direct repeats (DRs) leads to excision from the telomere and viral reactivation, but the expected "solo-DR scar" has not been described in vivo. Here we screened for integrated HHV-6 in 7,485 Japanese subjects using whole-genome sequencing (WGS). Integrated HHV-6 was associated with polymorphisms on chr22q. However, in contrast to prior work, we find that the reported MOV10L1 polymorphism is physically linked to an ancient endogenous HHV-6A variant integrated into the telomere of chr22q in East Asians. Unexpectedly, an HHV-6B variant has also endogenized in chr22q; two endogenous HHV-6 variants at this locus thus account for 72% of all integrated HHV-6 in Japan. We also report human genomes carrying only one portion of the HHV-6B genome, a solo-DR, supporting in vivo excision and possible viral reactivation. Together these results explain the recently-reported association between integrated HHV-6 and MOV10L1/piRNAs, suggest potential exaptation of HHV-6 in its coevolution with human chr22q, and clarify the evolution and risk of reactivation of the only intact (non-retro)viral genome known to be present in human germlines.

Evolutionary History of Endogenous Human Herpesvirus 6 Reflects Human Migration out of Africa.

Human herpesvirus 6A and 6B (HHV-6) can integrate into the germline, and as a result, ∼70 million people harbor the genome of one of these viruses in every cell of their body. Until now, it has been largely unknown if 1) these integrations are ancient, 2) if they still occur, and 3) whether circulating virus strains differ from integrated ones. Here, we used next-generation sequencing and mining of public human genome data sets to generate the largest and most diverse collection of circulating and integrated HHV-6 genomes studied to date. In genomes of geographically dispersed, only distantly related people, we identified clades of integrated viruses that originated from a single ancestral event, confirming this with fluorescent in situ hybridization to directly observe the integration locus. In contrast to HHV-6B, circulating and integrated HHV-6A sequences form distinct clades, arguing against ongoing integration of circulating HHV-6A or "reactivation" of integrated HHV-6A. Taken together, our study provides the first comprehensive picture of the evolution of HHV-6, and reveals that integration of heritable HHV-6 has occurred since the time of, if not before, human migrations out of Africa.

piRNA-Guided CRISPR-like Immunity in Eukaryotes.

Eukaryotic genomes contain virus-derived sequences called endogenous virus elements (EVEs). The majority of EVEs are related to retroviruses, which integrate into the host genome in order to replicate. Some retroviral EVEs encode a function; for example, some produce proteins that block infection by related viruses. EVEs derived from nonretroviral viruses - also recently found in many eukaryotic genomes - are more enigmatic. Here, we summarize the evidence that EVEs can act as templates to generate Piwi-interacting RNAs (piRNAs), whose canonical function is sequence-specific silencing of transposable elements (TEs) to maintain genomic integrity. We argue that EVEs may thus enable heritable, sequence-specific antiviral immune memory in eukaryotes - analogous to CRISPR-Cas immunity in prokaryotes.

Comprehensive discovery of CRISPR-targeted terminally redundant sequences in the human gut metagenome: Viruses, plasmids, and more.

Viruses are the most numerous biological entity, existing in all environments and infecting all cellular organisms. Compared with cellular life, the evolution and origin of viruses are poorly understood; viruses are enormously diverse, and most lack sequence similarity to cellular genes. To uncover viral sequences without relying on either reference viral sequences from databases or marker genes that characterize specific viral taxa, we developed an analysis pipeline for virus inference based on clustered regularly interspaced short palindromic repeats (CRISPR). CRISPR is a prokaryotic nucleic acid restriction system that stores the memory of previous exposure. Our protocol can infer CRISPR-targeted sequences, including viruses, plasmids, and previously uncharacterized elements, and predict their hosts using unassembled short-read metagenomic sequencing data. By analyzing human gut metagenomic data, we extracted 11,391 terminally redundant CRISPR-targeted sequences, which are likely complete circular genomes. The sequences included 2,154 tailed-phage genomes, together with 257 complete crAssphage genomes, 11 genomes larger than 200 kilobases, 766 genomes of Microviridae species, 56 genomes of Inoviridae species, and 95 previously uncharacterized circular small genomes that have no reliably predicted protein-coding gene. We predicted the host(s) of approximately 70% of the discovered genomes at the taxonomic level of phylum by linking protospacers to taxonomically assigned CRISPR direct repeats. These results demonstrate that our protocol is efficient for de novo inference of CRISPR-targeted sequences and their host prediction.

piRNAs derived from ancient viral processed pseudogenes as transgenerational sequence-specific immune memory in mammals.

Endogenous bornavirus-like nucleoprotein elements (EBLNs) are sequences within vertebrate genomes derived from reverse transcription and integration of ancient bornaviral nucleoprotein mRNA via the host retrotransposon machinery. While species with EBLNs appear relatively resistant to bornaviral disease, the nature of this association is unclear. We hypothesized that EBLNs could give rise to antiviral interfering RNA in the form of PIWI-interacting RNAs (piRNAs), a class of small RNA known to silence transposons but not exogenous viruses. We found that in both rodents and primates, which acquired their EBLNs independently some 25-40 million years ago, EBLNs are present within piRNA-generating regions of the genome far more often than expected by chance alone (ℙ = 8 × 10(-3)-6 × 10(-8)). Three of the seven human EBLNs fall within annotated piRNA clusters and two marmoset EBLNs give rise to bona fide piRNAs. In both rats and mice, at least two of the five EBLNs give rise to abundant piRNAs in the male gonad. While no EBLNs are syntenic between rodent and primate, some of the piRNA clusters containing EBLNs are; thus we deduce that EBLNs were integrated into existing piRNA clusters. All true piRNAs derived from EBLNs are antisense relative to the proposed ancient bornaviral nucleoprotein mRNA. These observations are consistent with a role for EBLN-derived piRNA-like RNAs in interfering with ancient bornaviral infection. They raise the hypothesis that retrotransposon-dependent virus-to-host gene flow could engender RNA-mediated, sequence-specific antiviral immune memory in metazoans analogous to the CRISPR/Cas system in prokaryotes.

Phenotypic properties of transmitted founder HIV-1.

Defining the virus-host interactions responsible for HIV-1 transmission, including the phenotypic requirements of viruses capable of establishing de novo infections, could be important for AIDS vaccine development. Previous analyses have failed to identify phenotypic properties other than chemokine receptor 5 (CCR5) and CD4+ T-cell tropism that are preferentially associated with viral transmission. However, most of these studies were limited to examining envelope (Env) function in the context of pseudoviruses. Here, we generated infectious molecular clones of transmitted founder (TF; n = 27) and chronic control (CC; n = 14) viruses of subtypes B (n = 18) and C (n = 23) and compared their phenotypic properties in assays specifically designed to probe the earliest stages of HIV-1 infection. We found that TF virions were 1.7-fold more infectious (P = 0.049) and contained 1.9-fold more Env per particle (P = 0.048) compared with CC viruses. TF viruses were also captured by monocyte-derived dendritic cells 1.7-fold more efficiently (P = 0.035) and more readily transferred to CD4+ T cells (P = 0.025). In primary CD4+ T cells, TF and CC viruses replicated with comparable kinetics; however, when propagated in the presence of IFN-α, TF viruses replicated to higher titers than CC viruses. This difference was significant for subtype B (P = 0.000013) but not subtype C (P = 0.53) viruses, possibly reflecting demographic differences of the respective patient cohorts. Together, these data indicate that TF viruses are enriched for higher Env content, enhanced cell-free infectivity, improved dendritic cell interaction, and relative IFN-α resistance. These viral properties, which likely act in concert, should be considered in the development and testing of AIDS vaccines.

Transmitted/founder and chronic HIV-1 envelope proteins are distinguished by differential utilization of CCR5.

Infection by HIV-1 most often results from the successful transmission and propagation of a single virus variant, termed the transmitted/founder (T/F) virus. Here, we compared the attachment and entry properties of envelope (Env) glycoproteins from T/F and chronic control (CC) viruses. Using a panel of 40 T/F and 47 CC Envs, all derived by single genome amplification, we found that 52% of clade C and B CC Envs exhibited partial resistance to the CCR5 antagonist maraviroc (MVC) on cells expressing high levels of CCR5, while only 15% of T/F Envs exhibited this same property. Moreover, subtle differences in the magnitude with which MVC inhibited infection on cells expressing low levels of CCR5, including primary CD4(+) T cells, were highly predictive of MVC resistance when CCR5 expression levels were high. These results are consistent with previous observations showing a greater sensitivity of T/F Envs to MVC inhibition on cells expressing very high levels of CCR5 and indicate that CC Envs are often capable of recognizing MVC-bound CCR5, albeit inefficiently on cells expressing physiologic levels of CCR5. When CCR5 expression levels are high, this phenotype becomes readily detectable. The utilization of drug-bound CCR5 conformations by many CC Envs was seen with other CCR5 antagonists, with replication-competent viruses, and did not obviously correlate with other phenotypic traits. The striking ability of clade C and B CC Envs to use MVC-bound CCR5 relative to T/F Envs argues that the more promiscuous use of CCR5 by these Env proteins is selected against at the level of virus transmission and is selected for during chronic infection.

Transmitted/founder and chronic subtype C HIV-1 use CD4 and CCR5 receptors with equal efficiency and are not inhibited by blocking the integrin α4ß7.

Sexual transmission of human immunodeficiency virus type 1 (HIV-1) most often results from productive infection by a single transmitted/founder (T/F) virus, indicating a stringent mucosal bottleneck. Understanding the viral traits that overcome this bottleneck could have important implications for HIV-1 vaccine design and other prevention strategies. Most T/F viruses use CCR5 to infect target cells and some encode envelope glycoproteins (Envs) that contain fewer potential N-linked glycosylation sites and shorter V1/V2 variable loops than Envs from chronic viruses. Moreover, it has been reported that the gp120 subunits of certain transmitted Envs bind to the gut-homing integrin α4ß7, possibly enhancing virus entry and cell-to-cell spread. Here we sought to determine whether subtype C T/F viruses, which are responsible for the majority of new HIV-1 infections worldwide, share biological properties that increase their transmission fitness, including preferential α4ß7 engagement. Using single genome amplification, we generated panels of both T/F (n = 20) and chronic (n = 20) Env constructs as well as full-length T/F (n = 6) and chronic (n = 4) infectious molecular clones (IMCs). We found that T/F and chronic control Envs were indistinguishable in the efficiency with which they used CD4 and CCR5. Both groups of Envs also exhibited the same CD4+ T cell subset tropism and showed similar sensitivity to neutralization by CD4 binding site (CD4bs) antibodies. Finally, saturating concentrations of anti-α4ß7 antibodies failed to inhibit infection and replication of T/F as well as chronic control viruses, although the growth of the tissue culture-adapted strain SF162 was modestly impaired. These results indicate that the population bottleneck associated with mucosal HIV-1 acquisition is not due to the selection of T/F viruses that use α4ß7, CD4 or CCR5 more efficiently.

Mucosal simian immunodeficiency virus transmission in African green monkeys: susceptibility to infection is proportional to target cell availability at mucosal sites.

African green monkeys (AGMs) are naturally infected with a simian immunodeficiency virus (SIVagm) that is nonpathogenic in its host. Although SIVagm is common and widespread, little is known about the mechanisms that govern its transmission. Since the earliest virus-host interactions may provide key insights into the nonpathogenic phenotype of SIVagm, we developed a mucosal transmission model for this virus. Using plasma from an acutely infected AGM as the virus inoculum, we exposed adult and juvenile AGMs, as well as pigtailed macaques (PTMs) as a nonnatural host control, by mucosal routes to increasing titers of virus and compared the doses needed to establish a productive infection. Four juvenile and four adult AGMs as well as two PTMs were intrarectally (IR) exposed, while two additional adult female AGMs were intravaginally (IVAG) exposed. No animal became infected following exposure to 10(5) RNA copies. Both PTMs but none of the AGMs became infected following exposure to 10(6) RNA copies. Finally, all adult AGMs and two of the four juvenile AGMs became infected following exposure to 10(7) RNA copies, acquiring either one (2 IR infected juveniles, 1 IR infected adult, 2 IVAG infected adults) or two (3 IR infected adults) transmitted founder viruses. These results were consistent with immunophenotypic data, which revealed a significant correlation between the percentage of CD4(+) T cells expressing CCR5 in the mucosa and the susceptibility to infection, in terms of both the viral dose and the numbers of transmitted founder viruses. Moreover, studies of uninfected AGMs showed that the fraction of CCR5-expressing CD4(+) T cells increased significantly with age. These results indicate that (i) AGMs are readily infected with SIVagm by both intrarectal and intravaginal routes, (ii) susceptibility to infection is proportional to the number of available CCR5(+) CD4(+) target cells in the mucosa, and (iii) the paucity of CCR5(+) CD4(+) target cells in infant and juvenile AGMs may explain the near absence of vertical transmission.

Bridging multiple dimensions: roles of transposable elements in higher-order genome regulation.

Transposable elements (TEs) such as endogenous retroviruses (ERVs), long interspersed nuclear elements (LINEs), and short interspersed nuclear elements (SINEs) occupy nearly half of typical mammalian genomes. Previous studies show that these parasitic elements, especially LINEs and ERVs, provide important activities promoting host germ cell and placental development, preimplantation embryogenesis, and maintenance of pluripotent stem cells. Despite being the most numerically abundant type of TEs in the genome, the consequences of SINEs on host genome regulation are less well characterized than those of ERVs and LINEs. Interestingly, recent findings reveal that SINEs recruit the key architectural protein CTCF (CCCTC-binding factor), indicating a role of these elements for 3D genome regulation. Higher-order nuclear structures are linked with important cellular functions such as gene regulation and DNA replication. SINEs and other TEs, therefore, may mediate distinct physiological processes with benefits to the host by modulating the 3D genome.

Prediction of the cell-type-specific transcription of non-coding RNAs from genome sequences via machine learning.

Gene transcription is regulated through complex mechanisms involving non-coding RNAs (ncRNAs). As the transcription of ncRNAs, especially of enhancer RNAs, is often low and cell type specific, how the levels of RNA transcription depend on genotype remains largely unexplored. Here we report the development and utility of a machine-learning model (MENTR) that reliably links genome sequence and ncRNA expression at the cell type level. Effects on ncRNA transcription predicted by the model were concordant with estimates from published studies in a cell-type-dependent manner, regardless of allele frequency and genetic linkage. Among 41,223 variants from genome-wide association studies, the model identified 7,775 enhancer RNAs and 3,548 long ncRNAs causally associated with complex traits across 348 major human primary cells and tissues, such as rare variants plausibly altering the transcription of enhancer RNAs to influence the risks of Crohn's disease and asthma. The model may aid the discovery of causal variants and the generation of testable hypotheses for biological mechanisms driving complex traits.

Mobile element variation contributes to population-specific genome diversification, gene regulation and disease risk.

Mobile genetic elements (MEs) are heritable mutagens that recursively generate structural variants (SVs). ME variants (MEVs) are difficult to genotype and integrate in statistical genetics, obscuring their impact on genome diversification and traits. We developed a tool that accurately genotypes MEVs using short-read whole-genome sequencing (WGS) and applied it to global human populations. We find unexpected population-specific MEV differences, including an Alu insertion distribution distinguishing Japanese from other populations. Integrating MEVs with expression quantitative trait loci (eQTL) maps shows that MEV classes regulate tissue-specific gene expression by shared mechanisms, including creating or attenuating enhancers and recruiting post-transcriptional regulators, supporting class-wide interpretability. MEVs more often associate with gene expression changes than SNVs, thus plausibly impacting traits. Performing genome-wide association study (GWAS) with MEVs pinpoints potential causes of disease risk, including a LINE-1 insertion associated with keloid and fasciitis. This work implicates MEVs as drivers of human divergence and disease risk.

Phenotypic and immunologic comparison of clade B transmitted/founder and chronic HIV-1 envelope glycoproteins.

Sexual transmission of human immunodeficiency virus type 1 (HIV-1) across mucosal barriers is responsible for the vast majority of new infections. This relatively inefficient process results in the transmission of a single transmitted/founder (T/F) virus, from a diverse viral swarm in the donor, in approximately 80% of cases. Here we compared the biological activities of 24 clade B T/F envelopes (Envs) with those from 17 chronic controls to determine whether the genetic bottleneck that occurs during transmission is linked to a particular Env phenotype. To maximize the likelihood of an intact mucosal barrier in the recipients and to enhance the sensitivity of detecting phenotypic differences, only T/F Envs from individuals infected with a single T/F variant were selected. Using pseudotyping to assess Env function in single-round infectivity assays, we compared coreceptor tropism, CCR5 utilization efficiencies, primary CD4(+) T cell subset tropism, dendritic cell trans-infections, fusion kinetics, and neutralization sensitivities. T/F and chronic Envs were phenotypically equivalent in most assays; however, T/F Envs were modestly more sensitive to CD4 binding site antibodies b12 and VRC01, as well as pooled human HIV Ig. This finding was independently validated with a panel of 14 additional chronic HIV-1 Env controls. Moreover, the enhanced neutralization sensitivity was associated with more efficient binding of b12 and VRC01 to T/F Env trimers. These data suggest that there are subtle but significant structural differences between T/F and chronic clade B Envs that may have implications for HIV-1 transmission and the design of effective vaccines.