Displayed in Figure are all genes that were
significantly upregulated or downregulated (at least 2-fold difference, P
< 0.001) in the strains that were highly sensitive to AIGN disease
(i.e., NZW, DBA/1, and 129/SvJ) compared with either of the control strains (C57BL/6
[B6] or BALB/c), following challenge with anti-GBM serum. Though several strain-specific
gene differences were also noted within this panel of genes, a subset of 50 genes within
this panel were consistently downregulated in all 3 of the highly disease-sensitive
strains compared with the control strains (shown enlarged on the right, Figure A). Intriguingly, 10 of these genes belonged to the
kallikrein (Klk) family, with the highest expression levels and
differences being noted in Klk1, Klk1b3,
Klk1b5, Klk1b26, and Klk1b27, as
summarized in Table . In addition to the
Klk genes displayed in Table (all of which were found to be different between the 2 sets of strains at
P < 0.001), a few other Klk genes, notably
Klk1b9 and Klk1b21, exhibited similar expression
differences between the disease-sensitive strains and the control strains, though these
differences did not reach statistical P values of 0.001 (data not
shown). Hence, in total, 12 Klk genes were underexpressed in the
kidneys of AIGN-sensitive strains, following anti-GBM challenge, compared with the
control strains.
| Table 1 Several kallikrein gene messages were significantly upregulated in the kidneys of
B6 and BALB/c mice compared with NZW, DBA1, and 129/SvJ mice, upon anti-GBM Ab
challenge |
Next, renal Klk gene expression was examined before and after induction
of AIGN, using real-time PCR as an orthogonal approach. Real-time PCR analyses validated
the above microarray results, indicating that whereas the B6 and BALB/c control strains
successfully upregulated Klk following anti-GBM Ab challenge, the
DBA/1, 129/SvJ, and NZW strains were ineffective at doing so (Figure B). In contrast, the basal, predisease levels of
renal Klk were similar in all strains (Figure B). These differences were also confirmed at the protein level by
Western blot analysis of renal cortex samples from these 5 strains, as illustrated for
Klk1 (Figure A and data not
shown). Parallel differences in Klk enzymatic activity were also noted in urine samples
from the same mice following anti-GBM Ab challenge (Figure B). Thus, whereas the BALB/c control strain exhibited a robust increase in
urinary Klk activity 10–14 days after challenge, in the AIGN-sensitive
strains DBA/1 and 129/SvJ, Klk was not significantly upregulated over the same period
(Figure B).
Given that Klk upregulation in the kidneys and urine of anti-GBM
Ab–challenged mice correlated well with the subdued nephritis noted in the
B6 and BALB/c control strains, we next asked whether Klk might have a disease-protective
role in immune nephritis. Kallikreins act through the generation of bradykinins (BKs),
which in turn exert their biological effects by binding BK (BK B1 and B2) receptors on
various cells (
25,
26). Selective receptor blockade using pharmacological inhibitors
further indicated that the biological effects of BK (and Klk) were mediated by the BK B2
receptor, since blocking this receptor aggravated proteinuria, azotemia, and GN
following anti-GBM Ab challenge in BALB/c mice (Figure , A–D). Importantly, the mice in which BK B2 receptor was blocked
exhibited significantly more severe GN compared with the other groups of mice (average
GN score of 2.1 versus 0.2,
P < 0.001, as partly illustrated in
Figure , C and D). Similar differences were noted
when the BK B2 receptor blockers were administered to anti-GBM–challenged B6
mice (data not shown). Conversely, the administration of BK dampened the severity of
anti-GBM disease in 129/SvJ mice, which otherwise develop severe AIGN following the
experimental insult (Figure E and data not shown).
With respect to the 3 AIGN-sensitive strains, DBA/1, 129/SvJ, and NZW, it is already
known that several genomic intervals in the latter 2 strains (including loci on
chromosomes 1 and 7) also contribute to spontaneous lupus nephritis (
13–
16). Indeed, the entire
Klk gene complex is encoded within a
lupus susceptibility interval on chromosome 7. The NZW-derived lupus susceptibility
locus on chromosome 7,
Sle3z (which includes the
Klk gene complex), has previously been introgressed onto the B6 genome
as a congenic interval, and this locus had already been shown to facilitate development
of spontaneous lupus nephritis (
17–
19). Through recursive
backcrossing of B6.
Sle3z congenic mice to B6 parents and
microsatellite-assisted selection, we generated
B6.
Sle3z157–158 recombinants
harboring the NZW-derived
Klk gene complex within a 4-Mb interval, with
termini at
D7mit157 and
D7mit158, as diagramed in
Figure A.
Importantly, B6.Sle3z mice and the newly generated
B6.Sle3z157–158 congenic
recombinants both displayed heightened sensitivity to AIGN, marked by elevated
proteinuria and severe nephritis (Figure ,
B–D). Moreover, the renal cortex of both these congenics failed to
efficiently upregulate Klk following anti-GBM Ab challenge, compared
with the B6 controls (Figure E), as assessed by
real-time PCR. Taken together with the functional data presented above, these findings
suggest that the z allele of Klk positioned within the
Sle3z157–158 subinterval may
harbor important culprit genes for the heightened experimental anti-GBM disease (and
spontaneous lupus nephritis) seen in NZW (and related strains of) mice.
Five of the most differentially expressed Klk genes,
Klk1, Klk1b3, Klk1b5,
Klk1b26, and Klk1b27 were sequenced (GenBank accession
numbers EU597301–EU597324). Klk4 and
Klk1b8, though differentially expressed, were not studied further
because of their relatively low expression levels in all strains (Table ). Several strain-specific differences were noted
in the promoter regions of the 5 sequenced genes, as summarized in Table . The B6 and BALB/c Klk genes were
almost identical, while the Klk genes from the AIGN-sensitive strains
NZW and DBA/1 were closely related to each other, as illustrated by the phylogenetic
trees in Figure A and the promoter region sequence
exemplified in Figure B. The Klk
genes from the 129/SvJ strain clustered more closely with the B6/BALB/c genes in some
cases, and with the NZW/DBA/1 genes in others; despite repeated attempts, we could not
amplify some of the critical regions of the 129/SvJ Klk1b26 gene (Table
). In contrast to the promoter regions, no
sequence differences were noted in the coding regions or the 3′ untranslated
regions of these Klk genes, when the 5 strains were compared (data not
shown).
| Table 2 Sequence polymorphisms in the promoter region of the mouse Klk
genes |
Some of the observed strain-specific promoter sequence differences in
Klk fell within putative transcription factor binding sites. In the
promoter region of Klk1b3, for example, SNPs were noted in potential
transcription factor binding sites for PBF, TFIID, HoxD9, TCF-4E, NF-S, NF-E, and LBP-1,
which distinguished the AIGN-sensitive strains from the controls. In particular, the
GGCTT[C→G]AAAAT SNP in the promoter region of Klk1b3 is
predicted to abrogate TFIID binding. Importantly, promoter-luciferase studies indicated
that sequence differences in the promoter region of Klk1b3 may
contribute in part to the reduced expression of this gene in the disease-sensitive
strains (Figure C). Given that more than 200
sequence variations have been noted in the regulatory regions of the
Klk genes (Table ), each
difference has to be systematically evaluated for its potential functional relevance.
As in murine lupus, the orthologous human interval encompassing the
KLK
genes on human chromosome
19q13 has also been implicated in human SLE
susceptibility in previous genome scans (
40–
43). To ascertain
whether
KLK might also be a culprit gene in human SLE, we examined
several SNPs in the renal-expressed
KLK genes encoded within this
interval, specifically human
KLK1,
KLK5,
KLK6, and
KLK7. We first genotyped a set of German SLE
patients (
n = 340) and a set of controls matched for ethnicity, age,
and sex (
n = 400). As detailed in Table , association to two
KLK1 SNPs,
rs1054713
(a synonymous coding SNP) and
rs2740502 (an intronic SNP), but not to
the nonsynonymous substitution
rs5517, was observed, with the strongest
associations being noted in SLE patients with nephritis (
KLK1 SNP,
rs2740502;
P = 0.007 compared with non-nephritic
SLE patients and
P = 0.01 compared with healthy controls). A weak
association was also observed for SNP
rs1897604 located in
KLK5; however, this SNP was not in Hardy-Weinberg equilibrium in the
controls (data not shown). No disease association was detected for any of the other
KLK genes analyzed in this preliminary study.
| Table 3 Association analysis of KLK SNPs in German SLE patients and controls |
To validate these preliminary associations, we genotyped 6 more cohorts of patients,
including additional European, European-American, and Korean patients with SLE, as
detailed in Tables and . For the KLK1 SNP rs2740502, we
could replicate the genetic association to lupus nephritis compared with controls in the
German and Italian samples but not in the remaining patient sets (Table ). Likewise, the German SLE patients exhibited
association to another KLK1 SNP, rs5516 (Table ). Due to likely heterogeneity in haplotype distribution and
potential differences in linkage disequilibrium between the analyzed SNPs and functional
SNPs in the different ethnic groups, we were unable to perform a Mantel-Haenszel
metaanalysis with these data sets.
| Table 4 Genotypic and allelic association of KLK1 rs5516 in 6 additional
sets of SLE patients with nephritis (cases) and healthy controls |
| Table 5 Genotypic and allelic association of KLK1 rs2740502 in 6
additional sets of SLE patients with nephritis (cases) and healthy controls |
For further independent validation, we gained access to the genotype data for SNPs of
the entire KLK region typed in 689 SLE patients and 3,718 healthy controls from the
International Consortium for Systemic Lupus Erythematosus Genetics (SLEGEN) and 595 SLE
patients from the UCSF Lupus Genetics Project cohort, both of which had been assembled
for genome-wide association studies. Although the sets of SNPs spanning the entire KLK
locus utilized in these studies were not completely overlapping, 56 SNPs were, and these
were examined further for disease association. The strongest association with SLE was
again noted to SNPs close to KLK1, in an intergenic region bordered at
the centromeric end by KLK1 and KLK15 and at the
telomeric end by KLK3 and KLK2 (Figure A and Supplemental Table 1; supplemental material
available online with this article; doi:10.1172/JCI36728DS1). Additional associations to
SNPs in KLK4 promoter, KLK5, KLK7,
KLK11, and KLK12 were also uncovered when cases
with nephritis were compared with cases without nephritis (Figure B and Supplemental Table 2). Next, we defined the haplotype blocks
across the KLK locus (Figure A)
and examined them for association to SLE. Again, the haplotype blocks harboring the
promoter region of KLK3 (block 1) displayed the strongest association,
followed by the haplotype containing KLK8 to KLK11
genes (block 8) (Figure C and
Supplemental Table 3). However, when only the patients with nephritis were considered as
cases, the association to the KLK3 promoter region was weaker; instead,
the haplotypes harboring the KLK4 promoter, KLK5, and
the KLK8–KLK11 block yielded stronger
association. Of note, the 2 strongest associations observed corresponded to haplotypes
with significantly lower frequencies in patients, reflecting a protective effect
conferred by the KLK polymorphism.
Overall, these 3 sets of independent findings support a likely genetic association of
the human KLK1 gene and the KLK3 promoter region with
lupus and lupus nephritis.