Host Transcriptome Research Articles

Heme Metabolism Mediates the Effects of Smoking on Gut Microbiome.

The number of smokers worldwide increased greatly during the past decades and reached 1.14 billion in 2019, becoming a leading risk factor for human health. Tobacco smoking has wide effects on human genetics, epigenetics, transcriptome, and gut microbiome. Although many studies have revealed effects of smoking on host transcriptome, research on the relationship between smoking, host gene expression, and the gut microbiome is limited. We first explored transcriptome and metagenome profile differences between smokers and nonsmokers. To evaluate the relationship between host gene expression and gut microbiome, we then applied bidirectional mediation analysis to infer causal relationships between smoking, gene expression, and gut microbes. Metagenome and transcriptome analyses revealed 71 differential species and 324 differential expressed genes between smokers and nonsmokers. With smoking as an exposure variable, we identified 272 significant causal relationships between gene expression and gut microbes, among which there were 247 genes that mediate the effect of smoking on gut microbes. Pathway-based enrichment analysis showed that these genes were significantly enriched in heme metabolic pathway, which mainly mediated the changes of Bacteroides finegoldii and Lachnospiraceae bacterium 9_1_43BFAA. Additionally, by performing metabolome data analysis in the Integrated Human Microbiome Project (iHMP) database, we verified the correlation between the intermediate products of the heme metabolism pathway (porphobilinogen, bilirubin, and biliverdin) and gut microbiome. By investigating the bidirectional interaction between smoking-related host gene expression and gut microbes, this study provided evidence for the mediation of smoking on gut microbes through co-involvement or interaction of heme metabolism. By comparing the metagenome and transcriptome sequencing profiles between 34 smokers and 33 age- and gender-matched nonsmokers, we are the first to reveal causal relationships among tobacco smoking, host gene expression, and gut microbes. These findings offer insight into how smoking affects gut microbes through host gene expression and metabolism, which highlights the importance of heme metabolism in modulating the effects of smoking on gut microbiome.

Liver microbial community and associated host transcriptome in calves with feed induced acidosis.

In the dairy industry, calves are typically fed diets rich in highly fermentable carbohydrates and low in fibrous feeds to maximize ruminal papillae and tissue development. Calves on such diets are vulnerable at developing ruminal acidosis. Prevalent in cattle, liver abscess (LA) is considered a sequela to ruminal acidosis. LAs can cause significant liver function condemnation and decreased growth and production. Currently, we know little about the liver microbiome in calves with feed-induced acidosis. Using our established model of ruminal acidosis, where young calves were fed an acidosis-inducing (AC) or -blunting (control) diet starting at birth until 17-week of age, we investigated microbial community changes in the liver resultant from ruminal acidosis. Eight calves were randomly assigned to each diet, with four animals per treatment. Rumen epithelium and liver tissues were collected at 17 weeks of age right after euthanasia. Total RNAs were extracted and followed by whole transcriptome sequencing. Microbial RNA reads were enriched bioinformatically and used for microbial taxonomy classification using Kraken2. AC Calves showed significantly less weight gain over the course of the experiment, in addition to significantly lower ruminal pH, and rumen degradation comparison to the control group (p < 0.05). In the liver, a total of 29 genera showed a significant (p < 0.05) abundance change (> 2-fold) between the treatments at 17-week of age. Among these, Fibrobacter, Treponema, Lactobacillus, and Olsenella have been reported in abscessed liver in cattle. Concurrent abundance changes in 9 of the genera were observed in both the liver and rumen tissues collected at 17-week of age, indicating potential crosstalk between the liver and rumen epithelial microbial communities. Significant association was identified between host liver gene and its embedded microbial taxa. Aside from identifying previously reported microbial taxa in cattle abscessed liver, new repertoire of actively transcribed microbial taxa was identified in this study. By employing metatranscriptome sequencing, our study painted a picture of liver microbiome in young calves with or without feed induced acidosis. Our study suggested that liver microbiome may have a critical impact on host liver physiology. Novel findings of this study emphasize the need for further in-depth analysis to uncover the functional roles of liver resident microbiome in liver metabolic acidosis resultant from feed-related ruminal acidosis.

