Investigating mobile element variations by statistical genetics.

The integration of structural variations (SVs) in statistical genetics provides an opportunity to understand the genetic factors influencing complex human traits and disease. Recent advances in long-read technology and variant calling methods for short reads have improved the accurate discovery and genotyping of SVs, enabling their use in expression quantitative trait loci (eQTL) analysis and genome-wide association studies (GWAS). Mobile elements are DNA sequences that insert themselves into various genome locations. Insertional polymorphisms of mobile elements between humans, called mobile element variations (MEVs), contribute to approximately 25% of human SVs. We recently developed a variant caller that can accurately identify and genotype MEVs from biobank-scale short-read whole-genome sequencing (WGS) datasets and integrate them into statistical genetics. The use of MEVs in eQTL analysis and GWAS has a minimal impact on the discovery of genome loci associated with gene expression and disease; most disease-associated haplotypes can be identified by single nucleotide variations (SNVs). On the other hand, it helps make hypotheses about causal variants or effector variants. Focusing on MEVs, we identified multiple MEVs that contribute to differential gene expression and one of them is a potential cause of skin disease, emphasizing the importance of the integration of MEVs in medical genetics. Here, I will provide an overview of MEVs, MEV calling from WGS, and the integration of MEVs in statistical genetics. Finally, I will discuss the unanswered questions about MEVs, such as rare variants.

Phylogenetic analysis of the promoter element 2 of paramyxo- and filoviruses.

Paramyxo- and filovirus genomes are equipped with bipartite promoters at their 3' ends to initiate RNA synthesis. The two elements, the primary promoter element 1 (PE1) and the secondary promoter element 2 (PE2), are separated by a spacer region that must be precisely a multiple of 6 nucleotides (nts), indicating these viruses adhere to the "rule of six." However, our knowledge of PE2 has been limited to a narrow spectrum of virus species. In this study, a comparative analysis of 1,647 paramyxoviral genomes from a public database revealed that the paramyxovirus PE2 can be clearly categorized into two distinct subcategories: one marked by C repeats at every six bases (exclusive to the subfamily Orthoparamyxovirinae) and another characterized by CG repeats every 6 nts (observed in the subfamilies Avulavirinae and Rubulavirinae). This unique pattern collectively mirrors the evolutionary lineage of these subfamilies. Furthermore, we showed that PE2 of the Rubulavirinae, with the exception of mumps virus, serves as part of the gene-coding region. This may be due to the fact that the Rubulavirinae are the only paramyxoviruses that cannot propagate without RNA editing. Filoviruses have three to eight consecutive uracil repeats every six bases (UN5) in PE2, which is located in the 3' end region of the genome. We obtained PE2 sequences from 2,195 filoviruses in a public database and analyzed the sequence conservation among virus species. Our results indicate that the continuity of UN5 hexamers is consistently maintained with a high degree of conservation across virus species. IMPORTANCE: The genomic intricacies of paramyxo- and filoviruses are highlighted by the bipartite promoters-promoter element 1 (PE1) and promoter element 2 (PE2)-at their 3' termini. The spacer region between these elements follows the "rule of six," crucial for genome replication. By a comprehensive analysis of paramyxoviral genome sequences, we identified distinct subcategories of PE2 based on C and CG repeats that were specific to Orthoparamyxovirinae and Avulavirinae/Rubulavirinae, respectively, mirroring their evolutionary lineages. Notably, the PE2 of Rubulavirinae is integrated into the gene-coding region, a unique trait potentially linked to its strict dependence on RNA editing for virus growth. This study also focused on the PE2 sequences in filovirus genomes. The strict conservation of the continuity of UN5 among virus species emphasizes its crucial role in viral genome replication.

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.

Variant spectrum of von Hippel-Lindau disease and its genomic heterogeneity in Japan.

