Ten years of bacterial genome sequencing: comparative-genomics-based discoveries

Recent Advances in Genomics-Based Approaches for the Development of Intracellular Bacterial Pathogen Vaccines

Infectious diseases continue to be a leading cause of morbidity and mortality worldwide. The majority of infectious diseases are caused by intracellular pathogenic bacteria (IPB). Historically, conventional vaccination drives have helped control the pathogenesis of intracellular bacteria and the emergence of antimicrobial resistance, saving millions of lives. However, in light of various limitations, many diseases that involve IPB still do not have adequate vaccines. In response to increasing demand for novel vaccine development strategies, a new area of vaccine research emerged following the advent of genomics technology, which changed the paradigm of vaccine development by utilizing the complete genomic data of microorganisms against them. It became possible to identify genes related to disease virulence, genetic patterns linked to disease virulence, as well as the genetic components that supported immunity and favorable vaccine responses. Complete genomic databases, and advancements in transcriptomics, metabolomics, structural genomics, proteomics, immunomics, pan-genomics, synthetic genomics, and population biology have allowed researchers to identify potential vaccine candidates and predict their effects in patients. New vaccines have been created against diseases for which previously there were no vaccines available, and existing vaccines have been improved. This review highlights the key issues and explores the evolution of vaccines. The increasing volume of IPB genomic data, and their application in novel genome-based techniques for vaccine development, were also examined, along with their characteristics, and the opportunities and obstacles involved. Critically, the application of genomics technology has helped researchers rapidly select and evaluate candidate antigens. Novel vaccines capable of addressing the limitations associated with conventional vaccines have been developed and pressing healthcare issues are being addressed.

The Prevalence of Aliarcobacter Species in the Fecal Microbiota of Farm Animals and Potential Effective Agents for Their Treatment: A Review of the Past Decade

There is an endless demand for livestock-originated food, so it is necessary to elucidate the hazard points for livestock breeding. Pathogens are one of the hazard points that threaten the biosecurity of farm-animal breeding and public health. As a potential foodborne pathogen, Aliarcobacter is a member of the intestinal microbiota of farm animals with and without diarrhea. Aliarcobacter spp. are capable of colonizing livestock intestines and are transmitted through the feces. Hence, they endanger slaughterhouses and milk products with fecal contamination. They also have other, rarer, vertical and horizontal transmission routes, including the offspring that abort in farm animals. Gastrointestinal symptoms and abort cases demonstrate potential financial losses to the industry. Viewed from this perspective, the global circulation of farm-animal products is a significant route for zoonotic agents, including Aliarcobacter. In the last decade, worldwide prevalence of Aliarcobacter in fecal samples has ranged from 0.8% in Italy to 100% in Turkey. Furthermore, antibiotic resistance is recognized as a new type of environmental pollutant and has become a hot topic in animal breeding and the food industry. Increasing antibiotic resistance has become a significant problem impacting productivity. The increase in antimicrobial resistance rates in Aliarcobacter is caused by the misuse of antimicrobial drugs in livestock animals, leading to the acquiring of resistance genes from other bacteria, as well as mutations in current resistance genes. The most resistant strains are A. butzleri, A. cryaerophilus, and A. skirrowii. This review analyzes recent findings from the past decade on the prevalence of Aliarcobacter in the intestinal microbiota and the current effective antibiotics against Aliarcobacter. The paper also highlights that A. cryaerophilus and A. skirrowii are found frequently in diarrheal feces, indicating that Aliarcobacter should be studied further in livestock diarrheal diseases. Moreover, Aliarcobacter-infected farm animals can be treated with only a limited number of antibiotics, such as enrofloxacin, doxycycline, oxytetracycline, and gentamicin.

Complete Genome Report of a Hydrocarbon-Degrading Sphingobium yanoikuyae S72

Sphingobium yanoikuyae S72 was isolated from the rhizosphere of sorghum plant in Mexico and we evaluated its survival and role in the degradation of some selected monoaromatic hydrocarbons and polycyclic aromatic hydrocarbons (PAHs) using minimal medium (Bushnell Hass medium (BH)) in which each of the hydrocarbons (naphthalene, phenanthrene, xylene, toluene, and biphenyl) served as sole carbon source. Gas column chromatography–mass spectrometry analysis was used to evaluate the effect of S72’s growth in the medium with the hydrocarbons. The genome of the S72 was sequenced to determine the genetic basis for the degradation of the selected hydrocarbon in S72. The genome was assembled de novo with Spades assembler and Velvet assembler and the obtained contigs were reduced to 1 manually using Consed software. Genome annotation was carried out Prokka version 1.12, and gene calling and further annotation was carried out with NCBI PGAAP. Pangenome analysis and COG annotation were done with bacteria pangenome analysis tool (BPGA) and with PATRIC online server, respectively. S72 grew effectively in the culture medium with the hydrocarbon with concentration ranging from 20–100 mg/mL for each hydrocarbon tested. S72 degraded biphenyl by 85%, phenanthrene by 93%, naphthalene by 81%, xylene by 19%, and toluene by 30%. The sequenced S72 genome was reduced to 1 contig and genome analysis revealed the presence of genes essential for the degradation of hydrocarbons in S72. A total of 126 unique genes in S72 are associated with the degradation of hydrocarbons and xenobiotics. S72 grew effectively in the tested hydrocarbon and shows good degradation efficiency. S72 will therefore be a good candidate for bioremediation of hydrocarbon contaminated soil.

Single-cell genomics for resolution of conserved bacterial genes and mobile genetic elements of the human intestinal microbiota using flow cytometry

As our understanding of the importance of the human microbiota in health and disease grows, so does our need to carefully resolve and delineate its genomic content. 16S rRNA gene-based analyses yield important insights into taxonomic composition, and metagenomics-based approaches reveal the functional potential of microbial communities. However, these methods generally fail to directly link genetic features, including bacterial genes and mobile genetic elements, to each other and to their source bacterial genomes. Further, they are inadequate to capture the microdiversity present within a genus, species, or strain of bacteria within these complex communities. Here, we present a method utilizing fluorescence-activated cell sorting for isolation of single bacterial cells, amplifying their genomes, screening them by 16S rRNA gene analysis, and selecting cells for genomic sequencing. We apply this method to both a cultured laboratory strain of Escherichia coli and human stool samples. Our analyses reveal the capacity of this method to provide nearly complete coverage of bacterial genomes when applied to isolates and partial genomes of bacterial species recovered from complex communities. Additionally, this method permits exploration and comparison of conserved and variable genomic features between individual cells. We generate assemblies of novel genomes within the Ruminococcaceae family and the Holdemanella genus by combining several 16S rRNA gene-matched single cells, and report novel prophages and conjugative transposons for both Bifidobacterium and Ruminococcaceae. Thus, we demonstrate an approach for flow cytometric separation and sequencing of single bacterial cells from the human microbiota, which yields a variety of critical insights into both the functional potential of individual microbes and the variation among those microbes. This method definitively links a variety of conserved and mobile genomic features, and can be extended to further resolve diverse elements present in the human microbiota.

