Challenges for the standardized reporting of NGS HLA genotyping: Surveying gaps between clinical and research laboratories.

Next generation sequencing (NGS) is being applied for HLA typing in research and clinical settings. NGS HLA typing has made it feasible to sequence exons, introns and untranslated regions simultaneously, with significantly reduced labor and reagent cost per sample, rapid turnaround time, and improved HLA genotype accuracy. NGS technologies bring challenges for cost-effective computation, data processing and exchange of NGS-based HLA data. To address these challenges, guidelines and specifications such as Genotype List (GL) String, Minimum Information for Reporting Immunogenomic NGS Genotyping (MIRING), and Histoimmunogenetics Markup Language (HML) were proposed to streamline and standardize reporting of HLA genotypes. As part of the 17th International HLA and Immunogenetics Workshop (IHIW), we implemented standards and systems for HLA genotype reporting that included GL String, MIRING and HML, and found that misunderstanding or misinterpretations of these standards led to inconsistencies in the reporting of NGS HLA genotyping results. This may be due in part to a historical lack of centralized data reporting standards in the histocompatibility and immunogenetics community. We have worked with software and database developers, clinicians and scientists to address these issues in a collaborative fashion as part of the Data Standard Hackathons (DaSH) for NGS. Here we report several categories of challenges to the consistent exchange of NGS HLA genotyping data we have observed. We hope to address these challenges in future DaSH for NGS efforts.

Frequencies and haplotype associations of non-expressed HLA alleles in ethnically diverse populations on the National Marrow Donor Program's Be The Match Registry.

HLA allele matching is critical to successful bone marrow transplantation between a patient and donor. Non-functional HLA alleles, so called 'null alleles', are not well described within a large population of well HLA-typed ethnically diverse individuals despite their impact on donor selection. A retrospective analysis was performed on 833,789 unrelated donors (URDs) in the National Marrow Donor Program's Be The Match Registry® typed for HLA-A, -B, -C, -DRB1, -DQB1, and -DPB1 by next-generation DNA sequencing. Results showed that null alleles occur in low frequency (2.30E-04) compared to expressed alleles. Their overall frequency ranged from 6.00E-07 to 9.25E-04 with a median of 1.20E-06. The expected allele associations were commonly observed for HLA-A*24:09N, HLA-B*51:11N, and HLA-C*04:09N; however, associations outside of the expected were also observed. Notably, 82% of the National Marrow Donor Program Registry URDs carrying HLA-A*24:11N showed a different HLA-C allele association, HLA-C*05:01, compared to the allele described by prior published work characterizing German donor populations, HLA-C*04:01. The frequencies of these observed null alleles and linkage disequilibrium information could be invaluable and helpful in guiding the HLA testing decisions.

Aflatoxicosis chemoprevention by probiotic Lactobacillius and lack of effect on the major histocompatibility complex.

Turkeys are extremely sensitive to aflatoxin B1 (AFB1) which causes decreased growth, immunosuppression and liver necrosis. The purpose of this study was to determine whether probiotic Lactobacillus, shown to be protective in animal and clinical studies, would likewise confer protection in turkeys, which were treated for 11 days with either AFB1 (AFB; 1 ppm in diet), probiotic (PB; 1 × 10(11) CFU/ml; oral, daily), probiotic + AFB1 (PBAFB), or PBS control (CNTL). The AFB1 induced drop in body and liver weights were restored to normal in CNTL and PBAFB groups. Hepatotoxicity markers were not significantly reduced by probiotic treatment. Major histocompatibility complex (MHC) genes BG1 and BG4, which are differentially expressed in liver and spleens, were not significantly affected by treatments. These data indicate modest protection, but the relatively high dietary AFB1 treatment, and the extreme sensitivity of this species may reveal limits of probiotic-based protection strategies.

Genetic variation at the MHC in a population of introduced wild turkeys.

Genetic variation in the major histocompatibility complex (MHC) is known to affect disease resistance in many species. Investigations of MHC diversity in populations of wild species have focused on the antigen presenting class IIß molecules due to the known polymorphic nature of these genes and the role these molecules play in pathogen recognition. Studies of MHC haplotype variation in the turkey ( Meleagris gallopavo ) are limited. This study was designed to examine MHC diversity in a group of Eastern wild turkeys ( Meleagris gallopavo silvestris ) collected during population expansion following reintroduction of the species in southern Wisconsin, USA. Southern blotting with BG and class IIß probes and single nucleotide polymorphism (SNP) genotyping was used to measure MHC variation. SNP analysis focused on single copy MHC genes flanking the highly polymorphic class IIß genes. Southern blotting identified 27 class IIß phenotypes, whereas SNP analysis identified 13 SNP haplotypes occurring in 28 combined genotypes. Results show that genetic diversity estimates based on RFLP (Southern blot) analysis underestimate the level of variation detected by SNP analysis. Sequence analysis of the mitochondrial D-loop identified 7 mitochondrial haplotypes (mitotypes) in the sampled birds. Results show that wild turkeys located in southern Wisconsin have a genetically diverse MHC and originate from several maternal lineages.

Defining the turkey MHC: identification of expressed class I- and class IIB-like genes independent of the MHC-B.

