Evidence supports a causal association between allele-specific vitamin D receptor binding and multiple sclerosis among Europeans.

Although evidence exists for a causal association between 25-hydroxyvitamin D (25(OH)D) serum levels, and multiple sclerosis (MS), the role of variation in vitamin D receptor (VDR) binding in MS is unknown. Here, we leveraged previously identified variants associated with allele imbalance in VDR binding (VDR-binding variant; VDR-BV) in ChIP-exo data from calcitriol-stimulated lymphoblastoid cell lines and 25(OH)D serum levels from genome-wide association studies to construct genetic instrumental variables (GIVs). GIVs are composed of one or more genetic variants that serve as proxies for exposures of interest. Here, GIVs for both VDR-BVs and 25(OH)D were used in a two-sample Mendelian Randomization study to investigate the relationship between VDR binding at a locus, 25(OH)D serum levels, and MS risk. Data for 13,598 MS cases and 38,887 controls of European ancestry from Kaiser Permanente Northern California, Swedish MS studies, and the UK Biobank were included. We estimated the association between each VDR-BV GIV and MS. Significant interaction between a VDR-BV GIV and a GIV for serum 25OH(D) was evidence for a causal association between VDR-BVs and MS unbiased by pleiotropy. We observed evidence for associations between two VDR-BVs (rs2881514, rs2531804) and MS after correction for multiple tests. There was evidence of interaction between rs2881514 and a 25(OH)D GIV, providing evidence of a causal association between rs2881514 and MS. This study is the first to demonstrate evidence that variation in VDR binding at a locus contributes to MS risk. Our results are relevant to other autoimmune diseases in which vitamin D plays a role.

Rare and low-frequency coding genetic variants contribute to pediatric-onset multiple sclerosis.

BACKGROUND: Rare genetic variants are emerging as important contributors to the heritability of multiple sclerosis (MS). Whether rare variants also contribute to pediatric-onset multiple sclerosis (POMS) is unknown. OBJECTIVE: To test whether genes harboring rare variants associated with adult-onset MS risk (PRF1, PRKRA, NLRP8, and HDAC7) and 52 major histocompatibility complex (MHC) genes are associated with POMS. METHODS: We analyzed DNA samples from 330 POMS cases and 306 controls from the US Network of Pediatric MS Centers and Kaiser Permanente Northern California for which Illumina ExomeChip genotypes were available. Using the gene-based method "SKAT-O," we tested the association between candidate genes and POMS risk. RESULTS: After correction for multiple comparisons, one adult-onset MS gene (PRF1, p = 2.70 × 10-3) and two MHC genes (BRD2, p = 5.89 × 10-5 and AGER, p = 7.96 × 10-5) were significantly associated with POMS. Results suggest these are independent of HLA-DRB1*1501. CONCLUSION: Findings support a role for rare coding variants in POMS susceptibility. In particular, rare minor alleles within PRF1 were more common among individuals with POMS compared to controls while the opposite was true for rare variants within significant MHC genes, BRD2 and AGER. These genes would not have been identified by common variant studies, emphasizing the merits of investigating rare genetic variation in complex diseases.

Mother-child histocompatibility and risk of rheumatoid arthritis and systemic lupus erythematosus among mothers.

The study objective was to test the hypothesis that having histocompatible children increases the risk of rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE), possibly by contributing to the persistence of fetal cells acquired during pregnancy. We conducted a case control study using data from the UC San Francisco Mother Child Immunogenetic Study and studies at the Inova Translational Medicine Institute. We imputed human leukocyte antigen (HLA) alleles and minor histocompatibility antigens (mHags). We created a variable of exposure to histocompatible children. We estimated an average sequence similarity matching (SSM) score for each mother based on discordant mother-child alleles as a measure of histocompatibility. We used logistic regression models to estimate odds ratios (ORs) and 95% confidence intervals. A total of 138 RA, 117 SLE, and 913 control mothers were analyzed. Increased risk of RA was associated with having any child compatible at HLA-B (OR 1.9; 1.2-3.1), DPB1 (OR 1.8; 1.2-2.6) or DQB1 (OR 1.8; 1.2-2.7). Compatibility at mHag ZAPHIR was associated with reduced risk of SLE among mothers carrying the HLA-restriction allele B*07:02 (n = 262; OR 0.4; 0.2-0.8). Our findings support the hypothesis that mother-child histocompatibility is associated with risk of RA and SLE.

