T Cells in Nonlymphoid Tissues Give Rise to Lymph-Node-Resident Memory T Cells.

Immunosurveillance of secondary lymphoid organs (SLO) is performed by central memory T cells that recirculate through blood. Resident memory T (Trm) cells remain parked in nonlymphoid tissues and often stably express CD69. We recently identified Trm cells within SLO, but the origin and phenotype of these cells remains unclear. Using parabiosis of "dirty" mice, we found that CD69 expression is insufficient to infer stable residence of SLO Trm cells. Restimulation of nonlymphoid memory CD8+ T cells within the skin or mucosa resulted in a substantial increase in bona fide Trm cells specifically within draining lymph nodes. SLO Trm cells derived from emigrants from nonlymphoid tissues and shared some transcriptional and phenotypic signatures associated with nonlymphoid Trm cells. These data indicate that nonlymphoid cells can give rise to SLO Trm cells and suggest vaccination strategies by which memory CD8+ T cell immunosurveillance can be regionalized to specific lymph nodes.

Allergen-induced DNA release by the airway epithelium amplifies type 2 immunity.

BACKGROUND: Alternaria alternata and house dust mite exposure evokes IL-33 secretion from the airway epithelium, which functions as an alarmin to stimulate type 2 immunity. Extracellular DNA (eDNA) is also an alarmin that intensifies inflammation in cystic fibrosis, chronic obstructive pulmonary disease, and asthma. OBJECTIVE: We investigated the mechanisms underlying allergen-evoked DNA mobilization and release from the airway epithelium and determined the role of eDNA in type 2 immunity. METHODS: Human bronchial epithelial (hBE) cells were used to characterize allergen-induced DNA mobilization and extracellular release using comet assays to measure DNA fragmentation, Qubit double-stranded DNA assays to measure DNA release, and DNA sequencing to determine eDNA composition. Mice were used to investigate the role of eDNA in type 2 immunity. RESULTS: Alternaria extract rapidly induces mitochondrial and nuclear DNA release from human bronchial epithelial cells, whereas house dust mite extract induces mitochondrial DNA release. Caspase-3 is responsible for nuclear DNA fragmentation and becomes activated after cleavage by furin. Analysis of secreted nuclear DNA showed disproportionally higher amounts of promotor and exon sequences and lower intron and intergenic regions compared to predictions of random DNA fragmentation. In mice, Alternaria-induced type 2 immune responses were blocked by pretreatment with a DNA scavenger. In caspase-3-deficient mice, Alternaria-induced DNA release was suppressed. Furthermore, intranasal administration of mouse genomic DNA with Alternaria amplified secretion of IL-5 and IL-13 into bronchoalveolar lavage fluid while DNA alone had no effect. CONCLUSION: These findings highlight a novel, allergen-induced mechanism of rapid DNA release that amplifies type 2 immunity in airways.

A De Novo Whole Genome Assembly and Annotation of Parelaphostrongylus tenuis.

Parelaphostrongylus tenuis causes ungulate morbidity and mortality in eastern and central North America, but no reference genome sequence exists to facilitate research. Here, we present a P. tenuis genome assembly and annotation, generated with PacBio and Illumina technologies. The assembly is 491 Mbp, with 7285 scaffolds and 185 kb N50.

Whole-genome assembly and annotation of northern wild rice, Zizania palustris L., supports a whole-genome duplication in the Zizania genus.

