Incorporated by reference in its entirety herein is a computer-readable nucleotide/amino acid sequence listing submitted concurrently herewith and identified as follows: One 29,679 byte XML file named “SeqList” created on May 23, 2023.
The invention relates to autoantibodies, methods, kits and their use to identify false positive antibody markers of cancer or autoimmune condition.
Autoantibodies have been reported in individuals with autoimmune disease and cancer. They were believed to be absent in healthy individuals due to the immune tolerance mechanism; however, some have been found frequently in healthy individuals, and these antibodies are called common autoantibodies. The presence of common autoantibodies can confound the search for disease-linked autoantibodies, and their documentation will simplify the identification of autoantibodies specific to certain diseases. Indeed, only a small fraction of the autoantibodies reported in the literature have been validated in independent cohorts, suggesting the classification performance for many reported autoantibodies requires further investigation.
A comprehensive documentation of common autoantibodies will facilitate the elucidation of the complex immunology underlying their elicitation. One class of common autoantibodies is referred to as natural antibodies (NAbs). Unlike adaptive antibodies, NAbs are synthesized by B1 lymphocytes (bearing CD20+CD27+CD43+CD70−) and marginal zone B cells, and do not undergo affinity maturation by antigen stimulation or extensive somatic mutation (Coutinho et al., 1995). Another class of common autoantibodies may arise from cross-reactive antibodies to infectious agent proteins when the similarity in foreign and self-peptides may activate self-reactive T or B cells. It has been experimentally demonstrated that patients with either measles virus or herpes simplex virus type 1 produce antibodies against a viral phosphoprotein that cross-react with an intermediate filament protein of human cells. Additionally, transgenic mice infected with lymphocytic choriomeningitis virus (LCMV) may develop chronic inflammation in the central nervous system (CNS) due to epitopes shared between LCMV proteins and CNS antigens. Several bioinformatics techniques have been developed to discover potential mimicry candidates.
The immunogenicity of a protein can be attributed to its intrinsic properties and extrinsic responses by the host. Biochemical and structural properties like flexibility, hydrophilicity and beta-turns can promote antigenicity while hydrophobicity, alpha-helices and beta-sheets can suppress antigenicity. That these common autoantibodies do not cause evident autoimmune disease is intriguing. The presence of autoantibodies in serum reflects leakiness of central and/or peripheral tolerance mechanisms. However, their presence does not guarantee a causal role in autoimmune disease development. For autoantibody-induced pathology, the autoantibody needs to bind to the autoantigen to form an immune complex. Sequestration of the autoantigen from autoantibodies can inhibit the autoantibody-induced pathology.
Described herein is a method of increasing the accuracy of a cancer or an autoimmune condition screening in a subject suspected of having cancer or an autoimmune condition by identifying at least one false positive indicator. The method comprises providing a biological sample from the subject; detecting the presence of one or more autoantibodies against at least one autoantigen listed in Table 1 in the biological sample and identifying at least one false positive indicator upon detecting the presence of the one or more autoantibodies. Also described is a method of verifying the accuracy of a diagnosis of a cancer or autoimmune condition in a subject. The method comprises providing a biological sample from the subject; detecting the presence of one or more autoantibodies against at least one autoantigen listed in Table 1 in the biological sample, an identifying the subject as falsely diagnosed with cancer or an autoimmune condition upon detecting the presence of the one or more autoantibodies.
In some implementations, the at least one autoantigen is selected from the group consisting of: AMY2A, SPAG8, GTSE1, MAK, RNF138, PMFBP1, CCDC34, ODF2, CSF3, CAPN3, MYLK2, TRIM29, SOX2, and STMN4. Where the biological sample is from the subject's pancreas, the at least one autoantigen is AMY2A. Where the biological sample is from the subject's lungs, the at least one autoantigen is CSF3. Where the biological sample is from the subject's salivary gland, nerve tibial, or brain, the at least one autoantigen is SOX2. Where the biological sample is from the subject's esophagus, skin, or vagina, the at least one autoantigen is TRIM29. Where the biological sample is from the subject's skeletal muscle, the at least one autoantigen is CAPN3 and/or MYLK2. Where the biological sample is from the subject's testis, the at least one autoantigen is selected from the group consisting of: SPAG8, GTSE1, MAK, RNF138, PMFBP1, CCDC34, and ODF2.
In other implementations, the at least one autoantigen is selected from the group consisting of: STMN4, ODF4, RBPJ, AMY2A, EPCAM, ZNF688, CSF3, RAD51AP1, PSKH1, LENG1, S1PR3, LYSMD1, FAM76A, CDR2L, and CCDC130. In some embodiments, where the at least one autoantigen is EPCAM, the method further comprises detecting the presence of at least one other autoantibody against EDG3 or CSF3. Detecting the presence of the at least one other autoantibody against EDG3 or CSF3 in addition to the presence of at least one autoantibody against EPCAM in the biological sample indicates that the subject is false positive for having cancer or an autoimmune condition.
In some aspects, the method comprises detecting at least two, at least three, at least four, or at least five autoantibodies against at least two, at least three, at least four, or at least five autoantigens selected from the group consisting of: STMN4, ODF4, RBPJ, AMY2A, EPCAM, ZNF688, CSF3, RAD51AP1, PSKH1, LENG1, S1PR3, LYSMD1, FAM76A, CDR2L, and CCDC130. Detecting the presence of the at least two, at least three, at least four, or at least five autoantibodies indicates the subject is false positive for having cancer or an autoimmune condition. In certain embodiments, wherein the method comprises detecting at least ten autoantibodies against at least ten selected from the group consisting of: STMN4, ODF4, RBPJ, AMY2A, EPCAM, ZNF688, CSF3, RAD51AP1, PSKH1, LENG1, S1PR3, LYSMD1, FAM76A, CDR2L, and CCDC130, detecting the presence of the at least ten autoantibodies indicate the subject is false positive for having cancer or an autoimmune condition. In a particular embodiment, the method comprises detecting the presence of at least one autoantibody against STMN4, ODF4, RBPJ, AMY2A, EPCAM, ZNF688, CSF3, RAD51AP1, PSKH1, LENG1, S1PR3, LYSMD1, FAM76A, CDR2L, and CCDC130 in the biological sample. Detecting the presence of at least one autoantibody against STMN4, ODF4, RBPJ, AMY2A, EPCAM, ZNF688, CSF3, RAD51AP1, PSKH1, LENG1, S1PR3, LYSMD1, FAM76A, CDR2L, and CCDC130 in the biological sample indicate the subject is false positive for having cancer or an autoimmune condition.
In certain implementations, the method comprises detecting the presence of at least one autoantibody against at least one autoantigen selected from the group consisting of: STMN4, ODF2, RBPJ, AMY2A, EPCAM, and ZNF688 in the biological sample, wherein the presence of the at least one autoantibody indicates the subject is false positive for having cancer or an autoimmune condition. In some aspects, the method comprises detecting the presence of at least two, at least three, at least four, or at least five autoantibodies against at least two, at least three, at least four, or at least five autoantigens selected from the group consisting of: STMN4, ODF2, RBPJ, AMY2A, EPCAM, and ZNF688 in the biological sample. Detecting the presence of the at least two, at least three, at least four, or at least five autoantibodies indicates the subject is false positive for having cancer or an autoimmune condition. In particular embodiments, the method comprises detecting the presence of at least one autoantibody against STMN4, ODF2, RBPJ, AMY2A, EPCAM, and ZNF688 in the biological sample. Detecting the presence of the at least one autoantibody against STMN4, ODF2, RBPJ, AMY2A, EPCAM, and ZNF688 in the biological sample indicates the subject is false positive for having cancer or an autoimmune condition.
A kit for diagnosing and verifying an autoimmune condition or cancer in a subject is further described. The kit comprises at least one antigen or antibody targeting a cancer marker, an autoimmune marker, or both and a false positive indicator, the false positive indicator comprising at least one autoantigen listed in Table 1. In some embodiments, the kit further comprises at least one labeled secondary antibody or antigen for detecting the binding of the at least one antigen or antibody with the cancer marker and/or the autoimmune marker; and at least one labeled secondary antibody or antigen detecting binding of at least one autoantibody to the false positive indicator. In some aspects, wherein the cancer marker is P53. In other aspects, the autoimmune marker is antinuclear antibodies.
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.
Detailed aspects and applications of the invention are described below in the drawings and detailed description of the invention. Unless specifically noted, it is intended that the words and phrases in the specification and the claims be given their plain, ordinary, and accustomed meaning to those of ordinary skill in the applicable arts.
In the following description, and for the purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the various aspects of the invention. It will be understood, however, by those skilled in the relevant arts, that the present invention may be practiced without these specific details. It should be noted that there are many different and alternative configurations, devices, and technologies to which the disclosed inventions may be applied. The full scope of the inventions is not limited to the examples that are described below.
The singular forms “a,” “an,” and “the” include plural referents unless the context clearly dictates otherwise. Thus, for example, reference to “a step” includes reference to one or more of such steps.
As used herein, the term “false positive” refers a test result which incorrectly indicates that a particular condition or attribute is present. Accordingly, the “false positive indicator”, in some aspects, refers to a biomarker that indicates the corresponding positive test result incorrectly indicates the presence of a particular condition or attribute, for example, cancer or an autoimmune condition.
Disclosed herein are a selection of common IgG autoantibodies that can be found in healthy individuals and thus are suitable target for identifying false positives in cancer or autoimmune condition diagnostics methods and assays.
Autoantibodies can be broadly divided into two types: 1) pathogenic autoantibodies which contribute to various immune-mediated diseases and 2) common autoantibodies which are found in apparently healthy individuals. While pathogenic autoantibodies can lead to autoimmune diseases, common autoantibodies can bind to a variety of microbial components, thereby providing first line of defense against infections (Elkon & Casali, 2008). They can also recognize self-antigens which help in B cell repertoire development and homeostasis of the immune system. Some of these common autoantibodies occur frequently enough to confound studies intended to find disease-related autoantibodies.
The number of unique IgG autoantibodies in healthy individuals increased with age from infancy to adolescence and then plateaued. This observation suggests that while response to infectious agents (and maybe vaccines) might contribute to autoantibodies through molecular mimicry, this mechanism does not appear to continue to accumulate autoantibodies throughout life. Gender did not appear to play a role in autoantibody production in healthy individuals. This stands in contrast to the observation that autoimmune diseases disproportionally affect females compared to males because male-predominant autoimmune disease is associated with acute inflammation, whereas female-predominant autoimmune disease is associated with antibody-mediated pathology. Several common autoantibodies co-occurred frequently. This could occur if the same antibody recognized two different proteins that share a common epitope. Other possibilities include sharing common HLA haplotypes or playing similar biological roles that lead to escape from tolerance. It is notable that the targets of several of the co-occurring antibodies play roles in stem cell proliferation and differentiation (EPCAM, EDG3 and CSF3) and two others play roles in DNA damage repair (PML and PSMD2). The meaning of this is not clear, but it occurred frequently enough (Phi correlation coefficient >0.6) that it is worth further investigation.
Viral proteins with sequences similar to a human protein may initiate cross-reactive antibodies leading to autoimmunity. There are around 20 autoimmune diseases reported in literature where autoantibodies are generated due to cross-reactivity to infectious agent proteins. Some of the common autoantibodies may be a result from cross-reactivity from anti-viral antibodies, albeit without causing subsequent pathology. The typical length of linear epitope of antibodies ranges from 7 to 9 amino acids and hence these specific matches have the potential to elicit cross-reactive antibodies. The fact that these matches occur significantly more frequently between viral proteins and common autoantigens but less frequently for unreactive proteins on the microarrays further suggests the role of molecular mimicry in common autoantibody elicitation.
The intrinsic properties of a protein, such as its chemical and structural complexity, can impact its antigenicity. Gene Set Enrichment Analysis (GSEA) revealed that common autoantigens tended to favor more basic, hydrophilic with fewer aromatic amino acids. In addition, common autoantigens were also found to be more flexible and have more beta-turns. Flexibility is a property that can help the polypeptide chain to bind easily to immunoglobulin compared to a stiff polypeptide chain. Also, beta-turns can be a potential site for antibody binding as the peptide chain reverses its direction at beta-turns with side chain projected outwards.
Accessibility of autoantigens to circulating autoantibodies is critical to autoimmune disease pathology. In systemic autoimmune diseases, a majority of the target antigens are intracellular molecules and therefore not normally accessible to the B cells or antibodies. Only after excessive cell death or ineffective clearance of apoptotic debris do these intracellular autoantigens become available for immune complex formation. In Wegener's granulomatosis, the autoantigen is an intracellular protease that becomes accessible to the autoantibodies only after an infection triggers translocation of the protease to the surface. Similarly, the autoantigen in Goodpasture's syndrome, normally ensconced in the basal membranes of alveolar capillaries, becomes accessible to the antibodies after an environmental insult to the capillaries, leading to pulmonary hemorrhage. A majority of the common autoantigens identified were located exclusively at intracellular sites, which make them inaccessible to circulating autoantibodies. Some of the common autoantigens are organ/tissue-specific and predominately expressed in the testis and brain, which are isolated from the immune system by the blood-testis or blood-brain barriers, respectively. No obvious form of sequestration was identified for the remaining autoantigens although this cannot be ruled out.
Thousands of studies over the past decade have investigated autoantibodies as potential biomarkers for disease risk assessment, diagnosis, and prognosis. Given the prevalence observed for these common autoantibodies in healthy individuals, in some cases exceeding a quarter of all individuals, they will be frequently encountered in such studies and may confound them as false positives. Table 1 lists 77 common autoantibodies. An examination of the AAgAtlas & PubMed revealed that 20 of our 77 common autoantibodies have been reported as disease-related biomarkers. These 20 common autoantibodies are STMN4, ODF2, RBPJ, AMY2A, EPCAM, ZNF688, CSF3, S1PR3, CDR2L, RASSF1, SPAG8, TRAP1, SART1, GATA2, PML, SOX2, PELI1, TAX1BP1, JUN, and PAK1.
As shown in the examples, a subset of common autoantibodies may be used as indicators of a false positive result in the cancer or autoimmunity screening. Thus, disclosed herein is a method of increasing the accuracy of a cancer or an autoimmune condition screening in a subject suspected of having cancer or an autoimmune condition by identifying at least one false positive indicator. The method comprises providing a biological sample from the subject; detecting the presence of one or more autoantibodies against at least one autoantigen listed in Table 1 in the biological sample, and identifying at least one false positive indicator upon detecting the presence of the one or more autoantibodies.
The biological sample from the subject can be any type of biological sample, for example, a tissue sample, a saliva sample, or a blood sample, including a fraction thereof such as a plasma sample or a serum sample. In certain embodiments, the biological sample is a serum sample or a plasma sample.
The step for detecting the presence of one or more autoantibodies against at least one autoantigen may be performed by methods well-established in the art, for example by immunoassay, where autoantibody capture molecules are immobilized on a solid support. For example, the presence of one or more autoantibodies against at least one autoantigen is detected using Enzyme Linked Immunosorbent Assay (ELISA), Protein microarrays, Western Blot or beads based immunoassay. In some implementation, the autoantigens are immobilized on a solid surface, in the immunoassay. In other implementation, peptide methods, such as peptide arrays and phage display, are used to detect the presence of one or more autoantibodies against at least one autoantigen.
In some implementations of the method of increasing the accuracy of a cancer or an autoimmune condition screening, the at least one autoantigen is selected from the group consisting of: AMY2A, SPAG8, GTSE1, MAK, RNF138, PMFBP1, CCDC34, ODF2, CSF3, CAPN3, MYLK2, TRIM29, SOX2, and STMN4. Where the biological sample is from the subject's pancreas, the at least one autoantigen is AMY2A. Where the biological sample is from the subject's lungs, the at least one autoantigen is CSF3. Where the biological sample is from the subject's salivary gland, nerve tibial, or brain, the at least one autoantigen is SOX2. Where the biological sample is from the subject's esophagus, skin, or vagina, the at least one autoantigen is TRIM29. Where the biological sample is from the subject's skeletal muscle, the at least one autoantigen is CAPN3 and/or MYLK2. Where the biological sample is from the subject's testis, the at least one autoantigen is selected from the group consisting of: SPAG8, GTSE1, MAK, RNF138, PMFBP1, CCDC34, and ODF2.
In other implementations of the method, the at least one autoantigen is selected from the group consisting of: STMN4, ODF4, RBPJ, AMY2A, EPCAM, ZNF688, CSF3, RAD51AP1, PSKH1, LENG1, S1PR3, LYSMD1, FAM76A, CDR2L, and CCDC130. In some embodiments, where the at least one autoantigen is EPCAM, the method further comprises detecting the presence of at least one other autoantibody against EDG3 or CSF3, wherein the presence of the at least one other autoantibody against EDG3 or CSF3 in the biological sample indicates the subject is false positive for having cancer or an autoimmune condition.
In some aspects, the method of increasing the accuracy of a cancer or an autoimmune condition screening comprises detecting at least two, at least three, at least four, or at least five autoantibodies against at least two, at least three, at least four, or at least five autoantigens selected from the group consisting of: STMN4, ODF4, RBPJ, AMY2A, EPCAM, ZNF688, CSF3, RAD51AP1, PSKH1, LENG1, S1PR3, LYSMD1, FAM76A, CDR2L, and CCDC130. Detecting the presence of the at least two, at least three, at least four, or at least five autoantibodies indicates that the subject is false positive for having cancer or an autoimmune condition. In certain embodiments, the method comprises detecting at least ten autoantibodies against at least ten selected from the group consisting of: STMN4, ODF4, RBPJ, AMY2A, EPCAM, ZNF688, CSF3, RAD51AP1, PSKH1, LENG1, S1PR3, LYSMD1, FAM76A, CDR2L, and CCDC130. Detecting the presence of the at least ten autoantibodies indicate that the subject is false positive for having cancer or an autoimmune condition. In a particular embodiment, wherein the method comprises detecting the presence of at least one autoantibody against STMN4, ODF4, RBPJ, AMY2A, EPCAM, ZNF688, CSF3, RAD51AP1, PSKH1, LENG1, S1PR3, LYSMD1, FAM76A, CDR2L, and CCDC130 in the biological sample, the presence of at least one autoantibody against STMN4, ODF4, RBPJ, AMY2A, EPCAM, ZNF688, CSF3, RAD51AP1, PSKH1, LENG1, S1PR3, LYSMD1, FAM76A, CDR2L, and CCDC130 in the biological sample indicate the subject is false positive for having cancer or an autoimmune condition.
In certain implementations, the method comprises detecting the presence of at least one autoantibody against at least one autoantigen selected from the group consisting of: STMN4, ODF2, RBPJ, AMY2A, EPCAM, and ZNF688 in the biological sample. Detecting the presence of the at least one autoantibody against STMN4, ODF2, RBPJ, AMY2A, EPCAM, or ZNF688 indicates the subject is false positive for having cancer or an autoimmune condition. In some aspects, the method comprises detecting the presence of at least two, at least three, at least four, or at least five autoantibodies against at least two, at least three, at least four, or at least five autoantigens selected from the group consisting of: STMN4, ODF2, RBPJ, AMY2A, EPCAM, and ZNF688 in the biological sample, wherein the presence of the at least two, at least three, at least four, or at least five autoantibodies indicates the subject is false positive for having cancer or an autoimmune condition. In particular embodiments, the method comprises detecting the presence of at least one autoantibody against STMN4, ODF2, RBPJ, AMY2A, EPCAM, and ZNF688 in the biological sample, wherein the presence of the at least one autoantibody against STMN4, ODF2, RBPJ, AMY2A, EPCAM, and ZNF688 in the biological sample indicates the subject is false positive for having cancer or an autoimmune condition.
Also described herein are a diagnostic kit and array for diagnosing and verifying an autoimmune condition or cancer in a subject. In some aspects, the kit and the array comprise at least one false positive indicator, wherein the false positive indicator comprising at least one autoantigen listed in Table 1. For example, the kit comprises at least one antigen or antibody targeting a cancer marker, an autoimmune marker, or both and a false positive indicator, the false positive indicator comprising at least one autoantigen listed in Table 1. In some embodiments, the kit further comprises at least one labeled secondary antibody or antigen for detecting the binding of the at least one antigen or antibody with the cancer marker and/or the autoimmune marker; and at least one labeled secondary antibody or antigen detecting binding of the false positive indicator to the at least one autoantigen. The diagnostic array comprises at least one antigen or antibody targeting a cancer marker, an autoimmune marker, or both and a false positive indicator, the false positive indicator being selected from at least one autoantigen listed in Table 1. In some aspects, wherein the cancer marker is P53. In other aspects, the autoimmune marker is antinuclear antibodies.
The present invention is further illustrated by the following examples that should not be construed as limiting. The contents of all references, patents, and published patent applications cited throughout this application, as well as the Figures, are incorporated herein by reference in their entirety for all purposes.
Autoantibody profiles for 272 healthy subjects from 9 case-control studies were compiled (Table 2). There were more females than males, 195 vs. 67, because several studies focused on female-specific diseases such as breast and ovarian cancers. These studies were diverse in terms of subject ages, ranging from infancy to adulthood, with most above 50 years old. Antibodies against 8,282 unique human proteins were studied; although, the number of proteins studied for each subject varied by study (
To minimize the effect of study heterogeneity, sample size based weighted prevalence was calculated as the sum of individual prevalence of antibody in each study multiplied by the sample size of the study. For the healthy subjects, 77 autoantibodies occurred frequently and had a weighted prevalence between 10%-47% (Table 1). These were termed as common autoantibodies. Antibodies against STMN4, ODF2, RBPJ, AMY2A, EPCAM, and ZNF688 showed the highest prevalence (Table 1).
160 healthy subjects from five studies that included age information (Studies I, II, IV, VI, VII, Table 2) were divided into five age groups based on human development stages to examine the time course of autoantibody development. The infant and early childhood age group (0-6 years) had the least number of autoantibodies. The number increased in the middle and late childhood age group (6-12 years) and then plateaued (
If these common autoantibodies observed in the healthy subjects were elicited through common non-pathogenic mechanisms, then they should also occur at similar frequencies in their matched disease cohorts. Indeed, the 77 common autoantibodies occurred at similar frequencies in diseased cohorts to those in healthy cohorts (
To answer whether any of these common autoantibodies were related to each other (in other words, was there any concordance among them or were their occurrences independent), the common autoantigens were analyzed pairwise to determine if any occur together in healthy individuals at frequencies greater than chance alone (
EDG3
EPCAM
0.83
0.67
6
<0.001
PML
PSMD2
0.73
0.63
3
<0.001
CSF3
EPCAM
0.67
0.65
6
<0.001
To understand the extent that common autoantibodies observed in our study resulted from cross-reactivity of antibodies induced by viral infection, the sequence similarities between viral proteins and common autoantigens was examined. As these autoantibodies developed early in age and did not change after adolescence, respiratory and common viruses found in children of the US were included in the analysis (Table 4). In order to avoid redundancy and false positives, duplicate proteins and consecutive amino acid repeats were removed from viral proteomes (
To determine whether any intrinsic biochemical and structural properties of the target antigens were responsible for common autoantibodies production, various properties were examined by comparing our list of common autoantigens to all 8,282 proteins using Gene Set Enrichment Analysis (GSEA). The 77 common autoantigens were significantly enriched with proteins having low aromaticity (NES or normalized enrichment score: −2.13, P<0.001), low hydrophobicity (NES: −2.01, P<0.001), high isoelectric point (NES: 1.58, P=0.018), high fraction of amino acids in β-turns (NES: 1.95, P=0.04), high Karplus & Schulz flexibility (NES: 4.40, P<0.001), high Parker hydrophilicity (NES: 2.33, P<0.001), and high Chou & Fasman β-turn score (NES: 2.61, P<0.001) (
The discovery of common autoantibodies in healthy individuals raised the question about why these antibodies do not lead to autoantibody-mediated pathology. A primary requirement for such pathology is the formation of immune complexes. The subcellular localization of the common autoantigens to see if they were antibody accessible. They were divided into three broad categories: “intracellular”, “cell membrane”, and “secreted” (Table 6). The localization of an autoantigen can belong to one or more of these 3 categories. 55 among 70 common autoantigens were located exclusively at intracellular sites. The percentage of common autoantigens with “intracellular only” subcellular localization was significantly higher than that for all the proteins studied on the microarrays (78% vs. 54%, P<0.001) (
Tissue-specific gene expression can impact autoantigen exposure to circulating autoantibodies and the potential to trigger autoimmune disease. The data from GTEx, which is a public resource portal for tissue-specific gene expression in multiple human tissues, was used. In the GTEx dataset, transcripts encoding for 14 common autoantigens were organ/tissue-specific (defined as having log2 ((organ expression)/(mean expression in all other organs)>3) (
a. Datasets
The healthy subjects included in this study were originally included in 9 different case-control studies (Table 2). These studies were all conducted in our lab; 5 of them were published (Table 2). The serum samples were collected from various parts of the USA and the UK. The goal of the original studies was to discover biomarkers of various cancers and autoimmune diseases by comparing the prevalence of antibodies present in diseased and healthy subjects. The presence of antibody was determined using protein microarrays that displayed thousands of human proteins as potential targets. Serum samples were probed on protein microarrays followed by a secondary antibody with a fluorophore tag specific for human IgG. Microarrays were scanned by a laser scanner. The microarray images from the 9 studies were qualitatively examined to identify protein targets that serum antibodies bound using Array-Pro Analyzer 6.3 (Media Cybernetics). All proteins were not probed by all samples included in our analysis (
b. Age and Gender Comparison
To understand the effect of age on autoantibody counts in healthy individuals, studies having both male and female subjects with age information were used (Studies I, II, IV, VI, VII, Table 2). A total of 160 subjects were divided into five age groups based on human development stages. The groups were 0 to 6 years old (infancy & early childhood), 6 to 12 years old (middle & late childhood), 12 to 18 years old (adolescence), 18 to 51 years old (early adulthood) and 51 to 84 years old (late adulthood). The number of autoantibodies in each subject was plotted using GraphPad Prism by age groups. To understand the effect of gender on autoantibody counts in healthy individuals, studies having both male and female subjects with matched age were used (Studies I, II, IV, VII, Table 2). The subjects were divided into male and female groups. The number of autoantibodies found in each subject was plotted using GraphPad Prism. The weighted prevalence of each autoantibody was calculated for male and female separately. The method of weighted prevalence calculation is described in the “Quantification and statistical analysis” subsection. Prevalence values for the 77 most common autoantibodies were plotted as a population pyramid using GraphPad Prism. A paired t-test was performed to determine the significance of the prevalence difference between genders. Pearson correlation of common autoantibodies frequency in diseased and healthy cohorts were plotted using Python seaborn package.
c. Correlation of Common Autoantibodies
As the presence of common autoantibodies were measured on a binary scale, a phi correlation coefficient was computed to measure associations between autoantibodies. Specifically, for each pair of antibodies, a phi correlation coefficient was computed for each study, and multiple phi correlation coefficients across different studies were combined into a single phi correlation coefficient using the R meta package. The R “pheatmap” package was then used to produce correlation heatmap plots for both healthy and diseased cohorts (
d. Sequence Similarity with Viral Proteins
The proteomes of respiratory and common viruses found in children of the US were downloaded from UniProt as a FASTA file. All the common human viruses were included except sexually transmitted ones as common autoantibodies that develop early in age and then plateau (Table 4). CD-HIT was employed to remove duplicate sequences in the file (sequence identity cut-off: 1) (Huang et al., 2010). The sequences were then segregated into 14-mer peptides using a Python script (sliding window: 1) and consecutive amino acid repeats (3 or more) were removed. The sequences of all the human proteins analyzed on microarrays were retrieved from DNASU (https://dnasu.org) and split into two sequence databases, “common autoantigens” and “unreactive proteins”. The “unreactive proteins” database comprises proteins from the microarrays without any autoantibody responses. Repeats and low-complexity regions were masked using BLAST+(Basic Local Alignment Search Tool, version 2.10.1) package “segmasker”. A protein-protein BLAST was run with the following parameters, “-ungapped, -db_hard_mask 21, -comp_based_stats F, -evalue 10”, between viral 14-mer peptides and “common autoantigens”. Another protein-protein BLAST was run between viral 14-mer peptides and “unreactive proteins” with similar parameters except adjusted “-evalue 593.89” to compensate for the bigger size of unreactive proteins database (Effective search space of “unreactive proteins” and “common autoantigens” databases were 15,970,464 and 268,912, respectively). The total number of amino acids matches higher than the threshold (7 ungapped amino acids match) was calculated for both databases and compared with the total number of amino acids in each database using a chi-square test (
e. Biochemical and Structural Properties
Biopython (version 1.75) module Bio.SeqUtils.ProtParam for Python (version 3.7.6) was used to calculate the values of aromaticity, isoelectric point, hydrophobicity, the fraction of amino acids in sheets and turns for each protein (Table 8). Secondary structure and antigenicity prediction methods from Immune Epitope Database (IEDB) were also used. Command-line tools from IEDB analysis resource were employed to calculate the values of Chou & Fasman beta-turn, Emini surface accessibility, Karplus & Schulz flexibility, and Parker hydrophilicity across the proteins, which were then averaged for each protein. The computed biochemical property values were used for the enrichment analysis on the identified common autoantigens using Gene Set Enrichment Analysis (GSEA) “GSEAPreranked” package (version 4.2) (Subramanian et al., 2005).
f. Subcellular Localization and Tissue Expression
All 8,282 proteins were used to query the UniProt database for subcellular localization (downloaded in December 2020), among which 6,875 proteins had subcellular localization data available in the database (Table 6). Some of the proteins were found simultaneously in more than one location, and hence, seven groups were created to segregate the proteins based on their subcellular localization profiles. Proteins that were found only in one subcellular location were put into “intracellular only”, “cell membrane only” and “secreted only” groups. Proteins that were found in two subcellular locations were put into “intracellular & cell membrane”, “cell membrane & secreted” and “secreted & intracellular” groups. Proteins that were found simultaneously inside the cell, in the cell membrane, and outside the cell were put into “intracellular, cell membrane and secreted” group. P value was calculated to assess the statistical significance of difference in fractions of “intracellular only” proteins for all proteins on the microarrays and for common autoantigens using the proportion test.
All 8,282 proteins were mapped to the Ensembl IDs using “BiomaRt” package available for R (version 3.5.0). The Ensembl IDs were used to identify the protein of interest in the Genotype-Tissue Expression (GTEx, version 8) dataset. The gene expression levels in 52 human tissue types, measured in transcripts per million (TPM), were downloaded from GTEx. Expression values for tissue types belonging to the same organ were averaged. Differentially expressed genes for each organ/tissue were identified using edgeR package for R (version 3.6.2) with a cutoff of Log2 (fold change)>3 to determine organ/tissue specificity, where the fold change for each gene was calculated by dividing the TPM value in a particular organ/tissue by the mean TPM values in all other organs/tissues. The Log2-scaled fold changes across the organs/tissues for each gene were standardized to the Z scores for data visualization. The Z score profiles were displayed in a heatmap with correlation-based average-linkage clustering by using the Python seaborn package.
g. Quantification and Statistical Analysis
i. Weighted Prevalence
Due to the heterogeneous number of proteins and subjects being analyzed in each study, we computed the weight for the jth antibody as {circumflex over (p)}j=Σi=1kwijpij/Σi=1kwij, where pij=xij/nij is the prevalence, xij is the total number of positive signals found for the jth antibody in the study i, and nij is the number of samples for the jth antibody in the study i, and k is the number of studies. Here, wij=(vij+τj2)−1 is the inverse variance-weighting which accounts for the heterogeneous effects between studies (Borenstein et al., 2010), where vij=nij/(pij(1−pij)), τj2=(Qj−k+1)/Uj if Qj>k−1 or τj2=0 otherwise, Qj=Σi=1kvij(pij−pj)2, Uj=(k−1)(
ii. Age and gender comparison
The significance of increase in the autoantibody counts among the five age groups was calculated using the Welch's t-test while the significance of difference in autoantibody counts between the male and female groups was calculated using a two-sample unpaired t-test.
iii. Biochemical and structural properties
The “GSEAPreranked” package available in GSEA software returned P values of “0.0” when the number of permutations was set to 1000, as it cannot calculate very small P values. Another R package named “fgsea” was used to calculate the very small P values with number of permutations set to 10,000 for more accurate calculation. To adjust multiple comparisons, we computed false discovery rate (FDR) adjusted P value using the “p.adjust” function in the R stats package.
This application claims the benefit of and priority to U.S. provisional patent application 63/347,560, filed May 31, 2022, the entirety of the disclosure of which is hereby incorporated by this reference thereto.
This invention was made with government support under U01 CA214201 awarded by the National Institutes of Health. The government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
63347560 | May 2022 | US |