Analysis of the B cell receptor repertoire in six immune-mediated diseases.

B cells are important in the pathogenesis of many, and perhaps all, immune-mediated diseases. Each B cell expresses a single B cell receptor (BCR)1, and the diverse range of BCRs expressed by the total B cell population of an individual is termed the 'BCR repertoire'. Our understanding of the BCR repertoire in the context of immune-mediated diseases is incomplete, and defining this could provide new insights into pathogenesis and therapy. Here, we compared the BCR repertoire in systemic lupus erythematosus, anti-neutrophil cytoplasmic antibody (ANCA)-associated vasculitis, Crohn's disease, Behçet's disease, eosinophilic granulomatosis with polyangiitis, and immunoglobulin A (IgA) vasculitis by analysing BCR clonality, use of immunoglobulin heavy-chain variable region (IGHV) genes and-in particular-isotype use. An increase in clonality in systemic lupus erythematosus and Crohn's disease that was dominated by the IgA isotype, together with skewed use of the IGHV genes in these and other diseases, suggested a microbial contribution to pathogenesis. Different immunosuppressive treatments had specific and distinct effects on the repertoire; B cells that persisted after treatment with rituximab were predominately isotype-switched and clonally expanded, whereas the inverse was true for B cells that persisted after treatment with mycophenolate mofetil. Our comparative analysis of the BCR repertoire in immune-mediated disease reveals a complex B cell architecture, providing a platform for understanding pathological mechanisms and designing treatment strategies.

Expression profiling of Sudanese visceral leishmaniasis patients pre- and post-treatment with sodium stibogluconate.

Visceral leishmaniasis (VL) in Sudan caused by Leishmania donovani is fatal in susceptible individuals if untreated. Treatment with sodium stibogluconate (SSG) leads to post-kala-azar dermal leishmaniasis (PKDL) in 58% of patients. Here, Affymetrix microarrays were used to identify genes differentially expressed in lymph nodes (N=9 paired samples) pre- and post-treatment with SSG. Using the Bioconductor package limma, 438 genes from 28 869 post-quality-control probe sets were differentially expressed (Pnominal ≤.02) post- vs pretreatment. Canonical pathway analysis using Ingenuity Pathway Analysis™ identified "role of nuclear factor of activated T-cell in regulation of immune response" (Pnominal =1.35×10-5 ; PBH-adjusted =4.79×10-3 ), "B-cell development" (Pnominal =2.04×10-4 ; PBH-adjusted =.024), "Fcγ receptor-mediated phagocytosis in macrophages and monocytes" (Pnominal =2.04×10-4 ; PBH-adjusted =.024) and "OX40 signalling" (Pnominal =2.82×10-4 ; PBH-adjusted =.025) as pathways differentially regulated post- vs pretreatment. Major network hub genes included TP53, FN1, MYC, BCL2, JUN, SYK, RUNX2, MMP1 and ACTA2. Top endogenous upstream regulators included IL-7 (P=2.28×10-6 ), TNF (P=4.26×10-6 ), Amyloid Precursor Protein (P=4.23×10-5 ) and SPI1/PI.1 (P=1.17×10-7 ). Top predicted chemical drug regulators included the flavonoid genistein (P=4.56×10-7 ) and the quinoline alkaloid camptothecin (P=5.14×10-5 ). These results contribute to our understanding of immunopathology associated with VL and response to SSG treatment. Further replication could identify novel therapeutic strategies that improve on SSG treatment and reduce the likelihood of progression to PKDL.

Global amplification of mRNA by template-switching PCR: linearity and application to microarray analysis.

Conventional approaches to target labelling for expression microarray analysis typically require relatively large amounts of total RNA, a serious limitation when the sample available is small. Here we explore the cycle-dependent amplification characteristics of Template-Switching PCR and validate its use for microarray target labelling. TS-PCR identifies up to 80% of the differentially expressed genes identified by direct labelling using 30-fold less input RNA for the amplification, with the equivalent of 1000-fold less starting material being used for each hybridisation. Moreover, the sensitivity of microarray experiments is increased considerably, allowing the identification of differentially expressed transcripts below the level of detection using targets prepared by direct labelling. We have also validated the fidelity of amplification and show that the amplified material faithfully represents the starting mRNA population. This method outperforms conventional labelling strategies, not only in terms of sensitivity and the identification of differentially expressed genes, but it is also faster and less labour intensive than other amplification protocols.

