Robert W. Reid

Association Between Composition of the Human Gastrointestinal Microbiome and Development of Fatty Liver With Choline Deficiency

BACKGROUND & AIMS Nonalcoholic fatty liver disease affects up to 30% of the US population, but the mechanisms underlying this condition are incompletely understood. We investigated how diet standardization and choline deficiency influence the composition of the microbial community in the human gastrointestinal tract and the development of fatty liver under conditions of choline deficiency. METHODS We performed a 2-month inpatient study of 15 female subjects who were placed on well-controlled diets in which choline levels were manipulated. We used 454-FLX pyrosequencing of 16S ribosomal RNA bacterial genes to characterize microbiota in stool samples collected over the course of the study. RESULTS The compositions of the gastrointestinal microbial communities changed with choline levels of diets; each individuals microbiome remained distinct for the duration of the experiment, even though all subjects were fed identical diets. Variations between subjects in levels of Gammaproteobacteria and Erysipelotrichi were directly associated with changes in liver fat in each subject during choline depletion. Levels of these bacteria, change in amount of liver fat, and a single nucleotide polymorphism that affects choline were combined into a model that accurately predicted the degree to which subjects developed fatty liver on a choline-deficient diet. CONCLUSIONS Host factors and gastrointestinal bacteria each respond to dietary choline deficiency, although the gut microbiota remains distinct in each individual. We identified bacterial biomarkers of fatty liver that result from choline deficiency, adding to the accumulating evidence that gastrointestinal microbes have a role in metabolic disorders.

RNA-Seq analysis and annotation of a draft blueberry genome assembly identifies candidate genes involved in fruit ripening, biosynthesis of bioactive compounds, and stage-specific alternative splicing.

BackgroundBlueberries are a rich source of antioxidants and other beneficial compounds that can protect against disease. Identifying genes involved in synthesis of bioactive compounds could enable the breeding of berry varieties with enhanced health benefits.ResultsToward this end, we annotated a previously sequenced draft blueberry genome assembly using RNA-Seq data from five stages of berry fruit development and ripening. Genome-guided assembly of RNA-Seq read alignments combined with output from ab initio gene finders produced around 60,000 gene models, of which more than half were similar to proteins from other species, typically the grape Vitis vinifera. Comparison of gene models to the PlantCyc database of metabolic pathway enzymes identified candidate genes involved in synthesis of bioactive compounds, including bixin, an apocarotenoid with potential disease-fighting properties, and defense-related cyanogenic glycosides, which are toxic. Cyanogenic glycoside (CG) biosynthetic enzymes were highly expressed in green fruit, and a candidate CG detoxification enzyme was up-regulated during fruit ripening. Candidate genes for ethylene, anthocyanin, and 400 other biosynthetic pathways were also identified. Homology-based annotation using Blast2GO and InterPro assigned Gene Ontology terms to around 15,000 genes. RNA-Seq expression profiling showed that blueberry growth, maturation, and ripening involve dynamic gene expression changes, including coordinated up- and down-regulation of metabolic pathway enzymes and transcriptional regulators. Analysis of RNA-seq alignments identified developmentally regulated alternative splicing, promoter use, and 3′ end formation.ConclusionsWe report genome sequence, gene models, functional annotations, and RNA-Seq expression data that provide an important new resource enabling high throughput studies in blueberry.

Early osteoblast responses to orthopedic implants: Synergy of surface roughness and chemistry of bioactive ceramic coating

Pro-osteogenic stimulation of bone cells by bioactive ceramic-coated orthopedic implants is influenced by both surface roughness and material chemistry; however, their concomitant impact on osteoblast behavior is not well understood. The aim of this study is to investigate the effects of nano-scale roughness and chemistry of bioactive silica-calcium phosphate nanocomposite (SCPC50) coated Ti-6Al-4V on modulating early bone cell responses. Cell attachment was higher on SCPC50-coated substrates compared to the uncoated controls; however, cells on the uncoated substrate exhibited greater spreading and superior quality of F-actin filaments than cells on the SCPC50-coated substrates. The poor F-actin filament organization on SCPC50-coated substrates is thought to be due to the enhanced calcium uptake by the ceramic surface. Dissolution analyses showed that an increase in surface roughness was accompanied by increased calcium uptake, and increased phosphorous and silicon release, all of which appear to interfere with F-actin assembly and osteoblast morphology. Moreover, cell attachment onto the SCPC50-coated substrates correlated with the known adsorption of fibronectin, and was independent of surface roughness. High-throughput genome sequencing showed enhanced expression of extracellular matrix and cell differentiation related genes. These results demonstrate a synergistic relationship between bioactive ceramic coating roughness and material chemistry resulting in a phenotype that leads to early osteoblast differentiation.

