Effective presentation of antigens by human leukocyte antigen (HLA) class I molecules to CD8+ T cells is required for viral elimination and generation of long-term immunological memory. In this study, we applied a single-cell, multiomic technology to generate a unified ex vivo characterization of the CD8+ T cell response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) across four major HLA class I alleles. We found that HLA genotype conditions key features of epitope specificity, TCRα/β sequence diversity, and the utilization of pre-existing SARS-CoV-2-reactive memory T cell pools. Single-cell transcriptomics revealed functionally diverse T cell phenotypes of SARS-CoV-2-reactive T cells, associated with both disease stage and epitope specificity. Our results show that HLA variations notably influence the CD8+ T cell repertoire shape and utilization of immune recall upon SARS-CoV-2 infection.

Effective presentation of antigens by HLA class I molecules to CD8 + T cells is required for viral elimination and generation of long-term immunological memory. In this study, we applied a single-cell, multi-omic technology to generate the first unified ex vivo characterization of the CD8 + T cell response to SARS-CoV-2 across 4 major HLA class I alleles. We found that HLA genotype conditions key features of epitope specificity, TCR α/β sequence diversity, and the utilization of pre-existing SARS-CoV-2 reactive memory T cell pools. Single-cell transcriptomics revealed functionally diverse T cell phenotypes of SARS-CoV-2-reactive T cells, associated with both disease stage and epitope specificity. Our results show that HLA variations influence pre-existing immunity to SARS-CoV-2 and shape the immune repertoire upon subsequent viral exposure. One-Sentence Summary We perform a unified, multi-omic characterization of the CD8 + T cell response to SARS-CoV-2, revealing pre-existing immunity conditioned by HLA genotype.

T cell-mediated immunity plays an important role in controlling SARS-CoV-2 infection, but the repertoire of naturally processed and presented viral epitopes on class I human leukocyte antigen (HLA-I) remains uncharacterized. Here, we report the first HLA-I immunopeptidome of SARS-CoV-2 in two cell lines at different times post infection using mass spectrometry. We found HLA-I peptides derived not only from canonical open reading frames (ORFs) but also from internal out-of-frame ORFs in spike and nucleocapsid not captured by current vaccines. Some peptides from out-of-frame ORFs elicited T cell responses in a humanized mouse model and individuals with COVID-19 that exceeded responses to canonical peptides, including some of the strongest epitopes reported to date. Whole-proteome analysis of infected cells revealed that early expressed viral proteins contribute more to HLA-I presentation and immunogenicity. These biological insights, as well as the discovery of out-of-frame ORF epitopes, will facilitate selection of peptides for immune monitoring and vaccine development.

The malaria mosquito, Anopheles stephensi, and other mosquitoes modulate their biology to match the time-of-day. In the present work, we used a non-hypothesis driven approach (untargeted proteomics) to identify proteins in mosquito tissue, and then quantified the relative abundance of the identified proteins from An. stephensi bodies. Using these quantified protein levels, we then analyzed the data for proteins that were only detectable at certain times-of-the day, highlighting the need to consider time-of-day in experimental design. Further, we extended our time-of-day analysis to look for proteins which cycle in a rhythmic 24-hour ("circadian") manner, identifying 31 rhythmic proteins. Finally, to maximize the utility of our data, we performed a proteogenomic analysis to improve the genome annotation of An. stephensi. We compare peptides that were detected using mass spectrometry but are 'missing' from the An. stephensi predicted proteome, to reference proteomes from 38 other primarily human disease vector species. We found 239 such peptide matches and reveal that genome annotation can be improved using proteogenomic analysis from taxonomically diverse reference proteomes. Examination of 'missing' peptides revealed reading frame errors, errors in gene-calling, overlapping gene models, and suspected gaps in the genome assembly.

Aim: The aim of this study was to determine the immunologic effects and safety of oral anti-CD3 in patients with ulcerative colitis (UC). Methods: An open-label pilot study of orally delivered anti-CD3 was performed in patients with moderate-to-severe UC. The primary end points were changes in immunologic parameters and evaluation for safety. Results: Six subjects received oral OKT3. Biologic effects of oral anti-CD3 included significantly increased proliferation in response to anti-CD3 and anti-inflammatory gene expression profile in peripheral blood mononuclear cells. No serious treatment-related adverse events occurred. Conclusion: Orally delivered anti-CD3 resulted in immunologic changes in patients with UC.

