{'fraction_non_alphanumeric': 0.04099883855981417, 'fraction_numerical': 0.02257065427797135, 'mean_word_length': 4.935431985294118, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Analysis of Dermal Papilla Cell Interactome Using STRING Database to Profile the ex Vivo Hair Growth Inhibition Effect of a Vinca Alkaloid Drug, Colchicine

Dermal papillae (DPs) control the formation of hair shafts. In clinical settings, colchicine (CLC) induces patients’ hair shedding. Compared to the control, the ex vivo hair fiber elongation of organ cultured vibrissa hair follicles (HFs) declined significantly after seven days of CLC treatment. The cultured DP cells (DPCs) were used as the experimental model to study the influence of CLC on the protein dynamics of DPs. CLC could alter the morphology and down-regulate the expression of alkaline phosphatase (ALP), the marker of DPC activity, and induce IκBα phosphorylation of DPCs. The proteomic results showed that CLC modulated the expression patterns (fold > 2) of 24 identified proteins, seven down-regulated and 17 up-regulated. Most of these proteins were presumably associated with protein turnover, metabolism, structure and signal transduction. Protein-protein interactions (PPI) among these proteins, established by Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database, revealed that they participate in protein metabolic process, translation, and energy production. Furthermore, ubiquitin C (UbC) was predicted to be the controlling hub, suggesting the involvement of ubiquitin-proteasome system in modulating the pathogenic effect of CLC on DPC.

Interactome (the network of interacting proteins) consists of the regulatory pathways to determine the cellular behavior. An advantage of PPI maps is that researchers can gain understanding about specific signaling pathways and the functions of individual proteins [21].
In our previous report, we demonstrated that CLC not only inhibited the proliferation of DPCs, but also stimulated DPCs to release histone H4 to retro-modulate themselves [22]. Due to the niche effect of DP on HF [11], we further examined the effect of CLC on HFs ex vivo and DPCs in vitro. Using a combined approach of 2-DE and nano-LC-ESI-MS/MS analysis, the protein patterns of DPCs in the presence and absence of CLC were profiled. The differentially expressed proteins were further analyzed using STRING database to derive a genome-scale PPI network and the potential signaling pathways, which may reveal the pharmacological mechanism of CLC on DPCs.

The Effect of Colchicine (CLC) on Hair Growth
To examine the effect of CLC on the hair shaft elongation, organ cultured Sprague-Dawley (SD) rat vibrissa HFs were exposed to CLC of different concentrations (0.1 and 1 µM) for seven days. As shown in Figure 1, CLC significantly inhibited the hair shaft elongation in a dose-dependent manner from 12.5% ± 2.3% by 0.1 µM to 58% ± 4.5% by 1 µM after seven days. The organ cultured rat vibrissa hair follicles (HFs) were treated with CLC (0.1 and 1 µM) for seven days and the length of HF fiber was determined. The vibrissa HFs with DMSO treatment was taken as control (100%). Data were expressed as the mean from three batch triplicate experiments; bars, ±SD. * p < 0.05, ** p < 0.01 vs. control group.

CLC Induced Morphological Changes of DPCs
After exposing DPCs to 2 µM CLC for 6 h, the morphologies of DPCs were altered significantly, from elongated fibroblast-like to discoid or polygonal epithelial-like cells (Figure 2A) without obvious DPC apoptosis as determined by trypan blue exclusion. Interfere with normal cytoskeletal dynamics by the addition of CLC brings great stresses on cells and causes the destabilization of polymerized microtubule and the loss of cell shape [23]. Accordingly, 2 µM CLC disrupted the microtubules of DPCs to cause the morphological changes of DPC even within 6 h. Our previous data have shown that the long-term (72 and 120 h) proliferation of DPCs was inhibited under the same treatment [22]. . In all WB analyses, the loading and transfer of equal amounts of total lysate (20 µg per lane) were confirmed by immuno-detection of Glyceraldehyde 3-phosphate dehydrogenase (GAPDH); (C) The expression level of ALP mRNA in DPCs was analyzed by qRT-PCR. The level of ALP mRNA was normalized to that of GAPDH. Data were expressed as the mean from triplicate experiments; bars, ±SD, ** p < 0.01 vs. control group; (D) the time course of IκBα phosphorylation of DPCs after CLC treatment at indicated times was analyzed by WB.

CLC Down-Regulated the Expression of ALP in DPCs
Since ALP activity is the widely used cell marker of DPCs [8], the ALP level of DPCs in the presence or absence of 2 µM CLC was determined by Western blot (WB) and real-time quantitative reverse transcription-PCR (qRT-PCR) analysis. CLC (2 µM) significantly down-regulated the expression of ALP in both protein and transcript levels of DPCs after 6 h ( Figure 2B,C) The Implication of ALP Down-Regulation on the Biological Function of DPCs DPs, which mainly consist of DPCs, play a key role in directing matrix epithelial cells to produce hair shafts in anagen [3,11]. CLC has been shown to induce animal and patients alopecia [17][18][19][20]. The ALP expression of DPCs is recognized both as the marker of DPC activity and an index for the degree of drug-induced damage to DPCs within HF dystrophy [10]. Cultured DPCs with strong ALP activity interact with epithelial cells to induce the neogenesis of HFs in vivo [9]. As shown in Figure 1, CLC inhibited the ex vivo elongation of rat vibrissa hair fiber in a dose-dependent manner. Although compared to non-proliferative DPCs, the active proliferating matrix cells are more sensitive to CLC, we deduce that the impact of CLC to the activity (ALP) of DPCs also participated in the inhibition of hair shaft formation resulted from CLC. Collectively, the CLC-treated DPCs with reduced ALP activity were suggested to be inactivated and have lost the ability to direct epithelial cells to produce hair shafts and new HFs.

CLC Induced IκBα Phosphorylation of DPCs
The time course of protein levels of phosphor-IκBα of DPC were monitored by WB following 2 µM CLC treatment ( Figure 2D). The phosphorylation levels of IκBα were increased in 5 min and the increment was about 3.5-fold at 15 min. Unto 30 min of treatment; levels of phosphor-IκBα were reduced to near undetectable level. The phosphorylation of IκBα promoted its ubiquitination and degradation in proteosome that led to release of NF-κB from the inhibitory IκBα-NF-κB complex. The released NF-κB then translocated into nucleus to regulate downstream genes [24].

Two-Dimensional (DE) Analysis for the Effect of CLC on DPCs
The representative 2-DE maps of DPCs treated with or without 2 µM CLC for 6 h were depicted and each gel resolved up to ~654 protein spots ( Figure 3A,B). The image analysis indicated that 26 protein spots changed their quantities greater than twofold after treating DPCs with CLC. Those protein spots were excised, in-gel trypsin digested and analyzed by the nano-LC-ESI-MS/MS to obtain the MS/MS peak lists. Twenty-six protein spots, nine down-regulated and 17 up-regulated, were identified to match to candidate proteins using Mascot online search engine (Tables 1 and 2). Both horizontal adjacent down-regulated protein spots 3 and 4 were identified to be protein disulfide-isomerase A6 (Pdia6). The other two horizontal adjacent down-regulated protein spots 7 and 8 were also both identified to be Annexin A5 (Anxa5). Therefore, seven proteins were proteomic identified to be down-regulated after treating DPCs with 2 µM CLC for 6 h. The pattern of protein changes of CLC-treated DPC, except Psmb7, is quite different from those of taxol-treated DPCs [29].  Table 1 lists the identities of the seven down-regulated proteins in DPCs treated by 2 µM CLC for 6 h. Among seven down-regulated proteins, heat shock cognate 71 kDa protein (Hspa8, Hsc70) belongs to the heat shock protein 70 family. Hsc70 is a constitutively expressed molecular chaperone and 85% identity with heat-shock protein 70. Hsc70 functions as an ATPase and facilitates correct folding of nascent polypeptides [30]. The reduction of Hsc70 was associated with induction of apoptosis and normal aging in primate retina [31]. Furthermore, Hsc70 is highly associated with vimentin (Vim), one of the down-regulated proteins, and suggested to participate in stabilizing the intermediate filaments [32]. Vimentin, a major constituent of the intermediate filament family of proteins, is ubiquitously expressed in normal mesenchymal cells and is known to maintain cellular integrity and provide resistance against stress. The expression of Vim is controlled by cell configuration, and down-regulation of Vim causes the impaired migration and adhesion of cancer cells [33,34]. ATP synthase subunit β, mitochondria 1 (Atp5b) is a subunit of mitochondrial ATP synthase. A decreased expression of Atp5b has been documented in malignant tumors when compared with its level in normal tissues [35,36]. Mitochondria aggregate as clusters and move on the cytoskeleton. CLC depolymerizes microtubules to disrupt mitochondria chains, which are required for mitochondria to respond to the alteration in local energy demands [37]. Protein disulfide-isomerase A6 (Pdia6) is a chaperone protein that inhibits aggregation of misfolded proteins. The ER resident Pdia6 is an attenuator of the unfolded protein response; it interacts covalently with activated IRE1 and facilitates the decay of active IRE1 [38]. Eif4a2 (eukaryotic initiation factor 4A2), a protein-synthesis initiation factor, belongs to the DEAD box helicase family, functions to resolve RNA secondary structures and participates in the binding of mRNA to the 40S ribosome during the initiation phase of eukaryotic translation [39]. The Eif4a2 mRNA synthesis and translation was found to associate preferentially with the growth-arrested (quiescent) state [40]. The 14-3-3 protein ε (Ywhae), a member of highly conserved 14-3-3 family proteins, binds to phosphoserine-containing proteins to mediate signal transduction, and heterodimerizes with 14-3-3 protein ξ (Ywhaz). The down-regulation of 14-3-3 proteins (β, σ, ξ) was found during the transition of anagen to catagen [41]. A proteomic study showed that Ywhae is critical for mediating the PPI in the mitogen-activated protein kinases (MAPK) signal module for modulating NF-κB translocation/ activity [42]. Down-regulation of Ywhae and annexin A5 was observed in lung adenocarcinoma treated with tobacco-specific nitrosamine 4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone [43]. Annexin A5 (anxA5) is a phospholipid binding protein that efficiently binds to phosphatidylserine and it was down-regulated in an epithelial cell treated with rifaximin [44].

The Seventeen Proteins Up-Regulated in DPCs by CLC
The expression of 17 proteins in DPCs was up-regulated by CLC (Table 2). Peroxiredoxin 6 (Prdx6), one of the highly conserved peroxidase family proteins, exhibits anti-oxidant character and prevents cell from peroxide-induced cytotoxicity. It has been shown that NF-κB acts as a stress-sensing molecule that determines optimal regulation of Prdx6 transcription for fine tuning the delicate redox balance necessary for the maintenance of cellular functions [45]. Protein disulfide-isomerase A3 (Pdia3) gene encodes a protein in the endoplasmic reticulum that interacts with lectin chaperones calreticulin and calnexin to modulate the folding of glycoproteins that are newly synthesized. The protein PDIA3 has been found to be active in several other locations and reactions, for example, interactions in the nucleus which involve DNA repair, DNA damage recognition, and apoptosis [46].  The 26S proteasome non-ATPase regulatory subunit 13 (Psmd13) is a regulatory subunit of the 26S proteasome, which involved in the ATP-dependent degradation of ubiquitinated proteins. Recent report showed that Psmd13 gene silencing suppressed the production of proinflammatory mediators by modulating ubiquitin-proteasome system-mediated neuro-inflammation via the down-regulation of IκBα degradation and NF-κB activation in LPS-stimulated BV2 microglia [47]. Proteasome subunit β type-7 (Psmb7) is also a component of the proteasome, a macromolecular machine integral to cellular proteolytic degradation capability [48]. Proteasome subunit α type-2 (Psma2) gene encodes a member of the peptidase T1A family that is a 20S core α subunit strongly interacts with other proteasomes. It may have a potential regulatory effect on other components of the proteasome complex through tyrosine phosphorylation [48].
The 60S acidic ribosomal protein P0 (Rplp0) is a component of the 60S subunit of ribosomes and ribosomes catalyze protein synthesis. Proteomic study showed that Rplp0 was up-regulated in Jurkat cells during heat stress-induced apoptosis [49]. Protein DJ-1 (Park7), a peptidase, acts as a redox-sensitive chaperone and plays an important role in attenuating oxidant stress and promoting cell survival [50]. Through binding and inhibiting a deubiquitination enzyme, Cezanne, Park7 enhanced NF-κB nuclear translocation and cell survival [51]. PDZ and LIM domain protein 1 (Pdlim1, Clim1), a cytoskeletal protein, involves in the assembly and maintenance of stress fibers [52]. [53]. It has been reported that DPCs exhibited HSC activity and repopulated the mouse hematopoietic system [54]. Latexin inhibits the self-renewal of HSCs by facilitating the lodgment of HSCs within a bone marrow niche to maintain HSC homeostasis [55]. There were two upstream NF-κB regulatory elements involved in the control of cell type-specific expression of rat latexin gene [56]. Transgelin (TAGLN) is an actin cross-linking/polymerization protein that belongs to the family of actin-associated proteins, and involved in the migration of epithelial cells by interacting with actin or promoting podosome formation [57]. Transgelin in the kidney was up-regulated in repopulating mesangial cells in vivo and supported their migratory and proliferative repair response after injury [58].

Latexin (Lxn) is an endogenous carboxypeptidase inhibitor (CPI), inhibitor for metallo-carboxypeptidases (MCPs) and a negative regulator of hematopoietic stem cells (HSCs)
Cytoplasmic malate dehydrogenase (Mdh1) is important in transporting NADH equivalents across the mitochondrial membrane, controlling tricarboxylic acid (TCA) cycle pool size [59]. A proteomic analysis indicated that the activation of the gluconeogenesis caused up-regulation of Mdh and annexin A1 in thick ascending limb of Henle's loop (TALH cells) in response to osmotic stress [60]. Annexin A1 (Anxa1) has been shown to promote metastasis through constitutively activate NF-κB by interaction with the IκB kinase (IKK) complex in breast cancer cells [61].
Mpst (3-mercaptopyruvate-sulfurtransferase), an antioxidant enzyme, is one of the three enzymes involved in synthesizing H2S, which is increased under inflammatory conditions or sepsis in mammalian tissues. Knockout mice study indicated a positive interlink between Anxa1 and H2S pathway in the control of inflammation [62]. Glutathione S-transferases (GSTs) play an important role in the protection of cells against xenobiotics and lipid hydroperoxides generated by oxidative stress. Glutathione S-transferase P1 (Gstp1) prevented LPS-induced excessive production of pro-inflammatory factors through inhibiting LPS-induced MAPK including ERK, JNK and p38 as well as NF-κB activation in RAW264.7 macrophage-like cells [63].
Recent proteomic study showed that human keratinocytes treated with glucose oxidase (source of oxidative stress) could up-regulate expression of several proteins involving in stress response, including Pdia3 and Protein DJ-1, and in cytoskeletal organization, including Annexin A2 and F-actin-capping protein subunit β (Capzb) [64]. F-actin-capping proteins bind in a calcium-independent manner to the fast growing ends of actin filaments (barbed end), thereby playing a role in the regulation of cell morphology and cytoskeletal organization [65]. Nevertheless, the detailed molecular mechanisms of many identified proteins affected by CLC in DPCs require further investigation.

Major Identified Proteins Associated with Protein Turnover and Metabolism
Base on Protein ANalysis THrough Evolutionary Relationships (PANTHER) classification system [66], the functional distribution of seven down-regulated and 17 up-regulated proteins were each classified into five groups (Tables 1 and 2), respectively. The major functions of these proteins were associated with protein turnover (37.5%), metabolism (16.7%), structure (16.7%), and signal transduction (16.7%) ( Figure 4A).  Tables 1 and 2 for description). The PPI were shown in the confidence view produced by STRING database. Thicker blue lines represented the stronger associations.

