BioMed Research International

BioMed Research International / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 3235786 |

Weiguo Lin, Chenhao Tong, Weiguang Zhang, Wenda Cen, Yali Wang, Jiandong Li, Zhiyang Zhu, Jianhua Yu, Baochun Lu, "Silencing ARAF Suppresses the Malignant Phenotypes of Gallbladder Cancer Cells", BioMed Research International, vol. 2020, Article ID 3235786, 8 pages, 2020.

Silencing ARAF Suppresses the Malignant Phenotypes of Gallbladder Cancer Cells

Academic Editor: Kui Li
Received06 May 2020
Revised15 Jul 2020
Accepted22 Jul 2020
Published19 Aug 2020


ARAF is a member of the RAF kinase family that is necessary for mitogen-activated protein kinase (MAPK) activation in various malignancies, including lung, colorectal, pancreatic, and breast cancers. As the most common biliary tract tumor, gallbladder cancer (GBC) seriously harms human health while the function of ARAF in GBC remains elusive. Here, we found that ARAF expression was upregulated in gallbladder cancer tissues. In vitro, ARAF silencing mediated by RNA interference effectively inhibited cell proliferation, colony formation, migration, and invasion of GBC cells. Moreover, knocking down ARAF suppressed tumor growth in vivo. Our results indicated that ARAF functions as an oncogene in GBC and, thus, could be a potential therapeutic target for GBC.

1. Introduction

Gallbladder cancer (GBC) is an aggressive malignancy of the biliary tract that originates from the gallbladder and cystic duct mucosal epithelia [1]. As the most common biliary tract cancer, GBC accounts for 80%-95% of all biliary malignancies and has a dismal prognosis [2, 3]. The recent optimization of medical auxiliary examinations and the widespread application of laparoscopic cholecystectomy have significantly increased the detection rate of gallbladder cancer; however, its prognosis has not improved because of late-stage diagnoses, high recurrence rates, and metastatic features [4]. Although surgical resection remains the most effective treatment for GBC, most patients are diagnosed with advanced-stage disease, meaning they are not candidates for surgery [5, 6]. What is worse, GBC has extremely poor sensitivity to radiotherapy and chemotherapy. Therefore, clearing the underlying molecular mechanisms of GBC tumorigenesis and metastasis will provide a theoretical basis for improving its diagnosis and treatment.

Located on human chromosome band Xp11.3, ARAF belongs to the serine/threonine protein kinase gene family [7]. Similar to other RAF family members, ARAF transduces mitogen-activated protein kinase (MAPK) signaling from RAS to MEK and ERK, thus promoting cell proliferation, differentiation, migration, and survival [8, 9]. The RAS-RAF-MEK-ERK cascade is considered to be a therapeutic target in various cancers [10, 11].

Early studies on the RAF family focused on B-Raf and C-Raf kinases, resulting in little understanding of the biological function of ARAF. Recent studies focused on the role of ARAF in tumor progression have made significant impact on the field. Early cancer sequencing studies identified high-copy number gains as well as oncogenic driver mutations in ARAF in lung cancer patients [12]. In 2014, a study demonstrated that ARAF was required for MAPK activation in a variety of cancer types (e.g., colorectal, pancreatic, and breast cancers) and further verified that ARAF enhanced the migration and invasive ability of these tumor cells [13]. Other studies reported that ARAF mutations could drive lung cancer and that the RAF-targeted kinase inhibitor sorafenib improved the prognosis of advanced lung cancer patients, thus providing a new opportunity for lung cancer treatment [14]. These findings suggested that ARAF could be a therapeutic target in numerous cancers. However, the functional role of ARAF in GBC has not been verified.

Here, we explored the functional roles of ARAF in relation to GBC tumorigenesis and progression. As shown in our results, both ARAF mRNA and its encoding protein were overexpressed in GBC compared with nontumoral tissues. After the expression level of ARAF gene was downregulated by RNA interference technology, the tumor phenotype of gallbladder cancer cells was considerably affected both in vivo and in vitro, which showed that the cell proliferation, metastasis, and other abilities were weakened. Therefore, we believe that ARAF promotes the development of GBC and regulates its growth and metastasis.

2. Materials and Methods

2.1. Clinical Tissue Samples