Assessment of the Degradation Potential and Genomic Insights towards Phenanthrene by Dietzia psychralcaliphila JI1D

Extreme marine environments are potential sources of novel microbial isolations with dynamic metabolic activity. Dietzia psychralcaliphila J1ID was isolated from sediments originated from Deception Island, Antarctica, grown over phenanthrene. This strain was also assessed for its emulsifying activity. In liquid media, Dietzia psychralcaliphila J1ID showed 84.66% degradation of phenanthrene examined with HPLC-PDA. The identification of metabolites by GC-MS combined with its whole genome analysis provided the pathway involved in the degradation process. Whole genome sequencing indicated a genome size of 4,216,480 bp with 3961 annotated genes. The presence of a wide range of monooxygenase and dioxygenase, as well as dehydrogenase catabolic genes provided the genomic basis for the biodegradation. The strain possesses the genetic compartments for a wide range of toxic aromatic compounds, which includes the benABCD and catABC clusters. COG2146, COG4638, and COG0654 through COG analysis confirmed the genes involved in the oxygenation reaction of the hydrocarbons by the strain. Insights into assessing the depletion of phenanthrene throughout the incubation process and the genetic components involved were obtained. This study indicates the degradation potential of the strain, which can also be further expanded to other model polyaromatic hydrocarbons.

Whole genome strategies and bioremediation insight into dehalogenase-producing bacteria

An integral approach to decoding both culturable and uncultured microorganisms' metabolic activity involves the whole genome sequencing (WGS) of individual/complex microbial communities. WGS of culturable microbes, amplicon sequencing, metagenomics, and single-cell genome analysis are selective techniques integrating genetic information and biochemical mechanisms. These approaches transform microbial biotechnology into a quick and high-throughput culture-independent evaluation and exploit pollutant-degrading microbes. They are windows into enzyme regulatory bioremediation pathways (i.e., dehalogenase) and the complete bioremediation process of organohalide pollutants. While the genome sequencing technique is gaining the scientific community's interest, it is still in its infancy in the field of pollutant bioremediation. The techniques are becoming increasingly helpful in unraveling and predicting the enzyme structure and explore metabolic and biodegradation capabilities.

Whole genome characterization and phenanthrene catabolic pathway of a biofilm forming marine bacterium Pseudomonas aeruginosa PFL-P1

Pseudomonas aeruginosa is a small rod shaped Gram-negative bacterium of Gammaproteobacteria class known for its metabolic versatility. P. aeruginosa PFL-P1 was isolated from Polycyclic Aromatic Hydrocarbons (PAHs) contaminated site of Paradip Port, Odisha Coast, India. The strain showed excellent biofilm formation and could retain its ability to form biofilm grown with different PAHs in monoculture as well as co-cultures. To explore mechanistic insights of PAHs metabolism, the whole genome of the strain was sequenced. Next generation sequencing unfolded a genome size of 6,333,060 bp encoding 5857 CDSs. Gene ontology distribution assigned to a total of 2862 genes, wherein 2235 genes were allocated to biological process, 1549 genes to cellular component and 2339 genes to molecular function. A total of 318 horizontally transferred genes were identified when the genome was compared with the reference genomes of P. aeruginosa PAO1 and P. aeruginosa DSM 50071. Further comparison of P. aeruginosa PFL-P1 genome with P. putida containing TOL plasmids revealed similarities in the meta cleavage pathway employed for degradation of aromatic compounds like xylene and toluene. Gene annotation and pathway analysis unveiled 145 genes involved in xenobiotic biodegradation and metabolism. The biofilm cultures of P. aeruginosa PFL-P1 could degrade ~74% phenanthrene within 120 h while degradation increased up to ~76% in co-culture condition. GC-MS analysis indicated presence of diverse metabolites indicating the involvement of multiple pathways for one of the PAHs (phenanthrene) degradation. The strain also possesses the genetic machinery to utilize diverse toxic aromatic compounds such as naphthalene, benzoate, aminobenzoate, fluorobenzoate, toluene, xylene, styrene, atrazine, caprolactam etc. Common catabolic gene clusters such as benABCD, xylXYZ and catAB were observed within the genome of P. aeruginosa PFL-P1 which play key roles in the degradation of various toxic aromatic compounds.

Comparative Genomics of Lactobacillus crispatus from the Gut and Vagina Reveals Genetic Diversity and Lifestyle Adaptation

Lactobacillus crispatus colonizes the human feces, human vagina, and the crops and ceca of chicken. To explore the genetic characteristics and evolutionary relationships of L. crispatus isolated from different niches, we selected 37 strains isolated from the human vagina (n = 17), human feces (n = 11), and chicken feces (n = 9), and used comparative genomics to explore the genetic information of L. crispatus from the feces and vagina. No significant difference was found in the three sources of genomic features such as genome size, GC content, and number of protein coding sequences (CDS). However, in a phylogenetic tree constructed based on core genes, vagina-derived L. crispatus and feces-derived strains were each clustered separately. Therefore, the niche exerted an important impact on the evolution of L. crispatus. According to gene annotation, the L. crispatus derived from the vagina possessed a high abundance of genes related to acid tolerance, redox reactions, pullulanase, and carbohydrate-binding modules (CBMs). These genes helped L. crispatus to better adapt to the acidic environment of the vagina and obtain more nutrients, maintaining its dominance in the vagina in competition with other strains. In feces-derived bacteria, more genes encoding CRISPR/Cas system, glycoside hydrolases (GHs) family, and tetracycline/lincomycin resistance genes were found to adapt to the complex intestinal environment. This study highlights the evolutionary relationship of L. crispatus strains isolated from the vagina and feces, and the adaptation of L. crispatus to the host environment.

Bacterial diversity and functional metagenomics expounding the diversity of xenobiotics, stress, defense and CRISPR gene ontology providing eco-efficiency to Himalayan Hot Springs

