TMT based deep proteome analysis of buffalo mammary epithelial cells and identification of novel protein signatures during lactogenic differentiation.

The lactating mammary gland harbours numerous matured alveoli with their lumen surrounded by differentiated mammary epithelial cells (MECs), which are exclusively involved in milk synthesis and secretion. Buffalo (Bubalus bubalis) is the second major milk-producing animal, and its physiology is different from cattle. The complete protein machinery involved in MECs differentiation is still not defined in ruminants, in particular, buffalo. Therefore, we have studied the differential expression of regulated proteins in the in vitro grown buffalo MECs (BuMECs) at different time points (on 3, 6, 12, and 15 days) of their differentiation in the presence of lactogenic hormones. TMT-based MS analysis identified 4,934 proteins; of them, 681 were differentially expressed proteins (DEPs). The principal component analysis suggested a highly heterogeneous expression of DEPs at the four-time points of hormone treatment, with most of them (307) attained the highest expression on 12 days. Bioinformatics analysis revealed the association of DEPs with 24 KEGG pathways. We observed few new proteins, namely ABCA13, IVL, VPS37, CZIB, RFX7, Rab5, TTLL12, SMEK1, GDI2, and TMEM131 in BuMECs. The function of one of the highly upregulated proteins, namely involucrin in the differentiation of BuMECs was confirmed based on biochemical inhibition assay. The results further conclude that the proteins with higher abundance can be considered as the potential biomarkers for differentiation, and they may have a significant association with the lactation process in buffalo too. The proteome dataset obtained can be used to understand the species-specific variations among other lactating animals.

Crystal structure of breast regression protein 39 (BRP39), a signaling glycoprotein expressed during mammary gland apoptosis, at 2.6 Å resolution.

Breast regression protein 39 (BRP39) is a 39 kDa protein that is a member of chitolectin class of glycosyl hydrolase family 18 (GH18). High expression levels of BRP39 have been detected in breast carcinoma. It helps in proliferation of cells during the progression of this disease and may act as a signaling factor. BRP39 may act as a potential candidate for rational structure-based drug design against breast carcinoma. In this study, we report the crystal structure of mouse recombinant BRP39 expressed in E. coli. The structure was solved by molecular replacement and refined to 2.6 Å resolution. The overall structure of BRP39 consisted of two globular domains: a large (ß/α)8 triosephosphate isomerase (TIM) barrel domain and a small (α + ß) domain. Three non-proline cis-peptides were detected in the sugar-binding cleft of BRP39, including Ser57-Phe58, Leu141-Tyr142, and Trp353-Ala354. The latter residues were conserved in other GH18 family members. It was notable that the conformation of critical Trp100 residue within the sugar-binding cleft was oriented away from the barrel. The side-chain conformation was found to be similar to that observed in chitinases, however, it was oriented into the barrel in other chitinase-like proteins (CLPs). The conformation of this critical residue may have significant implications in sugar binding. Further, two amino acid substitutions were observed in the sugar-binding groove of BRP39. The conserved Asn100 and Arg263 in Hcgp39 and other CLPs proteins (SPX-40 structures) were substituted by Lys101 and Lys264 in BRP39 which may have a significant impact on the sugar-binding properties.

New insights into the catalytic inactivity of mammary gland protein-40, a chitinase-like protein expressed during mammary gland involution.

MGP-40 is a mammary gland-specific glycoprotein which is expressed during involution and is an important marker for mammary gland apoptosis. It is an inactive chitinase-like protein belonging to Glycosyl Hydrolase family 18. The present study reports sequence characterization, tissue-specific expression analysis, production of recombinant MGP-40 and its mutant (A117D and L119E) in both E. coli and COS1 cells for their chitin-binding and chitinase activity analysis. The cDNA of buffalo MGP-40 was cloned and sequenced which corresponded to 1803 bp with an open reading frame of 1152 bp (361 aa), signal sequence of 63 bp (21 aa), 5' and 3' UTR of 144 bp and 507 bp, respectively. The 3' UTR analysis revealed potential sites for high level expression and stability during involution. The half-life of buffalo MGP-40 was found to be 11.7 h. MGP-40 was highly expressed in mammary gland followed by small intestine, spleen and mammary epithelial cells. The purified recombinant MGP-40 and its mutant expressed in E.coli were observed to bind chitin efficiently, however, no chitinase activity was observed. Further, chitinase activity was also not observed by expressing mutant recombinant MGP-40 in COS1 cells ruling out the possible role of post-translational modifications. Structure-based in-silico mutagenesis by FoldX algorithm showed a drastic decrease in overall fold stability which might be a possible reason for inability to recover its activity. Therefore, chitinase activity could not be restored in MGP-40 even after reverting back two critical residues in active site which may be due to detrimental effect of mutations on structural stability.

