Chromosome 10q24.32 Variants Associate With Brain Arterial Diameters in Diverse Populations: A Genome-Wide Association Study.

BACKGROUND: Brain arterial diameters (BADs) are novel imaging biomarkers of cerebrovascular disease, cognitive decline, and dementia. Traditional vascular risk factors have been associated with BADs, but whether there may be genetic determinants of BADs is unknown. METHODS AND RESULTS: The authors studied 4150 participants from 6 geographically diverse population-based cohorts (40% European, 14% African, 22% Hispanic, 24% Asian ancestries). Brain arterial diameters for 13 segments were measured and averaged to obtain a global measure of BADs as well as the posterior and anterior circulations. A genome-wide association study revealed 14 variants at one locus associated with global BAD at genome-wide significance (P<5×10-8) (top single-nucleotide polymorphism, rs7921574; ß=0.06 [P=1.54×10-8]). This locus mapped to an intron of CNNM2. A trans-ancestry genome-wide association study meta-analysis identified 2 more loci at NT5C2 (rs10748839; P=2.54×10-8) and AS3MT (rs10786721; P=4.97×10-8), associated with global BAD. In addition, 2 single-nucleotide polymorphisms colocalized with expression of CNNM2 (rs7897654; ß=0.12 [P=6.17×10-7]) and AL356608.1 (rs10786719; ß=-0.17 [P=6.60×10-6]) in brain tissue. For the posterior BAD, 2 variants at one locus mapped to an intron of TCF25 were identified (top single-nucleotide polymorphism, rs35994878; ß=0.11 [P=2.94×10-8]). For the anterior BAD, one locus at ADAP1 was identified in trans-ancestry genome-wide association analysis (rs34217249; P=3.11×10-8). CONCLUSIONS: The current study reveals 3 novel risk loci (CNNM2, NT5C2, and AS3MT) associated with BADs. These findings may help elucidate the mechanism by which BADs may influence cerebrovascular health.

Using source-associated mobile genetic elements to identify zoonotic extraintestinal E. coli infections.

A one-health perspective may provide new and actionable information about Escherichia coli transmission. E. coli colonizes a broad range of vertebrates, including humans and food-production animals, and is a leading cause of bladder, kidney, and bloodstream infections in humans. Substantial evidence supports foodborne transmission of pathogenic E. coli strains from food animals to humans. However, the relative contribution of foodborne zoonotic E. coli (FZEC) to the human extraintestinal disease burden and the distinguishing characteristics of such strains remain undefined. Using a comparative genomic analysis of a large collection of contemporaneous, geographically-matched clinical and meat-source E. coli isolates (n = 3111), we identified 17 source-associated mobile genetic elements - predominantly plasmids and bacteriophages - and integrated them into a novel Bayesian latent class model to predict the origins of clinical E. coli isolates. We estimated that approximately 8 % of human extraintestinal E. coli infections (mostly urinary tract infections) in our study population were caused by FZEC. FZEC strains were equally likely to cause symptomatic disease as non-FZEC strains. Two FZEC lineages, ST131-H22 and ST58, appeared to have particularly high virulence potential. Our findings imply that FZEC strains collectively cause more urinary tract infections than does any single non-E. coli uropathogenic species (e.g., Klebsiella pneumoniae). Our novel approach can be applied in other settings to identify the highest-risk FZEC strains, determine their sources, and inform new one-health strategies to decrease the heavy public health burden imposed by extraintestinal E. coli infections.

Chromosome 10q24.32 Variants Associate with Brain Arterial Diameters in Diverse Populations: A Genome-Wide Association Study.

