Abstract
It has been previously shown that anti-dengue virus (DENV) nonstructural protein NS1 antibodies could act as autoantibodies that direct against one or more of the host's own proteins, which has potential implications for dengue hemorrhagic fever pathogenesis. In the present study, we have employed suppression subtractive hybridization (SSH) to identify the differentially expressed genes from human microvascular endothelial cells (HMEC-1) in response to anti-dengue virus type 2 NS1 antibodies (anti-DENV2 NS1 Abs). A total of 35 clones from the SSH cDNA library were randomly selected for further analysis using bioinformatics tools after vector screening. After searching for sequence homology in NCBI GenBank database with BLASTN and BLASTX programs, 23 obtained sequences with significant matches (E-values <1×10−4) in the SSH library. The predicted genes in the subtracted library include immune response molecules (CD59 antigen preproprotein preproprotein, MURR1), signal transduction molecules (Nuclear casein kinase and cyclin-dependent kinase substrate 1), calcium-binding proteins (S100A6, Annexin A2 isoform 1/2), and cell-membrane component (Yip1 domain family). From these clones, 5 upregulated genes were selected for differential expression profiling by real-time RT-PCR to confirm their upregulated status. The results confirmed their differential upregulation, and thus verified the success of SSHs and the likely involvement of these genes in dengue pathogenesis.
Introduction
DF is an acute, self-limited febrile illness, whereas DHF and DSS are life-threatening and characterized by plasma leakage, low platelet counts, liver damage, and, in the most severe cases, death due to shock (21). Plasma leakage, also called the hemorrhagic manifestations, occurs once vessel endothelium is disrupted and is followed by the loss of its barrier function. Although dengue virus-infected endothelial cells in vitro had been reported (6,24), there is no direct evidence indicating DENV infection of endothelial cells in vivo (21). It is likely that different mechanisms are involved in the disease pathogenesis, including complement activation (4,32), lymphocytes and cytokines Tsunami (3,8), or cross-reactive antibodies to endothelial cells (2,12).
Antibodies generated by mice to the dengue virus nonstructural protein 1 (NS1) have been shown to cross-react with human fibrinogen, platelets, and endothelial cells (12,13). Serum samples from DHF/DSS patients show higher endothelium binding activity than do those of people with dengue fever (30). Anti-DENV NS1 antibodies act as autoantibodies that cross-react with host proteins and noninfected endothelial cells and trigger intracellular signaling leading endothelial dysfunction (2,29,31). However, the mechanism of vascular permeability and hemostatic disorders as a result of endothelium dysfunction caused by the autoimmune pathogenesis is still not fully understood.
In order to gain insight into the intracellular signaling triggered by the anti-endothelial cell autoantibodies, we employed suppression subtractive hybridization (SSH) to construct cDNA libraries for identifying differentially expressed genes from human microvascular endothelial cells (HMEC-1) in response to anti-dengue virus type 2 NS1 antibodies immunoglobulin G (IgG) (anti-DENV2 NS1 Abs). From genes in this SSH libraries, five upregulated genes were selected for differential expression profiling by real-time RT-PCR to confirm their upregulated status. The information provided from this study will enhance the understanding of anti-NS1 antibodies' role in dengue hemorrhagic fever pathogenesis.
Materials and Methods
Preparation of recombinant NS1 and generation of anti-NS1 Abs
The full-length cDNA of NS1 from DENV-2 virus New Guinea C strain (NGC) was cloned to expression vector pPICZαB (Invitrogen) to establish a pPICZαB-NS1 plasmid; the recombinant NS1 protein expression was induced by adding methanol as described before (49). Recombinant NS1 protein immunogen was prepared and antibodies against NS1 protein were generated by the immunization of New Zealand rabbits. Anti-NS1 polyclonal antibody was purified by caprylic acid-ammonium sulfate precipitation and protein A resin (GenScript, USA). Then the specificity, sensitivity, and valence of the antibody were detected by Western blot and ELISA. The antibodies were quantified by UV280, as we previously described (27).
Cell culture
Human microvascular endothelial cell line-1 (HMEC-1) was obtained from the Centers for Disease Control, Atlanta, GA (1) and passed in culture plates containing endothelial cell growth medium (Gibco, MCDB131, Cat. 10372-019) composed of 10% FBS (Gibco), 1 μg/mL hydrocortisone (Sigma-Aldrich), 10 ng/mL epidermal growth factor (Becton-Dickinson), 10 mM L-Glutamine (Gibco), and antibiotics. Cells were detached using 0.25% trypsin-EDTA (Gibco).
Anti-NS1 Abs cell binding assay
HMEC-1 cells were grown in culture flasks and incubated with 5 mL anti-DENV2 NS1 Abs (8 μg/mL) or control rabbit IgG at 37°C for 1 h. After washing three times with PBS, the cells were incubated with FITC-conjugated anti-rabbit IgG at 4°C for 40 min and analyzed by indirect immunofluresence. In addition, more than 1×106 HMEC-1 cells were collected for binding assay using flow cytometry (Coulter EPICS Elite Flow Cytometer Sorter).
RNA extraction and cDNA synthesis
RNA was extracted from around 2×106 HMEC-1 cells at eight time points (2, 4, 6, 8, 10, 12, 24, and 48 h) post incubation with anti-NS1 Abs using RNeasy Mini Kit (Qiagen). Normal rabbit IgG incubated cells were used as control. After the extraction, total RNA was visualized on a 1.0% agarose Tris–acetate–ethylenediaminetetraacetic acid (TAE) gel. Total RNA in each sample was quantified using a BioSpec-mini spectrophotometer (Shimadzu Biotech, Nakagyo-ku, Kyoto, Japan). For construction of the cDNA libraries for SSH, each sample was created by combining 2 μg RNA from each time point. The cDNA synthesis was performed using 3 μg of RNA from each sample using the SMART™ PCR cDNA Synthesis Kit (Clontech, CA), in accordance with the manufacturer's protocols.
Construction of suppression subtractive hybridization library and differential screening
A suppression subtractive hybridization (SSH) library was prepared with the PCR-Select™ cDNA subtraction kit (Clontech), according to the manufacturer's instructions. Briefly, the forward-subtracted library was constructed where the cDNAs generated from HMEC-1 cells cultured with anti-NS1 Abs served as the “tester”, while cDNAs generated from HMEC-1 cells cultured with control rabbit antibodies served as the “driver”, as well as the reverse-subtracted library where the samples exchanged their positions of tester and driver. The forward-subtracted and the reverse-subtracted cDNA library were cloned using the pGEM-T Easy Kit (Promega, WI, USA). Then, 5 μL of the ligation reaction was used to transform Trans10 competent cells (TransGen, Guangzhou, China) for subsequent ampicillin-resistant colony selection and blue-white screening. To construct probes for differential screening, the forward-subtraction and reverse-subtraction cDNAs as well as the unsubtracted cDNA libraries were labeled with digoxigenin, and Southern blot was performed for colony PCR positive clones to exclude the false positive ones using DIG High Prime DNA Labeling and Detection Starter Kit I, according to the manufacturer's instructions. Results analysis was done according to the manual of the PCR-Select™ differential screening kit (Clontech).
Sequencing and homology search
The differentially expressed clones were sequenced using an automated DNA sequencer (Invitrogen, Shanghai, China) and analyzed using the online BLAST program from NCBI's website (
Quantitative RT-PCR
For quantitative RT-PCR, the specific oligonucleotide primers to five selected genes including ADP-ribosylation factor-like 6 interacting protein “ARL6IP” (forward: 5′-ACCACAAGGTACAACATCAGTGCAATAA-3′ and reverse: 5′-TCTATGACTTTCACTCTAAGACTGGCAA-3′), aryl hydrocarbon receptor nuclear translocator isoform 1/3 “ARNT1/3” (forward: 5′-AATGTGCCTTATGTTTGTATGTCTGGAG-3′ and reverse: 5′-GTTACTACCCTCACCCTTTACCTTCCTC-3′), CD59 (forward:5′-GGCTCCAAAGCACTATGTCACCTACCTC-3′ and reverse: 5′-AATCAGTCCTGTCCTTGGCAAGCTCTTC-3′), nuclear ubiquitous casein kinase and cyclin-dependent kinase substrate 1 “NUCKS1” (forward: 5′-AACTTTGCTTTCCACAACACTTCATAAC-3′ and reverse: 5′-TTTATTGAAGGGTTTCCAGGACTGCGTG-3′), S100A6 (forward: 5′-TACAATTACCCACCACTGGATTTGACTC-3′ and reverse: 5′-ACCAGGAGGTGAACTTCCAGGAGTATGT-3′) and β-actin (forward: 5′-GAACCTCTCGTTGCCGATGGTG-3′ and reverse: 5′-GAAGCTGTGCTACGTGGCTCTG-3′) were used in the experiment. Two samples of RNA were prepared as mentioned before. Quantitative RT-PCR was performed to detect differences in transcript levels of the selected ESTs at different time points. Amplification was determined using SYBR Premix Ex Taq (Takara, Guangzhou, China) in the CFX96™ Real-Time PCR Detection System (Bio-Rad, CA, USA). Each real-time PCR reaction was done in triplicate. The results were presented as the relative expression ratios of the target genes, which were standardized to the expression level of the constitutive β-actin gene. The data were analyzed using SPSS statistics 17.0 (Chicago, USA). Statistical significance was calculated by paired-samples t-test and one-way ANOVA. Significance was accepted at p<0.05.
Results
The antibodies against DENV2 NGC strain NS1 protein were obtained by the immunization of New Zealand rabbits with the recombinant expressed NS1 protein. The antibody could specifically recognize native DENV2 NS1 protein analyzed by Western blot and the titer of the antiserum was up to 1:6000 analyzed by ELISA (27). The binding assay showed that the anti-DENV2 NS1 antibodies could cross-react with HMEC-1 cells by using immunofluorescence assay (data not shown) and flow cytometry (Fig. 1).

