PreImplantation Factor in endometriosis: A potential role in inducing immune privilege for ectopic endometrium

Endometriosis is a chronic, benign disease affecting 10% of women in their reproductive years and characterized by the presence of stromal and epithelial cells outside the uterine cavity [ 1 ]. Endometriosis lesions are found predominantly in the pelvis, potentially through the reflux of viable endometrial cells during retrograde menstruation [ 1 , 2 ]. However, recent evidence suggests that the pathogenesis of endometriosis requires more than retrograde menstruation, underlying the emerging role of stems cells and the immune response [ 3 – 5 ]. For example, the cyclic regeneration of a healthy endometrium depends on stem cells. Both endometrial- and bone marrow-derived stem cells may migrate to ectopic lesions and contribute to lesion growth [ 3 ]. Irrespective of their source however, once endometrial cells are present they secrete chemokines that stimulate immune cell infiltration [ 6 ]. Both the endometrial and infiltrating immune cells produce inflammatory cytokines, such as TNFα, that further stimulate a cascading inflammatory response. This creates a regulatory feed forward loop influencing both the progression and symptomology of the disease [ 7 , 8 ], resulting, in a unique microenvironment that contributes to these lesions being able to evade immune surveillance [ 4 ]. One of the key regulators of immune processes in endometriosis are regulatory T cells (Tregs) derived from CD4 lineage [ 9 ]. Tregs are produced naturally in the thymus and express the forkhead box P3 transcription factor (Foxp3 + ). Cytokine-induced increase of Foxp3 expression drives Tregs differentiation/activation and suppresses the response of effector T cells by inhibiting dendritic, or the other antigen presenting cells from triggering effector T cell proliferation [ 10 ]. Additionally, Tregs induce immune tolerance by production of IL-10, Transforming Growth Factor-β, and anti-inflammatory cytokines that inhibit T helper cell activation. Not surprisingly, this altered immunological microenvironment may lead to tumor growth [ 11 – 13 ]. On the contrary, the absence or depletion of Tregs lead to multi-systemic autoimmunity in mice and humans [ 14 ]. In line with a role for a modified immunity in the pathogenesis of endometriosis, CD4 + /FoxP3 + Tregs are present in endometriotic lesions [ 15 , 16 ]. PreImplantation Factor (PIF) is a fifteen amino acid linear peptide secreted by viable embryos [ 17 , 18 ] that targets immune cells directly [ 19 , 20 ]. PIF`s essential roles in pregnancy begins by priming the endometrium for implantation and continue by promoting trophoblast invasion through local and systemic immune modulation [ 21 – 25 ]. Not surprisingly, a synthetic analog of PIF (sPIF) was successfully tested in multiple non-pregnant and pregnant animal models including autoimmune, transplantation, radiation induced, and brain diseases [ 26 – 33 ]. For example, sPIF is neuroprotective in neonatal rats through modulation of microRNA let-7 and Protein Kinase A and Protein Kinase C signalling in neuronal and microglial cells [ 29 , 30 ]. The neuroprotective and immune modulatory effect was also observed in adult multiple sclerosis models [ 27 , 32 ] and the regulation of immune function was observed in graft-versus-host, inflammatory induced fetal loss, autoimmune diabetes, and radiation induced injury models [ 26 , 28 , 31 , 33 ]. Furthermore, sPIF is currently being tested in an ongoing University-sponsored, FDA Fast-Track approved clinical trial for an autoimmune disease manifested in liver inflammation ( NCT02239562 ). Given that PIF is an essential and safe pregnancy product and pregnancy itself suppresses endometriosis, we aimed to evaluate the role of PIF within the endometriotic environment.

To determine the role of PIF in endometriosis, tissue samples were obtained from women with or without endometriosis including the ectopic ovarian endometriomas and peritoneal implants during the proliferative and secretory phases (see S1 Table for characterization of the patients). We did not detect PIF positive cells in either the epithelial, or stromal cells of healthy controls in the proliferative and secretory phases ( Fig 1A and 1B ). However, we detected PIF positive cells in the epithelial cells of ectopic endometria, predominantly in the secretory phase documenting for the first time expression of PIF outside of pregnancy. We hypothesized that PIF expression may be induced as a protective mechanism so we tested the effect of sPIF using an in-vitro system. (A) PIF positive cells were not detected in tissues derived from healthy patients (upper panel). PIF positive cells were detected in epithelial compartment of ectopic endometrial tissues (lower panel). ( B) Semi-quantitative analysis of PIF staining. We detected increased number of PIF positive cells with high HSCORE in epithelial compartment of ectopic tissues only. ( C) and (D) Diverse influence of sPIF on cell survival. ( C) sPIF treatment of epithelial ectopic cells (EM ‘osis −) resulted in a significant decrease of cell viability. However, treatment of epithelial eutopic cells (EM E6/E7 --) resulted in a significant increase of cell viability. ( D) The viability of ectopic (−) and eutopic (--) stromal cells did not show a significant change in response to sPIF treatment. *p<0.05; **; p<0.01; ***p<0.001. Scale bar 20 μm. To elucidate the potential role of PIF in endometriosis we used epithelial cell lines isolated from the eutopic endometrium (EM E6/E7) [ 39 ] and ectopic endometriomas (EM’osis) [ 40 ] and primary endometrial stromal cells (ESC) from women with and without endometriosis. We tested sPIF treatment in ascending doses in epithelial cells first. sPIF treatment significantly decreased cell viability compared to control ( Fig 1C solid line). Furthermore, and supporting PIF`s role during embryo implantation, sPIF also increased the viability of eutopic endometrium cells ( Fig 1D dashed line) [ 20 ]. In line with our previous reports and observation in human tissues ( Fig 1A ), sPIF treatment did not alter stromal cell viability ( Fig 1D ) [ 22 ]. Together, both PIF expression and sPIF treatment influence epithelial, but not stromal endometriotic cells suggesting a novel pathogenesis of endometriosis. To further dissect the underlying mechanisms, we performed a global gene expression array from sPIF treated ectopic endometrial epithelial cells. Having a screening approach in mind, we performed a global gene array from sPIF`s treated epithelial ectopic cell lines ( Fig 1 ). In line with previous reports [ 23 ], sPIF treatment resulted in modulation of multiple signaling pathways as shown by string analysis (see Fig 2A and S2 Table for details). We detected significant changes in channels activity pathways with cation and potassium channel activity being crucial, suggestion they play a crucial role in sPIF function, especially as potassium channels have previously been identified as a sPIF target [ 44 ]. Multiple pathways involved in neuronal development, plasticity, and protection were influenced by PIF treatment, in line with previous reports of sPIF neuroprotective effects [ 29 , 30 , 45 ]. Interestingly, changes in T-cell receptor signalling were also detected (see Fig 2A green arrow and S3 Table for details) and T cell signalling has previously been related to endometriosis pathogenesis [ 5 , 46 ]. Therefore, we performed a detailed heat map analysis and detected multiple genes modulated by sPIF treatment ( Fig 2B ). Of special interest is the FoxP3 gene as FoxP3 + Tregs cells are not decreased in endometriotic lesions during the secretory phase, which may contribute to their survival [ 5 ]. Given the importance of FoxP3 in the pathogenesis of endometriosis [ 5 , 16 ], we aimed to confirm a potential PIF interaction with FoxP3 signaling. (A) Gene set network showing the relations between the significant gene sets after sPIF treatment of epithelial ectopic cells. Every node represents a gene set and every edge reflects the intersection between two sets. Node size reflects the number of genes in each gene set; node fill color reflects the p-value corrected for intersections with the color going from blue to red with decreasing p-value; node border reflects network prominence, how important a particular gene set is compared to others. Edge width reflects the relative size of the intersection between two gene sets. Edge direction, as indicated by the arrow heads, is towards to more significant gene set. Repeated arrows indicate subset relations. ( B) Detailed analysis of the T cell receptor signaling pathways showing multiple up (red) and down (blue) regulated genes (6 samples: 3 sPIF treated and 3 control samples). FoxP3 + cells were detected in both eutopic and ectopic tissues ( Fig 3A and 3B ) and in line with previous reports, we detected an increased number of FoxP3 + cells in the ectopic endometria, particularly during the secretory phase [ 15 , 16 ]. To further underline the potential interaction between PIF and FoxP3 + cells, we tested PIF`s effects in a well-defined system of peripheral non-pregnant blood mononuclear cells (PBMCs). Notably, we reported on FITC-PIF binding to CD14 + and CD3 + cells in both pregnant and non-pregnant patients previously [ 21 ]. In addition binding to activated CD4 + cells were documented. Using the same method of flow cytometry, we detected a dose dependent increase in the binding of FITC-PIF to CD4 + /CD25 + /FoxP3 + cells ( Fig 3C ). This observation is in line with close proximity of FoxP3 + cells (stromal compartment) and PIF + cells (epithelial compartment) of the ectopic endometria during the secretory phase ( Fig 3D ). Together, these observations suggest that the expression of PIF may contribute to the increased FoxP3 expression in these lesions. (A) FoxP3 positive cells were detected in both healthy (left panel) and endometriotic (right panel) tissues in stromal and epithelial compartments. ( B) Semi-quantitative analysis of FoxP3 staining demonstrates increased percentage of positive cells with higher HSCORE in stromal compartment of ectopic tissues. ( C) Dose dependent and specific PIF binding to Tregs. ( D) Epithelial PIF positive cells (green asterisk) in close proximity to FoxP3 positive (red asterisk) stromal cells in ectopic tissue. ( E) sPIF specific modulation of cell viability ( Fig 1C ) is abolished in the presence of TNFα. ( F) New Hypothesis of immune privileged environment in endometriosis. PIF re-expression leads to recruitment of FoxP3 positive cells into the lesion. The production of pro-inflammatory cytokines such as TNFα leads to a diverged PIF effects resulting in cell survival. Systemic sPIF application may lead to absents of TNFα production and FoxP3 recruitment leading to cell death due to PIF re-expression in the endometriotic lesions. Scale bar 20 μm. Finally, we aimed to detect sPIF effects in the inflammatory endometriotic environment. Notably, multiple chemokines affect FoxP3 + cell function and FoxP3 cells regulate the survival of ectopic endometrial cells [ 9 ]. We decided to test sPIF effect on cell viability in the presence of TNFα as TNFα was identified as one of the pivotal chemokines involved in endometriosis and cell death [ 4 ]. Indeed, in the presence of TNFα, sPIF`s specific effect on cell viability in epithelial ectopic cells ( Fig 1C ) was abolished suggesting the pivotal involvement of TNFα ( Fig 3E ) in sPIF`s mediated effects. This observation supports the notion of a disturbed balance between pro-inflammatory and anti-inflammatory factors in endometriosis affecting the cellular survival [ 47 ]. Combining this data together therefore, we hypothesize that PIF expresses in the epithelial compartment of ectopic endometria resulting in the recruitment of FoxP3 + cells into the stromal compartment. The additional influx of chemokines and pro-inflammatory factors such as TNFα result in divergent sPIF effect on these cells creating a positive feed-back loop. It may lead to cellular survival which again contribute to a significant control of the immune privilege environment in the peritoneal cavity ( Fig 3F ). Thus, we hypothesize that PIF may be a crucial factor contributing to the pathogenesis of endometriosis on the local level. However, the immune modulatory effect of sPIF applied peripherally as a therapeutic option may lead to a decreased inflammatory response on the global level, leading to the recovery of proper immune balance on the local level and finally cell death of ectopic endometria. Studies investigating this hypothesis are currently underway.