The Interactomic Analysis Derived a PPI Network with 3 Protein Clusters
Interactomics holds great promise in understanding the molecular mechanism of cells affected by medications. To examine the clusters of 24 differentially expressed proteins and their associated partners and their PPIs, the STRING database [67] was used to deduce enriched protein clusters and generated a PPI network ( Figure 4B). The computed predictions indicated that the PPI map contained 1 large (I), and 2 small (II, III) protein clusters.
The symmetric connected cluster I contains proteins of multicatalytic proteasome complex, which participates in protein metabolic process, including proteasome subunit α type-2 (Psma2) and other α-type subunits, proteasome subunit β type-7 (Psmb7) and other β-type subunits, as well as 26S proteasome non-ATPase regulatory subunits 10, 13 and 14 (Psmd10, Psmd13 and Psmd14). The α-type and β-type subunit complexes are the subunits of 20S core of proteasome, whereas Psmd10, Psmd13 and Psmd14 belong to the 19S regulator [48]. Ubiquitin-labeled proteins are degraded via the ubiquitin-proteasome pathway, which controls the distribution, stability and function of cellular proteins [48,68]. The MDA vincristine enhanced proteasome-mediated degradation and disappearance of tubulin in neural cells, which is inhibited by a proteasome inhibitor MG132 [69].

