Upstream Regulator Analysis of Wooden Breast Myopathy Proteomics in Commercial Broilers and Comparison to Feed Efficiency Proteomics in Pedigree Male Broilers.

In an effort to understand the apparent trade-off between the continual push for growth performance and the recent emergence of muscle pathologies, shotgun proteomics was conducted on breast muscle obtained at ~8 weeks from commercial broilers with wooden breast (WB) myopathy and compared with that in pedigree male (PedM) broilers exhibiting high feed efficiency (FE). Comparison of the two proteomic datasets was facilitated using the overlay function of Ingenuity Pathway Analysis (IPA) (Qiagen, CA, USA). We focused on upstream regulator analysis and disease-function analysis that provides predictions of activation or inhibition of molecules based on (a) expression of downstream target molecules, (b) the IPA scientific citation database. Angiopoeitin 2 (ANGPT2) exhibited the highest predicted activation Z-score of all molecules in the WB dataset, suggesting that the proteomic landscape of WB myopathy would promote vascularization. Overlaying the FE proteomics data on the WB ANGPT2 upstream regulator network presented no commonality of protein expression and no prediction of ANGPT2 activation. Peroxisome proliferator coactivator 1 alpha (PGC1α) was predicted to be inhibited, suggesting that mitochondrial biogenesis was suppressed in WB. PGC1α was predicted to be activated in high FE pedigree male broilers. Whereas RICTOR (rapamycin independent companion of mammalian target of rapamycin) was predicted to be inhibited in both WB and FE datasets, the predictions were based on different downstream molecules. Other transcription factors predicted to be activated in WB muscle included epidermal growth factor (EGFR), X box binding protein (XBP1), transforming growth factor beta 1 (TGFB1) and nuclear factor (erythroid-derived 2)-like 2 (NFE2L2). Inhibitions of aryl hydrocarbon receptor (AHR), AHR nuclear translocator (ARNT) and estrogen related receptor gamma (ESRRG) were also predicted in the WB muscle. These findings indicate that there are considerable differences in upstream regulators based on downstream protein expression observed in WB myopathy and in high FE PedM broilers that may provide additional insight into the etiology of WB myopathy.

The role of microRNA-33 as a key regulator in hepatic lipogenesis signaling and a potential serological biomarker for NAFLD with excessive dietary fructose consumption in C57BL/6N mice.

Limited studies reported mechanisms by which microRNAs (miRNA) are interlinked in the etiology of fructose-induced non-alcoholic fatty liver disease (NAFLD). Here, we aimed to investigate the significance of miRNAs in fructose-induced NAFLD pathogenesis through unbiased approaches. In experiment I, C57BL/6N mice were fed either water or 34% fructose for six weeks ad libitum. In experiment II, time course effects of fructose intervention were monitored using the same conditions; mice were killed at the baseline, fourth, and sixth weeks. Bioinformatic analyses for hepatic proteomics revealed that SREBP1 is the most significant upstream regulator influenced by fructose; miR-33-5p (miR-33) was identified as the key miRNA responsible for SREBP1 regulation upon fructose intake, which was validated by in vitro transfection assay. In experiment II, we confirmed that the longer mice consumed fructose, the more severe liver injury markers (e.g., serum AST) appeared. Moreover, hepatic Srebp1 mRNA expression was increased depending upon the duration of fructose consumption. Hepatic miR-33 was time-dependently decreased by fructose while serum miR-33 expression was increased; these observations indicated that miR-33 from the liver might be released upon cell damage. Finally we observed that fructose-induced ferroptosis might be a cause of liver toxicity, resulting from oxidative damage. Collectively, our findings suggest that fructose-induced oxidative damage induces ferroptosis, and miR-33 could be used as a serological biomarker of fructose-induced NAFLD.

Characterization of stress response involved in chicken myopathy.

