Multiple effects govern endogenous retrovirus survival patterns in human gene introns

Trans-species polymorphism, HLA-disease associations and the evolution of the MHC

Currently, the paradigm is that major histocompatibility complex (MHC) polymorphism is maintained by balancing selection on the immune genes. However, other evolutionary forces besides selection also play a role in the population genetics of this multigene family. van Oosterhout proposed a new theory of MHC evolution called associative balancing complex (ABC) evolution.1 This theory incorporates the effects of the evolutionary forces in the entire MHC region (peri-MHC), and it proposes that recessive deleterious mutations can accumulate in the peri-MHC in a process similar to Muller's ratchet.2 These mutations are not easily purged because epistasis and high gene diversity in the MHC reduce the efficacy of natural selection. Because natural selection is less efficient, it could also make the MHC prone to the onslaught of genomic parasites such as retroviruses and transposable elements (TEs). The accumulated genetic load has important consequences for the evolution of this immune gene family, and it can reinforce linkage disequilibria and help to maintain the MHC polymorphism. ABC evolution offers new insights into some of the most puzzling aspects of the MHC, including the occurrence of identical MHC sequences in diverged species (i.e., trans-species polymorphism). It may also explain why the large numbers of disease-associated mutations are not removed by natural selection, and why the genes that protect vertebrates against infectious diseases are associated to such a wide variety of genetic disorders.

Large-scale transcriptome data reveals transcriptional activity of fission yeast LTR retrotransposons

Background

Retrotransposons are transposable elements that proliferate within eukaryotic genomes through a process involving reverse transcription. The numbers of retrotransposons within genomes and differences between closely related species may yield insight into the evolutionary history of the elements. Less is known about the ongoing dynamics of retrotransposons, as analysis of genome sequences will only reveal insertions of retrotransposons that are fixed - or near fixation - in the population or strain from which genetic material has been extracted for sequencing. One pre-requisite for retrotransposition is transcription of the elements. Given their intrinsic sequence redundancy, transcriptome-level analyses of transposable elements are scarce. We have used recently published transcriptome data from the fission yeast Schizosaccharomyces pombe to assess the ability to detect and describe transcriptional activity from Long Terminal Repeat (LTR) retrotransposons. LTR retrotransposons are normally flanked by two LTR sequences. However, the majority of LTR sequences in S. pombe exist as solitary LTRs, i.e. as single terminal repeat sequences not flanking a retrotransposon. Transcriptional activity was analysed for both full-length LTR retrotransposons and solitary LTRs.

Results

Two independent sets of transcriptome data reveal the presence of full-length, polyadenylated transcripts from LTR retrotransposons in S. pombe during growth phase in rich medium. The redundancy of retrotransposon sequences makes it difficult to assess which elements are transcriptionally active, but data strongly indicates that only a subset of the LTR retrotransposons contribute significantly to the detected transcription. A considerable level of reverse strand transcription is also detected. Equal levels of transcriptional activity are observed from both strands of solitary LTR sequences. Transcriptome data collected during meiosis suggests that transcription of solitary LTRs is correlated with the transcription of nearby protein-coding genes.

Conclusions

Presumably, the host organism negatively regulates proliferation of LTR retrotransposons. The finding of considerable transcriptional activity of retrotransposons suggests that part of this regulation is likely to take place at a post-transcriptional level. Alternatively, the transcriptional activity may signify a hitherto unrecognized activity level of retrotransposon proliferation. Our findings underline the usefulness of transcriptome data in elucidating dynamics in retrotransposon transcription.

Extensive DRB region diversity in cynomolgus macaques: recombination as a driving force

The DR region of primate species is generally complex and displays diversity concerning the number and combination of distinct types of DRB genes present per region configuration. A highly variable short tandem repeat (STR) present in intron 2 of nearly all primate DRB genes can be utilized as a quick and accurate high through-put typing procedure. This approach resulted previously in the description of unique and haplotype-specific DRB-STR length patterns in humans, chimpanzees, and rhesus macaques. For the present study, a cohort of 230 cynomolgus monkeys, including self-sustaining breeding groups, has been examined. MtDNA analysis showed that most animals originated from the Indonesian islands, but some are derived from the mainland, south and north of the Isthmus of Kra. Haplotyping and subsequent sequencing resulted in the detection of 118 alleles, including 28 unreported ones. A total of 49 Mafa-DRB region configurations were detected, of which 28 have not yet been described. Humans and chimpanzees possess a low number of different DRB region configurations in concert with a high degree of allelic variation. In contrast, however, allelic heterogeneity within a given Mafa-DRB configuration is even less frequently observed than in rhesus macaques. Several of these region configurations appear to have been generated by recombination-like events, most probably propagated by a retroviral element mapping within DRB6 pseudogenes, which are present on the majority of haplotypes. This undocumented high level of DRB region configuration-associated diversity most likely represents a species-specific strategy to cope with various pathogens.

