#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Combined transcriptomics and proteomics forecast analysis for potential genes regulating the Columbian plumage color in chickens


Authors: Xinlei Wang aff001;  Donghua Li aff001;  Sufang Song aff002;  Yanhua Zhang aff001;  Yuanfang Li aff001;  Xiangnan Wang aff001;  Danli Liu aff001;  Chenxi Zhang aff001;  Yanfang Cao aff001;  Yawei Fu aff001;  Ruili Han aff001;  Wenting Li aff001;  Xiaojun Liu aff001;  Guirong Sun aff001;  Guoxi Li aff001;  Yadong Tian aff001;  Zhuanjian Li aff001;  Xiangtao Kang aff001
Authors place of work: College of Animal Science and Veterinary Medicine, Henan Agricultural University, Zhengzhou, Henan, China aff001;  College of Animal Science and Technology, Henan University of Animal Husbandry and Economy, Zhengzhou, Henan, China aff002
Published in the journal: PLoS ONE 14(11)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0210850

Summary

Background

Coloration is one of the most recognizable characteristics in chickens, and clarifying the coloration mechanisms will help us understand feather color formation. “Yufen I” is a commercial egg-laying chicken breed in China that was developed by a three-line cross using lines H, N and D. Columbian plumage is a typical feather character of the “Yufen I” H line. To elucidate the molecular mechanism underlying the pigmentation of Columbian plumage, this study utilizes high-throughput sequencing technology to compare the transcriptome and proteome differences in the follicular tissue of different feathers, including the dorsal neck with black and white striped feather follicles (Group A) and the ventral neck with white feather follicles (Group B) in the “Yufen I” H line.

Results

In this study, we identified a total of 21,306 genes and 5,203 proteins in chicken feather follicles. Among these, 209 genes and 382 proteins were differentially expressed in two locations, Group A and Group B, respectively. A total of 8 differentially expressed genes (DEGs) and 9 differentially expressed proteins (DEPs) were found to be involved in the melanogenesis pathway. Additionally, a specifically expressed MED23 gene and a differentially expressed GNAQ protein were involved in melanin synthesis. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis mapped 190 DEGs and 322 DEPs to 175 and 242 pathways, respectively, and there were 166 pathways correlated with both DEGs and DEPs. 49 DEPs/DEGs overlapped and were enriched for 12 pathways. Transcriptomic and proteomic analyses revealed that the following pathways were activated: melanogenesis, cardiomyocyte adrenergic, calcium and cGMP-PKG. The expression of DEGs was validated by real-time quantitative polymerase chain reaction (qRT-PCR) that produced results similar to those from RNA-seq. In addition, we found that the expression of the MED23, FZD10, WNT7B and WNT11 genes peaked at approximately 8 weeks in the “Yufen I” H line, which is consistent with the molting cycle. As both groups showed significant differences in terms of the expression of the studied genes, this work opens up avenues for research in the future to assess their exact function in determining plumage color.

Conclusion

Common DEGs and DEPs were enriched in the melanogenesis pathway. MED23 and GNAQ were also reported to play a crucial role in melanin synthesis. In addition, this study is the first to reveal gene and protein variations in in the “Yufen I” H line during Columbian feather color development and to discover principal genes and proteins that will aid in functional genomics studies in the future. The results of the present study provide a significant conceptual basis for the future breeding schemes with the “Yufen I” H line and provide a basis for research on the mechanisms of feather pigmentation.

Keywords:

Gene expression – Gene regulation – Proteomes – Transcriptome analysis – Melanin – RNA sequencing – Chickens – Feathers

Introduction

“Yufen I” is an egg-laying chicken breed in China, which was developed by a three-line cross using line H as the first male parent, line N as the maternal grandparent and line D as the final male parent. As authorized by the National Commission on Livestock and Poultry Genetic Resources in 2015, this breed has been bred true for at least six generations. However, a closed breeding method was used to develop the line H by crossing the barred plumaged-original Gushi chicken with an egg-laying grandparent line C, the brown-shelled Babcock B-380. Line H is a fast plumage line and is characterized by early maturity, high egg production and the Columbian plumage pattern. Columbian plumage is a character of feather color in the “Yufen I” H line in which the dorsal neck, tail and apex of the wing feathers have black and white stripes and other feathers are white. This pattern was also regarded as the main pigmentation character of the “Yufen I” H line.

Complex feather coloration is likely coordinated through multiple genes that regulate diverse mechanisms, with more than 200 genes involved in pigmentation that have been studied in mammals [1]. Many studies have identified key genes that determine the types of pigments (melanin) that are expressed by melanocytes [2]. Birds are among the most colorful vertebrates. Melanins, porphyrins, polyenes, carotenoids and structural colors have been discovered in feathers [3]. Feather color in chickens is a result of the melanin produced by the melanocytes of the feather follicles. Feathers and feather follicles are ideal tissues to explore the genetic mechanisms and complexity of color patterns in birds. As a derivative of chicken skin, the feather follicles give rise to the feathers and are capable of self-renewal, and their proliferation and differentiation result in feather formation [46]. In adult feather follicles, plumage pigmentation is mainly dependent on the interaction between feather follicle melanocytes and dermal papilla fibroblasts. Pigmentation activity occurs only during the growth period of feather follicles, and the transfer of melanin and pigment to keratinocytes depends on the activity of melanin precursors. The melanin or pigment is transferred to the skin and feather follicles through the regulation of signal transduction pathways [7,8]. Studies have revealed that many growth factors and receptors coordinate genes and that the environment and signaling pathways play an extremely important role during feather growth. To date, some of the confirmed pathways involved in pigmentation include cAMP pathway [9], SCF-KIT pathway [10], PI3K-Akt pathway [11], BMP signaling pathway [12], Notch pathway [13], ERK pathway [14], CREB/MITF/tyrosinase pathway [15], MCIR/(Gs-AC)/PKA pathway [1618], Wnt/β-catenin pathway [19,20] and MAPK pathway [21,22].

Uniformity in the appearance of birds is essential in the poultry industry [23]. Although plumage color is easily observed, the genetics behind feather pigmentation are governed by both qualitative and quantitative features [24]. Research reveals that the ratio of pheomelanin (yellow-red pigment) and eumelanin (black-brown pigment) pigmentation regulates feather color in chickens [25,26]. Eumelanin pigmentation requires that the melanoblasts migrate to the epidermis from the neural crest to ultimately reach the developing feather follicles. Another requirement is the proportion of the pigment subject to control by genes [27,28]. Melanosomes are responsible for synthesizing these pigments. These organelles are granule-like and develop within melanoblasts, which go through several steps of differentiation to form melanocytes. The preliminary steps for melanin production include the appearance of the neural crest, determination of melanoblasts, migration, proliferation and differentiation [2931]. Besides this, melanogenesis regulation after the melanocytes migrate into feather follicles is another important step. The variation in plumage may be caused by any mutations in relevant genes and changes in molecules (receptors of transcription factors on cells, structural proteins, enzymes and growth factors) involved in the aforementioned process [32,17,33]. Until now, several studies have focused on genes such as MC1R, ASIP, TYR, SLC24A5, KITLG, MITFCDKN2A/B, PMEL17 and DCT, which play a role in melanin proportion synthesis [34,25,3541].