Mapping by genetic interaction: high-resolution congenic mapping of the type 1 diabetes loci Idd10 and Idd18 in the NOD mouse.

As many of the linked chromosome regions that predispose to type 1 diabetes in the NOD mouse have been dissected, it has become apparent that the initially observed effect is in fact attributable to several loci. One such cluster of loci on distal chromosome 3, originally described as Idd10, is now known to comprise three separate loci, Idd10, Idd17, and Idd18. Although these loci have a significant combined effect on diabetes development, their individual effects are barely detectable when diabetes is used as a read-out, which makes fine-mapping them by use of a conventional congenic approach impractical. In this study, we demonstrate that it is possible to map loci, with modest effects, to regions small enough for systematic gene identification by capitalizing on the fact that the combined loci provide more profound, measurable protection. We have mapped the Idd10 and Idd18 loci to 1.3- and 2.0-cM intervals, respectively, by holding the Idd3 allele constant. In addition, we have excluded Csf1 and Nras as candidates for both loci.

Statistical modeling of interlocus interactions in a complex disease: rejection of the multiplicative model of epistasis in type 1 diabetes.

In general, common diseases do not follow a Mendelian inheritance pattern. To identify disease mechanisms and etiology, their genetic dissection may be assisted by evaluation of linkage in mouse models of human disease. Statistical modeling of multiple-locus linkage data from the nonobese diabetic (NOD) mouse model of type 1 diabetes has previously provided evidence for epistasis between alleles of several Idd (insulin-dependent diabetes) loci. The construction of NOD congenic strains containing selected segments of the diabetes-resistant strain genome allows analysis of the joint effects of alleles of different loci in isolation, without the complication of other segregating Idd loci. In this article, we analyze data from congenic strains carrying two chromosome intervals (a double congenic strain) for two pairs of loci: Idd3 and Idd10 and Idd3 and Idd5. The joint action of both pairs is consistent with models of additivity on either the log odds of the penetrance, or the liability scale, rather than with the previously proposed multiplicative model of epistasis. For Idd3 and Idd5 we would also not reject a model of additivity on the penetrance scale, which might indicate a disease model mediated by more than one pathway leading to beta-cell destruction and development of diabetes. However, there has been confusion between different definitions of interaction or epistasis as used in the biological, statistical, epidemiological, and quantitative and human genetics fields. The degree to which statistical analyses can elucidate underlying biologic mechanisms may be limited and may require prior knowledge of the underlying etiology.

NOD Idd5 locus controls insulitis and diabetes and overlaps the orthologous CTLA4/IDDM12 and NRAMP1 loci in humans.

A genome scan for B10-derived loci that reduce the frequency of diabetes and insulitis in NOD mice demonstrated a large region (34 cM) of linkage on the proximal end of chromosome 1. This locus was designated Idd5 and encompassed candidate genes including Il1r1, Il1r2, Stat1, Stat4, Nramp1, and Bcl2. In the current study, we have confirmed the existence of Idd5 by developing a series of congenic mouse strains that are resistant to diabetes and determined that Idd5 is actually two genes located within a 9.4-cM interval. Idd5.1 is in the proximal 1.5-cM portion of the interval and contains the candidates Casp8, Cflar (FLIP), Cd28, and Cd152 (CTLA4). Idd5.1 overlaps the orthologous CTLA4/IDDM12 locus in humans. Idd5.2 is in the distal 5.1-cM portion of the 9.4-cM interval and contains the candidates Nramp1, which has a functional polymorphism between NOD and B10, and Cmkar2 (CXCR2, interleukin [IL]-8 receptor alpha). Candidate genes eliminated by this analysis include Il1r1, Ilr2, Zap70, Orch5, Stat1, Stat4, Bcl2, Cmkar4 (CXCR4), and Il10. On its own, the Idd5 locus provides a significant amount of protection from diabetes (50% reduction from parental frequency) and when combined with another resistance locus (Idd3 on chromosome 3), provides nearly complete protection from diabetes and insulitis.

Analysis of the mouse CD30 gene: a candidate for the NOD mouse type 1 diabetes locus Idd9.2.

