Mechanistic Characterization of RASGRP1 Variants Identifies an hnRNP-K-Regulated Transcriptional Enhancer Contributing to SLE Susceptibility

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic component. We recently identified a novel SLE susceptibility locus near RASGRP1, which governs the ERK/MAPK kinase cascade and B-/T-cell differentiation and development. However, precise causal RASGRP1 functional variant(s) and their mechanisms of action in SLE pathogenesis remain undefined. Our goal was to fine-map this locus, prioritize genetic variants likely to be functional, experimentally validate their biochemical mechanisms, and determine the contribution of these SNPs to SLE risk. We performed a meta-analysis across six Asian and European cohorts (9,529 cases; 22,462 controls), followed by in silico bioinformatic and epigenetic analyses to prioritize potentially functional SNPs. We experimentally validated the functional significance and mechanism of action of three SNPs in cultured T-cells. Meta-analysis identified 18 genome-wide significant (p < 5 × 10⁻⁸) SNPs, mostly concentrated in two haplotype blocks, one intronic and the other intergenic. Epigenetic fine-mapping, allelic, eQTL, and imbalance analyses predicted three transcriptional regulatory regions with four SNPs (rs7170151, rs11631591-rs7173565, and rs9920715) prioritized for functional validation. Luciferase reporter assays indicated significant allele-specific enhancer activity for intronic rs7170151 and rs11631591-rs7173565 in T-lymphoid (Jurkat) cells, but not in HEK293 cells. Following up with EMSA, mass spectrometry, and ChIP-qPCR, we detected allele-dependent interactions between heterogeneous nuclear ribonucleoprotein K (hnRNP-K) and rs11631591. Furthermore, inhibition of hnRNP-K in Jurkat and primary T-cells downregulated RASGRP1 and ERK/MAPK signaling. Comprehensive association, bioinformatics, and epigenetic analyses yielded putative functional variants of RASGRP1, which were experimentally validated. Notably, intronic variant (rs11631591) is located in a cell type-specific enhancer sequence, where its risk allele binds to the hnRNP-K protein and modulates RASGRP1 expression in Jurkat and primary T-cells. As risk allele dosage of rs11631591 correlates with increased RASGRP1 expression and ERK activity, we suggest that this SNP may underlie SLE risk at this locus.

Amino acid signatures of HLA Class-I and II molecules are strongly associated with SLE susceptibility and autoantibody production in Eastern Asians

Human leukocyte antigen (HLA) is a key genetic factor conferring risk of systemic lupus erythematosus (SLE), but precise independent localization of HLA effects is extremely challenging. As a result, the contribution of specific HLA alleles and amino-acid residues to the overall risk of SLE and to risk of specific autoantibodies are far from completely understood. Here, we dissected (a) overall SLE association signals across HLA, (b) HLA-peptide interaction, and (c) residue-autoantibody association. Classical alleles, SNPs, and amino-acid residues of eight HLA genes were imputed across 4,915 SLE cases and 13,513 controls from Eastern Asia. We performed association followed by conditional analysis across HLA, assessing both overall SLE risk and risk of autoantibody production. DR15 alleles HLA-DRB1*15:01 (P = 1.4x10^-27, odds ratio (OR) = 1.57) and HLA-DQB1*06:02 (P = 7.4x10^-23, OR = 1.55) formed the most significant haplotype (OR = 2.33). Conditioned protein-residue signals were stronger than allele signals and mapped predominantly to HLA-DRB1 residue 13 (P = 2.2x10^-75) and its proxy position 11 (P = 1.1x10^-67), followed by HLA-DRB1-37 (P = 4.5x10^-24). After conditioning on HLA-DRB1, novel associations at HLA-A-70 (P = 1.4x10^-8), HLA-DPB1-35 (P = 9.0x10^-16), HLA-DQB1-37 (P = 2.7x10^-14), and HLA-B-9 (P = 6.5x10^-15) emerged. Together, these seven residues increased the proportion of explained heritability due to HLA to 2.6%. Risk residues for both overall disease and hallmark autoantibodies (i.e., nRNP: DRB1-11, P = 2.0x10^-14; DRB1-13, P = 2.9x10^-13; DRB1-30, P = 3.9x10^-14) localized to the peptide-binding groove of HLA-DRB1. Enrichment for specific amino-acid characteristics in the peptide-binding groove correlated with overall SLE risk and with autoantibody presence. Risk residues were in primarily negatively charged side-chains, in contrast with rheumatoid arthritis. We identified novel SLE signals in HLA Class I loci (HLA-A, HLA-B), and localized primary Class II signals to five residues in HLA-DRB1, HLA-DPB1, and HLA-DQB1. These findings provide insights about the mechanisms by which the risk residues interact with each other to produce autoantibodies and are involved in SLE pathophysiology.