Von Hippel-Lindau (VHL) disease is an autosomal dominant, inherited syndrome with variants in the VHL gene, causing predisposition to multi-organ neoplasms with vessel abnormality. Germline variants in VHL can be detected in 80-90% of patients clinically diagnosed with VHL disease. Here, we summarize the results of genetic tests for 206 Japanese VHL families, and elucidate the molecular mechanisms of VHL disease, especially in variant-negative unsolved cases. Of the 206 families, genetic diagnosis was positive in 175 families (85%), including 134 families (65%) diagnosed by exon sequencing (15 novel variants) and 41 (20%) diagnosed by multiplex ligation-dependent probe amplification (MLPA) (one novel variant). The deleterious variants were significantly enriched in VHL disease Type 1. Interestingly, five synonymous or non-synonymous variants within exon 2 caused exon 2 skipping, which is the first report of exon 2 skipping caused by several missense variants. Whole genome and target deep sequencing analysis were performed for 22 unsolved cases with no variant identified and found three cases with VHL mosaicism (variant allele frequency: 2.5-22%), one with mobile element insertion in the VHL promoter region, and two with a pathogenic variant of BAP1 or SDHB. The variants associated with VHL disease are heterogeneous, and for more accuracy of the genetic diagnosis of VHL disease, comprehensive genome and DNA/RNA analyses are required to detect VHL mosaicism, complicated structure variants and other related gene variants.

A comprehensive list of the Bunyavirales replication promoters reveals a unique promoter structure in Nairoviridae differing from other virus families.

Members of the order Bunyavirales infect a wide variety of host species, including plants, animals and humans, and pose a threat to public health. Major families in this order have tri-segmented negative-sense RNA genomes, the 5' and 3' ends of which form complementary strands that serve as a replication promoter. Elucidation of the mechanisms by which viral polymerases recognize the promoter to initiate RNA synthesis is important for understanding viral replication and pathogenesis, and developing antivirals. A list of replication promoter configuration patterns may provide details on the differences in the replication mechanisms among bunyaviruses. By using public sequence data of all known bunyavirus species, we constructed a comprehensive list of the replication promoters comprising 40 nucleotides in both the 5' and 3' ends of the genome that form a specific complementary strand. Among tri-segmented bunyaviruses, members of the family Nairoviridae, including the highly pathogenic Crimean-Congo hemorrhagic fever virus, have evolved a GC-rich promoter structure differing from that of other families. The unique promoter structure might be related to the large genome size of the family Nairoviridae among bunyaviruses, and the large genome architecture might confer pathogenic advantages. The promoter list provided in this report is useful for predicting the virus family-specific replication mechanisms of bunyaviruses.

Wall teichoic acid-dependent phagocytosis of intact cell walls of Lactiplantibacillus plantarum elicits IL-12 secretion from macrophages.

Selected lactic acid bacteria can stimulate macrophages and dendritic cells to secrete IL-12, which plays a key role in activating innate and cellular immunity. In this study, we investigated the roles of cell wall teichoic acids (WTAs) displayed on whole intact cell walls (ICWs) of Lactiplantibacillus plantarum in activation of mouse macrophages. ICWs were prepared from whole bacterial cells of several lactobacilli without physical disruption, and thus retaining the overall shapes of the bacteria. WTA-displaying ICWs of several L. plantarum strains, but not WTA-lacking ICWs of strains of other lactobacilli, elicited IL-12 secretion from mouse bone marrow-derived macrophages (BMMs) and mouse macrophage-like J774.1 cells. The ability of the ICWs of L. plantarum to induce IL-12 secretion was abolished by selective chemical elimination of WTAs from ICWs, but was preserved by selective removal of cell wall glycopolymers other than WTAs. BMMs prepared from TLR2- or TLR4-deficient mouse could secret IL-12 upon stimulation with ICWs of L. plantarum and a MyD88 dimerization inhibitor did not affect ICW-mediated IL-12 secretion. WTA-displaying ICWs, but not WTA-lacking ICWs, were ingested in the cells within 30 min. Treatment with inhibitors of actin polymerization abolished IL-12 secretion in response to ICW stimulation and diminished ingestion of ICWs. When overall shapes of ICWs of L. plantarum were physically disrupted, the disrupted ICWs (DCWs) failed to induce IL-12 secretion. However, DCWs and soluble WTAs inhibited ICW-mediated IL-12 secretion from macrophages. Taken together, these results show that WTA-displaying ICWs of L. plantarum can elicit IL-12 production from macrophages via actin-dependent phagocytosis but TLR2 signaling axis independent pathway. WTAs displayed on ICWs are key molecules in the elicitation of IL-12 secretion, and the sizes and shapes of the ICWs have an impact on actin remodeling and subsequent IL-12 production.

