Molecular characterization of KRTAP9.2
RNA Isolated from skin samples of healthy animals was used for cDNA synthesis. The later was amplified (Fig 1) and used for sequencing. PCR amplified cDNA of KRTAP9.2 gene was sequenced and aligned using the BioEdit program (BioEdit 5.0.9). KRTAP9.2 protein sequences from different species have been downloaded from Uniprot database to study evolutionary relationship of
Capra hircus KRT9 protein sequence (Fig 2) with other species Muscle (
Edgar, 2004) alignment tool has been used for alignment and phylogenetic tree was constructed using parsimony method (Fig 3). Sequence based analysis of KRTAP9.2 gene between pashmina goat and other species showed that the pashmina goat sequences exhibited 98% similarity with
Capra hircus KRTAP9.2 (
NM-001285537.2)
, Capra hircus KRTAP9.2 (
XM-018065082.1)
, Capra hircus KRTAP9.2 (
AY510124.1) and
Pantholops hodgsonii KRTAP9.2 (XM-005964828.1)
. The similarity was 97% with
Ovis aries musimon KRTAP9.9 (
XM-012152535.1)
, Capra hircus KRTAP9.2 (
XM-018065084.1), Ovis aries KRTAP9.9 (
XM-012186088.2) and
Pantholops hodgsonii KRTAP9.2 (
XM-005964822.1)
. Pashmina goat, goat and sheep were found to cluster together and form a group apart from other mammals. The relationship displayed in the phylogenetic tree was found to correspond to mammalian classification group. The homology indicates the relationship of KRTAP9.2 gene with fiber related traits.
Structural analysis of KRTAP 9.2
Due to lack of availability of template structures in Protein Data Bank (PDB) multi-threading and other method have been used for structure prediction for KRTAP9.2 protein. RaptorX (
Kallberg, 2012) predicts protein structure when close homologs are not available in the protein data bank (PDB). In our study, total of 5 structures were predicted by RaptorX using templates given in supplementary 1, best structure (Fig 4) is predicted using 3QWQ (p-value = 3.28e-04) 285 (77%) residues are modeled correctly by RaptorX. Ramachandran plot suggests 88.1% residues are in favored region and 8% residues are in allowed region of predicted model. However, QUARK (
Xu and Zhang 2012) and LOMET (
Wu and Zhang 2007) have also been used in structure prediction of
Capra hircus KRTAP9.2 protein. Quark is suitable for structure prediction without known homologs using amino acid sequence only. Lomet generates 3D model by identifying small high scoring targets on query sequence. Total of 20 structures were predicted by QUARK and LOMET server, we have used different parameters to check the validity of predicted structure like RMSD value, C-Score, TM-value and ramachandran values, all parameters suggest structure predicted by RaptorX is better in our case.
Our study identified 2 domains in KRT9.2 protein sequence, using RaptorX Binding server, a total of 13 binding sites, 5 in domain-1 and 8 in domain-2 were discovered, in addition possible ligands and binding residues are given in Table 2 below.
Expression analysis
The expression pattern of KRTAP9.2 gene at Anagen and Telogen stages of pashmina fiber development was determined by quantitative Real Time PCR. KRTAP9.2 gene was found to express remarkably at Telogen stage compared to Anagen stage in Changthangi goats of both the regions. In Changthangi goats of Ladakh region, the expression was found significantly higher at Telogen stage (70.875± 7.402 fold) compared to that at Anagen stage (2.039±0.904 fold) with P value less than 0.0001, likewise in Kashmir region, its expression was also higher at Telogen stage (58.125±54.556 fold) (Fig 5) compared to that at Anagen stage (1.935±1.513), however the difference was not statistically significant (P=0.149123) as described in Table 1. Comparative expression of KRTAP9.2 gene at Telogen phase was observed to be lower in goats maintained at MRSS and G Kashmir than the goats maintained at high altitude Ladakh region, this finding coincides with the delay in cashmere fiber development in goats at MRSS and G.
Hair follicle is a dynamic mini-organ of the body, formed due to interaction between ectoderm and mesoderm (
Botchkarev and Paus 2003;
Millar 2002;
Schneider et al., 2009; Hardy, 1992). It is broadly differentiated into three parts including the outer and inner root sheath, hair shaft and matrix. Formation of hair follicle occurs in a cyclical fashion, involving growth (Anagen), regression (Catagen) and resting phases (Telogen), followed by the shedding of the hair shaft. The process involves the expression of both hair keratin intermediate filament (KIFs) proteins and their keratin associated proteins (KRTAPs)
(Shimomura et al., 2002; Rogers et al., 2008; Rogers et al., 2002; Powell et al., 1995; Pruett et al., 2004). Both keratins (KRTs) and keratin associated proteins (KRTAPs) are the important structural proteins of the secondary hair follicle, whose content in hair follicle remarkably determines the quality and fineness of cashmere
(Jin et al., 2011). Being major constituent of hair proteome, KRTAPs have been long associated with phenotypic differences in hair and wool
(Khan et al., 2014). So far different KRTAPs have been implicated in development of cashmere fiber, most notably KRTAP6.1
(Liu et al., 2007), KRTAP7.1
(Jin et al., 2011), KRTAP8.2 (
Liu H et al., 2007) and KRTAP1.4
(Shah et al., 2013). In the present study, KRTAP9.2 variant of KRTAPs was sequence analysed to construct a phylogenetic tree which revealed the similarity between
Changthangi goat and
Capra hircus variants,
Pantholops hodgsonii and
Ovis aries musimon. Changthangi goat was found to cluster with goat (
Capra hircus) and sheep, forming a separate group, distinct from other mammals. The present study provides the first report about the phylogenetic analysis of KRTAP9.2 gene in Changthangi goat of Ladakh and Kahmir region of J&K, India. Such phylogeny based studies are beneficial for identifying the essential morphogenetic features in different species or sub species (
Avise, 2000). In fact phylogenetic reconstruction is one of the most important subjects in molecular evolution, in particular, they are important in understanding the molecular mechanism of evolutionary events that lead to genetic as well as species diversity. Besides, phylogenetics is viewed as an important tool for performing selective breeding between compatible species for better productivity and quality (
Tabbaa and Al-Atiyat, 2009).
Real time PCR based analysis of KRTAP9.2 gene revealed that it is expressed significantly at Telogen stage (70.875±7.402 folds), as compared to Anagen stage (2.039±0.904 folds) in the Chanthangi goats of J&K and Ladakh region (P<0.00001). These results are in consonance with those reported earlier by
Wang et al., (2014). It has been shown earlier that mRNAs of KRTAP variants like KRTAP7.1 and KRTAP8.2 do not get expressed in the tissues like kidney, heart, lung, liver and spleen, but are inimitably expressed in hair follicles participating in hair formation
(Jin et al., 2011). Rogers (2004) demonstrated that a large number of KRTAPs are expressed only in hair follicles and not in other tissues. From the present study, it may be inferred that KRTAP9.2 gene down regulates pashmina fiber growth owing to its high expression at telogen stage.
Wang et al., (2014) demonstrated similar findings while investigating expression pattern of KRTAP9.2 gene in
Shaanbei white cashmere goat. Also, our results indicate lower expression of KRTAP9.2 gene at Telogen phase in Changthangi goats maintained at MRSSandG Kashmir than the goats maintained at high altitude Ladakh region. Previously it has been established that duration of solar radiation has an important effect on the thermo-physiological and the growth parameters of goats (
Al-Tamimi, 2005). In particular, development of hair cycle is critically influenced by the photoperiod
(Liu et al., 2016). Thus photoperiod could be an important factor regulating the altered expression of KRTAP9.2 gene in pashmina/Changthangi goats. Secondly, since KRTAP9.2 is associated with fiber traits
(Wang et al., 2014), its differential expression in the Pashmina goats of Ladakh and Kashmir region might correspond to the difference in fiber traits between two geographical regions. Characteristics of pashmina (Cashmere) fiber change significantly in Cashmere goats belonging to different geographical regions or species probably due to differential expression of KRT/KRTAP genes (
Jeffrey and Warren 2000) in response to environmental conditions. KRT/KRTAPs being essential for the process of keratinisation are highly evolved in response to changing environmental conditions
(Khan et al., 2014). It is noteworthy to mention that keratinization is an essential process safeguarding animal body by forming a protective barrier between the body and external environment. Climatic factors like cold stress requires more adaptability towards environment, which may be accomplished by altered expression of specific genes for example, enhanced KRTs/KRTAPs expression causing more keratinisation. Changthangi goats from Ladakh region are naturally more adapted to extreme cold arid conditions, compared to those from Kashmir region which may correspond to the difference in their fiber quality. Certain species of goats like dairy goats do not express KRTAP9.2 at all
(Wang et al., 2012) thus pin pointing towards the possible key role of this gene in the regulation of cashmere growth in Changthangi goats. Furthermore, the data of changthangi goat maintained in the J&K had outlier that contributed to a marked standard error, making the difference nonsignificant. Such outliers cannote be under estimated.