Reproductive Biology and Endocrinology Analysis of Gene Expression Profiles in Hela Cells in Response to Overexpression or Sirna-mediated Depletion of Nasp

Background: NASP (Nuclear Autoantigenic Sperm Protein) is a linker histone chaperone required for normal cell division. Changes in NASP expression significantly affect cell growth and development; loss of gene function results in embryonic lethality. However, the mechanism by which NASP exerts its effects in the cell cycle is not understood. To understand the pathways and networks that may involve NASP function, we evaluated gene expression in HeLa cells in which NASP was either overexpressed or depleted by siRNA.


Background
Nuclear autoantigenic sperm protein (NASP) is highly expressed in all dividing cells including embryonic and malignant tissues as either tNASP (testicular/embryonic isoform) or sNASP (somatic/embryonic isoform). Human tNASP contains three functional histone binding sites; sNASP is the shorter version of tNASP with two deletions in the coding region arising from alternative splicing, resulting in two histone binding sites. In embryonic and transformed cell lines both sNASP and tNASP are expressed and tNASP is present in a variety of malignant tumors [1]. Different cancer types as well as different stages of the same cancer demonstrate specific NASP RNA expression profiles, for example expression levels are up regulated in grade 1 and 2 versus grade 3 in breast cancer [2], estrogen receptor positive versus negative tumor types, or sporadic versus BRCA1/BRCA2 mutation positive tumors [3]. Thus, the NASP expression profile could be used to establish the "poor prognosis signature" which consists of genes regulating cell cycle, invasion, metastasis and angiogenesis [4]. Indeed, NASP has been reported as a serologic marker for ovarian cancer, which could be suitable for clinical testing in high-risk populations [5].
NASP has interactions with a variety of chromatin remodeling proteins: 1) linker H1 histones [1], 2) CAF1(p150, p60, p48) and HIRA in DNA synthesis dependent or independent nucleosome assembly pathways promoted by histone H3.1 and H3.3 complexes [6,7], 3) proteins involved in DNA repair (NASP is phosphorylated after DNA damage by irradiation of U2OS cells [8] and binds to Ku70/Ku80 and DNA PK in HeLa cells [9]). As a linker histone H1 chaperone, NASP binds linker histones in the cytoplasm and transports them into the nucleus [10], where NASP facilitates the incorporation of linker histones onto nucleosome arrays [11]. The present list of NASP interacting proteins is far from complete because Sun et al. [12] reported 356 network connectivity episodes for NASP in developing embryonic stem cells, suggesting numerous direct protein-protein interactions.
NASP is a tightly regulated cell cycle protein because both increased levels of NASP by overexpression [13] and decreased levels induced by siRNA treatment [14] cause disruption of the cell cycle. NASP mRNA levels increase during S-phase and decline during G 2 concomitant with histone mRNA levels [1] and NASP is required for cell survival because the NASP -/null mutation causes embryonic lethality [14]. Because HeLa cells have well studied signaling pathways and are easy to transfect, we chose them to study the general effect of increased or decreased NASP levels. Those networks and pathways that are involved in reproductive physiology were identified and will be further studied in reproductive tissues.

Materials
All chemicals and reagents used in this study were molecular biology grade. Restriction enzymes were purchased from Roche Applied Science (Indianapolis, Indiana, USA). Purification of plasmid DNA and PCR products was carried out using QIAprep Miniprep and QIAquick PCR purification kits (Qiagen, Valencia, California, USA); sequencing was performed at the University of North Carolina at Chapel Hill automated sequencing facility. Goat antiserum to full-length human recombinant tNASP (GenBank: AAH10105) was made by Bethyl Laboratories (Montgomery, Texas, USA).

