Whole genome resequencing revealed genomic variants and functional pathways related to adaptation in Indian yak populations

The present study aims to identify genomic variants through a whole genome sequencing (WGS) approach and uncover biological pathways associated with adaptation and fitness in Indian yak populations. A total of 30 samples (10 from each population) were included from Arunachali, Himachali and Ladakhi yak populations. WGS analysis revealed a total of 32171644, 27260825, and 32632460 SNPs and 4865254, 4429941, and 4847513 Indels in the Arunachali, Himachali, and Ladakhi yaks, respectively. Genes such as RYR2, SYNE2, BOLA, HF1, and the novel transcript ENSBGRG00000011079 were found to have the maximum number of high impact variants in all three yak populations, and might play a major role in local adaptation. Functional enrichment analysis of genes harboring high impact SNPs revealed overrepresented pathways related to response to stress, immune system regulation, and high-altitude adaptation. This study provides comprehensive information about genomic variants and their annotation in Indian yak populations, thus would serve as a data resource for researchers working on the yaks. Furthermore, it could be well exploited for better yak conservation strategies by estimating population genetics parameters viz., effective population size, inbreeding, and observed and expected heterozygosity.

Molecular insights into Pashmina fiber production: comparative skin transcriptomic analysis of Changthangi goats and sheep

Ladakh, one of the highest inhabited regions globally, hosts the unique Changthangi goat, renowned for producing Pashmina, the world's most luxurious natural fiber. In comparison, the fiber derived from Changthangi sheep is considered next only to Pashmina. This research endeavors to compare the skin transcriptome profiles of Changthangi goats and Changthangi sheep, aiming to discern the molecular determinants behind the recognition of Changthangi goats as the source of Pashmina. Drawing upon previously conducted studies, a collective of 225 genes correlated with fiber characteristics were extracted from the differentially expressed genes noticed between the two species (p-value of ≤ 0.05 and a log2 fold change of ≥ 1.5). These genes were analyzed using DAVID software to understand their biological functions and to identify enriched KEGG and Reactome pathways. The protein-protein interaction networks were constructed using Cytoscape, cytoHubba, and STRING to focus on key genes and infer their biological significance. Comparative transcriptome analysis revealed significantly higher expression of genes involved in signaling pathways like Wnt, MAPK, PI3K-Akt, Hedgehog, associated with fiber development and quality in Changthangi goats. These pathways play crucial roles in hair follicle (HF) formation, maintenance of epidermal stem cells, and fiber characteristics. Findings also highlight the enrichment of cell adhesion molecules and ECM-receptor interaction, emphasizing their roles in HF structure, growth, and signaling. This investigation offers an in-depth understanding of the molecular intricacies governing Pashmina production in Changthangi goats, providing valuable insights into their unique genetic makeup and underlying mechanisms influencing the exceptional quality of Pashmina fibers.

Identification and validation of stable reference genes for expression profiling of target genes in diverse ovine tissues

Quantitative PCR (qPCR) is a widely-used technique for quantifying the expression of target genes across various tissues, as well as under different pathological and physiological conditions. One of the challenges associated with this method is the need to identify optimal reference genes (RGs) that maintain consistent expression levels under diverse experimental settings, thereby ensuring accurate biological interpretation. In this study, we conducted a thorough analysis of 18 candidate RGs (ACTB, BACH1, B2M, GAPDH, HMBS, HPRT1, PGK1, PPIA, PPIB, RPLP0, RPL19, RPS9, RPS15, RPS28, SDHA, TBP, UXT, and YWHAZ) across 10 ovine tissues (muscle, skin, kidney, liver, intestine, rumen, lung, testis, heart, and spleen) obtained from five individual sheep. We aimed to identify genes with stable expression across these tissues. A literature-based survey helped us shortlist candidate genes representing various functional classes from multiple livestock species. We employed four algorithms: geNorm, NormFinder, BestKeeper, and Delta Ct (ΔCt), to rank these genes based on their stability. A consistent trend in the rankings was observed across these different algorithms. RefFinder was then used for a comprehensive ranking, integrating the outputs from the various methods. ACTB, PPIB, BACH1, and B2M emerged as the most stable RGs, while RPS9, RPS15, and PGK1 displayed variable expression. We validated our findings through qPCR analysis of four target genes (ACTN2, CRYAB, DLK1, and TRIM54) in the skin samples from two different sheep breeds. Based on these results, we recommend ACTB, PPIB, BACH1, and B2M as reliable internal control genes for qPCR experiments involving diverse ovine tissues.

Comprehensive evaluation and validation of optimal reference genes for normalization of qPCR data in different caprine tissues

BACKGROUND

Quantitative real-time PCR (qPCR) is a highly reliable method for validating gene expression data in molecular studies due to its sensitivity, specificity, and efficiency. To ensure accurate qPCR results, it's essential to normalize the expression data using stable reference genes.

METHODS

This study aimed to identify suitable reference genes for qPCR studies in goats by evaluating 18 candidate reference genes (ACTB, BACH1, B2M, GAPDH, HMBS, HPRT1, PGK1, PPIA, PPIB, RPLP0, RPL19, RPS9, RPS15, RPS28, SDHA, TBP, UXT, and YWHAZ) in 10 different caprine tissues (heart, intestine, kidney, liver, lung, muscle, rumen, skin, spleen, and testis). An integrated tool called RefFinder, which incorporates various algorithms like NormFinder, GeNorm, BestKeeper, and ΔCt, was used to assess the stability of expression among these genes.

RESULTS

After thorough analysis, ACTB, PPIB, and B2M emerged as the most stable reference genes, while RPL19, RPS15, and RPS9 were found to be the least stable. The suitability of the selected internal control genes was further validated through target gene analysis, confirming their efficacy in ensuring accurate gene expression profiling in goats.

CONCLUSION

The study determined that the geometric average of ACTB, PPIB, and B2M creates an appropriate normalization factor for gene expression studies in goat tissues.

Genetic insights into fiber quality, coat color and adaptation in Changthangi and Muzzafarnagri sheep: A comparative skin transcriptome analysis

Changthangi sheep, which inhabit the high-altitude regions of Ladakh, are known for their fine fiber production and are characterized by grey skin and either black or white coats. In contrast, Muzzafarnagri sheep from the plains of Uttar Pradesh produce coarse wool and have white skin and coats. We conducted comparative global gene expression profiling on four biological replicates of skin from each breed. Notably, our analysis identified 149 up-regulated genes and 2,139 down-regulated genes in Changthangi sheep compared to Muzzafarnagri sheep, with a p-adjusted value (padj) of ≤0.05 and a Log2 fold change of ≥1.5. Gene Ontology analysis of the up-regulated genes revealed an enrichment of terms related to melanin biosynthesis and developmental pigmentation. Additionally, enriched KEGG pathways included tyrosine metabolism and metabolic pathways. Among the melanogenesis-related genes that exhibited higher expression in Changthangi sheep were TYR, TYRP1, DCT, SLC45A2, PMEL, MLANA, and OCA2. These findings confirm melanin's role in both the animals' black coat color and UV protection at high-altitude. Furthermore, we observed more pronounced expression of genes related to fiber quality, namely KRTAP6, KRTAP7, KRTAP13, and KRTAP2, in the fine wool-producing sheep from Ladakh. The results of the RNA sequencing were validated using real-time PCR on 10 genes governing fiber quality and coat color, with ACTB and PPIB serving as reference genes. In conclusion, our comparative skin transcriptome analysis of Changthangi and Muzzafarnagri sheep sheds light on the genetic differences associated with distinct phenotypic traits and environmental adaptability, offering valuable insights into the underlying mechanisms.