The Human leukocyte antigen (HLA) region is a key genetic factor conferring risk of systemic lupus erythematosus (SLE). In spite of multiple SLE association signals identified in the HLA region, only amino-acid residues within HLA-DRB1 have been specifically described previously. In this study, we performed an imputation-based analysis on individuals with East Asian ancestry, and characterized SLE risk within the HLA region for all involved independent genes (HLA-DRB1, HLA-DPB1, HLA-DQB1, HLA-A, and HLA-B). Furthermore, we identified a characteristic SLE risk residue signature as well as a pattern of specific nRNP and Ro/La autoantibody residues located in the peptide-binding grooves, suggesting their key involvement in autoantibody production.

Hormetic Effects of Thiophanate-Methyl in Multiple Isolates of Sclerotinia homoeocarpa

Twenty-eight isolates of Sclerotinia homoeocarpa, causal agent of dollar spot disease in turf, were assessed for fungicide hormesis at sublethal concentrations of thiophanate-methyl (T-methyl). Each isolate was grown in corn meal agar amended with 11 concentrations of T-methyl (30,500 to 0.047 µg/liter), and the area of mycelial growth was determined relative to the control. Three replicates were used per concentration, and the experiment was repeated three to five times for each isolate. Reference isolates (EC₅₀ > 20 µg/liter), with no prior history of T-methyl exposure, were highly sensitive and not stimulated by low doses. Likewise, no stimulation was observed in two highly sensitive isolates (EC₅₀ > 30 µg/liter) that had been preconditioned by exposure to T-methyl, or in four T-methyl-tolerant isolates. Seventeen (81%) preconditioned T-methyl-tolerant isolates (EC₅₀ = 294 to1,550 µg/liter) had statistically significant growth stimulation, in the range of 2.8 to 19.7% relative to the control. These results support that hormesis (low-dose stimulation, high-dose inhibition) is a common dose response in preconditioned S. homoeocarpa, particularly in response to subtoxic doses of T-methyl.

Confirmation of five novel susceptibility loci for systemic lupus erythematosus (SLE) and integrated network analysis of 82 SLE susceptibility loci

We recently identified ten novel SLE susceptibility loci in Asians and uncovered several additional suggestive loci requiring further validation. This study aimed to replicate five of these suggestive loci in a Han Chinese cohort from Hong Kong, followed by meta-analysis (11,656 cases and 23,968 controls) on previously reported Asian and European populations, and to perform bioinformatic analyses on all 82 reported SLE loci to identify shared regulatory signatures. We performed a battery of analyses for these five loci, as well as joint analyses on all 82 SLE loci. All five loci passed genome-wide significance: MYNN (rs10936599, Pmeta = 1.92 × 10-13, OR = 1.14), ATG16L2 (rs11235604, Pmeta = 8.87 × 10 -12, OR = 0.78), CCL22 (rs223881, Pmeta = 5.87 × 10-16, OR = 0.87), ANKS1A (rs2762340, Pmeta = 4.93 × 10-15, OR = 0.87) and RNASEH2C (rs1308020, Pmeta = 2.96 × 10-19, OR = 0.84) and co-located with annotated gene regulatory elements. The novel loci share genetic signatures with other reported SLE loci, including effects on gene expression, transcription factor binding, and epigenetic characteristics. Most (56%) of the correlated (r2 > 0.8) SNPs from the 82 SLE loci were implicated in differential expression (9.81 × 10-198 < P < 5 × 10-3) of cis-genes. Transcription factor binding sites for p53, MEF2A and E2F1 were significantly (P < 0.05) over-represented in SLE loci, consistent with apoptosis playing a critical role in SLE. Enrichment analysis revealed common pathways, gene ontology, protein domains, and cell type-specific expression. In summary, we provide evidence of five novel SLE susceptibility loci. Integrated bioinformatics using all 82 loci revealed that SLE susceptibility loci share many gene regulatory features, suggestive of conserved mechanisms of SLE etiopathogenesis.

