Mitochondrial Dysfunction and Fatigue in Sjögren's Disease.

Objectives: Sjögren's disease (SjD) is a common exocrine disorder typified by chronic inflammation and dryness, but also profound fatigue, suggesting a pathological basis in cellular bioenergetics. In healthy states, damaged or dysfunctional mitochondrial components are broken down and recycled by mitophagy, a specialized form of autophagy. In many autoimmune disorders, however, evidence suggests that dysfunctional mitophagy allows poorly functioning mitochondria to persist and contribute to a cellular milieu with elevated reactive oxygen species. We hypothesized that mitophagic processes are dysregulated in SjD and that dysfunctional mitochondria contribute to overall fatigue. We sought to link fatigue with mitochondrial dysfunction directly in SjD, heretofore unexamined, and further sought to assess the pathogenic extent and implications of dysregulated mitophagy in SjD. Methods: We isolated pan T cells via negative selection from the peripheral blood mononuclear cells of 17 SjD and 8 age-matched healthy subjects, all of whom completed fatigue questionnaires prior to phlebotomy. Isolated T cells were analyzed for mitochondrial oxygen consumption rate (OCR) and glycolysis using Seahorse, and linear correlations with fatigue measures were assessed. A mitophagy transcriptional signature in SjD was identified by reanalysis of whole-blood microarray data from 190 SjD and 32 healthy subjects. Differential expression analyses were performed by case/control and subgroup analyses comparing SjD patients by mitophagy transcriptional cluster against healthy subjects followed by bioinformatic interpretation using gene set enrichment analysis. Results: Basal OCR, ATP-linked respiration, maximal respiration, and reserve capacity were significantly lower in SjD compared to healthy subjects with no observed differences in non-mitochondrial respiration, basal glycolysis, or glycolytic stress. SjD lymphocytic mitochondria show structural alterations compared to healthy subjects. Fatigue scores related to pain/discomfort in SjD correlated with the altered OCR. Results from subgroup analyses by mitophagic SjD clusters revealed highly variable inter-cluster differentially expressed genes (DEGs) and expanded the number of SjD-associated gene targets by tenfold within the same dataset. Conclusion: Mitochondrial dysfunction, associated with fatigue, is a significant problem in SjD and warrants further investigation.

Genomic Characterization and Phylogenetic Analysis of Linezolid-Resistant Enterococcus from the Nostrils of Healthy Hosts Identifies Zoonotic Transmission.

Linezolid resistance in Enterococcus spp. is increasingly considered critically important and a public health threat which mandates the need to understand their genomic contents and dissemination patterns. Here, we used whole-genome sequencing to characterize the resistome, virulome and mobile genetic elements of nine linezolid-resistant (LZDR) enterococci (seven optrA-E. faecalis, one poxtA-E. faecium and one optrA-E. casseliflavus) previously obtained from the nares of healthy dogs, pigs, pig farmers and tracheal samples of nestling storks in Spain. Also, the relatedness of the isolates with publicly available genomes was accessed by core-genome single nucleotide polymorphism (SNP) analysis. The optrA gene of the E. faecalis and E. casseliflavus isolates was located downstream of the fexA gene. The optrA gene in the E. casseliflavus isolate was carried in a plasmid (pURX4962), while those in the seven E. faecalis isolates were chromosomally located. The OptrA proteins were mostly variants of wild type (DP-2: Y176D/T481P; RDK: I104R/Y176D/E256K; DD-3: Y176D/G393D; and EDD: K3E/Y176D/G393D), except two that were wild type (one E. faecalis and one E. casseliflavus). The poxtA gene in the E. faecium isolate was found alone within its contig. The cfrD was upstream of ermB gene in the E. casseliflavus isolate and flanked by ISNCY and IS1216. All the LZDR enterococci carried plasmid rep genes (2-3) containing tetracycline, chloramphenicol and aminoglycoside resistance genes. All isolates except E. casseliflavus carried at least one intact prophage, of which E. faecalis-ST330 (X4957) from a pig carried the highest (n = 5). Tn6260 was associated with lnuG in E. faecalis-ST330 while Tn554 was with fexA in E. feaecalis-ST59 isolates. All except E. casseliflavus (n = 0) carried at least two metal resistance genes (MRGs), of which poxtA-carrying E. faecium-ST1739 isolate contained the most (arsA, copA, fief, ziaA, znuA, zosA, zupT, and zur). SNP-based analyses identified closely related optrA-E. faecalis isolates from a pig and a pig farmer on the same farm (SNP = 4). Moreover, optrA- carrying E. faecalis-ST32, -ST59, and -ST474 isolates from pigs were related to those previously described from humans (sick and healthy) and cattle in Spain, Belgium, and Switzerland (SNP range 43-86). These findings strongly suggest the transmission of LZDR-E. faecalis between a pig and a pig farmer and potential inter-country dissemination. These highlight the need to strengthen molecular surveillance of LZDR enterococci in all ecological niches and body parts to direct appropriate control strategies.