Ubiquitin C (UbC) was Predicted as a Controlling Hub Protein
Controlling hubs are the principal proteins with a great number of interactions in a PPI network. We search STRING database for predicting controlling hubs against 24 identified proteins as seeds. The ubiquitin C (UbC) was predicted to be the computed major controlling hub; this protein acts as the connector of three clusters ( Figure 4B). Ubc is a protein-coding gene and encodes a polyubiquitin precursor. Free ubiquitin is created from UbC by deubiquitinases [70]. Conjugation of ubiquitin monomers or polymers can lead to various effects within a cell, depending on the residues to which ubiquitin is conjugated. Ubiquitination has been associated with protein degradation, DNA repair, cell cycle regulation, kinase modification, endocytosis, and regulation of other cell signaling pathways [71,72]. The ubiquitin system has been known to play diverse roles of in NF-κB activation [24]. Therefore, the network analysis reveals potential roles of ubiquitin-proteasome system in the pathogenic effect of CLC on PDCs.

Preparation of CLC Solution and Primary Culture of SD Rat Vibrissa DPCs
CLC was dissolved in DMSO to make a stock solution of 100 mM, and then stored in dark at 4 °C in aliquots. CLC and DMSO were purchased from Sigma (St. Louis, MO, USA). Primary culture of SD rat vibrissa DPCs was prepared as previous described [22,29].

