Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Host-Microbe Biology

Host Genetic Factors Associated with Vaginal Microbiome Composition in Kenyan Women

Supriya D. Mehta, Drew R. Nannini, Fredrick Otieno, Stefan J. Green, Walter Agingu, Alan Landay, Yinan Zheng, Lifang Hou
Anthony Fodor, Editor
Supriya D. Mehta
aDivision of Epidemiology & Biostatistics, University of Illinois at Chicago School of Public Health, Chicago, Illinois, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Supriya D. Mehta
Drew R. Nannini
bCenter for Global Oncology, Institute of Global Health, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
cDepartment of Preventive Medicine, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fredrick Otieno
dNyanza Reproductive Health Society, Kisumu, Kenya
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stefan J. Green
eGenome Research Core, University of Illinois at Chicago, College of Medicine, Chicago, Illinois, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Stefan J. Green
Walter Agingu
dNyanza Reproductive Health Society, Kisumu, Kenya
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alan Landay
fDepartment of Internal Medicine, Rush University College of Medicine, Chicago, Illinois, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yinan Zheng
bCenter for Global Oncology, Institute of Global Health, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
cDepartment of Preventive Medicine, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lifang Hou
bCenter for Global Oncology, Institute of Global Health, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
cDepartment of Preventive Medicine, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anthony Fodor
University of North Carolina at Charlotte
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSystems.00502-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Bacterial vaginosis (BV) affects 20% of women worldwide and is associated with adverse reproductive health outcomes and increased risk for HIV. Typically, BV represents a shift in the vaginal microbiome from one that is dominated by Lactobacillus to one that is diverse. Persistent racial differences in BV and diverse vaginal microbiome composition overlap with racial disparities in risks for HIV and sexually transmitted infection, especially among women of African descent. Risk factors for BV and nonoptimal vaginal microbiome include sexual practices, yet racial differences persist when adjusted for behavioral factors, suggesting a host genetic component. Here, we perform a genome-wide association study on vaginal microbiome traits in Kenyan women. Linear regression and logistic regression were performed, adjusting for age and principal components of genetic ancestry, to evaluate the association between Lactobacillus crispatus, Lactobacillus iners, Gardnerella vaginalis, Shannon diversity index, and community state type (CST) with host genetic single nucleotide polymorphisms (SNPs). We identified novel genomic loci associated with the vaginal microbiome traits, though no SNP reached genome-wide significance. During pathway enrichment analysis, Toll-like receptors (TLRs), cytokine production, and other components of innate immune response were associated with L. crispatus, L. iners, and CST. Multiple previously reported genomic loci were replicated, including IL-8 (Shannon, CST), TIRAP (L. iners, Shannon), TLR2 (Shannon, CST), MBL2 (L. iners, G. vaginalis, CST), and MYD88 (L. iners, Shannon). These genetic associations suggest a role for the innate immune system and cell signaling in vaginal microbiome composition and susceptibility to nonoptimal vaginal microbiome.

IMPORTANCE Globally, bacterial vaginosis (BV) is a common condition in women. BV is associated with poorer reproductive health outcomes and HIV risk. Typically, BV represents a shift in the vaginal microbiome from one that is dominated by Lactobacillus to one that is diverse. Despite many women having similar exposures, the prevalence of BV and nonoptimal vaginal microbiome is increased for women of African descent, suggesting a possible role for host genetics. We conducted a genome-wide association study of important vaginal microbiome traits in Kenyan women. We identified novel genetic loci and biological pathways related to mucosal immunity, cell signaling, and infection that were associated with vaginal microbiome traits; we replicated previously reported loci associated with mucosal immune response. These results provide insight into potential host genetic influences on vaginal microbiome composition and can guide larger longitudinal studies, with genetic and functional comparison across microbiome sites within individuals and across populations.

INTRODUCTION

Bacterial vaginosis (BV) is a condition of clinical and public health significance. BV is prevalent in 20 to 50% of women in sub-Saharan Africa (1). In a meta-analysis, Atashili et al. estimated the relative risk of HIV acquisition to be 1.6 times higher for women with BV, equating to a population attributable fraction of 15% due to the high prevalence in the population (2). BV is also associated with an increased likelihood of sexually transmitted infections (STIs) and adverse pregnancy outcomes (3, 4). BV represents a polymicrobial shift in the vaginal microbiome, often from a community dominated by Lactobacillus to one that is diverse, with multiple species of anaerobic bacteria (5). In particular, Lactobacillus crispatus enrichment has been shown to be protective against BV, HIV, and STIs (6).

BV is considered a sexually enhanced infection, and individuals with new or multiple sex partners, those who engage in unprotected vaginal sex, or those whose sex partner is uncircumcised are at an increased risk of BV (3, 7–9). In addition to sexual risk factors, nonsexual risk factors that are associated with BV include intravaginal and vaginal hygiene practices (10), cigarette smoking (11), and race/ethnicity (12). In numerous studies, women of African descent have increased risk of BV and nonoptimal vaginal microbiome composition (13). These associations with race persist even when controlling for sociodemographics and sexual practices (14, 15), leading to the hypothesis that genetic factors may influence the vaginal microbiome composition and subsequently BV (16).

To date, a limited number of studies have examined the host genetic contribution to the vaginal flora and BV. Primarily these studies have focused on host differences in candidate genes responsible for inflammatory mucosal immune response in other infectious conditions (16, 17). As reviewed by Turovskiy et al., studies have targeted single nucleotide polymorphisms (SNPs) in (i) Toll-like receptor (TLR) genes because of their role in recognizing potential threats and initiating inflammatory responses and activating other immune cells and (ii) proinflammatory and inflammatory cytokines (e.g., interleukin 1β [IL-1β], tumor necrosis factor alpha [TNF-α]) (17). Most of these studies, however, evaluated SNPs in relation to BV or to specific taxa associated with BV (e.g., Gardnerella vaginalis, Atopobium vaginae), and none of the studies measured the vaginal microflora or host genetic contribution comprehensively, such as with cultivation-independent molecular microbial characterization or a genome-wide association study (GWAS), respectively.

To address this gap, we conducted GWAS on vaginal microbiome traits among native Kenyan women. We hypothesized that by conducting GWAS, we would identify novel genomic loci associated with vaginal microbiome traits, furthering our understanding of the underlying genetic factors. Identifying genetic variants may aid in elucidating the biological mechanisms associated with these complex vaginal microbiome traits and associated diseases. Complementary to GWAS, we conducted replication analyses of SNPs previously reported to be associated with vaginal microflora or BV.

RESULTS

Study sample.Among 171 women included in this analysis, the prevalence of BV at baseline was 22% (Table 1). More than half (53%) of the women were herpes simplex virus 2 (HSV-2) seropositive, and 10% were HIV positive, which is in keeping with the prevalence of HIV among women in this age range in the general population in this region of western Kenya (18). Lactobacillus crispatus was present in 25% of samples (28% mean relative abundance), Lactobacillus iners was present in 83% of samples (mean relative abundance 45%), and Gardnerella vaginalis was present in 75% of samples (25% mean relative abundance).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Baseline participant characteristics and vaginal microbiome composition overall and stratified by community state typea

Genome-wide association results.The genomic control (GC) inflation factor, λ, values were 1.01, 1.00, 1.00, 1.00, and 1.00 for L. crispatus, L. iners, G. vaginalis, the Shannon diversity index, and CST, respectively, after adjusting for age and the first three principal components. Upon visual inspection of the quantile-quantile (Q-Q) plots (see Fig. S1 in the supplemental material), the observed P values do not deviate from the null, except at the extreme tail. Together with the genomic inflation factors, the Q-Q plots indicate proper control of population stratification.

FIG S1