Allelic heterogeneity in NCF2 associated with systemic lupus erythematosus (SLE) susceptibility across four ethnic populations

Recent reports have associated NCF2, encoding a core component of the multi-protein NADPH oxidase (NADPHO), with systemic lupus erythematosus (SLE) susceptibility in individuals of European ancestry. To identify ethnicity-specific and -robust variants within NCF2, we assessed 145 SNPs in and around the NCF2 gene in 5325 cases and 21 866 controls of European-American (EA), African-American (AA), Hispanic (HS) and Korean (KR) ancestry. Subsequent imputation, conditional, haplotype and bioinformatic analyses identified seven potentially functional SLE-predisposing variants. Association with non-synonymous rs17849502, previously reported in EA, was detected in EA, HS and AA (P(EA) = 1.01 × 10(-54), PHS = 3.68 × 10(-10), P(AA) = 0.03); synonymous rs17849501 was similarly significant. These SNPs were monomorphic in KR. Novel associations were detected with coding variants at rs35937854 in AA (PAA = 1.49 × 10(-9)), and rs13306575 in HS and KR (P(HS) = 7.04 × 10(-7), P(KR) = 3.30 × 10(-3)). In KR, a 3-SNP haplotype was significantly associated (P = 4.20 × 10(-7)), implying that SLE predisposing variants were tagged. Significant SNP-SNP interaction (P = 0.02) was detected between rs13306575 and rs17849502 in HS, and a dramatically increased risk (OR = 6.55) with a risk allele at each locus. Molecular modeling predicts that these non-synonymous mutations could disrupt NADPHO complex assembly. The risk allele of rs17849501, located in a conserved transcriptional regulatory region, increased reporter gene activity, suggesting in vivo enhancer function. Our results not only establish allelic heterogeneity within NCF2 associated with SLE, but also emphasize the utility of multi-ethnic cohorts to identify predisposing variants explaining additional phenotypic variance ('missing heritability') of complex diseases like SLE.

Predicting fusarium head blight epidemics with weather-driven pre- and post-anthesis logistic regression models

Our objective was to identify weather-based variables in pre- and post-anthesis time windows for predicting major Fusarium head blight (FHB) epidemics (defined as FHB severity ≥ 10%) in the United States. A binary indicator of major epidemics for 527 unique observations (31% of which were major epidemics) was linked to 380 predictor variables summarizing temperature, relative humidity, and rainfall in 5-, 7-, 10-, 14-, or 15-day-long windows either pre- or post-anthesis. Logistic regression models were built with a training data set (70% of the 527 observations) using the leaps-and-bounds algorithm, coupled with bootstrap variable and model selection methods. Misclassification rates were estimated on the training and remaining (test) data. The predictive performance of models with indicator variables for cultivar resistance, wheat type (spring or winter), and corn residue presence was improved by adding up to four weather-based predictors. Because weather variables were intercorrelated, no single model or subset of predictor variables was best based on accuracy, model fit, and complexity. Weather-based predictors in the 15 final empirical models selected were all derivatives of relative humidity or temperature, except for one rainfall-based predictor, suggesting that relative humidity was better at characterizing moisture effects on FHB than other variables. The average test misclassification rate of the final models was 19% lower than that of models currently used in a national FHB prediction system.

Admixture mapping in lupus identifies multiple functional variants within IFIH1 associated with apoptosis, inflammation, and autoantibody production

Systemic lupus erythematosus (SLE) is an inflammatory autoimmune disease with a strong genetic component. African-Americans (AA) are at increased risk of SLE, but the genetic basis of this risk is largely unknown. To identify causal variants in SLE loci in AA, we performed admixture mapping followed by fine mapping in AA and European-Americans (EA). Through genome-wide admixture mapping in AA, we identified a strong SLE susceptibility locus at 2q22–24 (LOD = 6.28), and the admixture signal is associated with the European ancestry (ancestry risk ratio ∼1.5). Large-scale genotypic analysis on 19,726 individuals of African and European ancestry revealed three independently associated variants in the IFIH1 gene: an intronic variant, rs13023380 [P_meta = 5.20×10⁻¹⁴; odds ratio, 95% confidence interval = 0.82 (0.78–0.87)], and two missense variants, rs1990760 (Ala946Thr) [P_meta = 3.08×10⁻⁷; 0.88 (0.84–0.93)] and rs10930046 (Arg460His) [P_dom = 1.16×10⁻⁸; 0.70 (0.62–0.79)]. Both missense variants produced dramatic phenotypic changes in apoptosis and inflammation-related gene expression. We experimentally validated function of the intronic SNP by DNA electrophoresis, protein identification, and in vitro protein binding assays. DNA carrying the intronic risk allele rs13023380 showed reduced binding efficiency to a cellular protein complex including nucleolin and lupus autoantigen Ku70/80, and showed reduced transcriptional activity in vivo. Thus, in SLE patients, genetic susceptibility could create a biochemical imbalance that dysregulates nucleolin, Ku70/80, or other nucleic acid regulatory proteins. This could promote antibody hypermutation and auto-antibody generation, further destabilizing the cellular network. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

