Dual RNA-seq Analysis of Patients' Cells and Viral Genome After Measles Infection.

During the infection of a host cell by an infectious agent, a series of gene expression changes occurs as a consequence of host-pathogen interactions. Unraveling this complex interplay is the key for understanding of microbial virulence and host response pathways, thus providing the basis for new molecular insights into the mechanisms of pathogenesis and the corresponding immune response. Dual RNA sequencing (dual RNA-seq) has been developed to simultaneously determine pathogen and host transcriptomes enabling both differential and coexpression analyses between the two partners as well as genome characterization in the case of RNA viruses. Here, we provide a detailed laboratory protocol and bioinformatics analysis guidelines for dual RNA-seq experiments focusing on - but not restricted to - measles virus (MeV) as a pathogen of interest. The application of dual RNA-seq technologies in MeV-infected patients can potentially provide valuable information on the structure of the viral RNA genome and on cellular innate immune responses and drive the discovery of new targets for antiviral therapy.

Nanopore-Sequencing Metabarcoding for Identification of Phytopathogenic and Endophytic Fungi in Olive (Olea europaea) Twigs.

Metabarcoding approaches for the identification of plant disease pathogens and characterization of plant microbial populations constitute a rapidly evolving research field. Fungal plant diseases are of major phytopathological concern; thus, the development of metabarcoding approaches for the detection of phytopathogenic fungi is becoming increasingly imperative in the context of plant disease prognosis. We developed a multiplex metabarcoding method for the identification of fungal phytopathogens and endophytes in olive young shoots, using the MinION sequencing platform (Oxford Nanopore Technologies). Selected fungal-specific primers were used to amplify three different genomic DNA loci (ITS, beta-tubulin, and 28S LSU) originating from olive twigs. A multiplex metabarcoding approach was initially evaluated using healthy olive twigs, and further assessed with naturally infected olive twig samples. Bioinformatic analysis of basecalled reads was carried out using MinKNOW, BLAST+ and R programming, and results were also evaluated using the BugSeq cloud platform. Data analysis highlighted the approaches based on ITS and their combination with beta-tubulin as the most informative ones according to diversity estimations. Subsequent implementation of the method on symptomatic samples identified major olive pathogens and endophytes including genera such as Cladosporium, Didymosphaeria, Paraconiothyrium, Penicillium, Phoma, Verticillium, and others.

A miR-124-mediated post-transcriptional mechanism controlling the cell fate switch of astrocytes to induced neurons.

The microRNA (miRNA) miR-124 has been employed supplementary to neurogenic transcription factors (TFs) and other miRNAs to enhance direct neurogenic conversion. The aim of this study was to investigate whether miR-124 is sufficient to drive direct reprogramming of astrocytes to induced neurons (iNs) on its own and elucidate its independent mechanism of reprogramming action. Our data show that miR-124 is a potent driver of the reprogramming switch of astrocytes toward an immature neuronal fate by directly targeting the RNA-binding protein Zfp36L1 implicated in ARE-mediated mRNA decay and subsequently derepressing Zfp36L1 neurogenic interactome. To this end, miR-124 contribution in iNs' production largely recapitulates endogenous neurogenesis pathways, being further enhanced upon addition of the neurogenic compound ISX9, which greatly improves iNs' differentiation and functional maturation. Importantly, miR-124 is potent in guiding direct conversion of reactive astrocytes to immature iNs in vivo following cortical trauma, while ISX9 supplementation confers a survival advantage to newly produced iNs.

Tracing the First Days of the SARS-CoV-2 Pandemic in Greece and the Role of the First Imported Group of Travelers.