Preparation of Non-overlapping Transposable Elements (TEs) Annotation by Interval Tree.

Transposable elements (TEs) are a major source of PIWI-interacting RNAs (piRNAs), therefore properly assigning piRNA library sequencing reads to the TEs from which they were derived is important for accurate assessment of piRNA biology. When calculating the abundance of small RNA-seq reads mapping to various TEs, a non-overlapping TE annotation is preferable because reads mapping to more than one genomic feature will often be excluded when counting reads. However, most unmodified TE annotations contain some degree of overlap between TE features. Here, I outline the principle and provide all scripts needed to resolve such overlapping regions of TE annotations to a single best TE annotation leveraging a computationally efficient tree algorithm. Non-overlapping annotations generated by this method can be directly used in commonly used read counting software.

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.

An endogenous bornavirus-like nucleoprotein in miniopterid bats retains the RNA-binding properties of the original viral protein.

Endogenous bornavirus-like nucleoprotein elements (EBLNs) are sequences derived from bornaviral N genes in vertebrate genomes. Some EBLNs have been suggested to encode functional proteins in host cells; however, little is known about their evolution and functional relationship to the viral genes from which EBLNs originate. Here, we predicted functionality of EBLNs based on the properties of N as an RNA-binding protein. We showed an EBLN in miniopterid bats (miEBLN-1) has evolved under purifying selection and encodes an RNA-binding protein (miEBLN-1p) with biochemical properties similar to bornaviral N. Furthermore, we revealed miEBLN-1p interacts with host RNA-binding proteins, such as MOV10. These data suggest that miEBLN-1p has been exapted as an RNA-binding protein with similar properties to exogenous bornaviral N in miniopterid bats.

Hidden Viral Sequences in Public Sequencing Data and Warning for Future Emerging Diseases.

RNA viruses cause numerous emerging diseases, mostly due to transmission from mammalian and avian reservoirs. Large-scale surveillance of RNA viral infections in these animals is a fundamental step for controlling viral infectious diseases. Metagenomic analysis is a powerful method for virus identification with low bias and has contributed substantially to the discovery of novel viruses. Deep-sequencing data have been collected from diverse animals and accumulated in public databases, which can be valuable resources for identifying unknown viral sequences. Here, we screened for infections of 33 RNA viral families in publicly available mammalian and avian sequencing data and found approximately 900 hidden viral infections. We also discovered six nearly complete viral genomes in livestock, wild, and experimental animals: hepatovirus in a goat, hepeviruses in blind mole-rats and a galago, astrovirus in macaque monkeys, parechovirus in a cow, and pegivirus in tree shrews. Some of these viruses were phylogenetically close to human-pathogenic viruses, suggesting the potential risk of causing disease in humans upon infection. Furthermore, infections of five novel viruses were identified in several different individuals, indicating that their infections may have already spread in the natural host population. Our findings demonstrate the reusability of public sequencing data for surveying viral infections and identifying novel viral sequences, presenting a warning about a new threat of viral infectious disease to public health. IMPORTANCE Monitoring the spread of viral infections and identifying novel viruses capable of infecting humans through animal reservoirs are necessary to control emerging viral diseases. Massive amounts of sequencing data collected from various animals are publicly available, and these data may contain sequences originating from a wide variety of viruses. Here, we analyzed more than 46,000 public sequencing data and identified approximately 900 hidden RNA viral infections in mammalian and avian samples. Some viruses discovered in this study were genetically similar to pathogens that cause hepatitis, diarrhea, or encephalitis in humans, suggesting the presence of new threats to public health. Our study demonstrates the effectiveness of reusing public sequencing data to identify known and unknown viral infections, indicating that future continuous monitoring of public sequencing data by metagenomic analyses would help prepare and mitigate future viral pandemics.