Background: Brain arterial diameters are novel imaging biomarkers of cerebrovascular disease, cognitive decline and dementia. Traditional vascular risk factors have been associated with brain arterial diameters but whether there may be genetic determinants of brain arterial diameters is unknown. Results: We studied 4150 participants from six geographically diverse population-based cohorts (40% European, 14% African, 22% Hispanic, 24% Asian ancestries). We measured brain arterial diameters for 13 segments and averaged them to obtain a global measure of brain arterial diameters as well as the posterior and anterior circulations. A genome-wide association study (GWAS) revealed 14 variants at one locus associated with global brain arterial diameter at genome-wide significance (P<5×10-8) (top SNP, rs7921574; ß =0.06, P=1.54×10-8). This locus mapped to an intron of CNNM2. A trans-ancestry GWAS meta-analysis identified two more loci at NT5C2 (rs10748839; P=2.54×10-8) and at AS3MT (rs10786721; P=4.97×10-8), associated with global brain arterial diameter. In addition, two SNPs co-localized with expression of CNNM2 (rs7897654, ß=0.12, P=6.17×10-7) and AL356608.1 (rs10786719, ß =-0.17, P=6.60×10-6) in brain tissue. For the posterior brain arterial diameter, two variants at one locus mapped to an intron of TCF25 were identified (top SNP, rs35994878; ß =0.11, P=2.94×10-8). For the anterior brain arterial diameter, one locus at ADAP1 was identified in trans-ancestry genome-wide association analysis (rs34217249; P=3.11×10-8). Conclusion: Our study reveals three novel risk loci (CNNM2, NT5C2 and AS3MT) associated with brain arterial diameters. Our finding may elucidate the mechanisms by which brain arterial diameters influence the risk of stroke and dementia.

Admixture mapping identifies novel Alzheimer's disease risk regions in African Americans.

BACKGROUND: This study used admixture mapping to prioritize the genetic regions associated with Alzheimer's disease (AD) in African American (AA) individuals, followed by ancestry-aware regression analysis to fine-map the prioritized regions. METHODS: We analyzed 10,271 individuals from 17 different AA datasets. We performed admixture mapping and meta-analyzed the results. We then used regression analysis, adjusting for local ancestry main effects and interactions with genotype, to refine the regions identified from admixture mapping. Finally, we leveraged in silico annotation and differential gene expression data to prioritize AD-related variants and genes. RESULTS: Admixture mapping identified two genome-wide significant loci on chromosomes 17p13.2 (p = 2.2 × 10-5 ) and 18q21.33 (p = 1.2 × 10-5 ). Our fine mapping of the chromosome 17p13.2 and 18q21.33 regions revealed several interesting genes such as the MINK1, KIF1C, and BCL2. DISCUSSION: Our ancestry-aware regression approach showed that AA individuals have a lower risk of AD if they inherited African ancestry admixture block at the 17p13.2 locus. HIGHLIGHTS: We identified two genome-wide significant admixture mapping signals: on chromosomes 17p13.2 and 18q21.33, which are novel in African American (AA) populations. Our ancestry-aware regression approach showed that AA individuals have a lower risk of Alzheimer's disease (AD) if they inherited African ancestry admixture block at the 17p13.2 locus. We found that the overall proportion of African ancestry does not differ between the cases and controls that suggest African genetic ancestry alone is not likely to explain the AD prevalence difference between AA and non-Hispanic White populations.

FMNL2 regulates gliovascular interactions and is associated with vascular risk factors and cerebrovascular pathology in Alzheimer's disease.

Alzheimer's disease (AD) has been associated with cardiovascular and cerebrovascular risk factors (CVRFs) during middle age and later and is frequently accompanied by cerebrovascular pathology at death. An interaction between CVRFs and genetic variants might explain the pathogenesis. Genome-wide, gene by CVRF interaction analyses for AD, in 6568 patients and 8101 controls identified FMNL2 (p = 6.6 × 10-7). A significant increase in FMNL2 expression was observed in the brains of patients with brain infarcts and AD pathology and was associated with amyloid and phosphorylated tau deposition. FMNL2 was also prominent in astroglia in AD among those with cerebrovascular pathology. Amyloid toxicity in zebrafish increased fmnl2a expression in astroglia with detachment of astroglial end feet from blood vessels. Knockdown of fmnl2a prevented gliovascular remodeling, reduced microglial activity and enhanced amyloidosis. APP/PS1dE9 AD mice also displayed increased Fmnl2 expression and reduced the gliovascular contacts independent of the gliotic response. Based on this work, we propose that FMNL2 regulates pathology-dependent plasticity of the blood-brain-barrier by controlling gliovascular interactions and stimulating the clearance of extracellular aggregates. Therefore, in AD cerebrovascular risk factors promote cerebrovascular pathology which in turn, interacts with FMNL2 altering the normal astroglial-vascular mechanisms underlying the clearance of amyloid and tau increasing their deposition in brain.

