Genome-wide Association Study Identifies HLA-DPB1 as a Significant Risk Factor for Severe Aplastic Anemia.

Severe aplastic anemia (SAA) is a rare disorder characterized by hypoplastic bone marrow and progressive pancytopenia. The etiology of acquired SAA is not understood but is likely related to abnormal immune responses and environmental exposures. We conducted a genome-wide association study of individuals with SAA genetically matched to healthy controls in discovery (359 cases, 1,396 controls) and validation sets (175 cases, 1,059 controls). Combined analyses identified linked SNPs in distinct blocks within the major histocompatibility complex on 6p21. The top SNP encodes p.Met76Val in the P4 binding pocket of the HLA class II gene HLA-DPB1 (rs1042151A>G, odds ratio [OR] 1.75, 95% confidence interval [CI] 1.50-2.03, p = 1.94 × 10-13) and was associated with HLA-DP cell surface expression in healthy individuals (p = 2.04 × 10-6). Phylogenetic analyses indicate that Val76 is not monophyletic and likely occurs in conjunction with different HLA-DP binding groove conformations. Imputation of HLA-DPB1 alleles revealed increased risk of SAA associated with Val76-encoding alleles DPB1∗03:01, (OR 1.66, p = 1.52 × 10-7), DPB1∗10:01 (OR 2.12, p = 0.0003), and DPB1∗01:01 (OR 1.60, p = 0.0008). A second SNP near HLA-B, rs28367832G>A, reached genome-wide significance (OR 1.49, 95% CI 1.22-1.78, p = 7.27 × 10-9) in combined analyses; the association remained significant after excluding cases with clonal copy-neutral loss-of-heterozygosity affecting class I HLA genes (8.6% of cases and 0% of controls). SNPs in the HLA class II gene HLA-DPB1 and possibly class I (HLA-B) are associated with SAA. The replacement of Met76 to Val76 in certain HLA-DPB1 alleles might influence risk of SAA through mechanisms involving DP peptide binding specificity, expression, and/or other factors affecting DP function.

Distinct frequency patterns of LILRB3 and LILRA6 allelic variants in Europeans.

The leukocyte immunoglobulin-like receptor (LILR)B3 and LILRA6 genes encode homologous myeloid inhibitory and activating orphan receptors, respectively. Both genes exhibit a strikingly high level of polymorphism at the amino acid level and LILRA6 (but not LILRB3) displays copy number variation (CNV). Although multiple alleles have been reported for both genes, limited data is available on frequencies of these alleles among humans. We have sequenced LILRB3/A6 exons encoding signal peptides and ectodomains in 91 healthy blood donors of European descent who carry one or two copies of LILRA6 per diploid genome. Analysis of haplotypes among individuals with two LILRA6 copies, representing the majority in this cohort (N = 86), shows that common LILRB3 and LILRA6 alleles encode some distinct amino acid sequences in homologous regions of the receptors, which could potentially impact their respective functions differentially. Comparison of sequences in individuals with one vs. two copies of LILRA6 supports non-allelic homologous recombination between LILRB3 and LILRA6 as a mechanism for generating LILRA6 CNV and LILRB3 diversity. These data characterize LILRB3/LILRA6 genetic variation in more detail than previously described and underscore the need to determine their ligands.

HLA tapasin independence: broader peptide repertoire and HIV control.

Human leukocyte antigen (HLA) class I allotypes vary in their ability to present peptides in the absence of tapasin, an essential component of the peptide loading complex. We quantified tapasin dependence of all allotypes that are common in European and African Americans (n = 97), which revealed a broad continuum of values. Ex vivo examination of cytotoxic T cell responses to the entire HIV-1 proteome from infected subjects indicates that tapasin-dependent allotypes present a more limited set of distinct peptides than do tapasin-independent allotypes, data supported by computational predictions. This suggests that variation in tapasin dependence may impact the strength of the immune responses by altering peptide repertoire size. In support of this model, we observed that individuals carrying HLA class I genotypes characterized by greater tapasin independence progress more slowly to AIDS and maintain lower viral loads, presumably due to increased breadth of peptide presentation. Thus, tapasin dependence level, like HLA zygosity, may serve as a means to restrict or expand breadth of the HLA-I peptide repertoire across humans, ultimately influencing immune responses to pathogens and vaccines.