Many genes and pathways have been shown to be associated with pigmentation. However, research on the genetics of Columbian plumage in chickens is lacking. Feathers have different forms in terms of color and morphology, not only among different bird species but also among different body regions of an individual bird [42,43]. Some studies used transcriptomics to analyze the skin and hair follicles of poultry and found differentially expressed genes (DEGs) related to feather color and morphology, which provided a theoretical basis for studying the formation, development and regeneration of feathers [44,45]. In this experiment, we collected the feather follicles in two parts: the dorsal feather follicles of the neck with black and white-striped feathers (Group A) and the ventral feather follicles of the neck with white feathers (Group B). This study aimed to discover the differentially expressed genes (DEGs), differentially expressed proteins (DEPs) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways by transcriptomic (RNA-seq) and proteomic (iTRAQ) analyses and to forecast potential candidate genes of Columbian plumage in order to determine which genes and pathways affect feather coloration.

Materials and methods

Disclosure of ethics

Experimental animals were maintained per the rules in National standards for the environment and facilities of experimental animals of China (GB14925-2010). All the chickens were healthy, with a coop size of 0.12m2 and 0.4m high for each. Adequate and clean drinking water and feed were provided. All protocols for animal experiments received approval from the Institutional Animal Care and Use Committee (IACUC) of Henan Agricultural University and from Henan Agricultural University’s Animal Care Committee, College of Animal Science and Veterinary Medicine, China (Permit Number: 17–0322).

Animals for experiments and tissue sources

The Animal Center of Henan Agricultural University provided us with chickens of the “Yufen I” H line breed for this study. Three 22-week-old female chickens were anaesthetized with “Su-Mian-Xin” (a type of anesthetic, Shengda Animal Pharmaceutical Co., Ltd, Dunhua, China), which is diluted with saline at a volume ratio of 1:2. The chickens were anesthetized by i.v. injection in the wing vein with a concentration of 0.2mL/kg for approximately 20 minutes, and then the neck feathers were plucked. New feathers emerged in the skin after two weeks. Next, the three chickens involved in this study were placed in an airtight box and humanely sacrificed by inhaling carbon dioxide in order to reduce their suffering. Six tissue samples including feather follicles and the circumjacent skin closely around the rachis of the feathers were collected from the dorsal and ventral areas, the locations at which black and white-striped feathers and white feathers occur, respectively (Fig 1). Three replicates were collected for each group (A1 to A3 for Group A and B1 to B3 for Group B). Approximately 20~30 feather follicles were placed into 2-mL tubes, and then the tubes were sealed, dipped in liquid nitrogen to freeze quickly and stored at -80°C storage isolation and sequencing of RNA and qRT-PCR analysis [46].

Fig. 1. “Yufen I” H line chicken and the feather bulbs used in this study.
“Yufen I” H line chicken and the feather bulbs used in this study.
(A). “Yufen I” H line chicken, (B). black and white striped-feathers in the dorsal neck, (C). white feathers in the ventral neck.

Extraction of total RNA and RNA-seq

Six frozen samples corresponding to the two feather follicle areas (Group A and Group B) were selected for isolation of total RNA. The feather follicles of the “Yufen I” H line were used to extract total RNA using TRIzol reagent (Invitrogen, USA) and used in library construction. The assessment of contamination and degradation of RNA were performed with standard denaturing agarose gel electrophoresis. The RNA integrity, concentration and purity were evaluated on an Agilent RNA6000 Nano Chip in Reagent Port 1 of the Bioanalyzer Agilent 2100 (Agilent Technologies, CA, USA). Sequencing of the libraries was performed on the Illumina HiSeq 4000 platform by BGI Co., Ltd. Data were deposeted in the NCBI Sequence Read Archive under Accession SRR7973871. To process the raw data, the FASTQ format was first processed through SOAPnuke (v1.5.2), and then reads carrying adapters, poly-N sequences and those of low quality were deleted from the raw data to obtain clean data. Additionally, we calculated the Q20 and Q30 at error rates of 1% and 0.1%, respectively, for the clean data. To ascertain whether resequencing was required, quality control (QC) for alignment was carried out. These high quality data were used for the analyses performed downstream.

Transcriptomic data processing

