Classical Alleles Research Articles

DMPPred: a tool for identification of antigenic regions responsible for inducing type 1 diabetes mellitus.

There are a number of antigens that induce autoimmune response against β-cells, leading to type 1 diabetes mellitus (T1DM). Recently, several antigen-specific immunotherapies have been developed to treat T1DM. Thus, identification of T1DM associated peptides with antigenic regions or epitopes is important for peptide based-therapeutics (e.g. immunotherapeutic). In this study, for the first time, an attempt has been made to develop a method for predicting, designing, and scanning of T1DM associated peptides with high precision. We analysed 815 T1DM associated peptides and observed that these peptides are not associated with a specific class of HLA alleles. Thus, HLA binder prediction methods are not suitable for predicting T1DM associated peptides. First, we developed a similarity/alignment based method using Basic Local Alignment Search Tool and achieved a high probability of correct hits with poor coverage. Second, we developed an alignment-free method using machine learning techniques and got a maximum AUROC of 0.89 using dipeptide composition. Finally, we developed a hybrid method that combines the strength of both alignment free and alignment-based methods and achieves maximum area under the receiver operating characteristic of 0.95 with Matthew's correlation coefficient of 0.81 on an independent dataset. We developed a web server 'DMPPred' and stand-alone server for predicting, designing and scanning T1DM associated peptides (https://webs.iiitd.edu.in/raghava/dmppred/).

Genome polarisation for detecting barriers to geneflow

Abstract Semi‐permeable barriers to geneflow in principle allow distantly related organisms to capture and exchange pre‐adapted genes potentially speeding adaptation. However, describing barriers to geneflow on a genomic scale is non‐trivial. We extend classic diagnostic allele counting measures of geneflow across a barrier to the case of genome‐scale data. Diagnostic index expectation maximisation (diem) polarises the labelling of bistate markers with respect to the sides of a barrier. An initial state of ignorance is enforced by starting with randomly generated marker polarisations. This means there is no prior on population or taxon membership of the genomes concerned. Using a deterministic data labelling, small numbers of classic diagnostic markers can be replaced by large numbers of markers, each with a diagnostic index. Individuals' hybrid indices (genome admixture proportions) are then calculated genome wide conditioned on marker diagnosticity; within diploid, haplodiploid and/or haploid genome compartments; or indeed over any subset of markers, allowing classical cline width/barrier strength comparisons along genomes. Along‐genome barrier strength heterogeneity allows for barrier regions to be identified. Furthermore, blocks of genetic material that have introgressed across a barrier are easily identified with high power. diem indicates panmixis among Myotis myotis bat genomes, with a barrier separating low data quality outliers. In a Mus musculus domesticus/Mus spretus system, diem adds multiple introgressions of olfactory (and vomeronasal) gene clusters in one direction to previous demonstrations of a pesticide resistance gene introgressing in the opposite direction across a strong species barrier. diem is a genome analysis solution, which scales over reduced representation genomics of thousands of markers to treatment of all variant sites in large genomes. While the method lends itself to visualisation, its output of markers with barrier‐informative annotation will fuel research in population genetics, phylogenetics and association studies. diem can equip such downstream applications with millions of informative markers.

Transgenerational genomic analyses reveal allelic oscillation and purifying selection in a gut parasite Nosema ceranae.

Standing genetic variation is the predominant source acted on by selection. Organisms with high genetic diversity generally show faster responses toward environmental change. Nosema ceranae is a microsporidian parasite of honey bees, infecting midgut epithelial cells. High genetic diversity has been found in this parasite, but the mechanism for the parasite to maintain this diversity remains unclear. This study involved continuous inoculation of N. ceranae to honey bees. We found that the parasites slowly increased genetic diversity over three continuous inoculations. The number of lost single nucleotide variants (SNVs) was balanced with novel SNVs, which were mainly embedded in coding regions. Classic allele frequency oscillation was found at the regional level along the genome, and the associated genes were enriched in apoptosis regulation and ATP binding. The ratio of synonymous and non-synonymous substitution suggests a purifying selection, and our results provide novel insights into the evolutionary dynamics in microsporidian parasites.

