Microbial Diversity in Bushmeat Samples Recovered from the Serengeti Ecosystem in Tanzania.

Bushmeat, the meat and organs derived from wildlife species, is a common source of animal protein in the diets of those living in sub-Saharan Africa and is frequently associated with zoonotic spillover of dangerous pathogens. Given the frequent consumption of bushmeat in this region and the lack of knowledge about the microbial communities associated with this meat, the microbiome of 56 fresh and processed bushmeat samples ascertained from three districts in the Western Serengeti ecosystem in Tanzania was characterized using 16S rRNA metagenomic sequencing. The results show that the most abundant phyla present in bushmeat samples include Firmicutes (67.8%), Proteobacteria (18.4%), Cyanobacteria (8.9%), and Bacteroidetes (3.1%). Regardless of wildlife species, sample condition, season, or region, the microbiome is diverse across all samples, with no significant difference in alpha or beta diversity. The findings also suggest the presence of DNA signatures of potentially dangerous zoonotic pathogens, including those from the genus Bacillus, Brucella, Coxiella, and others, in bushmeat. Together, this investigation provides a better understanding of the microbiome associated with this major food source in samples collected from the Western Serengeti in Tanzania and highlights a need for future investigations on the potential health risks associated with the harvesting, trade, and consumption of bushmeat in Sub-Saharan Africa.

iMAP: an integrated bioinformatics and visualization pipeline for microbiome data analysis.

BACKGROUND: One of the major challenges facing investigators in the microbiome field is turning large numbers of reads generated by next-generation sequencing (NGS) platforms into biological knowledge. Effective analytical workflows that guarantee reproducibility, repeatability, and result provenance are essential requirements of modern microbiome research. For nearly a decade, several state-of-the-art bioinformatics tools have been developed for understanding microbial communities living in a given sample. However, most of these tools are built with many functions that require an in-depth understanding of their implementation and the choice of additional tools for visualizing the final output. Furthermore, microbiome analysis can be time-consuming and may even require more advanced programming skills which some investigators may be lacking. RESULTS: We have developed a wrapper named iMAP (Integrated Microbiome Analysis Pipeline) to provide the microbiome research community with a user-friendly and portable tool that integrates bioinformatics analysis and data visualization. The iMAP tool wraps functionalities for metadata profiling, quality control of reads, sequence processing and classification, and diversity analysis of operational taxonomic units. This pipeline is also capable of generating web-based progress reports for enhancing an approach referred to as review-as-you-go (RAYG). For the most part, the profiling of microbial community is done using functionalities implemented in Mothur or QIIME2 platform. Also, it uses different R packages for graphics and R-markdown for generating progress reports. We have used a case study to demonstrate the application of the iMAP pipeline. CONCLUSIONS: The iMAP pipeline integrates several functionalities for better identification of microbial communities present in a given sample. The pipeline performs in-depth quality control that guarantees high-quality results and accurate conclusions. The vibrant visuals produced by the pipeline facilitate a better understanding of the complex and multidimensional microbiome data. The integrated RAYG approach enables the generation of web-based reports, which provides the investigators with the intermediate output that can be reviewed progressively. The intensively analyzed case study set a model for microbiome data analysis.

ERAIZDA: a model for holistic annotation of animal infectious and zoonotic diseases.

