Fine mapping the MHC region identified rs4997052 as a new variant associated with nonobstructive azoospermia in Han Chinese males

Partitioning genetic effects on birthweight at classical human leukocyte antigen loci into maternal and fetal components, using structural equation modelling

BACKGROUND

Single nucleotide polymorphisms in the human leukocyte antigen (HLA) region in both maternal and fetal genomes have been robustly associated with birthweight (BW) in previous genetic association studies. However, no study to date has partitioned the association between BW and classical HLA alleles into maternal and fetal components.

METHODS

We used structural equation modelling (SEM) to estimate the maternal and fetal effects of classical HLA alleles on BW. Our SEM leverages the data structure of the UK Biobank (UKB), which includes ∼270 000 participants' own BW and/or the BW of their firstborn child.

RESULTS

We show via simulation that our model yields asymptotically unbiased estimates of the maternal and fetal allelic effects on BW and appropriate type I error rates, in contrast to simple regression models. Asymptotic power calculations show that we have sufficient power to detect moderate-sized maternal or fetal allelic effects of common HLA alleles on BW in the UKB. Applying our SEM to imputed classical HLA alleles and own and offspring BW from the UKB replicated the previously reported association at the HLA-C locus and revealed strong evidence for maternal (HLA-A*03:01, B*35:01, B*39:06, P <0.001) and fetal allelic effects (HLA-B*39:06, P <0.001) of non-HLA-C alleles on BW.

CONCLUSIONS

Our model yields asymptotically unbiased estimates, appropriate type I error rates and appreciable power to estimate maternal and fetal effects on BW. These novel allelic associations between BW and classical HLA alleles provide insight into the immunogenetics of fetal growth in utero.

Immune and spermatogenesis-related loci are involved in the development of extreme patterns of male infertility

We conducted a genome-wide association study in a large population of infertile men due to unexplained spermatogenic failure (SPGF). More than seven million genetic variants were analysed in 1,274 SPGF cases and 1,951 unaffected controls from two independent European cohorts. Two genomic regions were associated with the most severe histological pattern of SPGF, defined by Sertoli cell-only (SCO) phenotype, namely the MHC class II gene HLA-DRB1 (rs1136759, P = 1.32E-08, OR = 1.80) and an upstream locus of VRK1 (rs115054029, P = 4.24E-08, OR = 3.14), which encodes a protein kinase involved in the regulation of spermatogenesis. The SCO-associated rs1136759 allele (G) determines a serine in the position 13 of the HLA-DRβ1 molecule located in the antigen-binding pocket. Overall, our data support the notion of unexplained SPGF as a complex trait influenced by common variation in the genome, with the SCO phenotype likely representing an immune-mediated condition.

A GWAS in a large case-control cohort of European ancestry identifies two genomic regions, the MHC class II gene HLA-DRB1 and an upstream locus of VRK1, that are associated with the most severe phenotype of spermatogenic failure.

Fine Mapping of the MHC Region Identifies Novel Variants Associated with HBV-Related Hepatocellular Carcinoma in Han Chinese

Introduction

Genome-wide association studies identified susceptibility loci in the major histocompatibility complex region for hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC). However, the causal variants underlying HBV-related HCC pathogenesis remain elusive.

Methods

With a total of 1,161 HBV-related HCC cases and 1,353 chronic HBV carriers without HCC, we imputed human leukocyte antigen (HLA) variants based on a Chinese HLA reference panel and evaluated the associations of these variants with the risk of HBV-related HCC. Conditional analyses were used to identify independent signals associated with the risk of HBV-related HCC (P false-discovery rate (FDR) <0.20). A total of 14,930 variants within the MHC region were genotyped or imputed.

Results

We identified two variants, rs114401688 (P = 1.05 × 10⁻⁶, P_FDR = 2.43 × 10⁻³) and rs115126566 (P = 9.04 × 10⁻⁵, P_FDR = 1.77 × 10⁻¹), that are independently associated with the risk of HBV-related HCC. Single nucleotide polymorphism (SNP) rs114401688 is in linkage disequilibrium with a previously reported SNP rs9275319. In the current study, we found that its association with HCC could be explained by HLA-DQB1*04 and HLA-DRB1*04. SNP rs115126566 is a novel risk variant and may function by regulating transcriptions of HLA-DPA1/DPB1 through enhancer-mediated mechanisms. HLA zygosity analysis showed that homozygosity at HLA-DQB1 gene is suggestively associated with a higher risk of HCC (P = 0.10) and the risk was more pronounced in the older age group (age ≥50, P = 0.03).

Discussion

Our findings further the understanding of the genetic basis for HBV-related HCC predisposition in chronic HBV carriers.

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 Landscape of Nonobstructive Azoospermia and New Perspectives for the Clinic

Nonobstructive azoospermia (NOA) represents the most severe expression of male infertility, involving around 1% of the male population and 10% of infertile men. This condition is characterised by the inability of the testis to produce sperm cells, and it is considered to have an important genetic component. During the last two decades, different genetic anomalies, including microdeletions of the Y chromosome, karyotype defects, and missense mutations in genes involved in the reproductive function, have been described as the primary cause of NOA in many infertile men. However, these alterations only explain around 25% of azoospermic cases, with the remaining patients showing an idiopathic origin. Recent studies clearly suggest that the so-called idiopathic NOA has a complex aetiology with a polygenic inheritance, which may alter the spermatogenic process. Although we are far from a complete understanding of the molecular mechanisms underlying NOA, the use of the new technologies for genetic analysis has enabled a considerable increase in knowledge during the last years. In this review, we will provide a comprehensive and updated overview of the genetic basis of NOA, with a special focus on the possible application of the recent insights in clinical practice.