Distinct contributions of cathelin-related antimicrobial peptide (CRAMP) derived from epithelial cells and macrophages to colon mucosal homeostasis.

The cathelin-related antimicrobial peptide CRAMP protects the mouse colon from inflammation, inflammation-associated carcinogenesis, and disrupted microbiome balance, as shown in systemic Cnlp-/- mice (also known as Camp-/- mice). However, the mechanistic basis for the role and the cellular source of CRAMP in colon pathophysiology are ill defined. This study, using either epithelial or myeloid conditional Cnlp-/- mice, demonstrated that epithelial cell-derived CRAMP played a major role in supporting normal development of colon crypts, mucus production, and repair of injured mucosa. On the other hand, myeloid cell-derived CRAMP potently supported colon epithelial resistance to bacterial invasion during acute inflammation with exacerbated mucosal damage and higher rate of mouse mortality. Therefore, a well concerted cooperation of epithelial- and myeloid-derived CRAMP is essential for colon mucosal homeostasis. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Population-specific diversity of the immunoglobulin constant heavy G chain (IGHG) genes.

Human immunoglobulin G (IgG) molecules, IgG1, IgG2 and IgG3, exhibit substantial inter-individual variation in their constant heavy chain regions, as discovered by serological methods. This polymorphism is encoded by the IGHG1, IGHG2, and IGHG3 genes and may influence antibody function. We sequenced the coding fragments of these genes in 95 European Americans, 94 African Americans, and 94 Black South Africans. Striking differences were observed between the population groups, including extremely low amino acid sequence variation in IGHG1 among South Africans, and higher IGHG2 and IGHG3 diversity in individuals of African descent compared to individuals of European descent. Molecular definition of the loci illustrates a greater level of allelic polymorphism than previously described, including the presence of common IGHG2 and IGHG3 variants that were indistinguishable serologically. Comparison of our data with the 1000 Genome Project sequences indicates overall agreement between the datasets, although some inaccuracies in the 1000 Genomes Project are likely. These data represent the most comprehensive analysis of IGHG polymorphisms across major populations, which can now be applied to deciphering their functional impact.

A Prime/Boost Vaccine Regimen Alters the Rectal Microbiome and Impacts Immune Responses and Viremia Control Post-Simian Immunodeficiency Virus Infection in Male and Female Rhesus Macaques.

An efficacious human immunodeficiency virus (HIV) vaccine will likely require induction of both mucosal and systemic immune responses. We compared the immunogenicity and protective efficacy of two mucosal/systemic vaccine regimens and investigated their effects on the rectal microbiome. Rhesus macaques were primed twice mucosally with replication-competent adenovirus type 5 host range mutant (Ad5hr)-simian immunodeficiency virus (SIV) recombinants and boosted twice intramuscularly with ALVAC-SIV recombinant plus SIV gp120 protein or with DNA for SIV genes and rhesus interleukin-12 plus SIV gp120 protein. Controls received empty Ad5hr vector and alum adjuvant only. Both regimens elicited strong, comparable mucosal and systemic cellular and humoral immunity. Prevaccination rectal microbiomes of males and females differed and significantly changed over the course of immunization, most strongly in females after Ad5hr immunizations. Following repeated low-dose intrarectal SIV challenges, both vaccine groups exhibited modestly but significantly reduced acute viremia. Male and female controls exhibited similar acute viral loads; however, vaccinated females, but not males, exhibited lower levels of acute viremia, compared to same-sex controls. Few differences in adaptive immune responses were observed between the sexes. Striking differences in correlations of the rectal microbiome of males and females with acute viremia and immune responses associated with protection were seen and point to effects of the microbiome on vaccine-induced immunity and viremia control. Our study clearly demonstrates direct effects of a mucosal SIV vaccine regimen on the rectal microbiome and validates our previously reported SIV vaccine-induced sex bias. Sex and the microbiome are critical factors that should not be overlooked in vaccine design and evaluation.IMPORTANCE Differences in HIV pathogenesis between males and females, including immunity postinfection, have been well documented, as have steroid hormone effects on the microbiome, which is known to influence mucosal immune responses. Few studies have applied this knowledge to vaccine trials. We investigated two SIV vaccine regimens combining mucosal priming immunizations and systemic protein boosting. We again report a vaccine-induced sex bias, with female rhesus macaques but not males displaying significantly reduced acute viremia. The vaccine regimens, especially the mucosal primes, significantly altered the rectal microbiome. The greatest effects were in females. Striking differences between female and male macaques in correlations of prevalent rectal bacteria with viral loads and potentially protective immune responses were observed. Effects of the microbiome on vaccine-induced immunity and viremia control require further study by microbiome transfer. However, the findings presented highlight the critical importance of considering effects of sex and the microbiome in vaccine design and evaluation.

