Immune gene expression changes more during a malaria transmission season than between consecutive seasons.

Plasmodium parasites, the causative organism of malaria, caused over 600,000 deaths in 2022. In Mali, Plasmodium falciparum causes the majority of malaria cases and deaths and is transmitted seasonally. Anti-malarial immunity develops slowly over repeated exposures to P. falciparum and some aspects of this immunity (e.g., antibody titers) wane during the non-transmission, dry season. Here, we sequenced RNA from 33 pediatric blood samples collected during P. falciparum infections at the beginning or end of a transmission season, and characterized the host and parasite gene expression profiles for paired, consecutive infections. We found that human gene expression changes more over the course of one transmission season than between seasons, with signatures of partial development of an adaptive immune response during one transmission season and stability in gene expression during the dry season. Additionally, we found that P. falciparum gene expression did not vary with timing during the season and remained stable both across and between seasons, despite varying human immune pressures. Our results provide insights into the dynamics of anti-malarial immune response development over short time frames that could be exploited by future vaccine and prevention efforts. IMPORTANCE: Our work seeks to understand how the immune response to Plasmodium falciparum malaria changes between infections that occur during low and high malaria transmission seasons, and highlights that immune gene expression changes more during the high transmission season. This provides important insight into the dynamics of the anti-malarial immune response that are important to characterize over these short time frames to better understand how to exploit this immune response with future vaccine efforts.

Children with hemoglobin C or S trait have low serologic responses to a subset of malaria variant surface antigens.

Children with hemoglobin AC or AS have decreased susceptibility to clinical malaria. Parasite variant surface antigen (VSA) presentation on the surface of infected erythrocytes is altered in erythrocytes with hemoglobin C (Hb AC) or sickle trait (Hb AS) mutations in vitro. The protective role of incomplete or altered VSA presentation against clinical malaria in individuals with Hb AC or AS is unclear. Using a high-throughput protein microarray, we sought to use serological responses to VSAs as a measure of host exposure to VSAs among Malian children with Hb AC, Hb AS, or wildtype hemoglobin (Hb AA). In uncomplicated malaria, when compared to Hb AA children, Hb AC children had significantly lower serological responses to extracellular Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) domains but did not differ in responses to intracellular PfEMP1 domains and other VSAs, including members of the repetitive interspersed family (RIFIN) and subtelomeric variable open reading frame (STEVOR) family. Healthy children with Hb AC and Hb AS genotypes recognized fewer extracellular PfEMP1s compared to children with Hb AA, especially CD36-binding PfEMP1s. These reduced serologic responses may reflect reduced VSA presentation or lower parasite exposure in children with Hb AC or AS and provide insights into mechanisms of protection.

Gene expression analyses reveal differences in children's response to malaria according to their age.

In Bandiagara, Mali, children experience on average two clinical malaria episodes per year. However, even in the same transmission area, the number of uncomplicated symptomatic infections, and their parasitemia, can vary dramatically among children. We simultaneously characterize host and parasite gene expression profiles from 136 Malian children with symptomatic falciparum malaria and examine differences in the relative proportion of immune cells and parasite stages, as well as in gene expression, associated with infection and or patient characteristics. Parasitemia explains much of the variation in host and parasite gene expression, and infections with higher parasitemia display proportionally more neutrophils and fewer T cells, suggesting parasitemia-dependent neutrophil recruitment and/or T cell extravasation to secondary lymphoid organs. The child's age also strongly correlates with variations in gene expression: Plasmodium falciparum genes associated with age suggest that older children carry more male gametocytes, while variations in host gene expression indicate a stronger innate response in younger children and stronger adaptive response in older children. These analyses highlight the variability in host responses and parasite regulation during P. falciparum symptomatic infections and emphasize the importance of considering the children's age when studying and treating malaria infections.

Gene expression analyses reveal differences in children's response to malaria according to their age.