Our study shows for the first time the presence of FoxP3 + cells in the endometriotic lesion in close proximity to PIF expressing epithelial cells ( Fig 3D ). The presence of Tregs in eutopic and ectopic endometrium of women with endometriosis has been reported previously [ 15 , 16 ]. FoxP3 mRNA is increased in ectopic endothelial tissue and the percentage of Tregs is significantly decreased in the peripheral blood of women with endometriosis, compared to healthy controls. However in the peritoneal fluid the Treg percentage is increased [ 48 ]. This discrepancy suggests a differential immune modulatory system at the local and global levels. Given that sPIF influences ectopic endometrial cells and potentially interacts with FoxP3 + cells (Figs 1 and 3 ), we speculate that PIF expression may mediate an immune privilege for endometriotic lesions. A similar role was previously reported for the Fas-FasLigand system [ 49 ]. Notably, in autoimmune disorders CD4 + /CD25 + /FoxP3 + dependent suppression of effector cells (macrophages, natural killer, dendritic, and cytotoxic T cells) was reported [ 5 , 46 ]. As FoxP3 + cells were present in ectopic and eutopic epithelial cells it is possible that FoxP3 contributes to epithelial transition and differentiation as a transcription factor rather than a reaction to the inflammatory process as in case of TNFα [ 4 ]. Recent evidence suggests apoptotic epithelial cells contribute to Treg survival and abundance [ 50 ]. A role for PIF in mediating this effect would be plausible as cellular expression increases close to implantation and thus is associated with the development of immune tolerance for the embryo [ 22 , 24 ]. The observation that PIF may interact with those cells ( Fig 3D ) combined with their close localization (Figs 1 and 3D ) support a potential regulatory function, as shown in both in vitro and in vivo setting in CD4 + cells. It also raises the possibility that a similar action takes place locally at the site of endometriosis. A number of preclinical models have documented that PIF has both a local and global protective effects [ 20 ] and a potential for a bi-directional communication between PIF expressing epithelial cells and Treg cells exists. PIF increased cell viability of epithelial cells derived from the eutopic endometrium and decreased cell viability of cells derived from ectopic tissue ( Fig 1C and 1D ). These results support a local paracrine effects of PIF on epithelial cells within the endometriotic lesion and this effect may be modulated by the inflammatory microenvironment via TNFα ( Fig 3E and 3F ). TNFα is a hierarchal cytokine that stimulates the expression of numerous other inflammatory mediators though the activation of the IKKβ/NFκB pathway, a pathway specifically activated in the epithelial cells of endometriotic lesions [ 51 ]. The whole genome transcriptome analysis ( Fig 2 ) showed a significant influence on the expression of T cell receptor signalling pathways raising the possibility that not only does the paracrine effects of PIF influence the survival of ectopic epithelial cells, but that it also influences gene expression and the ability of these epithelial cells to respond to the immune regulating effects of infiltrating Tregs. The variation in response to PIF between the eutopic and ectopic derived epithelial cells is also of interest (Figs 1 and 3 ). Although the pathogenesis of endometriosis is still not resolved, the notion that ectopic endometriotic cells have an inherent characteristic that leads to implantation is intriguing [ 3 ]. Such as a pathological alteration may be due to the ectopic environment or may have the origin in stem cells [ 3 ]. Since embryo development is largely dependent on PIF and the fact that sPIF promotes cultured embryo and stem cell development support such a premise [ 17 , 31 , 52 ]. Although hypothetical the idea that a peptide such as PIF may be the missing link between the role of stems cells and immune responses in the pathogenesis of endometriosis is intriguing [ 3 – 5 ]. Further investigation is needed. Lastly, endometriotic derived PIF may also influence the symptomology of endometriosis. sPIF is a well described neuroprotective compound [ 20 ]. Recent evidence suggest that the pain produced by endometriosis occurs through an interaction with endometriotic associated nerve fibers detected in close proximity to the lesions [ 53 , 54 ]. The neuroprotective effects of epithelia derived PIF may stimulate an enhanced nerve presence and influence the interpretation of pain in endometriosis. On the other hand peripheral sPIF and the ensuing reduced neuroinflammation may limit such neurotropic pain as well [ 30 , 45 ].