Integrated analysis of microbiome and host transcriptome reveals the damage/protective mechanism of corn oil and olive oil on the gut health of grouper (♀ Epinephelus fuscoguttatus × ♂ E. lanceolatu)

As digestive and immune organs of animals, the gut was frequently used to evaluate the health status of aquatic animals. In previous oil source alternatives study, corn oil (CO) had been found to induce gut inflammation, while olive oil (OO) had been found to be effective in protecting intestinal health. Three diets with different oil sources (fish oil, CO, OO) were formulated for an 8-week culture experiment, and it was proposed to combine 16S sequencing and transcriptome sequencing analysis to preliminarily elucidate the damage/protection mechanism of CO and OO on the gut health of grouper (♀ Epinephelus fuscoguttatus × ♂ E. lanceolatu). We found that CO indeed damaged to gut health and destroyed the gut structure, while OO had a positive outcome in protecting the gut structure, promoting digestibility and relieving enteritis. Photobacterium, Romboutsia and Epulopiscium were significantly enriched in OO group and Staphylococcus were significantly enriched in CO group. Transcriptome sequencing further revealed CO could activated Complement and coagulation cascades, Staphylococcus aureus infection, Systemic lupus erythematosus, and Tuberculosis pathways; conversely, OO activated B-cell signaling receptors, promoted B-cell proliferation and apoptosis, and thus activated B-cell signaling pathways to enhance immunity, whereas OO can regulate IL17 signaling pathway and TNF signaling pathway to inhibit NF-κB signaling pathway to reduce pro-inflammatory response. By integrating the microbiome and transcriptome, further identified all differential microorganisms were directly and significantly correlated with differential genes, and Clostridium_sensu_stricto_1, Romboutsia, Staphylococcus might as the core regulates the expression of differential gene in the organism. These results reveal that different oil sources alter gut gene expression mainly by modulating the composition and abundance of gut microbiota, further regulating the health status of the gut. Gut microbiota could be used as biomarkers to provide reference and solutions for the mitigation of inflammation in aquatic animals.

Multiomic profiling of cutaneous leishmaniasis infections reveals microbiota-driven mechanisms underlying disease severity.

Leishmania braziliensis is a parasitic infection that can result in inflammation and skin injury with highly variable and unpredictable clinical outcomes. Here, we investigated the potential impact of microbiota on infection-induced inflammatory responses and disease resolution by conducting an integrated analysis of the skin microbiome and host transcriptome on a cohort of 62 patients infected with L. braziliensis. We found that overall bacterial burden and microbiome configurations dominated with Staphylococcus spp. were associated with delayed healing and enhanced inflammatory responses, especially by IL-1 family members. Quantification of host and bacterial transcripts on human lesions revealed that high lesional S. aureus transcript abundance was associated with delayed healing and increased expression of IL-1β. This cytokine was critical for modulating disease outcomes in L. braziliensis-infected mice colonized with S. aureus, given that its neutralization reduced pathology and inflammation. These results highlight how the human microbiome can shape disease outcomes in cutaneous leishmaniasis and suggest pathways toward host-directed therapies to mitigate the inflammatory consequences.

Temporal and spatial variations in body mass and thermogenic capacity associated with alterations in the gut microbiota and host transcriptome in mammalian herbivores

Most wild animals follow Bergmann's rule and grow in body size as cold stress increases. However, the underlying thermogenic strategies and their relationship with the gut microbiota have not been comprehensively elucidated. Herein, we used the plateau pikas as a model to investigate body mass, thermogenic capacity, host transcriptome, gut microbiota and metabolites collected from seven sites ranging from 3100 to 4700 m on the Qinghai–Tibetan Plateau (QTP) in summer and winter to test the seasonal thermogenesis strategy in small herbivorous mammals. The results showed that the increase in pika body mass with altitude followed Bergmann's rule in summer and an inverted parabolic shape was observed in winter. However, physiological parameters and transcriptome profiles indicated that the thermogenic capacity of pikas increased with altitude in summer and decreased with altitude in winter. The abundance of Firmicutes declined, whereas that of Bacteroidetes significantly increased with altitude in summer. Phenylalanine, tyrosine, and proline were enriched in summer, whereas carnitine and succinate were enriched in winter. Spearman's correlation analysis revealed significant positive correlations between Prevotella, Bacteroides, Ruminococcus, Alistipes and Akkermansia and metabolites of amino acids, pika physiological parameters, and transcriptome profiles. Moreover, metabolites of amino acids further showed significant positive correlations with pika physiological parameters and transcriptome profiles. Our study highlights that the changes in body mass and thermogenic capacity with altitude distinctly differentiate small herbivorous mammals between summer and winter on the QTP, and that the gut microbiota may regulate host thermogenesis through its metabolites.

