Non-invasive prenatal testing of beta-hemoglobinopathies using next generation sequencing, in-silico sequence size selection, and haplotyping.

AIM: To develop a non-invasive prenatal test for beta-hemoglobinopathies based on analyzing maternal plasma by using next generation sequencing. METHODS: We applied next generation sequencing (NGS) of maternal plasma to the non-invasive prenatal testing (NIPT) of autosomal recessive diseases, sickle cell disease and beta-thalassemia. Using the Illumina MiSeq, we sequenced plasma libraries obtained via a Twist Bioscience probe capture panel covering 4 Kb of chromosome 11, including the beta-globin (HBB) gene and >450 genomic single-nucleotide polymorphisms (SNPs) used to estimate the fetal fraction (FF). The FF is estimated by counting paternally transmitted allelic sequence reads present in the plasma but absent in the mother. We inferred fetal beta-globin genotypes by comparing the observed mutation (Mut) and reference (Ref) read ratios to those expected for the three possible fetal genotypes (Mut/Mut; Mut/Ref; Ref/Ref), based on the FF. RESULTS: We bioinformatically enriched the FF by excluding reads over a specified length via in-silico size selection (ISS), favoring the shorter fetal reads, which increased fetal genotype prediction accuracy. Finally, we determined the parental HBB haplotypes, which allowed us to use the read ratios observed at linked SNPs to help predict the fetal genotype at the mutation site(s). We determined HBB haplotypes via Oxford Nanopore MinION sequencing of a 2.2 kb amplicon and aligned these sequences using Soft Genetics' NextGENe LR software. CONCLUSION: The combined use of ISS and HBB haplotypes enabled us to correctly predict fetal genotypes in cases where the prediction based on variant read ratios alone was incorrect.

High resolution HLA analysis reveals independent class I haplotypes and amino-acid motifs protective for multiple sclerosis.

We investigated association between HLA class I and class II alleles and haplotypes, and KIR loci and their HLA class I ligands, with multiple sclerosis (MS) in 412 European American MS patients and 419 ethnically matched controls, using next-generation sequencing. The DRB1*15:01~DQB1*06:02 haplotype was highly predisposing (odds ratio (OR) = 3.98; 95% confidence interval (CI) = 3-5.31; p-value (p) = 2.22E-16), as was DRB1*03:01~DQB1*02:01 (OR = 1.63; CI = 1.19-2.24; p = 1.41E-03). Hardy-Weinberg (HW) analysis in MS patients revealed a significant DRB1*03:01~DQB1*02:01 homozyote excess (15 observed; 8.6 expected; p = 0.016). The OR for this genotype (5.27; CI = 1.47-28.52; p = 0.0036) suggests a recessive MS risk model. Controls displayed no HW deviations. The C*03:04~B*40:01 haplotype (OR = 0.27; CI = 0.14-0.51; p = 6.76E-06) was highly protective for MS, especially in haplotypes with A*02:01 (OR = 0.15; CI = 0.04-0.45; p = 6.51E-05). By itself, A*02:01 is moderately protective, (OR = 0.69; CI = 0.54-0.87; p = 1.46E-03), and haplotypes of A*02:01 with the HLA-B Thr80 Bw4 variant (Bw4T) more so (OR = 0.53; CI = 0.35-0.78; p = 7.55E-04). Protective associations with the Bw4 KIR ligand resulted from linkage disequilibrium (LD) with DRB1*15:01, but the Bw4T variant was protective (OR = 0.64; CI = 0.49-0.82; p = 3.37-04) independent of LD with DRB1*15:01. The Bw4I variant was not associated with MS. Overall, we find specific class I HLA polymorphisms to be protective for MS, independent of the strong predisposition conferred by DRB1*15:01.

Correction: High resolution HLA analysis reveals independent class I haplotypes and amino-acid motifs protective for multiple sclerosis.

Since the publication of this article, the authors have found that the numbers of patients and controls were reversed. This study included 412 MS patients and 419 controls. This correction applies to the Abstract, the final paragraph of the Introduction, and the first paragraph of the Materials and Methods. This was entirely a reporting error and does not impact the Results or Conclusions.

