Combinatorial metabolomic and transcriptomic analysis of muscle growth in hybrid striped bass (female white bass Morone chrysops x male striped bass M. saxatilis).

BACKGROUND: Understanding growth regulatory pathways is important in aquaculture, fisheries, and vertebrate physiology generally. Machine learning pattern recognition and sensitivity analysis were employed to examine metabolomic small molecule profiles and transcriptomic gene expression data generated from liver and white skeletal muscle of hybrid striped bass (white bass Morone chrysops x striped bass M. saxatilis) representative of the top and bottom 10 % by body size of a production cohort. RESULTS: Larger fish (good-growth) had significantly greater weight, total length, hepatosomatic index, and specific growth rate compared to smaller fish (poor-growth) and also had significantly more muscle fibers of smaller diameter (≤ 20 µm diameter), indicating active hyperplasia. Differences in metabolomic pathways included enhanced energetics (glycolysis, citric acid cycle) and amino acid metabolism in good-growth fish, and enhanced stress, muscle inflammation (cortisol, eicosanoids) and dysfunctional liver cholesterol metabolism in poor-growth fish. The majority of gene transcripts identified as differentially expressed between groups were down-regulated in good-growth fish. Several molecules associated with important growth-regulatory pathways were up-regulated in muscle of fish that grew poorly: growth factors including agt and agtr2 (angiotensins), nicotinic acid (which stimulates growth hormone production), gadd45b, rgl1, zfp36, cebpb, and hmgb1; insulin-like growth factor signaling (igfbp1 and igf1); cytokine signaling (socs3, cxcr4); cell signaling (rgs13, rundc3a), and differentiation (rhou, mmp17, cd22, msi1); mitochondrial uncoupling proteins (ucp3, ucp2); and regulators of lipid metabolism (apoa1, ldlr). Growth factors pttg1, egfr, myc, notch1, and sirt1 were notably up-regulated in muscle of good-growing fish. CONCLUSION: A combinatorial pathway analysis using metabolomic and transcriptomic data collectively suggested promotion of cell signaling, proliferation, and differentiation in muscle of good-growth fish, whereas muscle inflammation and apoptosis was observed in poor-growth fish, along with elevated cortisol (an anti-inflammatory hormone), perhaps related to muscle wasting, hypertrophy, and inferior growth. These findings provide important biomarkers and mechanisms by which growth is regulated in fishes and other vertebrates as well.

Selection for high and low antibody responses to sheep red blood cells influences cytokine and chemokine expression in chicken peripheral blood leukocytes and splenic tissue.

White Leghorn chickens from a common founder population have been divergently selected for high (HAS) or low (LAS) antibody responses to sheep red blood cells (SRBC) for 49 generations resulting in 2 diverse lines for this trait. Much has been studied in these two lines; however, the impact of these selection pressures on cytokine and chemokine expression is not fully understood. The purpose of this study is to determine if selection for antibody response to SRBC impacts cytokine and chemokine expression in peripheral blood leukocytes (PBL) and spleen from HAS and LAS chickens. Total RNA was isolated from PBL and spleen after which mRNA expression of cytokines (IL4, IL6, IL10, TGF-ß4) and chemokines (CXCL8, CCL4) were determined by quantitative real-time RT-PCR (qRT-PCR). The data were analyzed using Student's t test comparing HAS and LAS (P < 0.05) and are reported as corrected 40-CT. PBL and spleen samples were analyzed separately. With respect to PBL, expression of IL6 was higher (P < 0.05) in PBL isolated from LAS chickens compared to those from the HAS line whereas there were no differences (P > 0.05) in IL4, IL10, CXCL8, CCL4, or TGF-ß4. The cytokine and chemokine mRNA expression profiles were different in the spleen between the two lines. IL4 and CXCL8 expression were higher (P < 0.05) in spleen samples from HAS chickens than LAS. The expression of IL6, IL10, CCL4, or TGF-ß4 in the spleens did not differ (P > 0.05) between the lines. The data indicate that selection for specific antibody responses to SRBC impacts the cytokine and chemokine expression profile in PBL and spleens but in different ways in HAS and LAS. These studies provide insight into the influence that selection pressures for antibody responses have on different immune response components, specifically cytokines and chemokines typically involved in the innate response.