Admixture Mapping of Alzheimer's disease in Caribbean Hispanics identifies a new locus on 22q13.1.

Late-onset Alzheimer's disease (LOAD) is significantly more frequent in Hispanics than in non-Hispanic Whites. Ancestry may explain these differences across ethnic groups. To this end, we studied a large cohort of Caribbean Hispanics (CH, N = 8813) and tested the association between Local Ancestry (LA) and LOAD ("admixture mapping") to identify LOAD-associated ancestral blocks, separately for ancestral components (European [EUR], African [AFR], Native American[NA]) and jointly (AFR + NA). Ancestral blocks significant after permutation were fine-mapped employing multi-ethnic whole-exome sequencing (WES) to identify rare variants associated with LOAD (SKAT-O) and replicated in the UK Biobank WES dataset. Candidate genes were validated studying (A) protein expression in human LOAD and control brains; (B) two animal AD models, Drosophila and Zebrafish. In the joint AFR + NA model, we identified four significant ancestral blocks located on chromosomes 1 (p value = 8.94E-05), 6 (p value = 8.63E-05), 21 (p value = 4.64E-05) and 22 (p value = 1.77E-05). Fine-mapping prioritized the GCAT gene on chromosome 22 (SKAT-O p value = 3.45E-05) and replicated in the UK Biobank (SKAT-O p value = 0.05). In LOAD brains, a decrease of 28% in GCAT protein expression was observed (p value = 0.038), and GCAT knockdown in Amyloid-ß42 Drosophila exacerbated rough eye phenotype (68% increase, p value = 4.84E-09). In zebrafish, gcat expression increased after acute amyloidosis (34%, p value = 0.0049), and decreased upon anti-inflammatory Interleukin-4 (39%, p value = 2.3E-05). Admixture mapping uncovered genomic regions harboring new LOAD-associated loci that might explain the observed different frequency of LOAD across ethnic groups. Our results suggest that the inflammation-related activity of GCAT is a response to amyloid toxicity, and reduced GCAT expression exacerbates AD pathology.

Genetic determinants of intracranial large artery stenosis in the northern Manhattan study.

BACKGROUND: Intracranial stenosis is one of the most common causes of stroke worldwide. Several single nucleotide polymorphisms have been associated with intracranial atherosclerosis, which is inferred to be the most common underlying cause of intracranial large artery stenosis (ILAS). We previously reviewed known genetic variants related to ILAS in predominantly Asian cohorts, but their prevalence and role in ILAS among western multiethnic populations are uncertain. METHODS: We leveraged existing imaging and genetic data from the Northern Manhattan Study, a multiethnic prospective cohort study. Based on literature review, we selected adiponectin Q (ADIPOQ) rs2241767 and rs182052, ring finger protein 213 (RNF213) rs112735431, apolipoprotein E (APOE) rs429358, phosphodiesterase 4D (PDE4D) rs2910829, lipoprotein lipase (LPL) rs320, and aldosterone synthase (CYP11B2) rs1799998 variants as candidates to explore. We defined ILAS as luminal stenosis >50% in any intracranial large artery using time-of-flight magnetic resonance angiography (MRA). RESULTS: We included 1109 participants (mean age 70 ± 9 years, 70% Hispanic, 60% women) in this study. ILAS was identified in 81 (7%) NOMAS participants. Logistic regression analyses adjusted for age, sex, principal components, and vascular risk factors showed ILAS prevalence associated with CYP11B2 rs1799998 under the dominant model (OR = 0.56, 95%CI: 0.35-0.89) and LPL rs320 heterozygote genotype (OR = 1.68, 95%CI: 1.05-2.71). The genotype distributions of ADIPOQ rs2241767 and rs182052, APOE rs429358 and CYP11B2 rs1799998 variants were significantly different among non-Hispanic white and Black, and Hispanic groups. When participants were further stratified by race/ethnicity, the estimates were consistent for CYP11B2 rs1799998 across race/ethnic groups but not for LPL rs320. CONCLUSION: The CYP11B2 rs1799998 variant may be a protective genetic factor for ILAS across race/ethnic groups, but the risk of ILAS associated with LPL rs320 varies by race/ethnic group. Further functional studies may help elucidate the role that these variants play in the pathophysiology of ILAS.

