Conventional Co-Housing Modulates Murine Gut Microbiota and Hematopoietic Gene Expression.

Specific-pathogen-free (SPF) mice have improved hematopoietic characteristics relative to germ-free mice, however, it is not clear whether improvements in hematopoietic traits will continue when the level of microorganism exposure is further increased. We co-housed SPF C57BL/6 mice in a conventional facility (CVT) and found a significant increase in gut microbiota diversity along with increased levels of myeloid cells and T cells, especially effector memory T cells. Through single cell RNA sequencing of sorted KL (c-Kit+Lin-) cells, we imputed a decline in long-term hematopoietic stem cells and an increase in granulocyte-monocyte progenitors in CVT mice with up-regulation of genes associated with cell survival. Bone marrow transplantation through competitive repopulation revealed a significant increase in KSL (c-Kit+Sca-1+Lin-) cell reconstitution in recipients of CVT donor cells which occurred when donors were co-housed for both one and twelve months. However, there was minimal to no gain in mature blood cell engraftment in recipients of CVT donor cells relative to those receiving SPF donor cells. We conclude that co-housing SPF mice with mice born in a conventional facility increased gut microbiota diversity, augmented myeloid cell production and T cell activation, stimulated KSL cell reconstitution, and altered hematopoietic gene expression.

TREML4 Promotes Inflammatory Programs in Human and Murine Macrophages and Alters Atherosclerosis Lesion Composition in the Apolipoprotein E Deficient Mouse.

The Triggering Receptor Expressed on Myeloid cells-like 4 (TREML4) is a member of the TREM receptor family, known modulators of inflammatory responses. We have previously found that TREML4 expression positively correlates with human coronary arterial calcification (CAC). However, the role of TREML4 in the pathogenesis of cardiovascular disease remains incompletely defined. Since macrophages play a key role in inflammatory conditions, we investigated if activated macrophages selectively expressed TREML4 and found that carriage of either one of the eQTL SNP's previously associated with increased TREML4 expression conferred higher expression in human inflammatory macrophages (M1) compared to alternatively activated macrophages (M2). Furthermore, we found that TREML4 expression in human M1 dysregulated several inflammatory pathways related to leukocyte activation, apoptosis and extracellular matrix degradation. Similarly, murine M1 expressed substantial levels of Treml4, as did oxLDL treated macrophages. Transcriptome analysis confirmed that murine Treml4 controls the expression of genes related to inflammation and lipid regulation pathways, suggesting a possible role in atherosclerosis. Analysis of Apoe-/-/Treml4-/- mice showed reduced plaque burden and lesion complexity as indicated by decreased stage scores, macrophage content and collagen deposition. Finally, transcriptome analysis of oxLDL-loaded murine macrophages showed that Treml4 represses a specific set of genes related to carbohydrate, ion and amino acid membrane transport. Metabolomic analysis confirmed that Treml4 deficiency may promote a beneficial relationship between iron homeostasis and glucose metabolism. Together, our results suggest that Treml4 plays a role in the development of cardiovascular disease, as indicated by Treml4-dependent dysregulation of macrophage inflammatory pathways, macrophage metabolism and promotion of vulnerability features in advanced lesions.

Attenuation of immune-mediated bone marrow damage in conventionally housed mice.

In humans, bone marrow (BM) failure syndromes, both constitutional and acquired, predispose to myeloid malignancies. We have modeled acquired immune aplastic anemia, the paradigmatic disease of these syndromes, in the mouse by infusing lymph node cells from specific pathogen-free (SPF) CD45.1 congenic C57BL/6 (B6) donors into hybrid CByB6F1 recipients housed either in conventional (CVB) or SPF facilities. The severity of BM damage was reduced in CVB recipients; they also had reduced levels of CD44+ CD62L- effector memory T cells, reduced numbers of donor-type CD44+ T cells, and reduced expansion of donor-type CD8 T cells carrying T-cell receptor ß-variable regions 07, 11, and 17. Analyses of fecal samples through 16S ribosomal RNA amplicon sequencing revealed greater gut microbial alpha diversity in CVB mice relative to that of SPF mice. Thus, the presence of a broader spectrum of gut microorganisms in CVB-housed CByB6F1 could have primed recipient animal's immune system leading to suppression of allogeneic donor T-cell activation and expansion and attenuation of host BM destruction. These results suggest the potential benefit of diverse gut microbiota in patients receiving BM transplants.

Gut microbiome-mediated bile acid metabolism regulates liver cancer via NKT cells.