NASP overexpression
The entire coding sequence of mouse tNASP (nucleotides 92-2405, GenBank: AF034610) was amplified from mouse testis Quick-clone cDNA (Clontech, Palo Alto, California, USA) using the Expand High Fidelity PCR System (Roche Applied Science, Indianapolis, Indiana, USA) and cloned into a Kpn1/BamH1 site in the pEGFP-N1 vector (Clontech, Palo Alto, California, USA) which contains the sequence for the green fluorescent protein. To prevent possible effects of the expressed GFP [15], the GFP sequence was removed. All constructs were sequenced to verify the correct reading frame. Plasmid-DNA complexes were transiently transfected into HeLa cells using Effectene transfection reagent (Qiagen, Valencia, California, USA) according to the manufacturer's recommendations. This method is based on a non-liposomal lipid formulation, and resulted in low cytotoxicity and high transfection efficiency (~97%) as determined by FACS analysis [13]. Control cells were transfected with the transfection reagent only (Effectene transfection reagent, Qiagen, Valencia, California, USA). This choice of control cell treatment was based on the data that "mock" transfection with lipid formulation transfection reagents could consistently affect gene expression [16]. Transfection efficiency was confirmed by Western blotting: lysates from HeLa cells overexpressing tNASP were separated by SDS-PAGE, blotted and probed with goat anti-NASP affinity purified antibody.

siRNA transfection
A series of siRNAs targeting the human NASP open reading frame were designed (Dharmacon, Lafayette, CO) and one (GCACAGUUCAGCAAAUCUAdTdT) was found to effectively deplete both tNASP and sNASP from HeLa cells [14]. Transfection with C2 siRNA, which had no cellular target, served as a negative control [17]. HeLa cells (8.5-10 × 10 5 cells per well in a 24-well plate) were transfected with NASP and C2 siRNA utilizing a two-hit siRNA transfection method with Lipofectamine™2000 for 18 hr as described [17]. Twenty four hours after the first transfection cells were trypsinized and split into 6-well plates.
Forty eight hours after the first transfection cells were retransfected. Ninety hours after the initial transfection cells were harvested for RNA purification.

RNA isolation and hybridization of RNA to oligonucleotide arrays
Total cellular RNA was purified from HeLa cells using RNeasy ® Mini Kit (Qiagen, Valencia, California, USA) according to the manufacturers' instructions. RNA samples representing four separate experiments from cells overexpressing tNASP and four experiments from cells treated with NASP siRNA, along with appropriate controls, were submitted for analysis. After the RNA Quality check was performed the double-stranded cDNA was synthesized from RNA via MMLV reverse transcriptase. Amplified labeled cRNA was created via T7 RNA polymerase, which simultaneously amplifies the target material and incorporates Cy3-or Cy5-labeled CTP with at least a 100-fold RNA amplification rate. cRNA from treated cells was amplified with incorporation of Cy5-CTP (fluorescent in the red region), while cRNA from control samples was labeled by Cy3-CTP (fluorescent in the green region) and purified. cDNA synthesis, cRNA synthesis, amplification and labeling were done using the Low RNA Input Linear Amplification Kit (Agilent Technologies, Wilmington, Delaware, USA). The labeled cRNA samples were then fragmented in fragmentation buffer at 60°C for 30 min before the microarray hybridization. Each sample was hybridized to a whole separate Human Genome (4 × 44K) microarray (Agilent Technologies, Wilmington, Delaware, USA) overnight at 65°C in a hybridization oven. The hybridization slides were washed, stabilized, dried, and immediately scanned by Agilent Technologies Microarray Scanner (Agilent Technologies, Wilmington, Delaware, USA). RNA hybridization was performed in the Genomics and Bioinformatics Core Facility (Lineberger Comprehensive Cancer Center, UNC-CH) according to the protocol suggested by Agilent (Agilent Technologies, Wilmington, Delaware, USA).