The first SARS-CoV-2 case in Greece was confirmed on February 26, 2020, and since then, multiple strains have circulated the country, leading to regional and country-wide outbreaks. Our aim is to enlighten the events that took place during the first days of the SARS-CoV-2 pandemic in Greece, focusing on the role of the first imported group of travelers. We used whole-genome SARS-CoV-2 sequences obtained from the infected travelers of the group as well as Greece-derived and globally subsampled sequences and applied dedicated phylogenetics and phylodynamics tools as well as in-house-developed bioinformatics pipelines. Our analyses reveal the genetic variants circulating in Greece during the first days of the pandemic and the role of the group's imported strains in the course of the first pandemic wave in Greece. The strain that dominated in Greece throughout the first wave, bearing the D614G mutation, was primarily imported from a certain group of travelers, while molecular and clinical data suggest that the infection of the travelers occurred in Egypt. Founder effects early in the pandemic are important for the success of certain strains, as those arriving early, several times, and to diverse locations lead to the formation of large transmission clusters that can be estimated using molecular epidemiology approaches and can be a useful surveillance tool for the prioritization of nonpharmaceutical interventions and combating present and future outbreaks. IMPORTANCE The strain that dominated in Greece during the first pandemic wave was primarily imported from a group of returning travelers in February 2020, while molecular and clinical data suggest that the origin of the transmission was Egypt. The observed molecular transmission clusters reflect the transmission dynamics of this particular strain bearing the D614G mutation while highlighting the necessity of their use as a surveillance tool for the prioritization of nonpharmaceutical interventions and combating present and future outbreaks.

Repeated out-of-Africa expansions of Helicobacter pylori driven by replacement of deleterious mutations.

Helicobacter pylori lives in the human stomach and has a population structure resembling that of its host. However, H. pylori from Europe and the Middle East trace substantially more ancestry from modern African populations than the humans that carry them. Here, we use a collection of Afro-Eurasian H. pylori genomes to show that this African ancestry is due to at least three distinct admixture events. H. pylori from East Asia, which have undergone little admixture, have accumulated many more non-synonymous mutations than African strains. European and Middle Eastern bacteria have elevated African ancestry at the sites of these mutations, implying selection to remove them during admixture. Simulations show that population fitness can be restored after bottlenecks by migration and subsequent admixture of small numbers of bacteria from non-bottlenecked populations. We conclude that recent spread of African DNA has been driven by deleterious mutations accumulated during the original out-of-Africa bottleneck.

Early Signs of Molecular Defects in iPSC-Derived Neural Stems Cells from Patients with Familial Parkinson's Disease.

Parkinson's disease (PD) is the second most common neurodegenerative disorder, classically associated with extensive loss of dopaminergic neurons of the substantia nigra pars compacta. The hallmark of the disease is the accumulation of pathogenic conformations of the presynaptic protein, α-synuclein (αSyn), and the formation of intraneuronal protein aggregate inclusions. Neurodegeneration of dopamine neurons leads to a prominent dopaminergic deficiency in the basal ganglia, responsible for motor disturbances. However, it is now recognized that the disease involves more widespread neuronal dysfunction, leading to early and late non-motor symptoms. The development of in vitro systems based on the differentiation of human-induced pluripotent stem cells provides us the unique opportunity to monitor alterations at the cellular and molecular level throughout the differentiation procedure and identify perturbations that occur early, even at the neuronal precursor stage. Here we aim to identify whether p.A53T-αSyn induced disturbances at the molecular level are already present in neural precursors. Towards this, we present data from transcriptomics analysis of control and p.A53T-αSyn NPCs showing altered expression in transcripts involved in axon guidance, adhesion, synaptogenesis, ion transport, and metabolism. The comparative analysis with the transcriptomics profile of p.A53T-αSyn neurons shows both distinct and overlapping pathways leading to neurodegeneration while meta-analysis with transcriptomics data from both neurodegenerative and neurodevelopmental disorders reveals that p.A53T-pathology has a significant overlap with the latter category. This is the first study showing that molecular dysregulation initiates early at the p.A53T-αSyn NPC level, suggesting that synucleinopathies may have a neurodevelopmental component.

