ERVK9, transposons and the evolution of MHC class I duplicons within the alpha-block of the human and chimpanzee

HLA-OLI: A new MHC class I pseudogene and HLA-Y are located on a 60 kb indel in the human MHC between HLA-W and HLA-J

Analysis of publicly available whole-genome sequence data from the Human Pangenome Project and the 1000 Genomes Project has identified a DNA segment of approximately 60 kb in the major histocompatibility complex (MHC) between HLA-W and HLA-J that is present in some MHC haplotypes but not others. This DNA segment is largely repeat element-rich but includes the pseudogene HLA-Y, thus pinpointing the location of this pseudogene, and a new HLA class I sequence we have called HLA-OLI. HLA-OLI clusters phylogenetically with the HLA class I pseudogenes, HLA-P and HLA-W, and appears to have a similar genetic structure. The availability of whole-genome sequence data from diverse populations enables a detailed characterization of the MHC at the population level and will have implications for understanding MHC disease associations and the non-HLA MHC factors that impact unrelated hematopoietic cell transplant outcomes.

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

Genetic variation and hitchhiking between structurally polymorphic Alu insertions and HLA-A, -B, and -C alleles and other retroelements within the MHC class I region

We investigated structurally polymorphic Alu insertions (POALINs) at five loci in the major histocompatibility complex (MHC) class I genomic region to determine their allele and haplotype frequencies and associations with the human leukocyte antigen (HLA)-A, -B, and -C genes in three populations, the Australian Caucasians, Japanese, and African Americans. The POALINs varied in allelic frequency between 0% and 42.3% with significant differences between populations at three of the five loci. The linkage disequilibrium (LD) between Alu insertions and the HLA-A, -B, or -C alleles and previously published polymorphic retroelements (four SVA and human endogenous retrovirus type 9 (HERVK9) loci) within the class I region of the MHC were calculated in pairwise analyses of haplotypes to show strong allelic associations and possible crossing-over events between some loci. Each POALIN was in significant LD with a variety of HLA-A, -B, or -C two-digit alleles probably as a result of hitchhiking. The POALINs helped to further stratify the HLA-A:B:C haplotypes into different POALIN:HLA-A:B:C haplotype frequencies. Of the multilocus haplotype analyses, the seven- and eight-locus haplotypes showed the largest number of differences between the populations, and fewer matched haplotypes between populations that ranged, for example, from 49% for HLA-B:HLA-A haplotypes to 7% for AluMICB:HLA-B:HLA-C:AluTF:AluHJ:HLA-A:AluHG:AluTF haplotypes in the Japanese. This comparative study of multilocus POALINs in the HLA class I region of three ethnic populations shows that POALINs alone or together with the HLA class I alleles and other retroelements are informative ancestral markers for assessing the interrelationship of HLA class I haplotype lineages, LD, and genetic diversity within the same and/or different populations.

Viral ancestors of antiviral systems

All life must survive their corresponding viruses. Thus antiviral systems are essential in all living organisms. Remnants of virus derived information are also found in all life forms but have historically been considered mostly as junk DNA. However, such virus derived information can strongly affect host susceptibility to viruses. In this review, I evaluate the role viruses have had in the origin and evolution of host antiviral systems. From Archaea through bacteria and from simple to complex eukaryotes I trace the viral components that became essential elements of antiviral immunity. I conclude with a reexamination of the 'Big Bang' theory for the emergence of the adaptive immune system in vertebrates by horizontal transfer and note how viruses could have and did provide crucial and coordinated features.

Viral ancestors of antiviral systems

All life must survive their corresponding viruses. Thus antiviral systems are essential in all living organisms. Remnants of virus derived information are also found in all life forms but have historically been considered mostly as junk DNA. However, such virus derived information can strongly affect host susceptibility to viruses. In this review, I evaluate the role viruses have had in the origin and evolution of host antiviral systems. From Archaea through bacteria and from simple to complex eukaryotes I trace the viral components that became essential elements of antiviral immunity. I conclude with a reexamination of the ‘Big Bang’ theory for the emergence of the adaptive immune system in vertebrates by horizontal transfer and note how viruses could have and did provide crucial and coordinated features.

HERV-mediated genomic rearrangement of EYA1 in an individual with branchio-oto-renal syndrome

