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.

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.

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.

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.

Hypomethylation of CYP2E1 and DUSP22 Promoters Associated With Disease Activity and Erosive Disease Among Rheumatoid Arthritis Patients

Objective

Epigenetic modifications have previously been associated with rheumatoid arthritis (RA). In this study, we aimed to determine whether differential DNA methylation in peripheral blood cell subpopulations is associated with any of 4 clinical outcomes among RA patients.

Methods

Peripheral blood samples were obtained from 63 patients in the University of California, San Francisco RA cohort (all satisfied the American College of Rheumatology classification criteria; 57 were seropositive for rheumatoid factor and/or anti‐cyclic citrullinated protein). Fluorescence‐activated cell sorting was used to separate the cells into 4 immune cell subpopulations (CD14+ monocytes, CD19+ B cells, CD4+ naive T cells, and CD4+ memory T cells) per individual, and 229 epigenome‐wide DNA methylation profiles were generated using Illumina HumanMethylation450 BeadChips. Differentially methylated positions and regions associated with the Clinical Disease Activity Index score, erosive disease, RA Articular Damage score, Sharp score, medication at time of blood draw, smoking status, and disease duration were identified using robust regression models and empirical Bayes variance estimators.

Results

Differential methylation of CpG sites associated with clinical outcomes was observed in all 4 cell types. Hypomethylated regions in the CYP2E1 and DUSP22 gene promoters were associated with active and erosive disease, respectively. Pathway analyses suggested that the biologic mechanisms underlying each clinical outcome are cell type–specific. Evidence of independent effects on DNA methylation from smoking, medication use, and disease duration were also identified.

Conclusion

Methylation signatures specific to RA clinical outcomes may have utility as biomarkers or predictors of exposure, disease progression, and disease severity.

Rheumatoid Arthritis Naive T Cells Share Hypermethylation Sites With Synoviocytes

Objective

To determine whether differentially methylated CpGs in synovium‐derived fibroblast‐like synoviocytes (FLS) of patients with rheumatoid arthritis (RA) were also differentially methylated in RA peripheral blood (PB) samples.

Methods

For this study, 371 genome‐wide DNA methylation profiles were measured using Illumina HumanMethylation450 BeadChips in PB samples from 63 patients with RA and 31 unaffected control subjects, specifically in the cell subsets of CD14+ monocytes, CD19+ B cells, CD4+ memory T cells, and CD4+ naive T cells.

Results

Of 5,532 hypermethylated FLS candidate CpGs, 1,056 were hypermethylated in CD4+ naive T cells from RA PB compared to control PB. In analyses of a second set of CpG candidates based on single‐nucleotide polymorphisms from a genome‐wide association study of RA, 1 significantly hypermethylated CpG in CD4+ memory T cells and 18 significant CpGs (6 hypomethylated, 12 hypermethylated) in CD4+ naive T cells were found. A prediction score based on the hypermethylated FLS candidates had an area under the curve of 0.73 for association with RA case status, which compared favorably to the association of RA with the HLA–DRB1 shared epitope risk allele and with a validated RA genetic risk score.

Conclusion

FLS‐representative DNA methylation signatures derived from the PB may prove to be valuable biomarkers for the risk of RA or for disease status.

Genome-wide profiling identifies associations between lupus nephritis and differential methylation of genes regulating tissue hypoxia and type 1 interferon responses

OBJECTIVE

Previous studies have shown that differential DNA methylation is associated with SLE susceptibility. How DNA methylation affects the clinical heterogeneity of SLE has not been fully defined. We conducted this study to identify differentially methylated CpG sites associated with nephritis among women with SLE.

METHODS

The methylation status of 428 229 CpG sites across the genome was characterised for peripheral blood cells from 322 women of European descent with SLE, 80 of whom had lupus nephritis, using the Illumina HumanMethylation450 BeadChip. Multivariable linear regression adjusting for population substructure and leucocyte cell proportions was used to identify differentially methylated sites associated with lupus nephritis. The influence of genetic variation on methylation status was investigated using data from a genome-wide association study of SLE. Pathway analyses were used to identify biological processes associated with lupus nephritis.

RESULTS

We identified differential methylation of 19 sites in 18 genomic regions that was associated with nephritis among patients with SLE (false discovery rate q<0.05). Associations for four sites in HIF3A, IFI44 and PRR4 were replicated when examining methylation data derived from CD4+ T cells collected from an independent set of patients with SLE. These associations were not driven by genetic variation within or around the genomic regions. In addition, genes associated with lupus nephritis in a prior genome-wide association study were not differentially methylated in this epigenome-wide study. Pathway analysis indicated that biological processes involving type 1 interferon responses and the development of the immune system were associated with nephritis in patients with SLE.

CONCLUSIONS