Vitamin D genes influence MS relapses in children.

OBJECTIVE: The aim of this study was to determine whether a vitamin D genetic risk score (vitDGRS) is associated with 25-hydroxyvitamin D (25(OH)D) level and multiple sclerosis (MS) relapses in children. METHODS: DNA samples were typed for single nucleotide polymorphisms (SNPs) from four genes previously identified to be associated with 25(OH)D levels. SNPs with strong associations with 25(OH)D after multiple comparison correction were used to create a genetic risk score (vitDGRS). Cox regression models tested associations of vitDGRS with relapse hazard. RESULTS: Two independent SNPs within or near GC and NADSYN1/DHCR7 genes were strongly associated with 25(OH)D levels in the discovery cohort (n = 182) after Bonferroni correction. The vitDGRS of these SNPs explained 4.5% of the variance of 25(OH)D level after adjustment for genetic ancestry. Having the highest versus lowest vitDGRS was associated with 11 ng/mL lower 25(OH)D level (95% confidence interval (CI) = -17.5, -4.5, p = 0.001) in the discovery cohort. Adjusting for ancestry, sex, disease-modifying therapy (DMT), and HLA-DRB1*15 carrier status, the highest versus lowest vitDGRS was associated with 2.6-fold (95% CI = 1.37, 5.03, p = 0.004) and 2.0-fold (95% CI = 0.75, 5.20, p = 0.16) higher relapse hazard in the discovery and replication cohorts, respectively. CONCLUSION: The vitDGRS identifies children at greater risk of relapse. These findings support a causal role for vitamin D in MS course.

Genetic risk factors for pediatric-onset multiple sclerosis.

BACKGROUND: Strong evidence supports the role of both genetic and environmental factors in pediatric-onset multiple sclerosis (POMS) etiology. OBJECTIVE: We comprehensively investigated the association between established major histocompatibility complex (MHC) and non-MHC adult multiple sclerosis (MS)-associated variants and susceptibility to POMS. METHODS: Cases with onset <18 years (n = 569) and controls (n = 16,251) were included from the United States and Sweden. Adjusted logistic regression and meta-analyses were performed for individual risk variants and a weighted genetic risk score (wGRS) for non-MHC variants. Results were compared to adult MS cases (n = 7588). RESULTS: HLA-DRB1*15:01 was strongly associated with POMS (odds ratio (OR)meta = 2.95, p < 2.0 × 10-16). Furthermore, 28 of 104 non-MHC variants studied (23%) were associated (p < 0.05); POMS cases carried, on average, a higher burden of these 28 variants compared to adults (ORavg = 1.24 vs 1.13, respectively), though the difference was not significant. The wGRS was strongly associated with POMS (ORmeta = 2.77, 95% confidence interval: 2.33, 3.32, p < 2.0 × 10-16) and higher, on average, when compared to adult cases. Additional class III risk variants in the MHC region associated with POMS were revealed after accounting for HLA-DRB1*15:01 and HLA-A*02. CONCLUSION: Pediatric and adult MS share many genetic variants suggesting similar biological processes are present. MHC variants beyond HLA-DRB1*15:01 and HLA-A*02 are also associated with POMS.

Causal Effect of Genetic Variants Associated With Body Mass Index on Multiple Sclerosis Susceptibility.