Statistical analysis
During the initial analysis at UNC Microarray Database (Genomics and Bioinformatics Core Facility, Lineberger Comprehensive Cancer Center, UNC-CH), all genes were retrieved, appropriately annotated, and filtered. Eventually, only genes with an absolute value of a Log 2 Red/ Green Lowess Normalized Ratio of at least 1 (doubled in intensity) for all 4 arrays were selected. The complete processed and raw data were deposited in Gene Expression Omnibus (GEO) and can be found as GSE14972 [18].
The UNC Microarray Database analysis generated a list of genes with an altered expression (at least 2 fold increased/ decreased) between overexpression/depletion and mock treated samples. To identify which of these genes were significantly differentially expressed (significant genes) we used a statistical technique called SAM (Significance Analysis of Microarrays; [19]). SAM assigns a score to each gene on the basis of a change in gene expression relative to the standard deviation of repeated measurements. For genes with scores greater than an adjustable threshold, SAM uses permutations of the repeated measurements to estimate the percentage of genes identified by chancethe false discovery rate (FDR). Analysis parameters (Delta) were set to result in zero FDR.
To provide a rapid biological interpretation (from "genes to themes") of the obtained data, significant genes from SAM analysis were analyzed by the Expression Analysis Systematic Explorer (EASE). EASE calculates over representation with respect to the total number of genes assayed and annotated within each system to allow comparisons of categories from categorization systems [20]. We used three main categories: biological processes, cellular components, and molecular function. For each classification within the category the Fisher exact probability of over representation was calculated. Presented in this study the EASE score serves as a p-value to the Fisher exact probability that weights significance in favor of themes supported by more genes [20].
Functional interpretation of significant genes in the context of gene ontology, molecular networks and relevance to canonical pathways was generated through the use of Ingenuity Pathways Analysis (IPA 6.5 software, Ingenuity Systems ® [21]). Gene ontology analysis was based on an approach similar to EASE analysis, but the main categories used were in connection to top biological functions: diseases and disorders, molecular and cellular functions, and physiological system development and function. The significant genes were categorized, compared to genetic categories in the IPA database, and ranked according to pvalues. P-values less than 0.05 indicate a statistically significant, non-random association between a set of significant genes and a set of all genes related to a given function in Ingenuity's knowledge base [22]. The IPA analysis determined the subcategories within each category supplied with an appropriate p-value and the number of genes identified.
The set of significant genes was used to find possible connections between genes/gene products and other genes based on interactions previously reported in the literature. Intermolecular connections were presented as molecular networks. Interacting genes (not found in SAM analysis but present in the IPA knowledge base) were used by IPA software to connect smaller groups of significant genes into a larger network. Since the size of the created network could potentially be enormous, the IPA software limited the number of molecules in the network to 35, leaving only the most important ones based on the number of connections for each focus gene (focus genes = a subset of uploaded significant genes having direct interactions with other genes in the database) to other significant genes [23]. Focus molecules and interacting molecules are presented separately. Networks are scored based on the number of focus molecules in the network, its size, the total number of focus molecules analyzed, and the total number of molecules in the knowledge database that could potentially be included in the networks [23]. The network score is the negative log of Fisher's Exact Test pvalue. Only networks with a score of at least 10 (p-value of 10 -10 ) were analyzed.
The final IPA analysis compared the list of significant genes with established pathways associated with metabolism and signaling (canonical pathway analysis). The results are presented in a diagram based on scoring and the ratio of significant genes present in the canonical pathway to the total number of molecules in the canonical pathway. The threshold level was set at p = 0.05. These analyses determined which pathways were involved under our experimental conditions.

Results and discussion
SAM analysis SAM analysis resulted in identification of groups of genes that were significantly differentially expressed after overexpression or depletion of NASP in HeLa cells at least two fold in all four experiments. Overexpression of tNASP in HeLa cells significantly affects the expression level of 54 genes (0.14% of 39,064 genes included in the microarray): 47 (0.12%) were up-regulated and 7 (0.02%) were down-regulated. Figure 1A shows a scatter plot summary of up-regulated and down-regulated genes from cells overexpressing NASP. Inhibition of NASP expression by siRNA treatment significantly affected the expression of 127 genes (0.32%): 56 (0.14%) were up-regulated and 71(0.18%) were down-regulated. Figure 1B shows a scatter plot summary of up-regulated and down-regulated genes from cells treated with NASP siRNA. The different ratios between up-regulated and down-regulated genes demonstrate that there was no bias between overexpression, siRNA treatment, and control groups due to dye affinity misbalance. The list of significant genes is presented in table 1.

