Search

Traumatic brain injury (TBI) is a leading cause of chronic brain impairment and results in a robust, but poorly understood, neuroinflammatory response that contributes to the long-term pathology. We used snRNA-seq to study transcriptomic changes in different cell populations in human brain tissue obtained acutely after severe, life-threatening TBI. This revealed a unique transcriptional response in oligodendrocyte precursors and mature oligodendrocytes, including the activation of a robust innate immune response, indicating an important role for oligodendroglia in the initiation of neuroinflammation. The activation of an innate immune response correlated with transcriptional upregulation of endogenous retroviruses in oligodendroglia. This observation was causally linkedin vitrousing human glial progenitors, implicating these ancient viral sequences in human neuroinflammation. In summary, this work provides a unique insight into the initiating events of the neuroinflammatory response in TBI, which has new therapeutic implications.

Retroviruses have left their legacy in host genomes over millions of years as endogenous retroviruses (ERVs), and their structure, diversity, and prevalence provide insights into the historical dynamics of retrovirus-host interactions. In bioinformatic analyses of koala (

A large portion of nuclear DNA is composed of transposable element (TE) sequences, whose transposition is controlled by diverse host defense strategies in order to maintain genomic integrity. One such strategy is the fungal-specific Repeat-Induced Point mutation (RIP) that hyper-mutates repetitive DNA sequences. While RIP is found across Fungi, it has been shown to vary in efficiency. The filamentous ascomycete Neurospora crassa has been a pioneer in the study of RIP, but data on TEs and RIP from other species in the genus is limited. In this study, we investigated 18 nearly gapless genome assemblies of ten Neurospora species, which diverged from a common ancestor about 7 MYA, to determine and compare genome-wide TE distribution and their associated RIP patterns. Four of these assemblies, generated by PacBio technology, represent new genomic datasets. We showed that the TE contents (8.7-18.9%) covary with genome sizes that range between 37.8-43.9 Mb. Degraded copies of Long Terminal Repeat (LTR) retrotransposons were abundant among the identified TEs, and these are distributed across the genome at varying frequencies. In all investigated Neurospora genomes, TE sequences had signs of numerous C-to-T substitutions, suggesting that RIP occurred in all species, and accordingly, RIP signatures correlated with TE-dense regions in all genomes. In conclusion, essentially gapless genome assemblies allowed us to identify TEs in Neurospora genomes, and reveal that TEs contribute to genome size variation in this group. Our study suggests that TEs and RIP are highly correlated in each examined Neurospora species, and hence, the pattern of interaction is conserved over the investigated evolutionary timescale. Finally, with our results, we verify that RIP signatures can be used to facilitate the identification of TE-rich regions in the genome. The comprehensive genomic dataset of Neurospora is a rich resource for further in- depth analyses of fungal genomes by the community.

Background Henny feathering in chickens is determined by a dominant mutation that transforms male-specific plumage to female-like plumage. Previous studies indicated that this phenotype is caused by ectopic expression in skin of CYP19A1 encoding aromatase that converts androgens to estrogen and thereby inhibits the development of male-specific plumage. A long terminal repeat (LTR) from an uncharacterized endogenous retrovirus (ERV) insertion was found in an isoform of the CYP19A1 transcript from henny feathering chicken. However, the complete sequence and the genomic position of the insertion were not determined. Results We used publicly available whole genome sequence data to determine the flanking sequences of the ERV, and then PCR amplified the entire insertion and sequenced it using Nanopore long reads and Sanger sequencing. The 7524 bp insertion contains an intact endogenous retrovirus that was not found in chickens representing 31 different breeds not showing henny feathering or in samples of the ancestral red junglefowl. The sequence shows over 99% sequence identity to the avian leukosis virus ev-1 and ev-21 strains, suggesting a recent integration. The ERV 3’LTR, containing a powerful transcriptional enhancer and core promoter with TATA box together with binding sites for EFIII and Ig/EBP inside the CYP19A1 5′ untranslated region, was detected partially in an aromatase transcript, which present a plausible explanation for ectopic expression of aromatase in non-ovarian tissues underlying the henny feathering phenotype. Conclusions We demonstrate that the henny feathering allele harbors an insertion of an intact avian leukosis virus at the 5’end of CYP19A1. The presence of this ERV showed complete concordance with the henny feathering phenotype both within a pedigree segregating for this phenotype and across breeds.

The underlying molecular mechanisms that determine long day versus short day breeders remain unknown in any organism. Atlantic herring provides a unique opportunity to examine the molecular mechanisms involved in reproduction timing, because both spring and autumn spawners exist within the same species. Although our previous whole genome comparisons revealed a strong association of TSHR alleles with spawning seasons, the functional consequences of these variants remain unknown. Here we examined the functional significance of six candidate TSHR mutations strongly associated with herring reproductive seasonality. We show that the L471M missense mutation in the spring-allele causes enhanced cAMP signaling. The best candidate non-coding mutation is a 5.2 kb retrotransposon insertion upstream of the TSHR transcription start site, near an open chromatin region, which is likely to affect TSHR expression. The insertion occurred prior to the split between Pacific and Atlantic herring and was lost in the autumn-allele. Our study shows that strongly associated coding and non-coding variants at the TSHR locus may both contribute to the regulation of seasonal reproduction in herring., This work was funded by Vetenskapsrådet (to L.A. and P.J.), the Knut and Alice Wallenberg foundation (to L.A.), Research Council of Norway project 254774, GENSINC (to L.A.) and Formas (to P.J.). Chunheng Mo was supported by the Overseas Exchange Program for Ph.D Students from China Scholarship Council (No. 201706240029). P.S. acknowledges funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through CRC 1423, project number 421152132, subprojects A01, A05, Z03; through CRC 1365, project number 394046635, subproject A03 (to G.K. and P.S.); through Germany´s Excellence Strategies—EXC 2008/1 (UniSysCat)—390540038. We thank the SNP&SEQ Technology Platform and UPPMAX in Uppsala for the provided service in high-throughput sequencing and computational infrastructure under the SNIC project 2020/15-137. These facilities are part of the National Genomics Infrastructure (NGI) of Sweden and Science for Life Laboratory. The SNP&SEQ Platform is also supported by the Swedish Research Council and the Knut and Alice Wallenberg Foundation. Open access funding provided by Uppsala University.

Endogenous retroviruses (ERVs) make up a large fraction of mammalian genomes and are thought to contribute to human disease, including brain disorders. In the brain, aberrant activation of ERVs is a potential trigger for an inflammatory response, but mechanistic insight into this phenomenon remains lacking. Using CRISPR/Cas9-based gene disruption of the epigenetic co-repressor protein Trim28, we found a dynamic H3K9me3-dependent regulation of ERVs in proliferating neural progenitor cells (NPCs), but not in adult neurons. In vivo deletion of Trim28 in cortical NPCs during mouse brain development resulted in viable offspring expressing high levels of ERVs in excitatory neurons in the adult brain. Neuronal ERV expression was linked to activated microglia and the presence of ERV-derived proteins in aggregate-like structures. This study demonstrates that brain development is a critical period for the silencing of ERVs and provides causal in vivo evidence demonstrating that transcriptional activation of ERV in neurons results in an inflammatory response.

Atlantic Halibut (Hippoglossus hippoglossus) has a X/Y genetic sex determination system, but the sex determining factor is not known. We produced a high-quality genome assembly and identified parts of chromosome 13 as the Y chromosome due to sequence divergence between sexes and segregation of sex genotypes in pedigrees. Linkage analysis revealed that all chromosomes exhibit heterochiasmy, i.e. male- and female restricted meiotic recombination intervals (MRR/FRR). We show that FRR/MRR intervals differ in nucleotide diversity and repeat class content and that this is true also for other Pleuronectidae species. We further show that remnants of a Gypsy-like transposable element insertion on chr13 promotes early male specific expression of gonadal somatic cell derived factor (gsdf). Less than 4 MYA, this male-determining element evolved on an autosomal FRR segment featuring pre-existing male meiotic recombination barriers, thereby creating a Y chromosome. We propose that heterochiasmy may facilitate the evolution of genetic sex determination systems.

BackgroundA large portion of nuclear DNA is composed of transposable element (TE) sequences, whose transposition is controlled by diverse host defense strategies in order to maintain genomic integrity. One such strategy is the fungal-specific Repeat-Induced Point (RIP) mutation that hyper-mutates repetitive DNA sequences. While RIP is found across Fungi, it has been shown to vary in efficiency. To date, detailed information on the TE landscapes and associated RIP patterns exist only in a few species belonging to highly divergent lineages.ResultWe investigated 18 nearly gapless genome assemblies of tenNeurosporaspecies, which diverged from a common ancestor about 7 MYA, to determine genome-wide TE distribution and their associated RIP patterns. We showed that the TE contents between 8.7-18.9% covary with genome sizes that range between 37.8-43.9 Mb. Degraded copies of Long Terminal Repeat (LTR) retrotransposons were abundant among the identified TEs, and these are distributed across the genome at varying frequencies. In all investigated genomes, TE sequences had signs of numerous C-to-T substitutions, suggesting that RIP occurred in all species. RIP signatures in all genomes correlated with TE-dense regions.ConclusionsEssentially gapless genome assemblies allowed us to identify TEs inNeurosporagenomes, and reveal that TEs contribute to genome size variation in this group. Our study suggests that TEs and RIP are highly correlated inNeurospora, and hence, the pattern of interaction is conserved over the investigated evolutionary timescale. We show that RIP signatures can be used to facilitate the identification of TE-rich region in the genome.

It is a broadly observed pattern that the non-recombining regions of sex-limited chromosomes (Y and W) accumulate more repeats than the rest of the genome, even in species like birds with a low genome-wide repeat content. Here, we show that in birds with highly heteromorphic sex chromosomes, the W chromosome has a transposable element (TE) density of greater than 55% compared to the genome-wide density of less than 10%, and contains over half of all full-length (thus potentially active) endogenous retroviruses (ERVs) of the entire genome. Using RNA-seq and protein mass spectrometry data, we were able to detect signatures of female-specific ERV expression. We hypothesize that the avian W chromosome acts as a refugium for active ERVs, probably leading to female-biased mutational load that may influence female physiology similar to the ‘toxic-Y’ effect in Drosophila males. Furthermore, Haldane's rule predicts that the heterogametic sex has reduced fertility in hybrids. We propose that the excess of W-linked active ERVs over the rest of the genome may be an additional explanatory variable for Haldane's rule, with consequences for genetic incompatibilities between species through TE/repressor mismatches in hybrids. Together, our results suggest that the sequence content of female-specific W chromosomes can have effects far beyond sex determination and gene dosage. This article is part of the theme issue ‘Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)’.