Branchio-oto-renal syndrome is characterized by branchial defects, hearing loss, preauricular pits, and renal anomalies. Mutations in EYA1 are the most common cause of branchio-oto-renal and branchio-otic syndromes. Large chromosomal aberrations of 8q13, including complex rearrangements occur in about 20% of these individuals. However, submicroscopic deletions and the molecular characterization of genomic rearrangements involving the EYA1 gene have rarely been reported. Using the array-comparative genomic hybridization, we identified non-recurrent genomic deletions including the EYA1 gene in three patients with branchio-oto-renal syndrome, short stature, and developmental delay. One of these deletions was mediated by two human endogenous retroviral sequence blocks, analogous to the AZFa microdeletion on Yq11, responsible for male infertility. This report describes the expanded phenotype of individuals, resulting from contiguous gene deletion involving the EYA1 gene and provides a molecular description of the genomic rearrangements involving this gene in branchio-oto-renal syndrome.

HLA-A allele associations with viral MER9-LTR nucleotide sequences at two distinct loci within the MHC alpha block

The study of the association of the Human Leukocyte Antigen (HLA) alleles and polymorphic retrotransposons such as Alu, HERV, and LTR at various loci within the Major Histocompatibility Complex allows for a better identification and stratification of disease associations and the origins of HLA haplotypes in different populations. This paper provides sequence and association data on two structurally polymorphic MER9-LTR retrotransposons that are located 54 kb apart and in close proximity to the multiallelic HLA-A gene involved in the regulation of the human immune system. Direct DNA sequencing and analysis of the PCR products identified DNA nucleotide variations between the MER9-LTR sequences at the two loci and their associations with HLA-A alleles as potential haplotype and evolutionary markers. All MER9-LTR sequences were haplotypic when associated with common HLA-A alleles. The number of SNP loci was 2.5 times greater for the solo LTR at the AK locus, which is located closer to the HLA-A gene than the solo or 3' LTR at the HG locus. Our study shows that the nucleotide variations of the MER9-LTR DNA sequences are additional informative markers in fine mapping HLA-A genomic haplotypes for future population, evolutionary, and disease studies.

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.

Compound evolutionary history of the rhesus macaque MHC class I B region revealed by microsatellite analysis and localization of retroviral sequences

In humans, the single polymorphic B locus of the major histocompatibility complex is linked to the microsatellite MIB. In rhesus macaques, however, haplotypes are characterized by the presence of unique combinations of multiple B genes, which may display different levels of polymorphism. The aim of the study was to shed light on the evolutionary history of this highly complex region. First, the robustness of the microsatellite MIB-linked to almost half of the B genes in rhesus macaques (Mamu-B)–for accurate B haplotyping was studied. Based on the physical map of an established haplotype comprising 7 MIB loci, each located next to a certain Mamu-B gene, two MIB loci, MIB1 and MIB6, were investigated in a panel of MHC homozygous monkeys. MIB1 revealed a complex genotyping pattern, whereas MIB6 analysis resulted in the detection of one or no amplicon. Both patterns are specific for a given B haplotype, show Mendelian segregation, and even allow a more precise haplotype definition than do traditional typing methods. Second, a search was performed for retroelements that may have played a role in duplication processes as observed in the macaque B region. This resulted in the description of two types of duplicons. One basic unit comprises an expressed Mamu-B gene, adjacent to an HERV16 copy closely linked to MIB. The second type of duplicon comprises a Mamu-B (pseudo)gene, linked to a truncated HERV16 structure lacking its MIB segment. Such truncation seems to coincide with the loss of B gene transcription. Subsequent to the duplication processes, recombination between MIB and Mamu-B loci appears to have occurred, resulting in a hyperplastic B region. Thus, analysis of MIB in addition to B loci allows deciphering of the compound evolutionary history of the class I B region in Old World monkeys.

Human endogenous retrovirus (HERVK9) structural polymorphism with haplotypic HLA-A allelic associations

The frequency and HLA-A allelic associations of a HERVK9 DNA structural polymorphism located in close proximity to the highly polymorphic HLA-A gene within the major histocompatibility complex (MHC) genomic region were determined in Japanese, African Americans, and Australian Caucasians to better understand its human population evolutionary history. The HERVK9 insertion or deletion was detected as a 3' LTR or a solo LTR, respectively, by separate PCR assays. The average insertion frequency of the HERVK9.HG was significantly different (P < 1.083e(-6)) between the Japanese (0.59) and the African Americans (0.34) or Australian Caucasians (0.37). LD analysis predicted a highly significant (P < 1.0e(-5)) linkage between the HLA-A and HERVK9 alleles, probably as a result of hitchhiking (linkage). Evolutionary time estimates of the solo, 5' and 3' LTR nucleotide sequence divergences suggest that the HERVK9 was inserted 17.3 MYA with the first structural deletion occurring 15.1 MYA. The LTR/HLA-A haplotypes appear to have been formed mostly during the past 3.9 MY. The HERVK9 insertion and deletion, detected by a simple and economical PCR method, is an informative genetic and evolutionary marker for the study of HLA-A haplotype variations, human migration, the origins of contemporary populations, and the possibility of disease associations.

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.