Differential methylation of genes regulating the response to tissue hypoxia and interferon-mediated signalling is associated with lupus nephritis among women with SLE. These findings have not been identified in genetic studies of lupus nephritis, suggesting that epigenome-wide association studies can help identify the genomic differences that underlie the clinical heterogeneity of SLE.

Genome-Wide Assessment of Differential DNA Methylation Associated with Autoantibody Production in Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is characterized by the development of autoantibodies associated with specific clinical manifestations. Previous studies have shown an association between differential DNA methylation and SLE susceptibility, but have not investigated SLE-related autoantibodies. Our goal was to determine whether DNA methylation is associated with production of clinically relevant SLE-related autoantibodies, with an emphasis on the anti-dsDNA autoantibody. In this study, we characterized the methylation status of 467,314 CpG sites in 326 women with SLE. Using a discovery and replication study design, we identified and replicated significant associations between anti-dsDNA autoantibody production and the methylation status of 16 CpG sites (p_discovery<1.07E-07 and p_replication<0.0029) in 11 genes. Associations were further investigated using multivariable regression to adjust for estimated leukocyte cell proportions and population substructure. The adjusted mean DNA methylation difference between anti-dsDNA positive and negative cases ranged from 1.2% to 19%, and the adjusted odds ratio for anti-dsDNA autoantibody production comparing the lowest and highest methylation tertiles ranged from 6.8 to 18.2. Differential methylation for these CpG sites was also associated with anti-SSA, anti-Sm, and anti-RNP autoantibody production. Overall, associated CpG sites were hypomethylated in autoantibody positive compared to autoantibody negative cases. Differential methylation of CpG sites within the major histocompatibility region was not strongly associated with autoantibody production. Genes with differentially methylated CpG sites represent multiple biologic pathways, and have not been associated with autoantibody production in genetic association studies. In conclusion, hypomethylation of CpG sites within genes from different pathways is associated with anti-dsDNA, anti-SSA, anti-Sm, and anti-RNP production in SLE, and these associations are not explained by genetic variation. Thus, studies of epigenetic mechanisms such as DNA methylation represent a complementary method to genetic association studies to identify biologic pathways that may contribute to the clinical heterogeneity of autoimmune diseases.

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×10⁻¹⁶). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53×10⁻¹²), *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×10⁻¹³. 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.

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and involvement of multiple organ systems. Although the cause of SLE remains unknown, several lines of evidence underscore the importance of genetic factors. As is true for most autoimmune diseases, a substantial genetic contribution to disease risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6. This region of the genome contains a large number of genes that participate in the immune response. However, the full contribution of this genomic region to SLE risk has not yet been defined. In the current study we characterize a large number of SLE patients and family members for approximately 2,000 MHC region variants to identify the specific genes that influence disease risk. Our results, for the first time, implicate four different MHC regions in 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.

AUXIN RESPONSE FACTOR 2 (ARF2): a pleiotropic developmental regulator

AUXIN RESPONSE FACTORS (ARFs) regulate auxin-mediated transcriptional activation/repression. They are encoded by a gene family in Arabidopsis, and each member is thought to play a central role in various auxin-mediated developmental processes. We have characterized three arf2 mutant alleles, arf2-6, arf2-7 and arf2-8. The mutants exhibit pleiotropic developmental phenotypes, including large, dark green rosette leaves, delayed flowering, thick and long inflorescence, abnormal flower morphology and sterility in early formed flowers, large organ size and delayed senescence and abscission, compared with wild-type plants. In addition, arf2 mutant seedlings have elongated hypocotyls with enlarged cotyledons under various light conditions. The transcription of ACS2, ACS6 and ACS8 genes is impaired in the developing siliques of arf2-6. The phenotypes of all three alleles are similar to those of the loss-of-function mutants obtained by RNA interference or co-suppression. There is no significant effect of the mutation on global auxin-regulated gene expression in young seedlings, suggesting that ARF2 does not participate in auxin signaling at that particular developmental stage of the plant life cycle. Because ARF2 is thought to function as a transcriptional repressor, the prospect arises that its pleiotropic effects may be mediated by negatively modulating the transcription of downstream genes in signaling pathways that are involved in cell growth and senescence.

Empirical analysis of transcriptional activity in the Arabidopsis genome

Functional analysis of a genome requires accurate gene structure information and a complete gene inventory. A dual experimental strategy was used to verify and correct the initial genome sequence annotation of the reference plant Arabidopsis. Sequencing full-length cDNAs and hybridizations using RNA populations from various tissues to a set of high-density oligonucleotide arrays spanning the entire genome allowed the accurate annotation of thousands of gene structures. We identified 5817 novel transcription units, including a substantial amount of antisense gene transcription, and 40 genes within the genetically defined centromeres. This approach resulted in completion of approximately 30% of the Arabidopsis ORFeome as a resource for global functional experimentation of the plant proteome.