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

Nutritional Strategies to Improve Meat Quality and Composition in the Challenging Conditions of Broiler Production: A Review

Poultry meat is becoming one of the most important animal protein sources for human beings in terms of health benefits, cost, and production efficiency. Effective genetic selection and nutritional programs have dramatically increased meat yield and broiler production efficiency. However, modern practices in broiler production result in unfavorable meat quality and body composition due to a diverse range of challenging conditions, including bacterial and parasitic infection, heat stress, and the consumption of mycotoxin and oxidized oils. Numerous studies have demonstrated that appropriate nutritional interventions have improved the meat quality and body composition of broiler chickens. Modulating nutritional composition [e.g., energy and crude protein (CP) levels] and amino acids (AA) levels has altered the meat quality and body composition of broiler chickens. The supplementation of bioactive compounds, such as vitamins, probiotics, prebiotics, exogenous enzymes, plant polyphenol compounds, and organic acids, has improved meat quality and changed the body composition of broiler chickens

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

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

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

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

Primers, probes, and conditions for RT-qPCR.

<p>For each <i>CPQ</i> locus: position is the [chromosome]:[chromosomal position] of the 5’ base according to the <i>Gallus_gallus-5</i> assembly; Primers are 5’-3’ for forward (F) and reverse (R); Probes are 5’-3’ with allele 1 (P1) labeled with FAM and allele 2 (P2) labeled with HEX. Annealing temperatures for qPCR are indicated in °C.</p

<i>CPQ</i> coding and 3’-UTR SNPs.

<p><i>CPQ</i> coding and 3’-UTR SNPs.</p

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 &alpha;-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 &beta;-oxidation process. Our unbiased transcriptomics approach and validation experiments suggested that IDH2 might play a key role in homeostasis of lipid metabolism

Scatterplot representing GBS data on chromosome 2.

<p>Difference in SNP frequencies (Y-axis) is plotted between the average non-reference SNP frequencies of the two replicate pools (resistant less susceptible) relative to chromosome 2 position (X-axis) in Mbp. Positive values indicate that non-reference SNP favors resistance, whereas negative values indicate the SNP favors susceptibility for ascites. <b>(Panel 1A):</b> Entire chromosome 2 for males. The black arrow indicates the peak associated with the <i>CPQ</i> gene locus. <b>(Panel 1B):</b> An expanded view of the peak region ((black arrow) from Panel 1A. <b>(Panel 1C):</b> The same region as Panel 1B only for females.</p