Functional analysis of recombinant buffalo lactoferrin and monoferric lobes and their cytotoxic effect on buffalo mammary epithelial cells.

Lactoferrin (Lf) has been involved in diverse type of cellular activities and its biochemical properties are species specific. Lf is a bilobal molecule in which each lobe binds with one Fe2+/Fe3+ ion. A lot of physiological effects of Lf are regulated by its iron binding and release properties; however these properties are species-specific. To understand the iron-binding, thermal stability and cytotoxic effect of buffalo Lf (buLf) and contribution of individual N- and C-terminal lobes therein, buLf and the truncated monoferric lobes were expressed in Kluyveromyces lactis or Pichia pastoris yeast expression systems. The iron-uptake/release behavior and thermal stability of recombinant buLf was observed similar to the Lf purified from buffalo milk. Supplementation of recombinant buLf to the buffalo mammary epithelial cells (BuMEC) culture decreased their proliferation and the cell viability in a dose dependent manner. The cell growth decreased by 37% at 1.0 mg/ml Lf. C-lobe decreased the viability of BuMEC by 15% at 1 mg/ml. The C-lobe showed greater cytotoxic effect against BuMEC in comparison to N-lobe. buLf caused a reduced expression of the casein in BuMEC. At 1.0 mg/ml of buLf, CSN2 transcript level was reduced by 74% and 78% in the normal and hormone free media, respectively. The expression of IL-1ß gene in BuMEC increased by 4-5 fold in the presence of 1.0 mg/ml of Lf. The effect was similar to that observed in the involutory mammary gland, suggesting the role of elevated level of Lf in remodeling of buffalo mammary tissue during involution.

Buffalo Leukemia Inhibitory Factor Induces Differentiation and Dome-Like Secondary Structures in COS-1 Cells.

This study aimed to understand the molecular characteristics of buffalo leukemia inhibitory factor (BuLIF) and the generation of a stably transfected COS-1_BuLIF cell line for its functional characterization. Cumulus cells, isolated from oocytes, were separated, and total cDNA was prepared. The BuLIF gene was ligated into the cloning vector pJET1.2/blunt and expression vector pAcGFP-N1 which was transfected into COS-1 cells and confirmed by qRT-PCR and Western blot. BuLIF was immunoprecipitated and evaluated through a MTT assay. qRT-PCR of STAT3 was performed. The multiple sequence alignment of BuLIF showed high similarity with sheep (98.77%) and cattle (96.62%) compared with other species. The BuLIF gene has an open reading frame of 609 nucleotides coding for 202 amino acids. BuLIF was integrated into the genome of COS-1 cells and resulted in the formation of dome-like secondary structures which are indicative of its functional role mediated through STAT3 proteins. In conclusion, this cell line is suitable for understanding LIF-mediated biological functions.

Identification of potential protein biomarkers for early detection of pregnancy in cow urine using 2D DIGE and label free quantitation.