Multiple sclerosis (MS) is an autoimmune disease with both genetic and environmental risk factors. Recent studies indicate that childhood and adolescent obesity double the risk of MS, but this association may reflect unmeasured confounders rather than causal effects of obesity. We used separate-sample Mendelian randomization to estimate the causal effect of body mass index (BMI) on susceptibility to MS. Using data from non-Hispanic white members of the Kaiser Permanente Medical Care Plan of Northern California (KPNC) (2006-2014; 1,104 cases of MS and 10,536 controls) and a replication data set from Sweden (the Epidemiological Investigation of MS (EIMS) and the Genes and Environment in MS (GEMS) studies, 2005-2013; 5,133 MS cases and 4,718 controls), we constructed a weighted genetic risk score using 97 variants previously established to predict BMI. Results were adjusted for birth year, sex, education, smoking status, ancestry, and genetic predictors of MS. Estimates in KPNC and Swedish data sets suggested that higher genetically induced BMI predicted greater susceptibility to MS (odds ratio = 1.13, 95% confidence interval: 1.04, 1.22 for the KPNC sample; odds ratio = 1.09, 95% confidence interval: 1.03, 1.15 for the Swedish sample). Although the mechanism remains unclear, to our knowledge, these findings support a causal effect of increased BMI on susceptibility to MS for the first time, and they suggest a role for inflammatory pathways that characterize both obesity and the MS disease process.

Increased risk of rheumatoid arthritis among mothers with children who carry DRB1 risk-associated alleles.

OBJECTIVE: To investigate whether a child's genotype affects a mother's risk of rheumatoid arthritis (RA) beyond the risk associated with her genotype and to test whether exposure to fetal alleles inherited from the father increases risk of RA among mothers without risk alleles. METHODS: A case-control study was conducted among 1165 mothers (170 cases/995 controls) and their respective 1482 children. We tested the association between having any child with alleles encoding amino acids (AAs) associated with RA including the 'shared epitope' (SE) and DERAA AA sequences at positions 70-74; AA valine, lysine and alanine at positions 11, 71 and 74 of HLA-DRB1; aspartic acid at position 9 of HLA-B and phenylalanine at position 9 of DPB1. We used logistic regression models to estimate OR and 95% CI for each group of alleles, adjusting for maternal genotype and number of live births. RESULTS: We found increased risk of RA among mothers who had any child with SE (OR 3.0; 95% CI 2.0 to 4.6); DERAA (OR 1.7; 95% CI 1.1 to 2.6); or valine (OR 2.3; 95% CI 1.6 to 3.5), lysine (OR 2.3; 95% CI 1.5 to 3.4) and alanine (OR 2.8; 95% CI 1.2 to 6.4) at DRB1 positions 11, 71 and 74, respectively. Among non-carrier mothers, increased risk of RA was associated with having children who carried DERAA (OR 1.7; 95% CI 1.0 to 2.7) and alleles encoding lysine at DRB1 position 71 (OR 2.3; 95% CI 1.5 to 4.8). CONCLUSION: Findings support the hypothesis that a child's genotype can contribute independently to risk of RA among mothers.

A Child's HLA-DRB1 genotype increases maternal risk of systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) disproportionately affects women of reproductive age. During pregnancy, women are exposed to various sources of fetal material possibly constituting a significant immunologic exposure relevant to the development of SLE. The objective of this study was to investigate whether having any children who carry DRB1 alleles associated with SLE increase the risk of maternal SLE. This case-control study is based on the University of California, San Francisco Mother-Child Immunogenetic Study and from studies at the Inova Translational Medicine Institute. Analyses were conducted using data for 1304 mothers (219 cases/1085 controls) and their respective 1664 children. We selected alleles based on their known association with risk of SLE (DRB1*03:01, *15:01, or *08:01) or Epstein-Barr virus (EBV) glycoproteins (*04:01) due to the established EBV association with SLE risk. We used logistic regression models to estimate odds ratios (OR) and 95% confidence intervals (CI) for each allele of interest, taking into account maternal genotype and number of live births. We found an increase in risk of maternal SLE associated with exposure to children who inherited DRB1*04:01 from their father (OR 1.9; 95% CI, 1.1-3.2), among *04:01 allele-negative mothers. Increased risk was only present among mothers who were positive for one or more SLE risk-associated alleles (*03:01, *15:01 and/or *08:01). We did not find increased risk of maternal SLE associated with any other tested allele. These findings support the hypothesis that a child's alleles inherited from the father influence a mother's subsequent risk of SLE.

Hypomethylation within gene promoter regions and type 1 diabetes in discordant monozygotic twins.