Myopathies (Woody Breast (WB) and White Striping (WS)) of broiler chickens have been correlated with fast growth. Recent studies reported that localized hypoxia and metabolic impairment may involve in these myopathies of birds. In order to better understand the stress response mechanisms affecting myopathies of broilers, the aim of this study was to examine effects of WB and both WB/WS on stress hormone corticosterone (CORT) levels and expressional changes of stress response genes including glucocorticoid (GC) receptor (GR), 11ß-Hydroxysteroid dehydrogenase type 1 (11ß-HSD1), DNA methylation regulators (DNMTs), and arginine vasotocin receptor 1a and 1b (V1aR, V1bR). Results of radioimmunoassay showed that CORT levels of WB and WB/WS birds were significantly higher compared to Con (p < 0.05), however, the combination of WB/WS was not significantly higher than WB birds, implying that the effects of WB and WS on CORT are not synergistic. Hepatic GR expression of both WB and WB/WS birds were significantly higher compared to Con (p < 0.05). However, GR expression levels in breast muscle of both WB and WB/WS birds were decreased compared to Con (p < 0.05). Hepatic 11ß-HSD1 expression was increased only in WB/WS birds compared to Con birds with no significant difference between Con and WB birds. 11ß-HSD1 expression was decreased and increased in WB and WB/WS birds compared to Con, respectively, in breast muscle (p < 0.05). DNMT1 expression was significantly decreased in both muscle and liver of WB birds, and in muscle of WB/WS birds, but not in liver of WB/WS birds, indicating differential effects of WS on the epigenetical stress response of muscle and liver compared to WB. V1aR expression was significantly increased in muscle of WB birds, and in liver of WB/WS birds compared to Con birds (p < 0.05). V1bR was not changed in muscle and liver of WB birds compared to Con birds. Taken together, results suggest that GC-induced myopathies occur in fast-growing broiler chickens and circulating CORT level might be a significant biochemical marker of myopathies (WB and WS) of birds. In addition, chronic stress responses in breast muscle and tissue-specific epigenetic changes of stress response genes by DNMTs may play a critical role in the occurrence of myopathies.

Acute alcohol consumption-induced let-7a inhibition exacerbates hepatic apoptosis by regulating Rb1 in mice.

Alcohol consumption is a critical risk factor for hepatic pathogenesis, including alcoholic liver diseases (ALD), but implications of alcohol-induced dysregulation of microRNA (miRNA) in ALD pathogenesis are not completely understood. In the present study, C57BL/6J male mice were treated with saline (CON; oral gavage; n = 8) or alcohol (EtOH; 3 g/kg body weight; oral gavage; n = 8) for 7 days. A total of 599 miRNAs and 158 key mRNAs related to fatty liver and hepatotoxicity pathways were assessed in mice liver tissues. The mRNA expression datasets were then utilized to predict interactions with miRNAs that were changed by alcohol consumption. Predicted miRNA-mRNA interactions were validated using in vitro miRNA transfection experiments. The results showed that let-7a was significantly decreased in the EtOH group and Rb1 mRNA was predicted as a target gene. This was further supported by an inverse correlation of RB1 and let-7a expression in mice liver tissue. Additionally, key protein expressions involved in RB1-apoptosis axis [i.e., p73, cleaved CASP-3 (cCASP-3), and cCASP-7] showed a trend of increase in the EtOH mice; this was also confirmed by capase-3 enzyme activity and terminal deoxynucleotidyl transferase-mediated dUTP nick-end labeling assay in livers of mice that had consumed alcohol. In line with our in vivo observations, alcohol treatment suppressed the let-7a expression and subsequently upregulated p73, cCASP-3, and cCASP-7 protein expressions in mice hepatocytes. Additional proteins in the apoptosis regulatory pathway (i.e., MDM2-p53 axis) were significantly changed in response to let-7a suppression in the cells. Taken together, the current study provides mechanistic evidence that alcohol consumption-induced let-7a suppression results in the upregulation of RB1, thereby promoting hepatic apoptosis through induction of pro-apoptotic proteins (e.g., p73), and by, at least in part, preventing MDM2-mediated p53 degradation.

Hepatic Transcriptomics Reveals that Lipogenesis Is a Key Signaling Pathway in Isocitrate Dehydrogenase 2 Deficient Mice.