Morphometric characteristics and microsatellite markers based diversity and differentiation recognizes the first prospective cattle breed from the Jharkhand state of India

India is bestowed with immense cattle biodiversity with 50 registered breeds. However, the majority (59.3%) is yet not characterized. Identification and characterization are the gateways to the management of prized indigenous resources. Present research described a unique cattle population of Jharkhand state, managed under a traditional low-input, low-output system. It was characterized by morphological traits, performance parameters, and management practices. Animals have the characteristic pre-scapular location of the hump. Genetic variation within this population and its differentiation with the six closely distributed cattle breeds were evaluated using FAO recommended microsatellite markers. Jharkhandi cattle have substantial genetic variation based on gene diversity (>0.6) and the average number of alleles per locus (>8). The population did not suffer from a genetic bottleneck in the recent past. Pairwise Nei's genetic distance, phylogenetic relationship, population differentiation, and the correct assignment of all the animals to self group substantiated its separate genetic identity. Since gene flow (Nm = 2.8-7.32) was identified and admixture was indicated by the Bayesian analysis there is a pressing need for scientific management of this population. Results endow authorities with critical information for registering a new Indian cattle breed (Medini) that contributes to the food security, livelihood, and economic sustainability of rural tribal households.

Tracing the genetic footprints: India's role as a gateway for pig migration and domestication across continents

This study explored the maternal genetic diversity in the pig genetic resources of India by analyzing a mitochondrial D-loop fragment and comparing it with the corresponding sequences of previously published studies involving domestic pigs and wild boars. Sequencing of 103 samples representing different domestic pig populations revealed existence of 32 maternal haplotypes. The indices of haplotype and nucleotide diversity in Indian domestic pigs were 0.9421 and 0.015, respectively. Median-Joining network revealed that Indian pigs belong to Clade A and show conformity to 6 haplogroups reported worldwide (D1a, D1a1, D1a2, D1e, D1h and D3a). Among these, D1e and D1a2 were shared with Asian wild boars too. Interestingly, haplotype sharing was evident between Indian pigs and samples from other countries representing Africa, Asia, Europe and Oceania. This study substantiates India's contribution as a possible pig domestication center and highlights the importance of the Indian subcontinent in dispersal of the species to other continents. Additionally, genetic evidence suggested the influence of trading routes and historical interactions in shaping pig genetic exchange. Overall, this investigation provides valuable insights into the genetic diversity, historical migration, and domestication of Indian domestic pigs, contributing to the broader understanding of global pig genetic resources and their evolutionary history.

Genetic diversity and differentiation of Thutho cattle from northeast India using microsatellite markers

Cattle are losing maximum breeds among the world's livestock. Genetic variability data is essentially required for conservation decision-making. Thutho is a recently registered Indian cattle breed (INDIA_CATTLE_1400_THUTHO_03047) from the northeast region (NE), a biodiversity hotspot. Genetic diversity in the Thutho population and its differentiation from the only other cattle breed of NE (Siri) and cattle (Bachaur) of the neighboring region was established using highly polymorphic, FAO-recommended microsatellite markers. Numerous alleles (253) were detected across the 25 loci. The mean observed and expected numbers of alleles in the population were 10.12 ± 0.5 and 4.5 ± 0.37, respectively. The observed heterozygosity (0.67 ± 0.04) was lower than the expected heterozygosity (0.73 ± 0.03) which indicated a departure from the Hardy-Weinberg equilibrium. A positive FIS value (0.097) confirmed the heterozygote deficiency in the Thutho population. Genetic distance, phylogenetic relationships, differentiation parameters, population assignment, and Bayesian analysis explicitly ascertained the unique genetic identity of the Thutho cattle. The population did not suffer any bottlenecks in the past. Thutho has minimum diversity among the three populations; hence, its scientific management needs to be initiated immediately. Interestingly, genetic variation is enough for formulating breeding programs for managing, improving, and conserving this precious indigenous cattle germplasm.

Expression profiling of cytokine genes in peripheral blood mononuclear cells from Anaplasma marginale infected and healthy cattle

In this study, changes in expression profiles of genes encoding 14 cytokines (IL1A, IL1B, IL2, IL4, IL6, IL8, IL10, IL12A, IL12B, IL16, IFNA, IFNB, TGFB1, and TNFA) were investigated amongst six Anaplasma marginale infected and six healthy crossbred cattle. Health status of the animals was determined based on clinical signs, blood smear examination and molecular detection using A. marginale-specific primers. Total RNA was isolated from the peripheral blood mononuclear cells of the infected animals as well as the healthy controls, which was further reverse transcribed to cDNA. Primers for real time PCR were designed using Primer3 software and the results were analyzed by the 2-ΔΔCt method with RPS15 and GAPDH as the reference genes. The expression levels of IL1A, IL1B, IL6, IL10, IL12A, IL12B, and TNFA varied significantly between the two groups, with higher expression in the infected cattle. The transcript abundance of IL4, IL16, and TGFB1 did not vary between the diseased and healthy animals. The expression of IL2 and IL8 was higher in the healthy animals, but the results were non-significant. Taken together, this study provides evidence for difference in expression of cytokine genes in response to anaplasmosis in crossbred cattle.

Assessment of genetic diversity of the fat-tailed Dumba sheep of India by mitochondrial and microsatellite markers

India has a centuries-old tradition of sheep production and breeding that accomplish economic, agricultural, and religious roles. In addition to the 44 registered sheep breeds, there is a fat-tailed sheep population referred to as Dumba. This study evaluated genetic variation in Dumba sheep and its differentiation from other Indian sheep breeds using mitochondrial DNA and genomic microsatellite loci. Haplotype and nucleotide diversity based on mitochondrial DNA analysis revealed substantially high maternal genetic diversity in Dumba sheep. Major ovine haplogroups A and B observed in sheep populations across the globe registered their presence in the Dumba sheep. The molecular genetic analysis using microsatellite markers also showed high measures of allele (10.125 ± 0.762) and gene diversity (0.749 ± 0.029). Results correspond to the non-bottleneck population that is near mutation-drift equilibrium despite some deficiency in the number of heterozygotes (FIS = 0.043 ± 0.059). Phylogenetic clustering confirmed Dumba to be a distinct population. Results of this study endow authorities with critical information imperative for sustainable utilization and conservation of Indian fat-tailed sheep, which is considered to be an untapped genetic resource contributing to the food security, livelihood, and economic sustainability of rural households in marginal areas of the country.

Revelation of genes associated with energy generating metabolic pathways in the fighter type Aseel chicken of India through skeletal muscle transcriptome sequencing

In this study, comparative analysis of skeletal muscle transcriptome was carried out for four biological replicates of Aseel, a fighter type breed and Punjab Brown, a meat type breed of India. The profusely expressed genes in both breeds were related to muscle contraction and motor activity. Differential expression analysis identified 961 up-regulated and 979 down-regulated genes in Aseel at a threshold of log2 fold change ≥ ±2.0 (padj<0.05). Significantly enriched KEGG pathways in Aseel included metabolic pathways and oxidative phosphorylation, with higher expression of genes associated with fatty acid beta-oxidation, formation of ATP by chemiosmotic coupling, response to oxidative stress, and muscle contraction. The highly connected hub genes identified through gene network analysis in the Aseel gamecocks were HNF4A, APOA2, APOB, APOC3, AMBP, and ACOT13, which are primarily associated with energy generating metabolic pathways. The up-regulated genes in Punjab Brown chicken were found to be related to muscle growth and differentiation. There was enrichment of pathways such as focal adhesion, insulin signaling pathway and ECM receptor interaction in these birds. The results presented in this study help to improve our understanding of the molecular mechanisms associated with fighting ability and muscle growth in Aseel and Punjab Brown chicken, respectively.