Daily turnover of active giant virus infection during algal blooms revealed by single-cell transcriptomics.

Giant viruses infect many unicellular eukaryotes, including algae that form massive oceanic blooms. Despite the major impact of viruses on the marine ecosystem, the ability to quantify and assess active viral infection in nature remains a major challenge. We applied single-cell RNA sequencing, to profile virus and host transcriptomes of 12,000 single algal cells from a coccolithophore bloom. Viral infection was detected already at early exponential bloom phase, negatively correlating with the bloom intensity. A consistent percent of infected coccolithophores displayed the early phase of viral replication for several consecutive days, indicating a daily turnover and continuous virocell-associated metabolite production, potentially affecting the surrounding microbiome. Linking single-cell infection state to host physiology revealed that infected cells remained calcified even in the late infection stage. These findings stress the importance of studying host-virus dynamics in natural populations, at single-cell resolution, to better understand virus life cycle and its impact on microbial food webs.

Mpox virus infection and drug treatment modelled in human skin organoids.

Mpox virus (MPXV) primarily infects human skin to cause lesions. Currently, robust models that recapitulate skin infection by MPXV are lacking. Here we demonstrate that human induced pluripotent stem cell-derived skin organoids are susceptible to MPXV infection and support infectious virus production. Keratinocytes, the predominant cell type of the skin epithelium, effectively support MPXV infection. Using transmission electron microscopy, we visualized the four stages of intracellular virus particle assembly: crescent formation, immature virions, mature virions and wrapped virions. Transcriptional analysis showed that MPXV infection rewires the host transcriptome and triggers abundant expression of viral transcripts. Early treatment with the antiviral drug tecovirimat effectively inhibits infectious virus production and prevents host transcriptome rewiring. Delayed treatment with tecovirimat also inhibits infectious MPXV particle production, albeit to a lesser extent. This study establishes human skin organoids as a robust experimental model for studying MPXV infection, mapping virus-host interactions and testing therapeutics.

How the luxR Gene Affects the Pathogenicity of Pseudomonas plecoglossicida and the Immune Response of Epinephelus coioides

This study aimed to investigate the effect of reduced expression of the luxR gene on the virulence of Pseudomonas plecoglossicida and the immune response of Epinephelus coioides. To achieve this, RNA interference (RNAi) was used to silence the luxR gene, and the pathogenicity of wild-type and luxR-RNAi strains of P. plecoglossicida, as well as the immune response of Epinephelus coioides to the infection of these two strains, were compared. The mutant strain with the highest silencing efficiency of 70.1% was selected for subsequent analysis. Silencing the luxR gene in the mutant strain resulted in a significant 30% reduction in mortality rates in artificially infected Epinephelus coioides compared to the wild-type strain. Transcriptome analysis revealed that the host transcriptome, particularly in the spleens of infected Epinephelus coioides, was markedly altered by the silencing of the luxR gene in the mutant strain. Tilapia infected with the luxR-RNAi strain exhibited altered immune defenses, with changes in gene expression primarily in the NOD-like receptor (NLR) signaling pathway. These results suggest that the luxR gene plays a crucial role in the host’s resistance to pathogen invasion, and reducing its expression could decrease quorum sensing (QS) signals while increasing the expression of the IL-1β gene in the host’s NLR pathway. This effect may lead to a pro-inflammatory response that enhances the immune response to infection. Further investigation of these mechanisms may lead to innovative approaches to treating bacterial infections.

The specific viral composition in triple-negative breast cancer tissue shapes the specific tumor microenvironment characterized on pathological images

