Interspecies comparison of the early transcriptomic changes associated with hepatitis B virus exposure in human and macaque immune cell populations.

Background and aims: Hepatitis B virus (HBV) infection affects 300 million individuals worldwide, representing a major factor for the development of hepatic complications. Although existing antivirals are effective in suppressing replication, eradication of HBV is not achieved. Therefore, a multi-faceted approach involving antivirals and immunomodulatory agents is required. Non-human primates are widely used in pre-clinical studies due to their close evolutionary relationship to humans. Nonetheless, it is fundamental to identify the differences in immune response between humans and these models. Thus, we performed a transcriptomic characterization and interspecies comparison of the early immune responses to HBV in human and cynomolgus macaques. Methods: We characterized early transcriptomic changes in human and cynomolgus B cells, T cells, myeloid and plasmacytoid dendritic cells (pDC) exposed to HBV ex vivo for 2 hours. Differentially-expressed genes were further compared to the profiles of HBV-infected patients using publicly-available single-cell data. Results: HBV induced a wide variety of transcriptional changes in all cell types, with common genes between species representing only a small proportion. In particular, interferon gamma signaling was repressed in human pDCs. At the gene level, interferon gamma inducible protein 16 (IFI16) was upregulated in macaque pDCs, while downregulated in humans. Moreover, IFI16 expression in pDCs from chronic HBV-infected patients anti-paralleled serum HBsAg levels. Conclusion: Our characterization of early transcriptomic changes induced by HBV in humans and cynomolgus macaques represents a useful resource for the identification of shared and divergent host responses, as well as potential immune targets against HBV.

Comparison of DNA extraction methods for 16S rRNA gene sequencing in the analysis of the human gut microbiome.

The gut microbiome is widely analyzed using high-throughput sequencing, such as 16S rRNA gene amplicon sequencing and shotgun metagenomic sequencing (SMS). DNA extraction is known to have a large impact on the metagenomic analyses. The aim of this study was to compare DNA extraction protocols for 16S sequencing. In that context, four commonly used DNA extraction methods were compared for the analysis of the gut microbiota. Commercial versions were evaluated against modified protocols using a stool preprocessing device (SPD, bioMérieux) upstream DNA extraction. Stool samples from nine healthy volunteers and nine patients with a Clostridium difficile infection were extracted with all protocols and 16S sequenced. Protocols were ranked using wet- and dry-lab criteria, including quality controls of the extracted genomic DNA, alpha-diversity, accuracy using a mock community of known composition and repeatability across technical replicates. SPD improved overall efficiency of three of the four tested protocols compared with their commercial version, in terms of DNA extraction yield, sample alpha-diversity, and recovery of Gram-positive bacteria. The best overall performance was obtained for the S-DQ protocol, SPD combined with the DNeasy PowerLyser PowerSoil protocol from QIAGEN. Based on this evaluation, we strongly believe that the use of such stool preprocessing device improves both the standardization and the quality of the DNA extraction in the human gut microbiome studies.

Immune Profiling Demonstrates a Common Immune Signature of Delayed Acquired Immunodeficiency in Patients With Various Etiologies of Severe Injury.