Zizania palustris L. (northern wild rice, NWR) is an aquatic grass native to North America that is notable for its nutritious grain. This is an important species with ecological, cultural and agricultural significance, specifically in the Great Lakes region of the USA. Using flow cytometry, we first estimated the NWR genome size to be 1.8 Gb. Using long- and short-range sequencing, Hi-C scaffolding and RNA-seq data from eight tissues, we generated an annotated whole-genome de novo assembly of NWR. The assembly was 1.29 Gb in length, highly repetitive (approx. 76.0%) and contained 46 421 putative protein-coding genes. The expansion of retrotransposons within the genome and a whole-genome duplication (WGD) after the Zizania-Oryza speciation event have both led to an increase in the genome size of NWR in comparison with Oryza sativa L. and Zizania latifolia. Both events depict a genome rapidly undergoing change over a short evolutionary time. Comparative analyses revealed the conservation of large syntenic blocks between NWR and O. sativa, which were used to identify putative seed-shattering genes. Estimates of divergence times revealed that the Zizania genus diverged from Oryza approximately 26-30 million years ago (26-30 MYA), whereas NWR and Z. latifolia diverged from one another approximately 6-8 MYA. Comparative genomics confirmed evidence of a WGD in the Zizania genus and provided support that the event occurred prior to the NWR-Z. latifolia speciation event. This genome assembly and annotation provides a valuable resource for comparative genomics in the Oryzeae tribe and provides an important resource for future conservation and breeding efforts of NWR.

Cell type- and replication stage-specific influenza virus responses in vivo.

Influenza A viruses (IAVs) remain a significant global health burden. Activation of the innate immune response is important for controlling early virus replication and spread. It is unclear how early IAV replication events contribute to immune detection. Additionally, while many cell types in the lung can be infected, it is not known if all cell types contribute equally to establish the antiviral state in the host. Here, we use single-cycle influenza A viruses (scIAVs) to characterize the early immune response to IAV in vitro and in vivo. We found that the magnitude of virus replication contributes to antiviral gene expression within infected cells prior to the induction of a global response. We also developed a scIAV that is only capable of undergoing primary transcription, the earliest stage of virus replication. Using this tool, we uncovered replication stage-specific responses in vitro and in vivo. Using several innate immune receptor knockout cell lines, we identify RIG-I as the predominant antiviral detector of primary virus transcription and amplified replication in vitro. Through a Cre-inducible reporter mouse, we used scIAVs expressing Cre-recombinase to characterize cell type-specific responses in vivo. Individual cell types upregulate unique sets of antiviral genes in response to both primary virus transcription and amplified replication. We also identified antiviral genes that are only upregulated in response to direct infection. Altogether, these data offer insight into the early mechanisms of antiviral gene activation during influenza A infection.

SQANTI: extensive characterization of long-read transcript sequences for quality control in full-length transcriptome identification and quantification.

High-throughput sequencing of full-length transcripts using long reads has paved the way for the discovery of thousands of novel transcripts, even in well-annotated mammalian species. The advances in sequencing technology have created a need for studies and tools that can characterize these novel variants. Here, we present SQANTI, an automated pipeline for the classification of long-read transcripts that can assess the quality of data and the preprocessing pipeline using 47 unique descriptors. We apply SQANTI to a neuronal mouse transcriptome using Pacific Biosciences (PacBio) long reads and illustrate how the tool is effective in characterizing and describing the composition of the full-length transcriptome. We perform extensive evaluation of ToFU PacBio transcripts by PCR to reveal that an important number of the novel transcripts are technical artifacts of the sequencing approach and that SQANTI quality descriptors can be used to engineer a filtering strategy to remove them. Most novel transcripts in this curated transcriptome are novel combinations of existing splice sites, resulting more frequently in novel ORFs than novel UTRs, and are enriched in both general metabolic and neural-specific functions. We show that these new transcripts have a major impact in the correct quantification of transcript levels by state-of-the-art short-read-based quantification algorithms. By comparing our iso-transcriptome with public proteomics databases, we find that alternative isoforms are elusive to proteogenomics detection. SQANTI allows the user to maximize the analytical outcome of long-read technologies by providing the tools to deliver quality-evaluated and curated full-length transcriptomes.

Identification of gut microbiota and microbial metabolites regulated by an antimicrobial peptide lipocalin 2 in high fat diet-induced obesity.