Mitochondrial nicotinamide adenine dinucleotide phosphate (NADP+)-dependent isocitrate dehydrogenase (IDH2) plays a key role in the intermediary metabolism and energy production via catalysing oxidative decarboxylation of isocitrate to α-ketoglutarate in the tricarboxylic acid (TCA) cycle. Despite studies reporting potential interlinks between IDH2 and various diseases, there is lack of effort to comprehensively characterize signature(s) of IDH2 knockout (IDH2 KO) mice. A total of 6583 transcripts were identified from both wild-type (WT) and IDH2 KO mice liver tissues. Afterwards, 167 differentially expressed genes in the IDH2 KO group were short-listed compared to the WT group based on our criteria. The online bioinformatic analyses indicated that lipid metabolism is the most significantly influenced metabolic process in IDH2 KO mice. Moreover, the TR/RXR activation pathway was predicted as the top canonical pathway significantly affected by IDH2 KO. The key transcripts found in the bioinformatic analyses were validated by qPCR analysis, corresponding to the transcriptomics results. Further, an additional qPCR analysis confirmed that IDH2 KO caused a decrease in hepatic de novo lipogenesis via the activation of the fatty acid ß-oxidation process. Our unbiased transcriptomics approach and validation experiments suggested that IDH2 might play a key role in homeostasis of lipid metabolism.

Hepatic transcriptomics analysis reveals that fructose intervention down-regulated xenobiotics-metabolising enzymes through aryl hydrocarbon receptor signalling suppression in C57BL/6N mice.

For decades, fructose intake has been recognised as an environmental risk for metabolic syndromes and diseases. Here we comprehensively examined the effects of fructose intake on mice liver transcriptomes. Fructose-supplemented water (34 %; w/v) was fed to both male and female C57BL/6N mice at their free will for 6 weeks, followed by hepatic transcriptomics analysis. Based on our criteria, differentially expressed genes (DEG) were selected and subjected to further computational analyses to predict key pathways and upstream regulator(s). Subsequently, predicted genes and pathways from the transcriptomics dataset were validated via quantitative RT-PCR analyses. As a result, we identified eighty-nine down-regulated and eighty-eight up-regulated mRNA in fructose-fed mice livers. These DEG were subjected to bioinformatics analysis tools in which DEG were mainly enriched in xenobiotic metabolic processes; further, in the Ingenuity Pathway Analysis software, it was suggested that the aryl hydrocarbon receptor (AhR) is an upstream regulator governing overall changes, while fructose suppresses the AhR signalling pathway. In our quantitative RT-PCR validation, we confirmed that fructose suppressed AhR signalling through modulating expressions of transcription factor (AhR nuclear translocator; Arnt) and upstream regulators (Ncor2, and Rb1). Altogether, we demonstrated that ad libitum fructose intake suppresses the canonical AhR signalling pathway in C57BL/6N mice liver. Based on our current observations, further studies are warranted, especially with regard to the effects of co-exposure to fructose on (1) other types of carcinogens and (2) inflammation-inducing agents (or even diets such as a high-fat diet), to find implications of fructose-induced AhR suppression.

Copy number variation study in Japanese quail associated with stress related traits using whole genome re-sequencing data.

Copy number variation (CNV) is a major driving factor for genetic variation and phenotypic diversity in animals. To detect CNVs and understand genetic components underlying stress related traits, we performed whole genome re-sequencing of pooled DNA samples of 20 birds each from High Stress (HS) and Low Stress (LS) Japanese quail lines using Illumina HiSeq 2×150 bp paired end method. Sequencing data were aligned to the quail genome and CNVnator was used to detect CNVs in the aligned data sets. The depth of coverage for the data reached to 41.4x and 42.6x for HS and LS birds, respectively. We identified 262 and 168 CNV regions affecting 1.6 and 1.9% of the reference genome that completely overlapped 454 and 493 unique genes in HS and LS birds, respectively. Ingenuity pathway analysis showed that the CNV genes were significantly enriched to phospholipase C signaling, neuregulin signaling, reelin signaling in neurons, endocrine and nervous development, humoral immune response, and carbohydrate and amino acid metabolisms in HS birds, whereas CNV genes in LS birds were enriched in cell-mediated immune response, and protein and lipid metabolisms. These findings suggest CNV genes identified in HS and LS birds could be candidate markers responsible for stress responses in birds.

Transcriptomic analysis of the development of skeletal muscle atrophy in cancer-cachexia in tumor-bearing mice.