Extended haplotype with rs41524547-G defines the ancestral origin of SCA10.

Spinocerebellar ataxia type 10 (SCA10) is a rare autosomal dominant ataxia caused by a large expansion of the (ATTCT)n repeat in ATXN10. SCA10 was described in Native American and Asian individuals which prompted a search for an expanded haplotype to confirm a common ancestral origin for the expansion event. All patients with SCA10 expansions in our cohort share a single haplotype defined at the 5'-end by the minor allele of rs41524547, located ~35 kb upstream of the SCA10 expansion. Intriguingly, rs41524547 is located within the miRNA gene, MIR4762, within its DROSHA cleavage site and just outside the seed sequence for mir4792-5p. The world-wide frequency of rs41524547-G is less than 5% and found almost exclusively in the Americas and East Asia-a geographic distribution that mirrors reported SCA10 cases. We identified rs41524547-G(+) DNA from the 1000 Genomes/International Genome Sample Resource and our own general population samples and identified SCA10 repeat expansions in up to 25% of these samples. The reduced penetrance of these SCA10 expansions may be explained by a young (pre-onset) age at sample collection, a small repeat size, purity of repeat units, or the disruption of miR4762-5p function. We conclude that rs41524547-G is the most robust at-risk SNP allele for SCA10, is useful for screening of SCA10 expansions in population genetics studies and provides the most compelling evidence to date for a single, prehistoric origin of SCA10 expansions sometime prior to or during the migration of individuals across the Bering Land Bridge into the Americas.

Towards facilitated interpretation of shotgun metagenomics long-read sequencing data analyzed with KMA for the detection of bacterial pathogens and their antimicrobial resistance genes.

Metagenomic sequencing is a promising method that has the potential to revolutionize the world of pathogen detection and antimicrobial resistance (AMR) surveillance in food-producing environments. However, the analysis of the huge amount of data obtained requires performant bioinformatics tools and databases, with intuitive and straightforward interpretation. In this study, based on long-read metagenomics data of chicken fecal samples with a spike-in mock community, we proposed confidence levels for taxonomic identification and AMR gene detection, with interpretation guidelines, to help with the analysis of the output data generated by KMA, a popular k-mer read alignment tool. Additionally, we demonstrated that the completeness and diversity of the genomes present in the reference databases are key parameters for accurate and easy interpretation of the sequencing data. Finally, we explored whether KMA, in a two-step procedure, can be used to link the detected AMR genes to their bacterial host chromosome, both detected within the same long-reads. The confidence levels were successfully tested on 28 metagenomics datasets which were obtained with sequencing of real and spiked samples from fecal (chicken, pig, and buffalo) or food (minced beef and food enzyme products) origin. The methodology proposed in this study will facilitate the analysis of metagenomics sequencing datasets for KMA users. Ultimately, this will contribute to improvements in the rapid diagnosis and surveillance of pathogens and AMR genes in food-producing environments, as prioritized by the EU.

The Sjögren's Working Group: The 2023 OMERACT meeting and provisional domain generation.

Sjögren's disease (SjD) is a systemic autoimmune exocrinopathy with key features of dryness, pain, and fatigue. SjD can affect any organ system with a variety of presentations across individuals. This heterogeneity is one of the major barriers for developing effective disease modifying treatments. Defining core disease domains comprising both specific clinical features and incorporating the patient experience is a critical first step to define this complex disease. The OMERACT SjD Working Group held its first international collaborative hybrid meeting in 2023, applying the OMERACT 2.2 filter toward identification of core domains. We accomplished our first goal, a scoping literature review that was presented at the Special Interest Group held in May 2023. Building on the domains identified in the scoping review, we uniquely deployed multidisciplinary experts as part of our collaborative team to generate a provisional domain list that captures SjD heterogeneity.

Exposure to air pollution as an environmental determinant of how Sjögren's disease is expressed at diagnosis.