The Antimicrobial Peptide CRAMP Is Essential for Colon Homeostasis by Maintaining Microbiota Balance.

Commensal bacteria are critical for physiological functions in the gut, and dysbiosis in the gut may cause diseases. In this article, we report that mice deficient in cathelin-related antimicrobial peptide (CRAMP) were defective in the development of colon mucosa and highly sensitive to dextran sulfate sodium (DSS)-elicited colitis, as well as azoxymethane-mediated carcinogenesis. Pretreatment of CRAMP-/- mice with antibiotics markedly reduced the severity of DSS-induced colitis, suggesting CRAMP as a limiting factor on dysbiosis in the colon. This was supported by observations that wild-type (WT) mice cohoused with CRAMP-/- mice became highly sensitive to DSS-induced colitis, and the composition of fecal microbiota was skewed by CRAMP deficiency. In particular, several bacterial species that are typically found in oral microbiota, such as Mogibacterium neglectum, Desulfovibrio piger, and Desulfomicrobium orale, were increased in feces of CRAMP-/- mice and were transferred to WT mice during cohousing. When littermates of CRAMP+/- parents were examined, the composition of the fecal microbiota of WT pups and heterozygous parents was similar. In contrast, although the difference in fecal microbiota between CRAMP-/- and WT pups was small early on after weaning and single mouse housing, there was an increasing divergence with prolonged single housing. These results indicate that CRAMP is critical in maintaining colon microbiota balance and supports mucosal homeostasis, anti-inflammatory responses, and protection from carcinogenesis.

Posttranscriptional Regulation of HLA-A Protein Expression by Alternative Polyadenylation Signals Involving the RNA-Binding Protein Syncrip.

Genomic variation in the untranslated region (UTR) has been shown to influence HLA class I expression level and associate with disease outcomes. Sequencing of the 3'UTR of common HLA-A alleles indicated the presence of two polyadenylation signals (PAS). The proximal PAS is conserved, whereas the distal PAS is disrupted within certain alleles by sequence variants. Using 3'RACE, we confirmed expression of two distinct forms of the HLA-A 3'UTR based on use of either the proximal or the distal PAS, which differ in length by 100 bp. Specific HLA-A alleles varied in the usage of the proximal versus distal PAS, with some alleles using only the proximal PAS, and others using both the proximal and distal PAS to differing degrees. We show that the short and the long 3'UTR produced similar mRNA expression levels. However, the long 3'UTR conferred lower luciferase activity as compared with the short form, indicating translation inhibition of the long 3'UTR. RNA affinity pull-down followed by mass spectrometry analysis as well as RNA coimmunoprecipitation indicated differential binding of Syncrip to the long versus short 3'UTR. Depletion of Syncrip by small interfering RNA increased surface expression of an HLA-A allotype that uses primarily the long 3'UTR, whereas an allotype expressing only the short form was unaffected. Furthermore, specific blocking of the proximal 3'UTR reduced surface expression without decreasing mRNA expression. These data demonstrate HLA-A allele-specific variation in PAS usage, which modulates their cell surface expression posttranscriptionally.

Sequence and Phylogenetic Analysis of the Untranslated Promoter Regions for HLA Class I Genes.

Polymorphisms located within the MHC have been linked to many disease outcomes by mechanisms not yet fully understood in most cases. Variants located within untranslated regions of HLA genes are involved in allele-specific expression and may therefore underlie some of these disease associations. We determined sequences extending nearly 2 kb upstream of the transcription start site for 68 alleles from 57 major lineages of classical HLA class I genes. The nucleotide diversity within this promoter segment roughly follows that seen within the coding regions, with HLA-B showing the highest (∼1.9%), followed by HLA-A (∼1.8%), and HLA-C showing the lowest diversity (∼0.9%). Despite its greater diversity, HLA-B mRNA expression levels determined in 178 European Americans do not vary in an allele- or lineage-specific manner, unlike the differential expression levels of HLA-A or HLA-C reported previously. Close proximity of promoter sequences in phylogenetic trees is roughly reflected by similarity of expression pattern for most HLA-A and -C loci. Although promoter sequence divergence might impact promoter activity, we observed no clear link between the phylogenetic structures as represented by pairwise nucleotide differences in the promoter regions with estimated differences in mRNA expression levels for the classical class I loci. Further, no pair of class I loci showed coordinated expression levels, suggesting that distinct mechanisms across loci determine their expression level under nonstimulated conditions. These data serve as a foundation for more in-depth analysis of the functional consequences of promoter region variation within the classical HLA class I loci.