Sikkim is one of the bio-diverse states of India, which harbors diverse alkaline and sulfur rich hot springs in its vicinity. However, there is a dearth of data present in terms of microbial and its functional diversity as only a few hot springs have been studied in this area. Thus, in this regard, microbial and functional diversity of two hot springs by NGS, PLFA, and culture-independent approaches were carried out. PLFA and culture-dependent analysis was complementary as the Gram-positive bacteria were abundant in both the hot springs with the dominance of phylum Firmicutes with Geobacillus. Metagenomic analysis revealed the abundance of Proteobacteria, Actinobacteria, and Firmicutes in both hot springs. Functional metagenomics suggested that both Yumthang and Reshi hot spring possess a diverse set of genes analogous to stress such as genes allied to osmotic, heat shock, and acid stresses; defense analogies such as multidrug resistance efflux pump, multidrug transport system, and β-lactamase; and CRISPR analogues such as related to Cas1, Cas2, Cas3, cmr1-5 proteins, CT1972, and CT1133 gene families. The xenobiotic analogues were found against benzoate, nitrotolune, xylene, DDT, and chlorocyclohexane/chlorobenzene degradation. Thus, these defensive mechanisms against environmental and anthropogenic hiccups and hindrances provide the eco-efficiency to such thermal habitats. The higher enzymatic, degradation, defense, stress potential and the lower percentage identity (< 95%) of isolates encourage the further exploration and exploitation of these habitats for industrial and biotechnological purposes.

Genomic comparison of Clostridium chauvoei isolates from classical and visceral clinical manifestation

Clostridium chauvoei is the etiological agent of blackleg, an infectious disease affecting cattle and small ruminants worldwide. This disease can manifest as classical blackleg, a condition in which skeletal muscles are affected and visceral blackleg, which affects the heart, sublingual muscles, and the diaphragm. The pathogenesis of the visceral form of the disease is poorly understood. The objective of this study is to determine and analyze complete genomic sequences of six C. chauvoei strains, five isolates from skeletal muscle and one isolate from a visceral case of blackleg in Brazil, to provide insights into the differences in pathogenic profiles of strains causing the different forms of disease. The full genomes of the six C. chauvoei strains were sequenced and comparative analyses were performed among these genomes and the C. chauvoei reference strain JF4335. The results of this study revealed that the genomes of the C. chauvoei strains analyzed are highly conserved; no particular differences were noted that could be associated with the two different clinical manifestations of the disease.

Genome sequence of a spore-laccase forming, BPA-degrading Bacillus sp. GZB isolated from an electronic-waste recycling site reveals insights into BPA degradation pathways

Bisphenol A (BPA) is a synthetic chemical with known deleterious effects on biota. A genome sequencing project is an important starting point for designing a suitable BPA bioremediation process, because it provides valuable genomic information about the physiological, metabolic, and genetic potential of the microbes used for the treatment. This study explored genomic insights provided by the BPA-degrading strain Bacillus sp. GZB, previously isolated from electronic-waste-dismantling site. The GZB genome is a circular chromosome, comprised of a total of 4,077,007 bp with G+C content comprising 46.2%. Genome contained 23 contigs encoded by 3881 protein-coding genes with nine rRNA and 53 tRNA genes. A comparative study demonstrated that strain GZB bloomed with some potential features as compared to other Bacillus species. In addition, strain GZB developed spore cells and displayed laccase activity while growing at elevated stress levels. Most importantly, strain GZB contained many protein-coding genes associated with BPA degradation, as well as the degradation of several other compounds. The protein-coding genes in the genome revealed the genetic mechanisms associated with the BPA degradation by strain GZB. This study predicts four possible degradation pathways for BPA, contributing to the possible use of strain GZB to remediate different polluted environments in the future.

Inhibition of Protein Secretion in Escherichia coli and Sub-MIC Effects of Arylomycin Antibiotics

At sufficient concentrations, antibiotics effectively eradicate many bacterial infections. However, during therapy, bacteria are unavoidably exposed to lower antibiotic concentrations, and sub-MIC exposure can result in a wide variety of other effects, including the induction of virulence, which can complicate therapy, or horizontal gene transfer (HGT), which can accelerate the spread of resistance genes. Bacterial type I signal peptidase (SPase) is an essential protein that acts at the final step of the general secretory pathway. This pathway is required for the secretion of many proteins, including many required for virulence, and the arylomycins are a class of natural product antibiotics that target SPase. Here, we investigated the consequences of exposing Escherichia coli cultures to sub-MIC levels of an arylomycin. Using multidimensional protein identification technology mass spectrometry, we found that arylomycin treatment inhibits the proper extracytoplasmic localization of many proteins, both those that appear to be SPase substrates and several that do not. The identified proteins are involved in a broad range of extracytoplasmic processes and include a number of virulence factors. The effects of arylomycin on several processes required for virulence were then individually examined, and we found that, at even sub-MIC levels, the arylomycins potently inhibit flagellation, motility, biofilm formation, and the dissemination of antibiotic resistance via HGT. Thus, we conclude that the arylomycins represent promising novel therapeutics with the potential to eradicate infections while simultaneously reducing virulence and the dissemination of resistance.

Genome Sequence and Metabolic Analysis of a Fluoranthene-Degrading Strain Pseudomonas aeruginosa DN1

Pseudomonas aeruginosa DN1, isolated from petroleum-contaminated soil, showed excellent degradation ability toward diverse polycyclic aromatic hydrocarbons (PAHs). Many studies have been done to improve its degradation ability. However, the molecular mechanisms of PAHs degradation in DN1 strain are unclear. In this study, the whole genome of DN1 strain was sequenced and analyzed. Its genome contains 6,641,902 bp and encodes 6,684 putative open reading frames (ORFs), which has the largest genome in almost all the comparative Pseudomonas strains. Results of gene annotation showed that this strain harbored over 100 candidate genes involved in PAHs degradation, including those encoding 25 dioxygenases, four ring-hydroxylating dioxygenases, five ring-cleaving dioxygenases, and various catabolic enzymes, transcriptional regulators, and transporters in the degradation pathways. In addition, gene knockout experiments revealed that the disruption of some key PAHs degradation genes in DN1 strain, such as catA, pcaG, pcaH, and rhdA, did not completely inhibit fluoranthene degradation, even though their degradative rate reduced to some extent. Three intermediate metabolites, including 9-hydroxyfluorene, 1-acenaphthenone, and 1, 8-naphthalic anhydride, were identified as the dominating intermediates in presence of 50 μg/mL fluoranthene as the sole carbon source according to gas chromatography mass spectrometry analysis. Taken together, the genomic and metabolic analysis indicated that the fluoranthene degradation by DN1 strain was initiated by dioxygenation at the C-1, 2-, C-2, 3-, and C-7, 8- positions. These results provide new insights into the genomic plasticity and environmental adaptation of DN1 strain.

The Xylella fastidosa RTX operons: evidence for the evolution of protein mosaics through novel genetic exchanges

BACKGROUND