OBJECTIVES: To analyse how the potential exposure to air pollutants can influence the key components at the time of diagnosis of Sjögren's phenotype (epidemiological profile, sicca symptoms, and systemic disease). METHODS: For the present study, the following variables were selected for harmonization and refinement: age, sex, country, fulfilment of 2002/2016 criteria items, dry eyes, dry mouth, and overall ESSDAI score. Air pollution indexes per country were defined according to the OECD (1990-2021), including emission data of nitrogen and sulphur oxides (NO/SO), particulate matter (PM2.5 and 1.0), carbon monoxide (CO) and volatile organic compounds (VOC) calculated per unit of GDP, Kg per 1000 USD. RESULTS: The results of the chi-square tests of independence for each air pollutant with the frequency of dry eyes at diagnosis showed that, except for one, all variables exhibited p-values <0.0001. The most pronounced disparities emerged in the dry eye prevalence among individuals inhabiting countries with the highest NO/SO exposure, a surge of 4.61 percentage points compared to other countries, followed by CO (3.59 points), non-methane (3.32 points), PM2.5 (3.30 points), and PM1.0 (1.60 points) exposures. Concerning dry mouth, individuals residing in countries with worse NO/SO exposures exhibited a heightened frequency of dry mouth by 2.05 percentage points (p<0.0001), followed by non-methane exposure (1.21 percentage points increase, p=0.007). Individuals inhabiting countries with the worst NO/SO, CO, and PM2.5 pollution levels had a higher mean global ESSDAI score than those in lower-risk nations (all p-values <0.0001). When systemic disease was stratified according to DAS into low, moderate, and high systemic activity levels, a heightened proportion of individuals manifesting moderate/severe systemic activity was observed in countries with worse exposures to NO/SO, CO, and PM2.5 pollutant levels. CONCLUSIONS: For the first time, we suggest that pollution levels could influence how SjD appears at diagnosis in a large international cohort of patients. The most notable relationships were found between symptoms (dryness and general body symptoms) and NO/SO, CO, and PM2.5 levels.

Influence of exposure to climate-related hazards in the phenotypic expression of primary Sjögren's syndrome.

OBJECTIVES: To analyse how the key components at the time of diagnosis of the Sjögren's phenotype (epidemiological profile, sicca symptoms, and systemic disease) can be influenced by the potential exposure to climate-related natural hazards. METHODS: For the present study, the following variables were selected for harmonisation and refinement: age, sex, country, fulfilment of 2002/2016 criteria items, dry eyes, dry mouth, and overall ESSDAI score. Climate-related hazards per country were defined according to the OECD and included seven climate-related hazard types: extreme temperature, extreme precipitation, drought, wildfire, wind threats, river flooding, and coastal flooding. Climatic variables were defined as dichotomous variables according to whether each country is ranked among the ten countries with the most significant exposure. RESULTS: After applying data-cleaning techniques and excluding people from countries not included in the OECD climate rankings, the database study analysed 16,042 patients from 23 countries. The disease was diagnosed between 1 and 3 years earlier in people living in countries included among the top 10 worst exposed to extreme precipitation, wildfire, wind threats, river flooding, and coastal flooding. A lower frequency of dry eyes was observed in people living in countries exposed to wind threats, river flooding, and coastal flooding, with a level of statistical association being classified as strong (p<0.0001 for the three variables). The frequency of dry mouth was significantly lower in people living in countries exposed to river flooding (p<0.0001) and coastal flooding (p<0.0001). People living in countries included in the worse climate scenarios for extreme temperature (p<0.0001) and river flooding (p<0.0001) showed a higher mean ESSDAI score in comparison with people living in no-risk countries. In contrast, those living in countries exposed to worse climate scenarios for wind threats (p<0.0001) and coastal flooding (p<0.0001) showed a lower mean ESSDAI score in comparison with people living in no-risk countries. CONCLUSIONS: Local exposure to extreme climate-related hazards plays a role in modulating the presentation of Sjögren across countries concerning the age at which the disease is diagnosed, the frequency of dryness, and the degree of systemic activity.

Mortality risk factors in primary Sjögren syndrome: a real-world, retrospective, cohort study.

