TLR2 Supports γδ T cell IL-17A Response to ocular surface commensals by Metabolic Reprogramming.

The ocular surface is a mucosal barrier tissue colonized by commensal microbes, which tune local immunity by eliciting IL-17 from conjunctival γδ T cells to prevent pathogenic infection. The commensal Corynebacterium mastitidis (C. mast) elicits protective IL-17 responses from conjunctival Vγ4 T cells through a combination of γδ TCR ligation and IL-1 signaling. Here, we identify Vγ6 T cells as a major C. mast-responsive subset in the conjunctiva and uncover its unique activation requirements. We demonstrate that Vγ6 cells require not only extrinsic (via dendritic cells) but also intrinsic TLR2 stimulation for optimal IL-17A response. Mechanistically, intrinsic TLR2 signaling was associated with epigenetic changes and enhanced expression of genes responsible for metabolic shift to fatty acid oxidation to support Il17a transcription. We identify one key transcription factor, IκBζ, which is upregulated by TLR2 stimulation and is essential for this program. Our study highlights the importance of intrinsic TLR2 signaling in driving metabolic reprogramming and production of IL-17A in microbiome-specific mucosal γδ T cells.

A CTCF-binding site in the Mdm1-Il22-Ifng locus shapes cytokine expression profiles and plays a critical role in early Th1 cell fate specification.

Cytokine expression during T cell differentiation is a highly regulated process that involves long-range promoter-enhancer and CTCF-CTCF contacts at cytokine loci. Here, we investigated the impact of dynamic chromatin loop formation within the topologically associating domain (TAD) in regulating the expression of interferon gamma (IFN-γ) and interleukin-22 (IL-22); these cytokine loci are closely located in the genome and are associated with complex enhancer landscapes, which are selectively active in type 1 and type 3 lymphocytes. In situ Hi-C analyses revealed inducible TADs that insulated Ifng and Il22 enhancers during Th1 cell differentiation. Targeted deletion of a 17 bp boundary motif of these TADs imbalanced Th1- and Th17-associated immunity, both in vitro and in vivo, upon Toxoplasma gondii infection. In contrast, this boundary element was dispensable for cytokine regulation in natural killer cells. Our findings suggest that precise cytokine regulation relies on lineage- and developmental stage-specific interactions of 3D chromatin architectures and enhancer landscapes.

Too Much of a Good Thing: Extended Duration of Gut Microbiota Depletion Reverses Protection From Experimental Autoimmune Uveitis.

Purpose: Using the model of experimental autoimmune uveitis (EAU) induced by immunization with a retinal antigen, two studies have reported contradictory results on disease development following oral antibiotic treatment (ABX). We showed that long-term ABX did not affect EAU, but another study showed that short-term ABX was protective. We therefore studied the effects of ABX on EAU, gut microbiota, and host immune responses as a function of treatment duration. Methods: EAU-susceptible mice were treated orally with broad-spectrum antibiotics starting at least 10 weeks (long-term) or 1 week (short-term) before immunization until termination of the experiment. Gut microbiota were characterized by 16S amplicon sequencing, and host gut immune elements were studied phenotypically and functionally. Results: Long-term ABX had no effect, whereas short-term ABX delayed EAU, as previously reported by us and others, respectively. Microbial sequencing revealed progressive reduction of gut microbiota that showed some differences in the two ABX groups. Interestingly, duration of ABX was associated with a gradual disappearance of the CD4+ and CD4+CD8+ subset of gut intraepithelial lymphocytes (IELs). This IEL subset is microbiota dependent and is absent in germ-free mice. Relative abundance of Lactobacillus reuteri correlated with the frequencies of CD4+CD8+ IELs. IELs suppressed antigen-specific activation of autoreactive T cells in culture. Conclusions: Gut microbiota may play dual roles in uveitis development: They promote EAU development but also help maintain IEL populations that have regulatory function against autoreactive T cells. We propose that the progressive loss of this population during long-term ABX reverses the EAU-ameliorating effects of microbiota depletion.

Single-cell multiomic analysis reveals the involvement of Type I interferon-responsive CD8+ T cells in amyloid beta-associated memory loss.