Structural brain changes in gray and white matter have been found in adults with Crohn’s disease (CD). In adults, increased disease activity has also been associated with decreased cognitive performance. We have shown similar effects for pediatric CD, particularly for gray matter and associated cognitive and emotional functions. Underlying inflammation and steroid therapy have emerged as contributors, however, associations between brain structure and immune pathways have not been studied. We here aim to replicate earlier findings of brain involvement in pediatric CD, including study of white matter structures, and to elucidate the impact of immune response on the central nervous system.

Vector-borne diseases are responsible for > 1 million deaths every year but genomic resources for most species responsible for their transmission are limited. This is true for neglected diseases such as sleeping sickness (Human African Trypanosomiasis), a disease caused by Trypanosoma parasites vectored by several species of tseste flies within the genus Glossina We describe an integrative approach that identifies statistical associations between trypanosome infection status of Glossina fuscipes fuscipes (Gff) flies from Uganda, for which functional studies are complicated because the species cannot be easily maintained in laboratory colonies, and ∼73,000 polymorphic sites distributed across the genome. Then, we identify candidate genes involved in Gff trypanosome susceptibility by taking advantage of genomic resources from a closely related species, G. morsitans morsitans (Gmm). We compiled a comprehensive transcript library from 72 published and unpublished RNAseq experiments of trypanosome-infected and uninfected Gmm flies, and improved the current Gmm transcriptome assembly. This new assembly was then used to enhance the functional annotations on the Gff genome. As a consequence, we identified 56 candidate genes in the vicinity of the 18 regions associated with Trypanosoma infection status in Gff Twenty-nine of these genes were differentially expressed (DE) among parasite-infected and uninfected Gmm, suggesting that their orthologs in Gff may correlate with disease transmission. These genes were involved in DNA regulation, neurophysiological functions, and immune responses. We highlight the power of integrating population and functional genomics from related species to enhance our understanding of the genetic basis of physiological traits, particularly in nonmodel organisms.

Background: Transcriptional profiling has been performed on biopsies from ulcerative colitis patients. Limitations in prior studies include the variability introduced by inflammation, anatomic site of biopsy, extent of disease, and medications. We sought to more globally understand the variability of gene expression from patients with ulcerative colitis to advance our understanding of its pathogenesis and to guide clinical study design. Methods: We performed transcriptional profiling on 13 subjects, including pediatric and adult patients from 2 hospital sites. For each patient, we collected 6 biopsies from macroscopically inflamed tissue and 4 biopsies from macroscopically healthy-appearing tissue. Isolated RNA was used for microarray gene expression analysis utilizing Affymetrix Human Primeview microarrays. Ingenuity pathway analysis was used to assess over-representation of gene ontology and biological pathways. RNAseq was also performed, and differential analysis was assessed to compare affected vs unaffected samples. Finally, we modeled the minimum number of biopsies required to reliably detect gene expression across different subject numbers. Results: Transcriptional profiles co-clustered independently of the hospital collection site, patient age, sex, and colonic location, which parallels prior gene expression findings. A small set of genes not previously described was identified. Our modeling analysis reveals the number of biopsies and patients per cohort to yield reliable results in clinical studies. Conclusions: Key findings include concordance, including some expansion, of previously published gene expression studies and similarity among different age groups. We also established a reliable statistical model for biopsy collection for future clinical studies.

Tsetse flies (genus Glossina) are the only vector for the parasitic trypanosomes responsible for sleeping sickness and nagana across sub-Saharan Africa. In Uganda, the tsetse fly Glossina fuscipes fuscipes is responsible for transmission of the parasite in 90% of sleeping sickness cases, and co-occurrence of both forms of human-infective trypanosomes makes vector control a priority. We use population genetic data from 38 samples from northern Uganda in a novel methodological pipeline that integrates genetic data, remotely sensed environmental data, and hundreds of field-survey observations. This methodological pipeline identifies isolated habitat by first identifying environmental parameters correlated with genetic differentiation, second, predicting spatial connectivity using field-survey observations and the most predictive environmental parameter(s), and third, overlaying the connectivity surface onto a habitat suitability map. Results from this pipeline indicated that net photosynthesis was the strongest predictor of genetic differentiation in G. f. fuscipes in northern Uganda. The resulting connectivity surface identified a large area of well-connected habitat in northwestern Uganda, and twenty-four isolated patches on the northeastern margin of the G. f. fuscipes distribution. We tested this novel methodological pipeline by completing an ad hoc sample and genetic screen of G. f. fuscipes samples from a model-predicted isolated patch, and evaluated whether the ad hoc sample was in fact as genetically isolated as predicted. Results indicated that genetic isolation of the ad hoc sample was as genetically isolated as predicted, with differentiation well above estimates made in samples from within well-connected habitat separated by similar geographic distances. This work has important practical implications for the control of tsetse and other disease vectors, because it provides a way to identify isolated populations where it will be safer and easier to implement vector control and that should be prioritized as study sites during the development and improvement of vector control methods.