MHC class II super-enhancer increases surface expression of HLA-DR and HLA-DQ and affects cytokine production in autoimmune vitiligo.

Genetic risk for autoimmunity in HLA genes is most often attributed to structural specificity resulting in presentation of self-antigens. Autoimmune vitiligo is strongly associated with the MHC class II region. Here, we fine-map vitiligo MHC class II genetic risk to three SNPs only 47 bp apart, located within a predicted super-enhancer in an intergenic region between HLA-DRB1 and HLA-DQA1, localized by a genome-wide association study of 2,853 Caucasian vitiligo patients. The super-enhancer corresponds to an expression quantitative trait locus for expression of HLA-DR and HLA-DQ RNA; we observed elevated surface expression of HLA-DR (P = 0.008) and HLA-DQ (P = 0.02) on monocytes from healthy subjects homozygous for the high-risk SNP haplotype. Unexpectedly, pathogen-stimulated peripheral blood mononuclear cells from subjects homozygous for the high-risk super-enhancer haplotype exhibited greater increase in production of IFN-γ and IL-1ß than cells from subjects homozygous for the low-risk haplotype. Specifically, production of IFN-γ on stimulation of dectin-1, mannose, and Toll-like receptors with Candida albicans and Staphylococcus epidermidis was 2.5- and 2.9-fold higher in high-risk subjects than in low-risk subjects, respectively (P = 0.007 and P = 0.01). Similarly, production of IL-1ß was fivefold higher in high-risk subjects than in low-risk subjects (P = 0.02). Increased production of immunostimulatory cytokines in subjects carrying the high-risk haplotype may act as an "adjuvant" during the presentation of autoantigens, tying together genetic variation in the MHC with the development of autoimmunity. This study demonstrates that for risk of autoimmune vitiligo, expression level of HLA class II molecules is as or more important than antigen specificity.

Association of Parkinson disease with structural and regulatory variants in the HLA region.

Historically, association of disease with the major histocompatibility complex (HLA) genes has been tested with HLA alleles that encode antigen-binding affinity. The association with Parkinson disease (PD), however, was discovered with noncoding SNPs in a genome-wide association study (GWAS). We show here that several HLA-region SNPs that have since been associated with PD form two blocks tagged by rs3129882 (p = 9 × 10(-11)) and by rs9268515 and/or rs2395163 (p = 3 × 10(-11)). We investigated whether these SNP-associations were driven by HLA-alleles at adjacent loci. We imputed class I and class II HLA-alleles for 2000 PD cases and 1986 controls from the NeuroGenetics Research Consortium GWAS and sequenced a subset of 194 cases and 204 controls. We were therefore able to assess accuracy of two imputation algorithms against next-generation-sequencing while taking advantage of the larger imputed data sets for disease study. Additionally, we imputed HLA alleles for 843 cases and 856 controls from another GWAS for replication. PD risk was positively associated with the B(∗)07:02_C(∗)07:02_DRB5(∗)01_DRB1(∗)15:01_DQA1(∗)01:02_DQB1(∗)06:02 haplotype and negatively associated with the C(∗)03:04, DRB1(∗)04:04 and DQA1(∗)03:01 alleles. The risk haplotype and DQA1(∗)03:01 lost significance when conditioned on the SNPs, but C(∗)03:04 (OR = 0.72, p = 8 × 10(-6)) and DRB1(∗)04:04 (OR = 0.65, p = 4 × 10(-5)) remained significant. Similarly, rs3129882 and the closely linked rs9268515 and rs2395163 remained significant irrespective of HLA alleles. rs3129882 and rs2395163 are expression quantitative trait loci (eQTLs) for HLA-DR and HLA-DQ (9 × 10(-5) ≥ PeQTL ≥ 2 × 10(-79)), suggesting that HLA gene expression might influence PD. Our data suggest that PD is associated with both structural and regulatory elements in HLA. Furthermore, our study demonstrates that noncoding SNPs in the HLA region can be associated with disease irrespective of HLA alleles, and that observed associations with HLA alleles can sometimes be secondary to a noncoding variant.

Human NK cells licensed by killer Ig receptor genes have an altered cytokine program that modifies CD4+ T cell function.