Endogenous retroviruses (ERVs) make up a large fraction of mammalian genomes and are thought to contribute to human disease, including brain disorders. In the brain, aberrant activation of ERVs is a potential trigger for neuroinflammation, but mechanistic insight into this phenomenon remains lacking. Using CRISPR/Cas9-based gene disruption of the epigenetic co-repressor protein Trim28, we found a dynamic H3K9me3-dependent regulation of ERVs in proliferating neural progenitor cells (NPCs), but not in adult neurons.In vivodeletion ofTrim28in cortical NPCs during mouse brain development resulted in viable offspring expressing high levels of ERV expression in excitatory neurons in the adult brain. Neuronal ERV expression was linked to inflammation, including activated microglia, and aggregates of ERV-derived proteins. This study demonstrates that brain development is a critical period for the silencing of ERVs and provides causalin vivoevidence demonstrating that transcriptional activation of ERV in neurons results in neuroinflammation.

More than half of the genomic landscape in humans and many other organisms is composed of repetitive DNA, which mostly derives from transposable elements (TEs) and viruses. Recent technological advances permit improved assessment of the repetitive content across genomes and newly developed molecular assays have revealed important roles of TEs and viruses in host genome evolution and organization. To update on our current understanding of TE biology and to promote new interdisciplinary strategies for the TE research community, leading experts gathered for the 2nd Uppsala Transposon Symposium on October 4–5, 2018 in Uppsala, Sweden. Using cutting-edge single-molecule and single-cell approaches, research on TEs and other repeats has entered a new era in biological and biomedical research.