Cancer-cachexia (CC) is a wasting condition directly responsible for 20-40% of cancer-related deaths. The mechanisms controlling development of CC-induced muscle wasting are not fully elucidated. Most investigations focus on the postcachectic state and do not examine progression of the condition. We recently demonstrated mitochondrial degenerations precede muscle wasting in time course progression of CC. However, the extent of muscle perturbations before wasting in CC is unknown. Therefore, we performed global gene expression analysis in CC-induced muscle wasting to enhance understanding of intramuscular perturbations across the development of CC. Lewis lung carcinoma (LLC) was injected into the hind-flank of C57BL6/J mice at 8 wk of age with tumor allowed to develop for 1, 2, 3, or 4 wk and compared with PBS-injected control. Muscle wasting was evident at 4 wk LLC. RNA sequencing of gastrocnemius muscle samples showed widespread alterations in LLC compared with PBS animals with largest differences seen in 4 wk LLC, suggesting extensive transcriptomic alterations concurrent to muscle wasting. Commonly altered pathways included: mitochondrial dysfunction and protein ubiquitination, along with other less studied processes in this condition regulating transcription/translation and cytoskeletal structure. Current findings present novel evidence of transcriptomic shifts and altered cellular pathways in CC-induced muscle wasting.

Whole Genome Resequencing of Arkansas Progressor and Regressor Line Chickens to Identify SNPs Associated with Tumor Regression.

Arkansas Regressor (AR) chickens, unlike Arkansas Progressor (AP) chickens, regress tumors induced by the v-src oncogene. To better understand the genetic factors responsible for this tumor regression property, whole genome resequencing was conducted using Illumina Hi-Seq 2 × 100 bp paired-end read method (San Diego, CA, USA) with AR (confirmed tumor regression property) and AP chickens. Sequence reads were aligned to the chicken reference genome (galgal5) and produced coverage of 11× and 14× in AR and AP, respectively. A total of 7.1 and 7.3 million single nucleotide polymorphisms (SNPs) were present in AR and AP genomes, respectively. Through a series of filtration processes, a total of 12,242 SNPs were identified in AR chickens that were associated with non-synonymous, frameshift, nonsense, no-start and no-stop mutations. Further filtering of SNPs based on read depth ≥ 10, SNP% ≥ 0.75, and non-synonymous mutations identified 63 reliable marker SNPs which were chosen for gene network analysis. The network analysis revealed that the candidate genes identified in AR chickens play roles in networks centered to ubiquitin C (UBC), phosphoinositide 3-kinases (PI3K), and nuclear factor kappa B (NF-kB) complexes suggesting that the tumor regression property in AR chickens might be associated with ubiquitylation, PI3K, and NF-kB signaling pathways. This study provides an insight into genetic factors that could be responsible for the tumor regression property.

MicroRNA profiling associated with muscle growth in modern broilers compared to an unselected chicken breed.

BACKGROUND: Genetically selected modern broiler chickens have acquired outstanding production efficiency through rapid growth and improved feed efficiency compared to unselected chicken breeds. Recently, we analyzed the transcriptome of breast muscle tissues obtained from modern pedigree male (PeM) broilers (rapid growth and higher efficiency) and foundational Barred Plymouth Rock (BPR) chickens (slow growth and poorer efficiency). This study was designed to investigate microRNAs that play role in rapid growth of the breast muscles in modern broiler chickens. RESULTS: In this study, differential abundance of microRNA (miRNA) was analyzed in breast muscle of PeM and BPR chickens and the results were integrated with differentially expressed (DE) mRNA in the same tissues. A total of 994 miRNA were identified in PeM and BPR chicken lines from the initial analysis of small RNA sequencing data. After filtering and statistical analyses, the results showed miR-2131-5p, miR-221-5p, miR-126-3p, miR-146b-5p, miR-10a-5p, let-7b, miR-125b-5p, and miR-146c-5p up-regulated whereas miR-206 down-regulated in PeM compared to BPR breast muscle. Based on inhibitory regulations of miRNAs on the mRNA abundance, our computational analysis using miRDB, an online software, predicated that 118 down-regulated mRNAs may be targeted by the up-regulated miRNAs, while 35 up-regulated mRNAs appear to be due to a down-regulated miRNA (i.e., miR-206). Functional network analyses of target genes of DE miRNAs showed their involvement in calcium signaling, axonal guidance signaling, and NRF2-mediated oxidative stress response pathways suggesting their involvement in breast muscle growth in chickens. CONCLUSION: From the integrated analyses of differentially expressed miRNA-mRNA data, we were able to identify breast muscle specific miRNAs and their target genes whose concerted actions can contribute to rapid growth and higher feed efficiency in modern broiler chickens. This study provides foundation data for elucidating molecular mechanisms that govern muscle growth in chickens.