There is an urgent need for a unified resource that integrates trans-disciplinary annotations of emerging and reemerging animal infectious and zoonotic diseases. Such data integration will provide wonderful opportunity for epidemiologists, researchers and health policy makers to make data-driven decisions designed to improve animal health. Integrating emerging and reemerging animal infectious and zoonotic disease data from a large variety of sources into a unified open-access resource provides more plausible arguments to achieve better understanding of infectious and zoonotic diseases. We have developed a model for interlinking annotations of these diseases. These diseases are of particular interest because of the threats they pose to animal health, human health and global health security. We demonstrated the application of this model using brucellosis, an infectious and zoonotic disease. Preliminary annotations were deposited into VetBioBase database (http://vetbiobase.igbb.msstate.edu). This database is associated with user-friendly tools to facilitate searching, retrieving and downloading of disease-related information. Database URL: http://vetbiobase.igbb.msstate.edu.

Transcriptomic dissection of the rice-Burkholderia glumae interaction.

BACKGROUND: Bacterial panicle blight caused by the bacterium Burkholderia glumae is an emerging disease of rice in the United States. Not much is known about this disease, the disease cycle or any source of disease resistance. To understand the interaction between rice and Burkholderia glumae, we used transcriptomics via next-generation sequencing (RNA-Seq) and bioinformatics to identify differentially expressed transcripts between resistant and susceptible interactions and formulate a model for rice resistance to the disease. RESULTS: Using inoculated young seedlings as sample tissues, we identified unique transcripts involved with resistance to bacterial panicle blight, including a PIF-like ORF1 and verified differential expression of some selected genes using qRT-PCR. These transcripts, which include resistance genes of the NBS-LRR type, kinases, transcription factors, transporters and expressed proteins with functions that are not known, have not been reported in other pathosystems including rice blast or bacterial blight. Further, functional annotation analysis reveals enrichment of defense response and programmed cell death (biological processes); ATP and protein binding (molecular functions); and mitochondrion-related (cell component) transcripts in the resistant interaction. CONCLUSION: Taken together, we formulated a model for rice resistance to bacterial panicle blight that involves an activation of previously unknown resistance genes and their activation partners upon challenge with B. glumae. Other interesting findings are that 1) though these resistance transcripts were up-regulated upon inoculation in the resistant interaction, some of them were already expressed in the water-inoculated control from the resistant genotype, but not in the water- and bacterium-inoculated samples from the susceptible genotype; 2) rice may have co-opted an ORF that was previously a part of a transposable element to aid in the resistance mechanism; and 3) resistance may have existed immediately prior to rice domestication.

Computational prediction of disease microRNAs in domestic animals.

BACKGROUND: The most important means of identifying diseases before symptoms appear is through the discovery of disease-associated biomarkers. Recently, microRNAs (miRNAs) have become highly useful biomarkers of infectious, genetic and metabolic diseases in human but they have not been well studied in domestic animals. It is probable that many of the animal homologs of human disease-associated miRNAs may be involved in domestic animal diseases. Here we describe a computational biology study in which human disease miRNAs were utilized to predict orthologous miRNAs in cow, chicken, pig, horse, and dog. RESULTS: We identified 287 human disease-associated miRNAs which had at least one 100% identical animal homolog. The 287 miRNAs were associated with 359 human diseases referenced in 2,863 Pubmed articles. Multiple sequence analysis indicated that over 60% of known horse mature miRNAs found perfect matches in human disease-associated miRNAs, followed by dog (50%). As expected, chicken had the least number of perfect matches (5%). Phylogenetic analysis of miRNA precursors indicated that 85% of human disease pre-miRNAs were highly conserved in animals, showing less than 5% nucleotide substitution rates over evolutionary time. As an example we demonstrated conservation of human hsa-miR-143-3p which is associated with type 2 diabetes and targets AKT1 gene which is highly conserved in pig, horse and dog. Functional analysis of AKT1 gene using Gene Ontology (GO) showed that it is involved in glucose homeostasis, positive regulation of glucose import, positive regulation of glycogen biosynthetic process, glucose transport and response to food. CONCLUSIONS: This data provides the animal and veterinary research community with a resource to assist in generating hypothesis-driven research for discovering animal disease-related miRNA from their datasets and expedite development of prophylactic and disease-treatment strategies and also influence research efforts to identify novel disease models in large animals. Integrated data is available for download at http://agbase.hpc.msstate.edu/cgi-bin/animal_mirna.cgi.

AgBase: supporting functional modeling in agricultural organisms.

