In-Depth Duodenal Transcriptome Survey in Chickens with Divergent Feed Efficiency Using RNA-Seq

<div><p>Since the feed cost is a major determinant of profitability in poultry industry, how to improve feed efficiency through genetic selection is an intriguing subject for breeders and producers. As a more suitable indicator assessing feed efficiency, residual feed intake (RFI) is defined as the difference between observed and expected feed intake based on maintenance and growth. However, the genetic mechanisms responsible for RFI in chickens are still less well appreciated. In this study, we investigated the duodenal transcriptome architecture of extreme RFI phenotypes in the six brown-egg dwarf hens (three per group) using RNA sequencing technology. Among all mapped reads, an average of 75.62% fell into annotated exons, 5.50% were located in introns, and the remaining 18.88% were assigned to intergenic regions. In total, we identified 41 promising candidate genes by differential expression analysis between the low and high RFI groups. Furthermore, qRT-PCR assays were designed for 10 randomly chosen genes, and nine (90.00%) were successfully validated. Functional annotation analyses revealed that these significant genes belong to several specific biological functions related to digestibility, metabolism and biosynthesis processes as well as energy homeostasis. We also predicted 253 intergenic coding transcripts, and these transcripts were mainly involved in fundamental biological regulation and metabolism processes. Our findings provided a pioneering exploration of biological basis underlying divergent RFI using RNA-Seq, which pinpoints promising candidate genes of functional relevance, is helpful to guide future breeding strategies to optimize feed efficiency and assists in improving the current gene annotation in chickens.</p></div

The percentage of reads mapped to exonic, intronic and intergenic regions.

<p>The percentage of reads mapped to exonic, intronic and intergenic regions.</p

Histogram presentation of gene ontology (GO) term for putative coding transcripts.

<p>The GO terms were classified into different categories at level 2.</p

Detailed information of differentially expressed genes responsible for divergent RFI.

<p>Abbreviations: Chr = chromosome; FC = fold change; FDR = false discovery rate; FCR = feed conversion ratio; RFI = residual feed intake; FI = feed intake.</p><p>^aIdentification of the gene according to Ensembl genes database 76</p><p>Detailed information of differentially expressed genes responsible for divergent RFI.</p

Summary statistics for sequence quality and alignment information of six samples.

<p>Q20 and Q30 represent the proportion of bases with a Phred quality score greater than 20 and 30, respectively.</p><p>Summary statistics for sequence quality and alignment information of six samples.</p

Volcano plot reporting <i>P</i> values against fold changes.

<p>The Volcano plot indicates-log₁₀ (<i>P</i>-value) for genome-wide genes (Y-axis) plotted against their respective log₂ (fold change) (X-axis). The red and blue dots represent significantly up- and down-regulated genes between the low and high residual feed intake groups respectively.</p

Sequence variants (non-reference alleles) in the coding region of <i>LRMDA</i>, <i>ZNF503</i>, and <i>COMTD1</i> detected by sequencing one homozygous <i>IG/IG</i> bird and the corresponding allele frequencies among non-IG <i>(N/-</i>) birds.

Sequence variants (non-reference alleles) in the coding region of LRMDA, ZNF503, and COMTD1 detected by sequencing one homozygous IG/IG bird and the corresponding allele frequencies among non-IG (N/-) birds.</p

Primer sequences used in this paper.

The biochemical pathway regulating the synthesis of yellow/red pheomelanin is less well characterized than the synthesis of black/brown eumelanin. Inhibitor of gold (IG phenotype) is a plumage colour variant in chicken that provides an opportunity to further explore this pathway since the recessive allele (IG) at this locus is associated with a defect in the production of pheomelanin. IG/IG homozygotes display a marked dilution of red pheomelanin pigmentation, whilst black pigmentation (eumelanin) is only slightly affected. Here we show that a 2-base pair insertion (frame-shift mutation) in the 5th exon of the Catechol-O-methyltransferase containing domain 1 gene (COMTD1), expected to cause a complete or partial loss-of-function of the COMTD1 enzyme, shows complete concordance with the IG phenotype within and across breeds. We show that the COMTD1 protein is localized to mitochondria in pigment cells. Knockout of Comtd1 in a mouse melanocytic cell line results in a reduction in pheomelanin metabolites and significant alterations in metabolites of glutamate/glutathione, riboflavin, and the tricarboxylic acid cycle. Furthermore, COMTD1 overexpression enhanced cellular proliferation following chemical-induced transfection, a potential inducer of oxidative stress. These observations suggest that COMTD1 plays a protective role for melanocytes against oxidative stress and that this supports their ability to produce pheomelanin.</div

HA-tagged COMTD1 localizes to mitochondria in immortalized mouse melanocytes.

(A-D) Immortalized melan-Ink4a cells from Ink4a-deficient C57BL/6J mice were transiently transfected to express COMTD1 fused with the HA11 epitope at either the N-terminus (HA-COMTD1; A, C) or C-terminus (COMTD1-HA; B, D). Two days later, cells were fixed and analyzed by bright field (BF) and immunofluorescence microscopy for HA and either the mitochondrial resident protein MAVS (A, B) or the ER resident protein calnexin (CNX; C, D). Individual images of labelled cells or the bright field image are shown in addition to an overlay of HA (green) with MAVS (red; HA/ MAVS), CNX (red; HA/ CNX), or the pseudocolored bright field image (magenta; HA/BF). Insets show a 5-fold magnified image of the boxed region to emphasize overlap or lack thereof. Main scale bar, 10 μm; inset scale bar, 2 μm. (E) Quantification of the degree of overlap of COMTD1-HA or HA-COMTD1, as indicated, with markers of the ER (CNX; N = 29 for COMTD1-HA, N = 17 for HA-COMTD1), mitochondria (MAVS; N = 25 for COMTD1-HA, N = 16 for HA-COMTD1), mature melanosomes (TYRP1; N = 16), immature melanosomes (PMEL; N = 17), late endosomes/ lysosomes (LAMP2; N = 15), or early endosomes (STX13; N = 21). Data from 4–5 individual experiments are presented as a box and whiskers plot in which the area of overlap is shown relative to the total area occupied by HA (e.g., CNX vs. HA) or by the indicated marker (e.g., HA vs. CNX). See S2 Fig for examples of the data for TYRP1, PMEL, LAMP2 and STX13. Statistical significance was determined by ordinary one-way ANOVA with Tukey’s tests for multiple comparisons; ****, P < 0.0001.</p

Genotype distribution of the 2-bp insertion in <i>COMTD1</i> associated with the Inhibitor of Gold (IG) phenotype in different populations sorted by phenotype.

Genotype distribution of the 2-bp insertion in COMTD1 associated with the Inhibitor of Gold (IG) phenotype in different populations sorted by phenotype.</p