In Bandiagara, Mali, children experience on average two clinical malaria episodes per season. However, even in the same transmission area, the number of uncomplicated symptomatic infections, and their parasitemia, vary dramatically among children. To examine the factors contributing to these variations, we simultaneously characterized the host and parasite gene expression profiles from 136 children with symptomatic falciparum malaria and analyzed the expression of 9,205 human and 2,484 Plasmodium genes. We used gene expression deconvolution to estimate the relative proportion of immune cells and parasite stages in each sample and to adjust the differential gene expression analyses. Parasitemia explained much of the variation in both host and parasite gene expression and revealed that infections with higher parasitemia had more neutrophils and fewer T cells, suggesting parasitemia-dependent neutrophil recruitment and/or T cell extravasation to secondary lymphoid organs. The child's age was also strongly correlated with gene expression variations. Plasmodium falciparum genes associated with age suggested that older children carried more male gametocytes, while host genes associated with age indicated a stronger innate response (through TLR and NLR signaling) in younger children and stronger adaptive immunity (through TCR and BCR signaling) in older children. These analyses highlight the variability in host responses and parasite regulation during P. falciparum symptomatic infections and emphasize the importance of considering the children's age when studying and treating malaria infections.

Gene expression analyses reveal differences in children's response to malaria according to their age.

In Bandiagara, Mali, children experience on average two clinical malaria episodes per season. However, even in the same transmission area, the number of uncomplicated symptomatic infections, and their parasitemia, vary dramatically among children. To examine the factors contributing to these variations, we simultaneously characterized the host and parasite gene expression profiles from 136 children with symptomatic falciparum malaria and analyzed the expression of 9,205 human and 2,484 Plasmodium genes. We used gene expression deconvolution to estimate the relative proportion of immune cells and parasite stages in each sample and to adjust the differential gene expression analyses. Parasitemia explained much of the variation in both host and parasite gene expression and revealed that infections with higher parasitemia had more neutrophils and fewer T cells, suggesting parasitemia-dependent neutrophil recruitment and/or T cell extravasation to secondary lymphoid organs. The child's age was also strongly correlated with gene expression variations. Plasmodium falciparum genes associated with age suggested that older children carried more male gametocytes, while host genes associated with age indicated a stronger innate response (through TLR and NLR signaling) in younger children and stronger adaptive immunity (through TCR and BCR signaling) in older children. These analyses highlight the variability in host responses and parasite regulation during P. falciparum symptomatic infections and emphasize the importance of considering the children's age when studying and treating malaria infections.

Natural immunity to malaria preferentially targets the endothelial protein C receptor-binding regions of PfEMP1s.

Antibody responses to variant surface antigens (VSAs) produced by the malaria parasite Plasmodium falciparum may contribute to age-related natural immunity to severe malaria. One VSA family, P. falciparum erythrocyte membrane protein-1 (PfEMP1), includes a subset of proteins that binds endothelial protein C receptor (EPCR) in human hosts and potentially disrupts the regulation of inflammatory responses, which may lead to the development of severe malaria. We probed peptide microarrays containing segments spanning five PfEMP1 EPCR-binding domain variants with sera from 10 Malian adults and 10 children to determine the differences between adult and pediatric immune responses. We defined serorecognized peptides and amino acid residues as those that elicited a significantly higher antibody response than malaria-naïve controls. We aimed to identify regions consistently serorecognized among adults but not among children across PfEMP1 variants, potentially indicating regions that drive the development of immunity to severe malaria. Adult sera consistently demonstrated broader and more intense serologic responses to constitutive PfEMP1 peptides than pediatric sera, including peptides in EPCR-binding domains. Both adults and children serorecognized a significantly higher proportion of EPCR-binding peptides than peptides that do not directly participate in receptor binding, indicating a preferential development of serologic responses at functional residues. Over the course of a single malaria transmission season, pediatric serological responses increased between the start and the peak of the season, but waned as the transmission season ended. IMPORTANCE Severe malaria and death related to malaria disproportionately affect sub-Saharan children under 5 years of age, commonly manifesting as cerebral malaria and/or severe malarial anemia. In contrast, adults in malaria-endemic regions tend to experience asymptomatic or mild disease. Our findings indicate that natural immunity to malaria targets specific regions within the EPCR-binding domain, particularly peptides containing EPCR-binding residues. Epitopes containing these residues may be promising targets for vaccines or therapeutics directed against severe malaria. Our approach provides insight into the development of natural immunity to a binding target linked to severe malaria by characterizing an "adult-like" response as recognizing a proportion of epitopes within the PfEMP1 protein, particularly regions that mediate EPCR binding. This "adult-like" response likely requires multiple years of malaria exposure, as increases in pediatric serologic response over a single malaria transmission season do not appear significant.