GBC and normal gallbladder tissues were obtained at Shaoxing People’s Hospital. All patients signed informed consent documents before inclusion in the study. Informed consent document and tissue acquisition protocol were approved by the Ethics Committee of Shaoxing People’s Hospital (Shaoxing, China). Cancer tissues were collected from GBC patients, while nontumoral tissues were harvested from patients with gallbladder polyps. Fresh tissues were stored in liquid nitrogen prior to RNA and protein extraction.

2.2. Cell Culture

The GBC cell line GBC-SD was purchased from the Chinese Academy of Sciences Shanghai Branch Cell Bank (Shanghai, China), and the SGC-996 cell line was obtained from Dr. Ying-Bin Liu’s lab at Xin Hua Hospital Affiliated to Shanghai Jiao Tong University School of Medicine, China. Both cell lines were cultured in RPMI-1640 medium (cat. no. GNM-31800-S; USEN Biological Technology Co., Ltd., Shanghai, China) with 10% fetal bovine serum (FBS; cat. no. 16140071; Gibco; Thermo Fisher Scientific, Inc., Waltham, MA, USA), 100 IU/ml penicillin, and 100 μg/ml streptomycin in a 37°C incubator with 5% CO2.

2.3. siRNA Transfection

To downregulate ARAF expression in GBC cell lines, ARAF-targeting siRNA (5-GGGATGGCATGAGTGTCTA-3) was purchased from RiboBio (Guangzhou, China). Control siRNA was also obtained from RiboBio, and the control sequence was not public. Before transfection, cells were seeded into dishes at 50%–60% confluence. Transfection was performed using Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocol.

2.4. RNA Extraction and RT-qPCR

TRIzol (Invitrogen) was used to isolate total RNA from tissues or cells. Then, RNA was reverse transcribed into cDNA using the PrimeScript Reagent Kit (Takara, Shiga, Japan) according to the manufacturer’s instruction. RT-qPCR was performed using the TB Green Kit (Takara) on an ABI 7500 Real-time PCR system (Applied Biosystems, Foster City, CA, USA). GAPDH was used as an endogenous control. Relative ARAF mRNA level was determined by the 2-ΔΔCt method [15]. The primer sequences are listed as follows: ARAF forward 5-CCCACATTCCAAGTCACCAGCA-3 and reverse 5-CCTCCCAGTAATAGCCTGAGTC-3 and GAPDH forward 5-GTCTCCTCTGACTTCAACAGCG-3 and reverse 5-ACCACCCTGTTGCTGTAGCCAA-3.

2.5. Western Blot Analysis

Western blotting was used to detect ARAF protein levels in tissues and cell lines. RIPA lysis buffer containing 1% phenylmethylsulfonyl fluoride (Beyotime Institute of Biotechnology, Nantong, China) was used to extract the total protein from tissues and cells. The protein concentration was determined using the BSA method (Beyotime Institute of Biotechnology). Briefly, 30 μg protein was loaded onto 10% SDS–PAGE gels, electrophoresed, and transferred onto to polyvinylidene fluoride membranes, which were blocked by 5% skim milk powder in TBST, incubated with primary antibodies overnight, and then with horseradish peroxidase-conjugated secondary antibodies (cat. nos. A0208 and A0216; Beyotime Institute of Biotechnology) at 1 : 10000 dilution for 2 h at room temperature. Primary antibodies against ARAF (dilution 1 : 1000; cat. no. 4432) and phospho-p44/42 MAPK (Erk1/2, Thr202/Tyr204, dilution 1 : 1000; cat. no. 4370) were purchased from Cell Signaling Technology (Danvers, MA, USA). Primary antibodies against PCNA (dilution 1 : 1000; cat. no. 60097-1-Ig), cyclin D1 (dilution 1 : 1000; cat. no. 26939-1-AP), and β-actin (dilution 1 : 1000; cat. no. 20536-1-AP) were purchased from Proteintech (Rosemont, IL, USA). β-Actin was used as the endogenous control. Immunoreactive bands were visualized using a chemiluminescence solution (Beyotime Institute of Biotechnology).

2.6. Cell Proliferation Assays