OBJECTIVES: The host response plays a central role in the pathophysiology of sepsis and severe injuries. So far, no study has comprehensively described the overtime changes of the injury-induced immune profile in a large cohort of critically ill patients with different etiologies. DESIGN: Prospective observational cohort study. SETTING: Adult ICU in a University Hospital in Lyon, France. PATIENTS: Three hundred fifty-three septic, trauma, and surgical patients and 175 healthy volunteers were included in the REAnimation Low Immune Status Marker study. INTERVENTIONS: None. MEASUREMENTS AND MAIN RESULTS: Extensive immune profiling was performed by assessing cellular phenotypes and functions, protein, and messenger RNA levels at days 1-2, 3-4, and 5-7 after inclusion using a panel of 30 standardized immune markers. Using this immunomonitoring panel, no specificity in the immune profile was observed among septic, trauma, and surgical patients. This common injury-induced immune response was characterized by an initial adaptive (i.e., physiologic) response engaging all constituents of the immune system (pro- and anti-inflammatory cytokine releases, and innate and adaptive immune responses) but not associated with increased risk of secondary infections. In contrary, the persistence in a subgroup of patients of profound immune alterations at the end of the first week after admission was associated with increased risk of secondary infections independently of exposure to invasive devices. The combined monitoring of markers of pro-/anti-inflammatory, innate, and adaptive immune responses allowed a better enrichment of patients with risk of secondary infections in the selected population. CONCLUSIONS: Using REAnimation Low Immune Status Marker immunomonitoring panel, we detected delayed injury-acquired immunodeficiency in a subgroup of severely injured patients independently of primary disease. Critically ill patients' immune status could be captured through the combined monitoring of a common panel of complementary markers of pro-/anti-inflammatory, innate, and adaptive immune responses. Such immune monitoring needs to be incorporated in larger study cohorts with more extensive immune surveillance to develop specific hypothesis allowing for identification of biological systems affecting altered immune function related to late infection in the setting of acute systemic injury.

Herpes DNAemia and TTV Viraemia in Intensive Care Unit Critically Ill Patients: A Single-Centre Prospective Longitudinal Study.

Introduction: We analysed blood DNAemia of TTV and four herpesviruses (CMV, EBV, HHV6, and HSV-1) in the REAnimation Low Immune Status Marker (REALISM) cohort of critically ill patients who had presented with either sepsis, burns, severe trauma, or major surgery. The aim was to identify common features related to virus and injury-associated pathologies and specific features linking one or several viruses to a particular pathological context. Methods: Overall and individual viral DNAemia were measured over a month using quantitative PCR assays from the 377 patients in the REALISM cohort. These patients were characterised by clinical outcomes [severity scores, mortality, Intensive Care Unit (ICU)-acquired infection (IAI)] and 48 parameters defining their host response after injury (cell populations, immune functional assays, and biomarkers). Association between viraemic event and clinical outcomes or immune markers was assessed using χ2-test or exact Fisher's test for qualitative variables and Wilcoxon test for continuous variables. Results: The cumulative incidence of viral DNAemia increased from below 4% at ICU admission to 35% for each herpesvirus during the first month. EBV, HSV1, HHV6, and CMV were detected in 18%, 12%, 10%, and 9% of patients, respectively. The incidence of high TTV viraemia (>10,000 copies/ml) increased from 11% to 15% during the same period. Herpesvirus viraemia was associated with severity at admission; CMV and HHV6 viraemia correlated with mortality during the first week and over the month. The presence of individual herpesvirus during the first month was significantly associated (p < 0.001) with the occurrence of IAI, whilst herpesvirus DNAemia coupled with high TTV viraemia during the very first week was associated with IAI. Herpesvirus viraemia was associated with a lasting exacerbated host immune response, with concurrent profound immune suppression and hyper inflammation, and delayed return to immune homeostasis. The percentage of patients presenting with herpesvirus DNAemia was significantly higher in sepsis than in all other groups. Primary infection in the hospital and high IL10 levels might favour EBV and CMV reactivation. Conclusion: In this cohort of ICU patients, phenotypic differences were observed between TTV and herpesviruses DNAemia. The higher prevalence of herpesvirus DNAemia in sepsis hints at further studies that may enable a better in vivo understanding of host determinants of herpesvirus viral reactivation. Furthermore, our data suggest that EBV and TTV may be useful as additional markers to predict clinical deterioration in ICU patients.

Blood transcriptomics reveal the evolution and resolution of the immune response in tuberculosis.

Blood transcriptomics have revealed major characteristics of the immune response in active TB, but the signature early after infection is unknown. In a unique clinically and temporally well-defined cohort of household contacts of active TB patients that progressed to TB, we define minimal changes in gene expression in incipient TB increasing in subclinical and clinical TB. While increasing with time, changes in gene expression were highest at 30 d before diagnosis, with heterogeneity in the response in household TB contacts and in a published cohort of TB progressors as they progressed to TB, at a bulk cohort level and in individual progressors. Blood signatures from patients before and during anti-TB treatment robustly monitored the treatment response distinguishing early and late responders. Blood transcriptomics thus reveal the evolution and resolution of the immune response in TB, which may help in clinical management of the disease.

