Robert W. Li

Analysis of copy number variations among diverse cattle breeds

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. Here, we describe the first systematic and genome-wide analysis of copy number variations (CNVs) in modern domesticated cattle using array comparative genomic hybridization (array CGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH). The array CGH panel included 90 animals from 11 Bos taurus, three Bos indicus, and three composite breeds for beef, dairy, or dual purpose. We identified over 200 candidate CNV regions (CNVRs) in total and 177 within known chromosomes, which harbor or are adjacent to gains or losses. These 177 high-confidence CNVRs cover 28.1 megabases or approximately 1.07% of the genome. Over 50% of the CNVRs (89/177) were found in multiple animals or breeds and analysis revealed breed-specific frequency differences and reflected aspects of the known ancestry of these cattle breeds. Selected CNVs were further validated by independent methods using qPCR and FISH. Approximately 67% of the CNVRs (119/177) completely or partially span cattle genes and 61% of the CNVRs (108/177) directly overlap with segmental duplications. The CNVRs span about 400 annotated cattle genes that are significantly enriched for specific biological functions, such as immunity, lactation, reproduction, and rumination. Multiple gene families, including ULBP, have gone through ruminant lineage-specific gene amplification. We detected and confirmed marked differences in their CNV frequencies across diverse breeds, indicating that some cattle CNVs are likely to arise independently in breeds and contribute to breed differences. Our results provide a valuable resource beyond microsatellites and single nucleotide polymorphisms to explore the full dimension of genetic variability for future cattle genomic research.

Characterization of the rumen microbiota of pre-ruminant calves using metagenomic tools

The temporal sequence of microbial establishment in the rumen of the neonatal ruminant has important ecological and pathophysiological implications. In this study, we characterized the rumen microbiota of pre-ruminant calves fed milk replacer using two approaches, pyrosequencing of hypervariable V3-V5 regions of the 16S rRNA gene and whole-genome shotgun approach. Fifteen bacterial phyla were identified in the microbiota of pre-ruminant calves. Bacteroidetes was the predominant phylum in the rumen microbiota of 42-day-old calves, representing 74.8% of the 16S sequences, followed by Firmicutes (12.0%), Proteobacteria (10.4%), Verrucomicrobia (1.2%) and Synergistetes (1.1%). However, the phylum-level composition of 14-day-old calves was distinctly different. A total of 170 bacterial genera were identified while the core microbiome of pre-ruminant calves included 45 genera. Rumen development seemingly had a significant impact on microbial diversity. The dazzling functional diversity of the rumen microbiota was reflected by identification of 8298 Pfam and 3670 COG protein families. The rumen microbiota of pre-ruminant calves displayed a considerable compositional heterogeneity during early development. This is evidenced by a profound difference in rumen microbial composition between the two age groups. However, all functional classes between the two age groups had a remarkably similar assignment, suggesting that rumen microbial communities of pre-ruminant calves maintained a stable function and metabolic potentials while their phylogenetic composition fluctuated greatly. The presence of all major types of rumen microorganisms suggests that the rumen of pre-ruminant calves may not be rudimentary. Our results provide insight into rumen microbiota dynamics and will facilitate efforts in formulating optimal early-weaning strategies.

Alterations in the Porcine Colon Microbiota Induced by the Gastrointestinal Nematode Trichuris suis

ABSTRACT Helminth parasites ensure their survival by regulating host immunity through mechanisms that dampen inflammation. These properties have recently been exploited therapeutically to treat human diseases. The biocomplexity of the intestinal lumen suggests that interactions between the parasite and the intestinal microbiota would also influence inflammation. In this study, we characterized the microbiota in the porcine proximal colon in response to Trichuris suis (whipworm) infection using 16S rRNA gene-based and whole-genome shotgun (WGS) sequencing. A 21-day T. suis infection in four pigs induced a significant change in the composition of the proximal colon microbiota compared to that of three parasite-naive pigs. Among the 15 phyla identified, the abundances of Proteobacteria and Deferribacteres were changed in infected pigs. The abundances of approximately 13% of genera were significantly altered by infection. Changes in relative abundances of Succinivibrio and Mucispirillum, for example, may relate to alterations in carbohydrate metabolism and niche disruptions in mucosal interfaces induced by parasitic infection, respectively. Of note, infection by T. suis led to a significant shift in the metabolic potential of the proximal colon microbiota, where 26% of all metabolic pathways identified were affected. Besides carbohydrate metabolism, lysine biosynthesis was repressed as well. A metabolomic analysis of volatile organic compounds (VOCs) in the luminal contents showed a relative absence in infected pigs of cofactors for carbohydrate and lysine biosynthesis, as well as an accumulation of oleic acid, suggesting altered fatty acid absorption contributing to local inflammation. Our findings should facilitate development of strategies for parasitic control in pigs and humans.

Butyrate induces profound changes in gene expression related to multiple signal pathways in bovine kidney epithelial cells