Background: Diagnostic evaluation of eosinophilic esophagitis (EoE) remains difficult, particularly the assessment of the patient's allergic status. Objective: This study sought to establish an automated medical algorithm to assist in the evaluation of EoE. Methods: Machine learning techniques were used to establish a diagnostic probability score for EoE, p(EoE), based on esophageal mRNA transcript patterns from biopsies of patients with EoE, gastroesophageal reflux disease and controls. Dimensionality reduction in the training set established weighted factors, which were confirmed by immunohistochemistry. Following weighted factor analysis, p(EoE) was determined by random forest classification. Accuracy was tested in an external test set, and predictive power was assessed with equivocal patients. Esophageal IgE production was quantified with epsilon germ line (IGHE) transcripts and correlated with serum IgE and the T(H)2-type mRNA profile to establish an IGHE score for tissue allergy. Results: In the primary analysis, a 3-class statistical model generated a p(EoE) score based on common characteristics of the inflammatory EoE profile. A p(EoE) >= 25 successfully identified EoE with high accuracy (sensitivity: 90.9%, specificity: 93.2%, area under the curve: 0.985) and improved diagnosis of equivocal cases by 84.6%. The p(EoE) changed in response to therapy. A secondary analysis loop in EoE patients defined an IGHE score of >= 37.5 for a patient subpopulation with increased esophageal allergic inflammation. Conclusions: The development of intelligent data analysis from a machine learning perspective provides exciting opportunities to improve diagnostic precision and improve patient care in EoE. The p(EoE) and the IGHE score are steps toward the development of decision trees to define EoE subpopulations and, consequently, will facilitate individualized therapy.

Background: Structural brain changes in gray and white matter have been previously found in adults with Crohn's disease (CD). We have recently shown similar effects for pediatric CD in 2 separate studies (Mrakotsky et al., 2012, 2013, 2015), particularly for cortical and subcortical brain regions important for cognition, memory and emotion. Our prior data also revealed serum markers of inflammation and steroid therapy to be negatively associated with cortical thickness, subcortical volume, cognitive and school function, however, associations between brain structure and more detailed inflammatory profiles have not been studied. This CCFA SRA aims to replicate prior studies on brain involvement in pediatric CD, and provide initial information on potential contributing immune pathways. Methods: Structural magnetic resonance imaging (MRI) was performed in 54 children age 10 to 14 years (CD patients: n = 37, age-and sex-matched healthy controls: n = 17) using a 3 Tesla Siemens Trio MRI system. Anatomical whole brain T1-weighted images were acquired at 1 mm3 voxel size. All subjects also underwent neuropsychological assessments (tests of attention, memory, IQ, emotional functions), disease severity ratings (PCDAI), and phlebotomy to obtain inflammatory markers. CD patients were scanned and tested during an active disease flare or in clinical remission. Cortical thickness and volume was analyzed using Freesurfer. Gene expression analysis of immunology-related genes was performed via multiplex array (Nanostring nCounter Human Immunology v2, nSolver v3 Advanced Analysis). Analyses were adjusted for age and sex, and included correlation, regression and group comparisons. Multiple comparisons were corrected with false discovery rate. Groups were comparable on demographics. Results: We confirmed our initial finding of reduced gray matter (cortical thickness in posterior brain regions (bilateral inferior parietal, superior parietal, supramarginal [all pFDR < 0.05]), subcortical volume (bilateral hippocampi and thalami [all pFDR < 0.05]), and poorer cognitive (P < 0.05) and memory function (P < 0.01) in CD patients compared to healthy children. We further found differences in mRNA gene expression in a range of immune-related genes in CD patients compared to controls (i.e., TIGIT, KLRG1, KLRK1 [downregulated], CD44, ITGA2B, NFKB [upregulated] [pFDR = 0.02–0.10]). Regression including both differentially expressed and a priori defined immune genes revealed some associations between cytokine (i.e., IL1R1, IL6R, IL18RAP), NK cells, NFKB1 gene expression with gray matter (supramarginal, inferoparietal thickness, thalamic volume), after controlling for steroids. Steroids remained a robust contributor to brain structural outcomes. Conclusions: Complete data from this CCFA SRA confirm our prior findings of reduced gray matter in pediatric CD, particularly in posterior cortical and subcortical limbic regions critical for higher cognition, emotional and pain processing. We found differentially expressed immune genes in CD patients that may play a role in inflammatory response, cell migration, and apoptosis and have also been implicated in brain development and injury. The preliminary association of immune gene expression with gray matter brain structures provides first insights in the molecular basis of brain involvement in IBD, and may further lend support for novel treatment targets.