Intracisternal A particle genes: Distribution in the mouse genome, active subtypes, and potential roles as species-specific mediators of susceptibility to cancer

Rodents, mice and rats in particular, are the species of choice for evaluating chemical carcinogenesis. However, different species and strains often respond very differently, undermining the logic of extrapolation of animal results to humans and complicating risk assessment. Intracisternal A particles (IAPs), endogenous retroviral sequences, are an important class of transposable elements that induce genomic mutations and cell transformation by disrupting gene expression. Several lines of evidence support a role of IAPs as mouse-specific genetic factors in responses to toxicity and expression of disease phenotypes. Since multiple subtypes and copies of IAPs are present in the mouse genome, their activity and locations relative to functional genes are of critical importance. This study identified the major "active" subtypes of IAPs (subtype 1/1a) that are responsible for newly transposed IAP insertions described in the literature, and confirmed that (1) polymorphisms for IAP insertions exist among different mouse strains and (2) promoter activity of the LTRs can be modulated by chemicals. This study further identified all the genes in the C57BL/6 mouse genome with IAP subtype 1 and 1a sequences inserted in their proximity, and the major biofunctional categories and cellular signaling networks of those genes. Since many "IAP-associated genes" play important roles in the regulation of cell proliferation, cell cycle, and cell death, the associated IAPs, upon activation, can affect cellular responses to xenobiotics and disease processes, especially carcinogenesis. This systemic analysis provides a solid foundation for further investigations of the role of IAPs as species- and strain-specific disease susceptibility factors.

Unique functions of repetitive transcriptomes

Repetitive sequences occupy a huge fraction of essentially every eukaryotic genome. Repetitive sequences cover more than 50% of mammalian genomic DNAs, whereas gene exons and protein-coding sequences occupy only ~3% and 1%, respectively. Numerous genomic repeats include genes themselves. They generally encode "selfish" proteins necessary for the proliferation of transposable elements (TEs) in the host genome. The major part of evolutionary "older" TEs accumulated mutations over time and fails to encode functional proteins. However, repeats have important functions also on the RNA level. Repetitive transcripts may serve as multifunctional RNAs by participating in the antisense regulation of gene activity and by competing with the host-encoded transcripts for cellular factors. In addition, genomic repeats include regulatory sequences like promoters, enhancers, splice sites, polyadenylation signals, and insulators, which actively reshape cellular transcriptomes. TE expression is tightly controlled by the host cells, and some mechanisms of this regulation were recently decoded. Finally, capacity of TEs to proliferate in the host genome led to the development of multiple biotechnological applications.

Large-scale discovery of insertion hotspots and preferential integration sites of human transposed elements

Throughout evolution, eukaryotic genomes have been invaded by transposable elements (TEs). Little is known about the factors leading to genomic proliferation of TEs, their preferred integration sites and the molecular mechanisms underlying their insertion. We analyzed hundreds of thousands nested TEs in the human genome, i.e. insertions of TEs into existing ones. We first discovered that most TEs insert within specific ‘hotspots’ along the targeted TE. In particular, retrotransposed Alu elements contain a non-canonical single nucleotide hotspot for insertion of other Alu sequences. We next devised a method for identification of integration sequence motifs of inserted TEs that are conserved within the targeted TEs. This method revealed novel sequences motifs characterizing insertions of various important TE families: Alu, hAT, ERV1 and MaLR. Finally, we performed a global assessment to determine the extent to which young TEs tend to nest within older transposed elements and identified a 4-fold higher tendency of TEs to insert into existing TEs than to insert within non-TE intergenic regions. Our analysis demonstrates that TEs are highly biased to insert within certain TEs, in specific orientations and within specific targeted TE positions. TE nesting events also reveal new characteristics of the molecular mechanisms underlying transposition.

Retroelements and their impact on genome evolution and functioning

Retroelements comprise a considerable fraction of eukaryotic genomes. Since their initial discovery by Barbara McClintock in maize DNA, retroelements have been found in genomes of almost all organisms. First considered as a "junk DNA" or genomic parasites, they were shown to influence genome functioning and to promote genetic innovations. For this reason, they were suggested as an important creative force in the genome evolution and adaptation of an organism to altered environmental conditions. In this review, we summarize the up-to-date knowledge of different ways of retroelement involvement in structural and functional evolution of genes and genomes, as well as the mechanisms generated by cells to control their retrotransposition.