SUMMARY TRIM28 is a corepressor that mediates transcriptional silencing by establishing local heterochromatin. Here, we show that deletion of TRIM28 in neural progenitor cells (NPCs) results in high-level expression of two groups of endogenous retroviruses (ERVs): IAP1 and MMERVK10C. We find that NPCs use TRIM28-mediated histone modifications to dynamically regulate transcription and silencing of ERVs, which is in contrast to other somatic cell types using DNA methylation. We also show that derepression of ERVs influences transcriptional dynamics in NPCs through the activation of nearby genes and the expression of long noncoding RNAs. These findings demonstrate a unique dynamic transcriptional regulation of ERVs in NPCs. Our results warrant future studies on the role of ERVs in the healthy and diseased brain., Graphical Abstract

Although extensive research has demonstrated host-retrovirus microevolutionary dynamics, it has been difficult to gain a deeper understanding of the macroevolutionary patterns of host-retrovirus interactions. Here we use recent technological advances to infer broad patterns in retroviral diversity, evolution, and host-virus relationships by using a large-scale phylogenomic approach using endogenous retroviruses (ERVs). Retroviruses insert a proviral DNA copy into the host cell genome to produce new viruses. ERVs are provirus insertions in germline cells that are inherited down the host lineage and consequently present a record of past host-viral associations. By mining ERVs from 65 host genomes sampled across vertebrate diversity, we uncover a great diversity of ERVs, indicating that retroviral sequences are much more prevalent and widespread across vertebrates than previously appreciated. The majority of ERV clades that we recover do not contain known retroviruses, implying either that retroviral lineages are highly transient over evolutionary time or that a considerable number of retroviruses remain to be identified. By characterizing the distribution of ERVs, we show that no major vertebrate lineage has escaped retroviral activity and that retroviruses are extreme host generalists, having an unprecedented ability for rampant host switching among distantly related vertebrates. In addition, we examine whether the distribution of ERVs can be explained by host factors predicted to influence viral transmission and find that internal fertilization has a pronounced effect on retroviral colonization of host genomes. By capturing the mode and pattern of retroviral evolution and contrasting ERV diversity with known retroviral diversity, our study provides a cohesive framework to understand host-virus coevolution better.