Members of the tumor necrosis factor receptor superfamily play an important role in the initiation, expansion, and termination of an immune response. It has recently been demonstrated that one member of this family, CD30, plays a central role in maintaining peripheral tolerance by controlling the expansion of autoreactive CD8+ T-cells. In the present study, Cd30 was mapped to a 5.6-cM interval on chromosome 4 containing the type 1 diabetes susceptibility locus Idd9.2. We determined the intron/exon structure of Cd30 and sequenced the exons, as well as 1.8 kb of the 5' putative promoter region, from 6 different mouse strains. Remarkably, 63 sequence variants, both coding and noncoding, were found. A total of 27 sequence variants, 4 of which were nonsynonymous, were found between the diabetes susceptible NOD strain and the resistant B10 strain. Of these sequence variants, 19 are within the promoter region. However, no difference between NOD and the congenic strain NOD.B10 Idd9R1, which has the B10 allele of Cd30, was observed in CD30 expression at either the mRNA or protein level. Given its role in protecting against autoimmunity, one or more of the coding variants within CD30 is a good candidate for the Idd9.2 etiological variant.

The NOD Idd9 genetic interval influences the pathogenicity of insulitis and contains molecular variants of Cd30, Tnfr2, and Cd137.

Previous analyses of NOD mice have shown that some genes control the development of both insulitis and diabetes, while other loci influence diabetes without reducing insulitis. Evidence for the existence of a gene only influencing diabetes, Idd9 on mouse chromosome 4, is provided here by the development of a novel congenic mouse strain, NOD.B10 Idd9. NOD.B10 Idd9 mice display profound resistance to diabetes even though nearly all develop insulitis. Subcongenic analysis has demonstrated that alleles of at least three B10 genes, Idd9.1, Idd9.2, and Idd9.3 are required to produce Idd9-mediated diabetes resistance. Candidate genes with amino acid differences between the NOD and B10 strains have been localized to the 5.6 cM Idd9.2 interval (Tnfr2, Cd30) and to the 2.0 cM Idd9.3 interval (Cd137).

Differential glycosylation of interleukin 2, the molecular basis for the NOD Idd3 type 1 diabetes gene?

The insulin-dependent diabetes (Idd) gene, Idd3, has been localised to a 0.35 cM region of chromosome 3 containing the structural gene for the cytokine interleukin 2 (IL-2). While variation of the N-terminal amino acid sequence of IL-2 has been shown to correlate with Idd3 allelic variation, differences in induction of proliferation by IL-2 allotypes have not been detected. In the current study, we examined the electrophoretic migration of IL-2 allotypes and have found two distinct patterns, consistent with differences in glycosylation, that correlate with diabetes-resistance and susceptibility. These findings strongly suggest that IL-2 variants may be functionally distinct.

Congenic mapping of the type 1 diabetes locus, Idd3, to a 780-kb region of mouse chromosome 3: identification of a candidate segment of ancestral DNA by haplotype mapping.

Type 1 diabetes in the nonobese diabetic (NOD) mouse arises as a consequence of T cell-mediated destruction of the insulin-producing beta cells of the pancreas. Although little is known of the events that initiate and subsequently drive beta-cell destruction it is clear that the entire process is under complex genetic control. At present 19 loci have been mapped that influence the development of diabetes either at the level of initiation of insulitis or at the level of progression from insulitis to overt diabetes, or both. Previously, we have mapped one of these loci, Idd3, to a 0.35-cM interval on proximal mouse chromosome 3. In the present study we have narrowed the map position of this locus to an interval of 0.15 cM by a combination of novel congenic strains and an ancestral haplotype analysis approach. We have constructed a physical contig in bacterial artificial chromosome (BAC) clones across the minimal interval. Restriction mapping of the BAC contig placed the maximum size of the Idd3 interval at 780 kb between the markers D3Nds36 and D3Nds76. To refine further the Idd3 interval we developed a series of novel single nucleotide polymorphisms (SNPs) and carried out haplotype analysis on DNA from mouse strains known to carry either Idd3 susceptibility or protective alleles. This haplotype analysis identified a 145-kb segment of ancestral DNA between the microsatellite marker D3Nds6 and the SNP 81.3. One haplotype of this ancestral segment of DNA is found in mouse strains carrying an Idd3 susceptibility allele and another is found in mouse strains carrying an Idd3 protective allelle. Within the 780-kb congenically defined interval this 145-kb segment represents the most likely location for Idd3. The Il2 gene, which encodes the cytokine interleukin 2 (IL2), maps to this interval and is a strong candidate for Idd3. To investigate whether sequence variation exists in the promoter region of the Il2 gene, which might alter its expression, we sequenced the promoter region of the Il2 gene from mouse strains carrying either an Idd3 susceptibility or resistance allele. Two sequence variants were identified, neither of which fell in known regulatory elements within the Il2 promoter. In agreement with this observation steady-state Il2 mRNA levels showed no variation between susceptible and resistant mouse strains. These data suggest that the profound protection from diabetes seen in congenic mice carrying an Idd3 protective allele is unlikely to be due to differences in the level of expression of the Il2 gene. Instead, all of the current data support our hypothesis that Idd3 corresponds to amino acid variation at the amino terminus of Il2.