EASE analysis
All genes reported in table 1 were subjected to Gene Ontology (GO) clustering by EASE software [20]. EASE analysis is presented in table 2 and table 3. Themes with a p value of < 0.05 and with at least three genes in the category are reported [24].
In HeLa cells overexpressing tNASP, EASE identified 39 up-regulated and 5 down-regulated GO categories (table  2). Fourteen percent of up-regulated genes are located on chromosome 4. In the group of up-regulated genes the highest EASE score as well as the highest number of genes was found within the categories represented by proteins participating in the immune response, the response to biotic and external stimuli, stress, and to pathogens. These results may indicate a rather non-specific reaction of HeLa cells to an excessive amount of expressed recombinant NASP. Down-regulated genes in NASP overexpressing cells represented proteins involved in transcription regulator activity, transcription factor complex structure and activity, and nucleoplasm structure.
Gene ontology categories affected in cells treated by siRNA are presented in table 3. Most of the up-regulated genes in this category are located on chromosome 1(17%), the same chromosome where the NASP gene is located. In cells with inhibited NASP expression, up-regulated categories (n = 13) represented gene products that were involved in morphogenesis and organogenesis, chromosome organization (chromosome, nucleic acid binding, chromatin, DNA packaging, DNA binding), development, and transporter activity. Down-regulated genes in NASP deficient cells represented a group of proteins involved in DNA packaging, nucleosome structure and assembly, establishment and maintenance of chromatin architecture, DNA metabolism, chromosome organization and biogenesis, chromatin assembly/disassembly, and chromatin/chromosome structure (a total of 16 categories). This result is consistent with earlier observations that histone H1 transfer between NASP and DNA affects chromatin structure [13].
Although EASE analysis identified general categories, it did not provide detailed subcategories or connect identified genes to discrete disorders and functions. Therefore additional functional GO analysis was carried out employing the IAP 6.5 software.

NASP overexpression
The highest scoring category in NASP overexpressing cells was "organismal injury", subcategory "fibrosis" (p-value 6.35 × 10 -14 , represented by 13 genes: ATF3, IFI6, IFI27, IFITM1, IL6, IL11, IRF7, IRF9, NT5E, OAS1, PTGS2, STAT1, TNF). A high number of up-regulated genes were functionally related to the "cell growth" subcategory, part of the category "cellular proliferation" (22 molecules out of 47 analyzed, including transcription factors ATF3, DDIT3, FOSL1, RELB, STAT1, cytokines CXCL1, CXCL2, CXCL3, TNF, IL6, IL8, IL11, IL1A, and growth factors EREG, GDF15). The "apoptosis of eukaryotic cell" subcategory of the "cell death" category was represented by an overlapping with the "cell growth" set of genes (total= 24 genes). The most over represented subcategories in the "cell cycle category" were "cell division process of cells "and "arrest in cell division process of cells" with 18 and 12 molecules respectively. The category "cancer" had the highest number of genes (30 out of 47 analyzed). The leading number of molecules (16) was in the subcategory "developmental process of tumor cell lines" and "apoptosis of tumor cell lines". "Immune response", "immunological disease", and "viral functions" categories were selected as a result of the increased expression of a group of cytokines: CXCL1, CXCL3, IL6, IL8, IL11, and TNF.
Down-regulated genes in NASP overexpressing cells were over represented in a group of interrelated functions (figure 2B) mostly as a result of down-regulation of 4 significant genes: transcription regulators EGR1, FOS, LMO1, and NFE2. Their changed expression affected subcategories related to differentiation, maturation, and development within such categories as cellular development (subcategory "developmental process of blood cells", pvalue 1.25 × 10 -4 ), hematological system (subcategory "differentiation of blood cells", p-value 2.67 × 10 -5 ), cell morphology (subcategory "morphology of tumor cells", p-value 3 × 10 -5 ), immune system (subcategory "development of macrophages", p-value 6.69 × 10 -5 ), gene expression (subcategory "transactivation", p-value 2.1 × 10 -3 ). All these categories are interrelated and the different level of significance is a reflection of the total number of genes in each category. The subcategory "arrest in G 0 /G 1 phase of eukaryotic cells" from category "cancer" had the highest significance (p-value of 1.82 × 10 -4 ) within a given category.

NASP inhibition
Although up-regulated genes after siRNA treatment identified a list of discrete and diverse functional categories (figure 2C), the number of molecules presented in each subcategory was low (1)(2)(3)(4). Therefore these categories and subcategories were assigned a low significance. In the "reproductive system disease" category the presence of transcription regulator HMGB2 determined the subcate- Transmembrane receptor TNFRSF11B determined the selection of multiple subcategories associated with bone mineral density (p-value 1.42 × 10 -4 ) within the "connective tissue function" category as well as myogenesis (pvalue 2.89 × 10 -2 ) within the "muscle development" category. Subcategories related to vascular development were determined by combined up-regulation of G-protein coupled receptor F2RL1 and growth factor receptor PDGFRB.