Uganda is the only country where the chronic and acute forms of human African Trypanosomiasis (HAT) or sleeping sickness both occur and are separated by < 100 km in areas north of Lake Kyoga. In Uganda, Glossina fuscipes fuscipes is the main vector of the Trypanosoma parasites responsible for these diseases as well for the animal African Trypanosomiasis (AAT), or Nagana. We used highly polymorphic microsatellite loci and a mitochondrial DNA (mtDNA) marker to provide fine scale spatial resolution of genetic structure of G. f. fuscipes from 42 sampling sites from the northern region of Uganda where a merger of the two disease belts is feared. Based on microsatellite analyses, we found that G. f. fuscipes in northern Uganda are structured into three distinct genetic clusters with varying degrees of interconnectivity among them. Based on genetic assignment and spatial location, we grouped the sampling sites into four genetic units corresponding to northwestern Uganda in the Albert Nile drainage, northeastern Uganda in the Lake Kyoga drainage, western Uganda in the Victoria Nile drainage, and a transition zone between the two northern genetic clusters characterized by high level of genetic admixture. An analysis using HYBRIDLAB supported a hybrid swarm model as most consistent with tsetse genotypes in these admixed samples. Results of mtDNA analyses revealed the presence of 30 haplotypes representing three main haplogroups, whose location broadly overlaps with the microsatellite defined clusters. Migration analyses based on microsatellites point to moderate migration among the northern units located in the Albert Nile, Achwa River, Okole River, and Lake Kyoga drainages, but not between the northern units and the Victoria Nile drainage in the west. Effective population size estimates were variable with low to moderate sizes in most populations and with evidence of recent population bottlenecks, especially in the northeast unit of the Lake Kyoga drainage. Our microsatellite and mtDNA based analyses indicate that G. f. fuscipes movement along the Achwa and Okole rivers may facilitate northwest expansion of the Rhodesiense disease belt in Uganda. We identified tsetse migration corridors and recommend a rolling carpet approach from south of Lake Kyoga northward to minimize disease dispersal and prevent vector re-colonization. Additionally, our findings highlight the need for continuing tsetse monitoring efforts during and after control.

The tsetse fly Glossina fuscipes fuscipes (Gff) is the insect vector of the two forms of Human African Trypanosomiasis (HAT) that exist in Uganda. Understanding Gff population dynamics, and the underlying genetics of epidemiologically relevant phenotypes is key to reducing disease transmission. Using ddRAD sequence technology, complemented with whole-genome sequencing, we developed a panel of ∼73,000 single-nucleotide polymorphisms (SNPs) distributed across the Gff genome that can be used for population genomics and to perform genome-wide-association studies. We used these markers to estimate genomic patterns of linkage disequilibrium (LD) in Gff, and used the information, in combination with outlier-locus detection tests, to identify candidate regions of the genome under selection. LD in individual populations decays to half of its maximum value (r2max/2) between 1359 and 2429 bp. The overall LD estimated for the species reaches r2max/2 at 708 bp, an order of magnitude slower than in Drosophila. Using 53 infected (Trypanosoma spp.) and uninfected flies from four genetically distinct Ugandan populations adapted to different environmental conditions, we were able to identify SNPs associated with the infection status of the fly and local environmental adaptation. The extent of LD in Gff likely facilitated the detection of loci under selection, despite the small sample size. Furthermore, it is probable that LD in the regions identified is much higher than the average genomic LD due to strong selection. Our results show that even modest sample sizes can reveal significant genetic associations in this species, which has implications for future studies given the difficulties of collecting field specimens with contrasting phenotypes for association analysis.

