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

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

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

HLA Class I and Genetic Susceptibility to Type 1 Diabetes: Results From the Type 1 Diabetes Genetics Consortium

OBJECTIVE-We report here genotyping data and type 1 diabetes association analyses for HLA class I loci (A, B, and C) on 1,753 multiplex pedigrees from the Type 1 Diabetes Genetics Consortium (T1DGC), a large international collaborative study. RESEARCH DESIGN AND METHODS-Complete eight-locus HLA genotyping data were generated. Expected patient class I (HLA-A, -B, and -C) allele frequencies were calculated, based on linkage disequilibrium (LD) patterns with observed HLA class II DRB1-DQA1-DQB1 haplotype frequencies. Expected frequencies were compared to observed allele frequencies in patients. RESULTS-Significant type 1 diabetes associations were observed at all class I HLA loci. After accounting for LD with HLA class II, the most significantly type 1 diabetes-associated alleles were B*5701 (odds ratio 0.19; P = 4 x 10(-11)) and B*3906 (10.31; P = 4 X 10(-10)). Other significantly type 1 diabetes-associated alleles included A*2402, A*0201, B*1801, and C*0501 (predisposing) and A*1101, A*3201, A*6601, B*0702, B*4403, B*3502, C*1601, and C*0401 (protective). Some alleles, notably B*3906, appear to modulate the risk of all DRB1-DQA1-DQB1 haplotypes on which they reside, suggesting a class I effect that is independent of class H. Other class I type 1 diabetes associations appear to be specific to individual class H haplotypes. Some apparent associations (e.g., C*1601) could be attributed to strong LD to another class I susceptibility locus (B*4403). CONCLUSIONS-These data indicate that HLA class I alleles, in addition to and independently from HLA class H alleles, are associated with type 1 diabetes. Diabetes 59:2972-2979, 201

Evaluation of the rheumatoid arthritis susceptibility loci HLA-DRB1, PTPN22, OLIG3/TNFAIP3, STAT4 and TRAF1/C5 in an inception cohort.

INTRODUCTION: This study investigated five confirmed rheumatoid arthritis (RA) susceptibility genes/loci (HLA-DRB1, PTPN22, STAT4, OLIG3/TNFAIP3 and TRAF1/C5) for association with susceptibility and severity in an inception cohort. METHODS: The magnitude of association for each genotype was assessed in 1,046 RA subjects from the Yorkshire Early RA cohort and in 5,968 healthy UK controls. Additional exploratory subanalyses were undertaken in subgroups defined by autoantibody status (rheumatoid factor and anti-cyclic citrullinated peptide) or disease severity (baseline articular erosions, Health Assessment Questionnaire (HAQ) score and swollen joint count (SJC)). RESULTS: In the total RA inception cohort, the HLA-DRB1 shared epitope (per-allele odds ratio (OR) = 2.1, trend P < 0.0001), PTPN22 (per-allele OR = 1.5, trend P < 0.0001), OLIG3/TNFAIP3 locus (per-allele OR = 1.2, trend P = 0.009) and TRAF1/C5 locus (per-allele OR = 1.1, trend P = 0.04) were associated with RA. The magnitude of association for these loci was increased in those patients who were autoantibody-positive. PTPN22 was associated with autoantibody-negative RA (per-allele OR = 1.3, trend P = 0.04). There was no evidence of association between these five genetic loci and baseline erosions or SJC in the total RA cohort, after adjustment for symptom duration. TRAF1/C5 was significantly associated with baseline HAQ, however, following adjustment for symptom duration (P trend = 0.03). CONCLUSIONS: These findings support the mounting evidence that different genetic loci are associated with autoantibody-positive and autoantibody-negative RA, possibly suggesting that many of the genes identified to date are associated with autoantibody production. Additional studies with a specific focus on autoantibody-negative RA will be needed to identify the genes predisposing to this RA subgroup. The TRAF1/C5 locus in particular warrants further investigation in RA as a potential disease severity locus