Network analysis NASP overexpression
Network analysis assembled three networks from up-regulated genes after NASP overexpression (table 4A). The network with the highest score (table 4A, network #1) included gene products associated with TNF, which is a multifunctional proinflammatory cytokine that belongs to the tumor necrosis factor (TNF) superfamily. This cytokine is involved in the regulation of a wide spectrum of biological processes including cell proliferation, and differentiation, and has been implicated in cancer and autoimmune diseases [25]. In response to TNF and growth factors, STAT1 protein (found up-regulated in NASP overexpressing cells), is phosphorylated and translocates to the nucleus where it acts as a transcription activator [26], mediating the expression of a variety of genes. In this network ATF3 (activating transcription factor) is connected with TNF, which is known to increase the expression of human ATF3 mRNA in LoVo cells [27].
The other two networks (table 4A, network #2 and #3) contained groups of interacting proteins that included GADD45A, interleukins IL6, IL1A, chemokines CXCL3, CXCL1, and NF-κB. NF-κB is a transcription factor that mediates the transcription of proteins involved in cell survival, proliferation, and inflammatory responses, and is the subject of active research for anti-cancer therapy [28]. The protein encoded by the GADD45A gene induces apoptosis and cell cycle arrest by maintaining p38 and c-JNK MAPK activation in keratinocytes. The absence of Gadd45a results in loss of sustained p38/JNK MAPK activity that leads to inadequate p53 activation and loss of normal activation of G 1 and G 2 checkpoints [29].

NASP inhibition
Distinct networks were assembled with up-regulated genes after NASP siRNA treatment (table 5). The highest scoring network (table 5, network #1) presents E2F8 interacting with TP53. E2F8, which is in the family of E2F transcription factors, is essential for orchestrating expression of genes required for cell cycle progression and proliferation [31]. In the assembled network E2F8 interacts with tumor protein p53 [32], which wasn't found as a focus gene in this study, but was added during network assemblage as an interacting gene. Protein p53 regulates target genes that induce cell cycle arrest, apoptosis, senescence, DNA repair and is postulated to function as a tumor suppressor [33]. Multiple focus genes TMSL8, OGG1, HMGB2, GSP, ARL5A, which were reported to be acted upon by p53 [34][35][36][37][38][39], were found to be up-regulated in this study. One of them, NAP1 (nucleosome assembly protein1) plays a role in chromatin maintenance by facilitating core histone exchange (by regulating the concentration of free histones) as well as nucleosome assembly and disassembly [40].
Network #2 ( Follistatin has been reported as a modulator of gonadal tumor progression in inhibin deficient mice [42]. Network #3 (table 5) of up-regulated genes after siRNA treatment identified TGFB1 as a core molecule indirectly affiliated with interacting proteins. TGFB is a multifunctional peptide that controls proliferation, differentiation, inducing transformation, and other functions in many cell types [43].
After siRNA treatment 71 down-regulated genes were almost evenly assembled into four networks of interacting proteins (table 6). The network with the highest score (table 6, network #1) included: HOXA10, FGF2, SRF, ITGA5, and ITGB3. HOXA10 is a DNA-binding transcription factor that may regulate gene expression, morphogenesis, and differentiation. More specifically, it may function in fertility and embryo viability [44]. FGF2 (fibroblast growth factor) and SRF (c-fos serum response element-binding transcription factor) participate in cell cycle regulation, apoptosis, cell growth and differentiation [45]. Interaction of SRF with other proteins, such as steroid hormone receptors, may contribute to regulation of muscle growth [46]. Network #2 of down regulated genes (table 6) includes TNF as a core protein along with several interacting proteins. Overlap of this network with another assembled network (table 6, network #4) revealed TRAF6 (TNF receptor associated factor), which mediates signal transduction from members of the TNF receptor superfamily  figure 3) include HIST1H1C [1] and HIST1H2BK [49], which implies that these histone gene products may be regulated in concert with NASP. None of the focus genes were found to be a hub in network #3 (table 6), therefore it was not investigated further.

Canonical pathway analysis NASP overexpression
We found that each treatment resulted in activation/inhibition of specific signaling pathways. A list of the significant genes present in canonical pathways after different treatments is presented in Additional file 5, Additional file 6, Additional file 7 and Additional file 8.
The highest activation level after NASP overexpression (figure 4A) was detected in pathways related to antiviral responses, activation of IRF (interferon activation factor), recognition of bacteria and viruses, interferon signaling and other interrelated pathways. Some pathways were inhibited by NASP overexpression ( figure 4B). The significance of all top 6 canonical pathways in this category barely exceeded the threshold level and all of them were selected as a result of down-regulation of only one gene, FOS, which is an important transcription regulator of cell proliferation, differentiation, and transformation [50]. The inhibition of a single gene makes understanding its role in any of those pathways difficult.

NASP inhibition
A similar result was observed with analysis of activated signaling pathways after siRNA treatment ( figure 4C). Only one pathway reached significance above the threshold level, the remaining five canonical pathways were below the threshold level. After siRNA treatment ( figure  4D) several signaling pathways were inhibited due to down-regulation of growth factor FGF2, transmembrane receptors ITGB3 and ITGA5 and other molecules such as ATCA1 and ATCA2. The two highest scores had canonical pathways related to caveolar-mediated (p-value 1.85 × 10 -4 , 4 out of 81 molecules were associated with a given pathway) and clathrin-mediated (p-value 2.39 × 10 -4 , 5/165) endocytosis. Selection of these pathways was based on down-regulation of both αand βintegrins (ITGA5 and ITGB3) which serve as transmembrane receptors. Activation of αand β-integrins by appropriate ECM proteins causes the activation of ERK/MAPK-mediated transcrip-tion, which eventually leads to cell proliferation [51]. Mitogen activated protein kinase MAP2K1 (synonym MEK1/2), an important link in both those pathways, was found to be down regulated and it makes NASP's effect on this pathway plausible. Down regulation of ACTA1 and ACTA2 (actin α-1 and actin α-2) along with MAP2K1, αintegrin, and FGF2 (fibroblast growth factor 2) indicates a high probability that there is an effect on the actin cytoskeleton signaling pathway (p-value 8.17 × 10 -4 , 5/ 222), which leads to actin reorganization and plays an important role in cell motility, cytokinesis and phagocytosis [52]. Decreased expression of kinase MAP2K1 and transcription factor SRF implies an effect on EGF (epidermal growth factor) signaling pathways (p-value 1.25 × 10 -2 , 2/47), which regulate cell growth and differentiation [53].

Conclusion
This study has demonstrated that NASP, a linker histone chaperone, belongs to a network of genes and gene functions that are critical for cell survival. We have confirmed the previously reported interactions between NASP and HSP90, HSP70, histone H1, histone H3, and TRAF6. Indeed, based on the number of critical pathways affected by the overexpression or inhibition of NASP expression, NASP may play a much wider role in gene expression events that require the participation of histones. Significantly this study found that during overexpression the network with the highest score of up-regulated genes included gene products associated with TNF and during inhibition network #2 of down regulated genes also included TNF. TNF receptor associated factor, TRAF6, has already been identified as interacting with NASP. Therefore these two identified networks may explain how the expression of NASP is modulated during cell proliferation and differentiation.
Of interest to reproductive biology, this study found a significant association of NASP with the FSH and follistatin gene pathways that are up-regulated after NASP siRNA treatment. This may have implications for the control of NASP expression during granulosa cell and endometrial cell proliferation.
Most of the pathways found in this study have not been previously reported to be connected to expression levels of NASP and demonstrate a rather complicated picture of changed gene expression in HeLa cells after NASP expression was increased or decreased. Gene ontology and protein network analysis identified general biological processes as well as individual genes/gene products and possible interaction networks. Some of these processes may relate to HeLa-type cells in tissue culture in which cell adhesion and migration are critical, while others may only be relevant in neoplasia. We found signaling path-Top canonical pathways affected by Figure 4 Top canonical pathways affected by. A. Up-regulated genes after NASP overexpression. B. Down-regulated genes after NASP overexpression. C. Up-regulated genes after NASP depletion. D. Down-regulated genes after NASP depletion. Bars represent -log (p-value) for disproportionate representation of affected genes in the selected pathway, yellow line represents the ratio of affected genes to the total number of genes in a pathway. Threshold (red line) denotes the p = 0.05 level.
ways which were affected as a result of changed NASP expression and despite some overlap each reactive response was associated with a unique gene signature. The results of this study have elucidated the changes that emerge from increased and decreased NASP expression and will help our understanding of the molecular mechanisms involved in NASP function. Confirmation of NASP's role in regulating the cell cycle may contribute to the development of new pharmaceutical approaches to control the relevant pathological conditions.