Epistatic Interaction between BANK1 and BLK in Rheumatoid Arthritis: Results from a Large Trans-Ethnic Meta-Analysis

Background

BANK1 and BLK belong to the pleiotropic autoimmune genes; recently, epistasis between BANK1 and BLK was detected in systemic lupus erythematosus. Although BLK has been reproducibly identified as a risk factor in rheumatoid arthritis (RA), reports are conflicting about the contribution of BANK1 to RA susceptibility. To ascertain the real impact of BANK1 on RA genetic susceptibility, we performed a large meta-analysis including our original data and tested for an epistatic interaction between BANK1 and BLK in RA susceptibility.

Patients and Methods

We investigated data for 1,915 RA patients and 1,915 ethnically matched healthy controls genotyped for BANK1 rs10516487 and rs3733197 and BLK rs13277113. The association of each SNP and RA was tested by logistic regression. Multivariate analysis was then used with an interaction term to test for an epistatic interaction between the SNPs in the 2 genes.

Results

None of the SNPs tested individually was significantly associated with RA in the genotyped samples. However, we detected an epistatic interaction between BANK1 rs3733197 and BLK rs13277113 (P_interaction = 0.037). In individuals carrying the BLK rs13277113 GG genotype, presence of the BANK1 rs3733197 G allele increased the risk of RA (odds ratio 1.21 [95% confidence interval 1.04–1.41], P = 0.015. Combining our results with those of all other studies in a large trans-ethnic meta-analysis revealed an association of the BANK1 rs3733197 G allele and RA (1.11 [1.02–1.21], P = 0.012).

Conclusion

This study confirms BANK1 as an RA susceptibility gene and for the first time provides evidence for epistasis between BANK1 and BLK in RA. Our results illustrate the concept of pleiotropic epistatic interaction, suggesting that BANK1 and BLK might play a role in RA pathogenesis.     

Relationship between Item Responses of Negative Affect Items and the Distribution of the Sum of the Item Scores in the General Population

BackgroundSeveral studies have shown that total depressive symptom scores in the general population approximate an exponential pattern, except for the lower end of the distribution. The Center for Epidemiologic Studies Depression Scale (CES-D) consists of 20 items, each of which may take on four scores: “rarely,” “some,” “occasionally,” and “most of the time.” Recently, we reported that the item responses for 16 negative affect items commonly exhibit exponential patterns, except for the level of “rarely,” leading us to hypothesize that the item responses at the level of “rarely” may be related to the non-exponential pattern typical of the lower end of the distribution. To verify this hypothesis, we investigated how the item responses contribute to the distribution of the sum of the item scores.MethodsData collected from 21,040 subjects who had completed the CES-D questionnaire as part of a Japanese national survey were analyzed. To assess the item responses of negative affect items, we used a parameter r, which denotes the ratio of “rarely” to “some” in each item response. The distributions of the sum of negative affect items in various combinations were analyzed using log-normal scales and curve fitting.ResultsThe sum of the item scores approximated an exponential pattern regardless of the combination of items, whereas, at the lower end of the distributions, there was a clear divergence between the actual data and the predicted exponential pattern. At the lower end of the distributions, the sum of the item scores with high values of r exhibited higher scores compared to those predicted from the exponential pattern, whereas the sum of the item scores with low values of r exhibited lower scores compared to those predicted.ConclusionsThe distributional pattern of the sum of the item scores could be predicted from the item responses of such items.     

Human Leukocyte Antigen and Systemic Sclerosis in Japanese: The Sign of the Four Independent Protective Alleles,DRB1*13:02, DRB1*14:06, DQB1*03:01, and DPB1*02:01

ObjectiveSeveral studies on associations between human leukocyte antigen (HLA) allele frequencies and susceptibility to systemic sclerosis (SSc) have been reported. Anti-centromere antibodies (ACA) and anti-topoisomerase I antibodies (ATA) are found in SSc patients. Here, we sought to identify HLA alleles associated with SSc in Japanese, and explored their associations with SSc phenotypes including the presence of autoantibodies.MethodsAssociations of HLA-DRB1, DQB1, and DPB1 were analyzed in 463 Japanese SSc patients and 413 controls.ResultsWe found that DRB1*13:02 (P = 0.0011, Pc = 0.0319, odds ratio [OR] 0.46, 95% confidence interval [CI] 0.29–0.73), DRB1*14:06 (P = 6.60X10^-5, Pc = 0.0020, OR 0.05, 95%CI 0.01–0.41), DQB1*03:01 (P = 0.0009, Pc = 0.0150, OR 0.56, 95%CI 0.40–0.79), and DPB1*02:01 (P = 5.16X10^-6, Pc = 8.77X10^-5, OR 0.52, 95%CI 0.39–0.69) were protectively associated with SSc. In addition, these four alleles seemed to be independently associated with the protection against the susceptibility of SSc. On the other hand, we could not find predisposing alleles for overall SSc. With respect to SSc subsets, a tendency for these four alleles to be protectively associated was observed. However, there was a significant association between DRB1*01:01, DRB1*10:01, DQB1*05:01, and DPB1*04:02 and the susceptibility to SSc with ACA. On the other hand, the presence of DRB1*15:02, DQB1*06:01, DPB1*03:01, and DPB1*09:01 was associated with SSc with ATA.ConclusionThus, the present study has identified protective associations of the four HLA class II alleles with overall Japanese SSc and predisposing associations of HLA class II alleles with Japanese SSc subsets.     

Heme Orientation of Cavity Mutant Hemoglobins (His F8 → Gly) in Either α or β Subunits: Circular Dichroism, 1H NMR,and Resonance Raman Studies

Native human adult hemoglobin (Hb A) has mostly normal orientation of heme, whereas recombinant Hb A (rHb A) expressed in E. coli contains both normal and reversed orientations of heme. Hb A with the normal heme exhibits positive circular dichroism (CD) bands at both the Soret and 260‐nm regions, while rHb A with the reversed heme shows a negative Soret and decreased 260‐nm CD bands. In order to examine involvement of the proximal histidine (His F8) of either α or β subunits in determining the heme orientation, we prepared two cavity mutant Hbs, rHb(αH87G) and rHb(βH92G), with substitution of glycine for His F8 in the presence of imidazole. CD spectra of both cavity mutant Hbs did not show a negative Soret band, but instead exhibited positive bands with strong intensity at the both Soret and 260‐nm regions, suggesting that the reversed heme scarcely exists in the cavity mutant Hbs. We confirmed by ¹H NMR and resonance Raman (RR) spectroscopies that the cavity mutant Hbs have mainly the normal heme orientation in both the mutated and native subunits. These results indicate that the heme Fe‐His F8 linkage in both α and β subunits influences the heme orientation, and that the heme orientation of one type of subunit is related to the heme orientation of the complementary subunits to be the same. The present study showed that CD and RR spectroscopies also provided powerful tools for the examination of the heme rotational disorder of Hb A, in addition to the usual ¹H NMR technique. Chirality 28:585–592, 2016. © 2016 Wiley Periodicals, Inc.     

Crystal structure of unautoprocessed precursor of subtilisin from a hyperthermophilic archaeon: evidence for Ca2+-induced folding

The crystal structure of an active site mutant of pro-Tk-subtilisin (pro-S324A) from the hyperthermophilic archaeon Thermococcus kodakaraensis was determined at 2.3 A resolution. The overall structure of this protein is similar to those of bacterial subtilisin-propeptide complexes, except that the peptide bond linking the propeptide and mature domain contacts with the active site, and the mature domain contains six Ca2+ binding sites. The Ca-1 site is conserved in bacterial subtilisins but is formed prior to autoprocessing, unlike the corresponding sites of bacterial subtilisins. All other Ca2+-binding sites are unique in the pro-S324A structure and are located at the surface loops. Four of them apparently contribute to the stability of the central alphabetaalpha substructure of the mature domain. The CD spectra, 1-anilino-8-naphthalenesulfonic acid fluorescence spectra, and sensitivities to chymotryptic digestion of this protein indicate that the conformation of pro-S324A is changed from an unstable molten globule-like structure to a stable native one upon Ca2+ binding. Another active site mutant, pro-S324C, was shown to be autoprocessed to form a propeptide-mature domain complex in the presence of Ca2+. The CD spectra of this protein indicate that the structure of pro-S324C is changed upon Ca2+ binding like pro-S324A but is not seriously changed upon subsequent autoprocessing. These results suggest that the maturation process of Tk-subtilisin is different from that of bacterial subtilisins in terms of the requirement of Ca2+ for folding of the mature domain and completion of the folding process prior to autoprocessing.     

Protein thermostabilization requires a fine-tuned placement of surface-charged residues

Using the information from the genome projects, recent comparative studies of thermostable proteins have revealed a certain trend of amino acid composition in which polar residues are scarce and charged residues are rich on the protein surface. To clarify experimentally the effect of the amino acid composition of surface residues on the thermostability of Escherichia coli Ribonuclease HI (RNase HI), we constructed six variants in which five to eleven polar residues were replaced by charged residues (5C, 7Ca, 7Cb, 9Ca, 9Cb and 11C). The thermal denaturation experiments indicated that all of the variant proteins are 3.2-10.1 degrees C in Tm less stable than the wild proteins. The crystal structures of resultant protein variants 7Ca, 7Cb, 9Ca and 11C closely resemble that of E. coli RNase HI in their global fold, and several different hydrogen bonding and ion-pair interactions are formed by the mutations. Comparison of the crystal structures of these variant proteins with that of E. coli RNase HI reveals that thermal destabilization is apparently related to electrostatic repulsion of the charged residues with neighbours. This result suggests that charged residues of natural thermostable proteins are strictly posted on the surface with optimal interactions and without repulsive interactions.     

Helicobacter pylori dampens gut epithelial self-renewal by inhibiting apoptosis, a bacterial strategy to enhance colonization of the stomach

Colonization of the gastric pits in the stomach by Helicobacter pylori (Hp) is a major risk factor for gastritis, gastric ulcers, and cancer. Normally, rapid self-renewal of gut epithelia, which occurs by a balance of progenitor proliferation and pit cell apoptosis, serves as a host defense mechanism to limit bacterial colonization. To investigate how Hp overcomes this host defense, we use the Mongolian gerbil model of Hp infection. Apoptotic loss of pit cells induced by a proapoptotic agent is suppressed by Hp. The ability of Hp to suppress apoptosis contributed to pit hyperplasia and persistent bacterial colonization of the stomach. Infection with WT Hp but not with a mutant in the virulence effector cagA increased levels of the prosurvival factor phospho-ERK and antiapoptotic protein MCL1 in the gastric pits. Thus, CagA activates host cell survival and antiapoptotic pathways to overcome self-renewal of the gastric epithelium and help sustain Hp infection.     

The genetic structure of endangered indigenous populations of the amago salmon, <Emphasis Type="Italic">Oncorhynchus</Emphasis><Emphasis Type="Italic">masou</Emphasis><Emphasis Type="Italic">ishikawae</Emphasis>, in Japan

The amago salmon, Oncorhynchus masou ishikawae, is an endemic subspecies of O. masou in Japan. Owing to the extensive stocking of hatchery fish throughout Japan, indigenous populations of O. m. ishikawae are now on the verge of extinction. We examined the genetic effects of stocking hatchery fish on wild populations in the River Koza, Japan, using microsatellite and mitochondrial DNA (mtDNA) markers. For mtDNA, haplotype mt1, which is common in wild populations, was present exclusively in isolated wild populations assumed to be unaffected by previous stocking, while it was never observed in hatchery fish. Genetic diversity was much higher in wild populations in the stocked area, which shared many mtDNA haplotypes with hatchery fish, than in isolated wild populations with haplotype mt1. Pairwise F _ST estimates based on microsatellites showed significant differentiation among the isolated populations with many microsatellite loci monomorphic. Significant deviation from Hardy–Weinberg equilibrium was observed in wild populations in the area subject to stocking, where a Bayesian-based assignment test showed a high level of introgression with hatchery fish. These results suggest that wild populations with haplotype mt1, which became isolated through anthropogenic environmental change in the 1950–1960s, represent indigenous populations of O. m. ishikawae in the River Koza. They have low genetic diversity, most likely caused by genetic bottlenecks following damming and environmental deterioration, while stocking of hatchery fish over the past 30 years apparently had a large impact on the genetic structure of wild populations in the main channel of the River Koza.