Alzheimer's disease (AD) is the leading cause of dementia worldwide, but there are limited therapeutic options and no current cure. While the involvement of microglia in AD has been highly appreciated, the role of other innate and adaptive immune cells remains largely unknown, partly due to their scarcity and heterogeneity. This study aimed to study non-microglial immune cells in wild type and AD-transgenic mouse brains across different ages. Our results uncovered the presence of a unique CD8+ T cell population that were selectively increased in aging AD mouse brains, here referred to as "disease-associated T cells (DATs)". These DATs were found to express an elevated tissue-resident memory and Type I interferon-responsive gene signature. Further analysis of aged AD mouse brains showed that these CD8+ T cells were not present in peripheral or meningeal tissues. Preventing CD8+ T cell development in AD-transgenic mice via genetic deletion of beta-2 microglobulin ( B2m ) led to a reduction of amyloid-ß plaque formation in aged mice, and improved memory in AD-transgenic mice as early as four months of age. The integration of transcriptomic and epigenomic profiles at the single-cell level revealed potential transcription factors that reshape the regulomes of CD8+ T cells. These findings highlight a critical role for DATs in the progression of AD and provide a new avenue for treatment.

Draft Reference Genome Sequence of Corynebacterium mastitidis RC, an Ocular Commensal, Isolated from Mouse Conjunctiva.

Here, we report the genome sequence of a protective commensal, Corynebacterium mastitidis RC, isolated from mouse conjunctiva. The C. mastitidis RC genome sequence is 2,153,054 bp in size and 96.95% complete, and we believe that it can contribute to the understanding of the functional immune attributes of the ocular commensal microbiome.

Profiling invasive Plasmodium falciparum merozoites using an integrated omics approach.

The symptoms of malaria are brought about by blood-stage parasites, which are established when merozoites invade human erythrocytes. Our understanding of the molecular events that underpin erythrocyte invasion remains hampered by the short-period of time that merozoites are invasive. To address this challenge, a Plasmodium falciparum gamma-irradiated long-lived merozoite (LLM) line was developed and investigated. Purified LLMs invaded erythrocytes by an increase of 10-300 fold compared to wild-type (WT) merozoites. Using an integrated omics approach, we investigated the basis for the phenotypic difference. Only a few single nucleotide polymorphisms within the P. falciparum genome were identified and only marginal differences were observed in the merozoite transcriptomes. By contrast, using label-free quantitative mass-spectrometry, a significant change in protein abundance was noted, of which 200 were proteins of unknown function. We determined the relative molar abundance of over 1100 proteins in LLMs and further characterized the major merozoite surface protein complex. A unique processed MSP1 intermediate was identified in LLM but not observed in WT suggesting that delayed processing may be important for the observed phenotype. This integrated approach has demonstrated the significant role of the merozoite proteome during erythrocyte invasion, while identifying numerous unknown proteins likely to be involved in invasion.

A Plasmodium yoelii HECT-like E3 ubiquitin ligase regulates parasite growth and virulence.

Infection of mice with strains of Plasmodium yoelii parasites can result in different pathology, but molecular mechanisms to explain this variation are unclear. Here we show that a P. yoelii gene encoding a HECT-like E3 ubiquitin ligase (Pyheul) influences parasitemia and host mortality. We genetically cross two lethal parasites with distinct disease phenotypes, and identify 43 genetically diverse progeny by typing with microsatellites and 9230 single-nucleotide polymorphisms. A genome-wide quantitative trait loci scan links parasite growth and host mortality to two major loci on chromosomes 1 and 7 with LOD (logarithm of the odds) scores = 6.1 and 8.1, respectively. Allelic exchange of partial sequences of Pyheul in the chromosome 7 locus and modification of the gene expression alter parasite growth and host mortality. This study identifies a gene that may have a function in parasite growth, virulence, and host-parasite interaction, and therefore could be a target for drug or vaccine development.Many strains of Plasmodium differ in virulence, but factors that control these distinctions are not known. Here the authors comparatively map virulence loci using the offspring from a P. yoelii YM and N67 genetic cross, and identify a putative HECT E3 ubiquitin ligase that may explain the variance.

Ethanol-induced stress response of Staphylococcus aureus.