Background: What baseline predictors would be involved in mortality in people with primary Sjögren syndrome (SjS) remains uncertain. This study aimed to investigate the baseline characteristics collected at the time of diagnosis of SjS associated with mortality and to identify mortality risk factors for all-cause death and deaths related to systemic SjS activity measured by the ESSDAI score. Methods: In this international, real-world, retrospective, cohort study, we retrospectively collected data from 27 countries on mortality and causes of death from the Big Data Sjögren Registry. Inclusion criteria consisted of fulfilling 2002/2016 SjS classification criteria, and exclusion criteria included chronic HCV/HIV infections and associated systemic autoimmune diseases. A statistical approach based on a directed acyclic graph was used, with all-cause and Sjögren-related mortality as primary endpoints. The key determinants that defined the disease phenotype at diagnosis (glandular, systemic, and immunological) were analysed as independent variables. Findings: Between January 1st, 2014 and December 31, 2023, data from 11,372 patients with primary SjS (93.5% women, 78.4% classified as White, mean age at diagnosis of 51.1 years) included in the Registry were analysed. 876 (7.7%) deaths were recorded after a mean follow-up of 8.6 years (SD 7.12). Univariate analysis of prognostic factors for all-cause death identified eight Sjögren-related variables (ocular and oral tests, salivary biopsy, ESSDAI, ANA, anti-Ro, anti-La, and cryoglobulins). The multivariate CPH model adjusted for these variables and the epidemiological features showed that DAS-ESSDAI (high vs no high: HR = 1.68; 95% CI, 1.27-2.22) and cryoglobulins (positive vs negative: HR = 1.72; 95% CI, 1.22-2.42) were independent predictors of all-cause death. Of the 640 deaths with available information detailing the specific cause of death, 14% were due to systemic SjS. Univariate analysis of prognostic factors for Sjögren-cause death identified five Sjögren-related variables (oral tests, clinESSDAI, DAS-ESSDAI, ANA, and cryoglobulins). The multivariate competing risks CPH model adjusted for these variables and the epidemiological features showed that oral tests (abnormal vs normal results: HR = 1.38; 95% CI, 1.01-1.87), DAS-ESSDAI (high vs no high: HR = 1.55; 95% CI, 1.22-1.96) and cryoglobulins (positive vs negative: HR = 1.52; 95% CI, 1.16-2) were independent predictors of SjS-related death. Interpretation: The key mortality risk factors at the time of SjS diagnosis were positive cryoglobulins and a high systemic activity scored using the ESSDAI, conferring a 2-times increased risk of all-cause and SjS-related death. ESSDAI measurement and cryoglobulin testing should be considered mandatory when an individual is diagnosed with SjS. Funding: Novartis.

Increase in invasive group A streptococcal infections and emergence of novel, rapidly expanding sub-lineage of the virulent Streptococcus pyogenes M1 clone, Denmark, 2023.

A highly virulent sub-lineage of the Streptococcus pyogenes M1 clone has been rapidly expanding throughout Denmark since late 2022 and now accounts for 30% of the new invasive group A streptococcal infections. We aimed to investigate whether a shift in variant composition can account for the high incidence rates observed over winter 2022/23, or if these are better explained by the impact of COVID-19-related restrictions on population immunity and carriage of group A Streptococcus.

Multiple Correspondence Analysis and HLA-Associations of Organ Involvement in a Large Cohort of African-American and European-American Patients with Sarcoidosis.

Sarcoidosis is a systemic granulomatous disease with predominant pulmonary involvement and vast heterogeneity of clinical manifestations and disease outcomes. African American (AA) patients suffer greater morbidity and mortality. Using Multiple Correspondence Analysis, we identified seven clusters of organ involvement in European American (EA; n = 385) patients which were similar to those previously described in a Pan-European (GenPhenReSa) and a Spanish cohort (SARCOGEAS). In contrast, AA (n = 987) had six, less well-defined and overlapping clusters with little similarity to the cluster identified in the EA cohort evaluated at the same U.S. institutions. Association of cluster membership with two-digit HLA-DRB1 alleles demonstrated ancestry-specific patterns of association and replicated known HLA effects.These results further support the notion that genetically influenced immune risk profiles, which differ based on ancestry, play a role in phenotypic heterogeneity. Dissecting such risk profiles will move us closer to personalized medicine for this complex disease.

Autoantibodies identify primary Sjögren's syndrome in patients lacking serum IgG specific for Ro/SS-A and La/SS-B.