(A to E) Q-Q plots for five vaginal microbiome traits measured in 171 Kenyan women. Inspection of Q-Q plots for each microbiome trait evaluated indicate proper control of population stratification. Download FIG S1, TIF file, 0.2 MB.
Copyright © 2020 Mehta et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Figure 1A to E display the Manhattan plots of the genome-wide P values for each vaginal microbiome trait. Table 2 summarizes the top SNPs (P < 1 × 10−5) for each vaginal microbiome trait and minor allele frequencies (MAF) derived from our data. Overall, no SNP reached genome-wide significance (0.05/336,151 = 1.49E−07). The SNP with the lowest adjusted P value associated with the presence of L. crispatus was rs73330467 (P = 4.79 × 10−6; Fig. 1A), located in an intergenic region between LOC101927488 and GRAMD2B on chromosome 5. The minor allele G (MAF = 0.07) was associated with an 11.6-fold increased odds of detecting L. crispatus. The most significant SNP associated with relative abundance of L. iners was rs527430 (P = 6.98 × 10−7; Fig. 1B) on chromosome 1, located in an intergenic region between FOXD2 and TRABD2B. The minor allele A (MAF = 0.06) was associated with an increase in relative abundance of L. iners (β = 1.05). The most significant SNP associated with G. vaginalis was rs1229660 (P = 4.65 × 10−6; Fig. 1C) located between SNX10-LOC441204 on chromosome 7. The minor allele C (MAF = 0.10) was associated with a decrease in G. vaginalis relative abundance (β = −0.99). The most significant SNP associated with the Shannon diversity index was rs972741 (P = 8.52 × 10−7; Fig. 1D) located between ZKSCAN2 and HS3ST4 on chromosome 16. The minor allele A (MAF = 0.19) was associated with an increase in the Shannon diversity index (β = 0.66). The most significant SNP associated with CST was rs2302902 (P = 3.09 × 10−6; Fig. 1E) located in ELK3 on chromosome 12. The minor allele T (MAF = 0.20) was associated with an increased likelihood of membership in lower-diversity CST-I and CST-II compared to CST-IV (β = 0.41). Sensitivity analyses were performed to further adjust for HIV and herpes simplex virus (HSV) status and yielded similar results and conclusions (see Table S1 in the supplemental material). Results for the comparison of the top SNPs between the vaginal microbiome traits are presented in Table S2. Overall, each microbiome trait was associated with at least one other top SNP, i.e., 1 for L. crispatus, 4 for L. iners, 4 for G. vaginalis, 3 for Shannon diversity index, and 5 for CST.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

(A to E) Manhattan plots for single nucleotide polymorphisms associated with vaginal microbiome traits. The x axis corresponds to the genomic position, and the y axis shows the −log10 of the P value. The horizontal red line corresponds to the genome-wide significance threshold line corresponds to the genome-wide significance threshold.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 2

Top single nucleotide polymorphisms associated with five vaginal microbiome traits

TABLE S1

Sensitivity analysis of the top GWAS SNPs adjusting for HIV and HSV-2 status. SNP, single nucleotide polymorphism; Chr, chromosome; A1/A2, allele 1/allele 2. aReported as odds ratio. Download Table S1, DOCX file, 0.01 MB.
Copyright © 2020 Mehta et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S2

Comparison of top SNPs across vaginal microbiome traits. P values corresponding to the top SNPs identified during single trait analysis are shown in italics. SNPs with P < 3.13 × 10−3 (0.05/16) are shown in boldface. Shannon, Shannon diversity index; CST, community state type. Download Table S2, DOCX file, 0.02 MB.
Copyright © 2020 Mehta et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Results from imputed SNPs.We performed genotype imputation in order to expand genomic coverage and interrogate additional SNPs not directly genotyped. After removing SNPs with low imputation quality and low frequency (R-squared [Rsq] < 0.30 and MAF < 1%), several imputed SNPs (but no additional genomic regions) were identified. Figure 2 presents the regional SNP association plots for each of the vaginal microbiome traits. Within the ZKSCAN2-HS3ST4 intergenic region on chromosome 16, several imputed SNPs were more significantly associated with the Shannon diversity index than those directly genotyped (Fig. 2E). The most significant SNP in this region was rs115869045 (P = 4.77 × 10−7; Rsq = 0.70).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

(A to E) Regional association plots of the top genomic loci associated with vaginal microbiome traits. P values (−log10) of the GWAS (solid circles) on the y axis are plotted against the genomic positions of each SNP on the x axis for each microbiome trait. Genes in the region are shown below. The linkage disequilibrium (LD) values (r2) between the lead SNP and the other SNPs are indicated in different colors.

Conditional analyses.We conducted conditional analyses for the genomic loci in which the most significant SNP was identified for each vaginal microbiome trait to determine whether additional SNPs (±400 kb) contribute to the observed associations. Figure S2 displays the regional SNP associations conditioning on the most significant directly genotyped SNP by including the SNP as a covariate in the regression model. As shown in each plot, all of the immediate SNP associations reduced toward the null, suggesting that the SNPs identified are the leading SNPs for the observed associations.

FIG S2

(A to E) Regional conditional association plots of the top genomic loci associated with vaginal microbiome traits. P values (−log10) of the GWAS (solid circles) on the y axis are plotted against the genomic positions of each SNP on the x axis for each microbiome trait. Genes in the region are shown below. The linkage disequilibrium (LD) values (r2) between the lead SNP and the other SNPs are indicated in different colors. Download FIG S2, TIF file, 0.9 MB.
Copyright © 2020 Mehta et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Analysis of previously reported loci.We investigated the associations of previously reported SNPs in relation to BV or vaginal microflora to determine whether these associations are relevant in a Kenyan population (19–26). Figure 3 and Table S3 summarize these findings. Overall, 49 autosomal SNPs were identified from previous studies, of which 45 SNPs exhibited high imputation quality (Rsq ≥ 0.30). When analyzed, 7 SNPs for L. iners, 8 SNPs for L. crispatus, 6 SNPs for G. vaginalis, 10 SNPs for the Shannon diversity index, and 8 SNPs for CST had a P < 0.05 in our sample. These findings, however, were not significant after applying a Bonferroni correction (0.05/49 = 0.00102).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Heatmap summarizing beta coefficients for the most significant SNPs within 100 kbp of previously reported SNPs. The heatmap on the left represents the direction and magnitude of the coefficients (beta) for the SNPs within 100 kb associated with five vaginal microbiome traits: (A) relative abundance of L. iners, (B) presence of L. crispatus, (C) relative abundance of G. vaginalis, (D) Shannon diversity index, (E) community state type (CST), with CST-IV as the reference. Negative coefficients are shaded in blue, and positive coefficients are shaded in red, with deeper intensity representing the magnitude of the coefficient. The heatmap on the right represents the P value of the corresponding coefficients for SNPs and vaginal microbiome traits (A to E). All associations shown have P values of <0.05. The Bonferroni cutoff for significance is <0.00102, and these P values are shaded yellow to purple. P values of >0.00102 to 0.05 are shaded gray.

TABLE S3

Analysis of previously reported loci. Note that this is an Excel file. Download Table S3, XLSX file, 0.03 MB.
Copyright © 2020 Mehta et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

To replicate genomic loci associated with the microbiome traits, we extracted the most significant imputed SNP within 100 kb of the previously reported SNPs for each microbiome trait in our study. This approach resulted in 6 genomic regions for L. iners, 1 genomic region for G. vaginalis, 7 genomic regions for Shannon diversity index, and 4 genomic regions for CST that were significant after Bonferroni correction (Fig. 3 and Table S3). No regions were significant after multiple testing correction for L. crispatus. The replicated genomic regions include IL1RN (CST), IL-5 (Shannon), IL5RA (L. iners, Shannon), IL-8 (Shannon, CST), TIRAP (L. iners, Shannon), TLR2 (Shannon, CST), TLR3 (L. iners, Shannon), TLR4 (L. iners), MBL2 (L. iners, G. vaginalis, CST), and MYD88 (L. iners, Shannon).

Gene, pathway, and disease/phenotype analyses.During gene-based analyses, no genes remained significant after multiple testing correction for any of the microbiome traits (Table S4). Table 3 presents the top five Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome pathways from WebGestalt for each vaginal microbiome trait. Four pathways were significantly associated with L. crispatus after Benjamini-Hochberg multiple testing adjustment, including G-protein-coupled receptor (GPCR) ligand binding and neutrophil degranulation. Five pathways were associated with L. iners: 4 related to TLR cascades, and one to the TLR upstream component (MyD88:Mal TIRAP [MyD88-adaptor-like TIR {Toll-interleukin 1 receptor} domain containing adaptor protein]). Two biological pathways related to GPCR signaling were associated with G. vaginalis, as well as “neutrophil degranulation,” though it did not meet Bonferroni significance. Similarly, two pathways were associated with CST, including the most significant pathway “Class B/2 (secretin family receptors),” and one pathway was associated with the Shannon diversity index, “Class B/2 (secretin family receptors).”

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 3