Transgenerational epigenetic inheritance and immunity in chickens that vary in Marek's disease resistance.

Marek's disease virus (MDV), a naturally oncogenic, highly contagious alpha herpesvirus, induces a T cell lymphoma in chickens that causes severe economic loss. Marek's disease (MD) outcome in an individual is attributed to genetic and environmental factors. Further investigation of the host-virus interaction mechanisms that impact MD resistance is needed to achieve greater MD control. This study analyzed genome-wide DNA methylation patterns in 2 highly inbred parental lines 63 and 72 and 5 recombinant congenic strains (RCS) C, L, M, N, and X strains from those parents. Lines 63 and 72, are MD resistant and susceptible, respectively, whereas the RCS have different combinations of 87.5% Line 63 and 12.5% Line 72. Our DNA methylation cluster showed a strong association with MD incidence. Differentially methylated regions (DMRs) between the parental lines and the 5 RCS were captured. MD-resistant and MD-susceptible markers of DNA methylation were identified as transgenerational epigenetic inheritable. In addition, the growth of v-src DNA tumors and antibody response against sheep red blood cells differed among the 2 parental lines and the RCS. Overall, our results provide very solid evidence that DNA methylation patterns are transgenerational epigenetic inheritance (TEI) in chickens and also play a vital role in MD tumorigenesis and other immune responses; the specific methylated regions may be important modulators of general immunity.

Combining supervised machine learning with statistics reveals differential gene expression patterns related to energy metabolism in the jejuna of chickens divergently selected for antibody response to sheep red blood cells.

Since the 1970s, 2 lines of White Leghorn chickens, HAS and LAS, have been continuously divergently selected for 5-day postinjection antibody titer to injection with sheep red blood cells (SRBC). Antibody response is a complex genetic trait and characterizing differences in gene expression could facilitate better understanding of physiological changes due to selection and antigen exposure. At 41 d of age, randomly selected HAS and LAS chickens, which had been coraised from hatch, were either injected with SRBC (HASI and LASI) or kept as the noninjected cohort (HASN and LASN). Five days later, all were euthanized, and samples collected from the jejunum for RNA isolation and sequencing. Resulting gene expression data were analyzed combining traditional statistics with machine learning to obtain signature gene lists for functional analysis. Differences in ATP production and cellular processes were observed in the jejunum between lines and following SRBC injection. HASN vs. LASN exhibited upregulation of ATP production, immune cell motility, and inflammation. LASI exhibits upregulation of ATP production and protein synthesis vs. LASN, reflective of what was observed in HASN vs. LASN. In contrast, no corresponding upregulation of ATP production was observed in HASI vs. HASN, and most other cellular processes appear inhibited. Without exposure to SRBC, gene expression in the jejunum indicates HAS generates more ATP than LAS, suggesting HAS maintains a "primed" system; and gene expression of HASI vs. HASN further suggests this basal ATP production is sufficient for robust antibody responses. Conversely, LASI vs. LASN jejunal gene expression implies a physiological need for increased ATP production with only minimal correlating antibody production. The results of this experiment provide insight into energetic resource needs and allocations in the jejunum in response to genetic selection and antigen exposure in HAS and LAS which may help explain phenotypic differences observed in antibody response.

Differences in the microbiome of the small intestine of Leghorn lines divergently selected for antibody titer to sheep erythrocytes suggest roles for commensals in host humoral response.