OBJECTIVE: Identify autoantibodies in anti-Ro/SS-A negative primary Sjögren's syndrome (SS). METHODS: This is a proof-of-concept, case-control study of SS, healthy (HC) and other disease (OD) controls. A discovery dataset of plasma samples (n=30 SS, n=15 HC) was tested on human proteome arrays containing 19 500 proteins. A validation dataset of plasma and stimulated parotid saliva from additional SS cases (n=46 anti-Ro+, n=50 anti-Ro-), HC (n=42) and OD (n=54) was tested on custom arrays containing 74 proteins. For each protein, the mean+3 SD of the HC value defined the positivity threshold. Differences from HC were determined by Fisher's exact test and random forest machine learning using 2/3 of the validation dataset for training and 1/3 for testing. Applicability of the results was explored in an independent rheumatology practice cohort (n=38 Ro+, n=36 Ro-, n=10 HC). Relationships among antigens were explored using Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) interactome analysis. RESULTS: Ro+ SS parotid saliva contained autoantibodies binding to Ro60, Ro52, La/SS-B and muscarinic receptor 5. SS plasma contained 12 novel autoantibody specificities, 11 of which were detected in both the discovery and validation datasets. Binding to ≥1 of the novel antigens identified 54% of Ro- SS and 37% of Ro+ SS cases, with 100% specificity in both groups. Machine learning identified 30 novel specificities showing receiver operating characteristic area under the curve of 0.79 (95% CI 0.64 to 0.93) for identifying Ro- SS. Sera from Ro- cases of an independent cohort bound 17 of the non-canonical antigens. Antigenic targets in both Ro+ and Ro- SS were part of leukaemia cell, ubiquitin conjugation and antiviral defence pathways. CONCLUSION: We identified antigenic targets of the autoantibody response in SS that may be useful for identifying up to half of Ro seronegative SS cases.

Anti-vimentin antibodies are associated with higher severity of Sjögren's disease.

Vimentin is a ubiquitously present Type III intermediate filament protein, often targeted by autoimmune responses in multiple rheumatic disorders. Although previous studies have reported anti-vimentin antibodies in Sjögren's disease (SjD) patients, the clinical significance of such antibodies is unknown. To address this issue, the presence of anti-vimentin antibodies was determined in serum samples from a well-characterized cohort of primary SjD patients, non-SjD Sicca, and healthy controls. The occurrence of anti-vimentin antibodies and their association with different clinical features of the disease were evaluated. Anti-vimentin antibodies were detected in 24% of primary SjD patients, compared to 4% in non-SjD sicca patients and 3% in healthy controls. In primary SjD patients, higher levels of anti-vimentin antibodies were significantly associated with reduced saliva and tear flow and severe ocular surface damage indicators. The anti-vimentin antibody levels did not show significant associations with the presence or absence of other autoantibodies like ANA, RF, and anti-Ro/La. Our data suggest that the anti-vimentin antibody specificity arises in a subset of primary SjD patients and is associated with oral and ocular features of the disease. Anti-vimentin can potentially serve as a novel biomarker for evaluating the severity of salivary and lacrimal gland dysfunction in primary SjD.

Dysregulated long non-coding RNA in Sjögren's disease impacts both interferon and adaptive immune responses.

OBJECTIVE: Sjögren's disease (SjD) is an autoimmune disease characterised by inflammatory destruction of exocrine glands. Patients with autoantibodies to Ro/SSA (SjDRo+) exhibit more severe disease. Long non-coding RNAs (lncRNAs) are a functionally diverse class of non-protein-coding RNAs whose role in autoimmune disease pathology has not been well characterised. METHODS: Whole blood RNA-sequencing (RNA-seq) was performed on SjD cases (n=23 Ro/SSA negative (SjDRo-); n=27 Ro/SSA positive (SjDRo+) and healthy controls (HCs; n=27). Bioinformatics and pathway analyses of differentially expressed (DE) transcripts (log2 fold change ≥2 or ≤0.5; padj<0.05) were used to predict lncRNA function. LINC01871 was characterised by RNA-seq analyses of HSB-2 cells with CRISPR-targeted LINC01871 deletion (LINC01871-/ -) and in vitro stimulation assays. RESULTS: Whole blood RNA-seq revealed autoantibody-specific transcription profiles and disproportionate downregulation of DE transcripts in SjD cases relative to HCs. Sixteen DE lncRNAs exhibited correlated expression with the interferon (IFN)-regulated gene, RSAD2, in SjDRo+ (r≥0.65 or ≤-0.6); four antisense lncRNAs exhibited IFN-regulated expression in immune cell lines. LINC01871 was upregulated in all SjD cases. RNA-seq and pathway analyses of LINC01871-/ - cells implicated roles in cytotoxic function, differentiation and IFNγ induction. LINC01871 was induced by IFNγ in a myeloid cell line and regulated by calcineurin/NFAT pathway and T cell receptor (TCR) signalling in primary human T cells. CONCLUSION: LINC01871 influences expression of many immune cell genes and growth factors, is IFNγ inducible, and regulated by calcineurin signalling and TCR ligand engagement. Altered LINC01871 expression may influence the dysregulated T cell inflammatory pathways implicated in SjD.