Both GBC cell lines were transfected with ARAF siRNA or siRNA control and, 6 h after transfection, were seeded into 96-well plates at 2000 cells per well. Every 24 h, cell growth was evaluated using Cell Counting Kit-8 (Beyotime Institute of Biotechnology). According to the manufacturer’s protocol, 10 μl of CCK-8 reagent was added to each well and incubated for 2 h. Then, cell viability was measured with an enzyme-labeling instrument (BioTek, Winooski, VT, USA) at 450 nm.

2.7. Colony Formation Assays

After transfection, 200 cells were seeded into 35 mm dishes and then cultured for 2 weeks. After fixation with 4% paraformaldehyde and staining with 0.1% crystal violet solution, colonies of >50 cells were counted.

2.8. Wound Healing Assays

Cells were seeded and grown to confluence on 35 mm cell dishes. Six hours posttransfection, a 10 μl pipette tip was used to scratch the confluent monolayers. Cells were then cultured in serum-free medium (inhibiting cell proliferation), and after 48 h, images of the wounds were captured at 100x magnification. Wound healing was quantified as the average linear speed of the wound edges. The scratch area was calculated by ImageJ software, and the cell mobility was calculated by the following formula: the .

2.9. Transwell Assays

Briefly, transfected cells per 100 μl in serum-free medium were added to the upper chamber of the insert (Corning Inc., Corning, NY, USA), while the lower chamber was filled with 0.5 ml of medium containing 20% FBS. Chambers were incubated for 24 h, and then, invasive cells were fixed with 4% paraformaldehyde and stained with 0.1% crystal violet. Four fields of view were taken under a 100x optical microscope, and the number of cells entering the lower chamber in each field was calculated and averaged.

2.10. Xenograft Formation Assays

All procedures were approved by the Ethics Committee of Shaoxing People’s Hospital and conformed to the Care and Use of Laboratory Animals guide published by the US National Institutes of Health (NIH Publication No. 85-23, revised 1996). Six-week-old athymic nude mice were supplied by Shanghai SLAC Laboratory Animal Co., Ltd. (Shanghai, China). The mice had ad libitum access to food and water and were maintained at 20°C, with 50% humidity under 12 : 12-h light-dark cycles. Then, cells were suspended in 0.2 ml PBS and subcutaneously injected into the back of nude mice. siRNA was injected intratumorally in a volume of 100 μl once every 3 days. At regular intervals, tumor sizes were measured, and tumor volumes were calculated according to the following formula: [16]. The mice with GBC-SD cells or SGC-996 cells were sacrificed via cervical dislocation under isoflurane anesthesia after 42 days or 24 days, respectively. Finally, all tumor specimens were collected and weighed.

2.11. Statistical Analysis

All experiments were repeated at least three times, and data are presented as the . Student’s -test was used to determine statistical significance between two groups. One-way ANOVA followed by the Tukey–Kramer adjustment was used to examine differences among multiple groups. All statistical analyses were conducted using SPSS v21.0 (IBM, Armonk, NY, USA), and was considered statistically significant.

3. Results

3.1. ARAF Expression Is Upregulated in GBC Tissues

To compare ARAF mRNA expression between GBC and nontumoral samples, RT-qPCR was performed, and the average ARAF mRNA expression of 12 nontumoral tissues was defined as the baseline expression of normal tissues. As shown in Figure 1(a), ARAF mRNA expression was significantly higher in GBC tissues than that in nontumoral tissues according to the results of RT-qPCR (Figure 1(a)). Additionally, ARAF protein was also significantly increased in GBC (Figure 1(b)).

3.2. Silencing ARAF Inhibits GBC Cell Proliferation

ARAF knockdown was employed to study its function in GBC cells. We examined ARAF expression in GBC-SD and SGC-996 cells transfected with ARAF siRNA or siRNA control. After siRNA transfection, ARAF mRNA and protein levels were both significantly lower than controls (Figures 2(a) and 2(b)). It was demonstrated that the ARAF siRNA successfully silenced endogenous ARAF in GBC cells. Inhibiting the rapid growing of cancer cells is an important way to treat cancers. As shown in Figure 2(c), the CCK-8 assays showed that ARAF siRNA inhibited the proliferation of different GBC cell lines, including GBC-SD cells and SGC-996 cells. Interestingly, PCNA, a reliable indicator of cell proliferation, was also higher in the control siRNA group compared with ARAF knockdown (Figure 2(d)). Furthermore, colony formation was significantly reduced in the ARAF siRNA group, compared with controls (Figures 2(e) and 2(f)).