Integration of GWAS and brain transcriptomic analyses in a multiethnic sample of 35,245 older adults identifies DCDC2 gene as predictor of episodic memory maintenance.

Identifying genes underlying memory function will help characterize cognitively resilient and high-risk declining subpopulations contributing to precision medicine strategies. We estimated episodic memory trajectories in 35,245 ethnically diverse older adults representing eight independent cohorts. We conducted apolipoprotein E (APOE)-stratified genome-wide association study (GWAS) analyses and combined individual cohorts' results via meta-analysis. Three independent transcriptomics datasets were used to further interpret GWAS signals. We identified DCDC2 gene significantly associated with episodic memory (Pmeta = 3.3 x 10-8 ) among non-carriers of APOE ε4 (N = 24,941). Brain transcriptomics revealed an association between episodic memory maintenance and (1) increased dorsolateral prefrontal cortex DCDC2 expression (P = 3.8 x 10-4 ) and (2) lower burden of pathological Alzheimer's disease (AD) hallmarks (paired helical fragment tau P = .003, and amyloid beta load P = .008). Additional transcriptomics results comparing AD and cognitively healthy brain samples showed a downregulation of DCDC2 levels in superior temporal gyrus (P = .007) and inferior frontal gyrus (P = .013). Our work identified DCDC2 gene as a novel predictor of memory maintenance.

KYNA/Ahr Signaling Suppresses Neural Stem Cell Plasticity and Neurogenesis in Adult Zebrafish Model of Alzheimer's Disease.

Neurogenesis decreases in Alzheimer's disease (AD) patients, suggesting that restoring the normal neurogenic response could be a disease modifying intervention. To study the mechanisms of pathology-induced neuro-regeneration in vertebrate brains, zebrafish is an excellent model due to its extensive neural regeneration capacity. Here, we report that Kynurenic acid (KYNA), a metabolite of the amino acid tryptophan, negatively regulates neural stem cell (NSC) plasticity in adult zebrafish brain through its receptor, aryl hydrocarbon receptor 2 (Ahr2). The production of KYNA is suppressed after amyloid-toxicity through reduction of the levels of Kynurenine amino transferase 2 (KAT2), the key enzyme producing KYNA. NSC proliferation is enhanced by an antagonist for Ahr2 and is reduced with Ahr2 agonists or KYNA. A subset of Ahr2-expressing zebrafish NSCs do not express other regulatory receptors such as il4r or ngfra, indicating that ahr2-positive NSCs constitute a new subset of neural progenitors that are responsive to amyloid-toxicity. By performing transcriptome-wide association studies (TWAS) in three late onset Alzheimer disease (LOAD) brain autopsy cohorts, we also found that several genes that are components of KYNA metabolism or AHR signaling are differentially expressed in LOAD, suggesting a strong link between KYNA/Ahr2 signaling axis to neurogenesis in LOAD.

Polygenic Risk Score for Alzheimer's Disease in Caribbean Hispanics.

