Whole-transcriptome sequencing of phagocytes reveals a ceRNA network contributing to natural resistance to tuberculosis infection

Tuberculosis (TB) is a major fatal infectious disease globally, exhibiting high morbidity rates and impacting public health and other socio-economic factors. However, some individuals are resistant to TB infection and are referred to as “Resisters”. Resisters remain uninfected even after exposure to high load of Mycobacterium tuberculosis (Mtb). To delineate this further, this study aimed to investigate the factors and mechanisms influencing the Mtb resistance phenotype. We assayed the phagocytic capacity of peripheral blood mononuclear cells (PBMCs) collected from Resisters, patients with latent TB infection (LTBI), and patients with active TB (ATB), following infection with fluorescent Mycobacterium bovis Bacillus Calmette-Guérin (BCG). Phagocytosis was stronger in PBMCs from ATB patients, and comparable in LTBI patients and Resisters. Subsequently, phagocytes were isolated and subjected to whole transcriptome sequencing and small RNA sequencing to analyze transcriptional expression profiles and identify potential targets associated with the resistance phenotype. The results revealed that a total of 277 mRNAs, 589 long non-coding RNAs, 523 circular RNAs, and 35 microRNAs were differentially expressed in Resisters and LTBI patients. Further, the endogenous competitive RNA (ceRNA) network was constructed from differentially expressed genes after screening. Bioinformatics, statistical analysis, and quantitative real-time polymerase chain reaction were used for the identification and validation of potential crucial targets in the ceRNA network. As a result, we obtained a ceRNA network that contributes to the resistance phenotype. TCONS_00034796-F3, ENST00000629441-DDX43, hsa-ATAD3A_0003-CYP17A1, and XR_932996.2-CERS1 may be crucial association pairs for resistance to TB infection. Overall, this study demonstrated that the phagocytic capacity of PBMCs was not a determinant of the resistance phenotype and that some non-coding RNAs could be involved in the natural resistance to TB infection through a ceRNA mechanism. •The PBMC phagocytosis of resisters is weaker than that of ATB patients and comparable to that of LTBI patients.•We obtained a ceRNA network involved in immune metabolism that contributes to the resistance phenotype.•The phagocytic capacity of PBMCs is not a determinant of the resistance phenotype, and some non-coding RNAs are involved in resistance to TB infection through the ceRNA mechanism.