Targeted Virome Sequencing Enhances Unbiased Detection and Genome Assembly of Known and Emerging Viruses-The Example of SARS-CoV-2.

Targeted virome enrichment and sequencing (VirCapSeq-VERT) utilizes a pool of oligos (baits) to enrich all known­up to 2015­vertebrate-infecting viruses, increasing their detection sensitivity. The hybridisation of the baits to the target sequences can be partial, thus enabling the detection and genomic reconstruction of novel pathogens with <40% genetic diversity compared to the strains used for the baits' design. In this study, we deploy this method in multiplexed mixes of viral extracts, and we assess its performance in the unbiased detection of DNA and RNA viruses after cDNA synthesis. We further assess its efficiency in depleting various background genomic material. Finally, as a proof-of-concept, we explore the potential usage of the method for the characterization of unknown, emerging human viruses, such as SARS-CoV-2, which may not be included in the baits' panel. We mixed positive samples of equimolar DNA/RNA viral extracts from SARS-CoV-2, coronavirus OC43, cytomegalovirus, influenza A virus H3N2, parvovirus B19, respiratory syncytial virus, adenovirus C and coxsackievirus A16. Targeted virome enrichment was performed on a dsDNA mix, followed by sequencing on the NextSeq500 (Illumina) and the portable MinION sequencer, to evaluate its usability as a point-of-care (PoC) application. Genome mapping assembly was performed using viral reference sequences. The untargeted libraries contained less than 1% of total reads mapped on most viral genomes, while RNA viruses remained undetected. In the targeted libraries, the percentage of viral-mapped reads were substantially increased, allowing full genome assembly in most cases. Targeted virome sequencing can enrich a broad range of viruses, potentially enabling the discovery of emerging viruses.

Molecular and Clinical Prognostic Biomarkers of COVID-19 Severity and Persistence.

The coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), poses several challenges to clinicians, due to its unpredictable clinical course. The identification of laboratory biomarkers, specific cellular, and molecular mediators of immune response could contribute to the prognosis and management of COVID-19 patients. Of utmost importance is also the detection of differentially expressed genes, which can serve as transcriptomic signatures, providing information valuable to stratify patients into groups, based on the severity of the disease. The role of biomarkers such as IL-6, procalcitonin, neutrophil-lymphocyte ratio, white blood cell counts, etc. has already been highlighted in recently published studies; however, there is a notable amount of new evidence that has not been summarized yet, especially regarding transcriptomic signatures. Hence, in this review, we assess the latest cellular and molecular data and determine the significance of abnormalities in potential biomarkers for COVID-19 severity and persistence. Furthermore, we applied Gene Ontology (GO) enrichment analysis using the genes reported as differentially expressed in the literature in order to investigate which biological pathways are significantly enriched. The analysis revealed a number of processes, such as inflammatory response, and monocyte and neutrophil chemotaxis, which occur as part of the complex immune response to SARS-CoV-2.

Hepatitis B Virus Integration into Transcriptionally Active Loci and HBV-Associated Hepatocellular Carcinoma.

Hepatitis B Virus (HBV) DNA integrations into the human genome are considered major causative factors to HBV-associated hepatocellular carcinoma development. In the present study, we investigated whether HBV preferentially integrates parts of its genome in specific genes and evaluated the contribution of the integrations in HCC development per gene. We applied dedicated in-house developed pipelines on all of the available HBV DNA integration data and performed a statistical analysis to identify genes that could be characterized as hotspots of integrations, along with the evaluation of their association with HBV-HCC. Our results suggest that 15 genes are recurrently affected by HBV integrations and they are significantly associated with HBV-HCC. Further studies that focus on HBV integrations disrupting these genes are mandatory in order to understand the role of HBV integrations in clonal advantage gain and oncogenesis promotion, as well as to determine whether inhibition of the HBV-disrupted genes can provide a therapy strategy for HBV-HCC.

A recurrent chromosomal inversion suffices for driving escape from oncogene-induced senescence via subTAD reorganization.