Tests for Genetic Interactions in Type 1 Diabetes: Linkage and Stratification Analyses of 4,422 Affected Sib-Pairs

OBJECTIVE - Interactions between genetic and environmental factors lead to immune dysregulation causing type 1 diabetes and other autoimmune disorders. Recently, many common genetic variants have been associated with type 1 diabetes risk, but each has modest individual effects. Familial clustering of type 1 diabetes has not been explained fully and could arise from many factors, including undetected genetic variation and gene interactions. RESEARCH DESIGN AND METHODS - To address this issue, the Type 1 Diabetes Genetics Consortium recruited 3,892 families, including 4,422 affected sib-pairs. After genotyping 6,090 markers, linkage analyses of these families were performed, using a novel method and taking into account factors such as genotype at known susceptibility loci. RESULTS - Evidence for linkage was robust at the HLA and INS loci, with logarithm of odds (LOD) scores of 398.6 and 5.5, respectively. There was suggestive support for five other loci. Stratification by other risk factors (including HLA and age at diagnosis) identified one convincing region on chromosome 6q14 showing linkage in male subjects (corrected LOD = 4.49; replication P = 0.0002), a locus on chromosome 19q in HLA identical siblings (replication P = 0.006), and four other suggestive loci. CONCLUSIONS - This is the largest linkage study reported for any disease. Our data indicate there are no major type 1 diabetes subtypes definable by linkage analyses; susceptibility is caused by actions of HLA and an apparently random selection from a large number of modest-effect loci; and apart from HLA and INS, there is no important susceptibility factor discoverable by linkage methods

HLA DPA1, DPB1 Alleles and Haplotypes Contribute to the Risk Associated With Type 1 Diabetes: Analysis of the Type 1 Diabetes Genetics Consortium Families

OBJECTIVE-To determine the relative risk associated with DPA1 and DPB1 alleles and haplotypes in type I diabetes. RESEARCH DESIGN AND METHODS-The frequency of DPA1 and DPB1 alleles and haplotypes in type I diabetic patients was compared to the family based control frequency in 1,771 families directly and conditional on FILA (B)-DRB1-DQA1-DQB1 linkage disequilibrium. A relative predispositional analysis (RPA) was performed in the presence or absence of the primary HLA DR-DQ associations and the contribution of DP haplotype to individual DR-DQ haplotype risks examined. RESULTS-Eight DPAI and thirty-eight DPB1 alleles forming seventy-four DPA1-DPB1 haplotypes were observed, nineteen DPB1 alleles were associated with multiple DPA1 alleles Following both analyses, type I diabetes susceptibility was significantly associated with DPB1*0301 (DPA1*0103-DPB1*0301) and protection with DPB1*0402 (DPA1*0103-DPB1*0402) and DPA1*0103-DPB1*0101 but not DPA1*0201-DPB1*0101. In addition, DPB1*0202 (DPA1*0103-DPB1*0202) and DPB1*0201 (DPA1*0103-DPB1*0201) were significantly associated with susceptibility in the presence of the high risk and protective DR-DQ haplotypes Three associations (DPB1*0301, *0402, and *0202) remained statistically significant when only the extended HLA-A1-B8-DR3 haplotype was considered, suggesting that DPB1 alone may delineate the risk associated with this otherwise conserved haplotype CONCLUSIONS-HLA DP allelic and haplotypic diversity contributes significantly to the risk for type I diabetes; DPB1*0301 (DPA1*0103-DPB1*0301) is associated with susceptibility and DPB1*0402 (DPA1*0103-DPB1*0402) and DPA1*0103-DPB1*0101 with protection Additional evidence is presented for the susceptibility association of DPB1*0202 (DPA1*0103-DPB1*0202) and for a contributory role of individual amino acids and DPA1 or a gene in linkage disequilibrium in DR3-DPB1*0101 positive haplotypes Diabetes 59:2055-2062, 201