NK cells are innate immune cells known for their cytolytic activities toward tumors and infections. They are capable of expressing diverse killer Ig-like receptors (KIRs), and KIRs are implicated in susceptibility to Crohn's disease (CD), a chronic intestinal inflammatory disease. However, the cellular mechanism of this genetic contribution is unknown. In this study, we show that the "licensing" of NK cells, determined by the presence of KIR2DL3 and homozygous HLA-C1 in host genome, results in their cytokine reprogramming, which permits them to promote CD4(+) T cell activation and Th17 differentiation ex vivo. Microfluidic analysis of thousands of NK single cells and bulk secretions established that licensed NK cells are more polarized to proinflammatory cytokine production than unlicensed NK cells, including production of IFN-γ, TNF-α, CCL-5, and MIP-1ß. Cytokines produced by licensed NK augmented CD4(+) T cell proliferation and IL-17A/IL-22 production. Ab blocking indicated a primary role for IFN-γ, TNF-α, and IL-6 in the augmented T cell-proliferative response. In conclusion, NK licensing mediated by KIR2DL2/3 and HLA-C1 elicits a novel NK cytokine program that activates and induces proinflammatory CD4(+) T cells, thereby providing a potential biologic mechanism for KIR-associated susceptibility to CD and other chronic inflammatory diseases.

A sequence-based approach demonstrates that balancing selection in classical human leukocyte antigen (HLA) loci is asymmetric.

Balancing selection has maintained human leukocyte antigen (HLA) allele diversity, but it is unclear whether this selection is symmetric (all heterozygotes are comparable and all homozygotes are comparable in terms of fitness) or asymmetric (distinct heterozygote genotypes display greater fitness than others). We tested the hypothesis that HLA is under asymmetric balancing selection in populations by estimating allelic branch lengths from genetic sequence data encoding peptide-binding domains. Significant deviations indicated changes in the ratio of terminal to internal branch lengths. Such deviations could arise even if no individual alleles present a strikingly altered branch length (e.g. if there is an overall distortion, with all or many terminal branches being longer than expected). DQ and DP loci were also analyzed as haplotypes. Using allele frequencies for 419 distinct populations in 10 geographical regions, we examined population differentiation in alleles within and between regions, and the relationship between allelic branch length and frequency. The strongest evidence for asymmetrical balancing selection was observed for HLA-DRB1, HLA-B and HLA-DPA1, with significant deviation (P ≤ 1.1 × 10(-4)) in about half of the populations. There were significant results at all loci except HLA-DQB1/DQA1. We observed moderate genetic variation within and between geographic regions, similar to the rest of the genome. Branch length was not correlated with allele frequency. In conclusion, sequence data suggest that balancing selection in HLA is asymmetric (some heterozygotes enjoy greater fitness than others). Because HLA polymorphism is crucial for pathogen resistance, this may manifest as a frequency-dependent selection with fluctuation in the fitness of specific heterozygotes over time.

HLA-DPB1*04:01 Protects Genetically Susceptible Children from Celiac Disease Autoimmunity in the TEDDY Study.

OBJECTIVES: Tissue transglutaminase autoantibodies (tTGAs) represent the first evidence of celiac disease (CD) development. Associations of HLA-DR3-DQA1*05:01-DQB1*02:01 (i.e., DR3-DQ2) and, to a lesser extent, DR4-DQA1*03:01-DQB1*03:02 (i.e., DR4-DQ8) with the risk of CD differ by country, consistent with additional genetic heterogeneity that further refines risk. Therefore, we examined human leukocyte antigen (HLA) factors other than DR3-DQ2 for their contribution to developing tTGAs. METHODS: The Environmental Determinants of Diabetes in the Young (TEDDY) study enrolled 8,676 infants at an increased HLA-DR-DQ risk for type 1 diabetes and CD into a 15-year prospective surveillance follow-up. Of those followed up, 21% (n=1,813) carried DR3-DQ2/DR3-DQ2, 39% (n=3,359) carried DR3-DQ2/DR4-DQ8, 20% (n=1701) carried DR4-DQ8/DR4-DQ8, and 17% (n=1,493) carried DR4-DQ8/DQ4. Within TEDDY, a nested case-control design of 248 children with CD autoimmunity (CDA) and 248 matched control children were genotyped for HLA-B, -DRB3, -DRB4, -DPA1, and -DPB1 genes, and the entire cohort was genotyped for single-nucleotide polymorphisms (SNPs) using the Illumina ImmunoChip. CDA was defined as a positive tTGA test at two consecutive clinic visits, whereas matching in those with no evidence of tTGAs was based on the presence of HLA-DQ2, country, and sex. RESULTS: After adjustment for DR3-DQ2 and restriction to allele frequency (AF) ≥5%, HLA-DPB1*04:01 was inversely associated with CDA by conditional logistic regression (AF=44%, odds ratio=0.71, 95% confidence interval (CI)=0.53-0.96, P=0.025). This association of time to CDA and HLA-DPB1*04:01 was replicated with statistical significance in the remainder of the cohort using imputation for specific HLA alleles based on SNP genotyping (hazard ratio=0.84, 95% CI=0.73-0.96, P=0.013). CONCLUSIONS: HLA-DPB1*04:01 may reduce the risk of tTGAs, an early marker of CD, among DR3-DQ2 children, confirming that additional variants in the HLA region influence the risk for CDA.