For forty generations, two lines of White Leghorn chickens have been selected for high (HAS) or low (LAS) antibody response to a low dose injection of sheep red blood cells (SRBCs). Their gut is home to billons of microorganisms and the largest number of immune cells in the body; therefore, the objective of this experiment was to gain understanding of the ways the microbiome may influence the differential antibody response observed in these lines. We achieved this by characterizing the small intestinal microbiome of HAS and LAS chickens, determining their functional microbiome profiles, and by using machine learning to identify microbes which best differentiate HAS from LAS and associating the abundance of those microbes with host gene expression. Microbiome sequencing revealed greater diversity in LAS but statistically higher abundance of several strains, particularly those of Lactobacillus, in HAS. Enrichment of microbial metabolites implicated in immune response such as lactic acid, short chain fatty acids, amino acids, and vitamins were different between HAS and LAS. The abundance of several microbial strains corresponds to enriched host gene expression pathways related to immune response. These data provide a compelling argument that the microbiome is both likely affected by host divergent genetic selection and that it exerts influence on host antibody response by various mechanisms.

Immunomodulatory effects of heat stress and lipopolysaccharide on the bursal transcriptome in two distinct chicken lines.

BACKGROUND: Exposure to heat stress suppresses poultry immune responses, which can increase susceptibility to infectious diseases and, thereby, intensify the negative effects of heat on poultry welfare and performance. Identifying genes and pathways that are affected by high temperatures, especially heat-induced changes in immune responses, could provide targets to improve disease resistance in chickens. This study utilized RNA-sequencing (RNA-seq) to investigate transcriptome responses in the bursa of Fabricius, a primary immune tissue, after exposure to acute heat stress and/or subcutaneous immune stimulation with lipopolysaccharide (LPS) in a 2 × 2 factorial design: Thermoneutral + Saline, Heat + Saline, Thermoneutral + LPS and Heat + LPS. All treatments were investigated in two chicken lines: a relatively heat- and disease-resistant Fayoumi line and a more susceptible broiler line. RESULTS: Differential expression analysis determined that Heat + Saline had limited impact on gene expression (N = 1 or 63 genes) in broiler or Fayoumi bursa. However, Thermoneutral + LPS and Heat + LPS generated many expression changes in Fayoumi bursa (N = 368 and 804 genes). Thermoneutral + LPS was predicted to increase immune-related cell signaling and cell migration, while Heat + LPS would activate mortality-related functions and decrease expression in WNT signaling pathways. Further inter-treatment comparisons in the Fayoumi line revealed that heat stress prevented many of the expression changes caused by LPS. Although fewer significant expression changes were observed in the broiler bursa after exposure to Thermoneutral + LPS (N = 59 genes) or to Heat + LPS (N = 146 genes), both treatments were predicted to increase cell migration. Direct comparison between lines (broiler to Fayoumi) confirmed that each line had distinct responses to treatment. CONCLUSIONS: Transcriptome analysis identified genes and pathways involved in bursal responses to heat stress and LPS and elucidated that these effects were greatest in the combined treatment. The interaction between heat and LPS was line dependent, with suppressive expression changes primarily in the Fayoumi line. Potential target genes, especially those involved in cell migration and immune signaling, can inform future research on heat stress in poultry and could prove useful for improving disease resistance.

Genome-wide standing variation facilitates long-term response to bidirectional selection for antibody response in chickens.

BACKGROUND: Long-term selection experiments provide a powerful approach to gain empirical insights into adaptation, allowing researchers to uncover the targets of selection and infer their contributions to the mode and tempo of adaptation. Here we implement a pooled genome re-sequencing approach to investigate the consequences of 39 generations of bidirectional selection in White Leghorn chickens on a humoral immune trait: antibody response to sheep red blood cells. RESULTS: We observed wide genome involvement in response to this selection regime. Many genomic regions were highly differentiated resulting from this experimental selection regime, an involvement of up to 20% of the chicken genome (208.8 Mb). While genetic drift has certainly contributed to this, we implement gene ontology, association analysis and population simulations to increase our confidence in candidate selective sweeps. Three strong candidate genes, MHC, SEMA5A and TGFBR2, are also presented. CONCLUSIONS: The extensive genomic changes highlight the polygenic genetic architecture of antibody response in these chicken populations, which are derived from a common founder population, demonstrating the extent of standing immunogenetic variation available at the onset of selection.