After read filtering, we used HISAT (v0.1.6-beta) to perform genome mapping. To map RNA-seq reads, a spliced alignment program, HISAT, is fast and sensitive with an accuracy that is equal to or better than that of other methods. The spliced mapping algorithm of HISAT has been applied for genome mapping of preprocessed reads [47,48]. Reads that passed the QC test were aligned to a reference genome assembly of chicken (Gallus gallus 5.0, https://www.ncbi.nlm.nih.gov/assembly/GCF_000002315.4/) from NCBI. Measurement of the abundance of expression for each assembled transcript was done using the Fragments per Kilobase of exon model per Million mapped reads (FPKM) values.

For the two groups, to analyze differential expression, the levels of gene and transcript expression were determined using RSEM software (V1.2.12). For each pair of samples, the DEGs were screened using a model derived from PossionDis software [49]. In this paper, “false discovery rate (FDR) ≤ 0.001 and the fold change ≥ 2 (absolute value of log2Ratio ≥ 1)” were set as the significant threshold values to ascertain the differences between the gene expression. The phyper function of R software was used to identify enriched KEGG (http://www.kegg.jp/) pathways (p ≤ 0.05) and Gene Ontology (GO: http://www.geneontology.org), respectively. Additionally, to identify the major metabolic and signal transduction pathways, KEGG enrichment analysis was performed on the DEGs. Likewise, GO enrichment analysis was performed to ascertain the main molecular functions, cellular components and biological processes associated with DEGs.

RNA validation through qRT-PCR

To validate the RNA-seq data, qRT-PCR was used, and seven DEGs were selected for the analysis. The same samples used for RNA sequencing were used for reverse transcription and synthesis of cDNA with the PrimeScript RT Reagent Kit with gDNA Eraser following the instructions of the manufacturer (TaKaRa, Dalian, China). To design the primers for qRT-PCR, the NCBI Primers-BLAST online program was used. Information regarding the primers of these genes can be found in S1 Table. Each 10 μL qRT-PCR reaction contained 1.0 μL of cDNA, 0.5 μL of each primer at 10 μM, 5.0 μL of 2×SYBR® Premix Ex Taq™ II (TaKaRa, Dalian, China), and 3 μL of deionized water. The reaction was carried out on a LightCycler® 96 Real-Time PCR system (Roche Applied Science, Indianapolis, USA).The internal control was the GAPDH gene, and the 2-ΔΔCT method was used to assess expression. The procedure for qRT-PCR amplification is shown: 95°C for 3 min; 35 cycles at 95°C for 30 s, 60°C for 30 s, and 72°C for 20 s; and a final extension for 10 min at 72°C. The data were statistically analyzed using SPSS V 21.0 (SPSS Inc., Chicago, IL, USA). One-way and repeated- analyses of variance followed by Dunnett’s test were carried out. The data are shown in the form of the mean ± SE with significance set at p≤ 0.05.

Protein extraction and iTRAQ reagent labeling

In the context of several experiments, an approach called isobaric tags for relative and absolute quantification (iTRAQ) has been successfully used. Here, the samples were used for RNA-seq as well as iTRAQ analysis. For protein extraction, the lysis of 2 g of each of the samples was carried out in lysis buffer 3 (TEAB with 1 mM PMSF, 8 M Urea, 10 mM DTT and 2 mM EDTA, pH 8.5) with 2 magnetic beads of 5 mm diameter. The samples were kept in a tissue lyser for 2 min at 50 Hz so that the proteins were released and then subjected to centrifugation at 25,000 x g for 20 min at 4 C, The the supernatant was transferred into a fresh tube, and 10 mM DTT (dithiothreitol) was added for reduction at 56°C for 60 mins, followed by alkylation using 55 mM IAM (iodoacetamide) at room temperature, in the dark for 45 min, and then centrifuged under the same conditions described above.

The concentration and quality of proteins were assessed with a Bradford assay of the supernatant and confirmed with 12% SDS-PAGE.Then, 100 mM TEAB was used to dilute 100μg of protein solution with 8M urea that was subjected to digestion using Trypsin Gold (40:1 protein: enzyme, Promega, Madison, WI, USA) at 37°C overnight. Subsequently, an iTRAQ Reagent 8-plex Kit was used to label samples, followed by combining the differently labeled peptides, desalting using a Strata X C18 column (Phenomenex), and finally vacuum drying in accordance with the instructions in the manufacturer's protocol. After this step, fractionation of these pooled mixtures was performed on a Shimadzu LC-20AB high pressure liquid chromatography (HPLC) pump system with a high pH RP column, followed by nanoelectrospray ionization and subsequent tandem mass spectrometry (MS/MS) on a Triple TOF 5600 system (SCIEX, Framingham, MA, USA). Triplicate experiments were carried out.

Proteomics database processing

The ProteoWizard tool was utilized to convert raw MS/MS data into “.mgf” format followed by searching these exported files with Mascot (v 2.3.02) in this project against a database (Gallus gallus 5.0, https://www.ncbi.nlm.nih.gov/assembly/GCF_000002315.4/). Identification required the presence of a minimum of one unique peptide. Parameters set included: MS/MS ion search; Trypsin enzyme; 0.1 Da was the mass tolerance of fragment; Monoisotopic mass values; Variable modifications oxidation (M) and iTRAQ8plex (Y); 0.05 Da was the tolerance for Peptide mass; Fixed modifications Carbamidomethyl (C), iTRAQ8plex (N-term), iTRAQ8plex (K); Database I-ZAwBa007 (50596 sequences); Database_info transcriptome. Identification was performed using peptides that reached a confidence level of 95%.

The iTRAQ peptides were analyzed in a quantitative format with IQuant automated software [50]. This software incorporates a method based on machine learning called Mascot Percolator [51], which can rescore results from databases to yield a significant scale of standards. A 1% FDR was used to pre-filter PSMs that were used to measure the confidence level. A set of confident proteins was assembled using the parsimony principle also termed as the “simple principle”. Using the picked protein FDR strategy as a basis [52], followed by protein inference with an FDR at a protein level lesser than or equal to 0.01, an FDR of 1% was measured in order to limit the number of false positives. The quantification of proteins involved identification of proteins, correction of tag impurities, normalization of data, imputation of missing values, calculation of protein ratio, and statistical approaches followed by presentation of the results.

The proteins that had a p-value lower than 0.05 and a fold change of 1.2 were classified as DEPs. Analysis of metabolic pathways was in accordance with KEGG, while analysis of other databases, COG and GO, was in lieu of earlier research [53]. The hypergeometric test was applied to analyze DEP enrichment analysis in KEGG and GO. The equation was as follows:

Where N represents the quantity of identified proteins that were linked to data from GO and KEGG analyses, n represents the measure of DEPs within N, M represents the quantity of proteins linked to a pathway or term of GO or KEGG, and m represents the quantity of DEPs linked to a pathway or term of GO or KEGG. An enrichment of differential proteins was considered when the p-value ≤ 0.05. The DEG analysis was the same as the DEP analysis.

Transcriptomics and proteomics: Association analysis

Genes are regulated at multiple levels during the expression process. At present, most studies have reported that the expression consistency between mRNAs and their corresponding proteins is not very high, so a combined analysis of the proteome and transcriptome is helpful for discovering the regulation of gene expression [54]. Cluster 3.0 software was utilized to analyze clusters for DEP expression with that of transcripts in order to recognize transcripts and DEPs with similarity across various tissues with a graphical output from Java Treeview software. All expression data related to proteomics and transcriptomics were analyzed, and the Spearman correlation coefficient was calculated using R software [55]. Associations between the expression of mRNA and proteins from these respective “omics” were quantified followed by analysis for GO as well as KEGG enrichmentin order to study the potential role of DEGs/DEPs in metabolism or signal transduction. The combined transcriptomic and proteomic analysis parameters are given in S2 Table.

Results

Analysis of data from RNA sequencing

The Illumina HiSeq 4000 platform was employed to sequence 6 preparations of cDNA library. The RNA-seq data showed high consistency among the libraries (Table 1). A total of 44.06 MB of clean reads and 45.05 MB of clean reads were obtained from the 6 cDNA libraries. The Q30 for each sample was more than 94%, and approximately 71% of the clean reads were mapped to the chicken genome (Gallus gallus 5.0).The data are indicative of data of high quality sequencing which can ensure the reliability of results for further investigation.

Tab. 1. Gene expression and clean reads analyses in “Yufen I” H line feather follicles.
Gene expression and clean reads analyses in “Yufen I” H line feather follicles.

Identification of DEGs

To identify genes involved in feather coloration, DEGs were detected in the “Yufen I” H line using PossionDis software [56]. Of a total of 21,306 identified genes, 597 novel genes and 12,687 novel isoforms were identified in chicken feather follicle libraries in this work. A total of 209 (83 upregulated and 126 downregulated) DEGs were detected at the two locations (the dorsal and ventral areas, where black and white-striped feathers and white feathers occur). All of the DEGs are illustrated in S3 Table. Earlier research showed several genes involved in pigmentation (e.g., ASIP, KITLG) [57,58].

In order to corroborate whether the RNA-seq gene expression data were accurate and reproducible, 7 genes were subjected to selection for qRT-PCR in the two groups. The expression of KITLG, FZD10, WNT7B, WNT9A, WNT11, PVALB and MED23 was validated by qRT-PCR (Fig 2). As the results from both analyses were consistent with each other, the reliability of the sequencing was indicated.

Fig. 2. Verification of DEGs via qRT-PCR.
Verification of DEGs via qRT-PCR.
The 2-ΔΔCT method was used for data analysis, and the housekeeping gene was GAPDH. Data shown on the vertical-axis represents the relative expression. Significant differences between two groups were determined by applying the unpaired Student’s t-test. * 0.01 < p < 0.05, ** p < 0.01. All dataare presented as means ± standard error (SE).

The downy feathers of chickens are reportedly replaced by the second generation at the age of 6 weeks. At the age of 8 weeks, the second generation of feathers begins to molt and a large number of third-generation feathers begin to appear. At this time, the feather color tends to be stable. We collected feather follicles at 0, 2, 4, 6, 8, 10 and 12 weeks from the dorsal neck of the “Yufen I” H line chicken. Using qRT-PCR, we found that MED23, FZD10, WNT7B and WNT11 gene expression peaked at approximately 8 weeks, which is consistent with the molting cycle (Fig 3). Furthermore, FZD10, WNT7B and WNT11 are located in the melanogenesis pathway, and MED23 is specifically expressed in the dorsal feather follicles of the neck with black and white-striped feathers.

Fig. 3. qRT-PCR validation of DEGs involved in the feather cycle.
qRT-PCR validation of DEGs involved in the feather cycle.
Changes in DEGs in feather follicles from the dorsal neck of the chicken at 0, 2, 4, 6, 8, 10 and 12 weeks.

Metabolic pathways and GO analysis of DEGs

KEGG pathway and GO analyses were performed in order to interpret the exact roles played by the DEGs that regulate feather pigmentation. Using Blast2GO, we classified the DEGs with GO terms and KEGG pathways according to their functions. With the KEGG and GO annotation results, official classification of DEGs was performed followed by functional enrichment study with these databases using phyper, an R software function.

Classification of 209 DEGs with GO terms was performed with the Blast2GO platform to identify functional roles with an abundance of the following terms: “biological processes,” “molecular functions” and “cellular components” (FDR ≤ 0.01). Many genes were associated with cell part, cell, organelle, and cellular process as well as binding. The categories with abundance were the processes associated with cell, single-organism, metabolic process and multicellular organisms, as well as regulation of biological process in the biological processes category. Cells, cell part and organelles showed maximum abundances in the category called cellular component. The majority of the unigenes could be classified into binding and catalytic activity functions in terms of function (Fig 4A).

Fig. 4. KEGG and GO analyses of DEGs.
KEGG and GO analyses of DEGs.
(A) GO functional annotation histogram of the DEGs. The vertical axis represents the three GO categories, the horizontal axis represents the gene number, and the number of genes is considered the difference in the proportion of the total. The GO annotations are classified in three basic categories, including cellular component, biological processes, and molecular function. (B) The degree of enrichment of the first 20 entries in the pathway. The pathway names are represented on the vertical axis, and the horizontal axis represents the pathways corresponding to the rich factor. The ratio of the number of DEGs and all annotated genes in the pathway is defined as the rich factor.

As a result, 190 DEGs were mapped to 175 pathways in KEGG (FDR ≤ 0.01); therein, 188 DEGs were enriched in environmental information processing, of which 81 DEGs were enriched at the level of signaling molecules and interaction and 107 DEGs showed enrichment at the level of signal transduction. The most enriched pathways of the DEGs were focal adhesion, metabolic pathways, PI3K-Akt, cGMP-PKG, calcium and adrenergic signaling pathway in cardiomyocytes showed significant enrichment in two groups (Fig 4B).

iTRAQ data analysis

iTRAQ analysis or proteomics was carried out to supplement the transcriptomics data from the same samples. The range of masses of identified proteins was 10 to 100 kDa, while the average coverage was challenged for groups with more than 100 kDa. These data are indicative of reliable proteomic analyses [59]. This technique aided the identification of 293,356 spectra in the samples which were subjected to data filtering to obtain 51,631 unique spectra that could be matched with 23,244 unique peptides, for the total identification of 5,203 proteins (Table 2).

Tab. 2. Information on identified proteins in chicken feather follicles.
Information on identified proteins in chicken feather follicles.

Functional classification and annotation of DEPs

The DEPs were identified with a p-value < 0.05 and fold change ≥ 1.2 between the dorsal and ventral feather follicles of the neck. In brief, 382 DEPs (160 upregulated and 222 downregulated) were detected in the feather follicles (S4 Table), followed by classification into functions of COG categories. Clearly, changes were observed in the levels of proteins involved in functions such as general function prediction only; posttranslational modification, protein turnover; signal transduction mechanisms; translation, ribosomal structure and biogenesis; transcription and cytoskeleton (Fig 5A). DEPs were subjected to enrichment and clustering using functional analyses of GO and KEGG. In the GO function analysis, cell part, intracellular part, and cytoplasm of cellular component; binding, catalytic activity and ion binding of molecular function; and biological regulation, response to stimulus and metabolic process of biological processes, were found to be different between the samples (Fig 5B). Regarding KEGG pathway functions associated with the feather follicles, significant differences were found in the following signaling pathways: calcium, PI3K-Akt, mTOR, cAMP, MAPK, Wnt, cGMP-PKG, Jak-STAT and adrenergic signaling in cardiomyocytes, melanogenesis, tyrosine metabolism, and melanoma. We displayed the top 20 enriched pathways in Fig 5C. At present, under feather regeneration treatment, segment polarity protein disheveled homolog DVL-3 isoform X6 protein (DVL3), 1-phosphatidylinositol 4,5-bisphosphate phosphodiesterase beta-1-like protein (LOC107052863), 1-phosphatidylinositol 4,5-bisphosphate phosphodiesterase beta-1 protein (PLCB1), GTPase KRas isoform X2 protein (KRAS), matrix-remodeling-associated protein 8 precursor protein (MXRA8), guanine nucleotide-binding protein G(q) subunit alpha protein (GNAQ), parvalbumin muscle isoform X1 protein (PVALB), calcium/calmodulin-dependent protein kinase type II subunit beta protein (CAMK2B) and calcium/calmodulin-dependent protein kinase type II subunit alpha isoform X3 protein (CAMK2A) were differentially expressed in the dorsal follicles of the neck compared with the ventral follicles of the neck and were mainly enriched for the melanogenesis pathway. The results indicated that upon plucking stimulation, the DEPs related to feather pigmentation were mainly concentrated in the melanogenesis pathway.

Fig. 5. DEP analysis in terms of GOG, GO and KEGG.
DEP analysis in terms of GOG, GO and KEGG.
(A) The horizontal axis represents the COG term, and the vertical axis represents the corresponding protein count illustrating the protein number of different function. (B) GO functional annotation histogram of the DEPs. The three GO categories are presented on the vertical axis under the GO term, the horizontal axis represents the gene number, and the number of genes is accounted for by differences in the proportion of the total. The GO annotations are classified in three basic categories including cellular component, biological processes, and molecular function. (C) The name of the pathway is mentioned on the vertical axis, and the pathway matching the rich factor are mentioned on the horizontal axis. The rich factor is the ratio of the number of DEPs in the pathway to the number of all annotated proteins in the pathway. After testing multiple hypotheses, Q values were completed with corrected P value in the range of 0–0.05. The enrichment was considered significant if the P-value was closer to zero.

Integrating transcriptomic and proteomic results

A majority of earlier reports are suggestive of a weak correlation between the expression of mRNA with protein attributed to posttranscriptional regulation or posttranslational modification or experimental errors [6062]. Nevertheless, the flow of information from RNA to proteins is the crux of the central dogma [63,64]. To allow for genes that are expressed differentially in terms of the transcriptome and proteome as well as singling out vital genes, we performed an integration of DEGs and DEPs. For multi-omics data, GO terms and KEGG pathways were enriched at the levels of transcriptome and proteome, respectively. Then, the data from the two groups were integrated and analyzed, which was conducive to the study of gene expression regulation at the level of gene set coexpression [6568].

To compare the proteomic and transcriptomic analyses, we compared the 382 DEPs with the 209 DEGs. According to the results, only 49 genes meeting the criteria overlapped (S5 Table). The changes at the transcript and protein levels showed a weak correlation for the proteins that were quantified. Biological pathways were elucidated with a change in the statistical reports at the protein level when the changes in mRNA were zero. To investigate the overall correlation between these transcripts and DEPs, all identified mRNAs with DEPs were matched, followed by transformation of DEP and transcript volume ratios into log2 forms. An investigation of changes at both the transcript and protein levels revealed a weak correlation (Spearman correlation coefficient, R = 0.0840) for all genes and proteins assessed. We then calculated the correlation between the 382 DEPs and the 209 DEGs, and a positive correlation of R = 0.3006 was obtained when all significantly changed proteins with a cognate mRNA were considered (Fig 6). This result showed a modest correlation between the mRNA and protein levels (between the proteome and transcriptome), which was consistent with previously reported results [69], and accounted for the complexity of the gene expression regulation mechanism [54,7072].

Fig. 6. Correlations between the expression of proteins and genes.
Correlations between the expression of proteins and genes.
The vertical-axis represents the protein expression level, and the horizontal-axis represents the genes expression level. (A) Scatter plots of the correlation between data sets of genes evaluated in both the proteomic and gene transcript analyses. (B) Scatter plots and coefficients of DEPs and DEGs correlation.

Forty-nine DEGs/DEPs were assigned GO terms to assess their functions that encompassed a vast range of cellular components, molecular functions and biological processes (S1 Fig). The GO analysis indicated that the DEGs/DEPs relevant to cellular, metabolic and biological regulation processes are possibly related to our study and showed that the DEGs/DEPs were relevant to molecular functions including catalysis, binding or structure.

Assignment to COG functional categories indicated that excluding carbohydrate transport and metabolism and the cytoskeleton, the DEGs/DEPs were classified into the categories of inorganic ion transport and metabolism and general function prediction.

Additionally, among the components of a single pathway the extent of comity between proteins and mRNA was studied. In total, 382 DEPs and 209 DEGs were mapped to 166 biological pathways, of which 12 pathways showed significant enrichment in both groups. This study involved the determination of vital DEPs and DEGs, involved in melanogenesis and the calcium, cGMP-PKG and adrenergic signaling in cardiomyocytes signaling pathway using pathway enrichment analysis (Fig 7).

Fig. 7. Correlation of KEGG enrichment between transcriptome and proteome.
Correlation of KEGG enrichment between transcriptome and proteome.
(A) Number of KEGG enrichment correlations between transcriptome and proteome. (B) The overview scatter diagram of KEGG enrichment correlations between the transcript levels and protein levels of genes.

Discussion

Feathers have evolved to have diverse coloration that exhibits a wide spectrum of colors arranged in noticeable patterns. Either physical or chemical factors or a combination can form these colors. Melanin is the major factor responsible for the different types of colors described in the existing literature.

DEGs and DEPs related to melanin synthesis

Pigmentation based on melanin is linked to a sequence of genes involved in the initiation and regulation of pigment production involving many signaling pathways as well as transcription factors such as the tyrosine kinase receptor KIT with SCF ligand, and MITF [73]. Many genes in the melanogenesis pathway are related to MITF, which is a sole transcription factor of the microphthalmia family that is involved in the regulation of melanocytes. MITF target genes regulate melanocyte pigmentation [74]. Previous studies revealed that MITF was involved in the development of melanocytes, along with reports of plumage color of Japanese quail and chicken and mutations in this gene [75].

In this study, a total of 8 DEGs (ASIP, KITLG, FZD10, WNT7B, WNT9A, WNT9B, WNT11, PVALB) and 9 DEPs (DVL3, KRAS, MXRA8, GNAQ, PVALB, CAMK2B, CAMK2A, PLCB1, LOC107052863) were found to be involved in the melanogenesis pathway (Fig 8). MED23 was specially expressed and significantly upregulated in dorsal follicles of the neck. A new study revealed that zebrafish pigmentation is regulated by MED23 via a modulation of the function of MITF as an enhancer [76]. GNAQ was significantly upregulated in the dorsal follicles of the neck in this pathway. Research showed that the GNAQ activates MITF via the MAPK pathway, and then affects melanin synthesis [77]. As a consequence, melanin synthesis can be regulated by influencing the expression of the MITF gene, which then regulates the development of feather coloration. In addition, we found that expression of the MED23, FZD10, WNT7B and WNT11 genes peaked at approximately 8 weeks in the “Yufen I” H line, which is consistent with the molting cycle. The results are indicative of an extensive and vital role of these genes in melanin synthesis, which is consistent with previous research results. Therefore, we speculated that the above genes were related to the formation of Columbian plumage in the “Yufen I” H line.

Fig. 8. Differentially expressed feather color genes and proteins identified in the analyzed chicken feather follicles and their involvement in the melanogenesis pathway.
Differentially expressed feather color genes and proteins identified in the analyzed chicken feather follicles and their involvement in the melanogenesis pathway.
DEGs have a blue frame, and the DEPs have a red frame.

Pathways related to melanin synthesis

To further understand the DEGs and DEPs in the pathways that regulate chicken feather color development, we performed KEGG pathway analysis of DEGs and DEPs (p < 0.05). Correlation analysis and integration were performed on DEPs and DEGs annotated in the same pathway. Combining these results with those of the transcriptome and proteome analyses showed thatthe melanogenesis, adrenergic signaling in cardiomyocytes, calcium and cGMP-PKG signaling pathways were highly prevalent at the locations where different feather colors occurred in “Yufen I” H line feather follicles. Therefore, we speculate that these four pathways may be related to the Columbian plumage of the “Yufen I” H line, although we primarily focused on the melanogenesis pathway in this study. We found that the DEGs were mainly concentrated upstream of this pathway, while the DEPs were concentrated downstream of the pathway (Fig 8).

Interestingly, 5 proteins (PLCB1, LOC107052863, PVALB, CAMK2A, CAMK2B) and the PAVLB gene regulated melamine synthesis in the melanogenesis pathway. Few studies have previously reported on the correlation between melamine and pigmentation. The effect of melamine on the skin color of darkbarbel catfish showed that the factors necessary for melanin formation may be suppressed by intake of melamine, resulting in a significant reduction in melanin in the skin of the dorsal surface [78]. Therefore, we speculate that the 5 genes that regulate melamine synthesis may also be related to feather color deposition by affecting melanin synthesis.

Conclusion

In summary, an incorporated and strong approach for analyzing the transcriptome and proteome was utilized to study the mechanism of melanin synthesis in feathers. Strikingly, this original report of a transcriptomics and proteomics analysis of the feather color in chicken follicles describes and reveals a set of DEGs that are putatively involved in determining feather color along with other physiological functions. MED23, FZD10, WNT7B, WNT11 and GNAQ expression were significantly different in feather follicles with black and white stripes versus those with white feather color, and elucidating the functional roles of these genes in the regulation of feather color will be of particular interest in future studies. These results provide a potential understanding of the mechanism underlying Columbian plumage in the “Yufen I” H line and solid genetic resources that allow for the selection of birds with uniform plumage for breeding.

Supporting information

S1 Fig [a]
GO enrichment correlation between transcriptome and proteome.

S1 Table [docx]
Information regarding the specific primers used for the qRT-PCR.

S2 Table [docx]
The combined transcriptomic and proteomic analysis parameters.

S3 Table [docx]
DEGs between the dorsal feather follicles and ventral feather follicles of the neck.

S4 Table [docx]
DEPs between the dorsal feather follicles and ventral feather follicles of the neck.

S5 Table [docx]
Correlation of DEPS and DEGs.


Zdroje

1. Keiji K, Toyoko A, Makoto M, Ai S, Akira I, Hassan YH, et al. Endothelin Receptor B2 (EDNRB2) Is Responsible for the Tyrosinase-Independent Recessive White (mow) and Mottled (mo) Plumage Phenotypes in the Chicken. Plos One. 2014; 9: e86361. doi: 10.1371/journal.pone.0086361 24466053

2. Vickrey AI, Rebecca B, Zev K, Emma M, J BR, T ME, et al. Introgression of regulatory alleles and a missense coding mutation drive plumage pattern diversity in the rock pigeon. Elife. 2018; 7: e34803–. doi: 10.7554/eLife.34803 30014848

3. Cooke TF, Fischer CR, Wu P, Jiang TX, Bustamante CD. Genetic Mapping and Biochemical Basis of Yellow Feather Pigmentation in Budgerigars. Cell. 2017; 171: 427–439. doi: 10.1016/j.cell.2017.08.016 28985565

4. Oro AE, Scott MP. Splitting hairs: dissecting roles of signaling systems in epidermal development. Cell. 1998; 95: 575–578. doi: 10.1016/s0092-8674(00)81624-4 9845357

5. Hardy MH. The secret life of the hair follicle. Trends in Genetics. 1992; 8: 55–61. doi: 10.1016/0168-9525(92)90350-d 1566372

6. Galbraith H, Galbraith H. Fundamental hair follicle biology and fine fibre production in animals. Animal. 2010; 4: 1490–1509. doi: 10.1017/S175173111000025X 22444696

7. Alexeev V, Igoucheva O, Domashenko A, Cotsarelis G, Yoon K. Localized in vivo genotypic and phenotypic correction of the albino mutationin skin by RNA-DNA oligonucleotide. Nature Biotechnology. 2000; 18: 43–47. doi: 10.1038/71901 10625389

8. Demehri S, Kopan R. Notch signaling in bulge stem cells is not required for selection of hair follicle fate. Development. 2009; 136: 891. doi: 10.1242/dev.030700 19211676

9. Englaro W, Rezzonico R, Durand-ClãMent M, Lallemand D, Ortonne JP, Ballotti R. Mitogen-activated protein kinase pathway and AP-1 are activated during cAMP-induced melanogenesis in B-16 melanoma cells. Journal of Biological Chemistry. 1995; 270: 24315–24320. doi: 10.1074/jbc.270.41.24315 7592642

10. Grichnik JM, Burch JA, Burchette J, Shea CR. The SCF/KIT pathway plays a critical role in the control of normal human melanocyte homeostasis. Journal of Dermatological Science. 1998; 16: 233–238.

11. Oka M, Nagai H, Ando H, Fukunaga M, Matsumura M, Araki K, et al. Regulation of melanogenesis through phosphatidylinositol 3-kinase-Akt pathway in human G361 melanoma cells. Journal of Investigative Dermatology. 2000; 115: 699–703. doi: 10.1046/j.1523-1747.2000.00095.x 10998146

12. Sharov AA, Fessing M, Atoyan R, Sharova TY, Haskellluevano C, Weiner L, et al. Bone morphogenetic protein (BMP) signaling controls hair pigmentation by means of cross-talk with the melanocortin receptor-1 pathway. Proceedings of the National Academy of Sciences of the United States of America. 2005; 102: 93–98. doi: 10.1073/pnas.0408455102 15618398

13. Schouwey K, Beermann F. The Notch pathway: hair graying and pigment cell homeostasis. Histology & Histopathology. 2008; 23: 609.

14. Li X, Guo L, Sun Y, Zhou J, Gu Y, Li Y. Baicalein inhibits melanogenesis through activation of the ERK signaling pathway. International Journal of Molecular Medicine. 2010; 25: 923–927. doi: 10.3892/ijmm_00000423 20428797

15. Chiang HM, Chien YC, Wu CH, Kuo YH, Wu WC, Pan YY, et al. Hydroalcoholic extract of Rhodiola rosea L. (Crassulaceae) and its hydrolysate inhibit melanogenesis in B16F0 cells by regulating the CREB/MITF/tyrosinase pathway. Food & Chemical Toxicology. 2014; 65: 129–139.

16. Yeo GSH, Farooqi IS, Challis BG, Jackson RS, O'Rahilly S. The role of melanocortin signalling in the control of body weight: evidence from human and murine genetic models. Qjm Monthly Journal of the Association of Physicians. 2000; 93: 7–14. doi: 10.1093/qjmed/93.1.7 10623776

17. Pape EL, Wakamatsu K, Ito S, Wolber R, Hearing VJ. Regulation of eumelanin / pheomelanin synthesis and visible pigmentation in melanocytes by ligands of the melanocortin 1 receptor. Pigment Cell & Melanoma Research. 2008; 21: 477–486.

18. V?Ge DI, Nieminen M, Anderson DG, R?Ed KH. Two missense mutations in melanocortin 1 receptor (MC1R) are strongly associated with dark ventral coat color in reindeer (Rangifer tarandus). Animal Genetics. 2015; 45: 750–753.

19. Liu F, Kohlmeier S, Wang CY. Wnt signaling and skeletal development. Cellular Signalling. 2008; 20: 999–1009. doi: 10.1016/j.cellsig.2007.11.011 18164181

20. Chen AE, Ginty DD, Fan CM. Protein kinase A signalling via CREB controls myogenesis induced by Wnt proteins. Nature. 2005; 433: 317–322. doi: 10.1038/nature03126 15568017

21. Wagner AJ, Fisher DE. Melanocyte signaling pathways and the etiology of melanoma. Drug Discovery Today Disease Mechanisms. 2006; 2: 179–183.

22. Russell M, Hoeffler JP. Signal transduction and the regulation of cell growth. The journal of investigative dermatology Symposium proceedings / the Society for Investigative Dermatology, Inc [and] European Society for Dermatological Research. 1996; 1: 119.

23. Hubbard JK, Uy JAC, Hauber ME, Hoekstra HE, Safran RJ. Vertebrate pigmentation: from underlying genes to adaptive function. Trends in Genetics. 2010; 26: 231–239. doi: 10.1016/j.tig.2010.02.002 20381892

24. Klungland H, Vage DI. Molecular Genetics of Pigmentation in Domestic Animals. Current Genomics. 2000; 1: -. doi: 10.2174/1389202003351535

25. Chang CM, Coville JL, Coquerelle G, Gourichon D, Oulmouden A, Tixier-Boichard M. Complete association between a retroviral insertion in the tyrosinase gene and the recessive white mutation in chickens. Bmc Genomics. 2006; 7: 19–19. doi: 10.1186/1471-2164-7-19 16457736

26. Sato S, Otake T, Suzuki C, Saburi J, Kobayashi E, Mapping of the recessive white locus and analysis of the tyrosinase gene in chickens. Poult Sci. 2007; 86: 2126–2133. doi: 10.1093/ps/86.10.2126 17878441

27. Yang L, Du X, Wei S, Gu L, Li N, Gong Y, et al. Genome-wide association analysis identifies potential regulatory genes for eumelanin pigmentation in chicken plumage. Animal Genetics. 2017; 48: 611. doi: 10.1111/age.12573 28639704

28. Slominski A, Wortsman J, Plonka PM, Schallreuter KU, Paus R, Tobin DJ. Hair Follicle Pigmentation. Journal of Investigative Dermatology. 2004; 124: 13–21.

29. Dupin E, Glavieux C, Vaigot P, Douarin NML. Endothelin 3 Induces the Reversion of Melanocytes to Glia through a Neural Crest-Derived Glial-Melanocytic Progenitor. Proceedings of the National Academy of Sciences of the United States of America. 2000; 97: 7882. doi: 10.1073/pnas.97.14.7882 10884419

30. Mills MG, Patterson LB. Not just black and white: Pigment pattern development and evolution in vertebrates. Seminars in Cell & Developmental Biology. 2009; 20: 72–81.

31. Li L, Li D, Liu L, Li S, Feng Y, Peng X, et al. Endothelin Receptor B2 (EDNRB2) Gene Is Associated with Spot Plumage Pattern in Domestic Ducks (Anas platyrhynchos). Plos One. 2015; 10: e0125883. doi: 10.1371/journal.pone.0125883 25955279

32. Gunnarsson U, Hellstrom A, Tixier-Boichard M, Minvielle F, Bed'Hom B, Ito S, et al. Mutations in SLC45A2 cause plumage color variation in chicken and Japanese quail. Genetics. 2007; 175: 867. doi: 10.1534/genetics.106.063107 17151254

33. Liu WB, Chen SR, Zheng JX, Qu LJ, Xu GY, Yang N. Developmental phenotypic-genotypic associations of tyrosinase and melanocortin 1 receptor genes with changing profiles in chicken plumage pigmentation. Poultry Science. 2010; 89: 1110. doi: 10.3382/ps.2010-00628 20460655

34. S K, P S, U G, H K, S B, R F, et al. The Dominant white, Dun and Smoky color variants in chicken are associated with insertion/deletion polymorphisms in the PMEL17 gene. Genetics. 2004; 168: 1507–1518. doi: 10.1534/genetics.104.027995 15579702

35. Yeo J, Lee Y, Hyeong KE, Ha J, Yi JK, Kim B, et al. Detection of exonic variants within the melanocortin 1 receptor (MC1R) gene in Black Silky, White Leghorn and Golden duckwing Araucana chicken. Molecular Biology Reports. 2014; 41: 4843–4846. doi: 10.1007/s11033-014-3394-0 24830563

36. Hoekstra HE. Genetics, development and evolution of adaptive pigmentation in vertebrates. Heredity. 2006; 97: 222–234. doi: 10.1038/sj.hdy.6800861 16823403

37. Z R, E B, A B, D W. From tyrosine to melanin: Signaling pathways and factors regulating melanogenesis. Postepy Hig Med Dosw. 2016; 70: 695–708.

38. Sultana H, Seo D, Choi NR, Msa B, Lee SH, Heo KN, et al. Identification of polymorphisms inMITFandDCTgenes and their associations with plumage colors in Asian duck breeds. Asian-Australas J Anim Sci. 2018; 31: 180–188. doi: 10.5713/ajas.17.0298 28823136

39. Lin SJ, Foley J, Jiang TX, Yeh CY, Wu P, Foley A, et al. Topology of feather melanocyte progenitor niche allows complex pigment patterns to emerge. Science. 2013; 340: 1442–1445. doi: 10.1126/science.1230374 23618762

40. Hoekstraweebers JEHM, Hubbard BK, Hauber ME, Uy JAC, Safran RJ. Vertebrate pigmentation: from underlying genes to adaptive function. Trends in Genetics. 2010; 26: 231–239. doi: 10.1016/j.tig.2010.02.002 20381892

41. Liu X, Zhou R, Peng Y, Zhang C, Li L, Lu C, et al. Hair follicles transcriptome profiles in Bashang long-tailed chickens with different plumage colors. Genes & Genomics. 2018:

42. Jing Y, Yanhua Q, Yuan H, Fumin L. Dynamic transcriptome profiling towards understanding the morphogenesis and development of diverse feather in domestic duck. Bmc Genomics. 2018; 19: 391. doi: 10.1186/s12864-018-4778-7 29793441

43. Chen Siang Ng, Chen C-K, Fan W-L, Wu P, Wu S-M. Transcriptomic analyses of regenerating adult feathers in chicken. BMC Genomics. 2015; 16: 756. doi: 10.1186/s12864-015-1966-6 26445093

44. Gong H, Wang H, Wang YX, Bai X, Zhang WG. Skin transcriptome reveals the dynamic changes in the Wnt pathway during integument morphogenesis of chick embryos. Plos One. 2018; 13: e0190933. doi: 10.1371/journal.pone.0190933 29351308

45. Li A, Figueroa S, Jiang T-X, Wu P, Widelitz R, Nie Q, et al. Diverse feather shape evolution enabled by coupling anisotropic signalling modules with self-organizing branching programme. Nature Communications. 2017; 8: ncomms14139. doi: 10.1038/ncomms14139 28106042

46. Zhang L, Nie Q, Su Y, Xie X, Luo W, Jia X, et al. MicroRNA profile analysis on duck feather follicle and skin with high-throughput sequencing technology. Gene. 2013; 519: 77–81. doi: 10.1016/j.gene.2013.01.043 23384715

47. Shen S, Park JW, Lu ZX, Lin L, Henry MD, Wu YN, et al. rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data. Proc Natl Acad Sci U S A. 2014; 111: 5593–5601.

48. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nature Methods. 2015; 12: 357–360. doi: 10.1038/nmeth.3317 25751142

49. S A, JM C. The significance of digital gene expression profiles. Genome Research. 1997; 7: 986–995. doi: 10.1101/gr.7.10.986 9331369

50. Wen B, Zhou R, Feng Q, Wang Q, Wang J, Liu S. IQuant: an automated pipeline for quantitative proteomics based upon isobaric tags. Proteomics. 2015; 14: 2280–2285.

51. Markus B, Lu Y, Tim H, Jyoti C. Accurate and sensitive peptide identification with Mascot Percolator. Journal of Proteome Research. 2009; 8: 3176–3181. doi: 10.1021/pr800982s 19338334

52. Savitski MM, Wilhelm M, Hahne H, Kuster B, Bantscheff M. A Scalable Approach for Protein False Discovery Rate Estimation in Large Proteomic Data Sets. Molecular & Cellular Proteomics Mcp. 2015; 14: 2394.

53. Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, Rao BS, et al. The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Research. 2001; 29: 22–28. doi: 10.1093/nar/29.1.22 11125040

54. Maier T, Güell M, Serrano L. Correlation of mRNA and protein in complex biological samples. Febs Letters. 2009; 583: 3966–3973. doi: 10.1016/j.febslet.2009.10.036 19850042

55. B F, GI L, P M, CS M, JI G. A sampling of the yeast proteome. Molecular and Cellular Biology. 1999; 19: 7357. doi: 10.1128/mcb.19.11.7357 10523624

56. Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Research. 1997; 7: 986–995. doi: 10.1101/gr.7.10.986 9331369

57. Jianqin Z, Fuzhu L, Junting C, Xiaolin L. Skin transcriptome profiles associated with skin color in chickens. Plos One. 2015; 10: e0127301. doi: 10.1371/journal.pone.0127301 26030885

58. Yu S, Wang G, Liao J, Tang M, Sun W. Transcriptome Profile Analysis of Mechanisms of Black and White Plumage Determination in Black-Bone Chicken. Cellular Physiology & Biochemistry. 2018; 46: 2373–2384.

59. Li JM, Huang XS, Li LT, Zheng DM, Xue C, Zhang SL, et al. Proteome analysis of pear reveals key genes associated with fruit development and quality. Planta. 2015; 241: 1363–1379. doi: 10.1007/s00425-015-2263-y 25682102

60. Greenbaum D, Jansen R, Gerstein M. Analysis of mRNA expression and protein abundance data: an approach for the comparison of the enrichment of features in the cellular population of proteins and transcripts. Bioinformatics. 2002; 18: 585–596. doi: 10.1093/bioinformatics/18.4.585 12016056

61. Zhang W, Culley DE, Scholten JCM, Hogan M, Vitiritti L, Brockman FJ. Global transcriptomic analysis of Desulfovibrio vulgaris on different electron donors. Antonie Van Leeuwenhoek. 2006; 89: 221. doi: 10.1007/s10482-005-9024-z 16710634

62. Culley DE. Integrative analysis of transcriptomic and proteomic data: challenges, solutions and applications. Critical Reviews in Biotechnology. 2007; 27: 63. doi: 10.1080/07388550701334212 17578703

63. Crick F. Central Dogma of Molecular Biology. Nature. 1970; 227: 561–563. doi: 10.1038/227561a0 4913914

64. Gygi SP, Rochon Y, Franza BR, Aebersold R. Correlation between Protein and mRNA Abundance in Yeast. Molecular & Cellular Biology. 1999; 19: 1720–1730.

65. Wang SH, You ZY, Ye LP, Che J, Qian Q, Nanjo Y, et al. Quantitative proteomic and transcriptomic analyses of molecular mechanisms associated with low silk production in silkworm Bombyx mori. Journal of Proteome Research. 2014; 13: 735–751. doi: 10.1021/pr4008333 24428189

66. Chen Q, Guo W, Feng L, Ye X, Xie W, Huang X, et al. Data for transcriptome and proteome analysis of Eucalyptus infected with Calonectria pseudoreteaudii. Data in Brief. 2015; 3: 24–28. doi: 10.1016/j.dib.2014.12.008 26217712

67. Trevisan S, Manoli A, Ravazzolo L, Botton A, Pivato M, Masi A, et al. Nitrate sensing by the maize root apex transition zone: a merged transcriptomic and proteomic survey. Journal of Experimental Botany. 2015; 66: 3699–3715. doi: 10.1093/jxb/erv165 25911739

68. Ghazalpour A, Bennett B, Petyuk VA, Orozco L, Hagopian R, Mungrue IN, et al. Comparative Analysis of Proteome and Transcriptome Variation in Mouse. Plos Genetics. 2011; 7: e1001393. doi: 10.1371/journal.pgen.1001393 21695224

69. Wu S, Zhu Y, He F. Progress in the comparison of transcriptome and proteome. Progress in Biochemistry & Biophysics. 2005; 32: 99–105.

70. Vogel C, Marcotte EM. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nature Reviews Genetics. 2012; 13: 227–232. doi: 10.1038/nrg3185 22411467

71. Bouchal P, Dvořáková M, Roumeliotis T, Bortlíček Z, Ihnatová I, Procházková I, et al. Combined Proteomics and Transcriptomics Identifies Carboxypeptidase B1 and Nuclear Factor κB (NF-κB) Associated Proteins as Putative Biomarkers of Metastasis in Low Grade Breast Cancer. European Journal of Cancer. 2015; 51: S24–S25.

72. Petersen HO, Höger SK, Looso M, Lengfeld T, Kuhn A, Warnken U, et al. A Comprehensive Transcriptomic and Proteomic Analysis of Hydra Head Regeneration. Molecular Biology & Evolution. 2015; 32: 1928–1947.

73. Hou L, Panthier JJ, Arnheiter H. Signaling and transcriptional regulation in the neural crest-derived melanocyte lineage: interactions between KIT and MITF. Development. 2000; 127: 5379. 11076759

74. D’Mello SAN, Finlay GJ, Baguley BC, AskarianAmiri ME. Signaling Pathways in Melanogenesis. International Journal of Molecular Sciences. 2016; 17: 1144.

75. Zhang J, Liu F, Cao J, Liu X. Skin Transcriptome Profiles Associated with Skin Color in Chickens. Plos One. 2015; 10: e0127301. doi: 10.1371/journal.pone.0127301 26030885

76. Xia M, Chen K, Yao X, Xu Y, Yao J, Yan J, et al. Mediator MED23 Links Pigmentation and DNA Repair through the Transcription Factor MITF. Cell Reports. 2017; 20: 1794. doi: 10.1016/j.celrep.2017.07.056 28834744

77. Gessi M, Hammes J, Lauriola L, Dörner E, Kirfel J, Kristiansen G, et al. GNA11 and N-RAS mutations: alternatives for MAPK pathway activating GNAQ mutations in primary melanocytic tumours of the central nervous system. Neuropathology & Applied Neurobiology. 2013; 39: 417.

78. Jipeng X, Effects of melamine, oxidized fish oil and lipid on growth and skin colour of darkbarbel catfish. Ph.D. Thesis, Ocean university of China, 2011.


Článok vyšiel v časopise

PLOS One


2019 Číslo 11
Najčítanejšie tento týždeň
Najčítanejšie v tomto čísle
Kurzy

Zvýšte si kvalifikáciu online z pohodlia domova

Aktuální možnosti diagnostiky a léčby litiáz
nový kurz
Autori: MUDr. Tomáš Ürge, PhD.

Všetky kurzy
Prihlásenie
Zabudnuté heslo

Zadajte e-mailovú adresu, s ktorou ste vytvárali účet. Budú Vám na ňu zasielané informácie k nastaveniu nového hesla.

Prihlásenie

Nemáte účet?  Registrujte sa

#ADS_BOTTOM_SCRIPTS#