Genomic data provide an excellent resource to improve understanding of retrovirus evolution and the complex relationships among viruses and their hosts. In conjunction with broad-scale in silico screening of vertebrate genomes, this resource offers an opportunity to complement data on the evolution and frequency of past retroviral spread and so evaluate future risks and limitations for horizontal transmission between different host species. Here, we develop a methodology for extracting phylogenetic signal from large endogenous retrovirus (ERV) datasets by collapsing information to facilitate broad-scale phylogenomics across a wide sample of hosts. Starting with nearly 90,000 ERVs from 60 vertebrate host genomes, we construct phylogenetic hypotheses and draw inferences regarding the designation, host distribution, origin, and transmission of the Gammaretrovirus genus and associated class I ERVs. Our results uncover remarkable depths in retroviral sequence diversity, supported within a phylogenetic context. This finding suggests that current infectious exogenous retrovirus diversity may be underestimated, adding credence to the possibility that many additional exogenous retroviruses may remain to be discovered in vertebrate taxa. We demonstrate a history of frequent horizontal interorder transmissions from a rodent reservoir and suggest that rats may have acted as important overlooked facilitators of gammaretrovirus spread across diverse mammalian hosts. Together, these results demonstrate the promise of the methodology used here to analyze large ERV datasets and improve understanding of retroviral evolution and diversity for utilization in wider applications.

Endogenous retroviruses (ERVs), which make up 8% of the human genome, have been proposed to participate in the control of gene regulatory networks. In this study, we find a region- and developmental stage-specific expression pattern of ERVs in the developing human brain, which is linked to a transcriptional network based on ERVs. We demonstrate that almost 10,000, primarily primate-specific, ERVs act as docking platforms for the co-repressor protein TRIM28 in human neural progenitor cells, which results in the establishment of local heterochromatin. Thereby, TRIM28 represses ERVs and consequently regulates the expression of neighboring genes. These results uncover a gene regulatory network based on ERVs that participates in control of gene expression of protein-coding transcripts important for brain development.

Ecological adaptation is of major relevance to speciation and sustainable population management, but the underlying genetic factors are typically hard to study in natural populations due to genetic differentiation caused by natural selection being confounded with genetic drift in subdivided populations. Here, we use whole genome population sequencing of Atlantic and Baltic herring to reveal the underlying genetic architecture at an unprecedented detailed resolution for both adaptation to a new niche environment and timing of reproduction. We identify almost 500 independent loci associated with a recent niche expansion from marine (Atlantic Ocean) to brackish waters (Baltic Sea), and more than 100 independent loci showing genetic differentiation between spring- and autumn-spawning populations irrespective of geographic origin. Our results show that both coding and non-coding changes contribute to adaptation. Haplotype blocks, often spanning multiple genes and maintained by selection, are associated with genetic differentiation. DOI: http://dx.doi.org/10.7554/eLife.12081.001, eLife digest The Atlantic herring is one of the most common fish in the world and has been a crucial food resource in northern Europe. One school of herring may comprise billions of fish, but previous studies had only revealed very few genetic differences in herring from different geographic regions. This was unexpected since Atlantic herring is one of the few marine species that can reproduce throughout the brackish Baltic Sea, which can be about a tenth as salty as the Atlantic Ocean. This unexpected finding could be explained in at least two different ways. Firstly, perhaps Atlantic herring are flexible enough to adapt to very different environments (i.e. high or low salinity) without much genetic change. Secondly, the previous studies only looked at a handful of sites in the Atlantic herring’s genome and so it is possible that genetic differences at other genes control this fish’s adaptation instead. Now, Martinez Barrio, Lamichhaney, Fan, Rafati et al. have sequenced entire genomes from groups of Atlantic herring and revealed hundreds of sites that are associated with adaptation to the Baltic Sea. The analysis also identified a number of genes that control when these fish reproduce by comparing herring that spawn in the autumn with those that spawn in spring. This is important because natural populations must carefully time when they reproduce to maximize the survival of their young. These new findings provide compelling evidence that changes in protein-coding genes and stretches of DNA that regulate the expression of other genes both contribute to adaptation in herrings. The analysis also clearly shows that variants of genes that contribute to adaptation were likely to evolve over time by accumulating multiple sequence changes affecting the same gene. Furthermore, these gene variants essentially form a rich “tool-box” that underlies the Atlantic herring’s adaptation to its environment, and different subpopulations of herring were found to have their own optimal sets of gene variants. For instance, autumn-spawning herring and spring-spawning herring from the Baltic Sea both have gene variants that favor adaptation to low salinity. However, autumn-spawning Baltic herring also share gene variants that favor spawning in the autumn with autumn-spawning herring from the North Sea, but not with spring-spawning Baltic herring. The next step will be to study how the 500 or so genes identified affect adaptation at the molecular level. This will likely involve experiments with other model fish such as zebrafish and sticklebacks. Finally, these new findings can be directly applied to monitor stocks of herring to make herring fisheries more sustainable. DOI: http://dx.doi.org/10.7554/eLife.12081.002