An integrated transcriptomics and metabolomics study of the immune response of newly hatched chicks to the cytosine-phosphate-guanine oligonucleotide stimulation.

The immunological immaturity of the innate immune system during the first-week post-hatch enables pathogens to infect chickens, leading to the death of the animals. Current preventive solutions to improve the resistance of chicks to infections include vaccination, breeding, and sanitation. Other prophylactic solutions have been investigated, such as the stimulation of animal health with immunostimulants. Recent studies showed that administration of immune-modulators to one-day-old chicks, or in ovo, significantly reduces mortality in experimental bacterial or viral infection challenge models. Owing to a lack of molecular biomarkers required to evaluate chicken immune responses and assess the efficacy of vaccines or immune-modulators, challenge models are still used. One way to reduce challenge experiments is to define molecular signatures through omics approaches, resulting in new methodologies to rapidly screen candidate molecules or vaccines. This study aims at identifying a dual transcriptomics and metabolomics blood signature after administration of CpG-ODN (cytosine-phosphate-guanine oligodeoxynucleotides), a reference immune-stimulatory molecule. A clinical study was conducted with chicks and transcriptomics and metabolomics analyses were performed on whole-blood and plasma samples, respectively. Differentially expressed genes and metabolites with different abundance were identified in chicks treated with CpG-ODN. The results showed that CpG-ODN activated the innate immune system, within hours after administration, and its effect lasted over time, as metabolomics and transcriptomics profiles still varied 6 D after administration. In conclusion, through an integrated clinical omics approach, we deciphered in part the mode of action of CpG-ODN in post-hatch chicks.

Caecal microbiota compositions from 7-day-old chicks reared in high-performance and low-performance industrial farms and systematic culturomics to select strains with anti-Campylobacter activity.

There is growing interest in exploring the chickens' intestinal microbiota and understanding its interactions with the host. The objective is to optimize this parameter in order to increase the productivity of farm animals. With the goal to isolate candidate probiotic strains, specific culturomic methods were used in our study to culture commensal bacteria from 7-days old chicks raised in two farms presenting long history of high performance. A total of 347 isolates were cultured, corresponding to at least 64 species. Among the isolates affiliated to the Firmicutes, 26 had less than 97% identity of their partial 16S sequence with that of the closest described species, while one presented less than 93% identity, thus revealing a significant potential for new species in this ecosystem. In parallel, and in order to better understand the differences between the microbiota of high-performing and low-performing animals, caecal contents of animals collected from these two farms and from a third farm with long history of low performance were collected and sequenced. This compositional analysis revealed an enrichment of Faecalibacterium-and Campylobacter-related sequences in lower-performing animals whereas there was a higher abundance of enterobacteria-related sequences in high-performing animals. We then investigated antibiosis activity against C. jejuni ATCC 700819 and C. jejuni field isolate as a first phenotypic trait to select probiotic candidates. Antibiosis was found to be limited to a few strains, including several lactic acid bacteria, a strain of Bacillus horneckiae and a strain of Escherichia coli. The antagonist activity depended on test conditions that mimicked the evolution of the intestinal environment of the chicken during its lifetime, i.e. temperature (37°C or 42°C) and oxygen levels (aerobic or anaerobic conditions). This should be taken into account according to the stage of development of the animal at which administration of the active strain is envisaged.

Evaluation of integrative clustering methods for the analysis of multi-omics data.