Numerous studies have shown that different subtypes of breast cancer (BC) have significant differences in terms of the tumor microbiome, host gene expression, and histopathological image, whereas the biological links between these cancer-associated indicators are still unknown. Here, we performed a comprehensive analysis with 610 patients of the four subtypes of BC with matched tissue microbiota, host transcriptome, and histopathological image samples. Correlation analysis showed that the composition of intratumoral viruses shaped the tumor microenvironment (TME) of patients with BC, and the TME was further reflected in the histopathological images. Of the four subtypes, patients with triple-negative breast cancer (TNBC) had unique intratumoral viral community composition, non-cancer cell infiltration in the TME, and histopathological image characteristics. Furthermore, we detected multiple virus-cell-image association axes in TNBC, in which tumor-associated macrophages (TAMs) have clinical prognostic implication. This study provides a comprehensive map of the associations between the intratumoral virome, TME, and histopathological image of TNBC, as well as insights into disease prognosis that can be crucial for precise therapeutic intervention strategies.

Characterizing the influence of various antimicrobials used for metaphylaxis against bovine respiratory disease on host transcriptome responses.

Currently, control against bovine respiratory disease (BRD) primarily consists of mass administration of an antimicrobial upon arrival to facility, termed "metaphylaxis." The objective of this study was to determine the influence of six different antimicrobials used as metaphylaxis on the whole blood host transcriptome in healthy steers upon and following arrival to the feedlot. One hundred and five steers were stratified by arrival body weight (BW = 247 ± 28 kg) and randomly and equally allocated to one of seven treatments: negative control (NC), ceftiofur (CEFT), enrofloxacin (ENRO), florfenicol (FLOR), oxytetracycline (OXYT), tildipirosin (TILD), or tulathromycin (TULA). On day 0, whole blood samples and BW were collected prior to a one-time administration of the assigned antimicrobial. Blood samples were collected again on days 3, 7, 14, 21, and 56. A subset of cattle (n = 6) per treatment group were selected randomly for RNA sequencing across all time points. Isolated RNA was sequenced (NovaSeq 6,000; ~35 M paired-end reads/sample), where sequenced reads were processed with ARS-UCD1.3 reference-guided assembly (HISAT2/StringTie2). Differential expression analysis comparing treatment groups to NC was performed with glmmSeq (FDR ≤ 0.05) and edgeR (FDR ≤ 0.1). Functional enrichment was performed with KOBAS-i (FDR ≤ 0.05). When compared only to NC, unique differentially expressed genes (DEGs) found within both edgeR and glmmSeq were identified for CEFT (n = 526), ENRO (n = 340), FLOR (n = 56), OXYT (n = 111), TILD (n = 3,001), and TULA (n = 87). At day 3, CEFT, TILD, and OXYT shared multiple functional enrichment pathways related to T-cell receptor signaling and FcεRI-mediated NF-kappa beta (kB) activation. On day 7, Class I major histocompatibility complex (MHC)-mediated antigen presentation pathways were enriched in ENRO and CEFT groups, and CEFT and FLOR had DEGs that affected IL-17 signaling pathways. There were no shared pathways or Gene Ontology (GO) terms among treatments at day 14, but TULA had 19 pathways and eight GO terms enriched related to NF- κβ activation, and interleukin/interferon signaling. Pathways related to cytokine signaling were enriched by TILD on day 21. Our research demonstrates immunomodulation and potential secondary therapeutic mechanisms induced by antimicrobials commonly used for metaphylaxis, providing insight into the beneficial anti-inflammatory properties antimicrobials possess.

INVADEseq to identify cell-adherent or invasive bacteria and the associated host transcriptome at single-cell-level resolution.

