Pedigree-Defined Haplotypes and Their Applications to Genetic Studies

Characterizing the diversity of MHC conserved extended haplotypes using families from the United Arab Emirates

Aside from its anthropological relevance, the characterization of the allele frequencies of genes in the human Major Histocompatibility Complex (MHC) and the combination of these alleles that make up MHC conserved extended haplotypes (CEHs) is necessary for histocompatibility matching in transplantation as well as mapping disease association loci. The structure and content of the MHC region in Middle Eastern populations remain poorly characterized, posing challenges when establishing disease association studies in ethnic groups that inhabit the region and reducing the capacity to translate genetic research into clinical practice. This study was conceived to address a gap of knowledge, aiming to characterize CEHs in the United Arab Emirates (UAE) population through segregation analysis of high-resolution, pedigree-phased, MHC haplotypes derived from 41 families. Twenty per cent (20.5%) of the total haplotype pool derived from this study cohort were identified as putative CEHs in the UAE population. These consisted of CEHs that have been previously detected in other ethnic groups, including the South Asian CEH 8.2 [HLA- C*07:02-B*08:01-DRB1*03:01-DQA1*05:01-DQB1*02:01 (H.F. 0.094)] and the common East Asian CEH 58.1 [HLA- C*03:02-B*58:01-DRB1*03:01- DQA1*05:01-DQB1*02:01 (H.F. 0.024)]. Additionally, three novel CEHs were identified in the current cohort, including HLA- C*15:02-B*40:06-DRB1*16:02-DQB1*05:02 (H.F. 0.035), HLA- C*16:02-B*51:01-DRB1*16:01-DQA1*01:02-DQB1*05:02 (H.F. 0.029), and HLA- C*03:02-B*58:01-DRB1*16:01-DQA1*01:02-DQB1*05:02 (H.F. 0.024). Overall, the results indicate a substantial gene flow with neighbouring ethnic groups in the contemporary UAE population including South Asian, East Asian, African, and European populations. Importantly, alleles and haplotypes that have been previously associated with autoimmune diseases (e.g., Type 1 Diabetes) were also present. In this regard, this study emphasizes that an appreciation for ethnic differences can provide insights into subpopulation-specific disease-related polymorphisms, which has remained a difficult endeavour.

New insights regarding origin of monosomy occurrence in early developing embryos as demonstrated in preimplantation genetic testing

Introduction

Analyses of miscarriage products indicate that the majority of aneuploidies in early developing embryos derive from errors occurring during maternal meiosis and the paternal contribution is less than 10%. Our aim was to assess the aneuploidy (mainly monosmies) frequencies at the earliest stages of embryo development, 3 days following fertilization during In vitro fertilization (IVF) treatments and to elucidate their parental origin. Later, we compared monosomies rates of day 3 to those of day 5 as demonstrated from Preimplantation Genetic Testing for Structural chromosomal Rearrangement (PGT-SR) results.

Methods

For a retrospective study, we collected data of 210 Preimplantation Genetic Testing for Monogenic Disorder (PGT-M) cycles performed between years 2008 and 2019.This study includes 2083 embryos, of 113 couples. It also included 432 embryos from 90 PGT-SR cycles of other 45 patients, carriers of balanced translocations. Defining the parental origin of aneuploidy in cleavage stage embryos was based on haplotypes analysis of at least six informative markers flanking the analyzed gene. For comprehensive chromosomal screening (CCS), chromosomal microarray (CMA) and next generation sequencing (NGS) was used.

Results

We inspected haplotype data of 40 genomic regions, flanking analyzed genes located on 9 different chromosomes.151 (7.2%) embryos presented numerical alterations in the tested chromosomes. We found similar paternal and maternal contribution to monosomy at cleavage stage. We demonstrated paternal origin in 51.5% of the monosomy, and maternal origin in 48.5% of the monosomies cases.

Conclusion

In our study, we found equal parental contribution to monosomies in cleavage-stage embryos. Comparison to CCS analyses of PGT-SR patients revealed a lower rate of monosomy per chromosome in embryos at day 5 of development. This is in contrast to the maternal dominancy described in studies of early miscarriage. Mitotic errors and paternal involvement in chemical pregnancies and IVF failure should be re-evaluated. Our results show monosomies are relatively common and may play a role in early development of ART embryos.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13039-022-00582-5.

The Path to Conserved Extended Haplotypes: Megabase-Length Haplotypes at High Population Frequency

This minireview describes the history of the conceptual development of conserved extended haplotypes (CEHs): megabase-length haplotypes that exist at high (≥0.5%) population frequency. My career began in internal medicine, shifted to pediatrics, and clinical practice changed to research. My research interest was initially in hematology: on plasma proteins, their metabolism, synthesis, and function. This narrowed to a focus on proteins of the human complement system, their role in immunity and their genetics, beginning with polymorphism and deficiency of C3. My group identified genetic polymorphisms and/or inherited deficiencies of C2, C4, C6, and C8. After defining glycine-rich beta glycoprotein as factor B (Bf) in the properdin system, we found that the genes for Bf (CFB), C2, C4A, and C4B were inherited as a single haplotypic unit which we named the "complotype." Complotypes are located within the major histocompatibility complex (MHC) between HLA-B and HLA-DRB1 and are designated (in arbitrary order) by their CFB, C2, C4A, and C4B types. Pedigree analysis revealed long stretches (several megabases) of apparently fixed DNA within the MHC that we referred to as "extended haplotypes" (later as "CEHs"). About 10 to 12 common CEHs constitute at least 25 - 30% of MHC haplotypes among European Caucasian populations. These CEHs contain virtually all the most common markers of MHC-associated diseases. In the case of type 1 diabetes, we have proposed a purely genetic and epigenetic model (with a small number of Mendelian recessive disease genes) that explains all the puzzling features of the disease, including its rising incidence.

