Skip to main content
  • Letter to the Editor
  • Open access
  • Published:

Response to comments on our article (Yin YL et al., Parasit Vectors, 10.1186/s13071-021-04739-w) by Yuqing Wang and colleagues


This letter responds to comments on our article (Yin YL et al., Parasit Vectors, 10.1186/s13071-021-04739-w) by Yuqing Wang and colleagues, who wrote a letter entitled “Microarray analysis of circular RNAs in HCT-8 cells infected with Cryptosporidium parvum” and discussed statistical procedures for microarray analysis during C. parvum infection. To further confirm our data, in this letter, a common R package for analyses of differentially expressed genes, namely DESeq2, with Benjamini-Hochberg correction, was used to analyze our microarray data and identified 26 significantly differentially expressed circRNAs using adjusted P value < 0.05 and | Log2 (fold change [FC]) | ≥ 1.0, including our circRNA ciRS-7 of interest. Therefore, the protocol for selecting circRNAs of interest for further study in our article is acceptable and did not affect the subsequent scientific findings in our article.

To the Editor

Cryptosporidiosis, caused by Cryptosporidium spp., is an important zoonotic parasitic disease, seriously threatening the health of humans and many animals [1, 2]. Considering that the severity of cryptosporidiosis is closely associated with host status, especially immunity, understanding the host response to infection is critical to effectively develop well-directed control strategies against cryptosporidiosis [3].

Microarray analysis is a sensitive tool to accurately investigate differentially expressed circRNAs during many disease processes [4,5,6,7,8]. It is based on hybridization and fluorescence detection and uses specific probes targeting the back-spliced junction of each circRNA. CapitalBio Technology Human CircRNA Array, which was used in our study, also has its own analyzing software, GeneSpring, and GeneSpring supports normalization processing and analysis of differential circRNAs of original data from microarray analysis. Therefore, we used this software to process our microarray data. Due to the generally mild impact of Cryptosporidium infection on host cell gene transcription compared to other pathogens [9,10,11,12], to obtain as many differentially expressed genes as possible, we chose fold changes > 2 combined with P ≤ 0.05 as the threshold to preliminarily screen differentially expressed circRNAs in the microarray analysis in our article [3]. The criterion of P <0.05 or P ≤ 0.05 was also used by recent studies [13,14,15,16]. From the subsequent verification of the experimental results, circRNAs of interest obtained using the threshold of P ≤ 0.05 are stable and accurate. Certainly, the unadjusted P-value will introduce false positives during multiple comparisons, as stated by Wang et al. [17], in their Letter to the Editor published commenting on our article. To reduce this possibility, expression of the circRNA of interest, ciRS-7, was further verified by using real-time PCR in multiple repetitions in our article. This protocol is also a standard method to select genes of interest for subsequent experiments [15, 16, 18, 19].

Nowadays, several software packages and models have been used to define significantly differentially expressed (DE) genes during multiple comparisons. Three R packages, namely limma, DESeq2 and edge R, are commonly used to analyze DE genes for microarray and RNA-seq data. Of these, both limma and DESeq2 are quite reliable and are not much different for multiple comparisons of these data, with > 90% of genes detected overlapping between the two methods [20]. It should be noted that limma can find the accurate DE genes better but obtains fewer significant DE genes because of its rigorous screening criteria, while DESeq2 is more suitable when more possible candidate DE genes are expected. However, both methods are preliminary screening tools that will inevitably introduce false-positive/negative results to a certain extent. Therefore, the convincing expression of candidate DE genes should be further confirmed by other more precise methods with multiple repetitions, such as the real-time PCR used in our article. Since Cryptosporidium infection has a mild impact on host cell gene transcription, here we further analyzed our microarray data by using DESeq2 for multiple comparisons using Benjamini-Hochberg correction and identified a total of 26 significantly DE circRNAs (Additional file 1: Table S1) using the standard of adjusted P value < 0.05 and | Log2 (fold change [FC]) | ≥ 1.0, including our circRNA of interest, ciRS-7 (P = 0.030272316 and Log2 FC = 2.497972579).

Collectively, the protocol for selecting circRNAs of interest for further study in our article is acceptable and did not affect the subsequent scientific findings. We welcome further discussions of our article.

Availability of data and materials