3.3. Silencing ARAF Inhibits the Migration and Invasion of GBC Cells

Metastasis is the most discouraging phenomenon in cancers and is also an important factor which leads the patients with the late stage of GBC lose the operation opportunity. The role of ARAF on the migration and invasion of GBC cells was further explored by wound healing and transwell assays. As shown in Figures 3(a) and 3(b), ARAF knockdown significantly attenuated cell migration compared with controls. Because we found that the invasiveness of SGC-996 cells is too poor to use for transwell assays, GBC-SD cells were employed to investigate the role of ARAF on cell invasion. We found that the invasion of GBC-SD cells was remarkably suppressed after ARAF was knockdown (Figure 3(c)). These results indicate that silencing ARAF inhibits the migration and invasion of GBC cells.

3.4. Silencing ARAF Suppresses Xenograft Tumor Growth In Vivo

Nude mouse xenograft formation assays were performed to investigate the biological significance of silencing ARAF in GBC by subcutaneously injecting GBC-SD and SGC-996 cells. Both the tumor volume and weight of nude mice in the ARAF silencing group were significantly reduced, compared with the control group (Figures 4(a)4(c)). It indicates that silencing ARAF also effectively suppresses tumor growth in vivo, consistent with the results in vitro.

4. Discussion

GBC is the seventh most common tumor worldwide and has a terrible prognosis [17]. The main reasons for its poor prognosis are late diagnosis, early metastasis, and limited therapeutic options, which make it urgently necessary to uncover the molecular mechanism of GBC.

Among patients with gallbladder cancer, it is clear that the proportion of women is significantly higher than that of men. Given the gender differences in gallbladder cancer prevalence, we tried to find out protooncogenes on the X chromosome or antioncogenes on the Y chromosome at the beginning of our study. After a literature search and preliminary experiments, the ARAF gene, located on the X chromosome, got our attention. Interestingly, the expression of ARAF was significantly higher in female patients, compared with that in male patients (Figure 4(d)). Given the potential association between the aberrant ARAF expression and the sexual dimorphism of GBC, we thought ARAF may be an oncogene which is more worth to be studied.

As a new star of the family of Raf kinases, ARAF plays an important role in the regulation of many cellular functions, including differentiation, cell proliferation, and transformation [18]. In the mouse experiment with gene knockout of ARAF, mouse embryonic fibroblasts delayed entering the S phase of cell cycle, indicating that ARAF maintained the progress of cell cycle [19]. ARAF also has been shown to play an important role in the proliferation of vascular smooth muscle while inhibiting the activity of Raf kinase could be used as a treatment for vascular hyperplastic diseases [20]. Many kinases which have regulating function during the process of embryonic developing always are potential protooncogenes, including ARAF.

Previous investigation of the protooncogene ARAF demonstrated that ARAF played an obligatory role in promoting MAPK activity as a kinase [13]. MAPK signaling, represented by the phosphorylation of Erk, plays a key role during the cell proliferation, migration, and invasion in various cancers. To clear whether silencing ARAF suppresses the malignant phenotypes of GBC cells through regulating MAPK signaling, phosphorylation of Erk was examined in our study. We found that the phosphorylation of Erk was significantly inhibited in GBC cells when ARAF was downregulated (Figure 4(e)). Our results showed that the silencing of ARAF could produce an inhibitory effect on GBC cell proliferation and colony formation. Coincidentally, a previous study on murine embryonic stem cells also revealed that ARAF is required for Erk activation and involved in the growth and colony formation [21]. More importantly, silencing ARAF limited the growth of xenograft tumors in nude mice. Cyclin D1 is induced by Raf/MAPK/ERK cascade and plays a key role during proliferation in various cancers [22, 23]. Just like the fact that the phosphorylation of Erk was inhibited, the same decreasing trend of cyclin D1 expression was observed after ARAF was knocked down (Figure 4(e)). Our results also demonstrated that the silencing ARAF impaired the migration and invasion of GBC cells. Interestingly, a previous work about trophoblasts reported that ARAF-mediated activation of the integrin/Erk signaling pathway promotes trophoblast migration and invasion [24]. Taken together, ARAF silencing suppresses the malignant phenotypes of gallbladder cancer cells, and the mechanism may be associated with regulating Erk/cyclin D1 axis.