EchinoDB, an application for comparative transcriptomics of deeply-sampled clades of echinoderms

BackgroundOne of our goals for the echinoderm tree of life project (http://echinotol.org) is to identify orthologs suitable for phylogenetic analysis from next-generation transcriptome data. The current dataset is the largest assembled for echinoderm phylogeny and transcriptomics. We used RNA-Seq to profile adult tissues from 42 echinoderm specimens from 24 orders and 37 families. In order to achieve sampling members of clades that span key evolutionary divergence, many of our exemplars were collected from deep and polar seas.DescriptionA small fraction of the transcriptome data we produced is being used for phylogenetic reconstruction. Thus to make a larger dataset available to researchers with a wide variety of interests, we made a web-based application, EchinoDB (http://echinodb.uncc.edu). EchinoDB is a repository of orthologous transcripts from echinoderms that is searchable via keywords and sequence similarity.ConclusionsFrom transcripts we identified 749,397 clusters of orthologous loci. We have developed the information technology to manage and search the loci their annotations with respect to the Sea Urchin (Strongylocentrotus purpuratus) genome. Several users have already taken advantage of these data for spin-off projects in developmental biology, gene family studies, and neuroscience. We hope others will search EchinoDB to discover datasets relevant to a variety of additional questions in comparative biology.

Determining gene expression on a single pair of microarrays

BackgroundIn microarray experiments the numbers of replicates are often limited due to factors such as cost, availability of sample or poor hybridization. There are currently few choices for the analysis of a pair of microarrays where N = 1 in each condition. In this paper, we demonstrate the effectiveness of a new algorithm called PINC (PINC is Not Cyber-T) that can analyze Affymetrix microarray experiments.ResultsPINC treats each pair of probes within a probeset as an independent measure of gene expression using the Bayesian framework of the Cyber-T algorithm and then assigns a corrected p-value for each gene comparison.The p-values generated by PINC accurately control False Discovery rate on Affymetrix control data sets, but are small enough that family-wise error rates (such as the Holms step down method) can be used as a conservative alternative to false discovery rate with little loss of sensitivity on control data sets.ConclusionPINC outperforms previously published methods for determining differentially expressed genes when comparing Affymetrix microarrays with N = 1 in each condition. When applied to biological samples, PINC can be used to assess the degree of variability observed among biological replicates in addition to analyzing isolated pairs of microarrays.

The phylogeny of extant starfish (Asteroidea: Echinodermata) including Xyloplax, based on comparative transcriptomics

Multi-locus phylogenetic studies of echinoderms based on Sanger and RNA-seq technologies and the fossil record have provided evidence for the Asterozoa-Echinozoa hypothesis. This hypothesis posits a sister relationship between asterozoan classes (Asteroidea and Ophiuroidea) and a similar relationship between echinozoan classes (Echinoidea and Holothuroidea). Despite this consensus around Asterozoa-Echinozoa, phylogenetic relationships within the class Asteroidea (sea stars or starfish) have been controversial for over a century. Open questions include relationships within asteroids and the status of the enigmatic taxon Xyloplax. Xyloplax is thought by some to represent a newly discovered sixth class of echinoderms - and by others to be an asteroid. To address these questions, we applied a novel workflow to a large RNA-seq dataset that encompassed a broad taxonomic and genomic sample. This study included 15 species sampled from all extant orders and 13 families, plus four ophiuroid species as an outgroup. To expand the taxonomic coverage, the study also incorporated five previously published transcriptomes and one previously published expressed sequence tags (EST) dataset. We developed and applied methods that used a range of alignment parameters with increasing permissiveness in terms of gap characters present within an alignment. This procedure facilitated the selection of phylogenomic data subsets from large amounts of transcriptome data. The results included 19 nested data subsets that ranged from 37 to 4,281loci. Tree searches on all data subsets reconstructed Xyloplax as a velatid asteroid rather than a new class. This result implies that asteroid morphology remains labile well beyond the establishment of the body plan of the group. In the phylogenetic tree with the highest average asteroid nodal support several monophyletic groups were recovered. In this tree, Forcipulatida and Velatida are monophyletic and form a clade that includes Brisingida as sister to Forcipulatida. Xyloplax is consistently recovered as sister to Pteraster. Paxillosida and Spinulosida are each monophyletic, with Notomyotida as sister to the Paxillosida. Valvatida is recovered as paraphyletic. The results from other data subsets are largely consistent with these results. Our results support the hypothesis that the earliest divergence event among extant asteroids separated Velatida and Forcipulatacea from Valvatacea and Spinulosida.

Phylotranscriptomic analysis uncovers a wealth of tissue inhibitor of metalloproteinases variants in echinoderms.

Tissue inhibitors of metalloproteinases (TIMPs) help regulate the extracellular matrix (ECM) in animals, mostly by inhibiting matrix metalloproteinases (MMPs). They are important activators of mutable collagenous tissue (MCT), which have been extensively studied in echinoderms, and the four TIMP copies in humans have been studied for their role in cancer. To understand the evolution of TIMPs, we combined 405 TIMPs from an echinoderm transcriptome dataset built from 41 specimens representing all five classes of echinoderms with variants from protostomes and chordates. We used multiple sequence alignment with various stringencies of alignment quality to cull highly divergent sequences and then conducted phylogenetic analyses using both nucleotide and amino acid sequences. Phylogenetic hypotheses consistently recovered TIMPs as diversifying in the ancestral deuterostome and these early lineages continuing to diversify in echinoderms. The four vertebrate TIMPs diversified from a single copy in the ancestral chordate, all other copies being lost. Consistent with greater MCT needs owing to body wall liquefaction, evisceration, autotomy and reproduction by fission, holothuroids had significantly more TIMPs and higher read depths per contig. Ten cysteine residues, an HPQ binding site and several other residues were conserved in at least 70% of all TIMPs. The conservation of binding sites and the placement of echinoderm TIMPs involved in MCT modification suggest that ECM regulation remains the primary function of TIMP genes, although within this role there are a large number of specialized copies.

Students' perspective on genomics: from sample to sequence using the case study of blueberry.

Advances in genomic sequencing technologies in the past decade have revolutionized the field of genomics, resulting in faster and less expensive sequencing. Holding back the potential for innovation, however, is a widespread lack of understanding of genomics and sequencing by the general public. In an attempt to remedy this problem, this paper presents an introduction to the fields of genomics, bioinformatics, and proteomics using the blueberry genome as a model case study of the plant genomics field. The blueberry (Vaccinium sect. Cyanococcus) is often cited as a “super food” in the media due to its nutritional benefits and global economic importance. There have been a number of related genomic publications in the past 20 years; however, a completed genome and a full analysis into the health-related pathways are still needed. As exemplified by this blueberry case study, there are opportunities for future genomic research into numerous beneficial plant species. The solid background presented in this paper provides future researchers the foundation to explore these uncharted areas.

LinkageMapView—rendering high-resolution linkage and QTL maps

Abstract Motivation Linkage and quantitative trait loci (QTL) maps are critical tools for the study of the genetic basis of complex traits. With the advances in sequencing technology over the past decade, linkage map densities have been increasing dramatically, while the visualization tools have not kept pace. LinkageMapView is a free add-on package written in R that produces high resolution, publication-ready visualizations of linkage and QTL maps. While there is software available to generate linkage map graphics, none are freely available, produce publication quality figures, are open source and can run on all platforms. LinkageMapView can be integrated into map building pipelines as it seamlessly incorporates output from R/qtl and also accepts simple text or comma delimited files. There are numerous options within the package to build highly customizable maps, allow for linkage group comparisons, and annotate QTL regions. Availability and implementation https://cran.r-project.org/web/packages/LinkageMapView/

Salt tolerance response revealed by RNA-Seq in a diploid halophytic wild relative of sweet potato.

Crop wild relatives harbor exotic and novel genetic resources, which hold great potential for crop improvement. Ipomoea imperati is a wild diploid relative of sweet potato with the capability of high salinity tolerance. We compared the transcriptomes of I. imperati under salt stress vs. control to identify candidate genes and pathways involved in salt response. De novo assembly produced 67,911 transcripts with a high depth of coverage. A total of 39,902 putative genes were assigned annotations, and 936 and 220 genes involved in salt response in roots and leaves, respectively. Functional analysis indicated a whole system response during salt stress in I. imperati, which included four metabolic processes: sensory initiation, transcriptional reprogramming, cellular protein component change, and cellular homeostasis regulation. We identified a number of candidate genes involved in the ABA signaling pathway, as well as transcription factors, transporters, antioxidant enzymes, and enzymes associated with metabolism of synthesis and catalysis. Furthermore, two membrane transporter genes, including vacuole cation/proton exchanger and inositol transporter, were considered to play important roles in salt tolerance. This study provided valuable information not only for understanding the genetic basis of ecological adaptation but also for future application in sweet potato and other crop improvements.