Analysis of mixtures using next generation sequencing of mitochondrial DNA hypervariable regions.

AIM: To apply massively parallel and clonal sequencing (next generation sequencing or NGS) to the analysis of forensic mixed samples. METHODS: A duplex polymerase chain reaction (PCR) assay targeting the mitochondrial DNA (mtDNA) hypervariable regions I/II (HVI/HVII) was developed for NGS analysis on the Roche 454 GS Junior instrument. Eight sets of multiplex identifier-tagged 454 fusion primers were used in a combinatorial approach for amplification and deep sequencing of up to 64 samples in parallel. RESULTS: This assay was shown to be highly sensitive for sequencing limited DNA amounts (~100 mtDNA copies) and analyzing contrived and biological mixtures with low level variants (~1%) as well as "complex" mixtures (≥3 contributors). PCR artifact "hybrid" sequences generated by jumping PCR or template switching were observed at a low level (<2%) in the analysis of mixed samples but could be eliminated by reducing the PCR cycle number. CONCLUSION: This study demonstrates the power of NGS technologies targeting the mtDNA HVI/HVII regions for analysis of challenging forensic samples, such as mixtures and specimens with limited DNA.

HLA-DP genetic variation, proxies for early life immune modulation and childhood acute lymphoblastic leukemia risk.

The human leukocyte antigen (HLA) genes are candidate genetic susceptibility loci for childhood acute lymphoblastic leukemia (ALL). We examined the effect of HLA-DP genetic variation on risk and evaluated its potential interaction with 4 proxies for early immune modulation, including measures of infectious exposures in infancy (presence of older siblings, daycare attendance, ear infections) and breastfeeding. A total of 585 ALL cases and 848 controls were genotyped at the HLA-DPA1 and DPB1 loci. Because of potential heterogeneity in effect by race/ethnicity, we included only non-Hispanic white (47%) and Hispanic (53%) children and considered these 2 groups separately in the analysis. Logistic regression analyses showed an increased risk of ALL associated with HLA-DPB1*01:01 (odds ratio [OR] = 1.43, 95% CI, 1.01-2.04) with no heterogeneity by Hispanic ethnicity (P = .969). Analyses of DPB1 supertypes showed a marked childhood ALL association with DP1, particularly for high-hyperdiploid ALL (OR = 1.83; 95% CI, 1.20-2.78). Evidence of interaction was found between DP1 and older sibling (P = .036), and between DP1 and breastfeeding (P = .094), with both showing statistically significant DP1 associations within the lower exposure categories only. These findings support an immune mechanism in the etiology of childhood ALL involving the HLA-DPB1 gene in the context of an insufficiently modulated immune system.

Development and characterization of the cobas human papillomavirus test.