For 10,000 years pigs and humans have shared a close and complex relationship. From domestication to modern breeding practices, humans have shaped the genomes of domestic pigs. Here we present the assembly and analysis of the genome sequence of a female domestic Duroc pig (Sus scrofa) and a comparison with the genomes of wild and domestic pigs from Europe and Asia. Wild pigs emerged in South East Asia and subsequently spread across Eurasia. Our results reveal a deep phylogenetic split between European and Asian wild boars ∼1 million years ago, and a selective sweep analysis indicates selection on genes involved in RNA processing and regulation. Genes associated with immune response and olfaction exhibit fast evolution. Pigs have the largest repertoire of functional olfactory receptor genes, reflecting the importance of smell in this scavenging animal. The pig genome sequence provides an important resource for further improvements of this important livestock species, and our identification of many putative disease-causing variants extends the potential of the pig as a biomedical model., The research leading to these results has received funding from the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement no. 249894 (SelSweep)

Retroviruses are the only group of viruses known to have left a fossil record, in the form of endogenous proviruses, and approximately 8% of the human genome is made up of these elements. Although many other viruses, including non-retroviral RNA viruses, are known to generate DNA forms of their own genomes during replication, none has been found as DNA in the germline of animals. Bornaviruses, a genus of non-segmented, negative-sense RNA virus, are unique among RNA viruses in that they establish persistent infection in the cell nucleus. Here we show that elements homologous to the nucleoprotein (N) gene of bornavirus exist in the genomes of several mammalian species, including humans, non-human primates, rodents and elephants. These sequences have been designated endogenous Borna-like N (EBLN) elements. Some of the primate EBLNs contain an intact open reading frame (ORF) and are expressed as mRNA. Phylogenetic analyses showed that EBLNs seem to have been generated by different insertional events in each specific animal family. Furthermore, the EBLN of a ground squirrel was formed by a recent integration event, whereas those in primates must have been formed more than 40 million years ago. We also show that the N mRNA of a current mammalian bornavirus, Borna disease virus (BDV), can form EBLN-like elements in the genomes of persistently infected cultured cells. Our results provide the first evidence for endogenization of non-retroviral virus-derived elements in mammalian genomes and give novel insights not only into generation of endogenous elements, but also into a role of bornavirus as a source of genetic novelty in its host.

For millions of years, retroviral infections have challenged vertebrates, occasionally leading to germline integration and inheritance as ERVs, genetic parasites whose remnants today constitute some 7% to 8% of the human genome. Although they have had significant evolutionary side effects, it is useful to view ERVs as fossil representatives of retroviruses extant at the time of their insertion into the germline and not as direct players in the evolutionary process itself. Expression of particular ERVs is associated with several positive physiological functions as well as certain diseases, although their roles in human disease as etiological agents, possible contributing factors, or disease markers—well demonstrated in animal models—remain to be established. Here we discuss ERV contributions to host genome structure and function, including their ability to mediate recombination, and physiological effects on the host transcriptome resulting from their integration, expression, and other events.