Genetics of neurosarcoidosis.

Sarcoidosis is a systemic, inflammatory, granulomatous disease characterized by great variability in organ involvement, clinical course, and severity. While pulmonary manifestations are almost universal, the central and peripheral nervous systems can also be affected. Neurosarcoidosis occurs in ∼5-15% of cases and is among the manifestations with the highest morbidity and mortality. It is known that sarcoidosis has genetic underpinnings and while multiple studies aimed at identifying associations to sarcoidosis susceptibility and prognosis, very few studies have focused on neurosarcoidosis. This review summarizes the genetic studies to date, compares and contrasts those findings with other genetic effects in sarcoidosis, and offers ideas for moving the field forward.

Genome-wide association study identifies Sjögren's risk loci with functional implications in immune and glandular cells.

Sjögren's disease is a complex autoimmune disease with twelve established susceptibility loci. This genome-wide association study (GWAS) identifies ten novel genome-wide significant (GWS) regions in Sjögren's cases of European ancestry: CD247, NAB1, PTTG1-MIR146A, PRDM1-ATG5, TNFAIP3, XKR6, MAPT-CRHR1, RPTOR-CHMP6-BAIAP6, TYK2, SYNGR1. Polygenic risk scores yield predictability (AUROC = 0.71) and relative risk of 12.08. Interrogation of bioinformatics databases refine the associations, define local regulatory networks of GWS SNPs from the 95% credible set, and expand the implicated gene list to >40. Many GWS SNPs are eQTLs for genes within topologically associated domains in immune cells and/or eQTLs in the main target tissue, salivary glands.

High Resolution Analysis of DMPK Hypermethylation and Repeat Interruptions in Myotonic Dystrophy Type 1.

Myotonic dystrophy type 1 (DM1) is a multisystemic neuromuscular disorder caused by the expansion of a CTG repeat in the 3'-UTR of DMPK, which is transcribed to a toxic gain-of-function RNA that affects splicing of a range of genes. The expanded repeat is unstable in both germline and somatic cells. The variable age at disease onset and severity of symptoms have been linked to the inherited CTG repeat length, non-CTG interruptions, and methylation levels flanking the repeat. In general, the genetic biomarkers are investigated separately with specific methods, making it tedious to obtain an overall characterisation of the repeat for a given individual. In the present study, we employed Oxford nanopore sequencing in a pilot study to simultaneously determine the repeat lengths, investigate the presence and nature of repeat interruptions, and quantify methylation levels in the regions flanking the CTG-repeats in four patients with DM1. We determined the repeat lengths, and in three patients, we observed interruptions which were not detected using repeat-primed PCR. Interruptions may thus be more common than previously anticipated and should be investigated in larger cohorts. Allele-specific analyses enabled characterisation of aberrant methylation levels specific to the expanded allele, which greatly increased the sensitivity and resolved cases where the methylation levels were ambiguous.

Novel HLA associations with outcomes of Mycobacterium tuberculosis exposure and sarcoidosis in individuals of African ancestry using nearest-neighbor feature selection.

Tuberculosis and sarcoidosis are inflammatory diseases characterized by granulomas that may occur in any organ but are often found in the lung. The panoply of classical human leukocyte antigen (HLA) alleles associated with occurrence and/or severity of both diseases varies considerably across studies. This heterogeneity of results, due to variation in factors like ancestry and disease subphenotype, as well as the use of simple modeling strategies to elucidate likely complex relationships, has made conclusions about underlying commonalities difficult. Here we perform HLA association analyses in individuals of African ancestry, using a greater resolution to include subphenotypes of disease and employing more comprehensive analytical techniques. Using a novel application of nearest-neighbor feature selection to score allelic importance, we investigated HLA allele association with Mycobacterium tuberculosis exposure outcomes in the first analysis of both latent Mycobacterium tuberculosis infection and active disease compared with those who, despite long-term exposure to active index cases, have neither positive diagnostic tests nor display clinical symptoms. We also compared persistent to resolved sarcoidosis. This led to the identification of novel HLA associations and evidence of main effects and interaction effects. We found strikingly similar main effects and interaction effects at HLA-DRB1, -DQB1, and -DPB1 in those resistant to tuberculosis (either latent or active) and persistent sarcoidosis.