Results of pathway analysis: top KEGG and Reactome pathways associated with vaginal microbiome traits

TABLE S4

Top genes from gene level analysis for each vaginal microbiome trait. Bonferroni corrected P value, 2.36 × 10−6 (0.05/21,213). Download Table S4, DOCX file, 0.01 MB.
Copyright © 2020 Mehta et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Table 4 summarizes the top five DisGeNET, GLAD4U, Online Mendelian Inheritance in Man (OMIM), and Human Phenotype Ontology terms from WebGestalt for each vaginal microbiome trait. Five phenotypes were significantly associated with L. crispatus after Benjamini-Hochberg multiple testing adjustment, including “Bacterial infections” (P = 0.011) as the most significant phenotype. Three phenotypes were associated with L. iners, with “Autosomal dominant inheritance” the most significant (P = 0.018). No phenotypes were significantly associated with G. vaginalis. One phenotype was associated with the Shannon diversity index, “Abnormality of the female genitalia” (P = 0.037). Five phenotypes were associated with CST, including “Abnormality of the integument” (P = 0.00023).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 4

Results of disease/phenotype analysis: top phenotypes associated with vaginal microbiome traits

DISCUSSION

We observed associations between genetic and biological pathways in native Kenyan women with biological relevance to host immunity, cell signaling, and infection, supporting the role of host genetics in interindividual variability in vaginal microbiome traits.

In GWAS, no SNP reached genome-wide significance. The most significant SNP associated with CST (rs2302902) is located in an intron of ELK3, a member of the ETS transcription factor family, which has not been previously associated with any vaginal microbiome traits or BV. Proteins in the ELK3 subfamily are recruited by serum response factor and participate in transcription regulation, and they have been associated with various cancers (27–29), including human papillomavirus (HPV)-positive tumors of oropharyngeal cancer (30) and HPV16 in cervical tumors (31). In meta-analysis, molecular and clinical measures of BV are associated with a twofold increase in risk of high-grade cervical epithelial neoplasia (CIN) or cancer (32). Nonoptimal vaginal microbiome leads to mucosal disruption and a proinflammatory environment that can facilitate HPV acquisition (33), but there may also be underlying shared or synergistic mechanism with ELK3. Other SNPs identified in GWAS occurred primarily in intergenic regions that have not previously been associated with BV or vaginal microbiome traits, and they should be assessed for replication in future studies.

BV was absent from CST-I (L. crispatus dominant) and infrequent (5.1%) in CST-III. Although CST-IV accounted for 89% of BV cases in this sample, 57% of women in CST-IV did not have BV. This is in keeping with many other studies: women with diverse CSTs are much more likely to have BV, but as many as 50% of women with these CSTs do not have BV (34). Host genetics may partially explain this variability. For example, women with similar vaginal bacterial colonization (such as those within CST-IV) may have different response to bacterial colonization depending on genetic traits. For women who mount an inflammatory response, this immune response may potentiate BV rather being solely a response to BV. In our study, restricting analyses to HIV uninfected women did not change our findings. This is not entirely unexpected, as associations between HIV and subsequent impact on the vaginal microbiome are variable, with some finding no effect (35, 36). Conversely, there are substantial data that demonstrate the temporal association between nonoptimal vaginal microbiome composition, mucosal inflammation, and subsequent HIV risk (37), and our results suggest that host genetic differences may underpin this. Well-powered longitudinal studies that examine how host genetics is associated with different CSTs, mucosal immune markers, and BV trajectories could be informative to understanding this variability in pathology within nonoptimal CSTs. Additionally, studies are needed to explore the potential long-range regulation of these SNPs on genes, as well as evaluate for expression or methylation quantitative trait locus associations.

Complex traits are polygenic, and our analysis of previously reported loci identified multiple SNPs across numerous genetic loci with Bonferroni corrected significance. Several SNPs in the MBL2 region were significantly associated with L. iners, G. vaginalis, and CST. Kalia et al. studied these previously reported MBL2 SNPs in association with BV in women in India due to associations between mannose binding lectin (MBL) insufficiency and increased susceptibility to infectious diseases (26). Among Brazilian women, a specific MBL2 polymorphism was associated with increased odds of recurrent BV compared to controls (38), though a study among Italian women that screened for three specific MBL2 SNPs found no differences for women with recurrent BV compared to controls (39). The relevance of this gene may be population specific or subject to interaction with other factors such as the prevalence of coinfections (e.g., HSV-2, HIV), sexual exposures (e.g., partner circumcision status, multiple sex partners, condom use), or nonsexual factors (e.g., cigarette smoking, intravaginal practices). Additionally, we identified SNPs in IL-1RN (CST), IL-5 (Shannon), and IL-5RA (L. iners, Shannon). IL-1RN (the gene for interleukin 1 receptor agonist [IL-1RA]) inhibits the inflammatory IL-1 cytokines, IL-1α and IL-1β. In a study of 62 Italian women, compared to women with Nugent scores of 0 to 3, women with BV (Nugent scores of 7 to 10) had increased IL-1RA, and increased inflammatory IL-5 was associated with decreased Lactobacillus (40). In candidate gene analysis, Si et al. targeted genetic variants of IL-5, finding association with increased abundances of various Prevotella (24). Our findings further these studies by analyzing genetic factors and vaginal microbiome traits in a novel population and by demonstrating trans-ethnic associations for these previously identified loci.

In pathway analysis, GPCR pathways were identified with L. crispatus, G. vaginalis, and CST, and class B/2 (secretin family receptors) within the GPCR family were identified with Shannon diversity index and CST. GPCRs play a significant role in intracellular signaling in response to a broad range of stimuli (e.g., hormones, neurotransmitters, proteins, etc.) and across a wide range of functions (e.g., growth, nutrition requirements, response to disease, etc.) (41). This likely includes vaginal microbiome modulation, and variation in GPCR signaling may be related to variations in innate immune response. LL-37 is elevated in women with BV, and in an ex vivo endocervical model, after application of a GPCR inhibitor, LL-37-mediated induction of IL-8 production was inhibited (42); this is relevant given that IL-8 is often elevated in women with BV (43). In analyses of previously reported SNPs, we found that SNPs in the IL-8 gene region were associated with Shannon diversity index and CST. In an experimental study, Mares et al. demonstrate that application of cervicovaginal lavage (CVL) from women with BV induced higher levels of IL-8 and NF-κB in human monocytes than CVL from women without BV (44). Among women in the CAPRISA 004 trial, those with vaginal CST-IV had increased mucosal inflammation (including elevated IL-8) and increased risk of HIV seroconversion (37). The effects of variation in GPCR signaling on innate immune response may extend to other outcomes. Among Dutch women, GWAS and pathway analysis found that genes encoding GPCR signaling were enriched for 71 Chlamydia trachomatis-seropositive women compared to 169 control women (45). GPCR signaling encompasses a broad range of actions, but in conjunction with our pathway results implicating IL-8 and TLRs and other studies showing potential roles for GPCR in cervicovaginal health, GPCR signaling merits further study as a potential pathway affecting acquisition or maintenance of vaginal L. crispatus, CST membership, and subsequent BV status.

Our pathway analysis additionally linked neutrophil degranulation with L. crispatus and G. vaginalis. Neutrophils are effector cells in innate immunity, and degranulation can have serious adverse consequences for the host due to release of damaging molecules. We are unaware of other studies linking neutrophil degranulation to the vaginal microbiome, BV, or other vaginal infections. However, myeloperoxidase (MPO) is a measure of neutrophils and is elevated in women with sexually transmitted infections (STIs) (46). Whether neutrophil degranulation is an important process in vaginal microbiome composition may be evaluated in subsequent studies, such as by targeting genes encoding molecules associated with this pathway (e.g., β-arrestins, Rac2) (47).

