Genomic and transcriptomic studies of expression quantitative trait loci (eQTL) revealed that SINE-VNTR-Alu (SVA) retrotransposon insertion polymorphisms (RIPs) within human genomes markedly affect the co-expression of many coding and noncoding genes by coordinated regulatory processes. This study examined the polymorphic SVA modulation of gene co-expression within the major histocompatibility complex (MHC) genomic region where more than 160 coding genes are involved in innate and adaptive immunity. We characterized the modulation of SVA RIPs utilizing the genomic and transcriptomic sequencing data obtained from whole blood of 1266 individuals in the Parkinson’s Progression Markers Initiative (PPMI) cohort that included an analysis of human leukocyte antigen (HLA)-A regulation in a subpopulation of the cohort. The regulatory properties of eight SVAs located within the class I and class II MHC regions were associated with differential co-expression of 71 different genes within and 75 genes outside the MHC region. Some of the same genes were affected by two or more different SVA. Five SVA are annotated in the human genomic reference sequence GRCh38.p14/hg38, whereas the other three were novel insertions within individuals. We also examined and found distinct structural effects (long and short variants and the CT internal variants) for one of the SVA (R_SVA_24) insertions on the differential expression of the HLA-A gene within a subpopulation (550 individuals) of the PPMI cohort. This is the first time that many HLA and non-HLA genes (multilocus expression units) and splicing mechanisms have been shown to be regulated by eight structurally polymorphic SVA within the MHC genomic region by applying precise statistical analysis of RNA data derived from the blood samples of a human cohort population. This study shows that SVA within the MHC region are important regulators or rheostats of gene co-expression that might have potential roles in diversity, health, and disease.
DawkinsRLeelayuwatCGaudieriSTayGHuiJCattleySMartinezPKulskiJ.Genomics of the major histocompatibility complex: haplotypes, duplication, retroviruses and disease.Immunol Rev1999;167:275–304
2.
VandiedonckCKnightJC.The human major histocompatibility complex as a paradigm in genomics research.Brief Funct Genomic Proteomic2009;8:379–94
3.
ShiinaTInokoHKulskiJK.An update of the HLA genomic region, locus information and disease associations: 2004.Tissue Antigens2004;64:631–49
4.
ShiinaTHosomichiKInokoHKulskiJK.The HLA genomic loci map: expression, interaction, diversity and disease.J Hum Genet2009;54:15–39
5.
KulskiJKSuzukiSShiinaT.Human leukocyte antigen super-locus: nexus of genomic supergenes, SNPs, indels, transcripts, and haplotypes.Hum Genome Var2022;9:49
6.
TaitBD.The importance of establishing genetic phase in clinical medicine.Int J Immunogenet2022;49:1–7
7.
BarkerDJMaccariGGeorgiouXCooperMAFlicekPRobinsonJMarshSGE. The IPD-IMGT/HLA database.Nucleic Acids Res2023;51:D1053–60
8.
HortonRGibsonRCoggillPMirettiMAllcockRJAlmeidaJForbesSGilbertJGHallsKHarrowJLHartEHoweKJacksonDKPalmerSRobertsANSimsSStewartCATraherneJATrevanionSWilmingLRogersJde JongPJElliottJFSawcerSToddJATrowsdaleJBeckS.Variation analysis and gene annotation of eight MHC haplotypes: the MHC haplotype project.Immunogenetics2008;60:1–18
9.
TrowsdaleJKnightJC.Major histocompatibility complex genomics and human disease.Annu Rev Genomics Hum Genet2013;14:301–23
10.
YunisEJLarsenCEFernandez-ViñaMAwdehZLRomeroTHansenJAAlperCA.Inheritable variable sizes of DNA stretches in the human MHC: conserved extended haplotypes and their fragments or blocks.Tissue Antigens2003;62:1–20
11.
GoodinDSKhankhanianPGourraudPAVinceN.Highly conserved extended haplotypes of the major histocompatibility complex and their relationship to multiple sclerosis susceptibility.PLoS ONE2018;13:e0190043
12.
MaiersMGragertLKlitzW.High-resolution HLA alleles and haplotypes in the United States population.Hum Immunol2007;68:779–88
13.
NevilleMJLeeWHumburgPWongDBarnardoMKarpeFKnightJC.High resolution HLA haplotyping by imputation for a British population bioresource.Hum Immunol2017;78:242–51
WissemannWTHill-BurnsEMZabetianCPFactorSAPatsopoulosNHoglundBHolcombCDonahueRJThomsonGErlichHPayamiH.Association of Parkinson disease with structural and regulatory variants in the HLA region.Am J Hum Genet2013;93:984–93
16.
NallsMABlauwendraatCVallergaCLHeilbronKBandres-CigaSChangDTanMKiaDANoyceAJXueABrasJYoungEvon CoellnRSimón-SánchezJSchulteCSharmaMKrohnLPihlstrømLSiitonenAIwakiHLeonardHFaghriFGibbsJRHernandezDGScholzSWBotiaJAMartinezMCorvolJCLesageSJankovicJShulmanLMSutherlandMTienariPMajamaaKToftMAndreassenOABangaleTBriceAYangJGan-OrZGasserTHeutinkPShulmanJMWoodNWHindsDAHardyJAMorrisHRGrattenJVisscherPMGrahamRRSingletonAB, 23andMe Research Team, System Genomics of Parkinson’s Disease Consortium, International Parkinson’s Disease Genomics Consortium. Identification of novel risk loci, causal insights, and heritable risk for Parkinson’s disease: a meta-analysis of genome-wide association studies.Lancet Neurol2019;18:1091–102
17.
YuEAmbatiAAndersenMSKrohnLEstiarMASainiPSenkevichKSoseroYLSreelathaAAKRuskeyJAAsayeshFSpiegelmanDToftMVikenMKSharmaMBlauwendraatCPihlstrømLMignotEGan-OrZ.Fine mapping of the HLA locus in Parkinson’s disease in Europeans.NPJ Park Dis2021;7:84
18.
LokkiMPaakkanenR.The complexity and diversity of major histocompatibility complex challenge diseasee association studies.HLA2018;93:3–15
19.
TraherneJA.Human MHC architecture and evolution: implications for disease association studies.Int J Immunogenet2008;35:179–92
20.
KulskiJKSuzukiSShiinaT.SNP-density crossover maps of polymorphic transposable elements and HLA genes within MHC class I haplotype blocks and junction.Front Genet2021;11:594318
21.
KulskiJKSuzukiSShiinaT.Haplotype shuffling and dimorphic transposable elements in the human extended major histocompatibility complex class II region.Front Genet2021;12:665899
22.
CunYShiLKulskiJKLiuSYangJTaoYZhangXShiLYaoY.Haplotypic associations and differentiation of MHC class II polymorphic alu insertions at five loci with HLA-DRB1 alleles in 12 minority ethnic populations in China.Front Genet2021;12:636236
23.
KulskiJKMawartAMarieKTayGKAlSafarHS.MHC class I polymorphic Alu insertion (POALIN) allele and haplotype frequencies in the Arabs of the United Arab Emirates and other world populations.Int J Immunogenet2019;46:247–62
24.
KulskiJKShigenariAInokoH.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.Tissue Antigens2011;78:359–77
25.
LamontagneMJoubertPTimensWPostmaDSHaoKNickleDSinDDParePDLavioletteMBosséY.Susceptibility genes for lung diseases in the major histocompatibility complex revealed by lung expression quantitative trait loci analysis.Eur Respir J2016;48:573–6
26.
LamTHShenMTayMZRenEC.Unique Allelic eQTL Clusters in Human MHC Haplotypes.G3 Genesgenomesgenetics2017;7:2595–604
27.
KoksSPfaffALBubbVJQuinnJP.Expression quantitative trait loci (eQTLs) associated with retrotransposons demonstrate their modulatory effect on the transcriptome.Int J Mol Sci2021;22:6319
28.
SmalleganMJShehataSSpradlinSFSwearingenAWheelerGDasACorbetGNebenfuehrBAhrensDTauberDLennonSChoiKHuynhTWieserTSchneiderKBradshawMBaskenJLaiMReadTHynes-GraceMTimmonsDDemasiJRinnJL.Genome-wide binding analysis of 195 DNA binding proteins reveals “reservoir” promoters and human specific SVA-repeat family regulation.PLoS ONE2021;16:e0237055
29.
KulskiJK.Long noncoding RNA HCP5, a Hybrid HLA class i endogenous retroviral gene: structure, expression, and disease associations.Cells2019;8:480
30.
FortVKhelifiGHusseinSMI. Long non-coding RNAs and transposable elements: a functional relationship.Biochim Biophys Acta Mol Cell Res2021;1868:118837
31.
WangLRishishwarLMariño-RamírezLJordanIK.Human population-specific gene expression and transcriptional network modification with polymorphic transposable elements.Nucleic Acids Res2016;45:gkw1286
32.
KwonYJChoiYEoJNohYNGimJAJungYDLeeJRKimHS.Structure and expression analyses of SVA elements in relation to functional genes.Genomics Inform2013;11:142–8
33.
VandiedonckCTaylorMSLockstoneHEPlantKTaylorJMDurrantCBroxholmeJFairfaxBPKnightJC.Pervasive haplotypic variation in the spliceo-transcriptome of the human major histocompatibility complex.Genome Res2011;21:1042–54
34.
SpiritoGMangoniDSangesRGustincichS.Impact of polymorphic transposable elements on transcription in lymphoblastoid cell lines from public data.BMC Bioinformatics2019;20:495
35.
MakałowskiWGoteaVPandeAMakałowskaI. Transposable elements: classification, identification, and their use as a tool for comparative genomics. In: AnisimovaM (ed.) Evolutionary genomics. New York: Springer, 2019, pp.177–207
36.
WangJSongLGroverDAzrakSBatzerMALiangP.dbRIP: a highly integrated database of retrotransposon insertion polymorphisms in humans.Hum Mutat2006;27:323–9
37.
AutioMIBin AminTPerrinAWongJYFooRS-YPrabhakarS.Transposable elements that have recently been mobile in the human genome.BMC Genomics2021;22:789
38.
SavageALBubbVJBreenGQuinnJP.Characterisation of the potential function of SVA retrotransposons to modulate gene expression patterns.BMC Evol Biol2013;13:101
39.
KulskiJKShigenariAInokoH.Polymorphic SVA retrotransposons at four loci and their association with classical HLA class I alleles in Japanese, Caucasians and African Americans.Immunogenetics2010;62:211–30
40.
PfaffALBubbVJQuinnJPKoksS.Reference SVA insertion polymorphisms are associated with Parkinson’s disease progression and differential gene expression.NPJ Park Dis2021;7:44
41.
MarekKChowdhurySSiderowfALaschSCoffeyCSCaspell-GarciaCSimuniTJenningsDTannerCMTrojanowskiJQShawLMSeibylJSchuffNSingletonAKieburtzKTogaAWMollenhauerBGalaskoDChahineLMWeintraubDForoudTTosun-TurgutDPostonKArnedoVFrasierMShererT, Parkinson’s Progression Markers Initiative. The Parkinson’s progression markers initiative (PPMI) – establishing a PD biomarker cohort.Ann Clin Transl Neurol2018;5:1460–77
42.
RauschTZichnerTSchlattlAStützAMBenesVKorbelJO.DELLY: structural variant discovery by integrated paired-end and split-read analysis.Bioinformatics2012;28:i333–9
43.
PurcellSNealeBTodd-BrownKThomasLFerreiraMABenderDMallerJSklarPde BakkerPIDalyMJShamPC.PLINK: a tool set for whole-genome association and population-based linkage analyses.Am J Hum Genet2007;81:559–75
44.
PatroRDuggalGLoveMIIrizarryRAKingsfordC.Salmon provides fast and bias-aware quantification of transcript expression.Nat Methods2017;14:417–9
45.
ShabalinAA.Matrix eQTL: ultra fast eQTL analysis via large matrix operations.Bioinformatics2012;28:1353–8
46.
Pérez-CaleroCBayona-FeliuAXueXBarrosoSIMuñozSGonzález-BasalloteVMSungPAguileraA.UAP56/DDX39B is a major cotranscriptional RNA–DNA helicase that unwinds harmful R loops genome-wide.Genes Dev2020;34:898–912
47.
SchottGGarcia-BlancoMA.MHC class III RNA binding proteins and immunity.RNA Biol2021;18:640–6
48.
ZhaoLLiuQHuangJLuYZhaoYPingJ.TREX (transcription/export)-NP complex exerts a dual effect on regulating polymerase activity and replication of influenza A virus.PLoS Pathog2022;18:e1010835
49.
SzymuraSJBernalGMWuLZhangZCrawleyCDVoceDJCampbellP-ARanoaDEWeichselbaumRRYaminiB.DDX39B interacts with the pattern recognition receptor pathway to inhibit NF-κB and sensitize to alkylating chemotherapy.BMC Biol2020;18:32
50.
MulcahyHO’RourkeKPAdamsCMolloyMGO’GaraF.LST1 and NCR3 expression in autoimmune inflammation and in response to IFN-γ, LPS and microbial infection.Immunogenetics2006;57:893–903
51.
FagerbergLHallströmBMOksvoldPKampfCDjureinovicDOdebergJHabukaMTahmasebpoorSDanielssonAEdlundKAsplundASjöstedtELundbergESzigyartoCASkogsMTakanenJOBerlingHTegelHMulderJNilssonPSchwenkJMLindskogCDanielssonFMardinogluASivertssonAvon FeilitzenKForsbergMZwahlenMOlssonINavaniSHussMNielsenJPontenFUhlénM.Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. Mol Cell Proteomics2014;13:397–406
52.
ChittendenTFlemingtonCHoughtonABEbbRGGalloGJElangovanBChinnaduraiGLutzRJ.A conserved domain in Bak, distinct from BH1 and BH2, mediates cell death and protein binding functions.EMBO J1995;14:5589–96
53.
KulskiJKGaudieriSMartinADawkinsRL.Coevolution of PERB11 (MIC) and HLA class I genes with HERV-16 and retroelements by extended genomic duplication.J Mol Evol1999;49:84–97
54.
GianfrancescoOBubbVJQuinnJP.SVA retrotransposons as potential modulators of neuropeptide gene expression.Neuropeptides2017;64:37
55.
GianfrancescoOGearyBSavageALBillingsleyKJBubbVJQuinnJP.The role of SINE-VNTR-Alu (SVA) retrotransposons in shaping the human genome.Int J Mol Sci2019;20:5977
56.
QuinnJPBubbVJ.SVA retrotransposons as modulators of gene expression.Mob Genet Elem2014;4:e32102
57.
KanevskiyLErokhinaSKobyzevaPStreltsovaMSapozhnikovAKovalenkoE.Dimorphism of HLA-E and its disease association.Int J Mol Sci2019;20:5496
58.
WangLNorrisETJordanIK.Human retrotransposon insertion polymorphisms are associated with health and disease via gene regulatory phenotypes.Front Microbiol2017;8:1418
59.
ShenLWuLSanliogluSChenRMendozaARDangelAWCarrollMCZipfWBYuCY.Structure and genetics of the partially duplicated gene RP located immediately upstream of the complement C4A and the C4B genes in the HLA class III region.J Biol Chem1994;269:8466–76
60.
ZhuZBHsiehS-LBentleyDRCampbellRDVolanakisJE.A variable number of tandem repeats locus within the human complement C2 gene is associated with a retroposon derived from a human endogenous retrovirus.J Exp Med1992;175:1783–7
61.
OnoMKawakamiMTakezawaT.A novel human nonviral retroposon derived from an endogenous retrovirus.Nucleic Acids Res1987;15:8725–37
62.
BenjaminiYHochbergY.Controlling the false discovery rate: a practical and powerful approach to multiple testing.J R Stat Soc Ser B Methodol1995;57:289–300
63.
KlenovaEScottACRobertsJShamsuddinSLovejoyEABergmannSBubbVJRoyerH-DQuinnJP.YB-1 and CTCF differentially regulate the 5-HTT polymorphic intron 2 enhancer which predisposes to a variety of neurological disorders.J Neurosci2004;24:5966–73
64.
PaciPFisconGConteFWangR-SFarinaLLoscalzoJ.Gene co-expression in the interactome: moving from correlation toward causation via an integrated approach to disease module discovery.NPJ Syst Biol Appl2021;7:3
65.
BarnadaSMIsopiATejada-MartinezDGoubertCPatooriSPagliaroliLTracewellMTrizzinoM.Genomic features underlie the co-option of SVA transposons as cis-regulatory elements in human pluripotent stem cells.PLoS Genet2022;18:e1010225
66.
PappalardoXGBarraV.Losing DNA methylation at repetitive elements and breaking bad.Epigenetics Chromatin2021;14:25
67.
TrigianteGBlanes RuizNCeraseA.Emerging roles of repetitive and repeat-containing RNA in nuclear and chromatin organization and gene expression.Front Cell Dev Biol2021;9:735527
68.
OttavianiDLeverEMitterRJonesTForshewTChristovaRTomazouEMRakyanVKKrawetzSAPlattsAESegaraneBBeckSSheerD.Reconfiguration of genomic anchors upon transcriptional activation of the human major histocompatibility complex.Genome Res2008;18:1778–86
69.
HancksDCEwingADChenJETokunagaKKazazianHHJr.Exon-trapping mediated by the human retrotransposon SVA.Genome Res2009;19:1983–91
70.
NadlerMJSChangWOzkaynakEHuoYNongYBoillotMJohnsonM, Moreno A, Anderson MP Hominoid SVA-lncRNA AK057321 targets human-specific SVA retrotransposons in SCN8A and CDK5RAP2 to initiate neuronal maturation.Commun Biol2023;6:347
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