All datasets generated for this study are included in the article and its additional file.


  1. Zhang G, Zhang Y, Niu Z, Wang C, Xie F, Li J, et al. Cryptosporidium parvum upregulates miR-942-5p expression in HCT-8 cells via TLR2/TLR4-NF-κB signaling. Parasit Vectors. 2020;13:435.

    Article  CAS  Google Scholar 

  2. Carter BL, Chalmers RM, Davies AP. Health sequelae of human cryptosporidiosis in industrialised countries: a systematic review. Parasit Vectors. 2020;13:443.

    Article  Google Scholar 

  3. Yin YL, Liu TL, Yao Q, Wang YX, Wu XM, Wang XT, et al. Circular RNA ciRS-7 affects the propagation of Cryptosporidium parvum in HCT-8 cells by sponging miR-1270 to activate the NF-κB signaling pathway. Parasit Vectors. 2021;14:238.

    Article  CAS  Google Scholar 

  4. Luo S, Deng M, Xie Z, Li X, Huang G, Zhou Z. Circulating circular RNAs profiles associated with type 1 diabetes. Diabetes Metab Res Rev. 2021;37:e3394.

    CAS  PubMed  Google Scholar 

  5. Kalifu B, Maitiseyiti A, Ge X, Chen X, Meng Y. Expression profile of circular RNAs in cystic echinococcosis pericystic tissue. J Clin Lab Anal. 2021;35:e23687.

    Article  CAS  Google Scholar 

  6. Li Y, Zhang Y, Zhang S, Huang D, Li B, Liang G, et al. circRNA circARNT2 suppressed the sensitivity of hepatocellular carcinoma cells to cisplatin by targeting the miR-155-5p/PDK1 Axis. Mol Ther Nucleic Acids. 2020;23:244–54.

    Article  Google Scholar 

  7. Yi F, Xin L, Feng L. Potential mechanism of circRNA_000585 in cholangiocarcinoma. J Int Med Res. 2021;49:3000605211024501.

    CAS  PubMed  Google Scholar 

  8. Peng D, Hou ZL, Zhang HX, Zhang S, Zhang SM, Lin RY, et al. Microarray expression profile and analysis of circular RNA regulatory network in pulmonary thromboembolism. Int J Gen Med. 2021;14:1239–49.

    Article  Google Scholar 

  9. Deng M, Lancto CA, Abrahamsen MS. Cryptosporidium parvum regulation of human epithelial cell gene expression. Int J Parasitol. 2004;3:73–82.

    Article  Google Scholar 

  10. Yang YL, Serrano MG, Sheoran AS, Manque PA, Buck GA, Widmer G. Over-expression and localization of a host protein on the membrane of Cryptosporidium parvum infected epithelial cells. Mol Biochem Parasitol. 2009;168:95–101.

    Article  CAS  Google Scholar 

  11. Zhou R, Hu G, Liu J, Gong AY, Drescher KM, Chen XM. NF-kappaB p65-dependent transactivation of miRNA genes following Cryptosporidium parvum infection stimulates epithelial cell immune responses. PLoS Pathog. 2009;5:e1000681.

    Article  Google Scholar 

  12. Zhao GH, Gong AY, Wang Y, Zhang XT, Li M, Mathy NW, et al. Nuclear delivery of parasite Cdg2_FLc_0220 RNA transcript to epithelial cells during Cryptosporidium parvum infection modulates host gene transcription. Vet Parasitol. 2018;251:27–33.

    Article  CAS  Google Scholar 

  13. Zhi F, Ding Y, Wang R, Yang Y, Luo K, Hua F. Exosomal hsa_circ_0006859 is a potential biomarker for postmenopausal osteoporosis and enhances adipogenic versus osteogenic differentiation in human bone marrow mesenchymal stem cells by sponging miR-431-5p. Stem Cell Res Ther. 2021;12:157.

    Article  Google Scholar 

  14. Gao YJ, Zhang RJ, Liu Q, Sun SG, Qi MY, Wang Y, et al. Functional predication of differentially expressed circRNAs/lncRNAs in the prefrontal cortex of Nrf2-knockout mice. Aging (Albany NY). 2021;13:8797–816.

    Article  Google Scholar 

  15. Li L, Wei H, Zhang H, Xu F, Che G. Circ_100565 promotes proliferation, migration and invasion in non-small cell lung cancer through upregulating HMGA2 via sponging miR-506-3p. Cancer Cell Int. 2020;20:160.

    Article  CAS  Google Scholar 

  16. Li Y, Shi P, Zheng T, Ying Z, Jiang D. Circular RNA hsa_circ_0131242 promotes triple-negative breast cancer progression by sponging hsa-miR-2682. Onco Targets Ther. 2020;13:4791–8.

    Article  CAS  Google Scholar 

  17. Wang Y, Zhao H, Zhang Y, Yan L. Microarray analysis of circular RNAs in HCT-8 cells infected with Cryptosporidium parvum. Parasit Vectors. 2021;12:99.

    CAS  Google Scholar 

  18. Xiang Q, Kang L, Wang J, Liao Z, Song Y, Zhao K, et al. CircRNA-CIDN mitigated compression loading-induced damage in human nucleus pulposus cells via miR-34a-5p/SIRT1 axis. EBioMedicine. 2020;53:102679.

    Article  Google Scholar 

  19. Shi Y, Fang N, Li Y, Guo Z, Jiang W, He Y, et al. Circular RNA LPAR3 sponges microRNA-198 to facilitate esophageal cancer migration, invasion, and metastasis. Cancer Sci. 2020;111:2824–36.

    Article  CAS  Google Scholar 

  20. Tong Y. The comparison of limma and DESeq2 in gene analysis. E3S Web Conf. 2021;271:03058.

    Article  Google Scholar 

Download references


This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Author information

Authors and Affiliations



GHZ and YLY conceived and prepared the first draft of the manuscript. GHZ, YLY and XY critically reviewed the draft. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Guang-Hui Zhao.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1. Significantly differentially expressed circRNAs between C. parvum-infected and -non-infected HCT-8 cells.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yin, YL., Yang, X. & Zhao, GH. Response to comments on our article (Yin YL et al., Parasit Vectors, 10.1186/s13071-021-04739-w) by Yuqing Wang and colleagues. Parasites Vectors 14, 484 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: