Animal & Avian Sciences Research Works

Permanent URI for this collectionhttp://hdl.handle.net/1903/1600

Browse

Now showing 1 - 20 of 89

CeHR Axenic Liquid Medium for C. elegans growth
(2005-02-24T14:24:16Z) Hamza, Iqbal
CeHR medium was originally formulated by Dr. Eric Clegg at the United States Army Center for Environmental Health Research to grow worms in an axenic liquid media. (http://usacehr.detrick.army.mil/CleggLab.html) We are thankful to Eric for his invaluable help and advice through all the modifications we have since made to optimize CeHR medium to suit our research goals. We modified CeHR medium (mCeHR) because of our interests in defining the individual micronutrient components in order to identify the cellular and molecular components of metal homeostasis using C. elegans as a model multicellular eukaryote.
Identification of QTL controlling meat quality traits in an F2 cross between two chicken lines selected for either low or high growth rate
(Springer Nature, 2007-06-08) Nadaf, Javad; Gilbert, Hélène; Pitel, Frédérique; Berri, Cécile M; Feve, Katia; Beaumont, Catherine; Duclos, Michel J; Vignal, Alain; Porter, Tom E; Simon, Jean; Aggrey, Samuel E; Cogburn, Larry A; Le Bihan-Duval, Elisabeth
Meat technological traits (i.e. meat pH, water retention and color) are important considerations for improving further processing of chicken meat. These quality traits were originally characterized in experimental lines selected for high (HG) and low (LG) growth. Presently, quantitative trait loci (QTL) for these traits were analyzed in an F2 population issued from the HG × LG cross. A total of 698 animals in 50 full-sib families were genotyped for 108 microsatellite markers covering 21 linkage groups. The HG and LG birds exhibit large differences in body weight and abdominal fat content. Several meat quality traits [pH at 15 min post-slaughter (pH15) and ultimate pH (pHu), breast color-redness (BCo-R) and breast color-yellowness (BCo-Y)] were lower in HG chickens. In contrast, meat color-lightness (BCo-L) was higher in HG chickens, whereas meat drip loss (DL) was similar in both lines. HG birds were more active on the shackle line. Association analyses were performed using maximum-likelihood interval mapping in QTLMAP. Five genome-wide significant QTLs were revealed: two for pH15 on GGA1 and GGA2, one for DL on GGA1, one for BCo-R and one for BCo-Y both on GGA11. In addition, four suggestive QTLs were identified by QTLMAP for BCo-Y, pHu, pH15 and DL on GGA1, GGA4, GGA12 and GGA14, respectively. The QTL effects, averaged on heterozygous families, ranged from 12 to 31% of the phenotypic variance. Further analyses with QTLExpress confirmed the two genome-wide QTLs for meat color on GGA11, failed to identify the genome-wide QTL for pH15 on GGA2, and revealed only suggestive QTLs for pH15 and DL on GGA1. However, QTLExpress qualified the QTL for pHu on GGA4 as genome-wide. The present study identified genome-wide significant QTLs for all meat technological traits presently assessed in these chickens, except for meat lightness. This study highlights the effects of divergent selection for growth rate on some behavioral traits, muscle biochemistry and ultimately meat quality traits. Several QTL regions were identified that are worthy of further characterization. Some QTLs may in fact co-localize, suggesting pleiotropic effects for some chromosomal regions.
Term-tissue specific models for prediction of gene ontology biological processes using transcriptional profiles of aging in drosophila melanogaster
(Springer Nature, 2008-02-28) Zhang, Wensheng; Zou, Sige; Song, Jiuzhou
Predictive classification on the base of gene expression profiles appeared recently as an attractive strategy for identifying the biological functions of genes. Gene Ontology (GO) provides a valuable source of knowledge for model training and validation. The increasing collection of microarray data represents a valuable source for generating functional hypotheses of uncharacterized genes. This study focused on using support vector machines (SVM) to predict GO biological processes from individual or multiple-tissue transcriptional profiles of aging in Drosophila melanogaster. Ten-fold cross validation was implemented to evaluate the prediction. One-tail Fisher's exact test was conducted on each cross validation and multiple testing was addressed using BH FDR procedure. The results showed that, of the 148 pursued GO biological processes, fifteen terms each had at least one model with FDR-adjusted p-value (Adj.p) <0.05 and six had the values between 0.05 and 0.25. Furthermore, all these models had the prediction sensitivity (SN) over 30% and specificity (SP) over 80%. We proposed the concept of term-tissue specific models indicating the fact that the major part of the optimized prediction models was trained from individual tissue data. Furthermore, we observed that the memberships of the genes involved in all the three pursued children biological processes on mitochondrial electron transport could be predicted from the transcriptional profiles of aging (Adj.p < 0.01). This finding may be important in biology because the genes of mitochondria play a critical role in the longevity of C. elegans and D. melanogaster.
Principal component tests: applied to temporal gene expression data
(Springer Nature, 2009-01-30) Zhang, Wensheng; Fang, Hong-Bin; Song, Jiuzhou
Clustering analysis is a common statistical tool for knowledge discovery. It is mainly conducted when a project still is in the exploratory phase without any priori hypotheses. However, the statistical significance testing between the clusters can be meaningful in helping the researchers to assess if the classification results from implementing a clustering algorithm need to be improved, even after the cluster number has been determined by a well-established criterion. This is important when we want to identify highly-specific patterns through classification. We proposed to use a principal component (PC) test, which is an implementation of an exact F statistic for the measures at multiple endpoints based on elliptical distribution theory, to assess the statistical significance between clusters. A challenge in the implementation is the choice of the number (q) of principal components to be considered, which can severely influence the statistical power of the method. We optimized the determination via validation according to a permutation test based on the clustering to be evaluated. The method was applied to a public dataset in classifying genes according to their temporal gene expression profiles. The results demonstrated that the PC testing were useful for determining the optimal number of clusters.
Gene expression responses in male fathead minnows exposed to binary mixtures of an estrogen and antiestrogen
(Springer Nature, 2009-07-13) Garcia-Reyero, Natàlia; Kroll, Kevin J; Liu, Li; Orlando, Edward F; Watanabe, Karen H; Sepúlveda, María S; Villeneuve, Daniel L; Perkins, Edward J; Ankley, Gerald T; Denslow, Nancy D
Aquatic organisms are continuously exposed to complex mixtures of chemicals, many of which can interfere with their endocrine system, resulting in impaired reproduction, development or survival, among others. In order to analyze the effects and mechanisms of action of estrogen/anti-estrogen mixtures, we exposed male fathead minnows (Pimephales promelas) for 48 hours via the water to 2, 5, 10, and 50 ng 17α-ethinylestradiol (EE2)/L, 100 ng ZM 189,154/L (a potent antiestrogen known to block activity of estrogen receptors) or mixtures of 5 or 50 ng EE2/L with 100 ng ZM 189,154/L. We analyzed gene expression changes in the gonad, as well as hormone and vitellogenin plasma levels. Steroidogenesis was down-regulated by EE2 as reflected by the reduced plasma levels of testosterone in the exposed fish and down-regulation of genes in the steroidogenic pathway. Microarray analysis of testis of fathead minnows treated with 5 ng EE2/L or with the mixture of 5 ng EE2/L and 100 ng ZM 189,154/L indicated that some of the genes whose expression was changed by EE2 were blocked by ZM 189,154, while others were either not blocked or enhanced by the mixture, generating two distinct expression patterns. Gene ontology and pathway analysis programs were used to determine categories of genes for each expression pattern. Our results suggest that response to estrogens occurs via multiple mechanisms, including canonical binding to soluble estrogen receptors, membrane estrogen receptors, and other mechanisms that are not blocked by pure antiestrogens.
Analysis of recent segmental duplications in the bovine genome
(Springer Nature, 2009-12-01) Liu, George E; Ventura, Mario; Cellamare, Angelo; Chen, Lin; Cheng, Ze; Zhu, Bin; Li, Congjun; Song, Jiuzhou; Eichler, Evan E
Duplicated sequences are an important source of gene innovation and structural variation within mammalian genomes. We performed the first systematic and genome-wide analysis of segmental duplications in the modern domesticated cattle (Bos taurus). Using two distinct computational analyses, we estimated that 3.1% (94.4 Mb) of the bovine genome consists of recently duplicated sequences (≥ 1 kb in length, ≥ 90% sequence identity). Similar to other mammalian draft assemblies, almost half (47% of 94.4 Mb) of these sequences have not been assigned to cattle chromosomes. In this study, we provide the first experimental validation large duplications and briefly compared their distribution on two independent bovine genome assemblies using fluorescent in situ hybridization (FISH). Our analyses suggest that the (75-90%) of segmental duplications are organized into local tandem duplication clusters. Along with rodents and carnivores, these results now confidently establish tandem duplications as the most likely mammalian archetypical organization, in contrast to humans and great ape species which show a preponderance of interspersed duplications. A cross-species survey of duplicated genes and gene families indicated that duplication, positive selection and gene conversion have shaped primates, rodents, carnivores and ruminants to different degrees for their speciation and adaptation. We identified that bovine segmental duplications corresponding to genes are significantly enriched for specific biological functions such as immunity, digestion, lactation and reproduction. Our results suggest that in most mammalian lineages segmental duplications are organized in a tandem configuration. Segmental duplications remain problematic for genome and assembly and we highlight genic regions that require higher quality sequence characterization. This study provides insights into mammalian genome evolution and generates a valuable resource for cattle genomics research.
Mapping main, epistatic and sex-specific QTL for body composition in a chicken population divergently selected for low or high growth rate
(Springer Nature, 2010-02-11) Ankra-Badu, Georgina A; Shriner, Daniel; Le Bihan-Duval, Elisabeth; Mignon-Grasteau, Sandrine; Pitel, Frédérique; Beaumont, Catherine; Duclos, Michel J; Simon, Jean; Porter, Tom E; Vignal, Alain; Cogburn, Larry A; Allison, David B; Yi, Nengjun; Aggrey, Samuel E
Delineating the genetic basis of body composition is important to agriculture and medicine. In addition, the incorporation of gene-gene interactions in the statistical model provides further insight into the genetic factors that underlie body composition traits. We used Bayesian model selection to comprehensively map main, epistatic and sex-specific QTL in an F2 reciprocal intercross between two chicken lines divergently selected for high or low growth rate. We identified 17 QTL with main effects across 13 chromosomes and several sex-specific and sex-antagonistic QTL for breast meat yield, thigh + drumstick yield and abdominal fatness. Different sets of QTL were found for both breast muscles [Pectoralis (P) major and P. minor], which suggests that they could be controlled by different regulatory mechanisms. Significant interactions of QTL by sex allowed detection of sex-specific and sex-antagonistic QTL for body composition and abdominal fat. We found several female-specific P. major QTL and sex-antagonistic P. minor and abdominal fatness QTL. Also, several QTL on different chromosomes interact with each other to affect body composition and abdominal fatness. The detection of main effects, epistasis and sex-dimorphic QTL suggest complex genetic regulation of somatic growth. An understanding of such regulatory mechanisms is key to mapping specific genes that underlie QTL controlling somatic growth in an avian model.
Transcriptional and pathway analysis in the hypothalamus of newly hatched chicks during fasting and delayed feeding
(2010-03-09) Higgins, Stacey E; Ellestad, Laura E; Trakooljul, Nares; McCarthy, Fiona; Saliba, Jason; Cogburn, Larry A; Porter, Tom E
Background: The hypothalamus plays a central role in regulating appetite and metabolism. However, the gene networks within the hypothalamus that regulate feed intake and metabolism, and the effects of fasting on those pathways are not completely understood in any species. The present experiment evaluated global hypothalamic gene expression in newly hatched chicks using microarray analysis to elucidate genes and pathways regulated by feeding, fasting, and delayed feeding. Ten groups of chicks were sampled over four days post-hatch, including fed, fasted, and 48 h fasted followed by access to feed for 4 h, 24 h, and 48 h. Hypothalamic samples were collected for microarray analysis (n = 4). Expression patterns of selected genes were confirmed by quantitative real-time PCR. Pathway analysis of the microarray results predicted a network of genes involved in neuropeptide or neurotransmitter signaling. To confirm the functionality of this predicted gene network, hypothalamic neurons from fed and fasted chicks were isolated and cultured in the presence of neuropeptide Y, somatostatin, amelanocyte stimulating hormone, norepinephrine, and L-phospho-serine. Results confirmed functional relationships among members of the predicted gene network. Moreover, the effects observed were dependant upon the nutritional state of the animals (fed vs. fasted). Results: Differences in gene expression (≥ 1.6 fold) were detected in 1,272 genes between treatments, and of those, 119 genes were significantly (P < 0.05) different. Pathway Miner analysis revealed that six genes (SSTR5, NPY5R, POMC, ADRB2, GRM8, and RLN3) were associated within a gene network. In vitro experiments with primary hypothalamic neurons confirmed that receptor agonists involved in this network regulated expression of other genes in the predicted network, and this regulation within the network was influenced by the nutritional status and age of the chick. Conclusions: Microarray analysis of the hypothalamus during different nutritional states revealed that many genes are differentially regulated. We found that functional interactions exist among six differentially regulated genes associated within a putative gene network from this experiment. Considering that POMC, an important gene in controlling metabolism, was central to this network, this gene network may play an important role in regulation of feeding and metabolism in birds.
Genomic characteristics of cattle copy number variations
(Springer Nature, 2011-02-23) Hou, Yali; Liu, George E; Bickhart, Derek M; Cardone, Maria Francesca; Wang, Kai; Kim, Eui-soo; Matukumalli, Lakshmi K; Ventura, Mario; Song, Jiuzhou; VanRaden, Paul M; Sonstegard, Tad S; Van Tassell, Curt P
Copy number variation (CNV) represents another important source of genetic variation complementary to single nucleotide polymorphism (SNP). High-density SNP array data have been routinely used to detect human CNVs, many of which have significant functional effects on gene expression and human diseases. In the dairy industry, a large quantity of SNP genotyping results are becoming available and can be used for CNV discovery to understand and accelerate genetic improvement for complex traits. We performed a systematic analysis of CNV using the Bovine HapMap SNP genotyping data, including 539 animals of 21 modern cattle breeds and 6 outgroups. After correcting genomic waves and considering the pedigree information, we identified 682 candidate CNV regions, which represent 139.8 megabases (~4.60%) of the genome. Selected CNVs were further experimentally validated and we found that copy number "gain" CNVs were predominantly clustered in tandem rather than existing as interspersed duplications. Many CNV regions (~56%) overlap with cattle genes (1,263), which are significantly enriched for immunity, lactation, reproduction and rumination. The overlap of this new dataset and other published CNV studies was less than 40%; however, our discovery of large, high frequency (> 5% of animals surveyed) CNV regions showed 90% agreement with other studies. These results highlight the differences and commonalities between technical platforms. We present a comprehensive genomic analysis of cattle CNVs derived from SNP data which will be a valuable genomic variation resource. Combined with SNP detection assays, gene-containing CNV regions may help identify genes undergoing artificial selection in domesticated animals.
A computational model of the hypothalamic - pituitary - gonadal axis in female fathead minnows (Pimephales promelas) exposed to 17α-ethynylestradiol and 17β-trenbolone
(Springer Nature, 2011-05-05) Li, Zhenhong; Kroll, Kevin J; Jensen, Kathleen M; Villeneuve, Daniel L; Ankley, Gerald T; Brian, Jayne V; Sepúlveda, María S; Orlando, Edward F; Lazorchak, James M; Kostich, Mitchell; Armstrong, Brandon; Denslow, Nancy D; Watanabe, Karen H
Endocrine disrupting chemicals (e.g., estrogens, androgens and their mimics) are known to affect reproduction in fish. 17α-ethynylestradiol is a synthetic estrogen used in birth control pills. 17β-trenbolone is a relatively stable metabolite of trenbolone acetate, a synthetic androgen used as a growth promoter in livestock. Both 17α-ethynylestradiol and 17β-trenbolone have been found in the aquatic environment and affect fish reproduction. In this study, we developed a physiologically-based computational model for female fathead minnows (FHM, Pimephales promelas), a small fish species used in ecotoxicology, to simulate how estrogens (i.e., 17α-ethynylestradiol) or androgens (i.e., 17β-trenbolone) affect reproductive endpoints such as plasma concentrations of steroid hormones (e.g., 17β-estradiol and testosterone) and vitellogenin (a precursor to egg yolk proteins). Using Markov Chain Monte Carlo simulations, the model was calibrated with data from unexposed, 17α-ethynylestradiol-exposed, and 17β-trenbolone-exposed FHMs. Four Markov chains were simulated, and the chains for each calibrated model parameter (26 in total) converged within 20,000 iterations. With the converged parameter values, we evaluated the model's predictive ability by simulating a variety of independent experimental data. The model predictions agreed with the experimental data well. The physiologically-based computational model represents the hypothalamic-pituitary-gonadal axis in adult female FHM robustly. The model is useful to estimate how estrogens (e.g., 17α-ethynylestradiol) or androgens (e.g., 17β-trenbolone) affect plasma concentrations of 17β-estradiol, testosterone and vitellogenin, which are important determinants of fecundity in fish.
Temporal transcriptome changes induced by MDV in marek's disease-resistant and -susceptible inbred chickens
(Springer Nature, 2011-10-12) Yu, Ying; Luo, Juan; Mitra, Apratim; Chang, Shuang; Tian, Fei; Zhang, Huanmin; Yuan, Ping; Zhou, Huaijun; Song, Jiuzhou
Marek's disease (MD) is a lymphoproliferative disease in chickens caused by Marek's disease virus (MDV) and characterized by T cell lymphoma and infiltration of lymphoid cells into various organs such as liver, spleen, peripheral nerves and muscle. Resistance to MD and disease risk have long been thought to be influenced both by genetic and environmental factors, the combination of which contributes to the observed outcome in an individual. We hypothesize that after MDV infection, genes related to MD-resistance or -susceptibility may exhibit different trends in transcriptional activity in chicken lines having a varying degree of resistance to MD. In order to study the mechanisms of resistance and susceptibility to MD, we performed genome-wide temporal expression analysis in spleen tissues from MD-resistant line 63, susceptible line 72 and recombinant congenic strain M (RCS-M) that has a phenotype intermediate between lines 63 and 72 after MDV infection. Three time points of the MDV life cycle in chicken were selected for study: 5 days post infection (dpi), 10dpi and 21dpi, representing the early cytolytic, latent and late cytolytic stages, respectively. We observed similar gene expression profiles at the three time points in line 63 and RCS-M chickens that are both different from line 72. Pathway analysis using Ingenuity Pathway Analysis (IPA) showed that MDV can broadly influence the chickens irrespective of whether they are resistant or susceptible to MD. However, some pathways like cardiac arrhythmia and cardiovascular disease were found to be affected only in line 72; while some networks related to cell-mediated immune response and antigen presentation were enriched only in line 63 and RCS-M. We identified 78 and 30 candidate genes associated with MD resistance, at 10 and 21dpi respectively, by considering genes having the same trend of expression change after MDV infection in lines 63 and RCS-M. On the other hand, by considering genes with the same trend of expression change after MDV infection in lines 72 and RCS-M, we identified 78 and 43 genes at 10 and 21dpi, respectively, which may be associated with MD-susceptibility. By testing temporal transcriptome changes using three representative chicken lines with different resistance to MD, we identified 108 candidate genes for MD-resistance and 121 candidate genes for MD-susceptibility over the three time points. Genes included in our resistance or susceptibility genes lists that are also involved in more than 5 biofunctions, such as CD8α, IL8, USP18, and CTLA4, are considered to be important genes involved in MD-resistance or -susceptibility. We were also able to identify several biofunctions related with immune response that we believe play an important role in MD-resistance.
miRNA-dysregulation associated with tenderness variation induced by acute stress in Angus cattle
(Springer Nature, 2012-06-01) Zhao, Chunping; Tian, Fei; Yu, Ying; Liu, George; Zan, Linsen; Updike, M Scott; Song, Jiuzhou
miRNAs are a class of small, single-stranded, non-coding RNAs that perform post-transcriptional repression of target genes by binding to 3’ untranslated regions. Research has found that miRNAs involved in the regulation of many metabolic processes. Here we uncovered that the beef quality of Angus cattle sharply diversified after acute stress. By performing miRNA microarray analysis, 13 miRNAs were significantly differentially expressed in stressed group compared to control group. Using a bioinformatics method, 135 protein-coding genes were predicted as the targets of significant differentially expressed miRNAs. Gene Ontology (GO) term and Ingenuity Pathway Analysis (IPA) mined that these target genes involved in some important pathways, which may have impact on meat quality and beef tenderness.
In vitro expansion of the mammary stem/progenitor cell population by xanthosine treatment
(Springer Nature, 2012-06-14) Choudhary, Ratan K; Capuco, Anthony V
Mammary stem cells are critical for growth and maintenance of the mammary gland and therefore are of considerable interest for improving productivity and efficiency of dairy animals. Xanthosine treatment has been demonstrated to promote expansion of putative mammary stem cells in vivo, and hepatic and hair follicle stem cells in vitro. In the latter, xanthosine promoted the symmetrical division of hepatic and hair follicle stem cells. The objective of this study was to determine if treating primary cultures of bovine mammary epithelial cells (MEC) with xanthosine increases the stem/progenitor cell population by promoting symmetrical division of mammary stem cells. In vitro treatment with xanthosine increased the population of MEC during the exponential phase of cell growth, reducing the doubling time from 86 h in control cultures to 60 h in xanthosine-treated cultures. The bromodeoxyuridine (BrdU) labeling index and the proportion of MEC in S-phase both were increased by xanthosine treatment, indicating that increased cell accretion was due to increased cell proliferation. Analysis of daughter-pairs indicated that xanthosine promoted a shift from asymmetric to symmetric cell division. Moreover, the 30 % increase in symmetric cell division was concomitant with an increase in the proportion of MEC that were positive for a putative stem cell marker (FNDC3B) and a trend toward increased telomerase activity. These results suggest that xanthosine treatment in vitro can increase cell proliferation, promote symmetric cell division and enhance stem/progenitor cell activity. Xanthosine treatment increased the proliferation rate of bovine MEC in vitro. This was likely to be mediated by an increase in the proportion of stem/progenitor cells in the MEC population due to promotion of symmetrical stem cell division by xanthosine.
Fine mapping of copy number variations on two cattle genome assemblies using high density SNP array
(Springer Nature, 2012-08-06) Hou, Yali; Bickhart, Derek M; Hvinden, Miranda L; Li, Congjun; Song, Jiuzhou; Boichard, Didier A; Fritz, Sébastien; Eggen, André; DeNise, Sue; Wiggans, George R; Sonstegard, Tad S; Van Tassell, Curtis P; Liu, George E
Btau_4.0 and UMD3.1 are two distinct cattle reference genome assemblies. In our previous study using the low density BovineSNP50 array, we reported a copy number variation (CNV) analysis on Btau_4.0 with 521 animals of 21 cattle breeds, yielding 682 CNV regions with a total length of 139.8 megabases. In this study using the high density BovineHD SNP array, we performed high resolution CNV analyses on both Btau_4.0 and UMD3.1 with 674 animals of 27 cattle breeds. We first compared CNV results derived from these two different SNP array platforms on Btau_4.0. With two thirds of the animals shared between studies, on Btau_4.0 we identified 3,346 candidate CNV regions representing 142.7 megabases (~4.70%) of the genome. With a similar total length but 5 times more event counts, the average CNVR length of current Btau_4.0 dataset is significantly shorter than the previous one (42.7 kb vs. 205 kb). Although subsets of these two results overlapped, 64% (91.6 megabases) of current dataset was not present in the previous study. We also performed similar analyses on UMD3.1 using these BovineHD SNP array results. Approximately 50% more and 20% longer CNVs were called on UMD3.1 as compared to those on Btau_4.0. However, a comparable result of CNVRs (3,438 regions with a total length 146.9 megabases) was obtained. We suspect that these results are due to the UMD3.1 assembly's efforts of placing unplaced contigs and removing unmerged alleles. Selected CNVs were further experimentally validated, achieving a 73% PCR validation rate, which is considerably higher than the previous validation rate. About 20-45% of CNV regions overlapped with cattle RefSeq genes and Ensembl genes. Panther and IPA analyses indicated that these genes provide a wide spectrum of biological processes involving immune system, lipid metabolism, cell, organism and system development. In this study using the high density BovineHD SNP array, we performed high resolution CNV analyses on both Btau_4.0 and UMD3.1 with 674 animals of 27 cattle breeds. We first compared CNV results derived from these two different SNP array platforms on Btau_4.0. With two thirds of the animals shared between studies, on Btau_4.0 we identified 3,346 candidate CNV regions representing 142.7 megabases (~4.70%) of the genome. With a similar total length but 5 times more event counts, the average CNVR length of current Btau_4.0 dataset is significantly shorter than the previous one (42.7 kb vs. 205 kb). Although subsets of these two results overlapped, 64% (91.6 megabases) of current dataset was not present in the previous study. We also performed similar analyses on UMD3.1 using these BovineHD SNP array results. Approximately 50% more and 20% longer CNVs were called on UMD3.1 as compared to those on Btau_4.0. However, a comparable result of CNVRs (3,438 regions with a total length 146.9 megabases) was obtained. We suspect that these results are due to the UMD3.1 assembly's efforts of placing unplaced contigs and removing unmerged alleles. Selected CNVs were further experimentally validated, achieving a 73% PCR validation rate, which is considerably higher than the previous validation rate. About 20-45% of CNV regions overlapped with cattle RefSeq genes and Ensembl genes. Panther and IPA analyses indicated that these genes provide a wide spectrum of biological processes involving immune system, lipid metabolism, cell, organism and system development. We present a comprehensive result of cattle CNVs at a higher resolution and sensitivity. We identified over 3,000 candidate CNV regions on both Btau_4.0 and UMD3.1, further compared current datasets with previous results, and examined the impacts of genome assemblies on CNV calling.
Application of orange essential oil as an antistaphylococcal agent in a dressing model
(Springer Nature, 2012-08-16) Muthaiyan, Arunachalam; Biswas, Debabrata; Crandall, Philip G; Wilkinson, Brian J; Ricke, Steven C
Staphylococcus aureus is the pathogen most often and prevalently involved in skin and soft tissue infections. In recent decades outbreaks of methicillin-resistant S. aureus (MRSA) have created major problems for skin therapy, and burn and wound care units. Topical antimicrobials are most important component of wound infection therapy. Alternative therapies are being sought for treatment of MRSA and one area of interest is the use of essential oils. With the increasing interest in the use and application of natural products, we screened the potential application of terpeneless cold pressed Valencia orange oil (CPV) for topical therapy against MRSA using an in vitro dressing model and skin keratinocyte cell culture model. The inhibitory effect of CPV was determined by disc diffusion vapor assay for MRSA and vancomycin intermediate-resistant S. aureus (VISA) strains. Antistaphylococcal effect of CPV in an in vitro dressing model was tested on S. aureus inoculated tryptic soya agar plate. Bactericidal effect of CPV on MRSA and VISA infected keratinocyte cells was examined by enumeration of extra- and intra-cellular bacterial cells at different treatment time points. Cytotoxic effects on human skin cells was tested by adding CPV to the keratinocyte (HEK001) cells grown in serum free KSFM media, and observed by phase-contrast microscope. CPV vapour effectively inhibited the MRSA and VISA strains in both disc diffusion vapour assay and in vitro dressing model. Compared to untreated control addition of 0.1% CPV to MRSA infected keratinocyte decreased the viable MRSA cells by 2 log CFU/mL in 1 h and in VISA strain 3 log CFU/mL reduction was observed in 1 h. After 3 h viable S. aureus cells were not detected in the 0.2% CPV treatment. Bactericidal concentration of CPV did not show any cytotoxic effect on the human skin keratinocyte cells in vitro. At lower concentration addition of CPV to keratinocytes infected with MRSA and VISA rapidly killed the bacterial cells without causing any toxic effect to the keratinocytes. Therefore, the results of this study warrant further in vivo study to evaluate the potential of CPV as a topical antistaphylococcal agent.
Marek’s disease virus infection induces widespread differential chromatin marks in inbred chicken lines
(Springer Nature, 2012-10-16) Mitra, Apratim; Luo, Juan; Zhang, Huanming; Cui, Kairong; Zhao, Keji; Song, Jiuzhou
Marek’s disease (MD) is a neoplastic disease in chickens caused by the MD virus (MDV). Successful vaccine development against MD has resulted in increased virulence of MDV and the understanding of genetic resistance to the disease is, therefore, crucial to long-term control strategies. Also, epigenetic factors are believed to be one of the major determinants of disease response. Here, we carried out comprehensive analyses of the epigenetic landscape induced by MDV, utilizing genome-wide histone H3 lysine 4 and lysine 27 trimethylation maps from chicken lines with varying resistance to MD. Differential chromatin marks were observed on genes previously implicated in the disease such as MX1 and CTLA-4 and also on genes reported in other cancers including IGF2BP1 and GAL. We detected bivalent domains on immune-related transcriptional regulators BCL6, CITED2 and EGR1, which underwent dynamic changes in both lines as a result of MDV infection. In addition, putative roles for GAL in the mechanism of MD progression were revealed. Our results confirm the presence of widespread epigenetic differences induced by MD in chicken lines with different levels of genetic resistance. A majority of observed epigenetic changes were indicative of increased levels of viral infection in the susceptible line symptomatic of lowered immunocompetence in these birds caused by early cytolytic infection. The GAL system that has known anti-proliferative effects in other cancers is also revealed to be potentially involved in MD progression. Our study provides further insight into the mechanisms of MD progression while revealing a complex landscape of epigenetic regulatory mechanisms that varies depending on host factors.
Transcriptional analysis of abdominal fat in genetically fat and lean chickens reveals adipokines, lipogenic genes and a link between hemostasis and leanness
(Springer Nature, 2013-08-16) Resnyk, Christopher W; Carré, Wilfrid; Wang, Xiaofei; Porter, Tom E; Simon, Jean; Le Bihan-Duval, Elisabeth; Duclos, Michael J; Aggrey, Sam E; Cogburn, Larry A
This descriptive study of the abdominal fat transcriptome takes advantage of two experimental lines of meat-type chickens (Gallus domesticus), which were selected over seven generations for a large difference in abdominal (visceral) fatness. At the age of selection (9 wk), the fat line (FL) and lean line (LL) chickens exhibit a 2.5-fold difference in abdominal fat weight, while their feed intake and body weight are similar. These unique avian models were originally created to unravel genetic and endocrine regulation of adiposity and lipogenesis in meat-type chickens. The Del-Mar 14K Chicken Integrated Systems microarray was used for a time-course analysis of gene expression in abdominal fat of FL and LL chickens during juvenile development (1–11 weeks of age). Microarray analysis of abdominal fat in FL and LL chickens revealed 131 differentially expressed (DE) genes (FDR≤0.05) as the main effect of genotype, 254 DE genes as an interaction of age and genotype and 3,195 DE genes (FDR≤0.01) as the main effect of age. The most notable discoveries in the abdominal fat transcriptome were higher expression of many genes involved in blood coagulation in the LL and up-regulation of numerous adipogenic and lipogenic genes in FL chickens. Many of these DE genes belong to pathways controlling the synthesis, metabolism and transport of lipids or endocrine signaling pathways activated by adipokines, retinoid and thyroid hormones. The present study provides a dynamic view of differential gene transcription in abdominal fat of chickens genetically selected for fatness (FL) or leanness (LL). Remarkably, the LL chickens over-express a large number of hemostatic genes that could be involved in proteolytic processing of adipokines and endocrine factors, which contribute to their higher lipolysis and export of stored lipids. Some of these changes are already present at 1 week of age before the divergence in fatness. In contrast, the FL chickens have enhanced expression of numerous lipogenic genes mainly after onset of divergence, presumably directed by multiple transcription factors. This transcriptional analysis shows that abdominal fat of the chicken serves a dual function as both an endocrine organ and an active metabolic tissue, which could play a more significant role in lipogenesis than previously thought.
Genome wide CNV analysis reveals additional variants associated with milk production traits in Holsteins
(Springer Nature, 2014-08-15) Xu, Lingyang; Cole, John B; Bickhart, Derek M; Hou, Yali; Song, Jiuzhou; VanRaden, Paul M; Sonstegard, Tad S; Van Tassell, Curtis P; Liu, George E
Milk production is an economically important sector of global agriculture. Much attention has been paid to the identification of quantitative trait loci (QTL) associated with milk, fat, and protein yield and the genetic and molecular mechanisms underlying them. Copy number variation (CNV) is an emerging class of variants which may be associated with complex traits. In this study, we performed a genome-wide association between CNVs and milk production traits in 26,362 Holstein bulls and cows. A total of 99 candidate CNVs were identified using Illumina BovineSNP50 array data, and association tests for each production trait were performed using a linear regression analysis with PCA correlation. A total of 34 CNVs on 22 chromosomes were significantly associated with at least one milk production trait after false discovery rate (FDR) correction. Some of those CNVs were located within or near known QTL for milk production traits. We further investigated the relationship between associated CNVs with neighboring SNPs. For all 82 combinations of traits and CNVs (less than 400 kb in length), we found 17 cases where CNVs directly overlapped with tag SNPs and 40 cases where CNVs were adjacent to tag SNPs. In 5 cases, CNVs located were in strong linkage disequilibrium with tag SNPs, either within or adjacent to the same haplotype block. There were an additional 20 cases where CNVs did not have a significant association with SNPs, suggesting that the effects of those CNVs were probably not captured by tag SNPs. We conclude that combining CNV with SNP analyses reveals more genetic variations underlying milk production traits than those revealed by SNPs alone.
Changes in various metabolic parameters in blood and milk during experimental Escherichia coli mastitis for primiparous Holstein dairy cows during early lactation
(Springer Nature, 2014-10-17) Moyes, Kasey M; Larsen, Torben; Sørensen, Peter; Ingvartsen, Klaus L
The objective of this study was to characterize the changes in various metabolic parameters in blood and milk during IMI challenge with Escherichia coli (E. coli) for dairy cows during early lactation. Thirty, healthy primiparous Holstein cows were infused (h = 0) with ~20-40 cfu of live E. coli into one front mammary quarter at ~4-6 wk in lactation. Daily feed intake and milk yield were recorded. At –12, 0, 3, 6, 12, 18, 24, 36, 48, 60, 72, 96, 108, 120, 132, 144, 156, 168, 180 and 192 h relative to challenge rectal temperatures were recorded and quarter foremilk was collected for analysis of shedding of E. coli. Composite milk samples were collected at -180, -132, -84, -36, -12, 12, 24, 36, 48, 60, 72, 84, 96, 132 and 180 h relative to challenge (h = 0) and analyzed for lactate dehydrogenase (LDH), somatic cell count, fat, protein, lactose, citrate, beta-hydroxybutyrate (BHBA), free glucose (fglu), and glucose-6-phosphate (G6P). Blood was collected at -12, 0, 3, 6, 12, 18, 24, 36, 60, 72, 84, 132 and 180 h relative to challenge and analyzed for plasma non-esterified fatty acids (NEFA), BHBA and glucose concentration. A generalized linear mixed model was used to determine the effect of IMI challenge on metabolic responses of cows during early lactation. By 12 h, E. coli was recovered from challenged quarters and shedding continued through 72 h. Rectal temperature peaked by 12 h post-challenge and returned to pre-challenge values by 36 h post-IMI challenge. Daily feed intake and milk yield decreased (P <0.05) by 1 and 2 d, respectively, after mastitis challenge. Plasma BHBA decreased (12 h; P <0.05) from 0.96 ± 1.1 at 0 h to 0.57 ± 0.64 mmol/L by 18 h whereas concentration of plasma NEFA (18 h) and glucose (24 h) were significantly greater, 11 and 27%, respectively, after challenge. In milk, fglu, lactose, citrate, fat and protein yield were lower whereas yield of BHBA and G6P were higher after challenge when compared to pre-challenge values. Changes in metabolites in blood and milk were most likely associated with drops in feed intake and milk yield. However, the early rise in plasma NEFA may also signify enhanced adipose tissue lipolysis. Lower concentrations of plasma BHBA may be attributed to an increase transfer into milk after IMI. Decreases in both milk lactose yield and % after challenge may be partly attributed to reduced conversion of fglu to lactose. Rises in G6P yield and concentration in milk after challenge (24 h) may signify increased conversion of fglu to G6P. Results identify changes in various metabolic parameters in blood and milk after IMI challenge with E. coli in dairy cows that may partly explain the partitioning of nutrients and changes in milk components after IMI for cows during early lactation.
X-inactivation informs variance-based testing for X-linked association of a quantitative trait
(Springer Nature, 2015) Ma, Li; Hoffman, Gabriel; Keinan, Alon
The X chromosome plays an important role in human diseases and traits. However, few X-linked associations have been reported in genome-wide association studies, partly due to analytical complications and low statistical power. In this study, we propose tests of X-linked association that capitalize on variance heterogeneity caused by various factors, predominantly the process of X-inactivation. In the presence of X-inactivation, the expression of one copy of the chromosome is randomly silenced. Due to the consequent elevated randomness of expressed variants, females that are heterozygotes for a quantitative trait locus might exhibit higher phenotypic variance for that trait. We propose three tests that build on this phenomenon: 1) A test for inflated variance in heterozygous females; 2) A weighted association test; and 3) A combined test. Test 1 captures the novel signal proposed herein by directly testing for higher phenotypic variance of heterozygous than homozygous females. As a test of variance it is generally less powerful than standard tests of association that consider means, which is supported by extensive simulations. Test 2 is similar to a standard association test in considering the phenotypic mean, but differs by accounting for (rather than testing) the variance heterogeneity. As expected in light of X-inactivation, this test is slightly more powerful than a standard association test. Finally, test 3 further improves power by combining the results of the first two tests. We applied the these tests to the ARIC cohort data and identified a novel X-linked association near gene AFF2 with blood pressure, which was not significant based on standard association testing of mean blood pressure. Variance-based tests examine overdispersion, thereby providing a complementary type of signal to a standard association test. Our results point to the potential to improve power of detecting X-linked associations in the presence of variance heterogeneity.