AgBase (http://www.agbase.msstate.edu/) provides resources to facilitate modeling of functional genomics data and structural and functional annotation of agriculturally important animal, plant, microbe and parasite genomes. The website is redesigned to improve accessibility and ease of use, including improved search capabilities. Expanded capabilities include new dedicated pages for horse, cat, dog, cotton, rice and soybean. We currently provide 590 240 Gene Ontology (GO) annotations to 105 454 gene products in 64 different species, including GO annotations linked to transcripts represented on agricultural microarrays. For many of these arrays, this provides the only functional annotation available. GO annotations are available for download and we provide comprehensive, species-specific GO annotation files for 18 different organisms. The tools available at AgBase have been expanded and several existing tools improved based upon user feedback. One of seven new tools available at AgBase, GOModeler, supports hypothesis testing from functional genomics data. We host several associated databases and provide genome browsers for three agricultural pathogens. Moreover, we provide comprehensive training resources (including worked examples and tutorials) via links to Educational Resources at the AgBase website.

Facilitating functional annotation of chicken microarray data.

BACKGROUND: Modeling results from chicken microarray studies is challenging for researchers due to little functional annotation associated with these arrays. The Affymetrix GenChip chicken genome array, one of the biggest arrays that serve as a key research tool for the study of chicken functional genomics, is among the few arrays that link gene products to Gene Ontology (GO). However the GO annotation data presented by Affymetrix is incomplete, for example, they do not show references linked to manually annotated functions. In addition, there is no tool that facilitates microarray researchers to directly retrieve functional annotations for their datasets from the annotated arrays. This costs researchers amount of time in searching multiple GO databases for functional information. RESULTS: We have improved the breadth of functional annotations of the gene products associated with probesets on the Affymetrix chicken genome array by 45% and the quality of annotation by 14%. We have also identified the most significant diseases and disorders, different types of genes, and known drug targets represented on Affymetrix chicken genome array. To facilitate functional annotation of other arrays and microarray experimental datasets we developed an Array GO Mapper (AGOM) tool to help researchers to quickly retrieve corresponding functional information for their dataset. CONCLUSION: Results from this study will directly facilitate annotation of other chicken arrays and microarray experimental datasets. Researchers will be able to quickly model their microarray dataset into more reliable biological functional information by using AGOM tool. The disease, disorders, gene types and drug targets revealed in the study will allow researchers to learn more about how genes function in complex biological systems and may lead to new drug discovery and development of therapies. The GO annotation data generated will be available for public use via AgBase website and will be updated on regular basis.

Gene Ontology annotation quality analysis in model eukaryotes.

Functional analysis using the Gene Ontology (GO) is crucial for array analysis, but it is often difficult for researchers to assess the amount and quality of GO annotations associated with different sets of gene products. In many cases the source of the GO annotations and the date the GO annotations were last updated is not apparent, further complicating a researchers' ability to assess the quality of the GO data provided. Moreover, GO biocurators need to ensure that the GO quality is maintained and optimal for the functional processes that are most relevant for their research community. We report the GO Annotation Quality (GAQ) score, a quantitative measure of GO quality that includes breadth of GO annotation, the level of detail of annotation and the type of evidence used to make the annotation. As a case study, we apply the GAQ scoring method to a set of diverse eukaryotes and demonstrate how the GAQ score can be used to track changes in GO annotations over time and to assess the quality of GO annotations available for specific biological processes. The GAQ score also allows researchers to quantitatively assess the functional data available for their experimental systems (arrays or databases).

Experimental-confirmation and functional-annotation of predicted proteins in the chicken genome.

BACKGROUND: The chicken genome was sequenced because of its phylogenetic position as a non-mammalian vertebrate, its use as a biomedical model especially to study embryology and development, its role as a source of human disease organisms and its importance as the major source of animal derived food protein. However, genomic sequence data is, in itself, of limited value; generally it is not equivalent to understanding biological function. The benefit of having a genome sequence is that it provides a basis for functional genomics. However, the sequence data currently available is poorly structurally and functionally annotated and many genes do not have standard nomenclature assigned. RESULTS: We analysed eight chicken tissues and improved the chicken genome structural annotation by providing experimental support for the in vivo expression of 7,809 computationally predicted proteins, including 30 chicken proteins that were only electronically predicted or hypothetical translations in human. To improve functional annotation (based on Gene Ontology), we mapped these identified proteins to their human and mouse orthologs and used this orthology to transfer Gene Ontology (GO) functional annotations to the chicken proteins. The 8,213 orthology-based GO annotations that we produced represent an 8% increase in currently available chicken GO annotations. Orthologous chicken products were also assigned standardized nomenclature based on current chicken nomenclature guidelines. CONCLUSION: We demonstrate the utility of high-throughput expression proteomics for rapid experimental structural annotation of a newly sequenced eukaryote genome. These experimentally-supported predicted proteins were further annotated by assigning the proteins with standardized nomenclature and functional annotation. This method is widely applicable to a diverse range of species. Moreover, information from one genome can be used to improve the annotation of other genomes and inform gene prediction algorithms.