The cobas human papillomavirus (HPV) test, approved by the FDA in April 2011, is a fully automated assay for the detection of 14 high-risk (hr) HPV genotypes from cervical specimens collected in liquid-based cytology medium using real-time PCR amplification of the L1 gene and TaqMan probes. Results are simultaneously reported as positive or negative for the pooled 12 oncogenic HPV types (HPV31, -33, -35, -39, -45, -51, -52, -56, -58, -59, -66, and -68) from channel 1, with HPV16 and HPV18 genotypes read individually from channels 2 and 3. A fourth channel detects the human ß-globin gene as a control for sample adequacy and assay inhibition. To optimize clinical sensitivity and specificity, cutoff values (cycle thresholds [C(T)]) were established for each channel based on the detection of cervical intraepithelial neoplasia grade 2 (CIN2) or greater (≥CIN2). For women aged ≥21 years with cytology results indicating atypical squamous cells of undetermined significance (ASC-US), CT values provided a sensitivity of 90% (95% confidence interval [CI], 81.5% to 94.8%) for the detection of ≥CIN2 and a specificity of 70.5% (95% CI, 68.1% to 72.7%). The analytic sensitivity (limit of detection) ranged from 150 to 2,400 copies/ml, depending on genotype. The analytic specificity, evaluated by comparing the HPV result with a combined comparator of Sanger sequencing and the Qiagen digene HC2 high-risk HPV DNA test (hc2), demonstrated overall positive agreement of 96.3% for 14 hrHPV types in women with ASC-US cytology results who were aged ≥21 years and 86.1% in women with NLIM (negative for intraepithelial neoplasia or malignancy) cytology who were aged ≥30 years. These and other performance validation studies demonstrate that the cobas HPV test is a fully automated and clinically validated robust test.

HLA-C dysregulation as a possible mechanism of immune evasion in SARS-CoV-2 and other RNA-virus infections.

One of the mechanisms by which viruses can evade the host's immune system is to modify the host's DNA methylation pattern. This work aims to investigate the DNA methylation and gene expression profile of COVID-19 patients, divided into symptomatic and asymptomatic, and healthy controls, focusing on genes involved in the immune response. In this study, changes in the methylome of COVID-19 patients' upper airways cells, the first barrier against respiratory infections and the first cells presenting viral antigens, are shown for the first time. Our results showed alterations in the methylation pattern of genes encoding proteins implicated in the response against pathogens, in particular the HLA-C gene, also important for the T-cell mediated memory response. HLA-C expression significantly decreases in COVID-19 patients, especially in those with a more severe prognosis and without other possibly confounding co-morbidities. Moreover, our bionformatic analysis revealed that the identified methylation alteration overlaps with enhancers regulating HLA-C expression, suggesting an additional mechanism exploited by SARS-CoV-2 to inhibit this fundamental player in the host's immune response. HLA-C could therefore represent both a prognostic marker and an excellent therapeutic target, also suggesting a preventive intervention that conjugate a virus-specific antigenic stimulation with an adjuvant increasing the T-cell mediated memory response.

Noninvasive Prenatal Test for ß-Thalassemia and Sickle Cell Disease Using Probe Capture Enrichment and Next-Generation Sequencing of DNA in Maternal Plasma.

BACKGROUND: Noninvasive prenatal testing (NIPT) of chromosomal aneuploidies based on next-generation sequencing (NGS) analysis of fetal DNA in maternal plasma is well established, but testing for autosomal recessive disorders remains challenging. NGS libraries prepared by probe capture facilitate the analysis of the short DNA fragments plasma. This system has been applied to the ß-hemoglobinopathies to reduce the risk to the fetus. METHOD: Our probe panel captures >4 kb of the HBB region and 435 single-nucleotide polymorphisms (SNPs) used to estimate fetal fraction. Contrived mixtures of DNA samples, plasma, and whole blood samples from 7 pregnant women with ß-thalassemia or sickle cell anemia mutations and samples from the father, sibling, and baby or chorionic villus were analyzed. The fetal genotypes, including point mutations and deletions, were inferred by comparing the observed and expected plasma sequence read ratios, based on fetal fraction, at the mutation site and linked SNPs. Accuracy was increased by removing PCR duplicates and by in silico size selection of plasma sequence reads. A probability was assigned to each of the potential fetal genotypes using a statistical model for the experimental variation, and thresholds were established for assigning clinical status. RESULTS: Using in silico size selection of plasma sequence files, the predicted clinical fetal genotype assignments were correct in 9 of 10 plasma libraries with maternal point mutations, with 1 inconclusive result. For 2 additional plasmas with deletions, the most probable fetal genotype was correct. The ß-globin haplotype determined from linked SNPs, when available, was used to infer the fetal genotype at the mutation site. CONCLUSION: This probe capture NGS assay demonstrates the potential of NIPT for ß-hemoglobinopathies.