Xylella fastidiosa (Xf) is a gram negative bacterium inhabiting the plant vascular system. In most species this bacterium lives as a benign symbiote, but in several agriculturally important plants (e.g. coffee, citrus, grapevine) Xf is pathogenic. Xf has four loci encoding homologues to hemolysin RTX proteins, virulence factors involved in a wide range of plant pathogen interactions.

RESULTS

We show that all four genes are expressed during pathogenesis in grapevine. The sequences from these four genes have a complex repetitive structure. At the C-termini, sequence diversity between strains is what would be expected from orthologous genes. However, within strains there is no N-terminal homology, indicating these loci encode RTXs of different functions and/or specificities. More striking is that many of the orthologous loci between strains share this extreme variation at the N-termini. Thus these RTX orthologues are most easily visualized as fusions between the orthologous C-termini and different N-termini. Further, the four genes are found in operons having a peculiar structure with an extensively duplicated module encoding a small protein with homology to the N-terminal region of the full length RTX. Surprisingly, some of these small peptides are most similar not to their corresponding full length RTX, but to the N-termini of RTXs from other Xf strains, and even other remotely related species.

CONCLUSIONS

These results demonstrate that these genes are expressed in planta during pathogenesis. Their structure suggests extensive evolutionary restructuring through horizontal gene transfers and heterologous recombination mechanisms. The sum of the evidence suggests these repetitive modules are a novel kind of mobile genetic element.

Genome sequence and comparative analysis of Bacillus cereus BC04, reveals genetic diversity and alterations for antimicrobial resistance

In this study, we delineated the genome sequence of a Bacillus cereus strain BC04 isolated from a stool sample in India. The draft genome is 5.1 Mb in size and consists of total 109 scaffolds, GC content is 35.2% with 5182 coding genes. The comparative analysis with other completely sequenced genomes highlights the unique presence of genomic islands, hemolysin, capsular synthetic protein, modifying enzymes accC7 and catA15, regulators of antibiotic resistance MarR and LysR with annotated functions related to virulence, stress response, and antimicrobial resistance. Overall, this study not only signifies the genetic diversity in gut isolate BC04 in particular, but also pinpoints the presence of unique genes possessed by B. cereus which can be pertinently exploited to design novel drugs and intervention strategies for the treatment of food borne diseases.

Exploitation of microbial forensics and nanotechnology for the monitoring of emerging pathogens

Emerging infectious diseases remain among the leading causes of global mortality. Traditional laboratory diagnostic approaches designed to detect and track infectious disease agents provide a framework for surveillance of bio threats. However, surveillance and outbreak investigations using such time-consuming approaches for early detection of pathogens remain the major pitfall. Hence, reasonable real-time surveillance systems to anticipate threats to public health and environment are critical for identifying specific aetiologies and preventing the global spread of infectious disease. The current review discusses the growing need for monitoring and surveillance of pathogens with the same zeal and approach as adopted by microbial forensics laboratories, and further strengthening it by integrating with the innovative nanotechnology for rapid detection of microbial pathogens. Such innovative diagnostics platforms will help to track pathogens from high risk areas and environment by pre-emptive approach that will minimize damages. The various scenarios with the examples are discussed where the high risk associated human pathogens in particular were successfully detected using various nanotechnology approaches with potential future prospects in the field of microbial forensics.

Characterization of the genome of a Nocardia strain isolated from soils in the Qinghai-Tibetan Plateau that specifically degrades crude oil and of this biodegradation

A strain of Nocardia isolated from crude oil-contaminated soils in the Qinghai-Tibetan Plateau degrades nearly all components of crude oil. This strain was identified as Nocardia soli Y48, and its growth conditions were determined. Complete genome sequencing showed that N. soli Y48 has a 7.3 Mb genome and many genes responsible for hydrocarbon degradation, biosurfactant synthesis, emulsification and other hydrocarbon degradation-related metabolisms. Analysis of the clusters of orthologous groups (COGs) and genomic islands (GIs) revealed that Y48 has undergone significant gene transfer events to adapt to changing environmental conditions (crude oil contamination). The structural features of the genome might provide a competitive edge for the survival of N. soli Y48 in oil-polluted environments and reflect the adaptation of coexisting bacteria to distinct nutritional niches.

First report of two complete Clostridium chauvoei genome sequences and detailed in silico genome analysis

Clostridium (C.) chauvoei is a Gram-positive, spore forming, anaerobic bacterium. It causes black leg in ruminants, a typically fatal histotoxic myonecrosis. High quality circular genome sequences were generated for the C. chauvoei type strain DSM 7528^T (ATCC 10092^T) and a field strain 12S0467 isolated in Germany. The origin of replication (oriC) was comparable to that of Bacillus subtilis in structure with two regions containing DnaA boxes. Similar prophages were identified in the genomes of both C. chauvoei strains which also harbored hemolysin and bacterial spore formation genes. A CRISPR type I-B system with limited variations in the repeat number was identified. Sporulation and germination process related genes were homologous to that of the Clostridia cluster I group but novel variations for regulatory genes were identified indicative for strain specific control of regulatory events. Phylogenomics showed a higher relatedness to C. septicum than to other so far sequenced genomes of species belonging to the genus Clostridium. Comparative genome analysis of three C. chauvoei circular genome sequences revealed the presence of few inversions and translocations in locally collinear blocks (LCBs). The species genome also shows a large number of genes involved in proteolysis, genes for glycosyl hydrolases and metal iron transportation genes which are presumably involved in virulence and survival in the host. Three conserved flagellar genes (fliC) were identified in each of the circular genomes. In conclusion this is the first comparative analysis of circular genomes for the species C. chauvoei, enabling insights into genome composition and virulence factor variation.

Antimicrobial resistance of Escherichia coli isolated in newly-hatched chickens and effect of amoxicillin treatment during their growth

The use of antimicrobials in food animals is the major determinant for the propagation of resistant bacteria in the animal reservoir. However, other factors may also play a part, and in particular vertical spread between the generations has been suggested to be an important transmission pathway. The objective of this paper was to determine the resistance patterns of Escherichia coli isolated from newly-hatched chickens as well as to study the antibiotic pressure effect when amoxicillin was administered during their growing period. With this aim, meconium from 22 one-day-old Ross chickens was analysed. In addition, during their growth period, amoxicillin treatments at days 7, 21 and 35 were carried out. Results showed a high number of E. coli-resistant strains were isolated from the treated one-day-old chickens, and were the highest for β-lactams group, followed by quinolone and tetracyclines. After treatment with amoxicillin, the highest percentage of resistances were detected for this antibiotic compared to the others analysed, with significant differences in resistance percentages between control and treated broilers detected in relation to ampicillin, cephalothin, streptomycin, kanamycin, gentamicin, chloramphenicol and tetracycline. Differences in resistances to ciprofloxacin and nalidixic acid between control and treated animals were not observed and there was lack of resistance for amikacin and ceftriaxone. These results suggest the possibility of vertical transmission of resistant strains to newly-hatched chicks from parent flocks, and seem to indicate that the treatment with amoxicillin increased the resistance of E. coli to other antibiotics.