Single-cell RNA sequencing (scRNAseq) technologies have been beneficial in revealing and describing cellular heterogeneity within mammalian tissues, including solid tumors. However, many of these techniques apply poly(A) selection of RNA, and thus have primarily focused on determining the gene signatures of eukaryotic cellular components of the tumor microenvironment. Microbiome analysis has revealed the presence of microbial ecosystems, including bacteria and fungi, within human tumor tissues from major cancer types. Imaging data have revealed that intratumoral bacteria may be located within epithelial and immune cell types. However, as bacterial RNA typically lacks a poly(A) tail, standard scRNAseq approaches have limited ability to capture this microbial component of the tumor microenvironment. To overcome this, we describe the invasion-adhesion-directed expression sequencing (INVADEseq) approach, whereby we adapt 10x Genomics 5' scRNAseq protocol by introducing a primer that targets a conserved region of the bacterial 16S ribosomal RNA gene in addition to the standard primer for eukaryotic poly(A) RNA selection. This 'add-on' approach enables the generation of eukaryotic and bacterial DNA libraries at eukaryotic single-cell level resolution, utilizing the 10x barcode to identify single cells with intracellular bacteria. The INVADEseq method takes 30 h to complete, including tissue processing, sequencing and computational analysis. As an output, INVADEseq has shown to be a reliable tool in human cancer cell lines and patient tumor specimens by detecting the proportion of human cells that harbor bacteria and the identities of human cells and intracellular bacteria, along with identifying host transcriptional programs that are modulated on the basis of associated bacteria.

The impact of Parkinson's disease-associated gut microbiota on the transcriptome in Drosophila.

Parkinson's disease (PD) is a common neurodegenerative disease in middle-aged and elderly people, and many studies have confirmed that the disorder of gut microbiota is involved in the pathophysiological process of PD. However, the molecular mechanism of gut microbiota in regulating the pathogenesis of PD is still lacking. In this study, to investigate the impact of PD-associated gut microbiota on host transcriptome, we established various PD models with fecal microbiota transplantation (FMT) in the model organism Drosophila followed by integrative data analysis of microbiome and transcriptome. We first constructed rotenone-induced PD models in Drosophila followed by FMT in different groups. Microbial analysis by 16S rDNA sequencing showed that gut microbiota from PD Drosophila could affect bacterial structure of normal Drosophila, and gut microbiota from normal Drosophila could affect bacterial structure of PD Drosophila. Transcriptome analysis revealed that PD-associated gut microbiota influenced expression patterns of genes enriched in neuroactive ligand-receptor interaction, lysosome, and diverse metabolic pathways. Importantly, to verify our findings, we transplanted Drosophila with fecal samples from clinical PD patients. Compared to the control, Drosophila transplanted with fecal samples from PD patients had reduced microbiota Acetobacter and Lactobacillus, and differentially expressed genes enriched in diverse metabolic pathways. In summary, our results reveal the influence of PD-associated gut microbiota on host gene expression, and this study can help better understand the link between gut microbiota and PD pathogenesis through gut-brain axis. IMPORTANCE Gut microbiota plays important roles in regulating host gene expression and physiology through complex mechanisms. Recently, it has been suggested that disorder of gut microbiota is involved in the pathophysiological process of Parkinson's disease (PD). However, the molecular mechanism of gut microbiota in regulating the pathogenesis of PD is still lacking. In this study, to investigate the impact of PD-associated gut microbiota on host transcriptome, we established various PD models with fecal microbiota transplantation in the model organism Drosophila followed by integrative data analysis of microbiome and transcriptome. We also verified our findings by transplanting Drosophila with fecal samples from clinical PD patients. Our results demonstrated that PD-associated gut microbiota can induce differentially expressed genes enriched in diverse metabolic pathways. This study can help better understand the link between gut microbiota and PD pathogenesis through gut-brain axis.

Multi-breed host rumen epithelium transcriptome and microbiome associations and their relationship with beef cattle feed efficiency

Understanding host-microbial interactions in the rumen and its influence on desirable production traits may lead to potential microbiota manipulation or genetic selection for improved cattle feed efficiency. This study investigated the host transcriptome and its correlation with the rumen archaea and bacteria differential abundance of two pure beef cattle breeds (Angus and Charolais) and one composite beef hybrid (Kinsella) divergent for residual feed intake (RFI; low-RFI vs. high-RFI). Using RNA-Sequencing of rumen tissue and 16S rRNA gene amplicon sequencing, differentially expressed genes (FDR ≤ 0.05, |log2(Fold-change) >|2) and differentially abundant (p-value < 0.05) archaea and bacteria amplicon sequence variants (ASV) were determined. Significant correlations between gene expression and ASVs (p-value < 0.05) were determine using Spearman correlation. Interesting associations with muscle contraction and the modulation of the immune system were observed for the genes correlated with bacterial ASVs. Potential functional candidate genes for feed efficiency status were identified for Angus (CCL17, CCR3, and CXCL10), Charolais (KCNK9, GGT1 and IL6), and Kinsella breed (ESR2). The results obtained here provide more insights regarding the applicability of target host and rumen microbial traits for the selection and breeding of more feed efficient beef cattle.