Transcriptional profiles of 2 unrelated clinical methicillin-resistant Staphylococcus aureus (MRSA) isolates were analyzed following 10% (v/v) ethanol challenge (15 min), which arrested growth but did not reduce viability. Ethanol-induced stress (EIS) resulted in differential gene expression of 1091 genes, 600 common to both strains, of which 291 were upregulated. With the exception of the downregulation of genes involved with osmotic stress functions, EIS resulted in the upregulation of genes that contribute to stress response networks, notably those altered by oxidative stress, protein quality control in general, and heat shock in particular. In addition, genes involved with transcription, translation, and nucleotide biosynthesis were downregulated. relP, which encodes a small alarmone synthetase (RelP), was highly upregulated in both MRSA strains following ethanol challenge, and relP inactivation experiments indicated that this gene contributed to EIS growth arrest. A number of persistence-associated genes were also upregulated during EIS, including those that encode toxin-antitoxin systems. Overall, transcriptional profiling indicated that the MRSA investigated responded to EIS by entering a state of dormancy and by altering the expression of elements from cross protective stress response systems in an effort to protect preexisting proteins.

MetaNetVar: Pipeline for applying network analysis tools for genomic variants analysis.

Network analysis can make variant analysis better. There are existing tools like HotNet2 and dmGWAS that can provide various analytical methods. We developed a prototype of a pipeline called MetaNetVar that allows execution of multiple tools. The code is published at https://github.com/NCBI-Hackathons/Network_SNPs. A working prototype is published as an Amazon Machine Image - ami-4510312f .

Comparative analyses of CTCF and BORIS occupancies uncover two distinct classes of CTCF binding genomic regions.

BACKGROUND: CTCF and BORIS (CTCFL), two paralogous mammalian proteins sharing nearly identical DNA binding domains, are thought to function in a mutually exclusive manner in DNA binding and transcriptional regulation. RESULTS: Here we show that these two proteins co-occupy a specific subset of regulatory elements consisting of clustered CTCF binding motifs (termed 2xCTSes). BORIS occupancy at 2xCTSes is largely invariant in BORIS-positive cancer cells, with the genomic pattern recapitulating the germline-specific BORIS binding to chromatin. In contrast to the single-motif CTCF target sites (1xCTSes), the 2xCTS elements are preferentially found at active promoters and enhancers, both in cancer and germ cells. 2xCTSes are also enriched in genomic regions that escape histone to protamine replacement in human and mouse sperm. Depletion of the BORIS gene leads to altered transcription of a large number of genes and the differentiation of K562 cells, while the ectopic expression of this CTCF paralog leads to specific changes in transcription in MCF7 cells. CONCLUSIONS: We discover two functionally and structurally different classes of CTCF binding regions, 2xCTSes and 1xCTSes, revealed by their predisposition to bind BORIS. We propose that 2xCTSes play key roles in the transcriptional program of cancer and germ cells.

PLZF regulates CCR6 and is critical for the acquisition and maintenance of the Th17 phenotype in human cells.

Th17 cells, which express the chemokine receptor CCR6, are implicated in many immune-mediated disorders, such as psoriasis and multiple sclerosis. We found that expression levels of CCR6 on human effector/memory CD4(+) T cells reflect a continuum of Th17 differentiation. By evaluating the transcriptome in cells with increasing CCR6, we detected progressive upregulation of ZBTB16, which encodes the broad complex, tramtrack, bric-à-brac-zinc finger transcription factor promyelocytic leukemia zinc finger protein (PLZF). Using chromatin immunoprecipitation for modified histones, p300, and PLZF, we identified enhancer-like sites at -9/-10 and -13/-14 kb from the upstream transcription start site of CCR6 that bind PLZF in CCR6(+) cells. For Th cells from adult blood, both in the CCR6(+) memory population and in naive cells activated ex vivo, knockdown of ZBTB16 downregulated CCR6 and other Th17-associated genes. ZBTB16 and RORC (which encodes the "master regulator" RORγt) cross-regulate each other, and PLZF binds at the RORC promoter in CCR6(+) cells. In naive Th cells from cord blood, ZBTB16 expression was confined to CD161(+) cells, which are Th17 cell precursors. ZBTB16 was not expressed in mouse Th17 cells, and Th17 cells could be made from luxoid mice, which harbor an inactivating mutation in Zbtb16. These studies demonstrate a role for PLZF as an activator of transcription important both for Th17 differentiation and the maintenance of the Th17 phenotype in human cells, expand the role of PLZF as a critical regulator in the human adaptive immune system, and identify a novel, essential element in a regulatory network that is of significant therapeutic interest.

SATRAT: Staphylococcus aureus transcript regulatory network analysis tool.