OBJECTIVE: Polygenic risk scores (PRSs) assess the individual genetic propensity to a condition by combining sparse information scattered across genetic loci, often displaying small effect sizes. Most PRSs are constructed in European-ancestry populations, limiting their use in other ethnicities. Here we constructed and validated a PRS for late-onset Alzheimer's Disease (LOAD) in Caribbean Hispanics (CH). METHODS: We used a CH discovery (n = 4,312) and independent validation sample (n = 1,850) to construct an ancestry-specific PRS ("CH-PRS") and evaluated its performance alone and with other predictors using the area under curve (AUC) and logistic regression (strength of association with LOAD and statistical significance). We tested if CH-PRS predicted conversion to LOAD in a subsample with longitudinal data (n = 1,239). We also tested the CH-PRS in an independent replication CH cohort (n = 200) and brain autopsy cohort (n = 33). Finally, we tested the effect of ancestry on PRS by using European and African American discovery cohorts to construct alternative PRSs ("EUR-PRS", "AA-PRS"). RESULTS: The full model (LOAD ~ CH-PRS + sex + age + APOE-ɛ4), achieved an AUC = 74% (ORCH-PRS  = 1.51 95%CI = 1.36-1.68), raising to >75% in APOE-ɛ4 non-carriers. CH-PRS alone achieved an AUC = 72% in the autopsy cohort, raising to AUC = 83% in full model. Higher CH-PRS was significantly associated with clinical LOAD in the replication CH cohort (OR = 1.61, 95%CI = 1.19-2.17) and significantly predicted conversion to LOAD (HR = 1.93, CI = 1.70-2.20) in the longitudinal subsample. EUR-PRS and AA-PRS reached lower prediction accuracy (AUC = 58% and 53%, respectively). INTERPRETATION: Enriching diversity in genetic studies is critical to provide an effective PRS in profiling LOAD risk across populations. ANN NEUROL 2021;90:366-376.

Association of Variants in PINX1 and TREM2 With Late-Onset Alzheimer Disease.

IMPORTANCE: Genetic causes of late-onset Alzheimer disease (LOAD) are not completely explained by known genetic loci. Whole-exome and whole-genome sequencing can improve the understanding of the causes of LOAD and provide initial steps required to identify potential therapeutic targets. OBJECTIVE: To identify the genetic loci for LOAD across different ethnic groups. DESIGN, SETTING, AND PARTICIPANTS: This multicenter cohort study was designed to analyze whole-exome sequencing data from a multiethnic cohort using a transethnic gene-kernel association test meta-analysis, adjusted for sex, age, and principal components, to identify genetic variants associated with LOAD. A meta-analysis was conducted on the results of 2 independent studies of whole-exome and whole-genome sequence data from individuals of European ancestry. This group of European American, African American, and Caribbean Hispanic individuals participating in an urban population-based study were the discovery cohort; the additional cohorts included affected individuals and control participants from 2 publicly available data sets. Replication was achieved using independent data sets from Caribbean Hispanic families with multiple family members affected by LOAD and the International Genetics of Alzheimer Project. MAIN OUTCOMES AND MEASURES: Late-onset Alzheimer disease. RESULTS: The discovery cohort included 3595 affected individuals, while the additional cohorts included 5931 individuals with LOAD and 5504 control participants. Of 3916 individuals in the discovery cohort, we included 3595 individuals (1397 with LOAD and 2198 cognitively healthy controls; 2451 [68.2%] women; mean [SD] age, 80.3 [6.83] years). Another 321 individuals (8.2%) were excluded because of non-LOAD diagnosis, age younger than 60 years, missing covariates, duplicate data, or genetic outlier status. Gene-based tests that compared affected individuals (n = 7328) and control participants (n = 7702) and included only rare and uncommon variants annotated as having moderate-high functional effect supported PINX1 (8p23.1) as a locus with gene-wide significance (P = 2.81 × 10-6) after meta-analysis across the 3 studies. The PINX1 finding was replicated using data from the family-based study and the International Genetics of Alzheimer Project. Full meta-analysis of discovery and replication cohorts reached a P value of 6.16 × 10-7 for PINX1 (in 7620 affected individuals vs 7768 control participants). We also identified TREM2 in an annotation model that prioritized highly deleterious variants with a combined annotation dependent depletion greater than 20 (P= 1.0 × 10-7). CONCLUSIONS AND RELEVANCE: This gene-based, transethnic approach identified PINX1, a gene involved in telomere integrity, and TREM2, a gene with a product of an immune receptor found in microglia, as associated with LOAD. Both genes have well-established roles in aging and neurodegeneration.