The Complex Relationship between Virulence and Antibiotic Resistance

Antibiotic resistance, prompted by the overuse of antimicrobial agents, may arise from a variety of mechanisms, particularly horizontal gene transfer of virulence and antibiotic resistance genes, which is often facilitated by biofilm formation. The importance of phenotypic changes seen in a biofilm, which lead to genotypic alterations, cannot be overstated. Irrespective of if the biofilm is single microbe or polymicrobial, bacteria, protected within a biofilm from the external environment, communicate through signal transduction pathways (e.g., quorum sensing or two-component systems), leading to global changes in gene expression, enhancing virulence, and expediting the acquisition of antibiotic resistance. Thus, one must examine a genetic change in virulence and resistance not only in the context of the biofilm but also as inextricably linked pathologies. Observationally, it is clear that increased virulence and the advent of antibiotic resistance often arise almost simultaneously; however, their genetic connection has been relatively ignored. Although the complexities of genetic regulation in a multispecies community may obscure a causative relationship, uncovering key genetic interactions between virulence and resistance in biofilm bacteria is essential to identifying new druggable targets, ultimately providing a drug discovery and development pathway to improve treatment options for chronic and recurring infection.

Surface proteome mining for identification of potential vaccine candidates against Campylobacter jejuni: an in silico approach

Campylobacter jejuni remains a major cause of human gastroenteritis with estimated annual incidence rate of 450 million infections worldwide. C. jejuni is a major burden to public health in both socioeconomically developing and industrialized nations. Virulence determinants involved in C. jejuni pathogenesis are multifactorial in nature and not yet fully understood. Despite the completion of the first C. jejuni genome project in 2000, there are currently no vaccines in the market against this pathogen. Traditional vaccinology approach is an arduous and time extensive task. Omics techniques coupled with sequencing data have engaged researcher's attention to reduce the time and resources applied in the process of vaccine development. Recently, there has been remarkable increase in development of in silico analysis tools for efficiently mining biological information obscured in the genome. In silico approaches have been crucial for combating infectious diseases by accelerating the pace of vaccine development. This study employed a range of bioinformatics approaches for proteome scale identification of peptide vaccine candidates. Whole proteome of C. jejuni was investigated for varied properties like antigenicity, allergenicity, major histocompatibility class (MHC)-peptide interaction, immune cell processivity, HLA distribution, conservancy, and population coverage. Predicted epitopes were further tested for binding in MHC groove using computational docking studies. The predicted epitopes were conserved; covered more than 80 % of the world population and were presented by MHC-I supertypes. We conclude by underscoring that the epitopes predicted are believed to expedite the development of successful vaccines to control or prevent C. jejuni infections albeit the results need to be experimentally validated.

Heterogeneity among Isolates Reveals that Fitness in Low Oxygen Correlates with Aspergillus fumigatus Virulence

Previous work has shown that environmental and clinical isolates of Aspergillus fumigatus represent a diverse population that occupies a variety of niches, has extensive genetic diversity, and exhibits virulence heterogeneity in a number of animal models of invasive pulmonary aspergillosis (IPA). However, mechanisms explaining differences in virulence among A. fumigatus isolates remain enigmatic. Here, we report a significant difference in virulence of two common lab strains, CEA10 and AF293, in the murine triamcinolone immunosuppression model of IPA, in which we previously identified severe low oxygen microenvironments surrounding fungal lesions. Therefore, we hypothesize that the ability to thrive within these lesions of low oxygen promotes virulence of A. fumigatus in this model. To test this hypothesis, we performed in vitro fitness and in vivo virulence analyses in the triamcinolone murine model of IPA with 14 environmental and clinical isolates of A. fumigatus Among these isolates, we observed a strong correlation between fitness in low oxygen in vitro and virulence. In further support of our hypothesis, experimental evolution of AF293, a strain that exhibits reduced fitness in low oxygen and reduced virulence in the triamcinolone model of IPA, results in a strain (EVOL20) that has increased hypoxia fitness and a corresponding increase in virulence. Thus, the ability to thrive in low oxygen correlates with virulence of A. fumigatus isolates in the context of steroid-mediated murine immunosuppression.

IMPORTANCE

Aspergillus fumigatus occupies multiple environmental niches, likely contributing to the genotypic and phenotypic heterogeneity among isolates. Despite reports of virulence heterogeneity, pathogenesis studies often utilize a single strain for the identification and characterization of virulence and immunity factors. Here, we describe significant variation between A. fumigatus isolates in hypoxia fitness and virulence, highlighting the advantage of including multiple strains in future studies. We also illustrate that hypoxia fitness correlates strongly with increased virulence exclusively in the nonleukopenic murine triamcinolone immunosuppression model of IPA. Through an experimental evolution experiment, we observe that chronic hypoxia exposure results in increased virulence of A. fumigatus We describe here the first observation of a model-specific virulence phenotype correlative with in vitro fitness in hypoxia and pave the way for identification of hypoxia-mediated mechanisms of virulence in the fungal pathogen A. fumigatus.

Impact of genomics on the understanding of microbial evolution and classification: the importance of Darwin's views on classification

Analyses of genome sequences, by some approaches, suggest that the widespread occurrence of horizontal gene transfers (HGTs) in prokaryotes disguises their evolutionary relationships and have led to questioning of the Darwinian model of evolution for prokaryotes. These inferences are critically examined in the light of comparative genome analysis, characteristic synapomorphies, phylogenetic trees and Darwin's views on examining evolutionary relationships. Genome sequences are enabling discovery of numerous molecular markers (synapomorphies) such as conserved signature indels (CSIs) and conserved signature proteins (CSPs), which are distinctive characteristics of different prokaryotic taxa. Based on these molecular markers, exhibiting high degree of specificity and predictive ability, numerous prokaryotic taxa of different ranks, currently identified based on the 16S rRNA gene trees, can now be reliably demarcated in molecular terms. Within all studied groups, multiple CSIs and CSPs have been identified for successive nested clades providing reliable information regarding their hierarchical relationships and these inferences are not affected by HGTs. These results strongly support Darwin's views on evolution and classification and supplement the current phylogenetic framework based on 16S rRNA in important respects. The identified molecular markers provide important means for developing novel diagnostics, therapeutics and for functional studies providing important insights regarding prokaryotic taxa.