Recent advances in sequencing, mass spectrometry and cytometry technologies have enabled researchers to collect large-scale omics data from the same set of biological samples. The joint analysis of multiple omics offers the opportunity to uncover coordinated cellular processes acting across different omic layers. In this work, we present a thorough comparison of a selection of recent integrative clustering approaches, including Bayesian (BCC and MDI) and matrix factorization approaches (iCluster, moCluster, JIVE and iNMF). Based on simulations, the methods were evaluated on their sensitivity and their ability to recover both the correct number of clusters and the simulated clustering at the common and data-specific levels. Standard non-integrative approaches were also included to quantify the added value of integrative methods. For most matrix factorization methods and one Bayesian approach (BCC), the shared and specific structures were successfully recovered with high and moderate accuracy, respectively. An opposite behavior was observed on non-integrative approaches, i.e. high performances on specific structures only. Finally, we applied the methods on the Cancer Genome Atlas breast cancer data set to check whether results based on experimental data were consistent with those obtained in the simulations.

Metagenomic Investigation of Torque Teno Mini Virus-SH in Hematological Patients.

A new member of Anelloviridae, named torque teno mini virus (TTMV)-SH, was recently identified in the serum of three Hodgkin's lymphoma patients suggesting that TTMV-SH may be associated with this type of hematological malignancy. We investigated by metagenomic analysis the presence of TTMV-SH-related viruses in plasma samples (n = 323) collected from patients with various hematological malignancies (multiple myeloma (MM, n = 256), non-Hodgkin's lymphoma (NHL, n = 20), acute myeloid leukemia (n = 10)) and from healthy donors (n = 37). TTMV-SH-related strains were identified in 24 samples corresponding to four MM and one NHL patients. Phylogenic analysis revealed that the 24 isolates were close to the TTMV-SH strains previously identified, sharing 79.6-86.7% ORF1 nucleotide sequence identity. These results suggest that TTMV-SH-related viruses might be found in hematological diseases other than Hodgkin's lymphoma. Due to the high genetic variability within Anelloviridae species, the association between a particular medical condition and a new genotype should be interpreted with caution.

Early herpes and TTV DNAemia in septic shock patients: a pilot study.

BACKGROUND: Septic shock patients exhibit an increased incidence of viral reactivation. Precise timing of such reactivation-as an early marker of immune suppression, or as a consequence of the later-is not known precisely. Here, using a fully designed nucleic acid extraction automated procedure together with tailored commercial PCR kits, we focused on the description of early reactivation within the first week of ICU admission of several herpes viruses and Torque Teno virus (TTV) in 98 septic shock patients. RESULTS: Most of septic shock patients had at least one viremia event during the first week (88%). TTV and herpesviruses were detected in 56% and 53% of septic shock patient, respectively. The two most frequent herpesviruses detected within the first week were EBV (35%) and HSV1 (26%). Different kinetic were observed among herpesviruses, faster for EBV and HSV1 than for CMV and HHV6. Although no association was found between herpes viremia and secondary infections, patients with herpesviridae-related viremia were more severe, e.g., higher SOFA scores and plasma lactate levels. While reactivating only 1 virus was not associated with mortality, patients with multiple viremia events had higher ICU mortality. Surprisingly, EBV + TTV early reactivation seemed associated with a lower D28 mortality. No clear association was observed between viremia and immune biomarkers. CONCLUSION: Applying a semi-automated process of viral DNAemia determination to this cohort of 98 patients with septic shock, we observed that the number of patients with positive viremia increased during the first week in the ICU. Of note, there was no improvement in predicting the outcome when using viremia status. Nevertheless, this pilot study, introducing standardized procedures from extraction to detection, provides the basis for future standardized diagnostic criteria. A prospective longitudinal clinical study using these procedures will enable determination of whether such viremia is due to a lack of a latent virus control by the immune system or a true clinical viral infection.

Metagenomic Next-Generation Sequencing Reveals Individual Composition and Dynamics of Anelloviruses during Autologous Stem Cell Transplant Recipient Management.

Over recent years, there has been increasing interest in the use of the anelloviruses, the major component of the human virome, for the prediction of post-transplant complications such as severe infections. Due to an important diversity, the comprehensive characterization of this viral family over time has been poorly studied. To overcome this challenge, we used a metagenomic next-generation sequencing (mNGS) approach with the aim of determining the individual anellovirus profile of autologous stem cell transplant (ASCT) patients. We conducted a prospective pilot study on a homogeneous patient cohort regarding the chemotherapy regimens that included 10 ASCT recipients. A validated viral mNGS workflow was used on 108 plasma samples collected at 11 time points from diagnosis to 90 days post-transplantation. A complex interindividual variability in terms of abundance and composition was noticed. In particular, a strong sex effect was found and confirmed using quantitative PCR targeting torque teno virus, the most abundant anellovirus. Interestingly, an important turnover in the anellovirus composition was observed during the course of the disease revealing a strong intra-individual variability. Although more studies are needed to better understand anellovirus dynamics, these findings are of prime importance for their future use as biomarkers of immune competence.

