A METAGENOMIC STUDY OF THE TICK MIDGUT

A Metagenomic Study of the Tick Midgut Daniel Yuan, B.S. Supervisory Professor : Steven J. Norris, Ph.D. Southern tick–associated rash illness (STARI) or Master’s disease is a Lyme-like illness that occurs following bites by Amblyomma americanum, the lone-star tick. Clinical symptoms include a bull’s eye rash similar to the erythema migrans lesions of Lyme disease, as well as fever and joint pains. Lyme disease is caused by Borrelia burgdorferi and related spirochetes. However, B. burgdorferi has not been detected in STARI patients, or in ticks in the South Central U.S. The causative agent of STARI has not been identified, although it was once thought to be caused by another Borrelia species, Borrelia lonestari. Furthermore, while adult A. americanum have up to a 5.6% Borrelia lonestari infection rate, the prevalence of all Borrelia species in Texas ticks as a whole is not known. Previous studies indicate that 6%-30% of Northern Ixodes scapularis ticks are infected by Borrelia burgdorferi while only 10% of Northern A. americanum and I. scapularis ticks are infected by Borrelia species. The first specific aim of this project was to determine the bacterial community that inhabits the midgut of Texas and Northeastern ticks by using high throughput metagenomic sequencing to sequence bacterial 16S rDNA. Through the use of massively parallel 454 sequencing, we were able to individually sequence hundreds of thousands of 16S rDNA regions of the bacterial flora from 133 ticks from the New York, Missouri and Texas. The presence of previously confirmed endosymbionts, specifically the Rickettsia spp. and Coxiella spp., that are commonly found in ticks were confirmed, as well as some highly prevalent genera that were previously undocumented. Furthermore, multiple pathogenic genera sequences were often found in the same tick, suggesting the possibility of co-infection of multiple pathogenic species. The second specific aim was to use Borrelia specific primers to screen 344 individual ticks from Missouri, Texas and the Northeast to determine the prevalence of Borrelia species in ticks. To screen for Borrelia species, two housekeeping genes, uvrA and recG, were selected as well as the 16S-23S rDNA intergenic spacer. Ticks from Missouri, Texas and New York were screened. None of the Missouri or Texas ticks tested positive for Borrelia spp. The rate of I. scapularis infection by B.burgdorferi is dependent on tick feeding activity as well as reservoir availability. B. burgdorferi is endemic in the Northeast, sometimes reported as highly present in over 50% of all I. scapularis ticks. 11.6% of all New York ticks were positive for a species of Borrelia, however only 6.9% of all New York ticks were positive for B. burgdorferi. Despite being significantly lower than 50%, the results still fall in line with previous reports of about the prevalence of B. burgdorferi. 1.5% of all Texas ticks were positive for a Borrelia species, specifically B. lonestari. While this study was unable to identify the causative agent for STARI, 454 sequencing was able to provide a tremendous insight into the bacterial flora and possible pathogenic species of both the I. scapularis and the A. americanum tick.

Whole genome sequencing on the identification of pathogens involved with high mortality in piglets.

2016

Parvovirus B19 1A complete genome from a fatal case in Brazil

Parvovirus B19 (B19V) infects individuals worldwide and is associated with an ample range of pathologies and clinical manifestations. B19V is classified into three distinct genotypes, all identified in Brazil. Here, we report a complete sequence of a B19V genotype 1A that was obtained by high-throughput metagenomic sequencing. This genome provides information that will contribute to the studies on B19V epidemiology and evolution.

Pro-inflammatory flagellin proteins of prevalent motile commensal bacteria are variably abundant in the intestinal microbiome of elderly humans

Some Eubacterium and Roseburia species are among the most prevalent motile bacteria present in the intestinal microbiota of healthy adults. These flagellate species contribute "cell motility" category genes to the intestinal microbiome and flagellin proteins to the intestinal proteome. We reviewed and revised the annotation of motility genes in the genomes of six Eubacterium and Roseburia species that occur in the human intestinal microbiota and examined their respective locus organization by comparative genomics. Motility gene order was generally conserved across these loci. Five of these species harbored multiple genes for predicted flagellins. Flagellin proteins were isolated from R. inulinivorans strain A2-194 and from E. rectale strains A1-86 and M104/1. The amino-termini sequences of the R. inulinivorans and E. rectale A1-86 proteins were almost identical. These protein preparations stimulated secretion of interleukin-8 (IL-8) from human intestinal epithelial cell lines, suggesting that these flagellins were pro-inflammatory. Flagellins from the other four species were predicted to be pro-inflammatory on the basis of alignment to the consensus sequence of pro-inflammatory flagellins from the beta- and gamma-proteobacteria. Many fliC genes were deduced to be under the control of sigma(28). The relative abundance of the target Eubacterium and Roseburia species varied across shotgun metagenomes from 27 elderly individuals. Genes involved in the flagellum biogenesis pathways of these species were variably abundant in these metagenomes, suggesting that the current depth of coverage used for metagenomic sequencing (3.13-4.79 Gb total sequence in our study) insufficiently captures the functional diversity of genomes present at low (<= 1%) relative abundance. E. rectale and R. inulinivorans thus appear to synthesize complex flagella composed of flagellin proteins that stimulate IL-8 production. A greater depth of sequencing, improved evenness of sequencing and improved metagenome assembly from short reads will be required to facilitate in silico analyses of complete complex biochemical pathways for low-abundance target species from shotgun metagenomes.