Design and validation of high-density SNP array of goats and population stratification of Indian goat breeds

Goats are the supporting pillars of rural economy contributing significantly to meat and milk production in India. It is a species targeted for fulfilling the interdependent goals of poverty reduction and creation of employment for supporting the rural income. The increased demand for goat products necessitates their genetic characterization and improvement to augment the production of native breeds. Bi-allelic, genome wide, densely placed single nucleotide polymorphism (SNP) markers are most suitable for this purpose. This paper describes the design and validation of an Affymetrix Axiom-based high-density (HD) SNP chip for goats. The array was designed using a panel of 225 samples from 15 diverse goat breeds of India. In total, more than 38 million high quality SNPs were subjected to stringent filtering and 626,975 SNPs were finally tiled on the array. The average coverage of SNPs in our chip is one SNP per four kilobase (kb), providing a denser coverage of the goat genome than previously available arrays. The HD chip (Axiom_Cahi) was validated by genotyping 443 samples from 26 indigenous goat breeds/populations. The results revealed 95.83% markers to be highly informative and polymorphic in Indian goats. Multivariate analysis indicated population structuring, as 15 breeds could be segregated using the designed array. Phylogenetic analysis suggested stratification of breeds by geographic proximity. This HD SNP chip for goats is a valuable resource for genomic selection, genome wide association as well as population genetic studies in goats.

Exploring the Transcriptome Dynamics of In Vivo Theileria annulata Infection in Crossbred Cattle

The molecular changes occurring in the host in response to in vivo Theileria annulata parasitic infection are not well understood. Therefore, the present study investigated the differential expression profiles of peripheral blood mononuclear cells (PBMCs) across Theileria annulata-infected and non-infected crossbred cows. The differential expression profiles from PBMCs of infected and non-infected crossbred cows were generated by RNA sequencing. A marked difference in the expression of genes associated with innate immunity (FTH1, ACTB, ISG15) was observed between the two groups. The over-represented pathways in Theileria annulata-infected cows were associated with the immune system and regulation of the mitotic cycle. Enriched genes and pathways in non-infected animals were associated with the maintenance of chromatin integrity and cell structure. The highly connected genes identified in this study form potential candidates for further investigation into host-parasite interactions in cattle. An improved understanding of the transcriptomic dynamics during theileriosis would lead to underpinning molecular level differences related to the health status of cattle.

Mitogenomic phylogeny reveals the predominance of the Nubian lineage of African wild ass in Indian donkeys

To contribute to the knowledge of maternal genetic diversity in domestic donkeys, this study investigated the mitochondrial DNA variations and analyzed the genetic structure in Indian donkeys based on 31 mitogenome sequences representing four breeds/populations (Agra, Halari, Kachchhi and Spiti). A total of 27 haplotypes with a haplotype diversity value of 0.989 were evident in the donkey genetic resources of India. The genetic differentiation between the investigated populations was evaluated using population pairwise F_ST values, which showed maximum differentiation between Kachchhi and Halari donkeys. The Neighbor-Joining (NJ) tree based on the whole mitogenome sequence and the Median-Joining (MJ) network for partial D-loop fragment showed clear demarcation of Indian donkeys into Nubian and Somali clades, substantiating African maternal origin of Indian domestic donkeys. The topology of the MJ network excluded the Asian wild asses as the possible progenitors of Indian donkeys. Halari and Agra donkeys showed conformity exclusively to the Nubian lineage of the African wild asses. However, representation of both the Nubian and Somali lineages was observed in Kachchhi and Spiti donkeys. Comprehensive analysis carried out by retrieving D-loop sequences from different countries representing Asia, Africa, Europe and South America revealed existence of shared haplotypes across geographically isolated regions of the globe. This observation is indicative of utility of donkeys as pack animals across inter-continental trading routes during development of human civilizations. Our results represent a valuable contribution to maternal genetic diversity of Indian donkeys and provide insights into the worldwide spread of the species following initial domestication in Africa.

Transcriptomics reveals key genes responsible for functional diversity in pectoralis major muscles of native black Kadaknath and broiler chicken

RNA sequencing-based expression profiles from pectoralis major muscles of black meat (Kadaknath) and white meat (broiler) chicken were compared to identify differentially expressed genes. A total of 156 genes with log2 fold change ≥  ± 2.0 showed higher expression in Kadaknath and 68 genes were expressed at a lower level in comparison to broiler. Significantly enriched biological functions of up-regulated genes in Kadaknath were skeletal muscle cell differentiation, regulation of response to reactive oxygen, positive regulation of fat cell differentiation and melanosome. Significant ontology terms up-regulated in broiler included DNA replication origin binding, G-protein coupled receptor signaling pathway and chemokine activity. Highly inter-connected differentially expressed genes in Kadaknath (ATFs, C/EPDs) were observed to be important regulators of cellular adaptive functions, while in broiler, the hub genes were involved in cell cycle progression and DNA replication. The study is an attempt to get an insight into the transcript diversity of pectoralis major muscles of Kadaknath and broiler chicken.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03682-0.

Functional Quality Characteristics of the Meat from a Dual-Purpose Poultry Crossbreed Suitable for Backyard Rearing in Comparison to Commercial Broilers

Backyard poultry farming contributes to food security, nutrition, and the regular income of rural farmers in India. Their products have a niche market here and fetch higher prices than those of commercial poultry. Improved varieties are being developed to overcome the slow growth, late sexual maturity, and low production of indigenous breeds, while retaining their positive attributes. A comprehensive study was conducted to analyze the functional attributes of meat from the Jabalpur color (JBC), a colored, improved dual-purpose synthetic line, developed by Nanaji Deshmukh Veterinary Science University, Jabalpur, India. The birds were managed in a deep litter system under a backyard type of housing (night shelter and free range). Primal meat cuts (breast and thigh) of the male birds (n = 20/group) were evaluated at the age of marketing. The corresponding attributes were compared with the results obtained for commercial Cobb (400) broilers. The protein concentration of JBC breast (25.65 ± 0.39 g/100 g of tissue) and thigh (19.04 ± 0.23 g/100 g of tissue) meat was superior (p ≤ 0.05) to that of Cobb broilers. Established assays (in vitro) identified a better (p ≤ 0.05) antioxidation capacity in the JBC meat. High-performance liquid chromatography confirmed a considerable quantity of functional biomolecules (carnosine, anserine, and creatine) in the JBC breast and thigh meat extracts. The average carnosine concentration (mg/g of tissue) was 2.66 ± 0.09 and 1.11 ± 0.04 in the JBC breast and thigh meat, respectively. The mRNA expression was quantified by qRT-PCR for the carnosine-related genes: β-alanine transporter (SLC36A1), carnosine-synthesizing enzyme (CARNS1), and carnosine-degrading enzyme (CNDP2); this explained the comparable carnosine in the JBC and Cobb meat. Meat extracts from both genetic groups (JBC and Cobb) had high anti-glycation potential. Higher protein content and antioxidant capacity, along with the bioactive dipeptides in the JBC meat, herald exciting research opportunities for its use in improving the traditional backyard poultry farming system.