Retrotransposons and non-protein coding RNAs

Retrotransposons constitute a significant fraction of mammalian genomes. Considering the finding of widespread transcriptional activity across entire genomes, it is not surprising that retrotransposons contribute to the collective RNA pool. However, the transcriptional output from retrotransposons does not merely represent spurious transcription. We review examples of functional RNAs transcribed from retrotransposons, and address the collection of non-protein coding RNAs derived from transposable element sequences, including numerous human microRNAs and the neuronal BC RNAs. Finally, we review the emerging understanding of how retrotransposons themselves are regulated by small RNAs.

The role of transposable elements in the regulation of IFN-lambda1 gene expression

IFNs lambda1, lambda2, and lambda3, or type III IFNs, are recently identified cytokines distantly related to type I IFNs. Despite an early evolutionary divergence, the 2 types of IFNs display similar antiviral activities, and both are produced primarily in dendritic cells. Although virus induction of the type I IFN-beta gene had served as a paradigm of gene regulation, relatively little is known about the regulation of IFN-lambda gene expression. Studies of virus induction of IFN-lambda1 identified an essential role of IFN regulatory factors (IRF) 3 and 7, which bind to a regulatory DNA sequence near the start site of transcription. Here, we report that the proximal promoter region of the IFN-lambda1 regulatory region is not sufficient for maximal gene induction in response to bacterial LPS, and we identify an essential cluster of homotypic NF-kappaB binding sites. Remarkably, these sites, which bind efficiently to NF-kappaB and function independently of the IRF3/7 binding sites, originate as transposable elements of the Alu and LTR families. We also show that depletion of the NF-kappaB RelA protein significantly reduces the level of the IFN-lambda1 gene expression. We conclude that IFN-lambda1 gene expression requires NF-kappaB, and we propose a model for IFN-lambda1 gene regulation, in which IRF and NF-kappaB activate gene expression independently via spatially separated promoter elements. These observations provide insights into the independent evolution of the IFN-lambda1 and IFN-beta promoters and directly implicate transposable elements in the regulation of the IFN-lambda1 gene by NF-kappaB.

Human-specific modulation of transcriptional activity provided by endogenous retroviral insertions

Many phenotypic differences exist between Homo sapiens and its closest relatives, chimpanzees, and these differences can arise as a result of variations in the regulation of certain genes common to these closely related species. Human-specific endogenous retroviruses (HERVs) and their solitary long terminal repeats (LTRs) are probable candidates for such a role due to the presence of regulatory elements, such as enhancers, promoters, splice sites, and polyadenylation signals. In this study we show for the first time that HERVs can participate in the specific antisense regulation of human gene expression owing to their LTR promoter activity. We found that two HERV LTRs situated in the introns of genes SLC4A8 (for sodium bicarbonate cotransporter) and IFT172 (for intraflagellar transport protein 172) in the antisense orientation serve in vivo as promoters for generating RNAs complementary to the exons of enclosing genes. The antisense transcripts formed from LTR promoter were shown to decrease the mRNA level of the corresponding genes. The human-specific regulation of these genes suggests their involvement in the evolutionary process.

Integration target site selection by a resurrected human endogenous retrovirus

At least 8% of the human genome was formed by integration of retroviral DNA sequences. Here we analyze the forces directing the accumulation of human endogenous retroviruses (HERVs) by comparing de novo HERV integration targeting with the distribution of fixed HERV elements in the human genome. All known genomic HERVs are inactive due to mutation, but we were able to study integration targeting using a reconstituted consensus HERV-K (designated HERV-K(Con)). We found that HERV-K(Con) integrated preferentially in transcription units, in gene-rich regions, and near features associated with active transcription units and associated regulatory regions. In contrast, genomic HERV-K proviruses are found preferentially outside transcription units. The minority of genomic HERVKs present inside transcription units are in opposite transcriptional orientation relative to the host gene, the orientation predicted to be minimally disruptive to host mRNA synthesis, but de novo HERV-K(Con) integration within transcription units showed no orientation bias. We also found that the youngest HERV-K elements in the human genome showed a distribution intermediate between de novo HERV-K(Con) integration sites and older fixed HERV-Ks. These findings indicate that accumulation of HERVs in the human germline is a two-step process: integration targeting biases direct initial accumulation, then purifying selection leads to loss of proviruses disrupting gene function.

Effects of retroviruses on host genome function

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.