Genomic perspectives on the evolution and spread of bacterial pathogens

Since the first complete sequencing of a free-living organism, Haemophilus influenzae, genomics has been used to probe both the biology of bacterial pathogens and their evolution. Single-genome approaches provided information on the repertoire of virulence determinants and host-interaction factors, and, along with comparative analyses, allowed the proposal of hypotheses to explain the evolution of many of these traits. These analyses suggested many bacterial pathogens to be of relatively recent origin and identified genome degradation as a key aspect of host adaptation. The advent of very-high-throughput sequencing has allowed for detailed phylogenetic analysis of many important pathogens, revealing patterns of global and local spread, and recent evolution in response to pressure from therapeutics and the human immune system. Such analyses have shown that bacteria can evolve and transmit very rapidly, with emerging clones showing adaptation and global spread over years or decades. The resolution achieved with whole-genome sequencing has shown considerable benefits in clinical microbiology, enabling accurate outbreak tracking within hospitals and across continents. Continued large-scale sequencing promises many further insights into genetic determinants of drug resistance, virulence and transmission in bacterial pathogens.

Molecular characterization of the human microbiome from a reproductive perspective

The process of reproduction inherently poses unique microbial challenges because it requires the transfer of gametes from one individual to the other, meanwhile preserving the integrity of the gametes and individuals from harmful microbes during the process. Advances in molecular biology techniques have expanded our understanding of the natural organisms living on and in our bodies, including those inhabiting the reproductive tract. Over the past two decades accumulating evidence has shown that the human microbiome is tightly related to health and disease states involving the different body systems, including the reproductive system. Here we introduce the science involved in the study of the human microbiome. We examine common methods currently used to characterize the human microbiome as an inseparable part of the reproductive system. Finally, we consider a few limitations, clinical implications, and the critical need for additional research in the field of human fertility.

Whole-genome comparative analysis of virulence genes unveils similarities and differences between endophytes and other symbiotic bacteria

Plant pathogens and endophytes co-exist and often interact with the host plant and within its microbial community. The outcome of these interactions may lead to healthy plants through beneficial interactions, or to disease through the inducible production of molecules known as virulence factors. Unravelling the role of virulence in endophytes may crucially improve our understanding of host-associated microbial communities and their correlation with host health. Virulence is the outcome of a complex network of interactions, and drawing the line between pathogens and endophytes has proven to be conflictive, as strain-level differences in niche overlapping, ecological interactions, state of the host's immune system and environmental factors are seldom taken into account. Defining genomic differences between endophytes and plant pathogens is decisive for understanding the boundaries between these two groups. Here we describe the major differences at the genomic level between seven grapevine endophytic test bacteria, and 12 reference strains. We describe the virulence factors detected in the genomes of the test group, as compared to endophytic and non-endophytic references, to better understand the distribution of these traits in endophytic genomes. To do this, we adopted a comparative whole-genome approach, encompassing BLAST-based searches through the GUI-based tools Mauve and BRIG as well as calculating the core and accessory genomes of three genera of enterobacteria. We outline divergences in metabolic pathways of these endophytes and reference strains, with the aid of the online platform RAST. We present a summary of the major differences that help in the drawing of the boundaries between harmless and harmful bacteria, in the spirit of contributing to a microbiological definition of endophyte.

Insights from 20 years of bacterial genome sequencing

Since the first two complete bacterial genome sequences were published in 1995, the science of bacteria has dramatically changed. Using third-generation DNA sequencing, it is possible to completely sequence a bacterial genome in a few hours and identify some types of methylation sites along the genome as well. Sequencing of bacterial genome sequences is now a standard procedure, and the information from tens of thousands of bacterial genomes has had a major impact on our views of the bacterial world. In this review, we explore a series of questions to highlight some insights that comparative genomics has produced. To date, there are genome sequences available from 50 different bacterial phyla and 11 different archaeal phyla. However, the distribution is quite skewed towards a few phyla that contain model organisms. But the breadth is continuing to improve, with projects dedicated to filling in less characterized taxonomic groups. The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas system provides bacteria with immunity against viruses, which outnumber bacteria by tenfold. How fast can we go? Second-generation sequencing has produced a large number of draft genomes (close to 90 % of bacterial genomes in GenBank are currently not complete); third-generation sequencing can potentially produce a finished genome in a few hours, and at the same time provide methlylation sites along the entire chromosome. The diversity of bacterial communities is extensive as is evident from the genome sequences available from 50 different bacterial phyla and 11 different archaeal phyla. Genome sequencing can help in classifying an organism, and in the case where multiple genomes of the same species are available, it is possible to calculate the pan- and core genomes; comparison of more than 2000 Escherichia coli genomes finds an E. coli core genome of about 3100 gene families and a total of about 89,000 different gene families. Why do we care about bacterial genome sequencing? There are many practical applications, such as genome-scale metabolic modeling, biosurveillance, bioforensics, and infectious disease epidemiology. In the near future, high-throughput sequencing of patient metagenomic samples could revolutionize medicine in terms of speed and accuracy of finding pathogens and knowing how to treat them.

Complete genome sequence analysis of Pseudomonas aeruginosa N002 reveals its genetic adaptation for crude oil degradation

The present research work reports the whole genome sequence analysis of Pseudomonas aeruginosa strain N002 isolated from crude oil contaminated soil of Assam, India having high crude oil degradation ability. The whole genome of the strain N002 was sequenced by shotgun sequencing using Ion Torrent method and complete genome sequence analysis was done. It was found that the strain N002 revealed versatility for degradation, emulsification and metabolizing of crude oil. Analysis of cluster of orthologous group (COG) revealed that N002 has significantly higher gene abundance for cell motility, lipid transport and metabolism, intracellular trafficking, secretion and vesicular transport, secondary metabolite biosynthesis, transport and catabolism, signal transduction mechanism and transcription than average levels found in other genome sequences of the same bacterial species. However, lower gene abundance for carbohydrate transport and metabolism, replication, recombination and repair, translation, ribosomal structure, biogenesis was observed in N002 than average levels of other bacterial species.

The role of prophage in plant-pathogenic bacteria