Staphylococcus aureus is a commensal organism that primarily colonizes the nose of healthy individuals. S. aureus causes a spectrum of infections that range from skin and soft-tissue infections to fatal invasive diseases. S. aureus uses a large number of virulence factors that are regulated in a coordinated fashion. The complex regulatory mechanisms have been investigated in numerous high-throughput experiments. Access to this data is critical to studying this pathogen. Previously, we developed a compilation of microarray experimental data to enable researchers to search, browse, compare, and contrast transcript profiles. We have substantially updated this database and have built a novel exploratory tool-SATRAT-the S. aureus transcript regulatory network analysis tool, based on the updated database. This tool is capable of performing deep searches using a query and generating an interactive regulatory network based on associations among the regulators of any query gene. We believe this integrated regulatory network analysis tool would help researchers explore the missing links and identify novel pathways that regulate virulence in S. aureus. Also, the data model and the network generation code used to build this resource is open sourced, enabling researchers to build similar resources for other bacterial systems.

Unipro UGENE NGS pipelines and components for variant calling, RNA-seq and ChIP-seq data analyses.

The advent of Next Generation Sequencing (NGS) technologies has opened new possibilities for researchers. However, the more biology becomes a data-intensive field, the more biologists have to learn how to process and analyze NGS data with complex computational tools. Even with the availability of common pipeline specifications, it is often a time-consuming and cumbersome task for a bench scientist to install and configure the pipeline tools. We believe that a unified, desktop and biologist-friendly front end to NGS data analysis tools will substantially improve productivity in this field. Here we present NGS pipelines "Variant Calling with SAMtools", "Tuxedo Pipeline for RNA-seq Data Analysis" and "Cistrome Pipeline for ChIP-seq Data Analysis" integrated into the Unipro UGENE desktop toolkit. We describe the available UGENE infrastructure that helps researchers run these pipelines on different datasets, store and investigate the results and re-run the pipelines with the same parameters. These pipeline tools are included in the UGENE NGS package. Individual blocks of these pipelines are also available for expert users to create their own advanced workflows.

Nucleotide excision repair/transcription gene defects in the fetus and impaired TFIIH-mediated function in transcription in placenta leading to preeclampsia.

BACKGROUND: Preeclampsia is a significant cause of maternal and fetal mortality and morbidity worldwide. We previously reported associations between trichothiodystrophy (TTD) nucleotide excision repair (NER) and transcription gene mutations in the fetus and the risk of gestational complications including preeclampsia. TTD NER/transcription genes, XPD, XPB and TTD-A, code for subunits of Transcription Factor (TF)IIH. Interpreting XPD mutations in the context of available biochemical data led us to propose adverse effects on CDK-activating kinase (CAK) subunit of TFIIH and TFIIH-mediated functions as a relevant mechanism in preeclampsia. In order to gain deeper insight into the underlying biologic mechanisms involving TFIIH-mediated functions in placenta, we analyzed NER/transcription and global gene expression profiles of normal and preeclamptic placentas and studied gene regulatory networks. RESULTS: We found high expression of TTD NER/transcription genes in normal human placenta, above the mean of their expression in all organs. XPD and XPB were consistently expressed from 14 to 40 weeks gestation while expression of TTD-A was strongly negatively correlated (r=-0.7, P<0.0001) with gestational age. Analysis of gene expression patterns of placentas from a case-control study of preeclampsia using Algorithm for Reconstruction of Accurate Cellular Networks (ARACNE) revealed GTF2E1, a component of TFIIE which modulates TFIIH, among major regulators of differentially-expressed genes in preeclampsia. The basal transcription pathway was among the largest dysregulated protein-protein interaction networks in this preeclampsia dataset. Within the basal transcription pathway, significantly down-regulated genes besides GTF2E1 included those coding for the CAK complex of TFIIH, namely CDK7, CCNH, and MNAT1. Analysis of other relevant gene expression and gene regulatory network data also underscored the involvement of transcription pathways and identified JUNB and JUND (components of transcription factor AP-1) as transcription regulators of the network involving the TTD genes, GTF2E1, and selected gene regulators implicated in preeclampsia. CONCLUSIONS: Our results indicate that TTD NER/transcription genes are expressed in placenta during gestational periods critical to preeclampsia development. Our overall findings suggest that impairment of TFIIH-mediated function in transcription in placenta is a likely mechanism leading to preeclampsia and provide etiologic clues which may be translated into therapeutic and preventive measures.

miR-155 activates cytokine gene expression in Th17 cells by regulating the DNA-binding protein Jarid2 to relieve polycomb-mediated repression.