Unraveling the genetic mechanisms governing the host response to bovine anaplasmosis

Bovine anaplasmosis caused by Anaplasma marginale is a tick-borne disease of livestock with widespread prevalence and huge economic implications. In order to get new insights into modulation of host gene expression in response to natural infections of anaplasmosis, this study is the first attempt that compared the transcriptome profiles of peripheral blood mononuclear cells (PBMCs) of A. marginale infected and healthy crossbred cattle. Transcriptome analysis identified shared as well as unique functional pathways in the two groups. Translation and structural constituent of ribosome were the important terms for the genes abundantly expressed in the infected as well as healthy animals. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the differentially expressed genes revealed that immunity and signal transduction related terms were enriched for the up-regulated genes in the infected animals. The over-represented pathways were cytokine-cytokine receptor interaction and signaling pathways involving chemokines, Interleukin 17 (IL17), Tumour Necrosis Factor (TNF), Nuclear Factor Kappa B (NFKB) etc. Interestingly, many genes previously associated with parasite-borne diseases such as amoebiasis, trypanosomiasis, toxoplasmosis, and leishmaniasis were profusely expressed in the dataset of the diseased animals. High expression was also evident for the genes for acute phase response proteins, anti-microbial peptides and many inflammatory cytokines. Role of cytokines in mediating communication between immune cells was the most conspicuous gene network identified through the Ingenuity Pathway Analysis. This study provides comprehensive information about the crosstalk of genes involved in host defense as well as parasite persistence in the host upon infection with A. marginale.

Validation of stable reference genes in peripheral blood mononuclear cells for expression studies involving vector-borne haemoparasitic diseases in bovines

Normalization of gene expression data using appropriate reference genes is critical to diminish any technical bias in an experiment involving quantitative real-time PCR (qPCR). To the best of our knowledge, this is the first report offering a systematic assessment of 14 potential reference genes (RPLP0, ACTB, RPS28, YWHAZ, SDHA, PPIA, RPS9, RPS15, UXT, GAPDH, B2M, BACH1, HMBS, and PPIB) for the identification of the most stable normalizers for qPCR of target genes in peripheral blood mononuclear cells (PBMCs) of bovines for vector-borne haemoparasitic diseases such as anaplasmosis, babesiosis, theileriosis, and trypanosomiasis. A total of 38 blood samples were collected from healthy as well as diseased cattle and buffaloes representing different haemoparasitic diseases. RNA isolated from the PBMCs was subjected to qPCR for the 14 prospective internal control genes. The comprehensive ranking of the genes was accomplished by the RefFinder tool that integrates the results of three algorithms (geNorm, NormFinder, and BestKeeper) and the comparative C_T method. RPS15, B2M, and GAPDH were ranked to be the most stable genes, whereas, PPIA and HMBS emerged to be the least suitable genes. Validation of the selected reference genes by the qPCR analysis of two immunity genes, ISG15 and GPX7 was congruent with the observations of this study. We recommend that a panel of three reference genes including RPS15, B2M, and GAPDH could prove useful in delineating the transcriptional landscape of PBMCs for vector-borne haemoparasitic diseases in bovines.

First report on delineation of differentially expressed genes and pathways in milk somatic cells of mastitic and healthy Murrah buffaloes

Despite immense contribution of buffaloes as dairy species, limited studies have addressed the bubaline mastitis as compared to cattle. This was the first differential transcriptomic study investigating the alterations induced by clinical mastitis in buffalo milk relative to healthy controls. Comparative gene expression profiling of three biological replicates of each group identified 1014 up-regulated and 999 down-regulated genes in the diseased buffaloes (Fold change > 2, FDR < 0.05). Activation of immune and inflammatory responses were the most enriched GO terms in the mastitic animals, with higher transcript abundance of many genes coding for anti-microbial proteins such as β-defensins, perforin, granzymes, granulysin, cathelicidins etc. Analysis of the gene regulatory interactions of the up-regulated DEGs identified many hub genes that govern the cellular and macromolecular metabolic processes (E2F4, E2F1, RBL2, FOXM1, IRF1 and MYB). This study contributes to an insightful understanding of molecular mechanisms governing immune response of buffaloes to mastitis.

Mitochondrial DNA diversity divulges high levels of haplotype diversity and lack of genetic structure in the Indian camels

Camels represent an important genetic resource of the desert ecosystems of India, with the dromedary and Bactrian camels inhabiting the hot and cold deserts, respectively. This study is the first attempt to investigate mitochondrial DNA based genetic diversity in the Indian camel populations and explores their relationship in the context of global genetic diversity of all the three large camel species (Camelus ferus, Camelus bactrianus and Camelus dromedaries). A mitochondrial DNA fragment encompassing part of cytochrome b gene, tRNA^Thr, tRNA^Pro and the beginning of the control region was amplified and analyzed in 72 dromedary and 8 Bactrian camels of India. Sequence analysis revealed that the haplotype and nucleotide diversity (Hd: 0.937 and π: 0.00431) in the Indian dromedaries was higher than the indices reported so far for the dromedary or Bactrian camels across the globe. The corresponding values in the Indian Bactrian camels were 1.000 and 0.00393, respectively. Signals of population expansion were evident in the dromedaries of India on the basis of mismatch analysis and Fu's Fs values. The analysis of molecular variance attributed most of the genetic variance (92.15%) between the dromedary, wild Bactrian and domestic Bactrian camels indicating separate maternal origins. The existence of three mitochondrial lineages in the old world camels (C. bactrianus: Lineage A; C. ferus: Lineage B and C. dromedarius: Lineage C) was also substantiated by the topology of the Median-Joining network.

Mining of diverse short non-coding RNAs from transcriptome of milk somatic cells of Murrah buffalo

The non-coding RNAs (ncRNA) are known to regulate expression of genes at the transcription, translation and processing levels. The present study was conducted to identify diverse short ncRNAs from milk somatic cells of lactating Murrah buffaloes. Elucidating the molecular drivers of lactation in dairy animals will help understand the process of lactation, eventually leading to improvement in milk production and quality. In order to discover the ncRNA, the transcriptome data of 12 samples of somatic cells from buffalo milk were analyzed. A web based pipeline, exceRpt was used to perform the analysis. The most abundant short ncRNA molecules discovered in buffalo milk were the miRNAs, followed by snRNAs. Least number of rRNAs was discovered in the investigated samples. The total number of rRNAs, tRNAs, snRNAs, snoRNAs and miRNAs were 12, 23, 72, 51 and 229 respectively, in the entire dataset. On matching with miRBase v22.1, a total of 1724, 897, 211 and 4 miRNAs were observed to be common to human, bovine, caprine and ovine genomes. The results provide information on the bioavailability of short ncRNAs in buffalo milk somatic cells, most of which are largely uncharacterized. The generated information is a step towards developing a database for ncRNAs in buffalo species.

Muscle transcriptome provides the first insight into the dynamics of gene expression with progression of age in sheep