Interaction between maternal killer immunoglobulin-like receptors and offspring HLAs and susceptibility of childhood ALL.

Acute lymphoblastic leukemia (ALL) in children is associated with a distinct neonatal cytokine profile. The basis of this neonatal immune phenotype is unknown but potentially related to maternal-fetal immune receptor interactions. We conducted a case-control study of 226 case child-mother pairs and 404 control child-mother pairs to evaluate the role of interaction between HLA genotypes in the offspring and maternal killer immunoglobulin-like receptor (KIR) genotypes in the etiology of childhood ALL, while considering potential mediation by neonatal cytokines and the immune-modulating enzyme arginase-II (ARG-II). We observed different associations between offspring HLA-maternal KIR activating profiles and the risk of ALL in different predicted genetic ancestry groups. For instance, in Latino subjects who experience the highest risk of childhood leukemia, activating profiles were significantly associated with a lower risk of childhood ALL (odds ratio [OR] = 0.59; 95% confidence interval [CI], 0.49-0.71) and a higher level of ARG-II at birth (coefficient = 0.13; 95% CI, 0.04-0.22). HLA-KIR activating profiles were also associated with a lower risk of ALL in non-Latino Asians (OR = 0.63; 95% CI, 0.38-1.01), although they had a lower tumor necrosis factor-α level (coefficient = -0.27; 95% CI, -0.49 to -0.06). Among non-Latino White subjects, no significant association was observed between offspring HLA-maternal KIR interaction and ALL risk or cytokine levels. The current study reports the association between offspring HLA-maternal KIR interaction and the development of childhood ALL with variation by predicted genetic ancestry. We also observed some associations between activating profiles and immune factors related to cytokine control; however, cytokines did not demonstrate causal mediation of the activating profiles on ALL risk.

The Environmental Determinants of Diabetes in the Young (TEDDY): genetic criteria and international diabetes risk screening of 421 000 infants.

AIMS: The Environmental Determinants of Diabetes in the Young (TEDDY) study seeks to identify environmental factors influencing the development of type 1 diabetes (T1D) using intensive follow-up of children at elevated genetic risk. This study requires a cost-effective yet accurate screening strategy to identify the high-risk cohort. METHODS: The TEDDY cohort was identified through newborn screening using human leukocyte antigen (HLA) class II genes based on criteria established with pre-TEDDY data. HLA typing was completed at six international centers using different genotyping methods that can achieve >98% accuracy. RESULTS: TEDDY developed separate inclusion criteria for the general population (GP) and first-degree relatives (FDRs) of T1D patients. The FDR eligibility includes nine haplogenotypes (DR3/4, DR4/4, DR4/8, DR3/3, DR4/4b, DR4/1, DR4/13, DR4/9, and DR3/9) for broad HLA diversity, whereas the GP eligibility includes only the first four haplogenotypes with DRB1*0403 as an exclusion allele. TEDDY has screened 414 714 GP infants, of which 19 906 (4.8%) were eligible, whereas 1415 of the 6333 screened FDR infants (22.2%) were eligible. High-resolution confirmation testing of the eligible subjects indicated that the low-cost and low-resolution genotyping techniques employed at the screening centers yielded an accuracy of 99%. There were considerable variations in eligibility rates among the centers for GP (3.5-7.4%) and FDR (19-32%) subjects. The eligibility rates among US ethnic groups were 0.9, 1.3, 5.0, and 6.9% for Asians, Black, Caucasians, and Hispanics, respectively. CONCLUSIONS: Different low-cost and low-resolution genotyping methods are useful for the efficient and accurate identification of a high-risk cohort for follow-up based on the TEDDY HLA inclusion criteria.

Juvenile idiopathic arthritis and HLA class I and class II interactions and age-at-onset effects.