Trans-natural antisense transcripts including noncoding RNAs in 10 species: implications for expression regulation

Natural antisense transcripts are at least partially complementary to their sense transcripts. Cis-Sense/Antisense pairs (cis-SAs) have been extensively characterized and known to play diverse regulatory roles, whereas trans-Sense/Antisense pairs (trans-SAs) in animals are poorly studied. We identified long trans-SAs in human and nine other animals, using ESTs to increase coverage significantly over previous studies. The percentage of transcriptional units (TUs) involved in trans-SAs among all TUs was as high as 4.13%. Particularly 2896 human TUs (or 2.89% of all human TUs) were involved in 3327 trans-SAs. Sequence complementarities over multiple segments with predicted RNA hybridization indicated that some trans-SAs might have sophisticated RNA-RNA pairing patterns. One-fourth of human trans-SAs involved noncoding TUs, suggesting that many noncoding RNAs may function by a trans-acting antisense mechanism. TUs in trans-SAs were statistically significantly enriched in nucleic acid binding, ion/protein binding and transport and signal transduction functions and pathways; a significant number of human trans-SAs showed concordant or reciprocal expression pattern; a significant number of human trans-SAs were conserved in mouse. This evidence suggests important regulatory functions of trans-SAs. In 30 cases, trans-SAs were related to cis-SAs through paralogues, suggesting a possible mechanism for the origin of trans-SAs. All trans-SAs are available at http://trans.cbi.pku.edu.cn/.

Impact of endogenous intronic retroviruses on major histocompatibility complex class II diversity and stability

The major histocompatibility complex (MHC) represents a multigene family that is known to display allelic and gene copy number variations. Primate species such as humans, chimpanzees (Pan troglodytes), and rhesus macaques (Macaca mulatta) show DRB region configuration polymorphism at the population level, meaning that the number and content of DRB loci may vary per haplotype. Introns of primate DRB alleles differ significantly in length due to insertions of transposable elements as long endogenous retrovirus (ERV) and human ERV (HERV) sequences in the DRB2, DRB6, and DRB7 pseudogenes. Although the integration of intronic HERVs resulted sooner or later in the inactivation of the targeted genes, the fixation of these endogenous retroviral segments over long time spans seems to have provided evolutionary advantage. Intronic HERVs may have integrated in a sense or an antisense manner. On the one hand, antisense-oriented retroelements such as HERV-K14I, observed in intron 2 of the DRB7 genes in humans and chimpanzees, seem to promote stability, as configurations/alleles containing these hits have experienced strong conservative selection during primate evolution. On the other hand, the HERVK3I present in intron 1 of all DRB2 and/or DRB6 alleles tested so far integrated in a sense orientation. The data suggest that multigenic regions in particular may benefit from sense introgressions by HERVs, as these elements seem to promote and maintain the generation of diversity, whereas these types of integrations may be lethal in monogenic systems, since they are known to influence transcript regulation negatively.

Genome-wide assessments reveal extremely high levels of polymorphism of two active families of mouse endogenous retroviral elements

Endogenous retroviral elements (ERVs) in mice are significant genomic mutagens, causing ∼10% of all reported spontaneous germ line mutations in laboratory strains. The majority of these mutations are due to insertions of two high copy ERV families, the IAP and ETn/MusD elements. This significant level of ongoing retrotranspositional activity suggests that inbred mice are highly variable in content of these two ERV groups. However, no comprehensive genome-wide studies have been performed to assess their level of polymorphism. Here we compared three test strains, for which sufficient genomic sequence is available, to each other and to the reference C57BL/6J genome and detected very high levels of insertional polymorphism for both ERV families, with an estimated false discovery rate of only 0.4%. Specifically, we found that at least 60% of IAP and 25% of ETn/MusD elements detected in any strain are absent in one or more of the other three strains. The polymorphic nature of a set of 40 ETn/MusD elements found within gene introns was confirmed using genomic PCR on DNA from a panel of mouse strains. For some cases, we detected gene-splicing abnormalities involving the ERV and obtained additional evidence for decreased gene expression in strains carrying the insertion. In total, we identified nearly 700 polymorphic IAP or ETn/MusD ERVs or solitary LTRs that reside in gene introns, providing potential candidates that may contribute to gene expression differences among strains. These extreme levels of polymorphism suggest that ERV insertions play a significant role in genetic drift of mouse lines.