Oncogene-induced senescence (OIS) is an inherent and important tumor suppressor mechanism. However, if not removed timely via immune surveillance, senescent cells also have detrimental effects. Although this has mostly been attributed to the senescence-associated secretory phenotype (SASP) of these cells, we recently proposed that "escape" from the senescent state is another unfavorable outcome. The mechanism underlying this phenomenon remains elusive. Here, we exploit genomic and functional data from a prototypical human epithelial cell model carrying an inducible CDC6 oncogene to identify an early-acquired recurrent chromosomal inversion that harbors a locus encoding the circadian transcription factor BHLHE40. This inversion alone suffices for BHLHE40 activation upon CDC6 induction and driving cell cycle re-entry of senescent cells, and malignant transformation. Ectopic overexpression of BHLHE40 prevented induction of CDC6-triggered senescence. We provide strong evidence in support of replication stress-induced genomic instability being a causative factor underlying "escape" from oncogene-induced senescence.

Upregulation of Human Endogenous Retroviruses in Bronchoalveolar Lavage Fluid of COVID-19 Patients.

Severe COVID-19 pneumonia has been associated with the development of intense inflammatory responses during the course of infections with SARS-CoV-2. Given that human endogenous retroviruses (HERVs) are known to be activated during and participate in inflammatory processes, we examined whether HERV dysregulation signatures are present in COVID-19 patients. By comparing transcriptomes of bronchoalveolar lavage fluid (BALF) of COVID-19 patients and healthy controls, and peripheral blood monocytes (PBMCs) from patients and controls, we have shown that HERVs are intensely dysregulated in BALF of COVID-19 patients compared to those in BALF of healthy control patients but not in PBMCs. In particular, upregulation in the expression of specific HERV families was detected in BALF samples of COVID-19 patients, with HERV-FRD being the most highly upregulated family among the families analyzed. In addition, we compared the expression of HERVs in human bronchial epithelial cells (HBECs) without and after senescence induction in an oncogene-induced senescence model in order to quantitatively measure changes in the expression of HERVs in bronchial cells during the process of cellular senescence. This apparent difference of HERV dysregulation between PBMCs and BALF warrants further studies in the involvement of HERVs in inflammatory pathogenetic mechanisms as well as exploration of HERVs as potential biomarkers for disease progression. Furthermore, the increase in the expression of HERVs in senescent HBECs in comparison to that in noninduced HBECs provides a potential link for increased COVID-19 severity and mortality in aged populations. IMPORTANCE SARS-CoV-2 emerged in late 2019 in China, causing a global pandemic. Severe COVID-19 is characterized by intensive inflammatory responses, and older age is an important risk factor for unfavorable outcomes. HERVs are remnants of ancient infections whose expression is upregulated in multiple conditions, including cancer and inflammation, and their expression is increased with increasing age. The significance of this work is that we were able to recognize dysregulated expression of endogenous retroviral elements in BALF samples but not in PBMCs of COVID-19 patients. At the same time, we were able to identify upregulated expression of multiple HERV families in senescence-induced HBECs in comparison to that in noninduced HBECs, a fact that could possibly explain the differences in disease severity among age groups. These results indicate that HERV expression might play a pathophysiological role in local inflammatory pathways in lungs afflicted by SARS-CoV-2 and their expression could be a potential therapeutic target.

Dual RNA-Seq Enables Full-Genome Assembly of Measles Virus and Characterization of Host-Pathogen Interactions.