BACKGROUND: An early, reliable and noninvasive method of early pregnancy diagnosis is prerequisite for efficient reproductive management in dairy industry. The early detection of pregnancy also help in to reduce the calving interval and rebreeding time which is beneficial for industries as well as farmers. The aim of this work is to identify potential biomarker for pregnancy detection at earlier stages (16-25 days). To achieve this goal we performed DIGE and LFQ for identification of protein which has significant differential expression during pregnancy. RESULTS: DIGE experiment revealed a total of eleven differentially expressed proteins out of which nine were up regulated having fold change ≥1.5 in all time points. The LFQ data analysis revealed 195 differentially expressed proteins (DEPs) out of 28 proteins were up-regulated and 40 down regulated having significant fold change ≥1.5 and ≤0.6 respectively. Bioinformatics analysis of DEPs showed that a majority of proteins were involved in regulation of leukocyte immunity, endopeptidase inhibitor activity, regulation of peptidase activity and polysaccharide binding. CONCLUSION: This is first report on differentially expressed protein during various time points of pregnancy in cow to our best knowledge. In our work, we identified few proteins such MBP, SERPIN, IGF which were differentially expressed and actively involved in various activities related to pregnancy such as embryo implantation, establishment and maintenance of pregnancy. Due to their involvement in these events, these can be considered as biomarker for pregnancy but further validation of is required.

Molecular cloning and production of caprine recombinant Oct4 protein for generation induced pluripotent stem cells.

Oct4, pluripotency marker and transcription factor, expresses in embryonic stem cells. It plays a pivotal role in determination of stem cells fate. Up and down regulation of Oct4 causes differentiation of embryonic stem cells. It is one of the main transcription factors which remained concerned in every study related to induced pluripotent stem cell. Here, we report the production of goat Oct4 protein using plasmid and lentiviral based vectors. Firstly, Oct4 ORF was cloned in pAcGFP1-N1 plasmid vector and positive clones were screened with colony PCR. Oct4 was over-expressed in CHO-K1 cell line and expression was confirmed by observing green florescent protein expression in CHO-K1 cells. Secondly, Oct4 lentiviral expression construct has been prepared using pLenti-gw vector. Oct4 ORF was cloned into pLenti4/V5-DEST vector and viral particles were produced in 293FT cells. Oct4 viral particles were used to infect goat fibroblast cells. Oct4 expression was observed and confirmed in transfected goat fibroblast cells using RT-PCR. Detection of Oct4 protein in western blotting assay affirmed the capacity of over-expression of our Oct4 lentiviral vector. The lentiviral expression construct and recombinant Oct4 protein may be used for reprogramming of somatic cell into induced pluripotent stem cell.

Molecular cloning, sequence characterization and recombinant expression of Nanog gene in goat fibroblast cells using lentiviral based expression system.

Nanog is a homeodomain containing protein which plays important roles in regulation of signaling pathways for maintenance and induction of pluripotency in stem cells. Because of its unique expression in stem cells it is also regarded as pluripotency marker. In this study goat Nanog (gNanog) gene has been amplified, cloned and characterized at sequence level with successful over-expression in CHO-K1 cell line using a lentiviral based system. gNanog ORF is 903 bp long which codes for Nanog protein of size 300 amino acids (aas). Complete nucleotide sequence shows some evolutionary mutation in goat in comparision to other species. Protein sequence of goat is highly similar to other species. Overall, gNanog nucleotide sequence and predicted protein sequence showed high similarity and minimum divergence with cattle (96 % identity/4 % divergence) and buffalo (94/5 %) while low similarity and high divergence with pig (84/15 %), human (81/23 %) and mouse (69/40 %) indicating evolutionary closeness of gNanog to cattle and buffalo. gNanog lentiviral expression construct was prepared for over-expression of Nanog gene in adult goat fibroblast cells. Lentiviral expression construct of Nanog enabled continuous protein expression for induction and maintenance of pluripotency. Western blotting revealed the expression of Nanog gene at protein level which supported that the lentiviral expression system is highly promising for Nanog protein expression in differentiated goat cell.

Stage specific expression of ATP-binding cassette and solute carrier superfamily of transporter genes in mammary gland of riverine buffalo (Bubalus bubalis).