Synthetic PIF 15 ( MVRIKPGSANKPSDD ) and a scrambled peptide sequence with the same amino acids in random order ( GRVDPSNKSMPKDIA ) were synthesized by solid-phase peptide synthesis (Peptide Synthesizer, Applied Biosystems) employing Fmoc (9-fluorenylmethoxycarbonyl) chemistry at Bio-Synthesis, Inc. (Lewisville, TX, USA). Final purification was carried out by reversed-phase HPLC and the peptide identity verified by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry and amino acid analysis at >95% purity. Fluorescein labeled FITC-PIF and scrambled PIF (FITC- PIFscr) were generated as previously reported [ 22 , 34 ] and the anti-PIF monoclonal antibody against MVRIKPGSANKPSDD was generated in (Genway, SanDiego, CA, USA). Tissue specimens were obtained from 25 women who underwent laparoscopic surgery for severe endometriosis according to the revised criteria of the American Society Reproductive Medicine [ 35 ] (see S1 Table for details). The surgical procedures were carried out in CERM-Hungaria Institute, Rome, Italy, from September 2014 through April 2015. The study was approved by Ethical Committee of CERM-Hungaria (Istitutional Review Board approved 10-02-2014) and conducted according to the Helsinky’s Declaration of Human Rights. Each patient undergoing surgical procedures signed a written informed consent. Samples were obtained from the ectopic endometrium, ovarian endometriomas and peritoneal implants. A total of 25 eutopic endometria, 25 ovarian endometriomas and 10 peritoneal implants were collected from patients. Furthermore, the endometria of 10 healthy women were used as controls. Biopsy samples were fixed in 4% neutral buffered formalin overnight with subsequently paraffin embedding. Before performing immunohistochemistry (IHC) tissue sections were stained with eosin and hematoxylin to select ectopic tissue containing epithelial cells. Serial sections of 5μm thick were used for IHC. Commercially available monoclonal antibodies were used for the detection of PIF (BioIncept, LLC, NJ, USA) and Foxp3 (number: sc-53876, Santa Cruz, CA, USA). IHC was performed according to manufacturer’s instructions. Briefly, tissue sections were dewaxed and re-hydrated and endogenous peroxidase activity quenched by incubation with 0.3% hydrogen peroxide in methanol for 30 minutes at room temperature. Sections were exposed to a non-immune block with normal horse serum for 30 minutes at room temperature. Incubation with the first antibody was carried out at 4° C overnight with a dilution of 1 to 100 for the monoclonal mouse anti-human PIF and with a dilution of 1 to 50 for the monoclonal mouse anti-human Foxp3. Thereafter tissue sections were labelled with the avidin-biotin-peroxidase detection system Vectastain (Vector Laboratories, Burlingthon VT, USA). Each step was followed by washing with PBS. Finally, 3,3'-diaminobenzidine (DAB) and/or 3-amino-9-ethylcarbazole (AEC) were used as chromogens for single or double staining. Counterstaining was performed with hematoxylin. The positive controls were previously confirmed PIF and Foxp3 positive tissue. Negative controls were performed by replacing the primary antibody with mouse immunoglobulin at the same concentration used previously. A semi-quantitative statistical analysis of specific staining was performed using an HSCORE system [ 36 ], calculated using the following equation: HSCORE = SPi(I+1), where i is the intensity of staining with a value of 1, 2, or 3 (weak, strong or very strong respectively) and Pi is the percentage of stained cells at each intensity, varying from 0% to 100%. For all samples, ten microscopic fields were counted by two observers blinded to group status and three slides of each sample were analysed. The HSCORE analysis was performed separately for epithelial and stromal compartment (each observer performed 4 different HSCOREs for each slide). The intra-observer and inter-observer coefficient of variation were 3.4% and 4.2% respectively. We reported on FITC-PIF binding in both pregnant and non-pregnant patients previously [ 21 ]. Briefly, at Millenova Immunology Laboratories non-pregnant, infertile and first trimester pregnant patients undergoing fertility treatment signed a standard informed consent (CARI, Institute, Chicago, USA). All experiments were performed in accordance with the guidelines and regulations of CARI, Institute, Chicago and with the approval from the Institutional Review Board of the University of Illinois, Chicago in March 2006, Dr. R. Roussev, PI. Blood was sampled as part of their work-up process and excess sample used. To test whether PIF targets regulatory T-cells we used specific anti-CD4 + , anti-CD25 + , and anti-FoxP3 + antibodies (BD, Pharmingen, USA). Binding was compared with scrambled-FITC-PIF used as a negative control. We isolated peripheral blood mononuclear cells (PBMC) following separation using Ficoll-Hypaque, Binding to Isotype control served as negative controls. Two- three color staining was performed using conventional techniques. Fluorescence measurements (20,000–50,000 gated events per sample) were performed in a Coulter ® Epics ® XL™ Flow Cytometer using System II software for data acquisition and analysis (Beckman Coulter, Inc., Miami, FL, USA). After ethical approval was obtained from the Bern Cantonal Ethical Committee (KEK) and written informed consent provided by participants endometrial biopsies were collected via soft curette (Pipelle de Cornier, Laboratorie CCD, France) and stored in RNAlater at -80°C from women undergoing laparoscopic surgery at the University Hospital Bern, Switzerland, as described previously [ 37 ]. The pelvic cavity was examined and any endometriotic lesions were removed and the presence of endometriosis confirmed via histological investigation. Endometrial biopsies were collected from women with (n = 4) and without (n = 5) endometriosis. All surgeries were performed during the proliferative phase of the menstrual cycle. Primary endometrial stromal cells (ESC) from women with and without endometriosis were prepared via collagenase digestion and size exclusion protocols as described previously [ 38 ]. Isolated ESC were maintained in Iscoves’s modified Eagle medium (IMEM) (Invitrogen) supplemented with 10% fetal calf serum (fcs) (Invitrogen) and 1% antibiotic/antimycotic (Invitrogen). Immortalized epithelial cell lines were kindly provided by Professor Kyo, Kanazawa, Japan and were isolated from eutopic endometrium, EM E6/E7 [ 39 ] and an ectopic endometriomas, EM’osis, [ 40 ]. Epithelial cells were cultured in Dulbecco’s modified Eagles medium (DMEM) (invitrogen) with 10% fcs and 1% antibiotic/antimycotic. To determine cell viability the CellTiter96 AQueous One Solution Cell Proliferation Assay (MTS) (Promega) was used. Cells were plated into 96 well plates at a density of approximately 6 x 10 3 /well. Sixteen hours prior to treatments cells were changed into a reduced serum media (0.5% FCS). sPIF was prepared by diluting into phosphate buffered saline (PBS) at a final concentration of 300nM. Subsequent 1:3 serial dilutions (111nM, 37nM, 12.3nM and 4.12nM) were prepared for a dose-response assay. Treatment was performed for a total of 48 hours with the treatment media replaced after the first 24 hours. Cell viability was measured via the addition of the tetrazolium compound provided with the kit and colorimetric change measured by spectrophotometer as per manufacturer’s instructions. A control (without sPIF) was included for each experiment and all subsequent values expressed as a percent of control. For assay including TNFα recombinant human TNFα (R&D, Cat No; 210-TA, Minneapolis, MN, USA) was diluted into PBS at 100ng/ml and included in treatment media. For whole transcriptome expression array analysis ectopic endometrioma (EM’osis) cells were seeded into 6 well plates at approximately 2x 10 5 /well and grown until approximately 80% confluent. Sixteen hours prior to treatment the cells were transferred to 0.5% FCS media to synchronize cell cycles. Treatment with 100nM PIF was performed for 48 hours with the treatment media replenished after 24 hours. After treatment period cells were lysed in Qiazol lysis buffer (Qiagen) and RNA isolated using the RNAeasy mini kit (Qiagen), as per the manufacturers instructions. RNA quantity was measured via the Nanodrop 2000 (Witec) and quality via the Bioanalyser 2000 (Agilent). RNA was considered of sufficient quality if RNA integrity number (RIN) was above 9.8. A final concentration for all samples was approximately 200ng/ml. 6 RNA samples were analysed using the Affymetrix platform according to the manual`s instructions (GeneChip® Human Transcriptome Array 2.0 and miRNA Array, Affymetrix). IHC data are reported as a mean ± standard deviation. Statistical analysis was performed by SPSS statistical package (Chicago, IL USA), using the Mann-Whitney Sum Rank test as appropriate. Analysis of cell viability in the in vitro cell models in response to PIF treatment was performed with a two-way analysis of variance (ANOVA) test with a post-hoc Sidak’s multiple comparison test to determine if cell viability was increased by treatment conditions from control. Significance was considered a value for p < 0.05 and analysis performed with Graphpad Prism 6. The raw microarray data was background-corrected, normalized using the RMA method as implemented in the R/Bioconductor package affy [ 41 ]. Probe sets where redefined using the alternative chip definition file [ 42 ]. Differential gene expression was calculated using the moderated t-test as described previously and implemented in the R/Bioconductor package limma [ 43 ]. For Pathway analysis the output of limma was used to perform gene set enrichment analysis (GSEA) using the SetRank method. The key principle of this algorithm is that it discards gene sets that have initially been flagged as significant, if their significance is only due to the overlap with another gene set. It calculates the p-value of a gene set using the ranking of its genes in the ordered list of p-values as calculated by limma. The following databases were searched for significant gene sets: BIOCYC, Gene Ontology, KEGG, Pathway Interaction Database, REACTOME, and WikiPathways.

(DOCX) Click here for additional data file. (XLSX) Click here for additional data file. (XML) Click here for additional data file.