The MHC of the turkey (Meleagris gallopavo) is divided into two genetically unlinked regions; the MHC-B and MHC-Y. Although previous studies found the turkey MHC-B to be highly similar to that of the chicken, little is known of the gene content and extent of the MHC-Y. This study describes two partially overlapping large-insert BAC clones that genetically and physically map to the turkey MHC chromosome (MGA18) but to a region that assorts independently of MHC-B. Within the sequence assembly, 14 genes were predicted including new class I- and class IIB-like loci. Additional unassembled sequences corresponded to multiple copies of the ribosomal RNA repeat unit (18S-5.8S-28S). Thus, this newly identified MHC region appears to represent a physical boundary of the turkey MHC-Y. High-resolution multi-color fluorescence in situ hybridization studies confirm rearrangement of MGA18 relative to the orthologous chicken chromosome (GGA16) in regard to chromosome architecture, but not gene order. The difference in centromere position between the species is indicative of multiple chromosome rearrangements or alternate events such as neocentromere formation/centromere inactivation in the evolution of the MHC chromosome. Comparative sequencing of commercial turkeys (six amplicons totaling 7.6 kb) identified 68 single nucleotide variants defining nine MHC-Y haplotypes. Sequences of the new class I- and class IIB-like genes are most similar to MHC-Y genes in the chicken. All three loci are expressed in the spleen. Differential transcription of the MHC-Y class IIB-like loci was evident as one class IIB-like locus was only expressed in some individuals.

Extended sequence of the turkey MHC B-locus and sequence variation in the highly polymorphic B-G loci.

Genetic variation in the major histocompatibility complex (MHC) is directly correlated to differences in disease resistance. Immunity is greatly dependent on highly polymorphic genes in the MHC, such as class I, class II, and class III complement genes. Preliminary studies of wild turkey populations show extreme polymorphisms in a family of genes exclusive to the avian MHC, the class IV or B-G genes. Significance of this variation is unclear as there are few and conflicting studies of the expression of these genes. Confounding understanding of B-G variation is the lack of a complete delineation of the number of loci in the turkey genome. Direct 454 sequencing of a clone from the CHORI-260 BAC library was used to extend the turkey MHC B-locus sequence, identifying five additional complete B-locus genes including two B-G loci. Sequences of the new B-G genes were compared with those of other turkey gene (BG1-3) and sequences available for other galliformes. Phylogenetic analysis shows species-specific gene evolution supporting a birth-death model of evolution for the B-G gene family. Analysis of variation within the signal peptide sequence (exon 1) found two clusters of polymorphism among the turkey B-G genes. Resequencing of exon 1 in a diverse sample including wild, heritage, and commercial turkeys confirmed multiple alleles at each B-G gene. Future studies aim to correlate B-G variation with group and individual immunological differences.

A candidate gene for choanal atresia in alpaca.

Choanal atresia (CA) is a common nasal craniofacial malformation in New World domestic camelids (alpaca and llama). CA results from abnormal development of the nasal passages and is especially debilitating to newborn crias. CA in camelids shares many of the clinical manifestations of a similar condition in humans (CHARGE syndrome). Herein we report on the regulatory gene CHD7 of alpaca, whose homologue in humans is most frequently associated with CHARGE. Sequence of the CHD7 coding region was obtained from a non-affected cria. The complete coding region was 9003 bp, corresponding to a translated amino acid sequence of 3000 aa. Additional genomic sequences corresponding to a significant portion of the CHD7 gene were identified and assembled from the 2x alpaca whole genome sequence, providing confirmatory sequence for much of the CHD7 coding region. The alpaca CHD7 mRNA sequence was 97.9% similar to the human sequence, with the greatest sequence difference being an insertion in exon 38 that results in a polyalanine repeat (A12). Polymorphism in this repeat was tested for association with CA in alpaca by cloning and sequencing the repeat from both affected and non-affected individuals. Variation in length of the poly-A repeat was not associated with CA. Complete sequencing of the CHD7 gene will be necessary to determine whether other mutations in CHD7 are the cause of CA in camelids.

Comparative genomics identifies new alpha class genes within the avian glutathione S-transferase gene cluster.

Glutathione S-transferases (GSTs: EC2.5.1.18) are a superfamily of multifunctional dimeric enzymes that catalyze the conjugation of glutathione (GSH) to electrophilic chemicals. In most animals and in humans, GSTs are the principal enzymes responsible for detoxifying the mycotoxin aflatoxin B(1) (AFB(1)) and GST dysfunction is a known risk factor for susceptibility towards AFB(1). Turkeys are one of the most susceptible animals known to AFB(1), which is a common contaminant of poultry feeds. The extreme susceptibility of turkeys is associated with hepatic GSTs unable to detoxify the highly reactive and electrophilic metabolite exo-AFB(1)-8,9-epoxide (AFBO). In this study, comparative genomic approaches were used to amplify and identify the alpha-class tGST genes (tGSTA1.1, tGSTA1.2, tGSTA1.3, tGSTA2, tGSTA3 and tGSTA4) from turkey liver. The conserved GST domains and four alpha-class signature motifs in turkey GSTs (with the exception of tGSTA1.1 which lacked one motif) confirm the presence of hepatic alpha-class GSTs in the turkey. Four signature motifs and conserved residues found in alpha-class tGSTs are (1) xMExxxWLLAAAGVE, (2) YGKDxKERAxIDMYVxG, (3) PVxEKVLKxHGxxxL and (4) PxIKKFLXPGSxxKPxxx. A BAC clone containing the alpha-class GST gene cluster was isolated and sequenced. The turkey alpha-class GTS genes genetically map to chromosome MGA2 with synteny between turkey and human alpha-class GSTs and flanking genes. This study identifies the alpha-class tGST gene cluster and genetic markers (SNPs, single nucleotide polymorphisms) that can be used to further examine AFB(1) susceptibility and resistance in turkeys. Functional characterization of heterologously expressed proteins from these genes is currently underway.