Quantitative trait loci identified for blood chemistry components of an advanced intercross line of chickens under heat stress.

BACKGROUND: Heat stress in poultry results in considerable economic losses and is a concern for both animal health and welfare. Physiological changes occur during periods of heat stress, including changes in blood chemistry components. A highly advanced intercross line, created from a broiler (heat susceptible) by Fayoumi (heat resistant) cross, was exposed to daily heat cycles for seven days starting at 22 days of age. Blood components measured pre-heat treatment and on the seventh day of heat treatment included pH, pCO2, pO2, base excess, HCO3, TCO2, K, Na, ionized Ca, hematocrit, hemoglobin, sO2, and glucose. A genome-wide association study (GWAS) for these traits and their calculated changes was conducted to identify quantitative trait loci (QTL) using a 600 K SNP panel. RESULTS: There were significant increases in pH, base excess, HCO3, TCO2, ionized Ca, hematocrit, hemoglobin, and sO2, and significant decreases in pCO2 and glucose after 7 days of heat treatment. Heritabilities ranged from 0.01-0.21 for pre-heat measurements, 0.01-0.23 for measurements taken during heat, and 0.00-0.10 for the calculated change due to heat treatment. All blood components were highly correlated within measurement days, but not correlated between measurement days. The GWAS revealed 61 QTL for all traits, located on GGA (Gallus gallus chromosome) 1, 3, 6, 9, 10, 12-14, 17, 18, 21-28, and Z. A functional analysis of the genes in these QTL regions identified the Angiopoietin pathway as significant. The QTL that co-localized for three or more traits were on GGA10, 22, 26, 28, and Z and revealed candidate genes for birds' response to heat stress. CONCLUSIONS: The results of this study contribute to our knowledge of levels and heritabilities of several blood components of chickens under thermoneutral and heat stress conditions. Most components responded to heat treatment. Mapped QTL may serve as markers for genomic selection to enhance heat tolerance in poultry. The Angiopoietin pathway is likely involved in the response to heat stress in chickens. Several candidate genes were identified, giving additional insight into potential mechanisms of physiologic response to high ambient temperatures.

Differences in mRNA populations of short and long bovine conceptuses on Day 15 of gestation.

The majority of pregnancy loss in cattle occurs between Days 8 and 16 of gestation, coincident with the initiation of conceptus elongation and the onset of maternal recognition of pregnancy. Differences in conceptus length on the same day of gestation may be related to an inherent lack of developmental competency or may simply be a consequence of asynchrony with the maternal environment. The objective of this work was to characterize differential patterns of mRNA expression between short and long bovine conceptuses recovered on Day 15 of gestation. Embryos were produced from super-ovulated Holstein donor cows, and groups of Grade-1 and Grade-3 compact morulas were transferred into recipient heifers at Day 6.5 of their cycle. Conceptuses were recovered at Day 15 of gestation, and measured to assess overall length and area. Total RNA was extracted and analyzed on individual GeneChip Bovine Genome Arrays (Affymetrix, Santa Clara, CA). Gene expression was compared between conceptuses derived from the transfer of Grade-1 versus Grade-3 embryos; no differences were identified in the profiles of Day-15 conceptuses of these different embryo grades. When gene expression was compared between conceptuses classified as either short (mean length of 4.2 ± 0.1 mm [standard error]) or long (24.7 ± 1.9 mm) upon recovery at Day 15 of gestation, a total of 348 genes were differentially expressed. Of these, 221 genes were up-regulated and 127 were down-regulated in long compared to short conceptuses. In summary, differences in gene expression were identified between conceptuses recovered on Day 15 of gestation, based on their length. These data may be used to identify genes and cellular pathways involved in enhanced conceptus elongation that could serve as markers of successful pregnancy. Mol. Reprod. Dev. 83: 424-441, 2016. © 2016 Wiley Periodicals, Inc.