Lipocalin 2 (Lcn2), as an antimicrobial peptide is expressed in intestine, and the upregulation of intestinal Lcn2 has been linked to inflammatory bowel disease. However, the role of Lcn2 in shaping gut microbiota during diet-induced obesity (DIO) remains unknown. We found that short-term high fat diet (HFD) feeding strongly stimulates intestinal Lcn2 expression and secretion into the gut lumen. As the HFD feeding prolongs, fecal Lcn2 levels turn to decrease. Lcn2 deficiency accelerates the development of HFD-induced intestinal inflammation and microbiota dysbiosis. Moreover, Lcn2 deficiency leads to the remodeling of microbiota-derived metabolome, including decreased production of short-chain fatty acids (SCFAs) and SCFA-producing microbes. Most importantly, we have identified Lcn2-targeted bacteria and microbiota-derived metabolites that potentially play roles in DIO and metabolic dysregulation. Correlation analyses suggest that Lcn2-targeted Dubosiella and Angelakisella have a novel role in regulating SCFAs production and obesity. Our results provide a novel mechanism involving Lcn2 as an antimicrobial host factor in the control of gut microbiota symbiosis during DIO.

The Impact of TCR Signal Strength on Resident Memory T Cell Formation during Influenza Virus Infection.

Resident memory T cells (TRM) in the lung are vital for heterologous protection against influenza A virus (IAV). Environmental factors are necessary to establish lung TRM; however, the role of T cell-intrinsic factors like TCR signal strength have not been elucidated. In this study, we investigated the impact of TCR signal strength on the generation and maintenance of lung TRM after IAV infection. We inserted high- and low-affinity OT-I epitopes into IAV and infected mice after transfer of OT-I T cells. We uncovered a bias in TRM formation in the lung elicited by lower affinity TCR stimulation. TCR affinity did not impact the overall phenotype or long-term maintenance of lung TRM Overall, these findings demonstrate that TRM formation is negatively correlated with increased TCR signal strength. Lower affinity cells may have an advantage in forming TRM to ensure diversity in the Ag-specific repertoire in tissues.

Distinct antiviral signatures revealed by the magnitude and round of influenza virus replication in vivo.

Influenza virus has a broad cellular tropism in the respiratory tract. Infected epithelial cells sense the infection and initiate an antiviral response. To define the antiviral response at the earliest stages of infection we used a series of single-cycle reporter viruses. These viral probes demonstrated cells in vivo harbor a range in magnitude of virus replication. Transcriptional profiling of cells supporting different levels of replication revealed tiers of IFN-stimulated gene expression. Uninfected cells and cells with blunted replication expressed a distinct and potentially protective antiviral signature, while cells with high replication expressed a unique reserve set of antiviral genes. Finally, we used these single-cycle reporter viruses to determine the antiviral landscape during virus spread, which unveiled disparate protection of epithelial cell subsets mediated by IFN in vivo. Together these results highlight the complexity of virus-host interactions within the infected lung and suggest that magnitude and round of replication tune the antiviral response.

Activated entomopathogenic nematode infective juveniles release lethal venom proteins.

Entomopathogenic nematodes (EPNs) are unique parasites due to their symbiosis with entomopathogenic bacteria and their ability to kill insect hosts quickly after infection. It is widely believed that EPNs rely on their bacterial partners for killing hosts. Here we disproved this theory by demonstrating that the in vitro activated infective juveniles (IJs) of Steinernema carpocapsae (a well-studied EPN species) release venom proteins that are lethal to several insects including Drosophila melanogaster. We confirmed that the in vitro activation is a good approximation of the in vivo process by comparing the transcriptomes of individual in vitro and in vivo activated IJs. We further analyzed the transcriptomes of non-activated and activated IJs and revealed a dramatic shift in gene expression during IJ activation. We also analyzed the venom proteome using mass spectrometry. Among the 472 venom proteins, proteases and protease inhibitors are especially abundant, and toxin-related proteins such as Shk domain-containing proteins and fatty acid- and retinol-binding proteins are also detected, which are potential candidates for suppressing the host immune system. Many of the venom proteins have conserved orthologs in vertebrate-parasitic nematodes and are differentially expressed during IJ activation, suggesting conserved functions in nematode parasitism. In summary, our findings strongly support a new model that S. carpocapsae and likely other Steinernema EPNs have a more active role in contributing to the pathogenicity of the nematode-bacterium complex than simply relying on their symbiotic bacteria. Furthermore, we propose that EPNs are a good model system for investigating vertebrate- and human-parasitic nematodes, especially regarding the function of excretory/secretory products.