OBJECTIVE: The aim of this study was to quantitate risk and to examine heterogeneity for HLA at high resolution in patients with the most common subtypes of juvenile idiopathic arthritis (JIA), IgM rheumatoid factor-negative polyarticular JIA and oligoarticular JIA. Use of 4-digit comprehensive HLA typing enabled great precision, and a large cohort allowed for consideration of both age at disease onset and disease subtype. METHODS: Polymerase chain reaction-based high-resolution HLA typing for class I and class II loci was accomplished for 820 patients with JIA and 273 control subjects. Specific HLA epitopes, potential interactions of alleles at specific loci and between loci (accounting for linkage disequilibrium and haplotypic associations), and an assessment of the current International League of Associations for Rheumatology classification criteria were considered. RESULTS: An HLA-DRB1/DQB1 effect was shown to be exclusively attributable to DRB1 and was similar between patients with oligoarticular JIA and a younger subgroup of patients with polyarticular JIA. Furthermore, patients with polyarticular JIA showed age-specific related effects, with disease susceptibility in the group older than age 6 years limited to an effect of the HLA-DRB1*08 haplotype, which is markedly different from the additional susceptibility haplotypes, HLA-DRB1*1103/1104, found in the group with oligoarticular JIA and the group of younger patients with polyarticular JIA. Also in contrast to findings for oligoarticular JIA, patients with polyarticular arthritis had no evidence of an HLA class I effect. Markers associated with a reduced risk of disease included DRB1*1501, DRB1*0401, and DRB1*0701. DRB1*1501 was shown to reduce risk across the whole cohort, whereas DRB1*0401 and DRB1*0701 were protective for selected JIA subtypes. Surprisingly, the disease predisposition mediated by DPB1*0201 in individuals without any disease-predisposing DRB1 alleles was great enough to overcome even the very strong protective effect observed for DRB1*1501. CONCLUSION: Inherited HLA factors in JIA show similarities overall as well as differences between JIA subtypes.

Advantage of rare HLA supertype in HIV disease progression.

The highly polymorphic human leukocyte antigen (HLA) class I molecules help to determine the specificity and repertoire of the immune response. The great diversity of these antigen-binding molecules confers differential advantages in responding to pathogens, but presents a major obstacle to distinguishing HLA allele-specific effects. HLA class I supertypes provide a functional classification for the many different HLA alleles that overlap in their peptide-binding specificities. We analyzed the association of these discrete HLA supertypes with HIV disease progression rates in a population of HIV-infected men. We found that HLA supertypes alone and in combination conferred a strong differential advantage in responding to HIV infection, independent of the contribution of single HLA alleles that associate with progression of the disease. The correlation of the frequency of the HLA supertypes with viral load suggests that HIV adapts to the most frequent alleles in the population, providing a selective advantage for those individuals who express rare alleles.

Resolution of mitochondrial DNA mixtures using a probe capture next generation sequencing system and phylogenetic-based software.

Interpreting mixtures with nuclear genetic markers remains one of the persisting challenges in forensic DNA analysis, particularly when the DNA is degraded or present in trace amounts. In these scenarios, analyzing mitochondrial (mt) DNA can be useful due to the higher copy number per cell compared to nuclear DNA. However, until the emergence of Next-Generation Sequencing (NGS) with its clonal sequencing capability, analysis of mtDNA mixtures was very challenging. A mitochondrial genome probe capture Next-Generation Sequencing (NGS) system was used to sequence complex mtDNA mixtures and two different software programs to analyze the sequence data. Analysis of contrived mixtures of two contributors in 50:50 and 95:5 ratios as well as three-person mixtures ranging from near equal proportions (~33:33:33 ratio) to low amounts of the minor contributors (e.g., a 90:5:5 ratio) is reported. This system was also applied to the analysis of mtDNA mixtures from forensically relevant samples. For data analysis, both the variant frequency-based software program GeneMarker®HTS and the phylogenetic-based software program Mixemt was used to de-convolute the mixtures. Using the massively parallel, clonal features of NGS, one can bioinformatically separate and count the individual sequence reads to calculate the proportions of individual contributors using phylogenetically informative polymorphisms. GeneMarker®HTS allows us to detect all mutations, including "private" mutations (non-phylogenetically informative polymorphisms) and assign them to individual contributors based on the frequency of the sequence reads, provided that the proportions of the various contributors are sufficiently different. Using a probe capture NGS system and both GeneMarker®HTS and Mixemt software programs, the interpretation of complex mixtures of equal proportion contributors, trace amount contributors, and more than two contributors in contrived mixtures, as well as interpretation of challenging forensic specimens is demonstrated.