Background: The extensive use of pyrethroids for control of malaria vectors, driven by their cost, efficacy and safety, has led to widespread resistance. To favor their sustainable use, the World Health Organization (WHO) formulated an insecticide resistance management plan, which includes the identification of the mechanisms of resistance and resistance surveillance. Recognized physiological mechanisms of resistance include target site mutations in the para voltage-gated sodium channel, metabolic detoxification and penetration resistance. Such understanding of resistance mechanisms has allowed the development of resistance monitoring tools, including genotyping of the kdr mutation L1014F/S in the para gene.Methods: The sequence-based technique RNA-seq was applied to study changes in the transcriptome of deltamethrin-resistant and -susceptible Anopheles gambiae mosquitoes from the Western Province of Kenya. The resulting gene expression profiles were compared to data in the most recent literature to derive a list of candidate resistance genes. RNA-seq data were analyzed also to identify sequence polymorphisms linked to resistance.Results: A total of five candidate-resistance genes (AGAP04177, AGAP004572, AGAP008840, AGAP007530 and AGAP013036) were identified with altered expression between resistant and susceptible mosquitoes from West and East Africa. A change from G to C at position 36043997 of chromosome 3R resulting in A101G of the sulfotransferase gene AGAP009551 was significantly associated with the resistance phenotype (odds ratio: 5.10). The kdr L1014S mutation was detected at similar frequencies in both phenotypically resistant and susceptible mosquitoes, suggesting it is no longer fully predictive of the resistant phenotype.Conclusions: Overall, these results support the conclusion that resistance to pyrethroids is a complex and evolving phenotype, dependent on multiple gene functions including, but not limited to, metabolic detoxification. Functional convergence among metabolic detoxification genes may exist, with the role of each gene being modulated by the life history and selection pressure on mosquito populations. As a consequence, biochemical assays that quantify overall enzyme activity may be a more suitable method for predicting metabolic resistance than gene-based assays.

Numerous studies have examined changes in transcript levels after Anopheles gambiae takes a blood meal. Marinotti et al. (2006) used microarrays and reported massive changes in transcript levels 3 h after feeding (BF3h) compared to non-blood fed (NBF). We were intrigued by the number of transcripts for structural cuticular proteins (CPs) that showed such major differences in levels and employed paired-end (50 bp) RNA-seq technology to compare whole body transcriptomes from 5-day-old females NBF and BF3h. We detected transcripts for the majority of CPs (164/243) but levels of only 12 were significantly altered by the blood meal. While relative transcript levels of NBF females were somewhat similar to the microarray data, there were major differences in BF3h animals, resulting in levels of many transcripts, both for CPs and other genes changing in the opposite direction. We compared our data also to other studies done with both microarrays and RNA-seq. Findings were consistent that a small number of CP genes have transcripts that persist even in 5-day-old adults. Some of these transcripts showed diurnal rhythms (Rund et al., 2013; Rinker et al., 2013). In situ hybridization revealed that transcripts for several of these CP genes were found exclusively or predominantly in the eye. Transcripts other than for CPs that changed in response to blood-feeding were predominantly expressed in midgut and Malpighian tubules. Even in these tissues, genes responsible for proteins with similar functions, such as immunity or digestion, responded differently, with transcript levels for some rising and others falling. These data demonstrate that genes coding for some CPs are dynamic in expression even in adults and that the response to a blood meal is rapid and precisely orchestrated.

Background Dengue is the most prevalent arboviral disease world-wide and its primary vector is the mosquito Aedes aegypti . The current lack of commercially-available vaccines makes control of vector populations the only effective strategy to prevent dengue transmission. Aedes aegypti geographic populations exhibit great variability in insecticide resistance and susceptibility to dengue infection. The characterization of single nucleotide polymorphisms (SNPs) as molecular markers to study quantitatively this variation is needed greatly because this species has a low abundance of microsatellite markers and limited known restriction fragments length polymorphisms (RFLPs) and single-strand conformation polymorphism (SSCP) markers. Results We used RNA-seq to characterize SNPs in three Ae. aegypti strains, including the Liverpool (LVP) strain, from which the current genome annotation is derived. We identified 131,764 unique genome locations with at least one alternative nucleotide to what is reported in the reference annotation. These comprised changes in both open-reading frames (ORFs) and untranslated regions (UTRs) of transcripts . An in depth-look at sequence variation in immunity genes revealed that those associated with autophagy, MD2-like receptors and Peptidoglycan Recognition Proteins had more sequence variation in their 3’UTRs than mutations associated with non-synonymous changes. This supports the conclusion that these genes had maintained their functional specificity while being adapted to different regulatory domains. In contrast, a number of peroxidases, serpins and Clip-domain serine proteases exhibited conservation of putative UTR regulatory sequences while displaying diversification of the ORFs. Transcriptome evidence also was found for ~2500 novel transcriptional units (NTUs) not annotated in the reference genome. Conclusions The transcriptome-wide assessment of within and inter-strain polymorphisms in Ae. aegypti adds considerably to the number of molecular markers available for genetic studies in this mosquito. Additionally, data supporting NTU discovery emphasizes the need for continuous amendments of the reference genome annotation.

This figure demonstrates the complexities of a typical RNA-seq experiment as analyzed with the "Tuxedo Protocol". It models a fairly simple two condition experiment with each condition having three replications. It is "abbreviated" in that it only displays the output files that will be used in the next step for all steps except the 'cuffdiff' step. The complete output of 'cuffdiff' is included to demonstrate the final challenge of integrating the data which exists in multiple cross-referenced files.