Genetic susceptibility to type 1 diabetes (T1D) is well supported by epidemiologic evidence; however, disease risk cannot be entirely explained by established genetic variants identified so far. This study addresses the question of whether epigenetic modification of the inherited DNA sequence may contribute to T1D susceptibility. Using the Infinium HumanMethylation450 BeadChip array (450k), a total of seven long-term disease-discordant monozygotic (MZ) twin pairs and five pairs of HLA-identical, disease-discordant non-twin siblings (NTS) were examined for associations between DNA methylation (DNAm) and T1D. Strong evidence for global hypomethylation of CpG sites within promoter regions in MZ twins with TID compared to twins without T1D was observed. DNA methylation data were then grouped into three categories of CpG sites for further analysis, including those within: 1) the major histocompatibility complex (MHC) region, 2) non-MHC genes with reported T1D association through genome wide association studies (GWAS), and 3) the epigenome, or remainder of sites that did not include MHC and T1D associated genes. Initial results showed modest methylation differences between discordant MZ twins for the MHC region and T1D-associated CpG sites, BACH2, INS-IGF2, and CLEC16A (DNAm difference range: 2.2%-5.0%). In the epigenome CpG set, the greatest methylation differences were observed in MAGI2, FANCC, and PCDHB16, (DNAm difference range: 6.9%-16.1%). These findings were not observed in the HLA-identical NTS pairs. Targeted pyrosequencing of five candidate CpG loci identified using the 450k array in the original discordant MZ twins produced similar results using control DNA samples, indicating strong agreement between the two DNA methylation profiling platforms. However, findings for the top five candidate CpG loci were not replicated in six additional T1D-discordant MZ twin pairs. Our results indicate global DNA hypomethylation within gene promoter regions may contribute to T1D; however, findings do not support the involvement of large DNAm differences at single CpG sites alone in T1D.

Genetic predictors of relapse rate in pediatric MS.

BACKGROUND: Genetic ancestry, sex, and individual alleles have been associated with multiple sclerosis (MS) susceptibility. OBJECTIVE: To determine whether established risk factors for disease onset are associated with relapse rate in pediatric MS. METHODS: Whole-genome genotyping was performed for 181 MS or high-risk clinically isolated syndrome patients from two pediatric MS centers. Relapses and disease-modifying therapies were recorded as part of continued follow-up. Participants were characterized for 25-hydroxyvitamin D serum status. Ancestral estimates (STRUCTURE v2.3.1), human leukocyte antigen (HLA)-DRB1*15 carrier status (direct sequencing), sex, and a genetic risk score (GRS) of 110 non-HLA susceptibility single-nucleotide polymorphisms (SNPs) were evaluated for association with relapse rate with Cox and negative binomial regression models. RESULTS: Over 622 patient-years, 408 relapses were captured. Girls had greater relapse rate than boys (incident rate ratio (IRR) = 1.40, 95% confidence interval (CI) = 1.04-1.87, p = 0.026). Participants were genetically diverse; ~40% (N = 75) had <50% European ancestry. HLA-DRB1*15 status modified the association of vitamin D status (pixn = 0.022) with relapse rate (per 10 ng/mL, in DRB1*15+ hazard ratio (HR) = 0.72, 95% CI = 0.58-0.88, p = 0.002; in DRB1*15- HR = 0.96, 95% CI = 0.83-1.12, p = 0.64). Neither European ancestry nor GRS was associated with relapse rate. CONCLUSION: We demonstrate that HLA-DRB1*15 modifies the association of vitamin D status with relapse rate. Our findings emphasize the need to pursue disease-modifying effects of MS genes in the context of environmental factors.

Smoking and risk of multiple sclerosis: evidence of modification by NAT1 variants.

