Development of a high-resolution NGS-based HLA-typing and analysis pipeline

The human leukocyte antigen (HLA) complex contains the most polymorphic genes in the human genome. The classical HLA class I and II genes define the specificity of adaptive immune responses. Genetic variation at the HLA genes is associated with susceptibility to autoimmune and infectious diseases and plays a major role in transplantation medicine and immunology. Currently, the HLA genes are characterized using Sanger- or next-generation sequencing (NGS) of a limited amplicon repertoire or labeled oligonucleotides for allele-specific sequences. High-quality NGS-based methods are in proprietary use and not publicly available. Here, we introduce the first highly automated open-kit/open-source HLA-typing method for NGS. The method employs in-solution targeted capturing of the classical class I (HLA-A, HLA-B, HLA-C) and class II HLA genes (HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, HLA-DPB1). The calling algorithm allows for highly confident allele-calling to three-field resolution (cDNA nucleotide variants). The method was validated on 357 commercially available DNA samples with known HLA alleles obtained by classical typing. Our results showed on average an accurate allele call rate of 0.99 in a fully automated manner, identifying also errors in the reference data. Finally, our method provides the flexibility to add further enrichment target regions

Towards the simplification of MHC typing protocols: targeting classical MHC class II genes in a passerine, the pied flycatcher Ficedula hypoleuca

<p>Abstract</p> <p>Background</p> <p>Major Histocompatibility Complex (MHC) has drawn the attention of evolutionary biologists due to its importance in crucial biological processes, such as sexual selection and immune response in jawed vertebrates. However, the characterization of classical MHC genes subjected to the effects of natural selection still remains elusive in many vertebrate groups. Here, we have tested the suitability of flanking intron sequences to guide the selective exploration of classical MHC genes driving the co-evolutionary dynamics between pathogens and their passerine (Aves, Order Passeriformes) hosts.</p> <p>Findings</p> <p>Intronic sequences flanking the usually polymorphic exon 2 were isolated from different species using primers sitting on conserved coding regions of MHC class II genes (β chain). Taking the pied flycatcher <it>Ficedula hypoleuca</it> as an example, we demonstrate that careful primer design can evade non-classical MHC gene and pseudogene amplification. At least four polymorphic and expressed loci were co-replicated using a single pair of primers in five non-related individuals (N = 28 alleles). The cross-amplification and preliminary inspection of similar MHC fragments in eight unrelated songbird taxa suggests that similar approaches can also be applied to other species.</p> <p>Conclusions</p> <p>Intron sequences flanking the usually polymorphic exon 2 may assist the specific investigation of classical MHC class II B genes in species characterized by extensive gene duplication and pseudogenization. Importantly, the evasion of non-classical MHC genes with a more specific function and non-functional pseudogenes may accelerate data collection and diminish lab costs. Comprehensive knowledge of gene structure, polymorphism and expression profiles may be useful not only for the selective examination of evolutionarily relevant genes but also to restrict chimera formation by minimizing the number of co-amplifying loci.</p

Evolutionary genetics of MHC class II beta genes in the brown hare, Lepus europaeus

The genes of the major histocompatibility complex (MHC) are attractive candidates for investigating the link between adaptive variation and individual fitness. High levels of diversity at the MHC are thought to be the result of parasite-mediated selection and there is growing evidence to support this theory. Most studies, however, target just a single gene within the MHC and infer any evidence of selection to be representative of the entire gene region. Here we present data from three MHC class II beta genes (DPB, DQB, and DRB) for brown hares in two geographic regions and compare them against previous results from a class II alpha-chain gene (DQA). We report moderate levels of diversity and high levels of population differentiation in the DQB and DRB genes (Na = 11, Dest = 0.071 and Na = 15, Dest = 0.409, respectively), but not for the DPB gene (Na = 4, Dest = 0.00). We also detected evidence of positive selection within the peptide binding region of the DQB and DRB genes (95% CI, ω > 1.0) but found no signature of selection for DPB. Mutation and recombination were both found to be important processes shaping the evolution of the class II genes. Our findings suggest that while diversifying selection is a significant contributor to the generally high levels of MHC diversity, it does not act in a uniform manner across the entire MHC class II region. The beta-chain genes that we have characterized provide a valuable set of MHC class II markers for future studies of the evolution of adaptive variation in Leporids

High major histocompatibility complex class I polymorphism despite bottlenecks in wild and domesticated populations of the zebra finch ()