African-Americans (AA) are at increased risk of systemic lupus erythematosus (SLE), but the genetic basis of this risk increase is largely unknown. We used admixture mapping to localize disease-causing genetic variants that differ in frequency across populations. This approach is advantageous for localizing susceptibility genes in recently admixed populations like AA. Our genome-wide admixture scan identified seven admixture signals, and we followed the best signal at 2q22–24 with fine-mapping, imputation-based association analysis and experimental validation. We identified two independent coding variants and a non-coding variant within the IFIH1 gene associated with SLE. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

Sublethal Doses of Mefenoxam Enhance Pythium Damping-off of Geranium

The effect of sublethal doses of fungicides on fungicide-resistant Pythium isolates is unknown but potentially relevant to disease management. Occasional grower reports of Pythium disease increases after fungicide applications and our observations of greater radial growth in vitro on fungicide-amended media than on nonamended media suggests that Pythium isolates may be stimulated by sublethal doses of fungicides. The objectives of this study were to determine whether Pythium isolates were stimulated by sublethal doses of mefenoxam in vitro and whether this stimulation had any influence on Pythium damping-off of geranium seedlings. A mefenoxam-resistant isolate of Pythium aphanidermatum displayed 10% mean radial growth increase in vitro with mefenoxam at 1 × 10^-10 μg/ml compared with growth on nonamended agar (nonsignificant). Geranium seedlings treated with one of eight mefenoxam concentrations were inoculated with 5-mm-diameter colonized agar plugs and evaluated for disease severity every 24 h. The area under the disease progress curve and the survival curve were estimated for each treatment and compared. Significant increases in damping-off were observed with mefenoxam at 1 × 10^-6 and 1 × 10^-10 μg/ml. Our data indicate that a Pythium isolate with resistance to mefenoxam can be stimulated by sublethal doses of the fungicide, and that this stimulation can result in significantly higher rates of Pythium damping-off of geranium seedlings.

Bean pod mottle virus Spread in Insect-Feeding-Resistant Soybean

Bean pod mottle virus (BPMV) infection reduces yield and seed quality in soybean. To test the hypothesis that virus incidence and movement within plots would be reduced in soybean with resistance to feeding by the virus' bean leaf beetle (Cerotoma trifurcata) vector, BPMV spread was evaluated in five soybean genotypes at two inoculum levels over 2 years at two locations in Ohio. Soybean genotypes included two insect-feeding-susceptible genotypes (Williams 82 and Resnik), two insect-feeding-resistant, semidwarf genotypes (HC95-15 and HC95-24), and an insect-feeding-susceptible, semidwarf genotype (Troll). BPMV incidence was assessed in individual plants at growth stages R5/R6 and R7/R8 using enzyme-linked immunosorbent assay. Beetle feeding was visually assessed in 2004. Data for infection of individual plants were analyzed using a generalized linear mixed model, with a binomial distribution and logit-link. Within plots, BPMV incidence was highest in Resnik and Williams 82 and significantly lower in Troll. Incidence in HC95-15 was not significantly different than in Williams 82 and Resnik but incidence in HC95-24 was lower than in Resnik. BPMV incidence was also significantly (P < 0.05) affected by year, location, inoculum level and sampling date, with increasing incidence over time and higher incidence at the higher inoculum level. Beetle feeding damage was affected by the interaction of location-genotype. Significant spatial aggregation of infected plants was found for most plots but aggregation was independent of host genotype and inoculum level. Although the results indicate that BPMV infection varied by genotype, they do not support the hypothesis that insect-feeding resistance is sufficient to reduce the incidence and spread of BPMV.