Organ Culture of SD Rat Vibrissa HF for Hair Fiber Elongation Analysis
Fresh intact vibrissa HFs of SD rats were isolated with a jewelry-maker forceps [73]. Isolated cleaned vibrissa HFs, which were first sterilized with an antibiotic wash medium (RPMI-1640 with 20% FCS, 2 mM L-glutamine, 50 IU/mL penicillin, and 50 µg/mL streptomycin, 10 µg/mL tetracycline, 150 µg/mL chloramphenicol, 10 µg/mL gentamycin sulfate, 4% fungizone) for 30 min in a 5% CO2 humidified atmosphere incubator at 37 °C. Sterilized HFs were maintained in William's E medium (containing 50 IU/mL penicillin, and 50 µg/mL streptomycin) in individual wells of 24-well plates for a 7-day culture period at 37 °C in a 5% CO2 humidified atmosphere incubator. The HF organ culture assay was used to evaluate the effect of CLC on the elongation of HF fiber. CLC was added to the HF culture on the first day at final concentrations of 0.1 and 1 µM, respectively. Each HF shaft length was determined under a binocular dissecting microscope, fitted with an eyepiece measuring graticule and measured after 7 days. The length of the HF fiber was defined as the distance between its base and cut tip of the fiber. Cell culture reagents were purchased from Hyclone (Logan, UT, USA). Eight-week-old SD rats were obtained from National laboratory animal center (Taipei, Taiwan).