Specification of the T helper 17 (Th17) cell lineage requires a well-defined set of transcription factors, but how these integrate with posttranscriptional and epigenetic programs to regulate gene expression is poorly understood. Here we found defective Th17 cell cytokine expression in miR-155-deficient CD4+ T cells in vitro and in vivo. Mir155 was bound by Th17 cell transcription factors and was highly expressed during Th17 cell differentiation. miR-155-deficient Th17 and T regulatory (Treg) cells expressed increased amounts of Jarid2, a DNA-binding protein that recruits the Polycomb Repressive Complex 2 (PRC2) to chromatin. PRC2 binding to chromatin and H3K27 histone methylation was increased in miR-155-deficient cells, coinciding with failure to express Il22, Il10, Il9, and Atf3. Defects in Th17 cell cytokine expression and Treg cell homeostasis in the absence of Mir155 could be partially suppressed by Jarid2 deletion. Thus, miR-155 contributes to Th17 cell function by suppressing the inhibitory effects of Jarid2.

Genome-wide polymorphisms and development of a microarray platform to detect genetic variations in Plasmodium yoelii.

The rodent malaria parasite Plasmodium yoelii is an important model for studying malaria immunity and pathogenesis. One approach for studying malaria disease phenotypes is genetic mapping, which requires typing a large number of genetic markers from multiple parasite strains and/or progeny from genetic crosses. Hundreds of microsatellite (MS) markers have been developed to genotype the P. yoelii genome; however, typing a large number of MS markers can be labor intensive, time consuming, and expensive. Thus, development of high-throughput genotyping tools such as DNA microarrays that enable rapid and accurate large-scale genotyping of the malaria parasite will be highly desirable. In this study, we sequenced the genomes of two P. yoelii strains (33X and N67) and obtained a large number of single nucleotide polymorphisms (SNPs). Based on the SNPs obtained, we designed sets of oligonucleotide probes to develop a microarray that could interrogate ∼11,000 SNPs across the 14 chromosomes of the parasite in a single hybridization. Results from hybridizations of DNA samples of five P. yoelii strains or cloned lines (17XNL, YM, 33X, N67 and N67C) and two progeny from a genetic cross (N67×17XNL) to the microarray showed that the array had a high call rate (∼97%) and accuracy (99.9%) in calling SNPs, providing a simple and reliable tool for typing the P. yoelii genome. Our data show that the P. yoelii genome is highly polymorphic, although isogenic pairs of parasites were also detected. Additionally, our results indicate that the 33X parasite is a progeny of 17XNL (or YM) and an unknown parasite. The highly accurate and reliable microarray developed in this study will greatly facilitate our ability to study the genetic basis of important traits and the disease it causes.

Development of genetic system to inactivate a Borrelia turicatae surface protein selectively produced within the salivary glands of the arthropod vector.

BACKGROUND: Borrelia turicatae, an agent of tick-borne relapsing fever, is an example of a pathogen that can adapt to disparate conditions found when colonizing the mammalian host and arthropod vector. However, little is known about the genetic factors necessary during the tick-mammalian infectious cycle, therefore we developed a genetic system to transform this species of spirochete. We also identified a plasmid gene that was up-regulated in vitro when B. turicatae was grown in conditions mimicking the tick environment. This 40 kilodalton protein was predicted to be surface localized and designated the Borrelia repeat protein A (brpA) due to the redundancy of the amino acid motif Gln-Gly-Asn-Val-Glu. METHODOLOGY/PRINCIPAL FINDINGS: Quantitative reverse-transcriptase polymerase chain reaction using RNA from B. turicatae infected ticks and mice indicated differential regulation of brpA during the tick-mammalian infectious cycle. The surface localization was determined, and production of the protein within the salivary glands of the tick was demonstrated. We then applied a novel genetic system for B. turicatae to inactivate brpA and examined the role of the gene product for vector colonization and the ability to establish murine infection. CONCLUSIONS/SIGNIFICANCE: These results demonstrate the complexity of protein production in a population of spirochetes within the tick. Additionally, the development of a genetic system is important for future studies to evaluate the requirement of specific B. turicatae genes for vector colonization and transmission.

Tea tree oil-induced transcriptional alterations in Staphylococcus aureus.