Host- and Helminth-Derived Endocannabinoids That Have Effects on Host Immunity Are Generated during Infection.

Helminths have coevolved with their hosts, resulting in the development of specialized host immune mechanisms and parasite-specific regulatory products. Identification of new pathways that regulate helminth infection could provide a better understanding of host-helminth interaction and may identify new therapeutic targets for helminth infection. Here we identify the endocannabinoid system as a new mechanism that influences host immunity to helminths. Endocannabinoids are lipid-derived signaling molecules that control important physiologic processes, such as feeding behavior and metabolism. Following murine infection with Nippostrongylus brasiliensis, an intestinal nematode with a life cycle similar to that of hookworms, we observed increased levels of endocannabinoids (2-arachidonoylglycerol [2-AG] or anandamide [AEA]) and the endocannabinoid-like molecule oleoylethanolamine (OEA) in infected lung and intestine. To investigate endocannabinoid function in helminth infection, we employed pharmacological inhibitors of cannabinoid subtype receptors 1 and 2 (CB1R and CB2R). Compared to findings for vehicle-treated mice, inhibition of CB1R but not CB2R resulted in increased N. brasiliensis worm burden and egg output, associated with significantly decreased expression of the T helper type 2 cytokine interleukin 5 (IL-5) in intestinal tissue and splenocyte cultures. Strikingly, bioinformatic analysis of genomic and transcriptome sequencing (RNA-seq) data sets identified putative genes encoding endocannabinoid biosynthetic and degradative enzymes in many parasitic nematodes. To test the novel hypothesis that helminth parasites produce their own endocannabinoids, we measured endocannabinoid levels in N. brasiliensis by mass spectrometry and quantitative PCR and found that N. brasiliensis parasites produced endocannabinoids, especially at the infectious larval stage. To our knowledge, this is the first report of helminth- and host-derived endocannabinoids that promote host immune responses and reduce parasite burden.

Phenotype and function of IL-10 producing NK cells in individuals with malaria experience.

Plasmodium falciparum infection can trigger high levels of inflammation that lead to fever and sometimes severe disease. People living in malaria-endemic areas gradually develop resistance to symptomatic malaria and control both parasite numbers and the inflammatory response. We previously found that adaptive natural killer (NK) cells correlate with reduced parasite load and protection from symptoms. We also previously found that murine NK cell production of IL-10 can protect mice from experimental cerebral malaria. Human NK cells can also secrete IL-10, but it was unknown what NK cell subsets produce IL-10 and if this is affected by malaria experience. We hypothesize that NK cell immunoregulation may lower inflammation and reduce fever induction. Here, we show that NK cells from subjects with malaria experience make significantly more IL-10 than subjects with no malaria experience. We then determined the proportions of NK cells that are cytotoxic and produce interferon gamma and/or IL-10 and identified a signature of adaptive and checkpoint molecules on IL-10-producing NK cells. Lastly, we find that co-culture with primary monocytes, Plasmodium -infected RBCs, and antibody induces IL-10 production by NK cells. These data suggest that NK cells may contribute to protection from malaria symptoms via IL-10 production.

Effects of the SLICK1 mutation in PRLR on regulation of core body temperature and global gene expression in liver in cattle.