WB and qRT-PCR Analysis
The protein concentration of DPC lysate was determined using the Bradford protein assay reagent. Samples (20 µg) were resolved in a 15% SDS-PAGE and proteins were electro-transferred onto an Immobilon-P membrane. The protein levels of ALP, phospho-IκBα and GADPH were detected by WB with anti-TNAP, anti-phospho-IκBα, and anti-GADPH, respectively, following the procedure described previously [22]. Immunoblots were scanned and quantified by ImageJ. Trizol reagent was used to isolate total DPC RNA. The qRT-PCR was performed to measure the levels of ALP and GAPDH mRNAs according to the previous described protocol [22].

Two-DE and Image Analysis
DPCs were treated or untreated with 2 µM CLC for 6 h. Trypan blue staining indicated the viability of DPC is greater than 95%. DPCs were washed with cold wash buffer (10 mM Tris-HCl, pH 8.0) and solubilized by sonication in lysis buffer (7 M urea, 2 M thiourea, 2% CHAPS, bromophenol blue, 2% DTT and 0.5% IPG buffer). A whole cell lysate (250 µg) was applied onto an IPG strip (130 × 3 × 0.5 mm, pH 4-7) which was rehydrated for 16 h at 20 °C. The IEF was performed on an Ettan IPGphor II (GE Healthcare, Pittsburgh, PA, USA) with 200 V for 1 h (200 V-h), 500 V for 1 h (500 V-h), 1000 V for 1 h (1000 V-h), a gradient of 8000 V (4000 V-h) and a "step-n-hold" of 8000 V for 6.75 h. After IEF, IPG strips were equilibrated with the equilibration buffer (6 M urea, 2% SDS, 0.05 M Tris-HCl, pH 8.8, 20% glycerol and 0.5% (w/v) DTT) for 15 min and blocked in the equilibrating buffer by substituting DTT with 4.5% (w/v) iodoacetamide for 15 min. Each focused IPG strip was separated by a 12.5% SDS-PAGE under a constant voltage of 300 V till finish. Sypro Ruby-stained 2-DE gels were scanned with Typhoon 8600 system (GE Healthcare, Pittsburgh, PA, USA) and analyzed using PDQuest 2-DE image analysis software (Bio-Rad, Hercules, CA, USA). To reduce variations among 2-DE gels resulted from the procedure of staining and destaining, a triplicate 2-DE analysis per sample was performed. Manual matching was done when necessary. 2-DE reagents were purchased from GE Healthcare (Pittsburgh, PA, USA).