In the present study, expression level of various ATP-binding cassette (ABC) viz., ABCA1, ABCA7, ABCG1, ABCG2, and ABCG5; associated transcription factors viz., SREBF1, LXRα (NR1H3), PPARA, and Solute Carriers (SLC); or Glucose transporters (GLUT) viz., SLC2A1(GLUT1), SLC2A4 (GLUT4), SLC2A8 (GLUT8), and SLC2A12 (GLUT12) superfamily of transporters were compared across physiological stages of buffalo mammary gland. The relative expression of ABCA1, and ABCG1 was significantly (p < 0.05) higher in mammary gland of heifer followed by involution and lactation stages. Similarly, ABCA7 gene expression was highest in heifer mammary gland followed by lactation and involution stages. ABCG2 gene expression was significantly (p < 0.05) high in lactating mammary gland in comparison to involution and heifer stages. On the other hand, ABCG5 gene expression was highest in involuting mammary gland followed by lactation and involution stages. Additionally, the expression of LXRα SREBF1, and PPARA which are known to regulate some of the ABC tranporters were also analyzed. The expression of LXRα gene was high in involuting as compared to lactating mammary gland. In contrast, SREBF1 and PPARA expression was significantly (p < 0.05) high in lactating mammary gland. Among the several SLC transporters studied, SLC2A1, SLC2A4, and SLC2A8 showed significant (p < 0.05) higher expression during lactation stage, whereas SLC2A12 expression was greater during heifer stage suggesting SLC2A1, SLC2A4, and SLC2A8 to be the major transporters associated with glucose uptake in buffalo mammary gland. The expression profile of (lactoferrin) LTF, known to be expressed at high level in mammary gland during involution was also studied. As expected, its expression was significantly (p < 0.05) higher during involution in comparison to lactating mammary gland.in buffaloes as well. The inclusion of LTF as a control gene further provided the confidence in the buffalo mammary gland expression data generated in the present study. This study thus helped to provide information about the distinct expression pattern of various transporters and their regulators in buffalo mammary gland during different physiological states.

Proteome analysis of functionally differentiated bovine (Bos indicus) mammary epithelial cells isolated from milk.

Mammary gland is made up of a branching network of ducts that end in alveoli. Terminally differentiated mammary epithelial cells (MECs) constitute the innermost layer of aveoli. They are milk-secreting cuboidal cells that secrete milk proteins during lactation. Little is known about the expression profile of proteins in the metabolically active MECs during lactation or their functional role in the lactation process. In the present investigation, we have reported the proteome map of MECs in lactating cows using 2DE MALDI-TOF/TOF MS and 1D-Gel-LC-MS/MS. MECs were isolated from milk using immunomagnetic beads and confirmed by RT-PCR and Western blotting. The 1D-Gel-LC-MS/MS and 2DE-MS/MS based approaches led to identification of 431 and 134 proteins, respectively, with a total of 497 unique proteins. Proteins identified in this study were clustered into functional groups using bioinformatics tools. Pathway analysis of the identified proteins revealed 28 pathways (p < 0.05) providing evidence for involvement of various proteins in lactation function. This study further provides experimental evidence for the presence of many proteins that have been predicted in annotated bovine genome. The data generated further provide a set of bovine MEC-specific proteins that will help the researchers to understand the molecular events taking place during lactation.

Postbiotic Lipoteichoic acid of probiotic Lactobacillus origin ameliorates inflammation in HT-29 cells and colitis mice.

Lipoteichoic acid (LTA) is a key surface component of probiotic lactobacilli that is involved in important cellular functions including cross talk with the host immune cells. In this study, the anti-inflammatory and ameliorative properties of LTA from probiotic lactobacilli strains were assessed in in vitro HT-29 cells and in vivo colitis mice. The LTA was extracted with n-butanol and its safety was confirmed based on its endotoxin content and cytotoxicity in HT-29 cells. In the Lipopolysaccharide stimulated HT-29 cells, the LTA from the test probiotics evoked a visible but non-significant increase in IL-10 and decrease in TNF-α levels. During the colitis mice study, probiotic LTA treated mice showed substantial improvement in external colitis symptoms, disease activity score and weight gain. The treated mice also showed improvements in key inflammatory markers such as the gut permeability, myeloperoxidase activity and histopathological damages in colon, although non-significant improvements were recorded for the inflammatory cytokines. Furthermore, structural studies by NMR and FTIR revealed increased level of D-alanine substitution in the LTA of LGG strain over MTCC5690. The present study demonstrates the ameliorative effect of LTA as a postbiotic component from probiotics which can be helpful in building effective strategies for combating gut inflammatory disorders.