The dynamic synergy of genes and pathways in muscles in relation to age affects the muscle characteristics. Investigating the temporal changes in gene expression will help illustrate the molecular mechanisms underlying muscle development. Here we report the gene expression changes in skeletal muscles through successive age groups in Bandur, a meat type sheep of India. RNA sequencing data was generated from the longissimus thoracis muscles from four age groups, ranging from lamb to adult. Analysis of 20 highest expressed genes common across the groups revealed muscle protein, phosphorylation, acetylation, metal binding and transport as significant functions. Maximum differentiation was observed after 2.5–3 years on transition from lambs to adult. Transcriptional regulation by the TFAP2 transcription factors, IL-6 signaling and PI3K/AKT signaling pathways were enriched in younger animals. The gene-protein network demarcated key interactive genes involved in muscle development and proliferation that can be used as candidates for future research on improvement of muscle characteristics.

Genetic differentiation of Indian dromedary and Bactrian camel populations based on mitochondrial ATP8 and ATP6 genes

Camelids are acknowledged worldwide to endure hostile conditions prevalent in the hot as well cold deserts across the globe. Adaptations to climatic extremes have been associated with mitochondrial protein variants such as ATP8 and ATP6 in different species. The camel genetic resources of India are represented by 9 breeds of dromedary camels which inhabit hot arid and semi-arid zones of the country and a small population of Bactrian camels found in the cold desert of Ladakh. In this study, within and between breed genetic diversity in Indian dromedaries and their divergence from Bactrian camels was investigated based on ATP8/6 genes. Sequence analysis of a mitochondrial DNA fragment encompassing ATP8 and ATP6 genes identified 15 haplotypes in the dromedaries of India and 3 haplotypes in Bactrian camels. The values of haplotype diversity and nucleotide diversity were 0.647 and 0.00187 in the former and 0.679 and 0.00098, respectively in the latter. AMOVA analysis revealed 97.81% variance between the two species. Median-Joining network delineated three distinct mitochondrial haplogroups for Camelus dromedarius, Camelus ferus and Camelus bactrianus. Clear demarcation of the old world (Dromedary and Bactrian camels) and new world camelids (Alpaca, llama, guanaco and vicugna) was evident through the phylogenetic analysis.

Transcriptomic diversity in longissimus thoracis muscles of Barbari and Changthangi goat breeds of India

The present study is an attempt to examine the differential expression of genes in longissimus thoracis muscles between meat and wool type Indian goat breeds. Barbari goat is considered the best meat breed while Changthangi is famous for its fine fibre quality. RNA sequencing data was generated from four biological replicates of longissimus thoracis muscles of Barbari and Changthangi goats. A clear demarcation could be observed between the breeds in terms of expression of genes associated with lipid metabolism (FASN, SCD, THRSP, DGAT2 and FABP3). Most significant genes with high connectivity identified by gene co-expression network analysis were associated with triacylglycerol biosynthesis pathway in Barbari goat. Highly interactive genes identified in Changthangi goat were mainly associated with muscle fibre type. This study provides an insight into the differential expression of genes in longissimus thoracis muscles between Barbari and Changthangi goats that are adapted to and reared in different agro-climatic regions.

Microsatellite and mitochondrial DNA analyses unveil the genetic structure of native sheep breeds from three major agro-ecological regions of India

Sheep farming has been fundamental to many civilizations in the world and is practiced in India since antiquity. Several thousand years of adaptation to local environmental conditions and selective breeding have evolved 44 sheep breeds in India. They are paramount in terms of economic, scientific, and cultural heritage. Genetic characterization information is imperative for sustainable utilization and conservation of ovine heritage. In this study, the genetic diversity, differentiation, and structure of 11 indigenous sheep breeds from three different agro-ecological zones of India were explored with genomic microsatellite loci and mitochondrial DNA (D loop). The estimated diversity parameters indicated that populations retained high levels of genetic diversity (Na = 8.27 ± 0.17; Ho = 0.65 ± 0.01), which provides an optimistic viewpoint for their survival. However, significant inbreeding was also observed in the nine populations. Moderate genetic differentiation existed among the groups (F_ST = 0.129 ± 0.012), and most likely clusters existing in the dataset are seven. Phylogenetic clustering was in line with the geographical locations of sheep populations. Mitochondrial sequences revealed high haplotype diversity with the existence of maternal haplogroups A, B, and C, and signals of population expansion. Decreased genetic diversity and unique maternal lineage (C) in endangered Tibetan and Bonpala sheep breed, warrant their immediate scientific management. Overall, the quantitative data reported here on the extant variability, and genetic relationships among the Indian sheep breeds, provide critically important inputs that will be valuable for the decision-making process on their management, both for the conservation of endangered breeds, and formulation of breeding programs to check genetic erosion.

Differential expression of miRNAs in skeletal muscles of Indian sheep with diverse carcass and muscle traits

The study presents the miRNA profiles of two Indian sheep populations with divergent carcass and muscle traits. The RNA sequencing of longissimus thoracis muscles from the two populations revealed a total of 400 known miRNAs. Myomirs or miRNAs specific to skeletal muscles identified in our data included oar-miR-1, oar-miR-133b, oar-miR-206 and oar-miR-486. Comparison of the two populations led to identification of 100 differentially expressed miRNAs (p < 0.05). A total of 45 miRNAs exhibited a log₂ fold change of ≥ ( ±) 3.0. Gene Ontology analysis revealed cell proliferation, epithelial to mesenchymal transition, apoptosis, immune response and cell differentiation as the most significant functions of the differentially expressed miRNAs. The differential expression of some miRNAs was validated by qRT-PCR analysis. Enriched pathways included metabolism of proteins and lipids, PI3K-Akt, EGFR and cellular response to stress. The microRNA-gene interaction network revealed miR-21, miR-155, miR-143, miR-221 and miR-23a as the nodal miRNAs, with multiple targets. MicroRNA-21 formed the focal point of the network with 42 interactions. The hub miRNAs identified in our study form putative regulatory candidates for future research on meat quality traits in Indian sheep. Our results provide insight into the biological pathways and regulatory molecules implicated in muscling traits of sheep.

Comparative gene expression profiling of milk somatic cells of Sahiwal cattle and Murrah buffaloes

India is the world's largest milk producing country because of massive contribution made by cattle and buffaloes. In the present investigation, comprehensive comparative profiling of transcriptomic landscape of milk somatic cells of Sahiwal cattle and Murrah buffaloes was carried out. Genes with highest transcript abundance in both species were enriched for biological processes such as lactation, immune response, cellular oxidant detoxification and response to hormones. Analysis of differential expression identified 377 significantly up-regulated and 847 significantly down-regulated genes with fold change >1.5 in Murrah buffaloes as compared to Sahiwal cattle (p_adj <0.05). Marked enrichment of innate and adaptive immune response related GO terms and higher expression of genes for various host defense peptides such as lysozyme, defensin β and granzymes were evident in buffaloes. Genes related to ECM-receptor interaction, complement and coagulation cascades, cytokine-cytokine receptor interaction and keratinization pathway showed more abundant expression in cattle. Network analysis of the up-regulated genes delineated highly connected genes representing immunity and haematopoietic cell lineage (CBL, CD28, CD247, PECAM1 and ITGA4). For the down-regulated dataset, genes with highest interactions were KRT18, FGFR1, GPR183, ITGB3 and DKK3. Our results lend support to more robust immune mechanisms in buffaloes, possibly explaining lower susceptibility to mammary infections as compared to cattle.

Gut-organ axis: a microbial outreach and networking