The SLICK1 mutation in bovine PRLR (c.1382del; rs517047387) is a deletion mutation resulting in a protein with a truncated intracellular domain. Cattle carrying at least one allele have a phenotype characterized by a short hair coat (slick phenotype) and increased resistance to heat stress. Given the pleiotropic nature of prolactin, the mutation may affect other physiological characteristics. The liver is one organ that could potentially be affected because of the expression of PRLR. The mutation is a dominant allele, and heterozygous animals have a similar hair coat to that of animals homozygous for the mutation. Present objectives were to determine whether inheritance of the SLICK1 mutation affects liver gene expression and if animals homozygous for the SLICK1 allele differ from heterozygotes in liver gene expression and regulation of body temperature during heat stress. In one experiment, rectal and ruminal temperatures were less for Holstein heifers that were heterozygous for the SLICK1 allele compared with wildtype heifers. There were 71 differentially expressed genes in liver, with 13 upregulated and 58 downregulated in SLICK1 heterozygotes. Among the ontologies characteristic of differentially expressed genes were those related to immune function and fatty acid and amino acid metabolism. In a prospective cohort study conducted with adult Senepol cattle, body temperature and hepatic gene expression were compared between animals heterozygous or homozygous for the SLICK1 mutation. There were no differences in ruminal temperatures between genotypes, rectal temperature was higher in animals homozygous for the SLICK1 mutation, and there was only one gene in liver that was differentially expressed. It was concluded that inheritance of the SLICK1 allele can exert functional changes beyond those related to hair growth although changes in liver gene expression were not extensive. Results are also consistent with the SLICK1 allele being dominant because there were few differences in phenotype between animals inheriting one or two copies of the allele.

Systemic lipolysis promotes physiological fitness in Drosophila melanogaster.

Since interventions such as caloric restriction or fasting robustly promote lipid catabolism and improve aging-related phenotypical markers, we investigated the direct effect of increased lipid catabolism via overexpression of bmm (brummer, FBgn0036449), the major triglyceride hydrolase in Drosophila, on lifespan and physiological fitness. Comprehensive characterization was carried out using RNA-seq, lipidomics and metabolomics analysis. Global overexpression of bmm strongly promoted numerous markers of physiological fitness, including increased female fecundity, fertility maintenance, preserved locomotion activity, increased mitochondrial biogenesis and oxidative metabolism. Increased bmm robustly upregulated the heat shock protein 70 (Hsp70) family of proteins, which equipped the flies with higher resistance to heat, cold, and ER stress via improved proteostasis. Despite improved physiological fitness, bmm overexpression did not extend lifespan. Taken together, these data show that bmm overexpression has broad beneficial effects on physiological fitness, but these effects did not impact lifespan.

Virus-induced transposable element expression up-regulation in human and mouse host cells.

Virus-host cell interactions initiate a host cell-defensive response during virus infection. How transposable elements in the host cell respond to viral stress at the molecular level remains largely unclear. By reanalyzing next generation sequencing data sets from dozens of virus infection studies from the Gene Expression Omnibus database, we found that genome-wide transposon expression up-regulation in host cells occurs near antiviral response genes and exists in all studies regardless of virus, species, and host cell tissue types. Some transposons were found to be up-regulated almost immediately upon infection and before increases in virus replication and significant increases in interferon ß expression. These findings indicate that transposon up-regulation is a common phenomenon during virus infection in human and mouse and that early up-regulated transposons are part of the first wave response during virus infection.

Hybrid Assembly of the Genome of the Entomopathogenic Nematode Steinernema carpocapsae Identifies the X-Chromosome.