Mouse chromosome 3.

> http://link.springer-ny. com/link/service/journals/00335/bibs/10n10p942.html

Localization of two insulin-dependent diabetes (Idd) genes to the Idd10 region on mouse chromosome 3.

Multiple genes control the development of autoimmune diabetes both in humans and in the nonobese diabetic (NOD) strain of mouse. Previously, three insulin-dependent diabetes (Idd) genes, Idd3, Idd10, and Idd17, were localized to mouse Chromosome (Chr) 3. The B10- or B6-derived resistance alleles at Idd10 and Idd3 together provide the NOD mouse with nearly complete protection from diabetes. In the present study, the 10.2-cM region encoding Idd10 was defined further with newly developed congenic strains. A locus, located in the centromeric 2.1 cM of the 10.2 cM region, contributed to the Idd10 trait. However, this locus did not account for the full effect of Idd10, suggesting the presence of a second gene in the distal portion of the 10.2-cM region. This second gene is designated as Idd18 and is localized to a 5.1-cM region. The resolution of the originally defined Idd3 locus into at least four separate loci, Idd3, Idd10, Idd17, and Idd18, illustrates the complex polygenic nature of diabetes.

Congenic mapping of the insulin-dependent diabetes (Idd) gene, Idd10, localizes two genes mediating the Idd10 effect and eliminates the candidate Fcgr1.

The development of autoimmune diabetes in the nonobese diabetic (NOD) mouse is under the control of multiple insulin-dependent diabetes (Idd) genes. The Idd3 gene, originally defined as a broad peak of linkage on mouse chromosome 3, was subsequently identified as two genes, Idd3 and Idd10, separated by at least 20 cM. The resistance alleles of Idd3 and Idd10 individually confer only partial protection from diabetes but, in combination, result in profound resistance to disease due to an epistatic genetic interaction. In this study, we used newly developed congenic strains to further localize Idd10. Surprisingly, we found that Idd10 itself comprises at least two linked loci: Idd10 and the newly designated Idd17. Idd17 was localized to a 1.1-cM region between D3Mit26 and D3Mit40, proximal to Fcgr1, a candidate gene encoding the high affinity Fc receptor for IgG. Idd10 was localized to a 10-cM region between D3Mit213 and D3Mit106, distal to Fcgr1. Thus, Fcgr1 was excluded as a candidate for either Idd10 or Idd17, despite the fact that the NOD strain expresses a mutant form of the receptor. Interestingly, although Idd10 and Idd17 participate in a genetic interaction with each other, Idd10 but not Idd17 participates in the genetic interaction with Idd3. Our study on chromosome 3 begins to reveal the extent of the polygenic nature of autoimmune diabetes, and demonstrates that the use of congenic strains is an effective mapping strategy, even in the dissection of multiple, linked genes with subtle effects.

A 2-Mb YAC contig and physical map of the natural killer gene complex on mouse chromosome 6.

We have constructed a physical map of a > 2-Mb region on mouse chromosome 6 that contains the natural killer gene complex (NKC). The map comprises a contig of 14 overlapping yeast artificial chromosomes onto which we positioned 25 NKC markers. NKC genetically linked genes encode > 17 proteins that directly control innate NK cell-mediated tumor lysis and disease resistance. Herein we show that Nkrp1 genes are clustered in a region flanked by A2m and Cd69 genes and that most Ly49 genes are clustered in a distal region -1 Mb distant. Importantly, syntenic intervals of mouse chromosome 6 and human chromosome 12p that include the NKC are conserved. NKC species conservation suggests that the human NKC may contain orthologues for the mouse viral disease resistance genes, Cmv1 and Rmp1. The high-resolution NKC map will facilitate investigation of NKC gene regulation and identification of phenotypically defined gene products that confer NK cell defense against viral pathogens.