Malian children infected with Plasmodium ovale and Plasmodium falciparum display very similar gene expression profiles.

Plasmodium parasites caused 241 million cases of malaria and over 600,000 deaths in 2020. Both P. falciparum and P. ovale are endemic to Mali and cause clinical malaria, with P. falciparum infections typically being more severe. Here, we sequenced RNA from nine pediatric blood samples collected during infections with either P. falciparum or P. ovale, and characterized the host and parasite gene expression profiles. We found that human gene expression varies more between individuals than according to the parasite species causing the infection, while parasite gene expression profiles cluster by species. Additionally, we characterized DNA polymorphisms of the parasites directly from the RNA-seq reads and found comparable levels of genetic diversity in both species, despite dramatic differences in prevalence. Our results provide unique insights into host-pathogen interactions during malaria infections and their variations according to the infecting Plasmodium species, which will be critical to develop better elimination strategies against all human Plasmodium parasites.

Protein Microarrays as a Tool to Analyze Antibody Responses to Variant Surface Antigens Expressed on the Surface of Plasmodium falciparum-Infected Erythrocytes.

Enzyme-linked immunosorbent assays (ELISAs) remain the gold standard for measuring antibodies, but are time-consuming and use significant amounts of precious sample and reagents. Protein microarrays represent an appealing alternative, particularly for studies focused on large gene families such as those encoding variant surface antigens in the malaria parasite Plasmodium falciparum. Such microarrays represent an ideal high-throughput platform to study antibody responses to hundreds of malaria parasite variant surface antigens at once, providing critical insights into the development of natural immunity to malaria. We describe the essential background and approach to run an assay using a P. falciparum microarray populated with variant surface antigens. This allows the user to define serologic profiles and identify serodominant antigens that represent promising targets for vaccine or therapeutic development.

Successful Profiling of Plasmodium falciparum var Gene Expression in Clinical Samples via a Custom Capture Array.

var genes encode Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) antigens. These highly diverse antigens are displayed on the surface of infected erythrocytes and play a critical role in immune evasion and sequestration of infected erythrocytes. Studies of var expression using non-leukocyte-depleted blood are challenging because of the predominance of host genetic material and lack of conserved var segments. Our goal was to enrich for parasite RNA, allowing de novo assembly of var genes and detection of expressed novel variants. We used two overall approaches: (i) enriching for total mRNA in the sequencing library preparations and (ii) enriching for parasite RNA with a custom capture array based on Roche's SeqCap EZ enrichment system. The capture array was designed with probes based on the whole 3D7 reference genome and an additional >4,000 full-length var gene sequences from other P. falciparum strains. We tested each method on RNA samples from Malian children with severe or uncomplicated malaria infections. All reads mapping to the human genome were removed, the remaining reads were assembled de novo into transcripts, and from these, var-like transcripts were identified and annotated. The capture array produced the longest maximum length and largest numbers of var gene transcripts in each sample, particularly in samples with low parasitemia. Identifying the most-expressed var gene sequences in whole-blood clinical samples without the need for extensive processing or generating sample-specific reference genome data is critical for understanding the role of PfEMP1s in malaria pathogenesis. IMPORTANCE Malaria parasites display antigens on the surface of infected red blood cells in the human host that facilitate attachment to blood vessels, contributing to the severity of infection. These antigens are highly variable, allowing the parasite to evade the immune system. Identifying these expressed antigens is critical to understanding the development of severe malarial disease. However, clinical samples contain limited amounts of parasite genetic material, a challenge for sequencing efforts further compounded by the extreme diversity of the parasite surface antigens. We present a method that enriches for these antigen sequences in clinical samples using a custom capture array, requiring minimal processing in the field. While our results are focused on the malaria parasite Plasmodium falciparum, this approach has broad applicability to other highly diverse antigens from other parasites and pathogens such as those that cause giardiasis and leishmaniasis.