Measles virus (MeV) has a negative-sense 15 kb long RNA genome, which is generally conserved. Recent advances in high-throughput sequencing (HTS) and Dual RNA-seq allow the analysis of viral RNA genomes and the discovery of viral infection biomarkers, via the simultaneous characterization of the host transcriptome. However, these host-pathogen interactions remain largely unexplored in MeV infections. We performed untargeted Dual RNA-seq in 6 pharyngeal and 6 peripheral blood mononuclear cell (PBMCs) specimens from patients with MeV infection, as confirmed via routine real-time PCR testing. Following optimised DNase treatment of total nucleic acids, we used the pharyngeal samples to build poly-A-enriched NGS libraries. We reconstructed the viral genomes using the pharyngeal datasets and we further conducted differential expression, gene-ontology and pathways enrichment analysis to compare both the pharyngeal and the peripheral blood transcriptomes of the MeV-infected patients vs. control groups of healthy individuals. We obtained 6 MeV genotype-B3 full-genome sequences. We minutely analyzed the transcriptome of the MeV-infected pharyngeal epithelium, detecting all known viral infection biomarkers, but also revealing a functional cluster of local antiviral and inflammatory immune responses, which differ substantially from those observed in the PBMCs transcriptome. The application of Dual RNA-seq technologies in MeV-infected patients can potentially provide valuable information on the virus genome structure and the cellular innate immune responses and drive the discovery of new targets for antiviral therapy.

Hepatitis B Virus DNA Integration, Chronic Infections and Hepatocellular Carcinoma.

Hepatitis B Virus (HBV) is an Old World virus with a high mutation rate, which puts its origins in Africa alongside the origins of Homo sapiens, and is a member of the Hepadnaviridae family that is characterized by a unique viral replication cycle. It targets human hepatocytes and can lead to chronic HBV infection either after acute infection via horizontal transmission usually during infancy or childhood or via maternal-fetal transmission. HBV has been found in ~85% of HBV-related Hepatocellular Carcinomas (HCC), and it can integrate the whole or part of its genome into the host genomic DNA. The molecular mechanisms involved in the HBV DNA integration is not yet clear; thus, multiple models have been described with respect to either the relaxed-circular DNA (rcDNA) or the double-stranded linear DNA (dslDNA) of HBV. Various genes have been found to be affected by HBV DNA integration, including cell-proliferation-related genes, oncogenes and long non-coding RNA genes (lincRNAs). The present review summarizes the advances in the research of HBV DNA integration, focusing on the evolutionary and molecular side of the integration events along with the arising clinical aspects in the light of WHO's commitment to eliminate HBV and viral hepatitis by 2030.

SARS-CoV-2 Molecular Transmission Clusters and Containment Measures in Ten European Regions during the First Pandemic Wave.

BACKGROUND: The spatiotemporal profiling of molecular transmission clusters (MTCs) using viral genomic data can effectively identify transmission networks in order to inform public health actions targeting SARS-CoV-2 spread. METHODS: We used whole genome SARS-CoV-2 sequences derived from ten European regions belonging to eight countries to perform phylogenetic and phylodynamic analysis. We developed dedicated bioinformatics pipelines to identify regional MTCs and to assess demographic factors potentially associated with their formation. RESULTS: The total number and the scale of MTCs varied from small household clusters identified in all regions, to a super-spreading event found in Uusimaa-FI. Specific age groups were more likely to belong to MTCs in different regions. The clustered sequences referring to the age groups 50-100 years old (y.o.) were increased in all regions two weeks after the establishment of the lockdown, while those referring to the age group 0-19 y.o. decreased only in those regions where schools' closure was combined with a lockdown. CONCLUSIONS: The spatiotemporal profiling of the SARS-CoV-2 MTCs can be a useful tool to monitor the effectiveness of the interventions and to reveal cryptic transmissions that have not been identified through contact tracing.

Impact of Helicobacter pylori Infection and Its Major Virulence Factor CagA on DNA Damage Repair.