Abstract PO-100: Fusobacterium nucleatum promotes epithelial-mesenchymal transition (EMT) of oral squamous cell carcinoma (OSCC) via INHBA-dependent SMAD/Laminin332/PI3K/AKT pathway

Abstract Introduction The enrichment of Fusobacterium nucleatum has been reported in both oral rinses and lesion samples from oral squamous cell carcinoma (OSCC) patients based on several case-control sequencing analyses, but studies on the potential role of F. nucleatum in OSCC tumorigenesis are still limited. In this study, we sought to elucidate the relationship between F. nucleatum and host transcriptome dysregulation in potentially carcinogenic pathways. Method The association between microbiota dysbiosis and host gene transcriptome was profiled by bacterial 16S rRNA V3-V4 amplicon sequencing and host gene long non-coding RNA sequencing (lncRNA-seq) for tumor and adjacent normal (AN) tissues from OSCC patients. Cell migration and invasion assays were performed using the in vitro bacteria-cell co-culture model. EMT-related molecular mechanisms induced by F. nucleatum were further validated by IncRNA-seq, real-time PCR, and western blot for bacteria-infected cells, respectively, at the genetic and proteomic levels. Result Significant enrichment of Fusobacterium was observed in OSCC tumor samples when compared to the paired AN samples. Correlation analysis between mucosal oral bacteria and host transcriptome in the OSCC tumor microenvironment elucidated the positive relationship between Fusobacterium and EMT pathways including significantly upregulated genes INHBA, SNAI2, and LAMC2. Using in-vitro cell models, we further observed that F. nucleatum significantly promoted cell migration and invasion, and significantly dysregulated the expression of several EMT factors including the increased inducer inhibin beta A (INHBA), core regulators SNAI1 (snai1) and SLUG (snai2), effectors vimentin (VIM) and N-cadherin (CDH2), and decreased E-cadherin (CDH1). Following the INHBA/SMAD pathway induced by F. nucleatum infection, the expression of an essential extracellular matrix glycoprotein, Laminin 332, was dramatically enhanced, which may trigger the phosphorylation of phosphatidylinositol 3-kinase and AKT kinase in the downstream PI3K/AKT pathway for the promoted cell migration and invasion ability. Conclusion Our study confirms the enrichment of Fusobacterium in the OSCC tumor microenvironment and reveals that F. nucleatum infection induces the INHBA/SAMD pathway and subsequently increases Laminin 332 expression, which may further promote the phosphorylation of PI3K/AKT and enhance OSCC migration and invasion. Citation Format: Hengyan Zhu, Liuyang Cai, Zigui Chen, Jason Y. K. Chan. Fusobacterium nucleatum promotes epithelial-mesenchymal transition (EMT) of oral squamous cell carcinoma (OSCC) via INHBA-dependent SMAD/Laminin332/PI3K/AKT pathway [abstract]. In: Proceedings of the AACR-AHNS Head and Neck Cancer Conference: Innovating through Basic, Clinical, and Translational Research; 2023 Jul 7-8; Montreal, QC, Canada. Philadelphia (PA): AACR; Clin Cancer Res 2023;29(18_Suppl):Abstract nr PO-100.

On a path toward a broad-spectrum anti-viral: inhibition of HIV-1 and coronavirus replication by SR kinase inhibitor harmine.

This study highlights the crucial role RNA processing plays in regulating viral gene expression and replication. By targeting SR kinases, we identified harmine as a potent inhibitor of HIV-1 as well as coronavirus (HCoV-229E and multiple SARS-CoV-2 variants) replication. Harmine inhibits HIV-1 protein expression and reduces accumulation of HIV-1 RNAs in both cell lines and primary CD4+ T cells. Harmine also suppresses coronavirus replication post-viral entry by preferentially reducing coronavirus sub-genomic RNA accumulation. By focusing on host factors rather than viral targets, our study offers a novel approach to combating viral infections that is effective against a range of unrelated viruses. Moreover, at doses required to inhibit virus replication, harmine had limited toxicity and minimal effect on the host transcriptome. These findings support the viability of targeting host cellular processes as a means of developing broad-spectrum anti-virals.