In conclusion, our results demonstrate that ARAF expression is highly expressed in human gallbladder cancer, and ARAF silencing has an inhibitory effect on various phenotypes of GBC. Based on these findings, ARAF should be regarded as oncogene in GBC progression. Targeting ARAF therefore represents a potential therapeutic target for GBC.

Data Availability

Data are available on request.

Conflicts of Interest

The authors confirm that there are no conflicts of interest.

Authors’ Contributions

Weiguo Lin and Chenhao Tong contributed equally to this work.


The work was sponsored by Zhejiang Provincial Natural Science Foundation of China under grant no. LY19H160016 and LQ14H160001; Zhejiang Province Public Welfare Technology Application Research Project under grant no. LGF18H030008; National Natural Science Foundation of China (NSFC) under grant no. 81602044; Zhejiang Provincial Medical and Health Science and Technology Project under grant no. 2020RC127, 2019ZD057, 2018KY836, and 2018RC077; and Public Welfare Technology Application Research Projects of Shaoxing under grant no. 2014B70073. The authors are grateful to Dr. Ying-Bin Liu for providing SGC-996 cell line.


  1. U. Dutta, “Gallbladder cancer: can newer insights improve the outcome?” Journal of Gastroenterology and Hepatology, vol. 27, no. 4, pp. 642–653, 2012. View at: Publisher Site | Google Scholar
  2. R. Hundal and E. A. Shaffer, “Gallbladder cancer: epidemiology and outcome,” Clinical Epidemiology, vol. 6, pp. 99–109, 2014. View at: Publisher Site | Google Scholar
  3. S. Zhao, Y. Cao, S. B. Liu et al., “The E545K mutation of PIK3CA promotes gallbladder carcinoma progression through enhanced binding to EGFR,” Journal of Experimental & Clinical Cancer Research, vol. 35, no. 1, p. 97, 2016. View at: Publisher Site | Google Scholar
  4. A. Sharma, K. L. Sharma, A. Gupta, A. Yadav, and A. Kumar, “Gallbladder cancer epidemiology, pathogenesis and molecular genetics: recent update,” World Journal of Gastroenterology, vol. 23, no. 22, pp. 3978–3998, 2017. View at: Publisher Site | Google Scholar
  5. M. Rakic, L. Patrlj, M. Kopljar et al., “Gallbladder cancer,” Hepatobiliary Surg Nutr, vol. 3, no. 5, pp. 221–226, 2014. View at: Publisher Site | Google Scholar
  6. P. Portincasa, A. di Ciaula, O. de Bari, G. Garruti, V. O. Palmieri, and D. Q. H. Wang, “Management of gallstones and its related complications,” Expert Review of Gastroenterology & Hepatology, vol. 10, no. 1, pp. 93–112, 2015. View at: Publisher Site | Google Scholar
  7. M. J. Robinson and M. H. Cobb, “Mitogen-activated protein kinase pathways,” Current Opinion in Cell Biology, vol. 9, no. 2, pp. 180–186, 1997. View at: Publisher Site | Google Scholar
  8. M. J. Garnett and R. Marais, “Guilty as charged: B-RAF is a human oncogene,” Cancer Cell, vol. 6, no. 4, pp. 313–319, 2004. View at: Publisher Site | Google Scholar
  9. A. P. Rebocho and R. Marais, “ARAF acts as a scaffold to stabilize BRAF:CRAF heterodimers,” Oncogene, vol. 32, no. 26, pp. 3207–3212, 2013. View at: Publisher Site | Google Scholar
  10. H. Yang, Z. Wan, C. Huang, H. Yin, and D. Song, “AMPH-1 is a tumor suppressor of lung cancer by inhibiting Ras-Raf-MEK-ERK signal pathway,” Lasers in Medical Science, vol. 34, no. 3, pp. 473–478, 2019. View at: Publisher Site | Google Scholar
  11. R. Roskoski Jr., “Targeting oncogenic Raf protein-serine/threonine kinases in human cancers,” Pharmacological Research, vol. 135, pp. 239–258, 2018. View at: Publisher Site | Google Scholar
  12. W. Lee, Z. Jiang, J. Liu et al., “The mutation spectrum revealed by paired genome sequences from a lung cancer patient,” Nature, vol. 465, no. 7297, pp. 473–477, 2010. View at: Publisher Site | Google Scholar
  13. J. Mooz, T. K. Oberoi-Khanuja, G. S. Harms et al., “Dimerization of the kinase ARAF promotes MAPK pathway activation and cell migration,” Science Signaling, vol. 7, no. 337, article ra73, 2014. View at: Publisher Site | Google Scholar
  14. M. Imielinski, H. Greulich, B. Kaplan et al., “Oncogenic and sorafenib-sensitive ARAF mutations in lung adenocarcinoma,” The Journal of Clinical Investigation, vol. 124, no. 4, pp. 1582–1586, 2014. View at: Publisher Site | Google Scholar
  15. K. J. Livak and T. D. Schmittgen, “Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method,” Methods, vol. 25, no. 4, pp. 402–408, 2001. View at: Publisher Site | Google Scholar
  16. S. Naito, A. von Eschenbach, R. Giavazzi, and I. J. Fidler, “Growth and metastasis of tumor cells isolated from a human renal cell carcinoma implanted into different organs of nude mice,” Cancer Research, vol. 46, no. 8, pp. 4109–4115, 1986. View at: Google Scholar
  17. Y. M. He, Z. L. Zhang, Q. Y. Liu et al., “Effect of CLIC1 gene silencing on proliferation, migration, invasion and apoptosis of human gallbladder cancer cells,” Journal of Cellular and Molecular Medicine, vol. 22, no. 5, pp. 2569–2579, 2018. View at: Publisher Site | Google Scholar
  18. S. An, Y. Yang, R. Ward, Y. Liu, X. X. Guo, and T. R. Xu, “A-Raf: a new star of the family of raf kinases,” Critical Reviews in Biochemistry and Molecular Biology, vol. 50, no. 6, pp. 520–531, 2015. View at: Publisher Site | Google Scholar
  19. K. Mercer, S. Giblett, A. Oakden, J. Brown, R. Marais, and C. Pritchard, “A-Raf and Raf-1 work together to influence transient ERK phosphorylation and Gl/S cell cycle progression,” Oncogene, vol. 24, no. 33, pp. 5207–5217, 2005. View at: Publisher Site | Google Scholar
  20. C. L. Cioffi, M. Garay, J. F. Johnston et al., “Selective inhibition of A-Raf and C-Raf mRNA expression by antisense oligodeoxynucleotides in rat vascular smooth muscle cells: role of A-Raf and C-Raf in serum-induced proliferation,” Molecular Pharmacology, vol. 51, no. 3, pp. 383–389, 1997. View at: Google Scholar
  21. W. Guo, B. Hao, Q. Wang, Y. Lu, and J. Yue, “Requirement of B-Raf, C-Raf, and A-Raf for the growth and survival of mouse embryonic stem cells,” Experimental Cell Research, vol. 319, no. 18, pp. 2801–2811, 2013. View at: Publisher Site | Google Scholar
  22. K. Roovers and R. K. Assoian, “Integrating the MAP kinase signal into the G1 phase cell cycle machinery,” BioEssays, vol. 22, no. 9, pp. 818–826, 2000. View at: Publisher Site | Google Scholar
  23. D. Xiao, D. Chinnappan, R. Pestell, C. Albanese, and H. C. Weber, “Bombesin regulates cyclin D1 expression through the early growth response protein Egr-1 in prostate cancer cells,” Cancer Research, vol. 65, no. 21, pp. 9934–9942, 2005. View at: Publisher Site | Google Scholar
  24. J. Zhang, H. Q. Mo, F. J. Tian et al., “EIF5A1 promotes trophoblast migration and invasion via ARAF-mediated activation of the integrin/ERK signaling pathway,” Cell Death & Disease, vol. 9, no. 9, p. 926, 2018. View at: Publisher Site | Google Scholar

Copyright © 2020 Weiguo Lin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.