Helicobacter pylori infection induces a plethora of DNA damages. Gastric epithelial cells, in order to maintain genomic integrity, require an integrous DNA damage repair (DDR) machinery, which, however, is reported to be modulated by the infection. CagA is a major H. pylori virulence factor, associated with increased risk for gastric carcinogenesis. Its pathogenic activity is partly regulated by phosphorylation on EPIYA motifs. Our aim was to identify effects of H. pylori infection and CagA on DDR, investigating the transcriptome of AGS cells, infected with wild-type, ΔCagA and EPIYA-phosphorylation-defective strains. Upon RNA-Seq-based transcriptomic analysis, we observed that a notable number of DDR genes were found deregulated during the infection, potentially resulting to base excision repair and mismatch repair compromise and an intricate deregulation of nucleotide excision repair, homologous recombination and non-homologous end-joining. Transcriptome observations were further investigated on the protein expression level, utilizing infections of AGS and GES-1 cells. We observed that CagA contributed to the downregulation of Nth Like DNA Glycosylase 1 (NTHL1), MutY DNA Glycosylase (MUTYH), Flap Structure-Specific Endonuclease 1 (FEN1), RAD51 Recombinase, DNA Polymerase Delta Catalytic Subunit (POLD1), and DNA Ligase 1 (LIG1) and, contrary to transcriptome results, Apurinic/Apyrimidinic Endodeoxyribonuclease 1 (APE1) upregulation. Our study accentuates the role of CagA as a significant contributor of H. pylori infection-mediated DDR modulation, potentially disrupting the balance between DNA damage and repair, thus favoring genomic instability and carcinogenesis.

Spatiotemporal Distribution and Genetic Characterization of Measles Strains Circulating in Greece during the 2017-2018 Outbreak.

Between May 2017 and November 2018, Greece has experienced a severe measles outbreak with a total of 3258 cases reported, after reaching its goal of eliminating measles since 2014-2015. In this study, we aimed to investigate the origin and the dispersal patterns of the measles strains that circulated in Greece during this outbreak and to identify possible transmission patterns of measles virus (MeV) in the country. Of the 832 measles suspect cases referred to the National Measles and Rubella Reference Laboratory for MeV RNA detection, 131 randomly selected positive samples, representative of the temporal and spatial distribution of the laboratory-confirmed measles cases in Greece, were processed for genotypic identification by an RT-PCR amplification of a 598 bp fragment containing the 450 bp hypervariable region of the measles virus N gene. Phylogenetic analysis was carried out by the approximate maximum likelihood method (ML) under the generalized time-reversible (GTR + cat) model. All samples analyzed were found to belong to genotype B3. Comparative analysis with other European and reference measles strains revealed three separate major clusters and other multiple viruses circulating simultaneously in Greece. They were all isolated from three main community groups, Greek-Roma children, non-minority Greek nationals and immigrants/refugees, a finding that is in accordance with what was also observed in the last two measles outbreaks in 2005-2006 and 2010-2011. Notably, for one of the three clusters, no similarity was detected with previously reported prototype strains. Our results indicate the need for a more intensive vaccination program against measles amongst minority populations and in refugee hot-spots as well as the importance of molecular surveillance as a tool for monitoring measles outbreaks.

SARS-CoV-2 exhibits intra-host genomic plasticity and low-frequency polymorphic quasispecies.

In December 2019, an outbreak of atypical pneumonia (Coronavirus disease 2019 -COVID-19) associated with a novel coronavirus (SARS-CoV-2) was reported in Wuhan city, Hubei province, China. The outbreak was traced to a seafood wholesale market and human to human transmission was confirmed. The rapid spread and the death toll of the new epidemic warrants immediate intervention. The intra-host genomic variability of SARS-CoV-2 plays a pivotal role in the development of effective antiviral agents and vaccines, as well as in the design of accurate diagnostics. We analyzed NGS data derived from clinical samples of three Chinese patients infected with SARS-CoV-2, in order to identify small- and large-scale intra-host variations in the viral genome. We identified tens of low- or higher- frequency single nucleotide variations (SNVs) with variable density across the viral genome, affecting 7 out of 10 protein-coding viral genes. The majority of these SNVs (72/104) corresponded to missense changes. The annotation of the identified SNVs but also of all currently circulating strain variations revealed colocalization of intra-host as well as strain specific SNVs with primers and probes currently used in molecular diagnostics assays. Moreover, we de-novo assembled the viral genome, in order to isolate and validate intra-host structural variations and recombination breakpoints. The bioinformatics analysis disclosed genomic rearrangements over poly-A / poly-U regions located in ORF1ab and spike (S) gene, including a potential recombination hot-spot within S gene. Our results highlight the intra-host genomic diversity and plasticity of SARS-CoV-2, pointing out genomic regions that are prone to alterations. The isolated SNVs and genomic rearrangements reflect the intra-patient capacity of the polymorphic quasispecies, which may arise rapidly during the outbreak, allowing immunological escape of the virus, offering resistance to anti-viral drugs and affecting the sensitivity of the molecular diagnostics assays.