BACKGROUND: Tobacco smoke is an established risk factor for multiple sclerosis (MS). We hypothesized that variation in genes involved in metabolism of tobacco smoke constituents may modify MS risk in smokers. METHODS: A three-stage gene-environment investigation was conducted for NAT1, NAT2, and GSTP1 variants. The discovery analysis was conducted among 1588 white MS cases and controls from the Kaiser Permanente Northern California Region HealthPlan (Kaiser). The replication analysis was carried out in 988 white MS cases and controls from Sweden. RESULTS: Tobacco smoke exposure at the age of 20 years was associated with greater MS risk in both data sets (in Kaiser, odds ratio [OR] = 1.51 [95% confidence interval (CI) = 1.17-1.93]; in Sweden, OR = 1.35 [1.04-1.74]). A total of 42 NAT1 variants showed evidence for interaction with tobacco smoke exposure (P(corrected) < 0.05). Genotypes for 41 NAT1 single nucleotide polymorphisms (SNPs) were studied in the replication data set. A variant (rs7388368C>A) within a dense transcription factor-binding region showed evidence for interaction (Kaiser, OR for interaction = 1.75 [95% CI = 1.19-2.56]; Sweden, OR = 1.62 [1.05-2.49]). Tobacco smoke exposure was associated with MS risk among rs7388368A carriers only; homozygote individuals had the highest risk (A/A, OR = 5.17 [95% CI = 2.17-12.33]). CONCLUSIONS: We conducted a three-stage analysis using two population-based case-control datasets that consisted of a discovery population, a replication population, and a pooled analysis. NAT1 emerged as a genetic effect modifier of tobacco smoke exposure in MS susceptibility.

Identification of Sjögren's syndrome patient subgroups by clustering of labial salivary gland DNA methylation profiles.

Heterogeneity in Sjögren's syndrome (SS), increasingly called Sjögren's disease, suggests the presence of disease subtypes, which poses a major challenge for the diagnosis, management, and treatment of this autoimmune disorder. Previous work distinguished patient subgroups based on clinical symptoms, but it is not clear to what extent symptoms reflect underlying pathobiology. The purpose of this study was to discover clinical meaningful subtypes of SS based on genome-wide DNA methylation data. We performed a cluster analysis of genome-wide DNA methylation data from labial salivary gland (LSG) tissue collected from 64 SS cases and 67 non-cases. Specifically, hierarchical clustering was performed on low dimensional embeddings of DNA methylation data extracted from a variational autoencoder to uncover unknown heterogeneity. Clustering revealed clinically severe and mild subgroups of SS. Differential methylation analysis revealed that hypomethylation at the MHC and hypermethylation at other genome regions characterize the epigenetic differences between these SS subgroups. Epigenetic profiling of LSGs in SS yields new insights into mechanisms underlying disease heterogeneity. The methylation patterns at differentially methylated CpGs are different in SS subgroups and support the role of epigenetic contributions to the heterogeneity in SS. Biomarker data derived from epigenetic profiling could be explored in future iterations of the classification criteria for defining SS subgroups.

Identification of Cell-Specific Differential DNA Methylation Associated With Methotrexate Treatment Response in Rheumatoid Arthritis.

OBJECTIVE: We undertook this study to estimate changes in cell-specific DNA methylation (DNAm) associated with methotrexate (MTX) response using whole blood samples collected from rheumatoid arthritis (RA) patients before and after initiation of MTX treatment. METHODS: Patients included in this study were from the Rheumatoid Arthritis Medication Study (n = 66) and the University of California San Francisco Rheumatoid Arthritis study (n = 11). All patients met the American College of Rheumatology RA classification criteria. Blood samples were collected at baseline and following treatment. Disease Activity Scores in 28 joints using the C-reactive protein level were collected at baseline and after 3-6 months of treatment with MTX. Methylation profiles were generated using the Illumina Infinium HumanMethylation450 and MethylationEPIC v1.0 BeadChip arrays using DNA from whole blood. MTX response was defined using the EULAR response criteria (responders showed good/moderate response; nonresponders showed no response). Differentially methylated positions were identified using the Limma software package and Tensor Composition Analysis, which is a method for identifying cell-specific differential DNAm at the CpG level from tissue-level ("bulk") data. Differentially methylated regions were identified using Comb-p software. RESULTS: We found evidence of differential global methylation between treatment response groups. Further, we found patterns of cell-specific differential global methylation associated with MTX response. After correction for multiple testing, 1 differentially methylated position was associated with differential DNAm between responders and nonresponders at baseline in CD4+ T cells, CD8+ T cells, and natural killer cells. Thirty-nine cell-specific differentially methylated regions associated with MTX treatment response were identified. There were no significant findings in analyses of whole blood samples. CONCLUSION: We identified cell-specific changes in DNAm that were associated with MTX treatment response in RA patients. Future studies of DNAm and MTX treatment response should include measurements of DNAm from sorted cells.