Malian adults maintain serologic responses to virulent PfEMP1s amid seasonal patterns of fluctuation.

Plasmodium falciparum erythrocyte membrane protein-1s (PfEMP1s), diverse malaria proteins expressed on the infected erythrocyte surface, play an important role in pathogenesis, mediating adhesion to host vascular endothelium. Antibodies to particular non-CD36-binding PfEMP1s are associated with protection against severe disease. We hypothesized that given lifelong P. falciparum exposure, Malian adults would have broad PfEMP1 serorecognition and high seroreactivity levels during follow-up, particularly to non-CD36-binding PfEMP1s such as those that attach to endothelial protein C receptor (EPCR) and intercellular adhesion molecule-1 (ICAM-1). Using a protein microarray, we determined serologic responses to 166 reference PfEMP1 fragments during a dry and subsequent malaria transmission season in Malian adults. Malian adult sera had PfEMP1 serologic responses throughout the year, with decreased reactivity to a small subset of PfEMP1 fragments during the dry season and increases in reactivity to a different subset of PfEMP1 fragments during the subsequent peak malaria transmission season, especially for intracellular PfEMP1 domains. For some individuals, PfEMP1 serologic responses increased after the dry season, suggesting antigenic switching during asymptomatic infection. Adults were more likely to experience variable serorecognition of CD36-binding PfEMP1s than non-CD36-binding PfEMP1s that bind EPCR or ICAM-1, which remained serorecognized throughout the year. Sustained seroreactivity to non-CD36-binding PfEMP1s throughout adulthood amid seasonal fluctuation patterns may reflect underlying protective severe malaria immunity and merits further investigation.

Burden of perianal Staphylococcus aureus colonization in nursing home residents increases transmission to healthcare worker gowns and gloves.

OBJECTIVE: To evaluate the effect of the burden of Staphylococcus aureus colonization of nursing home residents on the risk of S. aureus transmission to healthcare worker (HCW) gowns and gloves. DESIGN: Multicenter prospective cohort study. SETTING AND PARTICIPANTS: Residents and HCWs from 13 community-based nursing homes in Maryland and Michigan. METHODS: Residents were cultured for S. aureus at the anterior nares and perianal skin. The S. aureus burden was estimated by quantitative polymerase chain reaction detecting the nuc gene. HCWs wore gowns and gloves during usual care activities; gowns and gloves were swabbed and then cultured for the presence of S. aureus. RESULTS: In total, 403 residents were enrolled; 169 were colonized with methicillin-resistant S. aureus (MRSA) or methicillin-sensitive S. aureus (MSSA) and comprised the study population; 232 were not colonized and thus were excluded from this analysis; and 2 were withdrawn prior to being swabbed. After multivariable analysis, perianal colonization with S. aureus conferred the greatest odds for transmission to HCW gowns and gloves, and the odds increased with increasing burden of colonization: adjusted odds ratio (aOR), 2.1 (95% CI, 1.3-3.5) for low-level colonization and aOR 5.2 (95% CI, 3.1-8.7) for high level colonization. CONCLUSIONS: Among nursing home patients colonized with S. aureus, the risk of transmission to HCW gowns and gloves was greater from those colonized with greater quantities of S. aureus on the perianal skin. Our findings inform future infection control practices for both MRSA and MSSA in nursing homes.

Optimization of parasite DNA enrichment approaches to generate whole genome sequencing data for Plasmodium falciparum from low parasitaemia samples.