Whole genome resequencing identifies the CPQ gene as a determinant of ascites syndrome in broilers.

BACKGROUND: Ascites syndrome is the most severe manifestation of pulmonary hypertension in fast-growing broilers. The disease can be attributed to increased body weights of birds, where the higher metabolic load is not matched by sufficient oxygen supply to the cells and tissues. Although there are environmental components, the disease exhibits moderate to high heritability. The current study uses high throughput whole genome resequencing (WGR) to identify genes and chromosomal regions associated with ascites. RESULTS: The WGR data identified the CPQ gene on chromosome 2. The association was confirmed by genotyping a large collection of DNAs from phenotyped birds from three distinct broiler lines using SNPs in intron 6 and exon 8 of the CPQ gene. By combining the genotype data for these two SNP loci, we identified three different alleles segregating in the three broiler lines. Particular genotypes could be associated with resistance to ascites. We further determined that particular genotypes most associated with resistance overexpress CPQ mRNA in three tissues which might explain the role of these alleles in contributing to resistance. CONCLUSIONS: Our findings indicate CPQ is an important determinant of pulmonary hypertension syndrome leading to ascites in broilers. We identified particular SNPs that can be used for marker-assisted selection of broilers for resistance to the disease. Our findings validate WGR as a highly efficient approach to map determinants contributing to complex phenotypic or disease-related traits. The CPQ gene has been associated with pulmonary hypertension in genome-wide association studies in humans. Therefore, ascites investigations in broilers are likely to provide insights into some forms of hypertension in humans.

Identification and Differential Abundance of Mitochondrial Genome Encoding Small RNAs (mitosRNA) in Breast Muscles of Modern Broilers and Unselected Chicken Breed.

Background: Although small non-coding RNAs are mostly encoded by the nuclear genome, thousands of small non-coding RNAs encoded by the mitochondrial genome, termed as mitosRNAs were recently reported in human, mouse and trout. In this study, we first identified chicken mitosRNAs in breast muscle using small RNA sequencing method and the differential abundance was analyzed between modern pedigree male (PeM) broilers (characterized by rapid growth and large muscle mass) and the foundational Barred Plymouth Rock (BPR) chickens (characterized by slow growth and small muscle mass). Methods: Small RNA sequencing was performed with total RNAs extracted from breast muscles of PeM and BPR (n = 6 per group) using the 1 × 50 bp single end read method of Illumina sequencing. Raw reads were processed by quality assessment, adapter trimming, and alignment to the chicken mitochondrial genome (GenBank Accession: X52392.1) using the NGen program. Further statistical analyses were performed using the JMP Genomics 8. Differentially expressed (DE) mitosRNAs between PeM and BPR were confirmed by quantitative PCR. Results: Totals of 183,416 unique small RNA sequences were identified as potential chicken mitosRNAs. After stringent filtering processes, 117 mitosRNAs showing >100 raw read counts were abundantly produced from all 37 mitochondrial genes (except D-loop region) and the length of mitosRNAs ranged from 22 to 46 nucleotides. Of those, abundance of 44 mitosRNAs were significantly altered in breast muscles of PeM compared to those of BPR: all mitosRNAs were higher in PeM breast except those produced from 16S-rRNA gene. Possibly, the higher mitosRNAs abundance in PeM breast may be due to a higher mitochondrial content compared to BPR. Our data demonstrate that in addition to 37 known mitochondrial genes, the mitochondrial genome also encodes abundant mitosRNAs, that may play an important regulatory role in muscle growth via mitochondrial gene expression control.