Evaluation of Genetic Diversity and Parasite-Mediated Selection of MHC Class I Genes in Emberiza godlewskii (Passeriformes: Emberizidae)

The major histocompatibility complex (MHC) is a multi-copy immune gene family in vertebrates. Its genes are highly variable and code for antigen-presenting molecules. Characterization of MHC genes in different species and investigating the mechanisms that shape MHC diversity is an important goal in understanding the evolution of biological diversity. Here we developed a next-generation sequencing (NGS) protocol to genotype the MHC class I genes of 326 Godlewski’s buntings (Emberiza godlewskii) sampled in the Western mountain area of Beijing from 2014 to 2016. A total of 184 functional alleles were identified, including both non-classical and classical alleles, clustering into nine supertypes. Compared with other passerine birds, the number of MHC class I alleles per individual in Godlewski’s buntings is high (mean 16.1 ± 3.3, median 16). In addition, we demonstrated signatures of historical and contemporary selection on MHC genes. Reflecting historical selection, ten amino acid sites in the antigen-binding domain showed signatures of balancing selection, eight of which exhibit high amino acid polymorphism. In terms of contemporary selection, we found that specific MHC supertypes were nominally associated with the infection of two malaria parasite lineages. These findings indicate the action of historical and possibly also contemporary balancing selection and suggest negative frequency-dependent or fluctuating selection as possible selection mechanisms.

Complement component C4 structural variation and quantitative traits contribute to sex-biased vulnerability in systemic sclerosis

Copy number (CN) polymorphisms of complement C4 play distinct roles in many conditions, including immune-mediated diseases. We investigated the association of C4 CN with systemic sclerosis (SSc) risk. Imputed total C4, C4A, C4B, and HERV-K CN were analyzed in 26,633 individuals and validated in an independent cohort. Our results showed that higher C4 CN confers protection to SSc, and deviations from CN parity of C4A and C4B augmented risk. The protection contributed per copy of C4A and C4B differed by sex. Stronger protection was afforded by C4A in men and by C4B in women. C4 CN correlated well with its gene expression and serum protein levels, and less C4 was detected for both in SSc patients. Conditioned analysis suggests that C4 genetics strongly contributes to the SSc association within the major histocompatibility complex locus and highlights classical alleles and amino acid variants of HLA-DRB1 and HLA-DPB1 as C4-independent signals.

High frequency of CD8 escape mutations in elite controllers as new obstacle for HIV cure

ABSTRACT Accumulation of mutations in epitopes of cytolytic-T-lymphocytes immune response (CTL) in HIV-reservoir seems to be one of the reasons for shock-and-kill strategy failure. Ten non-controller patients on successful cART (TX) and seven elite controllers (EC) were included. HIV-Gag gene from purified resting memory CD4+ T-cells was sequenced by Next-Generation-Sequencing. HLA class-I alleles were typed to predict optimal HIV-Gag CTL epitopes. For each subject, the frequency of mutated epitopes in the HIV-Gag gene, the proportion of them considered as CTL-escape variants as well as their effect on antigen recognition by HLA were assessed. The proportion (%) of mutated HIV-Gag CTL epitopes in the reservoir was high and similar in EC and TX (86%[50–100] and 57%[48–82] respectively, p=0.315). Many of them were predicted to negatively impact antigen recognition. Moreover, the proportion of mutated epitopes considered to be CTL-escape variants was also similar in TX and EC (77%[49–92] vs. 50%[33–75] respectively, p=0.117). Thus, the most relevant finding of our study was the high and similar proportions of HIV-Gag CTL-escape mutations in the reservoir of both HIV-noncontroller patients with cART (TX) and patients with spontaneous HIV-control (EC). Our findings suggest that escape mutations of CTL-response may be another obstacle to eliminate the HIV reservoir and constitute a proof of concept that challenges HIV cure strategies focused on the reactivation of reservoirs. Due to the small sample size that could impact the robustness of the study, further studies with larger cohorts of elite controller patients are needed to confirm these results.