A diverse set of phage lineages is associated with the bacterial plant-pathogen genomes sequenced to date. Analysis of 37 genomes revealed 5,169 potential genes (approximately 4.3 Mbp) of phage origin, and at least 50% had no function assigned or are nonessential to phage biology. Some phytopathogens have transcriptionally active prophage genes under conditions that mimic plant infection, suggesting an association between plant disease and prophage transcriptional modulation. The role of prophages within genomes for cell biology varies. For pathogens such as Pectobacterium, Pseudomonas, Ralstonia, and Streptomyces, involvement of prophage in disease symptoms has been demonstrated. In Xylella and Xanthomonas, prophage activity is associated with genome rearrangements and strain differentiation. For other pathogens, prophage roles are yet to be established. This review integrates available information in a unique interface ( http://propnav.esalq.usp.br ) that may be assessed to improve research in prophage biology and its association with genome evolution and pathogenicity.

SPA: a short peptide assembler for metagenomic data

The metagenomic paradigm allows for an understanding of the metabolic and functional potential of microbes in a community via a study of their proteins. The substrate for protein identification is either the set of individual nucleotide reads generated from metagenomic samples or the set of contig sequences produced by assembling these reads. However, a read-based strategy using reads generated by next-generation sequencing (NGS) technologies, results in an overwhelming majority of partial-length protein predictions. A nucleotide assembly-based strategy does not fare much better, as metagenomic assemblies are typically fragmented and also leave a large fraction of reads unassembled. Here, we present a method for reconstructing complete protein sequences directly from NGS metagenomic data. Our framework is based on a novel short peptide assembler (SPA) that assembles protein sequences from their constituent peptide fragments identified on short reads. The SPA algorithm is based on informed traversals of a de Bruijn graph, defined on an amino acid alphabet, to identify probable paths that correspond to proteins. Using large simulated and real metagenomic data sets, we show that our method outperforms the alternate approach of identifying genes on nucleotide sequence assemblies and generates longer protein sequences that can be more effectively analysed.

The animal food supplement sepiolite promotes a direct horizontal transfer of antibiotic resistance plasmids between bacterial species

Animal fodder is routinely complemented with antibiotics together with other food supplements to improve growth. For instance, sepiolite is currently used as a dietary coadjuvant in animal feed, as it increases animal growth parameters and improves meat and derived final product quality. This type of food additive has so far been considered innocuous for the development and spread of antibiotic resistance. In this study, we demonstrate that sepiolite promotes the direct horizontal transfer of antibiotic resistance plasmids between bacterial species. The conditions needed for plasmid transfer (sepiolite and friction forces) occur in the digestive tracts of farm animals, which routinely receive sepiolite as a food additive. Furthermore, this effect may be aggravated by the use of antibiotics supplied as growth promoters.

Antibiotics and antibiotic resistance: a bitter fight against evolution

One of the most terrible consequences of Darwinian evolution is arguably the emergence and spread of antibiotic resistance, which is becoming a serious menace to modern societies. While spontaneous mutation, recombination and horizontal gene transfer are recognized as the main causes of this notorious phenomenon; recent research has raised awareness that sub-lethal concentrations of antibiotics can also foster resistance as an undesirable side-effect. They can produce genetic changes by different ways, including a raise of free radicals within the cell, induction of error-prone DNA-polymerases mediated by SOS response, imbalanced nucleotide metabolism or affect directly DNA. In addition to certain environmental conditions, subinhibitory concentrations of antimicrobials may increase, even more, the mutagenic effect of antibiotics. Here, we review the state of knowledge on antibiotics as promoters of antibiotic resistance.

Taxonomy becoming a driving force in genome sequencing projects

We studied the possible impact of genomic projects by comparing the number of published articles before and after the completion of the project. We found that for most species, there is no significant change in the number of citations. Also our study remarks the growing importance of taxonomy as main motivation for the sequencing of genomes.

Complete genome sequence of the novel Escherichia coli phage phAPEC8

Bacteriophage phAPEC8 is an Escherichia coli-infecting myovirus, isolated on an avian pathogenic Escherichia coli (APEC) strain. APEC strains cause colibacillosis in poultry, resulting in high mortality levels and important economic losses. Genomic analysis of the 147,737-bp double-stranded DNA phAPEC8 genome revealed that 53% of the 269 encoded proteins are unique to this phage. Its closest relatives include the Salmonella phage PVP-SE1 and the coliphage rv5, with 19% and 18% similar proteins, respectively. As such, phAPEC8 represents a novel, phylogenetically distinct clade within the Myoviridae, with molecular properties suitable for phage therapy applications.

Curli functional amyloid systems are phylogenetically widespread and display large diversity in operon and protein structure

Escherichia coli and a few other members of the Enterobacteriales can produce functional amyloids known as curli. These extracellular fibrils are involved in biofilm formation and studies have shown that they may act as virulence factors during infections. It is not known whether curli fibrils are restricted to the Enterobacteriales or if they are phylogenetically widespread. The growing number of genome-sequenced bacteria spanning many phylogenetic groups allows a reliable bioinformatic investigation of the phylogenetic diversity of the curli system. Here we show that the curli system is phylogenetically much more widespread than initially assumed, spanning at least four phyla. Curli fibrils may consequently be encountered frequently in environmental as well as pathogenic biofilms, which was supported by identification of curli genes in public metagenomes from a diverse range of habitats. Identification and comparison of curli subunit (CsgA/B) homologs show that these proteins allow a high degree of freedom in their primary protein structure, although a modular structure of tightly spaced repeat regions containing conserved glutamine, asparagine and glycine residues has to be preserved. In addition, a high degree of variability within the operon structure of curli subunits between bacterial taxa suggests that the curli fibrils might have evolved to fulfill specific functions. Variations in the genetic organization of curli genes are also seen among different bacterial genera. This suggests that some genera may utilize alternative regulatory pathways for curli expression. Comparison of phylogenetic trees of Csg proteins and the 16S rRNA genes of the corresponding bacteria showed remarkably similar overall topography, suggesting that horizontal gene transfer is a minor player in the spreading of the curli system.

Antimicrobials as promoters of genetic variation

The main causes of antibiotic resistance are the selection of naturally occurring resistant variants and horizontal gene transfer processes. In recent years, the implications of antibiotic contact or treatment in drug resistance acquisition by bacteria have been gradually more evident. The ultimate source of bacterial genetic alterations to face antibiotic toxicity is mutation. All evidence points to antibiotics, especially when present at sublethal concentrations, as responsible for increasing genetic variation and therefore participating in the emergence of antibiotic resistance. Antibiotics may cause genetic changes by means of different pathways involving an increase of free radicals inside the cell or oxidative stress, by inducing error-prone polymerases mediated by SOS response, misbalancing nucleotide metabolism or acting directly on DNA. In addition, the concerted action of certain environmental conditions with subinhibitory concentrations of antimicrobials may contribute to increasing the mutagenic effect of antibiotics even more. Here we review and discuss in detail the recent advances concerning these issues and their relevance in the field of antibiotic resistance.

Comparative analysis of the Hrp pathogenicity island of Rubus- and Spiraeoideae-infecting Erwinia amylovora strains identifies the IT region as a remnant of an integrative conjugative element

The Hrp pathogenicity island (hrpPAI) of Erwinia amylovora not only encodes a type III secretion system (T3SS) and other genes required for pathogenesis on host plants, but also includes the so-called island transfer (IT) region, a region that originates from an integrative conjugative element (ICE). Comparative genomic analysis of the IT regions of two Spiraeoideae- and three Rubus-infecting strains revealed that the regions in Spiraeoideae-infecting strains were syntenic and highly conserved in length and genetic information, but that the IT regions of the Rubus-infecting strains varied in gene content and length, showing a mosaic structure. None of the ICEs in E. amylovora strains were complete, as conserved ICE genes and the left border were missing, probably due to reductive genome evolution. Comparison of the hrpPAI region of E. amylovora strains to syntenic regions from other Erwinia spp. indicates that the hrpPAI and the IT regions are the result of several insertion and deletion events that have occurred within the ICE. It also suggests that the T3SS was present in a common ancestor of the pathoadapted Erwinia spp. and that insertion and deletion events in the IT region occurred during speciation.

Selection of marker genes using whole-genome DNA polymorphism analysis

Molecular markers serve to assign individual samples to specific groups. Such markers should be easily identified and have a high discrimination power, being highly conserved within groups while showing sufficient variability between the groups that are to be distinguished. The availability of a large number of complete genomic sequences now enables the informed selection of genes as molecular markers based on the observed patterns of variability. We derived a new scoring system based on observed DNA polymorphic differences, and which uses the Bayes theorem as adapted by Wilcox. For validation, we applied this system to the problem of identifying individual species within a prokaryotic (Vibrio) and a eukaryotic (Diphyllobothrium) genus for validation. Top-scoring candidates genes Chromosome segregation ATPase and ATPase-subunit 6 showed better discrimination power in Vibrio and Diphyllobothrium, respectively, as compared to standard molecular markers (recA, dnaJ and atpA for Vibrio, and 18s rRNA, ITS and COX1 for Diphyllobothrium).

Genomic variation in Salmonella enterica core genes for epidemiological typing

BACKGROUND

Technological advances in high throughput genome sequencing are making whole genome sequencing (WGS) available as a routine tool for bacterial typing. Standardized procedures for identification of relevant genes and of variation are needed to enable comparison between studies and over time. The core genes--the genes that are conserved in all (or most) members of a genus or species--are potentially good candidates for investigating genomic variation in phylogeny and epidemiology.

RESULTS

We identify a set of 2,882 core genes clusters based on 73 publicly available Salmonella enterica genomes and evaluate their value as typing targets, comparing whole genome typing and traditional methods such as 16S and MLST. A consensus tree based on variation of core genes gives much better resolution than 16S and MLST; the pan-genome family tree is similar to the consensus tree, but with higher confidence. The core genes can be divided into two categories: a few highly variable genes and a larger set of conserved core genes, with low variance. For the most variable core genes, the variance in amino acid sequences is higher than for the corresponding nucleotide sequences, suggesting that there is a positive selection towards mutations leading to amino acid changes.

CONCLUSIONS

Genomic variation within the core genome is useful for investigating molecular evolution and providing candidate genes for bacterial genome typing. Identification of genes with different degrees of variation is important especially in trend analysis.

Exploring bacterial insertion sequences with ISfinder: objectives, uses, and future developments

We describe here the use of the ISfinder database and its associated software. ISfinder was conceived -initially as a comprehensive database for prokaryotic insertion sequences (ISs). It now includes software for visualising complete and partial IS copies in whole genomes (ISbrowser) and for high-quality genome annotation (Insertion Sequence Semi-automatic Genome Annotation).

The Microbe: The Basics of Structure, Morphology, and Physiology as They Relate to Microbial Characterization and Attribution

This chapter is meant to (1) review classical methods used to characterize and classify microbes and (2) introduce new molecular methods used in microbial characterization. The fundamental composition of microbes is discussed as well as their importance in classification of microbes into genus and species. Classical microbiological methods in general seek to define the common features of specific bacterial groups as a means of classification and identification of microbes. Thus, the focus was to describe the common features which discriminated closely related groups of organisms. In contrast, the newer molecular methods often seek to expand the classification of microbes not only as a means to organize microbial phylogeny but also to differentiate signatures between microbes identified within a species in greater detail. Molecular biology tools are used both as an adjunct to established methods and as replacement for classical methods for detection, discrimination, or identification of bacterial and viral species.

A Partial Least Squares based algorithm for parsimonious variable selection

Background

In genomics, a commonly encountered problem is to extract a subset of variables out of a large set of explanatory variables associated with one or several quantitative or qualitative response variables. An example is to identify associations between codon-usage and phylogeny based definitions of taxonomic groups at different taxonomic levels. Maximum understandability with the smallest number of selected variables, consistency of the selected variables, as well as variation of model performance on test data, are issues to be addressed for such problems.

Results

We present an algorithm balancing the parsimony and the predictive performance of a model. The algorithm is based on variable selection using reduced-rank Partial Least Squares with a regularized elimination. Allowing a marginal decrease in model performance results in a substantial decrease in the number of selected variables. This significantly improves the understandability of the model. Within the approach we have tested and compared three different criteria commonly used in the Partial Least Square modeling paradigm for variable selection; loading weights, regression coefficients and variable importance on projections. The algorithm is applied to a problem of identifying codon variations discriminating different bacterial taxa, which is of particular interest in classifying metagenomics samples. The results are compared with a classical forward selection algorithm, the much used Lasso algorithm as well as Soft-threshold Partial Least Squares variable selection.

Conclusions

A regularized elimination algorithm based on Partial Least Squares produces results that increase understandability and consistency and reduces the classification error on test data compared to standard approaches.

Exploration of community traits as ecological markers in microbial metagenomes

The rate of information collection generated by metagenomics is uncoupled with its meaningful ecological interpretation. New analytical approaches based on functional trait-based ecology may help to bridge this gap and extend the trait approach to the community level in vast and complex environmental genetic data sets. Here, we explored a set of community traits that range from nucleotidic to genomic properties in 53 metagenomic aquatic samples from the Global Ocean Sampling (GOS) expedition. We found significant differences between the community profile derived from the commonly used 16S rRNA gene and from the functional trait set. The traits proved to be valuable ecological markers by discriminating between marine ecosystems (coastal vs. open ocean) and between oceans (Atlantic vs. Indian vs. Pacific). Intertrait relationships were also assessed, and we propose some that could be further used as habitat descriptors or indicators of artefacts during sample processing. Overall, the approach presented here may help to interpret metagenomics data to gain a full understanding of microbial community patterns in a rigorous ecological framework.