Eukaryotic genomes contain many endogenous retroviral sequences (ERVs). ERVs are often severely mutated, therefore difficult to detect. A platform independent (Java) program package, RetroTector© (ReTe), was constructed. It has three basic modules: (i) detection of candidate long terminal repeats (LTRs), (ii) detection of chains of conserved retroviral motifs fulfilling distance constraints and (iii) attempted reconstruction of original retroviral protein sequences, combining alignment, codon statistics and properties of protein ends. Other features are prediction of additional open reading frames, automated database collection, graphical presentation and automatic classification. ReTe favors elements >1000-bp long due to its dependence on order of and distances between retroviral fragments. It detects single or low-copy-number elements. ReTe assigned a ‘retroviral’ score of 890–2827 to 10 exogenous retroviruses from seven genera, and accurately predicted their genes. In a simulated model, ReTe was robust against mutational decay. The human genome was analyzed in 1–2 days on a LINUX cluster. Retroviral sequences were detected in divergent vertebrate genomes. Most ReTe detected chains were coincident with Repeatmasker output and the HERVd database. ReTe did not report most of the evolutionary old HERV-L related and MalR sequences, and is not yet tailored for single LTR detection. Nevertheless, ReTe rationally detects and annotates many retroviral sequences.

After performing de novo transcript assembly of >1 billion RNA-Sequencing reads obtained from 22 samples of different Norway spruce (Picea abies) tissues that were not surface sterilized, we found that assembled sequences captured a mix of plant, lichen, and fungal transcripts. The latter were likely expressed by endophytic and epiphytic symbionts, indicating that these organisms were present, alive, and metabolically active. Here, we show that these serendipitously sequenced transcripts need not be considered merely as contamination, as is common, but that they provide insight into the plant's phyllosphere. Notably, we could classify these transcripts as originating predominantly from Dothideomycetes and Leotiomycetes species, with functional annotation of gene families indicating active growth and metabolism, with particular regards to glucose intake and processing, as well as gene regulation.

The human genome is littered by endogenous retrovirus sequences (HERVs), which constitute up to 8% of the total genomic sequence. The sequencing of the human (Homo sapiens) and chimpanzee (Pan troglodytes) genomes has facilitated the evolutionary study of ERVs and related sequences. We screened both the human genome (version hg16) and the chimpanzee genome (version PanTro1) for ERVs and conducted a phylogenetic analysis of recent integrations. We found a number of recent integrations within both genomes. They segregated into four groups. Two larger gammaretrovirus-like groups (PtG1 and PtG2) occurred in chimpanzees but not in humans. The PtG sequences were most similar to two baboon ERVs and a macaque sequence but neither to other chimpanzee ERVs nor to any human gammaretrovirus-like ERVs. The pattern was consistent with cross-species transfer via predation. This appears to be an example of horizontal transfer of retroviruses with occasional fixation in the germ line. Human endogenous retrovirus sequences (HERVs), which have become trapped as Mendelian genes and fixed in the germ line (9), constitute 7 to 8% of the human genome (8, 26). The sequencing of the human (Homo sapiens) genome (26) and the draft sequence of the chimpanzee (Pan troglodytes) genome (16) have facilitated the evolutionary studies of ERVs and related sequences. The Homo-Pan speciation has been estimated to have occurred between 4.6 and 6.2 million years ago (15). The genomic nonidentity is approximately 1.2% (15, 17) up to 5% (14), depending on analysis methods. We used a newly developed bioinformatic tool, RetroTector (G. O. Sperber, T. Airola, P. Jern, A. Castell, and J. Blomberg, unpublished data), to screen both the human genome (version hg16) and the chimpanzee genome (version PanTro1) for ERVs. We have successfully used RetroTector-derived data in several papers; the program validation can be traced in these papers and their respective supplemental materials (1, 27–29, 42). The differences in ERV content serve as indicators of dynamic genomes. The proviral long terminal repeats (LTRs) are identical at the integration moment (50). We analyzed recent provirus integrations with less than 2% LTR difference. With a typical value of 0.2% substitutions per million years (34), this LTR difference served as an approximation to 5 million years and the Homo-Pan sp. split. In this paper, we refer to these proviruses as “recent,” knowing that postintegrational mutations can accumulate at different rates depending on the genetic environment (25, 33, 60). The potential effects of retroviral integrations have been studied in several papers. The effects include addition of promoters and enhancers (43) and

Endogenous retroviral sequences (ERVs) are dynamic genomic components with profound influences on gene expression and genomic structure. Their extent of expression is not well known. Several broadly targeted real-time reverse transcription PCR (QPCRs) systems for surveillance of RNA expression of the major groups of human gammaretroviral ERVs were constructed. The highly conserved reverse transcriptase (RT) and integrase (IN) domains of the pol gene were used as targets for the PCRs, which were both probe-based (TaqMan) and probe-less (SYBR Green). Different levels of primer and probe degeneracy, with or without inosine, were tested. Several of the PCRs had sensitivities of a few HERV nucleic acid copies per PCR reaction. Specificities were approximately as expected from the fit of primers and probes. Gammaretroviral HERV RNA expression was studied in different human tissues. Each HERV group had a specific pattern of expression. HERV-E was highly expressed in testis, HERV-I/T in brain and testis, HERV-H in brain and testis, while HERV-W was highly expressed in placenta. Gammaretroviral RNA was not detected in plasma from 50 blood donors in saliva from 20 persons. In conclusion, a set of tools for investigation of gammaretroviral HERV RNA expression was created.