Cross-Trait Mendelian Randomization Study to Investigate Whether Migraine Is a Risk Factor for Multiple Sclerosis.

BACKGROUND AND OBJECTIVES: Migraine is common among people with multiple sclerosis (MS), but the reasons for this are unknown. We tested 3 hypothesized mechanisms for this observed comorbidity, including migraine is a risk factor of MS, genetic variants are shared between the conditions, and migraine is because of MS. METHODS: Data were from 2 sources: publicly available summary statistics from genome-wide association studies of MS (N = 115,748) and migraine (N = 375,752 and N = 361,141) and a case-control study of MS recruited from the Kaiser Permanente Northern California Health Plan (N = 1,991). For the latter participants, migraine status was ascertained using a validated electronic health record migraine probability algorithm or self-report. Using the public summary statistics, we used 2-sample Mendelian randomization to test whether a migraine genetic instrumental variable was associated with MS. We used linkage disequilibrium score regression and LOGODetect to ascertain whether MS and migraine shared genetic variants across the genome and regionally. Using the Northern California MS cohort, we used logistic regression to identify whether people with both MS and migraine had different odds of clinical characteristics (e.g., age at MS onset, Perceived Deficits Questionnaire, and depression) or MS-specific risk factors (e.g., body mass index, smoking status, and infectious mononucleosis status) compared with people with MS without migraine. RESULTS: We did not find evidence supporting migraine as a causal risk factor of MS (p = 0.29). We did, however, identify 4 major histocompatibility complex (MHC) loci shared between MS and migraine. Among the Northern California MS cohort, 774 (39%) experienced migraine. People with both MS and migraine from this cohort were more likely to ever smoke (odds ratio [OR] = 1.30, 95% CI: 1.08-1.57), have worse self-reported cognitive deficits (OR = 1.04, 95% CI: 1.02-1.06), and ever experience depression (OR = 1.48, 95% CI: 1.22-1.80). DISCUSSION: Our findings do not support migraine as a causal risk factor of MS. Several genetic variants, particularly in the MHC, may account for some of the overlap. It seems likely that migraine within the context of MS is because of MS. Identifying what increases the risk of migraine within MS might lead to an improved treatment and quality of life.

High-density SNP screening of the major histocompatibility complex in systemic lupus erythematosus demonstrates strong evidence for independent susceptibility regions.

A substantial genetic contribution to systemic lupus erythematosus (SLE) risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6p21. Previous studies in SLE have lacked statistical power and genetic resolution to fully define MHC influences. We characterized 1,610 Caucasian SLE cases and 1,470 parents for 1,974 MHC SNPs, the highly polymorphic HLA-DRB1 locus, and a panel of ancestry informative markers. Single-marker analyses revealed strong signals for SNPs within several MHC regions, as well as with HLA-DRB1 (global p = 9.99 x 10(-16)). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53 x 10(-12)), *1401 (OR = 0.50, p = 0.0002), and *1501 (OR = 1.39, p = 0.0032). The MHC region SNP demonstrating the strongest evidence of association with SLE was rs3117103, with OR = 2.44 and p = 2.80 x 10(-13). Conditional haplotype and stepwise logistic regression analyses identified strong evidence for association between SLE and the extended class I, class I, class III, class II, and the extended class II MHC regions. Sequential removal of SLE-associated DRB1 haplotypes revealed independent effects due to variation within OR2H2 (extended class I, rs362521, p = 0.006), CREBL1 (class III, rs8283, p = 0.01), and DQB2 (class II, rs7769979, p = 0.003, and rs10947345, p = 0.0004). Further, conditional haplotype analyses demonstrated that variation within MICB (class I, rs3828903, p = 0.006) also contributes to SLE risk independent of HLA-DRB1*0301. Our results for the first time delineate with high resolution several MHC regions with independent contributions to SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.