Background Two subspecies of zebra finch, Taeniopygia guttata castanotis and T. g. guttata are native to Australia and the Lesser Sunda Islands, respectively. The Australian subspecies has been domesticated and is now an important model system for research. Both the Lesser Sundan subspecies and domesticated Australian zebra finches have undergone population bottlenecks in their history, and previous analyses using neutral markers have reported reduced neutral genetic diversity in these populations. Here we characterize patterns of variation in the third exon of the highly variable major histocompatibility complex (MHC) class I α chain. As a benchmark for neutral divergence, we also report the first mitochondrial NADH dehydrogenase 2 (ND2) sequences in this important model system. Results Despite natural and human-mediated population bottlenecks, we find that high MHC class I polymorphism persists across all populations. As expected, we find higher levels of nucleotide diversity in the MHC locus relative to neutral loci, and strong evidence of positive selection acting on important residues forming the peptide-binding region (PBR). Clear population differentiation of MHC allele frequencies is also evident, and this may be due to adaptation to new habitats and associated pathogens and/or genetic drift. Whereas the MHC Class I locus shows broad haplotype sharing across populations, ND2 is the first locus surveyed to date to show reciprocal monophyly of the two subspecies. Conclusions Despite genetic bottlenecks and genetic drift, all surveyed zebra finch populations have maintained high MHC Class I diversity. The diversity at the MHC Class I locus in the Lesser Sundan subspecies contrasts sharply with the lack of diversity in previously examined neutral loci, and may thus be a result of selection acting to maintain polymorphism. Given uncertainty in historical population demography, however, it is difficult to rule out neutral processes in maintaining the observed diversity. The surveyed populations also differ in MHC Class I allele frequencies, and future studies are needed to assess whether these changes result in functional immune differences

Age before beauty? Relationships between fertilization success and age-dependent ornaments in barn swallows

When males become more ornamented and reproduce more successfully as they grow older, phenotypic correlations between ornament exaggeration and reproductive success can be confounded with age effects in cross-sectional studies, and thus say relatively little about sexual selection on these traits. This is exemplified here in a correlative study of male fertilization success in a large colony of American barn swallows (Hirundo rustica erythrogaster). Previous studies of this species have indicated that two sexually dimorphic traits, tail length and ventral plumage coloration, are positively correlated with male fertilization success, and a mechanism of sexual selection by female choice has been invoked. However, these studies did not control for potential age-related variation in trait expression. Here, we show that male fertilization success was positively correlated with male tail length but not with plumage coloration. We also show that 1-year-old males had shorter tails and lower fertilization success than older males. This age effect accounted for much of the covariance between tail length and fertilization success. Still, there was a positive relationship between tail length and fertilization success among older males. But as this group consisted of males from different age classes, an age effect may be hidden in this relationship as well. Our data also revealed a longitudinal increase in both tail length and fertilization success for individual males. We argue that age-dependent ornament expression and reproductive performance in males complicate inferences about female preferences and sexual selection

Genotyping of black grouse MHC class II B using reference Strand-Mediated Conformational Analysis (RSCA)

<p>Abstract</p> <p>Background</p> <p>The Major Histocompatibility Complex (MHC) is a cluster of genes involved in the vertebrate immune system and includes loci with an extraordinary number of alleles. Due to the complex evolution of MHC genes, alleles from different loci within the same MHC class can be very similar and therefore difficult to assign to separate loci. Consequently, single locus amplification of MHC genes is hard to carry out in species with recently duplicated genes in the same MHC class, and multiple MHC loci have to be genotyped simultaneously. Since amplified alleles have the same length, accurate genotyping is difficult. Reference Strand-Mediated Conformational Analysis (RSCA), which is increasingly used in studies of natural populations with multiple MHC genes, is a genotyping method capable to provide high resolution and accuracy in such cases.</p> <p>Findings</p> <p>We adapted the RSCA method to genotype multiple MHC class II B (BLB) genes in black grouse (<it>Tetrao tetrix</it>), a non-model galliform bird species, using a 96-Capillary Array Electrophoresis, the MegaBACE™ 1000 DNA Analysing System (GE Healthcare). In this study we used fluorescently labelled reference strands from both black grouse and hazel grouse and observed good agreement between RSCA and cloning/sequencing since 71 alleles were observed by cloning/sequencing and 76 alleles by RSCA among the 24 individuals included in the comparison. At the individual level however, there was a trend towards more alleles scored with RSCA (1-6 per individual) than cloning/sequencing (1-4 per individual). In 63% of the pair-wise comparison, the identical allele was scored in RSCA as in cloning/sequencing. Nine out of 24 individuals had the same number of alleles in RSCA as in cloning/sequencing. Our RSCA protocol allows a faster RSCA genotyping than presented in many other RSCA studies.</p> <p>Conclusions</p> <p>In this study, we have developed the RSCA typing method further to work on a 96-Capillary Array Electrophoresis (MegaBACE™ 1000). Our RSCA protocol can be applied to fast and reliable screening of MHC class II B diversity of black grouse populations. This will facilitate future large-scale population studies of black grouse and other galliformes species with multiple inseparable MHC loci.</p