BACKGROUND: Owing to the large amount of host DNA in clinical samples, generation of high-quality Plasmodium falciparum whole genome sequencing (WGS) data requires enrichment for parasite DNA. Enrichment is often achieved by leukocyte depletion of infected blood prior to storage. However, leukocyte depletion is difficult in low-resource settings and limits analysis to prospectively-collected samples. As a result, approaches such as selective whole genome amplification (sWGA) are being used to enrich for parasite DNA. However, sWGA has had limited success in generating reliable sequencing data from low parasitaemia samples. In this study, enzymatic digestion with MspJI prior to sWGA and whole genome sequencing was evaluated to determine whether this approach improved genome coverage compared to sWGA alone. The potential of sWGA to cause amplification bias in polyclonal infections was also examined. METHODS: DNA extracted from laboratory-created dried blood spots was treated with a modification-dependent restriction endonuclease, MspJI, and filtered via vacuum filtration. Samples were then selectively amplified using a previously reported sWGA protocol and subjected to WGS. Genome coverage statistics were compared between the optimized sWGA approach and the previously reported sWGA approach performed in parallel. Differential amplification by sWGA was assessed by comparing WGS data generated from lab-created mixtures of parasite isolates, from the same geographical region, generated with or without sWGA. RESULTS: MspJI digestion did not enrich for parasite DNA. Samples that underwent vacuum filtration (without MspJI digestion) prior to sWGA had the highest parasite DNA concentration and displayed greater genome coverage compared to MspJI + sWGA and sWGA alone, particularly for low parasitaemia samples. The optimized sWGA (filtration + sWGA) approach was successfully used to generate WGS data from 218 non-leukocyte depleted field samples from Malawi. Sequences from lab-created mixtures of parasites did not show evidence of differential amplification of parasite strains compared to directly sequenced samples. CONCLUSION: This optimized sWGA approach is a reliable method to obtain WGS data from non-leukocyte depleted, low parasitaemia samples. The absence of amplification bias in data generated from mixtures of isolates from the same geographic region suggests that this approach can be appropriately used for molecular epidemiological studies.

Strains used in whole organism Plasmodium falciparum vaccine trials differ in genome structure, sequence, and immunogenic potential.

BACKGROUND: Plasmodium falciparum (Pf) whole-organism sporozoite vaccines have been shown to provide significant protection against controlled human malaria infection (CHMI) in clinical trials. Initial CHMI studies showed significantly higher durable protection against homologous than heterologous strains, suggesting the presence of strain-specific vaccine-induced protection. However, interpretation of these results and understanding of their relevance to vaccine efficacy have been hampered by the lack of knowledge on genetic differences between vaccine and CHMI strains, and how these strains are related to parasites in malaria endemic regions. METHODS: Whole genome sequencing using long-read (Pacific Biosciences) and short-read (Illumina) sequencing platforms was conducted to generate de novo genome assemblies for the vaccine strain, NF54, and for strains used in heterologous CHMI (7G8 from Brazil, NF166.C8 from Guinea, and NF135.C10 from Cambodia). The assemblies were used to characterize sequences in each strain relative to the reference 3D7 (a clone of NF54) genome. Strains were compared to each other and to a collection of clinical isolates (sequenced as part of this study or from public repositories) from South America, sub-Saharan Africa, and Southeast Asia. RESULTS: While few variants were detected between 3D7 and NF54, we identified tens of thousands of variants between NF54 and the three heterologous strains. These variants include SNPs, indels, and small structural variants that fall in regulatory and immunologically important regions, including transcription factors (such as PfAP2-L and PfAP2-G) and pre-erythrocytic antigens that may be key for sporozoite vaccine-induced protection. Additionally, these variants directly contributed to diversity in immunologically important regions of the genomes as detected through in silico CD8+ T cell epitope predictions. Of all heterologous strains, NF135.C10 had the highest number of unique predicted epitope sequences when compared to NF54. Comparison to global clinical isolates revealed that these four strains are representative of their geographic origin despite long-term culture adaptation; of note, NF135.C10 is from an admixed population, and not part of recently formed subpopulations resistant to artemisinin-based therapies present in the Greater Mekong Sub-region. CONCLUSIONS: These results will assist in the interpretation of vaccine efficacy of whole-organism vaccines against homologous and heterologous CHMI.