BackgroundGlobal gene expression profiles of bovine kidney epithelial cells regulated by sodium butyrate were investigated with high-density oligonucleotide microarrays. The bovine microarray with 86,191 distinct 60mer oligonucleotides, each with 4 replicates, was designed and produced with Maskless Array Synthesizer technology. These oligonucleotides represent approximately 45,383 unique cattle sequences.Results450 genes significantly regulated by butyrate with a median False Discovery Rate (FDR) = 0 % were identified. The majority of these genes were repressed by butyrate and associated with cell cycle control. The expression levels of 30 selected genes identified by the microarray were confirmed using real-time PCR. The results from real-time PCR positively correlated (R = 0.867) with the results from the microarray.ConclusionThis study presented the genes related to multiple signal pathways such as cell cycle control and apoptosis. The profound changes in gene expression elucidate the molecular basis for the pleiotropic effects of butyrate on biological processes. These findings enable better recognition of the full range of beneficial roles butyrate may play during cattle energy metabolism, cell growth and proliferation, and possibly in fighting gastrointestinal pathogens.

Worm Burden-Dependent Disruption of the Porcine Colon Microbiota by Trichuris suis Infection

Helminth infection in pigs serves as an excellent model for the study of the interaction between human malnutrition and parasitic infection and could have important implications in human health. We had observed that pigs infected with Trichuris suis for 21 days showed significant changes in the proximal colon microbiota. In this study, interactions between worm burden and severity of disruptions to the microbial composition and metabolic potentials in the porcine proximal colon microbiota were investigated using metagenomic tools. Pigs were infected by a single dose of T. suis eggs for 53 days. Among infected pigs, two cohorts were differentiated that either had adult worms or were worm-free. Infection resulted in a significant change in the abundance of approximately 13% of genera detected in the proximal colon microbiota regardless of worm status, suggesting a relatively persistent change over time in the microbiota due to the initial infection. A significant reduction in the abundance of Fibrobacter and Ruminococcus indicated a change in the fibrolytic capacity of the colon microbiota in T. suis infected pigs. In addition, ∼10% of identified KEGG pathways were affected by infection, including ABC transporters, peptidoglycan biosynthesis, and lipopolysaccharide biosynthesis as well as α-linolenic acid metabolism. Trichuris suis infection modulated host immunity to Campylobacter because there was a 3-fold increase in the relative abundance in the colon microbiota of infected pigs with worms compared to naïve controls, but a 3-fold reduction in worm-free infected pigs compared to controls. The level of pathology observed in infected pigs with worms compared to worm-free infected pigs may relate to the local host response because expression of several Th2-related genes were enhanced in infected pigs with worms versus those worm-free. Our findings provided insight into the dynamics of the proximal colon microbiota in pigs in response to T. suis infection.

Perturbation dynamics of the rumen microbiota in response to exogenous butyrate.

The capacity of the rumen microbiota to produce volatile fatty acids (VFAs) has important implications in animal well-being and production. We investigated temporal changes of the rumen microbiota in response to butyrate infusion using pyrosequencing of the 16S rRNA gene. Twenty one phyla were identified in the rumen microbiota of dairy cows. The rumen microbiota harbored 54.5±6.1 genera (mean ± SD) and 127.3±4.4 operational taxonomic units (OTUs), respectively. However, the core microbiome comprised of 26 genera and 82 OTUs. Butyrate infusion altered molar percentages of 3 major VFAs. Butyrate perturbation had a profound impact on the rumen microbial composition. A 72 h-infusion led to a significant change in the numbers of sequence reads derived from 4 phyla, including 2 most abundant phyla, Bacteroidetes and Firmicutes. As many as 19 genera and 43 OTUs were significantly impacted by butyrate infusion. Elevated butyrate levels in the rumen seemingly had a stimulating effect on butyrate-producing bacteria populations. The resilience of the rumen microbial ecosystem was evident as the abundance of the microorganisms returned to their pre-disturbed status after infusion withdrawal. Our findings provide insight into perturbation dynamics of the rumen microbial ecosystem and should guide efforts in formulating optimal uses of probiotic bacteria treating human diseases.

Gene expression in the digestive tissues of ruminants and their relationships with feeding and digestive processes