Mapping of the IDDM locus Idd3 to a 0.35-cM interval containing the interleukin-2 gene.

Currently, 16 loci that contribute to the development of IDDM in the NOD mouse have been mapped by linkage analysis. To fine map these loci, we used congenic mapping. Using this approach, we localized the Idd3 locus to a 0.35-cM interval on chromosome 3 containing the Il2 gene. Segregation analysis of the known variations within this interval indicated that only one variant, a serine-to-proline substitution at position 6 of the mature interleukin-2 (IL-2) protein, consistently segregates with IDDM in crosses between NOD and a series of nondiabetic mouse strains. These data, taken together with the immunomodulatory role of IL-2, provide circumstantial evidence in support of the hypothesis that Idd3 is an allelic variation of the Il2 gene, or a variant in strong linkage disequilibrium.

Effect of natural sequence variation at the H-2Ld-restricted CD8+ T cell epitope of the murine cytomegalovirus ie1-encoded pp89 on T cell recognition.

The amino acid sequence YPHFMPTNL of pp89, the ie1-encoded product of murine cytomegalovirus (MCMV; Smith strain), constitutes an immunodominant T cell epitope recognized in association with H-2Ld. Nucleotide sequencing of MCMV isolates derived from wild mice identified variation between amino acids 147-192 of pp89 in 19 of 27 isolates, including the region encompassing the CTL epitope (amino acid residues 168-176). Four groups of isolates with naturally occurring variant sequences for the CTL epitope were defined: (1) YPHFMPPNL; (2) YPHFMPPSL; (3) YPHFIPPSL; and (4) YLDFMPPNL. The remaining isolates, and the laboratory strains K181 and Vancouver, showed complete identity with the Smith strain. Polyclonal pp89 (Smith strain)-specific CTL only weakly recognized target cells infected with MCMV from most variant groups. No lysis of cells infected with isolate N1 from group 4 was detected. Analyses of cross-reactive recognition of YPHFMPTNL peptide-coated targets by CTL primed with variant MCMV isolates showed that the group 2 and 3 isolates, G4 and K6, respectively, but not the group 4 isolate N1, elicited CTL that exhibited a cross-reactive response. Furthermore, while the group 2 and 3 isolates G4 and K6 were able to prime CTL responses that displayed reactivity to homologous pp89 variant nonapeptides, the group 4 isolate N1 failed to do so. Finally, while immunization of mice with the nonapeptide YPHFMPTNL conferred significant protection against the laboratory strain K181 [correction of Kl81], no evidence of protection was observed for the group 2 and 4 variants G4 and N1, respectively. These observations raise the possibility that clinical isolates of HCMV may also differ in sequence from potential vaccine strains at immunodominant epitopes for CD8+ T cells thus reducing the efficacy of vaccination.

Assessment of antigenicity and genetic variation of glycoprotein B of murine cytomegalovirus.

An analysis of linear antibody-binding sites of the glycoprotein B (gB) molecule of murine cytomegalovirus (MCMV) and of genetic variation within these regions was performed. To achieve this, a series of overlapping fragments spanning the entire coding sequence of the gB gene of the K181 strain of MCMV was expressed in E. coli as fusion proteins with glutathione S-transferase (GST) using the pGEX expression system. Four antibody-binding regions were mapped to locations spanning amino acid residues 17-79 (BS), 155-278 (BE2), 809-926 (SS) and 347-508 (BB and EE), based on reactivity in Western blot analysis of GST-gB fusion proteins with murine polyclonal antiserum raised against MCMV. Only the antibody-binding region BE2 (155-278) elicited an antiserum that exhibited complement-dependent neutralizing activity, and immunization of mice with the fusion protein BE2 led to moderate but significant reductions in the level of MCMV replication in the spleen. Polyclonal antisera raised against the GST-gB fusion proteins detected purified virion proteins of 105 kDa (anti-BS and anti-BE2) and 52 kDa (anti-SS) and are therefore likely to recognize the N-terminal and C-terminal portions of the gB molecule, respectively. The antibody-binding region within amino acid residues 17-79 was found to be MCMV strain-specific, whereas antibody-binding regions within residues 155-278 and 809-926 were found to be conserved among MCMV field isolates. Comparative sequence analysis of the corresponding regions of MCMV gB revealed a level and extent of sequence of sequence heterogeneity consistent with these findings.