Development and Implementation of Dried Blood Spot-Based COVID-19 Serological Assays for Epidemiologic Studies.

Serological surveillance studies of infectious diseases provide population-level estimates of infection and antibody prevalence, generating crucial insight into population-level immunity, risk factors leading to infection, and effectiveness of public health measures. These studies traditionally rely on detection of pathogen-specific antibodies in samples derived from venipuncture, an expensive and logistically challenging aspect of serological surveillance. During the COVID-19 pandemic, guidelines implemented to prevent the spread of SARS-CoV-2 infection made collection of venous blood logistically difficult at a time when SARS-CoV-2 serosurveillance was urgently needed. Dried blood spots (DBS) have generated interest as an alternative to venous blood for SARS-CoV-2 serological applications due to their stability, low cost, and ease of collection; DBS samples can be self-generated via fingerprick by community members and mailed at ambient temperatures. Here, we detail the development of four DBS-based SARS-CoV-2 serological methods and demonstrate their implementation in a large serological survey of community members from 12 cities in the East Bay region of the San Francisco metropolitan area using at-home DBS collection. We find that DBS perform similarly to plasma/serum in enzyme-linked immunosorbent assays and commercial SARS-CoV-2 serological assays. In addition, we show that DBS samples can reliably detect antibody responses months postinfection and track antibody kinetics after vaccination. Implementation of DBS enabled collection of valuable serological data from our study population to investigate changes in seroprevalence over an 8-month period. Our work makes a strong argument for the implementation of DBS in serological studies, not just for SARS-CoV-2, but any situation where phlebotomy is inaccessible. IMPORTANCE Estimation of community-level antibody responses to SARS-CoV-2 from infection or vaccination is critical to inform public health responses. Traditional studies of antibodies rely on collection of blood via venipuncture, an invasive procedure not amenable to pandemic-related social-distancing measures. Dried blood spots (DBS) are an alternative to venipuncture, since they can be self-collected by study participants at home and do not require refrigeration for shipment or storage. However, DBS-based assays to measure antibody levels to SARS-CoV-2 have not been widely utilized. Here, we show that DBS are comparable to blood as a sampling method for antibody responses to SARS-CoV-2 infection and vaccination over time measured using four distinct serological assays. The DBS format enabled antibody surveillance in a longitudinal cohort where study participants self-collected samples, ensuring the participants' safety during an ongoing pandemic. Our work demonstrates that DBS are an excellent sampling method for measuring antibody responses whenever venipuncture is impractical.

Health inequities in SARS-CoV-2 infection, seroprevalence, and COVID-19 vaccination: Results from the East Bay COVID-19 study.

Comprehensive data on transmission mitigation behaviors and both SARS-CoV-2 infection and serostatus are needed from large, community-based cohorts to identify COVID-19 risk factors and the impact of public health measures. We conducted a longitudinal, population-based study in the East Bay Area of Northern California. From July 2020-March 2021, approximately 5,500 adults were recruited and followed over three data collection rounds to investigate the association between geographic and demographic characteristics and transmission mitigation behavior with SARS-CoV-2 prevalence. We estimated the populated-adjusted prevalence of antibodies from SARS-CoV-2 infection and COVID-19 vaccination, and self-reported COVID-19 test positivity. Population-adjusted SARS-CoV-2 seroprevalence was low, increasing from 1.03% (95% CI: 0.50-1.96) in Round 1 (July-September 2020), to 1.37% (95% CI: 0.75-2.39) in Round 2 (October-December 2020), to 2.18% (95% CI: 1.48-3.17) in Round 3 (February-March 2021). Population-adjusted seroprevalence of COVID-19 vaccination was 21.64% (95% CI: 19.20-24.34) in Round 3, with White individuals having 4.35% (95% CI: 0.35-8.32) higher COVID-19 vaccine seroprevalence than individuals identifying as African American or Black, American Indian or Alaskan Native, Asian, Hispanic, two or more races, or other. No evidence for an association between transmission mitigation behavior and seroprevalence was observed. Despite >99% of participants reporting wearing masks individuals identifying as African American or Black, American Indian or Alaskan Native, Asian, Hispanic, two or more races, or other, as well as those in lower-income households, and lower-educated individuals had the highest SARS-CoV-2 seroprevalence and lowest vaccination seroprevalence. Results demonstrate that more effective policies are needed to address these disparities and inequities.

