Jersey and Kashmiri cattle are important dairy breeds that contribute significantly to the total milk production of the Indian northern state of Jammu and Kashmir. The Kashmiri cattle germplasm has been extensively diluted through crossbreeding with Jersey cattle with the goal of enhancing its milk production ability. However, crossbred animals are prone to diseases resulting to unsustainable milk production. This study aimed to provide a comprehensive transcriptome profile of mammary gland epithelial cells at different stages of lactation and to find key differences in genes and pathways regulating milk traits between Jersey and Kashmiri cattle. Mammary epithelial cells (MEC) isolated from milk obtained from six lactating cows (three Jersey and three Kashmiri cattle) on day 15 (D15), D90 and D250 in milk, representing early, mid and late lactation, respectively were used. RNA isolated from MEC was subjected to next-generation RNA sequencing and bioinformatics processing. Casein and whey protein genes were found to be highly expressed throughout the lactation stages in both breeds. Largest differences in differentially expressed genes (DEG) were between D15 vs D90 (1,805 genes) in Kashmiri cattle and, D15 vs D250 (3,392 genes) in Jersey cattle. A total of 1,103, 1,356 and 1,397 genes were differentially expressed between Kashmiri and Jersey cattle on D15, D90 and D250, respectively. Antioxidant genes like RPLPO and RPS28 were highly expressed in Kashmiri cattle. Differentially expressed genes in both Kashmiri and Jersey were enriched for multicellular organismal process, receptor activity, catalytic activity, signal transducer activity, macromolecular complex and developmental process gene ontology terms. Whereas, biological regulation, endopeptidase activity and response to stimulus were enriched in Kashmiri cattle and, reproduction and immune system process were enriched in Jersey cattle. Most of the pathways responsible for regulation of milk production like JAK-STAT, p38 MAPK pathway, PI3 kinase pathway were enriched by DEG in Jersey cattle only. Although Kashmiri has poor milk production efficiency, the present study suggests possible physicochemical and antioxidant properties of Kashmiri cattle milk that needs to be further explored.

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/.

Pashmina goat (Capra hircus) is an economically important livestock species, which habitats the cold arid desert of the Ladakh region (India), and produces a princely animal fiber called Pashmina. The Pashmina goat has a double coat fleece as an adaptation to the very harsh cold winters the outer long coarse hair (guard hair) produced from primary hair follicles and the inner fine Pashmina fiber produced from secondary hair follicles. Pashmina fiber undergoes a circannual and synchronized growth cycle. In the present study, we analyzed transcriptome profiles from 10 different Pashmina goats during anagen and telogen to delineate genes and signaling pathways regulating active (anagen) and regressive (telogen) phases of the follicle growth. During anagen, 150 genes were expressed at significantly higher levels with log (FC) > 2 and p_adj < 0.05. The RNA seq results were subjected to qRT-PCR validation. Among the nine genes selected, the expression of HAS1, TRIB2, P2RX1. PRG4, CNR2, and MMP25 were significantly higher (p < 0.05) in the anagen phase, whereas MC4R, GIPC2, and CDO1 were significantly expressed (p < 0.05) in the telogen phase which supports and validates the gene expression pattern from the RNA-sequencing. Differentially expressed genes revealed that Pashmina fiber initiation is largely controlled by signaling pathways like Wnt, NF-Kappa, JAK-STAT, Hippo, MAPK, Calcium, and PI3K-Akt. Expression of genes from the Integrin family, Cell adhesion molecules, and ECM-receptors were observed to be at much higher levels during anagen. We identified key genes (IL36RN, IGF2, ITGAV, ITGA5, ITCCR7, CXCL5, C3, CCL19, and CXCR3) and a collagen cluster which might be tightly correlated with anagen-induction. The regulatory network suggests the potential role of RUNX3, NR2F1/2, and GATA family transcription factors in anagen-initiation and maintaining fiber quality in Pashmina goats.

Membrane proteins plays significant role in living cells. Transmembrane proteins are estimated to constitute approximately 30% of proteins at genomic scale. It has been a difficult task to develop specific alignment tools for transmembrane proteins due to limited number of experimentally validated protein structures. Alignment tools based on homology modeling provide fairly good result by recapitulating 70–80% residues in reference alignment provided all input sequences should have known template structures. However, homology modeling tools took substantial amount of time, thus aligning large numbers of sequences becomes computationally demanding. Here we present TM-Aligner, a new tool for transmembrane protein sequence alignment. TM-Aligner is based on Wu-Manber and dynamic string matching algorithm which has significantly improved its accuracy and speed of multiple sequence alignment. We compared TM-Aligner with prevailing other popular tools and performed benchmarking using three separate reference sets, BaliBASE3.0 reference set7 of alpha-helical transmembrane proteins, structure based alignment of transmembrane proteins from Pfam database and structure alignment from GPCRDB. Benchmarking against reference datasets indicated that TM-Aligner is more advanced method having least turnaround time with significant improvements over the most accurate methods such as PROMALS, MAFFT, TM-Coffee, Kalign, ClustalW, Muscle and PRALINE. TM-Aligner is freely available through http://lms.snu.edu.in/TM-Aligner/.

Changes in whole blood DNA methylation levels at several CpG sites have been associated with circulating blood lipids, specifically high-density lipoprotein and triglycerides. This study performs a discovery and validation epigenome-wide association study (EWAS) for circulating lipoprotein(a) [Lp(a)], an independent risk factor for cardiovascular diseases. Whole-blood DNA methylation profiles were assessed in a cohort of 1020 elderly individuals using the Illumina EPIC array and independent validation in 359 elderly males using the Illumina 450 k array. Plasma Lp(a) was measured using an apolipoprotein(a)-size-independent ELISA. Epigenome-wide rank regression analysis identified and validated a single CpG site, cg17028067 located in intron 1 of the LPA gene, that was significantly associated with plasma Lp(a) levels after correction for multiple testing. Genotyping of the site identified a relatively uncommon SNP (rs76735376, MAF <0.02) at the CpG site that largely explained the observed methylation effect. Rs76735376 is an expression quantitative trait loci for the LPA gene and could affect expression by altering enhancer activity. This EWAS for plasma Lp(a) identified a single CpG site within LPA. This association is due to an uncommon, but highly effective genetic variant, which was not in significant linkage disequilibrium with other variants known to influence Lp(a) levels or apo(a) isoform size. This study highlights the utility of CpG site methylation to identify potentially important genetic associations that would not be readily apparent in a comparable size genetic association study.