Advertisement
Research Article Free access | 10.1172/JCI29943
1Department of Immunology, Weizmann Institute of Science, Rehovot, Israel. 2Pediatric Department, Dana Children’s Hospital, Tel-Aviv Sourasky Medical Center, Tel-Aviv, Israel. 3Sackler Faculty of Medicine, Tel-Aviv University, Tel-Aviv, Israel.
Address correspondence to: Irun R. Cohen, Department of Immunology, Weizmann Institute of Science, Rehovot 76100, Israel. Phone: 972-8-934-2911; Fax: 972-8-934-4103; E-mail: irun.cohen@weizmann.ac.il.
Find articles by Merbl, Y. in: JCI | PubMed | Google Scholar
1Department of Immunology, Weizmann Institute of Science, Rehovot, Israel. 2Pediatric Department, Dana Children’s Hospital, Tel-Aviv Sourasky Medical Center, Tel-Aviv, Israel. 3Sackler Faculty of Medicine, Tel-Aviv University, Tel-Aviv, Israel.
Address correspondence to: Irun R. Cohen, Department of Immunology, Weizmann Institute of Science, Rehovot 76100, Israel. Phone: 972-8-934-2911; Fax: 972-8-934-4103; E-mail: irun.cohen@weizmann.ac.il.
Find articles by Zucker-Toledano, M. in: JCI | PubMed | Google Scholar
1Department of Immunology, Weizmann Institute of Science, Rehovot, Israel. 2Pediatric Department, Dana Children’s Hospital, Tel-Aviv Sourasky Medical Center, Tel-Aviv, Israel. 3Sackler Faculty of Medicine, Tel-Aviv University, Tel-Aviv, Israel.
Address correspondence to: Irun R. Cohen, Department of Immunology, Weizmann Institute of Science, Rehovot 76100, Israel. Phone: 972-8-934-2911; Fax: 972-8-934-4103; E-mail: irun.cohen@weizmann.ac.il.
Find articles by Quintana, F. in: JCI | PubMed | Google Scholar
1Department of Immunology, Weizmann Institute of Science, Rehovot, Israel. 2Pediatric Department, Dana Children’s Hospital, Tel-Aviv Sourasky Medical Center, Tel-Aviv, Israel. 3Sackler Faculty of Medicine, Tel-Aviv University, Tel-Aviv, Israel.
Address correspondence to: Irun R. Cohen, Department of Immunology, Weizmann Institute of Science, Rehovot 76100, Israel. Phone: 972-8-934-2911; Fax: 972-8-934-4103; E-mail: irun.cohen@weizmann.ac.il.
Find articles by Cohen, I. in: JCI | PubMed | Google Scholar
Published March 1, 2007 - More info
Natural antibodies are autoreactive/polyreactive antibodies believed to be secreted in the absence of xenoantigens. The origin and functional role of this limited and selective autoimmunity are not clear, nor is the specificity and range of autoantigens that drive the development of B cells producing natural antibodies. In this issue of the JCI, Merbl et al. report that in utero, humans generate natural IgM and IgA antibodies that recognize a uniform set of autoantigens (see the related article beginning on page 712), some of which are associated with autoimmune diseases. The authors postulate that this “autoimmunity” at birth favors the emergence of autoimmune diseases in later life. We present a molecular basis for the limited and common repertoire of antibodies produced by fetal B cells, which may be distinct from the abnormalities in B cell development described in patients with autoimmune diseases.
Eric Meffre, Jane E. Salmon
Autoimmune diseases are often marked by autoantibodies binding to self antigens. However, many healthy persons also manifest autoantibodies that bind to self antigens, known as natural autoantibodies. In order to characterize natural autoantibodies present at birth, we used an antigen microarray (antigen chip) to analyze informatically (with clustering algorithms and correlation mapping) the natural IgM, IgA, and IgG autoantibody repertoires present in 10 pairs of sera from healthy mothers and the cords of their newborn babies. These autoantibodies were found to bind to 305 different, mostly self, molecules. We report that in utero, humans develop IgM and IgA autoantibodies to relatively uniform sets of self molecules. The global patterns of maternal IgM autoantibodies significantly diverged from those at birth, although certain reactivities remained common to both maternal and cord samples. Because maternal IgG antibodies (unlike IgM and IgA) cross the placenta, maternal and cord IgG autoantibodies showed essentially identical reactivities. We found that some self antigens that bind cord autoantibodies were among the target self antigens associated with autoimmune diseases later in life. Thus, the obviously benign autoimmunity prevalent at birth may provide the basis for the emergence of some autoimmune diseases relatively prevalent later in life.
Natural antibodies are antibodies detected in the absence of known immunization (1, 2). Although autoimmunity is forbidden by the clonal selection theory (3), many natural antibodies are autoantibodies; they bind to self molecules. The functions of natural autoantibodies are not clear, but the specific self molecules recognized by these autoantibodies appear to form clinically defining signatures: some autoantibodies create a pattern that heralds susceptibility to a future autoimmune disease, while a different autoantibody pattern can mark resistance to the disease (4). Indeed, it has been proposed that natural autoantibodies and autoreactive T cells in healthy individuals may be directed to a specific and limited set of self molecules; this selective autoimmunity has been termed the immunological homunculus (5–7) or the immunculus (8) — the immune system’s internal representation of the body.
In order to characterize natural autoantibodies present at birth, their isotypes, and the self molecules they bind, we used an antigen microarray device to analyze informatically, with clustering algorithms and correlation mapping, the natural IgM, IgA, and IgG autoantibody repertoires present in 10 pairs of sera from healthy mothers and their newborn babies. These autoantibodies were found to bind to 305 different molecules, most of them self molecules. Because only maternal IgG antibodies, but not IgM or IgA antibodies, cross the placenta to the fetus (9–11), the IgM and IgA autoantibodies present in cord serum at birth would have had to arise as a consequence of prenatal immune activation in the baby. Thus, cord IgM and IgA antibodies originate within the developing baby; cord IgG comes from the mother. We now report that different babies manifested cord IgM autoantibodies binding to a highly correlated, relatively uniform set of self molecules and that cord and maternal IgM reactivities clustered separately. Thus, natural autoimmunity begins in utero in healthy humans, and the prenatal autoantibody repertoires are directed to a specific, standardized set of body molecules, the immunological homunculus (5). Many cord autoantibodies bound self molecules that are associated with major autoimmune diseases later in life. These findings relate to our understanding of both natural autoimmunity and autoimmune disease.
Analysis by microarray showed that certain autoantibodies were indeed quite prevalent at birth. Table 1 lists the self molecules bound by IgM, IgG, or IgA autoantibodies present in 8 or more of the 10 individual cord sera. Antibody binding to a molecule was scored as positive when the mean intensity of the laser signal was at least 2 standard deviations above the mean background control. Some of the biologic associations of these self antigens are also shown. Table 2 lists the self molecules bound by maternal and cord autoantibodies according to their mean reactivity index (MRI), rather than by their prevalence. The MRI denotes the fold increase (rounded off) above the control of the mean reactivity to the self antigen; the MRI values shown in Table 2 are at least 2-fold greater than the mean value from PBS incubation alone (see Methods) plus 2 standard deviations in the same sample. For reference, note that the MRI of natural antibodies to E. coli LPS in maternal sera manifested an MRI of 2 for IgG and IgA and 5 for IgM (Table 2); thus, some autoantibodies manifested a much higher degree of binding to self molecules than did antibodies to immunogenic foreign LPS.
In order to test whether the natural antibodies of different people are indeed directed to a common set of antigens, we used Pearson correlation coefficients to generate a correlation map of the global sets of positive IgM, IgA, and IgG reactivities for each of the 10 cord and maternal repertoires compared with the other 9 samples within each group (Figure 1). Correlations, which mathematically depict the degrees of relationship between individual whole repertoires, are represented by intensity of shading from white (correlation coefficient of 1, perfect correlation) to black (correlation coefficient of 0, no correlation). The range of lower and upper bounds for a 95% confidence interval for each group of samples are shown in Table 3. Figure 1 shows that individual maternal IgM, IgA, and IgG repertoires differed from one another; the white diagonal of self identity contrasts with the dark gray background of individual repertoire variability for each antibody isotype. The cords show a different picture: each cord differed from the others in IgG and IgA repertoires, but the cord IgM repertoires were remarkably uniform. The correlations between individual cord IgM repertoires were so high — so close to white — that the white diagonal of identity marking the perfect correlation of each serum’s global reactivity with itself is barely discernible. In other words, the global repertoire of IgM antibodies expressed in any cord is quite similar to that expressed by each of the other 9 cords. We conclude that different babies produce IgM antibodies in utero to a very similar set of self molecules (12, 13).
Correlation coefficient map comparing each sample with each of the other samples. The 10 cord and 10 maternal samples are shown, and the mean IgG, IgM, and IgA reactivities of the 305 tested antigens were compared to determine correlation. The scale shows the degree of correlation from 0 (black) to 1 (white).
Ranges of lower and upper bounds for a 95% confidence interval for each correlation coefficient
We used clustering algorithms to reveal relationships between defined subsets of the global data; membership in a cluster is determined automatically by an algorithm that mathematically measures shared similarities and dissimilarities. Figure 2 illustrates the results of hierarchical clustering of mother-cord sets (columns) according to each of their IgM and IgA antigen reactivities (rows); the relative degrees of antibody binding are denoted by color, ranging from dark blue to dark red. Note that the range of antibody binding intensity spans several orders of magnitude (see Methods), while the colors are relatively few; thus significantly different reactivities may appear to share the same color. Nevertheless, the representation by color does give some visual sense of the actual clustering. Also note that the positions from top to bottom of particular antigens in the clustering process (Figure 2, rows) were arranged automatically by the clustering algorithm according to relationships between the antigen reactivities; antigens that tend to manifest correlated antibody responses were automatically arranged in adjacent rows. Thus the order of the rows of antigens differs in Figure 2A and Figure 2B because the IgM and IgA repertories show correlated antibody reactivities to different sets of antigens. The clustering of maternal and cord samples shown in Figure 2 demonstrates that the IgM repertoires of the maternal and cord samples differ: the 10 cord IgM repertoires clustered separately from the 10 maternal IgM repertoires (Figure 2A). Thus, the global IgM repertoire of each cord is more similar to the global IgM repertoire of each of the other 9 cords than it is to any of the maternal global repertoires. Note that the maternal samples clustered together separately from the cord samples because the cord samples share a commonality that sets them apart from the maternal samples; the maternal repertoires do differ individually from one another (see Figure 1). In contrast, the IgA repertoires of cord and maternal samples could not be separated into clusters by the unsupervised clustering algorithm applied to the global reactivities to all 305 antigens (Figure 2B). However, we detected differences between the maternal and cord samples as separate groups by selecting a more limited set of 24 of the 305 arrayed antigens, all of which differed significantly between the maternal and cord samples by Wilcoxon rank-sum test (Table 4). Some of the antibodies manifested higher antigen reactivities in the maternal sera and some in the cords; the response to these antigens was not uniformly high or low between the groups. The highest ranking antigens that separated the IgM repertoires (which were significantly different between cord and maternal samples by Wilcoxon rank-sum test) are listed in Table 5. The IgG antibody repertoires of each mother and her newborn separated into essentially identical pairs upon hierarchical clustering (data not shown); this can be explained by the active transport of maternal IgG across the placenta (9–11).
Clustering of cord and maternal IgM and IgA. Hierarchical clusterings of IgM (A) and IgA (B) reactivities for the maternal and cord samples based on their reactivities to each of the 305 antigens in the array. A color code denoting maternal (blue) and cord (red) samples is shown at the bottom. The color scale shows the relative degree of antibody binding from low (dark blue) to high (dark red).
Obviously, the several hundred antigens used in the present study are not an exhaustive list of natural autoantibody self reactivity. Nevertheless, the sample of self molecules, however limited, is illuminating. We concluded that the cord IgM autoantibody repertoire was both selective (only some self molecules are recognized) and uniform (different individuals express a common repertoire). Moreover, the common cord IgM repertoire clustered separately from the maternal IgM repertoires (Figure 2A). Cord IgA was less uniform than was cord IgM, and cord IgA did not cluster separately from maternal IgA; cord and maternal samples, however, differed as groups in their IgA reactivities to a set of 24 antigens (Table 4). The organization of self antigens into correlated sets will be reported in a future study.
Mouthon and associates, using a blot technique to assay natural antibodies binding to undefined antigens in tissue extracts, could not detect the differences in IgM autoantibody repertoires between neonate and adult (demonstrated herein using the antigen chip), nor could their method name the antigens (12, 13). Mirilas and coworkers (14) could detect the partial evolution of the natural autoimmune repertoire based on testing a relatively limited number of self antigens.
The present antigen chip findings are compatible with the immunological homunculus hypothesis (5–8) — the idea that the healthy immune system includes autoreactivity to a selected set of particular self molecules. Indeed, the development of our antigen microarray chip (4) and the present work were undertaken to test the homunculus hypothesis. The fact that the primary autoantibody homunculus arises during the uterine life of the fetus, an environment normally free of foreign antigens, suggests that self molecules are likely to serve as the immunogenic stimulus inducing the IgM and IgA homuncular autoantibodies. How this happens is not yet known, but it seems that the development of B cells, like that of T cells, may involve positive selection for self reactivity (15).
The second disclosure of the present work is that some of the prevalent natural autoantibodies bind self antigens associated with prevalent autoimmune diseases. Tables 1 and 2 reveal that humans are born with autoantibodies to clinically important self antigens, including single-stranded and double-stranded DNA (16), targeted in SLE; glutamic acid decarboxylase (GAD), targeted in type 1 diabetes (17); myelin oligodendrocyte glycoprotein (MOG), targeted in multiple sclerosis (18); cytokeratin, targeted in rheumatoid arthritis (19); and thyroglobulin, targeted in thyroiditis (20). In patients with clinical autoimmune diseases, however, the amounts of the autoantibodies to these self molecules are usually, but not always, greater than in healthy individuals; moreover, the autoantibodies associated with autoimmune diseases tend to be predominantly of the IgG isotype, rather than of the IgM isotype.
Prevalent natural autoantibodies also bind to molecules not known to be associated with autoimmune diseases, such as ubiquitin, galectins, gelsolin, and annexin, which perform critical cellular functions. These molecules participate in cytoskeleton regulation, signal transduction, and proteosome activity, but they also affect immune regulation. Galectins, for example, are immune modulators (21), and gelsolin has been shown to bind the powerful immune activator endotoxin (22); the effects of natural autoantibodies on the functions of these molecules are currently unknown. Hsps also regulate immune reactivity, and T cell autoreactivity to Hsp molecules appears to downregulate inflammation and autoimmune diabetes and arthritis (23, 24). Some antibodies to LDL have been reported to protect against atherosclerosis (25). Thus, autoimmunity to some self molecules might actually serve to prevent autoimmune diseases. Note that much of this self reactivity is of the IgM isotope; it has been proposed that autoimmune diseases might be blocked by natural autoantibodies of the IgM isotype that mask self antigens from pathogenic autoimmune T cells (26).
The self molecule ubiquitin is essential to the destruction and processing of cell molecules; recently it has been reported that a domain of the autoimmune regulator (AIRE) molecule functions as a ligase for ubiquitin. The AIRE gene plays a central role in shaping the T cell repertoire (27); mutations in AIRE can lead to autoimmune disease (28). Table 2 shows that cord blood, but not maternal blood, contained a relatively high amount of IgM anti-ubiquitin. It remains to be seen how natural IgM autoantibodies to ubiquitin might influence the immune functions of AIRE in newborns. Other molecules such as plasmin and factor II were targeted by cord IgM, but not by maternal IgM; similarly, the maternal IgM autoimmune repertoire demonstrated reactivities absent from cord serum. Thus different natural autoantibodies appear to mark different stages of human growth and development. If we assume that the mothers began life with IgM autoantibody repertoires similar to those of the cords (Figures 1 and 2), then it would appear that the individuality of one’s autoimmune repertoire evolves through one’s postnatal immune experience.
The prevalence of cord autoantibodies to enzymes and plasma proteins (Tables 1 and 2) also is intriguing; β2-microglobulin, for example, takes part in antigen presentation to T cells (29). Some natural autoantibodies may have antitumor effects: tyrosinase, involved in the synthesis of melanin, is a melanoma-associated antigen (30), and carbonic anhydrase IX is a tumor-associated antigen in renal cell carcinoma (31).
Previously, monoclonal autoantibodies binding to 3 or more different self antigens, termed polyspecific autoantibodies, have been isolated from human fetuses and cord blood (32). These polyspecific autoantibodies were detected in ELISA assays by their binding to actin, tubulin, myosin, thyroglobulin, and/or transferrin and accounted for about 13%–16% of the clones that were isolated. Our microarray study did not test monoclonal antibodies, and so we could not determine the degree of polyspecificity of any of the autoantibodies we detected. However, autoantibodies binding to actin, tubulin or transferrin were not prevalent in our study. Thus the prevalent autoantibodies we detected with our antigen chip technology behaved differently from the classical polyspecific autoantibodies detected by other means.
Because the adaptive immune system has evolved to produce autoantibodies to certain self molecules even before birth, as we have shown here, it is reasonable to conclude that these autoantibodies must provide advantages that offset the occasional autoimmune disease associated with their target autoreactivities (5, 33). It is conceivable that natural autoantibodies, in recognizing specific body molecules, help the immune system gather and integrate essential information about the state of the body’s cells (34). In any case, the inclusion of some major disease–associated self antigens in the congenital autoantibody repertoire suggests that pathologic autoimmune disease could arise through a lapse in the regulation of natural, otherwise benign autoimmunity (5).
Serum samples. Blood samples were obtained by random availability from 10 healthy women at the onset of labor and from 10 serum samples of the cord blood of their newborns, in the course of normal procedures. All samples were collected with informed consent and approval of the Helsinki committee of the Tel-Aviv Sourasky Medical Center. The newborns were healthy at the term of pregnancy and appropriate in development and weight for gestational age. The blood samples were allowed to clot at room temperature. After centrifugation, sera were collected and stored at –20°C.
Antigens. We spotted 305 antigens on each microarray, as described previously (4). For the most part, we used the same antigens as in our previous studies of natural autoimmune repertoires (4); these included proteins, synthetic peptides from the sequences of key proteins, nucleotides, phospholipids, and other self and non-self molecules. See Supplemental Table 1 (supplemental material available online with this article; doi:10.1172/29943DS1) for the full list.
Antigen microarray. Antigen microarrays were prepared as described previously (4). Briefly, antigens were spotted in replicates of 4, and the microarrays were blocked for 1 hour at 37°C with 1% bovine serum albumin and incubated under a coverslip overnight at 4°C with a 1:5 dilution of the test serum in blocking buffer. The quantitative range of signal intensity of binding to each antigen spot was 0.01–65,000, and this range of detection made it possible to record reliable data with little dilution of test samples. The arrays were then washed and incubated for 1 hour at 37°C with a 1:500 dilution of detection antibodies. Three detection antibodies were used: a goat anti-human IgG Cy3-conjugated antibody, a goat anti-human IgM Cy5-conjugated antibody, and a goat anti-human IgA Cy5-conjugated antibody (all purchased from Jackson ImmunoResearch Laboratories Inc.). Each sample was analyzed in 2 microarray slides: one with the IgG and IgM fluorescence-labeled anti-isotypes and one with the IgG and IgA fluorescence-labeled anti-isotypes. Thus, the IgG repertoires — which showed no difference between maternal and cord samples because maternal IgG antibodies cross the placenta to the fetus (11) — served as controls for the maternal and cord IgM and IgA determinations measured simultaneously along with the IgG on each plate. Image acquisition by laser and quantification were performed as previously described (4).
Data preprocessing and background filtering. Antigen reactivity was defined by the mean intensity of the 4 replicates binding to that antigen on the microarray. Each chip was normalized by its median reactivity to account for differences in total protein concentrations that affect the background intensity level. We identified positive antibodies (Tables 1 and 2) in the following way. In order to establish the minimum level of significant antibody binding, we calculated the mean reactivity level of 32 spots incubated with PBS in place of an antigen on each microarray slide. A signal was scored as positive when it expressed intensity greater than the upper limit of the PBS control, which was defined as the mean intensity of the PBS spots plus twice the standard deviation. Intensity below this background was considered a negative result. The cutoff signal intensity level for each antibody was as follows: cord IgG, 4,830; cord IgM, 1,121; cord IgA, 863; maternal IgG, 4,333; maternal IgM, 5,944; maternal IgA, 2,204. The difference in background intensity between cord and maternal samples in IgM and IgA can be attributed to the relatively greater amounts of total IgM and IgA in maternal blood.
For global pattern analysis, such as in clustering correlation studies, we used the detected intensity of each spot without a threshold intensity cutoff; however, signal intensities that were lower than 0.01 were set to 0.01 and included in the data set.
Data analysis. The processed data set included 3 matrices of IgG, IgA, and IgM antibody reactivities, each consisting of 305 rows and 30 columns. Each column contained the reactivities measured on a given array (sample), and each row contained the reactivities measured for a given antigen over all arrays. Data were normalized by dividing each array by its median reactivity; this made it possible to transform different arrays to comparable reactivity levels.
Clustering. Raw data were normalized and analyzed using GeneSpring software (version 7; Silicon Genetics) and MATLAB (version 6.5.0.180913a, release 13; The MathWorks). Clustering analysis was performed for each isotype separately to test samples from the different groups (cord and mother). The Euclidean distance measure was used to detect related samples and related antigens, which are not discussed here. The results of clustering were robust in a leave-one-out test (35). In addition, we used a Wilcoxon rank-sum test (36) to identify antigens that separated each of the groups (see Tables 4 and 5). The Benjamini and Hochberg false discovery rate method (37) was applied, using a P value of 0.05 to determine significance.
Correlation mapping. Using the mean antibody reactivity value of each sample, we calculated the common Pearson correlation coefficient between that sample and each of the other samples within the cord sera and maternal sera groups. The correlation coefficients from 0 to 1 were expressed using a black and white color bar with gradations of gray (Figure 1). The range of lower and upper bounds for a 95% confidence interval for each of the correlation coefficients were evaluated and are shown in Table 3.
We thank Shirley Horn-Saban and Ron Ophir (Weizmann Institute of Science) for technical support. A grant from the Center for the Study of Emerging Diseases, Jerusalem, supported part of this study.
Address correspondence to: Irun R. Cohen, Department of Immunology, Weizmann Institute of Science, Rehovot 76100, Israel. Phone: 972-8-934-2911; Fax: 972-8-934-4103; E-mail: irun.cohen@weizmann.ac.il.
Nonstandard abbreviations used: AIRE, autoimmune regulator; GAD, glutamic acid decarboxylase; MOG, myelin oligodendrocyte glycoprotein; MRI, mean reactivity index.
Conflict of interest: The authors have declared that no conflict of interest exists.
Reference information: J. Clin. Invest.117:712–718 (2007). doi:10.1172/JCI29943
Yifat Merbl’s present address is: Systems Biology Department, Harvard Medical School, Boston, Massachusetts, USA.
Francisco J. Quintana’s present address is: Center for Neurologic Diseases, Harvard Medical School, Boston, Massachusetts, USA.
Yifat Merbl and Merav Zucker-Toledano contributed equally to this work.
See the related article at Autoantibody selection and production in early human life.