Integrative analysis of the Pekin duck (Anas anas) MicroRNAome during feather follicle development
© The Author(s). 2017
Received: 27 December 2016
Accepted: 5 July 2017
Published: 20 July 2017
The quality and yield of duck feathers are very important economic traits that might be controlled by miRNA regulation. The aim of the present study was to investigate the mechanism underlying the crosstalk between individual miRNAs and the activity of signaling pathways that control the growth of duck feathers during different periods. We therefore conducted a comprehensive investigation using Solexa sequencing technology on the Pekin duck microRNAome over six stages of feather development at days 11, 15, and 20 of embryonic development (during the hatching period), and at 1 day and 4 and 10 weeks posthatch.
There were a total of 354 known miRNAs and 129 novel candidate miRNAs found based on comparisons with known miRNAs in the Gallus gallus miRBase. The series of miRNAs related to feather follicle formation as summarized in the present study showed two expression patterns, with primary follicle developed during embryonic stage and secondary follicle developed mainly at early post hatch stage. Analysis of miRNA expression profiles identified 18 highly expressed miRNAs, which might be directly responsible for regulation of feather development. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis suggested that in addition to Wnt and transforming growth factor (TGFβ) signaling pathways, which were widely reported in response to follicle formation, another group of signaling pathways that regulate lipid synthesis and metabolism, such as the phosphatidylinositol signaling system and glycerolipid metabolism and signaling, are also responsible for follicle formation.
The highly expressed miRNAs provide a valuable reference for further investigation into the functional miRNAs important for feather development. Lipid synthesis and metabolism related signaling pathways might be responsible for lipid formation on the surface of feather, and should be paid much more attention for their relation to feather quality.
The structure of bird skin and feathers differ from that of similar structures in mammals . Flight feathers and down feathers serve different functions, and the quality and yield of down feathers are very important economic traits . The feathers of ducks are widely used for decoration and to line clothing, bedding, and so on . Fast-growing Pekin ducks usually reach market weight around 35–40 d of age [4, 5], when the plumage is still far from maturity . This inevitably reduces the value of the feathers at this stage. To maximize the economic potential of duck feathers, the diversity that is evident in feather morphology should be used as an advantage to gain deeper insight on the genomic and epigenetic levels. With recent advances in epigenetics, the present is the opportune time to develop our understanding of morphogen signals in feather follicles at the genomic level. Increasing evidence suggests that miRNAs serve as biological regulators that mediate gene expression .
Strands of miRNA exert observable impact on target genes of vertebrates and invertebrates and show hallmarks of incorporation into endogenous regulatory networks. Yuan et al. found 399 conserved miRNAs that played roles in the regulation of hair follicle cycling in cashmere goats (Capra hircus) . The most highly expressed miRNA in that study was quite different from that observed in duck feather regeneration . In addition, miRNA expression exhibits significant differences between embryonic stages and feather regeneration stages in the chicken and duck, respectfully. Bao et al. (2016) reported that 226 miRNA genes differentially expressed among three embryonic stages in the chicken in which 21 miRNAs targeted genes were feather related . However, none of the differentially expressed miRNAs were found during feather regeneration in duck .
In our previous study, genes that regulated feather regeneration differed significantly from those associated with similar functions in mammals . Thus, we assumed that feather development in avian species is regulated by miRNA, which might differ from that in mammalian skin and hair. Therefore, we selected skin samples of ducks from the early embryonic stages up until feathers were fully developed at 10 weeks of age. In addition, we investigated the miRNAome during feather development using Solexa sequencing. A better understanding of the complex mechanisms underlying feather development in ducks might be of benefit to the down feather industry and genes participating in the signaling pathways that regulate feather follicle formation might be used as candidate genes for molecular marker screening in the future.
All animal procedures were performed in accordance with guidelines developed by the China Council on Animal, Care and protocols were approved by the Animal Care and Use Committee of Anhui Agricultural University, China.
From our previous study on follicle development in broilers and geese, we deduced that follicle formation in the duck starts around day 11 of embryonic development (11EM), progresses so that the body becomes lightly covered with feathers by day 15 (15EM), following which feathers completely cover the body at day 20 (20EM) . A significant increase in secondary follicle growth has been observed post hatch . After 4–5 weeks of age, changes have been observed in duckling feathers, which become fully mature feathers at 10 weeks of age .
White Pekin duck embryos (Anhui Taiyang Poultry Co. Ltd., China) of similar genetic background were selected at days 11, 15, and 20 of embryonic development (11EM, 15EM, and 20EM). Ducks were reared under normal conditions of light and temperature and had free access to food and water at Anhui Taiyang Poultry Co. Ltd., China. From 1 day to 10 weeks of age, ducks were fed with commercial diets that met all National Research Council requirements (NRC, 1994). Ducks were slaughtered by sectioning the neck 2–4 h after the last meal at three stages of development: 1 d, 4 and 10 weeks of age (1 DB, 4WK, and 10WK) post-hatch. For each stage, six skin samples collected from an area just above the second to third thoracic vertebrae were used for total RNA extraction. Each of these samples was immediately stored in a tube with 1 mL RNALater (Qiagen) and stored at −20 °C. Another six skin samples from the same area were stored immediately in vials containing 10% buffered formalin for 24 h and transferred to 70% ethanol for storage until embedding and processing. Tissue sections were laterally cut for Masson trichrome staining . The analysis of the slides was performed according to Chen et al. .
Total RNA was extracted from the skin samples using a mirVanamiRNA Isolation Kit (Ambion, USA) according to Qin et al. (2013) . The purified RNA yield was determined by the absorbance at 260 nm with an ND-2000 NanoDrop spectrophotometer (Thermo Fisher, USA), and RNA quality was evaluated with the BioAnalyzer 2100 system (Agilent Technologies, USA) that was done at Shanghai Biotechnology Corporation.
Small RNA library preparation and deep sequencing
RNA smaller than 200 bp was enriched with the mirVanamiRNA isolation kit (Ambion). RNA was then precipitated with ethanol and dissolved in water. Small RNAs had linkers ligated to them and bar-coded cDNAs were prepared using a TruSeq Small RNA Sample Prep Kit (Illumina, USA) following the manufacturer’s instructions. Small RNA (1 μg) was ligated with adapters at the 3′ and 5′ ends. Adapter-ligated RNA was reverse-transcribed with SuperScript II Reverse Transcriptase (Invitrogen, USA), and then PCR-amplified (11 cycles). Individual libraries were analyzed on a BioAnalyzer (Agilent) for the presence of linked cDNA at the appropriate size (140–150 bp), as determined by the BioAnalyzer. The amplified cDNA constructs were then purified from agarose gel in preparation for sequencing analysis, which was conducted at the Shanghai Biotechnology Corporation, using the Illumina HiSeq 2500 System (Illumina, CA, USA) according to the manufacturer’s instructions.
Bioinformatic analyses of sequencing data
The small RNA sequence reads were pre-processed using the FASTX-Toolkit , excluding low-quality reads (ambiguous N, quality <10 nucleotides [nt], and length < 18 nt) as well as the 3′ and 5′ adapters and poly (A) sequences. Further annotation analyses were performed using the commercial software CLC Genomic Workbench 5.5. The resulting clean reads were aligned against various databases, including ncRNA, piRNA, and Rfam, allowing a maximum mismatch of 2 nt to remove noncoding RNA, such as rRNA, tRNA, snRNA, and snoRNA. The remaining sequences were analyzed by a BLAST Gallus gallus search against the Sanger miRBase (version 21.0). Reads that did not match any of the databases mentioned above were marked as non-annotated. Non-annotated sequences were searched against the Gallus gallus genome using the miRCat program included in the sRNAToolkit (http://srna-workbench.cmp.uea.ac.uk/tools/analysis-tools/mircat/). Using default settings, 100 nt flanking each side of the genomic sequences were extracted for prediction of RNA secondary structure using RNAfold . Only typical stem-loop hairpin structures with free energy lower than −20 kcal/mol were considered potential novel miRNAs. After the completion of all annotation steps, sequencing libraries were used for size distribution and saturation analysis. All sequence data have been submitted to the NCBI Sequence Read Archive (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE101542) under accession No. SRA073195.
Identification of miRNA differential expression
We first normalized miRNA-sequence data from 18 libraries as transcripts per million (TPM). The followingnormalization formula was used:
Normalized expression = Actual miRNA count/Total count of clean reads × 1000,000. If the normalized expression of a given miRNA was zero, its expression value was modified to 0.01. Normalized sequence counts were used to perform a one-way ANOVA to determine significant differences. The expression ofa specific miRNA was considered significantly different if the p-value and False Discovery Rate (FDR) value were both less than 0.05. FDR indicates the expected proportion of false positives among the results deemed significant .
Time series analysis
The Short Time-series Expression Miner (STEMv 1.1) program was used to cluster and visualize possible profiles and changes in expression over time (less than 8 time points) in differentially expressed (DE) miRNAs. The maximum unit change in model profiles between time points was adjusted to 1 and the maximum number of model profiles to 50. MicroRNA expression profiles were clustered according to correlation coefficients. Statistical significance of the number of genes assigned to each profile versus the numbers expected were computed by an algorithm and the default p-value was 1e-5 . Statistically significant model profiles that were similar to each other were grouped together to form clusters of profiles.
Prediction and analysis of miRNA target genes
As there is no appropriate database or method to predict duck miRNA target genes, miRanda version 3.1 (http://www.microrna.org/microrna/getMirnaForm) was used for target gene prediction. To completely assess functions of the differentially expressed miRNAs, we conducted gene ontology (GO) and KEGG pathway analyses enriched with predicted miRNA targets. The analyses were conducted using the fuzzy clustering algorithm” to reduce redundancy among functionally related pathways that share similar target genes. Terms with Benjamini-corrected enrichment p < 0.01 and FDR < 0.05 were considered. An association map was generated that summarizes the enriched pathways in a graphical representation of the relationships between terms based on the similarity of their target genes. All of the data analysis methods used in this experiment were the same as mentioned in Qin et al.  because we both performed the miRNA sequencing at Shanghai Biotechnology Corporation.
MiRNA validation via stem-loop RT-PCR
Small RNAs (< 200 nt) were isolated using the mirVanamiRNA Isolation Kit (Ambion, USA) following the manufacturer’s instructions. Quantitative real-time PCR (qRT-PCR) was carried out as previously described in an ABI PRISM 7500 Fast Realtime PCR System (Ambion, USA), using the SYBR Premix Ex Taq™ Kit (TaKaRa, Japan) . The reactions were carried out in a volume of 20 μL, containing 2 μL of diluted cDNA, 200 nM of each primer, and 16 μLPCR Master Mix under the following conditions: 95 °C for 30 s; 45 cycles at 95 °C for 5 s; 58 °C for 15 s; and 72 °C for 20 s. A thermal denaturing cycle then followed at 95 °C for 15 s and 60 °C for 1 min applied to determine the dissociation curves, which were used to verify the specificity of PCR amplifications. All reactions were run in triplicate for each sample. Ten miRNAs, including seven conserved and three novel miRNAs were validated, and 5.8S rRNA was selected as a reference gene for normalization (Additional file 1: Table S1). The experimental data was analyzed using the 2-∆∆CTmethod.
Microscopic observation of follicle characters
Overview of Solexa-sequencing of feather follicle small RNAs
Identification of potential novel miRNAs
To detect additional potential miRNAs, unannotated sequences within 18 to 35 nt were searched against chicken miRNA database and analyzed using the miRDeep software. Sequenced reads that did not match any known ncRNA, piRNA, and miRNAs from the Rfam database, and had typical stem-loop hairpin structures with free energy lower than −20 Kcal/mol were considered potential novel miRNA. There were totally 129 candidate novel miRNAs were predicted based on novel miRNA identification criterion and then named as gga-miR-new-N (1–129) (Additional file 4: Table S3). The most frequent sequence length was 18 nt (27.9%), followed by those that measured 22 nt (17.1%), 19 and 20 nt (15.5% in each), and 23 nt (9.3%).
Expression analysis of miRNAs and miRNA families in duck feather development
A total of 354 miRNAs were identified in duck feather follicles (Additional file 5: Table S4). All miRNAs were expressed at varying levels, ranging from 1 to more than 100,000,000 reads. Mullokandov et al. (2012) developed a Sensor-seq to show that over 60% of detected miRNAs had no discernible activity and only the most highly expressed miRNAs in a cell mediate target suppression . From the normalized miRNAs reads, there were 296 miRNAs with mean expression values below 1000 Transcript Per Million (TPM) (243 below 100 and 53 between 100 to 1000); these miRNAs account for almost 83.6% of the total number of miRNAs (68.6% and 15.0%, respectively) within the 18 libraries (Additional file 5: Table S4). The mean number of miRNAs with expression values over 1000 TPM was 58 (16.4%), among which 19 miRNAs (5.4%) expressed over 10,000 TPM, and two miRNAs expressed over 100,000 TPM, and the 19 top most expressed miRNAs accounted for 82.5% of the total unique miRNAs (Additional file 6: Figure S2 and Additional file 7: Table S5).
From the 129 novel miRNAs detected in 18 libraries, the total reads were at a relatively low level. There were 100 novel miRNAs expressed less than 100 times, and 17 were at 100–1000 TPM (of which four were sequenced more than 500 times), 12 were at 1000–10,000 times, and only gga-mir-new-12 was sequenced more than 10,000 times (Additional file 4: Table S3).
Identification of differentially expressed miRNAs at different developmental stages
Numbers of differentially expressed miRNAs between libraries
Numbers of differentially expressed miRNAs
Comparison between adjacent libraries
Total DE miRNAs
Total DE miRNAs-up
Total DE miRNAs-down
MicroRNA target predictions and KEGG Orthology analysis
Validation of sequencing data
Time points selection for miRNA libraries construction
Some miRNAs are involved in the development of hair follicles and skin, and are related to hair follicle cycles in mammals [8, 24, 25]. However, there have been few studies on poultry miRNAs and their relation to the development of feather growth. Avian skin and feathers are very different from mammalian skin and hair, because of the differences in function (despite similarities in structure) between flight and down feathers. Zhang et al. detected 96 miRNAs in duck skin during feather regeneration at 8 weeks of age . However, feather growth might already have been initiated during the early embryonic stages. Thus, we chose the six time points to elucidate the miRNA expression profiles of duck feather formation in this experiment.
MicroRNA database is reliable for analysis of miRNA expression profiles during duck feather growth
We identified 354 known gga-miRNAs and 129 potential novel miRNAs between the embryonic stages and 10 weeks of age. In viewing the sequenced reads, the libraries capacity approached saturation demonstrated that the deep sequencing data were able to represent the miRNA transcriptome profiles of duck feather growth. Among the annotated small RNAs, only one third were matched to chicken miRNA, suggesting that the miRNAs including very richful diversities, and the aforementioned results promote confidence in the deep sequencing data and the present data are reliable for analysis of miRNA expression profiles.
The highly expressed miRNA were mostly related to regulation of cell cycle
Skin development is governed by complex processes of gene activation and silencing, and entails miRNA-dependent modulation of gene expression. Many miRNAs that regulate hair follicle formation are also involved in the regulation and development of cancer and skin diseases in mice and humans, mainly because they regulate cell cycles. MiR-214 regulates skin morphogenesis and hair follicle cycling by targeting β-catenin [26, 27]. It can also be expressed in several human tumors, such as those of ovarian cancer and breast cancer, by targeting Bim to promote cell proliferation. However, few studies have focused on feather formation in birds, let alone in waterfowl.
There were only 19 miRNAs expressed over 10,000 TPM suggested that the majority of abundantly expressed miRNAs were fewer in number. We also found that more than half of the 19 miRNAs were related to regulation of the cell cycle. The most abundant miRNAs were gga-miR-10a and gga-miR-10b, which represented more than 2,000,000 TPM within 18 libraries. MiR-10a reportedly plays an important role in regulation of squamous cell growth in the neck area . The miRNAs miR-10a and miR-10b are close homologs, differing only by a single central nucleotide . In the mouse embryo, miR-10a is mainly expressed in a region of the posterior trunk . The miRNAs miR-10a and miR-10b are reportedly involved in cell proliferation and regulation of the cell cycle . Two other miRNAs, gga-miR-26a and gga-miR-181a, were highly expressed and ranked as the third and fourth most highly expressed miRNAs, respectively. They represented more than 1000,000 TPM in 18 libraries and have been identified to play a role in tissue regeneration via the BMP/SMAD1 signaling pathway [32, 33]. Frucht et al. reported that miR-181a plays a key role in regeneration of basilar papillae in the cochlea of the chicken . Leeper et al. also demonstrated that miR-26a could inhibit cellular differentiation and apoptosis by altering the TGFβ signaling pathway . It has been widely accepted that BMP and TGFβ signaling pathways play pivotal roles in feather follicle growth and development. MiR-205, the seventh most highly expressed miRNA, targets the lipid phosphatase SHIP2 in epithelial cells, whereas the corneal-specific miR-184 interferes with the ability of miR-205 to suppress SHIP2 levels . MiR-21, the 11th most highly expressed miRNA, negatively regulates BMP4 and has been reported in the epidermis and hair follicle epithelium of normal mouse skin . Other miRNAs, which were responsible for feather growth andregulation of gene expression, showed significantly lower levels of expression. Expression of miR-24 alters the normal process of hair keratinocyte differentiation, leading to altered expression of differentiation markers, by direct repression of the hair keratinocyte stemness regulator TCF-3, which sustains postnatal epidermal homeostasis . MiR-31 is responsible for the anagen phase of the follicle cycle and could alter hair shaft formation anddifferentiation in hair matrix keratinocytes . MiR-31 negatively regulates expression of Fgf10, components of the Wnt and BMP signaling pathways, including sclerostin, BAMBI, and Dlx3 transcription factor, as well as selected keratin genes. MiR-203 is a skin- and keratinocyte-specific miRNA that has been implicated in repression of ‘stemness’ in epidermal progenitors by targeting p63 and the suppressor of cytokine signaling-3 (SOCS-3) . Up-regulation of miR-203 is required for differentiation in human keratinocytes and is dependent on the activation of the Pkc/AP-1 pathway . Thus, it is clear that the 19 most highly expressed miRNAs might regulate feather generation via BMP and TGFβ pathways.
Comparing the 19 most highly expressed miRNAs with the total miRNA expression during the period of feather regeneration in ducks , only miR-30a, miR-126, and miR-181a were also highly expressed during feather regeneration in ducks . A single specific stage of feather growth cannot clearly illustrate miRNA function, because its expression pattern is both tissue- and developmental stage-specific. However, we propose that these three miRNAs are mainly responsible for feather growth, whereas other miRNAs might also share responsibility for feather and follicle growth.
Differentially expressed miRNA might suggest the different feather developmental stage
Based on the differentially expressed miRNAs at different developmental stages, and also combined with previous studies on feather growth , we speculated that the peak time for primary feather follicle formation might be during the incubation period, while the peak for second feather follicle formation was during the early growth period post hatch and before 4 weeks of age. However, we could not conclude the concise growth period for the secondary feather follicle growth since we did not construct the miRNA database from one day post hatch till 4 weeks of age. Only two genes detected in this experiment played a role in feather formation after 4 weeks of age. However, the feathers of ducklings are normally replaced around 4 weeks of age, suggesting no obvious changes in structure and morphology in the follicles of white Pekin duck after this period.
In birds, skin and feathers serve different functions, as the flight and down feathers are adapted for different roles. Follicle formation is initiated during the early embryonic stage . From the most highly expressed miRNAs, we found miR-10b and miR-205a were also found to be differentially expressed at the eighth embryonic day in chicken, which targeted genes ALDH1A3 and GALGA_5AKI_KRT19, respectively . However, feather growth differs from that of follicular formation . Of the most highly expressed miRNAs during follicle formation, we found miR-181a and miR-30a were also the most highly expressed during feather regeneration in duck . The involvement of mmu-miR-199a and miR-140 in hair follicle and skin development are related to hair follicle cycling in mammals [24, 43]. These two miRNAs were down-regulated from 11EM to 1 DB, suggesting their role in the regulation of feather generation during the early embryonic periods. Only two DE miRNAs were detected from 4WK to 10WK: miR2131–3p and miR-196-1-3p, which reportedly targets the 3′ untranslated region (UTR) of the rarab gene, and is responsible for pectoral fin bud initiation in zebrafish .
The KEGG signaling pathway analysis
Signal transduction pathways and molecules that control hair and follicle formation in mammals have been well understood over the past few decades. However, the signals and molecules that control avian feather growth and regeneration might be quite different because of the complexity of ectodermal organs with hierarchical branching patterns, that regulate endothermy, communication, and flight .
In order to evaluate the differentially expressed miRNAs, the KEGG signaling pathways were further analyzed. The major pathways predicted in the present study have been previously reported. Feather formation signal pathways comprised mainly the regulation of the actin cytoskeleton, phosphatidylinositol signaling system, insulin signaling pathway, focal adhesion, endocytosis, ubiquitin-mediated proteolysis, MAPK, ErbB, mTOR, and glycerolipid metabolism, Wnt, VEGF, Notch, the cell cycle, and TGF-β, among others. Most of these signal pathways have already implicated in other tissue regeneration processes, such as skin and hair follicle development in the mouse, human, and cashmere goats [8, 45–48]. However, it is worth mentioning, that with the exception of the signal pathways that are usually reported, there were other types of signal pathways that regulated lipid synthesis and metabolism, such as the phosphatidylinositol signaling system, glycerolipid metabolism, O-glycan biosynthesis, inositol phosphate metabolism, glycerophospholipid metabolism, and biosynthesis of unsaturated fatty acids. A layer of fat and oil typically coats the surface of feathers in waterfowl . A certain amount of fat and oil is necessary to maintain the water repellent properties of such feathers. Thus, we propose that genes that were involved in lipid synthesis and metabolism signal pathways could be candidate genes that are also responsible for feather quality.
We first sequenced the miRNAome in duck skin from its embryonic day till market size. The highly expressed miRNAs founded provide a valuable reference for further investigation into the functional miRNAs important for feather development. Feather follicle development could be divided into two distinct stages, one for primary follicle and the other for secondary follicle, which initiated different miRNAs expression. Lipid synthesis and metabolism related signaling pathways might responsible for feather quality.
This study was financially supported by the Natural Science Foundation of Anhui Province (1508085MC53),The Major Science and Technology project of Anhui Province (16030701067), key project in the National Science & Technology Pillar Program during the Twelfth Five-year plan period (2015BAD03B06). We thank the Shanghai Biotechnology Corporation for providing us with technical assistance in Solexa sequencing and bioinformatics analysis. We thank the Writing Center from Mississippi State University for grammar revision.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files].
XC, KG, MW, and CZ completed all experiments including sample collection, qRT-PCR verification, miRNA extraction and sequencing, and tissue staining, and participated in writing the manuscript and critically revised the manuscript. XC, MW, and ZG completed all bioinformatics analysis. ZG provided the experimental design and developed the methodology. ZG provided the financial supports for all sample sequencing and tissue histology, XC provided the financial support for gene quantification. All authors read and approved the final manuscript and agreed to be accountable for all aspects of the work.
Ethics approval and consent to participate
All animal procedures were performed in accordance with guidelines developed by the China Council on Animal, Care and protocols were approved by the Animal Care and Use Committee of Anhui Agricultural University, China.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Chen CF, Foley J, Tang PC, Li A, Jiang TX, Wu P, Widelitz RB, Chuong CM. Development, regeneration, and evolution of feathers. Annu Rev Anim Biosci. 2015;3:169–95.View ArticlePubMedGoogle Scholar
- Kozák J. An overview of feathers formation, moults and down production in geese. Asian-Aust J Anim Sci. 2011;24:881–7.View ArticleGoogle Scholar
- Liu Y, Chen X, Xin JH. Hydrophobic duck feathers and their simulation on textile substrates for water repellent treatment. Bioinspir Biomim. 2008;3:046007.View ArticlePubMedGoogle Scholar
- Fan HP, Xie M, Wang WW, Hou SS, Huang W. Effects of dietary energy on growth performance and carcass quality of white growing Pekin ducks from two to six weeks of age. Poult Sci. 2008;87:1162–4.View ArticlePubMedGoogle Scholar
- Jones TA, Dawkins MS. Environment and management factors affecting Pekin duck production and welfare on commercial farms in the UK. Br Poult Sci. 2010;51:12–21.View ArticlePubMedGoogle Scholar
- Johnson OW. Reproductive cycle of the mallard duck. Condor. 1961;63:351–64.View ArticleGoogle Scholar
- Shivdasani RA. MicroRNAs: regulators of gene expression and cell differentiation. Blood. 2006;108:3646–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Yuan C, Wang X, Geng R, He X, Qu L, Chen Y. Discovery of cashmere goat (Capra hircus) microRNAs in skin and hair follicles by Solexa sequencing. BMC Genomics. 2013;14:511.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang L, Xie X, Jia S, Xiao M, Lin S, An L, Luo W, Jia X, Nie Q, Zhang X. Characterization of microRNA species in Peking duck skin. J Int Agricul. 2013a;12:1614–9.View ArticleGoogle Scholar
- Bao W, Greenwold MJ, Sawyer RH. Expressed miRNAs target feather related mRNAs involved in cell signaling, cell adhesion and structure during chicken epidermal development. Gene. 2016;591:393–402.View ArticlePubMedGoogle Scholar
- Chen X, Xie S, Zhou L, Jiang R, Geng Z. Idnetification of differentially expressed genes in skin of Wanxi-white goose during regeneration of downy feather. Acta Veterinaria et Zootechnica Sinica. 2013;44:1030–6.Google Scholar
- Chen XY, Bai HC, Li L, Jia Ng RS, Zhang W, Geng ZY. Follicle characteristics and follicle developmental related Wnt6 polymorphism in Chinese indigenous Wanxi-white goose. Mol Biol Rep. 2012;39:9843–8.View ArticlePubMedGoogle Scholar
- Howell SNG, Corben C, Pyle P, Rogers DI. The first basic problem: a review of molt and plumage homologies. Condor. 2003;105:635–53.View ArticleGoogle Scholar
- Qin L, Chen Y, Liu X, Ye S, Yu K, Huang Z, Yu J, Zhou X, Chen H, Mo D. Integrative analysis of porcine microRNAome during skeletal muscle development. PLoS One. 2013;8:e72418.View ArticlePubMedPubMed CentralGoogle Scholar
- Roberts SB. QPX Genome Browser Feature Tracks; 2013. Database: figshare [Internet]. Accessed :http://figshare.com/articles/QPX_Genome_Browser_Feature_Tracks/701214.
- Hofacker IL. Vienna RNA secondary structure server. Nucleic Acids Res 2003; 31:3429-3431.Google Scholar
- Zhu E, Zhao FQ, Xu G, Hou HB, Zhou LL, Li XL, Sun Z, Wu J. mirTools: microRNA profiling and discovery based on high-throughput sequencing. Nucleic Acids Res. 2010;38:392–7.View ArticleGoogle Scholar
- Ernst J, Nau GJ, Bar-Joseph Z. Clustering short time series gene expression data. Bioinformatics. 2005;21:159–68.View ArticleGoogle Scholar
- Varkonyi-Gasic E, Wu R, Wood M, Walton EF, Hellens RP. Protocol: a highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods. 2007;3:12.View ArticlePubMedPubMed CentralGoogle Scholar
- Mullokandov G, Baccarini A, Ruzo A, Jayaprakash AD, Tung N, Israelow B, Evans MJ, Sachidanandam R, Brown BD. High-throughput assessment of microRNA activity and function using microRNA sensor and decoy libraries. Nat Methods. 2012;9:840–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Chu Q, Cai L, Fu Y, Chen X, Yan Z, Lin X, Zhou G, Han H, Widelitz RB, Chuong C, Wu W, Yue Z. Dkk2/Frzb in the dermal papillae regulates feather regeneration. Dev Biol. 2014;387:167–78.View ArticlePubMedPubMed CentralGoogle Scholar
- Widelitz RB, Jiang TX, Lu J, Chuong CM. β-catenin in epithelial morphogenesis: Coversion of part of avian foot scales into feather buds with a mutatedβ-catenin. Dev Biol. 2000;219:98–114.View ArticlePubMedGoogle Scholar
- Li A, Chen M, Jiang TX, Wu P, Nie Q, Widelitz R, Chuong CM. Shaping organs by a wingless-int/notch/nonmuscle myosin module which orients feather bud elongation. Proc Nat Acad Sci U S A. 2013;110:E1452–61.View ArticleGoogle Scholar
- Andl T, Murchison EP, Liu F, Zhang Y, Yunta-Gonzalez M, Tobias JW, Andl CD, Seykora JT, Hannon GJ, Millar SE. The miRNA-processing enzyme dicer is essential for the morphogenesis and maintenance of hair follicles. Curr Biol. 2006;16:1041–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Mardaryev AN, Ahmed MI, Vlahov NV, Fessing MY, Gill JH, Sharov AA, Botchkareva NV. Micro-RNA-31 controls hair cycle-associated changes in gene expression programs of the skin and hair follicle. FASEB J. 2010;24:3869–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Ahmed MI, Alam M, Emelianov VU, Poterlowicz K, Patel A, Sharov AA, Mardaryev AN, Botchkareva NV. MicroRNA-214 controls skin and hair follicle development by modulating the activity of the Wnt pathway. J Cell Biol. 2014;207:549–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Z, Li Y, Wang H, Fu S, Wang X, Zeng M, Zeng Y, Shao J. Knockdown of miR-214 promotes apoptosis and inhibits cell proliferation in nasopharyngeal carcinoma. PLoS One. 2014;9:e86149.View ArticlePubMedPubMed CentralGoogle Scholar
- Inoue N, Isomoto H, Matsushima K, Hayashi T, Kunizaki M, Hidaka S, Machida H, Mitsutake N, Nanashima A, Takeshima F, Nakayama T, Ohtsuru A, Nakashima M, Nagayasu T, Yamashita S, Nakao K, Kohno S. Down-regulation of microRNA10a expression in esophageal squamous cell carcinoma cells. Oncol Lett. 2010;1:527–31.PubMedPubMed CentralGoogle Scholar
- Ørom UA, Nielsen FC, Lund AH. MicroRNA-10a binds the 5’UTR of ribosomal protein mRNAs and enhances their translation. Mol Cell. 2008;30:460–71.View ArticlePubMedGoogle Scholar
- Jia H, Zhang Z, Zou D, Wang B, Yan Y, Luo M, Dong L, Yin H, Gong B, Li Z, Wang F, Song W, Liu C, Ma Y, Zhang J, Zhao H, Li J, Yu J. MicroRNA-10a is down-regulated by DNA methylation and functions as a tumor suppressor in gastric cancer cells. PLoS One. 2014;9:e88057.View ArticlePubMedPubMed CentralGoogle Scholar
- Dong CG, Wu WKK, Feng SY, Wang XJ, Shao JF, Qiao J. Co-inhibition of microRNA-10b and microRNA-21 exerts synergistic inhibition on the proliferation and invasion of human glioma cells. Int J Oncol. 2012;41:1005–12.PubMedGoogle Scholar
- Frucht CS, Santos-Sacchi J, Navaratnam DS. MicroRNA181a plays a key role in hair cell regeneration in the avian auditory epithelium. Neurosci Lett. 2011;493:44–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Icli B, Wara AKM, Moslehi J, Sun X, Plovie E, Cahill M, Marchini JF, Schissler A, Pader RF, Shi J, Cheng H, Raghuram S, Arany Z, Liao R, Croce K, MacRae C, Feinberg MW. MicroRNA-26a regulates pathological and physiological angiogenesis b targeting BMP/SMAD1 signaling. Circ Res. 2013;113:1231–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Leeper NJ, Raiesdana A, Kojima Y, Chun HJ, Azuma J, Maegdefessel L, Kundu RK, Quertermous T, Tsao PS, Spin JM. MicroRNA-26a is a novel regulator of vascular smooth muscle cell function. J Cell Physiol. 2011;226:1035–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu J, Ryan DG, Getsios S, Oliveira-Fernandes M, Fatima A, Lavker RM. MicroRNA-184 antagonizes microRNA-203 to maintain SHIP2 levels in epithelia. Proc Nat Acad Sci U S A. 2008;105:19300–5.View ArticleGoogle Scholar
- Ahmed MI, Mardaryev AN, Lewis CJ, Sharov AA, Botchkareva NV. MicroRNA-21 is an important downstream component of BMP signaling in epidermal keratinocytes. J Cell Sci. 2011;124:3399–404.View ArticlePubMedPubMed CentralGoogle Scholar
- Amelio I, Lena AM, Bonanno E, Melino G, Candi E. miR-24 affects hair follicle morphogenesis targeting Tcf-3. Cell Death Dis. 2013;4:e922.View ArticlePubMedPubMed CentralGoogle Scholar
- Wei T, Orfanidis K, Xu N, Janson P, Ståhle M, Pivarcsi A, Sonkoly E. The expression of microRNA-203 during human skin morphogenesis. Exp Dermatol. 2010;19:854–6.View ArticlePubMedGoogle Scholar
- Sonkoly E, Wei T, Loriè EP, Suzuki H, Kato M, Törmä H, Stahle M, Pivarcsi A. Protein kinase C-dependent upregulation of miR-203 induces the differentiation of human keratinocytes. J Invest Dermatol. 2010;130:124–34.View ArticlePubMedGoogle Scholar
- Zhang L, Nie Q, Su Y, Xie X, Luo W, Jia X, Zhang X. MicroRNA profile analysis on duck feather follicle and skin with high-throughput sequencing technology. Gene. 2013;519:77–81.View ArticlePubMedGoogle Scholar
- Eames BF, Schneider RA. Quail-duck chimeras reveal spatiotemporal plasticity in molecular and histogenic programs of cranial feather development. Development. 2005;132:1499–509.View ArticlePubMedPubMed CentralGoogle Scholar
- Alibardi L. Ultrastructure of the feather follicle in relation to the formation of the rachis in pennaceous feather. Anat Sci Int. 2010;85:79–91.View ArticlePubMedGoogle Scholar
- Guo L, Huang ZX, Chen XW, Deng QK, Yan W, Zhou MJ, Ou C, Ding Z. Differential expression profiles of microRNAs in NIH3T3 cells in response to UVB irradiation. Photochem Photobiol. 2009;85:765–73.View ArticlePubMedGoogle Scholar
- He X, Yan YL, Eberhart JK, Herpin A, Wagner TU, Schartl M, Postlethwait JH. miR-196 regulates axial patterning and pectoral appendage initiation. Dev Biol. 2011;357:463–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Teichert A, Elallieh H, Bikle D. Disruption of the hedgehog signaling pathway contributes to the hair follicle cycling deficiency in Vdr knockout mice. J Cell Physiol. 2010;225:482–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lei MX, Chuong CM, Widelitz RB. Tuning Wnt signals for more or fewer hairs. J Invest Dermatol. 2013;133:7–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Rishikaysh P, Dev K, Diaz D, Qureshi WMS, Filip S, Mokry J. Signaling involved in hair follicle morphogenesis and development. International J Mol Sci. 2014;15:1647–70.View ArticleGoogle Scholar
- Xiong Y, Liu Y, Song Z, Hao F, Yang X. Identification of Wnt/β-catenin signaling pathway in dermal papilla cells of human scalp hair follicles: TCF4 regulates the proliferation and secretory activity of dermal papilla cell. The J Dermatol. 2014;41:84–91.View ArticlePubMedGoogle Scholar
- Abdoli MA, Mohamadi F, Ghobadian B, Fayyazi E. Effective parameters on biodiesel production from feather fat oil as a cost-effective feedstock. Int J Environ Res. 2014;8:139–48.Google Scholar