Human gut microbiota (GM) includes a complex and dynamic population of microorganisms that are crucial for well-being and survival of the organism. It has been reported as diverse and relatively stable with shared core microbiota, including Bacteroidetes and Firmicutes as the major dominants. They are the key regulators of body homeostasis, involving both intestinal and extra-intestinal effects by influencing many physiological functions such as metabolism, maintenance of barrier homeostasis, inflammation and hematopoiesis. Any alteration in GM community structures not only trigger gut disorders but also influence other organs and cause associated diseases. In recent past, the GM has been defined as a 'vital organ' with its involvement with other organs; thus, establishing a link or a bi- or multidirectional communication axis between the organs via neural, endocrine, immune, humoral and metabolic pathways. Alterations in GM have been linked to several diseases known to humans; although the exact interaction mechanism between the gut and the organs is yet to be defined. In this review, the bidirectional relationship between the gut and the vital human organs was envisaged and discussed under several headings. Furthermore, several disease symptoms were also revisited to redefine the communication network between the gut microbes and the associated organs.

Identification of a new potential native Indian cattle breed by population differentiation based on microsatellite markers

India has a rich heritage of rearing cattle where farmers selected native cattle suitable to their local agro-ecological conditions for centuries. It is reflected in 50 indigenous breeds of cattle, besides many lesser known populations not explored so far. It is the need of the hour to characterize such populations to have prudent improvement and conservation options. Thus, present study was carried out to assess the genetic diversity and relationship between an unexplored local cattle population (Kathani) and four established cattle breeds of adjoining area (Gaolao, Kosali, Ongole and Motu) by using 20 FAO recommended microsatellite markers. High variability was recorded in the Kathani population with a total of 198 alleles that varied between 5 (ILSTS11, TGLA22, INRA05) and 17 (ILSTS34) with a mean of 9.9 ± 0.73. The average observed heterozygosity (Ho) was 0.658 ± 0.054. Heterozygote deficiency was not significant (F_IS = 0.029 ± 0.063) indicating random mating prevalent across this population. Mean estimates of observed number of alleles and heterozygosity over all the loci and five populations were 9.73 ± 0.421 and 0.617 ± 0.022, respectively. In the overall populations, the homozygote excess (F_IT) of 0.293 ± 0.032, was partly due to the homozygote excess within breeds (F_IS = 0.121 ± 0.025) and to a larger extent due to high (0.05 < F_ST < 0.15) genetic differentiation among them (F_ST = 0.195 ± 0.029). Substantial pairwise Nei's genetic distance and high population differentiation indicated towards separate genetic identity of Kathani cattle. The analysis of genetic structure based on Bayesian approach indicated that the most probable number of clusters is five confirming definitive genetic differentiation among all the popultions. Entire analysis showed that a significant amount of genetic variation is maintained in Kathani, a lesser known cattle population that is distinct from the recognized breeds in the proximity. As this autochthonous cattle plays role in the economic sustainability of a marginal and disadvantaged area, it is important to preserve and develop its breeding.

Comprehensive analysis of mitochondrial DNA based genetic diversity in Indian goats

Thirty four distinct breeds and many non-descript populations represent the caprine diversity of India. Genetic characterization of breeds is an essential element in designing breeding strategies and preserving genetic diversity. Considering the popularity of mitochondrial DNA for phylogeographical studies, this study involved an extensive analysis of population structure and genetic diversity of 28 defined breeds and 5 lesser known populations representing all four major agro-climatic zones of India using mitochondrial DNA markers. Analysis of hypervariable region 1 of mtDNA control region in 443 goats together with 22 reference sequences, delineated 341 distinct haplotypes belonging to four maternal haplogroups; A, B, C and D, with haplogroup A representing 90% of the individuals. The haplotype and nucleotide diversity indices of Indian goats were 0.998 ± 0.001 and 0.028 ± 0.001, respectively indicating abundant genetic variability. Estimates of population demographic parameters from mismatch analysis suggested a relatively good fit to the model of either spatial or demographic expansion of Indian goats. AMOVA analysis and topology of MJ network suggested lack of phylogeographic structure in domestic goats, which can be attributed to unstructured animal breeding, dwindling pastures and nomadic pastoralism. Genetic differentiation between goats from different agro-ecological regions was in accordance with their geographical propinquity.

A practical criterion for focusing of unstained cell samples using a digital holographic microscope

Digital holographic microscopy (DHM) is an important technique that may be used for quantitative phase imaging of unstained biological cell samples. Since the DHM technology is not commonly used in clinics or bioscience research labs, at present there is no well-accepted focusing criterion for unstained samples that users can follow while recording image plane digital holograms of cells. The usual sharpness metrics that are useful for auto-focusing of stained cells do not work well for unstained cells as there is no amplitude contrast. In this work, we report a practical method for estimating the best focus plane for unstained cells in the digital hologram domain. The method is based on an interesting observation that for the best focus plane the fringe pattern associated with individual unstained cells predominantly shows phase modulation effect in the form of bending of fringes and minimal amplitude modulation. This criterion when applied to unstained red blood cells shows that the central dip in the doughnut-like phase profile of cells is maximal in this plane. The proposed methodology is helpful for standardizing the usage of DHM technology across different users and application development efforts. LAY DESCRIPTION: Digital holographic microscopy (DHM) is slowly but steadily becoming an important microscopy modality and gaining acceptability for basic bio-science research as well as clinical usage. One of the important features of DHM is that it allows users to perform quantitative imaging of unstained transparent cells. Instead of using dyes or fluorescent labelling, DHM systems use quantitative phase as a contrast mechanism which depends on the natural refractive index variation within the cell samples. Since minimal wet lab processing is required in order to image cell samples with a DHM, cells can be imaged in their natural state. While DHM is gaining popularity among users, the imaging protocols across the labs or users need to be standardized in order to make sure that the same quantitative phase parameters are used for tasks such as quantitative phased based cell classification. One of the important operational tasks for any microscopy work is to focus the sample under study. While focusing comes naturally to users of brightfield microscopes based on image contrast, the focusing is not straightforward when samples are unstained so that they do not offer any amplitude contrast. When performing quantitative phase imaging, defocus can actually change the phase profile of the cell due to near-zone (Fresnel) diffraction effects. So unless a standardized focusing methodology is used, it will be difficult for multiple DHM users (potentially at different sites) to agree on quantitative results out of their phase images. DHM literature has prior works which perform numerical focusing of recovered complex wave-field in the hologram plane to find the best focus plane. However such methods are not user friendly and do not allow user the same focusing experience as in a brightfield microscope. The numerical focusing is therefore a reasonably good method for an optics researcher but not necessarily so for a microscopy technician looking at cell samples with a DHM system in a clinical setting. The present work provides a simple focusing criterion for unstained samples that works directly in the hologram domain. The technique is based on an interesting observation that the when an unstained cell sample is in the best-focus plane, its corresponding hologram (or fringe pattern) predominantly shows phase modulation manifested by bending of fringes at the location of the cell. This criterion can be converted into a simple numerical method as we have used to find the best-focus plane using a stack of through focus holograms. We believe that the technique can be used manually by visually observing the holograms or can be converted to an auto-focus algorithm for a motorized DHM system.

Identification of a new Indian camel germplasm by microsatellite markers based genetic diversity and population structure of three camel populations