The gastrointestinal tract (GIT) has multiple functions including digestion, nutrient absorption, secretion of hormones and excretion of wastes. In the ruminant animal, development of this organ system is more complex than that of the monogastric animal due to the necessity to establish a fully functional and differentiated rumen, in which a diverse microbial population of bacteria, fungi and protozoa support fermentation and digestion of dietary fiber. Central to the goal of animal scientists to enhance nutrient uptake and production efficiency of ruminants is the need for a comprehensive understanding of GIT development, as well as conditions that alter the digestion process. The relatively recent availability of genome sequence information has permitted physiological investigations related to the process of digestion for many agriculturally important species at the gene transcript level. For instance, numerous studies have evaluated the expression of ruminant GIT genes to gain insight into mechanisms involved in normal function, physiology and development, such as nutrient uptake and transport across the epithelial cell barrier throughout the alimentary canal, maintenance of rumen pH, and regulation of GIT motility and cell proliferation. Further, multiple studies have examined the effects of dietary modification, including feeding of supplemental fat, starch and protein, or a forage- v. concentrate-based diet on expression of critical gene pathways in the gut. In addition, the expression of genes in the GIT in response to disease, such as infection with gastrointestinal parasites, has been investigated. This review will summarize some of the recent scientific literature related to the gene expression in the GIT of ruminants, primarily cattle, sheep and goats, as it pertains to normal physiology, and dietary, developmental, and disease effects to provide an overview of critical proteins participating in the overall digestive processes, and their physiological functions. Recent findings from our laboratory will be highlighted also related to expression of the glucagon-like peptide two-hormone pathway in the GIT of dairy cattle during in various stages of the development and lactation, alterations in gene pathways associated with the rumen development and differentiation in the weaning calf, and genes of the GIT responding to Ostertagia, a common nematode infection of the cattle. Finally, prospective areas of investigation will be highlighted.

Pathway analysis identifies perturbation of genetic networks induced by butyrate in a bovine kidney epithelial cell line

Ruminant species have evolved to metabolize the short-chain volatile fatty acids (VFA), acetate, propionate, and butyrate, to fulfill up to 70% of their nutrient energy requirements. The inherent VFA dependence of ruminant cells was exploited to add a level of increased sensitivity to the study of the role of butyrate gene-response elements in regulatory biochemical pathways. Global gene expression profiles of the bovine kidney epithelial cells regulated by sodium butyrate were investigated with high-density oligonucleotide microarrays. The detailed mechanisms by which butyrate induces cell growth arrest and apoptosis were analyzed using the Ingenuity Pathways Knowledge Base. The functional category and pathway analyses of the microarray data revealed that four canonical pathways (Cell cycles: G2/M DNA damage checkpoint, and pyrimidine metabolism; G1/S checkpoint regulation and purine metabolism) were significantly perturbed. The biologically relevant networks and pathways of these genes were also identified. IGF2, TGFB1, TP53, E2F4, and CDC2 were established as being centered in these genomic networks. The present findings provide a basis for understanding the full range of the biological roles and the molecular mechanisms that butyrate may play in animal cell growth, proliferation, and energy metabolisms.

Quantification of Transcriptome Responses of the Rumen Epithelium to Butyrate Infusion using RNA-seq Technology

Short-chain fatty acids (SCFAs), such as butyrate, produced by gut microorganisms, play a critical role in energy metabolism and physiology of ruminants as well as in human health. In this study, the temporal effect of elevated butyrate concentrations on the transcriptome of the rumen epithelium was quantified via serial biopsy sampling using RNA-seq technology. The mean number of genes transcribed in the rumen epithelial transcriptome was 17,323.63 ± 277.20 (±SD; N = 24) while the core transcriptome consisted of 15,025 genes. Collectively, 80 genes were identified as being significantly impacted by butyrate infusion across all time points sampled. Maximal transcriptional effect of butyrate on the rumen epithelium was observed at the 72-h infusion when the abundance of 58 genes was altered. The initial reaction of the rumen epithelium to elevated exogenous butyrate may represent a stress response as Gene Ontology (GO) terms identified were predominantly related to responses to bacteria and biotic stimuli. An algorithm for the reconstruction of accurate cellular networks (ARACNE) inferred regulatory gene networks with 113,738 direct interactions in the butyrate-epithelium interactome using a combined cutoff of an error tolerance (ɛ = 0.10) and a stringent P-value threshold of mutual information (5.0 × 10−11). Several regulatory networks were controlled by transcription factors, such as CREBBP and TTF2, which were regulated by butyrate. Our findings provide insight into the regulation of butyrate transport and metabolism in the rumen epithelium, which will guide our future efforts in exploiting potential beneficial effect of butyrate in animal well-being and human health.

Detection of germline and somatic copy number variations in cattle.

As a complement to the Bovine HapMap Consortium project, we initiated a systematic study of the copy numbervariation (CNV) within the same cattle population using array comparative genomic hybridization (array CGH). Oligonucleotide CGH arrays were designed and fabricated to cover all chromosomes with an average interval of 6 kb using the latest bovine genome assembly. In the initial screening, three Holstein bulls were selected to represent major paternal lineages of the Holstein breed with some maternal linkages between these lines. Dual-label hybridizations were performed using either Hereford L1 Dominette 01449 or L1 Domino 99375 as reference. The CNVs were represented by gains and losses of normalized fluorescence intensities relative to the reference. The data presented here, for the first time, demonstrated that significant amounts of germline and fewer somatic CNVs exist in cattle, that many CNVs are common both across diverse cattle breeds and among individuals within a breed, and that array CGH is an effective tool to systematically detect bovine CNV. Selected CNVs have been confirmed by independent methods using real-time (RT) PCR. The strategy used in this study, based on genome higher-orderarchitecture variation, is a powerful approach to generating resources for the identification of novel genomic variation and candidate genes for economically important traits.