High HLA-DP Expression and Graft-versus-Host Disease.

BACKGROUND: Transplantation of hematopoietic cells from unrelated donors can cure blood disorders but carries a significant risk of acute graft-versus-host disease (GVHD). The risk is higher when the recipient and donor are HLA-DPB1-mismatched, but the mechanisms leading to GVHD are unknown. The HLA-DPB1 regulatory region variant rs9277534 is associated with HLA-DPB1 expression. We tested the hypothesis that the GVHD risk correlates with the rs9277534 allele linked to the mismatched HLA-DPB1 in the recipient. METHODS: We genotyped rs9277534 in 3505 persons to define rs9277534-DPB1 haplotypes. Among 1441 recipients of transplants from HLA-A,B,C,DRB1,DQB1-matched unrelated donors with only one HLA-DPB1 mismatch, linkage of the rs9277534 A and G alleles to the mismatched HLA-DPB1 was determined. HLA-DPB1 expression was assessed by means of a quantitative polymerase-chain-reaction assay. The risk of acute GVHD among recipients whose mismatched HLA-DPB1 allele was linked to rs9277534G (high expression) was compared with the risk among recipients whose mismatched HLA-DPB1 allele was linked to rs9277534A (low expression). RESULTS: The mean HLA-DPB1 expression was lower with rs9277534A than with rs9277534G. Among recipients of transplants from donors with rs9277534A-linked HLA-DPB1, the risk of acute GVHD was higher for recipients with rs9277534G-linked HLA-DPB1 mismatches than for recipients with rs9277534A-linked HLA-DPB1 mismatches (hazard ratio, 1.54; 95% confidence interval [CI], 1.25 to 1.89; P<0.001), as was the risk of death due to causes other than disease recurrence (hazard ratio, 1.25; 95% CI, 1.00 to 1.57; P=0.05). CONCLUSIONS: The risk of GVHD associated with HLA-DPB1 mismatching was influenced by the HLA-DPB1 rs9277534 expression marker. Among recipients of HLA-DPB1-mismatched transplants from donors with the low-expression allele, recipients with the high-expression allele had a high risk of GVHD. (Funded by the National Institutes of Health and others.).

Epigenetic regulation of differential HLA-A allelic expression levels.

MHC class I expression levels influence the strength of immune responses and represent another variable in determining outcome to disease beyond peptide binding alone. Identification of the HLA loci that vary in allelic expression levels and delineating the mechanism responsible for expression variation may provide the opportunity to modify their expression therapeutically. We have examined the expression levels of allelic lineages at the HLA-A locus in a sample of 216 European Americans using a real-time polymerase chain reaction assay, which amplifies all HLA-A lineages specifically with equal efficiency, and observed a gradient of expression that associates with HLA-A allelic lineage (R = 0.6, P = 5 × 10(-25)). DNA methylation of the HLA-A gene appears to contribute to the variation in HLA-A mRNA expression levels, as a significant inverse correlation was observed between HLA-A mRNA expression levels in untreated cells and the degree to which expression is increased after treatment of the cells with a DNA methyltransferase inhibitor (R = 0.6, P = 2.8 × 10(-6)). Further, deep-sequencing and immunoprecipitation assays revealed allelic lineage-specific methylation patterns within the HLA-A promoter region where increased DNA methylation levels correlated significantly with reduced HLA-A expression levels (R = 0.89, P = 3.7 × 10(-9)). These data demonstrate HLA-A allelic lineage-specific variation in expression levels, and DNA methylation as a likely factor in contributing to this variation.

Significance of positive selection and gene duplication in adaptive evolution: in memory of Austin L. Hughes.

In this minireview, we highlight the contributions of the late Austin L. Hughes to two areas of molecular evolution: the role of positive (Darwinian) selection, and the impact of gene duplications during genome evolution.

