Section 3. Medical science
https://doi.org/10.29013/AJT-20-5.6-14-21
Hruzevskyi Olexander, docent, department of microbiology, virology and immunology
Odessa National Medical University E-mail: [email protected]
A STATUS OF CELLULAR IMMUNITY IN BACTERIAL DYSBIOSIS AND BACTERIAL VAGINOSIS
Abstruct. Systemic immunodeficiency contributes to the development of bacterial dysbiosis and bacterial vaginosis (BV), therefore, it is relevant to study the systemic response of the immune system, particularly its cellular component, and to search for informative criteria of the pathological process development.
Purpose of the study was to determine the status of cellular immunity by populations of T-lymphocytes in normocenosis, bacterial dysbiosis and BV.
As result in dysbiosis and BV, CD3 and CD4 T-cell immunodeficiency was identified. IRI is a diagnostic factor both for the formation of dysbiosis and for the diagnostics of BV. The model proposed allows us to predict the probability of development of BV with a high level of confidence.
Keywords: bacterial vaginosis, T-cell immunodeficiency.
Bacterial vaginosis (BV) is a common but poorly BV is not usually associated with the redness, studied disease of the vagina, which promotes the de- swelling, or pain observed in classical inflammation, velopment of local and systemic immunodeficiency so it is called "vaginosis", not "vaginitis" [8, 555-63]. and facilitates the transmission of human immunode- However, BV is related to "subclinical" inflammation ficiency virus (HIV) [1, 219-228], papill omavirus and of the genitals, which is determined by the hyperre-cervical cancer [2, 9-18; 3, 553-558]. The problem of active response of the immune system [9, 965-76]. BV is very important for global health, because vaginal According to [10, 156-62], dysbiosis and BV are the microbiota with high levels ofstreptococcus or entero- main factors contributing to the persistent inflamma-bacteria, vaginal candidiasis and trichomoniasis causes tion of the genital organs and facilitating HIV infec-an increased risk of pelvic inflammatory disease, pre- tion among such women.
term birth and infections for mother and newly born It is the lack of local immunity that is an impor-
child [4, 859-864]. Over the last decade, an incidence tant factor for the occurrence of BV [11, 1399-1405; of BV has doubled and ranges from 26% to 40-45% 12, 103-7]. The reason may be the ability of bacteria [5, 472; 6]. For example, in the United States BV affects forming symbiote in BV to form an active biofilm, 29% ofwomen, and in sub-Saharan Africa, where HIV which causes suppression of the immune system, is wide spread - 52% ofwomen [7, 505-23]. chronicity of the process, and resistance to antibi-
otic therapy [13]. In general, dysbiosis and BV are characterized by a decrease of systemic and local inflammatory response [14, 1-5], which, as found in the study [15, 481-7], correlates with an increase of Gardnerella vaginalis and Mycoplasma hominis.
Therefore, according to the foregoing, it is relevant to study the systemic response of the immune system, in particular of its cellular component, and to search for informative indicators of the pathological process.
Purpose of the study was to determine the status of cellular immunity by populations of T-lympho-cytes in normocenosis, bacterial dysbiosis and BV.
Material and methods
This study examined 298 women aged from 16 to 64 who saw gynecologist for a preventive examination or with complaints of genital discomfort of various degrees of manifestation. Subsequent observation excluded patients who had at least one of the definitely pathogenic microorganisms (Trichomonas vaginalis, Neisseria gonorrhoeae, Chlamydia trachomatis and Herpes Simplex Virus 1, 2) in the material taken. Presence in the smear of more than 15-20 leukocytes, which indicated of an inflammatory reaction, was also the reason for exclusion from the number of patients.
During the examination, scraping of epithelium from the posterolateral vaginal paries was made using a urogenital probe. Molecular-genetic studies were performed using polymerase chain reaction (PCR) method. DNA was extracted using a kit of reagents "Proba-GS" ("DNK-Technologiia" LLC, RF). Amplification of tubes with the reaction mixture was performed in the amplifier "DTLite" ("DNK-Technologiia" LLC, RF) using the amplification programs recommended by the manufacturer of the reagent kit. Investigation ofvaginal biocenosis status was performed using the real-time PCR test system "Femoflor 16", which allowed to quantify the biota [16, 30] by the following indicators: total bacterial mass (TBM), normobiota (Lactobacillus spp.; NB), obligate anaerobes (Atopobium vaginalis, Eubacte-
rium spp., Gardnerella vaginalis, Prevotella bivia, Por-phyromonas spp., Lachnobacterium spp., Clostridium spp., Megasphaera spp., Veilonella spp., Dialister spp., Mobilunory spp., Mobiluncus spp. Peptostreptococ spp., Sneathia spp., Leptotrihia spp., Fusobacterium spp.), Facultative anaerobes (Enterobacteriaceae spp., Staphylococcus spp., Streptococcus spp.), mycoplasmas (Ureaplasma urealyticum + parvum, Mycoplasma hominis + genitalium) and yeast-like fungi (Candida spp.).
Criterion for the distribution of patients into groups was an index of conditionally pathogenic microflora (ICPM), which was calculated as the difference between the sum of all conditionally pathogenic microorganisms and the number of lactobacilli in lg GE/sample. In normocenosis ICPM was lower than -3 lg GE/sample (first group, n = 53); in grade I dysbiosis it was from -3 to -1 lg GE/sample (second group, n = 128); and in grade II dysbiosis (BV) it was more than -1 lg GE/sample (third group, n = 117) [17, 54-7; 18, 36-41]. In addition, groups with dysbiosis were subdivided into subgroups by normobiota indicator (NBI), which was calculated as the difference between TBM and lactobacilli (in lg GE/sample). In the 2nd group there are three subgroups: 1st subgroup - with NBI < 0.3 lg GE/sample (n = 23), 2nd subgroup - with NBI from 0.3 to 1.0 lg GE/sample (n = 83), and 3rd subgroup - with NBI > 1 lg GE/sample (n = 22). In the 3rd group there are two subgroups: 1st subgroup - with NBI < 1 lg GE/sample (n = 34), and 2nd subgroup - with NBI > 1 lg GE/ sample (n = 83). The maximum degree of dysbiosis was recognized in the 2nd subgroup of the 3rd group, which corresponded to the condition of BV [16].
Subpopulations of T-lymphocytes in the blood were quantified using the method of their visualization in the rosette formation reaction with monoclonal antibodies to T-lymphocytes (CD3), T-lym-phocytes helper (CD4), and T-lymphocytes killers/ suppressor (CD8) using erythrocyte diagnostic tubes manufactured by "Granum" NPP (Ukraine). Lymphocyte suspension was separated on a density gradient d = 1.077 ("Granum" NPP, Ukraine). The
ratio CD4/CD8 (immunoreactivity index - IRI) was calculated.
For descriptive statistics, arithmetic mean value (M) and standard error of mean (m), median (Me), 1st and 3rd quartiles (Q1; Q3) were used. Variation series were compared using the Student's test (t), single- and multi-factor analysis of variance (criterion F). Paired independent data samples were compared using Mann-Whitney (U) test. Influence of factor variables on the dependent indicators was investigated using linear regression analysis. Regres-
Table 1.- Indicators of
sion coefficients probability of their contrast to the main hypothesis, Wald-statistics value and the maximum probability coefficient for nonlinear models were calculated. Operational characteristics were evaluated using ROC-diagrams. Significance of all the variations was taken at p <0.05. Statistica 10 software package (StatSoft, Inc., USA) was used for statistical processing of the obtained data.
Results and Discussion. Indicators of cellular immunity in patients with normo- and dysbiosis are given in (Table 1).
cellular immunity (M ± m)
Group, subgroup Leucocytes,G/l CD3, G/l CD4, G/l CD8, G/l IPI, rel. units
1st (normocenosis). n=53 6.104±0.162 1.103±0.042 0.629±0.018 0.408±0.009 1.592±0.060
2nd (grade I dysbiosis). n=128 1st, n=23 6.478±0.296 1.147±0.063 0.681±0.037 0.428±0.015 1.631±0.104
2nd, n=83 7.142±0.180 1.165±0.031 0.690±0.019 0.414±0.006 1.695±0.054
3rd, n=22 7.476±0.337 1.159±0.058 0.469±0.020 0.422±0.015 1.140±0.061
3-H (grade II dysbiosis. n=117 1st, n=34 8.658±0.258 0.595±0.014 0.442±0.011 0.414±0.011 1.094±0.041
2nd, n=83 10.122±0.183 0.413±0.007 0.353±0.007 0.508±0.009 0.709±0.017
Statistical procedure of comparison of the results (p)
p(MW)1 0.282 0.527 0.190 0.218 0.709
p(MW)2 0.001 0.276 0.048 0.758 0.309
p(MW)3 0.001 0.446 <0.001 0.496 <0.001
p(MW)4 <0.001 <0.001 <0.001 0.811 <0.001
p(MW)5 <0.001 <0.001 <0.001 <0.001 <0.001
F 58.834 126.594 73.392 23.089 72.095
P <0.001 <0.001 <0.001 <0.001 <0.001
Notes: probability of discrepancies between the corresponding indicators in 1st and 2nd groups using MannWhitney test: p(MW)1 - in the 1st subgroup of the 2nd group, p(MW)2 - in the 2nd subgroup of the 2nd group,
nd
p(MW)3 - in the 3rd subgroup of the 2nd group, p(MW)4 - in the 1st subgroup of the 3rd group, p(MW)5 - in the 2 subgroup of the 3rd group; F - result, and p - probability of the dispersion analysis of variance of differences of the corresponding indicators among the subgroups
tory processes and/or infectious diseases, we considered these values to be the normal variant.
As dysbiosis progressed, progression ofleukocy-tosis was also observed (Fig. 1). Thus, the content of leukocytes in comparison with the 1st group was in-
Number of leukocytes in the blood under condition of normocenosis was slightly higher than the normal one (4.4-5.5 G/l [19, 960]). As the patients in 1st group did not experience inflamma-
creased by 1.2 times in the 2nd group (p = 0.001), by opment of BV (2nd subgroup of the 3rd group), leuko-
1.4 times in the 1st subgroup of the 3rd group and by cytosis acquired the maximum value and exceeded
1.7 times in 2nd subgroup of the 3rd group (p <0.001 such values in grades I and II dysbiosis (by 1.4 times
in both cases). Therefore, under conditions of devel- and 1.2 times respectively; p <0.05 in both cases).
G/l
Figure 1. The absolute number of leukocytes in the blood (G/l) depending on the grade of dysbiosis; statistical significance of the differences is shown in Table 1
Figure 2. Content of T-lymphocytes in the blood by fractions (CD3, CD4 and CD8) (G/l; right axis) and IRI value (rel. units; left axis) depending on the grade of dysbiosis; statistical significance of the differences is shown in Table 1
CD3 is a multiprotein membrane complex that lymphocytes. As we can see from Table 1, in grade is present on the surface of T-lymphocytes and I dysbiosis (2nd group) in comparison with the 1st
is the major coreceptor of T-cell receptor (TCR) group (normocenosis) CD3 content in the blood did
[20, 576]. Due to this, CD3 is a major marker of T- not actually change, whereas in grade II dysbiosis (3rd
group) it was reduced - by 1.8 times in the 1st subgroup and by 2.7 times in the 2nd subgroup (p <0.001 in both cases).
Such dynamics (Fig. 2) indicated the development of deep grade dysbiosis of CD3-lymphocy-topenia. Accentuated decrease in the number of CD3-lymphocytes in the blood was a characteristic of BV, when it was statistically significantly lower than in the case of grade II dysbiosis (by 1.4 times; p <0.001). It indicated the progression of T-cell immunodeficiency with increasing severity of dysbiosis.
CD4 is a monomeric transmembrane glycoprotein that belongs to the immunoglobulin superfamily and is a marker of T-helper cells [20]. It is expressed on thymocytes (80-90%), mature T-lymphocytes (65% of T-helper cells), monocytes, macrophages, Langerhans cells, dendritic cells.
In our study, we observed a significant decrease in the number of CD4-lymphocytes in the blood in severe dysbiosis, namely in the 3rd subgroup of the 2nd group and in the 3rd group, by 1.3-1.8 times (p <0.001 for all these observations). To a maximum extent (by 1.8 times; p <0.001), CD4-lymphocyte content was reduced in BV, when it was significantly lower than in dysbiosis (by 1.3 times; p <0.05). Also, inhibition of T-lymphocytic component in BV was identified [21], which the authors associate with a decrease in predominantly gamma-delta-TcR T-lym-phocytes, which recognize non-peptide microbial antigens, similar to NK-cells.
CD8 is transmembrane glycoprotein, which is a coreceptor of T-cell receptors (TCR). Like TCR, CD8 has the ability to bind to a class I major his-tocompatibility complex molecule [20]. Most of CD8 is represented by peripheral T-lymphocytes as a disulfide-linked alpha-chain homodimer, which is a marker of subpopulation of cytotoxic T-lymphocytes and T-suppressors. On the mature T-cells, either CD8 or CD4 is expressed.
In our study, in normocenosis and dysbiosis content of CD8-lymphocyte in the blood was almost the same (see Table 1 and Figure 2). 2nd subgroup
of the 3rd group was the exception - patients with BV showed a significant increase in content of CD8 lymphocytes in their blood (by 1.2 times compared to those in normocenosis; p <0.001).
Presence of pre-existing trends showing changes in CD3, CD4, CD8 lymphocyte content in the blood, and IRI values in the development of dysbiosis and BV, was confirmed by analysis of variance estimates among the subgroups (see Table 1): F value was from 23.1 to 126.6 (p <0.001 for all the indicators).
Thus, in dysbiosis, there was a gradual, associated with an increase of severity, increase of CD3- and CD4-lymphocytes in the blood, which was expressed to a maximum extent in case of BV. This allowed us to establish the progression of T-cell immunodeficiency in the development of dysbiosis. In contrast, the response of cytotoxic T-lymphocytes and T-suppressors (increase in CD8 content in the blood) was observed only in case of BV, which evidenced the development of compensatory activation of the cellular component of immune system.
The revealed pattern was confirmed and quantified by the dynamics of IRI (see Table 1 and Figure 2). In the 3rd subgroup of the 2nd group and in the 3rd group, IRI was significantly reduced in comparison with normocenosis (by 1.4-2.2 times; p < 0.001). Moreover, in case of BV, IRI was reduced to the maximum extent - by 2.2 times in comparison with normocenosis, and in by 1.5-1.6 times in comparison with the indicators in dysbiosis (p < 0.001 in all cases).
The analysis ofvariance (Table 2) confirmed the revealed differences of IRI among the groups of patients with normocenosis (1st group) in comparison with dysbiosis (2nd and 3rd groups). When analyzing the quantitative data, it was found that the value of IRI can be considered a diagnostic marker of development of dysbiosis. IRI is less than Q1 = 1.340 rel. units was the limit below which the dysbiosis could be statistically diagnosed (F = 20.47; p < 0.001).
Table 2.- Values of IRI in normocenosis and dysbiosis
Indicators Groups F P
1st group (normocenosis), n=53 2nd and 3rd groups (grade I and II dysbiosis), n=245
M±m 1.592±0.060 1.222±0.036 20.47 < 0.001
Me (Q1; Q3) 1.543 (1.340; 1.805) 1.081 (0.758; 1.335)
Notes: M ± m - mean value and standard error of the mean; Me (Q1; Q3) - median, 1st and 3rd quartiles; F -Fisher's criterion for analysis of variance; p - statistical significance of differences among groups (accepted at p <0,05)
Analysis of the effect of IRI on the development in all the subgroups with dysbiosis with the 2nd sub-of BV compared to dysbiosis was performed when group of the 3rd group (Table 3). comparing the data ofvariance analysis of IRI values
Table 3.- IRI values in dysbiosis and bacterial vaginosis
Indicators Groups F p
2nd group and 1st subgroup of the 2nd group (dysbiosis), n=162 2nd subgroup of the 3rd group (BV), n=83
M±m 1.484±0.040 0.709±0.017 186.032 < 0.001
Me (Q1; Q3) 1.401 (1.081; 1.805) 0.700 (0.595; 0.800)
Notes: M ± m - mean value and standard error of the mean; Me (Q1; Q3) - median, 1st and 3rd quartiles; F -Fisher's criterion for analysis ofvariance; p - statistical significance of differences among groups (accepted at p <0,05)
It is established that Q1-Q3 intervals for IRI in ematical models of immunity disorders is justified by
groups ofpatients with dysbiosis and BV do not inter- [22]. According to [23], the need to find optimal and
sect. In this regard, IRI less than 0.800 rel. units may be simple prognostic models is an important modern
considered as a margin, below which there is a prob- problem, especially in the context of the tendency
ability of BV development (F = 186.032; p < 0.001). of BV to become chronic.
The fundamental possibility of predicting the development of dysbacteriosis and BV using math-
Table 4.- Statistical significance of indicators, logistic regression coefficients and their probability for the dependent variable PBV
Indicators Wald ^-coef. ±SEp BI±95% t p
IPI 46.764 -10.546 1.542 -(13.568-7.523) -6.826 < 0.001
Independent indicator 46.655 9.305 72.922 6.635-11.976 6.820 < 0.001
Notes: Wald - resulting value of Wald-statistics; ft - regression coefficient; ± SEft - regression coefficient error; BI ± 95% - 95% confidence interval; t - Student's coefficient; p - significance of differences compared to the main hypothesis (accepted at p <0,05)
Logistic regression analysis was performed to cal- was the resulting feature (categorical values "YES" culate the probability of BV development. IRI values and "NO" were assigned indicator values: "1" and were used as a factor trait. Presence or absence of BV "0" respectively). When developing analysis model,
245 patients with dysbiosis were engaged, 83 patients among them had BV. To evaluate the contribution of IRI values to the probability of BV prediction (PBV), analysis of the resulting Wald-statistics value was conducted. (^-coefficients of the regression equation, their standard errors, confidence intervals, and significance of differences compared to the main hypothesis were calculated (Table 4).
Calculated ^-coefficients of the regression equation, their standard errors and confidence intervals indicate of the presence of adverse relation of IRI with the probability of BV development (t = -6,826; p < 0.001 compared with the main hypothesis). Operational characteristics of the developed model, calculated using ROC-analysis, showed the satisfactory quality: AUC = 0.961 ± 0.011 (CI ± 95% 0.940-0.983; p < 0.001). Assessment of the model from the point of view of conformity with the used regression construction method also had satisfactory parameters: -2log = = 115,426 (x2 = 198,281; p < 0.001; df = 1). Accuracy of the error-free detection of BV was 91.0%.
Thus, it was determined that IRI is a diagnostic factor both for formation of dysbiosis and for diagnostics of BV. The model proposed meets the quality criteria and allows us to predict the probability of BV development with a high degree of certainty.
Conclusions
1. In dysbiosis on the background of general leukocytosis, there was a gradual, related to the progression of severity of dysbiosis, decrease in the content of blood of CD3- and CD4-lymphocytes, which was expressed to a maximum extent in BV. This allowed us to establish the progression of T-cell immunodeficiency in the development of dysbiosis.
2. The response of cytotoxic T-lymphocytes and T-suppressors (CD8) was observed only in BV. IRI was a diagnostic factor both for the formation of dysbiosis (IRI < 1.340 units) and for the diagnostics of BV (IRI < 0.800 units).
3. Regression model developed meets the required criteria and allows us to predict the probability of development of BV with an accuracy of 91.0%.
References:
1. McKinnon L. R., Achilles S. L., Bradshaw C. S., Burgener A., Crucitti T., Fredricks D. N., Jaspan H. B., Kaul R. et al. The evolving facets of bacterial vaginosis: implications for HIV transmission. AIDS Res Hum Retroviruses. 2019. Mar,- 35(3).- P. 219-228. Doi: 10.1089/AID.2018.0304.
2. Brusselaers N., Shrestha S., van de Wijgert J., Verstraelen H. Vaginal dysbiosis and the risk of human papillomavirus and cervical cancer: systematic review and meta-analysis. Am J Obstet Gynecol. 2019. Jul, - 221(1).- P. 9-18.e8. Doi: 10.1016/j.ajog.2018.12.011.
3. Reid G. Is bacterial vaginosis a disease? Appl Microbiol Biotechnol. 2018. Jan,- 102(2).- P. 553-558. Doi: 10.1007/s00253-017-8659-9.
4. Van de Wijgert JHHM, Jespers V2. The global health impact of vaginal dysbiosis.Res Microbiol. 2017. Nov-Dec,- 168(9-10).- P. 859-864. Doi: 10.1016/j.resmic.2017.02.003.
5. Kira E. F. [Bacterial vaginosis].- Moscow: Medical Information Agency, 2012.- 472 p. [in Russian].
6. Mark K. S., Tenorio B., Stennett C. A., Ghanem K. G., Brotman R. M. Bacterial vaginosis diagnosis and treatment in postmenopausal women: a survey of clinician practices. Menopause. 2020. Mar, 2. Doi: 10.1097/GME.0000000000001515.
7. Kenyon C., Colebunders R., Crucitti T. The global epidemiology of bacterial vaginosis: A systematic review. Am J Obstet Gynecol 2013.- 209.- P. 505-23. Doi: 10.1016/j.ajog.2013.05.006.
8. Mitchell C., MarrazzoJ. Bacterial vaginosis and the cervicovaginal immune response. Am J Reprod Immunol. 2014.- 71.- P. 555-63. Doi: 10.1111/aji.12264
9. Anahtar M. N., Byrne E. H., Doherty K. E., Bowman B. A., Yamamoto H. S., Soumillon M., Padavattan N., Ismail N. et al. Cervicovaginal bacteria are a major modulator of host inflammatory responses in the female genital tract. Immunity. 2015. May, 19.- 42(5).- P. 965-76. Doi: 10.1016/j.immuni.2015.04.019.
10. Passmore J. A.,Jaspan H. B., Masson L. Genital inflammation, immune activation and risk ofsexual HIV acquisition. Curr Opin HIVAIDS. 2016. Mar, - 11(2). - P. 156-62. Doi: 10.1097/C0H.0000000000000232.
11. Muzny C. A., Taylor C. M., Swords W. E., Tamhane A., Chattopadhyay D., Cerca N., Schwebke J. R. An updated conceptual model on the pathogenesis of bacterial vaginosis. J Infect Dis. 2019. Sep 26.- 220 (9).- P. 1399-1405. Doi: 10.1093/infdis/jiz342.
12. Gruzevskyy O. A. [Colonization resistance in vaginal dysbiosis: the state of humoral and cellular links]. Bul marine med. 2017.- 4(77).- P. 103-7. [in Russian].
13. Vestby L. K., Gronseth T., Simm R., Nesse L. L. Bacterial biofilm and its role in the pathogenesis of disease. Antibiotics (Basel). 2020. Feb, - 3; 9(2). pii: E59. doi: 10.3390/antibiotics9020059.
14. Muzny C. A., Schwebke J. R. Pathogenesis ofbacterial vaginosis: discussion of current hypotheses. J Infect Dis. 2016. Aug, 15; 214 Suppl 1.- P. 1-5. Doi: 10.1093/infdis/jiw121.
15. Cox C., Watt A. P., Mc Kenna J. P., Coyle P. V. Mycoplasma hominis and Gardnerella vaginalis display a significant synergistic relationship in bacterial vaginosis. Eur J Clin Microbiol Infect Dis. 2016. MapT.-35 (3).- P. 481-7. Doi: 10.1007 / s10096-015-2564-x.
16. Lipova E. V., Boldyreva M. N., Trofimov D. Yu., Vitvitskaya Yu.G. [Femoflor. Urogenital infections caused by conditionally biota in women of reproductive age (clinical and laboratory diagnostics). Manual for doctors].- Moscow: DNA technology. 2015.- 30 p. [in Russian].
17. Hruzevskyi O. A., Vladymirova M. P. [Results of a complex bacteriological study of vaginal contents under the conditions of bacterial vaginosis]. Ach biol and med. 2014.- 2.- P. 54-7. [in Ukrainian].
18. Hruzevskyi O. A. [Normocenosis of the vagina: qualitative and quantitative characteristics] Odessa Med J. 2015.- 1(147).- P. 36-41. [in Ukrainian].
19. Tits N. U. Encyclopedia of clinical laboratory tests.- Moscow: Labinform. 1997.- 960 p. [in Russian].
20. Delves P. J., Martin S. J., Burton D. R., Roitt I. M. Roitt's Essential Immunology, 13th Edition. 2016. Wiley-Blackwell,- 576 p.
21. Alcaide M. L., Strbo N., Romero L., Jones D. L., Rodriguez V.J., Arheart K., Martinez O., Bolivar H., Podack E. R., Fischl M. A. Bacterial vaginosis is associated with loss of gamma delta T cells in the female reproductive tract in women in the Miami Women Interagency HIV study (WIHS): A cross sectional study. PLoS One. 2016. Apr 14.- 11(4): e0153045. Doi: 10.1371/journal.pone.0153045.
22. Fichorova R. N., Morrison C. S., Chen P. L., Yamamoto H. S., Govender Y., Junaid D., Ryan S., Kwok C., Chipato T., Salata R. A., Doncel G. F. Aberrant cervical innate immunity predicts onset of dysbiosis and sexually transmitted infections in women of reproductive age. PLoS One. 2020. Jan 8.- 15(1).-e0224359. Doi: 10.1371/journal.pone.0224359.
23. Sobel J. D., Kaur N., Woznicki N. A., Boikov D., Aguin T., Gill G., Akins R. A. Prognostic indicators of recurrence of bacterial vaginosis. J Clin Microbiol. 2019. Apr 26. - 57(5).- pii: e00227-19. Doi: 10.1128/ JCM.00227-19.