Unique Spatial Immune Profiling in Pancreatic Ductal Adenocarcinoma with Enrichment of Exhausted and Senescent T Cells and Diffused CD47-SIRPα Expression.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is resistant to single-agent immunotherapies. To understand the mechanisms leading to the poor response to this treatment, a better understanding of the PDAC immune landscape is required. The present work aims to study the immune profile in PDAC in relationship to spatial heterogeneity of the tissue microenvironment (TME) in intact tissues. METHODS: Serial section and multiplex in situ analysis were performed in 42 PDAC samples to assess gene and protein expression at single-cell resolution in the: (a) tumor center (TC), (b) invasive front (IF), (c) normal parenchyma adjacent to the tumor, and (d) tumor positive and negative draining lymph nodes (LNs). RESULTS: We observed: (a) enrichment of T cell subpopulations with exhausted and senescent phenotype in the TC, IF and tumor positive LNs; (b) a dominant type 2 immune response in the TME, which is more pronounced in the TC; (c) an emerging role of CD47-SIRPα axis; and (d) a similar immune cell topography independently of the neoadjuvant chemotherapy. CONCLUSION: This study reveals the existence of dysfunctional T lymphocytes with specific spatial distribution, thus opening a new dimension both conceptually and mechanistically in tumor-stroma interaction in PDAC with potential impact on the efficacy of immune-regulatory therapeutic modalities.

Impact of smoking cessation, coffee and bread consumption on the intestinal microbial composition among Saudis: A cross-sectional study.

The gut microbiota is often affected by the dietary and lifestyle habits of the host, resulting in a better efficacy that favors energy harvesting from the consumed food. Our objective was to characterize the composition of gut microbiota in adult Saudis and investigate possible association with lifestyle and dietary practices. Feces from 104 Saudi volunteers (48% males) were tested for microbiota by sequencing the V3-V4 region of bacterial 16S ribosomal RNA (rRNA). For all participants, data were collected related to their lifestyle habits and dietary practices. The relative abundance (RA) of Fusobacteria was significantly higher in normal weight Saudis (P = 0.005, false discovery rate-FDR = 0.014). Individuals who consumed more coffee presented marginally significant more RA of Fusobacteria (P = 0.02, FDR = 0.20) in their gut microbiota compared to those reporting low or no coffee intake, but the RA of Fusobacteria was significantly higher in smokers compared to non-smokers (P = 0.009, FDR = 0.027). The RA of Fusobacteria was also significantly higher in those reporting daily consumption of bread (P = 0.005, FDR = 0.015). At the species level, the gut microbiota of people who consumed coffee was dominated by Bacteroides thetaiotaomicron followed by Phascolarctobacterium faecium and Eubacterium rectale. Similarly, the gut microbiota of smokers was also enriched by B. thetaiotaomicron and Lactobacillus amylovorus. Smoking cessation, bread and coffee consumption induce changes in the intestinal microbial composition of Saudis. This indicates the significance of diet and lifestyle practices in the determination of the composition of the gut microbiota, which could possibly lead later to changes in metabolic profile and weight.