The effects of maternal care on the developmental transcriptome and metatranscriptome of a wild bee

Maternal care acts as a strong environmental stimulus that can induce phenotypic plasticity in animals and may also alter their microbial communities through development. Here, we characterize the developmental metatranscriptome of the small carpenter bee, Ceratina calcarata, across developmental stages and in the presence or absence of mothers. Maternal care had the most influence during early development, with the greatest number and magnitude of differentially expressed genes between maternal care treatments, and enrichment for transcription factors regulating immune response in motherless early larvae. Metatranscriptomic data revealed fungi to be the most abundant group in the microbiome, with Aspergillus the most abundant in early larvae raised without mothers. Finally, integrative analysis between host transcriptome and metatranscriptome highlights several fungi correlating with developmental and immunity genes. Our results provide characterizations of the influence of maternal care on gene expression and the microbiome through development in a wild bee.

Multi-omics reveal mechanisms underlying chronic kidney disease of unknown etiology (CKDu) pathogenesis using zebrafish

Chronic kidney disease of unknown etiology (CKDu) is an endemic disease in the dry zone of farming communities, Sri Lanka. The drinking water in a CKDu prevalent area contains a high concentration of F−, hardness and other environmental pollutants, including heavy metals and microcystin, which are considered possible etiology of CKDu in these areas. Here, multi-omics data with host transcriptome, metabolome and gut microbiomes were obtained using simulated local drinking water of Sri Lanka after their exposure to adult zebrafish. Based on an integrated multi-omics analysis in the context of host physiology in the kidney injury samples with different pathologic grades, two common pathways necroptosis and purine metabolism were identified as potentially important pathways that affect kidney injury. The key metabolite acetyl adenylate in the purine metabolism pathway was significantly positively correlated with Comamonas (rho = 0.72) and significantly negatively correlated with Plesiomonas (rho = −0.58). This crucial metabolite and two key gut bacteria genera may not only be potential markers but also potential therapeutic targets in the uric acid metabolic pathway, which is an important factor in the pathogenesis of acute kidney injury (AKI) in general, as well as of chronic kidney disease (CKD). Based on this, we revealed the urea metabolism pathway of kidney injury in zebrafish and provided a new avenue for the treatment of CKDu in Sri Lanka.

Systems analysis reveals differential expression of endocervical genes in African women randomized to DMPA-IM, LNG implant or cu-IUD

Although effective contraceptives are crucial for preventing unintended pregnancies, evidence suggests that their use may perturb the female genital tract (FGT). A comparative analysis of the effects of the most common contraceptives on the FGT have not been evaluated in a randomized clinical trial setting. Here, we evaluated the effect of three long-acting contraceptive methods: depot medroxyprogesterone acetate(DMPA-IM), levonorgestrel(LNG) implant, and a copper intrauterine device (Cu-IUD), on the endocervical host transcriptome in 188 women from the Evidence for Contraceptive Options and HIV Outcomes Trial (ECHO) trial. Cu-IUD usage showed the most extensive transcriptomic changes, and was associated with inflammatory and anti-viral host responses. DMPA-IM usage was enriched for pathways associated with T cell responses. LNG implant had the mildest effect on endocervical gene expression, and was associated with growth factor signaling. These data provide a mechanistic basis for the diverse influence that varying contraceptives have on the FGT.

Multi-omics analysis revealing the interplay between gut microbiome and the host following opioid use