Binding of anti-DENV2 NS1 antibodies and HMEC-1 cells analyzed using flow cytometry. HMEC-1 cells were incubated with anti-DENV2 NS1 antibodies (black), rabbit IgG (gray), or no antibodies (solid). The results is shown as a representative of triplicate experiments.
A subtracted cDNA library enriched in differentially expressed genes in HMEC-1 response of anti-DENV2 NS1 Abs was generated using SSH. The subtraction efficiency was evaluated by comparing the transcript abundance of G3PDH before and after subtraction. The result showed that G3PDH was not detectable in the subtracted cDNA samples, whilst the presence of G3PDH was discovered in the nonsubtracted cDNA samples whenever the subtraction was forward or reverse (data not shown), implying that the SSHs were successful.
Sequence analysis of subtracted cDNA library
A total of 35 clones from the SSH cDNA library were randomly selected and sequenced using the M13 reverse primer. After searching for sequence homology in NCBI GenBank database with BLASTN and BLASTX programs, 23 obtained sequences with significant matches (E-values <1×10−4) in the SSH library. Known genes were categorized according to their putative functions (Table 1). The predicted genes in the subtracted library include immune response molecules (CD59 antigen preproprotein preproprotein, MURR1), signal transduction molecules (Nuclear casein kinase and cyclin-dependent kinase substrate 1), calcium-binding proteins (S100A6, Annexin A2 isoform 1/2), and cell-membrane component (Yip1 domain family).
Determination of differentially expressed cDNAs by real-time qPCR
To confirm the differential expression of genes identified from the SSH library and to understand gene expression variations coupled with HMEC-1 response of anti-DENV2 NS1 Abs, nuclear ubiquitous casein kinase and cyclin-dependent kinase substrate 1 (NUCKS1), CD59 preproprotein, S100A6, aryl hydrocarbon receptor nuclear translocator isoform 1/3 (ARNT1/3), and ADP-ribosylation factor-like 6 interacting protein (ARL6IP) were selected for real-time PCR based evaluation of their relative expression levels (Fig. 2). The expression of NUCKS1 was upregulated from 6 h up to 48 h (except for at 24 h) post anti-DENV2 NS1 Abs incubation. The transcripts of NUCKS1 gene showed much higher in the anti-DENV2 NS1 Abs-treated group than control Abs-treated group at 6 h, 8 h, and 10 h, but their transcripts levels were almost the same at 12 h and 48 h. Strangely, the expression levels of NUCKS1 were much lower in the anti-DENV2 NS1 Abs-treated group than control Abs-treated group at 24 h. The expression of CD59 was also upregulated after 6 h post anti-DENV2 NS1 Abs incubation, except for a relatively lower expression level at 8 h and 24 h, respectively. The transcript levels were considerably higher in the anti-DENV2 NS1 Abs-treated group than control Abs-treated group at 6 h and 10 h. However, the expression of S100A6 increased only at 6 h and 24 h when compared with control Abs treated group. The expression pattern of ARNT was almost the same as NUCKS1 and CD59: the transcript levels increased after 6 h of anti-DENV2 NS1 Abs incubation except for a relatively lower expression level at 8 h and 24 h, respectively. The transcript levels were consistently higher in the anti-DENV2 NS1 Abs-treated group than control Abs-treated group at 6 h and 10 h, while they were almost the same at 12 h and 48 h. The expression of ARL6IP was upregulated after 10 h up to 48 h post anti-DENV2 NS1 Abs incubation except 12 h. The transcript levels were consistently higher in the anti-DENV2 NS1 Abs-treated group than control Abs-treated group at 10 h and 24 h.