Last, pathway analysis identified “MyD88:MAL (TIRAP) cascade” with L. iners, and we identified SNPs within the MyD88 gene in relation to L. iners and Shannon diversity index. MyD88:Mal (TIRAP) is an adaptor to MyD88, the first downstream component of TLR4 and TLR2, and part of the larger overall process for inflammatory cytokine effects (48). In pathway analysis, several TLR signaling associations were associated with L. iners relative abundance, and in analysis of previously reported SNPs, TLR2 and TLR3 had associations with Shannon diversity index, CST, and L. iners. As summarized in Table S3 in the supplemental material, several investigators have found polymorphisms in TLR genes related to altered risk of BV and enrichment of BV-related bacteria. TLRs, especially TLR2 alone or in combination with TLR1/TLR6 and TLR4, have been recognized for their association with BV and BV-related bacteria, as these TLRs precipitate release of proinflammatory cytokines and recruitment of inflammatory cells (20–23, 25). As summarized by Taylor et al., TLR gene variants can lead to inadequate or excess inflammatory immune response, which can affect disease susceptibility or progression (25). Additionally, studies are needed to evaluate the potential biological role TLRs have on vaginal microbiome composition. Given that the microbiome is a complex trait and complex traits are pleiotropic, this may provide understanding of how TLR variants contribute to risk of BV, response to BV treatment, or other outcomes of nonoptimal vaginal microbiome composition, such as adverse pregnancy outcomes.

We hypothesized that host genetic factors might be associated with vaginal microbiome traits, as has been found in the gut and airway, including associations with host genetic components encoding variants in the innate immune system (49, 50). As summarized by Thaiss et al., the innate immune system senses bacteria and their metabolic products and translates this to host physiologic responses to support organismal homeostasis (51). For example, microbiota influence the function of interleukins, potentially driving inflammation, and MyD88 expression at the epithelial barrier supports intestinal homeostasis (51). Taken together, our findings demonstrate numerous biologically relevant pathways associated with vaginal microbiome traits. However, such associations do not imply causality, and additional study is necessary to address unknown gene expression, pathways with multiple functions, and nonspecificity of gene polymorphisms or pathways with the vaginal microbiome, with larger study in multiple ethnic groups.

Limitations.The power of this study is limited by small sample size: in our GWAS, no SNPs reached genome-wide significance after multiple correction, and we were unable to assess genes by environment interactions. We calculated the number of associations passing a range of P value threshold counts (Table S5), which shows that we would not have detected significant SNPs at the genome-wide significance level (5 × 10−8) or the Bonferroni significance level of our study (1.49 × 10−7). It is highly likely we have false-negative findings due to small sample size. To understand the desired sample size for a study replicating our approach, we estimated minimum sample sizes necessary to reach the Bonferroni significance level (Table S6): GWAS study of these microbiome traits in a sample similar to ours will need a sample size of at least 188 to 376 women, and sample size varies based on trait and measure of association.

TABLE S5

Count of top single nucleotide polymorphisms (SNPs) by P value threshold for each vaginal microbiome trait. *5.00 × 10−8 is the conventional genome-wide significance level used in genome-wide association studies. 1.49 × 10−7 is the Bonferroni significance level in this study. Download Table S5, DOCX file, 0.01 MB.
Copyright © 2020 Mehta et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S6

Power table representing the minimum sample sizes needed to reach Bonferroni significance for each vaginal microbiome trait, expressed at three statistical power settings: 80%, 90%, and 95% power. *Bonferroni significance P value in our study, P = 1.49 × 10−7. **“N Extra” represents the number of samples needed in addition to the current samples in the study. SNP, single nucleotide polymorphism; Chr, chromosome. Download Table S6, DOCX file, 0.02 MB.
Copyright © 2020 Mehta et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

We found substantial functional overlap in several biologically relevant innate immunity pathways. This may in part be driven by the compositionality of individual taxa. To mitigate this, we limited examination to a few traits of known importance in conjunction with comprehensive community characterization. We analyzed CST linearly, though multinomial or quasi-Poisson may have been more applicable given the likely unequal effects across levels; however, these approaches are not supported by genetic statistical analysis software. We also confirmed numerous SNPs in close proximity to genomic loci associated with innate immunity reported in previous studies of host genetic traits associated with BV and vaginal taxa. Although we identified novel candidate regions associated with vaginal microbiome traits, additional studies are needed to replicate our findings. Despite the limitations of small sample size, the advantage of our unbiased discovery approach is that it can more easily be combined with other future studies. Evaluation of functional impact with multiomic studies (such as gene expression, epigenetics, and fine mapping) would allow mechanistic inference and aid in elucidating regulatory elements, which could lead to translational advances in optimizing the vaginal microbiome and preventing BV and associated consequences.

Conclusions.Given the multiple associations of important microbiome traits with various immune signaling pathways, the composition of the vaginal microbiome (and subsequently, BV risk) may be influenced by host genetics through immune phenotypes. Replication studies and functional studies are needed to confirm and further elucidate understanding of how vaginal bacteria interact with the genetic makeup of the host to affect nonoptimal vaginal microbiome and related outcomes. Given the cooccurrence of risk for BV, nonoptimal vaginal microbiome, and HIV and STI among African and African American women, well-powered cross-national studies are needed to answer whether common variants are replicated across populations and to identify important genes by environment interactions that inform how much risk for nonoptimal vaginal microbiome may be modifiable, such as through behavioral modification or prevention or treatment of coinfections (e.g., HSV-2, HIV).

MATERIALS AND METHODS

This study was approved by the Ethical Review Committee of Maseno University (Kisumu, Kenya; MSU/DRPC/MUERC/00054/13; 13 January 2014), and the Institutional Review Board of the University of Illinois at Chicago (USA; 2013-0511; 12 February 2014).

Study design and participants.Subjects in this analysis were enrolled in Afya Jozi, Afya Jamii (Kiswahili for “Healthy Pair, Healthy Community”), a prospective cohort study of heterosexual couples in Kisumu, Kenya. Recruitment and eligibility criteria have been previously published (52). Briefly, we included men aged 18 to 35 years old and their female partners aged 16 years and older, who independently confirmed a sexual relationship of at least 6 months duration. Men or women reporting antibiotic use within the past 60 days, men who were circumcised within the past 6 months, and women who reported vaginal douching within the past 7 days were temporarily excluded until sufficient time had passed. After the baseline visit, couples were scheduled for follow-up at 1 month, 6 months, and 12 months. Couples were enrolled from 1 April 1 2014 through 22 June 2016, and 12-month follow-up was completed 21 June 2017.

Specimen collection.At each scheduled study visit, participants underwent a standardized medical history and physical examination, followed by a personal interview to obtain sociodemographic and behavioral information. At baseline and each follow-up visit, clinicians collected vaginal swabs for assessment of BV. Gram-stained slides were evaluated according to Nugent’s criteria, and a score of 7 to 10 was defined as BV following clinical recommendations (53). Vaginal swab collection was followed by collection of CVL in 10 ml sterile water. Aliquots of 2 ml were stored at −80°C until shipment to Chicago, IL, for processing.

Vaginal microbiome characterization.We used the vaginal microbiome as measured from the baseline visit for this analysis, because antimicrobial treatment was provided for BV and other cervicovaginal infections which may confound subsequent vaginal microbiome observations. DNA extraction, library preparation, and sequencing were performed at the University of Illinois at Chicago Sequencing Core (UICSQC). DNA was extracted from CVL specimens using a PowerFecal DNA kit (Qiagen). A two-stage PCR protocol was used to amplify the V3-V4 variable region of bacterial 16S rRNA genes (54). Illumina MiSeq was used to sequence the amplicons after pooling, implementing V3 chemistry (600 cycles). Quality control and taxonomic annotation were conducted by University of Maryland Institute for Genomic Science (UMD IGS) following previously published protocols (55). Community state types (CSTs) were determined by the VALENCIA algorithm implemented by UMD IGS, which uses a distance-based metric to classify each sample to a CST based on the similarity of that sample to the centroid of CSTs identified in a reference set (56). Briefly, dominant taxa by CST are as follows: CST-I, L. crispatus; CST-II, Lactobacillus gasseri; CST-III, L. iners; CST-IV, diverse; CST-V, Lactobacillus jensenii (12).