Entomopathogenic nematodes from the genus Steinernema are lethal insect parasites that quickly kill their insect hosts with the help of their symbiotic bacteria. Steinernema carpocapsae is one of the most studied entomopathogens due to its broad lethality to diverse insect species and its effective commercial use as a biological control agent for insect pests, as well as a genetic model for studying parasitism, pathogenesis, and symbiosis. In this study, we used long-reads from the Pacific Biosciences platform and BioNano Genomics Irys system to assemble the most complete genome of the S. carpocapsae ALL strain to date, comprising 84.5 Mb in 16 scaffolds, with an N50 of 7.36 Mb. The largest scaffold, with 20.9 Mb, was identified as chromosome X based on sex-specific genome sequencing. The high level of contiguity allowed us to characterize gene density, repeat content, and GC content. RNA-seq data from 17 developmental stages, spanning from embryo to adult, were used to predict 30,957 gene models. Using this improved genome, we performed a macrosyntenic analysis to Caenorhabditis elegans and Pristionchus pacificus and found S. carpocapsae's chromosome X to be primarily orthologous to C. elegans' and P. pacificus' chromosome II and IV. We also investigated the expansion of protein families and gene expression differences between adult male and female stage nematodes. This new genome and more accurate set of annotations provide a foundation for additional comparative genomic and gene expression studies within the Steinernema clade and across the Nematoda phylum.

Adapting the Smart-seq2 Protocol for Robust Single Worm RNA-seq.

Most nematodes are small worms that lack enough RNA for regular RNA-seq protocols without pooling hundred to thousand of individuals. We have adapted the Smart-seq2 protocol in order to sequence the transcriptome of an individual worm. While developed for individual Steinernema carpocapsae and Caenorhabditis elegans larvae as well as embryos, the protocol should be adaptable for other nematode species and small invertebrates. In addition, we describe how to analyze the RNA-seq results using the Galaxy online environment. We expect that this method will be useful for the studying gene expression variances of individual nematodes in wild type and mutant backgrounds.

Comparative Transcriptomics of Steinernema and Caenorhabditis Single Embryos Reveals Orthologous Gene Expression Convergence during Late Embryogenesis.

Cells express distinct sets of genes in a precise spatio-temporal manner during embryonic development. There is a wealth of information on the deterministic embryonic development of Caenorhabditis elegans, but much less is known about embryonic development in nematodes from other taxa, especially at the molecular level. We are interested in insect pathogenic nematodes from the genus Steinernema as models of parasitism and symbiosis as well as a satellite model for evolution in comparison to C. elegans. To explore gene expression differences across taxa, we sequenced the transcriptomes of single embryos of two Steinernema species and two Caenorhabditis species at 11 stages during embryonic development and found several interesting features. Our findings show that zygotic transcription initiates at different developmental stages in each species, with the Steinernema species initiating transcription earlier than Caenorhabditis. We found that ortholog expression conservation during development is higher at the later embryonic stages than at the earlier ones. The surprisingly higher conservation of orthologous gene expression in later embryonic stages strongly suggests a funnel-shaped model of embryonic developmental gene expression divergence in nematodes. This work provides novel insight into embryonic development across distantly related nematode species and demonstrates that the mechanisms controlling early development are more diverse than previously thought at the transcriptional level.

Comparative genomics of Steinernema reveals deeply conserved gene regulatory networks.

BACKGROUND: Parasitism is a major ecological niche for a variety of nematodes. Multiple nematode lineages have specialized as pathogens, including deadly parasites of insects that are used in biological control. We have sequenced and analyzed the draft genomes and transcriptomes of the entomopathogenic nematode Steinernema carpocapsae and four congeners (S. scapterisci, S. monticolum, S. feltiae, and S. glaseri). RESULTS: We used these genomes to establish phylogenetic relationships, explore gene conservation across species, and identify genes uniquely expanded in insect parasites. Protein domain analysis in Steinernema revealed a striking expansion of numerous putative parasitism genes, including certain protease and protease inhibitor families, as well as fatty acid- and retinol-binding proteins. Stage-specific gene expression of some of these expanded families further supports the notion that they are involved in insect parasitism by Steinernema. We show that sets of novel conserved non-coding regulatory motifs are associated with orthologous genes in Steinernema and Caenorhabditis. CONCLUSIONS: We have identified a set of expanded gene families that are likely to be involved in parasitism. We have also identified a set of non-coding motifs associated with groups of orthologous genes in Steinernema and Caenorhabditis involved in neurogenesis and embryonic development that are likely part of conserved protein-DNA relationships shared between these two genera.