BUD23-TRMT112 interacts with the L protein of Borna disease virus and mediates the chromosomal tethering of viral ribonucleoproteins.

Persistent intranuclear infection is an uncommon infection strategy among RNA viruses. However, Borna disease virus 1 (BoDV-1), a nonsegmented, negative-strand RNA virus, maintains viral infection in the cell nucleus by forming structured aggregates of viral ribonucleoproteins (vRNPs), and by tethering these vRNPs onto the host chromosomes. To better understand the nuclear infection strategy of BoDV-1, we determined the host protein interactors of the BoDV-1 large (L) protein. By proximity-dependent biotinylation, we identified several nuclear host proteins interacting with BoDV-1 L, one of which is TRMT112, a partner of several methyltransferases (MTases). TRMT112 binds with BoDV-1 L at the RNA-dependent RNA polymerase domain, together with BUD23, an 18S ribosomal RNA MTase and 40S ribosomal maturation factor. We then discovered that BUD23-TRMT112 mediates the chromosomal tethering of BoDV-1 vRNPs, and that the MTase activity is necessary in the tethering process. These findings provide us a better understanding on how nuclear host proteins assist the chromosomal tethering of BoDV-1, as well as new prospects of host-viral interactions for intranuclear infection strategy of orthobornaviruses.

A Human Endogenous Bornavirus-Like Nucleoprotein Encodes a Mitochondrial Protein Associated with Cell Viability.

Endogenous retroviruses (ERVs) are sequences in animal genomes that originated from ancient retrovirus infections; they provide genetic novelty in hosts by being coopted as functional genes or elements during evolution. Recently, we demonstrated that endogenous elements from not only from retroviruses but also nonretroviral RNA viruses are a possible source of functional genes in host animals. The remnants of ancient bornavirus infections, called endogenous bornavirus-like elements (EBLs), are present in the genomes of a wide variety of vertebrate species, and some express functional products in host cells. Previous studies have predicted that the human EBL locus derived from bornavirus nucleoprotein, termed hsEBLN-2, expresses mRNA encoding a protein, suggesting that hsEBLN-2 has acquired a cellular function during evolution. However, the detailed function of the hsEBLN-2-derived product remains to be elucidated. In this study, we show that the hsEBLN-2-derived protein E2 acts as a mitochondrial protein that interacts with mitochondrial host factors associated with apoptosis, such as HAX-1. We also demonstrate that knockdown of hsEBLN-2-derived RNA increased the levels of PARP and caspase-3 cleavage and markedly decreased cell viability. In contrast, overexpression of E2 enhanced cell viability, as well as the intracellular stability of HAX-1, under stress conditions. Our results suggest that hsEBLN-2 has been coopted as a host gene, the product of which is involved in cell viability by interacting with mitochondrial proteins. IMPORTANCE Our genomes contain molecular fossils of ancient viruses, called endogenous virus elements (EVEs). Mounting evidence suggests that EVEs derived from nonretroviral RNA viruses have acquired functions in host cells during evolution. Previous studies have revealed that a locus encoding a bornavirus-derived EVE, hsEBLN-2, which was generated approximately 43 million years ago in a human ancestor, may be linked to the development of some tumors. However, the function of hsEBLN-2 has not been determined. In this study, we found that the E2 protein, an expression product of hsEBLN-2, interacts with apoptosis-related host proteins as a mitochondrial protein and affects cell viability. This study suggests that nonretroviral RNA viral EVEs have been coopted by hosts with more diverse functions than previously thought, showing a pivotal role for RNA virus infection in evolution.