Differential proteomic analysis of virus-enriched fractions obtained from plasma pools of patients with dengue fever or severe dengue.

BACKGROUND: Dengue is the most widespread mosquito-borne viral disease of public health concern. In some patients, endothelial cell and platelet dysfunction lead to life-threatening hemorrhagic dengue fever or dengue shock syndrome. Prognostication of disease severity is urgently required to improve patient management. The pathogenesis of severe dengue has not been fully elucidated, and the role of host proteins associated with viral particles has received little exploration. METHODS: The proteomes of virion-enriched fractions purified from plasma pools of patients with dengue fever or severe dengue were compared. Virions were purified by ultracentrifugation combined with a water-insoluble polyelectrolyte-based technique. Following in-gel hydrolysis, peptides were analyzed by nano-liquid chromatography coupled to ion trap mass spectrometry and identified using data libraries. RESULTS: Both dengue fever and severe dengue viral-enriched fractions contained identifiable viral envelope proteins and host cellular proteins. Canonical pathway analysis revealed the identified host proteins are mainly involved in the coagulation cascade, complement pathway or acute phase response signaling pathway. Some host proteins were over- or under-represented in plasma from patients with severe dengue compared to patients with dengue fever. ELISAs were used to validate differential expression of a selection of identified host proteins in individual plasma samples of patients with dengue fever compared to patients with severe dengue. Among 22 host proteins tested, two could differentiate between dengue fever and severe dengue in two independent cohorts (olfactomedin-4: area under the curve (AUC), 0.958; and platelet factor-4: AUC, 0.836). CONCLUSION: A novel technique of virion-enrichment from plasma has allowed to identify two host proteins that have prognostic value for classifying patients with acute dengue who are more likely to develop a severe dengue. The impact of these host proteins on pathogenicity and disease outcome are discussed.

Gene expression profiles in alveolar macrophages induced by lipopolysaccharide in humans.

Lipopolysaccharide (LPS) is ubiquitous in the environment. Inhalation of LPS has been implicated in the pathogenesis and/or severity of several lung diseases, including pneumonia, chronic obstructive pulmonary disease and asthma. Alveolar macrophages are the main resident leukocytes exposed to inhaled antigens. To obtain insight into which innate immune pathways become activated within human alveolar macrophages upon exposure to LPS in vivo, we conducted a study in eight healthy humans, in which we instilled sterile saline into a lung segment by bronchoscope, followed by instillation of LPS into the contralateral lung. Six hours later, a bilateral bronchoalveolar lavage was performed and whole-genome transcriptional profiling was done on purified alveolar macrophages, comparing cells exposed to saline or LPS from the same individuals. LPS induced differential expression of 2,932 genes in alveolar macrophages; 1,520 genes were upregulated, whereas 1,440 genes were downregulated. A total of 26 biological functions were overrepresented in LPS-exposed macrophages; 44 canonical pathways affected by LPS were identified, among which the genes associated with the role of pattern recognition receptors in recognition of bacteria and viruses represented the top pathway. Other pathways included cellular immune response, signaling by tumor necrosis factor (receptor) family members, cytokine signaling and glucocorticoid receptor signaling. These results reveal for the first time a large number of functional pathways influenced by the biologically relevant challenge provided by LPS administered into the airways. These data can assist in identifying novel targets for therapeutic intervention in pulmonary diseases associated with LPS exposure, including pneumonia, asthma and chronic obstructive pulmonary disease.

Importance of correlation between gene expression levels: application to the type I interferon signature in rheumatoid arthritis.