HLA-C expression levels define permissible mismatches in hematopoietic cell transplantation.

Life-threatening graft-versus-host disease (GVHD) limits the use of HLA-C-mismatched unrelated donors in transplantation. Clinicians lack criteria for donor selection when HLA-C-mismatched donors are a patient's only option for cure. We examined the role for HLA-C expression levels to identify permissible HLA-C mismatches. The median fluorescence intensity, a proxy of HLA-C expression, was assigned to each HLA-C allotype in 1975 patients and their HLA-C-mismatched unrelated transplant donors. The association of outcome with the level of expression of patients' and donors' HLA-C allotypes was evaluated in multivariable models. Increasing expression level of the patient's mismatched HLA-C allotype was associated with increased risks of grades III to IV acute GVHD, nonrelapse mortality, and mortality. Increasing expression level among HLA-C mismatches with residue 116 or residue 77/80 mismatching was associated with increased nonrelapse mortality. The immunogenicity of HLA-C mismatches in unrelated donor transplantation is influenced by the expression level of the patient's mismatched HLA-C allotype. HLA-C expression levels provide new information on mismatches that should be avoided and extend understanding of HLA-C-mediated immune responses in human disease.

Genetic interplay between HLA-C and MIR148A in HIV control and Crohn disease.

Variation in the 3' untranslated region (3'UTR) of the HLA-C locus determines binding of the microRNA Hsa-miR-148a, resulting in lower cell surface expression of alleles that bind miR-148a relative to those alleles that escape its binding. The HLA-C 3'UTR variant was shown to associate with HIV control, but like the vast majority of disease associations in a region dense with causal candidates, a direct effect of HLA-C expression level on HIV control was not proven. We demonstrate that a MIR148A insertion/deletion polymorphism associates with its own expression levels, affecting the extent to which HLA-C is down-regulated, the level of HIV control, and the risk of Crohn disease only among those carrying an intact miR-148a binding site in the HLA-C 3'UTR. These data illustrate a direct effect of HLA-C expression level on HIV control that cannot be attributed to other HLA loci in linkage disequilibrium with HLA-C and highlight the rich complexity of genetic interactions in human disease.

Sequencing rare and common APOL1 coding variants to determine kidney disease risk.

A third of African Americans with sporadic focal segmental glomerulosclerosis (FSGS) or HIV-associated nephropathy (HIVAN) do not carry APOL1 renal risk genotypes. This raises the possibility that other APOL1 variants may contribute to kidney disease. To address this question, we sequenced all APOL1 exons in 1437 Americans of African and European descent, including 464 patients with biopsy-proven FSGS/HIVAN. Testing for association with 33 common and rare variants with FSGS/HIVAN revealed no association independent of strong recessive G1 and G2 effects. Seeking additional variants that might have been under selection by pathogens and could represent candidates for kidney disease risk, we also sequenced an additional 1112 individuals representing 53 global populations. Except for G1 and G2, none of the 7 common codon-altering variants showed evidence of selection or could restore lysis against trypanosomes causing human African trypanosomiasis. Thus, only APOL1 G1 and G2 confer renal risk, and other common and rare APOL1 missense variants, including the archaic G3 haplotype, do not contribute to sporadic FSGS and HIVAN in the US population. Hence, in most potential clinical or screening applications, our study suggests that sequencing APOL1 exons is unlikely to bring additional information compared to genotyping only APOL1 G1 and G2 risk alleles.

Genetic variation in IL28B and spontaneous clearance of hepatitis C virus.

Hepatitis C virus (HCV) infection is the most common blood-borne infection in the United States, with estimates of 4 million HCV-infected individuals in the United States and 170 million worldwide. Most (70-80%) HCV infections persist and about 30% of individuals with persistent infection develop chronic liver disease, including cirrhosis and hepatocellular carcinoma. Epidemiological, viral and host factors have been associated with the differences in HCV clearance or persistence, and studies have demonstrated that a strong host immune response against HCV favours viral clearance. Thus, variation in genes involved in the immune response may contribute to the ability to clear the virus. In a recent genome-wide association study, a single nucleotide polymorphism (rs12979860) 3 kilobases upstream of the IL28B gene, which encodes the type III interferon IFN-3, was shown to associate strongly with more than a twofold difference in response to HCV drug treatment. To determine the potential effect of rs12979860 variation on outcome to HCV infection in a natural history setting, we genotyped this variant in HCV cohorts comprised of individuals who spontaneously cleared the virus (n = 388) or had persistent infection (n = 620). We show that the C/C genotype strongly enhances resolution of HCV infection among individuals of both European and African ancestry. To our knowledge, this is the strongest and most significant genetic effect associated with natural clearance of HCV, and these results implicate a primary role for IL28B in resolution of HCV infection.