100-My history of bornavirus infections hidden in vertebrate genomes.

Although viruses have threatened our ancestors for millions of years, prehistoric epidemics of viruses are largely unknown. Endogenous bornavirus-like elements (EBLs) are ancient bornavirus sequences derived from the viral messenger RNAs that were reverse transcribed and inserted into animal genomes, most likely by retrotransposons. These elements can be used as molecular fossil records to trace past bornaviral infections. In this study, we systematically identified EBLs in vertebrate genomes and revealed the history of bornavirus infections over nearly 100 My. We confirmed that ancient bornaviral infections have occurred in diverse vertebrate lineages, especially in primate ancestors. Phylogenetic analyses indicated that primate ancestors were infected with various bornaviral lineages during evolution. EBLs in primate genomes formed clades according to their integration ages, suggesting that bornavirus lineages infected with primate ancestors had changed chronologically. However, some bornaviral lineages may have coexisted with primate ancestors and underwent repeated endogenizations for tens of millions of years. Moreover, a bornaviral lineage that coexisted with primate ancestors also endogenized in the genomes of some ancestral bats. The habitats of these bat ancestors have been reported to overlap with the migration route of primate ancestors. These results suggest that long-term virus-host coexistence expanded the geographic distributions of the bornaviral lineage along with primate migration and may have spread their infections to these bat ancestors. Our findings provide insight into the history of bornavirus infections over geological timescales that cannot be deduced from research using extant viruses alone, thus broadening our perspective on virus-host coevolution.

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.

Identification of novel avian and mammalian deltaviruses provides new insights into deltavirus evolution.

Hepatitis delta virus (HDV) is a satellite virus that requires hepadnavirus envelope proteins for its transmission. Although recent studies identified HDV-related deltaviruses in certain animals, the evolution of deltaviruses, such as the origin of HDV and the mechanism of its coevolution with its helper viruses, is unknown, mainly because of the phylogenetic gaps among deltaviruses. Here, we identified novel deltaviruses of passerine birds, woodchucks, and white-tailed deer by extensive database searches and molecular surveillance. Phylogenetic and molecular epidemiological analyses suggest that HDV originated from mammalian deltaviruses and the past interspecies transmission of mammalian and passerine deltaviruses. Further, metaviromic and experimental analyses suggest that the satellite-helper relationship between HDV and hepadnavirus was established after the divergence of the HDV lineage from non-HDV mammalian deltaviruses. Our findings enhance our understanding of deltavirus evolution, diversity, and transmission, indicating the importance of further surveillance for deltaviruses.

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.

Effect of Spacer Length in Pyrene-Modified-Phenylboronic Acid Probe/CyD Complexes on Fluorescence-based Recognition of Monosaccharides in Aqueous Solution.

The chemical sensing of saccharides is of importance for the diagnosis of diabetes. Various enzymatic sensors have been developed, but their heat and pH instability issues need to be resolved. In this regard, the development of artificial saccharide sensors with high stability is attracting attention. We have designed a heat- and pH-stable supramolecular inclusion complex system composed of cyclodextrin (CyD) as a host and a phenylboronic acid (PB) probe possessing pyrene as a fluorescent guest. Several probes possessing alkyl spacers having various lengths between the PB and the pyrene moiety, Cn-APB (n = 1 - 4), were newly synthesized and evaluated with respect to their monosaccharide recognition ability on the basis of the fluorescence response through the cyclic esterification of monosaccharide and PB. These Cn-APB/CyD supramolecular inclusion complexes have exhibited a selective fluorescence response towards fructose in aqueous solution based on the photo-induced electron transfer mechanism. The spacer length of the alkyl group in Cn-APB significantly affects the affinity for saccharides. With respect to the complex between C4-APB and PB-modified CyD (3-PB-γ-CyD), it was found that the supramolecular inclusion complexes had high selectivity for glucose with significant fluorescence enhancement. These results indicate that the lengths of the alkyl spacers in the probe molecules are important to control the recognition of saccharides in aqueous solution.