Developmental profile of claudin-3, -5, and -16 proteins in the epithelium of chick intestine.

Proteins in the claudin family are a main component of tight junctions and form a seal that modulates paracellular transport in intestinal epithelium. This research tests the hypothesis that claudins 3, 5, and 16 will appear in the epithelium of embryonic intestine during functional differentiation. Immunohistochemistry is utilized to explore the developmental patterns of claudin-3, -5, and -16 proteins in the epithelium of embryonic chick intestine from 9 days prior to hatching through the early post-hatch period. These claudin proteins either changed their cellular localization or first appeared around the time of hatching. After hatching, claudin-3 expression was prominent in basal-lateral regions of the epithelium along the entire villus, but was absent from crypts. Claudin-5 was expressed most strongly in the crypt and lower villus epithelium within junctional complexes, whereas immunostaining of claudin-16 was localized within goblet cells of the upper villus region. The relative mRNA levels of claudin-3, -5, and -16 showed similar patterns; transcript levels rose between 18 and 20 days of development, then dropped by 2 days post-hatch. Results of this work indicate that the claudin proteins assume their final locations within the epithelium around the time of hatching, suggesting that in addition to their known barrier and fence functions within tight junctions, these claudins may have additional roles in the differentiation and/or physiological function of chick intestine. The localization of claudin-16 to goblet cells and its distribution in the more mature cells of the upper villus region suggest an unexpected role in goblet cell maturation and mucus secretion.

Molecular genotype identification of the Gallus gallus major histocompatibility complex.

The chicken major histocompatibility complex (MHC) is commonly defined by serologic reactions of erythrocytes with antibodies specific to the highly polymorphic MHC class I (BF) and MHC class IV (BG) antigens. The microsatellite marker LEI0258 is known to be physically located within the MHC, between the BG and BF regions. DNA from various serologically defined MHC haplotypes was amplified by polymerase chain reaction with primers surrounding this marker. Twenty-six distinctive allele sizes were identified. Some serologically well-defined MHC haplotypes shared a common LEI0258 allele size but could be distinguished either by the addition of information from another nearby marker (MCW0371) or by small indels or single nucleotide polymorphism (SNP) differences between the alleles. The association between LEI0258 allele and serologically defined MHC haplotype was very consistent for the same haplotype from multiple sources. Sequence information for the region defined by LEI0258 was obtained for 51 different haplotypes. Two internal repeats whose lengths were 13 and 12 bp, respectively, are the primary basis for allelic variability. Allele size variation ranges from 182 to 552 bp. Four indels and five SNPs in the surrounding sequence provide additional means for distinguishing alleles. Typing with LEI0258 and MCW0371 will be useful in identifying MHC haplotypes in outbred populations of chickens particularly for the initial development of serological reagents.

Expressed sequence tag analysis of Eimeria-stimulated intestinal intraepithelial lymphocytes in chickens.

Intraepithelial lymphocytes (IELs) play a critical role in protective immune response to intestinal pathogens such as Eimeria, the etiologic agent of avian coccidiosis. A list of genes expressed by intestinal IELs of Eimeria-infected chickens was compiled using the expressed sequence tag (EST) strategy. The 14,409 ESTs consisted of 1851 clusters and 7595 singletons, which revealed 9446 unique genes in the data set. Comparison of the sequence data with chicken DNA sequences in GenBank identified 125 novel clones. This EST library will provide a valuable resource for profiling global gene expression in normal and pathogen-infected chickens and identifying additional unique immune-related genes.

Plasma clearance and tissue distribution of radiolabeled leptin in the chicken.

