We generated linear regression models that incorporated total bacterial load, Nugent BV status, and vaginal CST, either alone or in combination with one another, to predict vaginal levels of sE-cad, IL-1α, and IP-10 since these immune factors were most strongly associated with vaginal microbiota composition and total bacterial load. Here, IL1A is linked to bacterial vaginosis.