Serologic responses to the PfEMP1 DBL-CIDR head structure may be a better indicator of malaria exposure than those to the DBL-α tag.

BACKGROUND: Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) antigens play a critical role in host immune evasion. Serologic responses to these antigens have been associated with protection from clinical malaria, suggesting that antibodies to PfEMP1 antigens may contribute to natural immunity. The first N-terminal constitutive domain in a PfEMP1 is the Duffy binding-like alpha (DBL-α) domain, which contains a 300 to 400 base pair region unique to each particular protein (the DBL-α "tag"). This DBL-α tag has been used as a marker of PfEMP1 diversity and serologic responses in malaria-exposed populations. In this study, using sera from a malaria-endemic region, responses to DBL-α tags were compared to responses to the corresponding entire DBL-α domain (or "parent" domain) coupled with the succeeding cysteine-rich interdomain region (CIDR). METHODS: A protein microarray populated with DBL-α tags, the parent DBL-CIDR head structures, and downstream PfEMP1 protein fragments was probed with sera from Malian children (aged 1 to 6 years) and adults from the control arms of apical membrane antigen 1 (AMA1) vaccine clinical trials before and during a malaria transmission season. Serological responses to the DBL-α tag and the DBL-CIDR head structure were measured and compared in children and adults, and throughout the season. RESULTS: Malian serologic responses to a PfEMP1's DBL-α tag region did not correlate with seasonal malaria exposure, or with responses to the parent DBL-CIDR head structure in either children or adults. Parent DBL-CIDR head structures were better indicators of malaria exposure. CONCLUSIONS: Larger PfEMP1 domains may be better indicators of malaria exposure than short, variable PfEMP1 fragments such as DBL-α tags. PfEMP1 head structures that include conserved sequences appear particularly well suited for study as serologic predictors of malaria exposure.

Antibodies to Peptides in Semiconserved Domains of RIFINs and STEVORs Correlate with Malaria Exposure.

The repetitive interspersed family (RIFIN) and the subtelomeric variable open reading frame (STEVOR) family represent two of three major Plasmodium falciparum variant surface antigen families involved in malaria pathogenesis and immune evasion and are potential targets in the development of natural immunity. Protein and peptide microarrays populated with RIFINs and STEVORs associated with severe malaria vulnerability in Malian children were probed with adult and pediatric sera to identify epitopes that reflect malaria exposure. Adult sera recognized and reacted with greater intensity to all STEVOR proteins than pediatric sera did. Serorecognition of and seroreactivity to peptides within the semiconserved domain of STEVORs increased with age and seasonal malaria exposure, while serorecognition and seroreactivity increased for the semiconserved and second hypervariable domains of RIFINs only with age. Serologic responses to RIFIN and STEVOR peptides within the semiconserved domains may play a role in natural immunity to severe malaria.IMPORTANCE Malaria, an infectious disease caused by the parasite Plasmodium falciparum, causes nearly 435,000 deaths annually worldwide. RIFINs and STEVORs are two variant surface antigen families that are involved in malaria pathogenesis and immune evasion. Recent work has shown that a lack of humoral immunity to these proteins is associated with severe malaria vulnerability in Malian children. This is the first study to have compared serologic responses of children and adults to RIFINs and STEVORs in settings of malaria endemicity and to examine such serologic responses before and after a clinical malaria episode. Using microarrays, we determined that the semiconserved domains in these two parasite variant surface antigen families harbor peptides whose seroreactivity reflects malaria exposure. A similar approach has the potential to illuminate the role of variant surface antigens in the development of natural immunity to clinical malaria. Potential vaccines for severe malaria should include consideration of peptides within the semiconserved domains of RIFINs and STEVORs.