Genotyping and quality control.Genomic DNA was obtained from 200 Kenyan women using oral swabs and was extracted using QIAamp DNA blood minikit (Qiagen). DNA extraction and isolation were performed at the Center for Population Epigenetics at Northwestern University. Study participants were genotyped using the Infinium Global Screening Array (∼640,000 markers; Illumina, Inc., San Diego, CA). The software Illumina GenomeStudio (v2.0.4; Illumina, Inc., San Diego, CA) was used to call single nucleotide polymorphisms (SNPs) at the Genomic Facility, University of Chicago. Study participants with a genotyping call rate of less than 95% were excluded from analysis as were study participants with missing vaginal microbiome traits. Due to cryptic relatedness identified during quality control, we restricted all analyses to the maximum number of unrelated individuals. Implementation of these exclusion criteria resulted in 171 study participants for analysis. SNPs with a call rate of less than 95%, a minor allele frequency of less than 1%, or a Hardy-Weinberg equilibrium P value of less than 10−6 were omitted. After these exclusion criteria, 336,151 remained for further analysis.

Genotype imputation.We performed genotype imputation to identify additional SNPs associated with the vaginal microbiome traits that were not directly genotyped using the Luhya in Webuye, Kenya, reference panel from the 1000 Genomes Project (57). The program Shapeit2 (58) was used to phase genotypes, and imputation was performed using Minimac3 (59). Imputed genotypes were coded as allelic dosages (estimated counts ranging from 0 to 2). Low-quality imputed SNPs, i.e., Rsq < 0.30 and MAF < 1%, were excluded from further analyses.

Statistical analysis.We evaluated the association between single nucleotide polymorphisms (SNPs) with five vaginal microbiome traits: (i) presence of L. crispatus, (ii) relative abundance of L. iners, (iii) relative abundance of G. vaginalis, (iv) Shannon diversity index, and (v) CST. We chose L. crispatus for its known importance in vaginal microbiome health, L. iners for its predominance among Lactobacillus species in our sample and other African samples, and G. vaginalis for its inclusion as a morphotype in Nugent scoring and determination of BV. For alpha diversity measures, we selected the Shannon diversity index, as all alpha diversity measures were highly correlated (Spearman correlations, 0.77 to 0.98). Prior to analyses, the relative abundance of L. iners and Shannon diversity index were inverse log transformed to approximate a more normal distribution (see Fig. S1 in the supplemental material). As transformation did not lead to normal approximation, the relative abundance of G. vaginalis was analyzed by quartile. Given the infrequency of L. crispatus (present in 24% of samples), we analyzed it as presence versus absence. Of 171 subjects in the analysis, 168 (98%) baseline vaginal microbial communities were classified as CST-I (n = 14; 0% with BV), CST-III (n = 77; 5.2% with BV), or CST-IV (n = 77; 43% with BV). There were three outliers (two with CST-II and one with CST-V), and these observations were excluded from CST analyses. BV was not analyzed as a separate outcome because it was colinear with CST. We analyzed CST as a linear outcome (i.e., equal distance between CST-I, CST-II, and CST-III) given the ordinal nature of BV risk and to meet parameters for statistical software.

Principal components of genetic ancestry were generated using the program EIGENSOFT (60), and the first three principal components were retained and included as covariates. Quantile-quantile (Q-Q) plots were generated for each trait to visualize test statistics, and the genomic control inflation factors were calculated. We used PLINK (v1.90) (61) to test the associations between L. iners, G. vaginalis, Shannon diversity index, and CST with SNPs using linear regression, and association between L. crispatus and SNPs using logistic regression, adjusting for age and principal components. We assumed an additive genetic effect model for all analyses. To account for genotype uncertainty, the allelic dosages from imputed SNPs were analyzed using Mach2QTL for linear regression and Mach2DAT for logistic regression (see “Web resources” below). SNPs with a P value of less than 1.49 × 10−7 (0.05/336,151) were declared genome-wide significant. The software R (3.5.0) (62) and LocusZoom (63) (hg19/1000 Genomes AFR) were used for graphing. We calculated the count of SNPs meeting a range of P value thresholds (see Table S5 in the supplemental material). Additionally, we estimated the minimum sample sizes needed to reach Bonferroni level significance (1.49 × 10−7) based on linear regression of each trait and SNPs at 80%, 90%, and 95% power thresholds (Table S6). We were unable to generate a reliable sample size estimate for L. crispatus due to the odds ratio estimate and low prevalence of the trait.

Gene, pathway, and disease/phenotype analyses.We performed gene, pathway, and disease/phenotype analyses using SKAT-O (64) and WebGestalt (65) to examine the combined genetic effects on the underlying biological mechanisms influencing the vaginal microbiome traits. We mapped all directly genotyped SNPs to autosomal genes based on the genomic positions of the GRCh37/hg19 assembly and included a ±20-kb gene boundary to capture proximal regulatory and functional elements that contribute to gene regulation, resulting in 21,213 genes for analysis. Gene set associations were performed using SKAT-O, adjusting for age and the first three principal components. Overrepresentation analyses of pathways (Kyoto Encyclopedia of Genes and Genomes [KEGG] and Reactome) and diseases/phenotypes (DisGeNET, GLAD4U, Online Mendelian Inheritance in Man [OMIM], and Human Phenotype Ontology) were analyzed using WebGestalt. Genes with a P value of less than 2.36 × 10−6 (0.05/21,213) were declared significant. Pathways, diseases, and phenotypes with a false-discovery rate (FDR) P value of <0.05 after Benjamini-Hochberg multiple testing adjustment were declared significant. The top five genes, pathways, and diseases/phenotypes from SKAT-O and WebGestalt are reported.