The molecular origin and consequences of escape from miRNA regulation by HLA-C alleles.

Differential expression of human leukocyte antigen C (HLA-C) allotypes is mediated by the binding of a microRNA, miR-148a, to the 3' untranslated region of some, but not all, HLA-C alleles. The binding results in lower levels of HLA-C expression, which is associated with higher levels of HIV-1 viral load among infected individuals. The alternative set of HLA-C alleles has several substitutions in the miR-148a binding site that prevent binding and HLA-C downregulation; these high-expression alleles associate with control of HIV-1 viral load. We show that the common ancestor of all extant HLA-C alleles was suppressed by miR-148a. Substitutions that prevent miR-148a binding arose by a sequence exchange event between an HLA-C allele and an HLA-B (MIM 142830) allele of a B(∗)07-like lineage. The event occurred 3-5 million years ago, resulting in an HLA-C variant that escape from miR-148a downregulation. We present evidence suggesting that selection played a role in the successful spread of the HLA-C escape alleles, giving rise to 7 of the 14 extant HLA-C lineages. Notably, critical peptide and KIR binding residues of the escape variants have selectively converged to resemble the sequence of their inhibited counterparts, such that the inhibited and escape groupings differ primarily by their levels of expression.

Impact of HLA class I functional divergence on HIV control.

Heterozygosity of Human leukocyte antigen (HLA) class I genes is linked to beneficial outcomes after HIV infection, presumably through greater breadth of HIV epitope presentation and cytotoxic T cell response. Distinct allotype pairs, however, differ in the extent to which they bind shared sets of peptides. We developed a functional divergence metric that measures pairwise complementarity of allotype-associated peptide binding profiles. Greater functional divergence for pairs of HLA-A and/or HLA-B allotypes was associated with slower AIDS progression and independently with enhanced viral load control. The metric predicts immune breadth at the peptide level rather than gene level and redefines HLA heterozygosity as a continuum differentially affecting disease outcome. Functional divergence may affect response to additional infections, vaccination, immunotherapy, and other diseases where HLA heterozygote advantage occurs.

The beneficial effects of commensal E. coli for colon epithelial cell recovery are related with Formyl peptide receptor 2 (Fpr2) in epithelial cells.

BACKGROUND: Formyl peptide receptor 2 (Fpr2) plays a crucial role in colon homeostasis and microbiota balance. Commensal E. coli is known to promote the regeneration of damaged colon epithelial cells. The aim of the study was to investigate the connection between E. coli and Fpr2 in the recovery of colon epithelial cells. RESULTS: The deficiency of Fpr2 was associated with impaired integrity of the colon mucosa and an imbalance of microbiota, characterized by the enrichment of Proteobacteria in the colon. Two serotypes of E. coli, O22:H8 and O91:H21, were identified in the mouse colon through complete genome sequencing. E. coli O22:H8 was found to be prevalent in the gut of mice and exhibited lower virulence compared to O91:H21. Germ-free (GF) mice that were pre-orally inoculated with E. coli O22:H8 showed reduced susceptibility to chemically induced colitis, increased proliferation of epithelial cells, and improved mouse survival. Following infection with E. coli O22:H8, the expression of Fpr2 in colon epithelial cells was upregulated, and the products derived from E. coli O22:H8 induced migration and proliferation of colon epithelial cells through Fpr2. Fpr2 deficiency increased susceptibility to chemically induced colitis, delayed the repair of damaged colon epithelial cells, and heightened inflammatory responses. Additionally, the population of E. coli was observed to increase in the colons of Fpr2-/- mice with colitis. CONCLUSION: Commensal E. coli O22:H8 stimulated the upregulation of Fpr2 expression in colon epithelial cells, and the products from E. coli induced migration and proliferation of colon epithelial cells through Fpr2. Fpr2 deficiency led to an increased E. coli population in the colon and delayed recovery of damaged colon epithelial cells in mice with colitis. Therefore, Fpr2 is essential for the effects of commensal E. coli on colon epithelial cell recovery.