Rare Variants Imputation in Admixed Populations: Comparison Across Reference Panels and Bioinformatics Tools.

BACKGROUND: Imputation has become a standard approach in genome-wide association studies (GWAS) to infer in silico untyped markers. Although feasibility for common variants imputation is well established, we aimed to assess rare and ultra-rare variants' imputation in an admixed Caribbean Hispanic population (CH). METHODS: We evaluated imputation accuracy in CH (N = 1,000), focusing on rare (0.1% ≤ minor allele frequency (MAF) ≤ 1%) and ultra-rare (MAF < 0.1%) variants. We used two reference panels, the Haplotype Reference Consortium (HRC; N = 27,165) and 1000 Genome Project (1000G phase 3; N = 2,504) and multiple phasing (SHAPEIT, Eagle2) and imputation algorithms (IMPUTE2, MACH-Admix). To assess imputation quality, we reported: (a) high-quality variant counts according to imputation tools' internal indexes (e.g., IMPUTE2 "Info" ≥ 80%). (b) Wilcoxon Signed-Rank Test comparing imputation quality for genotyped variants that were masked and imputed; (c) Cohen's kappa coefficient to test agreement between imputed and whole-exome sequencing (WES) variants; (d) imputation of G206A mutation in the PSEN1 (ultra-rare in the general population an more frequent in CH) followed by confirmation genotyping. We also tested ancestry proportion (European, African and Native American) against WES-imputation mismatches in a Poisson regression fashion. RESULTS: SHAPEIT2 retrieved higher percentage of imputed high-quality variants than Eagle2 (rare: 51.02% vs. 48.60%; ultra-rare 0.66% vs. 0.65%, Wilcoxon p-value < 0.001). SHAPEIT-IMPUTE2 employing HRC outperformed 1000G (64.50% vs. 59.17%; 1.69% vs. 0.75% for high-quality rare and ultra-rare variants, respectively, Wilcoxon p-value < 0.001). SHAPEIT-IMPUTE2 outperformed MaCH-Admix. Compared to 1000G, HRC-imputation retrieved a higher number of high-quality rare and ultra-rare variants, despite showing lower agreement between imputed and WES variants (e.g., rare: 98.86% for HRC vs. 99.02% for 1000G). High Kappa (K = 0.99) was observed for both reference panels. Twelve G206A mutation carriers were imputed and all validated by confirmation genotyping. African ancestry was associated with higher imputation errors for uncommon and rare variants (p-value < 1e-05). CONCLUSION: Reference panels with larger numbers of haplotypes can improve imputation quality for rare and ultra-rare variants in admixed populations such as CH. Ethnic composition is an important predictor of imputation accuracy, with higher African ancestry associated with poorer imputation accuracy.

Escherichia coli ST131-H22 as a Foodborne Uropathogen.