Web resources.The Web resources are as follows: Mach2DAT (http://csg.sph.umich.edu/abecasis/MACH/download/) and Mach2QTL (http://csg.sph.umich.edu/abecasis/MACH/download/).

Data availability.Raw sequence data files are available in the Sequence Read Archive (National Center for Biotechnology Information; BioProject identifier PRJNA516684).

DATA SET S1

Final analytic data set GWAS VMB (n = 171). Note that this is an Excel file. Download Data Set S1, XLS file, 0.05 MB.
Copyright © 2020 Mehta et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

ACKNOWLEDGMENTS

This study was supported by grant R01-AI110369 (principal investigator [PI], S. D. Mehta) from the National Institutes of Health, National Institute of Allergy and Infectious Diseases, Division of Microbiology, and by the National Cancer Institute of the National Institutes of Health under award U54CA221205 (multiple principal investigators [MPI], L. Hou, R. Murphy).

The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

No conflicts of interest have been declared.

FOOTNOTES

    • Received June 5, 2020.
    • Accepted July 14, 2020.
  • Copyright © 2020 Mehta et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Torrone EA,
    2. Morrison CS,
    3. Chen PL,
    4. Kwok C,
    5. Francis SC,
    6. Hayes RJ,
    7. Looker KJ,
    8. McCormack S,
    9. McGrath N,
    10. van de Wijgert J,
    11. Watson-Jones D,
    12. Low N,
    13. Gottlieb SL
    , STIMA Working Group. 2018. Prevalence of sexually transmitted infections and bacterial vaginosis among women in sub-Saharan Africa: an individual participant data meta-analysis of 18 HIV prevention studies. PLoS Med 15:e1002511. doi:10.1371/journal.pmed.1002511.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Atashili J,
    2. Poole C,
    3. Ndumbe PM,
    4. Adimora AA,
    5. Smith JS
    . 2008. Bacterial vaginosis and HIV acquisition: a meta-analysis of published studies. AIDS 22:1493–1501. doi:10.1097/QAD.0b013e3283021a37.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    1. Lewis FM,
    2. Bernstein KT,
    3. Aral SO
    . 2017. Vaginal microbiome and its relationship to behavior, sexual health, and sexually transmitted diseases. Obstet Gynecol 129:643–654. doi:10.1097/AOG.0000000000001932.
    OpenUrlCrossRef
  4. 4.↵
    1. van de Wijgert J,
    2. Jespers V
    . 2017. The global health impact of vaginal dysbiosis. Res Microbiol 168:859–864. doi:10.1016/j.resmic.2017.02.003.
    OpenUrlCrossRef
  5. 5.↵
    1. Fredricks DN
    . 2011. Molecular methods to describe the spectrum and dynamics of the vaginal microbiota. Anaerobe 17:191–195. doi:10.1016/j.anaerobe.2011.01.001.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. van de Wijgert JH,
    2. Borgdorff H,
    3. Verhelst R,
    4. Crucitti T,
    5. Francis S,
    6. Verstraelen H,
    7. Jespers V
    . 2014. The vaginal microbiota: what have we learned after a decade of molecular characterization? PLoS One 9:e105998. doi:10.1371/journal.pone.0105998.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Fethers KA,
    2. Fairley CK,
    3. Hocking JS,
    4. Gurrin LC,
    5. Bradshaw CS
    . 2008. Sexual risk factors and bacterial vaginosis: a systematic review and meta-analysis. Clin Infect Dis 47:1426–1435. doi:10.1086/592974.
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    1. Gray RH,
    2. Kigozi G,
    3. Serwadda D,
    4. Makumbi F,
    5. Nalugoda F,
    6. Watya S,
    7. Moulton L,
    8. Chen MZ,
    9. Sewankambo NK,
    10. Kiwanuka N,
    11. Sempijja V,
    12. Lutalo T,
    13. Kagayii J,
    14. Wabwire-Mangen F,
    15. Ridzon R,
    16. Bacon M,
    17. Wawer MJ
    . 2009. The effects of male circumcision on female partners’ genital tract symptoms and vaginal infections in a randomized trial in Rakai, Uganda. Am J Obstet Gynecol 200:42.e1. doi:10.1016/j.ajog.2008.07.069.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Verstraelen H,
    2. Verhelst R,
    3. Vaneechoutte M,
    4. Temmerman M
    . 2010. The epidemiology of bacterial vaginosis in relation to sexual behaviour. BMC Infect Dis 10:81. doi:10.1186/1471-2334-10-81.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Low N,
    2. Chersich MF,
    3. Schmidlin K,
    4. Egger M,
    5. Francis SC,
    6. van de Wijgert JH,
    7. Hayes RJ,
    8. Baeten JM,
    9. Brown J,
    10. Delany-Moretlwe S,
    11. Kaul R,
    12. McGrath N,
    13. Morrison C,
    14. Myer L,
    15. Temmerman M,
    16. van der Straten A,
    17. Watson-Jones D,
    18. Zwahlen M,
    19. Hilber AM
    . 2011. Intravaginal practices, bacterial vaginosis, and HIV infection in women: individual participant data meta-analysis. PLoS Med 8:e1000416. doi:10.1371/journal.pmed.1000416.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Nelson TM,
    2. Borgogna JC,
    3. Michalek RD,
    4. Roberts DW,
    5. Rath JM,
    6. Glover ED,
    7. Ravel J,
    8. Shardell MD,
    9. Yeoman CJ,
    10. Brotman RM
    . 2018. Cigarette smoking is associated with an altered vaginal tract metabolomic profile. Sci Rep 8:852. doi:10.1038/s41598-017-14943-3.
    OpenUrlCrossRef
  12. 12.↵
    1. Ravel J,
    2. Gajer P,
    3. Abdo Z,
    4. Schneider GM,
    5. Koenig SS,
    6. McCulle SL,
    7. Karlebach S,
    8. Gorle R,
    9. Russell J,
    10. Tacket CO,
    11. Brotman RM,
    12. Davis CC,
    13. Ault K,
    14. Peralta L,
    15. Forney LJ
    . 2011. Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci U S A 108:4680–4687. doi:10.1073/pnas.1002611107.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Peebles K,
    2. Velloza J,
    3. Balkus JE,
    4. McClelland RS,
    5. Barnabas RV
    . 2019. High global burden and costs of bacterial vaginosis: a systematic review and meta-analysis. Sex Transm Dis 46:304–311. doi:10.1097/OLQ.0000000000000972.
    OpenUrlCrossRef
  14. 14.↵
    1. Alcendor DJ
    . 2016. Evaluation of health disparity in bacterial vaginosis and the implications for HIV-1 acquisition in African American women. Am J Reprod Immunol 76:99–107. doi:10.1111/aji.12497.
    OpenUrlCrossRef
  15. 15.↵
    1. Fettweis JM,
    2. Brooks JP,
    3. Serrano MG,
    4. Sheth NU,
    5. Girerd PH,
    6. Edwards DJ,
    7. Strauss JF, The Vaginal Microbiome Consortium,
    8. Jefferson KK,
    9. Buck GA
    . 2014. Differences in vaginal microbiome in African American women versus women of European ancestry. Microbiology 160:2272–2282. doi:10.1099/mic.0.081034-0.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Martin DH,
    2. Marrazzo JM
    . 2016. The vaginal microbiome: current understanding and future directions. J Infect Dis 214(Suppl 1):S36–S41. doi:10.1093/infdis/jiw184.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Turovskiy Y,
    2. Sutyak Noll K,
    3. Chikindas ML
    . 2011. The aetiology of bacterial vaginosis. J Appl Microbiol 110:1105–1128. doi:10.1111/j.1365-2672.2011.04977.x.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    National AIDS Control Council. 2016. Kenya HIV county profiles. National AIDS Control Council, Nairobi, Kenya. http://nacc.or.ke/wp-content/uploads/2016/12/Kenya-HIV-County-Profiles-2016.pdf.
  19. 19.↵
    1. Goepfert AR,
    2. Varner M,
    3. Ward K,
    4. Macpherson C,
    5. Klebanoff M,
    6. Goldenberg RL,
    7. Mercer B,
    8. Meis P,
    9. Iams J,
    10. Moawad A,
    11. Carey JC,
    12. Leveno K,
    13. Wapner R,
    14. Caritis SN,
    15. Miodovnik M,
    16. Sorokin Y,
    17. O’Sullivan MJ,
    18. Van Dorsten JP,
    19. Langer O
    , NICHD Maternal-Fetal Medicine Units Network. 2005. Differences in inflammatory cytokine and Toll-like receptor genes and bacterial vaginosis in pregnancy. Am J Obstet Gynecol 193:1478–1485. doi:10.1016/j.ajog.2005.03.053.
    OpenUrlCrossRefPubMedWeb of Science
  20. 20.↵
    1. Ryckman KK,
    2. Williams SM,
    3. Krohn MA,
    4. Simhan HN
    . 2009. Genetic association of Toll-like receptor 4 with cervical cytokine concentrations during pregnancy. Genes Immun 10:636–640. doi:10.1038/gene.2009.47.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    1. Royse KE,
    2. Kempf MC,
    3. McGwin G, Jr,
    4. Wilson CM,
    5. Tang J,
    6. Shrestha S
    . 2012. Toll-like receptor gene variants associated with bacterial vaginosis among HIV-1 infected adolescents. J Reprod Immunol 96:84–89. doi:10.1016/j.jri.2012.08.002.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Taylor BD,
    2. Darville T,
    3. Ferrell RE,
    4. Ness RB,
    5. Kelsey SF,
    6. Haggerty CL
    . 2014. Cross-sectional analysis of Toll-like receptor variants and bacterial vaginosis in African-American women with pelvic inflammatory disease. Sex Transm Infect 90:563–566. doi:10.1136/sextrans-2014-051524.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Mackelprang RD,
    2. Scoville CW,
    3. Cohen CR,
    4. Ondondo RO,
    5. Bigham AW,
    6. Celum C,
    7. Campbell MS,
    8. Essex M,
    9. Wald A,
    10. Kiarie J,
    11. Ronald A,
    12. Gray G,
    13. Lingappa JR
    , Partners in Prevention HSV/HIV Transmission Study Team. 2015. Toll-like receptor gene variants and bacterial vaginosis among HIV-1 infected and uninfected African women. Genes Immun 16:362–365. doi:10.1038/gene.2015.13.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Si J,
    2. You HJ,
    3. Yu J,
    4. Sung J,
    5. Ko G
    . 2017. Prevotella as a hub for vaginal microbiota under the influence of host genetics and their association with obesity. Cell Host Microbe 21:97–105. doi:10.1016/j.chom.2016.11.010.
    OpenUrlCrossRef
  25. 25.↵
    1. Taylor BD,
    2. Totten PA,
    3. Astete SG,
    4. Ferris MJ,
    5. Martin DH,
    6. Ness RB,
    7. Haggerty CL
    . 2018. Toll-like receptor variants and cervical Atopobium vaginae infection in women with pelvic inflammatory disease. Am J Reprod Immunol 79:e12804. doi:10.1111/aji.12804.
    OpenUrlCrossRef
  26. 26.↵
    1. Kalia N,
    2. Singh J,
    3. Sharma S,
    4. Kaur M
    . 2019. Impact of SNPs interplay across the locus of MBL2, between MBL and Dectin-1 gene, on women’s risk of developing recurrent vulvovaginal infections. Cell Biosci 9:35. doi:10.1186/s13578-019-0300-4.
    OpenUrlCrossRef
  27. 27.↵
    1. Yang H,
    2. Schramek D,
    3. Adam RC,
    4. Keyes BE,
    5. Wang P,
    6. Zheng D,
    7. Fuchs E
    . 2015. ETS family transcriptional regulators drive chromatin dynamics and malignancy in squamous cell carcinomas. Elife 4:e10870. doi:10.7554/eLife.10870.
    OpenUrlCrossRef
  28. 28.↵
    1. Heo SH,
    2. Lee JY,
    3. Yang KM,
    4. Park KS
    . 2015. ELK3 expression correlates with cell migration, invasion, and membrane type 1-matrix metalloproteinase expression in MDA-MB-231 breast cancer cells. Gene Expr 16:197–203. doi:10.3727/105221615X14399878166276.
    OpenUrlCrossRef
  29. 29.↵
    1. Lee JH,
    2. Hur W,
    3. Hong SW,
    4. Kim JH,
    5. Kim SM,
    6. Lee EB,
    7. Yoon SK
    . 2017. ELK3 promotes the migration and invasion of liver cancer stem cells by targeting HIF-1alpha. Oncol Rep 37:813–822. doi:10.3892/or.2016.5293.
    OpenUrlCrossRef
  30. 30.↵
    1. Suarez E,
    2. Gonzalez L,
    3. Perez-Mitchell C,
    4. Ortiz AP,
    5. Ramirez-Sola M,
    6. Acosta J,
    7. Bernabe-Dones RD,
    8. Gonzalez-Aquino C,
    9. Montes-Rodriguez I,
    10. Cadilla CL
    . 2016. Pathway analysis using gene-expression profiles of HPV-positive and HPV-negative oropharyngeal cancer patients in a Hispanic population: methodological procedures. P R Health Sci J 35:3–8.
    OpenUrl
  31. 31.↵
    1. Thorland EC,
    2. Myers SL,
    3. Gostout BS,
    4. Smith DI
    . 2003. Common fragile sites are preferential targets for HPV16 integrations in cervical tumors. Oncogene 22:1225–1237. doi:10.1038/sj.onc.1206170.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    1. Brusselaers N,
    2. Shrestha S,
    3. van de Wijgert J,
    4. Verstraelen H
    . 2019. Vaginal dysbiosis and the risk of human papillomavirus and cervical cancer: systematic review and meta-analysis. Am J Obstet Gynecol 221:9–18.e8. doi:10.1016/j.ajog.2018.12.011.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Mitra A,
    2. MacIntyre DA,
    3. Marchesi JR,
    4. Lee YS,
    5. Bennett PR,
    6. Kyrgiou M
    . 2016. The vaginal microbiota, human papillomavirus infection and cervical intraepithelial neoplasia: what do we know and where are we going next? Microbiome 4:58. doi:10.1186/s40168-016-0203-0.
    OpenUrlCrossRef
  34. 34.↵
    1. Ravel J,
    2. Brotman RM
    . 2016. Translating the vaginal microbiome: gaps and challenges. Genome Med 8:35. doi:10.1186/s13073-016-0291-2.
    OpenUrlCrossRef
  35. 35.↵
    1. Reimers LL,
    2. Mehta SD,
    3. Massad LS,
    4. Burk RD,
    5. Xie X,
    6. Ravel J,
    7. Cohen MH,
    8. Palefsky JM,
    9. Weber KM,
    10. Xue X,
    11. Anastos K,
    12. Minkoff H,
    13. Atrio J,
    14. D’Souza G,
    15. Ye Q,
    16. Colie C,
    17. Zolnik CP,
    18. Spear GT,
    19. Strickler HD
    . 2016. The cervicovaginal microbiota and its associations with human papillomavirus detection in HIV-infected and HIV-uninfected women. J Infect Dis 214:1361–1369. doi:10.1093/infdis/jiw374.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Watts DH,
    2. Springer G,
    3. Minkoff H,
    4. Hillier SL,
    5. Jacobson L,
    6. Moxley M,
    7. Justman J,
    8. Cejtin H,
    9. O’Connell C,
    10. Greenblatt RM
    . 2006. The occurrence of vaginal infections among HIV-infected and high-risk HIV-uninfected women: longitudinal findings of the women's interagency HIV study. J Acquir Immune Defic Syndr 43:161–168. doi:10.1097/01.qai.0000242448.90026.13.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.↵
    1. Eastment MC,
    2. McClelland RS
    . 2018. Vaginal microbiota and susceptibility to HIV. AIDS 32:687–698. doi:10.1097/QAD.0000000000001768.
    OpenUrlCrossRef
  38. 38.↵
    1. Giraldo PC,
    2. Babula O,
    3. Goncalves AK,
    4. Linhares IM,
    5. Amaral RL,
    6. Ledger WJ,
    7. Witkin SS
    . 2007. Mannose-binding lectin gene polymorphism, vulvovaginal candidiasis, and bacterial vaginosis. Obstet Gynecol 109:1123–1128. doi:10.1097/01.AOG.0000260386.17555.a5.
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    1. Milanese M,
    2. Segat L,
    3. De Seta F,
    4. Pirulli D,
    5. Fabris A,
    6. Morgutti M,
    7. Crovella S
    . 2008. MBL2 genetic screening in patients with recurrent vaginal infections. Am J Reprod Immunol 59:146–151. doi:10.1111/j.1600-0897.2007.00549.x.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Campisciano G,
    2. Zanotta N,
    3. Licastro D,
    4. De Seta F,
    5. Comar M
    . 2018. In vivo microbiome and associated immune markers: new insights into the pathogenesis of vaginal dysbiosis. Sci Rep 8:2307. doi:10.1038/s41598-018-20649-x.
    OpenUrlCrossRef
  41. 41.↵
    1. Hilger D,
    2. Masureel M,
    3. Kobilka BK
    . 2018. Structure and dynamics of GPCR signaling complexes. Nat Struct Mol Biol 25:4–12. doi:10.1038/s41594-017-0011-7.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Frew L,
    2. Makieva S,
    3. McKinlay AT,
    4. McHugh BJ,
    5. Doust A,
    6. Norman JE,
    7. Davidson DJ,
    8. Stock SJ
    . 2014. Human cathelicidin production by the cervix. PLoS One 9:e103434. doi:10.1371/journal.pone.0103434.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. St John E,
    2. Mares D,
    3. Spear GT
    . 2007. Bacterial vaginosis and host immunity. Curr HIV/AIDS Rep 4:22–28. doi:10.1007/s11904-007-0004-y.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Mares D,
    2. Simoes JA,
    3. Novak RM,
    4. Spear GT
    . 2008. TLR2-mediated cell stimulation in bacterial vaginosis. J Reprod Immunol 77:91–99. doi:10.1016/j.jri.2007.04.004.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Roberts CH,
    2. Ouburg S,
    3. Preston MD,
    4. de Vries HJC,
    5. Holland MJ,
    6. Morré SA
    . 2018. Pathway-wide genetic risks in chlamydial infections overlap between tissue tropisms: a genome-wide association scan. Mediators Inflamm 2018:3434101. doi:10.1155/2018/3434101.
    OpenUrlCrossRef
  46. 46.↵
    1. Spear GT,
    2. Kendrick SR,
    3. Chen HY,
    4. Thomas TT,
    5. Bahk M,
    6. Balderas R,
    7. Ghosh S,
    8. Weinberg A,
    9. Landay AL
    . 2011. Multiplex immunoassay of lower genital tract mucosal fluid from women attending an urban STD clinic shows broadly increased IL1ss and lactoferrin. PLoS One 6:e19560. doi:10.1371/journal.pone.0019560.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Lacy P
    . 2006. Mechanisms of degranulation in neutrophils. Allergy Asthma Clin Immunol 2:98–108. doi:10.1186/1710-1492-2-3-98.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Gangloff M,
    2. Gay NJ
    . 2004. MD-2: the Toll 'gatekeeper' in endotoxin signalling. Trends Biochem Sci 29:294–300. doi:10.1016/j.tibs.2004.04.008.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    1. Goodrich JK,
    2. Davenport ER,
    3. Clark AG,
    4. Ley RE
    . 2017. The relationship between the human genome and microbiome comes into view. Annu Rev Genet 51:413–433. doi:10.1146/annurev-genet-110711-155532.
    OpenUrlCrossRef
  50. 50.↵
    1. Igartua C,
    2. Davenport ER,
    3. Gilad Y,
    4. Nicolae DL,
    5. Pinto J,
    6. Ober C
    . 2017. Host genetic variation in mucosal immunity pathways influences the upper airway microbiome. Microbiome 5:16. doi:10.1186/s40168-016-0227-5.
    OpenUrlCrossRef
  51. 51.↵
    1. Thaiss CA,
    2. Zmora N,
    3. Levy M,
    4. Elinav E
    . 2016. The microbiome and innate immunity. Nature 535:65–74. doi:10.1038/nature18847.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Mehta SD,
    2. Nordgren RK,
    3. Agingu W,
    4. Otieno F,
    5. Odongo W,
    6. Odhiambo F,
    7. Bailey RC
    . 2018. Sexual quality of life and association with HIV and sexually transmitted infections among a cohort of heterosexual couples in Kenya. J Sex Med 15:1446–1455. doi:10.1016/j.jsxm.2018.08.007.
    OpenUrlCrossRef
  53. 53.↵
    1. Nugent RP,
    2. Krohn MA,
    3. Hillier SL
    . 1991. Reliability of diagnosing bacterial vaginosis is improved by a standardized method of gram stain interpretation. J Clin Microbiol 29:297–301. doi:10.1128/JCM.29.2.297-301.1991.
    OpenUrlAbstract/FREE Full Text
  54. 54.↵
    1. Naqib A,
    2. Poggi S,
    3. Wang W,
    4. Hyde M,
    5. Kunstman K,
    6. Green SJ
    . 2018. Making and sequencing heavily multiplexed, high-throughput 16S ribosomal RNA gene amplicon libraries using a flexible, two-stage PCR protocol. Methods Mol Biol 1783:149–169. doi:10.1007/978-1-4939-7834-2_7.
    OpenUrlCrossRef
  55. 55.↵
    1. Holm JB,
    2. Humphrys MS,
    3. Robinson CK,
    4. Settles ML,
    5. Ott S,
    6. Fu L,
    7. Yang H,
    8. Gajer P,
    9. He X,
    10. McComb E,
    11. Gravitt PE,
    12. Ghanem KG,
    13. Brotman RM,
    14. Ravel J
    . 2019. Ultrahigh-throughput multiplexing and sequencing of >500-base-pair amplicon regions on the Illumina HiSeq 2500 platform. mSystems 4:e00029-19. doi:10.1128/mSystems.00029-19.
    OpenUrlAbstract/FREE Full Text
  56. 56.↵
    1. France M,
    2. Ma B,
    3. Gajer P,
    4. Brown S,
    5. Humphrys M,
    6. Holm J,
    7. Brotman R,
    8. Ravel J
    . 2020. VALENCIA: a nearest centroid classification method for vaginal microbial communities based on composition. Research Square doi:10.21203/rs.2.24139/v1.
  57. 57.↵
    1000 Genomes Project Consortium, Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, Marchini JL, McCarthy S, McVean GA, Abecasis GR. 2015. A global reference for human genetic variation. Nature 526:68–74. doi:10.1038/nature15393.
    OpenUrlCrossRefPubMed
  58. 58.↵
    1. Delaneau O,
    2. Marchini J
    , 1000 Genomes Project Consortium. 2014. Integrating sequence and array data to create an improved 1000 Genomes Project haplotype reference panel. Nat Commun 5:3934. doi:10.1038/ncomms4934.
    OpenUrlCrossRefPubMed
  59. 59.↵
    1. Das S,
    2. Forer L,
    3. Schonherr S,
    4. Sidore C,
    5. Locke AE,
    6. Kwong A,
    7. Vrieze SI,
    8. Chew EY,
    9. Levy S,
    10. McGue M,
    11. Schlessinger D,
    12. Stambolian D,
    13. Loh PR,
    14. Iacono WG,
    15. Swaroop A,
    16. Scott LJ,
    17. Cucca F,
    18. Kronenberg F,
    19. Boehnke M,
    20. Abecasis GR,
    21. Fuchsberger C
    . 2016. Next-generation genotype imputation service and methods. Nat Genet 48:1284–1287. doi:10.1038/ng.3656.
    OpenUrlCrossRefPubMed
  60. 60.↵
    1. Patterson N,
    2. Price AL,
    3. Reich D
    . 2006. Population structure and eigenanalysis. PLoS Genet 2:e190. doi:10.1371/journal.pgen.0020190.
    OpenUrlCrossRefPubMed
  61. 61.↵
    1. Chang CC,
    2. Chow CC,
    3. Tellier LC,
    4. Vattikuti S,
    5. Purcell SM,
    6. Lee JJ
    . 2015. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4:7. doi:10.1186/s13742-015-0047-8.
    OpenUrlCrossRefPubMed
  62. 62.↵
    R Development Core Team. 2009. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
  63. 63.↵
    1. Pruim RJ,
    2. Welch RP,
    3. Sanna S,
    4. Teslovich TM,
    5. Chines PS,
    6. Gliedt TP,
    7. Boehnke M,
    8. Abecasis GR,
    9. Willer CJ
    . 2010. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26:2336–2337. doi:10.1093/bioinformatics/btq419.
    OpenUrlCrossRefPubMedWeb of Science
  64. 64.↵
    1. Lee S,
    2. Emond MJ,
    3. Bamshad MJ,
    4. Barnes KC,
    5. Rieder MJ,
    6. Nickerson DA, NHLBI GO Exome Sequencing Project–ESP Lung Project Team,
    7. Christiani DC,
    8. Wurfel MM,
    9. Lin X
    . 2012. Optimal unified approach for rare-variant association testing with application to small-sample case-control whole-exome sequencing studies. Am J Hum Genet 91:224–237. doi:10.1016/j.ajhg.2012.06.007.
    OpenUrlCrossRefPubMed
  65. 65.↵
    1. Liao Y,
    2. Wang J,
    3. Jaehnig EJ,
    4. Shi Z,
    5. Zhang B
    . 2019. WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res 47:W199–W205. doi:10.1093/nar/gkz401.
    OpenUrlCrossRef
PreviousNext
Back to top
Download PDF
Citation Tools
Host Genetic Factors Associated with Vaginal Microbiome Composition in Kenyan Women
Supriya D. Mehta, Drew R. Nannini, Fredrick Otieno, Stefan J. Green, Walter Agingu, Alan Landay, Yinan Zheng, Lifang Hou
mSystems Jul 2020, 5 (4) e00502-20; DOI: 10.1128/mSystems.00502-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSystems article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Host Genetic Factors Associated with Vaginal Microbiome Composition in Kenyan Women
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Host Genetic Factors Associated with Vaginal Microbiome Composition in Kenyan Women
Supriya D. Mehta, Drew R. Nannini, Fredrick Otieno, Stefan J. Green, Walter Agingu, Alan Landay, Yinan Zheng, Lifang Hou
mSystems Jul 2020, 5 (4) e00502-20; DOI: 10.1128/mSystems.00502-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

vaginal microbiome
vaginal microbiota
bacterial vaginosis
Gardnerella vaginalis
Lactobacillus crispatus
Lactobacillus iners
L. crispatus
L. iners
G. vaginalis
Shannon diversity index
community state type
genome wide association study
pathway analysis
Toll-like receptors
innate immune response
Kenya
community state type
GWAS

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSystemsJ

@ASMicrobiology

       

 

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5077