Camel invokes fascinating chapter of Indian desert history and is integral component of its ecosystem. Camel population has reached a crisis point after three decades of decline (75%) causing major concern to the policy makers. >28% of Indian camel is not yet characterized. It is imperative to describe country’s camel germplasm and its existing diversity for designing conservation plan. One such population is Sindhi, distributed along border with Pakistan. Twenty five microsatellite markers being valuable tool for estimating genetic diversity were selected to elucidate genetic variability and relationship of Sindhi with two registered camel breeds of India- Marwari and Kharai. The standard metrics of genomic diversity detected moderate variability in all the three populations. A total of 303 alleles with a mean of 8.116 ± 0.587 alleles per locus were found in total of 143 animals. Sindhi population had intermediate allelic diversity with 8.522 ± 1.063 alleles per locus. Corresponding values in Marwari and Kharai were 8.783 ± 0.962 and 7.043 ± 1.030, respectively. Genetic variability within the breeds was moderate as evidenced by the mean observed heterozygosity of 0.556 ± 0.025. Sindhi camel population harbors higher genetic variability (Ho = 0.594) as compared to the two registered camel breeds (Marwari, 0.543 and Kharai, 0.531). Mean expected heterozygosity under Hardy-Weinberg equilibrium was higher than the observed values across the three camel groups, indicating deviations from assumptions of this model. In fact, average positive F value of 0.084 to 0.206 reflected heterozygote deficiency in these populations. These Indian camel populations have not experienced serious demographic bottlenecks in the recent past. Differences among populations were medium and accounted for 7.3% of total genetic variability. Distinctness of three camel populations was supported by all the approaches utilized to study genetic relationships such as genetic distances, phylogenetic relationship, correspondence analysis, clustering method based on Bayesian approach and individual assignment. Sindhi camel population was clearly separated from two registered breeds of Indian camel. Results conclude Sindhi to be a separate genepool. Moderate genetic diversity provides an optimistic viewpoint for the survival of severely declining indigenous camel populations with appropriate planning strategies for conserving the existing genetic variation and to avoid any escalation of inbreeding.

Skin transcriptome profiling of Changthangi goats highlights the relevance of genes involved in Pashmina production

Pashmina, the world's finest natural fiber is derived from secondary hair follicles of Changthangi goats which are domesticated in Ladakh region of Jammu and Kashmir by nomadic pastoralists. Complex epithelial-mesenchymal interactions involving numerous signal molecules and signaling pathways govern hair follicle morphogenesis and mitosis across different species. The present study involved transcriptome profiling of skin from fiber type Changthangi goats and meat type Barbari goats to unravel gene networks and metabolic pathways that might contribute to Pashmina development. In Changthangi goats, 525 genes were expressed at significantly higher levels and 54 at significantly lower levels with fold change >2 (p_adj < 0.05). Functional annotation and enrichment analysis identified significantly enriched pathways to be formation of the cornified envelope, keratinization and developmental biology. Expression of genes for keratins (KRTs) and keratin-associated proteins (KRTAPs) was observed to be much higher in Changthangi goats. A host of transcriptional regulator genes for hair follicle keratin synthesis such as GPRC5D, PADI3, HOXC13, FOXN1, LEF1 and ELF5 showed higher transcript abundance in Pashmina producing goats. Positive regulation of Wnt signaling pathway and negative regulation of Oncostatin M signaling pathway may be speculated to be important contributors to hair follicle development and hair shaft differentiation in Changthangi goats.

Transcriptome profiling of longissimus thoracis muscles identifies highly connected differentially expressed genes in meat type sheep of India

This study describes the muscle transcriptome profile of Bandur breed, a consumer favoured, meat type sheep of India. The transcriptome was compared to the less desirable, unregistered local sheep population, in order to understand the molecular factors related to muscle traits in Indian sheep breeds. Bandur sheep have tender muscles and higher backfat thickness than local sheep. The longissimus thoracis transcriptome profiles of Bandur and local sheep were obtained using RNA sequencing (RNA Seq). The animals were male, non-castrated, with uniform age and reared under similar environment, as well as management conditions. We could identify 568 significantly up-regulated and 538 significantly down-regulated genes in Bandur sheep (p≤0.05). Among these, 181 up-regulated and 142 down-regulated genes in Bandur sheep, with a fold change ≥1.5, were considered for further analysis. Significant Gene Ontology terms for the up-regulated dataset in Bandur sheep included transporter activity, substrate specific transmembrane, lipid and fatty acid binding. The down-regulated activities in Bandur sheep were mainly related to RNA degradation, regulation of ERK1 and ERK2 cascades and innate immune response. The MAPK signaling pathway, Adipocytokine signaling pathway and PPAR signaling pathway were enriched for Bandur sheep. The highly connected genes identified by network analysis were CNOT2, CNOT6, HSPB1, HSPA6, MAP3K14 and PPARD, which may be important regulators of energy metabolism, cellular stress and fatty acid metabolism in the skeletal muscles. These key genes affect the CCR4-NOT complex, PPAR and MAPK signaling pathways. The highly connected genes identified in this study, form interesting candidates for further research on muscle traits in Bandur sheep.

Buffalo milk transcriptome: A comparative analysis of early, mid and late lactation

The expression of genes and their regulation during lactation in buffaloes remains less understood. To understand the interplay of various genes and pathways, the milk transcriptome from three lactation stages of Murrah buffalo was analyzed by RNA sequencing. The filtered reads were mapped to the Bubalus bubalis as well as Bos taurus reference assemblies. The average mapping rate to water buffalo and Btau 4.6 reference sequence, was 75.5% and 75.7% respectively. Highly expressed genes (RPKM > 3000), throughout lactation included CSN2, CSN1S1, CSN3, LALBA, SPP1 and TPT1. A total of 12833 transcripts were common across all the stages, while 271, 205 and 418 were unique to early, mid and late lactation respectively. Majority of the genes throughout lactation were linked to biological functions like protein metabolism, transport and immune response. A discernible shift from metabolism in early stage to metabolism and immune response in mid stage, and an increase in immune response functions in late lactation was observed. The results provide information of candidate genes and pathways involved in the different stages of lactation in buffalo. The study also identified 14 differentially expressed and highly connected genes across the three lactation stages, which can be used as candidates for future research.

Evaluation of therapeutic potential of recombinant buffalo lactoferrin N-lobe expressed in E. coli

Lactoferrin (Lf) is a multifunctional bi-lobate iron-binding glycoprotein belonging to transferrin family with a mass of approximately 80 kD. Being ubiquitously present in almost all biological secretions, it performs important biological functions. One of the earliest and very well-documented functions of Lf is the antibacterial effect against broad spectrum Gram-negative and Gram-positive bacteria. In this study, buffalo Lf N-lobe cDNA was amplified, cloned and expressed as a fusion protein in Escherichia coli cells using pQE30 expression vector. After post-induction confirmation of expressed protein by SDS-PAGE, purification of recombinant protein using Ni-NTA was attempted and the yield of recombinant buffalo N-lobe Lf was estimated to be 1 mg/ml. Antibacterial activity of recombinant buffalo Lf N-lobe was assessed on pathogenic E. coli and Staphylococcus aureus strains. Peptic digest of recombinant N-lobe buffalo Lf showed antibacterial activity comparable to commercially available bovine Lf. The successful expression and characterization of functional recombinant N-lobe of buffalo Lf expressed in E. coli opens new vistas for developing alternate therapeutics, particularly against the diseases caused by Gram-negative microbes such as septicemia and diarrhea in newborn calves and mastitis in dairy animals.

Promoter methylation and expression analysis of Bvh gene in bulls with varying semen motility parameters