Current status and emerging role of glutathione in food grade lactic acid bacteria.

Lactic acid bacteria (LAB) have taken centre stage in perspectives of modern fermented food industry and probiotic based therapeutics. These bacteria encounter various stress conditions during industrial processing or in the gastrointestinal environment. Such conditions are overcome by complex molecular assemblies capable of synthesizing and/or metabolizing molecules that play a specific role in stress adaptation. Thiols are important class of molecules which contribute towards stress management in cell. Glutathione, a low molecular weight thiol antioxidant distributed widely in eukaryotes and Gram negative organisms, is present sporadically in Gram positive bacteria. However, new insights on its occurrence and role in the latter group are coming to light. Some LAB and closely related Gram positive organisms are proposed to possess glutathione synthesis and/or utilization machinery. Also, supplementation of glutathione in food grade LAB is gaining attention for its role in stress protection and as a nutrient and sulfur source. Owing to the immense benefits of glutathione, its release by probiotic bacteria could also find important applications in health improvement. This review presents our current understanding about the status of glutathione and its role as an exogenously added molecule in food grade LAB and closely related organisms.

Molecular cloning, sequence characterization and heterologous expression of buffalo (Bubalus bubalis) oviduct-specific glycoprotein in E. coli.

Oviductin is a high molecular weight oviduct-specific glycoprotein secreted by the non-ciliated epithelial cells of oviduct during estrous cycle and early pregnancy. It plays an important role during fertilization and early embryonic development. The oviductin gene from oviductal tissues of buffalo was successfully cloned and sequenced. The sequence analysis revealed that buffalo and cattle oviductin share very high homology between their cDNA sequences. The predicted amino acid sequences of the buffalo oviductin exhibited the highest percent of identity of 97 % with bovine followed by 94 % with goat, 93 % with sheep, 78 % with porcine, 72 % with human, 67 % with hamster and rabbit and 65 % with mouse. Oviductin was also observed to share high similarity with the mammalian chitinase, however oviductins do not show chitinase activity due to Glu→Ile mutation in the active site responsible for chitinase activity. The phylogenetic tree based on amino acid sequences of oviductin indicated that buffalo oviductin was closely related to its cattle counterpart, and this clustering is in accordance with the classic taxonomic relationship. Tissue specific expression of the transcripts for buffalo oviductin revealed a high level expression in oviduct and ovary followed by testis, mammary gland, kidney, while in mammary epithelial cells and liver its expression was very low. The full length matured oviductin and its domains constituting chitinase-like domain and mucin-like domain were cloned into pET and pGEX series of expression vectors and over expressed in E. coli. The soluble recombinant oviductin was successfully purified to homogeneity. Full length recombinant oviductin was expressed partially in soluble form, where as the chitinase-like and mucin-like domains of oviductin were expressed in insoluble form and aggregating to form inclusion bodies at both 37 and 16 °C induction temperatures.

Ameliorative Effect of Surface Proteins of Probiotic Lactobacilli in Colitis Mouse Models.

The increase in concern from viable cells of probiotics specifically in acute inflammatory conditions has led to the emergence of the concept of postbiotics as a safer alternative therapy in the field of health and wellness. The aim of the present study was to evaluate the efficacy of surface proteins from three probiotic strains in dextran sodium sulfate and trinitrobenzenesulphonic acid = induced colitis mouse models. The molecular weight of total surface proteins extracted from the three probiotic strains ranged from ∼25 to ∼250 kDa with the presence of negligible levels of endotoxins. Surface layer proteins (SLPs) (∼45 kDa) were found to be present only in the Lactobacillus acidophilus NCFM strain. In the in vivo study, significant differences were not observed in the weight loss and general appetite, however, the decrease in colon length was apparent in TNBS colitis control mice. Further, the administration of these surface proteins significantly reversed the histopathological damages induced by the colitogens and improved the overall histological score. The oral ingestion of these surface proteins also led to a decrease in myeloperoxidase activity and TNF-α expression while the IL-10 levels significantly increased for the strain NCFM followed by MTCC 5690 and MTCC 5689. Overall, the present study signifies the ameliorative role of probiotic surface proteins in colitis mice, thereby, offering a potential and safer alternative for the management of inflammatory bowel disorders.