Microbiome and resitome and foods of animal origin

To ensure food safety and to prevent food-borne illnesses, rapid and accurate detection of pathogenic agents is essential. It has already been demonstrated that shotgun metagenomic sequencing can be used to detect pathogens and their antibiotic resistance genes in food. In the studies presented in this thesis, the application shotgun metagenomic sequencing has been applied to investigate both the microbiome and resistome of foods of animal origin in order to assess advantages and disadvantages of shotgun metagenomic sequencing in comparison to the cultural methods. In the first study, it has been shown that shotgun metagenomics can be applied to detect microorganisms experimentally spiked in cold-smoked salmon. Nevertheless, a direct correlation between cell concentration of each spiked microorganism and number of corresponding reads cannot be established yet. In the second and third studies, the microbiomes and resistomes characterizing caeca and the corresponding carcasses of the birds reared in the conventional and antibiotic free farms were compared. The results highlighted the need to reduce sources of microbial contamination and antimicrobial resistance not only at the farm level but also at the post-harvest one. In the fourth study, it has been demonstrated that testing a single aliquot of a food homogenate is representative of the whole homogenate because biological replicates displayed overlapping taxonomic and functional composition. All in all, the results obtained confirmed that the application of shotgun metagenomic sequencing represents a powerful tool that can be used in the identification of both spoilage and pathogenic microorganism, and their resistome in foods of animal origin. However, a robust relationship between sequence read abundance and concentration of colony-forming unit must be still established.

New Fefe-hydrogenase Genes Identified In A Metagenomic Fosmid Library From A Municipal Wastewater Treatment Plant As Revealed By High-throughput Sequencing.

A fosmid metagenomic library was constructed with total community DNA obtained from a municipal wastewater treatment plant (MWWTP), with the aim of identifying new FeFe-hydrogenase genes encoding the enzymes most important for hydrogen metabolism. The dataset generated by pyrosequencing of a fosmid library was mined to identify environmental gene tags (EGTs) assigned to FeFe-hydrogenase. The majority of EGTs representing FeFe-hydrogenase genes were affiliated with the class Clostridia, suggesting that this group is the main hydrogen producer in the MWWTP analyzed. Based on assembled sequences, three FeFe-hydrogenase genes were predicted based on detection of the L2 motif (MPCxxKxxE) in the encoded gene product, confirming true FeFe-hydrogenase sequences. These sequences were used to design specific primers to detect fosmids encoding FeFe-hydrogenase genes predicted from the dataset. Three identified fosmids were completely sequenced. The cloned genomic fragments within these fosmids are closely related to members of the Spirochaetaceae, Bacteroidales and Firmicutes, and their FeFe-hydrogenase sequences are characterized by the structure type M3, which is common to clostridial enzymes. FeFe-hydrogenase sequences found in this study represent hitherto undetected sequences, indicating the high genetic diversity regarding these enzymes in MWWTP. Results suggest that MWWTP have to be considered as reservoirs for new FeFe-hydrogenase genes.

Bacterial diversity of soil under eucalyptus assessed by 16S rDNA sequencing analysis

Studies on the impact of Eucalyptus spp. on Brazilian soils have focused on soil chemical properties and isolating interesting microbial organisms. Few studies have focused on microbial diversity and ecology in Brazil due to limited coverage of traditional cultivation and isolation methods. Molecular microbial ecology methods based on PCR amplified 16S rDNA have enriched the knowledge of soils microbial biodiversity. The objective of this work was to compare and estimate the bacterial diversity of sympatric communities within soils from two areas, a native forest (NFA) and an eucalyptus arboretum (EAA). PCR primers, whose target soil metagenomic 16S rDNA were used to amplify soil DNA, were cloned using pGEM-T and sequenced to determine bacterial diversity. From the NFA soil 134 clones were analyzed, while 116 clones were analyzed from the EAA soil samples. The sequences were compared with those online at the GenBank. Phylogenetic analyses revealed differences between the soil types and high diversity in both communities. Soil from the Eucalyptus spp. arboretum was found to have a greater bacterial diversity than the soil investigated from the native forest area.