Relative fold expression of five different target genes, nuclear ubiquitous casein kinase and cyclin-dependent kinase substrate 1 (NUCKS1), CD59, S100A6, aryl hydrocarbon receptor nuclear translocator isoform 1/3 (ARNT1/3), and ADP-ribosylation factor-like 6 interacting protein (ARL6IP), expressed in the HMEC-1 cells incubated with anti-NS1 Abs (solid bar) compared to normal rabbit IgG incubated cells (open bar) and standardized against β-actin as the internal reference, were analyzed by real-time RT-PCR at 2, 4, 6, 8, 10, 12, 24, and 48 hours post incubation. Samples from each time point were tested in triplicate and the data represent the mean (±1 standard deviation) relative expression level was used for analysis that the target gene is downregulated, the same, or upregulated, respectively, compared with the control (*p<0.05; **p<0.01).
Discussion
Plasma leakage and hemorrhagic are the hallmarks of dengue hemorrhagic fever. Although studies showed Abs directed against DENV NS1 would cross-react with endothelial cells and cause damage, the pathogenic mechanism, however, has not been well understood (12,30). To develop a better understanding of the intracellular signaling triggered by the anti-endothelial cell autoantibodies, SSH is a powerful approach adopted herein to identify endothelial cell-responsive genes upon anti-DENV NS1 Abs incubation. Suppression subtractive hybridization (SSH), which is a technology that combines normalization and subtraction in a single procedure based primarily on suppression polymerase chain reaction (PCR), allows the isolation of differentially expressed cDNAs with specific genomic sequences (11). SSH could be employed to construct differentially expressed gene libraries between anti-DENV NS1 Abs treated and normal endothelial cells. From our SSH library, we found certain genes involved in various cellular functions that related to immune response regulation, cell growth, and transcriptional regulation.
Nuclear ubiquitous casein and cyclin-dependent kinases substrate (NUCKS) is a nuclear DNA binding protein occurring in almost all types of human cells, adult and fetal tissues (18,41). It has all the features of being a housekeeping gene (18). The abundance of NUCKS was found in rapidly growing cells, as well as in certain cancers, suggests that it may be involved in facilitating and maintaining activity of transcription of some genes during rapid proliferation and in cancer (39). The importance of NUCKS remains still unknown and disputable. Overexpression of NUCKS is probably related to high levels of transcription. Future studies should also pay additional attention in elucidating potential functions of specific post-translational modifications of NUCKS in etiology and progression of this disease.
Complement plays a well-known important role as a multi-functional complex system comprising more than 30 proteins, in host defense against infectious agents, in the removal of apoptotic cells and immune complexes, and in the modulation of the adaptive immune system (36). CD59 is the key regulator of the terminal pathway which inhibits the formation of membrane attack complex (MAC), that limits host cell damage by binding to C5b-8 or C5b-9 and preventing leakage or lysis of the host cell membrane (23,33). Cells with normal levels of the complement regulatory proteins DAF, MCP, and CD59 defend themselves against immediate lysis, but the attack persists as long as more antigen is expressed, more antibody is available, and complement levels remain sufficient to maintain the attack. Further study is needed to reveal the potential role of CD59 in response against leakage of endothelial cells triggered by anti-DENV NS1 Abs.
The family of S100 proteins is one of the largest subfamilies of EF-hand proteins. The S100A6 (calcyclin) protein, belongs to the S100 protein family, is a 10.5 kDa Ca2+-binding protein (38,47). The function of S100A6 is not clear at present, but it has been suggested that it may be involved in cell proliferation, cytoskeletal dynamics, and tumorigenesis. S100A6 expression at the protein and/or mRNA level can be upregulated by multiple factors, including platelet-derived growth factor (PDGF), epidermal growth factor (EGF), serum, TNFα, retinoic acid, estrogen, palmitate, vasopressin, and gastrin (5,10,16,25,28,42,43). The interaction of S100A6 and annexins, tropomyosin, caldesmon, and calponin suggest that S100A6 may play a regulatory role in the actin cytoskeleton (14,16,17,46,48). In this research, we found that S100A6 was upregulated following anti-DENV NS1 Abs treating; as a result, we suggest that S100A6 might be implicated in the increasing of endothelial cells permeability response to the antibody, thus it is attractive to examine the function of S100A6 for being a signal molecules and its interaction with the corresponding proteins.
The AhR nuclear translocator (ARNT; also known as HIF-1β) is the dimerization partner of the aryl hydrocarbon receptor (AhR) in the nucleus. The AhR-ARNT heterodimer has been shown to regulate the expression of several batteries of genes involved in diverse signaling pathways, including phase I/II metabolism, inflammation, and cell cycling (7,15). The AhR is a member of the PAS (Per-ARNT-Sim) superfamily of proteins. Physiologically, many of these proteins act by sensing molecules and stimuli from the cellular environment, and initiating signaling cascades to elicit the appropriate cellular responses. Some experimental results suggest that AhR plays an important role during development and maintenance of the hematopoietic and immune system (22,35,40,44). Our SSH library also revealed an upregulated state of this gene, so it merits further functional investigation of ARNT in endothelial cells response.
ADP-ribosylation factor-like 6 interacting protein 1 (ARL6IP1) is an apoptotic regulator, as well as involved in protein transport, membrane trafficking, or cell signaling during hematopoietic maturation (37,45). Mouse ARL6IP1 was confirmed to physically interact with mouse ARL6, which belongs to ARF/ARL family in the Ras superfamily (9,26). However, the biological functions of ARL6IP1 have not been fully elucidated. It raises the interesting question for further study of whether ARL6IP1 might be also implicated in the endothelial cells response to anti-DENV NS1 Abs.
In conclusion, subtractive hybridization libraries of HMEC-1 response of anti-DENV2 NS1 Abs were successfully constructed and some differentially expressed genes were identified. The information indicated there were many genes involved in the pathogenic mechanism directly or indirectly through signaling molecules. This baseline study can serve as a starting point for further work on the cell-autoantibody cross-reaction mechanism and increase our understanding of dengue pathogenesis.
Acknowledgments
This study was supported by Fund of Technology Program of Guangdong Province of P.R. China (2011B031800360).
Author Disclosure Statement
No competing financial interests exist.