316.11: HLA-DPB1 Molecular Mismatches Are Associated With Acute Rejection and Long-term Graft Function in First Kidney Transplants

Introduction: It is generally accepted that HLA-DP mismatches (MM) do not impact the survival of first transplants but that they may exert a deleterious effect in retransplants. This study aimed to investigate the impact of HLA-DPB1 allele and molecular MM (mMM) on the occurrence of acute rejection (AR) and low 5-year graft function (5-yGF) in first kidney transplants (KT). Method: HLA-DPB1 allele and mMM were determined in 130 first KT from deceased donors, performed between 2014 and 2016. The following mMM models were investigated: expression MM, with the high expression G allele in the donor, T cell epitope (TCE) MM, classified as permissive and nonpermissive MM, epitope MM (EMM), considering all six HLA-DPB1 hypervariable regions (EMM-ABCDEF HVR) or only ABEF regions (EMM-ABEF HVR), Eplet MM (EpMM), considering all or only antibody-verified (AbVer) eplets, and solvent accessible amino acid MM (SAMM). EMM, EpletMM and SAMM loads were categorized using thresholds determined by ROC analysis for each outcome. The outcomes were AR and 5-yGF (estimated by CKD-EPI). Logistic and Cox regressions were used to investigate variables associated with 5-yGF and AR, respectively. Results: 119 (91.5%) KT were performed with one or two HLA-DPB1 allele MM. The overall incidence of AR was 17.7% and the median 5-yGF was 40.3 (27.0 - 51.7) mL/min/1.73m2. No association between HLA-DPB1 allele MM and 5-yGF or rejection was observed. Considering only recipients of standard criteria donor, five logistic models for 5-yGF were built including variables that reached a p-value <0.10 in the univariable analysis (age and donor sex, cause of brain death, ABDR MM and AR) and the mMM (TCE, SAMM, all and AbVer EpMM, EMM-ABCDEF and ABEF HVR). In the model 2, variables associated with 5-yGF≤ 40 were donor age (OR=1.08; CI95%=1.02-1.14; p=0.005) and SAMM≥7 (OR=4.40; CI95%=1.56-12.36, p<0.005), while in the model 4, donor age (OR= 1.09; CI95%=1.03-1.15; p=0.003) and EpMM(AbVer≥2 (OR=3.39; CI95%=1.2-9.45, p=0.02) were associated with low GF. The AUC-ROC for predicting 5-yGF was 0.83 (0.75-0.92) for the model 2, and 0.68 (0.57-0.79) for the model 4. Lastly, for the outcome AR in all recipients, beyond the final donor creatinine, the following mMM models were statistically significant: TCE (HR=3.01; CI95%=1.33-6.84; p=0.008), SAMM≥5 (HR=2.69; CI95%=1.10-6.56; p=0.03), AbVerEpMM ≥3 (HR=3.38; 95%CI=1.32-8.70; p=0.01), and EMM-ABCDEF HVR ≥ 6 (HR=5.08; 95%CI=2.19-11.82; p<0.001). Conclusion: To the best of our knowledge, this is the first study that shows that some mMM, in contrast to classical allele HLA-DPB1 MM, are associated with acute rejection and long-term graft function in a population of exclusively first kidney transplants.

Classical and Non-Classical HLA Alleles as Supplementary Markers in Indirect Kinship Parentage Testing

A civil paternity investigation involving the parents of the deceased alleged father in order to establish a family relationship is presented. On the basis of the 23 autosomal short tandem repeat (aSTR) genotyping results, conclusive proof of paternity was not achieved, as the probability of paternity (W) was calculated to 0.99988. Additional genetic data of 17 classical and non-classical human leukocyte alleles (HLA) typing by next-generation sequencing (NGS) at a high-resolution level supported the hypothesis of grandpaternity over the hypothesis of coincidental paternal obligate allele (POA) sharing (total WaSTR&amp;HLA = 0.9999998). The present study demonstrates the utility of 17 HLA genetic markers-typing in the solution of deficiency cases of disputed parentage.

Snowflake: A deep learning-based human leukocyte antigen matching algorithm considering allele-specific surface accessibility.

Histocompatibility in solid-organ transplantation has a strong impact on long-term graft survival. Although recent advances in matching of both B-cell epitopes and T-cell epitopes have improved understanding of allorecognition, the immunogenic determinants are still not fully understood. We hypothesized that HLA solvent accessibility is allele-specific, thus supporting refinement of HLA B-cell epitope prediction. We developed a computational pipeline named Snowflake to calculate solvent accessibility of HLA Class I proteins for deposited HLA crystal structures, supplemented by constructed HLA structures through the AlphaFold protein folding predictor and peptide binding predictions of the APE-Gen docking framework. This dataset trained a four-layer long short-term memory bidirectional recurrent neural network, which in turn inferred solvent accessibility of all known HLA Class I proteins. We extracted 676 HLA Class-I experimental structures from the Protein Data Bank and supplemented it by 37 Class-I alleles for which structures were predicted. For each of the predicted structures, 10 known binding peptides as reported by the Immune Epitope DataBase were rendered into the binding groove. Although HLA Class I proteins predominantly are folded similarly, we found higher variation in root mean square difference of solvent accessibility between experimental structures of different HLAs compared to structures with identical amino acid sequence, suggesting HLA’s solvent accessible surface is protein specific. Hence, residues may be surface-accessible on e.g. HLA-A*02:01, but not on HLA-A*01:01. Mapping these data to antibody-verified epitopes as defined by the HLA Epitope Registry reveals patterns of (1) consistently accessible residues, (2) only subsets of an epitope’s residues being consistently accessible and (3) varying surface accessibility of residues of epitopes. Our data suggest B-cell epitope definitions can be refined by considering allele-specific solvent-accessibility, rather than aggregating HLA protein surface maps by HLA class or locus. To support studies on epitope analyses in organ transplantation, the calculation of donor-allele-specific solvent-accessible amino acid mismatches was implemented as a cloud-based web service.

Genetic analysis of Caenorhabditis elegans Haspin-like genes shows that hasp-1 plays multiple roles in the germline.

ABSTRACTHaspin is a histone kinase that promotes error-free chromosome segregation by recruiting the chromosomal passenger complex (CPC) to mitotic and meiotic chromosomes. Haspin remains less well studied than other M-phase kinases, and the models explaining Haspin function have been developed primarily in mitotic cells. Here, we generate strains containing new conditional or nonsense mutations in the Caenorhabditis elegans Haspin homologs hasp-1 and hasp-2 and characterize their phenotypes. We show that hasp-1 is responsible for all predicted functions of Haspin and that loss of function of hasp-1 using classical and conditional alleles produces defects in germline stem cell proliferation and spermatogenesis, and confirms its role in oocyte meiosis. Genetic analysis suggests that hasp-1 acts downstream of the Polo-like kinase plk-2 and shows synthetic interactions between hasp-1 and two genes expected to promote recruitment of the CPC by a parallel pathway that depends on the kinase Bub1. This work adds to the growing understanding of Haspin function by characterizing a variety of roles in an intact animal.

Human Leukocyte Antigen Fine-Mapping and Correlation Analysis of Han and Minority Leprosy Patients in Southern China.

Backround: Leprosy is very prevalent in many populations around the world, which is well known that both alleles for human leukocyte antigen (HLA) as well as single nucleotide polymorphisms (SNPs) in the HLA region are common in leprosy patients. Previous studies have identified leprosy-associated susceptibility genes that explain only part of disease risk and heritability. In view of the complicated characteristics of the major histocompatibility complex (MHC) region, this study aimed to explore the development and variation of HLA in leprosy and its possible mechanism. Methods: Previous genome-wide association data were extracted from Han and minority populations in southern China for HLA fine-mapping studies. Insertion and deletion (INDEL), SNP, and copy number variation (CNV) imputation were determined by using the Thousand People Database (1KGP Phase 3 Dataset) as a reference panel. The HAN-MHC database was used to input the HLA classical alleles and amino acids in the MHC region, and further step-regression analysis was performed to analyze independent variation signals associated with leprosy. Results: The most significant locus rs75324027 (the same locus as rs602875 in the HLA-DR region) [p = 7.49E-09, OR= 0.62, 95%,CI: 0.52–0.73] in the intergene region between HLA-DQA1 and HLA-DRB1 was related with leprosy in M-S(Han leprosy patients in south China)disease. In M-SM (Leprosy patients of ethnic minorities in south China)disease, one of the most significant loci of the HLA-DQB1 gene was 6-32626438-A-T (p = 4.49E-08, OR = 0.36, 95%,CI: 0.25–0.52). Therefore, rs75324027 is a locus in M-S disease, and 6-32626438-a-T may be a new locus in M-SM disease. The interaction between 6 and 32626438-A-T and RS75324027 was analyzed, and A significant interaction relationship was found. In the optimal model, the accuracy of prediction was 0.5974, cross-validation Consistency:10, p = 0.0107. Conclusion: In conclusion, this study is the first to assess the association between HLA and leprosy susceptibility in Han and other minority populations in southern China using the Thousand Population database and the Han MHC database. In addition, our analysis validated the previously reported locus rs602875 in the HLA-DR region and for the first time identified an unreported independent locus in leprosy among ethnic minorities in southern China.

Evolution of an Apomixis-Specific Allele Class in Supernumerary Chromatin of Apomictic Boechera.

Asexual reproduction through seeds in plants (i.e., apomixis) is a heritable trait, and apomixis- linked loci have been identified in multiple species. However, direct identification of genomic elements is typically hindered as apomixis-linked loci and are commonly found in recombination-suppressed and repetitive regions. Heterochromatinized elements, such as B chromosomes and other supernumerary chromosomal DNA fragments have long been known to be associated with asexuality in both plants and animals and are prime candidate regions for the evolution of multiple apomixis factors controlling the individual elements of apomixis. Here, we examined molecular evolution, gene regulation, and chromosomal location of a male apomeiosis factor (UPG2), a long noncoding RNA gene, in sexual and apomictic Boechera with and without male apomeiosis (i.e., balanced and unbalanced apomicts). We revealed the origin of the gene in the apomixis genome on an apomixis-specific, supernumerary heterochromatic Boechera chromosome (Boe1). The UPG2 is active in the tapetum at male meiosis. We found allele classes specific to apomictic and sexual Boechera accessions and a third class that shares the features of both and points to a convergent transition state. Sex alleles are found only in some of the sexual accessions and have higher nucleotide divergence and lower transcriptional activity compared to apo alleles. These data demonstrate selective pressure to maintain the function of UPG2 for unreduced pollen formation in apomicts as the occasional transmission of the allele from unbalanced apomicts into sexual organisms that lead to pseudogenization and functional decay of copies in sexual organisms.

MO206: HLA Frequency in A Spanish Transplant Patients Cohort With ANCA-Associated Vasculitis

Abstract BACKGROUND AND AIMS ANCA-associated vasculitis (AAV) is an autoimmune disease that causes end-stage renal disease in ∼20%–25% of patients who develop kidney failure with progression to end-stage kidney disease with requirement of renal replacement therapy. Previous studies have found a genetic association between some major histocompatibility complex (MHC) loci, like the locus DPB1*0401, and AAV. We conducted a retrospective observational unicentric study with the objective to identify other possible loci related to AAV. METHOD As MHC determination is not usually done in patients with AAV, we looked for patients with AAV included in the kidney transplant list with a complete classic HLA allele determination registered (A, B, C, DR, and DQ). Only 14 patients with diagnosis of AAV included in the kidney transplant list from August 1993 to December 2020 in our hospital centre fulfilled this inclusion criterion. To the statistical analysis, we compared the observed proportion to the theoretical one using a normal test, since the groups in our sample included at least five subjects. We also use a continuity correction factor to approach to a binomial test. RESULTS The highest allele frequency in our cohort were the HLA B*51and the haplotype HLA DRB1*04-DQB1*03:02. We identified a B*51 frequency of 42.8% (6/14 patients) in our AAV register, while its population frequency in Spain is 15%, difference = 27.86% [Confidence interval (CI) 1.93%–53.78%] P = .01. We detected a DRB1*04 frequency of 57.1% (8/14 patients) in our cohort, while the population frequency of DRB1*04 and its associated loci DQB1*03:02 is 22%, difference = 35.14% (CI 9.22%–61.07%) P &amp;lt; .05. CONCLUSION HLA B*51 and the haplotype DRB1*04-DQB1*03:02 are alleles more frequent in patients with AAV than general population. This could translate a genetic predisposition of people with these alleles to develop AAV.

Large-scale fungal strain sequencing unravels the molecular diversity in mating loci maintained by long-term balancing selection.

Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of strains. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci define mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty strains of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. The sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. 17,550 mating types were predicted based on the allelic classes. In vitro crosses allowed us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. We conclude that the genetic diversity of mating loci in Trichaptum is due to long-term balancing selection, with limited recombination and duplication activity. The large number of sequenced strains highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Trans-ancestral fine-mapping of MHC reveals key amino acids associated with spontaneous clearance of hepatitis C in HLA-DQβ1

Spontaneous clearance of acute hepatitis C virus (HCV) infection is associated with single nucleotide polymorphisms (SNPs) on the MHC class II. We fine-mapped the MHC region in European (n = 1,600; 594 HCV clearance/1,006 HCV persistence) and African (n = 1,869; 340 HCV clearance/1,529 HCV persistence) ancestry individuals and evaluated HCV peptide binding affinity of classical alleles. In both populations, HLA-DQβ1Leu26 (p valueMeta = 1.24 × 10−14) located in pocket 4 was negatively associated with HCV spontaneous clearance and HLA-DQβ1Pro55 (p valueMeta = 8.23 × 10−11) located in the peptide binding region was positively associated, independently of HLA-DQβ1Leu26. These two amino acids are not in linkage disequilibrium (r2 < 0.1) and explain the SNPs and classical allele associations represented by rs2647011, rs9274711, HLA-DQB1∗03:01, and HLA-DRB1∗01:01. Additionally, HCV persistence classical alleles tagged by HLA-DQβ1Leu26 had fewer HCV binding epitopes and lower predicted binding affinities compared to clearance alleles (geometric mean of combined IC50 nM of persistence versus clearance; 2,321 nM versus 761.7 nM, p value = 1.35 × 10−38). In summary, MHC class II fine-mapping revealed key amino acids in HLA-DQβ1 explaining allelic and SNP associations with HCV outcomes. This mechanistic advance in understanding of natural recovery and immunogenetics of HCV might set the stage for much needed enhancement and design of vaccine to promote spontaneous clearance of HCV infection.

Association with HLA-DRβ1 position 37 distinguishes juvenile dermatomyositis from adult-onset myositis.

Juvenile dermatomyositis (JDM) is a rare, severe autoimmune disease and the most common idiopathic inflammatory myopathy of children. JDM and adult-onset dermatomyositis (DM) have similar clinical, biological and serological features, although these features differ in prevalence between childhood-onset and adult-onset disease, suggesting that age of disease onset may influence pathogenesis. Therefore, a JDM-focused genetic analysis was performed using the largest collection of JDM samples to date. Caucasian JDM samples (n = 952) obtained via international collaboration were genotyped using the Illumina HumanCoreExome chip. Additional non-assayed human leukocyte antigen (HLA) loci and genome-wide single-nucleotide polymorphisms (SNPs) were imputed. HLA-DRB1*03:01 was confirmed as the classical HLA allele most strongly associated with JDM [odds ratio (OR) 1.66; 95% confidence interval (CI) 1.46, 1.89; P = 1.4 × 10−14], with an independent association at HLA-C*02:02 (OR = 1.74; 95% CI 1.42, 2.13, P = 7.13 × 10−8). Analyses of amino acid positions within HLA-DRB1 indicated that the strongest association was at position 37 (omnibus P = 3.3 × 10−19), with suggestive evidence this association was independent of position 74 (omnibus P = 5.1 × 10−5), the position most strongly associated with adult-onset DM. Conditional analyses also suggested that the association at position 37 of HLA-DRB1 was independent of some alleles of the Caucasian HLA 8.1 ancestral haplotype (AH8.1) such as HLA-DQB1*02:01 (OR = 1.62; 95% CI 1.36, 1.93; P = 8.70 × 10−8), but not HLA-DRB1*03:01 (OR = 1.49; 95% CR 1.24, 1.80; P = 2.24 × 10−5). No associations outside the HLA region were identified. Our findings confirm previous associations with AH8.1 and HLA-DRB1*03:01, HLA-C*02:02 and identify a novel association with amino acid position 37 within HLA-DRB1, which may distinguish JDM from adult DM.

A genotype-phenotype correlation matrix for ABCA4 disease based on long-term prognostic outcomes.

BackgroundMore than 1500 variants in the ATP-binding cassette, sub-family A, member 4 (ABCA4), locus underlie a heterogeneous spectrum of retinal disorders ranging from aggressive childhood-onset chorioretinopathy to milder late-onset macular disease. Genotype-phenotype correlation studies have been limited in clinical applicability as patient cohorts are typically small and seldom capture the full natural history of individual genotypes. To overcome these limitations, we constructed a genotype-phenotype correlation matrix that provides quantifiable probabilities of long-term disease outcomes associated with specific ABCA4 genotypes from a large, age-restricted patient cohort.MethodsThe study included 112 unrelated patients at least 50 years of age in whom 2 pathogenic variants were identified after sequencing of the ABCA4 locus. Clinical characterization was performed using the results of best corrected visual acuity, retinal imaging, and full-field electroretinogram testing.ResultsFour distinct prognostic groups were defined according to the spatial severity of disease features across the fundus. Recurring genotypes were observed in milder prognoses, including a newly defined class of rare hypomorphic alleles. PVS1 (predicted null) variants were enriched in the most severe prognoses; however, missense variants were present in a larger-than-expected fraction of these patients. Analysis of allele combinations and their respective prognostic severity showed that certain variants, such as p.(Gly1961Glu), and both rare and frequent hypomorphic alleles, were “clinically dominant” with respect to patient phenotypes irrespective of the allele in trans.ConclusionThese results provide much-needed structure to the complex genetic and clinical landscape of ABCA4 disease and add a tool to the clinical repertoire to quantitatively assess individual genotype-specific prognoses in patients.FUNDINGNational Eye Institute, NIH, grants R01 EY028203, R01 EY028954, R01 EY029315, P30 19007 (Core Grant for Vision Research); the Foundation Fighting Blindness USA, grant no. PPA-1218-0751-COLU; and Research to Prevent Blindness.

Fine-Mapping of the Major Histocompatibility Complex Region Linked to Leprosy in Northern China.

Background: Leprosy is a chronic infectious skin and neurological disease, and genetic background is considered to be one of the major factors of risk. The major histocompatibility complex (MHC) region not only affects susceptibility to leprosy but also its development and outcome. Given the complex traits of the MHC region, variants and the potential mechanism by which HLA influences leprosy development need to be further explored. Methods: We extracted previous genome-wide association study data from the Northern Han Chinese population to perform HLA fine-mapping. Using the 1,000 Genome Project Phase 3 dataset as the reference panel, single-nucleotide polymorphisms (SNP), insertion and deletion (INDEL) and copy number variant (CNV) imputation were carried out. HLA classical alleles and amino acids in the MHC region were imputed using the HAN-MHC database. Further stepwise regression analysis was conducted to analyze independent signals of variants related to leprosy. Results: We identified four independent variants: esv3608598, rs7754498, rs3130781 and rs144388449. Among them, esv3608598 is a CNV and the first HLA CNV associated with leprosy risk. SNP annotation using RegulomeDB, HaploReg, and rVarBase showed that three SNPs are likely to affect the pathogenesis of leprosy. Conclusion: In summary, this is the first study to assess the association between HLA CNV and leprosy susceptibility in a Northern Han Chinese population. By fine mapping of the MHC region in this population, our findings provide evidence for the contribution of HLA to leprosy susceptibility.

Structure of the translating Neurospora ribosome arrested by cycloheximide

Ribosomes translate RNA into proteins. The protein synthesis inhibitor cycloheximide (CHX) is widely used to inhibit eukaryotic ribosomes engaged in translation elongation. However, the lack of structural data for actively translating polyribosomes stalled by CHX leaves unanswered the question of which elongation step is inhibited. We elucidated CHX's mechanism of action based on the cryo-electron microscopy structure of actively translating Neurospora crassa ribosomes bound with CHX at 2.7-Å resolution. The ribosome structure from this filamentous fungus contains clearly resolved ribosomal protein eL28, like higher eukaryotes but unlike budding yeast, which lacks eL28. Despite some differences in overall structures, the ribosomes from Neurospora, yeast, and humans all contain a highly conserved CHX binding site. We also sequenced classic Neurospora CHX-resistant alleles. These mutations, including one at a residue not previously observed to affect CHX resistance in eukaryotes, were in the large subunit proteins uL15 and eL42 that are part of the CHX-binding pocket. In addition to A-site transfer RNA (tRNA), P-site tRNA, messenger RNA, and CHX that are associated with the translating N. crassa ribosome, spermidine is present near the CHX binding site close to the E site on the large subunit. The tRNAs in the peptidyl transferase center are in the A/A site and the P/P site. The nascent peptide is attached to the A-site tRNA and not to the P-site tRNA. The structural and functional data obtained show that CHX arrests the ribosome in the classical PRE translocation state and does not interfere with A-site reactivity.

Ssl2/TFIIH function in transcription start site scanning by RNA polymerase II in Saccharomyces cerevisiae.

In Saccharomyces cerevisiae, RNA polymerase II (Pol II) selects transcription start sites (TSSs) by a unidirectional scanning process. During scanning, a preinitiation complex (PIC) assembled at an upstream core promoter initiates at select positions within a window ~40-120 bp downstream. Several lines of evidence indicate that Ssl2, the yeast homolog of XPB and an essential and conserved subunit of the general transcription factor (GTF) TFIIH, drives scanning through its DNA-dependent ATPase activity, therefore potentially controlling both scanning rate and scanning extent (processivity). To address questions of how Ssl2 functions in promoter scanning and interacts with other initiation activities, we leveraged distinct initiation-sensitive reporters to identify novel ssl2 alleles. These ssl2 alleles, many of which alter residues conserved from yeast to human, confer either upstream or downstream TSS shifts at the model promoter ADH1 and genome-wide. Specifically, tested ssl2 alleles alter TSS selection by increasing or narrowing the distribution of TSSs used at individual promoters. Genetic interactions of ssl2 alleles with other initiation factors are consistent with ssl2 allele classes functioning through increasing or decreasing scanning processivity but not necessarily scanning rate. These alleles underpin a residue interaction network that likely modulates Ssl2 activity and TFIIH function in promoter scanning. We propose that the outcome of promoter scanning is determined by two functional networks, the first being Pol II activity and factors that modulate it to determine initiation efficiency within a scanning window, and the second being Ssl2/TFIIH and factors that modulate scanning processivity to determine the width of the scanning widow.