The Widespread Evolutionary Significance of Viruses

In the last 30 years, the study of virus evolution has undergone a transformation. Originally concerned with disease and its emergence, virus evolution had not been well integrated into the general study of evolution. This chapter reviews the developments that have brought us to this new appreciation for the general significance of virus evolution to all life. We now know that viruses numerically dominate all habitats of life, especially the oceans. Theoretical developments in the 1970s regarding quasispecies, error rates, and error thresholds have yielded many practical insights into virus–host dynamics. The human diseases of HIV-1 and hepatitis C virus cannot be understood without this evolutionary framework. Yet recent developments with poliovirus demonstrate that viral fitness can be the result of a consortia, not one fittest type, a basic Darwinian concept in evolutionary biology. Darwinian principles do apply to viruses, such as with Fisher population genetics, but other features, such as reticulated and quasispecies-based evolution distinguish virus evolution from classical studies. The available phylogenetic tools have greatly aided our analysis of virus evolution, but these methods struggle to characterize the role of virus populations. Missing from many of these considerations has been the major role played by persisting viruses in stable virus evolution and disease emergence. In many cases, extreme stability is seen with persisting RNA viruses. Indeed, examples are known in which it is the persistently infected host that has better survival. We have also recently come to appreciate the vast diversity of phage (DNA viruses) of prokaryotes as a system that evolves by genetic exchanges across vast populations (Chapter 10). This has been proposed to be the “big bang” of biological evolution. In the large DNA viruses of aquatic microbes we see surprisingly large, complex and diverse viruses. With both prokaryotic and eukaryotic DNA viruses, recombination is the main engine of virus evolution, and virus host co-evolution is common, although not uniform. Viral emergence appears to be an unending phenomenon and we can currently witness a selective sweep by retroviruses that infect and become endogenized in koala bears.

Class I genes have split from the MHC in the tammar wallaby

Genes within the Major Histocompatibility Complex (MHC) are critical to the immune response and immunoregulation. Comparative studies have revealed that the MHC has undergone many changes throughout evolution yet in tetrapods the three different classes of MHC genes have maintained linkage, suggesting that there may be some functional advantage obtained by maintaining this clustering of MHC genes. Here we present data showing that class II and III genes, the antigen processing gene TAP2, and MHC framework genes are found together in the tammar wallaby on chromosome 2. Surprisingly class I loci were not found on chromosome 2 but were mapped to ten different locations spread across six chromosomes. This distribution of class I loci in the wallaby on nearly all autosomes is not a characteristic of all marsupials and may be a relatively recent phenomenon. It highlights the need for the inclusion of more than one marsupial species in comparative studies and raises questions regarding the functional significance of the clustering of MHC genes.

Polymorphic Alu insertions within the Major Histocompatibility Complex class I genomic region: a brief review

Most polymorphic Alu insertions (POALINs) belong to a subgroup of the Alu multicopy retrotransposon family of short interspersed nucleotide elements (SINEs) that are categorized as AluYb8 and AluYa5. The number of AluYb8/AluYa5 members (approximately 4,492 copies) is significantly less than the approximately one million fixed Alu copies per human genome. We have studied the presence of POALINs within the Major Histocompatibility Complex (MHC) class I region on the short arm of chromosome 6 (6p21.3) because this region has a high gene density, many genes with immune system functions, large sequence variations and diversity, duplications and redundancy, and a strong association with more than 100 different diseases. Since little is known about POALINs within the MHC genomic region, we undertook to identify some of the members of the AluYb8/AluYa5 subfamily and to study their frequency of distribution and genetic characteristics in different populations. As a result of our comparative genomic analyses, we identified the insertion sites for five POALINs distributed within the MHC class I region. This brief review outlines the locations of the insertions and sequence features of the five MHC POALINs, their single site and haplotype frequencies in different geographic populations, and their association with different HLA class I genes and disease. We show that the MHC POALINs have a potential value as lineage and linkage markers for the study of human population genetics, disease associations, genomic diversity and evolution.