Exploring Bacterial Diversity of Endodontic Microbiota by Cloning and Sequencing 16S rRNA

Introduction: The characterization of microbial communities infecting the endodontic system in each clinical condition may help on the establishment of a correct prognosis and distinct strategies of treatment. The purpose of this study was to determine the bacterial diversity in primary endodontic infections by 16S ribosomal-RNA (rRNA) sequence analysis. Methods: Samples from root canals of untreated asymptomatic teeth (n = 12) exhibiting periapical lesions were obtained, 165 rRNA bacterial genomic libraries were constructed and sequenced, and bacterial diversity was estimated. Results: A total of 489 clones were analyzed (mean, 40.7 +/- 8.0 clones per sample). Seventy phylotypes were identified of which six were novel phylotypes belonging to the family Ruminococcaceae. The mean number of taxa per canal was 10.0, ranging from 3 to 21 per sample; 65.7% of the cloned sequences represented phylotypes for which no cultivated isolates have been reported. The most prevalent taxa were Atopobium rimae (50.0%), Dialister invisus, Pre-votella oris, Pseudoramibacter alactolyticus, and Tannerella forsythia (33.3%). Conclusions: Although several key species predominate in endodontic samples of asymptomatic cases with periapical lesions, the primary endodontic infection is characterized by a wide bacterial diversity, which is mostly represented by members of the phylum Firmicutes belonging to the class Clostridia followed by the phylum Bacteroidetes. (J Ended 2011;37:922-926)

Bacterial diversity of soil under eucalyptus assessed by 16S rDNA sequencing analysis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Algebraic approach to optimal clone selection applied in metagenomic projects

Due to the wide diversity of unknown organisms in the environment, 99% of them cannot be grown in traditional culture medium in laboratories. Therefore, metagenomics projects are proposed to study microbial communities present in the environment, from molecular techniques, especially the sequencing. Thereby, for the coming years it is expected an accumulation of sequences produced by these projects. Thus, the sequences produced by genomics and metagenomics projects present several challenges for the treatment, storing and analysis such as: the search for clones containing genes of interest. This work presents the OCI Metagenomics, which allows defines and manages dynamically the rules of clone selection in metagenomic libraries, thought an algebraic approach based on process algebra. Furthermore, a web interface was developed to allow researchers to easily create and execute their own rules to select clones in genomic sequence database. This software has been tested in metagenomic cosmid library and it was able to select clones containing genes of interest. Copyright 2010 ACM.

Metagenomic assembly and draft genome sequence of an Uncharacterized prevotella sp. from Nelore rumen

Prevotella is one of the most abundant genera in bovine rumen, although no genome has yet been assembled by a metagenomics approach applied to Brazilian Nelore. We report the draft genome sequence of Prevotella sp., comprising 2,971,040 bp, obtained using the Illumina sequencing platform. This genome includes 127 contigs and presents a low 48% GC.

Metagenomic trajectory of gut microbiome in the human lifespan

Co-evolving with the human host, gut microbiota establishes configurations, which vary under the pressure of inflammation, disease, ageing, diet and lifestyle. In order to describe the multi-stability of the microbiome-host relationship, we studied specific tracts of the bacterial trajectory during the human lifespan and we characterized peculiar deviations from the hypothetical development, caused by disease, using molecular techniques, such as phylogenetic microarray and next-generation sequencing. Firstly, we characterized the enterocyte-associated microbiota in breast-fed infants and adults, describing remarkable differences between the two groups of subjects. Successively, we investigated the impact of atopy on the development of the microbiome in Italian childrens, highlithing conspicuous deviations from the child-type microbiota of the Italian controls. To explore variation in the gut microbiota depending on geographical origins, which reflect different lifestyles, we compared the phylogenetic diversity of the intestinal microbiota of the Hadza hunter-gatherers of Tanzania and Italian adults. Additionally, we characterized the aged-type microbiome, describing the changes occurred in the metabolic potential of the gut microbiota of centenarians with respect to younger individuals, as a part of the pathophysiolology of the ageing process. Finally, we evaluated the impact of a probiotics intervention on the intestinal microbiota of elderly people, showing the repair of some age-related dysbioses. These studies contribute to elucidate several aspects of the intestinal microbiome development during the human lifespan, depicting the microbiota as an extremely plastic entity, capable of being reconfigured in response to different environmental factors and/or stressors of endogenous origin.