Crossbreds of low-producing indigenous cattle and high-producing exotic dairy bulls (Holstein Friesian and Jersey) have contributed in ensuring that India continues to be the world's top milk-producing country. However, subfertility observed in crossbred male progenies has been a major obstacle in exploitation of heterosis due to crossbreeding. There is sufficient scientific evidence in support of genetic and epigenetic regulation of key physiological processes including spermatogenesis. Bovine Vasa Homology (Bvh) is considered a molecular marker for the study of gametogenesis. Significant negative correlation between DNA methylation and gene expression has been reported in cattle-yaks hybrids and their parents. The present study analyzed promoter methylation status and expression profile of Bvh gene in spermatozoa from exotic Holstein Friesian cattle, indigenous Sahiwal cattle and their crossbreds with varying semen motility parameters. The degree of methylation of the Bvh promoter region was significantly higher in poor motility crossbred bulls (13.3%) as compared to good motility crossbreds (5.3%), Sahiwal (3%) and Holstein Friesian bulls (1%) (P < 0.05). Gene expression analysis revealed significantly higher mRNA abundance of Bvh in purebreds (Holstein Friesian and Sahiwal) as compared to crossbred counterparts (P < 0.001). Inverse correlation observed in this study between promoter methylation and gene expression of Bvh gene in spermatozoa from crossbred bulls with poor motility phenotype as compared to purebred parents provides an important insight into understanding the graded fertility of crossbred bulls.

Differential promoter methylation of DAZL gene in bulls with varying seminal parameters

In India, cross-breeding of indigenous cattle with exotic cattle such as Holstein Friesian and Jersey has been going on since last four decades to improve milk production. Although it has led to increased milk yield, the subfertility in male cross-bred progeny has remained a significant problem. Epigenetic modifications (DNA methylation, histone modifications and chromatin remodelling) are regarded as key players influencing gene expression. DAZL gene plays an important role in germline development and gametogenesis. The methylation and mRNA expression level of this gene have been significantly negatively correlated in the testes of cattle-yak hybrids and their parents. This study analysed the methylation profile of DAZL gene promoter in bull spermatozoa in an attempt to speculate its role in cross-bred cattle subfertility. Semen samples from Sahiwal, Holstein Friesian and Frieswal bulls (Sahiwal X Holstein Friesian) with varying semen motility parameters were collected, and DNA was isolated. Methylation-specific primers were used to amplify part of promoter and exon 1 of DAZL gene using bisulphite-converted DNA. The amplified products were sequenced after cloning in pTZ57R/T vector. Sequence analysis revealed significantly higher DNA methylation of DAZL gene in Frieswal bulls with poor motility (28.26%) as compared to medium (15.21%) and high motility phenotype (6.52%). In pure-bred counterparts, Sahiwal and Holstein Friesian, epigenetic marks were more in the former (15.21%) than the latter (4.34%), but in both cases, the values were lower as compared to the poor motility Frieswal bulls. This suggests that differential hypermethylation of the CpG islands could possibly influence reproductive parameters in bovines.

MDCT of hand and wrist infections: emphasis on compartmental anatomy

Hand and wrist infections can present with a spectrum of manifestations ranging from cellulitis to deep-space collections. The various infectious processes can be categorised as superficial or deep infections based on their respective locations relative to the tendons. Superficial hand infections are located superficial to the tendons and are comprised of cellulitis, lymphangitis, paronychia, pulp-space infections, herpetic whitlow, and include volar as well as dorsal subcutaneous abscesses. Deep hand infections are located deep to the tendon sheaths and include synovial space infections, such as infectious tenosynovitis, deep fascial space infections, septic arthritis, necrotising fasciitis, and osteomyelitis. Knowledge of hand and wrist compartmental anatomy is essential for the accurate diagnosis and management of hand infections. Although early and superficial infections of the hand may respond to non-surgical management, most hand infections are surgical emergencies. Multidetector computed tomography (MDCT), with its muliplanar reformation (MPR) and three-dimensional (3D) capabilities, is a powerful tool in the emergency setting for the evaluation of acute hand and wrist pathology. The clinical and imaging features of hand and wrist infections as evident on MDCT will be reviewed with emphasis on contiguous and closed synovial and deep fascial spaces. Knowledge of hand compartmental anatomy enables accurate characterisation of the infectious process and localise the extent of disease in the acute setting.

Evidence of positive selection and concerted evolution in the rapidly evolving PRDM9 zinc finger domain in goats and sheep

Meiotic recombination contributes to augmentation of genetic diversity, exclusion of deleterious alleles and proper segregation of chromatids. PRDM9 has been identified as the gene responsible for specifying the location of recombination hotspots during meiosis and is also the only known vertebrate gene associated with reproductive isolation between species. PRDM9 encodes a protein with a highly variable zinc finger (ZF) domain that varies between as well as within species. In the present study, the ZF domain of PRDM9 on chromosome 1 was characterized for the first time in 15 goat breeds and 25 sheep breeds of India. A remarkable variation in the number and sequence of ZF domains was observed. The number of ZF repeats in the ZF array varied from eight to 12 yielding five homozygous and 10 heterozygous genotypes. The number of different ZF domains was 84 and 52 producing 36 and 26 unique alleles in goats and sheep respectively. The posterior mean of dN/dS or omega values were calculated using the codeml tool of pamlx to identify amino acids that are evolving positively in goats and sheep, as positions -1, +3 and +6 in the ZF domain have been reported to experience strong positive selection across different lineages. Our study identified sites -5, -1, +3, +4 and +6 to be experiencing positive selection. Small ruminant zinc fingers were also found to be evolving under concerted evolution. Our results demonstrate the existence of a vast diversity of PRDM9 in goats and sheep, which is in concert with reports in many metazoans.

Genetic diversity estimates point to immediate efforts for conserving the endangered Tibetan sheep of India

Tibetan is a valuable Himalayan sheep breed classified as endangered. Knowledge of the level and distribution of genetic diversity in Tibetan sheep is important for designing conservation strategies for their sustainable survival and to preserve their evolutionary potential. Thus, for the first time, genetic variability in the Tibetan population was accessed with twenty five inter-simple sequence repeat markers. All the microsatellites were polymorphic and a total of 148 alleles were detected across these loci. The observed number of alleles across all the loci was more than the effective number of alleles and ranged from 3 (BM6506) to 11 (BM6526) with 5.920 ± 0.387 mean number of alleles per locus. The average observed heterozygosity was less than the expected heterozygosity. The observed and expected heterozygosity values ranged from 0.150 (BM1314) to 0.9 (OarCP20) with an overall mean of 0.473 ± 0.044 and from 0.329 (BM8125) to 0.885 (BM6526) with an overall mean 0.672 ± 0.030, respectively. The lower heterozygosity pointed towards diminished genetic diversity in the population. Thirteen microsatellite loci exhibited significant (P < 0.05) departures from the Hardy-Weinberg proportions in the population. The estimate of heterozygote deficiency varied from - 0.443 (OarCP20) to 0.668 (OarFCB128) with a mean positive value of 0.302 ± 0.057. A normal 'L' shaped distribution of mode-shift test and non-significant heterozygote excess on the basis of different models suggested absence of recent bottleneck in the existing Tibetan population. In view of the declining population of Tibetan sheep (less than 250) in the breeding tract, need of the hour is immediate scientific management of the population so as to increase the population hand in hand with retaining the founder alleles to the maximum possible extent.