Dataset Information

Genomic structure of the horse major histocompatibility complex class II region resolved using PacBio long-read sequencing technology.

ABSTRACT: The mammalian Major Histocompatibility Complex (MHC) region contains several gene families characterized by highly polymorphic loci with extensive nucleotide diversity, copy number variation of paralogous genes, and long repetitive sequences. This structural complexity has made it difficult to construct a reliable reference sequence of the horse MHC region. In this study, we used long-read single molecule, real-time (SMRT) sequencing technology from Pacific Biosciences (PacBio) to sequence eight Bacterial Artificial Chromosome (BAC) clones spanning the horse MHC class II region. The final assembly resulted in a 1,165,328?bp continuous gap free sequence with 35 manually curated genomic loci of which 23 were considered to be functional and 12 to be pseudogenes. In comparison to the MHC class II region in other mammals, the corresponding region in horse shows extraordinary copy number variation and different relative location and directionality of the Eqca-DRB, -DQA, -DQB and -DOB loci. This is the first long-read sequence assembly of the horse MHC class II region with rigorous manual gene annotation, and it will serve as an important resource for association studies of immune-mediated equine diseases and for evolutionary analysis of genetic diversity in this region.

SUBMITTER: Viluma A

PROVIDER: S-EPMC5374520 | biostudies-literature | 2017 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Genomic structure of the horse major histocompatibility complex class II region resolved using PacBio long-read sequencing technology.

Viļuma Agnese A Mikko Sofia S Hahn Daniela D Skow Loren L Andersson Göran G Bergström Tomas F TF

Scientific reports 20170331

The mammalian Major Histocompatibility Complex (MHC) region contains several gene families characterized by highly polymorphic loci with extensive nucleotide diversity, copy number variation of paralogous genes, and long repetitive sequences. This structural complexity has made it difficult to construct a reliable reference sequence of the horse MHC region. In this study, we used long-read single molecule, real-time (SMRT) sequencing technology from Pacific Biosciences (PacBio) to sequence eight ...[more]

PMID: 28361880

Similar Datasets

Project description:BACKGROUND: The quail and chicken major histocompatibility complex (Mhc) genomic regions have a similar overall organization but differ markedly in that the quail has an expanded number of duplicated class I, class IIB, natural killer (NK)-receptor-like, lectin-like and BG genes. Therefore, the elucidation of genetic factors that contribute to the greater Mhc diversity in the quail would help to establish it as a model experimental animal in the investigation of avian Mhc associated diseases. AIMS AND APPROACHES: The main aim here was to characterize the genetic and genomic features of the transcribed major quail MhcIIB (CojaIIB) region that is located between the Tapasin and BRD2 genes, and to compare our findings to the available information for the chicken MhcIIB (BLB). We used four approaches in the study of the quail MhcIIB region, (1) haplotype analyses with polymorphic loci, (2) cloning and sequencing of the RT-PCR CojaIIB products from individuals with different haplotypes, (3) genomic sequencing of the CojaIIB region from the individuals with the different haplotypes, and (4) phylogenetic and duplication analysis to explain the variability of the region between the quail and the chicken. RESULTS: Our results show that the Tapasin-BRD2 segment of the quail Mhc is highly variable in length and in gene transcription intensity and content. Haplotypic sequences were found to vary in length between 4 to 11 kb. Tapasin-BRD2 segments contain one or two major transcribed CojaIIBs that were probably generated by segmental duplications involving c-type lectin-like genes and NK receptor-like genes, gene fusions between two CojaIIBs and transpositions between the major and minor CojaIIB segments. The relative evolutionary speed for generating the MhcIIBs genomic structures from the ancestral BLB2 was estimated to be two times faster in the quail than in the chicken after their separation from a common ancestor. Four types of genomic rearrangement elements (GRE), composed of simple tandem repeats (STR), were identified in the MhcIIB genomic segment located between the Tapasin-BRD2 genes. The GREs have many more STR numbers in the quail than in the chicken that displays strong linkage disequilibrium. CONCLUSION: This study suggests that the Mhc classIIB region has a flexible genomic structure generated by rearrangement elements and rapid SNP accumulation probably as a consequence of the quail adapting to environmental conditions and pathogens during its migratory history after its divergence from the chicken.

Project description:BackgroundWe have previously identified associations between major histocompatibility complex (MHC) class I and resistance towards bacterial and viral pathogens in Atlantic salmon. To evaluate if only MHC or also closely linked genes contributed to the observed resistance we ventured into sequencing of the duplicated MHC class I regions of Atlantic salmon.ResultsNine BACs covering more than 500 kb of the two duplicated MHC class I regions of Atlantic salmon were sequenced and the gene organizations characterized. Both regions contained the proteasome components PSMB8, PSMB9, PSMB9-like and PSMB10 in addition to the transporter for antigen processing TAP2, as well as genes for KIFC1, ZBTB22, DAXX, TAPBP, BRD2, COL11A2, RXRB and SLC39A7. The IA region contained the recently reported MHC class I Sasa-ULA locus residing approximately 50 kb upstream of the major Sasa-UBA locus. The duplicated class IB region contained an MHC class I locus resembling the rainbow trout UCA locus, but although transcribed it was a pseudogene. No other MHC class I-like genes were detected in the two duplicated regions. Two allelic BACs spanning the UBA locus had 99.2% identity over 125 kb, while the IA region showed 82.5% identity over 136 kb to the IB region. The Atlantic salmon IB region had an insert of 220 kb in comparison to the IA region containing three chitin synthase genes.ConclusionWe have characterized the gene organization of more than 500 kb of the two duplicated MHC class I regions in Atlantic salmon. Although Atlantic salmon and rainbow trout are closely related, the gene organization of their IB region has undergone extensive gene rearrangements. The Atlantic salmon has only one class I UCA pseudogene in the IB region while trout contains the four MHC UCA, UDA, UEA and UFA class I loci. The large differences in gene content and most likely function of the salmon and trout class IB region clearly argues that sequencing of salmon will not necessarily provide information relevant for trout and vice versa.

Project description: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.

Dataset Information

Genomic structure of the horse major histocompatibility complex class II region resolved using PacBio long-read sequencing technology.

Publications

Genomic structure of the horse major histocompatibility complex class II region resolved using PacBio long-read sequencing technology.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets