- Research article
Study on LOC426217 as a candidate gene for beak deformity in chicken
BMC Geneticsvolume 17, Article number: 44 (2016)
The beak deformity (crossed beaks) was found in some indigenous chickens of China, such as Beijing-You (BJY), Qingyuan Partridge, and Huxu Chickens. Birds with deformed beaks have reduced feed intake and drinking, impeded growth rate, and poor production performance. Beak deformity reduces the economy of poultry industry and affects animal welfare as well. The genetic basis of this malformation remains incompletely understood. LOC426217, also named claw keratin-like, was the most up-regulated gene in the deformed beaks from a previous digital gene expression (DGE) analysis and was selected as an important candidate gene for further analysis.
In the present study, quantitative real-time PCR (qRT-PCR) was firstly performed to determine the expression pattern of LOC426217 gene in deformed and normal beaks to verify the DGE results. Tissue-specific expression profile of this gene in 14 tissues was also determined using qRT-PCR. The LOC426217 was amplified from the genomic DNA of 171 deformed and 164 normal beaks, and sequenced to detect the single nucleotide polymorphisms (SNPs). The results showed that LOC426217 was significantly high-expressed in the deformed beaks, which was in good agreement with the DGE results. This gene was specifically high-expressed in beaks than other tissues. Eight SNPs were detected in LOC426217: -62G > T, 24 T > C, 36G > C, 192A > T, 204C > T, 222 T > C, 285G > T, and 363 T > C. Genotype frequency of G-62 T, T24C, G36C, T222C, and T363C loci was significant different between deformed and normal beaks. Haplotype analysis revealed one block with SNPs T24C and G36C, and one block with SNPs A192T, C204T, T222C, and G285T in normal birds, while the block with SNPs G36C and A192T in deformed ones.
It was concluded from these results that the over-expression of LOC426217 in the beak maybe related to the malformation. The polymorphisms of LOC426217 gene were associated with the beak deformity trait where the SNPs of G-62 T, T24C, G36C, T222C, and T363C loci maybe used as markers. The specific haplotype block in deformed birds may be a potential linkage marker for this trait.
The beak is an external structure of birds, consisting of the upper and lower mandibles covered with a thin keratinized layer of epidermis . It is used for many important activities such as feeding, drinking, fighting, and preening. In addition to striking morphological differences between species, beak deformities of different forms (noticeably elongated, crossed, bent at right angles) have been documented in many wild birds [2–7]. Frequencies of 1 % to 3 % of beak deformity (normally a crossed beak) were found in various indigenous chickens of China, such as Beijing-You (BJY) (studied here), Silkies, Qingyuan Partridge, and Huxu Chickens. Chickens with deformed beaks have reduced feed intake and growth rate. Therefore, beak deformity represents an economic as well as an animal welfare problem in poultry industry. According to our observations in a BJY population, in the absence of known environmental factors contributing to the malformation, birds with deformed beaks present consistently in each generation and cannot be eliminated from a population simply on the basis of the phenotype. This indicated the genetic effects underlying this trait. Studies have been performed to identify the teratogenic genes or molecular genetic background of beak deformity. Previously recognized genetic factors associated with beak deformity include some knwon genes such as fibroblast growth factor 8 (FGF8) , bone morphogenetic protein 4 (BMP4) [9–11], calmodulin (CaM) , and ALX homeobox 1 (ALX1) . The over-expression of homeobox A1 (HOXA1) and homeobox D3 (HOXD3) may result in beak deformity in chicks .
Sets of differently expressed genes in the deformed and normal beaks have been detected using digital gene expression (DGE) profile analysis based on high-sequencing technology. Of these genes, LOC426217, also known as claw keratin-like gene, was the most up-regulated in the deformed beak (log2-Ratio (deformed/normal) = 10.91) . Located on GGA 25, LOC426217 is a member of the keratin family, containing 417 base pairs with only one exon (Fig. 1). Keratin is crucial for maintaining normal cell morphology involved in the cytoskeleton remodeling keratin filaments and cytoskeletal signaling pathways. Change of its structure results in dysmorphic cells . The cytoskeleton is a complex of intracellular proteins that contribute to shape, support, and movement of cells . Up to now, less study was reported about this gene in chickens. Although highlighted in the DGE analysis, further study of this gene is still needed to verify its roles in beak malformation.
In the present study, qRT-PCR was used to detect the relative expression of LOC426217 gene in the deformed and normal beaks, to verify the results of DGE profiling. Tissue expression profile of this gene was also determined in 14 tissues of the birds. Eventually, LOC426217 was amplified and sequenced to seeking the SNPs and haplotypes related with the beak malformation.
Animals and samples collection
The Institutional Animal Care and Use Committee at Institute of Animal Science, Chinese Academy of Agricultural Sciences (IAS, CAAS) approved all procedures involving the use of animals. All efforts were made to minimize the suffering of animals following the animal care guidelines . The animals used in this study came from a pure-line stock of a local breed (Beijing-You) kept by IAS, CAAS (Beijing, China). They were incubated contemporaneously and housed under the same conditions.
The lower mandibles of the beaks were collected from 18 BJY chickens of 56 days of age: 9 with crossed beaks and 9 with normal beaks. Total RNA of the lower mandibles of the crossed and normal beaks above was collected for the verification of DGE profiling results using quantitative real-time PCR (qRT-PCR).
Three normal birds of 56 days of age were killed by stunning and exsanguination. Tissues samples including bursa of fabricius, beak, brain, breast, feather, heart, kidney, thigh, liver, lung, skin, small intestine, stomach, and testicle (50–100 mg) were rapidly collected and snap-froze in liquid nitrogen and storage at -80 °C. The RNA of these samples was used to determine the tissue expression profile of LOC426217.
Blood samples were collected from the brachial vein by venipuncture from 171 beak- deformed birds (deformed) and 164 normal ones (control). Based on the case-control study design, we selected these birds according to the phenotype of the birds without family structure. The beak-deformity birds were collected from two generations. The normal birds were selected randomly from the same generation. DNA was isolated from the blood samples and stored at -20 °C for the detection of SNPs located in LOC426217 gene.
DNA and RNA extraction and reverse transcription (RT)
Genomic DNA (gDNA) was extracted from blood samples using phenol-chloroform. Total RNA was isolated at 4 °C using the Trizol reagent (Invitrogen, USA). Any residual gDNA and protein were removed with Dnase I (TaKaRa, Japan) and RNA clean kit (TIANGEN, China). The purified RNA was dissolved (200-400 ng/mL, OD260/OD280 = 1.8–2.0), and stored at -80 °C. Total RNA was used for RT (in 20 μL final volume) following the manufacturer's instruction (Promega, USA). The cDNA was stored at -80 °C for subsequent qRT-PCR.
PCR amplification and qRT-PCR
PCR amplification was performed using PCR Gene Amplifier (Bio-Rad, USA) in a total volume of 25 μL which contains 12 μL of 2 × Taq PCR StarMix (GenStar, China), 1 μL (10 pmol) of each primer (Table 1), 1 μL of gDNA (50 ng) and 10 μL of ddH2O. After an initial denaturing for 2 min at 95 °C, there were 35 cycles of amplification (94 °C for 30 s, 60 °C for 30 s, and 72 °C for 90 s) and extension for 5 min at last. PCR products were detected by 1 % agarose gel electrophoresis for 15 min 120 V, stained with ethidium bromide, examined under UV light, and photographed. The PCR products were then sequenced by BGI (Beijing, China).
To validate the DGE results, qRT-PCR was performed to determine the expression of LOC426217 in the 9 deformed and 9 normal beaks using the ABI 7500 Real-time Detection System (Applied Biosystems, USA) and TaKaRa DRR018A reagents. Each 20 μL PCR mixture contained10 μL of SYBR Premix Ex Taq™ II, 0.8 μL (10 pM) of each primer (Table 2), 0.4 μL of ROX Reference Dye II (50×), 2 μL of cDNA (100 ng) and 6 μL of ddH2O. After an initial denaturing for 30 s at 95 °C, there were 40 cycles of amplification (95 °C for 5 s and 60 °C for 32 s), followed by thermal denaturing to generate melting curves to verify amplification specificity. β-actin was amplified in the same plates as endogenous control. Samples were assayed in triplicate for standard curves. PCR efficiency of the LOC426217 gene and β-actin was consistent. cDNA from normal beaks served as a standard control for tissue-specific expression profile study. The amplification efficiency of transcripts of interest and the internal standard (β-actin) were consistent. Dissociation curves verified that amplification was specific.
SNPs filter and genotyping
PCR amplification product of LOC426217 gene was then directly sequenced by BGI company (Beijing, China) using Sanger sequencing methods . The SNPs and amino acids were determined and filtered by the software DNAStar (Version 5.01).
The relative abundance of transcripts was calculated from 2−ΔΔCT . All data presented graphically are means ± SEM. The significance level was P < 0.05 or P < 0.01. Student’s t-tests were used to evaluate the relative expression differences of LOC426217 between the RNA samples of deformed and normal beaks. The ANOVA procedure of SAS 8.0 was used to assess the differences expression of LOC426217 in all the tissues. All SNPs were checked for Hardy-Weinberg Equilibrium (HWE) in both groups (P > 0.05 means equilibrium). Allele frequency, genotype frequency, and polymorphism information content (PIC, PIC < 0.25: low polymorphism; 0.25 < PIC < 0.5: moderate polymorphism; PIC > 0.5: high polymorphism) were calculated by PopGene (Version 1.31). Chi-square tests were used to evaluate the genotype frequency differences between deformed and normal beaks. Benjamini & Hochberg method was used for the Bonferroni correction . Linkage Disequilibrium (LD) pattern for the SNPs genotyped was plotted using Haploview (Version 4.2). The sliding window method was used to generate different haplotypes between two groups [22, 23].
Verification for LOC426217 gene of DGE results
To verify the previous DGE analysis, where LOC426217 was the most up-regulated in the deformed beaks, qRT-PCR was used to estimate the expression of this gene in 9 deformed and 9 normal beaks. As shown in Fig. 2, the relative expression of LOC426217 in deformed beaks was significantly higher than that in normal ones (P < 0.05). This was in good agreement with the DGE analysis.
Tissue expression profile of LOC426217
As shown in Fig. 3, LOC426217 gene was hardly expressed in brain, heart, bursa, or small intestine. The relative expression in beak was significantly higher than that in other tissues (P < 0.01). The results revealed that this gene maybe specifically expressed in beak tissue. In order to seek the structural changes of this gene, DNA sequencing was carried out subsequently.
DNA sequencing and SNPs detection
LOC426217 gene of 171 beak-deformed and 164 normal birds was amplified and sequenced. Eight SNPs were detected, including one locus in promoter region: -62G > T, and seven loci in coding region: 24 T > C, 36G > C, 192A > T, 204C > T, 222 T > C, 285G > T, and 363 T > C (Fig. 4). The loci in coding region were synonymous mutations resulting with no amino acid changing. These SNPs were presented in both beak-deformed and normal birds. The 18 birds used in the qRT-PCR analysis for DGE verification were also sequenced and their genotypes were shown in Table 3 and Additional file 1: Table S1.
Genetic diversity analysis
Genetic diversity was analyzed by PopGene (Version 1.31). All SNPs were firstly checked for HWE. In the normal birds (n = 164), T24C, G36C, A192T, T222C, G285T, and T363C were not in agreement with the HWE (P < 0.05), while in the birds with a deformed beak (n = 171), T24C, G285T, and T363C were not in agreement with the HWE. Two loci (G-62 T and C204T) of normal beaks and two loci (G-62 T and T363C) of deformed ones were low polymorphism (PIC < 0.25) (Tables 4 and 5).
Genotype frequency differences in the beak-deformed and normal birds
Chi-square tests were used to evaluate the genotype frequency differences of SNP loci in LOC426217 gene between two groups with a Bonferroni correction for the P-values. As shown in Table 6, the genotype frequency of G-62 T locus showed significant difference between two groups (P < 0.05), while the genotype frequencies of T24C, G36C, T222C and T363C loci showed highly significant differences (P < 0.01). There was no significant difference for the rest SNPs (P > 0.05).
In the normal birds, two haplotype blocks of LOC426217 were identified: one block with SNPs T24C and G36C, and one block with SNPs A192T, C204T, T222C, and G285T. In the birds with a deformed beak, one block with SNP G36C and A192T was identified (Table 7 and Fig. 5). The haplotype blocks of deformed birds were different from those of the normal ones.
The molecular genetic mechanism underlying beak deformity is likely to be very complex. Beak deformities of different forms (noticeably elongated, crossed, bent at right angles) have been documented in many wild birds. The molecular mechanism of beak deformity trait is not clear yet. For the wild birds, it is difficult to obtain the individuals for genetic study. Beak deformity was also found in various indigenous chickens and it is easy to collect individuals. This made the chicken a perfect model for the genetic study of this defect. Based on the previous DGE profiling and bioinformatics analyses, we identified a cluster of differentially expressed genes (DEGs) in the deformed and normal beaks. Some of the DEGs were quite extreme, especially the LOC426217 gene in the present study.
Validation and tissue expression profile
In order to validate the reliability and accuracy of DGE results, qRT-PCR was performed. The results showed that the expression of LOC426217 in the deformed beaks was significantly higher than that in the normal beaks, which was in great agreement with DGE analysis. The result of tissue expression profile also revealed that this gene was specifically expressed in beak tissue. LOC426217 is a member of the keratin family . Keratin is a key gene family for maintaining normal cell morphology . Variation of keratin structure can lead to beak deformity . It is also an intermediate filament protein that has essential functions in maintaining the structural integrity of epidermis and its appendages , presumably including the beak. In addition, keratin is the main composition of the chicken beak. This may be the main reason for its high expression in beak tissue. According to our early observation of beak anatomy, the lower mandibles of the beaks were abnormal/asymmetry. It was reported that avian keratin disorder could result in gross over growth of the rhamphotheca . The beak deformity caused by the excessive growth of one side of the lower mandibles maybe a result of abnormal high expression of LOC426217.
SNPs and haplotypes associated with beak deformity
Based on the case-control study design, the investigated birds used for sequencing were selected from the pure line BJY population according to the beak phenotype of the birds with no family structures. Since the percentage of beak-deformity birds in the chicken population was 1 %-3 %, it was not easy to collect enough birds for gene sequencing. The beak-deformity birds were collected from two generations (about 5000 birds in each generation) and the normal birds were selected randomly from the same generation.
Using direct sequencing, eight SNPs located in the LOC426217 gene were detected, including one locus in promoter region and seven loci in coding region. These loci were found in both birds with deformed and normal beaks. Four loci (T24C, G36C, T222C, and T363C) were previously released on the database of NCBI (http://www.ncbi.nlm.nih.gov) while the rest four SNPs found in this study were new. Furthermore, five SNPs, individuals with different genotypes showed significant differences between beak-deformed and normal birds (P < 0.05). Previous research demonstrated that promoter mutation is a kind of the mutations which can enhance or reduce the expression level of a gene and lead to different phenotypes [27, 28]. In this current study, one SNP (-G62T) was detected in the promoter region and the genotype frequency was significant different between two groups. This mutation could possibly affect the start, expression time and level of this gene. As shown in Table 3, the genotype of 9 normal birds were all TT, while 3 out of the 9 beak-deformed were TG. Further validation study with larger sample size is needed to validate the mRNA expression profile of different genotypes. Similarly, the mutations in coding region are also very important. As we all know, exon is a nucleotide sequence in DNA that carries the code for the final mRNA molecule and thus defines the amino acid sequence during protein synthesis. If the base pairs changed in coding region, the structure and function of protein could possibly change, such as missense mutation or silent mutation [29–32]. In this study, all the loci in coding region were synonymous mutations with no amino acid and protein changing. Synonymous mutations, which do not alter the protein sequence, have also been shown to affect protein function [33, 34] and play key roles in human diseases . Therefore, the mutations detected here may also affect this gene function. Sometimes, the phenotype altering or diseases occurring was not associated with one single SNP changing but related to several SNPs. Global patterns of DNA sequence variation (haplotypes) defined by common SNPs have important implications for identifying disease and traits [36, 37]. A previous study showed that the ALX1 haplotype has contributed to diversification of beak shapes among the Darwin’s finches and, thereby, to an expanded utilization of food resources . In this present study, different haplotypes were analyzed between two groups. It indicated that the beak deformity trait might be related to the specific haplotype block which was only existed in beak-deformed chickens. To sum up, SNPs and haplotypes described here were interesting and worthy of further study.
To the best of our knowledge, this is the first time that LOC426217 was studied as an important candidate gene for beak deformity in birds. The over-expression of LOC426217 may be the cause of beaks malformation. The genotype frequency of SNPs at G-62 T, T24C, G36C, T222C, and T363C loci showed significant differences between deformed and normal birds, and might be used as candidate SNP markers for this trait. The specific haplotype block in the deformed group could be served as a potential linkage marker for this trait. Further functional verification studies like over-expression or RNA interfere of LOC426217 during the embryonic development are required to reveal its roles in the beak malformation.
Availability of supporting data
All the supporting data are included in the manuscript as well as additional files in the supplementary section.
- ALX1 :
ALX homeobox 1
- BMP4 :
Bone morphogenetic protein 4
- CaM :
Differentially expressed genes
Digital gene expression
- FGF8 :
Fibroblast growth factor 8
- HOXA1 :
- HOXD3 :
- PIC :
Polymorphism information content
Quantitative real-time PCR
Single nucleotide polymorphisms
Seki Y, Bodde SG, Meyers MA. Toucan and hornbill beaks: a comparative study. Acta Biomater. 2010;6(2):331–43.
Rintoul DA. Beak deformity in a brown-headed cowbird, with notes on causes of beak deformities in birds. Kansas Ornithol Soc Bull. 2005;56:29–32.
Clabaut C, Herrel A, Sanger TJ, Smith TB, Abzhanov A. Development of beak polymorphism in the African seedcracker, Pyrenestes ostrinus. Evol Dev. 2009;11(6):636–46.
Handel CM, Pajot LM, Matsuoka SM, Hemert CV, Terenzi J, Talbot SL, et al. Epizootic of beak deformities among wild birds in Alaska: an emerging disease in North America? Auk. 2010;127(4):882–98.
Hemert CV, Handel CM. Beak deformities in Northwestern Crows: evidence of a multispecies epizootic. Auk. 2010;127(4):746–51.
Demery ZP, Chappell J, Martin GR. Vision, touch and object manipulation in Senegal parrots Poicephalus senegalus. Proc Biol Sci. 2011;278:3687–93.
Van Hemert C, Handel CM, O'Hara TM. Evidence of accelerated beak growth associated with avian keratin disorder in black-capped chickadees (Poecile atricapillus). J Wildlife Dis. 2012;48(3):686–94.
MacDonald ME, Abbott UK, Richman JM. Upper beak truncation in chicken embryos with the cleft primary palate mutation is due to an epithelial defect in the frontonasal mass. Dev Dynam. 2004;230(2):335–49.
Abzhanov A, Protas M, Grant BR, Grant PR, Tabin CJ. Bmp4 and morphological variation of beaks in Darwin’s finches. Science. 2004;305(5689):1462–5.
Wu P, Jiang T, Suksaweang S, Widelitz RB, Chuong C. Molecular shaping of the beak. Science. 2004;305(5689):1465–6.
Wu P, Jiang TX, Shen JY, Widelitz RB, Chuong CM. Morphoregulation of avian beaks: comparative mapping of growth zone activities and morphological evolution. Dev Dynam. 2006;235(5):1400–12.
Abzhanov A, Kuo WP, Hartmann C, Grant BR, Grant PR, Tabin CJ. The calmodulin pathway and evolution of elongated beak morphology in Darwin’s finches. Nature. 2006;442(7102):563–7.
Lamichhaney S, Berglund J, Almén MS, Maqbool K, Grabherr M, Martinez-Barrio A, et al. Evolution of Darwin/’s finches and their beaks revealed by genome sequencing. Nature. 2015;518(7539):371–5.
Jaszczak K, Malewski T, Parada R, Malec H. Expression of Hoxal and Hoxd3 genes in chicken embryos with exencephaly. J Anim Feed Sci. 2006;15(3):463.
Bai H, Zhu J, Sun Y, Liu RR, Liu N, Li DL, et al. Identification of genes related to beak deformity of chickens using digital gene expression profiling. PLoS One. 2014;9:e107050.
Steinert PM. Structure, function, and dynamics of keratin intermediate filaments. J Invest Dermatol. 1993;100(6):729–34.
Frixione E. Recurring views on the structure and function of the cytoskeleton: a 300-year epic. Cell Motil Cytoskel. 2000;46(2):73–94.
Institute for Laboratory Animal Research (US), Division on Earth and Life Studies (US), National Research Council (US). Wang JF, Zhou Y, Liu JH, Lu JQ, Ge LJ (Translation). Guide for the Care and Use of Laboratory Animals: Eighth Edition. Shanghai Scientific & Technical Publishers (SSTP), 2012.
Sanger F, Nicklen S, Coulson AR. DNA sequencing with chain-terminating inhibitors. Proc Natl Acad Sci U S A. 1977;74(12):5463–7.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2 − ΔΔCT method. Methods. 2001;25(4):402–8.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B. 1995;57:289–300.
Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–5.
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, et al. The structure of haplotype blocks in the human genome. Science. 2002;296(5576):2225–9.
Eckhart L, Dalla VL, Jaeger K, Ballaun C, Szabo S, Nardi A, et al. Identification of reptilian genes encoding hair keratin-like proteins suggests a new scenario for the evolutionary origin of hair. Proc Natl Acad Sci U S A. 2008;105(47):18419–23.
D'Alba L, Saranathan V, Clarke JA, Vinther JA, Prum RO, Shawkey MD. Colour-producing β-keratin nanofibres in blue penguin (Eudyptula minor) feathers. Biol Lett-UK. 2011;7(4):543–6.
Törmä H. Regulation of keratin expression by retinoids. Dermato-endocrinol. 2011;3(3):136–40.
Chang GB, Liu XP, Ma T, Xu L, Wang HZ, Li ZT, et al. A mutation in the NLRC5 promoter limits NF-κB signaling after Salmonella Enteritidis infection in the spleen of young chickens. Gene. 2015;568:117–23.
Mirhafez SR, Avan A, Pasdar A, Kazemi E, Ghasemi F, Tajbakhsh A, et al. Association of tumor necrosis factor-α promoter G-308A gene polymorphism with increased triglyceride level of subjects with metabolic syndrome. Gene. 2015;568:81–4.
Stanford PM, Halliday GM, Brooks WS, Kwok JB, Storey CE, Creasey H, et al. Progressive supranuclear palsy pathology caused by a novel silent mutation in exon 10 of the tau gene. Brain. 2000;123(5):880–93.
Houston RD, Cameron ND, Rance KA. A melanocortin-4 receptor (MC4R) polymorphism is associated with performance traits in divergently selected large white pig populations. Anim Genet. 2004;35(5):386–90.
Hou JX, An XP, Han P, Peng JY, Cao BY. Two missense mutations in exon 9 of caprine PRLR gene were associated with litter size. Anim Genet. 2015;46(1):87–90.
Salomons GS, Bok LA, Struys EA, Pope LL, Darmin PS, Mills PB, et al. An intriguing “silent” mutation and a founder effect in antiquitin (ALDH7A1). Ann Neurol. 2007;62(4):414–8.
Ramamurthi KS, Schneewind O. A synonymous mutation in Yersinia enterocolitica yopE affects the function of the YopE type III secretion signal. J Bacteriol. 2005;187(2):707–15.
Sauna ZE, Kimchi-Sarfaty C. Understanding the contribution of synonymous mutations to human disease. Nat Rev Genet. 2011;12(10):683–91.
Gartner JJ, Parker SC, Prickett TD, Dutton-Regester K, Stitzel ML, Lin JC, et al. Whole-genome sequencing identifies a recurrent functional synonymous mutation in melanoma. Proc Natl Acad Sci U S A. 2013;110(33):13481–6.
Patil N, Berno AJ, Hinds DA, Barrett WA, Doshi JM, Hacker CR, et al. Blocks of limited haplotype diversity revealed by high-resolution scanning of human chromosome 21. Science. 2001;294(5547):1719–23.
Salman M, Dasgupta S, Cholendra A, Venugopal PN, Lakshmi GL, Xaviour D, et al. MTNR1B gene polymorphisms and susceptibility to Type 2 Diabetes: a pilot study in South Indians. Gene. 2015;566(2):189–93.
This work was supported by 1. National Natural Science Foundation of China (31501949) (http://www.nsfc.gov.cn/); 2. Special fund for basic scientific research business of central public research institutes (2014ywf-yb-3). 3. National High Technology Research and Development Program (863 Program) (2011AA100305); and 4. Agricultural Science and Technology Innovation Program (ASTIP-IAS04), Chinese Academy of Agricultural Sciences, China (http://www.caas.cn/en). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors declare that they have no competing interests.
JLC, HB and JZ conceived and designed the experiments. HB, JZ and NL performed the experiments. DLL, FGX and YLL participated in data collection and sequence alignment. HB and YS analyzed the data. HB and YS drafted the manuscript. JLC, YS and HB revised the manuscript critically. All authors read and approved the final manuscript.