Tea tree oil (TTO) is a steam distillate of Melaleuca alternifolia that demonstrates broad-spectrum antibacterial activity. This study was designed to document how TTO challenge influences the Staphylococcus aureus transcriptome. Overall, bioinformatic analyses (S. aureus microarray meta-database) revealed that both ethanol and TTO induce related transcriptional alterations. TTO challenge led to the down-regulation of genes involved with energy-intensive transcription and translation, and altered the regulation of genes involved with heat shock (e.g. clpC, clpL, ctsR, dnaK, groES, groEL, grpE and hrcA) and cell wall metabolism (e.g. cwrA, isaA, sle1, vraSR and vraX). Inactivation of the heat shock gene dnaK or vraSR which encodes a two-component regulatory system that responds to peptidoglycan biosynthesis inhibition led to an increase in TTO susceptibility which demonstrates a protective role for these genes in the S. aureus TTO response. A gene (mmpL) encoding a putative resistance, nodulation and cell division efflux pump was also highly induced by TTO. The principal antimicrobial TTO terpene, terpinen-4-ol, altered ten genes in a transcriptional direction analogous to TTO. Collectively, this study provides additional insight into the response of a bacterial pathogen to the antimicrobial terpene mixture TTO.

Genome-wide SNP analysis of the Systemic Capillary Leak Syndrome (Clarkson disease).

The Systemic Capillary Leak Syndrome (SCLS) is an extremely rare, orphan disease that resembles, and is frequently erroneously diagnosed as, systemic anaphylaxis. The disorder is characterized by repeated, transient, and seemingly unprovoked episodes of hypotensive shock and peripheral edema due to transient endothelial hyperpermeability. SCLS is often accompanied by a monoclonal gammopathy of unknown significance (MGUS). Using Affymetrix Single Nucleotide Polymorphism (SNP) microarrays, we performed the first genome-wide SNP analysis of SCLS in a cohort of 12 disease subjects and 18 controls. Exome capture sequencing was performed on genomic DNA from nine of these patients as validation for the SNP-chip discoveries and de novo data generation. We identified candidate susceptibility loci for SCLS, which included a region flanking CAV3 (3p25.3) as well as SNP clusters in PON1 (7q21.3), PSORS1C1 (6p21.3), and CHCHD3 (7q33). Among the most highly ranked discoveries were gene-associated SNPs in the uncharacterized LOC100130480 gene (rs6417039, rs2004296). Top case-associated SNPs were observed in BTRC (rs12355803, 3rs4436485), ARHGEF18 (rs11668246), CDH13 (rs4782779), and EDG2 (rs12552348), which encode proteins with known or suspected roles in B cell function and/or vascular integrity. 61 SNPs that were significantly associated with SCLS by microarray analysis were also detected and validated by exome deep sequencing. Functional annotation of highly ranked SNPs revealed enrichment of cell projections, cell junctions and adhesion, and molecules containing pleckstrin homology, Ras/Rho regulatory, and immunoglobulin Ig-like C2/fibronectin type III domains, all of which involve mechanistic functions that correlate with the SCLS phenotype. These results highlight SNPs with potential relevance to SCLS.

The fat body transcriptomes of the yellow fever mosquito Aedes aegypti, pre- and post- blood meal.

BACKGROUND: The fat body is the main organ of intermediary metabolism in insects and the principal source of hemolymph proteins. As part of our ongoing efforts to understand mosquito fat body physiology and to identify novel targets for insect control, we have conducted a transcriptome analysis of the fat body of Aedes aegypti before and in response to blood feeding. RESULTS: We created two fat body non-normalized EST libraries, one from mosquito fat bodies non-blood fed (NBF) and another from mosquitoes 24 hrs post-blood meal (PBM). 454 pyrosequencing of the non-normalized libraries resulted in 204,578 useable reads from the NBF sample and 323,474 useable reads from the PBM sample. Alignment of reads to the existing reference Ae. aegypti transcript libraries for analysis of differential expression between NBF and PBM samples revealed 116,912 and 115,051 matches, respectively. De novo assembly of the reads from the NBF sample resulted in 15,456 contigs, and assembly of the reads from the PBM sample resulted in 15,010 contigs. Collectively, 123 novel transcripts were identified within these contigs. Prominently expressed transcripts in the NBF fat body library were represented by transcripts encoding ribosomal proteins. Thirty-five point four percent of all reads in the PBM library were represented by transcripts that encode yolk proteins. The most highly expressed were transcripts encoding members of the cathepsin b, vitellogenin, vitellogenic carboxypeptidase, and vitelline membrane protein families. CONCLUSION: The two fat body transcriptomes were considerably different from each other in terms of transcript expression in terms of abundances of transcripts and genes expressed. They reflect the physiological shift of the pre-feeding fat body from a resting state to vitellogenic gene expression after feeding.