We defined the abundant human endogenous retrovirus group HERV-H based on pol similarity. Among 3661 pol-containing elements, 1124 integrations were similar to HERV-H RGH2 pol using translated pol sequences. A clustering procedure lessened these to 234 representatives, amenable to detailed study. Among the 1124, 926 clustered into HERV-H and 106 into adjacent HERV-H-like, the remainder being more distant to HERV-H. The HERV-H group was divided into RTVLH2-like (705 elements) and RGH2-like (77 elements) subgroups. Among 926 HERV-H, LTR differences were 1–33%, 10% had env, 78% had gag, 66% had a histidine primer binding site (PBS), and 3% (both subgroups) had a phenylalanine PBS. Allelic differences in env were studied using a convenient temperature gradient gel electrophoresis (TGGE) method and a genomic single nucleotide polymorphism (SNP) search. A pattern of abundant defective elements and less abundant less defective ones led us to formulate a “midwife” master model where more complete elements help the others in trans to transpose.

Summary: High-throughput genotyping and sequencing technologies facilitate studies of complex genetic traits and provide new research opportunities. The increasing popularity of genome-wide association studies (GWAS) leads to the discovery of new associated loci and a better understanding of the genetic architecture underlying not only diseases, but also other monogenic and complex phenotypes. Several softwares are available for performing GWAS analyses, R environment being one of them. Results: We present cgmisc, an R package that enables enhanced data analysis and visualization of results from GWAS. The package contains several utilities and modules that complement and enhance the functionality of the existing software. It also provides several tools for advanced visualization of genomic data and utilizes the power of the R language to aid in preparation of publication-quality figures. Some of the package functions are specific for the domestic dog (Canis familiaris) data. Availability and implementation: The package is operating system-independent and is available from: https://github.com/cgmisc-team/cgmisc Contact: marcin.kierczak@imbim.uu.se Supplementary information: Supplementary data are available at Bioinformatics online.

Background Phenomena such as incomplete lineage sorting, horizontal gene transfer, gene duplication and subsequent sub- and neo-functionalisation can result in distinct local phylogenetic relationships that are discordant with species phylogeny. In order to assess the possible biological roles for these subdivisions, they must first be identified and characterised, preferably on a large scale and in an automated fashion. Results We developed Saguaro, a combination of a Hidden Markov Model (HMM) and a Self Organising Map (SOM), to characterise local phylogenetic relationships among aligned sequences using cacti, matrices of pair-wise distance measures. While the HMM determines the genomic boundaries from aligned sequences, the SOM hypothesises new cacti in an unsupervised and iterative fashion based on the regions that were modelled least well by existing cacti. After testing the software on simulated data, we demonstrate the utility of Saguaro by testing two different data sets: (i) 181 Dengue virus strains, and (ii) 5 primate genomes. Saguaro identifies regions under lineage-specific constraint for the first set, and genomic segments that we attribute to incomplete lineage sorting in the second dataset. Intriguingly for the primate data, Saguaro also classified an additional ~3% of the genome as most incompatible with the expected species phylogeny. A substantial fraction of these regions was found to overlap genes associated with both the innate and adaptive immune systems. Conclusions Saguaro detects distinct cacti describing local phylogenetic relationships without requiring any a priori hypotheses. We have successfully demonstrated Saguaro’s utility with two contrasting data sets, one containing many members with short sequences (Dengue viral strains: n = 181, genome size = 10,700 nt), and the other with few members but complex genomes (related primate species: n = 5, genome size = 3 Gb), suggesting that the software is applicable to a wide variety of experimental populations. Saguaro is written in C++, runs on the Linux operating system, and can be downloaded from http://saguarogw.sourceforge.net/.