Segregation Analysis of Genotyped and Family-Phased, Long Range MHC Classical Class I and Class II Haplotypes in 5 Families With Type 1 Diabetes Proband in the United Arab Emirates

The classical Human Leucocyte Antigen (HLA) class II haplotypes of the Major Histocompatibility Complex (MHC) that are associated with type 1 diabetes (T1D) were identified in five families from the United Arab Emirates (UAE). Segregation analyses were performed on these 5 families with the disease, 3 with one child and 2 with 2 children diagnosed with T1D. Three HLA-DR4 haplotypes were identified: HLA- DRB1^∗04:01:01-DQB1^∗03:02:01:01; HLA- DRB1^∗04:02:01- DQB1^∗03:02:01; and HLA -DRB1^∗04:05:01-DQB1^∗02:02:01:02. All have previously been identified to be associated with T1D in studies of the Arabian population. In the 10 parents from the 5 families, 9 had at least one HLA-DR4 and HLA-DR3 haplotype which potentially increases the risk of T1D. Of these 9 parents, 3 were heterozygous for HLA-DR4/HLA-DR3 and one was homozygous for HLA-DR3. Two haplotypes that were identified here extend to the HLA class I region were previously designated AH8.2 (HLA -A^∗26-B^∗08-DRB1^∗03) and AH50.2 (HLA -C^∗06-B^∗50-DRB1^∗03:01-DQ^∗02) and associated with diabetes in neighboring North Indian populations. This study provides examples of MHC haplotype analysis in pedigrees to improve our understanding of the genetics of T1D in the understudied population of the UAE.

Haplotype Shuffling and Dimorphic Transposable Elements in the Human Extended Major Histocompatibility Complex Class II Region

The major histocompatibility complex (MHC) on chromosome 6p21 is one of the most single-nucleotide polymorphism (SNP)-dense regions of the human genome and a prime model for the study and understanding of conserved sequence polymorphisms and structural diversity of ancestral haplotypes/conserved extended haplotypes. This study aimed to follow up on a previous analysis of the MHC class I region by using the same set of 95 MHC haplotype sequences downloaded from a publicly available BioProject database at the National Center for Biotechnology Information to identify and characterize the polymorphic human leukocyte antigen (HLA)-class II genes, the MTCO3P1 pseudogene alleles, the indels of transposable elements as haplotypic lineage markers, and SNP-density crossover (XO) loci at haplotype junctions in DNA sequence alignments of different haplotypes across the extended class II region (∼1 Mb) from the telomeric PRRT1 gene in class III to the COL11A2 gene at the centromeric end of class II. We identified 42 haplotypic indels (20 Alu, 7 SVA, 13 LTR or MERs, and 2 indels composed of a mosaic of different transposable elements) linked to particular HLA-class II alleles. Comparative sequence analyses of 136 haplotype pairs revealed 98 unique XO sites between SNP-poor and SNP-rich genomic segments with considerable haplotype shuffling located in the proximity of putative recombination hotspots. The majority of XO sites occurred across various regions including in the vicinity of MTCO3P1 between HLA-DQB1 and HLA-DQB3, between HLA-DQB2 and HLA-DOB, between DOB and TAP2, and between HLA-DOA and HLA-DPA1, where most XOs were within a HERVK22 sequence. We also determined the genomic positions of the PRDM9-recombination suppression sequence motif ATCCATG/CATGGAT and the PRDM9 recombination activation partial binding motif CCTCCCCT/AGGGGAG in the class II region of the human reference genome (NC_ 000006) relative to published meiotic recombination positions. Both the recombination and anti-recombination PRDM9 binding motifs were widely distributed throughout the class II genomic regions with 50% or more found within repeat elements; the anti-recombination motifs were found mostly in L1 fragmented repeats. This study shows substantial haplotype shuffling between different polymorphic blocks and confirms the presence of numerous putative ancestral recombination sites across the class II region between various HLA class II genes.

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.

MHC Genomics and Disease: Looking Back to Go Forward

Ancestral haplotypes are conserved but extremely polymorphic kilobase sequences, which have been faithfully inherited over at least hundreds of generations in spite of migration and admixture. They carry susceptibility and resistance to diverse diseases, including deficiencies of CYP21 hydroxylase (47.1) and complement components (18.1), as well as numerous autoimmune diseases (8.1). The haplotypes are detected by segregation within ethnic groups rather than by SNPs and GWAS. Susceptibility to some other diseases is carried by specific alleles shared by multiple ancestral haplotypes, e.g., ankylosing spondylitis and narcolepsy. The difference between these two types of association may explain the disappointment with many GWAS. Here we propose a pathway for combining the two different approaches. SNP typing is most useful after the conserved ancestral haplotypes have been defined by other methods.