Endogenous retroviruses drive species-specific germline transcriptomes in mammals.

Gene regulation in the germline ensures the production of high-quality gametes, long-term maintenance of the species and speciation. Male germline transcriptomes undergo dynamic changes after the mitosis-to-meiosis transition and have been subject to evolutionary divergence among mammals. However, the mechanisms underlying germline regulatory divergence remain undetermined. Here, we show that endogenous retroviruses (ERVs) influence species-specific germline transcriptomes. After the mitosis-to-meiosis transition in male mice, specific ERVs function as active enhancers to drive germline genes, including a mouse-specific gene set, and bear binding motifs for critical regulators of spermatogenesis, such as A-MYB. This raises the possibility that a genome-wide transposition of ERVs rewired germline gene expression in a species-specific manner. Of note, independently evolved ERVs are associated with the expression of human-specific germline genes, demonstrating the prevalence of ERV-driven mechanisms in mammals. Together, we propose that ERVs fine-tune species-specific transcriptomes in the mammalian germline.

ADAR2 Is Involved in Self and Nonself Recognition of Borna Disease Virus Genomic RNA in the Nucleus.

Cells sense pathogen-derived double-stranded RNA (dsRNA) as nonself. To avoid autoimmune activation by self dsRNA, cells utilize A-to-I editing by adenosine deaminase acting on RNA 1 (ADAR1) to disrupt dsRNA structures. Considering that viruses have evolved to exploit host machinery, A-to-I editing could benefit innate immune evasion by viruses. Borna disease virus (BoDV), a nuclear-replicating RNA virus, may require escape from nonself RNA-sensing and immune responses to establish persistent infection in the nucleus; however, the strategy by which BoDV evades nonself recognition is unclear. Here, we evaluated the involvement of ADARs in BoDV infection. The infection efficiency of BoDV was markedly decreased in both ADAR1 and ADAR2 knockdown cells at the early phase of infection. Microarray analysis using ADAR2 knockdown cells revealed that ADAR2 reduces immune responses even in the absence of infection. Knockdown of ADAR2 but not ADAR1 significantly reduced the spread and titer of BoDV in infected cells. Furthermore, ADAR2 knockout decreased the infection efficiency of BoDV, and overexpression of ADAR2 rescued the reduced infectivity in ADAR2 knockdown cells. However, the growth of influenza A virus, which causes acute infection in the nucleus, was not affected by ADAR2 knockdown. Moreover, ADAR2 bound to BoDV genomic RNA and induced A-to-G mutations in the genomes of persistently infected cells. We finally demonstrated that BoDV produced in ADAR2 knockdown cells induces stronger innate immune responses than those produced in wild-type cells. Taken together, our results suggest that BoDV utilizes ADAR2 to edit its genome to appear as "self" RNA in order to maintain persistent infection in the nucleus.IMPORTANCE Cells use the editing activity of adenosine deaminase acting on RNA proteins (ADARs) to prevent autoimmune responses induced by self dsRNA, but viruses can exploit this process to their advantage. Borna disease virus (BoDV), a nuclear-replicating RNA virus, must escape nonself RNA sensing by the host to establish persistent infection in the nucleus. We evaluated whether BoDV utilizes ADARs to prevent innate immune induction. ADAR2 plays a key role throughout the BoDV life cycle. ADAR2 knockdown reduced A-to-I editing of BoDV genomic RNA, leading to the induction of a strong innate immune response. These data suggest that BoDV exploits ADAR2 to edit nonself genomic RNA to appear as self RNA for innate immune evasion and establishment of persistent infection.