A Neanderthal haplotype introgressed into the human genome confers protection against membranous nephropathy

Class 2 HLA and PLA2R1 alleles are exceptionally strong genetic risk factors for membranous nephropathy (MN), leading, through an unknown mechanism, to a targeted autoimmune response. Introgressed archaic haplotypes (introduced from an archaic human genome into the modern human genome) might influence phenotypes through gene dysregulation. Here, we investigated the genomic region surrounding the PLA2R1 gene. We reconstructed the phylogeny of Neanderthal and modern haplotypes in this region and calculated the probability of the observed clustering being the result of introgression or common descent. We imputed variants for the participants in our previous genome-wide association study and we compared the distribution of Neanderthal variants between MN cases and controls. The region associated with the lead MN risk locus in the PLA2R1 gene was confirmed and showed that, within a 507 kb region enriched in introgressed sequence, a stringently defined 105 kb haplotype, intersecting the coding regions for PLA2R1 and ITGB6, is inherited from Neanderthals. Thus, introgressed Neanderthal haplotypes overlapping PLA2R1 are differentially represented in MN cases and controls, with enrichment In controls suggesting a protective effect.