BACKGROUND: The analysis of gene expression data shows that many genes display similarity in their expression profiles suggesting some co-regulation. Here, we investigated the co-expression patterns in gene expression data and proposed a correlation-based research method to stratify individuals. METHODOLOGY/PRINCIPAL FINDINGS: Using blood from rheumatoid arthritis (RA) patients, we investigated the gene expression profiles from whole blood using Affymetrix microarray technology. Co-expressed genes were analyzed by a biclustering method, followed by gene ontology analysis of the relevant biclusters. Taking the type I interferon (IFN) pathway as an example, a classification algorithm was developed from the 102 RA patients and extended to 10 systemic lupus erythematosus (SLE) patients and 100 healthy volunteers to further characterize individuals. We developed a correlation-based algorithm referred to as Classification Algorithm Based on a Biological Signature (CABS), an alternative to other approaches focused specifically on the expression levels. This algorithm applied to the expression of 35 IFN-related genes showed that the IFN signature presented a heterogeneous expression between RA, SLE and healthy controls which could reflect the level of global IFN signature activation. Moreover, the monitoring of the IFN-related genes during the anti-TNF treatment identified changes in type I IFN gene activity induced in RA patients. CONCLUSIONS: In conclusion, we have proposed an original method to analyze genes sharing an expression pattern and a biological function showing that the activation levels of a biological signature could be characterized by its overall state of correlation.

Lack of association of IL-10 promoter gene variants with type 1 diabetes in a French population.

Although genetic predisposition to type 1 diabetes shows a strong association with human leukocyte antigen (HLA) class II alleles, additional genes may influence the immune process and the progression of beta cell loss. Preliminary reports suggested that IL-10 gene polymorphisms contribute to susceptibility to type 1 diabetes. We analyzed the frequencies of three main variants of the promoter region of the IL-10 gene at the positions -1082, -819, and -592 in a cohort of 358 type 1 diabetic patients representing the same regional population pattern and 519 controls from the same region using an enzyme-linked oligonucleotide sorbent assay. We did not find any statistical association in the entire cohort or after stratification for high-risk HLA-DQ alleles. However, the IL-10 -1082 polymorphism was significantly associated with GAD and IA-2 antibodies at clinical onset. Such polymorphism is known to be associated with the reduction of secreted IL-10 which may support the concept of accelerated Th-1 T-cell reactivity. In conclusion, IL-10 promoter gene variants may contribute, but to a minor extent, to disease susceptibility in juvenile type 1 diabetes and should not be included in the routine genetic screening of high-risk individuals.

Morphological convergence in Hippidion and Equus (Amerhippus) South American equids elucidated by ancient DNA analysis.

Unusual equids named hippidions inhabited South America for more than 2 MY (million years). Like many other animals they succumbed to the worldwide climatic change that occurred 10 KY (thousand years) ago and completely disappeared during the great late Pleistocene megafaunal extinction. According to fossil records and numerous dental, cranial, and postcranial characters, Hippidion and Equus lineages are known to have diverged prior to 10 MY. Some equid bones from Rio Verde and Ultima Esperanza (Patagonia, Chile) dating back to the late Pleistocene period (8-13 KY) have been identified as Hippidion saldiasi, while a few teeth have been assigned to Equus. Six samples of those remains have been obtained from the Zoological Museum of Amsterdam for ancient DNA analysis to try to place Hippidion in the evolutive tree of Perissodactyla. Two samples of Hippidion and one sample of Equus yielded 241-394 bp of the mtDNA control region and 172-296 bp of the cytochrome b gene. Unexpectedly, all the sequences clustered deep inside the Equus genus, casting doubt on the initial identification of the bones. For paleontologists, one of the striking and classical diagnostic characters of Hippidion is their extremely short and massive metapodials, a probable locomotory adaptation to the Andine steep slopes. However, our DNA analysis reveals that a very Hippidion-like metapod might also have been possessed by another South American equid, i.e., Equus (Amerhippus), an interpretation supported by complementary anatomical observations. This adaptive convergence between members of the two South American equid genera may lead paleontologists to limb bone misidentification.