In-depth proteome analysis of more than 12,500 proteins in buffalo mammary epithelial cell line identifies protein signatures for active proliferation and lactation.

The mature mammary gland is made up of a network of ducts that terminates in alveoli. The innermost layer of alveoli is surrounded by the differentiated mammary epithelial cells (MECs), which are responsible for milk synthesis and secretion during lactation. However, the MECs are in a state of active proliferation during pregnancy, when they give rise to network like structures in the mammary gland. Buffalo (Bubalus bubalis) constitute a major source of milk for human consumption, and the MECs are the major precursor cells which are mainly responsible for their lactation potential. The proteome of MECs defines their functional state and suggests their role in various cellular activities such as proliferation and lactation. To date, the proteome profile of MECs from buffalo origin is not available. In the present study, we have profiled in-depth proteome of in vitro cultured buffalo MECs (BuMECs) during active proliferation using high throughput tandem mass spectrometry (MS). MS analysis identified a total of 8330, 5970, 5289, 4818 proteins in four sub-cellular fractions (SCFs) that included cytosolic (SCF-I), membranous and membranous organelle's (SCF-II), nuclear (SCF-III), and cytoskeletal (SCF-IV). However, 792 proteins were identified in the conditioned media, which represented the secretome. Altogether, combined analysis of all the five fractions (SCFs- I to IV, and secretome) revealed a total of 12,609 non-redundant proteins. The KEGG analysis suggested that these proteins were associated with 325 molecular pathways. Some of the highly enriched molecular pathways observed were metabolic, MAPK, PI3-AKT, insulin, estrogen, and cGMP-PKG signalling pathway. The newly identified proteins in this study are reported to be involved in NOTCH signalling, transport and secretion processes.

Cloning, expression, characterization and crystallization of BRP39, a signalling glycoprotein expressed during mammary gland apoptosis.

Breast regression protein (BRP39) is a glycoprotein, which is expressed during mammary gland involution in mouse. The physiological function of BRP39 is not known. High levels of expression of BRP39 have also been associated with breast cancer development. In the present investigation a cDNA encoding rBRP39 (recombinant BRP39) was cloned by PCR techniques. It consists of 1,143 nucleotides and encodes an open reading frame of 381 amino acid residues including a signal sequence of 21 amino acids. Recombinant BRP39 was produced in E. coli in a soluble form at low temperature (15 degrees C). Expression and purification of rBRP39 was confirmed by western blot analysis. Purified rBRP39 showed high chitin-binding activity but no chitinase activity. The lack of chitinase activity may be attributed to the mutation of critical active site residue Glu120 to Leu120 and Asp118 to Ala118 in BRP39. However, a mutant in which the residue was reverted back to Glu, by site directed mutagenesis, displayed no chitinase activity. Purified recombinant BRP39 was crystallized and the crystals diffracted X-rays to 2.8A resolution. The crystals belonged to the space group C2 with unit cell parameters a=130.4A, b=81.3A, c=229.2A, beta=105.9 degrees. The structure refinement is in progress.

Structural and functional characterization of buffalo oviduct-specific glycoprotein (OVGP1) expressed during estrous cycle.