ABSTRACT Opioid crisis is an ongoing epidemic since the past several decades in the United States. Opioid use-associated microbial dysbiosis is emerging as a key regulator of intestinal homeostasis and behavioral responses to opioid. However, the mechanistic insight into the role of microbial community in modulating host response is unavailable. To uncover the role of opioid-induced dysbiosis in disrupting intestinal homeostasis we utilized whole genome sequencing, untargeted metabolomics, and mRNA sequencing to identify changes in microbiome, metabolome, and host transcriptome respectively. Morphine treatment resulted in significant expansion of Parasuterella excrementihominis, Burkholderiales bacterium 1_1_47, Enterococcus faecalis, Enterorhabdus caecimuris and depletion of Lactobacillus johnsonii. These changes correlated with alterations in lipid metabolites and flavonoids. Significant alteration in microbial metabolism (metabolism of lipids, amino acids, vitamins and cofactors) and increased expression of virulence factors and biosynthesis of lipopolysaccharides (LPS) and lipoteichoic acid (LTA) were observed in microbiome of morphine-treated animals. In concurrence with changes in microbiome and metabolome extensive changes in innate and adaptive immune response, lipid metabolism, and gut barrier dysfunction were observed in the host transcriptome. Microbiome depleted mice displayed lower levels of inflammation, immune response and tissue destruction compared to mice harboring a dysbiotic microbiome in response to morphine treatment, thus establishing dysbiotic microbiome as mediator of morphine gut pathophysiology. Integrative analysis of multi-omics data highlighted the associations between Parasutterella excrementihominis, Burkholderiales bacterium 1_1_47, Enterococcus faecalis, Enterorhabdus caecimuris and altered levels of riboflavin, flavonoids, and lipid metabolites including phosphocholines, carnitines, bile acids, and ethanolamines with host gene expression changes involved in inflammation and barrier integrity of intestine. Omic analysis also highlighted the role of probiotic bacteria Lactobacillus johnsonii, metabolites flavonoids and riboflavin that were depleted with morphine as important factors for intestinal homeostasis. This study presents for the first time ever an interactive view of morphine-induced changes in microbial metabolism, strain level gut microbiome analysis and comprehensive view of changes in gut transcriptome. We also identified areas of potential therapeutic interventions to limit microbial dysbiosis and present a unique resource to the opioid research community.

SARS-CoV-2 infection severity and mortality is modulated by repeat-mediated regulation of alternative splicing.

Like single-stranded RNA viruses, SARS-CoV-2 hijacks the host transcriptional machinery for its own replication. Numerous traditional differential gene expression-based investigations have examined the diverse clinical symptoms caused by SARS-CoV-2 infection. The virus, on the other hand, also affects the host splicing machinery, causing host transcriptional dysregulation, which can lead to diverse clinical outcomes. Hence, in this study, we performed host transcriptome sequencing of 125 hospital-admitted COVID-19 patients to understand the transcriptomic differences between the severity sub-phenotypes of mild, moderate, severe, and mortality. We performed transcript-level differential expression analysis, investigated differential isoform usage, looked at the splicing patterns within the differentially expressed transcripts (DET), and elucidated the possible genome regulatory features. Our DTE analysis showed evidence of diminished transcript length and diversity as well as altered promoter site usage in the differentially expressed protein-coding transcripts in the COVID-19 mortality patients. We also investigated the potential mechanisms driving the alternate splicing and discovered a compelling differential enrichment of repeats in the promoter region and a specific enrichment of SINE (Alu) near the splicing sites of differentially expressed transcripts. These findings suggested a repeat-mediated plausible regulation of alternative splicing as a potential modulator of COVID-19 disease severity. In this work, we emphasize the role of scarcely elucidated functional role of alternative splicing in influencing COVID-19 disease severity sub-phenotypes, clinical outcomes, and its putative mechanism. IMPORTANCE The wide range of clinical symptoms reported during the COVID-19 pandemic inherently highlights the numerous factors that influence the progression and prognosis of SARS-CoV-2 infection. While several studies have investigated the host response and discovered immunological dysregulation during severe infection, most of them have the common theme of focusing only up to the gene level. Viruses, especially RNA viruses, are renowned for hijacking the host splicing machinery for their own proliferation, which inadvertently puts pressure on the host transcriptome, exposing another side of the host response to the pathogen challenge. Therefore, in this study, we examine host response at the transcript-level to discover a transcriptional difference that culminates in differential gene-level expression. Importantly, this study highlights diminished transcript diversity and possible regulation of transcription by differentially abundant repeat elements near the promoter region and splicing sites in COVID-19 mortality patients, which together with differentially expressed isoforms hold the potential to elaborate disease severity and outcome.