HLA genotyping in the international Type 1 Diabetes Genetics Consortium.

BACKGROUND: Although human leukocyte antigen (HLA) DQ and DR loci appear to confer the strongest genetic risk for type 1 diabetes, more detailed information is required for other loci within the HLA region to understand causality and stratify additional risk factors. The Type 1 Diabetes Genetics Consortium (T1DGC) study design included high-resolution genotyping of HLA-A, B, C, DRB1, DQ, and DP loci in all affected sibling pair and trio families, and cases and controls, recruited from four networks worldwide, for analysis with clinical phenotypes and immunological markers. PURPOSE: In this article, we present the operational strategy of training, classification, reporting, and quality control of HLA genotyping in four laboratories on three continents over nearly 5 years. METHODS: Methods to standardize HLA genotyping at eight loci included: central training and initial certification testing; the use of uniform reagents, protocols, instrumentation, and software versions; an automated data transfer; and the use of standardized nomenclature and allele databases. We implemented a rigorous and consistent quality control process, reinforced by repeated workshops, yearly meetings, and telephone conferences. RESULTS: A total of 15,246 samples have been HLA genotyped at eight loci to four-digit resolution; an additional 6797 samples have been HLA genotyped at two loci. The genotyping repeat rate decreased significantly over time, with an estimated unresolved Mendelian inconsistency rate of 0.21%. Annual quality control exercises tested 2192 genotypes (4384 alleles) and achieved 99.82% intra-laboratory and 99.68% inter-laboratory concordances. LIMITATIONS: The chosen genotyping platform was unable to distinguish many allele combinations, which would require further multiple stepwise testing to resolve. For these combinations, a standard allele assignment was agreed upon, allowing further analysis if required. CONCLUSIONS: High-resolution HLA genotyping can be performed in multiple laboratories using standard equipment, reagents, protocols, software, and communication to produce consistent and reproducible data with minimal systematic error. Many of the strategies used in this study are generally applicable to other large multi-center studies.

A meta-analysis of candidate gene polymorphisms and ischemic stroke in 6 study populations: association of lymphotoxin-alpha in nonhypertensive patients.

BACKGROUND AND PURPOSE: Ischemic stroke is a multifactorial disease with a strong genetic component. Pathways, including lipid metabolism, systemic chronic inflammation, coagulation, blood pressure regulation, and cellular adhesion, have been implicated in stroke pathophysiology, and candidate gene polymorphisms in these pathways have been proposed as genetic risk factors. METHODS: We genotyped 105 simple deletions and single nucleotide polymorphisms from 64 candidate genes in 3550 patients and 6560 control subjects from 6 case-control association studies conducted in the United States, Europe, and China. Genotyping was performed using the same immobilized probe typing system and meta-analyses were based on summary logistic regressions for each study. The primary analyses were fixed-effects meta-analyses adjusting for age and sex with additive, dominant, and recessive models of inheritance. RESULTS: Although 7 polymorphisms showed a nominal additive association, none remained statistically significant after adjustment for multiple comparisons. In contrast, after stratification for hypertension, 2 lymphotoxin-alpha polymorphisms, which are in strong linkage disequilibrium, were significantly associated among nonhypertensive individuals: LTA 252A>G (additive model; OR, 1.41 with 95% CI, 1.20 to 1.65; P=0.00002) and LTA 26Thr>Asn (OR, 1.19 with 95% CI, 1.06 to 1.34; P=0.003). LTA 252A>G remained significant after adjustment for multiple testing using either the false discovery rate or by permutation testing. The 2 single nucleotide polymorphisms showed no association in hypertensive subjects (eg, LTA 252A>G, OR, 0.93; 95% CI, 0.84 to 1.03; P=0.17). CONCLUSIONS: These observations may indicate an important role of LTA-mediated inflammatory processes in the pathogenesis of ischemic stroke.