Primary liver tumors and liver metastasis currently represent the leading cause of cancer-related death. Commensal bacteria are important regulators of antitumor immunity, and although the liver is exposed to gut bacteria, their role in antitumor surveillance of liver tumors is poorly understood. We found that altering commensal gut bacteria in mice induced a liver-selective antitumor effect, with an increase of hepatic CXCR6+ natural killer T (NKT) cells and heightened interferon-γ production upon antigen stimulation. In vivo functional studies showed that NKT cells mediated liver-selective tumor inhibition. NKT cell accumulation was regulated by CXCL16 expression of liver sinusoidal endothelial cells, which was controlled by gut microbiome-mediated primary-to-secondary bile acid conversion. Our study suggests a link between gut bacteria-controlled bile acid metabolism and liver antitumor immunosurveillance.

Integrative DNA, RNA, and protein evidence connects TREML4 to coronary artery calcification.

Coronary artery calcification (CAC) is a heritable and definitive morphologic marker of atherosclerosis that strongly predicts risk for future cardiovascular events. To search for genes involved in CAC, we used an integrative transcriptomic, genomic, and protein expression strategy by using next-generation DNA sequencing in the discovery phase with follow-up studies using traditional molecular biology and histopathology techniques. RNA sequencing of peripheral blood from a discovery set of CAC cases and controls was used to identify dysregulated genes, which were validated by ClinSeq and Framingham Heart Study data. Only a single gene, TREML4, was upregulated in CAC cases in both studies. Further examination showed that rs2803496 was a TREML4 cis-eQTL and that the minor allele at this locus conferred up to a 6.5-fold increased relative risk of CAC. We characterized human TREML4 and demonstrated by immunohistochemical techniques that it is localized in macrophages surrounding the necrotic core of coronary plaques complicated by calcification (but not in arteries with less advanced disease). Finally, we determined by von Kossa staining that TREML4 colocalizes with areas of microcalcification within coronary plaques. Overall, we present integrative RNA, DNA, and protein evidence implicating TREML4 in coronary artery calcification. Our findings connect multimodal genomics data with a commonly used clinical marker of cardiovascular disease.

Gene expression signatures of coronary heart disease.

OBJECTIVE: To identify transcriptomic biomarkers of coronary heart disease (CHD) in 188 cases with CHD and 188 age- and sex-matched controls who were participants in the Framingham Heart Study. APPROACH AND RESULTS: A total of 35 genes were differentially expressed in cases with CHD versus controls at false discovery rate<0.5, including GZMB, TMEM56, and GUK1. Cluster analysis revealed 3 gene clusters associated with CHD, 2 linked to increased erythrocyte production and a third to reduced natural killer and T cell activity in cases with CHD. Exon-level results corroborated and extended the gene-level results. Alternative splicing analysis suggested that GUK1 and 38 other genes were differentially spliced in cases with CHD versus controls. Gene Ontology analysis linked ubiquitination and T-cell-related pathways with CHD. CONCLUSIONS: Two bioinformatically defined groups of genes show consistent associations with CHD. Our findings are consistent with the hypotheses that hematopoesis is upregulated in CHD, possibly reflecting a compensatory mechanism, and that innate immune activity is disrupted in CHD or altered by its treatment. Transcriptomic signatures may be useful in identifying pathways associated with CHD and point toward novel therapeutic targets for its treatment and prevention.

Synthetic spike-in standards improve run-specific systematic error analysis for DNA and RNA sequencing.

While the importance of random sequencing errors decreases at higher DNA or RNA sequencing depths, systematic sequencing errors (SSEs) dominate at high sequencing depths and can be difficult to distinguish from biological variants. These SSEs can cause base quality scores to underestimate the probability of error at certain genomic positions, resulting in false positive variant calls, particularly in mixtures such as samples with RNA editing, tumors, circulating tumor cells, bacteria, mitochondrial heteroplasmy, or pooled DNA. Most algorithms proposed for correction of SSEs require a data set used to calculate association of SSEs with various features in the reads and sequence context. This data set is typically either from a part of the data set being "recalibrated" (Genome Analysis ToolKit, or GATK) or from a separate data set with special characteristics (SysCall). Here, we combine the advantages of these approaches by adding synthetic RNA spike-in standards to human RNA, and use GATK to recalibrate base quality scores with reads mapped to the spike-in standards. Compared to conventional GATK recalibration that uses reads mapped to the genome, spike-ins improve the accuracy of Illumina base quality scores by a mean of 5 Phred-scaled quality score units, and by as much as 13 units at CpG sites. In addition, since the spike-in data used for recalibration are independent of the genome being sequenced, our method allows run-specific recalibration even for the many species without a comprehensive and accurate SNP database. We also use GATK with the spike-in standards to demonstrate that the Illumina RNA sequencing runs overestimate quality scores for AC, CC, GC, GG, and TC dinucleotides, while SOLiD has less dinucleotide SSEs but more SSEs for certain cycles. We conclude that using these DNA and RNA spike-in standards with GATK improves base quality score recalibration.