Leptin is an adipose and liver tissue-derived secreted protein in chickens that has been implicated in the regulation of food intake and whole-body energy balance. In this study, the metabolic clearance and tissue uptake of leptin were examined in the chicken (Gallus gallus). Four-week-old broiler males were infused with (125)I-labeled mouse leptin. Chromatography of radiolabeled leptin in plasma produced two peaks, one at 16 kDa (free leptin) and a free iodine peak. No leptin binding protein in blood was detected. Leptin was cleared with a half-life estimate of 23 min. In order to investigate the tissue distribution and uptake of radiolabeled leptin, multiple tissues were removed from infused birds at 15 and 240 min post-infusion, and trichloroacetic acid (TCA)-precipitable radioactivity was determined. The amounts of radioactivity at 15 min post-infusion in the tissues in rank order were: kidney, testis, lung, spleen, heart, liver, small and large intestine, gizzard, pancreas, bursa, leg and breast muscle, adrenals, and brain. A slightly different pattern of distribution was observed at 240 min post-infusion. We conclude from these studies that unlike mammals, no circulating leptin binding protein is present in chickens. Leptin is metabolized and cleared very rapidly from blood by the kidney.

Feed restriction significantly alters lipogenic gene expression in broiler breeder chickens.

Broiler breeder pullets were divided into two groups at 21 wk of age. One group was given free access to feed (ad libitum) and the other fed a limited amount of feed (restricted). At 22 wk, all birds were photostimulated and maintained throughout an egg-laying cycle ending at 36 wk. Samples of liver and abdominal fat pad were collected just before photostimulation (prelight), after photostimulation at first egg and at peak egg production (plateau). Hepatic expression of sterol regulatory element binding protein-1, ATP-citrate lyase, fatty acid synthase, malic enzyme, acetyl-CoA carboxylase and stearoyl-CoA (Delta9) desaturase 1 genes in ad libitum birds declined from their highest levels just before photostimulation as the birds came into and maintained egg production. In contrast, the restricted birds had significant (P < 0.05) increases in the expression of these genes after photostimulation at first egg with a subsequent decline as they reached peak egg production. Hepatic expression of fatty acid binding protein, VLDL apolipoprotein (apoVLDL-II) and apoB genes increased significantly (P < 0.05) in both ad libitum and restricted breeders after photostimulation, whereas apoA1 gene expression declined during this time. Abdominal fat pad weights were significantly (P < 0.05) higher in the ad libitum compared with restricted birds after photostimulation. Lipoprotein lipase in this tissue showed a pattern of expression similar to that observed for the hepatic lipogenic enzyme genes. In conclusion, feed restriction during the pullet-to-breeder transition period significantly (P < 0.05) altered hepatic lipogenic gene expression in broiler breeders.

Expression of an uncoupling protein gene homolog in chickens.

An avian uncoupling protein (UCP) gene homolog was recently sequenced from skeletal muscle and was proposed to have a role in thermogenesis in chickens, ducks and hummingbirds. Since mammalian UCP 2 and UCP 3 also appear to have functions associated with energy and substrate partitioning and body weight regulation, the purpose of this study was to further characterize chicken UCP under conditions of nutritional stress and/or leptin administration. Male 3-week-old chickens were starved for 24 or 48 h and then half of each group was refed for an additional 24 h. In a follow-up experiment, chickens were fed or starved for 48 h with or without leptin administration. Feed deprivation increased UCP mRNA expression in skeletal muscle by up to 260% (P<0.001), and in a time-dependent manner in pectoralis muscle. Refeeding for 24 h normalized muscle UCP mRNA levels. Leptin administration had no effect on muscle UCP. Chicken muscle UCP mRNA levels were highly correlated with plasma triglyceride and non-esterified fatty acid (NEFA) concentrations, and with circulating levels of insulin, insulin-like growth factor (IGF)-I and IGF-II. These results suggest that, as in mammals, avian UCP is up-regulated during feed deprivation and is highly correlated with increased fatty acid oxidation and flux into skeletal muscle.