In-Gel Protein Digestion and Protein Identification by Nano-LC-ESI-MS/MS
Protein spots of interest were recovered from gels and dried in vacuum. In-gel tryptic digestion was performed and the resulting peptide mixture was extracted in 50% ACN containing 0.1% (w/v) TFA. Desalted extracted trypsin digested peptides were separated and analyzed by nano-LC-ESI-MS/MS as previous described [22]. The identification of protein was based on MS/MS peak lists, which were generated by the Mascot Distiller program (v2.3.2, Matrix Science, London, UK) using default parameters with Mascot online search engine (v2.2.0.6, Matrix Science, London, UK) against SwissProt 2011 database (524420 sequences; Rattus norvegicus taxonomy). All the searching parameters have been described previously [22].

PANTHER and STRING Database Analysis
Accession numbers of identified proteins were analyzed for their biological process and interactions using PANTHER classification system (available online: http://www.pantherdb.org) and STRING database version 9.1 (available online: http://string.embl.de) [67], respectively.

Statistical Data Analysis
Data were processed in a model of Student's paired t-tests to evaluate the significance of differences using SPSS analytical software. The value of p for each data is * p < 0.05, ** p < 0.01 vs. control group.

Conclusions
The present study demonstrates the CLC-modulated DPC interactome provided by STRING database and suggests following conclusions. First, a CLC-induced multi-factor pathogenesis on DPCs is derived from this proteomic data. Second, a sub-lethal dose of CLC reduced the ALP expression of DPC to cause inactivation of DPCs. Third, proteins with functions of protein turnover, metabolism, structure, and signal transduction are greatly affected by CLC. Forth, ubiquitin C (UbC) is predicted by STRING database as a major controlling hub, implying the involvement of ubiquitin-proteasome system in mediating CLC effect on DPCs. However, the detailed molecular mechanism of the pathogenic effect of CLC on DPCs will be further investigated. Finally, the interactomic research can be applied to discover proteins in DPCs related to drug-induced effluvium.

Author Contributions
The work presented here was carried out as collaboration between all of the authors. Ching-Wu Hsia performed the experiment, participated in the design of the study and drafted the manuscript. Ming-Yi Ho carried out Nano-LC-ESI-MS/MS for protein identification. Hao-Ai Shui participated in analyzing data and provided critical regents and instruments. Chong-Bin Tsai and Min-Jen Tseng defined the research theme, participated in its design and coordination, analyzed and interpreted the data and helped to draft and revised the manuscript.