ZBED genes originate from domesticated hAT DNA transposons and encode regulatory proteins of diverse function in vertebrates. Here we reveal the evolutionary relationship between ZBED genes and demonstrate that they are derived from at least two independent domestication events in jawed vertebrate ancestors. We show that ZBEDs form two monophyletic clades, one of which has expanded through several independent duplications in host lineages. Subsequent diversification of ZBED genes has facilitated regulation of multiple diverse fundamental functions. In contrast to known examples of transposable element exaptation, our results demonstrate a novel unprecedented capacity for the repeated utilization of a family of transposable element-derived protein domains sequestered as regulators during the evolution of diverse host gene functions in vertebrates. Specifically, ZBEDs have contributed to vertebrate regulatory innovation through the donation of modular DNA and protein interacting domains. We identify that C7ORF29, ZBED2, 3, 4, and ZBEDX form a monophyletic group together with ZBED6, that is distinct from ZBED1 genes. Furthermore, we show that ZBED5 is related to Buster DNA transposons and is phylogenetically separate from other ZBEDs. Our results offer new insights into the evolution of regulatory pathways, and suggest that DNA transposons have contributed to regulatory complexity during genome evolution in vertebrates.

Domestication of wild boar ( Sus scrofa ) and subsequent selection have resulted in dramatic phenotypic changes in domestic pigs for a number of traits, including behavior, body composition, reproduction, and coat color. Here we have used whole-genome resequencing to reveal some of the loci that underlie phenotypic evolution in European domestic pigs. Selective sweep analyses revealed strong signatures of selection at three loci harboring quantitative trait loci that explain a considerable part of one of the most characteristic morphological changes in the domestic pig—the elongation of the back and an increased number of vertebrae. The three loci were associated with the NR6A1, PLAG1 , and LCORL genes. The latter two have repeatedly been associated with loci controlling stature in other domestic animals and in humans. Most European domestic pigs are homozygous for the same haplotype at these three loci. We found an excess of derived nonsynonymous substitutions in domestic pigs, most likely reflecting both positive selection and relaxed purifying selection after domestication. Our analysis of structural variation revealed four duplications at the KIT locus that were exclusively present in white or white-spotted pigs, carrying the Dominant white , Patch , or Belt alleles. This discovery illustrates how structural changes have contributed to rapid phenotypic evolution in domestic animals and how alleles in domestic animals may evolve by the accumulation of multiple causative mutations as a response to strong directional selection.

APOBEC3 proteins function to restrict the replication of retroviruses. One mechanism of this restriction is deamination of cytidines to uridines in (-) strand DNA, resulting in hypermutation of guanosines to adenosines in viral (+) strands. However, Moloney murine leukemia virus (MoMLV) is partially resistant to restriction by mouse APOBEC3 (mA3) and virtually completely resistant to mA3-induced hypermutation. In contrast, the sequences of MLV genomes that are in mouse DNA suggest that they were susceptible to mA3-induced deamination when they infected the mouse germline. We tested the possibility that sensitivity to mA3 restriction and to deamination resides in the viral gag gene. We generated a chimeric MLV in which the gag gene was from an endogenous MLV in the mouse germline, while the remainder of the viral genome was from MoMLV. This chimera was fully infectious but its response to mA3 was indistinguishable from that of MoMLV. Thus, the Gag protein does not seem to control the sensitivity of MLVs to mA3. We also found that MLVs inactivated by mA3 do not synthesize viral DNA upon infection; thus mA3 restriction of MLV occurs before or at reverse transcription. In contrast, HIV-1 restricted by mA3 and MLVs restricted by human APOBEC3G do synthesize DNA; these DNAs exhibit APOBEC3-induced hypermutation.

Host-retrovirus interactions influence the genomic landscape and have contributed substantially to mammalian genome evolution. To gain further insights, we analyzed a female boxer (Canis familiaris) genome for complexity and integration pattern of canine endogenous retroviruses (CfERV). Intriguingly, the first such in-depth analysis of a carnivore species identified 407 CfERV proviruses that represent only 0.15% of the dog genome. In comparison, the same detection criteria identified about six times more HERV proviruses in the human genome that has been estimated to contain a total of 8% retroviral DNA including solitary LTRs. These observed differences in man and dog are likely due to different mechanisms to purge, restrict and protect their genomes against retroviruses. A novel group of gammaretrovirus-like CfERV with high similarity to HERV-Fc1 was found to have potential for active retrotransposition and possibly lateral transmissions between dog and human as a result of close interactions during at least 10.000 years. The CfERV integration landscape showed a non-uniform intra- and inter-chromosomal distribution. Like in other species, different densities of ERVs were observed. Some chromosomal regions were essentially devoid of CfERVs whereas other regions had large numbers of integrations in agreement with distinct selective pressures at different loci. Most CfERVs were integrated in antisense orientation within 100 kb from annotated protein-coding genes. This integration pattern provides evidence for selection against CfERVs in sense orientation relative to chromosomal genes. In conclusion, this ERV analysis of the first carnivorous species supports the notion that different mammals interact distinctively with endogenous retroviruses and suggests that retroviral lateral transmissions between dog and human may have occurred.