Escherichia coli sequence type 131 (ST131) has emerged rapidly to become the most prevalent extraintestinal pathogenic E. coli clones in circulation today. Previous investigations appeared to exonerate retail meat as a source of human exposure to ST131; however, these studies focused mainly on extensively multidrug-resistant ST131 strains, which typically carry allele 30 of the fimH type 1 fimbrial adhesin gene (ST131-H30). To estimate the frequency of extraintestinal human infections arising from foodborne ST131 strains without bias toward particular sublineages or phenotypes, we conducted a 1-year prospective study of E. coli from meat products and clinical cultures in Flagstaff, Arizona. We characterized all isolates by multilocus sequence typing, fimH typing, and core genome phylogenetic analyses, and we screened isolates for avian-associated ColV plasmids as an indication of poultry adaptation. E. coli was isolated from 79.8% of the 2,452 meat samples and 72.4% of the 1,735 culture-positive clinical samples. Twenty-seven meat isolates were ST131 and belonged almost exclusively (n = 25) to the ST131-H22 lineage. All but 1 of the 25 H22 meat isolates were from poultry products, and all but 2 carried poultry-associated ColV plasmids. Of the 1,188 contemporaneous human clinical E. coli isolates, 24 were ST131-H22, one-quarter of which occurred in the same high-resolution phylogenetic clades as the ST131-H22 meat isolates and carried ColV plasmids. Molecular clock analysis of an international ST131-H22 genome collection suggested that ColV plasmids have been acquired at least six times since the 1940s and that poultry-to-human transmission is not limited to the United States.IMPORTANCEE. coli ST131 is an important extraintestinal pathogen that can colonize the gastrointestinal tracts of humans and food animals. Here, we combined detection of accessory traits associated with avian adaptation (ColV plasmids) with high-resolution phylogenetics to quantify the portion of human infections caused by ST131 strains of food animal origin. Our results suggest that one ST131 sublineage-ST131-H22-has become established in poultry populations around the world and that meat may serve as a vehicle for human exposure and infection. ST131-H22 is just one of many E. coli lineages that may be transmitted from food animals to humans. Additional studies that combine detection of host-associated accessory elements with phylogenetics may allow us to quantify the total fraction of human extraintestinal infections attributable to food animal E. coli strains.

Penile Anaerobic Dysbiosis as a Risk Factor for HIV Infection.

Sexual transmission of HIV requires exposure to the virus and infection of activated mucosal immune cells, specifically CD4+ T cells or dendritic cells. The foreskin is a major site of viral entry in heterosexual transmission of HIV. Although the probability of acquiring HIV from a sexual encounter is low, the risk varies even after adjusting for known HIV risk factors. The genital microbiome may account for some of the variability in risk by interacting with the host immune system to trigger inflammatory responses that mediate the infection of mucosal immune cells. We conducted a case-control study of uncircumcised participants nested within a randomized-controlled trial of male circumcision in Rakai, Uganda. Using penile (coronal sulcus) swabs collected by study personnel at trial enrollment, we characterized the penile microbiome by sequencing and real-time PCR and cytokine levels by electrochemiluminescence assays. The absolute abundances of penile anaerobes at enrollment were associated with later risk of HIV seroconversion, with a 10-fold increase in Prevotella, Dialister, Finegoldia, and Peptoniphilus increasing the odds of HIV acquisition by 54 to 63%, after controlling for other known HIV risk factors. Increased abundances of anaerobic bacteria were also correlated with increased cytokines, including interleukin-8, which can trigger an inflammatory response that recruits susceptible immune cells, suggesting a mechanism underlying the increased risk. These same anaerobic genera can be shared between heterosexual partners and are associated with increased HIV acquisition in women, pointing to anaerobic dysbiosis in the genital microbiome and an accompanying inflammatory response as a novel, independent, and transmissible risk factor for HIV infection.IMPORTANCE We found that uncircumcised men who became infected by HIV during a 2-year clinical trial had higher levels of penile anaerobes than uncircumcised men who remained HIV negative. We also found that having higher levels of penile anaerobes was also associated with higher production of immune factors that recruit HIV target cells to the foreskin, suggesting that anaerobes may modify HIV risk by triggering inflammation. These anaerobes are known to be shared by heterosexual partners and are associated with HIV risk in women. Therefore, penile anaerobes may be a sexually transmissible risk factor for HIV, and modifying the penile microbiome could potentially reduce HIV acquisition in both men and women.