An Analytical Workflow for Metagenomic Data and its Application to the Study of Chronic Obstructive Pulmonary Disease (COPD)

Metagenomics is the culture-independent study of genetic material obtained directly from environmental samples. It has become a realistic approach to understanding microbial communities thanks to advances in high-throughput DNA sequencing technologies over the past decade. Current research has shown that different sites of the human body house varied bacterial communities. There is a strong correlation between an individual’s microbial community profile at a given site and disease. Metagenomics is being applied more often as a means of comparing microbial profiles in biomedical studies. The analysis of the data collected using metagenomics can be quite challenging and there exist a plethora of tools for interpreting the results. An automatic analytical workflow for metagenomic analyses has been implemented and tested using synthetic datasets of varying quality. It is able to accurately classify bacteria by taxa and correctly estimate the richness and diversity of each set. The workflow was then applied to the study of the airways microbiome in Chronic Obstructive Pulmonary Disease (COPD). COPD is a progressive lung disease resulting in narrowing of the airways and restricted airflow. Despite being the third leading cause of death in the United States, little is known about the differences in the lung microbial community profiles of healthy individuals and COPD patients. Bronchoalveolar lavage (BAL) samples were collected from COPD patients, active or ex-smokers, and never smokers and sequenced by 454 pyrosequencing. A total of 56 individuals were recruited for the study. Substantial colonization of the lungs was found in all subjects and differentially abundant genera in each group were identified. These discoveries are promising and may further our understanding of how the structure of the lung microbiome is modified as COPD progresses. It is also anticipated that the results will eventually lead to improved treatments for COPD.

A metagenomic comparison of endemic viruses from broiler chickens with runting stunting syndrome and from normal birds

Runting-stunting syndrome (RSS) in broiler chickens is an enteric disease that causes significant economic losses to poultry producers worldwide due to elevated feed conversion ratios, decreased body weight during growth, and excessive culling. Of specific interest are the viral agents associated with RSS which have been difficult to fully characterise to date. Past research into the aetiology of RSS has implicated a wide variety of RNA and DNA viruses however, to date, no individual virus has been identified as the main agent of RSS and the current opinion is that it may be caused by a community of viruses, collectively known as the virome. This paper attempts to characterise the viral pathogens associated with 2 – 3 week old RSS-affected and unaffected broiler chickens using next-generation sequencing and comparative metagenomics. Analysis of the viromes identified a total of 20 DNA & RNA viral families, along with 2 unidentified categories, comprised of 31 distinct viral genera and 7 unclassified genera. The most abundant viral families identified in this study were the Astroviridae, Caliciviridae, Picornaviridae, Parvoviridae, Coronaviridae, Siphoviridae, and Myoviridae. This study has identified historically significant viruses associated with the disease such as chicken astrovirus, avian nephritis virus, chicken parvovirus, and chicken calicivirus along with relatively novel viruses such as chicken megrivirus and sicinivirus 1 and will help expand the knowledge related to enteric disease in broiler chickens, provide insights into the viral constituents of a healthy avian gut, and identify a variety of enteric viruses and viral communities appropriate for further study.

Genometa - A Fast and Accurate Classifier for Short Metagenomic Shotgun Reads

Metagenomic studies use high-throughput sequence data to investigate microbial communities in situ. However, considerable challenges remain in the analysis of these data, particularly with regard to speed and reliable analysis of microbial species as opposed to higher level taxa such as phyla. We here present Genometa, a computationally undemanding graphical user interface program that enables identification of bacterial species and gene content from datasets generated by inexpensive high-throughput short read sequencing technologies. Our approach was first verified on two simulated metagenomic short read datasets, detecting 100% and 94% of the bacterial species included with few false positives or false negatives. Subsequent comparative benchmarking analysis against three popular metagenomic algorithms on an Illumina human gut dataset revealed Genometa to attribute the most reads to bacteria at species level (i.e. including all strains of that species) and demonstrate similar or better accuracy than the other programs. Lastly, speed was demonstrated to be many times that of BLAST due to the use of modern short read aligners. Our method is highly accurate if bacteria in the sample are represented by genomes in the reference sequence but cannot find species absent from the reference. This method is one of the most user-friendly and resource efficient approaches and is thus feasible for rapidly analysing millions of short reads on a personal computer.