Correction: Hypomethylation mediates genetic association with the major histocompatibility complex genes in Sjögren's syndrome.

[This corrects the article DOI: 10.1371/journal.pone.0248429.].

Hypomethylation mediates genetic association with the major histocompatibility complex genes in Sjögren's syndrome.

Differential methylation of immune genes has been a consistent theme observed in Sjögren's syndrome (SS) in CD4+ T cells, CD19+ B cells, whole blood, and labial salivary glands (LSGs). Multiple studies have found associations supporting genetic control of DNA methylation in SS, which in the absence of reverse causation, has positive implications for the potential of epigenetic therapy. However, a formal study of the causal relationship between genetic variation, DNA methylation, and disease status is lacking. We performed a causal mediation analysis of DNA methylation as a mediator of nearby genetic association with SS using LSGs and genotype data collected from 131 female members of the Sjögren's International Collaborative Clinical Alliance registry, comprising of 64 SS cases and 67 non-cases. Bumphunter was used to first identify differentially-methylated regions (DMRs), then the causal inference test (CIT) was applied to identify DMRs mediating the association of nearby methylation quantitative trait loci (MeQTL) with SS. Bumphunter discovered 215 DMRs, with the majority located in the major histocompatibility complex (MHC) on chromosome 6p21.3. Consistent with previous findings, regions hypomethylated in SS cases were enriched for gene sets associated with immune processes. Using the CIT, we observed a total of 19 DMR-MeQTL pairs that exhibited strong evidence for a causal mediation relationship. Close to half of these DMRs reside in the MHC and their corresponding meQTLs are in the region spanning the HLA-DQA1, HLA-DQB1, and HLA-DQA2 loci. The risk of SS conferred by these corresponding MeQTLs in the MHC was further substantiated by previous genome-wide association study results, with modest evidence for independent effects. By validating the presence of causal mediation, our findings suggest both genetic and epigenetic factors contribute to disease susceptibility, and inform the development of targeted epigenetic modification as a therapeutic approach for SS.

Global expression and CpG methylation analysis of primary endothelial cells before and after TNFa stimulation reveals gene modules enriched in inflammatory and infectious diseases and associated DMRs.

Endothelial cells are a primary site of leukocyte recruitment during inflammation. An increase in tumor necrosis factor-alpha (TNFa) levels as a result of infection or some autoimmune diseases can trigger this process. Several autoimmune diseases are now treated with TNFa inhibitors. However, genomic alterations that occur as a result of TNF-mediated inflammation are not well understood. To investigate molecular targets and networks resulting from increased TNFa, we measured DNA methylation and gene expression in 40 human umbilical vein endothelial cell primary cell lines before and 24 hours after stimulation with TNFa via microarray. Weighted gene co-expression network analysis identified 15 gene groups (modules) with similar expression correlation patterns; four modules showed a strong association with TNFa treatment. Genes in the top TNFa-associated module were all up-regulated, had the highest proportion of hypomethylated regions, and were associated with 136 Disease Ontology terms, including autoimmune/inflammatory, infectious and cardiovascular diseases, and cancers. They included chemokines CXCL1, CXCL10 and CXCL8, and genes associated with autoimmune diseases including HLA-C, DDX58, IL4, NFKBIA and TNFAIP3. Cardiovascular and metabolic disease genes, including APOC1, ACLY, ELOVL6, FASN and SCD, were overrepresented in a module that was not associated with TNFa treatment. Of 223 hypomethylated regions identified, several were in promoters of autoimmune disease GWAS loci (ARID5B, CD69, HDAC9, IL7R, TNIP1 and TRAF1). Results reveal specific gene groups acting in concert in endothelial cells, delineate those driven by TNFa, and establish their relationship to DNA methylation changes, which has strong implications for understanding disease etiology and precision medicine approaches.