The laboratory mouse is the most widely used mammal for biological research. Hundreds of inbred mouse strains have been developed that vary in characteristics such as susceptibility to cancer or other diseases. There is much interest in uncovering differences between strains that result in different traits and, to aid this effort, millions of single nucleotide differences or polymorphisms between strains have been cataloged. To date, there has been less emphasis placed on other sources of genetic variation. In this study, we have conducted a genome-wide analysis to examine the level of polymorphism of mouse endogenous retroviral sequences (ERVs). ERVs are derived from infectious retroviruses that now exist in the genome and are inherited as part of chromosomes. Unlike in humans, genomic insertions of ERVs cause many new mutations in mice but their extent of variation between strains has been difficult to study because of their high copy numbers. By comparing genomic sequences of four common mouse strains, we found very high levels of polymorphism for two large active families of ERVs. Moreover, we documented nearly 700 polymorphic ERVs located within gene introns and found evidence that some of these affect gene transcript levels. This study demonstrates that ERV polymorphisms are a major source of genetic variability among mouse strains and likely contribute to strain-specific traits.

Species-specific endogenous retroviruses shape the transcriptional network of the human tumor suppressor protein p53

The evolutionary forces that establish and hone target gene networks of transcription factors are largely unknown. Transposition of retroelements may play a role, but its global importance, beyond a few well described examples for isolated genes, is not clear. We report that LTR class I endogenous retrovirus (ERV) retroelements impact considerably the transcriptional network of human tumor suppressor protein p53. A total of 1,509 of approximately 319,000 human ERV LTR regions have a near-perfect p53 DNA binding site. The LTR10 and MER61 families are particularly enriched for copies with a p53 site. These ERV families are primate-specific and transposed actively near the time when the New World and Old World monkey lineages split. Other mammalian species lack these p53 response elements. Analysis of published genomewide ChIP data for p53 indicates that more than one-third of identified p53 binding sites are accounted for by ERV copies with a p53 site. ChIP and expression studies for individual genes indicate that human ERV p53 sites are likely part of the p53 transcriptional program and direct regulation of p53 target genes. These results demonstrate how retroelements can significantly shape the regulatory network of a transcription factor in a species-specific manner.

Widely variable endogenous retroviral methylation levels in human placenta

It is generally assumed that transposable elements, including endogenous retroviruses (ERVs), are silenced by DNA methylation/chromatin structure in mammalian cells. However, there have been very few experimental studies to examine the methylation status of human ERVs. In this study, we determined and compared the methylation status of the 5′ long terminal repeats (LTRs) of different copies of the human endogenous retrovirus (HERV) family HERV-E, which are inserted in various genomic contexts. We found that three HERV-E LTRs which function as alternative gene promoters in placenta are unmethylated in that tissue but heavily methylated in blood cells, where these LTRs are not active promoters. This difference is not solely due to global hypomethylation in placenta, since two general measures of methylation levels of HERV-E and HERV-K LTRs suggest only 10–15% lower overall HERV methylation in placenta compared to blood. Comparisons between methylation levels of the LTR-derived gene promoters and six random HERV-E LTRs in placenta showed that the former display significantly lower methylation levels than random LTRs. Moreover, the differences in methylation between LTRs cannot always be explained by their genomic environment, since methylation of flanking sequences can be very different from methylation of the LTR itself.

Retrotransposons influence the mouse transcriptome: implication for the divergence of genetic traits

Massive accumulation of retrotransposons, comprising >40% of human and mouse genomes, is one of the major events in the evolution of the genome. However, most retrotransposons have lost retrotransposition competency, which makes studying their role in genome evolution elusive. Intracisternal A-particle (IAP) elements are long terminal repeat (LTR)-type mouse retrotransposons consisting of full-length and internally deleted types. Some are retrotransposition competent and their upregulated activity has been reported in mutant mice deficient in genome defense systems, suggesting that IAP elements provide a unique platform for studying the interaction between retrotransposons and mammalian genomes. Using the IAP element as a model case, here we show that mobilization of retrotransposons alters the mouse transcriptome. Retrotransposition assay in cultured cells demonstrated that a subset of internally deleted IAP elements, called IDelta1 type, retrotranspose efficiently when supplied with functional IAP proteins. Furthermore, the IDelta1 type IAP element exhibited substantial transcription-inducing activity in the flanking region. Genomewide transcript analysis of embryonic stem (ES) cells identified IAP-induced transcripts, including fusion transcripts between IAP sequence and endogenous genes. Unexpectedly, nearly half of these IAP elements obtained from ES cells derived from 129 mouse strain were absent in the C57BL/6 genome, suggesting that IAP-driven transcription contributes to the unique trait of the individual mouse strain. On the basis of these data, we propose that retrotransposons are one of the drivers that shape the mammalian transcriptome.