Oviduct-specific glycoprotein (OVGP1) is a high molecular weight chitinase-like protein belonging to GH18 family. It is secreted by non-ciliated epithelial cells of oviduct during estrous cycle providing an essential milieu for fertilization and embryo development. The present study reports the characterization of buffalo OVGP1 through structural modeling, carbohydrate-binding properties and evolutionary analysis. Structural model displayed the typical fold of GH18 family members till the boundary of chitinase-like domain further consisting of a large (ß/α)8 TIM barrel sub-domain and a small (α+ß) sub-domain. Two critical catalytic residues were found substituted in the catalytic centre (Asp to Phe118, Glu to Leu120) compared with the active chitinase. The carbohydrate-binding groove in TIM barrel was lined with various conserved aromatic residues. Molecular docking with different sugars revealed the involvement of various residues in hydrogen-bonding and non-bonded contacts. Most of the substrate-binding residues were conserved except for a few replacements (Ser13, Lys48, Asp49, Pro50, Asp167, Glu199, Gln272 and Phe275) in comparison with other GH18 members. The residues Trp10, Trp79, Asn80, Gln272, Phe275 and Trp334 were involved in recognition of all six ligands. The α+ß sub-domain participated in sugar-binding through Thr270, Gln272, Tyr242 and Phe275. The binding assays revealed significant sugar-binding with purified native and recombinant OVGP1. Phylogenetic analysis revealed that OVGP1 was closely related to AMCases followed by other CLPs and evolution of OVGP1 occurred through several gene duplications. This is the first study describing the structural characteristics of OVGP1 that will further help to understand its interaction with gametes to perform crucial reproductive functions.

Comparative transcriptome analysis reveals the genetic basis of coat color variation in Pashmina goat.

The genetics of coat color variation remains a classic area. Earlier studies have focused on a limited number of genes involved in color determination; however, the complete set of trait determinants are still not well known. In this study, we used high-throughput sequencing technology to identify and characterize intricate interactions between genes that cause complex coat color variation in Changthangi Pashmina goats, producer of finest and costly commercial animal fiber. We systematically identified differentially expressed mRNAs and lncRNAs from black, brown and white Pashmina goat skin samples by using RNA-sequencing technique. A pairwise comparison of black, white and brown skin samples yielded 2479 significantly dysregulated genes (2422 mRNA and 57 lncRNAs). Differentially expressed genes were enriched in melanin biosynthesis, melanocyte differentiation, developmental pigmentation, melanosome transport activities GO terms. Our analysis suggested the potential role of lncRNAs on color coding mRNAs in cis and trans configuration. We have also developed online data repository as a component of the study to provide a central location for data access, visualization and interpretation accessible through http://pcd.skuastk.org/ .

Primary structures of different isoforms of buffalo pregnancy-associated glycoproteins (BuPAGs) during early pregnancy and elucidation of the 3-dimensional structure of the most abundant isoform BuPAG 7.

Pregnancy-associated glycoproteins (PAGs) are expressed during pregnancy by the trophoectodermal cells of fetus. Presence of PAGs in dam's circulation has been widely used in pregnancy diagnosis. The present study reports the identification and characterization of different PAG isoforms in buffalo during early stages of pregnancy. The PAG mRNAs isolated from fetal cotyledons (Pregnancy stages: 45, 75 and 90 days) were successfully cloned in pJET1.2 vector and transformed in E. coli. A total of 360 random clones were sequenced and correlated with their stages of expression. A total of 12 isoforms namely, BuPAG 1, 2, 4, 6, 7, 8, 9, 13, 15, 16, 18 and one new isoform were identified. BuPAG 7 was found as the most abundant isoform in all three stages followed by BuPAG 18. Further, a large number of variants were found for most of these isoforms. Phylogenetic relationship of identified BuPAGs showed that BuPAG 2 belonged to an ancient group while other members clustered with modern group. Three-dimensional (3D) structure of BuPAG 7 was determined by homology modeling and molecular dynamic (MD) simulations which displayed a typical fold represented by other aspartic proteinase (AP) family members. Molecular docking of Pepstatin inhibitor with BuPAG 7 revealed to interact through various hydrogen bonding and hydrophobic interactions. Various amino acid substitutions were observed in peptide-binding cleft of BuPAG 7. Superimposition of BuPAG 7 with homologous structures revealed the presence of a 35-41 amino acid long insertion (alpha helix connected by two loops) near the N- terminus which seems to be a unique feature of BuPAG 7 in AP family. This is the first report on identification and sequence characterization of PAG isoforms in buffalo with unique finding that these isoforms represent many transcript variants. We also report 3D structure of the most abundant isoform BuPAG 7 for the first time.