Genetic variants at the 16p13 locus confer risk for eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is a chronic inflammatory disease of the esophagus triggered by immune hypersensitivity to food. Herein, we tested whether genetic risk factors for known, non-allergic, immune-mediated diseases, particularly those involving autoimmunity, were associated with EoE risk. We used the high-density Immunochip platform, encoding 200,000 genetic variants for major auto-immune disease. Accordingly, 1214 subjects with EoE of European ancestry and 3734 population controls were genotyped and assessed using data directly generated or imputed from the previously published GWAS. We found lack of association of EoE with the genetic variants in the major histocompatibility complex (MHC) class I, II, and III genes and nearly all other loci using a highly powered study design with dense genotyping throughout the locus. Importantly, we identified an EoE risk locus at 16p13 with genome-wide significance (Pcombined=2.05 × 10-9, odds ratio = 0.76-0.81). This region is known to encode for the genes CLEC16A, DEXI, and CIITI, which are expressed in immune cells and esophageal epithelial cells. Suggestive EoE risk were also seen 5q23 (intergenic) and 7p15 (JAZF1). Overall, we have identified an additional EoE risk locus at 16p13 and highlight a shared and unique genetic etiology of EoE with a spectrum of immune-associated diseases.

Whole-exome sequencing uncovers oxidoreductases DHTKD1 and OGDHL as linkers between mitochondrial dysfunction and eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is an allergic inflammatory esophageal disorder with a complex underlying genetic etiology often associated with other comorbidities. Using whole-exome sequencing (WES) of 63 patients with EoE and 60 unaffected family members and family-based trio analysis, we sought to uncover rare coding variants. WES analysis identified 5 rare, damaging variants in dehydrogenase E1 and transketolase domain-containing 1 (DHTKD1). Rare variant burden analysis revealed an overabundance of putative, potentially damaging DHTKD1 mutations in EoE (P = 0.01). Interestingly, we also identified 7 variants in the DHTKD1 homolog oxoglutarate dehydrogenase-like (OGDHL). Using shRNA-transduced esophageal epithelial cells and/or patient fibroblasts, we further showed that disruption of normal DHTKD1 or OGDHL expression blunts mitochondrial function. Finally, we demonstrated that the loss of DHTKD1 expression increased ROS production and induced the expression of viperin, a gene previously shown to be involved in production of Th2 cytokines in T cells. Viperin had increased expression in esophageal biopsies of EoE patients compared with control individuals and was upregulated by IL-13 in esophageal epithelial cells. These data identify a series of rare genetic variants implicating DHTKD1 and OGDHL in the genetic etiology of EoE and underscore a potential pathogenic role for mitochondrial dysfunction in EoE.

Paired Ig-like Receptor B Inhibits IL-13-Driven Eosinophil Accumulation and Activation in the Esophagus.

Eosinophilic esophagitis (EoE) is a Th2 cytokine-associated disease characterized by eosinophil infiltration, epithelial cell hyperplasia, and tissue remodeling. Recent studies highlighted a major contribution for IL-13 in EoE pathogenesis. Paired Ig-like receptor B is a cell surface immune-inhibitory receptor that is expressed by eosinophils and postulated to regulate eosinophil development and migration. We report that Pirb is upregulated in the esophagus after inducible overexpression of IL-13 (CC10-Il13(Tg) mice) and is overexpressed by esophageal eosinophils. CC10-Il13(Tg)/Pirb(-/-) mice displayed increased esophageal eosinophilia and EoE pathology, including epithelial cell thickening, fibrosis, and angiogenesis, compared with CC10-Il13(Tg)/Pirb(+/+) mice. Transcriptome analysis of primary Pirb(+/+) and Pirb(-/-) esophageal eosinophils revealed increased expression of transcripts associated with promoting tissue remodeling in Pirb(-/-) eosinophils, including profibrotic genes, genes promoting epithelial-to-mesenchymal transition, and genes associated with epithelial growth. These data identify paired Ig-like receptor B as a molecular checkpoint in IL-13-induced eosinophil accumulation and activation, which may serve as a novel target for future therapy in EoE.