101
|
Vázquez-Castellanos JF, García-López R, Pérez-Brocal V, Pignatelli M, Moya A. Comparison of different assembly and annotation tools on analysis of simulated viral metagenomic communities in the gut. BMC Genomics 2014; 15:37. [PMID: 24438450 PMCID: PMC3901335 DOI: 10.1186/1471-2164-15-37] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2013] [Accepted: 01/16/2014] [Indexed: 02/07/2023] Open
Abstract
BACKGROUND The main limitations in the analysis of viral metagenomes are perhaps the high genetic variability and the lack of information in extant databases. To address these issues, several bioinformatic tools have been specifically designed or adapted for metagenomics by improving read assembly and creating more sensitive methods for homology detection. This study compares the performance of different available assemblers and taxonomic annotation software using simulated viral-metagenomic data. RESULTS We simulated two 454 viral metagenomes using genomes from NCBI's RefSeq database based on the list of actual viruses found in previously published metagenomes. Three different assembly strategies, spanning six assemblers, were tested for performance: overlap-layout-consensus algorithms Newbler, Celera and Minimo; de Bruijn graphs algorithms Velvet and MetaVelvet; and read probabilistic model Genovo. The performance of the assemblies was measured by the length of resulting contigs (using N50), the percentage of reads assembled and the overall accuracy when comparing against corresponding reference genomes. Additionally, the number of chimeras per contig and the lowest common ancestor were estimated in order to assess the effect of assembling on taxonomic and functional annotation. The functional classification of the reads was evaluated by counting the reads that correctly matched the functional data previously reported for the original genomes and calculating the number of over-represented functional categories in chimeric contigs. The sensitivity and specificity of tBLASTx, PhymmBL and the k-mer frequencies were measured by accurate predictions when comparing simulated reads against the NCBI Virus genomes RefSeq database. CONCLUSIONS Assembling improves functional annotation by increasing accurate assignations and decreasing ambiguous hits between viruses and bacteria. However, the success is limited by the chimeric contigs occurring at all taxonomic levels. The assembler and its parameters should be selected based on the focus of each study. Minimo's non-chimeric contigs and Genovo's long contigs excelled in taxonomy assignation and functional annotation, respectively.tBLASTx stood out as the best approach for taxonomic annotation for virus identification. PhymmBL proved useful in datasets in which no related sequences are present as it uses genomic features that may help identify distant taxa. The k-frequencies underperformed in all viral datasets.
Collapse
|
102
|
Džunková M, Garcia-Garcerà M, Martínez-Priego L, D’Auria G, Calafell F, Moya A. Direct squencing from the minimal number of DNA molecules needed to fill a 454 picotiterplate. PLoS One 2014; 9:e97379. [PMID: 24887077 PMCID: PMC4041646 DOI: 10.1371/journal.pone.0097379] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2014] [Accepted: 04/17/2014] [Indexed: 02/07/2023] Open
Abstract
The large amount of DNA needed to prepare a library in next generation sequencing protocols hinders direct sequencing of small DNA samples. This limitation is usually overcome by the enrichment of such samples with whole genome amplification (WGA), mostly by multiple displacement amplification (MDA) based on φ29 polymerase. However, this technique can be biased by the GC content of the sample and is prone to the development of chimeras as well as contamination during enrichment, which contributes to undesired noise during sequence data analysis, and also hampers the proper functional and/or taxonomic assignments. An alternative to MDA is direct DNA sequencing (DS), which represents the theoretical gold standard in genome sequencing. In this work, we explore the possibility of sequencing the genome of Escherichia coli fs 24 from the minimum number of DNA molecules required for pyrosequencing, according to the notion of one-bead-one-molecule. Using an optimized protocol for DS, we constructed a shotgun library containing the minimum number of DNA molecules needed to fill a selected region of a picotiterplate. We gathered most of the reference genome extension with uniform coverage. We compared the DS method with MDA applied to the same amount of starting DNA. As expected, MDA yielded a sparse and biased read distribution, with a very high amount of unassigned and unspecific DNA amplifications. The optimized DS protocol allows unbiased sequencing to be performed from samples with a very small amount of DNA.
Collapse
|
103
|
Pérez-Cobas AE, Artacho A, Ott SJ, Moya A, Gosalbes MJ, Latorre A. Structural and functional changes in the gut microbiota associated to Clostridium difficile infection. Front Microbiol 2014; 5:335. [PMID: 25309515 PMCID: PMC4163665 DOI: 10.3389/fmicb.2014.00335] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2014] [Accepted: 06/16/2014] [Indexed: 02/05/2023] Open
Abstract
Antibiotic therapy is a causative agent of severe disturbances in microbial communities. In healthy individuals, the gut microbiota prevents infection by harmful microorganisms through direct inhibition (releasing antimicrobial compounds), competition, or stimulation of the host's immune defenses. However, widespread antibiotic use has resulted in short- and long-term shifts in the gut microbiota structure, leading to a loss in colonization resistance in some cases. Consequently, some patients develop Clostridium difficile infection (CDI) after taking an antibiotic (AB) and, at present, this opportunistic pathogen is one of the main causes of antibiotic-associated diarrhea in hospitalized patients. Here, we analyze the composition and functional differences in the gut microbiota of C. difficile infected (CDI) vs. non-infected patients, both patient groups having been treated with AB therapy. To do so we used 16S rRNA gene and metagenomic 454-based pyrosequencing approaches. Samples were taken before, during and after AB treatment and were checked for the presence of the pathogen. We performed different analyses and comparisons between infected (CD+) vs. non-infected (CD-) samples, allowing proposing putative candidate taxa and functions that might protect against C. difficile colonization. Most of these potentially protective taxa belonged to the Firmicutes phylum, mainly to the order Clostridiales, while some candidate protective functions were related to aromatic amino acid biosynthesis and stress response mechanisms. We also found that CDI patients showed, in general, lower diversity and richness than non-infected, as well as an overrepresentation of members of the families Bacteroidaceae, Enterococcaceae, Lactobacillaceae and Clostridium clusters XI and XIVa. Regarding metabolic functions, we detected higher abundance of genes involved in the transport and binding of carbohydrates, ions, and others compounds as a response to an antibiotic environment.
Collapse
|
104
|
Martínez-Cano DJ, Reyes-Prieto M, Martínez-Romero E, Partida-Martínez LP, Latorre A, Moya A, Delaye L. Evolution of small prokaryotic genomes. Front Microbiol 2014; 5:742. [PMID: 25610432 PMCID: PMC4285135 DOI: 10.3389/fmicb.2014.00742] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2014] [Accepted: 12/07/2014] [Indexed: 02/05/2023] Open
Abstract
As revealed by genome sequencing, the biology of prokaryotes with reduced genomes is strikingly diverse. These include free-living prokaryotes with ∼800 genes as well as endosymbiotic bacteria with as few as ∼140 genes. Comparative genomics is revealing the evolutionary mechanisms that led to these small genomes. In the case of free-living prokaryotes, natural selection directly favored genome reduction, while in the case of endosymbiotic prokaryotes neutral processes played a more prominent role. However, new experimental data suggest that selective processes may be at operation as well for endosymbiotic prokaryotes at least during the first stages of genome reduction. Endosymbiotic prokaryotes have evolved diverse strategies for living with reduced gene sets inside a host-defined medium. These include utilization of host-encoded functions (some of them coded by genes acquired by gene transfer from the endosymbiont and/or other bacteria); metabolic complementation between co-symbionts; and forming consortiums with other bacteria within the host. Recent genome sequencing projects of intracellular mutualistic bacteria showed that previously believed universal evolutionary trends like reduced G+C content and conservation of genome synteny are not always present in highly reduced genomes. Finally, the simplified molecular machinery of some of these organisms with small genomes may be used to aid in the design of artificial minimal cells. Here we review recent genomic discoveries of the biology of prokaryotes endowed with small gene sets and discuss the evolutionary mechanisms that have been proposed to explain their peculiar nature.
Collapse
|
105
|
Oakeson KF, Gil R, Clayton AL, Dunn DM, von Niederhausern AC, Hamil C, Aoyagi A, Duval B, Baca A, Silva FJ, Vallier A, Jackson DG, Latorre A, Weiss RB, Heddi A, Moya A, Dale C. Genome degeneration and adaptation in a nascent stage of symbiosis. Genome Biol Evol 2014; 6:76-93. [PMID: 24407854 PMCID: PMC3914690 DOI: 10.1093/gbe/evt210] [Citation(s) in RCA: 134] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Symbiotic associations between animals and microbes are ubiquitous in nature, with an estimated 15% of all insect species harboring intracellular bacterial symbionts. Most bacterial symbionts share many genomic features including small genomes, nucleotide composition bias, high coding density, and a paucity of mobile DNA, consistent with long-term host association. In this study, we focus on the early stages of genome degeneration in a recently derived insect-bacterial mutualistic intracellular association. We present the complete genome sequence and annotation of Sitophilus oryzae primary endosymbiont (SOPE). We also present the finished genome sequence and annotation of strain HS, a close free-living relative of SOPE and other insect symbionts of the Sodalis-allied clade, whose gene inventory is expected to closely resemble the putative ancestor of this group. Structural, functional, and evolutionary analyses indicate that SOPE has undergone extensive adaptation toward an insect-associated lifestyle in a very short time period. The genome of SOPE is large in size when compared with many ancient bacterial symbionts; however, almost half of the protein-coding genes in SOPE are pseudogenes. There is also evidence for relaxed selection on the remaining intact protein-coding genes. Comparative analyses of the whole-genome sequence of strain HS and SOPE highlight numerous genomic rearrangements, duplications, and deletions facilitated by a recent expansion of insertions sequence elements, some of which appear to have catalyzed adaptive changes. Functional metabolic predictions suggest that SOPE has lost the ability to synthesize several essential amino acids and vitamins. Analyses of the bacterial cell envelope and genes encoding secretion systems suggest that these structures and elements have become simplified in the transition to a mutualistic association.
Collapse
|
106
|
López-Madrigal S, Beltrà A, Resurrección S, Soto A, Latorre A, Moya A, Gil R. Molecular evidence for ongoing complementarity and horizontal gene transfer in endosymbiotic systems of mealybugs. Front Microbiol 2014; 5:449. [PMID: 25206351 PMCID: PMC4144094 DOI: 10.3389/fmicb.2014.00449] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2014] [Accepted: 08/06/2014] [Indexed: 02/05/2023] Open
Abstract
Intracellular bacterial supply of essential amino acids is common among sap-feeding insects, thus complementing the scarcity of nitrogenous compounds in plant phloem. This is also the role of the two mealybug endosymbiotic systems whose genomes have been sequenced. In the nested endosymbiotic system from Planococcus citri (Pseudococcinae), "Candidatus Tremblaya princeps" and "Candidatus Moranella endobia" cooperate to synthesize essential amino acids, while in Phenacoccus avenae (Phenacoccinae) this function is performed by its single endosymbiont "Candidatus Tremblaya phenacola." However, little is known regarding the evolution of essential amino acid supplementation strategies in other mealybug systems. To address this knowledge gap, we screened for the presence of six selected loci involved in essential amino acid biosynthesis in five additional mealybug species. We found evidence of ongoing complementarity among endosymbionts from insects of subfamily Pseudococcinae, as well as horizontal gene transfer affecting endosymbionts from insects of family Phenacoccinae, providing a more comprehensive picture of the evolutionary history of these endosymbiotic systems. Additionally, we report two diagnostic motifs to help identify invasive mealybug species.
Collapse
|
107
|
D'Auria G, Peris-Bondia F, Džunková M, Mira A, Collado MC, Latorre A, Moya A. Active and secreted IgA-coated bacterial fractions from the human gut reveal an under-represented microbiota core. Sci Rep 2013; 3:3515. [PMID: 24343271 PMCID: PMC3865468 DOI: 10.1038/srep03515] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2013] [Accepted: 11/27/2013] [Indexed: 02/07/2023] Open
Abstract
Host-associated microbiota varies in distribution depending on the body area inhabited. Gut microbes are known to interact with the human immune system, maintaining gut homoeostasis. Thus, we studied whether secreted-IgA (S-IgA) coat specific microbial taxa without inducing strong immune responses. To do so, we fractionated gut microbiota by flow cytometry. We found that active and S-IgA-coated bacterial fractions were characterized by a higher diversity than those observed in raw faecal suspensions. A long-tail effect was observed in family distribution, revealing that rare bacteria represent up to 20% of total diversity. While Firmicutes was the most abundant phylum, the majority of its sequences were not assigned at the genus level. Finally, the single-cell-based approach enabled us to focus on active and S-IgA-coated bacteria. Thus, we revealed a microbiota core common to the healthy volunteers participating in the study. Interestingly, this core was composed mainly of low frequency taxa (e.g. Sphingomonadaceae).
Collapse
|
108
|
Durbán A, Abellán JJ, Jiménez-Hernández N, Artacho A, Garrigues V, Ortiz V, Ponce J, Latorre A, Moya A. Instability of the faecal microbiota in diarrhoea-predominant irritable bowel syndrome. FEMS Microbiol Ecol 2013; 86:581-9. [PMID: 23889283 DOI: 10.1111/1574-6941.12184] [Citation(s) in RCA: 62] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2013] [Revised: 06/20/2013] [Accepted: 07/18/2013] [Indexed: 02/05/2023] Open
Abstract
The irritable bowel syndrome (IBS) is a functional gastrointestinal disorder with a largely unknown aetiology and a wide range of symptoms. Most cross-sectional studies carried out so far suggest subtle alterations in the structure of the intestinal microbiota that are barely reproduced, partly because of the high inter-subject variation in the community composition and disorder-specific features. We performed a longitudinal study to explore the within-subject variation in the faecal microbiota in two patients with IBS classified into the diarrhoea subtype and the healthy spouse of one of them. Faecal communities were monitored over 6-8 weeks and analysed through metagenomic and metatranscriptomic approaches. We found a higher temporal instability in the fraction of active microbiota related to the IBS condition and fluctuating symptoms. Strong and quick shifts in the distribution of the active microbiota and changes in the global pattern of gene expression were detected in association with acute diarrhoea, whereas microbial composition and encoded functions were more stable. The specific alterations in the microbiota were barely reproduced within and between patients. Further research is needed to assess whether these changes are a consequence of the abnormal gut function in acute diarrhoeic episodes and the potential usefulness of tackling them.
Collapse
|
109
|
Pérez-Cobas AE, Gosalbes MJ, Friedrichs A, Knecht H, Artacho A, Eismann K, Otto W, Rojo D, Bargiela R, von Bergen M, Neulinger SC, Däumer C, Heinsen FA, Latorre A, Barbas C, Seifert J, dos Santos VM, Ott SJ, Ferrer M, Moya A. Gut microbiota disturbance during antibiotic therapy: a multi-omic approach. Gut 2013; 62:1591-601. [PMID: 23236009 PMCID: PMC3812899 DOI: 10.1136/gutjnl-2012-303184] [Citation(s) in RCA: 378] [Impact Index Per Article: 34.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
OBJECTIVE Antibiotic (AB) usage strongly affects microbial intestinal metabolism and thereby impacts human health. Understanding this process and the underlying mechanisms remains a major research goal. Accordingly, we conducted the first comparative omic investigation of gut microbial communities in faecal samples taken at multiple time points from an individual subjected to β-lactam therapy. METHODS The total (16S rDNA) and active (16S rRNA) microbiota, metagenome, metatranscriptome (mRNAs), metametabolome (high-performance liquid chromatography coupled to electrospray ionisation and quadrupole time-of-flight mass spectrometry) and metaproteome (ultra high performing liquid chromatography coupled to an Orbitrap MS(2) instrument [UPLC-LTQ Orbitrap-MS/MS]) of a patient undergoing AB therapy for 14 days were evaluated. RESULTS Apparently oscillatory population dynamics were observed, with an early reduction in Gram-negative organisms (day 6) and an overall collapse in diversity and possible further colonisation by 'presumptive' naturally resistant bacteria (day 11), followed by the re-growth of Gram-positive species (day 14). During this process, the maximum imbalance in the active microbial fraction occurred later (day 14) than the greatest change in the total microbial fraction, which reached a minimum biodiversity and richness on day 11; additionally, major metabolic changes occurred at day 6. Gut bacteria respond to ABs early by activating systems to avoid the antimicrobial effects of the drugs, while 'presumptively' attenuating their overall energetic metabolic status and the capacity to transport and metabolise bile acid, cholesterol, hormones and vitamins; host-microbial interactions significantly improved after treatment cessation. CONCLUSIONS This proof-of-concept study provides an extensive description of gut microbiota responses to follow-up β-lactam therapy. The results demonstrate that ABs targeting specific pathogenic infections and diseases may alter gut microbial ecology and interactions with host metabolism at a much higher level than previously assumed.
Collapse
|
110
|
Moya A. Simbiosis y parasitismo. ARBOR-CIENCIA PENSAMIENTO Y CULTURA 2013. [DOI: 10.3989/arbor.1997.i623-624.1797] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
111
|
González-Candelas F, Bracho MA, Wróbel B, Moya A. Molecular evolution in court: analysis of a large hepatitis C virus outbreak from an evolving source. BMC Biol 2013; 11:76. [PMID: 23870105 PMCID: PMC3717074 DOI: 10.1186/1741-7007-11-76] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2012] [Accepted: 05/24/2013] [Indexed: 02/07/2023] Open
Abstract
BACKGROUND Molecular phylogenetic analyses are used increasingly in the epidemiological investigation of outbreaks and transmission cases involving rapidly evolving RNA viruses. Here, we present the results of such an analysis that contributed to the conviction of an anesthetist as being responsible for the infection of 275 of his patients with hepatitis C virus. RESULTS We obtained sequences of the NS5B and E1-E2 regions in the viral genome for 322 patients suspected to have been infected by the doctor, and for 44 local, unrelated controls. The analysis of 4,184 cloned sequences of the E1-E2 region allowed us to exclude 47 patients from the outbreak. A subset of patients had known dates of infection. We used these data to calibrate a relaxed molecular clock and to determine a rough estimate of the time of infection for each patient. A similar analysis led to an estimate for the time of infection of the source. The date turned out to be 10 years before the detection of the outbreak. The number of patients infected was small at first, but it increased substantially in the months before the detection of the outbreak. CONCLUSIONS We have developed a procedure to integrate molecular phylogenetic reconstructions of rapidly evolving viral populations into a forensic setting adequate for molecular epidemiological analysis of outbreaks and transmission events. We applied this procedure to a large outbreak of hepatitis C virus caused by a single source and the results obtained played a key role in the trial that led to the conviction of the suspected source.
Collapse
|
112
|
Hernández E, Bargiela R, Diez MS, Friedrichs A, Pérez-Cobas AE, Gosalbes MJ, Knecht H, Martínez-Martínez M, Seifert J, von Bergen M, Artacho A, Ruiz A, Campoy C, Latorre A, Ott SJ, Moya A, Suárez A, Martins dos Santos VA, Ferrer M. Functional consequences of microbial shifts in the human gastrointestinal tract linked to antibiotic treatment and obesity. Gut Microbes 2013; 4:306-15. [PMID: 23782552 PMCID: PMC3744515 DOI: 10.4161/gmic.25321] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
The microbiomes in the gastrointestinal tract (GIT) of individuals receiving antibiotics and those in obese subjects undergo compositional shifts, the metabolic effects and linkages of which are not clearly understood. Herein, we set to gain insight into these effects, particularly with regard to carbohydrate metabolism, and to contribute to unravel the underlying mechanisms and consequences for health conditions. We measured the activity level of GIT carbohydrate-active enzymes toward 23 distinct sugars in adults patients (n = 2) receiving 14-d β-lactam therapy and in obese (n = 7) and lean (n = 5) adolescents. We observed that both 14 d antibiotic-treated and obese subjects showed higher and less balanced sugar anabolic capacities, with 40% carbohydrates being preferentially processed as compared with non-treated and lean patients. Metaproteome-wide metabolic reconstructions confirmed that the impaired utilization of sugars propagated throughout the pentose phosphate metabolism, which had adverse consequences for the metabolic status of the GIT microbiota. The results point to an age-independent positive association between GIT glycosidase activity and the body mass index, fasting blood glucose and insulin resistance (r ( 2) ≥ 0.95). Moreover, antibiotics altered the active fraction of enzymes controlling the thickness, composition and consistency of the mucin glycans. Our data and analyses provide biochemical insights into the effects of antibiotic usage on the dynamics of the GIT microbiota and pin-point presumptive links to obesity. The knowledge and the hypotheses generated herein lay a foundation for subsequent, systematic research that will be paramount for the design of "smart" dietary and therapeutic interventions to modulate host-microbe metabolic co-regulation in intestinal homeostasis.
Collapse
|
113
|
Pérez-Brocal V, García-López R, Vázquez-Castellanos JF, Nos P, Beltrán B, Latorre A, Moya A. Study of the viral and microbial communities associated with Crohn's disease: a metagenomic approach. Clin Transl Gastroenterol 2013; 4:e36. [PMID: 23760301 PMCID: PMC3696940 DOI: 10.1038/ctg.2013.9] [Citation(s) in RCA: 95] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
OBJECTIVES This study aimed to analyze and compare the diversity and structure of the viral and microbial communities in fecal samples from a control group of healthy volunteers and from patients affected by Crohn's disease (CD). METHODS Healthy adult controls (n=8) and patients affected by ileocolic CD (n=11) were examined for the viral and microbial communities in their feces and, in one additional case, in the intestinal tissue. Using two different approaches, we compared the viral and microbial communities in several ways: by group (patients vs. controls), entity (viruses vs. bacteria), read assembly (unassembled vs. assembled reads), and methodology (our approach vs. an existing pipeline). Differences in the viral and microbial composition, and abundance between the two groups were analyzed to identify taxa that are under- or over-represented. RESULTS A lower diversity but more variability between the CD samples in both virome and microbiome was found, with a clear distinction between groups based on the microbiome. Only ≈5% of the differential viral biomarkers are more represented in the CD group (Synechococcus phage S CBS1 and Retroviridae family viruses), compared with 95% in the control group. Unrelated patterns of bacteria and bacteriophages were observed. CONCLUSIONS Our use of an extensive database is critical to retrieve more viral hits than in previous approaches. Unrelated patterns of bacteria and bacteriophages may be due to uneven representation of certain viruses in databases, among other factors. Further characterization of Retroviridae viruses in the CD group could be of interest, given their links with immunodeficiency and the immune responses. To conclude, some methodological considerations underlying the analysis of the viral community composition and abundance are discussed.
Collapse
|
114
|
López-Madrigal S, Latorre A, Porcar M, Moya A, Gil R. Mealybugs nested endosymbiosis: going into the 'matryoshka' system in Planococcus citri in depth. BMC Microbiol 2013; 13:74. [PMID: 23548081 PMCID: PMC3620526 DOI: 10.1186/1471-2180-13-74] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2013] [Accepted: 03/25/2013] [Indexed: 02/07/2023] Open
Abstract
BACKGROUND In all branches of life there are plenty of symbiotic associations. Insects are particularly well suited to establishing intracellular symbiosis with bacteria, providing them with metabolic capabilities they lack. Essential primary endosymbionts can coexist with facultative secondary symbionts which can, eventually, establish metabolic complementation with the primary endosymbiont, becoming a co-primary. Usually, both endosymbionts maintain their cellular identity. An exception is the endosymbiosis found in mealybugs of the subfamily Pseudoccinae, such as Planococcus citri, with Moranella endobia located inside Tremblaya princeps. RESULTS We report the genome sequencing of M. endobia str. PCVAL and the comparative genomic analyses of the genomes of strains PCVAL and PCIT of both consortium partners. A comprehensive analysis of their functional capabilities and interactions reveals their functional coupling, with many cases of metabolic and informational complementation. Using comparative genomics, we confirm that both genomes have undergone a reductive evolution, although with some unusual genomic features as a consequence of coevolving in an exceptional compartmentalized organization. CONCLUSIONS M. endobia seems to be responsible for the biosynthesis of most cellular components and energy provision, and controls most informational processes for the consortium, while T. princeps appears to be a mere factory for amino acid synthesis, and translating proteins, using the precursors provided by M. endobia. In this scenario, we propose that both entities should be considered part of a composite organism whose compartmentalized scheme (somehow) resembles a eukaryotic cell.
Collapse
|
115
|
Gosalbes MJ, Llop S, Vallès Y, Moya A, Ballester F, Francino MP. Meconium microbiota types dominated by lactic acid or enteric bacteria are differentially associated with maternal eczema and respiratory problems in infants. Clin Exp Allergy 2013; 43:198-211. [PMID: 23331561 DOI: 10.1111/cea.12063] [Citation(s) in RCA: 236] [Impact Index Per Article: 21.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2012] [Revised: 10/09/2012] [Accepted: 11/01/2012] [Indexed: 02/05/2023]
Abstract
BACKGROUND Culture-dependent methods have shown that meconium, the newborn's first intestinal discharge, is not sterile, but the diversity of bacteria present in this material needs to be further characterized by means of more sensitive molecular techniques. OBJECTIVE Our aims were to characterize molecularly the meconium microbiota in term infants, to assess whether it contributes to the future microbiota of the infants' gastrointestinal tract, and to evaluate how it relates to lifestyle variables and atopy-related conditions. METHODS We applied high-throughput pyrosequencing of the 16S rRNA gene to study the meconium microbiota in twenty term newborns from a Spanish birth cohort. For comparison, we characterized the microbiota in fecal samples from seven pregnant women days before delivery and in two series of infant samples spanning the first seven months of life. We also compared our data with vaginal and skin microbiota characterized in independent studies. Different types of meconium microbiota were defined based on taxonomic composition and abundance and their associations with different factors were statistically evaluated. RESULTS The meconium microbiota differs from those in adult feces, vagina and skin, but resembles that of fecal samples from young infants. Meconium samples clustered into two types with different bacterial diversity, richness and composition. One of the types was less diverse, dominated by enteric bacteria and associated with a history of atopic eczema in the mother (P = 0.038), whereas the second type was dominated by lactic acid bacteria and associated with respiratory problems in the infant (P = 0.040). CONCLUSIONS & CLINICAL RELEVANCE Our findings suggest that the meconium microbiota has an intrauterine origin and participates in gut colonization. Although based on a small population sample, our association analyses also suggest that the type of bacteria detected in meconium is influenced by maternal factors and may have consequences for childhood health.
Collapse
|
116
|
Rodríguez-Barreiro R, Abendroth C, Vilanova C, Moya A, Porcar M. Towards a microbial thermoelectric cell. PLoS One 2013; 8:e56358. [PMID: 23468862 PMCID: PMC3582603 DOI: 10.1371/journal.pone.0056358] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2012] [Accepted: 01/14/2013] [Indexed: 02/07/2023] Open
Abstract
Microbial growth is an exothermic process. Biotechnological industries produce large amounts of heat, usually considered an undesirable by-product. In this work, we report the construction and characterization of the first microbial thermoelectric cell (MTC), in which the metabolic heat produced by a thermally insulated microbial culture is partially converted into electricity through a thermoelectric device optimized for low ΔT values. A temperature of 41°C and net electric voltage of around 250-600 mV was achieved with 1.7 L baker's yeast culture. This is the first time microbial metabolic energy has been converted into electricity with an ad hoc thermoelectric device. These results might contribute towards developing a novel strategy to harvest excess heat in the biotechnology industry, in processes such as ethanol fermentation, auto thermal aerobic digestion (ATAD) or bioremediation, which could be coupled with MTCs in a single unit to produce electricity as a valuable by-product of the primary biotechnological product. Additionally, we propose that small portable MTCs could be conceived and inoculated with suitable thermophilic of hyperthermophilic starter cultures and used for powering small electric devices.
Collapse
|
117
|
Patiño-Navarrete R, Moya A, Latorre A, Peretó J. Comparative genomics of Blattabacterium cuenoti: the frozen legacy of an ancient endosymbiont genome. Genome Biol Evol 2013; 5:351-61. [PMID: 23355305 PMCID: PMC3590773 DOI: 10.1093/gbe/evt011] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Many insect species have established long-term symbiotic relationships with intracellular bacteria. Symbiosis with bacteria has provided insects with novel ecological capabilities, which have allowed them colonize previously unexplored niches. Despite its importance to the understanding of the emergence of biological complexity, the evolution of symbiotic relationships remains hitherto a mystery in evolutionary biology. In this study, we contribute to the investigation of the evolutionary leaps enabled by mutualistic symbioses by sequencing the genome of Blattabacterium cuenoti, primary endosymbiont of the omnivorous cockroach Blatta orientalis, and one of the most ancient symbiotic associations. We perform comparative analyses between the Blattabacterium cuenoti genome and that of previously sequenced endosymbionts, namely those from the omnivorous hosts the Blattella germanica (Blattelidae) and Periplaneta americana (Blattidae), and the endosymbionts harbored by two wood-feeding hosts, the subsocial cockroach Cryptocercus punctulatus (Cryptocercidae) and the termite Mastotermes darwiniensis (Termitidae). Our study shows a remarkable evolutionary stasis of this symbiotic system throughout the evolutionary history of cockroaches and the deepest branching termite M. darwiniensis, in terms of not only chromosome architecture but also gene content, as revealed by the striking conservation of the Blattabacterium core genome. Importantly, the architecture of central metabolic network inferred from the endosymbiont genomes was established very early in Blattabacterium evolutionary history and could be an outcome of the essential role played by this endosymbiont in the host's nitrogen economy.
Collapse
|
118
|
Matas M, Picornell A, Cifuentes C, Payeras A, Bassa A, Homar F, González-Candelas F, López-Labrador FX, Moya A, Ramon MM, Castro JA. Generalized Linear Model (GLM) framework for the association of host variables and viral strains with liver fibrosis in HCV/HIV coinfected patients. INFECTION, GENETICS AND EVOLUTION : JOURNAL OF MOLECULAR EPIDEMIOLOGY AND EVOLUTIONARY GENETICS IN INFECTIOUS DISEASES 2013; 13:284-91. [PMID: 23174528 DOI: 10.1016/j.meegid.2012.10.020] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/02/2012] [Revised: 10/16/2012] [Accepted: 10/21/2012] [Indexed: 02/07/2023]
Abstract
Chronic hepatitis C virus (HCV) infection is the main cause of advanced and end-stage liver disease world-wide, and an important factor of morbidity and mortality in Human Immunodeficiency virus-1 (HIV-1) co-infected individuals. Whereas the genetic variability of HCV has been studied extensively in monoinfected patients, comprehensive analyses of both patient and virus characteristics are still scarce in HCV/HIV co-infection. In order to find correlates for liver damage, we sought to analyze demographic, epidemiological and clinical features of HCV/HIV co-infected patients along with the genetic makeup of HCV (viral subtypes and lineage studied by nucleotide sequencing and phylogenetic analysis of the NS5B region). We used the Generalized Linear Model (GLM) methodology in order to integrate data from the virus and the infected host to find predictors for liver damage. The degree of liver disease was evaluated indirectly by means of two indexes (APRI and FIB-4) and accounting for the time since infection, to estimate fibrosis progression rates. Our analyses identified a reduced number of variables (both from the virus and the host) implicated in liver damage, which included the stage of HIV infection, levels of gamma-glutamil transferase and cholesterol, and some distinct HCV phylogenetic clades.
Collapse
|
119
|
López-Madrigal S, Balmand S, Latorre A, Heddi A, Moya A, Gil R. How does Tremblaya princeps get essential proteins from its nested partner Moranella endobia in the Mealybug Planoccocus citri? PLoS One 2013; 8:e77307. [PMID: 24204799 PMCID: PMC3804617 DOI: 10.1371/journal.pone.0077307] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2013] [Accepted: 09/09/2013] [Indexed: 02/07/2023] Open
Abstract
Many insects maintain intracellular mutualistic symbiosis with a wide range of bacteria which are considered essential for their survival (primary or P-endosymbiont) and typically suffer drastic genome degradation. Progressive loss of P-endosymbiont metabolic capabilities could lead to the recruitment of co-existent facultative endosymbiont (secondary or S-endosymbiont), thus adding more complexity to the symbiotic system. Planococcus citri, among other mealybug species, harbors an unconventional nested endosymbiotic system where every Tremblaya princeps cell (β-proteobacterium) harbors many Moranella endobia cells (γ-proteobacterium). In this system, T. princeps possess one of the smallest prokaryote genome known so far. This extreme genome reduction suggests the supply of many metabolites and essential gene products by M. endobia. Although sporadic cell lysis is plausible, the bacterial participation on the regulation of the predicted molecular exchange (at least to some extent) cannot be excluded. Although the comprehensive analysis of the protein translocation ability of M. endobia PCVAL rules out the existence of specific mechanisms for the exportation of proteins from M. endobia to T. princeps, immunolocation of two M. endobia proteins points towards a non-massive but controlled protein provision. We propose a sporadic pattern for the predicted protein exportation events, which could be putatively controlled by the host and/or mediated by local osmotic stress.
Collapse
|
120
|
Pérez-Cobas AE, Artacho A, Knecht H, Ferrús ML, Friedrichs A, Ott SJ, Moya A, Latorre A, Gosalbes MJ. Differential effects of antibiotic therapy on the structure and function of human gut microbiota. PLoS One 2013; 8:e80201. [PMID: 24282523 PMCID: PMC3839934 DOI: 10.1371/journal.pone.0080201] [Citation(s) in RCA: 135] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2013] [Accepted: 09/28/2013] [Indexed: 02/07/2023] Open
Abstract
The human intestinal microbiota performs many essential functions for the host. Antimicrobial agents, such as antibiotics (AB), are also known to disturb microbial community equilibrium, thereby having an impact on human physiology. While an increasing number of studies investigate the effects of AB usage on changes in human gut microbiota biodiversity, its functional effects are still poorly understood. We performed a follow-up study to explore the effect of ABs with different modes of action on human gut microbiota composition and function. Four individuals were treated with different antibiotics and samples were taken before, during and after the AB course for all of them. Changes in the total and in the active (growing) microbiota as well as the functional changes were addressed by 16S rRNA gene and metagenomic 454-based pyrosequencing approaches. We have found that the class of antibiotic, particularly its antimicrobial effect and mode of action, played an important role in modulating the gut microbiota composition and function. Furthermore, analysis of the resistome suggested that oscillatory dynamics are not only due to antibiotic-target resistance, but also to fluctuations in the surviving bacterial community. Our results indicated that the effect of AB on the human gut microbiota relates to the interaction of several factors, principally the properties of the antimicrobial agent, and the structure, functions and resistance genes of the microbial community.
Collapse
|
121
|
Santos-Garcia D, Farnier PA, Beitia F, Zchori-Fein E, Vavre F, Mouton L, Moya A, Latorre A, Silva FJ. Complete genome sequence of "Candidatus Portiera aleyrodidarum" BT-QVLC, an obligate symbiont that supplies amino acids and carotenoids to Bemisia tabaci. J Bacteriol 2012; 194:6654-5. [PMID: 23144402 PMCID: PMC3497482 DOI: 10.1128/jb.01793-12] [Citation(s) in RCA: 58] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
The genome of "Candidatus Portiera aleyrodidarum," the primary endosymbiont of the whitefly Bemisia tabaci (Mediterranean species), is reported. It presents a reduced genome (357 kb) encoding the capability to synthetize, or participate in the synthesis of, several amino acids and carotenoids, being the first insect endosymbiont capable of supplying carotenoids.
Collapse
|
122
|
García-López R, Pérez-Brocal V, Diez-Domingo J, Moya A. Gut microbiota in children vaccinated with rotavirus vaccine. Pediatr Infect Dis J 2012; 31:1300-2. [PMID: 22828641 DOI: 10.1097/inf.0b013e318269e3ec] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
To assess the effect that the rotavirus vaccine RotaTeq may have on the gut microbiota, this study searched for differences in intestinal bacterial composition between vaccinated and unvaccinated children. Bacterial diversity in fecal samples was evaluated by pyrosequencing of the 16S rRNA gene and taxonomic analyses using bioinformatics tools. No evidence of such differences was observed.
Collapse
|
123
|
Krüger B, Liang C, Prell F, Fieselmann A, Moya A, Schuster S, Völker U, Dandekar T. Metabolic adaptation and protein complexes in prokaryotes. Metabolites 2012; 2:940-58. [PMID: 24957769 PMCID: PMC3901225 DOI: 10.3390/metabo2040940] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2012] [Revised: 11/10/2012] [Accepted: 11/12/2012] [Indexed: 02/07/2023] Open
Abstract
Protein complexes are classified and have been charted in several large-scale screening studies in prokaryotes. These complexes are organized in a factory-like fashion to optimize protein production and metabolism. Central components are conserved between different prokaryotes; major complexes involve carbohydrate, amino acid, fatty acid and nucleotide metabolism. Metabolic adaptation changes protein complexes according to environmental conditions. Protein modification depends on specific modifying enzymes. Proteins such as trigger enzymes display condition-dependent adaptation to different functions by participating in several complexes. Several bacterial pathogens adapt rapidly to intracellular survival with concomitant changes in protein complexes in central metabolism and optimize utilization of their favorite available nutrient source. Regulation optimizes protein costs. Master regulators lead to up- and downregulation in specific subnetworks and all involved complexes. Long protein half-life and low level expression detaches protein levels from gene expression levels. However, under optimal growth conditions, metabolite fluxes through central carbohydrate pathways correlate well with gene expression. In a system-wide view, major metabolic changes lead to rapid adaptation of complexes and feedback or feedforward regulation. Finally, prokaryotic enzyme complexes are involved in crowding and substrate channeling. This depends on detailed structural interactions and is verified for specific effects by experiments and simulations.
Collapse
|
124
|
Furió V, Garijo R, Durán M, Moya A, Bell JC, Sanjuán R. Relationship between within-host fitness and virulence in the vesicular stomatitis virus: correlation with partial decoupling. J Virol 2012; 86:12228-36. [PMID: 22951843 PMCID: PMC3486475 DOI: 10.1128/jvi.00755-12] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Given the parasitic nature of viruses, it is sometimes assumed that rates of viral replication and dissemination within hosts (within-host fitness) correlate with virulence. However, there is currently little empirical evidence supporting this principle. To test this, we quantified the fitness and virulence of 21 single- or double-nucleotide mutants of the vesicular stomatitis virus in baby hamster kidney cells (BHK-21). We found that, overall, these two traits correlated positively, but significant outliers were identified. Particularly, a single mutation in the conserved C terminus of the N nucleocapsid (U1323A) had a strongly deleterious fitness effect but did not alter or even slightly increased virulence. We also found a double mutant of the M matrix protein and G glycoprotein (U2617G/A3802G mutant) with high fitness yet low virulence. We further characterized these mutants in primary cultures from mouse brain cells and in vivo and found that their relative fitness values were similar to those observed in BHK-21 cells. The mutations had weak effects on the virus-induced death rate of total brain cells, although they specifically reduced neuron death rates. Furthermore, increased apoptosis levels were detected in neurons infected with the U2617G/A3802G mutant, consistent with its known inability to block interferon secretion. In vivo, this mutant had reduced virulence and, despite its low brain titer, it retained a relatively high fitness value owing to its ability to suppress competitor viruses. Overall, our results are in broad agreement with the notion that viral fitness and virulence should be positively correlated but show that certain mutations can break this association and that the fitness-virulence relationship can depend on complex virus-host and virus-virus interactions.
Collapse
|
125
|
Cabrera-Rubio R, Garcia-Núñez M, Setó L, Antó JM, Moya A, Monsó E, Mira A. Microbiome diversity in the bronchial tracts of patients with chronic obstructive pulmonary disease. J Clin Microbiol 2012; 50:3562-8. [PMID: 22915614 PMCID: PMC3486223 DOI: 10.1128/jcm.00767-12] [Citation(s) in RCA: 151] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Culture of bacteria from bronchial secretions in respiratory patients has low sensitivity and does not allow for complete assessment of microbial diversity across different bronchial compartments. In addition, a significant number of clinical studies are based on sputum samples, and it is not known to what extent they describe the real diversity of the mucosa. In order to identify previously unrecognized lower airway bacteria and to investigate the complexity and distribution of microbiota in patients with chronic obstructive pulmonary disease (COPD), we performed PCR amplification and pyrosequencing of the 16S rRNA gene in patients not showing signs or symptoms of infection. Four types of respiratory samples (sputum, bronchial aspirate, bronchoalveolar lavage, and bronchial mucosa) were taken from each individual, obtaining on average >1,000 16S rRNA sequences per sample. The total number of genera per patient was >100, showing a high diversity, with Streptococcus, Prevotella, Moraxella, Haemophilus, Acinetobacter, Fusobacterium, and Neisseria being the most commonly identified. Sputum samples showed significantly lower diversity than the other three sample types. Lower-bronchial-tree samples, i.e., bronchoalveolar lavage and bronchial mucosa, showed a very similar bacterial compositions in contrast to sputum and bronchial aspirate samples. Thus, sputum and bronchial aspirate samples are upper bronchial tree samples that are not representative of the lower bronchial mucosa flora, and bronchoalveolar lavage samples showed the results closest to those for the bronchial mucosa. Our data confirm that the bronchial tree is not sterile in COPD patients and support the existence a different microbiota in the upper and lower compartments.
Collapse
|